Question

您好，请问您在13年回答的问题“matlab拟合出椭圆方程”，参照的公式是什么？

题目：给出一些椭圆上离散的点的横纵坐标，怎么用matlab拟合出椭圆方程
链接：http://zhidao.baidu.com/question/571561726.html?from=pubpage&msgtype=2
请问您参照的公式或文献是什么，我拟合长轴与x轴交角时，误差较大。我不擅长数学推导，向您请教。

• 你的回答被采纳后将获得：
• 系统奖励15（财富值+成长值）+难题奖励30（财富值+成长值）

Answer 1

clear, clc
x=[-1.5484, -1.3843, -1.1175, -0.76781, -0.36121, 0.072186, 0.50023, 0.89117, 1.216, 1.4507, 1.5777, 1.5878, 1.4801, 1.2626, 0.9515, 0.56982, 0.14587, -0.2889, -0.70224, -1.0635, -1.3459, -1.5284, -1.5977]';
y=[3.0182, 1.3891, -0.34312, -2.0499, -3.6046, -4.8919, -5.8165, -6.3097, -6.3349, -5.8903, -5.0088, -3.7559, -2.2244, -0.52788, 1.2077, 2.8538, 4.2882, 5.4045, 6.1201, 6.3817, 6.17, 5.5008, 4.4235]';
fx1=@(b,x)b(1)*x +sqrt (b(2) + b(3)*x.^2);
fx2=@(b,x)b(1)*x - sqrt(b(2) + b(3)*x.^2);
b=[-2.7007649, 22.106671, -8.6561792];
figure(1),clf
plot(x,y,'o','markerfacecolor','k')
rg=range(x);
x1=min(x)-rg/500:rg/1550:max(x)+rg/300;
y1=fx1(b,x1);
y2=fx2(b,x1);
hold on
plot(x1,y1,'r-',x1,y2,'r-','linewidth',2)
legend('data','fit','location','best')