matlab与化工计算,matlab在化工计算中的迭代运算

1.jpg (63.03 KB, 下载次数: 7)

2010-11-13 13:16 上传

2.jpg (66.42 KB, 下载次数: 13)

2010-11-13 13:16 上传

3.jpg (54.06 KB, 下载次数: 3)

2010-11-13 13:16 上传

4.jpg (68.78 KB, 下载次数: 3)

2010-11-13 13:16 上传

5.jpg (62.6 KB, 下载次数: 1)

2010-11-13 13:16 上传

6.jpg (55.83 KB, 下载次数: 5)

2010-11-13 13:16 上传

7.jpg (56.3 KB, 下载次数: 5)

2010-11-13 13:16 上传

8.jpg (53.35 KB, 下载次数: 3)

2010-11-13 13:16 上传

9.jpg (31.05 KB, 下载次数: 5)

2010-11-13 13:16 上传

第一问的用matlab编程过程，已经出来正确结果，这个应该没什么问题。

%-- 10-11-12 上午9:45 --%

t=70                                                 %温度

y=[0.5 0.3 0.2]                                 %气相各组分的摩尔分数

A=[10.0311 9.6672 14.2043]

B=[2940.46 2945.47 2972.64]

C=[237.22 224 209]

P=exp(A-(B./(t+C)))                         % 各组分的饱和蒸汽压

p=1/sum(y./P)                                 % 总压

x=p*y./P                                          % 液相各组分的摩尔分数

第二问的用matlab编程过程，从设取温度初值到第一次迭代结束，程序都已经编好，也正确，但是从第二次迭代开始就不知道该如何继续往下编程，

请大家一起集思广益，帮忙解决一下问题。谢谢了！下面是我编了一部分的程序。

%-- 10-11-12上午9:45 --%p=80                                       %溶液总压强

x=[0.3 0.45 0.25]                     %液相各组分的摩尔分数

A=[10.0311 9.6672 14.2043]

B=[2940.46 2945.47 2972.64]

C=[237.22 224 209]t=(B./(A-log(p)))-C                  %假设各组分的饱和温度t0=sum(x.*t)                           %温度初值

%取k=3P=exp(A-(B./(t0+C)))             %各组分的饱和蒸汽压P30=2550.9                           %组分3的饱和蒸汽压P3S=p./(sum(x.*P/P30))         %任选组分3的饱和蒸汽压t3=(2972.64./(14.2043-log(P3S)))-209      %假设各组分的饱和温度

abs(t3-t0)

[本帖最后由 李佩珍 于 2010-11-13 13:27 编辑]