实验代数插值实验().docVIP

实验6 代数插值（） 姓名 报告日期 5.3 . 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：进一步熟练掌握Lagrange插值算法，提高编程能力和解决插值问题的实践技能。 实验内容：依照Lagrange算法编写基于N+1个点的 的N次Lagrange插值多项式系数计算的程序，以p163例4.7为例进行调试并画出y=f(x)=cox(x)及插值多项式P2(x)和P3(x)的图像。 实验原理 已知N+1个点的 计算Lagrange插值基函数 (k=0,…,N) 计算 N次Lagrange插值多项式函数 P(x)= 实验说明 需要建立两个数组X，Y分别存放（，…，，（，…， 需要考虑多项式乘积运算的算法或者函数。 实验步骤 1 要求上机实验前先编写出程序代码 2 编辑录入程序 3 调试程序并记录调试过程中出现的问题及修改程序的过程 4 经反复调试后，运行程序并验证程序运行是否正确。 5 记录运行时的输入和输出。 实验总结 实验报告：根据实验情况和结果撰写并递交实验报告。 参考程序（MATLAB程序） 分析需要解决的问题，编写代码。 建立zhangyuelagran.m文件，写入程序代码如下： function [C,L]=zhangyuelagran(X,Y) disp(If there is any problem,please consult ZhangYue for help!) w=length(X); n=w-1; L=zeros(w,w); for k=1:n+1 V=1; for j=1:n+1 if k~=j V=conv(V,poly(X(j)))/(X(k)-X(j)); end end L(k,:)=V; end C=Y*L; ③ 在command window里按照如下操作，得到二次插值的结果如下： X=[0.0 0.6 1.2]; Y=[cos(0.0) cos(0.6) cos(1.2)]; [C,L]=zhangyuelagran(X,Y) If there is any problem,please consult ZhangYue for help! C = -0.4004 -0.0508 1.0000 L = 1.3889 -2.5000 1.0000 -2.7778 3.3333 0 1.3889 -0.8333 0 得到插值多项式为 ④在command window里按照如下操作，得到三次插值的结果如下： X=[0.0 0.4 0.8 1.2]; Y=[cos(0.0) cos(0.4) cos(0.8) cos(1.2)]; [C,L]=zhangyuelagran(X,Y) If there is any problem,please consult ZhangYue for help! C = 0.0922 -0.5651 0.0139 1.0000 L = -2.6042 6.2500 -4.5833 1.0000 7.8125 -15.6250 7.5000 0 -7.8125 12.5000 -3.7500 0 2.6042 -3.1250 0.8333 0 得到插值多项式为 ⑤在command window里按照如下操作，分别画出y=cox(x)的图像，插值多项式P2(x)的图像，插值多项式P3(x)的图像，比较三个图像的差距。 ezplot(cos(x),[0.0,1.2]) ezplot(-0.4004*x*x-0.0508*x+1,[0.0,1.2]) ezplot(0.0922*x*x*x-0.5651*x*x+0.0139*x+1,[0.0,1.2]) 实验总结：通过实验可知，通过二次插值，三次插值都能较好的解得答案，但是两种计算方式的误差还有待讨论，可以通过编程实现进行比较，而且不同的插值法通过计算机求解的速度有所不同。通过实验的验证，可以更深刻地理解插值法在求近似值中的重要性，同时其基本原理在计算机计算方程近似解的应用，也更深刻地理解了插值法的计算原理。我们在实验过程中必须不怕困难，勇于创新，细致耐心。