非线性方程组数值上实机验报告.docxVIP

非线性方程组数值上机实验报告练习6.23. 用Newton 方法求解方程 的根，并绘制误差下降曲线。然后，试用割线法重复上述工作。实验步骤：1.编写解非线性方程的Newton迭代法和割线法程序。2.作出误差下降曲线。实验图像：实验图像分析：通过简单计算可知,方程 的三个根为：2（二重根）、-3。1.牛顿法在单根处为二阶收敛，在重根处仅为一阶收敛。取初值为1时，收敛到二重根2，取初值为-2时，收敛到单根-3.从上图中Newton法解方程的误差下降曲线,可以看出取初值为-2时，收敛到-3的速率远大于取初值为1时，收敛到2的速率。这与理论分析相符。2.割线法的启动需要两个初值，此处是由Newton法迭代一步算出另一初值后再进行割线法。从上图亦可以看出：割线法同样对单根和重根的收敛速度有差别。对单根的收敛速度较快，对重根的收敛速度较慢。练习6.24编制Newton方法和Broyden 方法求解非线性方程组取相同的初值，比较两者的迭代次数；若初值为呢？观察Broyden中的是否收敛到Jacobi矩阵呢？若非线性方程组为呢？初值取为，观察 Broyden 中的是否收敛到Jacobi矩阵呢？实验数据：1.取初值为（-0.15,1.4），设置相同的误差容限为，两者的迭代次数分别为：5和7。迭代7次后Broyden算法中：真解的Jacobi矩阵为，故可以认为收敛到。下图显示了两种方法的误差下降曲线。2.当取初值为时，同样设置相同的误差容限为,此处由于即为其真解，牛顿法需迭代一次，而Broyden法需迭代两次，即可达到所需量级。此时,可认为其收敛到Jacobi矩阵。3.对于新方程组，其有两组解和，取初值为时，收敛到解，牛顿法需要迭代7次，Broyden法需要迭代9次。此时Broyden法中矩阵而真解处对应的Jacobi矩阵为，注意到矩阵的和位置元素与0和6相差较远，故可认为并不收敛到Jacobi矩阵。所得误差下降曲线如下：实验数据分析：1.从上面两图中可以看出，Broyden 方法的收敛速度稍慢于Newton方法。这与理论相符。2.从上面结果可以看出，Broyden方法中也不一定收敛。事实上，我们只需要拟牛顿方法中的修正方向和原始Newton迭代中的修正方向趋于一致，即可保证数值解趋于真解，无需收敛于对应的Jacobi矩阵。练习6.25 编制修正Newton法（与m的关系 ）、离散Newton法和两点序列割线法求解上述非线性方程组。实验数据与分析：对于上述非线性方程组，其有两组解和.①修正的Newton法：为了充分体现修正的Newton法与m的关系，在迭代时取较大的误差容限为,取初值为，取,恰好均收敛到解,得到相应的迭代次数如下：M123456789迭代次数：765444333由于修正的Newton法在一定条件下至少为m+1阶收敛，m=1,则对应2阶收敛，依次类推。从上表中可以看出，m=1时迭代7次，相当于m=7时迭代3次下图为误差下降曲线图：此处为了近似地每一步的收敛速度，图像中横坐标步长与m对应，对应m=1:4,相当于每一步迭代了m次。从图中可以看出，就收敛速度而言，m=1到m=4依次递减。这是合理的，因为m=1,相当于每一步都要形成Jacobi矩阵，而m=4 相当于四步形成一次Jacobi矩阵，故收敛较慢。但在实际应用中，为了在不增加求逆次数的情况下提高收敛速度，利用修正的Newton法更好。②离散Newton法（单点牛顿法）：其思想为，其中为的一种差商逼近。这里取（其中 来近似地代替 。③两点序列割线法：其迭代格式为：，这里取教材上的第二种两点序列割线法，即用差商，来近似地代替对于离散Newton法和两点序列割线法，均取初始值为,设置误差容限为,得到图像如下：从上图可以看出，两点序列割线法稍快于离散牛顿法，这与理论分析相符。练习6.26.非线性方程组的应用：Newton法可应用于特征值问题的求解，方式如下：取T为以前的三对角矩阵，阶数为；取任意的初值和，做Newton迭代。将得到的结果与幂法作比较。实验数据与分析：1.当取初始向量为任意值时，经试验可通过取不同的初值来使得上述Newton算法给出的全部特征值。2.为了便于与乘幂法比较，取初始向量为，此时Newton方法恰好收敛到 的最大特征值。对乘幂法和上述Newton法设置相同的相邻误差限为,实验结果为：乘幂法迭代44次终止，而如上所述的Newton方法仅迭代15次便终止。下图为两者的误差下降曲线图。从下图中可以看出Newton法的收敛速度远大于乘幂法。综合比较乘幂法和如上的牛顿方法：1.乘幂法只可以求出第一主特征值，而上述的Newton方法可以通过取不同初值求出的所有特征值。2.乘幂法的收敛速度依赖于第二主特征值与第一主特征值的比值大小，而Newton法的收敛速度在单根时为二阶局部收敛