非线性方程组求解的牛顿迭代法用MATLAB实现培训讲学.docVIP

1 1. 二元函数的newton迭代法理论分析 设在点的某一邻域内连续且有直到2阶的连续偏导数，为该邻域内任意一点，则有 其中 ， 于是方程可近似表示为 即 同理，设在点的某一邻域内连续且有直到2阶的连续偏导数，为该邻域内任意一点，亦有 其中， 于是方程可近似表示为 即 于是得到方程组 求解这个方程组，当 时 从而 (1) 记符号 于是(1)式可改写为 (2) 迭代公式为: (3) 通过迭代公式(3)可以迭代出当时，的值，当 (为给定的误差控制项)时，原方程组的根即为。 newton迭代法求解给定的线性方程组 方程组 其中 求解过程如下 于是迭代公式为 为了解出正负轴的两个解，需要对函数f进行变形。 3.MATLAB编程实现过程 先画出函数图像找出大概位置 ezplot(exp(x^(-2)+y^(-2))=4,[-6,6,-6,6])%画出函数g y=2.4:0.001:3.8; x=(tan(1)+4-y.^(3/2)).^3; hold on plot(x,y)%画出函数f 将图放大观察 由图可以看出两个交点的大概位置是（-1,3.4）和（1，2.6）。 所以将这两个点作为初始值进行迭代计算，MATLAB编程如下： for m=1:2;