從空間維數(shù)的本質(zhì)探討多體問題下衛(wèi)星實(shí)際運(yùn)動的軌道

周新力

(湖南邵陽學(xué)院 城建系，湖南 邵陽 422000)

1 引言

迄今為止，國內(nèi)外對衛(wèi)星運(yùn)動的基本規(guī)律的研究都主要是基于物理學(xué)中牛頓的萬有引力定律和愛因斯坦的相對論的理論框架所展開的[1-7]，并且對于衛(wèi)星繞地球的運(yùn)行狀態(tài)也主要是歸因于其受諸多引力作用的結(jié)果。在此基礎(chǔ)上，對于只考慮地心引力作用的“二體問題”，已基本得到衛(wèi)星運(yùn)動的嚴(yán)密分析解；而對于同時考慮其他天體引力作用的“多體問題”，卻至今還未得到有效地解決[8-12]。當(dāng)然，最近幾年來國內(nèi)外一些專家和學(xué)者就“多體問題”的算法[13-16]和衛(wèi)星軌道的測控技術(shù)[17-20]這兩方面提出了比較好的處理方案。但不管是“二體問題”還是“多體問題”下的衛(wèi)星運(yùn)動，其得出的解所描述的只能是一個更接近于衛(wèi)星實(shí)際軌道的運(yùn)動軌道。本文則是從有別于現(xiàn)有理論框架的另一個新的角度即空間幾何的本質(zhì)來推導(dǎo)的衛(wèi)星運(yùn)動的實(shí)際的而非近似的軌道方程，供學(xué)界參考。

2 問題的提出

從空間的幾何性質(zhì)看，假設(shè)朝1維平直方向無限延伸的空間點(diǎn)集為N(N=1,2,…,∞)，那么，同時朝2維或3維平直方向無限延伸的空間點(diǎn)集應(yīng)為N2或N3。相比較可知，盡管在1維、2維和3維等平直空間里的點(diǎn)集都是無限的，但N2的無窮大量總是要比N高一階，而N3則又是比N2更高一階的無窮大量。這意味著，假如數(shù)量為N2的空間點(diǎn)集不能向2維平直空間無限延伸而只能選擇朝某1維平直空間無限延伸，如圖1(a)所示，該空間的點(diǎn)集實(shí)質(zhì)上將由N×N變成為1×N2。同樣，假如數(shù)量為N3的空間點(diǎn)集只能選擇朝某2維或1維平直空間無限延伸而不能向3維平直空間無限延伸，如圖1(b)所示，該空間的點(diǎn)集實(shí)質(zhì)上也將由N·N·N變成為1·N·N2或1·1·N3。

圖1 高維空間的點(diǎn)集到低維空間的幾何映射

