CN108507569A - 用于星光/惯性复合制导的弹载恒星星库快速生成方法 - Google Patents

Abstract

本发明公开了一种用于星光/惯性复合制导的弹载恒星星库快速生成方法，充分考虑弹载计算机的存储和计算速度条件，按照星云/星团、惯性平台限制、地平规避、强光源规避的次序生成弹载星库，利用天极转换法实现地平规避，占用内存少，计算速度快。本发明星云/星团的规避计算在将恒星数据库装载到导弹上之前完成，减小了弹上计算量；弹上生成星库时，按照惯性平台限制、地平规避、强光源规避的次序生成弹载星库，最大限度地降低了弹载计算机的计算量，降低了对计算机的存储限制；利用天极转换法实现地平规避，降低了已有算法的计算复杂度，且提高了计算精确度。

Description

用于星光/惯性复合制导的弹载恒星星库快速生成方法

技术领域

本发明涉及制导控制技术领域，可应用于远程弹道导弹、运载火箭的主动段星光/惯性复合制导，尤其是应用于复合制导中星光制导的弹载恒星星库的生成。

背景技术

现有的远程弹道导弹、运载火箭的主动段一般采用惯性制导技术，即通过加速度计、陀螺等惯性器件测量载体的视加速度、角速度等信息，经导航积分获得载体的实时运行状态，并传递给制导系统。制导系统计算得到制导指令，控制载体的运动状态达到期望值。在惯性制导过程中，不可避免地存在初始定位定向误差、初始对准误差、惯性器件漂移误差等误差因素的影响，使惯性制导的精度降低，为此提出星光/惯性复合制导方法。

星光/惯性复合制导，是一种星光制导和惯性制导相结合的复合制导方式，更准确地说，是在纯惯性制导的基础上辅以星光制导。它利用在弹上测量恒星所获得的准确的惯性空间方位基准，来估计惯性导航平台(物理平台或数学平台)与发射惯性坐标系之间的失准角，并根据失准角估计值来修正惯性导航误差造成的落点偏差，以达到综合利用惯性导航信息和星光信息来提高导弹制导精度的目的。星光/惯性复合制导需要在飞行过程中通过星敏感器测量导航星的方位，而导航星是在射前从弹载恒星星库中选定的。考虑到星敏感器测量的限制，生成弹载星库时需要考虑惯性平台框架角限制，星云/星团规避，地平规避，太阳、月球、大行星等强光源的规避；同时还要考虑弹载计算机的存储和运算性能限制，要求生成弹载恒星星库占用的存储量尽可能小，运算量尽可能低。

参考文献：

[1]陈克俊，刘鲁华，孟云鹤编著.远程火箭飞行动力学与制导[M].北京：国防工业出版社，2014.

[2]王宏力,陆敬辉，崔祥祥著.大视场星敏感器星光制导技术及应用[M].北京：国防工业出版社，2015.

[3]刘朝山，刘光斌，王新国，李爱华编著.弹载星敏感器原理及系统应用[M].北京：国防工业出版社，2010.

发明内容

本发明所要解决的技术问题是，针对现有技术不足，提供一种用于星光/惯性复合制导的弹载恒星星库快速生成方法，考虑星云/星团遮挡、惯性平台框架角限制、地球遮挡和强光源规避问题，占用内存少，计算速度快。

为解决上述技术问题，本发明所采用的技术方案是：一种用于星光/惯性复合制导的弹载恒星星库快速生成方法，包括以下步骤：

1)获取原始的恒星数据库，所述恒星数据库内包括恒星的编号ID、视星等M_V、赤经α₀、赤纬δ₀、赤经自行dα/dt、赤纬自行dδ/dt，所述恒星数据库的技术文档包括恒星数据库所采用的基准历元MJD₀、基准惯性坐标系ICRS、以及ICRS与J2000.0地心协议惯性坐标系的转换关系矩阵获取原始的星云/星团数据库，其中包括星云/星团亮度，以及星云/星团在J2000.0坐标系中的赤经α_N、赤纬δ_N信息；

