计算机在化工中的应用 结课作业

n=size(p); x2=1-x1; y2e=1-y1e;

ps=(exp([A(1)-A(2)./(T+A(3));A(4)-A(5)./(T+A(6))]))';%饱和蒸汽压 ps1=ps(:,1)'; ps2=ps(:,2)';

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% for i=1:n(2)

y1c(i)=exp(a(1)*x2(i)/(x1(i)+a(1)*x2(i)))*ps1(i)*x1(i)/...

( p(i)*(x1(i)+a(1)*x2(i))*exp((a(2)*x2(i)/(x2(i)+a(2)*x1(i))))); y2c(i)=exp(a(2)*x1(i)/(x2(i)+a(2)*x1(i)))*ps2(i)*x2(i)/...

( p(i)*(x2(i)+a(2)*x1(i))*exp((a(1)*x1(i)/(x1(i)+a(1)*x2(i)))));

z(i)=(y1e(i)-y1c(i))^2+(y2e(i)-y2c(i))^2;%以两个组分计算值与实验值之差的平方和作为目标函数

% z(i)=y1e(i)-y1c(i);%以其中一个组分的计算值与实验值之差的平方和作为目标函数 % z(i)=y2e(i)-y2c(i); end

四、模拟结果

1. 安托因常数 物质 甲醇 糠醛

2. 活度系数参数Λ12=2.8807,Λ21=-0.0184 3. 模拟得到的气相组成 T/K 434.18 380.45 358.15 351.35 347.95 347.15 345.15 343.60 341.15 339.95 337.85

yCH3OH 0 0.9540 0.9754 0.9832 0.9867 0.9922 0.9911 0.9941 0.9968 0.9982 1.0000

y1c 0 1.0780 0.9539 0.8131 0.8967 0.9637 0.9631 0.9128 0.8451 1.0341 0.7037

Er 0 -0.1240 0.0215 0.1699 0.0900 0.0285 0.0280 0.0813 0.1517 -0.0359 0.2963

A 19.2 16.6

B 1172.1 1089.9

C -191.8 -222.1

五、心得与体会

自己刚刚接触matlab有半个学期的时间，说实话我现在对MATLAB还是摸不着头脑，一方面是自己接触的时间太短，另一方面，就是自己在上机方面投入的时间有限，实践比较少。现在，我对MATLAB的印象仅仅在解决习题和绘制图形上，但是我很喜欢MATLAB的简单的语法,易于绘制图形,编程也非常容易, 并且具有功能强大的开放式的toolbox。 因此,尽管我一直没有这方面的应用,但是我还是对 它非常感兴趣，自己正打算暑假好好研究研究MATLAB。

六、参考文献

[1] 刘卫国. MATLAB程序设计教程. 北京：中国水利水电出版社，2005. [2] 陈小华.《化学工程》. 1987年第03期.