(无奖求助)二元非线性方程求根问题

944f4fe875fb39d

那么，你原来给出的 1.0082 和 1.6719 是什么？

gamma_1 = 1.031

gamma_2 = 1.667

把这两个数据代进去试试，看看对不对？

醉乡常客

1.0082和1.6719是我自己用Antoine方程算出的活度系数。

教材上给出的1.013和1.667代回去不对。

Antoine方程（压力单位是bar）算得的64.3℃饱和蒸汽压
甲醇 1.004817bar
甲基乙基酮 0.609528bar

我觉得书上给出的0.982×10^5、0.6067×10^5两个数有问题，此后的活度系数也就有点问题。
（这是陈钟秀书上的例7-3，我帖的图的上半截就是书上的解题过程。）

944f4fe875fb39d

按1.013和1.667，得到 Lambda12=1.08182135，Lambda21=0.377766147

是不是 1.0818、0.3778？用它们计算 t 和 y1 结果是否更好些？

另外，你再查一下安东尼常数给得对不对。

有时候，人们在处理实验数据时用一套数据得到了较好的结果，写书时就忘了。

醉乡常客

浙大陈钟秀教材上的甲醇和甲乙酮数据和例题一致。

天大马沛生教材上Antoine方程是分段函数，且无甲乙酮数据。

louselice^_^

醉兄，我从移动硬盘里找到软件了，需要哪个，我给你传^_^

944f4fe875fb39d

醉兄，我试算了你的题目，几乎可以肯定，题目给出的数据不是他原来计算时用的数据。

他原来的计算结果是合理的。

我取 gamma_1 = 1.031，gamma_2 = 1.667，得到 Lambda12 = 1.0818，Lambda21 = 0.3778，用此数据和你所给的安东尼常数计算得到：

x1  y1  y1*  t (°C)  t* (°C)
0.076  0.193  0.188  75.3  75.5
0.197  0.377  0.381  70.7  71.1
0.356  0.528  0.541  67.5  67.6
0.498  0.622  0.641  65.9  65.7
0.622  0.695  0.714  65.1  64.6
0.747  0.777  0.784  64.4  64.0
0.844  0.832  0.846  64.3  63.8
0.936  0.926  0.922  64.4  64.0

*号为计算值。

我觉得问题在安东尼常数和恒沸点的x1的不一致上，按该组常数和恒沸点的x1，使得Lambda12 和 Lambda21 无正确的解。

按最可能犯的错误判断，他们当初计算该例题时用的是另外的安东尼常数，或者是单位搞错了，或者是单位转换时什么地方搞混淆了……

醉乡常客

x1	T	y1	y1 实验值	温度 实验值
0	352.75	0.00000
0.05	349.5	0.14303
0.076	348.12	0.20220	0.193	348.45
0.1	346.94	0.24963
0.197	343.37	0.39687	0.337	343.85
0.2	343.3	0.40086
0.3	340.83	0.50452
0.356	339.78	0.54999	0.528	340.65
0.4	339.1	0.58165
0.498	337.95	0.64329	0.622	339.05
0.5	337.93	0.64442
0.6	337.13	0.69859
0.622	337	0.71004	0.695	338.25
0.7	336.65	0.75052
0.747	336.49	0.77484	0.777	337.55
0.8	336.42	0.80538
0.829	336.41	0.82330	0.932	337.45
0.842	336.43	0.83230
0.9	336.59	0.87696
0.936	336.8	0.91129	0.926	337.55
0.95	336.95	0.92798
1	337.65	0.99978

醉乡常客

0.161737  1.26246

x1	T	y1	y1实验值	温度实验值
0	352.75	0
0.05	346.27	0.222321
0.076	344.73	0.273694	0.193	348.45
0.1	343.73	0.30786
0.197	341.39	0.395499	0.337	343.85
0.2	341.34	0.397617
0.3	340	0.46128
0.356	339.44	0.494634	0.528	340.65
0.4	339.07	0.520999
0.498	338.39	0.581303	0.622	339.05
0.5	338.38	0.582616
0.6	337.88	0.648874
0.622	337.8	0.664437	0.695	338.25
0.7	337.57	0.722145
0.747	337.49	0.75959	0.777	337.55
0.8	337.43	0.803842
0.829	337.42	0.829335	0.932	337.45
0.842	337.42	0.84105
0.9	337.45	0.895483
0.936	337.5	0.931404	0.926	337.55
0.95	337.55	0.946743
1	337.65	0.999775

ereree

不要用matlab那个巨无霸，用Scilab吧免费的。

brook

难得一见的理工讨论贴，可为论坛理工讨论的模板，经典！

不知醉版后事如何？

醉乡常客

引用第49楼brook于2007-09-22 23:31发表的 :
难得一见的理工讨论贴，可为论坛理工讨论的模板，经典！

不知醉版后事如何？

我的结论和bookish老师说的一样。该书此题未验算！

当然，我不敢肯定。毕竟这本书算得上是一本非常优秀的热力学教材。

ly188

额 做了个验算的工具
怎么解我真的忘了 - -+
这里高手真多