数值计算期末报告书.docVIP

数值计算期末报告书 学 院 专 业 年级班别 1 学 号 学生姓名 指导教师 2013 年 11 月 1 日 编写Gauss-Seidel迭代法求解线性方程组的程序，要求附有算例。 答： （1）问题描述：由Jacobi迭代法中，每一次的迭代只用到前一次的迭代值，每一次迭代充分利用当前最新的迭代值，即在计算第个分量时，用最新分量，代替旧分量，。这样就是解方程组的Gauss-Seidel迭代法。 (2)算例：5X1+2X2+X3=-12 -X1+4X2+2X3=20 2x1-3X2+10X3=3 (3)计算程序: %Gauss_Seidel.m % A为方程的增广矩阵 clc; format long; A=[5 2 1 -12;-1 4 2 20;2 -3 10 3]; [n,m]=size(A); %最多进行50次迭代 Maxtime=100; %控制误差 Eps=10E-5; %初始值迭代值 x=zeros(1,n); disp(x=); %迭代次数小于最小最大迭代次数，进入迭代 for k=1:Maxtime disp(x); for i=1:n s=0.0; for j=1:n if i~=j s=s+A(i,j)*x(j);%计算和 end end x(i)=(A(i,n+1)-s)/A(i,i);%求出此时迭代的值 end %因为方程的精确解为整数，所以这里将迭代结果向整数靠近的误差作为判断迭代是否停止的条件 if sum((x-floor(x)).^2)Eps break; end; end; X=x; disp(迭代结果:); X format short; (4)迭代结果：X = -3.99997381777344 3.00007473745117 2.00001718479004 编写lagrange多项式插值的程序，要求附有算例。 (1)问题描述： 定义在一定区间上的实函数上的有限个节点里求出lagrange插值基函数，然后再构造插值多项式。 (2)计算程序： % M-file, cyLagrange.m function cyLagrange(x,x0,y) %x是指插值点 %x0是给定函数的区域 %y是插值函数 n=length(x0); m=length(x); for i=1:m %分别对每一点进行插值 t=0.0; for j=1:n u=1.0; %进入迭代计算过程 for k=1:n, if k~=j, u=u*(x(i)-x0(k))/(x0(j)-x0(k)); end end t=t+u*y(j); end s=t; end S=s; disp(插值结果:); S (3)输入 x0=[2 5 7]; x =0; fun =@(x) sqrt(x); y=fun(x0); cyLagrange(x,x0,y) (4)计算结果： 插值结果: S = 0.7281 编写分段梯形积分方法的程序，要求附有算例。 （1）问题描述： 梯形积分法是将被积函数近似看成由无数多个小梯形组成，即将无数小梯形的面积之和近似成被积函数的积分值。梯形法数值积分采用的梯形公式是最简单的数值积分公式，函数在区间[a,b]上计算梯形法数值积分表达式为 ..由于用梯形公式来求积分十分粗糙，误差也比较大，后来改进后提出了复合梯形公式：,其中，n为积分区间划分的个数；h为积分步长。所以给出被积函数后要知道积分上下限、知道要将积分区间分割的段数就可以利用梯形积分法求出被积函数的积分。 （2）计算程序： %用来就数值积分--Trapezia(a,b,fx,E,step);分段函数积分 function Trapezia(a,b,fx,E,step); %fx是由syms定义的函数表达式 %a b 是被积函数的上下限 %E是精度 %step是积分区间的分割线