回 帖 发 新 帖 刷新版面

主题:求助-二次多项式拟合

x: 0.1   0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1
  y: 1.98  3.28 6.16 7.34 7.66 9.58 9.56 9.48 9.30 11.2
 分别用解超定方程和多项式拟合命令进行求解~!!!!  
    各位大虾门 俺的幸福交到你们的手里了!

回复列表 (共12个回复)

11 楼


对于一个系数矩阵是m×n的线性方程组A x = b来说,如果m>n,也就是说方程的个数多于
未知数,则称为超定方程组。
   书上是这么说的 呵呵

12 楼

哦,那就很简单
若用m次多项式拟合:
b0+b1X1+b2X1^2+bmX1^m=Y1
b0+b1X2+b2X2^2+bmX2^m=Y2
......
b0+b2Xn+b2Xn^2+bmXn^m=Yn
用矩阵形式表示即为:
Xb=Y,其中X是n*(m+1)结构矩阵,第1列全是1,第二列则为原来的X数据,第3列是X的平方,第m+1列是X的m次方;Y是Y数据向量。当m<n时,需用最小二乘法求解b,即b=(X'*X)^(-1)*X'*Y,在matlat下
>b=X'*X\X'*Y即可
又及,本例以m=4较为合适,结果将与用多项式拟合(polyfit)的完全一样。

我来回复

您尚未登录,请登录后再回复。点此登录或注册