一般橢圓公式：
ax^2+2bxy+cy^2+2dx+2ey+1=0
C=[a b;b c];（表示第一行是a b；第二行是b c）
D=[a b d;b c e;d e 1];（表示第一行是a b d；等等）
其中記矩陣C的特征值為h1,h2
那么 長半軸=（D的行列式除以h1與C的行列式的乘積）的絕對值；
短半軸=（D的行列式除以h2與C的行列式的乘積）的絕對值；
請問以上求長、短半軸的公式是怎么得到的?
強(qiáng)行計(jì)算吧
設(shè)
C=T'*diag{h1,h2}*T,其中T正交.
由于原方程為
[x,y]C[x,y]'+2[d,e][x,y]'+1=0
做坐標(biāo)變換[nx,ny]=[x,y]T'
于是方程變成
[nx,ny]*diag{h1,h2}*[nx,ny]'+2[d,e]T[nx,ny]'+1=0
設(shè)[nd,ne]=[d,e]T
即
h1*nx^2+h2*ny^2+2*nd*nx+2*ne*ny+1=0
h1*(nx+nd/h1)^2+h2*(ny+ne/h2)^2+1-nd^2/h1-ne^2/h2=0
于是得到兩個半軸長度為
sqrt((nd^2/h1+ne^2/h2-1)/h1)
sqrt((nd^2/h1+ne^2/h2-1)/h2)
現(xiàn)在計(jì)算
|D|=|h1,0,nd;0,h2,ne;nd,ne,1|=(h2 - ne^2)*h1 - h2*nd^2
代入即可