CN110765404B - 一种机弹干扰辨识算法 - Google Patents
一种机弹干扰辨识算法 Download PDFInfo
- Publication number
- CN110765404B CN110765404B CN201910514373.3A CN201910514373A CN110765404B CN 110765404 B CN110765404 B CN 110765404B CN 201910514373 A CN201910514373 A CN 201910514373A CN 110765404 B CN110765404 B CN 110765404B
- Authority
- CN
- China
- Prior art keywords
- missile
- moment
- interference
- projectile
- bullet
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/02—Wind tunnels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
本发明公开了一种机弹干扰辨识算法,属于导弹飞行试验数据分析领域,其步骤包括:步骤一、对惯性测量装置获取的角速度信息进行预处理,包含单位变换及滤除高频噪声;步骤二、根据处理后的角速度信息,选取合适的时间间隔,解算角加速度信息;步骤三、根据角加速度信息和导弹转动惯量,求解弹体系下各轴总的力矩分量;步骤四、根据马赫数、合成攻角、气流滚转角、质心、舵偏角等信息,求解导弹自由流下弹体系下各轴的力矩分量;步骤五、根据弹体系下各轴总的力矩减去自由流下的力矩,得到弹体系下各轴机弹干扰力矩分量。本发明针对性地滤除高频信号,并尽量减小相位滞后,确保处理后的角速度信息可用于提取机弹干扰。
Description
技术领域
本发明涉及一种机载导弹飞行试验数据分析领域,具体涉及一种机弹干扰辨识算法,特别涉及内埋或外挂弹射发射机载导弹弹射分离段气动干扰辨识方法。
背景技术
为满足作战飞机挂点多样性需求和隐身性能需求,机载导弹的发射方式由传统的滑轨式向弹射式发展。由于导弹与载机之间距离较近,导弹与载机之间将会产生较强的干扰力和干扰力矩,形成较为复杂的流场环境。获取精确的机弹干扰,对实现导弹安全弹射分离尤为重要。
目前,主要通过CTS试验模拟机弹分离段流场环境,获取机弹干扰,用于机弹分离段稳定控制系统设计和验证。通过飞行试验数据分析,获取弹射分离段机弹干扰,有利于验证CTS试验气动数据有效性,并可用于弹射分离段稳定控制系统优化设计。但是,由于导弹弹射机构的作用,弹射分离段容易激发导弹一阶、二阶弹性振动,从而导致通过惯性测量组合得到的角速度信息包含了高频噪声信号,无法直接用于提取角加速度信息,造成无法从飞行试验数据中精确辨识机弹干扰。
发明内容
本发明的技术解决的问题是从飞行试验惯性测量装置中得到的包含高频噪声的角速度信息中,经过数据预处理后,用于辨识机弹干扰,从而获取机弹干扰。
本发明解决其技术问题所采用的技术方案是,涉及一种机弹干扰辨识算法,其步骤包括:
(1)对惯性测量装置获取的角速度信息进行预处理,包含单位变换及滤除高频噪声;
(2)根据步骤(1)处理后的角速度信息,选取合适的时间间隔,解算角加速度信息;
(3)根据步骤(2)得到的角加速度信息和导弹转动惯量,求解弹体系下各轴总的力矩分量;
(4)根据马赫数、合成攻角、气流滚转角、质心、舵偏角等信息,求解导弹自由流下弹体系下各轴的力矩分量;
(5)根据步骤(3)得到的总的力矩减去步骤(4)得到的自由流下的力矩,得到弹体系下各轴机弹干扰力矩分量。
进一步地,所述步骤(1)中,首先对将飞行试验数据中惯性测量装置测量的弹体系角速度除以57.3,将单位从°/s转换为rad/s。
进一步地,所述步骤(1)中,根据傅里叶变换方法得到的幅频特性曲线的特征,选用相应的滤波器滤除高频信号,具体为:如果有明显较强的特定频率的高频信号,可串联多个陷波滤波器(形如ω0为陷波滤波器中心频率,Q为品质因素,Q越大,滤波宽度越小,可串联多个陷波滤波器),滤除相应的高频信号;反之,可采用低通滤波器(形如/>T1为时间常数),滤除高频噪声。
进一步地,所述步骤(2)中,解算角加速度信息,所选取的时间间隔应不小于原始数据采样步长:
式中,Δt为飞行试验数据的采样步长,n≥1,一般n取3~6次,wx为步骤(1)预处理后弹体系滚转角速度,wy为步骤(1)预处理后弹体系偏航角速度,wz为步骤(1)预处理后弹体系俯仰角速度。为弹体系滚转角加速度,/>为弹体系偏航角加速度,/>弹体系俯仰角加速度。k为正整数,wx(k)、wy(k)、wz(k)表示第k步的姿态角速度。
进一步地,所述步骤(3)中,根据步骤(2)得到的角加速度信息和导弹转动惯量,求解弹体系下各轴总的力矩分量;
式中,Ix,Iy,Iz为导弹转动惯量。Mx、My、Mz为作用在导弹上的所有外力对质心的力矩在弹体系各轴上的分量。
进一步地,所述步骤(4)中,采用气动参数计算模型求解的力矩系数应坐标变换到弹体系:
式中,my_z1,mz_z1为作用在导弹上的不含干扰的力对质心的力矩系数在弹体旋转坐标系各轴的分量。my_z,mz_z为作用在导弹上的不含干扰的力对质心的力矩系数在弹体坐标系各轴的分量。为气流滚转角。
本发明的有益效果是:
(1)本发明根据惯性测量组合得到的角速度数据的幅频特性,选用陷波滤波器或低通滤波器,针对性地滤除高频信号,并尽量减小相位滞后,确保处理后的角速度信息可用于提取机弹干扰;
(2)本发明通过合理选取解算角加速度的时间间隔,提高干扰辨识精度。
附图说明
以下结合附图和实施例对本发明作进一步说明。
图1是本发明实施例的一种机弹干扰辨识流程图。
图2是坐标系示意图。
具体实施方式
以下对本发明提供的机弹干扰辨识算法的具体实施方式进行阐述:
首先,结合附图介绍涉及的坐标系和角度定义。
1、弹体坐标系oxyz
坐标原点o位于弹体质心,ox轴沿弹体纵轴,由尾部指向头部为正,oy轴在包含ox轴的导弹纵向对称面内且与ox轴垂直,指向上方为正,oz轴与ox轴和oy轴组成右手坐标系。
2、弹体旋转坐标系ox1y1z1
坐标原点o位于弹体质心,ox1轴沿弹体纵轴,由尾部指向头部为正,oy1轴与来流方向在弹体坐标系oyz平面内的投影重合,oz1轴与ox1轴和oy1轴组成右手坐标系。
3、合成攻角:速度矢量和ox轴的夹角;
4、气流滚转角Φ:oy1轴与oy轴夹角,沿弹体头部向尾部看,逆时针旋转为正。
首先说明利用飞行试验数据开展机弹干扰辨识所需的飞行试验数据和导弹相关信息:
(1)飞行试验数据:飞行时间、采样时间间隔、弹体系角速度(含俯仰、偏航、滚转角速度)、动压、马赫数、合成攻角、气流滚转角、舵偏角;
(2)导弹参数:导弹特征长度、导弹特征面积、导弹质心、导弹转动惯量;
(3)气动专业根据风洞试验建立的气动参数计算模型。
以下介绍具体实施步骤:
步骤一、对惯性测量装置获取的角速度信息进行预处理,包含单位变换及滤除高频噪声;
步骤一的具体实施过程为:
(1)将飞行试验数据中惯性测量装置测量的弹体系角速度除以57.3,将单位从°/s转换为rad/s。
(2)采用傅里叶变换方法,对弹射分离段角速度信息中包含的高频分量的频率进行辨识。(频谱分析为本领域公知常识,在此不再详述)
(3)根据傅里叶变换得到的幅频特性曲线,选取不同的滤波器滤除高频噪声:如果有明显较强的特定频率的高频信号(一般为弹体一阶、二阶弹性频率),采用陷波滤波器(形如ω0为陷波滤波器中心频率,Q为品质因素,Q越大,滤波宽度越小,可串联多个陷波滤波器),滤除相应的高频信号;反之,可采用低通滤波器(形如T1为时间常数),滤除高频噪声。
步骤二、根据步骤一处理后的角速度信息,选取合适的时间间隔(nΔt,Δt为采样间隔时间,n≥1,一般n取3~6次),解算角加速度信息;
式中,Δt为飞行试验数据的采样步长。wx为步骤(1)预处理后弹体系滚转角速度,wy为步骤(1)预处理后弹体系偏航角速度,wz为步骤(1)预处理后弹体系俯仰角速度。为弹体系滚转角加速度,/>为弹体系偏航角加速度,/>弹体系俯仰角加速度。k为正整数,wx(k)、wy(k)、wz(k)表示第k步的姿态角速度。
步骤三、根据步骤二得到的角加速度信息和导弹转动惯量,求解弹体系下各轴总的力矩分量;
式中,Ix,Iy,Iz为导弹转动惯量。Mx、My、Mz为作用在导弹上的所有外力对质心的力矩在弹体系各轴上的分量。
步骤四、根据马赫数、合成攻角、气流滚转角、质心、舵偏角等信息,求解导弹自由流下弹体系下各轴的力矩分量;
步骤四的计算需用到气动专业根据风洞试验建立的气动参数计算模型,气动参数计算模型的使用为本领域公知常识,在此不再详述,由于气动参数计算模型输出的力矩系数在弹体旋转坐标系下,需坐标变换到弹体系,按如下公式:
式中,my_z1,mz_z1为作用在导弹上的不含干扰的力对质心的力矩系数在弹体旋转坐标系各轴的分量。my_z,mz_z为作用在导弹上的不含干扰的力对质心的力矩系数在弹体坐标系各轴的分量。为气流滚转角。
根据自由流下的力矩系数和动压、特征长度、特征面积,计算导弹自由流下的力矩。
Mx_z=mx_zqSL
My_z=my_zqSL
Mz_z=mz_zqSL
式中:q为动压,S为导弹特征面积,L为导弹特征长度。Mx_z,My_z,Mz_z为作用在导弹上的不含干扰的力对质心的力矩在弹体坐标系各轴的分量。
步骤五、根据总的力矩减去自由流下的力矩,得到弹体系下各轴机弹干扰力矩分量。
Mx_n=Mx-Mx_z
My_n=My-My_z
Mz_n=Mz-Mz_z
式中,Mx_n,My_n,Mz_n为作用在导弹上的干扰力对质心的力矩在弹体坐标系各轴的分量。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。
Claims (3)
1.一种机弹干扰辨识方法,其特征在于:包含如下步骤:
(1)对惯性测量装置获取的角速度信息进行预处理,包含单位变换及滤除高频噪声;
(2)根据处理后的角速度信息,选取合适的时间间隔,解算角加速度信息;
(3)根据角加速度信息和导弹转动惯量,求解弹体系下各轴总的力矩分量;
(4)根据马赫数、合成攻角、气流滚转角、质心、舵偏角等信息,求解导弹自由流下弹体系下各轴的力矩分量;
(5)根据弹体系下各轴总的力矩减去自由流下的力矩,得到弹体系下各轴机弹干扰力矩分量;
所述步骤(1)中,根据傅里叶变换方法得到的幅频特性曲线的特征,选用相应的滤波器滤除高频信号,具体为:如果有明显较强的特定频率的高频信号,可串联多个陷波滤波器:ω0为陷波滤波器中心频率,Q为品质因素,Q越大,滤波宽度越小,可串联多个陷波滤波器),滤除相应的高频信号;反之,可采用低通滤波器:/>T1为时间常数,滤除高频噪声;
所述步骤(2)中,解算角加速度信息,所选取的时间间隔应不小于原始数据采样步长:
式中,Δt为飞行试验数据的采样步长,n≥1,n取3~6次;
wx为步骤(1)预处理后弹体系滚转角速度,wy为步骤(1)预处理后弹体系偏航角速度,wz为步骤(1)预处理后弹体系俯仰角速度;
为弹体系滚转角加速度,/>为弹体系偏航角加速度,/>弹体系俯仰角加速度;k为正整数,wx(k)、wy(k)、wz(k)表示第k步的姿态角速度。
2.根据权利要求1所述的一种机弹干扰辨识方法,其特征在于:所述步骤(3)中,根据步骤(2)得到的角加速度信息和导弹转动惯量,求解弹体系下各轴总的力矩分量;
式中,Ix,Iy,Iz为导弹转动惯量;Mx、My、Mz为作用在导弹上的所有外力对质心的力矩在弹体系各轴上的分量。
3.根据权利要求2所述的一种机弹干扰辨识方法,其特征在于:所述步骤(4)中,采用气动参数计算模型求解的力矩系数应坐标变换到弹体系:
式中,my_z1,mz_z1为作用在导弹上的不含干扰的力对质心的力矩系数在弹体旋转坐标系各轴的分量;my_z,mz_z为作用在导弹上的不含干扰的力对质心的力矩系数在弹体坐标系各轴的分量;为气流滚转角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910514373.3A CN110765404B (zh) | 2019-06-14 | 2019-06-14 | 一种机弹干扰辨识算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910514373.3A CN110765404B (zh) | 2019-06-14 | 2019-06-14 | 一种机弹干扰辨识算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110765404A CN110765404A (zh) | 2020-02-07 |
CN110765404B true CN110765404B (zh) | 2023-08-25 |
Family
ID=69329033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910514373.3A Active CN110765404B (zh) | 2019-06-14 | 2019-06-14 | 一种机弹干扰辨识算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110765404B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113609581B (zh) * | 2021-07-30 | 2024-04-19 | 航天科工火箭技术有限公司 | 运载火箭弹性频率在线辨识的方法及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105423823A (zh) * | 2015-10-29 | 2016-03-23 | 中国科学院长春光学精密机械与物理研究所 | 激光制导导弹的注入式半实物仿真系统 |
CN106125571A (zh) * | 2016-03-18 | 2016-11-16 | 中国运载火箭技术研究院 | 一种巡航导弹的速度控制建模方法 |
CN106643710A (zh) * | 2016-11-25 | 2017-05-10 | 北京航天自动控制研究所 | 基于动力学方程的飞行器惯性参数和测压孔压力预示方法 |
CN109612676A (zh) * | 2018-12-12 | 2019-04-12 | 北京空天技术研究所 | 基于飞行试验数据的气动参数反算方法 |
-
2019
- 2019-06-14 CN CN201910514373.3A patent/CN110765404B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105423823A (zh) * | 2015-10-29 | 2016-03-23 | 中国科学院长春光学精密机械与物理研究所 | 激光制导导弹的注入式半实物仿真系统 |
CN106125571A (zh) * | 2016-03-18 | 2016-11-16 | 中国运载火箭技术研究院 | 一种巡航导弹的速度控制建模方法 |
CN106643710A (zh) * | 2016-11-25 | 2017-05-10 | 北京航天自动控制研究所 | 基于动力学方程的飞行器惯性参数和测压孔压力预示方法 |
CN109612676A (zh) * | 2018-12-12 | 2019-04-12 | 北京空天技术研究所 | 基于飞行试验数据的气动参数反算方法 |
Non-Patent Citations (1)
Title |
---|
飞行器动力学信息辅助MEMS惯导系统;李烨等;《系统工程与电子技术》;20160602(第08期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110765404A (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105202972B (zh) | 一种基于模型预测控制技术的多导弹协同作战制导方法 | |
CN103471803B (zh) | 一种模型自由飞试验的气动参数确定方法 | |
CN102591212B (zh) | 一种时变测量延迟输出信号飞行器纵向运动状态观测方法 | |
CN107844128B (zh) | 一种基于复合比例导引的高超声速飞行器巡航段制导方法 | |
CN112836581B (zh) | 一种基于相关性分析的敏感故障特征提取方法及装置 | |
CN112836383B (zh) | 一种基于靶式数据修正的武器效能建模方法 | |
US8296053B1 (en) | System and method for determining relative motion between ship combat system elements | |
CN102680201A (zh) | 基于视频测量的抖振风洞试验方法 | |
CN109612676B (zh) | 基于飞行试验数据的气动参数反算方法 | |
CN111964862A (zh) | 一种机弹分离动力学相似风洞试验方法 | |
CN111141313B (zh) | 一种提高机载局部相对姿态匹配传递对准精度的方法 | |
CN110765404B (zh) | 一种机弹干扰辨识算法 | |
CN105300382B (zh) | 一种基于惯性测量单元的大角度船体变形测量方法 | |
CN107367941B (zh) | 高超声速飞行器攻角观测方法 | |
CN105180728A (zh) | 基于前数据的旋转制导炮弹快速空中对准方法 | |
CN112257186A (zh) | 针对小型四旋翼飞行器气动参数的时域辨识方法 | |
CN108318219B (zh) | 自由飞试验相似律重力补偿设计方法 | |
CN108225375A (zh) | 一种基于中值滤波的抗外速度野值的优化粗对准方法 | |
CN110487300A (zh) | 减振装置对惯性导航系统的性能影响测试方法 | |
CN110398342B (zh) | 飞机挂载物分离风洞自由飞试验相似律设计方法 | |
RU2408031C2 (ru) | Способ сопровождения пилотируемой воздушной цели | |
CN113359192A (zh) | 一种微弱磁异常目标检测及定位方法 | |
CN110533005B (zh) | 一种复杂海况下船体形变测量方法 | |
CN107991054B (zh) | 一种捕获轨迹实验方法 | |
CN106950982B (zh) | 再入飞行器姿控动力系统高空力矩特性辨识方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |