| 查看: 1614 | 回复: 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 ] |
回复此楼» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
5楼2009-08-03 20:28:16
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
