4-3.m

来自「MATLAB科学计算与工程分析源代码源程序3」· M 代码 · 共 8 行

x=0:0.1:pi;
y=sin(x);
[p,mu]=polyfit(x,y,9)

x1=0:0.1:2*pi;
y1=sin(x1);
y2=polyval(p,x1);
plot(x1,y,'k*',x1,y1,'k-')