一般多项式拟合,最小平方误差要怎么求?要源代码

一般多项式拟合,最小平方误差要怎么求?要源代码
一般多项式拟合,最小平方误差要怎么求?要源代码

一般多项式拟合,最小平方误差要怎么求?要源代码
>>x0=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24];
>>y0=[8.28,8.76,9.29,10.73,10.88,11.34,11.97,12.02,12.16,12.83,13.90,14.71,16.11,16.40,17.07,6.96,16.88,17.20,19.87,20.19,20.00,19.81,19.40,20.48];
>> plot(x0,y0,x0,y0,'*')
>> p2=polyfit(x0,y0,2);
>>ye=y0-polyval(p2,x0);
>> ye2s=sum(ye.^2);
>> disp(sprintf('误差的平方和=%d',ye2s));
误差的平方和=9.468029e+000
>> plot(x0,y0,x0,y0,'*',x0,polyval(p2,x0),'--')