找回密码
 欢迎注册
查看: 21603|回复: 11

[原创] 动椭圆的运动区域的面积

[复制链接]
发表于 2014-2-22 11:53:34 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
一个给定大小的恒过两点的动椭圆的运动区域的面积
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-3-2 11:28:53 | 显示全部楼层
这个问题计算量很大:

可以设椭圆的一般方程为

\[\frac{\Big((x-s)\sin(\theta)+(y-t)\cos(\theta)\Big)^2}{m^2}+\frac{\Big(-(x-s)\cos(\theta)+(y-t)\sin(\theta)\Big)^2}{n^2}=1\]

由于椭圆过定点\((x_1,y_1),(x_2,y_2)\)代入方程

\[\left\{ \begin{split}\frac{\Big((x_1-s)\sin(\theta)+(y_1-t)\cos(\theta)\Big)^2}{m^2}+\frac{\Big(-(x_1-s)\cos(\theta)+(y_1-t)\sin(\theta)\Big)^2}{n^2}&=1\\
\frac{\Big((x_2-s)\sin(\theta)+(y_2-t)\cos(\theta)\Big)^2}{m^2}+\frac{\Big(-(x_2-s)\cos(\theta)+(y_2-t)\sin(\theta)\Big)^2}{n^2}&=1\end{split} \right.\]

然后作代换\(\sin(\theta)=\dfrac{2r}{1+r^2},\cos(\theta)=\dfrac{1-r^2}{1+r^2}\)

消元\({s,t}\),再对\(r\)求导,再消元\(r\)即可得到椭圆扫过的边界方程

例如取\(m=5,n=3,x_1=1,y_1=1,x_2=-1,y_2=-1\)容易算得

\(10625r^{12}x^4+7650r^{12}x^2y^2+1377r^{12}y^4+20000r^{11}x^4+54400r^{11}x^3y+14400r^{11}x^2y^2+19584r^{11}xy^3+2592r^{11}y^4-114850r^{12}x^2+171900r^{12}xy-96354r^{12}y^2+9350r^{10}x^4+102400r^{10}x^3y+150348r^{10}x^2y^2+36864r^{10}xy^3+27846r^{10}y^4-80000r^{11}x^2-131584r^{11}xy-10368r^{11}y^2-42400r^9x^4+23936r^9x^3y+239808r^9x^2y^2+198016r^9xy^3+44640r^9y^4+19652r^{12}-717516r^{10}x^2+752872r^{10}xy-688972r^{10}y^2+11407r^8x^4-262144r^8x^3y+114750r^8x^2y^2+262144r^8xy^3+168623r^8y^4+110976r^{11}+38528r^9x^2-656896r^9xy-309632r^9y^2+68672r^7x^4-30464r^7x^3y-561024r^7x^2y^2+178432r^7xy^3+173120r^7y^4+326808r^{10}-1164670r^8x^2+2578500r^8xy-2003390r^8y^2+25364r^6x^4+319488r^6x^3y-55896r^6x^2y^2-598016r^6xy^3+284308r^6y^4+464000r^9+118528r^7x^2+523264r^7xy-299264r^7y^2-68672r^5x^4+30464r^5x^3y+561024r^5x^2y^2-178432r^5xy^3-173120r^5y^4+294780r^8-2172584r^6x^2+3995056r^6xy-1772968r^6y^2+11407r^4x^4-262144r^4x^3y+114750r^4x^2y^2+262144r^4xy^3+168623r^4y^4-171264r^7-118528r^5x^2-523264r^5xy+299264r^5y^2+42400r^3x^4-23936r^3x^3y-239808r^3x^2y^2-198016r^3xy^3-44640r^3y^4-24752r^6-1164670r^4x^2+2578500r^4xy-2003390r^4y^2+9350r^2x^4+102400r^2x^3y+150348r^2x^2y^2+36864r^2xy^3+27846r^2y^4+171264r^5-38528r^3x^2+656896r^3xy+309632r^3y^2-20000rx^4-54400rx^3y-14400rx^2y^2-19584rxy^3-2592ry^4+294780r^4-717516r^2x^2+752872r^2xy-688972r^2y^2+10625x^4+7650x^2y^2+1377y^4-464000r^3+80000rx^2+131584rxy+10368ry^2+326808r^2-114850x^2+171900xy-96354y^2-110976r+19652=0\)

然后对\(r\)求导并消元,其实表达式也惊人的长

我们可以画出轨迹图如下:

椭圆过两定点轨迹.gif
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-3-2 15:09:28 | 显示全部楼层
从2#的图上看,星空有增解需要舍弃。

如果椭圆只是靠着两定点滑动,扫掠区域只有上面所画的一半(如下图所示)。
滑动的椭圆.png
星空将扫掠区域的镜像(以两定点连线为镜像线)混进去了。
扫掠区域是关于两定点连线的中垂线轴对称的,星空描绘的区域不全。

点评

呵,是取参数范围太小,重新将参数范围调整了一下,并且椭圆运动区域有两个,当然这两个区域不是连续的而已。  发表于 2014-3-2 19:48

评分

参与人数 1经验 +6 收起 理由
数学星空 + 6

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-3-2 19:49:02 | 显示全部楼层
椭圆过两定点轨迹2.gif
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-3-2 19:51:09 | 显示全部楼层
给定大小是什么意思

点评

算算椭圆x^2/4+y^2/3=1在x轴上滚动的包络线方程  发表于 2014-3-2 20:26
就是椭圆的长轴与短轴是恒定的,即椭圆只能平移和旋转。  发表于 2014-3-2 20:07
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-3-2 21:51:11 | 显示全部楼层
对于椭圆中心在x轴上的运动轨迹问题:

可以设椭圆的方程为

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

且 \[x_0=a {\Large\int}_0^t{\sqrt{1-\frac{a^2-b^2}{a^2}\sin(\theta)^2}\dif\theta}\]

则包络曲线就是对\(t\)求导后,消元\(t\)后的表达式(不能用初等函数表达)

例如:取\(a=5,b=3\),我们可以得到如下轨迹图

椭圆中心在X轴的运动轨迹.gif
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-3-3 00:01:48 | 显示全部楼层
对于椭圆在x轴上的运动轨迹问题:

可以设椭圆的方程为

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

且 \[x_0=a {\Large\int}_0^t{\sqrt{1-\frac{a^2-b^2}{a^2}\sin(\theta)^2}\dif\theta}\]

当 \(x=x_0,y=0\) 时,\(y'=0\)

则包络曲线就是对\(t\)求导后,消元\(t\)后的表达式(不能用初等函数表达)

例如:取\(a=5,\ b=3\),我们可以得到如下轨迹图

椭圆在X轴的运动轨迹(平逆).gif

椭圆在X轴的运动轨迹(平顺).gif
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-3-3 21:40:51 | 显示全部楼层
QQ截图20140303214326.png QQ截图20140303214353.png
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2014-6-23 20:45:55 | 显示全部楼层
sdg.gif
给个动态图!
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-24 04:44 , Processed in 0.028076 second(s), 21 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表