| 查看: 1348 | 回复: 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 ] |
回复此楼» 猜你喜欢
机械专硕270求调剂，接受跨专业
已经有6人回复
-
材料与化工300求调剂
已经有27人回复
-
一志愿西北工业大学289 085602
已经有10人回复
-
求调剂 材料与工程 324分 专硕
已经有6人回复
-
调剂
已经有6人回复
-
求调剂
已经有5人回复
-
一志愿华东师范生物学326分，求调剂
已经有5人回复
-
还有化工二轮调剂的学校吗
已经有40人回复
-
080100力学316求调剂
已经有5人回复
-
一志愿矿大，材料工程专硕314分，0856可调都可以
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求高手帮忙看一下我的结构式画的有没有问题？
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
上海交通大学化学化工学院膜分离-纳流控课题组诚聘博士后2-3名
+1/88
-
邀您投稿 Polymers 特刊-高分子材料加工主题
+1/87
-
085404 计算机技术 真诚求老师带走 杭电 336 六级已过
+1/85
-
科研小助理招聘~~~开工啦
+1/38
-
福建医科大学生信系饶世涛老师2026年招收调剂优秀硕士生1名-生物医学信息学
+1/37
-
喀什大学环境生态修复功能材料与技术团队接收2026年调剂研究生
+1/34
-
岭南师范学院生物与医药专业（生物技术与工程领域）2026年调剂公告
+1/21
-
广东工业大学生态环境与海洋学院 张轩教授 招收膜分离方向 2026年9月入学 博士生1人
+1/17
-
齐鲁工业大学（山东省科学院）海洋流体机械团队 2026 年博士研究生招生
+1/17
-
香港科技大学（广州）杨其纯课题组招收水文学方向博士后研究人员两名
+1/13
-
上海应用技术大学 化工学院 硕士招生 明天面试 化学、化工、材料、轻工、环境等专业
+1/12
-
做流体实验还在凑设备？这套高精密流场可视化方案真的适合高校吗？
+1/11
-
2026年西华大学材料学院-材料近净成形与表面技术研究团队-招收调剂研究生
+1/10
-
福建农林大学资环学院李香真-姚敏杰教授课题组招收硕士调剂生
+1/7
-
贵州财经大学新一代电子信息技术（含量子技术等）硕士专业（085401）接收调剂
+1/5
-
【调剂上岸】985联合培养！1007/1055/08/07等，不限专业，想联培的速来！
+1/4
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/4
-
【调剂上岸】985联合培养！1007/1055/08/07等，不限专业，想联培的速来！
+1/2
-
丽水学院材能学院-先进能源存储与器件课题组-26届研究生招生-（新能源电池方向）
+1/2
-
长江师范学院 材料工程 招收调剂学生
+1/1
3楼2010-12-07 10:40:43
2楼2010-12-07 09:18:15
