一种折射星的模拟方法技术

本发明专利技术涉及一种折射星的模拟方法，其特征在于：步骤1：扫描星表，判断当前恒星是否为折射星；若判断为折射星,则进入步骤2,否则继续扫描；步骤2：计算星光折射角；步骤3：求解折射星视位置的赤经、赤纬；步骤4：判断折射星的折射星光是否被星敏感器捕捉到；若判断折射星的折射星光被星敏感器捕捉到,则进入步骤5,否则返回步骤1；步骤5：进行星点成像中心计算并生成模拟折射星图。

【技术实现步骤摘要】
一种折射星的模拟方法
本专利技术涉及一种折射星的模拟方法。
技术介绍
随着载人航天、深空探测等空间技术的发展，各类航天器对自主导航能力的要求越来越高。基于星敏感器的天文导航方法由于导航精度高、误差不随时间积累等特点成为了空间应用中的有效自主导航手段。在影响天文导航精度的因素中，地平的敏感精度是很重要的一点，根据星敏感器敏感地平方式的不同，可分为直接敏感地平和利用星光折射间接敏感地平两种方法。直接敏感地平法受地平仪精度的制约，导航精度较低。利用星光折射间接敏感地平的方法是一种低成本、高精度的自主导航方法。该方法的研究工作可以追溯到20世纪60年代，美国在实施Apollo计划的过程中，就对利用天体掩星、星光在大气中的折射、星光穿越大气时的衰减等实现自主导航的方案进行了研究.1975年由美国海军研究局和美国国防部高级研究计划局共同投资，麻省理工学院Draper实验室对星光折射/星光色散自主导航方案进行了研究和论证，结果表明在一个轨道周期可观测40颗折射星的理想条件下，导航精度可以达到100m。90年代初投入使用的MADAN(multi-missionattitudedeterminationandautonomousnavigation)导航系统(多任务姿态确定和自主导航系统)便利用了星光折射原理。二十世纪80年代初期，法国也进行了星光折射法自主导航的研究。1985年和1986年，CNES多次释放平流层气球对星光折射进行了实际测量，在此基础上，对大气折射的精确模型、测量方案、自然环境对系统观测的约束、误差分配和系统性能优化等方面进行了深入的分析和仿真试验，当时预计该系统导航精度为300m。星光折射间接敏感地平法精度高且成本低，是一种很有发展前途的导航方法，因此吸引了大量的学者对其进行研究。然而，由于航天实验费用巨大，往往不可能进行实时的星空拍摄，因此目前的折射星导航的研究，大多数都是基础理论的研究，实验验证比较困难。还有极少一部分研究成果是结合星光折射的几何原理和大气折射模型直接计算星光折射角，然后使用星光折射角来研究星光折射导航。这种方法直接使用了计算得到的星光折射角，没有对折射星图进行模拟，忽略了环境噪声、星敏感器噪声、星点位置噪声等众多噪声的影响，与真实情况差距较大，因此这些研究成果对于实际应用还有一定的距离。鉴于此，本专利技术提供了一种折射星的模拟方法，而且能够根据实际情况考虑星敏感器噪声、环境噪声的影响，为折射星导航的研究提供了一种更加符合实际情况且操作简单，成本低廉的方法。
技术实现思路
本专利技术目的在于提供一种不但能够模拟折射星在星敏感器像平面的位置，而且能够考虑星敏感器噪声，且操作简单，成本低廉的折射星的模拟方法。实现本专利技术目的技术方案：步骤1：扫描星表，判断当前恒星是否为折射星；若判断为折射星,则进入步骤2,否则继续扫描；步骤2：计算星光折射角；步骤3：求解折射星视位置的赤经、赤纬；步骤4：判断折射星的折射星光是否被星敏感器捕捉到；若判断折射星的折射星光被星敏感器捕捉到,则进入步骤5,否则返回步骤1；步骤5：进行星点成像中心计算并生成模拟折射星图。步骤1中,符合下列方程的恒星判断为折射星，式中，s为星表中记录的恒星星光的方向矢量；θR为折射高度为20km时对应的星光折射角，可通过大气折射模型计算得到；r＝[xyz]为载体在地心赤道惯性坐标系的位置矢量，为其模长，代表载体到地心的距离；α和β可表示为，式中，ha＝20km，hb＝50km分别为根据平流层厚度所选择的最小折射高度和最大折射高度；Re为地球半径。步骤2中，星光折射角通过解如下方程组获得，式中，ha为折射视高度；u＝|r·u|＝|xsx+ysy+zsz|，u＝[sxsysz]T为折射前的星光矢量，sx、sy、sz为恒星在天球坐标系的方向矢量，可通过查找星表得到；R为星光折射角。步骤3中，恒星视位置的赤经、赤纬通过解如下方程组获得，式中，P＝[cosγcosδsinγcosδsinδ]T为折射星光在地心赤道坐标系的方向矢量，γ和δ分别为折射星视位置的赤经和赤纬。步骤5中，具体包括如下步骤，步骤5.1：星等到灰度的模拟；将恒星在星敏感器中成像的亮度用灰度表示，在基本星表中星等的范围是-1至11之间的浮点型数，取0到6星等的恒星，令d为观测星的星等，g为所成像的灰度值，则模拟恒星在像平面的灰度为：g＝100+20·(6-d)步骤5.2：星点成像中心模拟；通过如下公式计算折射星成像中心位置(xs,ys)，式中，Nx、Ny分别为星敏感器像平面两个方向的像素数目；θFOV为星敏感器的视场角；γ0、δ0分别星敏感器的视轴指向所对应的赤经、赤纬；γ和δ分别为折射星视位置的赤经和赤纬；步骤5.3：星象点模拟；星像能量分布公式如下，式中，σ代表星点扩散半径表示星点弥散斑的大小，(x,y)为像平面的像元坐标，f(x,y)为像元(x,y)对应的灰度值；步骤5.4：星图生成；根据星敏感器的参数设定好星点扩散半径、背景灰度和背景噪声后生成的数字星图可表示为，P(x,y)＝f(x,y)+B+NB(x,y)+N(x,y)式中，B为背景的灰度值，NB(x,y)为背景噪声，设置为均值为0方差为σB的高斯白噪声；N(x,y)为随机噪声，设置为均值为0，方差为σN的高斯白噪声。本专利技术具有的有益效果：本专利技术首先扫描星表，判断当前恒星是否为折射星；若判断为折射星,则计算星光折射角；求解折射星视位置的赤经、赤纬；然后，判断折射星的折射星光是否被星敏感器捕捉到；若判断折射星的折射星光被星敏感器捕捉到,则进行星点成像中心计算并生成模拟折射星图。本专利技术不但能够模拟折射星在星敏感器像平面的位置，而且能够根据实际情况考虑星敏感器噪声、环境噪声的影响，使得折射星的模拟和折射角的计算更加符合实际，并且操作简单，成本低廉。附图说明图1是折射星的判断原理图；图2是星光折射原理图；图3是恒星矢量与主轴平面夹角示意图。具体实施方式步骤1：扫描星表，判断当前恒星是否为折射星；若判断为折射星,则进入步骤2,否则继续扫描；图1中：ha-最小折射高度20km，hb-最大折射高度50km，r-载体位置矢量，Re-地球半径s1、s2-恒星星光的方向矢量，θR-折射高度为20km时的星光折射角。在使用折射星导航时，根据平流层的厚度折射高度一般选取20km-50km，即图1中ha＝20km，hb＝50km；假设某恒星星光的方向矢量为s，则由图1得符合下列方程的恒星可以被选为折射星其中，θR为折射高度为20km时对应的星光折射角，可通过大气折射模型计算得到；r＝[xyz]为载体在地心赤道惯性坐标系的位置矢量，为其模长，代表载体到地心的距离；α和β如图1所示，可表示为步骤2：计算星光折射角；图2中，R-星光折射角，ha-折射视高度，hg-折射高度，O-地心，u-星光折射前的方向矢量若某颗恒星利用式(1)判断为折射星之后，就计算该恒星的星光折射角R；如图2所示，根据折射高度与星光折射角之间的几何关系可得如下方程其中，u＝|r·u|＝|xsx+ysy+zsz|，u＝[sxsysz]T为折射前的星光矢量,sx、sy、sz为恒星在天球坐标系的方向矢量，可通过查找星表得到；同时根据大气折射模型也可以得到折射视高度与星光折射角之间本文档来自技高网...

