![卫星轨道力学算法](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=1739092078-orLVkZdiMArAX6M5ygG05TDGXlI8deni-0-9a5d140f187e6fff67734b2af08bbf75)
的一个解,它由下列级数表达:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-4.jpg?sign=1739092078-ViVfWVHRUzYi79AvTYFsRBBtF9342oR5-0-be75453e0f82db310c2b4e0d8bc34868)
其中n为整数(n=0,1,2,…),x为任意实数,而k!由下式定义
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-5.jpg?sign=1739092078-rkfWp1cdDpiZm3Yi7twoC46InE1CAkVV-0-944e627a82d575c280b027a198c07f3b)
Jn(x)又是展开式的函数,即
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-7.jpg?sign=1739092078-BVs5KkSvCizd0ZdjIeTxM09IaulyGbTt-0-06189372d61cef2f6d56758908f32ead)
其中e是自然对数的底,而z可以是复变量。由此可给出Jn(x)的积分表达式,即
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-8.jpg?sign=1739092078-IpHqgNzq3sFjI6kmI5pPxIm0yzwyLZ8H-0-e183605109cf5af0a5c8e122c126d56e)
根据Jn(x)的定义不难得出下列一些重要性质:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-61-9.jpg?sign=1739092078-Ac2gTxkeLcHeEFnYpaBwOkOP68Sv1nL2-0-14ed545df24ef84dd9a7ac19f8b72249)
超几何函数F(a,b,c;x)是如下二阶线性常微分方程
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-1.jpg?sign=1739092078-G6kVJjRzsm8orGVLaUdsciOqeuTDU44S-0-25c8e9952bcf77966b84df7993b6274b)
的一个解,其中a,b,c是常数,解的形式如下:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-2.jpg?sign=1739092078-WGQ5vJVR5phI0Ifiy6pnBubOO33AeEEI-0-e51b9485cc83ca6c562cfe8ed8073f01)
2.3.1 sinkE和coskE的展开式
这里直接列出展开结果,它们在本章参考文献[1]~[2]中有详细的推导。对k>1有
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-3.jpg?sign=1739092078-4qpp4Flu3z5xaQgkPdKM02MssDndUjCR-0-4dc3cfb95d554329e65a0369d8cbcf95)
对k=1有
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-4.jpg?sign=1739092078-upMdzwNnia64EhLxBVfmbbuFxX3yadbC-0-2fbf17ab09611dc215a51464db184fa7)
2.3.2
和
的展开式
由
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-7.jpg?sign=1739092078-5jONkuAtD1CEWbRHZ7QEOu1NAtonDyvx-0-d1c02a24f6e004b53582a17fc3c5864d)
可给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-62-8.jpg?sign=1739092078-kL03IGljHdWdkbqhGzQJmETNULoPmGWa-0-001264c4004ecba90d8dac23c4e8e6cf)
2.3.3 sinf和cosf的展开式
利用偏导数关系式(2.61)可得
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-1.jpg?sign=1739092078-umThscYSpjAkv2ItxtFF5h6ZS4ILUBZg-0-b2d78c28665cbf3f990b3df47826478e)
由此给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-2.jpg?sign=1739092078-6wK7xT5Ws2WjAhSNJPIPkJi5pAX370sh-0-6b787361c0b8f53e7fc4ad6cb32ee460)
由轨道方程(2.16)给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-3.jpg?sign=1739092078-ZyHZBRK7jPe2RsZ69IX8d3NHd47OF06q-0-ee376225cdbf97ff70dfb541e09093ce)
2.3.4 f的展开式
利用sinf和cosf的展开式,取到e4项有
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-4.jpg?sign=1739092078-urWPYeqAwvfyBwq1dAb6fYe1cSHqEQ7m-0-540afc197ab67f555f97cd8e5404b160)
2.3.5
和
的展开式
这里n和m均为任意整数(包括零)。若仅用上述基本展开式,要给出这两个函数对M的三角级数(特别是一般表达式),那是相当困难的,下面就对这两个函数直接进行傅立叶(Fourier)展开。函数F(f)展成傅立叶级数的基本形式为
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-7.jpg?sign=1739092078-S2lQFDtr6VFVFHTpoQJUtf1H2aB0n96x-0-86dd7ed3807b0bd6e89e6b6795d64293)
是偶函数,bp=0,且
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-9.jpg?sign=1739092078-W1c9V1aZ7UgNMYbGeGEq3vMUR1zPGKZW-0-a658fdd71efd36030df51241b6bcdb5d)
对于被积函数的第二部分,可令p=-p,对应p=-1,-2,…,-∞,由此给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-10.jpg?sign=1739092078-RhSxLQIXHSG1nu0mHPMDazdjemrb3pbX-0-c80fb171031b073ed96a0bfa50660d36)
其中
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-63-11.jpg?sign=1739092078-eHzdUAQBZgVYaFCbxuJWucbK1UzjQR7y-0-54ffacba2c1b4dad9ec8ebacb5c8cd8b)
是奇函数,ap=0,bp的计算公式为
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-2.jpg?sign=1739092078-D2ZBC64ouVNeByQFvYxniIURsy1XagfI-0-b665a4b62a75a7f4cd658b7b6ec18f6f)
对被积函数第二部分的处理同上,结果为
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-3.jpg?sign=1739092078-Inm2KtnG0TJ4H22L4KNgE0DCFmFomqSx-0-9ad386edb1cd50eb3fb2eb24560af9b9)
由于上述两个函数的展开式系数相同,可用指数形式将它们表达成统一形式,即
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-4.jpg?sign=1739092078-fXwugITPHMrvH6LJWgStkcDa9RpajcIZ-0-900e247c88bb75c6a7356d68b9c0d118)
其中是虚数单位。因
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-6.jpg?sign=1739092078-6BrG4bSgSQRB9M4n4ZzHHLuUpPBLwL5n-0-313964ffc9d508bec6d5653f4920860a)
(2.99)式中的就是由(2.96)式表达的
,称为汉森(Hansen)系数,它是偏心率e的函数,无法用初等函数来表达它的具体形式,只能引用贝塞耳函数和超几何函数,详细推导见本章参考文献[7],这里直接列出展开结果。
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-9.jpg?sign=1739092078-Sur49apT9SWJB6UX91Ra3VkZF38vQGhe-0-5c0b880ee6d55a6a8b97ec8e5bfcc81d)
由(2.96)式即可给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-10.jpg?sign=1739092078-9pngGvaRPyc9vD8SF207ZWrCHoYJjcsy-0-a03052bc8ef807778e6435c7b179e6ab)
又根据可知
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-12.jpg?sign=1739092078-LfJ87HEtGZm2v3dCUUudrnb8ypjNhkSL-0-6798e98f5969da24f3630ea668fe65d0)
由上述展开式可以看出,要具体给出和
的展开式,是较麻烦的。为此,针对实际应用状况,作者给出了精确到O(e4)的
表达式[8],形式如下:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-64-16.jpg?sign=1739092078-HZyLFpLcWH6iEPqPVVjDwgIxpGJBIs7P-0-bda230441ff339f4e739ed36ef292af4)
以上各展开式的系数都是关于偏心率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=1739092078-zHh4IFdLob57zJQgbqZSqyqMfQlXP1vQ-0-b218317934540be285b80266366bb711)
其中p为正、负整数,β的意义同前,见(2.101)式,Tn(p,q)由超几何函数定义[2],即
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-1.jpg?sign=1739092078-YORIUEOEugbQzOStpfVnBgKOrftlGehr-0-c30bc5e9aa2a4c671ebf3066de586435)
当p=-1,-2时,有
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-2.jpg?sign=1739092078-aVXfIaybfhLACv8TRZX8pCqn9dWXI0DX-0-332b9769b2099a6f489e144116821917)
由上述一般表达式可给出如下两个具体展开式:
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-3.jpg?sign=1739092078-DmNYcDIv5j4PdhkjOLMS1SWgZayY2GqD-0-41a92d0b019d7332b712c2937d3d99d3)
利用这两个展开式,由
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-4.jpg?sign=1739092078-OjTMHCxcXTtFZxEX87LNqpnBFxFYSWeq-0-c5feda67f5ed8d8af09c40ac8f165484)
积分即给出
![](https://epubservercos.yuewen.com/A499B2/17770365207625606/epubprivate/OEBPS/Images/image-66-5.jpg?sign=1739092078-xNPMqqKK4tTQ8sCEEObqTV4A7aRbQdNZ-0-b3f4ea5b4c3f5db8b80782c2a77dade6)