数值逼近课程设计.docVIP

第一题 龙格现象 定义在区间[-1,1]上。这是一个很光滑的函数，它的任意阶导数都存在。对它在[-1,1]上做等距节点插值时，会发现插值多项式在变量靠近-1和1时余项会随着n的增大而很大；另一个现象就是插值多项式随节点增多而振动更多。这种插值多项式当节点增多时反而不能更好的接近被插函数的现象称为Runge现象。 function f=runge(n) syms t z z1; x=-1:2/n:1; y=1./(1+25*x.^2); z=0; for i=1:(n+1) z1=y(1,i); for j=1:(n+1) if(i~=j) z1=z1.*(t-x(1,j))./(x(1,i)-x(1,j)); end end z=z+z1; end f=expand(z); 在命令窗口输入命令： y1=runge(4) y1=1250/377*t^4-3225/754*t^2+1; y2=runge(8) y2=1+228601250/3725137*t^4-383000000/3725137*t^6+200000000/3725137*t^87450274*t^2; y3=runge(12) y3=1+367051586875/1847048164*t^4-107641853578125/112669938004*t^6+62017871484375/28167484501*t^8+25628906250000/28167484501*t^12-551599221900/28167484501*t^2-65809335937500/28167484501*t^10; x=-1:0.01:1; y=(1+25*x.^2).^(-1); t=-1:0.01:1; y1=1250/377.*t.^4-3225/754.*t.^2+1; y2=1+228601250/3725137.*t.^4-383000000/3725137.*t.^6+200000000/3725137.*t.^87450274.*t.^2; y3=1+367051586875/1847048164.*t.^4-107641853578125/112669938004.*t.^6+62017871484375/28167484501.*t.^8+25628906250000/28167484501.*t.^12-551599221900/28167484501.*t.^2-65809335937500/28167484501.*t.^10; plot(x,y,t,y1,:,t,y2,-.,t,y3,*) legend(f(x),n=4,n=8,n=12) 结论： 由上图可以看出用高次插值多项式是不妥当的，误差会随着阶数的提高越来越严重。因此，为提高插值的精度，可以采用分段插值。 第二题 Chebyshev多项式 画出Chebyshev多项式在n=2，3，4，5时的图像 x=-1:0.1:1; T2=cos(2*acos(x)); T3=cos(3*acos(x)); T4=cos(4*acos(x)); T5=cos(5*acos(x)); plot(x,T2,x,T3,-.,x,T4,:,x,T5,-) legend(n=2,n=3,n=4,n=5) 第三题 Remez算法 1.求函数f(x)=在[-1,1]上的二次多项式逼近。 编写M文件 function remezs(F,x1,x2,x3,x4,e,min,max) format long; syms x ; y1=F; xA=5; xB=5; xC=5; xD=5; m=0; n=0; X1(1)=x1;X2(1)=x2;X3(1)=x3;X4(1)=x4;P1(1)=subs(F,x,x1);P2(1)=subs(F,x,x2);P3(1)=subs(F,x,x3);P4(1)=subs(F,x,x4); d=2; while(abs(abs(subs(y1,x,xA))-abs(m))e | abs(abs(subs(y1,x,xB))-abs(m))e | abs(abs(subs(y1,x,xC))-abs(m))e | abs(abs(subs(y1,x,xD))-abs(m))e) a=[1,1,x1,x1^2;1,-1,x2,x2^2;1,1,x3,x3^2;1,-1,x4,x4^2]; z1=subs(F,x,x1); z2=su