版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

cheerfulchen

铁杆木虫 (小有名气)

应助: 6 (幼儿园)
金币: 4934.6
红花: 2
帖子: 127
在线: 62.4小时
虫号: 946644
注册: 2010-01-22
性别: GG
专业: 分离过程

[求助] 用lsqnonlin拟合 NRTL三元方程回归参数出现错误求解脱

function NRTLBT

a0=[1 1 1 1 1 1 1 1 1 ];
x=Lsqnonlin(@csjs,a0)
function z=csjs
x1 =[0.0022 0.0358 0.0134 0.0097 0.0073 0.0035];
x2 =[0.0224 0.339 0.1227 0.0857 0.0721 0.032];
x3=1-x1-x2;
y1e = [0.0098 0.0379 0.0329 0.037 0.0316 0.0168];
y2e = [0.0098 0.0379 0.0329 0.037 0.0316 0.0168];
y3e=1-y1e-y2e;
p =[2561 3112 3454 3600 3600 3710];
ps1 = [2530.011324 3731.190276 3731.190276 3731.190276 3731.190276 3731.190276];
ps2 = [38.27422323 91.99677979 91.99677979 91.99677979 91.99677979 91.99677979];
ps3 = [19.78428868 47.08467815 47.08467815 47.08467815 47.08467815 47.08467815];
T=[333 353 353 353 353 353];
R=8.314;
for i = 1:6
%G12=exp(-a(1)*a(4)/R/T(i));tao12=a(4)/R/T(i);%a12=a(1),a13=a(2),a23=a(3),g12-g22=a(4),g13-g33=a(5),g23-g33=a(6)
%G13=exp(-a(2)*a(5)/R/T(i));tao13=a(5)/R/T(i);
%G23=exp(-a(3)*a(6)/R/T(i));tao23=a(6)/R/T(i);
a=zero[];
tao12=a(1)/R/T(i);tao21=a(2)/R/T(i);
tao13=a(3)/R/T(i);tao31=a(4)/R/T(i);
tao23=a(5)/R/T(i);tao32=a(6)/R/T(i);
G12=exp(-a(7)*a(1)/R/T(i));G21=exp(-a(7)*a(2)/R/T(i));
G13=exp(-a(8)*a(3)/R/T(i));G31=exp(-a(8)*a(4)/R/T(i));
G23=exp(-a(9)*a(5)/R/T(i));G32=exp(-a(9)*a(6)/R/T(i));
y1c(i) = x1(i)*ps1(i)/p(i)*exp((G21*x2(i)+G31*x3(i))*(tao21*G21*x2(i)+tao31*G31*x3(i))/(x1(i)+x2(i)*G21+x3(i)*G31)^2)...
+(tao12*G12*x2(i)^2+G12*G32*x2(i)*x3(i)*(tao12-tao32))/(x1(i)*G12+x2(i)+x3(i)*G32)^2+(tao13*G13*x3(i)^2+G13*G23*x2(i)*x3(i)*(tao13-tao23))/(x1(i)*G13+x2(i)*G23+x3(i))^2;
y2c(i) = x2(i)*ps2(i)/p(i)*exp((G12*x1(i)+G32*x3(i))*(tao12*G12*x1(i)+tao32*G32*x3(i))/(x1(i)*G12+x2(i)+x3(i)*G32)^2)...
+(tao23*G23*x3(i)^2+G13*G23*x1(i)*x3(i)*(tao23-tao13))/(x1(i)*G13+x2(i)*G23+x3(i))^2+(tao21*G21*x1(i)^2+G21*G31*x1(i)*x3(i)*(tao21-tao31))/(x1(i)+x2(i)*G21+x3(i)*G31)^2;
y3c(i) = x3(i)*ps3(i)/p(i)*exp((G13*x1(i)+G23*x2(i))*(tao13*G13*x1(i)+tao23*G23*x2(i))/(x1(i)*G13+x2(i)*G23+x3(i))^2)...
+(tao21*G31*x1(i)^2+G21*G31*x1(i)*x2(i)*(tao31-tao21))/(x1(i)+x2(i)*G21+x3(i)*G31)^2+(tao32*G32*x2(i)^2+G12*G32*x1(i)*x2(i)*(tao32-tao12))/(x1(i)*G12+x3(i)*G32+x2(i))^2;

z(i)=(y1e(i)-y1c(i))^2 + (y2e(i) - y2c(i))^2+(y3e(i) - y3c(i))^2;

end

%Matlab函数：lsqnonlin('NRTL',x0)

[ Last edited by cheerfulchen on 2012-4-6 at 09:43 ]

回复此楼

» 猜你喜欢

如何从铁电相到顺电相。已经有1人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有252人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有42人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有24人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
PVA溶解已经有2人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-04-05 20:57:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chaiqing123

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 679.7
红花: 1
帖子: 71
在线: 65.2小时
虫号: 1704607
注册: 2012-03-20
性别: GG
专业: 化工热力学和基础数据

目标函数错了本身就包括求和累加

赞一下