>> % 多项式拟合
>> % 给定数据对
>> x0=0:0.1:1;
>> y0=[-0.447,1.978,3.11,5.25,5.02,4.66,4.01,4.58,3.45,5.35,9.22];
>> % 求拟合多项式
>> n=3;
>> p=polyfit(x0,y0,n)
p =
56.6915 -87.1174 40.0070 -0.9043
>> % 图示拟合情况
>> xx=0:0.01:1;
>> yy=polyval(p,xx);
>> plot(xx,yy,'-b',x0,y0,'.r','MarkerSize',20)
>> legend('拟合曲线','原数据','location','southeast')
>> xlabel('x')