插值MATLAB源程序

1、Lagrange插值MATLAB源程序 function y=lagrange(x0,y0,x); % x0自變量取值向量已知 y0為已知對(duì)應(yīng)x0的函數(shù)取值，x為要求插值點(diǎn)坐標(biāo)n=length(x0); m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j=k p=p*(z-x0(j)/(x0(k)-x0(j);  %插值基函數(shù) end end s=p*y0(k)+s; %lagrange插值多項(xiàng)式 end y(i)=s;end 測(cè)試：x0=0:2 y0=2 3 5 x=0.5Lagran

2、ge(x0,y0,x)x0=0:2 y0=2 3 5 x=0:0.01:2Lagrange(x0,y0,x)  Newton插值MATLAB源程序function f=Newton(x,y,x0) syms t; n = length(x); c(1:n) = 0.0; f = y(1); y1 = 0; l = 1; for(i=1:n-1) for(j=i+1:n) y1(j) = (y(j)-y(i)/(x(j)-x(i); end c(i) = y1(i+1); l = l*(t-x(i); f = f + c(i)*l; simplify(f); y = y1; if(i=

3、n-1) if(nargin = 3) f = subs(f,'t',x0); else f = collect(f); end end endtest: x=1 -1 2 y=0 -3 4 x0=-1:0.1:2Newton(x,y,x0)分段插值MATLAB源程序function y = div_linear(x0,y0,x,n) for i = 1:n-1 if (x >= x0(i) && (x <= x0(i+1) y = (x - x0(i+1)/(x0(i) - x0(i+1)*y0(i) + ( x - x0(i)/(x0(i+1)

4、-x0(i)*y0(i+1); else continue; end endtest：（例題3）x0=-1:0.2:1;y0= 1./(25*x0.2+1); y=interp1(x0,y0,x0,'linear') plot(x0,y0,x0,y,'p');Newton插值MATLAB源程序(2)function = newton(x,y,v)x=input(“X數(shù)組=”); y=input(“Y數(shù)組=”);v=input(“插值點(diǎn)數(shù)值=”); n=length(x); t=zeros(n,n);u=0; for i=1:n t(i,1)=y(i); end

5、for j=2:n for i=2:n if i>=j t(i,j)=(t(i,j-1)-t(i-1,j-1)/(x(i)-x(i-j+1); end end end for k=1:n s=1; m=1; for j=1:k if j<k s=s*(v-x(j); end end m=s*t(k,k); u=u+m; end disp(“插值結(jié)果=”); disp(u); end>>newton Newton插值MATLAB源程序  function yi=newton(x,y,xi) m=length(x);n=length(y); i

6、f m=n error('x and y must same');end f=zeros(n+1,1); k=2; f(1)=y(1) while k=n+1 f(1)=y(k);k,x(k) for i=1:k-1 if i=k-1 f(i+1)=(f(i)-y(i)/(x(k)-x(i); end end cs(i)=f(i+1); y(k)=f(k); k=k+1;end cfwh=0; for i=1:n-2 w=1; for j=1:i w=w*(xi-x(j); end cfwh=cfwh+cs(i)*w; end yi=y(1)+cfwh;x=0:2 y=2 3

7、5 xi=0.5newton(x,y,xi)Newton2)function f = Newton(x,y,x0) syms t;  if(length(x) = length(y)       n = length(x);  c(1:n) = 0.0;  else       disp('x和y的維

8、數(shù)不相等！');      return;  end  f = y(1);  y1 = 0;  l  = 1;  for(i=1:n-1)        for(j=i+1:n)        

9、0;  y1(j) = (y(j)-y(i)/(x(j)-x(i);      end       c(i) = y1(i+1);           l = l*(t-x(i);        f =

10、0;f + c(i)*l;  simplify(f);      y = y1;      if(i=n-1)           if(nargin = 3)           

11、0;   f = subs(f,'t',x0);          else             f = collect(f);              &#

12、160; %將插值多項(xiàng)式展開              f = vpa(f, 6);          end     end endx=0:2 y=2 3 5 x0=0.5Newton(x,y,x0) function  = Newton

13、(x,y,v) x=input('X數(shù)組=：'); y=input('Y數(shù)組=：'); v=input('插值點(diǎn)數(shù)值=：'); n=length(x); t=zeros(n,n); u=0;   for i=1:n       t(i,1)=y(i);   end   for j=2:n  

14、   for i=2:n         if i>=j       t(i,j)=(t(i,j-1)-t(i-1,j-1)/(x(i)-x(i-j+1);         end end   end    

15、60;for k=1:n     s=1;     m=1;      for j=1:k        if j<k        s=s*(v-x(j);        end