2)根据星光/惯性复合制导中对导航星星等M_VN的要求，从原始的恒星数据库中剔除掉星等M_V大于M_VN的恒星，设剩余恒星的数量为N₀；设导弹上所用星敏感器可测量恒星视星等的上限为M_VM，则剔除掉星云/星团数据库中星等高于M_VM的星云/星团，设剩余星云/星团的个数为N_Ne；

3)根据恒星的赤经α₀、赤纬δ₀、赤经自行dα/dt、赤纬自行dδ/dt、基准历元MJD₀，计算历元J2000.0时刻的赤经α和赤纬δ，根据历元J2000.0时刻的赤经α和赤纬δ计算恒星在ICRS坐标系中的单位方向矢量i_ICRS，根据恒星在ICRS坐标系中的单位方向矢量i_ICRS计算恒星在J2000.0坐标系中的单位方向矢量i_I和发射惯性系与J2000.0坐标系的方向余弦阵根据i_I和计算恒星在发射惯性系中的单位方向矢量、高低角e和方位角σ；计算星云/星团在J2000.0坐标系中的单位方向矢量i_Ne，对恒星数据库中的任意一颗恒星P，设该恒星P在J2000.0坐标系中的单位方向矢量为i_I,P；对星云/星团数据库中的任意一个星云/星团Q，设其在J2000.0坐标系中的单位方向矢量为i_Ne,Q；若两个单位方向矢量满足下式，则从经步骤2)处理后得到的恒星数据库中剔除该恒星：i_I,P·i_Ne,Q≥cosα_Ne,0＝常数；其中，α_Ne,0是根据星敏感器性能给定的星云/星团规避角；

4)根据下式将高低角和方位角不满足要求的恒星从步骤3)处理后得到的恒星数据库中剔除(下式确定的为不满足要求的恒星)：

|e|＞e_max,|σ|＞σ_max；

5)对步骤4)得到的恒星数据库中的每一颗恒星，依次判断下式是否成立：

其中，R_E为地球平均半径，h为大气层高度，为测星点到地心的距离；

若上式成立，说明恒星被地平遮挡，则从恒星数据库中剔除此颗恒星；其中，Δθ为距离地球大气上界的附加规避角；

对恒星数据库中的每一颗恒星，依次判断下式是否成立：

i_I·i_sun≥cosα_sun,0＝常数；

如果上式成立，说明恒星与太阳的距离不满足规避角要求，则从恒星数据库中剔除此颗恒星；其中，α_sun、δ_sun分别为太阳的赤经和赤纬；α_sun,0为太阳的规避角；

对恒星数据库中的每一颗恒星，依次判断下式是否成立：

i_I·i_moon≥cosα_moon,0＝常数；

如果上式成立，说明恒星与月球的距离不满足规避角要求，则从恒星数据库中剔除此颗恒星；其中，α_moon、δ_moon分别为月球的赤经和赤纬；α_moon,0为月球的规避角；

对恒星数据库中的每一颗恒星，依次判断下式是否成立：

i_I·i_planet≥cosα_planet,0＝常数；

如果上式成立，说明恒星与大行星的距离不满足规避角要求，则从恒星数据库中剔除此颗恒星；其中，α_planet、δ_planet分别为大行星的赤经和赤纬；

6)将经步骤5)处理后得到的恒星数据库中各恒星的编号、高低角、方位角作为基本信息，生成供弹上存储用的弹载恒星星库。

所述历元J2000.0时刻的赤经α和赤纬δ的计算公式为：

所述恒星在ICRS坐标系中的单位方向矢量i_ICRS的计算公式为：

所述恒星在J2000.0坐标系中的单位方向矢量i_I的计算公式为：

所述发射惯性系与J2000.0坐标系的方向余弦阵的计算公式为：

其中，A₀为发射方位角，λ₀、φ₀为发射点的经度、纬度，Ω_G为发射时刻格林尼治天文台的恒星时。

所述恒星在发射惯性系中的单位方向矢量i_L的计算公式为：

