CN111913491A - 一种基于视线角非线性抗饱和与不确定性补偿的制导方法 - Google Patents
一种基于视线角非线性抗饱和与不确定性补偿的制导方法 Download PDFInfo
- Publication number
- CN111913491A CN111913491A CN202011003323.8A CN202011003323A CN111913491A CN 111913491 A CN111913491 A CN 111913491A CN 202011003323 A CN202011003323 A CN 202011003323A CN 111913491 A CN111913491 A CN 111913491A
- Authority
- CN
- China
- Prior art keywords
- signal
- line
- sight angle
- sight
- saturation
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000009466 transformation Effects 0.000 claims abstract description 10
- 230000006641 stabilisation Effects 0.000 claims description 8
- 238000011105 stabilization Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000003786 synthesis reaction Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 4
- 238000013016 damping Methods 0.000 abstract description 2
- 230000001965 increasing effect Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
- G05D1/0816—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
- G05D1/0825—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Combustion & Propulsion (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种基于视线角非线性抗饱和与不确定性补偿的制导方法,属于飞行器飞行制导技术领域,其特点在于仅通过测量飞行器的视线角信息,通过构造三阶惯性非线性滤波器获得近似差分信号为飞行器制导提供所需的阻尼。然后通过非线性抗饱和变换得到相应的非线性信号。最后,通过构造制导干扰与不确定补偿器,得到制导系统的不确定性干扰的补偿信号,从而为提高整个制导设计的精度与适应度打下了良好的基础。最终综合了上述补偿信号、视线角信号以及抗饱和非线性变换信号与积分信号形成了制导综合信号。该基于干扰估计补偿的新型制导方法具有测量信息简单、制导精度高,尤其是鲁棒性好与适应度广泛的优点。
Description
技术领域
本发明属于飞行器飞行制导领域,尤其是涉及一种飞行器视线角非线性抗饱和与不确定性补偿的制导方法。
背景技术
传统的飞行器导引法一般是采用导引头测量视线角速率,但该角速度的准确测量并不容易。而且视线角速率到了目标末端,会出现增大,而使得视线角速率出现饱和情况。如比例导引,采用视线角速率形成最终的导引信号时在末端就会出现信号过早进入饱和而导致脱靶量增加的问题。
同时,由于高速运动的飞行器动力学模型非常复杂,因此在设计中采用的风洞数据模型尽管以及尽量精确,但仍然不可避免地存在模型误差、干扰等不确定性。为了综合考虑该模型误差等不确定性对导引精度的影响,对系统不确定性进行一定程度的预估补偿,是一种比传统导引方法更为高级的设计思想,其也是更深程度的一种反馈,从而能大大增加导引精度,尤其是增大初始条件变换情况下的导引律适应度,也就是增大了导引律的适应度。
基于以上背景原因,本发明提出一种抗饱和的非线性变换与不确定性补偿的制导方法,解决了视线角与视线角差分信号过早进入饱和问题,同时不确定性预估补偿的方法也大大提高了导引精度与鲁棒性。
需要说明的是,在上述背景技术部分发明的信息仅用于加强对本发明的背景的理解,因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。
发明内容
本发明的目的在于提供一种基于视线角非线性抗饱和与不确定性补偿的制导方法,进而至少在一定程度上克服传统导引方法的鲁棒性不强与适应性不强的问题。
本发明提供了一种基于视线角非线性抗饱和与不确定性补偿的制导方法,包括以下步骤:
步骤S10:采用导引头测量飞行器与目标的视线角信号,或者由地面站提供飞行器与目标的视线角信息;
步骤S20:根据所述的视线角信号,构造三阶惯性非线性滤波器,求解视线角信号的近似差分信号;
步骤S30:根据所述的视线角信号与视线角近似差分信号,分别进行非线性变换,得到视线角信号的非线性抗饱和信号与视线角近似差分信号的非线性抗饱和信号;
步骤S40:根据所述的视线角信号、视线角近似差分信号、视线角信号的非线性抗饱和信号、以及视线角近似差分信号的非线性抗饱和信号构造制导干扰与不确定补偿器,得到不确定补偿信号;
步骤S50:根据所述的不确定补偿信号、视线角信号、视线角近似差分信号、视线角信号的非线性抗饱和信号、以及视线角近似差分信号的非线性抗饱和信号进行线性叠加,组成最终的导引综合信号,输送给飞行器姿态稳定回路,实现对目标的精确导引。
在本发明的一种示例实施例中,根据所述的视线角信号,构造三阶惯性非线性滤波器,求解视线角信号的近似差分信号包括:
z0(1)=0、z2(1)=0、z3(1)=0;
z4=z0-qr;
其中z0、z2、z3为三阶惯性非线性滤波器的状态变量,其初始值均设置为0。qr为视线角信号,也是三阶惯性非线性滤波器的输入信号,z4为输入状态偏差信号,z1、zd2、zd3为三阶惯性非线性滤波器的中间变量,k1、k2、k3、k4、k5、k6、k7、ε1、ε2、ε3为常值参数,其详细设计见后文案例实施。z0(n)、z2(n)、z3(n)为状态z0、z2、z3的第n个数据。T为数据间的时间间隔,其详细设计见后文案例实施。z3即为最终所求的视线角信号的近似差分信号。
在本发明的一种示例实施例中,根据所述的视线角信号与视线角近似差分信号,分别进行非线性变换,得到相应的非线性抗饱和信号包括:
其中qr为视线角信号,z3为视线角近似差分信号,ε4、ε5为常值参数。qrb为视线角信号的非线性抗饱和信号,z3b为视线角近似差分信号的非线性抗饱和信号。
在本发明的一种示例实施例中,根据所述的视线角信号、视线角近似差分信号、视线角信号的非线性抗饱和信号、以及视线角近似差分信号的非线性抗饱和信号构造制导干扰与不确定补偿器,得到不确定补偿信号包括:
y=x+qr;
xa=-ka1qr-ka2z3-ka3qrb-ka4z3b-y-ua;
x(n+1)=x(n)+xaT;
其中x为干扰与不确定补偿器的状态信号,设置其初始状态为0,即x(0)=0。x(n)、x(n+1)分别为状态x的第n个数据与第n+1个数据。ka1、ka2、ka3、ka4为常值参数,T为数据间的时间间隔。qr为视线角信号、z3为视线角近似差分信号、qrb视线角信号的非线性抗饱和信号、z3b为视线角近似差分信号的非线性抗饱和信号,ua为导引综合信号,其详细设计见下文。而y即为本步最终所求的不确定补偿信号,其由状态x与视线角信号叠加而成。
在本发明的一种示例实施例中,根据所述的不确定补偿信号组成最终的导引综合信号包括:
ua=kb1qr+kb2z3+kb3qrb+kb4z3b+kb5∫qrdt-y;
其中dt表示对时间信号积分,kb1、kb2、kb3、kb4、kb5为常值参数,qr为视线角信号、z3为视线角近似差分信号、qrb视线角信号的非线性抗饱和信号、z3b为视线角近似差分信号的非线性抗饱和信号,y为不确定补偿信号。ua即为本步最终所求导引综合信号。
最终,将上述导引综合信号提供给飞行器姿态控制系统,即可实现飞行器对目标的精确导引。
本发明提供一种基于视线角非线性抗饱和与不确定性补偿的制导方法,该方法的特点有三,其一是提出了一种三阶惯性非线性滤波器的方法来求解视线角信号的近似差分,从而为导引提供阻尼信号,减少了测量元器件对速率信号的测量,从而节省了经济成本。其二在于提出了一类非线性抗饱和方法,减慢了信号进入饱和的时间,从而提高了导引精度。其三在于,采用了模型不确定干扰预估补偿的思想,增大了反馈的深度,从而大大提高了导引方法对不同导引初始条件下的适应度与鲁棒性。因此,本发明具有很高的工程应用价值。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本发明。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的一种基于视线角非线性抗饱和与不确定性补偿的制导方法的设计实施流程图;
图2是本发明实施例所提供方法的飞行器与目标视线角信号曲线(单位:度);
图3是本发明实施例所提供方法的视线角信号的近似差分信号信号曲线(无单位);
图4是本发明实施例所提供方法的视线角信号的非线性抗饱和信号曲线(单位:度);
图5是本发明实施例所提供方法的视线角近似差分信号的非线性抗饱和信号曲线(无单位);
图6是本发明实施例所提供方法的不确定补偿信号(无单位);
图7是本发明实施例所提供方法的导引综合信号(无单位);
图8本发明实施例所提供方法的飞行器偏航舵偏角曲线(单位:度);
图9是本发明实施例所提供方法的飞行器侧滑角曲线(单位:度);
图10本发明实施例所提供方法的飞行器偏航角曲线(单位:度);
图11本发明实施例所提供方法的飞行器与目标运动轨迹(单位:米);
图12本发明实施例所提供方法的器目距离曲线(单位:米)。
具体实施方式
现在将参考附图更全面地描述示例实施方式。然而,示例实施方式能够以多种形式实施,且不应被理解为限于在此阐述的范例;相反,提供这些实施方式使得本发明将更加全面和完整,并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。在下面的描述中,提供许多具体细节从而给出对本发明的实施方式的充分理解。然而,本领域技术人员将意识到,可以实践本发明的技术方案而省略所述特定细节中的一个或更多,或者可以采用其它的方法、组元、装置、步骤等。在其它情况下,不详细示出或描述公知技术方案以避免喧宾夺主而使得本发明的各方面变得模糊。
本发明一种基于视线角非线性抗饱和与不确定性补偿的制导方法,仅需要采用导引头或其它设备测量飞行器与目标的视线角,然后构造三阶惯性非线性滤波器,求解其近似差分信号。再进行非线性抗饱和变换得到非线性抗饱和信号。最后,提出了一种采用视线角与近似差分信号构造制导系统不确定性观测预估的方法,对综合制导信号进行补偿,从而提高制导律的精度与适应度以及鲁棒性。
以下,将结合附图对本发明实例实施例中提及的一种基于视线角非线性抗饱和与不确定性补偿的制导方法进行解释以及说明。参考图1所示,一种基于视线角非线性抗饱和与不确定性补偿的制导方法可以包括以下步骤:
步骤S10:采用导引头测量飞行器与目标的视线角信号,或者由地面站提供飞行器与目标的视线角信息。
具体的,如果飞行器上能够安装导引头等测量设备,则采用导引头设备测量飞行器与目标的视线角,记作qr。如果飞行器上空间不允许,则可以由地面设备进行目标与飞行器间相对运动信息测量,通过转换计算得到视线角信息qr。
步骤S20:根据所述的视线角信号,构造三阶惯性非线性滤波器,求解视线角信号的近似差分信号。
具体的,首先定义z0、z2、z3为三阶惯性非线性滤波器的状态变量,其初始值均设置为0,即z0(1)=0、z2(1)=0、z3(1)=0。
其次,设置视线角为三阶惯性非线性滤波器的输入信号,求解状态变量z0与输入信号的偏差。定义该偏差为输入状态偏差信号,记作z4,其计算方式如下:
z4=z0-qr;
再次,设计三阶惯性非线性滤波器的状态调节规律如下所示:
其中z1、zd2、zd3为三阶惯性非线性滤波器的中间变量,k1、k2、k3、k4、k5、k6、k7、ε1、ε2、ε3为常值参数,其详细设计见后文案例实施。z0(n)、z2(n)、z3(n)为状态z0、z2、z3的第n个数据。T为数据间的时间间隔,其详细设计见后文案例实施。z3即为最终所求的视线角信号的近似差分信号。
步骤S30:根据所述的视线角信号与视线角近似差分信号,分别进行非线性变换,得到视线角信号的非线性抗饱和信号与视线角近似差分信号的非线性抗饱和信号。
具体的,首先针对所述的视线角信号与视线角近似差分信号,进行非线性抗饱和变换,得到视线角信号的非线性抗饱和信号,记作qrb,其变换方式如下:
其中qr为视线角信号,z3为视线角近似差分信号,ε4为常值参数,其详细设计见后文案例实施。
其次,对所述的所述的视线角信号与视线角近似差分信号,进行非线性抗饱和变换,得到视线角近似差分信号的非线性抗饱和信号,记作z3b,其计算方式如下:
其中qr为视线角信号,z3为视线角近似差分信号,ε5为常值参数,其详细设计见后文案例实施。
步骤S40:根据所述的视线角信号、视线角近似差分信号、视线角信号的非线性抗饱和信号、以及视线角近似差分信号的非线性抗饱和信号构造制导干扰与不确定补偿器,得到不确定补偿信号。
具体的,首先设置干扰与不确定补偿器的状态信号为x,设置其初始状态为0,即x(0)=0。
其次,初步计算不确定补偿信号,将其定义为不确定性补偿器状态x与视线角信号qr的线性叠加,其中不确定补偿信号记作y,其计算方式为:
y=x+qr;
最后,进行制导干扰与不确定补偿器的状态调节律设计如下:
xa=-ka1qr-ka2z3-ka3qrb-ka4z3b-y-ua;
x(n+1)=x(n)+xaT;
其中x(n)、x(n+1)分别为状态x的第n个数据与第n+1个数据。ka1、ka2、ka3、ka4为常值参数,T为数据间的时间间隔,其详细设计见后文案例实施。qr为视线角信号、z3为视线角近似差分信号、qrb视线角信号的非线性抗饱和信号、z3b为视线角近似差分信号的非线性抗饱和信号,ua为导引综合信号,其详细设计见后文。而y即为本步最终所求的不确定补偿信号,其由状态x与视线角信号叠加而成。
步骤S50:根据所述的不确定补偿信号、视线角信号、视线角近似差分信号、视线角信号的非线性抗饱和信号、以及视线角近似差分信号的非线性抗饱和信号进行线性叠加,组成最终的导引综合信号,输送给飞行器姿态稳定回路,实现对目标的精确导引。
具体的,首先,将导引综合信号记作ua,其计算如下:
ua=kb1qr+kb2z3+kb3qrb+kb4z3b+kb5∫qrdt-y;
其中dt表示对时间信号积分,kb1、kb2、kb3、kb4、kb5为常值参数,其详细设计见后文案例实施。qr为视线角信号、z3为视线角近似差分信号、qrb视线角信号的非线性抗饱和信号、z3b为视线角近似差分信号的非线性抗饱和信号,y为不确定补偿信号。
其次,将导引综合信号输送给飞行器姿态稳定系统,使得飞行器偏航角跟踪该综合导引信号,即可实现飞行器对目标的精确导引。在此不再对姿态稳定系统进行详细说明,主要是由于姿态稳定系统的设计方法比较成熟,而且非本发明重点保护的内容,因此仅在案例实施中选取一类PID算法实现姿态稳定控制,来实现本发明例的导引方法设计,从而验证本发明方法的正确性。
案例实施与计算机仿真模拟结果分析
为验证本发明例所提供方法的有效性,进行以下的案例仿真。首先设置飞行器初始位置为三维坐标为(0,0,0),目标的初始位置三维坐标为(5800,1,-750)。其中1为目标的高度,-750为目标的初始偏差,5800为飞行器距离目标的初始距离。而目标速度设定为在水平面内运动,速度大小为26m/s,方向为与x轴方向成42度。
在步骤S10中,采用导引头测量飞行器与目标的视线角信号如图2所示。
在步骤S20中,选取k1=8、k2=10、k3=5、k4=2、k5=4、k6=3、k7=30、ε1=2、ε2=4、ε3=1、T=0.001,根据所述的视线角信号,构造三阶惯性非线性滤波器,求解视线角信号的近似差分信号如图3所示。
在步骤S30中,选取ε4=3;ε5=3,根据所述的视线角信号与视线角近似差分信号,分别进行非线性变换,得到视线角信号的非线性抗饱和信号与视线角近似差分信号的非线性抗饱和信号如图4与图5所示。
在步骤S40中,选取ka1=-3、ka2=-0.1、ka3=-2、ka4=-0.1,得到不确定补偿信号如图6所示。
在步骤S50中,选取kb1=2、kb2=1、kb3=0.5、kb4=2、kb5=0.1,得到最终的导引综合信号如图7所示,输送给飞行器姿态稳定回路,实现对目标的精确导引。最终得到的偏航舵偏角如图8所示,飞行器侧滑角如图9所示,飞行器的偏航角角如图10所示,飞行器与目标的相对运动曲线如图11所示,飞行器与目标的相对距离曲线如图12所示。
由图8可以看出,飞行器的舵偏角小于8度,从而没有超出舵偏角的限制,而由图10可以看出,飞行器的偏航角也变化平稳,没有剧烈颤振。而由图6与图7可以看出,系统不确定性估计值在最终的导引律中占有很大的比例,从而可以看出不确定性的估计在导引律对干扰、初始位置变化等起到了很大的稳定调节作用。图12的脱靶量曲线以及后续的目标位置改变的仿真实验也验证了这一点。这是其它方法所难以保证的。
由图11与图12可以看出,飞行器能够精准命中目标,脱靶量小于0.6米。而且跟换目标初始位置如增大5000米,减少2000米等,以及改变目标速度、速度方向后,脱靶量均小于1米,可见本发明所提供方法具有很高的精度,而且具有很好的适应性,从而具有很高的工程应用价值。
本领域技术人员在考虑说明书及实践这里发明的发明后,将容易想到本发明的其他实施例。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未发明的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由权利要求指出。
Claims (5)
1.一种基于视线角非线性抗饱和与不确定性补偿的制导方法,其特征在于,包括以下步骤:
步骤S10:采用导引头测量飞行器与目标的视线角信号,或者由地面站提供飞行器与目标的视线角信息;
步骤S20:根据所述的视线角信号,构造三阶惯性非线性滤波器,求解视线角信号的近似差分信号;
步骤S30:根据所述的视线角信号与视线角近似差分信号,分别进行非线性变换,得到视线角信号的非线性抗饱和信号与视线角近似差分信号的非线性抗饱和信号;
步骤S40:根据所述的视线角信号、视线角近似差分信号、视线角信号的非线性抗饱和信号、以及视线角近似差分信号的非线性抗饱和信号构造制导干扰与不确定补偿器,得到不确定补偿信号;
步骤S50:根据所述的不确定补偿信号、视线角信号、视线角近似差分信号、视线角信号的非线性抗饱和信号、以及视线角近似差分信号的非线性抗饱和信号进行线性叠加,组成最终的导引综合信号,输送给飞行器姿态稳定回路,实现对目标的精确导引。
2.根据权利要求1所述的所述的一种基于视线角非线性抗饱和与不确定性补偿的制导方法,其特征在于,根据所述的视线角信号,构造三阶惯性非线性滤波器,求解视线角信号的近似差分信号包括:
z0(1)=0、z2(1)=0、z3(1)=0;
z4=z0-qr;
z0(n+1)=z0(n)+z1T
z2(n+1)=z2(n)+zd2T;
z3(n+1)=z3(n)+zd3T
其中z0、z2、z3为三阶惯性非线性滤波器的状态变量,其初始值均设置为0。qr为视线角信号,也是三阶惯性非线性滤波器的输入信号,z4为输入状态偏差信号,z1、zd2、zd3为三阶惯性非线性滤波器的中间变量,k1、k2、k3、k4、k5、k6、k7、ε1、ε2、ε3为常值参数,其详细设计见后文案例实施。z0(n)、z2(n)、z3(n)为状态z0、z2、z3的第n个数据。T为数据间的时间间隔,其详细设计见后文案例实施。z3即为最终所求的视线角信号的近似差分信号。
4.根据权利要求1所述的所述的一种基于视线角非线性抗饱和与不确定性补偿的制导方法,其特征在于,根据所述的视线角信号、视线角近似差分信号、视线角信号的非线性抗饱和信号、以及视线角近似差分信号的非线性抗饱和信号构造制导干扰与不确定补偿器,得到不确定补偿信号包括:
y=x+qr;
xa=-ka1qr-ka2z3-ka3qrb-ka4z3b-y-ua;
x(n+1)=x(n)+xaT;
其中x为干扰与不确定补偿器的状态信号,设置其初始状态为0,即x(0)=0。x(n)、x(n+1)分别为状态x的第n个数据与第n+1个数据。ka1、ka2、ka3、ka4为常值参数,T为数据间的时间间隔。qr为视线角信号、z3为视线角近似差分信号、qrb视线角信号的非线性抗饱和信号、z3b为视线角近似差分信号的非线性抗饱和信号,ua为导引综合信号,其详细设计见下文。而y即为本步最终所求的不确定补偿信号,其由状态x与视线角信号叠加而成。
5.根据权利要求4所述的所述的一种基于视线角非线性抗饱和与不确定性补偿的制导方法,其特征在于,根据所述的不确定补偿信号组成最终的导引综合信号包括:
ua=kb1qr+kb2z3+kb3qrb+kb4z3b+kb5∫qrdt-y;
其中dt表示对时间信号积分,kb1、kb2、kb3、kb4、kb5为常值参数,qr为视线角信号、z3为视线角近似差分信号、qrb视线角信号的非线性抗饱和信号、z3b为视线角近似差分信号的非线性抗饱和信号,y为不确定补偿信号。ua即为本步最终所求导引综合信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011003323.8A CN111913491B (zh) | 2020-09-22 | 2020-09-22 | 一种基于视线角非线性抗饱和与不确定性补偿的制导方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011003323.8A CN111913491B (zh) | 2020-09-22 | 2020-09-22 | 一种基于视线角非线性抗饱和与不确定性补偿的制导方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111913491A true CN111913491A (zh) | 2020-11-10 |
CN111913491B CN111913491B (zh) | 2022-04-01 |
Family
ID=73265381
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011003323.8A Active CN111913491B (zh) | 2020-09-22 | 2020-09-22 | 一种基于视线角非线性抗饱和与不确定性补偿的制导方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111913491B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114153144A (zh) * | 2021-11-01 | 2022-03-08 | 天津大学 | 一种输入受限和输入扰动的弹性高超声速飞行器控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IT8249227A0 (it) * | 1981-10-08 | 1982-10-07 | Saab Scania Ab | Procedimento ed apparecchio per il pilotaggio di un corpo aerodinamico dotato di radioguida |
CN107390523A (zh) * | 2017-07-13 | 2017-11-24 | 西北工业大学 | 空间绳系复合体系统的自适应神经网络动态面控制器 |
CN111290415A (zh) * | 2019-12-04 | 2020-06-16 | 中国人民解放军海军航空大学 | 一种基于近似差分的飞行器综合前置导引方法 |
CN111399529A (zh) * | 2020-04-02 | 2020-07-10 | 上海交通大学 | 一种基于非线性滑模与前置的飞行器复合导引方法 |
-
2020
- 2020-09-22 CN CN202011003323.8A patent/CN111913491B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IT8249227A0 (it) * | 1981-10-08 | 1982-10-07 | Saab Scania Ab | Procedimento ed apparecchio per il pilotaggio di un corpo aerodinamico dotato di radioguida |
IT1203644B (it) * | 1981-10-08 | 1989-02-15 | Saab Scania Ab | Procedimento ed apparecchio per il pilotaggio di un corpo aerodinamico dotato di radioguida |
CN107390523A (zh) * | 2017-07-13 | 2017-11-24 | 西北工业大学 | 空间绳系复合体系统的自适应神经网络动态面控制器 |
CN111290415A (zh) * | 2019-12-04 | 2020-06-16 | 中国人民解放军海军航空大学 | 一种基于近似差分的飞行器综合前置导引方法 |
CN111399529A (zh) * | 2020-04-02 | 2020-07-10 | 上海交通大学 | 一种基于非线性滑模与前置的飞行器复合导引方法 |
Non-Patent Citations (6)
Title |
---|
WANG KAI,SUN CHUNZHEN,JIANG YI: "Research on Adaptive Guidance Technology of UAV Ship Landing System Based on Net Recovery", 《PROCEDIA ENGINEERING》 * |
WANG ZHUANG: "Design of Agent Training Environment for Aircraft Landing Guidance Based on Deep Reinforcement Learning", 《2018 11TH INTERNATIONAL SYMPOSIUM ON COMPUTATIONAL INTELLIGENCE AND DESIGN (ISCID), VOL 2》 * |
史小平等: "Hamilton系统的鲁棒自适应控制及在航天器中的应用", 《四川大学学报(工程科学版)》 * |
孙经广,孟庆鹏,李传明: "带有落角约束的反舰导弹自适应饱和制导律设计", 《指挥控制与仿真》 * |
林海兵: "可重复使用运载器约束制导与鲁棒容错控制", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
雷军委,肖支才,白祥军: "基于柔化与抗饱和函数的高超声速飞行器自适应滑模控制研究", 《战术导弹技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114153144A (zh) * | 2021-11-01 | 2022-03-08 | 天津大学 | 一种输入受限和输入扰动的弹性高超声速飞行器控制方法 |
CN114153144B (zh) * | 2021-11-01 | 2023-08-22 | 天津大学 | 一种输入受限和输入扰动的弹性高超声速飞行器控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111913491B (zh) | 2022-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111399529B (zh) | 一种基于非线性滑模与前置的飞行器复合导引方法 | |
CN103512426B (zh) | 一种次优的带末角约束制导方法 | |
CN108153323B (zh) | 一种高空无人飞行器高精度再入制导方法 | |
CN111538236B (zh) | 基于分数阶近似微分实现阻尼的飞行器纵向过载控制方法 | |
CN105157705A (zh) | 一种半捷联雷达导引头视线角速度提取方法 | |
CN111309040B (zh) | 一种采用简化分数阶微分的飞行器纵向俯仰角控制方法 | |
CN102425980A (zh) | 利用加速度计实现过载驾驶仪的控制方法 | |
CN105486308A (zh) | 估计弹目视线角速率的快收敛Kalman滤波器的设计方法 | |
CN111399530A (zh) | 一种基于反解传函的小型飞行器攻角滑模跟踪方法 | |
CN111367307A (zh) | 一种用校正网络代替角加速度计的飞行器侧向过载跟踪方法 | |
CN111290415A (zh) | 一种基于近似差分的飞行器综合前置导引方法 | |
CN111913491B (zh) | 一种基于视线角非线性抗饱和与不确定性补偿的制导方法 | |
CN110703793B (zh) | 采用姿态角测量的飞行器积分比例导引攻击机动目标的方法 | |
CN111102886B (zh) | 一种小微飞行器无陀螺精确制导方法 | |
CN111708382B (zh) | 一种基于非线性比例积分的飞行器导引方法 | |
CN111324136B (zh) | 一种微型飞行器位置与距离复合作用导引方法 | |
CN111650947B (zh) | 一种平流层飞艇高度非线性控制方法 | |
CN108762528B (zh) | 适用于空中飞鼠的姿态解算方法 | |
CN111734536B (zh) | 一种采用自适应分数阶微分的发动机稳定供油方法 | |
CN111290418B (zh) | 一种小微飞行器无稳定回路精确微分导引方法 | |
CN111309049B (zh) | 一种微型飞行器的虚拟目标五米小幅高精度导引方法 | |
CN112034889B (zh) | 一种采用非线性超前网络的飞行器过载控制方法 | |
CN112027117B (zh) | 一种基于姿态测量的飞行器侧滑滚转复合转弯控制方法 | |
CN111766776B (zh) | 一种采用非线性比例与积分型补偿预估器的前置导引方法 | |
CN112097765B (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 |