一种近地空间短波红外星敏感器流场气动热效应分析方法

本发明属于气动光学领域，涉及星敏感器在气动条件下的辐射分析，具体是一种近地空间短波红外星敏感器流场气动热效应分析方法。

背景技术：

星敏感器是以天球坐标系为参照系，以恒星为探测目标的高精度姿态测量仪器，主要由光学系统、图像传感器电路和控制与数据处理电路构成。通过图像传感器拍摄恒星图像，由数据处理电路对所拍摄的图像进行处理，提取出所需的星点位置和亮度信息，进行星图匹配，计算出星敏感器的三轴姿态，完成飞行器在空间惯性坐标系的姿态测量，为卫星、深空探测器等各类航天器提供高精度的姿态信息。由于作为光测设备的星敏感器受背景杂光的影响较大，所以以往星敏感器仅用于卫星、宇宙飞船、火箭等空间航天器的导航。传统星敏感器一般为可见光波段，在近地空间工作时受天空杂散光影响，难以满足全天时的应用需求，因此，越来越多的人开始研究全天时短波红外星敏感器。

目前的全天时短波红外星敏感器波长比传统的红外星敏感器更长，但是更易受热辐射影响，当应用于飞机、导弹等高速运动平台时，流场温度迅速升高，产生强烈的流场气动热辐射效应，会影响短波红外星敏感器正常工作。目前的星敏感器气动光学效应分析大多是针对可见光波段，可见光波段星敏感器波段短，受流场气动热辐射影响较小，但是分析气动光学效应时大多未考虑流场热效应的影响，导致分析星敏感器是否能探测到恒星时的准确率低。

技术实现要素：

本发明目的是为了解决目前的星敏感器气动光学效应分析未考虑流场气动热效应对星敏感器的影响导致的分析星敏感器是否能探测到恒星时的准确率低的问题，而提出了一种近地空间短波红外星敏感器流场气动热效应分析方法，具体过程为：

步骤一、根据恒星的坐标判断恒星是否在短波红外星敏感器视场范围内；

步骤二、获取在短波红外星敏感器视场范围内的恒星在2mass星表内的星等信息；

步骤三、利用短波红外色指数计算视场范围内恒星光谱；

步骤四、根据流场物理特性逐层计算流场吸收系数；

步骤五、根据流场吸收系数逐层迭代计算衰减的星光光谱辐射亮度和流场背景光谱辐射亮度；

步骤六、根据步骤三获得恒星光谱和步骤五获得的光谱辐射亮度和流场背景光谱辐射亮度计算视场内恒星信噪比；

步骤七、将获取的信噪比与短波红外星敏感器的探测极限信噪比进行比较，若获取的信噪比大于探测极限信噪比，则能够探测到该恒星，否则不能探测到该恒星。

本发明的有益效果为：

本发明考虑了气动光学效应时大多未考虑流场热效应对星敏感器的成像影响，基于2mass星表相关信息，结合逐线积分法计算流场光谱特性，并逐层迭代计算星光衰减能量与背景热辐射能量，进而分析全天时短波红外星敏感器在气动条件下的探测能力，增加了分析星敏感器是否能探测到恒星的准确率。

附图说明

图1为一种近地空间短波红外星敏感器流场气动热效应分析方法的流程图；

图2为2mass星表光谱响应曲线；

图3为建立的恒星短波红外色指数与恒星有效温度关系曲线；

图4为逐线积分法示意图。

具体实施方式

具体实施方式一：本实施方式一种近地空间短波红外星敏感器流场气动热效应分析方法具体过程为：

步骤一、根据恒星的坐标判断恒星是否在短波红外星敏感器视场范围内：

αi∈(αg-fovx/2cosδg,αg+fovx/2cosδg)

δi∈(δg-fovy/2,δg+fovy/2)

其中，(αi,δi)为第i颗恒星在天体坐标系下的坐标，(αg,δg)为星敏感器光轴指向；fovx和fovy分别表示x轴和y轴方向视场大小；