所述恒星在发射惯性系的高低角e和方位角σ的计算公式为：

其中，i_lx、i_ly、i_lz分别表示恒星在发射惯性系中的单位方向矢量i_L在发射惯性系x、y、z三轴上的分量。

本发明中，0＜e_max≤90°，0＜σ_max≤180°。

与现有技术相比，本发明所具有的有益效果为：本发明充分考虑弹载计算机的存储和计算速度条件，按照星云/星团、惯性平台限制、地平规避、强光源规避的次序生成弹载星库，利用天极转换法实现地平规避，占用内存少，计算速度快。本发明优点在于：(1)星云/星团的规避计算在将恒星数据库装载到导弹上之前完成，减小了弹上计算量；(2)弹上生成星库时，按照惯性平台限制、地平规避、强光源规避(太阳、月球、大行星)的次序生成弹载星库，最大限度地降低了弹载计算机的计算量，降低了对计算机的存储限制；(3)利用天极转换法实现地平规避，降低了已有算法的计算复杂度，且提高了计算精确度。

附图说明

图1地平规避的几何关系；

图2为弹载恒星星库生成流程图；

图3为弹载恒星星库的方位角、高低角分布(横轴为方位角，纵轴为高低角)。

具体实施方式

如图2所示，本发明实现步骤如下：

1)获取原始的恒星数据库。该数据库应包括恒星的编号ID、视星等M_V、赤经α₀、赤纬δ₀、赤经自行dα/dt、赤纬自行dδ/dt等信息，数据库的技术文档应包括恒星数据库所采用的基准历元MJD₀、基准惯性坐标系ICRS、以及ICRS与J2000.0地心协议惯性坐标系(以下简称J2000.0坐标系)的转换关系

2)根据星光/惯性复合制导中对导航星星等M_VN的要求，从原始数据库中剔除掉星等M_V大于M_VN的恒星，设剩余恒星的数量为N₀。

3)恒星信息预处理。根据恒星的赤经α₀、赤纬δ₀、赤经自行dα/dt、赤纬自行dδ/dt、基准历元MJD₀，计算历元J2000.0时刻的赤经α和赤纬δ

计算恒星在ICRS坐标系中的单位方向矢量

计算恒星在J2000.0坐标系中的单位方向矢量

计算发射惯性系与J2000.0坐标系的方向余弦阵

式中，A₀为发射方位角，λ₀、φ₀为发射点的经度、纬度，Ω_G为发射时刻格林尼治天文台的恒星时。

计算恒星在发射惯性系中的单位方向矢量

计算恒星在发射惯性系的高低角e和方位角σ

4)星云/星团规避。星云/星团作为面光源，生成星库过程中需要对其进行规避。首先获取原始的星云/星团数据库，其中包括星云/星团亮度，以及在J2000.0坐标系中的赤经α_N、赤纬δ_N信息。

设导弹上所用星敏感器可测量恒星视星等的上限为M_VM，则剔除掉星云/星团数据库中星等高于M_VM的星云/星团，设剩余星云/星团的个数为N_Ne。

计算星云/星团在J2000.0坐标系中的单位方向矢量

对恒星数据库中的任意一颗恒星P，设其在J2000.0坐标系中的单位方向矢量为i_I,_P；对星云/星团数据库中的任意一个星云/星团Q，设其在J2000.0坐标系中的单位方向矢量为i_Ne,Q；若两个单位方向矢量满足下式，则从恒星数据库中剔除该恒星

i_I,P·i_Ne,Q≥cosα_Ne,0＝常数

式中，α_Ne,0是根据星敏感器性能给定的星云/星团规避角。

5)惯性平台框架角限制。惯性平台框架角对恒星的高低角和方位角有限制，要求高低角不高于e_max，方位角不大于σ_max。根据3)中计算的恒星的高低角和方位角，根据下式将高低角和方位角不满足要求的恒星从星库中剔除

|e|＞e_max,|σ|＞σ_max

