椭圆参数的推导过程

一、用一般式方程的系数计算椭圆参数

二次曲线的一般式方程为Ax^2+Bxy+Cy^2+Dx+Ey+F=0,当B^2-4AC<0时所表示的图形为椭圆。

设椭圆的圆心为(x_0,y_0),长半轴长度为 a,短半轴长度为 b,长半轴与 x 轴的夹角为 θ。

通过将坐标做平移和旋转变化,可以得到标准形式的椭圆方程:

\frac{[(x-x_0)\cosθ+(y-y_0)\sinθ]^2}{a^2}+\frac{[-(x-x_0)\sinθ+(y-y_0)\cosθ]^2}{b^2}=1

展开后,与一般式做对比,并且由于一般式方程乘以任意的非零常数 k 仍成立,可得:

\begin{cases} a^2\sin^2θ+b^2\cos^2θ=kA \\ 2(b^2-a^2)\sinθ\cosθ=kB \\ a^2\cos^2θ+b^2\sin^2θ=kC \\ -2Ax_0-By_0=D \\ -Bx_0-2Cy_0=E \\ k(Ax_0^2+Bx_0y_0+Cy_0^2)-a^2b^2=kF \end{cases}

可以先解出:

\begin{cases} x_0 = \cfrac{2CD-BE}{B^2-4AC} \\\\ y_0 = \cfrac{2AE-BD}{B^2-4AC} \end{cases}

利用\sin2θ=2\sinθ\cosθ,\cos2θ=\cos^2θ-\sin^2θ可得:

\begin{cases} (b^2-a^2)\sin2θ=kB \\ (b^2-a^2)\cos2θ=k(A-C) \\ \end{cases}

所以θ = \frac{1}{2}\arctan(\frac{B}{A-C}),利用\tan(\frac{1}{2}\arctan{x})= \frac{\sqrt{x^2+1}-1}{x}可得:
θ = \arctan\frac{\pm\sqrt{(A-C)^2+B^2}-(A-C)}{B}

利用\sin^2θ+\cos^2θ=1,\sin(\arctan{x})=\dfrac{x}{\sqrt{x^2+1}}可得:

\begin{cases} a^2+b^2=k(A+C) \\ b^2-a^2=\pm k\sqrt{(A-C)^2+B^2} \\ a^2b^2=k(Ax_0^2+Bx_0y_0+Cy_0^2-F)=k(\frac{AE^2-BDE+CD^2}{4AC-B^2}-F) \end{cases}

假设 A>0,否则所有系数乘以-1,最后解出:

\begin{cases} a^2 = \cfrac{2(\frac{AE^2-BDE+CD^2}{4AC-B^2}-F)}{A+C-\sqrt{(A-C)^2+B^2}} \\\\ b^2 = \cfrac{2(\frac{AE^2-BDE+CD^2}{4AC-B^2}-F)}{A+C+\sqrt{(A-C)^2+B^2}} \\ \end{cases}

长半轴对应的倾角为θ = \arctan \frac{C-A-\sqrt{(A-C)^2+B^2}}{B}

二、坐标变换的方法可用参数方程来推导

把标准椭圆的参数方程

\begin{cases} x = a\cos t \\ y = b\sin t \end{cases}

旋转 θ 角,再平移到(x_0,y_0),得到的参数方程为:

\begin{cases} x = a\cos t \cosθ - b\sin t \sinθ + x_0 \\\\ y = a\cos t \sinθ + b\sin t \cosθ + y_0 \end{cases}

解出:

\begin{cases} \cos t = \frac{1}{a} [(x-x_0) \cosθ + (y-y_0) \sinθ] \\\\ \sin t = \frac{1}{b} [-(x-x_0) \sinθ + (y-y_0) \cosθ] \end{cases}

\sin^2t+\cos^2t=1可得:

\frac{[(x-x_0)\cosθ+(y-y_0)\sinθ]^2}{a^2}+\frac{[-(x-x_0)\sinθ+(y-y_0)\cosθ]^2}{b^2}=1

三、用矩阵形式推导椭圆参数

形如ax^2+2bxy+cy^2+2dx+2ey+f=0的二次曲线方程可以写成矩阵形式:

𝘅^TQ𝘅=0

其中𝘅=\begin{bmatrix}x\\y\\1\end{bmatrix},称作(x,y)的齐次坐标,Q=\begin{bmatrix}a&b&d\\b&c&e\\d&e&f\end{bmatrix}为系数矩阵。当\det Q > 0时,Q 表示的二次曲线为椭圆。

如果允许𝘅=\begin{bmatrix}x\\y\\0\end{bmatrix},𝘅 的定义域从二维欧氏空间扩展到二维投影空间。

l=Q𝘅称作极点 𝘅 关于椭圆 Q 的极线。当 𝘅 在 Q 上时,l 经过 𝘅,l^T𝘅=0,且 l 是 Q 在 𝘅 处的切线。当 𝘅 椭圆外面时,𝘅 离椭圆越远,l 越靠近椭圆的圆心。

无穷远点\begin{bmatrix}1\\0\\0\end{bmatrix}的极线ax+by+d=0和无穷远点\begin{bmatrix}0\\1\\0\end{bmatrix}的极线bx+cy+e=0都过圆心。

联立\begin{cases}ax+by+d=0\\\\ bx+cy+e=0\end{cases}可解出椭圆的圆心坐标为:

\begin{cases} x_0 = \cfrac{be-cd}{ac-b^2} \\\\ y_0 = \cfrac{bd-ae}{ac-b^2} \end{cases}

Q 可以分解为:

Q=\begin{bmatrix}1&0&0\\0&1&0\\-x_0&-y_0&1\end{bmatrix}\begin{bmatrix}a&b&0\\b&c&0\\0&0&g\end{bmatrix}\begin{bmatrix}1&0&-x_0\\0&1&-y_0\\0&0&1\end{bmatrix}

