CN112034870B - 一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法 - Google Patents

一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法 Download PDF

Info

Publication number
CN112034870B
CN112034870B CN202010836796.XA CN202010836796A CN112034870B CN 112034870 B CN112034870 B CN 112034870B CN 202010836796 A CN202010836796 A CN 202010836796A CN 112034870 B CN112034870 B CN 112034870B
Authority
CN
China
Prior art keywords
gliding
guided projectile
representing
guided
control
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
CN202010836796.XA
Other languages
English (en)
Other versions
CN112034870A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN202010836796.XA priority Critical patent/CN112034870B/zh
Publication of CN112034870A publication Critical patent/CN112034870A/zh
Application granted granted Critical
Publication of CN112034870B publication Critical patent/CN112034870B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • G05D1/0816Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
    • G05D1/0825Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/107Simultaneous control of position or course in three dimensions specially adapted for missiles
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

本发明公开了一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法。该方法为:建立滑翔制导炮弹的姿态运动学模型;设计滑翔制导炮弹控制模型的辅助系统;测量与更新滑翔制导炮弹运动参数,并更新滑翔制导炮弹控制模型的辅助变量;选择状态跟踪误差变量;设计自动驾驶控制律;输入指令到滑翔制导炮弹鸭舵执行器,控制弹体的飞行;重复进行测量与更新运动参数和姿态控制,直至滑翔制导炮弹自动驾驶结束。本发明结构简单,鲁棒性强,便于实际应用。

Description