6)地平规避。地球中心(以下简称地心)用O_E表示，发射点用O表示，测星点用R表示。发射点与地心的距离为测星点与发射点的距离为l_OR，测星点与地心的距离为测星点在发射惯性坐标系中的位置矢量为

计算角度α与β

通过下式所述转换矩阵，将天极坐标系转换为以地心为新的天极的新天极坐标系

新的天极坐标系中，星光的单位方向矢量为

计算角度θ

式中，R_E为地球平均半径，h为大气厚度，Δθ为距离地球大气上界的附加规避角。

对恒星星库中的每一颗恒星，依次判断下式是否成立

z_p≥cosθ＝cos(θ₀+Δθ)＝常数

如果上式成立，说明恒星被地平遮挡，则从星库中剔除此颗恒星。

7)强光源规避。星光/惯性复合制导中，需要考虑的强光源规避包括太阳、月球和大行星。设太阳、月球、大行星的规避角分别为α_sun,0、α_moon,0和α_planet,0。从天文年历中获得太阳的赤经α_sun和赤纬δ_sun，计算太阳在J2000.0坐标系中的单位方向矢量

对恒星星库中的每一颗恒星，依次判断下式是否成立

i_I·i_sun≥cosα_sun,0＝常数

如果上式成立，说明恒星与太阳的距离不满足规避角要求，则从星库中剔除此颗恒星。

从天文年历中获得月球的赤经α_moon和赤纬δ_moon，计算月球在J2000.0坐标系中的单位方向矢量

对恒星星库中的每一颗恒星，依次判断下式是否成立

i_I·i_moon≥cosα_moon,0＝常数

如果上式成立，说明恒星与月球的距离不满足规避角要求，则从星库中剔除此颗恒星。

对水星、金星、火星、木星、土星、天王星六大行星，依次从天文年历中获得其赤经α_planet和赤纬δ_planet，计算大行星在J2000.0坐标系中的单位方向矢量

对恒星星库中的每一颗恒星，依次判断下式是否成立

i_I·i_planet≥cosα_planet,0＝常数

如果上式成立，说明恒星与大行星的距离不满足规避角要求，则从星库中剔除此颗恒星。

8)存储弹载恒星星库。对完成星云/星团、地平、强光源规避的恒星星库，将各恒星的编号、高低角、方位角作为基本信息，生成供弹上存储用的弹载恒星星库。以生成某远程弹道导弹星光/惯性复合制导系统的弹载恒星数据库为例，说明本发明的实施方式。

设发射时间为2018年1月1日12时15分30.0秒，发射点的纬度为20度，经度为100度，发射点高程为200m，发射方位角为45度。

S1.以依巴谷(Hipparchus)星表作为原始的恒星数据库，该星表可由因特网上的中国天文数据中心(CAsDC)获得。星表中包括118218颗恒星，包含了恒星的依巴谷编号(HIPID)、赤经、赤纬、赤经自行、赤纬自行等信息。说明文档中指明，恒星数据库采用的基准历元为J1991.25，即MJD₀＝48348.5625。基准惯性坐标系与J2000.0协议地心惯性坐标系的方向余弦阵为

式中

M_x、M_y、M_z为基础方向余弦阵(参考文献[1])。

S2.设星光/惯性复合制导中要求导航星星等M_VN≤5.5等。从原始恒星数据库中剔除掉星等M_V大于5.5等的恒星，剩余恒星数量2226颗。

S3.恒星信息预处理。根据依巴谷星表中提供的恒星赤经α₀、赤纬δ₀、赤经自行dα/dt、赤纬自行dδ/dt、基准历元MJD₀＝48348.5625等信息，根据公式(2)计算2226颗恒星在历元J2000.0时刻的赤经α和赤纬δ

根据公式(3)依次计算2226颗恒星在ICRS坐标系中的单位方向矢量

根据公式(4)依次计算2226颗恒星在J2000.0坐标系中的单位方向矢量

根据给定的发射时刻，由中国天文年历可查得发射时格林尼治天文台的恒星时为Ω_G＝285.2598°。由此可根据公式(5)计算得到发射惯性系与J2000.0坐标系的方向余弦阵

