发明名称 一种船载无线电测量设备天线变形测量方法
摘要 一种船载无线电测量设备天线变形测量方法,涉及航天器姿态控制地面应用领域,解决现有技术无法实现对船载无线电测量设备自身变形测量的问题,本发明方法为:一、在船载无线电设备基座和三轴中心处分别安装星敏感器SS1,星敏感器SS2;二、船进坞坐墩时,无线电天线指向A=0°、E=90°并锁定,标定两星敏感器的安装矩阵;三、船航行时,锁定天线指向A=0°、E=90°,两星敏感器同时测星,并分别解算基座处和三轴中心处的船体姿态角;四、解算天线变形矩阵,求解天线等效变形角。本发明方法对惯导地平系下星敏感器姿态角进行了蒙气差修正,提高了变形角计算的精度。
申请公布号 CN103322969B 申请公布日期 2016.01.20
申请号 CN201310207175.5 申请日期 2013.05.30
申请人 中国科学院长春光学精密机械与物理研究所 发明人 郭敬明;何昕;魏仲慧;张同双;张磊;王子辰
分类号 G01B21/32(2006.01)I;G01C1/00(2006.01)I 主分类号 G01B21/32(2006.01)I
代理机构 长春菁华专利商标代理事务所 22210 代理人 陶尊新
主权项 一种船载无线电测量设备天线变形测量方法,其特征是,该方法包括以下步骤:步骤一、在船载无线电设备基座和三轴中心处分别安装星敏感器SS1和星敏感器SS2;步骤二、船进坞坐墩时,船载无线电设备天线指向方位角A=0°、俯仰角E=90°并锁定,船载标校经纬仪标定零位差和照准差,所述星敏感器SS1和星敏感器SS2分别标定主点、焦距以及光学畸变参数;船载标校经纬仪通过测星或者瞄方位标的方法确定测量船航向角K<sub>真</sub>,采用高精度水平仪标定全船水平基准纵摇角Ψ<sub>真</sub>、横摇角θ<sub>真</sub>,获得惯导地平系到甲板坐标系转换矩阵R<sub>0</sub>;两台星敏感器同时测星,分别获得星敏感器在J2000.0坐标系下姿态矩阵,对所述星敏感器SS1和星敏感器SS2在J2000.0坐标系下姿态矩阵经岁差、章动、自转、极移、船位修正,获得星敏感器在惯导地平系下姿态矩阵<img file="FDA0000761594910000011.GIF" wi="122" he="95" />对所述的姿态矩阵<img file="FDA0000761594910000012.GIF" wi="100" he="95" />经蒙气差修正,获得重构后的星敏感器在惯导地平系下姿态矩阵<img file="FDA0000761594910000013.GIF" wi="143" he="96" />计算两星敏感器的安装矩阵<img file="FDA0000761594910000014.GIF" wi="116" he="89" />所述两星敏感器的安装矩阵<img file="FDA0000761594910000015.GIF" wi="90" he="92" />为星敏感器坐标系到甲板坐标系下的转换矩阵:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>s</mi><mi>i</mi></mrow><mrow><mi>b</mi><mi>i</mi></mrow></msubsup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>M</mi><mrow><mi>g</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mi>s</mi><mi>i</mi></mrow></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>*</mo><msub><mi>R</mi><mn>0</mn></msub><mo>,</mo></mrow>]]></math><img file="FDA0000761594910000016.GIF" wi="446" he="94" /></maths>所述i=1,2;步骤三、船航行时,锁定船载无线电设备天线指向方位角A=0°、俯仰角E=90°,两星敏感器同时测星,获得船航行时星敏感器在J2000.0坐标系下姿态矩阵,对所述姿态矩阵经岁差、章动、自转、极移、船位修正,获得船航行时星敏感器在惯导地平系下的姿态矩阵<img file="FDA0000761594910000017.GIF" wi="110" he="98" />然后对姿态矩阵<img file="FDA0000761594910000018.GIF" wi="84" he="95" />经蒙气差修正,获得重构后的船航行时星敏感器在惯导地平系下姿态矩阵<img file="FDA0000761594910000019.GIF" wi="131" he="97" />根据步骤二获得的安装矩阵<img file="FDA00007615949100000110.GIF" wi="90" he="87" />和重构后的船航行时星敏感器在惯导地平系下姿态矩阵<img file="FDA00007615949100000111.GIF" wi="125" he="94" />计算船航行时船体姿态矩阵<img file="FDA00007615949100000112.GIF" wi="115" he="95" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>g</mi><mi>i</mi></mrow><mrow><mi>b</mi><mi>i</mi></mrow></msubsup><mo>=</mo><msubsup><mi>C</mi><mrow><mi>g</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mi>s</mi><mi>i</mi></mrow></msubsup><mo>*</mo><msubsup><mi>R</mi><mrow><mi>s</mi><mi>i</mi></mrow><mrow><mi>b</mi><mi>i</mi></mrow></msubsup><mo>,</mo></mrow>]]></math><img file="FDA0000761594910000021.GIF" wi="379" he="100" /></maths>其中,g为惯导地平系,b为甲板坐标系,s为星敏感器坐标系;步骤四、计算天线变形矩阵,获得天线等效变形角,实现对天线变形的测量;设定三轴中心处相对基座存在变形B,根据步骤三计算的星敏感器SS1和星敏感器SS2测量船体姿态矩阵分别为<img file="FDA0000761594910000022.GIF" wi="263" he="103" />具体关系为:<img file="FDA0000761594910000023.GIF" wi="259" he="100" />式中,B为等效变形矩阵:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><msubsup><mi>R</mi><msub><mi>g</mi><mn>1</mn></msub><msub><mi>b</mi><mn>1</mn></msub></msubsup><msup><mrow><mo>(</mo><msubsup><mi>R</mi><msub><mi>g</mi><mn>2</mn></msub><msub><mi>b</mi><mn>2</mn></msub></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>b</mi><mn>11</mn></msub></mtd><mtd><msub><mi>b</mi><mn>12</mn></msub></mtd><mtd><msub><mi>b</mi><mn>13</mn></msub></mtd></mtr><mtr><mtd><msub><mi>b</mi><mn>21</mn></msub></mtd><mtd><msub><mi>b</mi><mn>22</mn></msub></mtd><mtd><msub><mi>b</mi><mn>23</mn></msub></mtd></mtr><mtr><mtd><msub><mi>b</mi><mn>31</mn></msub></mtd><mtd><msub><mi>b</mi><mn>32</mn></msub></mtd><mtd><msub><mi>b</mi><mn>33</mn></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000761594910000024.GIF" wi="880" he="289" /></maths>则由变形矩阵B解算变形等效欧拉角:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><msub><mi>k</mi><mrow><mi>b</mi><mn>0</mn></mrow></msub><mo>=</mo><mi>t</mi><mi>a</mi><msup><mi>n</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>(</mo><msub><mi>b</mi><mn>31</mn></msub><mo>/</mo><msub><mi>b</mi><mn>11</mn></msub><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>&psi;</mi><mrow><mi>b</mi><mn>0</mn></mrow></msub><mo>=</mo><msup><mi>sin</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>(</mo><msub><mi>b</mi><mn>21</mn></msub><mo>)</mo></mtd></mtr><mtr><mtd><msub><mi>&theta;</mi><mrow><mi>b</mi><mn>0</mn></mrow></msub><mo>=</mo><mi>t</mi><mi>a</mi><msup><mi>n</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>(</mo><mo>-</mo><msub><mi>b</mi><mn>23</mn></msub><mo>/</mo><msub><mi>b</mi><mn>22</mn></msub><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000761594910000025.GIF" wi="485" he="254" /></maths>式中,k<sub>b0</sub>、ψ<sub>b0</sub>及θ<sub>b0</sub>分别为天线三轴中心处的变形角中的等效艏挠角、等效纵挠角和等效横扭角;天线三轴中心相对天线基座的变形角的计算公式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><msub><mi>k</mi><mi>b</mi></msub><mo>=</mo><mo>-</mo><msub><mi>k</mi><mrow><mi>b</mi><mn>0</mn></mrow></msub></mtd></mtr><mtr><mtd><mrow><msub><mi>&psi;</mi><mi>b</mi></msub><mo>=</mo><msub><mi>&psi;</mi><mrow><mi>b</mi><mn>0</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&theta;</mi><mi>b</mi></msub><mo>=</mo><msub><mi>&theta;</mi><mrow><mi>b</mi><mn>0</mn></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000761594910000026.GIF" wi="233" he="235" /></maths>式中,k<sub>b</sub>、ψ<sub>b</sub>和θ<sub>b</sub>分别为天线基座变形角中的艏挠角、纵挠角和横扭角。
地址 130033 吉林省长春市东南湖大路3888号