楕円のプログラム(一般→標準形）について

質問

こんばんわ、質問させて頂きます。
この掲示板の「楕円のプログラム」を参照して楕円のプログラムを組みました。（EXCELを用いて近いうちにCで組む予定です。)

ある計測点１２点より
楕円の一般方程式に当てはめ、
Ax^2+By^2+Cx+Dy+E=1
として各値A,B,C,D,Eを転置、逆行列、積を用いて解きました。
そして、この値は
A=1/a^2=-3.16665e-7
B=1/b^2=1.4901e-8
C=-2c/a^2=0
D=-2d/b^2=-2.32831e-10
E=c^2/a^2+d^2/b^2=1
となりました。ここで、小文字a,b,c,d,eは楕円の方程式の標準形
(x-c)^2/a^2＋(y-d)^2/b^2=1
となりますが、
Aの値は負となります。つまりa^2は虚数解を持ちます。
これはどういう意味を持っているのでしょうか？
ご教授よろしくお願いします。

用いた計測座標（x,y)
0.004343 　　　　　0.013968823
0.00314534 0.01454172
0.005537 　　　　　0.016603472
0.00640595 0.019503824
0.00554108 0.022148798
0.00296192 0.024042826
0.000376075 0.024724047
-0.002755　 0.024067637
-0.0051356 0.022425276
-0.0061896 0.019473675
-0.0051498 0.016384854
-0.0028986 0.014560819

回答No.1

Tacosan

なんか, 「きちんと計算できているとは思えない」結果なんですけど....
とりあえず, 全ての数値を 1000倍くらいして計算し直してみるとか....

補足

返信ありがとうございます。
値をすべて1000倍にすると
A=4.76e-13
B=-2.13e-14
C=-3.55e-15
D=2.27374e-13
E=1
となりました。

計算内容の方は
「測定点の楕円近似」の質問欄にある式をそのまま使いました。
p[1]=x^2,p[2]=y^2,p[3]=x,p[4]=y,p[5]=1,
q[1]=A,q[2]=B,q[3]=C,q[4]=D,q[5]=E, r=1
ここで、
=mmult(minverse(mmult(transpose(P),P),mmult(transpose(P),r))
としました。上記の計算内容が違うのでしょうか？