对x、y进行n阶幂函数的拟合 p ( x ) = p 1 x n + p 2 x n − 1 + ⋯ + p n x + p n + 1 p(x)=p_1x^n+p_2x^{n-1}+ \dots+p_nx+p_{n+1} p(x)=p1xn+p2xn−1+⋯+pnx+pn+1
[p,S]=polyfit(x,y,n)
结构体S可用作polyval的输入来获取误差估计值
[p,S,mu]=polyfit(x,y,n)
返回的mu(1)为mean(x);mu(2)为std(x)
举例
x =[10111213141516171819202122232425];
y =[19241825202837403020152623171417];plot(x,y,'*');
hold on
y1 =polyfit(x,y,1);
f1 =poly2sym(y1);
y2 =polyfit(x,y,2);
f2 =poly2sym(y2);
y3 =polyfit(x,y,3);
f3 =poly2sym(y3);fplot(f1,'r');fplot(f2,'b');fplot(f3,'k');axis([10251045]);