0

我有一个二维椭圆,由正定矩阵 X 定义如下:如果 x'*X*x <= 1,点 x 在椭圆中。如何在 matlab 中绘制这个椭圆?我做了一些搜索,但发现却很少。

实际上想出了答案:我会将此作为答案发布,但它不会让我(新用户):

经过一番折腾后想通了。基本上,我们将椭圆边界 (x'*X*x = 1) 上的点表示为 X 的特征向量的加权组合,这使得一些数学运算更容易找到这些点。我们可以只写 (au+bv)' X (au+bv)=1 并计算出 a,b 之间的关系。Matlab 代码如下(对不起,它很乱,只是使用了与笔/纸相同的符号):

function plot_ellipse(X, varargin)
% Plots an ellipse of the form x'*X*x <= 1

% plot vectors of the form a*u + b*v where u,v are eigenvectors of X

[V,D] = eig(X);
u = V(:,1);
v = V(:,2);
l1 = D(1,1);
l2 = D(2,2);

pts = [];

delta = .1;

for alpha = -1/sqrt(l1)-delta:delta:1/sqrt(l1)+delta
    beta = sqrt((1 - alpha^2 * l1)/l2);
    pts(:,end+1) = alpha*u + beta*v;
end
for alpha = 1/sqrt(l1)+delta:-delta:-1/sqrt(l1)-delta
    beta = -sqrt((1 - alpha^2 * l1)/l2);
    pts(:,end+1) = alpha*u + beta*v;
end

plot(pts(1,:), pts(2,:), varargin{:})
4

1 回答 1

4

我在搜索这个主题时偶然发现了这篇文章,即使它已经解决了,我想如果矩阵是对称的,我可能会提供另一个更简单的解决方案。

另一种方法是使用在 Matlab 中实现的半定正矩阵 E 的 Cholesky 分解作为chol函数。它计算一个上三角矩阵 R,使得 X = R' * R。使用这个,x'*X*x = (R*x)'*(R*x) = z'*z,如果我们将 z 定义为 R *X。

因此,要绘制的曲线变成了 z'*z=1,这是一个圆。因此,一个简单的解决方案是 z = (cos(t), sin(t)),对于 0<=t<=2 pi。然后乘以 R 的倒数得到椭圆。

翻译成以下代码非常简单:

function plot_ellipse(E)
 % plots an ellipse of the form xEx = 1
 R = chol(E);
 t = linspace(0, 2*pi, 100); % or any high number to make curve smooth
 z = [cos(t); sin(t)];
 ellipse = inv(R) * z;
 plot(ellipse(1,:), ellipse(2,:))
end 

希望这可能会有所帮助!

于 2015-12-16T09:42:46.610 回答