CN114398755B - 一种弹性滤波器设计方法 - Google Patents
一种弹性滤波器设计方法 Download PDFInfo
- Publication number
- CN114398755B CN114398755B CN202111499268.0A CN202111499268A CN114398755B CN 114398755 B CN114398755 B CN 114398755B CN 202111499268 A CN202111499268 A CN 202111499268A CN 114398755 B CN114398755 B CN 114398755B
- Authority
- CN
- China
- Prior art keywords
- filter
- elastic
- transfer function
- designed
- design
- 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
- 238000013461 design Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000012546 transfer Methods 0.000 claims abstract description 37
- 238000001914 filtration Methods 0.000 claims abstract description 23
- 238000012937 correction Methods 0.000 claims abstract description 6
- 101150091511 glb-1 gene Proteins 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 9
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000012938 design process Methods 0.000 abstract description 2
- 230000007246 mechanism Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 2
- 229920001971 elastomer Polymers 0.000 description 2
- 239000000806 elastomer Substances 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种弹性滤波器设计方法,包括步骤:计算箭体传递函数Gjt、弹性振动传递函数Gtx、执行机构的传递函数Gzx、以及设计的校正网络传递函数Gjz;计算控制系统在固定弹性频率点omga0下的幅值裕度;固定弹性频率点omga0下的要求的设计幅值裕度;确定需要设计的滤波深度;计算此滤波器在控制系统下的滤波宽度;设计此弹性滤波器为连续域滤波器;计算控制回路中弹性频率点omga0处的幅值裕度;在箭体上下拉偏工况下,计算判断频率点的幅值裕度是否满足设计要求。本发明的设计方法通过实际应用,在滤波器设计过程中,知道具体某个频率点的滤波深度和滤波宽度的情况下,最多通过一次或者两次就可以设计出来。
Description
技术领域
本发明属于运载火箭姿态控制系统设计技术领域,具体为一种弹性滤波器设计方法。
背景技术
在运载火箭姿态控制系统中,稳定和抑制箭体的弹性振动是极其重要的一个任务。因为如果箭体的弹性振动不能很快衰减,而是等幅振荡甚至发散,这不仅会使弹体结构强度受到破坏,也对刚性弹体姿态运动的稳定性极为不利。如果在控制系统回路中引入惯组敏感到的弹性角速率信息后,经导航、制导计算得到的姿态角偏差信号中包含明显的的弹性分量,经校正网络计算得到的控制指令中亦包含高频弹性频率信息。执行机构具有饱和特性,尤其伺服机构的速度特性容易进入饱和状态,伺服机构进入非线性饱和工作状态后,会使刚体控制信号减小,严重时控制系统将不能有足够的偏转角、角速度来提供需要的控制力,随之导致姿态角偏差增大以至失稳。
目前解决箭体弹性振动的途径主要有两种:一是运用鲁棒控制理论设计出统一的控制器,二是采用滤波器实现对弹性振动进行抑制;在实际的工程应用中上具有广泛的形式就是设计一种陷波滤波器,稳定和抑制箭体的弹性振动。但是在实际进行陷波器设计时,根据相关专业提供的弹性参数,进行点对点的设计,如果弹性参数有上下拉偏的情况下,很难确定陷波器明确陷波器的深度和宽度,需要根据通过蒙特卡洛打靶仿真进行测试得出设计结果是否满足要求,反复设计,不断迭代,使得工程任务量增大,延缓甚至严重耽误设计进度。
发明内容
针对现有技术的缺陷,本发明提出一种弹性滤波器设计方法,设计方式简单快速,相当程度上减少设计时间从而加快研制设计进度。
为实现上述目的,本发明提供一种弹性滤波器设计方法,应用于运载火箭姿态控制系统中,稳定和抑制箭体的弹性振动,所述运载火箭是轴对称形状,俯仰通道和偏航通道的姿态控制方式和弹性滤波器设计方式相同,其特征在于包括如下步骤
步骤一:挑选需要设计的飞行时间点T和所述飞行时间点T下箭体的弹性振动频率omga0(rad/s),根据确定的所述飞行时间点T和所述飞行时间点T对应的总体参数数据,计算箭体传递函数Gjt、弹性振动传递函数Gtx、执行机构的传递函数Gzx、以及设计的校正网络传递函数Gjz;
步骤二:计算系统的开环传递函数:Gkh=(Gjt+Gtx)×Gzx×Gjz。通过开环传递函数,计算控制系统在固定弹性频率点omga0下的幅值裕度,记为GmMargin0(dB);
步骤三:根据弹性频率拉偏范围,上下限分别记为omga1(rad/s)和omga2(rad/s),固定弹性频率点omga0下的要求的设计幅值裕度为GmMargin1(dB);
步骤四:确定需要设计的滤波深度为,Rs=|GmMargin1-GmMargin0|;
步骤五:设定控制系统的控制步长为Tfk,根据弹性频率上下拉偏范围,计算此滤波器在控制系统下的滤波宽度omgaP,
omgaP=[omga0/2π、omga1/2π]×Tfk×2
步骤六:通过步骤四和步骤五得到的滤波深度和宽度,设计此弹性滤波器为连续域滤波器Glb2(s);
步骤七:将连续域中的滤波器传递函数Glb2带入系统的开环控制函数Gkh_All=(Gjt+Gtx)×Gzx×Gjz×Glb2,计算控制回路中弹性频率点omga0处的幅值裕度;
步骤八:利用步骤七设计好的连续域滤波器Glb2(s),分别在箭体上下拉偏工况下,计算判断频率点的幅值裕度是否满足设计要求;若存在上下拉偏工况下不满足系统幅值裕度要求的情况时,回到步骤四和步骤五,调整增大滤波深度和宽度,进行后续的操作设计和判断。
进一步地,所述步骤六中首先设计的弹性滤波器为二阶的切比雪夫带阻滤波器,通过设置带阻滤波器的滤波宽度和滤波深度,设计出该离散域滤波器Glb1(z)的分子和分母的系数,所述滤波器的离散域滤波器传递函数为:
其中A1、A2、A3、A4、A5和B1、B2、B3、B4、B5分别为离散域滤波器传递函数Glb1(z)分子、分母的系数,Z表示离散的采样点;*号为乘号;
将离散域滤波器Glb1(z)转换为连续域滤波器Glb2(s),Z表示连续的采样点。
进一步地,所述双线性变化的方法具体为:
将离散域滤波器传递函数Glb1(z)中的Z用替换进行计算,转换成连续域下的传递函数Glb2(s);其中T为采样时间,大小与系统的控制步长Tfk相同;*号为乘号。
进一步地,所述步骤八中适当调整滤波深度和宽度的原则和方法是:
对于运载火箭而言,控制系统在额定工况状态下的幅值裕度大于等于6dB,上下拉偏工况下幅值裕度大于等于3dB;当设计出来的滤波器在弹性频率点处的幅值裕度小于指标要求时,增大增大滤波深度和宽度即可完成目标要求。
本发明与传统的设计方法比,优点如下:
传统的设计方法,在频率点有上下拉偏时,滤波器设计需要多个陷波滤波器相乘的形式来实现,每一个陷波滤波器在设计时需要多次试凑和迭代过程才可实现;
本发明的设计方法通过实际应用,在滤波器设计过程中,知道具体某个频率点的滤波深度和滤波宽度的情况下,最多通过一次或者两次就可以设计出来。
总之,本发明的设计方法适应能力强,设计方式简单快速,明显减少设计时间并且加快了工程设计进度。
附图说明
图1为运载火箭姿态控制系统原理示意图。
具体实施方式
本发明提供一种弹性滤波器设计方法,应用于运载火箭姿态控制系统中,稳定和抑制箭体的弹性振动,所述运载火箭是轴对称形状,俯仰通道和偏航通道的姿态控制方式和弹性滤波器设计方式相同,下面以俯仰通道的设计方法为例,分步骤描述此弹性滤波器的设计方法。
首先介绍箭体姿态控制系统原理:将传感器惯侧组合的输出数据进行相应的数据处理以及姿态角偏差计算,姿态角偏差信号经本发明的弹性滤波器和校正网络计算之后,输出姿态稳定控制量,控制执行机构进行动作,产生需要的控制力和控制力矩,对箭体姿态进行控制。如图1所示。
本发明的弹性滤波器设计方法包括如下步骤
步骤一:挑选需要设计的飞行时间点T和所述飞行时间点T下箭体的弹性振动频率omga0(rad/s),根据确定的所述飞行时间点T和所述飞行时间点T对应的总体参数数据,计算箭体传递函数Gjt、弹性振动传递函数Gtx、执行机构的传递函数Gzx、以及设计的校正网络传递函数Gjz;
步骤二:计算系统的开环传递函数:Gkh=(Gjt+Gtx)×Gzx×Gjz。通过开环传递函数,计算控制系统在固定弹性频率点omga0下的幅值裕度,记为GmMargin0(dB);
步骤三:根据弹性频率拉偏范围,上下限分别记为omga1(rad/s)和omga2(rad/s),固定弹性频率点omga0下的要求的设计幅值裕度为GmMargin1(dB);
步骤四:确定需要设计的滤波深度为,Rs=|GmMargin1-GmMargin0|;
步骤五:设定控制系统的控制步长为Tfk,根据弹性频率上下拉偏范围,计算此滤波器在控制系统下的滤波宽度omgaP,
omgaP=[omga0/2π、omga1/2π]×Tfk×2
步骤六:通过步骤四和步骤五得到的滤波深度和宽度,设计此弹性滤波器为连续域滤波器Glb2(s);
所述步骤六中首先设计的弹性滤波器为二阶的切比雪夫带阻滤波器,通过设置带阻滤波器的滤波宽度和滤波深度,设计出该离散域滤波器Glb1(z)的分子和分母的系数,所述滤波器的离散域滤波器传递函数为:
其中A1、A2、A3、A4、A5和B1、B2、B3、B4、B5分别为离散域滤波器传递函数Glb1(z)分子、分母的系数,Z表示离散的采样点;*号为乘号;
将离散域滤波器Glb1(z)转换为连续域滤波器Glb2(s),Z表示连续的采样点。
所述步骤六中通过双线性变化将离散域转换为连续域滤波器Glb2(s),所述双线性变化的方法具体为:
将离散域滤波器传递函数Glb1(z)中的Z用替换进行计算,转换成连续域下的传递函数Glb2(s);其中T为采样时间,大小与系统的控制步长Tfk相同;*号为乘号。
步骤七:将连续域中的滤波器传递函数Glb2带入系统的开环控制函数Gkh_All=(Gjt+Gtx)×Gzx×Gjz×Glb2,计算控制回路中弹性频率点omga0处的幅值裕度;
步骤八:利用步骤七设计好的连续域滤波器Glb2(s),分别在箭体上下拉偏工况下,计算判断频率点的幅值裕度是否满足设计要求;若存在上下拉偏工况下不满足系统幅值裕度要求的情况时,回到步骤四和步骤五,调整增大滤波深度和宽度,进行后续的操作设计和判断。
所述步骤八中适当调整滤波深度和宽度的原则和方法是:
对于运载火箭而言,控制系统在额定工况状态下的幅值裕度大于等于6dB,上下拉偏工况下幅值裕度大于等于3dB;当设计出来的滤波器在弹性频率点处的幅值裕度小于指标要求时,增大滤波深度和宽度即可完成目标要求。
因为俯仰通道和偏航通道的姿态控制方式和弹性滤波器设计方式相同,在实际设计工作中需要设计出两个弹性滤波器,然后在控制系统中分别使用。
Claims (4)
1.一种弹性滤波器设计方法,应用于运载火箭姿态控制系统中,稳定和抑制箭体的弹性振动,所述运载火箭是轴对称形状,俯仰通道和偏航通道的姿态控制方式和弹性滤波器设计方式相同,其特征在于包括如下步骤
步骤一:挑选需要设计的飞行时间点T和所述飞行时间点T下箭体的弹性振动频率omga0(rad/s),根据确定的所述飞行时间点T和所述飞行时间点T对应的总体参数数据,计算箭体传递函数Gjt、弹性振动传递函数Gtx、执行机构的传递函数Gzx、以及设计的校正网络传递函数Gjz;
步骤二:计算系统的开环传递函数:Gkh=(Gjt+Gtx)×Gzx×Gjz;通过开环传递函数,计算控制系统在固定弹性频率点omga0下的幅值裕度,记为GmMargin0(dB);
步骤三:根据弹性频率拉偏范围,上下限分别记为omga1(rad/s)和omga2(rad/s),固定弹性频率点omga0下的要求的设计幅值裕度为GmMargin1(dB);
步骤四:确定需要设计的滤波深度为,Rs=|GmMargin1-GmMargin0|;
步骤五:设定控制系统的控制步长为Tfk,根据弹性频率上下拉偏范围,计算此滤波器在控制系统下的滤波宽度omgaP,
omgaP=[omga0/2π、omga1/2π]×Tfk×2
步骤六:通过步骤四和步骤五得到的滤波深度和宽度,设计此弹性滤波器为连续域滤波器Glb2(s);
步骤七:将连续域中的滤波器传递函数Glb2带入系统的开环控制函数Gkh_All=(Gjt+Gtx)×Gzx×Gjz×Glb2,计算控制回路中弹性频率点omga0处的幅值裕度;
步骤八:利用步骤七设计好的连续域滤波器Glb2(s),分别在箭体上下拉偏工况下,计算判断频率点的幅值裕度是否满足设计要求;若存在上下拉偏工况下不满足系统幅值裕度要求的情况时,回到步骤四和步骤五,适当调整滤波深度和宽度,进行后续的操作设计和判断。
2.根据权利要求1所述的弹性滤波器设计方法,其特征在于,所述步骤六中首先设计的弹性滤波器为二阶的切比雪夫带阻滤波器,通过设置带阻滤波器的滤波宽度和滤波深度,设计出该离散域滤波器Glb1(z)的分子和分母的系数,所述滤波器的离散域滤波器传递函数为:
其中A1、A2、A3、A4、A5和B1、B2、B3、B4、B5分别为离散域滤波器传递函数Glb1(z)分子、分母的系数,Z表示离散的采样点;*号为乘号;
将离散域滤波器Glb1(z)转换为连续域滤波器Glb2(s),Z表示连续的采样点。
3.根据权利要求1所述的弹性滤波器设计方法,其特征在于,所述步骤六中通过双线性变化将离散域转换为连续域滤波器Glb2(s),所述双线性变化的方法具体为:
将离散域滤波器传递函数Glb1(z)中的Z用替换进行计算,转换成连续域下的传递函数Glb2(s);其中T为采样时间,大小与系统的控制步长Tfk相同;*号为乘号。
4.根据权利要求1所述的弹性滤波器设计方法,其特征在于,所述步骤八中适当调整滤波深度和宽度的原则和方法是:
对于运载火箭而言,控制系统在额定工况状态下的幅值裕度大于等于6dB,上下拉偏工况下幅值裕度大于等于3dB;当设计出来的滤波器在弹性频率点处的幅值裕度小于指标要求时,增大滤波深度和宽度即可完成目标要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111499268.0A CN114398755B (zh) | 2021-12-09 | 2021-12-09 | 一种弹性滤波器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111499268.0A CN114398755B (zh) | 2021-12-09 | 2021-12-09 | 一种弹性滤波器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114398755A CN114398755A (zh) | 2022-04-26 |
CN114398755B true CN114398755B (zh) | 2024-04-26 |
Family
ID=81226363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111499268.0A Active CN114398755B (zh) | 2021-12-09 | 2021-12-09 | 一种弹性滤波器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114398755B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116643482B (zh) * | 2023-07-27 | 2023-10-20 | 航天科工火箭技术有限公司 | 一种运载火箭侧喷流姿态冗余控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107515612A (zh) * | 2017-10-20 | 2017-12-26 | 湖北航天技术研究院总体设计所 | 基于侧喷流控制的弹性振动抑制方法 |
CN113609581A (zh) * | 2021-07-30 | 2021-11-05 | 航天科工火箭技术有限公司 | 运载火箭弹性频率在线辨识的方法及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2759401T3 (es) * | 2016-09-23 | 2020-05-08 | Airbus Defence & Space Sau | Método para cancelar de manera adaptativa en tiempo real modos elásticos en señales discretas en el tiempo |
US10386858B2 (en) * | 2017-05-12 | 2019-08-20 | Bell Helicopter Textron Inc. | System and method for rotorcraft active lateral shuffle filter |
-
2021
- 2021-12-09 CN CN202111499268.0A patent/CN114398755B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107515612A (zh) * | 2017-10-20 | 2017-12-26 | 湖北航天技术研究院总体设计所 | 基于侧喷流控制的弹性振动抑制方法 |
CN113609581A (zh) * | 2021-07-30 | 2021-11-05 | 航天科工火箭技术有限公司 | 运载火箭弹性频率在线辨识的方法及存储介质 |
Non-Patent Citations (3)
Title |
---|
基于广义最小二乘法的空空导弹幅值滤波器设计方法;李良;李友年;陈星阳;;航空兵器;20200215(01);全文 * |
导弹姿态控制系统弹性振动抑制问题研究;宿常鹏;王雪梅;许哲;李骏霄;王兴龙;;弹箭与制导学报;20191231(05);全文 * |
雷达回波模拟器通道幅相误差校正方法仿真;吴一飙;章锦文;徐若曦;;计算机仿真;20200715(07);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114398755A (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xia et al. | Back-stepping sliding mode control for missile systems based on an extended state observer | |
CN107807657B (zh) | 一种基于路径规划的挠性航天器姿态自适应控制方法 | |
CN106078742B (zh) | 一种针对带有输出约束的柔性机械臂的振动控制方法 | |
CN103207568B (zh) | 一种抗舵机饱和的船舶航向自适应控制方法 | |
CN110716566B (zh) | 一种欠驱动无人艇的有限时间轨迹跟踪控制方法 | |
CN110687800B (zh) | 一种数据驱动自适应抗干扰控制器及其估计方法 | |
CN114398755B (zh) | 一种弹性滤波器设计方法 | |
CN109709978B (zh) | 一种高超声速飞行器制导控制一体化设计方法 | |
CN108762083B (zh) | 一种基于加速度观测器的自动控制系统 | |
CN107515612B (zh) | 基于侧喷流控制的弹性振动抑制方法 | |
CN108646570A (zh) | 一种改进极点配置的混沌轨迹跟踪方法 | |
CN109298636B (zh) | 一种改进的积分滑模控制方法 | |
JP7412799B2 (ja) | ロボットの末端経路輪郭誤差の制御方法 | |
CN109143866A (zh) | 一种考虑执行器受限问题的刚性飞行器自适应固定时间姿态跟踪控制方法 | |
CN113110048B (zh) | 采用hosm观测器的非线性系统输出反馈自适应控制系统和方法 | |
Theis et al. | Observer-based LPV control with anti-windup compensation: A flight control example | |
CN110308726B (zh) | 基于非线性反步法的欠驱动船舶航向控制方法 | |
CN108803325B (zh) | 永磁同步电机伺服系统鲁棒有限时间控制方法 | |
CN113419550B (zh) | 飞行器执行机构输出力矩饱和抗干扰姿态控制方法及系统 | |
CN109814583A (zh) | 基于自主水下机器人航向运动的动态状态反馈控制方法 | |
CN115793449B (zh) | 直升机飞行姿态控制器设计方法及装置、设备、存储介质 | |
CN111752158A (zh) | 一种有限时间收敛的二阶滑模控制方法 | |
Fang et al. | Performance identification and compensation of simulator motion cueing delays | |
CN110209197B (zh) | 一种飞行器控制系统设计方法 | |
CN115248554A (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 |