所述恒星的位置通过查阅星表获得；

步骤二、获取在短波红外星敏感器视场范围内的恒星在2mass星表内的星等信息：

步骤三、利用短波红外色指数计算视场范围内恒星光谱，具体过程为：

步骤三一、将恒星短波红外色指数信息转化为恒星有效温度，恒星短波红外色指数和恒星有效温度函数关系为：

其中，

其中，m1、m2为两个不同波段恒星星等，m1-m2是短波红外色指数，b(λ,t')为黑体光谱辐射亮度，t0_vega为织女星恒星有效温度，tstar恒星有效温度，s(λ)为2mass星表响应曲线，λ1是恒星波段上限，λ2是恒星波段下限，k取1、2，分别对应m1、m2星等所在波段，lk(t)为黑体光谱辐射亮度b(λ,t)与2mass星响应曲线s(λ)乘积对波长的积分，t'取t0_vega或tstar；

步骤三二、基于普朗克黑体辐射公式计算恒星光谱：

其中，em(λ)为该星的光谱辐射照度，m是星等，e0(λ)是0星等时光谱辐照度，m是恒星星等；

其中，

其中，m(λ,tstar)为恒星光谱辐射出射度，λ1波段上限，λ2波段下限，β是对应波段0星等恒星2mass星表的辐照度；

其中，

其中，h为普朗克常数，c为光速，k为玻尔兹曼常数，tstar为恒星有效温度，λ是恒星波段；

步骤四、根据流场物理特性逐层计算流场吸收系数，具体过程为：

步骤四一、建立近地空间全天时短波红外星敏感器工作平台的有限元分析模型，输入气动条件下流场的物理特性，获取每层流场温度场t、压强场信息p；

步骤四二、按视轴方向将流场划分为n层，每层长度为δs，根据每层流场温度t，计算温度t下的谱线强度积分：

si,j(t)为第i层介质第j条谱线分子的谱线积分强度，h为普朗克常数，k为玻尔兹曼常数，η0为谱线中心位置，e"为低态谱项，qv(t)为振动配分函数，qr(t)为转动配分函数，s(t0)是p0＝1.01325×10⁵pa、t0＝296k条件下的谱线强度，qv(t0)是p0＝1.01325×10⁵pa、t0＝296k条件下的振动配分函数，qr(t0)是p0＝1.01325×10⁵pa、t0＝296k条件下的转动配分函数，c为光速，p0是标准态的压强，t0是标准态的温度。

步骤四三、根据谱线积分强度逐线积分，计算水和二氧化碳吸收系数：

式中，κi,η为第i层介质吸收系数，f(η-ηj)为谱线线型函数，ηj为第j条谱线中心波数，n为分子数密度：

其中，

步骤五、逐层迭代计算衰减的星光光谱辐射亮度和流场背景光谱辐射亮度，具体过程为：

li,λ＝τi,λli-1,λ+(1-τi,λ)bi,λ

其中，li,λ是第i层流场背景光谱辐射亮度，是星光光谱辐射亮度，bi,λ为黑体辐射亮度，τi,λ为第i层流场气体透过率，li-1,λ第i-1层流场背景光谱辐射亮度，第i-1层星光光谱辐射亮度。

其中，τi,λ＝exp(-δsκi,η)

其中，δs为该层流场长度。

步骤六、计算视场内恒星信噪比大小：

其中，k0为星点中心像元占星点总能量的比值；为读出噪声；idark为暗电流；tint为曝光时间；sm是恒星星信号电子数，sbgk是背景辐射电子数。

其中，

式中，d为光学系统孔径；τopt(λ)为光学系统透过率；tint为曝光时间；wph为单个光子能量；q0为量子效率，ω是单位项元所占立体角，是第n层星光光谱辐射亮度，ln,λ是第n层流场背景光谱辐射亮度。

步骤七、将获取的信噪比与短波红外星敏感器的探测极限信噪比进行比较，若获取的信噪比大于探测极限信噪比，则能够探测到该恒星，否则不能探测到该恒星。