任意椭圆和直线求交点方程组解析

以下为标准椭圆方程与一般直线方程的交点解析：
$$\{\begin{array}{c}y=kx+m\text{直线方程}\\ \frac{x^2}{a^2}+\frac{y^2}{b^2}=1\text{椭圆方程}\end{array}$$
将直线方程带入椭圆方程得：
$$\Rightarrow \frac{x^2}{a^2}+\frac{{\left(kx+m\right)}^2}{b^2}=1$$

$$\Rightarrow \frac{x^2}{a^2}+\frac{{\left(kx+m\right)}^2}{b^2}=1$$

$$\Rightarrow b^2{x}^2+a^2{\left(kx+m\right)}^2=a^2{b}^2$$

$$\Rightarrow b^2{x}^2+a^2\left(k^2{x}^2+2kmx+m^2\right)=a^2{b}^2$$

$$\Rightarrow b^2{x}^2+a^2{k}^2{x}^2+2a^{2}kmx+a^2{m}^2-a^2{b}^2=0$$

$$\Rightarrow \left(b^2+a^2{k}^2\right)x^2+2a^{2}kmx+a^2\left(m^2-b^2\right)=0$$

$$\Rightarrow \text{设}\{\begin{array}{l}A=\left(b^2+a^2{k}^2\right)\\ B=2a^{2}km\\ C=a^2\left(m^2-b^2\right)\end{array}$$

$$\Rightarrow \text{则方程式为： }A{x}^2+Bx+C=0$$

---

对于如下形式的一元二次方程：

$$a{x}^2+bx+c=0$$

方程的两个根为：

$$x_1=\frac{-b+\sqrt{b^2-4ac}}{2a}\text{；}{x}_2=\frac{-b-\sqrt{b^2-4ac}}{2a}$$

$$\text{一元二次方程求根公式推导过程：}\{\begin{array}{l}a{x}^2+bx+c=0\\ \Rightarrow x^2+\frac{b}{a}x+\frac{c}{a}=0\\ \Rightarrow x^2+\frac{b}{a}x=-\frac{c}{a}\\ \Rightarrow x^2+2\frac{b}{2a}x=-\frac{c}{a}\\ \Rightarrow x^2+2\frac{b}{2a}x+{\left(\frac{b}{2a}\right)}^2=-\frac{c}{a}+{\left(\frac{b}{2a}\right)}^2\\ \Rightarrow x^2+2\frac{b}{2a}x+\frac{b^2}{4a^2}=-\frac{c}{a}+\frac{b^2}{4a^2}\\ \Rightarrow {\left(x+\frac{b}{2a}\right)}^2=-\frac{c}{a}+\frac{b^2}{4a^2}\\ \Rightarrow {\left(x+\frac{b}{2a}\right)}^2=-\frac{4ac}{4a^2}+\frac{b^2}{4a^2}\\ \Rightarrow {\left(x+\frac{b}{2a}\right)}^2=\frac{b^2-4ac}{4a^2}\\ \Rightarrow x+\frac{b}{2a}=\frac{\pm \sqrt{\left(b^2-4ac\right)}}{2a}\\ \Rightarrow x=\frac{-b\pm \sqrt{\left(b^2-4ac\right)}}{2a}\end{array}$$

---

以下为通用椭圆方程与通用直线方程的交点解析：
参考：标准椭圆和任意椭圆方程之间的变换公式推导
$$kx+my+n=0通用直线方程$$
将直线方程化简：
$$y=-\frac{k}{m}x-\frac{n}{m}$$

$$\Rightarrow \text{设}\{\begin{array}{l}K=-\frac{k}{m}\\ M=-\frac{n}{m}\end{array}$$
将K，M带入：
$$\Rightarrow 得：y=Kx-M$$

---

我们在高中数学中就学习过标准的椭圆方程如下：

$$\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$$

这个方程表示一个中心在坐标原点、长轴沿$x$轴d 长度为$a$，短轴沿$y$轴的轴长度为$b$的正椭圆（$a$ > $b$ > 0 的条件下），如下图所示：

一般来说，椭圆可以以任何一点为中心，也可以有与坐标轴不平行的轴。这样的椭圆总是可以从标准位置的椭圆开始，然后进行旋转和/或平移得到。对于一般性的公式，我们可以包括通过一个角度为0的旋转(即根本没有旋转)和通过零向量的平移(根本没有平移)来进行变换。也就是说，每一个椭圆都可以通过在标准位置上旋转和平移得到。因此，对椭圆的标准方程进行旋转和平移，可以得到任意椭圆的方程。

是旋转然后平移，还是相反，这是一个需要选择的问题。为了理解这个，让R表示一个旋转，考虑点$x=(x,y)$，会发生什么如果我们先平移向量$v$，然后应用 $R$。因为$R$是线性的，变换之后的结果为$R(x+v)=Rx+Rv$，然而这和先旋转$x$，然后用$Rv$平移是一样的。这表明，每一个椭圆都可以从一个标准位置的椭圆中得到，要么是旋转后再平移，要么是平移后再旋转。在推导椭圆的一般方程时，我们将使用先旋转后平移的方法。

参考旋转矩阵的两种用法一文，我们使用极坐标工具，绕原点逆时针旋转，通过一个角$\alpha$，很容易可以将$(x,y)$变成$(xcos\alpha -ysin\alpha ,ycos\alpha +xsin\alpha )$。写成矩阵的形式为：
$$\left[\begin{array}{c}u\\ v\end{array}\right]=\left[\begin{array}{cc}cos\alpha & -sin\alpha \\ sin\alpha & cos\alpha \end{array}\right]\left[\begin{array}{c}x\\ y\end{array}\right]$$
它的逆运算可以通过旋转$2\pi -\alpha$得到，因此将$(u,v)$变换为$(xcos\alpha +ysin\alpha ,ycos\alpha -xsin\alpha )$。即：
$$\left[\begin{array}{c}u\\ v\end{array}\right]=\left[\begin{array}{cc}cos\alpha & sin\alpha \\ -sin\alpha & cos\alpha \end{array}\right]\left[\begin{array}{c}x\\ y\end{array}\right]$$
事实上，这表明了一个定理，即对于旋转矩阵：
$$R^{-1}=R^T$$

应用变换椭圆方程的方法，可以得到以下方程，对于一个经过 $\theta$角旋转的标准椭圆：
$$\frac{{\left(xcos\theta +ysin\theta \right)}^2}{a^2}+\frac{{\left(ycos\theta -xsin\theta \right)}^2}{b^2}=1$$

$$\frac{x^{2}co{s}^2\theta +2xycos\theta sin\theta +y^{2}si{n}^2\theta }{a^2}+\frac{y^{2}co{s}^2\theta -2xysin\theta cos\theta +x^{2}si{n}^2\theta }{b^2}=1$$

$$\frac{co{s}^2\theta }{a^2}{x}^2+\frac{si{n}^2\theta }{b^2}{x}^2+\frac{2cos\theta sin\theta }{a^2}xy-\frac{2sin\theta cos\theta }{b^2}xy+\frac{si{n}^2\theta }{a^2}{y}^2+\frac{co{s}^2\theta }{b^2}{y}^2=1$$

$$\left(\frac{co{s}^2\theta }{a^2}+\frac{si{n}^2\theta }{b^2}\right)x^2+2sin\theta cos\theta \left(\frac{1}{a^2}-\frac{1}{b^2}\right)xy+\left(\frac{si{n}^2\theta }{a^2}+\frac{co{s}^2\theta }{b^2}\right)y^2=1$$

$$\left(\frac{co{s}^2\theta }{a^2}+\frac{si{n}^2\theta }{b^2}\right)x^2+2sin\theta cos\theta \left(\frac{1}{a^2}-\frac{1}{b^2}\right)xy+\left(\frac{si{n}^2\theta }{a^2}+\frac{co{s}^2\theta }{b^2}\right)y^2=1$$

$$\Rightarrow \text{设}\{\begin{array}{l}A=\frac{co{s}^2\theta }{a^2}+\frac{si{n}^2\theta }{b^2}\\ B=2sin\theta cos\theta \left(\frac{1}{a^2}-\frac{1}{b^2}\right)\\ C=\frac{si{n}^2\theta }{a^2}+\frac{co{s}^2\theta }{b^2}\end{array}$$

则中心点在原点$(0,0)$做任意旋转的椭圆方程式可简化为：
$$A{x}^2+Bxy+C{y}^2=1$$

---

$$A-C可得：\left(\frac{co{s}^2\theta }{a^2}+\frac{si{n}^2\theta }{b^2}\right)-\left(\frac{si{n}^2\theta }{a^2}+\frac{co{s}^2\theta }{b^2}\right)$$

$$=\frac{b^{2}co{s}^2\theta }{a^2{b}^2}+\frac{a^{2}si{n}^2\theta }{a^2{b}^2}-\frac{b^{2}si{n}^2\theta }{a^2{b}^2}-\frac{a^{2}co{s}^2\theta }{a^2{b}^2}$$

$$=\frac{b^{2}co{s}^2\theta }{a^2{b}^2}-\frac{a^{2}co{s}^2\theta }{a^2{b}^2}+\frac{a^{2}si{n}^2\theta }{a^2{b}^2}-\frac{b^{2}si{n}^2\theta }{a^2{b}^2}$$

$$=\left(\frac{b^{2}co{s}^2\theta }{a^2{b}^2}-\frac{a^{2}co{s}^2\theta }{a^2{b}^2}\right)-\left(\frac{b^{2}si{n}^2\theta }{a^2{b}^2}-\frac{a^{2}si{n}^2\theta }{a^2{b}^2}\right)$$

$$=\frac{co{s}^2\theta \left(b^2-a^2\right)-si{n}^2\theta \left(b^2-a^2\right)}{a^2{b}^2}$$

$$=\frac{\left(co{s}^2\theta -si{n}^2\theta \right)\left(b^2-a^2\right)}{a^2{b}^2}$$

$$B通分后可得：\frac{2sin\theta cos\theta \left(b^2-a^2\right)}{a^2{b}^2}$$

$$\frac{A-C}{B}为：\frac{\left(co{s}^2\theta -si{n}^2\theta \right)\left(b^2-a^2\right)}{a^2{b}^2}\times \frac{a^2{b}^2}{2sin\theta cos\theta \left(b^2-a^2\right)}=\frac{co{s}^2\theta -si{n}^2\theta }{2sin\theta cos\theta }$$

根据三角函数二倍角公式：
$$\{\begin{array}{l}\sin 2\theta =2sin\theta cos\theta \\ \cos 2\theta =co{s}^2\theta -si{n}^2\theta \end{array}$$
可得：
$$cot\left(2\theta \right)=\frac{A-C}{B}$$

然后，再次使用坐标旋转变换公式，即可得到椭圆的长短轴信息，参考Rotated Conic Section Identifying这个视频；

---

通过上述求证，则可以得到，将椭圆中心点平移到$\left(x^{\prime },y^{\prime }\right)$，并且以$\left(x^{\prime },y^{\prime }\right)$为中心的旋转椭圆的一般方程为：
$$A{\left(x-x^{\prime }\right)}^2+B\left(x-x^{\prime }\right)\left(y-y^{\prime }\right)+C{\left(y-y^{\prime }\right)}^2=1$$
将方程式进一步拆解：

$$A\left(x^2-2x{x}^{\prime }+x^{\prime 2}\right)+B\left(xy-x{y}^{\prime }-x^{\prime }y+x^{\prime }{y}^{\prime }\right)+C\left(y^2-2y{y}^{\prime }+y^{\prime 2}\right)=1$$

$$A{x}^2-2Ax{x}^{\prime }+A{x}^{\prime 2}+Bxy-Bx{y}^{\prime }-B{x}^{\prime }y+B{x}^{\prime }{y}^{\prime }+C{y}^2-2Cy{y}^{\prime }+C{y}^{\prime 2}=1$$

$$A{x}^2+Bxy+C{y}^2-\left(2A{x}^{\prime }+B{y}^{\prime }\right)x-\left(2C{y}^{\prime }+B{x}^{\prime }\right)y+\left(A{x}^{\prime 2}+B{x}^{\prime }{y}^{\prime }+C{y}^{\prime 2}\right)=1$$
最终椭圆的通用方程式为：
$$x^2+\frac{B}{A}xy+\frac{C}{A}{y}^2-\left(2x^{\prime }+\frac{B}{A}{y}^{\prime }\right)x-\left(2\frac{C}{A}{y}^{\prime }+\frac{B}{A}{x}^{\prime }\right)y+\left(x^{\prime 2}+\frac{B}{A}{x}^{\prime }{y}^{\prime }+\frac{C}{A}{y}^{\prime 2}-\frac{1}{A}\right)=0$$

---

将直线方程式与通用椭圆方程式联立求解：
$$\{\begin{array}{l}y=Kx+M\\ x^2+\frac{B}{A}xy+\frac{C}{A}{y}^2-\left(2x^{\prime }+\frac{B}{A}{y}^{\prime }\right)x-\left(2\frac{C}{A}{y}^{\prime }+\frac{B}{A}{x}^{\prime }\right)y+\left(x^{\prime 2}+\frac{B}{A}{x}^{\prime }{y}^{\prime }+\frac{C}{A}{y}^{\prime 2}-\frac{1}{A}\right)=0\end{array}$$

$$x^2+\frac{B}{A}x\left(Kx+M\right)+\frac{C}{A}{\left(Kx+M\right)}^2-\left(2x^{\prime }+\frac{B}{A}{y}^{\prime }\right)x-\left(2\frac{C}{A}{y}^{\prime }+\frac{B}{A}{x}^{\prime }\right)\left(Kx+M\right)+\left(x^{\prime 2}+\frac{B}{A}{x}^{\prime }{y}^{\prime }+\frac{C}{A}{y}^{\prime 2}-\frac{1}{A}\right)=0$$

$$x^2+\frac{B}{A}K{x}^2+\frac{B}{A}Mx+\frac{C}{A}\left(K^2{x}^2+2KMx+M^2\right)-\left(2x^{\prime }+\frac{B}{A}{y}^{\prime }\right)x-\left(2\frac{C}{A}Kx{y}^{\prime }+2\frac{C}{A}M{y}^{\prime }+\frac{B}{A}Kx{x}^{\prime }+\frac{B}{A}M{x}^{\prime }\right)+\left(x^{\prime 2}+\frac{B}{A}{x}^{\prime }{y}^{\prime }+\frac{C}{A}{y}^{\prime 2}-\frac{1}{A}\right)=0$$

$$x^2+\frac{B}{A}K{x}^2+\frac{B}{A}Mx+\frac{C}{A}{K}^2{x}^2+2\frac{C}{A}KMx+\frac{C}{A}{M}^2-\left(2x^{\prime }+\frac{B}{A}{y}^{\prime }\right)x-2\frac{C}{A}K{y}^{\prime }x-2\frac{C}{A}M{y}^{\prime }-\frac{B}{A}K{x}^{\prime }x-\frac{B}{A}M{x}^{\prime }+\left(x^{\prime 2}+\frac{B}{A}{x}^{\prime }{y}^{\prime }+\frac{C}{A}{y}^{\prime 2}-\frac{1}{A}\right)=0$$

$$\left(1+\frac{B}{A}K+\frac{C}{A}{K}^2\right)x^2+\left(\frac{B}{A}M+2\frac{C}{A}KM-2x^{\prime }-\frac{B}{A}{y}^{\prime }-2\frac{C}{A}K{y}^{\prime }-\frac{B}{A}K{x}^{\prime }\right)x+\left(\frac{C}{A}{M}^2-2\frac{C}{A}M{y}^{\prime }-\frac{B}{A}M{x}^{\prime }+x^{\prime 2}+\frac{B}{A}{x}^{\prime }{y}^{\prime }+\frac{C}{A}{y}^{\prime 2}-\frac{1}{A}\right)=0$$

最终得到一元二次方程式为：
$$\left(A+BK+C{K}^2\right)x^2+\left(BM+2CKM-2A{x}^{\prime }-B{y}^{\prime }-2CK{y}^{\prime }-BK{x}^{\prime }\right)x+\left(C{M}^2-2CM{y}^{\prime }-BM{x}^{\prime }+A{x}^{\prime 2}+B{x}^{\prime }{y}^{\prime }+C{y}^{\prime 2}-1\right)=0$$

---

综合前期
$$\Rightarrow \text{记}\{\begin{array}{l}K=-\frac{k}{m}\\ M=-\frac{n}{m}\end{array}$$

$$\Rightarrow \text{记}\{\begin{array}{l}A=\frac{co{s}^2\theta }{a^2}+\frac{si{n}^2\theta }{b^2}\\ B=2sin\theta cos\theta \left(\frac{1}{a^2}-\frac{1}{b^2}\right)\\ C=\frac{si{n}^2\theta }{a^2}+\frac{co{s}^2\theta }{b^2}\end{array}$$

椭圆中心点坐标$\left(x^{\prime },y^{\prime }\right)$

以上所述以及结合韦达定理，可得直线与任意椭圆交点：

$$韦达定理\{\begin{array}{l}a{x}^2+bx+c=0\left(a,b,c\in R,a\ne 0\right)\\ x_1=\frac{-b+\sqrt{b^2-4ac}}{2a}\text{；}{x}_2=\frac{-b-\sqrt{b^2-4ac}}{2a}\\ \Delta =b^2-4ac\\ \{\begin{array}{ll}\Delta >0& \text{方程有两个不相等的实根}\\ \Delta =0& \text{方程有两个相等的实根}\\ \Delta <0& \text{方程有两个复根}\end{array}\\ x_1+x_2=\frac{-b+\sqrt{b^2-4ac}}{2a}+\frac{-b-\sqrt{b^2-4ac}}{2a}=-\frac{b}{a}\\ x_1\cdot x_2=\frac{-b+\sqrt{b^2-4ac}}{2a}\times \frac{-b-\sqrt{b^2-4ac}}{2a}=\frac{c}{a}\end{array}$$

2021-7-10 完结！