matlab散点拟合椭圆,根据离散点的坐标拟合椭圆

clear,clc

fx1=@(b,x)b(1)+b(2)*x+sqrt(b(3)+b(4)*x+b(5)*x.^2);

fx2=@(b,x)b(1)+b(2)*x-sqrt(b(3)+b(4)*x+b(5)*x.^2);

figure(1),clf

b=[304.896155  -0.0468025122  -25124.8078  330.246104  -1.05243991];

X=xlsread('tynh','a1:b156');

x=X(:,1);y=X(:,2);

rg=range(x);

x1=linspace(min(x)-rg/100,max(x)+rg/100,200);

plot(x,y,'o','markersize',10,'markerfacecolor','k')

hold on

y1=fx1(b,x1);y2=fx2(b,x1);

plot(x1,y1,'r-',x1,y2,'r-','linewidth',2.5)

axis tight

figure(2),clf

b=[372.4589157  -0.05152819875  -25265.89178  321.9513439  -0.9973063478];

X=xlsread('tynh','d1:e269');

x=X(:,1);y=X(:,2);

rg=range(x);

x1=linsp