一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法
技术领域
本发明涉及滑翔制导炮弹技术领域,特别是一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法。
背景技术
现代导弹防御技术的快速发展,很大程度上限制了传统弹道导弹和飞弹的穿透能力,在这种背景下,滑翔制导炮弹被研制出来。作为一种新型杀伤性武器,滑翔制导炮弹结合了弹道导弹和飞弹的优点。但是由于滑翔制导炮弹的旋转特性,及其他例如舵机系统的响应滞后等因素,使得俯仰通道和偏航通道之间存在不可忽视的交叉耦合。更进一步的,滑翔制导炮弹的控制模型的建立通常需要基于一些假设条件,且在飞行中会受到强烈的内外扰动,因此滑翔制导炮弹的控制是一项极富挑战的任务。国内外学者就此提出了一些设计方法,例如MortezaShahravi等(MortezaShahravi,Mansour Kabganian,AriaAlasty.Adaptive robust attitude control of a flexible spacecraft[J].International Journal of Robust&Nonlinear Control,2010,16(6):287-302.)基于滑模控制提出了一种姿态控制方法,这种方法可以使扰动和不确定的影响降到最低。但是,这种方法仅仅考虑了外部扰动的影响,实际工程中很难获得精确的控制模型,且由于飞行条件的复杂性,气动参数等是时变的,这种内部扰动会给控制系统带来不可忽视的影响。Yangrong-jun等(Yang rong-jun,et al.Nonlinear Discrete-time Sliding ModeController Design for Extended Range Guided Projectiles[J].Journal of NanjingUniversity of Science and Technology,2012,36(1):137-141.)使用输入输出反馈线性化理论把控制系统解耦为两个独立的子系统,并基于滑模控制设计了相应的姿态控制律,为了压制子系统中不确定性的影响,这种方法需要假定扰动是存在上限的,且这种上限是已知的。应用中由于很难精确获知扰动的上限,需要对此进行估计,而不合理的假设会不可避免的降低系统的控制性能。Guo Z等(Guo Z,Zhou J,Guo J,et al.CouplingCharacterization-Based Robust Attitude Control Scheme for Hypersonic Vehicles[J].IEEE Transactions on Industrial Electronics,2017:1-1.)使用扰动观测器来估计总扰动,这种方法需要已知总扰动的二阶导数值,且所设计的观测器具有较多调试参数,不便于实际工程应用。
发明内容
本发明的目的在于提供一种结构简单、鲁棒性强的应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法,且便于实际工程应用。
实现本发明目的的技术解决方案为:一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法,包括以下步骤:
步骤1,建立滑翔制导炮弹的姿态运动学模型;
步骤2,设计滑翔制导炮弹控制模型的辅助系统;
步骤3,测量与更新滑翔制导炮弹运动参数,并更新滑翔制导炮弹控制模型的辅助变量;
步骤4,选择状态跟踪误差变量;
步骤5,设计自动驾驶控制律;
步骤6,输入指令到滑翔制导炮弹鸭舵执行器,控制弹体的飞行;
步骤7,重复步骤3-步骤6,直至滑翔制导炮弹自动驾驶结束。
进一步地,步骤1所述的建立滑翔制导炮弹的姿态运动学模型,具体如下:
滑翔制导炮弹的姿态运动学模型建立为:
Figure GDA0003627600380000021
其中α为攻角,β为侧滑角;θ和
Figure GDA0003627600380000022
分别为轨迹倾角和轨迹偏离角;δz和δy分别为偏航通道的等效舵偏角和俯仰通道的等效多偏角;ai和bi为动力学系数:
Figure GDA0003627600380000023
Figure GDA0003627600380000024
Figure GDA0003627600380000025
其中m为质量,v表示速度,D表示弹径,S表示特征面积,l表示特征长度,Q表示动压,
Figure GDA0003627600380000026
表示弹体转速,E和C分别为赤道阻尼系数和极阻尼系数,
Figure GDA0003627600380000027
Figure GDA0003627600380000028
为升力系统导数,C”μ表示马格努斯力矩导数,m'z表示静力矩导数,m'zz表示赤道阻尼力矩导数,m'σ表示鸭舵控制力矩导数,m”y表示力矩系数的导数;
舵机系统的输入输出模型建立为
Figure GDA0003627600380000029
其中,
Figure GDA0003627600380000031
δyc和δzc表示等效多偏角控制指令,γd表示总滞后角,Ts表示舵机系统时间,μs表示舵机系统阻尼系数,τ表示滞后时间;
假设
Figure GDA0003627600380000032
并使用
Figure GDA0003627600380000033
弹体姿态运动学模型进一步表示为
Figure GDA0003627600380000034
其中,k1=b3+a1,k2=b5+a2,k3=b1-b3a1+b5a2,k4=b2-b3a2-b5a1,k5=b4-b3a3,k6=b5a3,k7=a3为中间变量,无具体含义;w1和w2表示模型简化误差和外部扰动。
进一步地,步骤2所述的设计滑翔制导炮弹控制模型的辅助系统,具体如下:
Figure GDA0003627600380000037
Figure GDA0003627600380000038
其中zp和zq表示辅助变量,
Figure GDA0003627600380000035
表式参考系统状态变量,Apm,Bpm,Aqm,Bqm为选定的参考模型的状态矩阵和控制矩阵,βc和αc表示姿态参考指令。
进一步地,步骤3所述的测量与更新滑翔制导炮弹运动参数,并更新滑翔制导炮弹控制模型的辅助变量,具体如下:
通过包括陀螺仪、加速度计的在板测量系统更新运动参数:α,β,v,
Figure GDA0003627600380000036
并更新气动参数;使用数值积分方法更新辅助变量zp和zq,zp(0)=-Bp Txp(0)和zq(0)=-Bq Txq(0)。
进一步地,步骤4所述的选择状态跟踪误差变量,具体如下:
选择状态跟踪误差变量σp和σq为:
σp=Bp Txp+zp (19)
σq=Bq Txq+zq (20)
其中,Bp=Bq=[0,1]T
进一步地,步骤5所述的设计自动驾驶控制律,具体如下:
Figure GDA0003627600380000041
其中,矩阵Lp,Mp,Lq和Mq满足:
Ap-Apm=BpLp,Bpm=BpMp (22)
Aq-Aqm=BqLq,Bqm=BqMq (23)
Figure GDA0003627600380000042
其中,Kp>0,Kq>0,ξp>0和ξq>0为常数;
中间变量κ5和κ6满足:
Figure GDA0003627600380000043
本发明与现有技术相比,具有以下显著优点:(1)得益于引入的滑翔制导炮弹控制模型的辅助系统,无需获知系统精确模型,结构简单;(2)采用鲁棒控制技术,能够对抗强内外扰动,鲁棒性强;(3)无需已知扰动上限,便于实际应用。
附图说明
图1为本发明应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法的流程图。
图2为仿真过程中侧滑角指令跟踪轨迹图。
图3为仿真过程中攻角指令跟踪轨迹图。
图4为仿真过程中多偏角δy的轨迹图。
图5为仿真过程中多偏角δz的轨迹图。
具体实施方式
本发明应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法,包括以下步骤:
步骤1,建立滑翔制导炮弹的姿态运动学模型;
步骤2,设计滑翔制导炮弹控制模型的辅助系统;
步骤3,测量与更新滑翔制导炮弹运动参数,并更新滑翔制导炮弹控制模型的辅助变量;
步骤4,选择状态跟踪误差变量;
步骤5,设计自动驾驶控制律;
步骤6,输入指令到滑翔制导炮弹鸭舵执行器,控制弹体的飞行;
步骤7,重复步骤3-步骤6,直至滑翔制导炮弹自动驾驶结束。
进一步地,步骤1所述的建立滑翔制导炮弹的姿态运动学模型,具体如下:
滑翔制导炮弹的姿态运动学模型建立为:
Figure GDA0003627600380000051
其中α为攻角,β为侧滑角;θ和
Figure GDA0003627600380000052
分别为轨迹倾角和轨迹偏离角;δz和δy分别为偏航通道的等效舵偏角和俯仰通道的等效多偏角;ai和bi为动力学系数:
Figure GDA0003627600380000053
Figure GDA0003627600380000054
Figure GDA0003627600380000055
其中m为质量,v表示速度,D表示弹径,S表示特征面积,l表示特征长度,Q表示动压,
Figure GDA0003627600380000056
表示弹体转速,E和C分别为赤道阻尼系数和极阻尼系数,
Figure GDA0003627600380000057
Figure GDA0003627600380000058
为升力系统导数,C”μ表示马格努斯力矩导数,m'z表示静力矩导数,m'zz表示赤道阻尼力矩导数,m'σ表示鸭舵控制力矩导数,m”y表示力矩系数的导数;
舵机系统的输入输出模型建立为
Figure GDA0003627600380000059
其中,
Figure GDA00036276003800000510
δyc和δzc表示等效多偏角控制指令,γd表示总滞后角,Ts表示舵机系统时间,μs表示舵机系统阻尼系数,τ表示滞后时间;
假设
Figure GDA00036276003800000511
并使用
Figure GDA00036276003800000512
弹体姿态运动学模型进一步表示为
Figure GDA0003627600380000061
其中,k1=b3+a1,k2=b5+a2,k3=b1-b3a1+b5a2,k4=b2-b3a2-b5a1,k5=b4-b3a3,k6=b5a3,k7=a3为中间变量,无具体含义;w1和w2表示模型简化误差和外部扰动。
进一步地,步骤2所述的设计滑翔制导炮弹控制模型的辅助系统,具体如下:
Figure GDA0003627600380000062
Figure GDA0003627600380000063
其中zp和zq表示辅助变量,
Figure GDA0003627600380000064
表式参考系统状态变量,Apm,Bpm,Aqm,Bqm为选定的参考模型的状态矩阵和控制矩阵,βc和αc表示姿态参考指令。
进一步地,步骤3所述的测量与更新滑翔制导炮弹运动参数,并更新滑翔制导炮弹控制模型的辅助变量,具体如下:
通过包括陀螺仪、加速度计的在板测量系统更新运动参数:α,β,v,
Figure GDA0003627600380000065
并更新气动参数;使用数值积分方法更新辅助变量zp和zq,zp(0)=-Bp Txp(0)和zq(0)=-Bq Txq(0)。
进一步地,步骤4所述的选择状态跟踪误差变量,具体如下:
选择状态跟踪误差变量σp和σq为:
σp=Bp Txp+zp (31)
σq=Bq Txq+zq (32)
其中,Bp=Bq=[0,1]T
进一步地,步骤5所述的设计自动驾驶控制律,具体如下:
Figure GDA0003627600380000066
其中,矩阵Lp,Mp,Lq和Mq满足:
Ap-Apm=BpLp,Bpm=BpMp (34)
Aq-Aqm=BqLq,Bqm=BqMq (35)
Figure GDA0003627600380000071
其中,Kp>0,Kq>0,ξp>0和ξq>0为常数;
中间变量κ5和κ6满足:
Figure GDA0003627600380000072
下面结合附图及具体实施例对本发明作进一步详细说明。
实施例
结合图1,本实施例应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法,包括以下步骤:
步骤1,建立滑翔制导炮弹的姿态运动学模型,具体如下:
滑翔制导炮弹的姿态运动学模型可以建立为:
Figure GDA0003627600380000073
其中α为攻角,β为侧滑角;θ和
Figure GDA0003627600380000074
分别为轨迹倾角和轨迹偏离角;δz和δy分别为偏航通道的等效舵偏角和俯仰通道的等效多偏角;ai和bi为动力学系数:
Figure GDA0003627600380000075
Figure GDA0003627600380000076
Figure GDA0003627600380000077
其中m为质量,v表示速度,D表示弹径,S表示特征面积,l表示特征长度,Q表示动压,
Figure GDA0003627600380000078
表示弹体转速,E和C分别为赤道阻尼系数和极阻尼系数,
Figure GDA0003627600380000079
Figure GDA00036276003800000710
为升力系统导数,C”μ表示马格努斯力矩导数,m'z表示静力矩导数,m'zz表示赤道阻尼力矩导数,m'σ表示鸭舵控制力矩导数,m”y表示
舵机系统的的输入输出模型可以建立为
Figure GDA00036276003800000711
其中,
Figure GDA0003627600380000081
δyc和δzc表示等效多偏角控制指令,γd表示总滞后角,Ts表示舵机系统时间,μs表示舵机系统阻尼系数,τ表示滞后时间。
假设
Figure GDA0003627600380000082
并使用
Figure GDA0003627600380000083
弹体姿态运动学模型可进一步表示为
Figure GDA0003627600380000084
其中,k1=b3+a1,k2=b5+a2,k3=b1-b3a1+b5a2,k4=b2-b3a2-b5a1,k5=b4-b3a3,k6=b5a3,k7=a3为中间变量,无具体含义。w1和w2表示模型简化误差和外部扰动。
步骤2,设计滑翔制导炮弹控制模型的辅助系统,具体如下:
Figure GDA0003627600380000085
Figure GDA0003627600380000086
其中zp和zq表示辅助变量,
Figure GDA0003627600380000087
表式参考系统状态变量,Apm,Bpm,Aqm,Bqm为选定的参考模型的状态矩阵和控制矩阵,βc和αc表示姿态参考指令。
步骤3,更新运动参数与求解滑翔制导炮弹控制模型的辅助系统,具体如下:
通过陀螺仪、加速度计等在板测量系统更新运动参数:α,β,v,
Figure GDA0003627600380000088
等,并更新气动参数;使用数值积分方法更新辅助变量zp和zq,注意zp(0)=-Bp Txp(0)和zq(0)=-Bq Txq(0)。
步骤4,选择状态跟踪误差变量,具体如下:
选择状态跟踪误差变量为
σp=Bp Txp+zp (43)
σq=Bq Txq+zq (44)
其中,Bp=Bq=[0,1]T
步骤5,设计自动驾驶控制律,具体如下:
Figure GDA0003627600380000091
其中,矩阵Lp,Mp,Lq和Mq满足:
Ap-Apm=BpLp,Bpm=BpMp (46)
Aq-Aqm=BqLq,Bqm=BqMq (47)
Figure GDA0003627600380000092
其中,Kp>0,Kq>0,ξp>0和ξq>0为设计常数;
中间变量κ5和κ6满足:
Figure GDA0003627600380000093
步骤6,输入指令到滑翔制导炮弹的鸭舵执行器,控制弹体飞行;
步骤7,重复步骤3-步骤6,直至自动驾驶结束。
利用本实施例应用于滑翔制导炮弹的鲁棒自动驾驶方法,设计姿态控制仿真实验。仿真中所使用的某型滑翔制导炮弹的参数在表1中给出。控制参数设计为:Kp=4,ξp=0.11,Kq=1,ξq=0.1,
Figure GDA0003627600380000094
初始条件设置为:x(0)=[0.01 0 0.01 0]T
表1某型滑翔制导炮弹的相关参数
Figure GDA0003627600380000095
仿真中跟踪指令为:
Figure GDA0003627600380000101
图2为仿真过程中侧滑角指令跟踪轨迹图,从图2中可以看出,滑翔制导炮弹的侧滑角能够很好地跟踪阶跃指令和曲线型指令。
图3为仿真过程中攻角指令跟踪轨迹图,从图2中可以看出,滑翔制导炮弹的攻角能够很好地跟踪曲线型指令和阶跃型指令。
图4为仿真过程中多偏角δy的轨迹图,从图4中可以看出,仿真过程中舵偏角|δy|<20°,且较为平滑。
图5为仿真过程中多偏角δz的轨迹图,从图5中可以看出,仿真过程中舵偏角|δz|<10°,且较为平滑。
综上所述,本发明得益于引入的滑翔制导炮弹控制模型的辅助系统,无需获知系统精确模型,结构简单;采用鲁棒控制技术,能够对抗强内外扰动,鲁棒性强;无需已知扰动上限,便于实际应用。

