大工矩阵与数值分析实验报告.docVIP

矩阵与数值分析实验报告 姓 名： 学 号： 院 系：机械工程学院 班 级：003班 老 师： 2014年12月20日 第 PAGE \* MERGEFORMAT 29页，共29页 方程在x=3.0附近有根，试写出其三种不同的等价形式以构成两种不同的迭代格式，再用这两种迭代求根，并绘制误差下降曲线，观察这两种迭代是否收敛及收敛的快慢。 解：三种不同的等价形式： 两种不同的迭代方法： 方法一：简单迭代法 选择作为简单迭代法的迭代格式，时，迭代终止，这里取δ=1E-5，简单迭代法的程序如下： 程序1.1：简单迭代法： %第一题简单迭代法程序 clc clear del=1E-5; x=3; e=1E10; n=0; while edel n=n+1; t=sqrt((2*x^3-19*x+42)/5); e=abs(t-x); E(n)=e; x=t; end X,n; plot(E,g,linewidth,3) title(误差下降曲线图) xlabel(迭代次数) ylabel({迭,代,误,差},rotation,0)计算结果1.1： x = 2.0000 n = 12结果分析1.1： 这里的迭代结果是收敛的的，但是我也试过第3个迭代方程，发现是发散的，由计算结果可知，迭代法的收敛性与构造的迭代函数是直接相关的，我们应该根据定理4.5来构造迭代函数，这样可以得到收敛的结果。 方法2：Newton迭代法 当时，迭代终止，x0=3,δ=1e-5. 程序1.2：Newton迭代法： %第一题Newton迭代法程序 clc clear %Newton迭代法 del=1E-5; x=3; e=1E10;%条件判断值，足够大即可 n=0; while edel n=n+1; t=x-(2*x^3-5*x^2-19*x+42)/(6*x^2-10*x-19); e=abs(t-x); E(n)=e; x=t; end x n plot(E,g,linewidth,3) title(误差下降曲线图) xlabel(迭代次数) ylabel({迭,代,误,差},rotation,0)计算结果1.2： x = 3.5000 n =6 结果分析1.2： Newton迭代法是二阶收敛，弦截法是在前者基础上变化来的，收敛阶1.618，在该题中使用简单迭代法的迭代步数为12，Newton迭代法为6，很明显，Newton迭代法收敛比较快。 2.用复化梯形公式、复化辛普森公式、龙贝格公式求下列定积分，要求绝对误差为，并将计算结果与精确解进行比较： （1） （2） 解：程序和结果如下： 程序2.1.1：使用复化梯形公式解第一个等式： function s=tixing1(a,b,n) format long a=1; b=2; n=150000; eror=100; h=(b-a)/(2*n); index1=(a+h):(2*h):(b-h); index2=(a+2*h):(2*h):(b-2*h); s1=sum(subs(fun2,index1)); s2=sum(subs(fun2,index2)); s=h*(subs(fun2,a)+subs(fun2,b)+2*s2) c=exp(4) error=abs(c-s) 计算结果2.1.1： s = 54.598150039042402 c = 54.598150033144236 error = 5.898165511553088e-009程序2.1.2：使用复化Simpson公式解第一个等式： function s=simpon(a,b,n) format long a=1; b=2; n=370; eror=100; h=(b-a)/(2*n); index1=(a+h):(2*h):(b-h); index2=(a+2*h):(2*h):(b-2*h); s1=sum(subs(fun2,index1)); s2=sum(subs(fun2,index2)); s=h*(subs(fun2,a)+subs(fun2,b)+4*s1+2*s2)/3 c=exp(4) error=abs(c-s) 计算结果2.1.2： s = 54.598150034153207 c