《数值分析》实验报告书.docVIP

《数值分析》实验报告 实验一、误差分析 误差问题是数值分析的基础，又是数值分析中一个困难的课题。在实际计算中，如果选用了不同的算法，由于舍入误差的影响，将会得到截然不同的结果。因此，选取算法时注重分析舍入误差的影响，在实际计算中是十分重要的。同时，由于在数值求解过程中用有限的过程代替无限的过程会产生截断误差，因此算法的好坏会影响到数值结果的精度。 一、实验目的 通过上机编程，复习巩固以前所学程序设计语言及上机操作指令； 通过上机计算，了解误差、绝对误差、误差界、相对误差界的有关概念； 通过上机计算，了解舍入误差所引起的数值不稳定性。 二、实验任务 对，计算定积分 . 算法1：利用递推公式 , , 取 . 算法2：利用递推公式 . 注意到 , 取 . 思考：从计算结果看，哪个算法是不稳定的，哪个算法是稳定的。 算法1： t=log(6.0)-log(5.0); n=0; y=zeros(1,21); y(1)=t; for k=2:21 y(k)=1/k-5*y(k-1); n=n+1; end y(1:6) y(7:21) 运行结果： ans = 0.1823 -0.4116 2.3914 -11.7069 58.7343 -293.5049 算法2： y=zeros(21,1); n=1; y1=(1/105+1/126)/20; for k=21:-1:2 y(k-1)=1/(5*k)-y(k)/5; n=n+1; end 运行结果： y =0.0884 0.0580 0.0431 0.0343 0.0285 0.0243 0.0212 0.0188 0.0169 0.0154 0.0141 0.0130 0.0120 0.0112 0.0105 0.0099 0.0093 0.0089 0.0081 0.0095 0 由数据对比可知，算法2较为稳定。 实验二、插值法 插值法是函数逼近的一种重要方法，它是数值积分、微分方程数值解等数值计算的基础与工具，其中多项式插值是最常用和最基本的方法。拉格朗日插值多项式的优点是表达式简单明确，形式对称，便于记忆，它的缺点是如果想要增加插值节点，公式必须整个改变，这就增加了计算工作量。而牛顿插值多项式对此做了改进，当增加一个节点时只需在原牛顿插值多项式基础上增加一项，此时原有的项无需改变，从而达到节省计算次数、节约存储单元、应用较少节点达到应有精度的目的。 一、实验目的 1、理解插值的基本概念，掌握各种插值方法，包括拉格朗日插值和牛顿插值等，注意其不同特点； 2、通过实验进一步理解并掌握各种插值的基本算法。 二、实验任务 已知函数表 0.56160 0.56280 0.56401 0.56521 0.82741 0.82659 0.82577 0.82495 用二次拉格朗日插值多项式求时的函数近似值。 已知函数表 0.4 0.55 0.65 0.8 0.9 0.41075 0.57815 0.69675 0.88811 1.02652 用牛顿插值多项式求和。 1. function [y,R]=lagranzi(X,Y,x,M) x=0.5635; M=2; X=[0.56160,0.56280,0.56401,0.56521]; Y=[0.82741,0.82659,0.82577,0.82495]; 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~=k p=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; end R=M.*q1./c1; 运行结果： ans = 0.8261 2. N3(0.596) function [y,R]= newcz(X,Y,x,M) x=0.596; M=3; X=[0.4,0.65,0.9]; Y=[0.41075,0.69675,1.02652]; n=length(X); m=length(x);