Claims (1)

1.一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法,其特征在于,包括以下步骤:
步骤1,建立滑翔制导炮弹的姿态运动学模型;
步骤2,设计滑翔制导炮弹控制模型的辅助系统;
步骤3,测量与更新滑翔制导炮弹运动参数,并更新滑翔制导炮弹控制模型的辅助变量;
步骤4,选择状态跟踪误差变量;
步骤5,设计自动驾驶控制律;
步骤6,输入指令到滑翔制导炮弹鸭舵执行器,控制弹体的飞行;
步骤7,重复步骤3-步骤6,直至滑翔制导炮弹自动驾驶结束;
步骤1所述的建立滑翔制导炮弹的姿态运动学模型,具体如下:
滑翔制导炮弹的姿态运动学模型建立为:
Figure FDA0003627600370000011
其中α为攻角,β为侧滑角;θ和
Figure FDA0003627600370000012
分别为轨迹倾角和轨迹偏离角;δz和δy分别为偏航通道的等效舵偏角和俯仰通道的等效多偏角;ai和bi为动力学系数:
Figure FDA0003627600370000013
Figure FDA0003627600370000014
Figure FDA0003627600370000015
其中m为质量,v表示速度,D表示弹径,S表示特征面积,l表示特征长度,Q表示动压,
Figure FDA0003627600370000016
表示弹体转速,E和C分别为赤道阻尼系数和极阻尼系数,
Figure FDA0003627600370000017
Figure FDA0003627600370000018
为升力系统导数,C″μ表示马格努斯力矩导数,m′z表示静力矩导数,m′zz表示赤道阻尼力矩导数,m′σ表示鸭舵控制力矩导数,m″y表示力矩系数的导数;
舵机系统的输入输出模型建立为
Figure FDA0003627600370000019
其中,
Figure FDA0003627600370000021
δyc和δzc表示等效多偏角控制指令,γd表示总滞后角,Ts表示舵机系统时间,μs表示舵机系统阻尼系数,τ表示滞后时间;
假设
Figure FDA0003627600370000022
并使用
Figure FDA0003627600370000023
弹体姿态运动学模型进一步表示为
Figure FDA0003627600370000024
其中,k1=b3+a1,k2=b5+a2,k3=b1-b3a1+b5a2,k4=b2-b3a2-b5a1,k5=b4-b3a3,k6=b5a3,k7=a3为中间变量,无具体含义;w1和w2表示模型简化误差和外部扰动;
步骤2所述的设计滑翔制导炮弹控制模型的辅助系统,具体如下:
Figure FDA0003627600370000025
Figure FDA0003627600370000026
其中zp和zq表示辅助变量,
Figure FDA0003627600370000027
表式参考系统状态变量,Apm,Bpm,Aqm,Bqm为选定的参考模型的状态矩阵和控制矩阵,βc和αc表示姿态参考指令;
步骤3所述的测量与更新滑翔制导炮弹运动参数,并更新滑翔制导炮弹控制模型的辅助变量,具体如下:
通过包括陀螺仪、加速度计的在板测量系统更新运动参数:α,β,v,
Figure FDA0003627600370000028
并更新气动参数;使用数值积分方法更新辅助变量zp和zq,zp(0)=-Bp Txp(0)和zq(0)=-Bq Txq(0);
步骤4所述的选择状态跟踪误差变量,具体如下:
选择状态跟踪误差变量σp和σq为:
σp=Bp Txp+zp (7)
σq=Bq Txq+zq (8)
其中,Bp=Bq=[0,1]T
步骤5所述的设计自动驾驶控制律,具体如下:
Figure FDA0003627600370000031
其中,矩阵Lp,Mp,Lq和Mq满足:
Ap-Apm=BpLp,Bpm=BpMp (10)
Aq-Aqm=BqLq,Bqm=BqMq (11)
Figure FDA0003627600370000032
其中,Kp>0,Kq>0,ξp>0和ξq>0为常数;
中间变量κ5和κ6满足:
Figure FDA0003627600370000033
CN202010836796.XA 2020-08-19 2020-08-19 一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法 Active CN112034870B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010836796.XA CN112034870B (zh) 2020-08-19 2020-08-19 一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010836796.XA CN112034870B (zh) 2020-08-19 2020-08-19 一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法

