发明名称 基于差分卫星导航测量的跟踪转台定姿与指令补偿方法
摘要 本发明涉及一种基于卫星导航系统的跟踪转台快速定姿与旋转指令生成的方法,该方法解决了复杂室外环境下当地基准难以确定的问题。在转台上安装两根垂直的长杆,利用长杆上安装的卫星导航系统接收天线进行转台旋转中心的位置确定和转台平台的方向确定。结合目标的状态信息,分别用位置速度转换法和坐标转换矩阵法计算转台旋转指令,位置速度转换方法可以计算转台指令角速度信息,可用于跟踪动目标。坐标转换矩阵方法无法计算转台指令角速度信息,可用于瞄准静目标。
申请公布号 CN104596503A 申请公布日期 2015.05.06
申请号 CN201510036930.7 申请日期 2015.01.26
申请人 中国人民解放军国防科学技术大学 发明人 陈磊;白显宗;张国军;柴华
分类号 G01C21/00(2006.01)I;G01C25/00(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 湖南省国防科技工业局专利中心 43102 代理人 冯青
主权项 一种基于差分卫星导航测量的跟踪转台定姿与指令补偿方法,其特征在于,在转台上安装两根相互垂直的长杆,在长杆上安装卫星导航系统接收天线进行转台旋转中心的位置确定和转台平台的方向确定,分别用位置速度转换法和坐标转换矩阵法计算转台旋转指令,具体过程如下:一、平台定姿(一)长杆的安装在转台安装面上从方位旋转中心沿俯仰轴方向和垂直轴方向各安装一根长杆,定义这两个相互垂直的长杆的交点A和两个端点分别为B、C点,A点与转台的旋转中心<img file="697450dest_path_image001.GIF" wi="24" he="25" />不重合,有一个偏移量<img file="321067dest_path_image002.GIF" wi="15" he="20" />,两个长杆固连在转台的安装平面上,当转台的转动方位角和俯仰角分别为0时,这两个杆的方向确定平台的指向,A、B、C三点的位置确定转台转动中心<img file="617925dest_path_image001.GIF" wi="24" he="25" />的位置;(二)平台定向在A、B、C三点安装卫星导航系统载波相位测量系统,任意时刻测得三点的大地经度L、大地纬度B和大地高H数据,通过下式得到A、B、C三点地心地固坐标系ECF的位置矢量<img file="740733dest_path_image003.GIF" wi="26" he="26" />、<img file="185359dest_path_image004.GIF" wi="26" he="26" />和<img file="556209dest_path_image005.GIF" wi="26" he="26" />:<img file="645257dest_path_image006.GIF" wi="185" he="94" />(1)式中:<i>N</i>为地球椭球的卯酉圈半径,<i>e</i>为椭球的第一偏心率;<img file="375184dest_path_image007.GIF" wi="221" he="49" />(2)其中<i>a</i>和<i>b</i>分别为椭球的长半径和短半径;计算ECF坐标系A点到B、C点的相对位置矢量<img file="542992dest_path_image008.GIF" wi="33" he="26" />和<img file="649357dest_path_image009.GIF" wi="33" he="26" /><img file="711685dest_path_image010.GIF" wi="232" he="26" />(3)以A点为原点建立北东地坐标系<i>G</i><sub>A</sub>,将A到B、C的相对位置矢量<img file="612514dest_path_image008.GIF" wi="33" he="26" />和<img file="595514dest_path_image009.GIF" wi="33" he="26" />转换到<i>G</i><sub>A</sub>坐标系,得<img file="239990dest_path_image011.GIF" wi="38" he="26" />和<img file="100368dest_path_image012.GIF" wi="38" he="26" />:<img file="63776dest_path_image013.GIF" wi="225" he="63" />(4)由<img file="585936dest_path_image011.GIF" wi="38" he="26" />和<img file="158738dest_path_image012.GIF" wi="38" he="26" />得北东地坐标系中转台坐标系三个轴的单位矢量为<img file="437403dest_path_image014.GIF" wi="317" he="53" />(5)得北东地坐标系到转台坐标系的方向余弦阵,也即坐标转换矩阵为<img file="945614dest_path_image015.GIF" wi="162" he="30" />(6)而由北东地坐标系到转台坐标系的转换矩阵,写为北东地坐标系到平台坐标系、再由平台坐标系到转台坐标系的转换矩阵的乘积<img file="11528dest_path_image016.GIF" wi="149" he="26" />(7)因此北东地坐标系到平台坐标系的转换矩阵写为<img file="890202dest_path_image017.GIF" wi="249" he="28" />(8)而平台坐标系到转台坐标系的转换矩阵为<img file="662855dest_path_image018.GIF" wi="168" he="54" />(9)将式(9)代入式(8),得<img file="968066dest_path_image019.GIF" wi="218" he="28" />(10)至此,得到北东地坐标系到平台坐标系的转换矩阵,也即方向余弦阵,将矩阵写为元素的形式<img file="255696dest_path_image020.GIF" wi="166" he="76" />(11)北东地坐标系到平台坐标系的转换矩阵<img file="107984dest_path_image021.GIF" wi="49" he="26" />与平台方向角:方位角<img file="626821dest_path_image022.GIF" wi="22" he="25" />、俯仰角<img file="279431dest_path_image023.GIF" wi="22" he="25" />、滚动角<img file="319937dest_path_image024.GIF" wi="22" he="25" />的关系为<img file="211801dest_path_image025.GIF" wi="553" he="96" />(12)由式(11)和式(12)得平台俯仰角为<img file="83680dest_path_image026.GIF" wi="229" he="46" />(13)定义函数<img file="901332dest_path_image027.GIF" wi="247" he="241" />(14)函数<img file="992916dest_path_image028.GIF" wi="53" he="28" />的值域为<img file="118829dest_path_image029.GIF" wi="58" he="25" />;如果得到的平台俯仰角<img file="408996dest_path_image030.GIF" wi="121" he="26" />,则平台方位角和滚动角为<img file="600812dest_path_image031.GIF" wi="337" he="52" />(15)如果得到的平台俯仰角<img file="412648dest_path_image032.GIF" wi="77" he="26" />,则有<img file="82533dest_path_image033.GIF" wi="177" he="28" />(16)在已经得知转台此时转动角度<img file="872546dest_path_image034.GIF" wi="17" he="15" />和<img file="595783dest_path_image035.GIF" wi="17" he="22" />的情况下,求解得到平台定向参数<img file="832598dest_path_image022.GIF" wi="22" he="25" />、<img file="368490dest_path_image023.GIF" wi="22" he="25" />和<img file="508616dest_path_image024.GIF" wi="22" he="25" />;(三)平台定位方法长杆的A点与转台转动中心<img file="104551dest_path_image001.GIF" wi="24" he="25" />点不重合,二者之间的距离<img file="908557dest_path_image002.GIF" wi="15" he="20" />由转台、长杆、天线的外形尺寸决定,通过图纸或测量得到,为已知量,转台坐标系中A点到<img file="484026dest_path_image001.GIF" wi="24" he="25" />点的相对位置矢量为<img file="977193dest_path_image036.GIF" wi="33" he="26" />表示为<img file="744030dest_path_image037.GIF" wi="125" he="30" />,由<img file="955437dest_path_image038.GIF" wi="49" he="26" />和A点的大地坐标得到转台坐标系到地心地固系的坐标转换矩阵为<img file="396914dest_path_image039.GIF" wi="145" he="26" />(17)因此在地心地固系中A点到<img file="750447dest_path_image001.GIF" wi="24" he="25" />点的相对位置矢量为<img file="688185dest_path_image040.GIF" wi="34" he="26" />表示为<img file="888353dest_path_image041.GIF" wi="118" he="26" />(18)故地心地固系中<img file="304160dest_path_image001.GIF" wi="24" he="25" />点的位置矢量为<img file="834236dest_path_image042.GIF" wi="109" he="26" />(19)将<img file="874699dest_path_image043.GIF" wi="28" he="26" />用大地坐标表示,即得到了转台转动中心的位置,完成了平台的定位;(四)多次定位定向方法由于A、B、C三点的位置测量存在误差,因此由前述得到的平台定位定向参数也存在误差,进行多次测量,得到多组数据,计算多组参数值,最终参数值取各次参数的平均值;二、转台旋转指令计算转台的基本功能是通过转动两个角度:方位角<img file="296584dest_path_image034.GIF" wi="17" he="15" />和俯仰角<img file="843978dest_path_image035.GIF" wi="17" he="22" />,使得安装在转台上的有效载荷在当地北东地坐标系下指向所需的方位角<img file="166244dest_path_image044.GIF" wi="22" he="25" />和俯仰角<img file="9566dest_path_image045.GIF" wi="24" he="25" />,北东地坐标系下的目标状态、转台旋转中心位置参数、平台定向参数均为已知量,所求的未知量是转台的旋转角度与旋转角速度,转台指令参数的计算方法有两种,分别为位置速度转换法和坐标转换矩阵法。
地址 410073 湖南省长沙市砚瓦池正街47号