CN105203113B - 一种线阵红外地球敏感器的两通道姿态解算方法 - Google Patents
一种线阵红外地球敏感器的两通道姿态解算方法 Download PDFInfo
- Publication number
- CN105203113B CN105203113B CN201510608148.8A CN201510608148A CN105203113B CN 105203113 B CN105203113 B CN 105203113B CN 201510608148 A CN201510608148 A CN 201510608148A CN 105203113 B CN105203113 B CN 105203113B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- msup
- mfrac
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
Abstract
本发明公开了一种线阵红外地球敏感器的两通道姿态解算方法，首先对相邻辐射数据进行判断，获得穿越位置的整数部分，进而计算其小数部分，然后通过单镜头线性扫描标定试验，线性拟合线性运动轨迹，将误差项正弦拟合以及二次拟合得到辐射量校正公式，进行各通道穿越位置的修正。继而通过获得轨道零位标定表，得到姿态角零位，在任意轨道、未知轨道应用相应模式解算姿态角。本发明在应用于不同轨道姿态测量背景下，基于提高线阵红外地球敏感器测量精度、降低测量偏差的算法上，损失部分精度的条件下，实现了硬件层面的冗余备份，提高在轨使用的可靠性。在精度要求不高的使用环境下，将镜头数量减至最少，最大程度降低设备体积和质量。
Description
技术领域
本发明涉及一种星载红外地球敏感器的姿态解算方法，可用于空间低轨、变轨航天器的在轨使用。
背景技术
红外地球敏感器，是基于地球红外辐射敏感原理的卫星姿态光学敏感器，可用于航天器相对于地球局地垂线的俯仰、滚动姿态角信号的测量、初始状态时航天器对地球的捕获和稳态运行时航天器的姿态控制。
根据红外地球敏感器内部是否含机械扫描机构，可分为扫描式和静态两类：其中扫描式又可分为圆锥扫描式(单圆锥、双圆锥)和摆动扫描式两种，而静态则分为线阵和面阵两种。如今，国内外已研发出多种类型的扫描式红外地平仪，并广泛用于空间，其精度已可达到0.07°(3σ)。
近年来，随着探测器的发展，我国已利用线阵和面阵焦平面红外探测器研制出两类静态红外地球敏感器，具有体积小、重量轻、无扫描机构等优点，并分别在小卫星、高轨卫星上得到应用。其中，面阵红外地球敏感器具有精度高的优势，可达到0.06°(3σ)，但相对研发成本较高，且主要适用于高轨卫星的；而已在轨应用的线阵红外地球敏感器虽然成本较低，但其精度低，测量理论偏差最大将达到0.6°(3σ＝0.5°)。
线阵红外地球敏感器中探测器位于光学系统的焦平面上，属于凝视型结构。当航天器运行于地球上空时，从太空航天器上观察地球时，得到相当于在4K冷背景中的一个平均亮温约为220K～240K的圆盘，圆盘的边缘称为地平圆。航天器运行于地球上空时，红外地球敏感器通过线列红外探测器检测地平圆的4个方位上14μm～16.25μm波段的地球红外辐射能量，确定线列阵红外探测器对应地平圆4个点的方位角位置，根据之间的几何关系，实现对卫星姿态的测量，得到航天器相对于地球当地垂线的俯仰角和滚动角。一般采用典型设计，按“X”结构对称排列四个探头(光学系统和探测器组成，探测器位于光学系统焦平面上)，滚动轴与星体飞行方向一致，而俯仰轴垂直与轨道面。A、B、C、D四个探头与滚动轴和俯仰轴成45°分布，相邻两个探头光轴夹角为90°。
目前，已有人提出提高其精度的方法，确实能够在一定程度上提高精度(3σ)，但未改变测量的绝对偏差。且因为航天型号产品需要高可靠性，产品内部一般采用反熔丝的FPGA芯片进行处理，使得这种方法在应用上存在较大难度：FPGA难以进行如此复杂的乘除运算，且将远远超出芯片使用容量。
同时，随着应用领域的进一步拓宽，航天器的飞行情况要求红外地球敏感器在轨能够适应于不同轨道、甚至未知轨道下的姿态测量。
发明内容
本发明的目的在于在应用于不同轨道姿态测量背景下，基于提高线阵红外地球敏感器测量精度、降低测量偏差的算法上，损失部分精度的条件下，实现了硬件层面的冗余备份，提高在轨使用的可靠性。该方法能够提高线阵列红外地球敏感器的在轨使用寿命，在精度要求不高的使用环境下，将镜头数量减至最少，最大程度降低设备体积和质量。
一种线阵红外地球敏感器的两通道姿态解算方法的处理步骤为：
(1)、获得为线阵列探测器A通道从当前第i元到第i-5元的辐射量数据灰度值，依次记为D6、D5、D4、D3、D2、D1，其中i是0-N的整数，N为线阵红外地球敏感器所用线阵列探测器的元素，N为大于6的整数。当前5元不足(i＜5)时，由上一周期最末像元依次补齐。
(2)、根据步骤(1)获得的当前位置i，按以下公式计算相邻辐射量差值diff，除法运算分子da和整数部分z。
当i＝0时，分别计算D4-D3、D5-D4、D6-D5：当D4-D3最大时，diff＝D4-D3，da＝D1+D6-2D3，z＝N-4；当D5-D4最大时，diff＝D5-D4，da＝D2+D6-2D4，z＝N-3；当D6-D5最大时，diff＝D6-D5，da＝D2+D6-2D4，z＝N-2；
当i＝1,2,3时，diff＝0，da＝0，z＝0；
当i＝4时，diff＝D4-D3，da＝D6-D3，z＝0；
当i＝5时，diff＝D4-D3，da＝D2+D6-2D3，z＝1；
当4＜i＜N时，diff＝D4-D3，da＝D1+D6-2D3，z＝i-4。
(3)、根据步骤(2)的计算结果，判断相邻辐射量差值diff首次出现最大值的情况，使用该情况下的相邻辐射量差值diff、除法运算分子da和整数部分z的数据。
(4)、根据步骤(3)的数据，按以下方法获得A通道穿越位置的整数部分Za：
若步骤(2)计算所得的除法运算分子da大于2倍的相邻辐射量差值diff，则穿越位置的整数部分Za的值为z-1，同时修正除法运算分子da为除法运算分子da与相邻辐射量差值diff的差。否则，穿越位置的整数部分Za的值即为z。
(5)、根据步骤(3)和(4)的计算结果，计算A通道穿越位置的小数部分Xa：小数计算的分子为除法运算分子da，分母为2倍的相邻辐射量差值diff，均采用M位有效数字，M为不大于32的正整数，并将分子、分母等倍扩大，使得分母最高位为1。
循环采用二分逼近的判断：如果分子值大于分母值的二分之一，则结果左移一位，舍最高位，并添最低位为1；否则，则结果左移一位，舍最高位，并添最低位为0。如此循环Y次进行二分逼近判断的方法，计算可得具有Y位有效数字的小数部分Xa，其中Y为不大于16的正整数。
(6)、针对四通道，进行单镜头线性扫描试验，获得各通道穿越位置的线性变化曲线图。
(7)、将步骤(6)中的曲线进行线性拟合，A通道姿态线性校准参数Ka的值即为其线性函数斜率均值的倒数；同时获得拟合值与实际值的差指，再作差值与实际值的变化曲线。
(8)、将步骤(7)中的曲线进行正弦函数拟合，拟合后再进行二次拟合修正。
(9)、由步骤(8)得到拟合函数fA(x)＝a1ax2+a2ax+a3a+b1a*sin(2π(x-b2a))，a1a、a2a、a3a分别为二次拟合的二次项系数、一次项系数和常数项系数，b1a、b2a分别为正弦函数拟合的幅度系数、相位系数，即为A通道辐射量校正公式。
(10)、由步骤(4)获得的整数部分满足大于1且小于(N-2)时，代入由步骤(9)获得辐射量校正公式，按公式计算，得到A通道穿越位置A’。
(11)、按步骤(1)-(10)，依次获得B通道穿越位置B’、C通道穿越位置C’、D通道穿越位置D’。
(12)、将线阵红外地球敏感器置于地球模拟系统中，设置地球模拟系统在线阵红外地球敏感器的工作轨道高度H下保持俯仰角和滚动角均为0，由步骤(1)-(11)获得的四通道穿越位置，依次对应即为各通道在当前轨道高度H下的穿越位置零位，依次对应记作A0、B0、C0、D0。
(13)、按步骤(12)，线阵红外地球敏感器的工作轨道要求，设置不同轨道高度，形成轨道零位标定表，该表包括不同轨道高度下，A、B、C、D通道穿越位置零位的值。
(14)、在任意轨道高度下，根据公式，计算零姿态穿越位置的理论值L，其值为l，其中h为线阵红外地球敏感器工作的轨道高度，单位为km，θ为地球敏感器每个像元的视场角，H0为线阵红外地球敏感器设计的标称轨道高度，单位为km。
(15)、根据步骤(13)获得的轨道零位标定表按以下方法进行，获得在理论值为l时，A、B、C、D通道穿越位置标定零位的值a0、b0、c0、d0。线阵红外地球敏感器的轨道零位标定表的内容为不同理论值L对应的A通道穿越位置标定零位A0、B通道穿越位置标定零位B0、C通道穿越位置标定零位C0和D通道穿越位置标定零位D0：
若由步骤(14)获得的l在该表中可查，则直接从表中获得；
若由步骤(14)获得的l在表中不可查，则选取不大于l的最小组数据，这组数据的理论值L和A、B、C、D通道穿越位置标定零位A0、B0、C0、D0的值分别记为l1、a1、b1、c1、d1，和不小于l的最大组数据，这组数据的理论值L和A、B、C、D通道穿越位置标定零位A0、B0、C0、D0的值分别记为l2、a2、b2、c2、d2。按以下公式，计算获得：
(16)、根据以下公式解算卫星姿态：
当设计只使用A、B通道时，采用2/4-AB模式，公式为
其中P为卫星姿态俯仰角，R为卫星姿态滚动角，θ为每个像元的视场角。
当设计只使用B、C通道时，采用2/4-BC模式，公式为
当设计只使用C、D通道时，采用2/4-CD模式，公式为
当设计只使用A、D通道时，采用2/4-AD模式，公式为
本专利的优点：
本专利在应用于不同轨道姿态测量背景下，基于提高线阵红外地球敏感器测量精度、降低测量偏差的算法上，损失部分精度的条件下，实现了硬件层面的冗余备份，提高在轨使用的可靠性。并可在精度要求不高的使用环境下，将镜头数量减至最少，最大程度降低设备体积和质量。
附图说明
图1为整体处理方法流程图。
具体实施方式
按照本发明所述方法，在某型号飞行器红外地球敏感器做了验证性实验，并结合地球模拟器进行误差评估。
验证性实验中，FPGA芯片采用A54SX72，其内容量为72000门，算法占资源为54.3％，各参量的取值如下：
| 文中及公式中的标号 | 取值或参数 |
| N | 16元 |
| θ | 2° |
| H0 | 500km |
| Y | 8位 |
| M | 12位 |
测试结果：
标定参数如下：
1.四通道辐射量校正公式
fA＝0.0030x2+0.9561x+0.1343-0.070sin(2π(x-0.42))；
fB＝0.0018x2+0.9739x+0.0772-0.070sin(2π(x-0.42))；
fC＝0.0028x2+0.9598x+0.1212-0.070sin(2π(x-0.42))；
fD＝0.0021x2+0.9697x+0.0924-0.078sin(2π(x-0.42))。
2.轨道零位标定表
3.其他参数
A通道姿态线性校准参数Ka＝1.038，
B通道姿态线性校准参数Kb＝1.038，
C通道姿态线性校准参数Kc＝1.038，
D通道姿态线性校准参数Kd＝1.038。
Claims (1)
1.一种线阵红外地球敏感器的两通道姿态解算方法,其特征在于包括以下步骤：
(1)、获得为线阵列探测器A通道从当前第i元到第i-5元的辐射量数据灰度值，依次记为D6、D5、D4、D3、D2、D1，其中i是0-N的整数，N为线阵红外地球敏感器所用线阵列探测器的元素，N为大于6的整数；当前5元不足，i＜5时，由上一周期最末像元依次补齐；
(2)、根据步骤(1)获得的当前位置i，按以下公式计算相邻辐射量差值diff，除法运算分子da和整数部分z；
当i＝0时，分别计算D4-D3、D5-D4、D6-D5：当D4-D3最大时，diff＝D4-D3，da＝D1+D6-2D3，z＝N-4；当D5-D4最大时，diff＝D5-D4，da＝D2+D6-2D4，z＝N-3；当D6-D5最大时，diff＝D6-D5，da＝D2+D6-2D4，z＝N-2；
当i＝1,2,3时，diff＝0，da＝0，z＝0；
当i＝4时，diff＝D4-D3，da＝D6-D3，z＝0；
当i＝5时，diff＝D4-D3，da＝D2+D6-2D3，z＝1；
当4＜i＜N时，diff＝D4-D3，da＝D1+D6-2D3，z＝i-4；
(3)、根据步骤(2)的计算结果，判断相邻辐射量差值diff首次出现最大值的情况，使用该情况下的相邻辐射量差值diff、除法运算分子da和整数部分z的数据；
(4)、根据步骤(3)的数据，按以下方法获得A通道穿越位置的整数部分Za：
若步骤(2)计算所得的除法运算分子da大于2倍的相邻辐射量差值diff，则穿越位置的整数部分Za的值为z-1，同时修正除法运算分子da为除法运算分子da与相邻辐射量差值diff的差；否则，穿越位置的整数部分Za的值即为z；
(5)、根据步骤(3)和(4)的计算结果，计算A通道穿越位置的小数部分Xa：小数计算的分子为除法运算分子da，分母为2倍的相邻辐射量差值diff，均采用M位有效数字，M为不大于32的正整数，并将分子、分母等倍扩大，使得分母最高位为1；
循环采用二分逼近的判断：如果分子值大于分母值的二分之一，则结果左移一位，舍最高位，并添最低位为1；否则，则结果左移一位，舍最高位，并添最低位为0；如此循环Y次进行二分逼近判断的方法，计算可得具有Y位有效数字的小数部分Xa，其中Y为不大于16的正整数；
(6)、针对四通道，进行单镜头线性扫描试验，获得各通道穿越位置的线性变化曲线图；
(7)、将步骤(6)中的曲线进行线性拟合，A通道姿态线性校准参数Ka的值即为其线性函数斜率均值的倒数；同时获得拟合值与实际值的差指，再作差值与实际值的变化曲线；
(8)、将步骤(7)中的曲线进行正弦函数拟合，拟合后再进行二次拟合修正；
(9)、由步骤(8)得到拟合函数fA(x)＝a1ax2+a2ax+a3a+b1a*sin(2π(x-b2a))，a1a、a2a、a3a分别为二次拟合的二次项系数、一次项系数和常数项系数，b1a、b2a分别为正弦函数拟合的幅度系数、相位系数，即为A通道辐射量校正公式；
(10)、由步骤(4)获得的整数部分满足大于1且小于(N-2)时，代入由步骤(9)获得辐射量校正公式，按公式计算，得到A通道穿越位置A’；
(11)、按步骤(1)-(10)，依次获得B通道穿越位置B’、C通道穿越位置C’、D通道穿越位置D’；
(12)、将线阵红外地球敏感器置于地球模拟系统中，设置地球模拟系统在线阵红外地球敏感器的工作轨道高度H下保持俯仰角和滚动角均为0，由步骤(1)-(11)获得的四通道穿越位置，依次对应即为各通道在当前轨道高度H下的穿越位置零位，依次对应记作A0、B0、C0、D0；
(13)、按步骤(12)，线阵红外地球敏感器的工作轨道要求，设置不同轨道高度，形成轨道零位标定表，该表包括不同轨道高度下，A、B、C、D通道穿越位置零位的值；
(14)、在任意轨道高度下，根据公式，计算零姿态穿越位置的理论值L，其值为l，其中h为线阵红外地球敏感器工作的轨道高度，单位为km，θ为地球敏感器每个像元的视场角，H0为线阵红外地球敏感器设计的标称轨道高度，单位为km；
(15)、根据步骤(13)获得的轨道零位标定表按以下方法进行，获得在理论值为l时，A、B、C、D通道穿越位置标定零位的值a0、b0、c0、d0；线阵红外地球敏感器的轨道零位标定表的内容为不同理论值L对应的A通道穿越位置标定零位A0、B通道穿越位置标定零位B0、C通道穿越位置标定零位C0和D通道穿越位置标定零位D0：
若由步骤(14)获得的l在该表中可查，则直接从表中获得；
若由步骤(14)获得的l在表中不可查，则选取不大于l的最小组数据，这组数据的理论值L和A、B、C、D通道穿越位置标定零位A0、B0、C0、D0的值分别记为l1、a1、b1、c1、d1，和不小于l的最大组数据，这组数据的理论值L和A、B、C、D通道穿越位置标定零位A0、B0、C0、D0的值分别记为l2、a2、b2、c2、d2，按以下公式，计算获得：
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>l</mi>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>l</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
(16)、根据以下公式解算卫星姿态：
当设计只使用A、B通道时，采用2/4-AB模式，公式为
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<msqrt>
<mn>2</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<mi>&theta;</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>b</mi>
</msub>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>R</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>2</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<mi>&theta;</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>b</mi>
</msub>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
其中P为卫星姿态俯仰角，R为卫星姿态滚动角，θ为每个像元的视场角，Kb为B通道姿态线性校准参数，Kc为C通道姿态线性校准参数，Kd为D通道姿态线性校准参数；
当设计只使用B、C通道时，采用2/4-BC模式，公式为：
<mrow>
<mi>P</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>2</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<mi>&theta;</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>b</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>C</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>c</mi>
</msub>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>R</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>2</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<mi>&theta;</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>C</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>c</mi>
</msub>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
当设计只使用C、D通道时，采用2/4-CD模式，公式为：
<mrow>
<mi>P</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>2</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<mi>&theta;</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>C</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>D</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>R</mi>
<mo>=</mo>
<mfrac>
<msqrt>
<mn>2</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<mi>&theta;</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>C</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>D</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
当设计只使用A、D通道时，采用2/4-AD模式，公式为：
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<msqrt>
<mn>2</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<mi>&theta;</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>D</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>R</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msqrt>
<mn>2</mn>
</msqrt>
<mn>2</mn>
</mfrac>
<mi>&theta;</mi>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>A</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>D</mi>
<mo>,</mo>
</msup>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mo>&rsqb;</mo>
<mo>.</mo>
</mrow>
3
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510608148.8A CN105203113B (zh) | 2015-09-22 | 2015-09-22 | 一种线阵红外地球敏感器的两通道姿态解算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510608148.8A CN105203113B (zh) | 2015-09-22 | 2015-09-22 | 一种线阵红外地球敏感器的两通道姿态解算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105203113A CN105203113A (zh) | 2015-12-30 |
| CN105203113B true CN105203113B (zh) | 2017-10-13 |
Family
ID=54950924
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510608148.8A Active CN105203113B (zh) | 2015-09-22 | 2015-09-22 | 一种线阵红外地球敏感器的两通道姿态解算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105203113B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106500730B (zh) * | 2016-11-30 | 2019-09-27 | 中国科学院上海技术物理研究所 | 基于线阵地球敏感器可变轨道地球模拟器的姿态模拟方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6026337A (en) * | 1997-09-12 | 2000-02-15 | Lockheed Martin Corporation | Microbolometer earth sensor assembly |
| CN101462599A (zh) * | 2008-12-15 | 2009-06-24 | 中国科学院上海技术物理研究所 | 静态红外地平仪地面检测用的新型地球模拟器 |
| CN102175247A (zh) * | 2011-01-14 | 2011-09-07 | 中国科学院上海技术物理研究所 | 提高线阵列静态红外地平仪姿态测量精度的方法 |
-
2015
- 2015-09-22 CN CN201510608148.8A patent/CN105203113B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6026337A (en) * | 1997-09-12 | 2000-02-15 | Lockheed Martin Corporation | Microbolometer earth sensor assembly |
| CN101462599A (zh) * | 2008-12-15 | 2009-06-24 | 中国科学院上海技术物理研究所 | 静态红外地平仪地面检测用的新型地球模拟器 |
| CN102175247A (zh) * | 2011-01-14 | 2011-09-07 | 中国科学院上海技术物理研究所 | 提高线阵列静态红外地平仪姿态测量精度的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 红外地球敏感器修正的卫星姿态角动量反馈控制;苏威等;《西安交通大学学报》;20140430;第48卷(第4期);第115-118页 * |
| 面向微小卫星的红外静态焦平面地球敏感器设计;沈国权等;《传感技术学报》;20120531;第25卷(第5期);第571-576页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105203113A (zh) | 2015-12-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105136171B (zh) | 基于线阵红外地球敏感器电信号产生装置的模拟方法 | |
| CN102175247B (zh) | 提高线阵列静态红外地平仪姿态测量精度的方法 | |
| CN105300385B (zh) | 一种线阵红外地球敏感器的四通道姿态解算方法 | |
| CN104807477B (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
| CN106840212A (zh) | 基于地面激光光斑质心位置的星载激光在轨几何检校方法 | |
| CN104459728B (zh) | 一种基于gnss定位的磁偏角校准方法 | |
| CN110220491A (zh) | 一种无人机的光学吊舱安装误差角估算方法 | |
| CN105446346B (zh) | 遥感卫星对月相对定标姿态调整方法 | |
| CN103983224A (zh) | 一种大尺度部件实测位姿拟合方法 | |
| CN105222790B (zh) | 一种线阵红外地球敏感器的三通道姿态解算方法 | |
| CN105444778A (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
| Swaczyna et al. | Distance to the IBEX ribbon source inferred from parallax | |
| CN105486315B (zh) | 遥感卫星对月绝对定标姿态调整方法 | |
| CN203551198U (zh) | 一种大型结构体转动惯量测量系统 | |
| CN105023281B (zh) | 基于点扩散函数波前修正的星点像质心计算方法 | |
| CN106873004A (zh) | 基于太阳高度角自适应的高轨面阵相机在轨几何检校方法 | |
| CN105203113B (zh) | 一种线阵红外地球敏感器的两通道姿态解算方法 | |
| CN103699121B (zh) | 一种基于解析冗余关系的卫星控制系统敏感器可诊断性确定方法 | |
| CN101251384B (zh) | 一种利用三轴转台测试紫外敏感器的方法 | |
| CN102607463A (zh) | 一种旋转反射体面天线经纬仪激光靶标测量方法 | |
| CN105842710A (zh) | 一种基于vrs差分原理的低轨双星时差频差精密修正方法 | |
| CN104976991A (zh) | 一种考虑卫星姿态变化的三线阵影像像方偏差的获取方法 | |
| CN109670204B (zh) | 运载火箭视频图像判读误差修正方法 | |
| CN106647704A (zh) | 一种geo轨道数传天线跟踪精度确定方法 | |
| CN104391311B (zh) | 基于gps广播数据的星上无源定位方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |