參數(shù)說(shuō)明：
$M:地球總質(zhì)量$
$dm:地球上的質(zhì)量微元$
$m^{\prime }:地球外部空間上的一個(gè)質(zhì)點(diǎn)$
$p:m與dm之間的距離$
萬(wàn)有引力大小為：
$f=\frac{Gdm?m^{\prime }}{p^2}$

引力位函數(shù)

$若外有引力做工(從無(wú)窮遠(yuǎn)移動(dòng)到半徑p處)為:A={\int }_{\infty }^p?\frac{Gdm?m}{p^2}=\frac{Gdm?m^{\prime }}{p}$
根據(jù)能量守恒，引力做工必然等于$m^{\prime }$減少的位能（勢(shì)能），若單位質(zhì)量的點(diǎn)(從無(wú)窮遠(yuǎn)移動(dòng)到半徑p處)，勢(shì)能必然減少，將減少量定義為位函數(shù):
$$dV=\frac{Gdm}{p}$$

地球外部空間的調(diào)和函數(shù)

地球產(chǎn)生的各質(zhì)量微元函數(shù)之和為：
$$V={\int }_{M}dV=G{\int }_M\frac{dm}{p}（公式1）$$
又因?yàn)?#xff1a;
$$▽V=\frac{{?}^{2}V}{?x^2}+\frac{{?}^{2}V}{?y^2}+\frac{{?}^{2}V}{?z^2}=0(證略)$$
所以$V$是在地球質(zhì)量M的外部空間上是調(diào)和函數(shù)

勒讓德多項(xiàng)式的生成函數(shù)

參數(shù)說(shuō)明：
$m^{\prime }：的球坐標(biāo)為(r,\theta ,\lambda )$
$dm：的球坐標(biāo)為(r,{\theta }^{\prime },{\lambda }^{\prime })$
$向量R與r夾角為\phi$
$m{0}^{\prime }是m^{\prime }在球面上的投影$
在,$o,dm,m^{\prime }$構(gòu)成的三角形中，根據(jù)余弦定理,可以得到：
$$p^2=r^2+R^2?2Rrcos\phi =r^2(1?2ax+a^2)(a=\frac{R}{r},x=cos\phi )$$
再開(kāi)平方，取倒數(shù)得到：
$$\frac{1}{p}=\frac{1}{r}(1?2ax+a^2{)}^{?\frac{1}{2}}（公式2）$$
又因?yàn)?#xff1a;
$$(1?2ax+a^2{)}^{?\frac{1}{2}}=\sum _{n=0}^{\infty }{P}_n(x)a^n(勒讓德多項(xiàng)式的生成函數(shù))（公式3）$$
所以$(1?2ax+a^2{)}^{?\frac{1}{2}}$為勒讓德多項(xiàng)式的生成函數(shù)

球函數(shù)的加法公式

根據(jù)球面的三角余弦定理有:$cos\phi =cos\theta cos{\theta }^{\prime }+sin\theta sin{\theta }^{\prime }cos(\lambda ?{\lambda }^{\prime })$
帶入$$P_n(cos\phi )=\sum _{k=0}^n\frac{2(n?k)!}{(1+{\delta }_k)(n+k)!}{P}_n^k(cos\theta )P_n^k(cos{\theta }^{\prime })cosk{\lambda }^{\prime }+P_n^{k}cos\theta sin(k\lambda )P_n^k(cos{\theta }^{\prime })(sink{\lambda }^{\prime })(公式4)$$

利用德讓勒多項(xiàng)式求解V（調(diào)和函數(shù)）

