| 查看: 465 | 回复: 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])  |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab中如何自己定义一个新函数，然后计算时可以调用？
已经有12人回复
- Matlab中有这两个现成的函数吗？没有的话，哪位朋友能帮忙写一写。。。急用
已经有3人回复
- 用matlab 的 fminu 函数做非线性最优化 出现错误， 求高手帮忙
已经有4人回复
- 求助：能不能帮忙写一个用matlab求一个函数的傅里叶变换的程序
已经有16人回复
- 求助会matlab的大侠帮我做下这个复杂函数求导
已经有4人回复
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
2楼2013-01-29 20:40:58
greetrunners
金虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 487.8
- 散金: 544
- 红花: 6
- 帖子: 1134
- 在线: 190.1小时
- 虫号: 1632555
- 注册: 2012-02-21
- 性别: GG
- 专业: 分离过程
3楼2013-01-30 10:43:58
