cn.m

matlab的拟合工具包

% 数据点
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