| 查看: 1347 | 回复: 2 | |||
[交流]
【求助】求高手帮看一下程序（泡点温度计算）
|
|
使用SRK方程和wilson方程计算泡点温度和组成 m文件中，y（1）是组分1 的气相组成，y（2）是温度 m文件如下 function f=myfun(y,x1) p=101325; T1c=579.35;V1c=219E-06;Z1c=0.259;p1c=56.9E+05;w1=0.193; ANTA1=36.6016;ANTB1=-2979.4;ANTC1=-10.104;ANTD1=1.1445E-09;ANTE1=3.2472E-06; T2c=591.79;V2c=315.8e-06;Z2c=0.264;p2c=41.09e+05;w2=0.264; ANTA2=34.0755;ANTB2=-3.0379e+03;ANTC2=-9.1635;ANTD2=1.0289E-11;ANTE2=2.7035E-06; b=[4198.6 -3001.6]; x2=1-x1; p1s=10^(ANTA1+ANTB1/y(2)+ANTC1*log10(y(2))+ANTD1*y(2)+ANTE1*y(2)* y(2))*133.322; p2s=10^(ANTA2+ANTB2/y(2)+ANTC2*log10(y(2))+ANTD2*y(2)+ANTE2*y(2)* y(2))*133.322; SRKa1=0.42748*8.31441*8.31441*T1c*T1c/p1c*(1+(0.480+1.574*w1-0.176*w1^2)*(1-(y(2)/T1c)^0.5))^2; SRKb1=0.08664*8.31441*T1c/p1c; SRKk1=-8.31441*y(2)/p1s; SRKm1=SRKa1/p1s-SRKb1*8.31441*y(2)/p1s-(SRKb1)^2; SRKn1=-SRKa1*SRKb1/p1s; SRKp1=[1 SRKk1 SRKm1 SRKn1]; SRKV1=roots(SRKp1); V1=max(SRKV1); Z1=p1s*V1/8.31441/y(2); fai1s=exp((Z1-1)-log(p1s*(V1-SRKb1)/8.31441/y(2))-SRKa1*log(1+SRKb1/SRKV1)/SRKb1/8.31441/y(2)); SRKa2=0.42748*8.31441*8.31441*T2c*T2c/p2c*(1+(0.480+1.574*w2-0.176*w2^2)*(1-(y(2)/T2c)^0.5))^2; SRKb2=0.08664*8.31441*T2c/p2c; SRKk2=-8.31441*y(2)/p2s; SRKm2=SRKa2/p2s-SRKb2*8.31441*y(2)/p2s-(SRKb2)^2; SRKn2=-SRKa2*SRKb2/p2s; SRKp2=[1 SRKk2 SRKm2 SRKn2]; SRKV2=roots(SRKp2); V2=max(SRKV2); Z2=p2s*V2/8.31441/y(2); fai2s=exp((Z2-1)-log(p2s*(V2-SRKb2)/8.31441/y(2))-SRKa2*log(1+SRKb2/V2)/SRKb2/8.31441/y(2)); Tc12=(T1c*T2c)^0.5; Vc12=((nthroot(V1c,3)+nthroot(V2c,3))/2)^3; Zc12=(Z1c+Z2c)/2; pc12=Zc12*8.31441*Tc12/Vc12; w12=(w1+w2)/2; SRKa12=(SRKa1*SRKa2)^0.5; SRKam=y(1)*y(1)*SRKa1+2*y(1)*(1-y(1))*SRKa12+(1-y(1))*(1-y(1))*SRKa2; SRKbm=y(1)*SRKb1+(1-y(1))*SRKb2; SRKkm=-8.31441*y(2)/p; SRKmm=SRKam/p-SRKbm*8.31441*y(2)/p-(SRKbm)^2; SRKnm=-SRKam*SRKbm/p; SRKpm=[1 SRKkm SRKmm SRKnm]; SRKVm=roots(SRKpm); Vm=max(SRKVm); Zm=p*Vm/8.31441/y(2); fai1=exp(SRKb1/SRKbm*(Zm-1)-log(p*(Vm-SRKbm)/8.31441/y(2))+SRKam*(SRKb1/SRKbm-2*(y(1)*SRKa1+(1-y(1))*SRKa12)/SRKam)*log(1+SRKbm/SRKVm)/SRKbm/8.31441/y(2)); fai2=exp(SRKb2/SRKbm*(Zm-1)-log(p*(Vm-SRKbm)/8.31441/y(2))+SRKam*(SRKb2/SRKbm-2*(y(1)*SRKa12+(1-y(1))*SRKa2)/SRKam)*log(1+SRKbm/SRKVm)/SRKbm/8.31441/y(2)); V1L=V1c*Z1c^((1-y(2)/T1c)^(2/7)); V2L=V2c*Z2c^((1-y(2)/T2c)^(2/7)); A12=(V2L/V1L)*exp(-b(1)/8.314/y(2)); A21=(V1L/V2L)*exp(-b(2)/8.314/y(2)); r1cal=exp(-log(x1+A12*x2)+x2*(A12/(x1+A12*x2)-A21/(x2+A21*x1))); r2cal=exp(-log(x2+A21*x1)-x1*(A12/(x1+A12*x2)-A21/(x2+A21*x1))); f=abs(y(1)*p*fai1-r1cal*p1s*fai1s*x1*exp(V1L*(p-p1s)/8.31441/y(2)))+abs((1-y(1))*p*fai2-r2cal*p2s*fai2s*x2*exp(V2L*(p-p2s)/8.31441/y(2))); main文件如下 x1=[0.9643 0.9165 0.8322 0.7520 0.6696 0.6362 0.5576 0.4416 0.3864 0.3407 0.2988 0.2628 0.2271 0.1877 0.1482 0.0933 0.0586 0.0375]; n=18; for i=1:n x0=[1;355]; [x fval]=fminsearch(@(x) myfun(x,x1(i)),x0) x1(i)=x(1); x2(i)=x(2); end; x1; x2; 运行提示Subscripted assignment dimension mismatch. Error in ==> fminsearch at 205 fv(:,1) = funfcn(x,varargin{:}); [ Last edited by flying_fay on 2010-12-1 at 11:26 ] |
回复此楼» 猜你喜欢
288求调剂
已经有14人回复
-
一志愿厦大生物学332求调剂
已经有7人回复
-
085600，321分求调剂
已经有17人回复
-
一志愿鲁东大学071000生物学学硕初试分数276求调剂
已经有16人回复
-
材料化工总分334求调剂
已经有11人回复
-
288求调剂，一志愿华南理工大学071005
已经有13人回复
-
复试调剂，一志愿郑州大学材料与化工289分
已经有20人回复
-
一志愿中科大070300化学，314分求调剂
已经有4人回复
-
285求调剂
已经有9人回复
-
材料调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求高手帮忙看一下我的结构式画的有没有问题？
已经有16人回复
- 求matlab编程代码，求帮助，谢谢！
已经有11人回复
- 跪求高手帮忙看看这组数据如何进行拟合
已经有7人回复
- 很简单的Fortran赋值出错啦，求大神们帮忙找bug...
已经有11人回复
- 求高手帮忙看一段英文，涉及到数据分析问题，如能给出说得通的全文翻译另有50BB相赠
已经有5人回复
- 高手 帮忙看看这两数据
已经有5人回复
- 求高手分析一下DSC 曲线
已经有10人回复
- 请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
已经有9人回复
- 麻烦高手帮忙解一下三元二次方程组！谢谢谢谢！
已经有4人回复
- 高手帮我看看，这种示意图用什么软件画比较方便？
已经有7人回复
- matlab高手帮忙做个自定义复杂函数拟合，含有复数，
已经有25人回复
- 硫酸铵沉淀的一些问题，求高手帮忙看一下！
已经有12人回复
- 求SWMM软件高手指导
已经有15人回复
- 帮忙计算半衰期
已经有13人回复
- 急！！！！液质图谱超诡异~~~不知道是不是被污染了 ~~~求高手帮忙看一下
已经有21人回复
- 求有机合成高人帮在下看一下这个嘧啶成环反应能发生吗？？？？？？？
已经有18人回复
- 求高人帮忙看一下我这个Fortran程序的错误出在哪里？
已经有12人回复
- 求助 哪位傅里叶红外高人帮忙看一下我的红外图谱
已经有18人回复
- 请各位高人看一下这够不够JBC的要求？
已经有4人回复
- 运行一个python程序，出现问题，请高手赐教
已经有9人回复
- 刚投就被拒稿，请高手帮忙分析一下
已经有18人回复
- LABVIEW四阶龙格库塔法求解一阶微分方程组
已经有12人回复
» 抢金币啦！回帖就可以得到:
-
浙江师范大学数理医学招聘专任教师
+5/325
-
★★爆★★0860生物与医药调剂考生看过来，河北大学化学与材料科学学院
+1/218
-
课题组(化学、化工、环境、能源、材料相关方向)欢迎你的加入
+1/187
-
★★精品★★0860生物与医药还有调剂指标，最后的机会，捡漏上岸
+1/132
-
香港科技大学(广州) CCUS及吸附分离方向全奖博士生，年津贴18万
+1/91
-
烟台大学环境与材料工程学院先进储能材料与器件课题组招收材料、化学、环境方向硕士
+1/87
-
B区学生调剂-兰州交通大学材料科学与工程学院
+1/45
-
哈尔滨工业大学化工学院中俄联合办学硕士生招生
+1/21
-
武汉工程大学化学与环境工程学院2026年硕士研究生调剂公告
+1/13
-
山东大学谢玉生课题组-诚聘博士后/科研助理（长期有效）
+1/8
-
中国农业大学资源与环境学院 - 招聘杰出人才（教授-专技四级）、青年研究员（副教授-�
+1/6
-
【#上海调剂急录#】能接受985联合培养的速来！带你发一区文章！
+1/6
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/6
-
0854,0858调剂招生
+1/5
-
贵州财经大学新一代电子信息技术（含量子技术等）硕士专业（085401）接收调剂
+1/5
-
找好工作来我这-广东唯一石化院校——资源与环境专硕招生
+1/5
-
哈尔滨医科大学李老师课题组招收生物信息学方向博士、硕士研究生
+1/4
-
军事医学研究员毒物药物研究所药理学生物技术药物代谢方向招收调剂生 4名
+1/2
-
丽水学院材能学院-先进能源存储与器件课题组-26届研究生招生-（新能源电池方向）
+1/2
-
一志愿211电子信息347求调剂
+1/1
2楼2010-12-07 09:18:15
3楼2010-12-07 10:40:43