其中g=\frac{\det Q}{ac-b^2}= \frac{(ac-b^2)f-(ae^2-2bde+cd^2)}{ac-b^2}

形如E=\begin{bmatrix}a&b&0\\b&c&0\\0&0&g\end{bmatrix}的系数矩阵可表示圆心在原点上的椭圆。

E的本征值和本征向量,Ev=λv,左侧的v是坐标点,可以看作是从原点到极点的向量;右侧的v是极线的参数,即v^T𝘅=0v作为向量则表示垂直于极线的方向。所以本征向量v是作为坐标点的v对应的极线垂直于作为向量的v的向量。当坐标点v在椭圆上时,该点的切线垂直于该点到原点的连线。

仅当极点在轴线上时,对应的极线垂直于轴线。仅当轴线过原点时轴线的方向与轴线上点的坐标值一致。当两条轴线都过原点时,椭圆的圆心在原点上。

因为v是齐次坐标,只有前两个维度表示坐标值或轴线方向,所以有

\begin{bmatrix}a&b\\b&c\end{bmatrix}v=λv

解出来的本征向量就是椭圆轴线的方向:

\begin{cases} v_1 = (a-c - \sqrt{(a-c)^2 + 4 b^2},\ 2b) \\\\ v_2 = (a-c + \sqrt{(a-c)^2 + 4 b^2},\ 2b) \end{cases}

本征值为:

\begin{cases} λ_1 = \frac{1}{2}(a+c - \sqrt{(a-c)^2 + 4 b^2}) \\\\ λ_2 = \frac{1}{2}(a+c + \sqrt{(a-c)^2 + 4 b^2}) \end{cases}

两个正交的轴线方向组成一个旋转矩阵,坐标变换后椭圆的系数矩阵变为\begin{bmatrix}λ_1&0&0\\0&λ_2&0\\0&0&g\end{bmatrix},对比椭圆的标准型\begin{bmatrix}\frac{1}{r_1^2}&0&0\\ 0&\frac{1}{r_2^2}&0\\0&0&-1\end{bmatrix},可知\sqrt{\frac{-g}{λ}}为椭圆半轴的长度。本征值与半轴的平方成反比,所以绝对值较小的本征值对应椭圆的长半轴。当 a 和 c 为正数时0<λ_1<λ_2v_1为长轴方向;当 a 和 c 为负数时λ_1<λ_2<0v_2为长轴方向。

半轴长度的平方为:

\begin{cases} r_1^2 = \frac{2(ae^2-2bde+cd^2-(ac-b^2)f)}{(ac-b^2)(a+c - \sqrt{(a-c)^2 + 4 b^2})} \\\\ r_2^2 = \frac{2(ae^2-2bde+cd^2-(ac-b^2)f)}{(ac-b^2)(a+c + \sqrt{(a-c)^2 + 4 b^2})} \end{cases}

r_1与 x 轴的夹角为:

θ = \arctan\frac{2 b}{a-c - \sqrt{(a-c)^2 + 4 b^2}}=\arctan\frac{c-a - \sqrt{(a-c)^2 + 4 b^2}}{2 b}

四、由椭圆参数计算椭圆方程的系数

x_0,y_0,a,b,θ 已知时,根据前面的推导:

\begin{cases} a^2\sin^2θ+b^2\cos^2θ=kA \\ 2(b^2-a^2)\sinθ\cosθ=kB \\ a^2\cos^2θ+b^2\sin^2θ=kC \\ -2Ax_0-By_0=D \\ -Bx_0-2Cy_0=E \\ k(Ax_0^2+Bx_0y_0+Cy_0^2)-a^2b^2=kF \end{cases}

k=1,得:

\begin{cases} A=a^2\sin^2θ+b^2\cos^2θ \\ B=2(b^2-a^2)\sinθ\cosθ \\ C=a^2\cos^2θ+b^2\sin^2θ \\ D=-2Ax_0-By_0 \\ E=-Bx_0-2Cy_0 \\ F=-\frac{1}{2}(Dx_0+Ey_0)-a^2b^2\end{cases}

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 159,458评论 4 363
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 67,454评论 1 294
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 109,171评论 0 243
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 44,062评论 0 207
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,440评论 3 287
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,661评论 1 219
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,906评论 2 313
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,609评论 0 200
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,379评论 1 246
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,600评论 2 246
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,085评论 1 261
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,409评论 2 254
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,072评论 3 237
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,088评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,860评论 0 195
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,704评论 2 276
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,608评论 2 270

推荐阅读更多精彩内容

  • 为了简便起见,这里的椭球指的是椭球面,如果还要包含椭球内部,则将等于号改为小于等于即可。 标准椭球方程的几何意义 ...
    Azur_wxj阅读 14,953评论 0 12
  • 一、选择题(本大题共4小题,共20.0分) 双曲线的两顶点为,,虚轴两端点为,,两焦点为,,若以为直径的圆内切于菱...
    7300T阅读 766评论 0 13
  • 1.7、坐标系与参数方程 一、坐标系 1、直角坐标系 建立坐标系必须满足的条件任意一点都有确定的坐标与之对应;反之...
    椰子数学阅读 8,777评论 0 7
  • 椭圆的几何参数包括:中心点坐标 、半长轴 、半短轴 、长轴与 X 轴夹角 椭圆方程为: 给定椭圆几何参数,求其方程...
    victorbian阅读 3,874评论 0 1
  • 处理DXF图形时,解析得到的椭圆弧 ellipse 的数据中,起始参数和结束参数是椭圆参数方程中的参数。与圆弧不同...
    师照照阅读 3,109评论 0 0