用的x、y的散点坐标,去拟合椭圆的参数,即长半轴a、短半轴b,中心坐标(h,k)。要拟合这些参数,可以按下列步骤来实现:
首先,我们先要自定义椭圆的一般式函数,即
func=inline('(X(:,1)-a(1)).^2/a(2).^2+(X(:,2)-a(3)).^2/a(4).^2-1','a','X');
其二,根据已知的x、y的散点坐标,利用nlinfit
()非线性回归函数,拟合出a、b、h、k数值。
然后,利用得到a、b、h、k计算出拟合后的y值。
最后,用polt()函数在同一窗口绘制原始数据曲线与拟合数据曲线
运行程序后可得到如下结果。
![](https://iknow-pic.cdn.bcebos.com/7c1ed21b0ef41bd50855fc905eda81cb39db3d4b?x-bce-process=image%2Fresize%2Cm_lfit%2Cw_600%2Ch_800%2Climit_1%2Fquality%2Cq_85%2Fformat%2Cf_auto)