1、函数插值
一维插值: interp1(x,y,cx,’method’)
一维插值: interp1(x,y,z,cx,cy,’method’)
method:nearest、linear、spline、cubic
例:
clear
echo on
x=-2:0.4:2;
y=[2.8 2.96 2.54 3.44 3.56 5.4
6.0 8.7 10.1 13.3 14.0];
t=-2:0.01:2;
nst=interp1(x,y,t,'nearest');
plot(x,y,'r*',t,nst)
title('最临近点插值')
lnr=interp1(x,y,t,'linear');
figure(2)
plot(x,y,'r*',t,lnr,'b:')
title('线性插值')
spl=interp1(x,y,t,'spline');
figure(3)
plot(x,y,'r*',t,spl)
title('样条插值')
cbc=interp1(x,y,t,'cubic');
figure(4)
plot(x,y,'r*',t,cbc,'k-')
title('三次插值')
2、曲线拟合
多项式拟合:polyfit(x,y,m) 线性:m=1, 二次:m=2, …
例:
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];
A=polyfit(x,y,2)
Z=polyval(A,x);
Plot(x,y,’r*’,x,z,’b’)
你所做的就是选择合适的差值方式,画图,或者插值,再画图 |