根据公式(6)依次计算2226颗恒星在发射惯性系中的单位方向矢量

根据公式(7)依次计算2226颗恒星在发射惯性系的高低角e和方位角σ

S4.星云/星团规避。星云/星团作为面光源，生成星库过程中需要对其进行规避。原始的星云/星团数据库采用修订的星云星团总表(RNGC)，可以从因特网上的中国天文数据中心(CAsDC)获得。RNGC包括星云/星团在J2000.0坐标系中的赤经α_N、赤纬δ_N信息，以及星云/星团的亮度信息。

设导弹上所用星敏感器可测量恒星视星等的上限为M_VM＝6.5等，则剔除掉星云/星团数据库中星等高于M_VM＝6.5等的星云/星团，剩余星云/星团的个数为N_Ne＝114个。

根据公式(8)依次计算114个星云/星团在J2000.0坐标系中的单位方向矢量

设根据星敏感器性能，给定星云/星团的规避角α_Ne,0＝3°，则

cosα_Ne,0＝0.998629

对恒星数据库中2226颗恒星中的任意一颗P，设其在J2000.0坐标系中的单位方向矢量为i_I,P；对星云/星团数据库中114个星云/星团中的任意一个Q，设其在J2000.0坐标系中的单位方向矢量为i_Ne,Q；若两者满足公式(9)，则从恒星数据库中剔除该恒星

i_I,P·i_Ne,Q≥cosα_Ne,0＝0.998629 (9)

经计算，共有203颗恒星满足公式(9)。故考虑星云/星团规避后，剩余2023颗恒星。星云/星团的规避在装载到导弹上之前完成，故装载到导弹上的初始恒星数据库共有恒星2023颗。

S5.惯性平台框架角限制。惯性平台框架角对恒星的高低角和方位角有限制。设要求高低角不高于e_max＝75°，方位角不大于σ_max＝75°。根据S3中计算的恒星的高低角和方位角，根据下式将高低角和方位角不满足要求的恒星从星库中剔除

|e|＞e_max＝75°,|σ|＞σ_max＝75° (10)

经计算，高低角限制共剔除恒星55颗，方位角共剔除恒星1158颗，剩余恒星810。

S6.地平规避。地球中心(以下简称地心)用O_E表示，发射点用O表示，测星点用R表示。

设测星点在发射惯性坐标系中的位置矢量为

根据给定的发射点位置，可计算得到发射点与地心的距离为测星点与发射点的距离为l_OR＝824.621km，测星点与地心的距离为计算公式计算角度α与β

将天极坐标系转换为以地心为新的天极的新天极坐标系，转换的方向余弦阵由公式(13)计算

根据公式(14)，依次计算810颗恒星在新的天极坐标系中的单位方向矢量为

地球平均半径为6371km，取地球大气厚度h＝90km，距离地球大气上界的附加规避角Δθ＝10°，则根据公式(15)计算角度θ

则

cosθ＝cos(θ₀+Δθ)＝0.047973

对恒星星库中剩余810颗恒星的任一颗，依次判断下式是否成立

z_p≥cosθ＝cos(θ₀+Δθ)＝常数 (15)

如果上式成立，说明恒星被地平遮挡，则从星库中剔除此颗恒星。

经计算，地平共遮挡恒星287颗，故剩余恒星523颗。

S7.强光源规避。星光/惯性复合制导中，需要考虑的强光源规避包括太阳、月球和大行星。设太阳、月球、大行星的规避角分别为α_sun,0＝20°、α_moon,0＝10°和α_planet,0＝3°。故有

cosα_sun,0＝0.939692,cosα_moon,0＝0.984808,cosα_planet,0＝0.998629S7.1从中国天文年历中，可以获得发射时刻太阳的赤经α_sun＝-77.999°和赤纬δ_sun＝-22.978°。根据公式(16)计算太阳在J2000.0坐标系中的单位方向矢量

