化工数值计算-5探究.ppt

? 为例。设方程组解的一组初始近似值为x0,y0,f1和f2对x和y的一阶偏导数存在且连续,则将它们在近似值附近泰勒展开并忽略高于一次项,可得 ?(5-26) 令 ? 这样,非线性方程的线性化近似后就得到一新的线性方程组: ? 若方程组系数矩阵 ? 非奇异,则可解出Δx和Δy,即得 ? 以此作为非线性方程组解的新近似值,重复上述过程,就可得到一系列的近似解,直到满足下列条件终止迭代过程。最后的迭代值可作为待求的解。 ? 更普遍的是,设X(0)是非线性方程组(5-23)的近似解,且F(X)在解的附近一阶可微,则利用泰勒展开可得 ? 即 ? 式中雅可比矩阵 ? 若矩阵DF(X)非奇异,则方程组(5-32)有唯一解,由此可得: ? 这就是解非线性方程组的牛顿-拉夫森迭代公式。显然,上式是一元非线性方程牛顿公式的直接推广,其中雅可比矩阵的逆矩阵?为迭代矩阵。牛顿-拉夫森法具有良 好的收敛性,缺点是要在每一步中计算雅可比矩阵的逆矩阵。 实际计算时,可将牛顿-拉夫森迭代公式(5-34)改写为 ? 5.3.2　程序框图 图5-6为牛顿-拉夫森法解非线性方程组的通用计算程序框图,其中线性方程组DF(X(k)) ΔX(k) = -F(X(k))的求解采用列主元消去法,参见图5-2和图5-3。 程序框图中主要变量的含义如下: N　　　变量的个数 E　　　允许误差 X　　　一维数组,