发明名称 一种捷联寻的导引头的制导信息提取方法
摘要 本发明提供了一种捷联寻的导引头的制导信息提取方法,包括如下步骤:(1)获取载体坐标系o-x<sub>b</sub>y<sub>b</sub>z<sub>b</sub>下导引头的体视线方位角和体视线高低角;(2)处理捷联导引头测量信息;(3)计算载体坐标系下弹目相对距离矢量以及目标相对飞行器的运动速度在载体坐标系的分量;(4)计算载体坐标系下当前时刻飞行器的惯性视线方位角速率和当前时刻飞行器的惯性视线高低角速率;(5)将步骤(4)获得的结果作为制导信息送至捷联寻的导引头。本发明针对现有技术的不足,直接使用体视线角测量信息和陀螺角速度测量信息提取惯性视线角速度,提高了制导精度,降低制导系统的设计难度,能够广泛应用于各种捷联寻的制导武器中。
申请公布号 CN105021092A 申请公布日期 2015.11.04
申请号 CN201510373110.7 申请日期 2015.06.30
申请人 北京航天长征飞行器研究所;中国运载火箭技术研究院 发明人 田源;黄朝东;方海红;吴昭辉;赵洋;李涛;鞠晓燕;孟刚;刘佳琪;刘志轩;苏晓东;闫新峰;秦雪;孙月光;赵春明;孙忠旭
分类号 F42B15/01(2006.01)I;G01C21/20(2006.01)I;F41G9/00(2006.01)I;F41G11/00(2006.01)I 主分类号 F42B15/01(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 陈鹏
主权项 一种捷联寻的导引头的制导信息提取方法,其特征在于,包括如下步骤:(1)获取载体坐标系o‑x<sub>b</sub>y<sub>b</sub>z<sub>b</sub>下导引头的体视线方位角q<sub>bh</sub>和体视线高低角q<sub>bv</sub>;所述导引头的体视线方位角q<sub>bh</sub>为弹目相对距离矢量<img file="FDA0000748967900000013.GIF" wi="47" he="65" />在ox<sub>b</sub>z<sub>b</sub>平面内的投影与ox<sub>b</sub>轴的夹角;所述导引头的体视线高低角q<sub>bv</sub>为弹目相对距离矢量<img file="FDA0000748967900000015.GIF" wi="48" he="65" />在ox<sub>b</sub>z<sub>b</sub>平面内的投影与弹目相对距离矢量<img file="FDA0000748967900000014.GIF" wi="47" he="69" />的夹角;所述的载体坐标系o‑x<sub>b</sub>y<sub>b</sub>z<sub>b</sub>的原点o位于导引头所在飞行器的质心,ox<sub>b</sub>轴沿飞行器纵向,oy<sub>b</sub>轴在飞行器纵向对称面内指向上方,oz<sub>b</sub>轴与ox<sub>b</sub>轴和oy<sub>b</sub>轴构成右手坐标系;所述的弹目相对距离矢量<img file="FDA00007489679000000110.GIF" wi="48" he="65" />为虚拟目标点与载体坐标系原点o的连线;所述的虚拟目标点为目标和飞行器质心之间的连线与以载体坐标系原点为中心的单位球面的交点;(2)对体视线方位角q<sub>bh</sub>及体视线高低角q<sub>bv</sub>进行滤波处理,获得导引头测量的当前时刻体视线方位角q<sub>bh,k</sub>的滤波值q<sub>bh,kf</sub>及导引头测量的当前时刻体视线高低角q<sub>bv,k</sub>的滤波值q<sub>bh,kf</sub>;计算当前时刻体视线方位角q<sub>bh,k</sub>及当前时刻体视线高低角q<sub>bv,k</sub>的导数并对计算结果进行滤波处理,获得当前时刻体视线方位角q<sub>bh,k</sub>的导数滤波值<img file="FDA0000748967900000016.GIF" wi="100" he="60" />及当前时刻体视线高低角q<sub>bv,k</sub>的导数滤波值<img file="FDA0000748967900000017.GIF" wi="105" he="70" />;(3)计算载体坐标系下弹目相对距离矢量<img file="FDA0000748967900000018.GIF" wi="48" he="65" />以及目标相对飞行器的运动速度在载体坐标系的分量;<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>r</mi><mi>x</mi></msub><mo>=</mo><mi>cos</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>h</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>r</mi><mi>y</mi></msub><mo>=</mo><mi>sin</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>r</mi><mi>z</mi></msub><mo>=</mo><mo>-</mo><mi>cos</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>h</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000748967900000011.GIF" wi="592" he="255" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>v</mi><mi>x</mi></msub><mo>=</mo><msub><mi>&omega;</mi><mrow><mi>y</mi><mn>1</mn></mrow></msub><msub><mi>r</mi><mi>z</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mrow><mi>z</mi><mn>1</mn></mrow></msub><msub><mi>r</mi><mi>y</mi></msub><mo>-</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mi>sin</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>h</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mi>b</mi><mi>h</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>h</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>v</mi><mi>y</mi></msub><mo>=</mo><msub><mi>&omega;</mi><mrow><mi>z</mi><mn>1</mn></mrow></msub><msub><mi>r</mi><mi>x</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mrow><mi>x</mi><mn>1</mn></mrow></msub><msub><mi>r</mi><mi>z</mi></msub><mo>+</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>v</mi><mi>z</mi></msub><mo>=</mo><msub><mi>&omega;</mi><mrow><mi>x</mi><mn>1</mn></mrow></msub><msub><mi>r</mi><mi>y</mi></msub><mo>-</mo><msub><mi>&omega;</mi><mrow><mi>y</mi><mn>1</mn></mrow></msub><msub><mi>r</mi><mi>x</mi></msub><mo>+</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mi>sin</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>h</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mi>b</mi><mi>h</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>v</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>b</mi><mi>h</mi><mo>,</mo><mi>k</mi><mi>f</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000748967900000012.GIF" wi="1510" he="260" /></maths>其中,r<sub>x</sub>,r<sub>y</sub>,r<sub>z</sub>为弹目相对距离矢量<img file="FDA0000748967900000019.GIF" wi="52" he="65" />在载体坐标系的位置分量,v<sub>x</sub>,v<sub>y</sub>,v<sub>z</sub>为目标相对飞行器的运动速度在载体坐标系的分量,ω<sub>x1,k</sub>,ω<sub>y1,k</sub>,ω<sub>z1,k</sub>为当前时刻飞行器的姿态角速度;(4)计算载体坐标系下当前时刻飞行器的惯性视线方位角速率<img file="FDA0000748967900000021.GIF" wi="76" he="73" />和当前时刻飞行器的惯性视线高低角速率<img file="FDA0000748967900000022.GIF" wi="109" he="74" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mi>h</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><msub><mi>r</mi><mi>z</mi></msub><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mi>r</mi><mi>x</mi></msub><msub><mi>v</mi><mi>z</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mi>v</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><msub><mi>r</mi><mi>x</mi></msub><msub><mi>v</mi><mi>y</mi></msub><mo>-</mo><msub><mi>r</mi><mi>y</mi></msub><msub><mi>v</mi><mi>x</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000748967900000023.GIF" wi="371" he="170" /></maths>(5)将步骤(4)获得的结果作为制导信息送至捷联寻的导引头。
地址 100076 北京市丰台区北京9200信箱76分箱6号