| 查看: 1617 | 回复: 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 ] |
回复此楼» 猜你喜欢
一志愿华东师范大学有机化学专业，初试351分，复试被刷求调剂!
已经有4人回复
-
298求调剂
已经有4人回复
-
材料学硕333求调剂
已经有11人回复
-
0856求调剂
已经有7人回复
-
356求调剂
已经有3人回复
-
332求调剂
已经有7人回复
-
348求调剂
已经有4人回复
-
305求调剂
已经有8人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有13人回复
-
调剂310
已经有9人回复
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
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
5楼2009-08-03 20:28:16
spc08
荣誉版主 (文学泰斗)
- 应助: 11 (小学生)
- 贵宾: 4.043
- 金币: 230245
- 散金: 955
- 红花: 57
- 沙发: 11
- 帖子: 181659
- 在线: 2753.4小时
- 虫号: 712283
- 注册: 2009-03-01
- 专业: 无机材料化学
- 管辖: 仿真模拟
★ ★ ★ ★ ★
cooler8395(金币+5,VIP+0): 8-25 21:23
cooler8395(金币+5,VIP+0): 8-25 21:23
|
楼主弄成附件，想帮忙的人还得花金币下载，这次还得我来替楼主贴出来，下次不管了啊， 楼主的程序： % 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)%最大偏差 |
6楼2009-08-03 20:28:49
7楼2009-08-03 20:40:45
