CN109543256B - 一种过载响应时间计算方法和系统 - Google Patents
一种过载响应时间计算方法和系统 Download PDFInfo
- Publication number
- CN109543256B CN109543256B CN201811322929.0A CN201811322929A CN109543256B CN 109543256 B CN109543256 B CN 109543256B CN 201811322929 A CN201811322929 A CN 201811322929A CN 109543256 B CN109543256 B CN 109543256B
- Authority
- CN
- China
- Prior art keywords
- sliding
- overload
- sliding mode
- time
- model
- 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
Images
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
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
本申请实施例提供了一种复合控制过载响应时间计算方法,采用滑模控制开关逻辑方法,首先计算到达滑模面所用时间,其次计算到达稳定所用时间,完成推导过载控制回路响应时间解析表达式,从而为设计合理的控制参数提供依据。
Description
技术领域
本申请涉及导弹飞行参数计算领域,特别涉及一种复合控制过载响应时间计算方法和系统。
背景技术
导弹在飞行过程中,由于高空气动力作用减弱,过载响应慢,为了减小中末制导交班误差,需要采用直接力与气动力复合控制,可以实现快速建立攻角,大幅提高过载响应速度。但是,对于导弹飞行过程中过载响应时间的分析与计算并没有合适的计算手段。
发明内容
为解决上述问题之一,本申请提供了一种过载响应时间计算方法和系统。
根据本申请实施例的第一个方面,提供了一种过载响应时间计算方法,方法的步骤包括:
步骤S3、基于从系统初始状态收敛到滑模面所用时间t1和在滑模面上滑动到达系统稳定点所用时间t2,获得过载响应总时间t。。
根据本申请实施例的第一个方面,提供了一种过载响应时间计算系统,该系统包括:
过载时间计算模块,基于从系统初始状态收敛到滑模面所用时间t1和在滑模面上滑动到达系统稳定点所用时间t2,获得过载响应总时间t。
本申请所述技术方案通过计算到达滑模面所用时间和到达稳定所用时间,完成推导过载控制回路响应时间解析表达式,从而为设计合理的控制参数提供依据。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1示出本方案所述过载响应时间计算方法的示意图;
图2示出本方案所述理论结果与仿真计算结果进行比对图。
具体实施方式
为了使本申请实施例中的技术方案及优点更加清楚明白,以下结合附图对本申请的示例性实施例进行进一步详细的说明,显然,所描述的实施例仅是本申请的一部分实施例,而不是所有实施例的穷举。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
本方案的核心思路是给出直气复合控制下过载响应时间解析表达式,通过该方法实现对过载回路控制参数的评价,将过载响应时间理论计算结果与仿真计算结果进行对比,证明过载响应时间计算方法的合理性。
如图1所示,本方案公开了一种复合控制过载响应时间计算方法的具体步骤为:
第一步到达滑模面所用时间
滑动模态公式:
s=(Nyc-Ny)+k1×(ωzc-ωz) (1)
式中,Nyc为过载指令,Ny为过载响应,k1为滑模控制参数,ωzc为俯仰角速度指令,ωz为俯仰角速度。
趋近律公式:
式中,k2为趋近律控制参数,T为轨控推力常数,sign(·)为符号函数。
设定在轨控直接力开机的条件下推导,即认为sign(s)=1,求解方程(2)可以得出滑动模态表达式:
式中,c为方程常值。
根据方程初值可以解出c值大小:
当t=0时,
式中,s0为滑动模态初值。
设定滑模边界层为0,当s到达0时,则
则根据上式得出到达滑模面所用时间t1:
第二步到达稳定所用时间
根据滑动模态公式,
s=(Nyc-Ny)+k1×(ωzc-ωz) (6)
假设(Nyc-Ny)服从以下关系:
式中,V为导弹速度,g为重力加速度,a4为气动力系数,αc为攻角指令,α为当前时刻攻角,s3为过载偏差。
近似认为
解得,
根据Ny=kNyc,解出:
式中,k5为常值增益,k为到达滑模面的过载指令倍数。
式中,t2为滑模到达稳定所用时间。
过载响应总时间:
t=t1+t2
通过该方法实现对过载控制参数的评价,要求过载响应时间理论计算值满足系统要求,可根据要求的过载响应时间反推出系统控制参数。
本方案进一步公开了一种过载响应时间计算系统,该系统包括:第一时间计算模块,基于滑动模态模型和趋近律模型,计算从系统初始状态收敛到滑模面所用时间t1:其中,k1为滑模控制参数,s0为滑动模态初值,k2为趋近律控制参数,T为轨控推力常数;
过载时间计算模块,基于从系统初始状态收敛到滑模面所用时间t1和在滑模面上滑动到达系统稳定点所用时间t2,获得过载响应总时间t。
本方案中,所述过载响应时间计算也可以通过电子设备实现其计算功能,所述电子设备包括:存储器,一个或多个处理器;存储器与处理器通过通信总线相连;处理器被配置为执行存储器中的指令;所述存储介质中存储有用于执行如上所述方法中各个步骤的指令。
本方案中,所述过载响应时间计算也可以记载于计算机可读存储介质中,通过计算机可读存储介质上存储有计算机程序实现计算功能,该程序被处理器执行时实现如上所述方法的步骤。
按照上述过载响应时间计算方法,给出过载响应时间的理论计算结果,设定参数如下:
根据公式计算得出基于从系统初始状态收敛到滑模面所用时间t1为0.0633s,在滑模面上滑动到达系统稳定点所用时间t2为0.0854s,从而得到过载响应总时间为0.1726s。
按照上述理论计算结果,对控制参数进行设计,将理论结果与仿真计算结果进行比对,根据图2,可得过载响应时间为0.175s,理论计算得到的响应时间为0.1726s,两者之间相对误差为1.37%,误差在允许范围之内,得到的仿真计算的过载响应时间与理论方法计算出来的结果吻合度较好,从而证明过载响应时间计算方法的合理性。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上仅为本发明的实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均包含在申请待批的本发明的权利要求范围之内。
Claims (7)
1.一种过载响应时间计算方法,其特征在于,该方法的步骤包括:
步骤S3、基于从系统初始状态收敛到滑模面所用时间t1和在滑模面上滑动到达系统稳定点所用时间t2,获得过载响应总时间t;
所述步骤S1包括:
预先构建滑动模态模型:s=(Nyc-Ny)+k1×(ωzc-ωz),其中,Nyc为过载指令,Ny为过载响应,k1为滑模控制参数,ωzc为俯仰角速度指令,ωz为俯仰角速度,s为滑动模态,以及,趋近律模型:其中,k2为趋近律控制参数,T为轨控推力常数,sign(·)为符号函数;
3.根据权利要求1所述的计算方法,其特征在于,所述过载响应总时间t:t=t1+t2。
6.根据权利要求4所述的过载响应时间计算系统,其特征在于,所述第二时间计算模块具体执行如下步骤:
预先构建滑动模态模型:s=(Nyc-Ny)+k1×(ωzc-ωz),其中,Nyc为过载指令,Ny为过载响应,k1为滑模控制参数,ωzc为俯仰角速度指令,ωz为俯仰角速度,s为滑动模态;
基于上述公式,得到在滑模面上滑动到达系统稳定点所用时间t2。
7.根据权利要求4所述的过载响应时间计算系统,其特征在于,所述过载响应总时间t为t=t1+t2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811322929.0A CN109543256B (zh) | 2018-11-08 | 2018-11-08 | 一种过载响应时间计算方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811322929.0A CN109543256B (zh) | 2018-11-08 | 2018-11-08 | 一种过载响应时间计算方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109543256A CN109543256A (zh) | 2019-03-29 |
CN109543256B true CN109543256B (zh) | 2023-01-20 |
Family
ID=65845364
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811322929.0A Active CN109543256B (zh) | 2018-11-08 | 2018-11-08 | 一种过载响应时间计算方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109543256B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104808492A (zh) * | 2015-03-23 | 2015-07-29 | 北京航天自动控制研究所 | 一种升力式飞行器的攻角指令生成方法 |
CN108363305A (zh) * | 2018-03-14 | 2018-08-03 | 中国人民解放军国防科技大学 | 基于主动干扰补偿的战术导弹鲁棒过载自驾仪设计方法 |
CN108536020A (zh) * | 2018-07-17 | 2018-09-14 | 哈尔滨工业大学 | 一种针对垂直起降重复使用运载器的模型参考自适应滑模控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106427957B (zh) * | 2015-08-11 | 2019-03-29 | 比亚迪股份有限公司 | 基于四轮驱动的电动汽车稳定控制系统及方法及电动汽车 |
-
2018
- 2018-11-08 CN CN201811322929.0A patent/CN109543256B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104808492A (zh) * | 2015-03-23 | 2015-07-29 | 北京航天自动控制研究所 | 一种升力式飞行器的攻角指令生成方法 |
CN108363305A (zh) * | 2018-03-14 | 2018-08-03 | 中国人民解放军国防科技大学 | 基于主动干扰补偿的战术导弹鲁棒过载自驾仪设计方法 |
CN108536020A (zh) * | 2018-07-17 | 2018-09-14 | 哈尔滨工业大学 | 一种针对垂直起降重复使用运载器的模型参考自适应滑模控制方法 |
Non-Patent Citations (1)
Title |
---|
飞行器制导与控制系统优化设计及弹道仿真;马娜娜;《中国优秀硕士学位论文全文数据库》;20180615;C032-66 * |
Also Published As
Publication number | Publication date |
---|---|
CN109543256A (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Robust attitude controller design for miniature quadrotors | |
Zuo et al. | Adaptive trajectory tracking control of output constrained multi‐rotors systems | |
Lu et al. | Finite‐time fault‐tolerant control for rigid spacecraft with actuator saturations | |
Yang et al. | Self-tuning PID control design for quadrotor UAV based on adaptive pole placement control | |
CN108037662A (zh) | 一种基于积分滑模障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法 | |
You et al. | Comprehensive design of uniform robust exact disturbance observer and fixed‐time controller for reusable launch vehicles | |
Lungu et al. | Landing auto‐pilots for aircraft motion in longitudinal plane using adaptive control laws based on neural networks and dynamic inversion | |
Norton et al. | Adaptive fuzzy multi‐surface sliding control of multiple‐input and multiple‐output autonomous flight systems | |
CN112925221B (zh) | 基于数据回注的辅助驾驶闭环测试方法 | |
CN112180957A (zh) | 无人机的三余度仲裁切换方法、装置及计算机设备 | |
CN109543256B (zh) | 一种过载响应时间计算方法和系统 | |
CN110825115A (zh) | 一种飞机迎角和过载的极限限制控制方法 | |
CN110058601B (zh) | 无人直升机起降快速判定的方法、装置及电子设备 | |
Gu et al. | Agile flight for a quadrotor via robust geometry control: Theory and experimental verification | |
CN107942672B (zh) | 一种基于对称时不变障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法 | |
Li et al. | A finite time composite control method for quadrotor UAV with wind disturbance rejection | |
CN112148035B (zh) | 多无人机轨迹优化方法、装置、存储介质及计算机设备 | |
CN110362108B (zh) | 无人机接地控制方法、装置,存储介质及电子设备 | |
CN108107726A (zh) | 一种基于对称时变障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法 | |
CN104699918A (zh) | 一种计算舵偏角的方法及系统 | |
CN113485393B (zh) | 飞行设备的控制方法、装置、存储介质及飞行设备 | |
Li et al. | On Anti-Periodic Solutions for High-Order Cohen-Grossberg Neural Networks with Bounded Delays | |
CN105759602A (zh) | 基于硬件模拟运算的闭环控制系统及其应用 | |
CN117452831B (zh) | 一种四旋翼无人机控制方法、装置、系统及存储介质 | |
CN108107900B (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 |