Publications (2)

Publication Number Publication Date
CN112034870A CN112034870A (zh) 2020-12-04
CN112034870B true CN112034870B (zh) 2022-09-06

Family

ID=73578573

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010836796.XA Active CN112034870B (zh) 2020-08-19 2020-08-19 一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法

Country Status (1)

Country Link
CN (1) CN112034870B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5348249A (en) * 1993-01-11 1994-09-20 Hughes Missile Systems Company Retro reflection guidance and control apparatus and method
CN103353759A (zh) * 2013-06-25 2013-10-16 西安交通大学 一种基于cdm的导弹自动驾驶仪设计方法
CN106054612A (zh) * 2016-06-29 2016-10-26 河南科技大学 一种btt导弹飞行轨迹自动控制方法
CN108363305A (zh) * 2018-03-14 2018-08-03 中国人民解放军国防科技大学 基于主动干扰补偿的战术导弹鲁棒过载自驾仪设计方法
CN110162067A (zh) * 2018-12-25 2019-08-23 福建工程学院 一种无人机非奇异自适应非光滑姿态跟踪控制方法
CN110220416A (zh) * 2019-05-15 2019-09-10 南京理工大学 一种自适应快速弹道跟踪制导方法
CN110471456A (zh) * 2019-08-22 2019-11-19 中国人民解放军国防科技大学 高超声速飞行器俯冲段制导、姿控、变形一体化控制方法
CN111103890A (zh) * 2019-12-17 2020-05-05 西北工业大学 一种高精度强鲁棒的进场着陆引导控制方法
CN111324142A (zh) * 2020-01-07 2020-06-23 湖北航天技术研究院总体设计所 一种导弹驾驶仪扰动补偿控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109144084B (zh) * 2018-07-11 2019-08-06 哈尔滨工业大学 一种基于固定时间收敛观测器的垂直起降重复使用运载器姿态跟踪控制方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5348249A (en) * 1993-01-11 1994-09-20 Hughes Missile Systems Company Retro reflection guidance and control apparatus and method
CN103353759A (zh) * 2013-06-25 2013-10-16 西安交通大学 一种基于cdm的导弹自动驾驶仪设计方法
CN106054612A (zh) * 2016-06-29 2016-10-26 河南科技大学 一种btt导弹飞行轨迹自动控制方法
CN108363305A (zh) * 2018-03-14 2018-08-03 中国人民解放军国防科技大学 基于主动干扰补偿的战术导弹鲁棒过载自驾仪设计方法
CN110162067A (zh) * 2018-12-25 2019-08-23 福建工程学院 一种无人机非奇异自适应非光滑姿态跟踪控制方法
CN110220416A (zh) * 2019-05-15 2019-09-10 南京理工大学 一种自适应快速弹道跟踪制导方法
CN110471456A (zh) * 2019-08-22 2019-11-19 中国人民解放军国防科技大学 高超声速飞行器俯冲段制导、姿控、变形一体化控制方法
CN111103890A (zh) * 2019-12-17 2020-05-05 西北工业大学 一种高精度强鲁棒的进场着陆引导控制方法
CN111324142A (zh) * 2020-01-07 2020-06-23 湖北航天技术研究院总体设计所 一种导弹驾驶仪扰动补偿控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
导弹末制导过程自抗扰连续 2 阶滑模导引律设计;曲悦;《兵工学报》;20200430;第14-24页 *

