本帖最后由 胜源毅 于 2013-2-6 11:04 编辑
刚刚读了论坛里分享的化工流程模拟实训第十三章收敛部分,书上提到“Aspenplus 中有7中收敛方法: “Wegstein、Direct、Secant、Broyden、Newton、Complex、SQP”。这不是数据处理方法吗? 翻阅了一下相关资料,自己尝试了一下,发现有些方法不太好收敛吧。
先试一下Wegstein,比如公式 'x=-0.9x^2+1.7x+2.5',初始值为1和5
结果如下所示
收敛的不错,但尝试了一下 x=-x^2,就不太好收敛了。
Direct方法
Aspen官方用户指南上说”对于直接迭代 撕裂流股变量的新值是由以前的流程计算循环公式导出来的值,X_k+1 = G(X_k)“ ,然后说”使用直接迭代这种方法 收敛较慢但肯定会收敛“。这个比较纳闷,就像上面那个公式 'x=-0.9x^2+1.7x+2.5' ,随便带入个1进去,值岂不是越来越小,最后小到无穷,怎么收敛。
Secant方法和牛顿法倒是挺容易收敛的,牛顿法f(x)积分不容易求,secant可以转换
只知道这些方法的数据处理,应用到Aspen就不太懂了{:1106_371:}
祝大家新年快乐
|