如下:
function fit_tulun
%多项式拟合方式
x=linspace(0,4*pi,6);
y_jingque=sin(x); %决定稀疏样本点数据
p_poly=polyfit(x,y_jingque,5);
x_poly_fit=linspace(0,4*pi,100);
y_poly_fit=polyval(p_poly,x_poly_fit);
%三次样条拟合方式
sp=csapi(x,y_jingque);
%求三次样条函数的导数.
s_diff=fnder(sp,1);
plot(x_poly_fit,y_poly_fit,'ko',x_poly_fit,y_poly_fit,'b:')
%plot(x_poly_fit,y_poly_fit,'b:')
hold on
fnplt(sp,'r')
fnplt(s_diff,'c')
x1=linspace(0,4*pi,200);
plot(x1,sin(x1),'m','linewidth',1.8)
legend('多项式拟合样本点','多项式拟合曲线','三次样条拟合曲线','三次样条导数曲线','正弦曲线精确图形').
x = 0:10;
y = sin(x);% 已知数据
xx = 0:.25:10; %插值点
yy = spline(x,y,xx); %三次样条插值
plot(x,y,'o',xx,yy) %绘图
圆圈代表原来的数据,曲线代表插值之后的图形曲线。