值得注意的問題是，由于朝任何一個平直的方向無限延伸的宇宙空間其相互垂直的平直的幾何維度最多只有3個。假設(shè)在宇宙空間里建立一個3維直角坐標(biāo)系o-xyz，對于數(shù)量為Nn(n≥4)的空間點(diǎn)集，如圖2(a)所示，盡管它既可以選擇全部都在1維的平直空間o-y內(nèi)無限延伸，又可以選擇部分在o-y之外沿2維平直空間o-yz的z方向分布的其他1維平直空間zi-y(i=1,2,3,…)內(nèi)無限延伸，或者還可以選擇部分在o-yz之外沿3維平直空間o-xyz的x方向分布的其他2維平直空間xi-yz(i=1,2,3,…)內(nèi)無限延伸；但事實(shí)上卻不能進(jìn)一步選擇向o-xyz之外沿4維或4維以上方向分布的另外3維平直空間內(nèi)無限延伸，如圖2(b)所示。這就決定了數(shù)量為Nn(n≥4)的宇宙空間里所有點(diǎn)集在本質(zhì)上為3維的平直的無限延伸的立體幾何空間里分布時，只能重疊在3維空間的各個點(diǎn)上而發(fā)生變化，即由N·N·N·N·…·N變成為Ni·Nj·Nn-i-j·1·…·1(i,j

圖2 4維或4維以上空間與3維空間的幾何性質(zhì)的差異

3 理論依據(jù)的分析

3.1 空間維數(shù)的本質(zhì)

若把靜止于3維直角坐標(biāo)軸上的點(diǎn)集定為0階，則在該3維直角坐標(biāo)系上任意標(biāo)準(zhǔn)時刻t的瞬間某一個區(qū)域空間V內(nèi)的點(diǎn)集實(shí)質(zhì)上應(yīng)包括0-n階。假設(shè)用(x,y,z)t代表空間位置點(diǎn)0階的3維直角坐標(biāo)值，然后用(x(1),y(1),z(1))t、(x(2),y(2),z(2))t、……、(x(n),y(n),z(n))t等分別表示空間位置點(diǎn)的勻速變化、勻加速變化、……、勻變變加速變化的各階3維直角坐標(biāo)值，那么，在任意時刻t該區(qū)域空間V的體積大小實(shí)質(zhì)上應(yīng)為

(1)

3.2 空間距離的本質(zhì)

圖3 空間維數(shù)的幾何本質(zhì)

(2)

從式(2)可以看出，只有當(dāng)A點(diǎn)的空間位置變化的速度、加速度和變加速度等各階坐標(biāo)值在任意時刻均為零時，o點(diǎn)到A點(diǎn)的空間距離soA才是3維直角坐標(biāo)系o-xyz空間里的兩個坐標(biāo)位置點(diǎn)之間的幾何空間距離loA，該值始終是一個固定不變的定值。而當(dāng)A點(diǎn)的任意一階位置變化的坐標(biāo)值在某一時刻t1不為零時，soA實(shí)際上是一段4維或4維以上的空間里兩個點(diǎn)之間的幾何空間距離，該值是隨時間而不斷變化的。宇宙空間內(nèi)任意兩個質(zhì)點(diǎn)，只有當(dāng)兩者相對靜止時，才共同處在一個3維或3維以下的空間里；而當(dāng)它們相對運(yùn)動時，實(shí)質(zhì)上是處在一個4維或4維以上的高維空間里。這也意味著，當(dāng)用電磁波測量o點(diǎn)與A點(diǎn)之間的空間距離時，有一點(diǎn)可以肯定的是，假設(shè)以坐標(biāo)系o-xyz空間的時間為標(biāo)準(zhǔn)，那么，電磁波在t1時刻無論是由o點(diǎn)發(fā)出還是由A點(diǎn)發(fā)出，在t2時刻到達(dá)對方時，它實(shí)際所覆蓋的空間距離的理論值應(yīng)該為soA而絕不是loA。

4 探討的結(jié)果

4.1 衛(wèi)星在太空中的運(yùn)動軌道方程

空間維數(shù)的本質(zhì)和空間距離的本質(zhì)說明，對于宇宙空間里任意一個物體或質(zhì)點(diǎn)，其靜止時的位置和運(yùn)動時的軌道，都可以直接使用3維位置坐標(biāo)和各階位置變化的運(yùn)動坐標(biāo)來正確地描述之。假設(shè)以地球為靜止參照系建立一個以地球的質(zhì)心Q為原點(diǎn)、以地球自轉(zhuǎn)軸為Z軸、固聯(lián)在地球上并隨地球自轉(zhuǎn)的3維直角坐標(biāo)系Q-XYZ，并設(shè)t時刻衛(wèi)星j在該坐標(biāo)系上的位置坐標(biāo)為(xj,yj,zj)t，由于衛(wèi)星和地球、月球等都在同一個地球公轉(zhuǎn)的軌道面上，即它們在太陽系這一階空間的階數(shù)(假設(shè)為n)是一致的，則t時刻該衛(wèi)星的位置變化的速度、加速度和變加速度等各階坐標(biāo)值C(1)～C(n)在Q-XYZ上的3個分量分別為

對于衛(wèi)星在太空中運(yùn)動的實(shí)際軌道實(shí)質(zhì)上可以通過直接使用3維位置坐標(biāo)和各階位置變化的運(yùn)動坐標(biāo)來建立一個類似泰勒展開式的運(yùn)動方程予以正確地描述，即

(3)

其中

(4)

衛(wèi)星中心與地心之間的空間距離的理論值為

(5)

4.2 方程的意義

根據(jù)以上的衛(wèi)星運(yùn)動方程可以得出這樣一個觀點(diǎn)：即假設(shè)衛(wèi)星僅僅受地球引力的作用而作圓周運(yùn)動時的運(yùn)動方程為

x=Ssinωt，y=Scosωt，z=0

若將此方程化為式(3)的形式，即

則該方程意味著，通常所認(rèn)為的“衛(wèi)星圍繞地球作圓周運(yùn)動是地球引力作用的結(jié)果”實(shí)質(zhì)上可以理解為衛(wèi)星的運(yùn)動并不是什么引力的作用，而僅僅是因為它處在一個以地球為中心的高維空間里的某一個點(diǎn)上，該點(diǎn)在以地心為原點(diǎn)的坐標(biāo)系上的空間坐標(biāo)分別為

同樣，對于所謂的“衛(wèi)星在多體問題下的自由運(yùn)動(即無人為作用的運(yùn)動)”，實(shí)質(zhì)上也可以理解為“衛(wèi)星并沒有受到地球、太陽和月亮等諸多引力的作用而運(yùn)動，只是因為它處在地球圍繞太陽公轉(zhuǎn)的這一階高維空間里的某一個點(diǎn)上，該點(diǎn)在以地心為原點(diǎn)的坐標(biāo)系上的空間坐標(biāo)分別為

5 結(jié)束語

(6)

圖4 衛(wèi)星的幾何定軌原理——距離交會法

當(dāng)將式(6)的計算結(jié)果代入式(3)和式(4)的衛(wèi)星運(yùn)動方程后，則進(jìn)一步可以解算出這個衛(wèi)星在任意ti時刻的空間位置3維坐標(biāo)值(x,y,z)ti和空間位置變化的k階3維坐標(biāo)值(x(1),y(1),z(1),x(2),y(2),

[1] 劉林.人造地球衛(wèi)星軌道力學(xué)[M].北京：高等教育出版社，1992.

[2] 愛因斯坦.狹義和廣義相對論淺說[M].楊潤殷，譯.北京：北京大學(xué)出版社，2006.

[3] 吳連大，王昌彬，童傅.人造衛(wèi)星二階攝動理論的半分析半數(shù)值方法[J].天文學(xué)報，1978，19(2)：131-150.

[4] 柳響林，PAVEL D.基于B-spline和正則化算法的低軌衛(wèi)星軌道平滑[J].地球物理學(xué)報，2006，49(1)：99-104.

[5] RUSSEL R P，LARA M.Long-lifetime Lunar Repeat Ground Track Orbits[J].Journal of Guidance，Control and Dynamics，2007，30(4)：982- 993.

[6] 毛悅，宋小勇，賈小林.GEO/IGSO/MEO衛(wèi)星軌道根數(shù)演化分析[J].測繪科學(xué)，2009，34(1)：119-121.

[7] 羅志才，鐘波，寧津生，等.GOCE衛(wèi)星軌道攝動的數(shù)值模擬與分析[J].武漢大學(xué)學(xué)報：信息科學(xué)版，2009，34(7)：757-760.

[8] 馮來平，賈小林，吳顯兵，等.一種基于抗差M估計和動力學(xué)平滑的衛(wèi)星軌道綜合方法[J].中國科學(xué)G輯，2010，40(5)：603-607.

[9] 馮晶瑯，袁建平，陳記爭.不同高度范圍月球衛(wèi)星軌道受攝分析[J].西北工業(yè)大學(xué)學(xué)報，2011，29(2)：194-197.

[10] 童科偉，劉偉，高朝輝，等.月球衛(wèi)星軌道攝動及凍結(jié)軌道研究[J].航天器工程，2012，21(2)：12-17.

[11] 劉大杰，施一民，過靜珺.全球定位系統(tǒng)(GPS)的原理與數(shù)據(jù)處理[M].上海：同濟(jì)大學(xué)出版社，2001：97-102.

[12] MARCHAL C，孫義燧.多體問題中慣量矩的演變[J].中國科學(xué)A輯，1985(6)：557-565.

[13] LUK’YANOVL G L，NASONOVA P，SHIRMIN G I.The Lagrange-Jacobi Equation in the Finite-size Many-body Problem[J].Astronomy Letters，2003，29(9)：635-639.

[14] GUILLOT A.The Painlevé Property for Quasihomogenous Systems and a Many-body Problem in the Plane[J].Communications in Mathematical Physics，2005，260(1)：181-194.

[15] 梁雨婷，譚毅，程楠，等.一種混合算法求解多體問題下行星際轉(zhuǎn)移軌道[J].計算機(jī)與數(shù)字工程，2010，38(4)：15-29.

[16] 李琪剛，柴亞輝，徐煒民，等.多體問題FMM 算法在加速部件FPGA 研究與實(shí)現(xiàn)[J].計算機(jī)工程與設(shè)計，2011，32(10)：3391-3394.

[17] 劉麗麗，文浩，金棟平，等.繩系衛(wèi)星軌道轉(zhuǎn)移的最優(yōu)控制[J].航空學(xué)報，2009，30(2)：332-336.

[18] 李恒年，李濟(jì)生，黃永宣.動力學(xué)補(bǔ)償受控衛(wèi)星軌道確定算法[J].宇航學(xué)報，2010，31(10)：2269-2275.

[19] 耿濤，劉經(jīng)南，趙齊樂，等.星地監(jiān)測網(wǎng)下的北斗導(dǎo)航衛(wèi)星軌道確定[J].測繪學(xué)報，2011(增刊)：46-51.

[20] 陳明，張宇，曹建峰，等.嫦娥二號衛(wèi)星軌道確定與測軌技術(shù)[J].科學(xué)通報，2012，57(9)：689-696.