| 查看: 1803 | 回复: 6 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】MATLAB回归wilson方程参数
|
|||
|
求助高手，回归wilson方程参数的matlab程序。本人编的误差比较大，不知哪里有问题。 % wilson 回归参数 调用格式-----[x,fval]=fminsearch(@wilson,[100 100])单纯形法求最小值 function f = wilson(x) A12=x(1); A21=x(2); YY=0; YYY=0; FFF=0; X1=[0.0593 0.1306 0.2660 0.3842 0.4708 0.6256 0.7178 0.8214 0.9037 0.9397]; Y1=[0.3285 0.4902 0.6234 0.6831 0.7125 0.7524 0.7832 0.8349 0.8864 0.9215]; for i=1:length(X1) i x1=X1(i);%提取本次试验点的实验值 x2=1-x1; % x1,x2 试验值 y1=Y1(i); y2=1-y1; % y1,y2 试验值 A12; A21; Lnr1=-log(x1+A12*x2)+x2*(A12/(x1+A12*x2)-A21/(x2+A21*x1)); Lnr2=-log(x2+A21*x1)+x1*(A21/(x2+A21*x1)-A12/(x1+A12*x2)); r1=exp(Lnr1); r2=exp(Lnr2); A1=4.1973;A2=3.4337; B1=1575.0;B2=1413.0; C1=-34.29;C2=-44.25; t=70; T=t+273.15; P10=exp(A1-B1/(C1+t)); P20=exp(A2-B2/(C2+t)); P=x1*r1*P10+x2*r2*P20 % P-总压；P10-物质1的饱和蒸汽压；P20-物质2的饱和蒸汽压. 单位：Kpa y1cal=x1*r1*P10/P y2cal=x2*r2*P20/P %------------------------------ Y=abs(y1-y1cal); %Y-本次计算点的绝对误差 YY=YY+Y; %YY-已计算点的绝对误差之和 %------------------------------- MM=abs((y1-y1cal)./y1); YYY=YYY+MM; %------------------ m=abs(y1cal-y1); M(1,i)=m; %----------------------- FF=(y1cal-y1).^2+(y2cal-y2).^2; FFF=FFF+FF;%目标函数 %a=log(r1/r2); %b=log(r1)/x2^2; % c=-log(r2)/x1^2; end f=FFF; fata=YY/i%平均偏差 sigma=YYY/i%平均相对偏差 sigmaM=max(M)%最大偏差 [ Last edited by cooler8395 on 2009-8-3 at 20:43 ] |
回复此楼» 猜你喜欢
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
40+中年人失业后，还能找到好工作吗？
已经有7人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
7楼2009-08-03 20:40:45
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
2楼2009-08-03 18:20:56
3楼2009-08-03 18:42:21
4楼2009-08-03 20:24:01
200