对恒星星库中剩余的523颗恒星中的任一颗，依次判断下式是否成立

i_I·i_sun≥cosα_sun,0＝0.939692

如果上式成立，说明恒星与太阳的距离不满足规避角要求，则从星库中剔除此颗恒星。经计算，共剔除恒星0颗，故太阳没有遮挡星库中的恒星，剩余恒星523颗。

S7.2从天文年历中获得月球的赤经α_moon＝92.757°和赤纬δ_moon＝19.882°，根据公式(17)计算月球在J2000.0坐标系中的单位方向矢量

对恒星星库中剩余的523颗恒星中的任一颗，依次判断下式是否成立

i_I·i_moon≥cosα_moon,0＝0.984808

如果上式成立，说明恒星与月球的距离不满足规避角要求，则从星库中剔除此颗恒星。经计算，共剔除恒星12颗，故太阳没有遮挡星库中的恒星，剩余恒星511颗。

S7.3对水星、金星、火星、木星、土星、天王星六大行星，依次从天文年历中获得其赤经α_planet和赤纬δ_planet。以金星为例，其赤经α_venus＝-80.201°，赤纬δ_venus＝22.642°，根据公式(18)计算金星在J2000.0坐标系中的单位方向矢量

对恒星星库中剩余的511颗恒星中的任一颗，依次判断下式是否成立

i_I·i_venus≥cosα_planet,0＝0.998629

如果上式成立，说明恒星与金星的距离不满足规避角要求，则从星库中剔除此颗恒星。

经计算，六颗大行星都不遮挡任何恒星，故最终剩余恒星511颗。图3给出了这511颗恒星的高低角、方位角的分布情况。

S8.存储弹载恒星星库。对完成星云/星团、地平、强光源规避的恒星星库，将剩余的511颗恒星的编号(HIP ID)、高低角、方位角作为基本信息，生成供弹上存储用的弹载恒星星库。

Claims (8)

1.一种用于星光/惯性复合制导的弹载恒星星库快速生成方法，其特征在于，包括以下步骤：

1)获取原始的恒星数据库，所述恒星数据库内包括恒星的编号ID、视星等M_V、赤经α₀、赤纬δ₀、赤经自行dα/dt、赤纬自行dδ/dt，所述恒星数据库的技术文档包括恒星数据库所采用的基准历元MJD₀、基准惯性坐标系ICRS、以及ICRS与J2000.0地心协议惯性坐标系的转换关系矩阵获取原始的星云/星团数据库，其中包括星云/星团亮度，以及星云/星团在J2000.0坐标系中的赤经α_N、赤纬δ_N信息；

2)根据星光/惯性复合制导中对导航星星等M_VN的要求，从原始的恒星数据库中剔除掉星等M_V大于M_VN的恒星，设剩余恒星的数量为N₀；设导弹上所用星敏感器可测量恒星视星等的上限为M_VM，则剔除掉星云/星团数据库中星等高于M_VM的星云/星团，设剩余星云/星团的个数为N_Ne；

3)根据恒星的赤经α₀、赤纬δ₀、赤经自行dα/dt、赤纬自行dδ/dt、基准历元MJD₀，计算历元J2000.0时刻的赤经α和赤纬δ，根据历元J2000.0时刻的赤经α和赤纬δ计算恒星在ICRS坐标系中的单位方向矢量i_ICRS，根据恒星在ICRS坐标系中的单位方向矢量i_ICRS计算恒星在J2000.0坐标系中的单位方向矢量i_I和发射惯性系与J2000.0坐标系的方向余弦阵根据i_I和计算恒星在发射惯性系中的单位方向矢量i_L、高低角e和方位角σ；计算星云/星团在J2000.0坐标系中的单位方向矢量i_Ne，对恒星数据库中的任意一颗恒星P，设该恒星P在J2000.0坐标系中的单位方向矢量为i_I,P；对星云/星团数据库中的任意一个星云/星团Q，设其在J2000.0坐标系中的单位方向矢量为i_Ne,Q；若两个单位方向矢量满足下式，则从经步骤2)处理后得到的恒星数据库中剔除该恒星：i_I,P·i_Ne,Q≥cosα_Ne,0＝常数；其中，α_Ne,0是根据星敏感器性能给定的星云/星团规避角；

