cn.asv

来自「matlab的拟合工具包」· ASV 代码 · 共 83 行

% 数据点
x=[1,2,4,7,9,12,13,15,17];
f=[1.5,3.9,6,11.7,12.6,18.8,20.3,20.6,21.1];
% 设定
axis([0,18, 1,22])
xlabel('x')
ylabel('f')
% 插值点
y=1:0.1:17;      
画出已知数据点、
figure(1)
plot(x,f,'o')
gtext('已知数据点')

hold on;

pause

bb1=interp1(x,f,y,'nearest')
plot(y,bb1)
gtext('最近邻点插值')

hold on;

pause

a=polyfit(x,f,3)
aa=polyval(a,y)
plot(y,aa)
gtext('三次曲线拟合')

hold off;

pause

figure(2)
plot(x,f,'o')
gtext('已知数据点')

hold on;

pause

bb1=interp1(x,f,y,'linest')
plot(y,bb1)
gtext('直线插值')

hold on;
pause

a=polyfit(x,f,3)
aa=polyval(a,y)
plot(y,aa)
gtext('三次曲线拟合')

hold off;

pause

figure(3)

plot(x,f,'o')
gtext('已知数据点')
hold on;

pause

bb1=interp1(x,f,y,'spline')
plot(y,bb1)
gtext('样条插值')

hold on;

pause

a=polyfit(x,f,3)
aa=polyval(a,y)
plot(y,aa)
gtext('三次曲线拟合')

hold off;

pause