- 注册时间
- 2009-2-12
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 24148
- 在线时间
- 小时
|
发表于 2020-8-12 09:10:58
|
显示全部楼层
圆锥曲线族有一个定义就是 到定点的距离与到定直线的距离的比值为定值。所以拿这个定义 就能一举证明所有的情况
设 定直线为准线$x=-p = -a^2/c$,那么定点是焦点,接下来我们用离心率$e=c/a$和$p=a^2/c$来表达方程的标准形式(中心在原点,对称轴是坐标轴),
对于椭圆,双曲线,焦点就是就是$(-pe^2, 0)$,方程是 $(x +p e^2)^2 + y^2 - e^2 (x+ p)^2=0$
对于抛物线,焦点就是$(p,0)$,方程就是$(x-p)^2+y^2-(x+p)^2=0$
##############################################
也可以设准线$x=p$,那么椭圆双曲线的焦点就是$(pe^2, 0)$,但是抛物线的焦点就是$(2p,0)$,这个时候的抛物线不是标准的开口向右的方程了,准线设成$x=-p$,那么就能直接得到抛物线方程$ y^2=4 p x$
##############################################
三个点的坐标用极坐标表达,幅角顺次增加$\frac{2\pi}{3}$,消元
- f[x_,y_]:=(x - p e^2)^2 + y^2 - e^2 (x - p)^2;
- eq=Eliminate[{f[x+r Cos[\[Theta]],y+r Sin[\[Theta]]]==0,f[x+r Cos[\[Theta]+2Pi/3],y+r Sin[\[Theta]+2Pi/3]]==0,f[x+r Cos[\[Theta]-2Pi/3],y+r Sin[\[Theta]-2Pi/3]]==0},r]/.And->List;
- ans=Table[{i,TrigFactor[(eq[[i]]/.Equal->List).{1,-1}]},{i,Length[eq]}];
复制代码
这个时候消元 消不动,请教chyanog之后,把三角函数方程换成有理表达,并且用Groebner Basis消元,就可以瞬间解决掉了
- f[x_,y_]:=(x-p)^2+y^2-(x+p)^2;
- f[x_,y_]:=x^2+y^2-e^2(x+P)^2;
- f[x_,y_]:=(x+p e^2)^2+y^2- e^2(x+p)^2;
- n=3;Table[f[x+r Cos[\[Theta]+(2k \[Pi])/n],y+r Sin[\[Theta]+(2k \[Pi])/n]]==0,{k,n}]/.\[Theta]->2ArcTan[t]//TrigExpand//Simplify
- GroebnerBasis[%,{},{r,t},MonomialOrder->EliminationOrder]//Factor
复制代码
得到方程$-e^6 p^2+e^8 p^2+16 x^2-40 e^2 x^2+33 e^4 x^2-9 e^6 x^2+16 y^2-8 e^2 y^2+e^4 y^2$,增根是原曲线$-e^2 p^2+e^4 p^2+x^2-e^2 x^2+y^2=0$,舍去, 稍微调整一下就是 $(1-e^2) (3 e^2-4)^2 (x/p)^2+ (e^2-4)^2 (y/p)^2 = e^6 (1-e^2) $,
- Collect[-e^6 p^2 + e^8 p^2 + 16 x^2 - 40 e^2 x^2 + 33 e^4 x^2 - 9 e^6 x^2 + 16 y^2 - 8 e^2 y^2 + e^4 y^2, {x, y}, Factor]
复制代码
对于抛物线,得到方程$32 p^2 - 4 p x + 9 y^2=0$,增根是原曲线$4 p x - y^2=0$,舍去。
##############################################
还有一种是非标准方程的设定,好处是一个参数设定全部囊括了椭圆双曲线和抛物线,得到一个大一统的方程,直接设焦点是原点$(0,0)$,准线是$x=-P$,那么圆锥曲线方程就是$x ^2 + y^2 - e^2 (x + P)^2 = 0$ ,得到,
\[9 e^6 P^2-16 e^4 P^2+18 e^6 P x-48 e^4 P x+32 e^2 P x+9 e^6 x^2-33 e^4 x^2+40 e^2 x^2-e^4 y^2+8 e^2 y^2-16 x^2-16 y^2=0\]
经过hujunhua提醒,配方会失去抛物线,所以,配方前先分析一下:
对于抛物线,代进$e=1$,得到轨迹方程是 $-7 P^2+2 P x-9 y^2=0$
- a (x-b)^2+c+d y^2/.Factor[SolveAlways[-16 e^4 P^2+9 e^6 P^2+32 e^2 P x-48 e^4 P x+18 e^6 P x-16 x^2+40 e^2 x^2-33 e^4 x^2+9 e^6 x^2-16 y^2+8 e^2 y^2-e^4 y^2==a (x-b)^2+c+d y^2&&a!=0,{x,y}]]
复制代码
对于椭圆和双曲线,$e!=1$,对表达式进行配方,得到:
\[\left(3 e^2-4\right)^2 \left(e^2 P+\left(e^2-1\right) x\right)^2-\left(e^2-1\right) \left(e^2-4\right)^2 y^2=e^6 P^2\]
|
评分
-
参与人数 1 | 威望 +12 |
金币 +12 |
贡献 +12 |
经验 +12 |
鲜花 +12 |
收起
理由
|
数学星空
| + 12 |
+ 12 |
+ 12 |
+ 12 |
+ 12 |
赞一个! |
查看全部评分
|