ex1801.m

来自「MATLAB7.0课本所有程序清单,共37章,特别适合初学者,可以帮助节省大量时」· M 代码 · 共 16 行

%例18-1：产生一条正弦曲线，然后用三次样条插值进行拟合，并对比拟合结合和原数据。
x = 0:10
y = sin(x);
xx = 0:.25:10;
yy = spline(x,y,xx);  %三次样条插值
plot(x,y,'o',xx,yy)   %将实际曲线与样条断点在同一图上显示
%...
pp = spline(x,y);
pp
%...
[breaks,coefs,nploys,ncofs,dim] = unmkpp(pp);    %分解pp形式
%...
pp1 = mkpp(breaks, coefs);    %重构pp形式
pp1