牛顿迭代法和雅克比迭代法和高斯赛德尔迭代逐次超松弛迭代.docVIP

实验三 非线性方程的牛顿法和线性方程组的迭代数值求解 信息与计算科学金融 崔振威 201002034031 实验目的： 设计牛顿迭代算法和线性方程组的常用迭代算法 实验内容：p69.3、p129.1 实验要求： 根据题目要求构造迭代格式 对线性方程组的迭代求解要求构造三种迭代格式 试比较牛顿迭代法和不动点迭代法的优劣（p69.3） 试比较线性方程组三种迭代的优劣（收敛和收敛速度） 主程序： a、牛顿迭代法程序： function [x,k]=newton(f,p0,e) %求f(x)=0在给定p0的根。 %f为所求的函数f(x)，p0为迭代初始值，e为迭代精度。k为迭代次数 %diff(f)为对函数求导，subs是赋值函数，用数值替代符号变量替换函数例 利用牛顿迭代法，在matlab窗口输入下面命令，得出结果： format long syms x f=(x^2-8) f = x^2-8 [x,k]=newton(f,3,10^-10) x = 2.82842712474619 k = 4 由输出结果可以看出：该方程迭代次数为4，值趋近2.82842712474619 (b)p0=10，,求的近似值 解：令f(x)=x^2-91, 在matlab的命令窗口中输入命令： fplot([x^2-91],[-15,15]);grid 可得下图： 利用牛顿迭代法，在matlab窗口输入下面命令，得出结果： format long syms x f=(x^2-91) f = x^2-91 [x,k]=newton(f,10,10^-10) x = 9.53939201416946 k = 4 由输出结果可以看出：该方程迭代次数为4，值趋近9.53939201416946 (c)、p0=-3，,求-的近似值 解：令f(x)=x^2-8, 在matlab的命令窗口中输入命令： fplot([x^2-8],[-3,3]);grid 可得下图： 利用牛顿迭代法，在matlab窗口输入下面命令，得出结果： format long syms x f=(x^2-8) f = x^2-8 [x,k]=newton(f,-3,10^-10) x = -2.82842712474619 k = 4 由输出结果可以看出：该方程迭代次数为4，值趋近-2.82842712474619 P129.1 1、4x-y=15 x+5y=9 解：（1）、利用雅克比迭代法，在matlab窗口输入下面命令，得出结果： format long A=[4 -1;1 5] A = 4 -1 1 5 b=[15 9] b = 15 9 x0=ones(2,1) x0 = 1 1 [x,n]=Jacobi_Solve(A,b,x0) x = 4.00000000000000 0.99999981250000 n = 11 由输出结果可以看出：该方程迭代次数为11，值x1、x2趋近4.00000000000000、0.99999981250000 、利用高斯_赛德尔迭代法，在matlab窗口输入下面命令，得出结果： format long A=[4 -1;1 5] A = 4 -1 1 5 b=[15 9] b = 15 9 x0=zeros(2,1) x0 = 0 0 [x,n]=gaussseide(A,b,x0) x = 3.99999999609375 1.00000000078125 n = 7 由输出结果可以看出：该方程迭代次数为7，值x1、x2趋近3.99999999609375、 1.00000000078125 、利用逐次超松弛迭代法，在matlab窗口输入下面命令，得出结果： format long A=[4 -1;1 5] A = 4 -1 1 5 b=[15 9] b = 15 9 x0=[0 0] x0 = 0 0 [x,n]=SOR_Solve(A,b,x0,1.25) x = 3.99999988432715 1.00000006817230 n = 18 由输出结果可以看出：该方程迭代次数为18，值x1、x2趋近3.99999988432715、1.00000006817230 2、8x-3y=10 -1x+4y=6 解：（1）、利用雅克比迭代法