数值分析教学课件shiyan.pptVIP

下面6张是第1题 下面6张是第2题 下面5张是第三题 第二章 插值与拟合 　　数值试验题2 2.1 编制分段线性插值和分段三次Hermite插值程序，对被插函数 插值区间 ，分成10等分，求分段线性插值函数在各节点间中点处的值，并画出分段线性插值函数和 　　的图形。 2.2 给定数据如下表： 编制程序求三次样条插值函数在插值中点的样条函数值，并作点集 和样条插值函数的图形，满足的边界条件为 2.3 对下列数据作三次多项式拟合，取权数 ，给出拟合多项式的系数、平方误差并作离散数据 　和拟合多项式的图形。 x=[-5,-4,-3,-2,-1,0,1,2,3,4,5];y=1./(1+x.^2) y = 0.0385 0.0588 0.1000 0.2000 0.5000 1.0000 0.5000 0.2000 0.1000 0.0588 0.0385 首先算出f(x)的值如下： 分段线性插值 function I=f(x) a=[-5,-4,-3,-2,-1,0,1,2,3,4,5]; b=[ 0.0385,0.0588,0.1000,0.2000,0.5000,1.0000,0.5000,0.2000,0.1000,0.0588,0.0385]; n=length(x); for i=1:n-1 I(i)=b(i)*(x(i)-a(i+1))/(a(i)-a(i+1))+b(i+1)*(x(i)-a(i))/(a(i+1)-a(i)); end x=[-4.5,-3.5,-2.5,-1.5,-0.5,0.5,1.5,2.5,3.5,4.5];I=f(x) I = 0.0486 0.0794 0.1500 0.3500 0.7500 0.7500 0.3500 0.1500 0.0794 0.0486 运行结果如下： function I=f(x) a=[-5,-4,-3,-2,-1,0,1,2,3,4,5];m=[0.0148,0.0277,0.0600,0.1600,0.5000,0, -0.5000,-0.1600,-0.0600,-0.0277,-0.0148] b=[ 0.0385,0.0588,0.1000,0.2000,0.5000,1.0000,0.5000,0.2000,0.1000,0.0588,0.0385]; n=length(x); for i=1:n-1 l(i)=(1+2*(x(i)-a(i+1))/(a(i)-a(i+1))*((x(i)-a(i))/(a(i+1)-a(i))).^2; n(i)=(x(i)-a(i+1))*((x(i)-a(i))/(a(i+1)-a(i))).^2; I(i)=b(i)*a(i)+m(i)*n(i) end Hermite 插值 x=[-4.5,-3.5,-2.5,-1.5,-0.5,0.5,1.5,2.5,3.5,4.5];f(x) ans = Columns 1 through 7 0.0486 0.0794 0.1500 0.3500 0.7500 0.7500 0.3500 Columns 8 through 10 0.1500 0.0794 0.0486 得f(x)图形如下： 输入： fplot(1./(1+x.^2),[-5,5]) x=[-5,-4,-3,-2,-1,0,1,2,3,4,5];y=[ 0.0385,0.0588,0.1000,0.2000,0.5000,1.0000,0.5000,0.2000,0.1000,0.0588,0.0385];plot(x,y) 画图命令如下： 图形如下： function Sanciyangtiao( X, Y, a, b ) % a = f(1), b = f(n) n = length( X ); for i = 1 : n-1 h(i) = X(i+1)-X(i); end for j = 2 : n-1 u(j) = h(j-1)/( h(j-1) + h(j) ); t(j) = 1 - u(j); d(j) = 6*( (Y(j+1) - Y(j))/h(j) - (Y(j) -