how to find volume of the ellipse not in the origin by using the eigenvalues

Question

nadia nadi 2017-4-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/336963-how-to-find-volume-of-the-ellipse-not-in-the-origin-by-using-the-eigenvalues

编辑： nadia nadi 2017-4-25

dear all,

i have matrix A has the information of an ellipse not in the origin. i need to find the volume of the ellipse using the eigenvalues but i got wrong volume. in 2-dimension i have points x enclosed by the ellipse that its information is A and b.

x = [ 0.55  0.0  ;
      0.25   0.35  
      -0.2   0.2     
      -0.25 -0.1   
      -0.0  -0.3  
      0.4  -0.2 ]';
  A=[ 2.4282   -0.0557
     -0.0557    2.9680]; 
  b=[-0.3393; -0.0562];
ei=eig(A)
vol=pi*sqrt(ei(1))*sqrt(ei(2));

the eigenvalues are ei=[ 2.4225, 2.9737 ]; the volume is 8.4320 which i think is wrong. can anyone help me to move the ellipse to the center to find the eigenvalues then the volume please. i tried to do that but its not correct.

many thanks,

Nadia,

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2017-4-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/336963-how-to-find-volume-of-the-ellipse-not-in-the-origin-by-using-the-eigenvalues#answer_264209

编辑：Torsten 2017-4-24

vol=pi/sqrt(det(A))

Best wishes

Torsten.

5 个评论
显示 3更早的评论隐藏 3更早的评论

nadia nadi 2017-4-24

编辑：nadia nadi 2017-4-24

在 MATLAB Online 中打开

Dear Torsten,

thanks for the link, it seems the same matrix A but i have different center. this matrix result from the problem of finding the minimum volume ellipsoid (described by matrix A and vector b) that covers all of the points x by solving the optimization problem

 [n,m]= size(x)
  % Create and solve the model
  cvx_begin
      variable A(n,n) symmetric
      variable b(n)
      maximize( det_rootn( A ) ); %%or (n*log(det_rootn(A)))% minimize( - log_det( A ) )
      subject to
          norms( A * x + b * ones( 1, m ), 2 ) <= 1;
  cvx_end
  scatter ( x(1,:), x(2,:) ); 
    hold on
    k = convhull ( x(1,:), x(2,:) );
    plot ( x(1,k), x(2,k),'k' );
    noangles = 200;
    angles   = linspace( 0, 2 * pi, noangles );
    ellipse  = A \ [ cos(angles) - b(1) ; sin(angles) - b(2) ];
    plot( x(1,:), x(2,:), 'ko', ellipse(1,:), ellipse(2,:), 'b-' );
    axis on
     hold off

you cant run the cvx because it is a package work with Matlab but you can run the plot and see that it is different from the plot in that link.

thank you.

Torsten 2017-4-25

In your problem formulation, you'll have to work with A^2 instead of A. If you do this, you get an approximate volume of 0.436.

Best wishes

Torsten.

nadia nadi 2017-4-25

编辑：nadia nadi 2017-4-25

在 MATLAB Online 中打开

Dear Torsten,

that's great thanks, it gives good approximation. I also used this

ei=eig(A);
VOL=pi*(1/(ei(1)))*(1/(ei(2)))

and it gives me the same from your formula, now for three dimension I will write the two formulas like this

vol=(4/3)*pi/sqrt(det(A^2))
VOL=4/3*pi*(1/(ei(1)))*(1/(ei(2)))*(1/(ei(3)))

and for high dimension I need to rewrite this formula so i write this. do you think is correct?.

m=>4
mm=((m/2)+1)-1;%or(m/2)
mm=ceil(mm);
vol=1/sqrt(det(A^2))*((pi^(m/2))/factorial(mm))

many thanks,

Nadia

请先登录，再进行评论。

how to find volume of the ellipse not in the origin by using the eigenvalues

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

5 个评论
显示 3更早的评论隐藏 3更早的评论

另请参阅

类别

标签

Community Treasure Hunt

how to find volume of the ellipse not in the origin by using the eigenvalues

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

5 个评论 显示 3更早的评论隐藏 3更早的评论

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论