CN104270120B - 一种利用双惯组进行弹性振荡抑制的方法及系统 - Google Patents

一种利用双惯组进行弹性振荡抑制的方法及系统 Download PDF

Info

Publication number
CN104270120B
CN104270120B CN201410452874.0A CN201410452874A CN104270120B CN 104270120 B CN104270120 B CN 104270120B CN 201410452874 A CN201410452874 A CN 201410452874A CN 104270120 B CN104270120 B CN 104270120B
Authority
CN
China
Prior art keywords
omg0
omg1
omg
used group
function
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
Application number
CN201410452874.0A
Other languages
English (en)
Other versions
CN104270120A (zh
Inventor
朱伟
涂正光
罗志伟
范威
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Designing Institute of Hubei Space Technology Academy
Original Assignee
General Designing Institute of Hubei Space Technology Academy
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by General Designing Institute of Hubei Space Technology Academy filed Critical General Designing Institute of Hubei Space Technology Academy
Priority to CN201410452874.0A priority Critical patent/CN104270120B/zh
Publication of CN104270120A publication Critical patent/CN104270120A/zh
Application granted granted Critical
Publication of CN104270120B publication Critical patent/CN104270120B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

本发明公开了一种利用双惯组进行弹性振荡抑制的方法和系统，方法包括以下步骤：S1、双惯组分别采集运载器角速率信息；S2、根据获取的运载器角速率信息获取线性加权系数K；S3、根据运载器角速率信息，进行加权系数K的线性叠加，从而获取阻断了弹性振荡后的真实弹体角速率。本发明还提供了实现上述方法的系统。本发明方法和系统利用双惯组采集的不同信息，自动选取加权系数进行线性叠加，从而输出只有微弱弹性信息的惯组数据，进而阻断了弹性振荡对姿态控制系统的不利影响。本发明方法和系统具有算法简单，易于工程实现的优点。

Description

一种利用双惯组进行弹性振荡抑制的方法及系统
技术领域
本发明属于弹性振荡抑制技术领域，具体涉及到一种利用双惯组进行弹性振荡抑制的方法及系统。
背景技术
国内目前运载器弹性振荡抑制的方法一般均采用滤波器的措施，常使用的滤波器为低通滤波器和陷波滤波器。低通滤波器对硬件的相位特性要求很高，而陷波滤波器则要求弹性频率大于稳定系统剪切频率5倍以上。在此背景下，提出一种弹性振荡自适应控制方法，然而现有弹性振荡自适应控制方法在弹性信息过大的状态下弹性信息抑制能力较差，甚至会出现姿态控制系统发散的情况。
因此需要提出一种改进的弹性振荡抑制方法，在保持对弹性振荡阻断的能力的基础上，避免现有技术在弹性信息过大的状态下发散的缺陷。
发明内容
针对现有技术的缺陷和技术需求，本发明提供了一种利用双惯组进行弹性振荡抑制的方法及系统，可解决现有弹性振荡自适应控制方法在弹性信息过大的状态下弹性信息抑制能力差、甚至出现姿态控制系统发散的缺陷。
为实现上述目的，按照本发明，提供了一种利用双惯组进行弹性振荡抑制的方法，所述方法包括以下步骤：
S1、双惯组中第一惯组和第二惯组分别采集运载器角速率信息omg0和omg1；
S2、根据采集的运载器角速率信息omg0和omg1获取线性加权系数K：
K＝K0+int[C*sign(omg_in1)*omg_in2]，
其中，K0＝k00/(k00-k10)，k00为运载器发动机点火前第一惯组所处位置上的振型斜率；k10为发动机点火前第二惯组所处位置上的振型斜率；
函数int(x)表示对自变量x的积分，初始值为0；
函数sign(x)表示对自变量x取符号，当x为正值时，sign(x)＝1；当x为负值时，sign(x)＝-1；
C为收敛因子；
omg_in1＝omg0-omg1，
omg_in2＝omg0-K(omg0-omg1)-notch(omg0)；
函数notch(x)为陷波滤波器，其传递函数为
notch＝(S2+2*0.02*wn*S+wn2)/(S2+2*0.96*wn*S+wn2),
wn为运载器的固有频率；
S3、根据线性加权系数K获取阻断了弹性振荡后的真实弹体角速率omg＝omg0-(omg0-omg1)K。
相应地，本发明还提供了一种利用双惯组进行弹性振荡抑制的系统，所述系统包括双惯组、线性加权模块和真实弹体角速率获取模块：
所述双惯组中第一惯组和第二惯组分别用于采集运载器角速率信息omg0和omg1；
所述线性加权模块用于根据采集的运载器角速率信息omg0和omg1获取线性加权系数K：
K＝K0+int[C*sign(omg_in1)*omg_in2]，
其中，K0＝k00/(k00-k10)，k00为运载器发动机点火前第一惯组所处位置上的振型斜率；k10为发动机点火前第二惯组所处位置上的振型斜率；
函数int(x)表示对自变量x的积分，初始值为0；
函数sign(x)表示对自变量x取符号，当x为正值时，sign(x)＝1；当x为负值时，sign(x)＝-1；
C为收敛因子；
omg_in1＝omg0-omg1，
omg_in2＝omg0-K(omg0-omg1)-notch(omg0)；
函数notch(x)为陷波滤波器，其传递函数为
notch＝(S2+2*0.02*wn*S+wn2)/(S2+2*0.96*wn*S+wn2),
wn为运载器的固有频率；
所述真实弹体角速率获取模块用于根据线性加权系数K获取阻断了弹性振荡后的真实弹体角速率omg＝omg0-(omg0-omg1)K。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，主要具备以下的技术优点：本发明方法和系统利用双惯组采集的不同信息，进行加权系数的线性叠加从而输出只有微弱弹性信息的惯组数据，阻断了弹性振荡对飞行器姿态控制系统的不利影响。本发明方法具有算法简单，易于工程实现的优点。
附图说明
图1为本发明方法工作流程图；
图2为本发明方法控制原理框图；
图3(a)为利用本发明方法进行仿真得到的加权系数效果图；
图3(b)为利用本发明方法进行仿真得到的输出角速率效果图；
图3(c)为利用本发明方法进行仿真得到的残余弹性角速率效果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明通过利用双惯组进行弹性振荡抑制。根据双惯组采集的不同信息，进行不同加权系数的线性叠加从而输出只包含有微弱弹性信息的真实弹体角速率，从而阻断弹性振荡对姿态控制系统的不利影响。
如图1和图2所示，本发明弹性振荡抑制方法包括以下步骤：
S1、双惯组中惯组0和惯组1分别采集运载器角速率信息omg0和omg1；
S2、获取线性加权系数K：
K＝K0+int[C*sign(omg_in1)*omg_in2]，
其中，K0＝k00/(k00-k10)，k00是在运载器发动机点火前，惯组0所处位置上的振型斜率；k10是在发动机点火前，惯组1所处位置上的振型斜率；
函数int(x)表示对自变量x的积分，初始值为0；函数sign(x)表示对自变量x取符号，当x为正值时，sign(x)＝1；当x为负值时，sign(x)＝-1；
C为关于K的收敛因子，C越大，K收敛速度越快。C可根据不同的状态需要进行设计(如惯组的物理特性以及其在运载器中的安装位置)，根据经验值调整确定；
omg_in1＝omg0-omg1；
omg_in2＝omg0-K(omg0-omg1)-notch(omg0)；
函数notch()为陷波滤波器，其传递函数为
notch＝(S2+2*0.02*wn*S+wn2)/(S2+2*0.96*wn*S+wn2),
wn为运载器的固有频率；
S3、获取阻断了弹性振荡后的真实弹体角速率omg＝omg0-(omg0-omg1)K。
上述步骤S2中，加权系数K的初始设置值K0＝k00/(k00-k10)，k00和k10可以在地面直接利用仪器测得。陷波滤波器在使用前需要根据姿态控制系统的采样频率进行双线性变化。
运载器实际飞行过程中，令真实弹体角速率为omgb，弹性角速率为omgt，则其满足如下关系：omg0＝omgb+k0*omgt，omg1＝omgb+k1*omgt；k0和k1分别是惯组0和惯组1所处位置上的实时振型斜率，其中在实际飞行过程中k0和k1无法实时测量。本发明方法将omg0和omg1进行关于加权系数K的线性叠加，根据步骤S1-S3可得到阻断了弹性振荡后的真实弹体角速率omg。
以下结合一个具体实施例对本发明方法作进一步说明。
假设某运载器飞行过程中，相关仿真条件设定为k0＝0.1、k1＝-0.25、K0＝0、C＝0.5。弹性角速率omgt幅值为100deg/s，频率10Hz；真实弹体角速率omgb幅值为50deg/s，频率1Hz。
由图2可以看出，仿真实验中，惯组0采集运载角速率信息omg0＝50sin(2πt)+0.1*100sin(20πt),惯组1采集运载角速率信息omg1＝50sin(2πt)+(-0.25)*100sin(20πt),根据上述步骤S1-S3中公式可计算得到如图3所示的加权系数、输出角速率，并计算得到残余弹性角速率。
由图3(a)-3(b)可知，本发明方法输出的角速率与理论真实弹体角速率非常接近，且残余弹性角速率很小，因此可很好地对弹性振荡进行抑制。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种利用双惯组进行弹性振荡抑制的方法，其特征在于，所述方法包括以下步骤：
S1、双惯组中第一惯组和第二惯组分别采集运载器角速率信息omg0和omg1；
S2、根据采集的运载器角速率信息omg0和omg1获取线性加权系数K：
K＝K0+int[C*sign(omg_in1)*omg_in2]，
其中，K0＝k00/(k00-k10)，k00为运载器发动机点火前第一惯组所处位置上的振型斜率；k10为发动机点火前第二惯组所处位置上的振型斜率；
函数int(x)表示对自变量x的积分，初始值为0；
函数sign(x)表示对自变量x取符号，当x为正值时，sign(x)＝1；当x为负值时，sign(x)＝-1；
C为收敛因子；
omg_in1＝omg0-omg1，
omg_in2＝omg0-K(omg0-omg1)-notch(omg0)；
函数notch(x)为陷波滤波器，其传递函数为
notch＝(S2+2*0.02*wn*S+wn2)/(S2+2*0.96*wn*S+wn2),
wn为运载器的固有频率；
S3、根据线性加权系数K获取阻断了弹性振荡后的真实弹体角速率omg＝omg0-(omg0-omg1)K。
2.一种利用双惯组进行弹性振荡抑制的系统，其特征在于，所述系统包括双惯组、线性加权模块和真实弹体角速率获取模块：
所述双惯组中第一惯组和第二惯组分别用于采集运载器角速率信息omg0和omg1；
所述线性加权模块用于根据采集的运载器角速率信息omg0和omg1获 取线性加权系数K：
K＝K0+int[C*sign(omg_in1)*omg_in2]，
其中，K0＝k00/(k00-k10)，k00为运载器发动机点火前第一惯组所处位置上的振型斜率；k10为发动机点火前第二惯组所处位置上的振型斜率；
函数int(x)表示对自变量x的积分，初始值为0；
函数sign(x)表示对自变量x取符号，当x为正值时，sign(x)＝1；当x为负值时，sign(x)＝-1；
C为收敛因子；
omg_in1＝omg0-omg1，
omg_in2＝omg0-K(omg0-omg1)-notch(omg0)；
函数notch(x)为陷波滤波器，其传递函数为
notch＝(S2+2*0.02*wn*S+wn2)/(S2+2*0.96*wn*S+wn2),
wn为运载器的固有频率；
所述真实弹体角速率获取模块用于根据线性加权系数K获取阻断了弹性振荡后的真实弹体角速率omg＝omg0-(omg0-omg1)K。
CN201410452874.0A 2014-09-05 2014-09-05 一种利用双惯组进行弹性振荡抑制的方法及系统 Active CN104270120B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410452874.0A CN104270120B (zh) 2014-09-05 2014-09-05 一种利用双惯组进行弹性振荡抑制的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410452874.0A CN104270120B (zh) 2014-09-05 2014-09-05 一种利用双惯组进行弹性振荡抑制的方法及系统

Publications (2)

Publication Number Publication Date
CN104270120A CN104270120A (zh) 2015-01-07
CN104270120B true CN104270120B (zh) 2017-04-05

Family

ID=52161621

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410452874.0A Active CN104270120B (zh) 2014-09-05 2014-09-05 一种利用双惯组进行弹性振荡抑制的方法及系统

Country Status (1)

Country Link
CN (1) CN104270120B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101629969A (zh) * 2009-08-20 2010-01-20 北京航空航天大学 一种低精度光纤惯组输出误差的标定补偿和测试方法及装置
CN103925917A (zh) * 2014-05-05 2014-07-16 上海新跃仪表厂 一种运载火箭姿态角速率信号的测量系统及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2005045829A1 (ja) * 2003-11-11 2008-06-12 松下電器産業株式会社 フィルタ係数調整回路

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101629969A (zh) * 2009-08-20 2010-01-20 北京航空航天大学 一种低精度光纤惯组输出误差的标定补偿和测试方法及装置
CN103925917A (zh) * 2014-05-05 2014-07-16 上海新跃仪表厂 一种运载火箭姿态角速率信号的测量系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
光纤陀螺捷联惯组姿态算法改进研究;王修等;《科学技术与工程》;20140228;第14卷(第4期);第98-102页 *

Also Published As

Publication number Publication date
CN104270120A (zh) 2015-01-07

Similar Documents

Publication Publication Date Title
CN106067783B (zh) 基于粒子群算法的fir滤波器设计方法
CN107515612B (zh) 基于侧喷流控制的弹性振动抑制方法
CN107750027A (zh) 场景自适应主动降噪耳机及场景自适应主动降噪设计方法
CN104965760B (zh) 一种管理软件功能模块生命周期的方法和装置
CN108549106A (zh) 混叠噪声压制方法及装置
Chi et al. Privacy partition: A privacy-preserving framework for deep neural networks in edge networks
CN106341101B (zh) 群时延滤波器系数的计算方法以及滤波器的模拟系统
CN104270120B (zh) 一种利用双惯组进行弹性振荡抑制的方法及系统
CN110034961B (zh) 以ooda链为元体的渗率计算方法
CN105006233A (zh) 窄带前馈主动噪声控制系统及目标噪声抑制方法
Goodman et al. Delays and growth rates of multiple TEOAE components
CN111207915A (zh) 一种避共振调优的测试方法及系统
CN105187029B (zh) 一种基于ifx‑lms自适应算法的控制方法及装置
CN104410388A (zh) 基于非线性变换的自适应时间延迟估计方法
CN109389216A (zh) 神经网络的动态裁剪方法、装置及存储介质
WO2020057286A1 (zh) 波场正演模拟方法及装置
CN109543225B (zh) 车辆的控制方案生成方法、装置、存储介质和电子设备
CN106021972A (zh) 一种基于非参数核密度估计的风电功率数据修正方法
CN110363288A (zh) 一种神经网络的输入图像生成方法和系统
CN110798176B (zh) 一种任意波宽带信号预失真滤波器构建方法及数字滤波器
CN109150136A (zh) 一种功率信号中的有色噪声滤除方法及系统
CN108182933A (zh) 汽车降噪方法及汽车
CN109951174B (zh) 一种快速傅里叶变换滤波方法
CN108072899B (zh) 间断Galerkin有限元地震数值模拟算法的自适应实现方法
CN106130508A (zh) 基于fir滤波器的数字万用表降噪方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant