exa4_6.m

来自「MATLAB 的一些编程,基本的程序,为出学MATLAB的初学者」· M 代码 · 共 17 行

clear all;
global x y yy
t=[1 2 3 9 10];
x=0:0.1:1;y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2];
figure;
xx=linspace(0,1);
subplot(2,3,1);plot(x,y,'o');
for ii=2:6
    p=polyfit(x,y,t(ii-1));yy=polyval(p,xx);
    subplot(2,3,ii);plot(x,y,'o',xx,yy,'b-');
end
figure(2);
y1=interp1(x,y,xx,'linear');subplot(1,2,1);plot(x,y,'o',xx,y1,'r-');
y2=interp1(x,y,xx,'spline');subplot(1,2,2);plot(x,y,'o',xx,y2,'r-');