Question

知道一组数据，如何用matlab进行三次样条函数拟合？

Answer 1

如下：

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('多项式拟合样本点','多项式拟合曲线','三次样条拟合曲线','三次样条导数曲线','正弦曲线精确图形').

Answer 2

x = 0:10; 

y = sin(x);% 已知数据

xx = 0:.25:10; %插值点

yy = spline(x,y,xx); %三次样条插值

plot(x,y,'o',xx,yy) %绘图

圆圈代表原来的数据，曲线代表插值之后的图形曲线。