发明名称 消除全捷联导引头制导回路刻度尺系数误差的装置和方法
摘要 本发明公开了一种消除全捷联导引头制导回路刻度尺系数误差的装置和方法,所述装置包括输入模块、滤波模块、评估模块和补偿模块,通过输入模块在过载指令信息上叠加强迫抖动信号,弹体前向通道和姿态陀螺补偿通道在经过滤波模块过滤后的信号中提取视线角速率信号弹体姿态角速率信号,通过评估模块建立状态方程和量测方程以及卡尔曼滤波方程,通过卡尔曼滤波方程评估计算刻度尺系数误差的估计值,通过补偿模块将刻度尺系数误差的估计值和弹体姿态角速率信息相乘,并将乘得的补偿信息输入至含弹体扰动的惯性系下弹目视线角速率信息中,从而完成刻度尺系数误差的补偿,消除刻度尺系数误差对制导回路稳定性的影响,并显著提高制导精度,减小脱靶量。
申请公布号 CN104111078B 申请公布日期 2017.02.08
申请号 CN201410176948.2 申请日期 2014.04.29
申请人 北京理工大学;中国北方工业公司 发明人 王辉;杨哲;唐义平;郑多;田松;吴骏雄;李斌
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 北京康思博达知识产权代理事务所(普通合伙) 11426 代理人 路永斌;余光军
主权项 一种消除全捷联导引头制导回路刻度尺系数误差的装置,所述全捷联导引头制导回路包括:二次积分模块,其对接收到的目标加速度信息进行二次积分处理,并输出导弹与目标位置偏差信息;数学模块,其用于接收导弹与目标位置偏差信息,并输出弹目视线角信息;微分模块,其对弹目视线角信息进行微分计算得到弹目视线角速率信息;比例导引律模块,其接收弹目视线角速率信息,并将计算得到的过载指令信息输送至驾驶仪动力学模块;驾驶仪动力学模块,其输出导弹加速度响应信息;角速率陀螺动力学模块,其根据导弹加速度响应信息测量并输出弹体角速率信息;积分模块,其接收弹体角速率信息并对接收到的信息进行积分处理,得到弹体姿态角信息;其中,通过弹体前向通道提取弹目视线角速率信息,通过姿态陀螺补偿通道提取弹体姿态角速率信息;其特征在于,该装置包括:输入模块、滤波模块、评估模块和补偿模块,其中:输入模块用于输出强迫抖动信号,并将该强迫抖动信号叠加在由比例导引律模块输出的过载指令信息上,所述的强迫抖动信号为固定幅值和固定频率的正弦信号,其中,幅值大于或等于弹体需用过载,频率为15rad/s以上;滤波模块包括带通滤波器,带通滤波器设置于弹体前向通道和姿态陀螺补偿通道上,带通滤波器用于过滤进入弹体前向通道的信号和进入姿态陀螺补偿通道的信号,弹体前向通道和姿态陀螺补偿通道从经过带通滤波器过滤的信号中提取由强迫抖动信号引起的视线角速率信号和由强迫抖动信号引起的弹体姿态角速率信号;评估模块用于接收弹体前向通道和姿态陀螺补偿通道提取出的由强迫抖动信号引起的弹目视线角速率信息和由强迫抖动信号引起的弹体姿态角速率信息,跟据接收到的信息建立状态方程和量测方程,并根据状态方程和量测方程建立卡尔曼滤波方程,通过卡尔曼滤波方程计算得到刻度尺系数误差的估计值,把刻度尺系数误差的估计值输入至补偿模块,所述状态方程如下<maths num="0001"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mtd></mtr><mtr><mtd><mover><mi>q</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover></mtd></mtr><mtr><mtd><mrow><mover><mi>C</mi><mo>&CenterDot;</mo></mover></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>q</mi></mtd></mtr><mtr><mtd><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mtd></mtr><mtr><mtd><mi>C</mi></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>&mu;</mi><mi>s</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000498244930000021.GIF" wi="587" he="234" /></maths>量测方程如下<maths num="0002"><math><![CDATA[<mrow><msup><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>&theta;</mi></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>q</mi></mtd></mtr><mtr><mtd><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mtd></mtr><mtr><mtd><mi>C</mi></mtd></mtr></mtable></mfenced><mo>+</mo><msub><mi>&mu;</mi><mi>n</mi></msub></mrow>]]></math><img file="FDA0000498244930000022.GIF" wi="489" he="232" /></maths>卡尔曼滤波方程如下<maths num="0003"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mtd></mtr><mtr><mtd><mover><mi>q</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover></mtd></mtr><mtr><mtd><mover><mi>C</mi><mover><mo>^</mo><mo>&CenterDot;</mo></mover></mover></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>q</mi></mtd></mtr><mtr><mtd><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mtd></mtr><mtr><mtd><mover><mi>C</mi><mo>^</mo></mover></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>K</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>K</mi><mn>3</mn></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mo>*</mo></msup><mo>-</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>&theta;</mi></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>q</mi></mtd></mtr><mtr><mtd><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mtd></mtr><mtr><mtd><mover><mi>C</mi><mo>^</mo></mover></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000498244930000023.GIF" wi="1018" he="270" /></maths>其中,q为惯性系下弹目视线角,<img file="FDA0000498244930000024.GIF" wi="38" he="91" />为惯性系下弹目视线角速率,C为寄生回路刻度尺系数误差,μ<sub>s</sub>为系统过程噪声,<img file="FDA0000498244930000025.GIF" wi="41" he="96" />为含弹体扰动的惯性系下弹目视线角速率,θ为弹体姿态角,<img file="FDA0000498244930000026.GIF" wi="54" he="82" />为刻度尺系数误差的估计值,μ<sub>n</sub>为系统量测噪声,K<sub>1</sub>、K<sub>2</sub>、K<sub>3</sub>为卡尔曼滤波方程增益;补偿模块用于接收评估模块输出的刻度尺系数误差的估计值,并将刻度尺系数误差的估计值乘以弹体姿态角速率得到补偿信息,补偿模块将其计算得到的补偿信息反馈至含弹体扰动的惯性系下弹目视线角速率信息中。
地址 100081 北京市海淀区中关村南大街5号