资源描述：
《数值分析报告某实验报告材料》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实用标准文案实验一、误差分析一、实验目的1．通过上机编程，复习巩固以前所学程序设计语言及上机操作指令；2．通过上机计算，了解误差、绝对误差、误差界、相对误差界的有关概念；3．通过上机计算，了解舍入误差所引起的数值不稳定性。二．实验原理误差问题是数值分析的基础，又是数值分析中一个困难的课题。在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。因此，选取算法时注重分析舍入误差的影响，在实际计算中是十分重要的。同时，由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差，因此算
2、法的好坏会影响到数值结果的精度。三．实验内容对，计算定积分.算法1：利用递推公式,,取.算法2：利用递推公式.注意到,取.:四.实验程序及运行结果程序一：t=log(6)-log(5);精彩文档实用标准文案n=1;y(1)=t;fork=2:1:20y(k)=1/k-5*y(k-1);n=n+1;endyy=0.0884y=0.0581y=0.0431y=0.0346y=0.0271y=0.0313y=-0.0134y=0.1920y=-0.8487y=4.3436y=-21.6268y=108.2176
3、y=-541.0110y=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;fork=20:-1:2y(k-1)=1/(5*k)-(1/5)*y(k);n=n+1;end运行结果：y=0.08840.05800.04310.03430.02850.02120.01880.01
4、69精彩文档实用标准文案0.01540.01410.01300.01200.01120.01050.00990.00930.0089精彩文档实用标准文案实验二、插值法一、实验目的1、理解插值的基本概念，掌握各种插值方法，包括拉格朗日插值和牛顿插值等，注意其不同特点；2、通过实验进一步理解并掌握各种插值的基本算法。二．实验原理插值法是函数逼近的一种重要方法，它是数值积分、微分方程数值解等数值计算的基础与工具，其中多项式插值是最常用和最基本的方法。拉格朗日插值多项式的优点是表达式简单明确，形式对称，便于记忆，
5、它的缺点是如果想要增加插值节点，公式必须整个改变，这就增加了计算工作量。而牛顿插值多项式对此做了改进，当增加一个节点时只需在原牛顿插值多项式基础上增加一项，此时原有的项无需改变，从而达到节省计算次数、节约存储单元、应用较少节点达到应有精度的目的。三．实验内容1、已知函数表0.561600.562800.564010.565210.827410.826590.825770.82495用二次拉格朗日插值多项式求时的函数近似值。2、已知函数表0.40.550.650.80.90.410750.578150.69
6、6750.888111.02652用牛顿插值多项式求和。四．实验程序及运行结果1.程序如下：function[y,R]=lagranzi(X,Y,x,M)n=length(X);m=length(x);fori=1:mz=x(i);s=0.0;fork=1:np=1.0;q1=1.0;c1=1.0;forj=1:nifj~=kp=p*(z-X(j))/(X(k)-X(j));end精彩文档实用标准文案q1=abs(q1*(z-X(j)));c1=c1*j;ends=p*Y(k)+s;endy(i)=s;e
7、ndR=M*q1/c1;在MATLAB工作窗口输入程序：>>x=0.5635;>>M=1;>>X=[0.56280,0.56401];>>Y=[0.82659,0.82577];>>[y,R]=lagranzi(X,Y,x,M)实验结果如下：y=0.8261R=1.7850e-0071.程序如下：function[y,R]=newcz(X,Y,x,M)n=length(X);m=length(x);fort=1:mz=x(t);A=zeros(n,n);A(:,1)=Y';s=0.0;p=1.0;q1=1
8、.0;c1=1.0;forj=2:nfori=j:nA(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1));endq1=abs(q1*(z-X(j-1)));c1=c1*j;endC=A(n,n);q1=abs(q1*(z-X(n)));fork=(n-1):-1:1C=conv(C,poly(X(k)));d=length(C);C(d)=C(d)+A(k,k);endy(k)=polyva
