（MATLAB）散点椭圆拟合与绘制代码

by 今天不飞了

二维平面散点的单椭圆拟合

散点椭圆拟合

1. 拟合

function W = fitellipse(x,y)

% 构造矩阵
D = [x.*x, x.*y, y.*y, x, y,ones(size(x))];
S = D'*D;
G = zeros(6);
G(1,3) = 2; G(3,1) = 2; G(2,2) = -1;

% 求解
[vec, val] = eig(S\G);
[~, idx] = find(val>0&~isinf(val));
W = vec(:,idx);
W = sqrt(1/(W'*S*W))*W;
end

2. 参数计算

中心，长轴，短轴，倾角

function [Center,Axis,Theta] = calellipseparams(W)

a = W(1); 
b = W(2); 
c = W(3); 
d = W(4); 
e = W(5); 
f = W(6);

% 中心
cx