数值分析实验报告--插值与拟合.docVIP

插值与拟合 1. 观察直接利用拉格朗日插值多项式的病态性 实验目的：如果直接利用拉格朗日插值多项式的的定义，即 确定系数，就必须解方程 该方程的系数矩阵为范德蒙（Vandermonde）矩阵，它是严重病态的，因此这种方法不能使用。本实验的目的是观察这种病态性。 实验内容：取不同的n ，并在[0，1]区间上取n个等分点，计算上述系数矩阵A的条件数，画出㏑(cond(A)）～n之间的曲线。（可利用MatLab的库函数cond完成） 实验算法：（源程序） function y=tiaojianshu(x) for i=1:x n=i; B=[]; A=ones(n+1,1); for j=1:n+1 a(j)=(j-1)/n; end for i=1:n+1 B=[B,a(i)]; end B=B; for k=1:n A=[A,B.*A(:,k)]; end A y=log(cond(A)); plot(n,y,.r); xlabel(n); ylabel(1n(cond(A))); hold on end 实验结果： ㏑(cond(A)）～n之间的曲线如下： 实验说明：随着n的值增大，㏑(cond(A)）也线性增大，即cond(A)呈现指数上升趋势。 当n取一个较大的值时，cond(A)会相当大，方程组呈现变态。 2. 牛顿插值法 实验