| 查看: 534 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
greetrunners金虫 (著名写手)
|
[求助]
这是我编的目标函数MATLAB代码，总是出现错误，请大侠帮忙。
|
|
|
这是我编的目标函数MATLAB代码，总是出现错误，请大侠帮忙。 function obj=lsqfun1(t12,t21,t13,t31,t23,t32) %t12,t21,t13,t31,t23,t32，是所需要关联的NRTL方程参数； %G12，G21，G13，G31，G23，G32，也是NRTL方程参数； %a1,a2,a3,是组分1,2,3，液相平衡数据； %y1,y2,y3,是组分1,2,3，气相平衡数据； %ps1，ps2，ps3，是组分1,2,3的饱和蒸汽压； %p==101.35; p=101.35; G12=exp(-0.3*t12); G21=exp(-0.3*t21);%G21是方程参数 G13=exp(-0.3*t13); G31=exp(-0.3*t31); G23=exp(-0.3*t23); G32=exp(-0.3*t32); k1=(a2.^2*t21*G21^2+a3.^2*t31*G31^2+a2.*a3*G21*G31*(t21+t31))./((a1+a2.*G21+a3.*G31).^2); k2=(a2.^2*t12*G12+a2.*a3*G12*G32*(t12-t32))./((a1*G12+a2+a3*G32).^2); k3=(a3.^2*t13*G13+a2.*a3*G13*G23*(t13-t23))./((a1*G13+a2*G23+a3).^2); k4=(a1.^2*t21*G21+a1.*a3*G21*G31*(t21-t31))./((a1+a2*G21+a3*G31).^2); k5=(a1.^2*t12*G12^2+a3.^2*t32*G32^2+a1.*a3*G12*G32*(t12+t32))./((a1*G12+a2+a3*G32).^2); k6=(a3.^2*t23*G23+a1.*a3*G13*G23*(t23-t13))./((a1*G13+a2*G23+a3).^2); k7=(a1.^2*t31*G31+a1.*a2*G21*G31*(t31-t21))./((a1+a2*G21+a3*G31).^2); k8=(a2.^2*t32*G32^2+a1.*a2*G12*G32*(t32-t12))./((a1*G12+a2+a3*G32).^2); k9=(a1.^2*t13*G13^2+a2.^2*t23*G23^2+a1.*a2*G13*G23*(t13+t23))./((a1*G13+a2*G23+a3).^2); gama1=exp(k1+k2+k3); gama2=exp(k4+k5+k6); gama3=exp(k7+k8+k9); y1c=ps1.*gama1.*a1./p; y2c=ps2*gama2*a2./p; y3c=ps3*gama3*a3./p; obj=sum(sum((y1c-y1).^2+(y2c-y2).^2+(y3c-y3).^2)); >> y1=[0.8907 0.8924 0.733 0.7731 0.698 0.7778 0.7622 0.7796 0.7934]; >> y2=[0.1082 0.1067 0.6626 0.5623 0.732 0.4326 0.3978 0.3766 0.3473]; >> y3=[0 0 0 0 0 0 0 0 0]; >> a1=[0.8918 0.8933 0.3374 0.4377 0.268 0.5674 0.6022 0.6234 0.6527]; >> a2=[0.10723 0.10649 0.62613 0.53391 0.70928 0.42518 0.39054 0.37007 0.34364]; >> a3=[0.00097 0.00021 0.03647 0.02839 0.02272 0.00742 0.00726 0.00653 0.00366]; [t12 t21 t13 t31 t23 t32]=lsqnonlin(@lsqfun1,[1 1 1 1 1 1])  |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab中如何自己定义一个新函数，然后计算时可以调用？
已经有12人回复
- Matlab中有这两个现成的函数吗？没有的话，哪位朋友能帮忙写一写。。。急用
已经有3人回复
- 用matlab 的 fminu 函数做非线性最优化 出现错误， 求高手帮忙
已经有4人回复
- 求助：能不能帮忙写一个用matlab求一个函数的傅里叶变换的程序
已经有16人回复
- 求助会matlab的大侠帮我做下这个复杂函数求导
已经有4人回复
greetrunners
金虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 487.8
- 散金: 544
- 红花: 6
- 帖子: 1134
- 在线: 190.1小时
- 虫号: 1632555
- 注册: 2012-02-21
- 性别: GG
- 专业: 分离过程
3楼2013-01-30 10:43:58
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
2楼2013-01-29 20:40:58
20