Also Published As

Publication number Publication date
CN112034870A (zh) 2020-12-04

Similar Documents

Publication Publication Date Title
CN106681348B (zh) 考虑全捷联导引头视场约束的制导控制一体化设计方法
CN108803649B (zh) 一种垂直起降重复使用运载器自抗扰滑模控制方法
CN108363305B (zh) 基于主动干扰补偿的战术导弹鲁棒过载自驾仪设计方法
CN110316358A (zh) 基于动态逆的战斗机大迎角控制方法
CN111399531A (zh) 一种高超声速飞行器滑翔段制导与姿态控制一体化设计方法
CN109703768B (zh) 一种基于姿态/轨迹复合控制的软式空中加油对接方法
CN111367182A (zh) 考虑输入受限的高超声速飞行器抗干扰反步控制方法
CN107957686B (zh) 基于预见控制的无人直升机自动着舰控制系统
CN111610794A (zh) 基于滑模干扰观测器的战斗机大迎角动态逆控制方法
CN111045440B (zh) 一种高超声速飞行器俯冲段快速滚转控制方法
CN104597911A (zh) 空中加油受油机自适应最优对接轨迹跟踪飞行控制方法
CN114815888B (zh) 一种仿射形式的制导控制一体化控制方法
CN117471952A (zh) 一种飞行器反步超螺旋滑模制导控制一体化方法
CN109703769A (zh) 一种基于预瞄策略的空中加油对接控制方法
CN111240204A (zh) 一种基于模型参考滑模变结构控制的巡飞弹控制方法
CN113485395A (zh) 一种误差约束下固定时间相对姿轨跟踪控制方法
CN112034870B (zh) 一种应用于滑翔制导炮弹的鲁棒姿态自动驾驶方法
CN113485396A (zh) 一种航天器交会对接最后逼近段相对轨道与姿态跟踪控制方法
CN112327629A (zh) 基于动态补偿的小型无人直升机自适应容错控制方法
CN116795126A (zh) 一种输入饱和与输出受限的变形飞行器控制方法
Yao et al. Trajectory tracking controller based on PID-NLADRC
CN116360258A (zh) 基于固定时间收敛的高超声速变形飞行器抗干扰控制方法
Elbaioumy et al. Modelling and Simulation of Surface to Surface Missile General Platform
CN114265420B (zh) 适于高动态、控制慢响应的制导控制一体化设计方法
CN116736692A (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