于是过冷水就去查了一下解非线性方程组的算法,觉得Newton-Raphson method算法针对我们的问题比较合适,本期过冷水就给大家讲讲该算法思路
已知方程f(x)=0有近似根xk将函数f(x)在xk...这是个线性方程,记其根为xk+1,则xk+1的计算公式为:
?
这就是解一元非线性方程的牛顿迭代法公式,我们的问题是非线性方程组,需要把一元扩展到二元。...,对优化算法有了解的读者就知道初始值对优化算法影响是很大的,针对上述特定问题过冷水就想出了一种特殊判断初始值的方法。...x轴为B12,y轴为B21,Z轴是f(g21,g21,T),在f(g21,g21,T)这个曲面上找出满足所有f(g21,g21,T1)=0的[g21,g21] 可知其为一条线。用等高线1表示。...然后再找出满足所有f(g21,g21,T1)=0的[g21,g21],可知其为另外一条等高线2线。两条两条线的交点就是该方程组的解。如图。
?
图像代码如下: