《数值分析报告某实验报告材料》由会员分享，可在线阅读，更多相关《数值分析报告某实验报告材料（14页珍藏版）》请在金锄头文库上搜索。
1、word实验一、误差分析一、实验目的1通过上机编程，复习巩固以前所学程序设计语言与上机操作指令；2通过上机计算，了解误差、绝对误差、误差界、相对误差界的有关概念；3通过上机计算，了解舍入误差所引起的数值不稳定性。二实验原理误差问题是数值分析的根底，又是数值分析中一个困难的课题。在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。因此，选取算法时注重分析舍入误差的影响，在实际计算中是十分重要的。同时，由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差，因此算法的好坏会影响到数值结果的精度。三实验内容对，计算定积分 .算法1：利用递推公式 , ,取 .算法2：
2、利用递推公式.注意到,取 .:程序一：t=log(6)-log(5);n=1;y(1)=t;for k=2:1:20 y(k)=1/k-5*y(k-1); n=n+1;endyy =2.7051e+003y =-1.3526e+004y =6.7628e+004y =-3.3814e+005y =1.6907e+006y =-8.4535e+006y =4.2267e+007程序2：y=zeros(20,1);n=1;y1=(1/105+1/126)/2;y(20)=y1;for k=20:-1:2 y(k-1)=1/(5*k)-(1/5)*y(k); n=n+1;end运行结果：y =实验二
3、、插值法一、实验目的1、理解插值的根本概念，掌握各种插值方法，包括拉格朗日插值和牛顿插值等，注意其不同特点；2、通过实验进一步理解并掌握各种插值的根本算法。二实验原理插值法是函数逼近的一种重要方法，它是数值积分、微分方程数值解等数值计算的根底与工具，其中多项式插值是最常用和最根本的方法。拉格朗日插值多项式的优点是表达式简单明确，形式对称，便于记忆，它的缺点是如果想要增加插值节点，公式必须整个改变，这就增加了计算工作量。而牛顿插值多项式对此做了改良，当增加一个节点时只需在原牛顿插值多项式根底上增加一项，此时原有的项无需改变，从而达到节省计算次数、节约存储单元、应用较少节点达到应有精度的目的。三实
4、验内容1、 函数表用二次拉格朗日插值多项式求时的函数近似值。2、 函数表用牛顿插值多项式求和。四实验程序与运行结果1. 程序如下：function y,R=lagranzi(X,Y,x,M)n=length(X); m=length(x);for i=1:m z=x(i);s=0.0; for k=1:n p=1.0; q1=1.0; c1=1.0;for j=1:n if j=kp=p*(z-X(j)/(X(k)-X(j); end q1=abs(q1*(z-X(j);c1=c1*j; end s=p*Y(k)+s; end y(i)=s;endR=M*q1/c1;在MATLAB工作窗口输入
5、程序： x=0.5635; M=1; X=0.56280,0.56401; Y=0.82659,0.82577; y,R=lagranzi(X,Y,x,M)实验结果如下：2. 程序如下： function y,R= newcz(X,Y,x,M)n=length(X); m=length(x);for t=1:m z=x(t); A=zeros(n,n);A(:,1)=Y; s=0.0; p=1.0; q1=1.0; c1=1.0; for j=2:n for i=j:n A(i,j)=(A(i,j-1)-A(i-1,j-1)/(X(i)-X(i-j+1); end q1=abs(q1*(z-X
6、(j-1);c1=c1*j; end C=A(n,n);q1=abs(q1*(z-X(n);for k=(n-1):-1:1C=conv(C,poly(X(k);d=length(C); C(d)=C(d)+A(k,k);end y(k)= polyval(C, z);endR=M*q1/c1;在MATLAB工作窗口输入程序：三次牛顿差值程序如下： X=0.4,0.55,0.6; Y=0.41075,0.57815,0.69675; y,R= newcz(X,Y,x,M)实验结果如下：四次牛顿差值程序如下： x=0.895; X=0.55,0.65,0.8,0.9; Y=0.57815,0.6
7、9675,0.88811,1.02652;实验结果如下：实验三、解线性方程组的直接法一、实验目的1了解求线性方程组的直接法的有关理论和方法；2会编制列主元消去法、LU分解法的程序；二实验原理解线性方程组的直接法是指经过有限步运算后能求得方程组准确解的方法。但由于实际计算中舍入误差是客观存在的，因而使用这类方法也只能得到近似解。目前较实用的直接法是古老的高斯消去法的变形，即主元素消去法与矩阵的三角分解法。引进选主元的技巧是为了控制计算过程中舍入误差的增长，减少舍入误差的影响。一般说来，列主元消去法与列主元三角分解法是数值稳定的算法，它具有准确度较高、计算量不大和算法组织容易等优点，是目前计算机上
8、解中、小型稠密矩阵方程组可靠而有效的常用方法。三实验内容1.用列主元高斯消去法求解方程组.2用矩阵直接三角分解法求解方程组，其中， .四实验程序与运行结果function RA,RB,n,X=liezhu(A,b)B=A b; n=length(b); RA=rank(A); RB=rank(B);zhica=RB-RA;if zhica0,disp(请注意：因为RA=RB，所以此方程组无解.)returnendif RA=RB if RA=ndisp(请注意：因为RA=RB=n，所以此方程组有唯一解.) X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1Y,
9、j=max(abs(B(p:n,p); C=B(p,:);B(p,:)= B(j+p-1,:); B(j+p-1,:)=C;for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);endend b=B(1:n,n+1);A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)/A(q,q); endelse disp(请注意：因为RA=RB A=0.101 2.304 3.555;-1.347 3.712 4.
10、623;-2.835 1.072 5.643;b=1.183;2.137;3.035;RA,RB,n,x=liezhu(A,b)运行结果RA = 3RB = 3n = 3x =function X=LUjfcz(A,b)n,n =size(A); X=zeros(n,1); Y=zeros(n,1); C=zeros(1,n);r=1:n;for p=1:n-1max1,j=max(abs(A(p:n,p);C=A(p,:); A(p,:)= A(j+p-1,:); A(j+p-1,:)=C;g=r(p); r(p)= r(j+p-1);r(j+p-1)=g;if A(p,p)=0 disp(
11、A是奇异阵,方程组无唯一解); break;endfor k=p+1:nH= A(k,p)/A(p,p); A(k,p) = H; A(k,p+1:n)=A(k,p+1:n)- H* A(p,p+1:n);endendY(1)=b(r(1);for k=2:nY(k)= b(r(k)- A(k,1:k-1)* Y(1:k-1); endX(n)= Y(n)/ A(n,n);for i=n-1:-1:1 X(i)= (Y(i)- A(i, i+1:n) * X (i+1:n)/ A(i,i);end在MATLAB工作窗口输入程序： A=1 2 -12 8;5 4 7 -2;-3 7 9 5;6
12、-12 -8 3;b=27;4;11;49;X=LUjfcz(A,b)运行结果X =实验四、解线性方程组的迭代法一、实验目的 1、熟悉迭代法的有关理论和方法；2、会编制雅可比迭代法、高斯-塞德尔迭代法的程序；3、注意所用方法的收敛性与其收敛速度问题。二实验原理解线性方程组的迭代法是用某种极限过程去逐步逼近线性方程组准确解的方法，即是从一个初始向量 出发，按照一定的迭代格式产生一个向量序列 ，使其收敛到方程组 的解。迭代法的优点是所需计算机存储单元少，程序设计简单，原始系数矩阵在计算过程中始终不变等。但迭代法存在收敛性与收敛速度问题。迭代法是解大型稀疏矩阵方程组的重要方法。.2用龙格-库塔方法求解function X=jacdd(A,b,X0,P,wucha,max1)n m=size(A);for j=1:ma(j)=sum(abs(A(:,j)-2*(abs(A(j,j);endfor k=1:max1k for j=1:mX(j)=(b(j)-A(j,1:j-1,j+1:m)*X0(1:j-1,j+1:m)/A(j,j);endX,djwcX=norm(X-X0,P); xdwcX=djwcX/(norm(X,P)+eps); X0=X;X1=Ab;if (djwcXwucha)&(xdwcXwucha)&(xdwcXwu