將勒讓德多項(xiàng)式的生成函數(shù)（公式(3)），將帶入公式（2），再將公式2帶入地球外部空間的調(diào)和函數(shù)公式（1）。
得到:
$R_e:旋轉(zhuǎn)橢圓球的長(zhǎng)半徑$
$u=GM:地球引力常數(shù)$
${\delta }_k=\{\begin{array}{ll}1& k=0\\ 0& k=1\end{array}$
$$V=\frac{G}{r}\sum _{n=0}^{\infty }(\frac{R_e}{r}{)}^n{\int }_M(\frac{R}{R_e}{)}^n{P}_n(cos\phi )dm$$
將公式（4）帶入:
得到:
$$V=\frac{u}{r}\sum _{n=0}^{\infty }(\frac{R_e}{r}{)}^n\sum _{k=0}^n(C_n^{k}cosk\lambda +S_n^{k}sink\lambda )P_n^k(cos\theta )$$
$$C_n^k=\frac{2(n?k)!}{M(1+{\delta }_k)(n+k)!}{\int }_M(\frac{R}{R_e}{)}^n{P}_n^k(cos{\theta }^{\prime })cosk{\lambda }^{\prime }dm$$
$$S_n^k=\frac{2(n?k)!}{M(1+{\delta }_k)(n+k)!}{\int }_M(\frac{R}{R_e}{)}^n{P}_n^k(cos{\theta }^{\prime })sink{\lambda }^{\prime }dm$$
已知直角坐標(biāo)系到球體坐標(biāo)系的轉(zhuǎn)換關(guān)系如下：
$$?\begin{array}{l}x=Rsin{\theta }^{\prime }cos{\lambda }^{\prime }\\ y=Rsin{\theta }^{\prime }sin{\lambda }^{\prime }\\ z=Rcos{\theta }^{\prime }\end{array}$$
剛體的張量定義如下:剛體的慣性張量及其物理意義
$$I=?\begin{array}{ccc}{I}_x& ?I_{xy}& ?I_{xz}\\ ?I_{xy}& I_y& ?I_{yz}\\ ?I_{xz}& ?I_{yz}& I_{xz}\end{array}?={\int }_M?\begin{array}{ccc}{y}^2+z^2& ?xy& ?xz\\ ?xy& x^2+z^2& ?yz\\ ?xz& ?yz& x^2+y^2\end{array}?dm$$
可以得到：
$C_0^0=1$
$C_0^1=\frac{1}{R_e}(\frac{1}{M}{\int }_{M}zdm)$
$C_1^1=\frac{1}{R_e}(\frac{1}{M}{\int }_{M}xdm)$
$S_1^1=\frac{1}{R_e}(\frac{1}{M}{\int }_{M}ydm)$
$C_2^0=\frac{1}{M{R}_e^2}({\int }_M{z}^2?\frac{x^2+y^2}{2}dm)=\frac{1}{M{R}_e^2}(\frac{I_x+I_y}{2}?I_z)$
$C_2^1=\frac{1}{M{R}_e^2}{\int }_{M}xzdm=\frac{I_{xz}}{M{R}_e^2}$
$C_2^2=\frac{1}{4M{R}_e^2}({\int }_M{x}^2?y^{2}dm)=\frac{I_y?I_x}{4M{R}_e^2}$
$S_2^1=\frac{1}{M{R}_e^2}({\int }_{M}yzdm)=\frac{I_{yz}}{M{R}_e^2}$
$S_2^2=\frac{1}{4M{R}_e^2}({\int }_{M}xydm)=\frac{I_{xy}}{2M{R}_e^2}$
如果定義直角坐標(biāo)與地球慣性慣性主軸重合，則：$I_{xy}=I_{yz}=I_{xz}=0$

對(duì)于實(shí)際應(yīng)用時(shí)，地球坐標(biāo)一邊選擇坐標(biāo)原點(diǎn)與地球質(zhì)心重合，oz與地球自轉(zhuǎn)平行，ox軸在赤道面上且指向0度經(jīng)線，這時(shí)坐標(biāo)軸往往與慣性主軸不重合。
利用三角函數(shù)恒等式：
得到:
$$V=\frac{u}{r}\{1?\sum _{n=1}^{\infty }(\frac{R_e}{r}{)}^n[J_n{P}_n(cos\theta )+\sum _{k=1}^n{J}_n^k{P}_n^k(cos\theta )cosk(\lambda +{\lambda }_n^k)]\}$$
$$J_n=?C_n^0,J_n^k=\sqrt{(C_n^k{)}^2+(S_n^k{)}^2},{\lambda }_n^k=?arctan(\frac{S_n^k}{C_n^k})/k$$
$J_2=C_2^0=\frac{1}{M{R}_e^2}(\frac{I_x+I_y}{2}?I_z)(動(dòng)力扁率：反映赤道與極軸上轉(zhuǎn)動(dòng)慣量的差別)$
$\frac{u}{r}:球形地球引起的引力位$
在實(shí)際應(yīng)用中因?yàn)槠渌禂?shù)比$J_2$小三個(gè)數(shù)量級(jí)，所以只需考慮$\frac{u}{r}和J_2$的影響，可以將$V$近似看為：
$$V=\frac{u}{r}[1?\frac{J_2{R}_e^2}{2r^2}(3co{s}^2\theta ?1)]（只考慮主諧項(xiàng)）$$

衛(wèi)星在慣性坐標(biāo)系下的運(yùn)動(dòng)方程（只考慮主諧項(xiàng)）

已知
$位移矢量：r={\left[\begin{array}{ccc}x& y& z\end{array}\right]}^T$
$引力：f={\left[\begin{array}{ccc}{f}_x& f_y& f_z\end{array}\right]}^T$
$cos\theta =\frac{z}{r}$
對(duì)x,y,z求2階片導(dǎo)，求得加速度方程為：
$$\ddot{r}=?\frac{u}{r^3}[I+\frac{3}{2}{J}_2(\frac{R_e}{r}{)}^2(D?5(u_r?u_p)I)]r$$
其中:$D=diag\left(\begin{array}{ccc}1& 1& 3\end{array}\right)$
$u_r表示r上的單位矢量$
$u_p表示自轉(zhuǎn)軸上的的單位矢量$
$u_r?u_p表示u_p與u_r之間夾角的余弦值$

總結(jié)

以上是生活随笔為你收集整理的捷联惯导系统学习3.3(引力位函数)的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。