Runge函数的插值方法对比.docxVIP

Runge函数的插值方法对比 一、对Runge函数的牛顿插值 等距牛顿插值图像 编制程序计算均差 clear,clc a=-l; b=l; n=20; x=a:(b-a)/n:b; %插值节点 %y=sin(x); y=l./(l+25*x.A2); plot(x,y/b*) %用蓝色线作被插函数图象 hold on z=a:(b-a)/(2*n):b; n=length(x); for j=2:n for i=n:-l:j y(i)=(y(i)?y(i?l))/(x(i)?x(i?j+l));% 计算差商 end end u=y(n); m=lcngth(z); for j=l:m for i=n-l:-l:l u二y(i)+u*(z(j)?x(i)); %计算牛顿插值多项式的值 V(j)=U； end u=y(n); end plot(z, V,巧 ％用红色线作牛顿插值多项式图象 hold off 得20次多项式为 Ln(x)= 1.-.101 e7*xA 18+.260e6*xA20-. 144e7*xA 14+.164e7*xA 16+47 l.*xA4-.612e4*xA6 +.493e5*xA8-.245e6*xA 10+.757e6*xA 12-24.1 *xA2 Chebyshev多项式零点插值 由于算式结果比较复杂，Subs命令不太好使，得修改程序如下 clear,clc close all t=cos((2*(0:l:20)+l)/42*pi); %t=-l:().l:l; Rx=l./(l+25*t.A2); %Rx=sin(t); n=length(t)-l; k=l; for x=-l:0.01:l for i=l:n+l temp= 1; for j=l:n+l ifQ~=i) temp=l(i); end end end Ln=Rx*l; yi(k)=Ln; k=k+1; end % Ln=simple(Ln); % cof=sym2poly(Ln); %vpa(Ln,3); plot(t Rx,To) hold on % yi = interp 1 (x^Rx^xi/spline*); xi=-1:0.01:1 plot(xi,yi,b*) hold off 插值结果 等距分段线性插值 这个可以直接连线接起来就可以了 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 -1 -0.8 -0.6 * * * 4- * * -0.4 -0.2 0 0.2 0.4 0.6 0.8 f 三次样条插值 function yi=cubic_spline(x,y,ydot,xi) %三次样条插值公式(三弯矩方程，笫一类边界条件)，-其中, %x为向量，全部的插值节点； %y为向虽，插值节点处的函数值； %ydot为向屋，端点处的导数值， %如果此处值缺省，则用均差代替导数， %xi为标量，自变量x; %yi为xi处的两数估计值。 n=length(x); ny=length(y); %输入的插值点与它的函数值应冇相同的个数。 if n ?二 ny errorCThe lengths of X and Y must be equal!1); return; end %如果没冇给出y的导数值，则用均差代替导数 if iscmpty(ydot)==l ydot=[(y(2)-y( 1 ))/(x(2)-x( 1)) (y(n)-y(n-l))/(x(n)-x(n-l))]; end h=zeros( 1 ,n);lambda=ones(l ,n);mu=ones(l ,n); M=zeros(n,l );d=zeros(n, 1); for k=2:n h(k)=x(k)-x(k?l); %插值点必须互异 if abs(h(k))eps crror(thc DATA is error!); return; end end for k=2:n-l lambda(k)=h(k+ l)/(h(k)+h(k+1)); mu(k)= 1 -lambda(k); d(k)=6/(h(k)+h(k+l))... *((y(k+ l)-y(k))/h(k+ l)-(y(k)-y(k-1 ))/h(k)); end d( 1 )=6/h(2)*((y(2)-y(l ))/h(2)-ydot( 1 ))； d(n)=6/h(n)*(ydot(2)-(y(n)-y(n-l ))/h(n)); A=diag(2*ones( 1 ,n)); for i=l :n-l A(i,i+ l)=lambda(i);A(i+1 ,i)=mu(i+1); end M=A\d; for k=2:n if x(k-l)=xi xi=x(k)