![卫星轨道力学算法](https://wfqqreader-1252317822.image.myqcloud.com/cover/375/33239375/b_33239375.jpg)
2.3 椭圆运动的展开式
在很多问题中,需要将有关量通过平近点角M表示成时间t的显函数,但由Kepler方程可知,这将涉及超越函数关系,无法直接达到上述要求。因此,必须将等量展开成M的三角级数,而在这些展开式中又要用到两个特殊函数:第一类贝塞耳(Bessel)函数和超几何函数(或称超几何级数),为了读者引用方便,首先简单地介绍一下这两个函数的有关知识,详细内容请翻阅特殊函数一类书籍。
第一类贝塞耳函数Jn(x)是二阶线性常微分方程
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-3.jpg?sign=1739156280-T8nnM12F0J0icOvpnTy46NsySsLokj1V-0-cd2aef944b3a72233e74ca4492a24c99)
的一个解,它由下列级数表达:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-4.jpg?sign=1739156280-Z6H91kyGPivbc8NrJKXZ9ACwIvvjonl6-0-7fb9153dcfeb7ed8023a006d7b1d61ae)
其中n为整数(n=0,1,2,…),x为任意实数,而k!由下式定义
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-5.jpg?sign=1739156280-BLi3dXx1e8NCdMrw7UBCOXkWiqgyWiQC-0-649d9c47d483bbaa1e67e3f711a9dfc8)
Jn(x)又是展开式的函数,即
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-7.jpg?sign=1739156280-YBuKtQd4UG3q2mYrHLY4nzX2jw6vT383-0-c1654dbbe2346ddcd5a4ddc036932336)
其中e是自然对数的底,而z可以是复变量。由此可给出Jn(x)的积分表达式,即
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-8.jpg?sign=1739156280-HC7IL7y2VTzp67XO8Y1HG2tPRvXKHA2V-0-3250e0f6e3b6d1315f30a122178aa6a4)
根据Jn(x)的定义不难得出下列一些重要性质:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-9.jpg?sign=1739156280-4MB7RU18W0cbGvwR6EdQPaQrhCXZXfZp-0-5d0a613a04b6f775b13f210391bd91aa)
超几何函数F(a,b,c;x)是如下二阶线性常微分方程
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-1.jpg?sign=1739156280-KK6h7F0tXMPq2mckEIzMUxPCmvNyBvFs-0-a6b9b05d93fa3447b724137f27d1f954)
的一个解,其中a,b,c是常数,解的形式如下:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-2.jpg?sign=1739156280-AiPoY37XHbzzefmEweegGCZksga7hCL3-0-78c03ca2c7d4feb1799b5cc1e46b2d83)
2.3.1 sinkE和coskE的展开式
这里直接列出展开结果,它们在本章参考文献[1]~[2]中有详细的推导。对k>1有
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-3.jpg?sign=1739156280-tvgnH1LR2SnKLVridVwWbSpZbu0BsVKb-0-0b25b8123bf1bb96fa91c93ad6e4449e)
对k=1有
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-4.jpg?sign=1739156280-xCyZVRJh4lqDdLCUThWWbGPt8HBwhkJ3-0-c0248f4eceb91c4052bde1673562ea43)
2.3.2
和
的展开式
由
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-7.jpg?sign=1739156280-fpazrQFsG0f3Qss65aFPtqLnd4kxcDRN-0-72e77ebd28450d1e1cc80a88014fa912)
可给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-8.jpg?sign=1739156280-UIEz87UyMKf4eiKnck05PKOofPJ4lrha-0-0b0137f02457a5de04a356bd8ce0e9f5)
2.3.3 sinf和cosf的展开式
利用偏导数关系式(2.61)可得
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-1.jpg?sign=1739156280-mqb8FoAGMALt7m6B7BggwgSk7tIMt5F6-0-a494d20718884e718b0260f149111c7b)
由此给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-2.jpg?sign=1739156280-o61Y8CyYEzNjzFl1AYDvLZtalP1iP3rC-0-d6e5b662b0af661cc9c53c6e464f4e57)
由轨道方程(2.16)给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-3.jpg?sign=1739156280-soDHKx5JW1aSWPW8F25IKGcoTwDDqf2W-0-0e45ea4e0aa6af3ef3d30c08767a5348)
2.3.4 f的展开式
利用sinf和cosf的展开式,取到e4项有
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-4.jpg?sign=1739156280-DT7AejT3RUBRryxrQhWG8S3uqCFGjWM7-0-9c38f5cbeb4afd9cfb06d97fc0bb5162)
2.3.5
和
的展开式
这里n和m均为任意整数(包括零)。若仅用上述基本展开式,要给出这两个函数对M的三角级数(特别是一般表达式),那是相当困难的,下面就对这两个函数直接进行傅立叶(Fourier)展开。函数F(f)展成傅立叶级数的基本形式为
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-7.jpg?sign=1739156280-8i2COzxnqc2Ot0PT8dd4SVs2IKh7jZ5E-0-f3e04fe2457f693d806045efc8c7aa06)
是偶函数,bp=0,且
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-9.jpg?sign=1739156280-4VlWZQJbd0xKi6bZOXtdOITGUeNi8ikm-0-c0482471ab782b081b1ba796c0968ff1)
对于被积函数的第二部分,可令p=-p,对应p=-1,-2,…,-∞,由此给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-10.jpg?sign=1739156280-hnhqBdfUnNJIPxV0MxTTzInMpj3BMzQJ-0-896f21c46023de9a34f9cdfab6e38de7)
其中
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-11.jpg?sign=1739156280-d8hNFSJJNw4kb9nH9nCJQVS29361LmRl-0-9ff307adcee49378e0c9aff448c08d0f)
是奇函数,ap=0,bp的计算公式为
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-2.jpg?sign=1739156280-fBkX0tvbKDQ4nvLQLEQ8JjWNzM7MC6Ae-0-c6b00daa8d023967f1018147cf93474b)
对被积函数第二部分的处理同上,结果为
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-3.jpg?sign=1739156280-LFMknx3jqf9uUduJcf7PPxpUb9jNdIgH-0-bbccc3c7778b89714c80c1daaab17c94)
由于上述两个函数的展开式系数相同,可用指数形式将它们表达成统一形式,即
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-4.jpg?sign=1739156280-8wkvJ7U2lFuLAJ63xuZhos8Dk65wN6ik-0-593565e8a25234ef79a0c38a6b002748)
其中是虚数单位。因
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-6.jpg?sign=1739156280-HMzmOOJGJES2J3ix5dQeIEnXK80HOsP8-0-baecb9f14bc9fd12ffbef806203af25e)
(2.99)式中的就是由(2.96)式表达的
,称为汉森(Hansen)系数,它是偏心率e的函数,无法用初等函数来表达它的具体形式,只能引用贝塞耳函数和超几何函数,详细推导见本章参考文献[7],这里直接列出展开结果。
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-9.jpg?sign=1739156280-pwbNexpO39mvp8vZe0H5V2rrGy6TJB3F-0-2049b93c0cee192529f7b19bab8badf9)
由(2.96)式即可给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-10.jpg?sign=1739156280-FMe6ppOguohFOKRQ84ofnU6vHDg0ycj1-0-393d692bc74225f4b6a14e6c33a7f1db)
又根据可知
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-12.jpg?sign=1739156280-BQ2jCdZV6WOHoabpFh923RVDF4FetlSa-0-d38033dc35e75cd05f83327f9774db63)
由上述展开式可以看出,要具体给出和
的展开式,是较麻烦的。为此,针对实际应用状况,作者给出了精确到O(e4)的
表达式[8],形式如下:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-16.jpg?sign=1739156280-NrouCx8hxco4iIWzNY5kloPFqCd2JYqi-0-365044a940ef6f98615ce119fe794fdb)
以上各展开式的系数都是关于偏心率e的无穷级数,只有当e<e1=0.6627…时才收敛,e1就称为拉普拉斯(Laplace)极限。
除上述展开式外,有些问题还需要其他类型的展开式,下面给出。
2.3.6
对f的展开式
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-18.jpg?sign=1739156280-kbFhrD8wiYoxkvBLkxnBjrVPcAg4Gk1f-0-c58beaa82db0e3663482041acab5a6b5)
其中p为正、负整数,β的意义同前,见(2.101)式,Tn(p,q)由超几何函数定义[2],即
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-1.jpg?sign=1739156280-SkiZQQKhfaVnWsLheMWdt9Xr5sQuDijt-0-f80ebc6ba2ee69564b88e3a812973ee7)
当p=-1,-2时,有
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-2.jpg?sign=1739156280-GptVnmGoslXUNEG5ahQElsRtVEeaff2s-0-12a7bcd93809b8de087f9891672b8ae6)
由上述一般表达式可给出如下两个具体展开式:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-3.jpg?sign=1739156280-oXjmlWQM6JuXV5OcDK2hIHmAzLboi39b-0-052c23854a67156fd6d4dea4e65095d8)
利用这两个展开式,由
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-4.jpg?sign=1739156280-aRUvGXgONyle0prJjgLgDILhFQ6O5Fhl-0-bf49041e82fd0e35a6e62b15a8bbc952)
积分即给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-5.jpg?sign=1739156280-iwzWTUoZwcq8h28doAKCoi9vleZ30lXk-0-f3ce1fc55244f2361c68977aacfa8f35)