| 查看: 646 | 回复: 1 | ||
cheerfulchen铁杆木虫 (小有名气)
|
[求助]
用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 ] |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有170人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ABAQUS安装成功后，运行出现问题，请大家帮忙看一下
已经有14人回复
- 添加聚乙烯醇后 出现 团聚，求教原因
已经有10人回复
- 材料内出现纳米晶
已经有6人回复
- 辰华-恒电流计时电位法测催化剂耐久性如何设参数？
已经有7人回复
- 求高手看下这个DNA单链的结构错误何在？
已经有13人回复
- 氘代水中能出现活泼氢么？
已经有5人回复
- 这是什么错误？
已经有3人回复
- R语言图书-求解差分方程和并行计算！
已经有43人回复
- 化药6类四种介质溶出度拟合时用做方法学吗？
已经有8人回复
- 求助这个非线性方程改怎么解啊？？
已经有4人回复
- 求助：origin高斯拟合
已经有10人回复
- 拟合过过弹性模量的没？
已经有43人回复
- 求绳子的曲线方程
已经有6人回复
- 求matlab 解方程的问题
已经有5人回复
- 对数据进行曲线拟合后如何得到拟合误差
已经有30人回复
- 脉冲微分方程内涵
已经有3人回复
- 重金加人民币求拟合wilson方程和NRTL方程参数的程序
已经有17人回复
chaiqing123
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 679.7
- 红花: 1
- 帖子: 71
- 在线: 65.2小时
- 虫号: 1704607
- 注册: 2012-03-20
- 性别: GG
- 专业: 化工热力学和基础数据
2楼2012-04-26 09:01:58
