牛顿迭代公式
设待求解的方程为f(x)=0。一个根的初始猜测值记为X0,那么在X0点的切线可以表示为
公式1
这条切线与x轴的交点记为X1,则
公式2
再过X1做切线,同理得到X2,X3...
递推表达式为
公式3
由此不断更新的Xn会趋于精确解。
收敛条件:一般|Xn+1-Xn|≤0.0001,作为判断条件。(PS:有时候值本身很小,就需要处理数据,如取对数等)
应用
使用牛顿迭代法求PR方程解。
压强P的显函数,很容易求解。
求摩尔体积V工程上最常用的情况,因为P,T易测。用迭代法求解。
步骤
构造函数
求导:
迭代公式:
迭代起点:
迭代终点:
用EXCEL编写PR方程计算表,结果如下:
在ASPEN软件中使用PR方程结果如下:
-END-
扫码关注我们
查看更多精彩内容