【技术保护点】
一种折射星的模拟方法，其特征在于:步骤1：扫描星表，判断当前恒星是否为折射星；若判断为折射星,则进入步骤2,否则继续扫描；步骤2：计算星光折射角；步骤3：求解折射星视位置的赤经、赤纬；步骤4：判断折射星的折射星光是否被星敏感器捕捉到；若判断折射星的折射星光被星敏感器捕捉到,则进入步骤5,否则返回步骤1；步骤5：进行星点成像中心计算并生成模拟折射星图。

【技术特征摘要】
1.一种折射星的模拟方法，其特征在于:步骤1：扫描星表，判断当前恒星是否为折射星；若判断为折射星,则进入步骤2,否则继续扫描；步骤2：计算星光折射角；星光折射角通过解如下方程组获得，式中，ha为折射视高度；u＝|r·u|＝|xsx+ysy+zsz|，u＝[sxsysz]T为折射前的星光矢量，sx、sy、sz为恒星在天球坐标系的方向矢量，可通过查找星表得到；R为星光折射角；步骤3：求解折射星视位置的赤经、赤纬；步骤4：判断折射星的折射星光是否被星敏感器捕捉到；若判断折射星的折射星光被星敏感器捕捉到,则进入步骤5,否则返回步骤1；步骤5：进行星点成像中心计算并生成模拟折射星图；具体包括如下步骤，步骤5.1：星等到灰度的模拟；将恒星在星敏感器中成像的亮度用灰度表示，在基本星表中星等的范围是-1至11之间的浮点型数，取0到6星等的恒星，令d为观测星的星等，g为所成像的灰度值，则模拟恒星在像平面的灰度为：g＝100+20·(6-d)步骤5.2：星点成像中心模拟；通过如下公式计算折射星成像中心位置(xs,ys)，式中，Nx、Ny分别为星敏感器像平面两个方向的像素数目；θFOV为星敏感器的视场角；γ0...

【专利技术属性】
技术研发人员：钱华明，孙龙，钱林琛，蔡佳楠，黄蔚，沈忱，吴永慧，李浩，李通，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江;23