亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        高斯無奇點(diǎn)衛(wèi)星運(yùn)動(dòng)方程的數(shù)學(xué)推導(dǎo)
        ——謹(jǐn)以本文恭賀師兄歐吉坤教授七十壽辰

        2015-05-08 00:37:03許國(guó)昌沈云中蔣春華
        導(dǎo)航定位學(xué)報(bào) 2015年3期
        關(guān)鍵詞:許國(guó)將式拉格朗

        許國(guó)昌,陳 武,沈云中,江 楠,蔣春華

        高斯無奇點(diǎn)衛(wèi)星運(yùn)動(dòng)方程的數(shù)學(xué)推導(dǎo)
        ——謹(jǐn)以本文恭賀師兄歐吉坤教授七十壽辰

        許國(guó)昌1,2,陳 武3,沈云中4,江 楠5,1,蔣春華1

        (1.山東大學(xué) 空間科學(xué)研究院,山東 威海 264209;2.中國(guó)空間技術(shù)研究院 中國(guó)航天錢學(xué)森實(shí)驗(yàn)室,北京 100094; 3.香港理工大學(xué),中國(guó) 香港 999077;4.同濟(jì)大學(xué) 測(cè)繪與地理信息學(xué)院,上海 200092;5.柏林工業(yè)大學(xué),德國(guó) 柏林 10623)

        本文給出高斯無奇點(diǎn)運(yùn)動(dòng)方程的嚴(yán)格數(shù)學(xué)推導(dǎo)。該方程的無奇點(diǎn)特性使軌道力學(xué)中的奇點(diǎn)問題獲得了完全的解決。

        天體力學(xué);解析方法;奇點(diǎn)問題;高斯無奇點(diǎn)運(yùn)動(dòng)方程;拉格朗日無奇點(diǎn)運(yùn)動(dòng)方程

        0 引言

        牛頓在1687年首先在其著作數(shù)學(xué)原理中提出萬有引力和運(yùn)動(dòng)方程[1-5]。其后,天體力學(xué)和N體問題的理論研究主要基于拉格朗日方程和源于拉格朗日方程的高斯方程[6-10]。奇點(diǎn)問題一直是天體力學(xué)和N體問題研究中的核心問題之一[11-16]。人們?yōu)槔碚撗芯窟M(jìn)展作了很多努力[17-24]。衛(wèi)星軌道理論研究主要基于二體問題的攝動(dòng)理論[31-39]。二階攝動(dòng)理論解研究近年來有了系統(tǒng)的進(jìn)展[26-28,30],所以奇點(diǎn)問題的解決成為迫切的需要。無奇點(diǎn)理論2012年以不定積分的解的形式提出[29]并于2013年有了以微分方程描述的版本[25]。但其公式推導(dǎo)主要是邏輯推理,2014年有了拉格朗日無奇點(diǎn)運(yùn)動(dòng)方程純數(shù)學(xué)的嚴(yán)格具體的證明推導(dǎo)[40]。本文給予高斯無奇點(diǎn)運(yùn)動(dòng)方程純數(shù)學(xué)的嚴(yán)格具體的證明推導(dǎo)。

        1 高斯無奇點(diǎn)運(yùn)動(dòng)方程的數(shù)學(xué)推導(dǎo)

        拉格朗日無奇點(diǎn)運(yùn)動(dòng)方程的數(shù)學(xué)推導(dǎo)已在文獻(xiàn)[40]中詳細(xì)給出。高斯運(yùn)動(dòng)方程的推導(dǎo)的基礎(chǔ)是拉格朗日運(yùn)動(dòng)方程,其推導(dǎo)可以在很多文獻(xiàn)[24-25]中找到,為了獲得優(yōu)化的推導(dǎo)公式,直接從拉格朗日無奇點(diǎn)運(yùn)動(dòng)方程開始。

        拉格朗日無奇點(diǎn)運(yùn)動(dòng)方程[40]即Lagrange-Xu運(yùn)動(dòng)方程為

        (1)

        (2)

        (3)

        (4)

        (5)

        (6)

        R對(duì)各軌道根數(shù)的偏導(dǎo)數(shù)為[24]:

        (7)

        (8)

        (9)

        (10)

        (11)

        (12)

        式(7)~式(12)中fa、fh、fr分別表示軌道坐標(biāo)系的攝動(dòng)力分量。軌道參數(shù)之間還有關(guān)系

        (13)

        rcosf=a(cosE-e)

        (14)

        (15)

        式(13)~式(15)中E為偏近點(diǎn)角。下面我們分別將式(7)~式(15)代入無奇點(diǎn)拉格朗日運(yùn)動(dòng)方程式(1)~式(6)來推導(dǎo)高斯無奇點(diǎn)方程。

        將式(12)代入式(1),可得

        (16)

        將式(11)、式(12)代入式(2)得

        (17)

        (18)

        (19)

        否則,

        (20)

        或者表示成

        (21)

        將式(10)、式(11)代入式(4)得:

        (22)

        將式(9)代入式(5),同理

        (23)

        (24)

        (25)

        (26)

        (27)

        (28)

        (29)

        (30)

        (31)

        (32)

        式(27)~式(32)即為無奇點(diǎn)的高斯運(yùn)動(dòng)方程,亦稱無奇點(diǎn)的Gauss-Xu方程。上述推導(dǎo)數(shù)學(xué)上是嚴(yán)格的。

        2 與原始高斯方程的比較

        與如下原始高斯衛(wèi)星運(yùn)動(dòng)方程相比[19-20,32,35,37,40]

        (33)

        很容易注意到原高斯方程僅適用于非奇異情形。而無奇點(diǎn)的高斯方程即Gauss-Xu方程式(27)~式(32)則一般地成立,且不存在奇異的可能性。式(27)與式(33) 的第一個(gè)方程相同,因其不會(huì)奇異。式(28)與式(30)實(shí)際上可以直接從式(33)中的第二個(gè)和第四個(gè)方程通過把奇異因子e和 sini從右側(cè)移到左側(cè)而獲得。

        3 與文獻(xiàn)[25]的比較

        高斯運(yùn)動(dòng)方程可作如下變形:

        (34)

        文獻(xiàn)[25]將高斯方程式(34)改寫成[25,29]

        (35)

        式(35)中,δ及其下標(biāo)表示式(34)的右邊各項(xiàng)。式(35)中的第二和第六個(gè)方程乘以e,第四個(gè)和第五個(gè)方程乘以sini,第三個(gè)方程乘以esini,有

        (36)

        式(36)是無奇點(diǎn)的 (n在第六個(gè)方程中被忽略)。解 (即式 (35) 的不定積分) 可用類似于[21,26-28]的方法獲得且有形式

        (37)

        式(37)中,右側(cè)Δy表示y的不定積分,左側(cè)Δx表示dx/dt的不定積分,x表示開普勒軌道根數(shù)。

        則無奇點(diǎn)的高斯方程的解為

        Δa=Δδa

        Δcosi=-siniΔδi

        (38)

        基本高斯運(yùn)動(dòng)方程被定義為[26,40]

        (39)

        由式(38) 和式(39)可得

        (40)

        同上,分別對(duì)式 (40) 條件做進(jìn)一步的推導(dǎo),據(jù)定義[26]

        (41)

        由式(41) 和 (34),條件if(Δδω2,sini)則等價(jià)于if(δω2,sini),

        (42)

        故if(Δδω2,e)的等價(jià)條件為:if(fh,sini)。

        條件if(Δδω1,e)則等價(jià)于if(δω1,e),

        (43)

        故if(Δδω1,e)的等價(jià)條件為:if(fa,e)∧if(fr,e)。

        若if(ΔδΩ,sini)則等價(jià)于if(δΩ,sini),

        (44)

        故if(ΔδΩ,sini)的等價(jià)條件為:if(fh,sini)。

        若if(ΔδM2,e)則等價(jià)于if(δM2,e)

        (45)

        故if(ΔδM2,e)的等價(jià)條件為:if(fa,e)∧if(fr,e)

        則式(40)可以寫成

        (46)

        式(46)與高斯無奇點(diǎn)運(yùn)動(dòng)方程式(27)到式(32)除微小差別(式(46)中最后一式n被忽略)外,完全一致。

        4 結(jié)束語

        衛(wèi)星軌道無奇點(diǎn)運(yùn)動(dòng)方程的正確性通過嚴(yán)格

        的數(shù)學(xué)證明予以證實(shí)。同時(shí)高斯—許方程被以嚴(yán)格的數(shù)學(xué)方法獲得再次推導(dǎo)。

        致謝:本研究受山東大學(xué)和中國(guó)空間技術(shù)研究院支持并部分地獲得山東大學(xué)學(xué)科建設(shè)經(jīng)費(fèi)和國(guó)家自然科學(xué)基金的資助。第二、三作者參與了本文有關(guān)的有價(jià)值的討論。第四、五作者分別就理論公式做了推導(dǎo)和驗(yàn)證。

        [1] NEWTON I.Philosophiae Naturalis Principia Mathematica[M].London:[s.n.],1687.

        [2] GOLDSTEIN H.Classical Mechanics[M].2nd ed.New York:Addison-Wesley,1980.

        [3] LANDAU L D,LIFSHITZ E M.MECHANICS[M].3rd ed.New York:Pergamon Press,1976.

        [4] EULER L.Proprietates Triangulorum,Quorum Anguli Certam inter se Rationem Tenent[J].Novi Commentarii Academiae Scientiarum Imperialis Petropolitanae,1767,10:207-242.

        [5] LAGRANGE J L.Motion for Three Bodies[J].Miscellanea Taurinensia,1772,4:118-243.

        [6] TISSERAND F F.Mecanique Celeste,Tome III[M].Paris:[s.n.],1894:27.

        [7] HAGIHARA Y.Celestial Mechanics[M].[ S.l.]:MIT Press,1970.

        [8] MITTAG-LEFFLER G.The N-body Problem (Price Announcement)[M].[ S.l.]:Acta Matematica,1885/1886:7.

        [9] POINCARE H.New Methods of Celestial Mechanics[M].[ S.l.]:AIP,1992.

        [10] BARROW-GREEN J.Poincare and the Three Body Problem[M].[ S.l.]:American Mathematical Society,1996.

        [11] XIA Zhi-ling.The Existence of Noncollision Singularities in Newtonian System[J].Annals Math.1992,135(3) 411-468.

        [12] SUNDMAN K F.Mémoire sur le Problème des Trois Corps[J].Acta Mathematica,1913,36(1):105-179.

        [13] DIACU F.The Solution of the N-body Problem[J].The Mathematical Intelligencer.1996,18(3):66-70.

        [14] DIACU F,HOLMES P.Celestial Encounters:The Origins or Chaos and Stability[M].NJ:Princeton University Press,1996.

        [15] WANG Qing-de.The Global Solution of The N-body Problem[J].Celestial Mechanics and Dynamical Astronomy,1991,50(1):73-88.

        [16] DIACU F.Singularities of the N-body Problem[M].Montreal:Les Publications CRM,1992.

        [17] AARSETH S J,TOUT C A,MARDLING R A.The Cambridge N-body Lectures[M].Berlin:Springer Heidelberg,2008.

        [18] HAVEL K.N-body Gravitational Problem:Unrestricted Solution[M].Brampton:Grevyt Press,2008.

        [19] DINO B,GIUSEPPE P.Theory of Orbits:Volume 1:Integrable Systems and Non-perturbative Methods[M].Berlin:Springer-Verlag,1998.

        [20] BATTIN R H.An Introduction to the Mathematics and Methods of Astrodynamics[M/OL].[2015-04-20].http://www.gravityassist.com/IAF-3.3%202010/Ref.%203-227.pdf.

        [21] KAULA W.Theory of Satellite Geodesy[M].New York:Blaisdell Publishing Company,Dover Publications,2000.

        [22] 劉林,趙德滋.人造地球衛(wèi)星運(yùn)動(dòng)理論[M].南京:南京大學(xué)出版社,1979.

        [23] AARSETH S J.Gravitational N-body Simulations Tools and Algorithms[M].London:Cambridge University Press,2003.

        [24] XU Guo-chang.Orbits[M].Berlin:Springer Heidelberg,2008.

        [25] XU Guo-chang,XU Jia.Orbits-2nd Order Singularity-Free Solutions[M].2nd ed.Berlin:Springer Heidelberg,2013.

        [26] XU Guo-chang,XU Tian-he,CHEN Wu,et al.Analytic Solution of Satellite Orbit Disturbed by Lunar and Solar Gravitation[J].Monthly Notices of the Royal Astronomical Society,2010,410(1):645-653.

        [27] XU Guo-chang,XU Tian-he,YEH T K,et al.Analytic Solution of Satellite Orbit Disturbed by Atmospheric Drag[J].Monthly Notices of the Royal Astronomical Society,2010,410(1):654-662.

        [28] XU Yan,YANG Yuan-xi,ZHANG Qin,et al.Solar Oblateness and Mercury’s Perihelion Precession[J].Monthly Notices of the Royal Astronomical Society,2011,415(4):3335-3343.

        [29] XU Guo-chang,XU Jia.On The Singularity Problem in Orbital Mechanics[J].Monthly Notices of the Royal Astronomical Society,2013,429(2):1139-1148.

        [30] XU Guo-chang,XU Jia.On Orbital Disturbing Effects of The Solar Radiation[J].Monthly Notices of the Royal Astronomical Society,2013,432(1):584-588.

        [31] XU Guo-chang.Analytic Orbit Theory[M].Berlin:Springer Heidelberg,2010.

        [32] BROUWER D,CLEMENCE G M.Methods of Celestial Mechanics[M].New York:Academic Pres,1961.

        [33] CHOBOTOV V A.Orbital Mechanics[M].Washington,D.C:AIAA,1991.

        [34] CUI Chun-fang.Die Bewegung künstlicher Satelliten im anisotropen Gravitationsfeld einer gleichm?ssig rotierenden starren Modellerde:-eine analytische L?sung 2[M].München:Verlag der Bayer Akademie der Wissenschaften,1990.

        [35] HERRICK S.Astrodynamics[M].London:Van Nostrand Reinhold,1972.

        [36] SHAPIRO I I.The Prediction of Satellite Orbits[C]// ROY M.Dynamics of Satellites/Dynamique des Satellites.Paris:IUTAM Symposium,1962.

        [37] VALLADO D A.Fundamentals of Astrodynamics and Applications[M].New York:Springer-Verlag,2007.

        [38] WNUK E.Tesseral Harmonic Perturbations in the Keplerian Orbital Elements[J].Acta,Astronomica,1990,40(1/2):191-198.

        [39] DVORAK R,LHOTKA C.Celestial Dynamics-Chaoticity and Dynamics of Celestial Systems[M].Weinheim:Wiley,2013.

        [40] 許國(guó)昌,呂志平,沈云中.拉格朗日無奇點(diǎn)衛(wèi)星運(yùn)動(dòng)方程的數(shù)學(xué)推導(dǎo)[C]//孫和平,熊熊,王勇.大地測(cè)量與地球動(dòng)力學(xué)進(jìn)展(第2輯).武漢:湖北科學(xué)技術(shù)出版社,2014:1134-1146.

        A Mathematical Derivation of Singularity-free Gaussian Equations of Planetary Motion

        XUGuo-chang1,2,CHENWu3,SHENYun-zhong4,JIANGNan5,1,JIANGChun-hua1

        (1.Institute of Space of Sciences Shandong University,Weihai 264209,China; 2.Qian Xuesen Laboratory,China Academy of Space Technology,Beijing 100094,China; 3.The Hong Kong Polytechnic University,Hong Kong 999077,China; 4.College of Surveying and Geo-informatics,Tongji University,Shanghai 200092,China; 5.Technical University Berlin,Berlin 10623,Germany)

        An exact mathematical derivation is given for the singularity-free Gaussian equations of planetary motion.The singularity-free property of the derived equations solved the singularity problem of the orbital theory completely.

        celestial mechanics;analytical method;singularity problem;singularity-free Gaussian equations of motion;singularity-free Lagrange equations of motion

        許國(guó)昌,陳武, 沈云中,等.高斯無奇點(diǎn)衛(wèi)星運(yùn)動(dòng)方程的數(shù)學(xué)推導(dǎo)[J].導(dǎo)航定位學(xué)報(bào),2015,3(3):05-12.(XU Guo-chang,CHEN Wu,SHEN Yun-zhong,et al.A Mathematical Derivation of Singularity-free Gaussian Equations of Planetary Motion[J].Journal of Navigation and Positioning,2015,3(3):05-12.)

        10.16547/j.cnki.10-1096.20150302.

        2015-05-18

        國(guó)家自然科學(xué)基金(41274042)。

        許國(guó)昌(1953—),福建莆田人,博士,教授;山東大學(xué)國(guó)家特聘教授;長(zhǎng)安大學(xué)、解放軍信息工程大學(xué)兼職教授;西南交通大學(xué)榮譽(yù)教授;中國(guó)空間技術(shù)研究院國(guó)家千人特聘專家;中國(guó)航天集團(tuán)錢學(xué)森實(shí)驗(yàn)室首席科學(xué)家;德波茨坦地球科學(xué)中心資深研究員(1993年~2014年);Springer Sciences of Geodesy主編;教育部長(zhǎng)江學(xué)者計(jì)劃海外通訊評(píng)審專家;Porto大學(xué)、柏林工業(yè)大學(xué)、德波茨坦地球科學(xué)中心及香港理工大學(xué)兼職博導(dǎo)。

        P236

        A

        2095-4999(2015)-03-0005-08

        猜你喜歡
        許國(guó)將式拉格朗
        AKNS方程的三線性型及周期孤立波解
        因子von Neumann代數(shù)上非線性*-Lie導(dǎo)子的刻畫
        精準(zhǔn)扶貧
        鄧稼先 此身許國(guó) 再難許卿
        海峽姐妹(2019年10期)2019-11-23 08:20:26
        單自由度系統(tǒng)
        Nearly Kaehler流形S3×S3上的切觸拉格朗日子流形
        拉格朗日代數(shù)方程求解中的置換思想
        基于拉格朗日的IGS精密星歷和鐘差插值分析
        阻尼系統(tǒng)的特征
        拉格朗日點(diǎn)
        太空探索(2014年3期)2014-07-10 14:59:39
        日本一区二区在线播放| 国产v片在线播放免费无码| 加勒比一本heyzo高清视频| 一本色道久久99一综合| 国产精品高潮呻吟av久久无吗| 97精品国产高清自在线看超 | 伊人久久中文大香线蕉综合| 呦泬泬精品导航| 午夜精品人妻中字字幕| 亚洲福利二区三区四区| 天堂一区二区三区在线观看视频| 小鲜肉自慰网站| 婷婷综合缴情亚洲| 国产美女在线一区二区三区| 亚洲视频精品一区二区三区 | 亚洲 欧美 国产 制服 动漫| 无码午夜成人1000部免费视频| 中文字幕久无码免费久久| 蜜臀aⅴ永久无码一区二区| 白白色福利视频在线观看| 本道天堂成在人线av无码免费| 亚洲色成人网站www永久| 国产精品久久久av久久久 | 激情内射亚洲一区二区| 国产一区二区三区四区在线视频| 国产精品久久久亚洲| 小12箩利洗澡无码视频网站| 久久国产品野战| 日本女优禁断视频中文字幕| 成年免费a级毛片免费看无码| 熟妇人妻无码中文字幕| 永久免费的拍拍拍网站| 青青草成人在线播放视频 | 丰满五十六十老熟女hd| 视频精品熟女一区二区三区| 日韩三级一区二区三区| 日本50岁丰满熟妇xxxx| 国产目拍亚洲精品一区二区| 国产av无毛无遮挡网站| 亚洲中字慕日产2020| 国产精品乱一区二区三区|