4)根据下式将高低角和方位角不满足要求的恒星从步骤3)处理后得到的恒星数据库中剔除：

|e|＞e_max,|σ|＞σ_max；

5)对步骤4)得到的恒星数据库中的每一颗恒星，依次判断下式是否成立：

z_p≥cosθ＝cos(θ₀+Δθ)＝常数；

其中，R_E为地球平均半径，h为大气层高度，为测星点到地心的距离；θ＝θ₀+Δθ，Δθ为距离地球大气上界的附加规避角；z_p是星光在新天极坐标系中单位方向矢量的z向分量，根据天极坐标系与新天极坐标系间的方向余弦阵和恒星在发射惯性系中的单位方向矢量i_L计算；方向余弦阵是发射点和测星点位置的函数；若上式成立，说明恒星被地平遮挡，则从恒星数据库中剔除此颗恒星；对恒星数据库中的每一颗恒星，依次判断下式是否成立：

i_I·i_sun≥cosα_sun,0＝常数；

如果上式成立，说明恒星与太阳的距离不满足规避角要求，则从恒星数据库中剔除此颗恒星；其中，α_sun、δ_sun分别为太阳的赤经和赤纬；α_sun,0为太阳的规避角；

对恒星数据库中的每一颗恒星，依次判断下式是否成立：

i_I·i_moon≥cosα_moon,0＝常数；

如果上式成立，说明恒星与月球的距离不满足规避角要求，则从恒星数据库中剔除此颗恒星；其中，α_moon、δ_moon分别为月球的赤经和赤纬；α_moon,0为月球的规避角；对恒星数据库中的每一颗恒星，依次判断下式是否成立：

i_I·i_planet≥cosα_planet,0＝常数；

α_planet,0为大行星规避角；如果上式成立，说明恒星与大行星的距离不满足规避角要求，则从恒星数据库中剔除此颗恒星；其中，α_planet、δ_planet分别为大行星的赤经和赤纬；

6)将经步骤5)处理后得到的恒星数据库中各恒星的编号、高低角、方位角作为基本信息，生成供弹上存储用的弹载恒星星库。

2.根据权利要求1所述的用于星光/惯性复合制导的弹载恒星星库快速生成方法，其特征在于，所述历元J2000.0时刻的赤经α和赤纬δ的计算公式为：

3.根据权利要求2所述的用于星光/惯性复合制导的弹载恒星星库快速生成方法，其特征在于，所述恒星在ICRS坐标系中的单位方向矢量i_ICRS的计算公式为：

4.根据权利要求3所述的用于星光/惯性复合制导的弹载恒星星库快速生成方法，其特征在于，所述恒星在J2000.0坐标系中的单位方向矢量i_I的计算公式为：

5.根据权利要求4所述的用于星光/惯性复合制导的弹载恒星星库快速生成方法，其特征在于，所述发射惯性系与J2000.0坐标系的方向余弦阵的计算公式为：

其中，A₀为发射方位角，λ₀、φ₀为发射点的经度、纬度，Ω_G为发射时刻格林尼治天文台的恒星时。

6.根据权利要求5所述的用于星光/惯性复合制导的弹载恒星星库快速生成方法，其特征在于，所述恒星在发射惯性系中的单位方向矢量i_L的计算公式为：

7.根据权利要求6所述的用于星光/惯性复合制导的弹载恒星星库快速生成方法，其特征在于，所述恒星在发射惯性系的高低角e和方位角σ的计算公式为：

其中，i_lx、i_ly、i_lz分别表示恒星在发射惯性系中的单位方向矢量i_L在发射惯性系x、y、z三轴上的分量。

8.根据权利要求1所述的用于星光/惯性复合制导的弹载恒星星库快速生成方法，其特征在于，0＜e_max≤90°，0＜σ_max≤180°。