发明内容
本发明的目的是为了解决上述问题,提出一种无人机纵向控制律平滑切换方法,通过该方法可确保无人机的舵面在纵向控制律切换时刻不发生瞬变,同时还可确定控制律切换后飞行控制系统的稳定裕度,从而提高飞行安全。
一种无人机纵向控制律平滑切换方法,针对制导回路控制律中有校正网络,姿控回路控制律中有校正网络的情况,包括以下几个步骤:
步骤一:将校正网络变换成新的表达形式,分离出隐含的积分器;
校正网络和校正网络分别属于制导回路和姿态回路,其中s为拉普拉斯算子,a、b、c、d、a1、b1、c1、d1为校正网络的系数,将上述两个校正网络变换成以下形式:
1)制导回路中的校正网络:
令
将Φ(s)视为单位反馈闭环传递函数,则其开环传递函数为
2)姿态回路中的校正网络:
令
将Φ1(s)视为单位反馈闭环传递函数,则其开环传递函数为
步骤二:定义积分器初值变量和中间变量;
1)定义积分器的初值变量
(1)制导回路PID控制器的积分器的初值变量为x0;
(2)制导回路中校正网络的积分器的初值变量为x1;
(3)姿控回路中校正网络的积分器的初值变量为x2;
2)定义中间变量
(1)制导回路中校正网络的输入信号为xin1;
(2)姿控回路中校正网络的输入信号为xin2;
(3)姿控回路的输入信号为;
(4)升降舵舵回路的输入信号为δzch;
步骤三:计算积分器初值;
控制律在ts时刻进行切换,已知由前一组控制律在ts时刻解算出的升降舵舵偏指令为δzcq,通过对新的控制律中的积分器赋初值x0、x1、x2,使得由新的控制律在ts+τ时刻解算出的升降舵舵偏指令δzcq等于δzch;
以下按照由内到外的顺序分别计算积分器初值x2、x1、x0:
1)姿控回路
在控制律切换时刻,使得以下等式成立:
则姿控回路中校正网络的输入信号xin2和积分器的初值x2分别为
由于
其中为俯仰角,ωz为俯仰角速度,为俯仰角反馈增益系数,为俯仰角速度反馈增益系数,
则姿控回路的输入信号为
2)制导回路
在控制律切换时刻,使得以下等式成立:
则制导回路校正网络的输入信号xin1和积分器的初值x1分别为
由于
xin1=kp·(hg-h)+x0+kd·(vyg-vy)
其中hg为高度给定指令,h为高度,vyg为升降速度给定指令,vy为升降速度,kp为高度反馈增益系数,kd为升降速度反馈增益系数,
则制导回路PID控制器中积分器的初值x0为
所以积分器初值x0、x1、x2的表达式分别为
其中
步骤四:将步骤三中计算的积分器初值x0、x1、x2代入需要切换的控制律中的积分器,最终实现两组纵向控制律之间的平稳切换。
本发明的优点在于:
(1)本发明提出的纵向控制律平滑切换方法,仅需一个运行周期即可实现控制律切换,无需通过仿真调整淡入淡出参数,控制律切换过程不依赖工程经验;
(2)本发明提出的纵向控制律平滑切换方法,在控制律切换时刻,不仅保证了由前后两组控制律解算的舵偏角相同,还保证了舵偏速率相同;
(3)本发明提出的纵向控制律平滑切换方法,控制律切换后,立即运算新的控制律,飞控系统的稳定裕度可根据新的控制律进行计算。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明的一种无人机纵向控制律平滑切换方法的构思是:
第一,明确控制回路中所有的积分器。针对高度跟踪控制,为了精确跟踪高度给定指令,通常采用PID控制,引入积分控制用于消除跟踪静差。有时因常规的PID控制无法满足系统的控制需求,通常会引入滞后超前校正网络串联入控制回路,来提高系统的稳定性能。其实,用于串联校正的校正网络中也隐含有积分器,可经过适当形式变换得到。可见,制导回路PID控制器和校正网络中均含有积分器。
第二,对积分器赋初值。在控制律切换时刻,已知当前的舵偏指令、新的控制律参数以及无人机的飞行状态,其中飞行状态包括高度、升降速度、俯仰角及俯仰角速率。可通过合适的计算方式,得到新的控制律中积分器的初值。
第三,将积分器的初值代入新的控制律中的积分器。由新的控制律在控制律切换后的第一个运行周期与由前一组控制律在控制律切换时刻分别解算得到的舵偏指令相等,从而实现控制律之间的平滑切换。
根据以上构思,按照制导回路控制律和姿控回路控制律中是否存在校正网络,对纵向控制律的平滑切换方式进行划分,存在以下四种情况:
1)制导回路控制律和姿控回路控制律中均没有校正网络;
2)制导回路控制律中有校正网络,姿控回路控制律中没有校正网络;
3)制导回路控制律中没有校正网络,姿控回路控制律中有校正网络;
4)制导回路控制律中有校正网络,姿控回路控制律中有校正网络。
针对前三种情况,在没有校正网络的控制回路中,可将比例系数1视为该回路的校正网络,可见前三种情况是第四种情况的特例,第四种情况更具有代表性、一般性和普遍性。以下针对第四种情况具体说明本发明纵向控制律平滑切换方法的实施步骤:
步骤一:将校正网络变换成新的表达形式,分离出隐含的积分器;
校正网络和校正网络分别属于制导回路和姿态回路,其中s为拉普拉斯算子,a、b、c、d、a1、b1、c1、d1为校正网络的系数,以下将上述两个校正网络变换成以下形式:
1)制导回路中的校正网络:
令
将Φ(s)视为单位反馈闭环传递函数,则其开环传递函数为
2)姿态回路中的校正网络:
令
将Φ1(s)视为单位反馈闭环传递函数,则其开环传递函数为
步骤二:定义积分器初值变量和中间变量;
由于整个控制回路中的信号节点较多,为避免计算积分器初值时造成混乱,在对控制律中的积分器定义初值变量的同时,还需定义一些中间变量。
1)定义积分器的初值变量
(1)制导回路PID控制器的积分器的初值变量为x0;
(2)制导回路中校正网络1的积分器的初值变量为x1;
(3)姿控回路中校正网络2的积分器的初值变量为x2;
2)定义中间变量
(1)制导回路中校正网络1的输入信号为xin1;
(2)姿控回路中校正网络2的输入信号为xin2;
(3)姿控回路的输入信号(制导回路的输出信号)为;
(4)升降舵舵回路的输入信号(姿控回路的输出信号)为δzch;
步骤三:计算积分器初值;
控制律在ts时刻进行切换,已知由前一组控制律在ts时刻解算出的升降舵舵偏指令为δzcq,通过对新的控制律中的积分器赋初值x0、x1、x2,使得由新的控制律在ts+τ时刻解算出的升降舵舵偏指令δzcq等于δzcq。
以下按照由内到外的顺序分别计算积分器初值x2、x1、x0:
1)姿控回路
在控制律切换时刻,使得以下等式成立:
则校正网络2的输入信号xin2和积分器的初值x2分别为
由于
其中为俯仰角,ωz为俯仰角速度,为俯仰角反馈增益系数,为俯仰角速度反馈增益系数,
则姿控回路的输入信号为
2)制导回路
在控制律切换时刻,使得以下等式成立:
则校正网络1的输入信号xin1和积分器的初值x1分别为
由于
xin1=kp·(hg-h)+x0+kd·(vyg-vy)
其中hg为高度给定指令,h为高度,vyg为升降速度给定指令,vy为升降速度,kp为高度反馈增益系数,kd为升降速度反馈增益系数,
则制导回路PID控制器中积分器的初值x0为
所以积分器初值x0、x1、x2的表达式分别为
其中
步骤四:将步骤三中计算的积分器初值x0、x1、x2代入需要切换的控制律中的积分器,最终实现两组纵向控制律之间的平稳切换。
本发明涉及一种无人机纵向控制律平滑切换方法。该方法通过对积分器赋初值,来保证控制律切换时前后两组控制律运算得出的舵偏角指令相同,实现控制律之间的平滑切换。
考虑某无人机沿五边定高飞行,无人机高度为100m,空速为38m/s,运行周期τ为40毫秒。在五边飞行过程中需放下起落架,并使襟翼偏至进场着陆所需要的度数20°。将起落架放下前定义为五边飞行前段,将起落架放下后定义为五边飞行后段。由于起落架收放前后无人机气动数据变化较大,为了保证整个飞行过程系统均具有较好的稳定性能,分别设计五边飞行前段和五边飞行后段的纵向控制律,图3给出了五边飞行前段的升降舵通道控制结构,图中kpq、kiq分别为高度跟踪误差比例增益系数和积分增益系数,kdq为升降速度跟踪误差增益系数,为俯仰角跟踪误差比例增益系数,为俯仰角速度反馈增益系数,δzcq为升降舵舵偏指令,图4给出了五边飞行后段的升降舵通道控制结构,图中kph、kih分别为高度跟踪误差比例增益系数和积分增益系数,kdh为升降速度跟踪误差增益系数,ah、bh、ch、dh为校正网络1的系数,为俯仰角跟踪误差比例增益系数,为俯仰角速度反馈增益系数,a1h、b1h、c1h、d1h为校正网络2的系数,δzch为升降舵舵回路的输入信号。利用经典控制理论分别设计五边飞行前段的控制律参数kpq、kiq、kdq、、及五边飞行后段的控制律参数kph、kih、kdh、、和校正网络参数ah、bh、ch、dh、a1h、b1h、c1h、d1h,使其满足控制性能指标。以下结合该实例给出本发明提出的一种无人机纵向控制律平滑切换方法的具体实施步骤:
步骤一:将校正网络变换成新的表达形式,分离出隐含的积分器;
校正网络和校正网络分别属于五边飞行后段的制导回路和姿态回路,其中s为拉普拉斯算子,ah、bh、ch、dh、a1h、b1h、c1h、d1h为校正网络的系数,以下将上述两个校正网络变换成以下形式:
1)制导回路中的校正网络:
令
将Φ(s)视为单位反馈闭环传递函数,则其开环传递函数为
制导回路中校正网络1拆分后的结构图如图5所示,图中含有积分器
2)姿态回路中的校正网络:
令
将Φ1h(s)视为单位反馈闭环传递函数,则其开环传递函数为
姿态回路中校正网络2拆分后的结构图如图6所示,图中含有积分器
步骤二:定义积分器初值变量和中间变量;
由于整个控制回路中的信号节点较多,为避免计算积分器初值时造成混乱,在对控制律中的积分器定义初值变量的同时,还需定义一些中间变量。
1)定义积分器的初值变量
(1)制导回路PID控制器的积分器的初值变量为x0;
(2)制导回路中校正网络1的积分器的初值变量为x1;
(3)姿控回路中校正网络2的积分器的初值变量为x2;
2)定义中间变量
(1)制导回路中校正网络1的输入信号为xin1;
(2)姿控回路中校正网络2的输入信号为xin2;
(3)姿控回路的输入信号(制导回路的输出信号)为;
(4)升降舵舵回路的输入信号(姿控回路的输出信号)为δzch;
定义积分器初值变量及中间变量后的控制结构图如图7所示。
步骤三:计算积分器初值;
控制律在ts=10s时刻进行切换,已知由五边飞行前段控制律在ts时刻解算出的升降舵舵偏指令为δzcq,通过对五边飞行后段控制律中的积分器赋初值x0、x1、x2,使得由五边飞行后段的控制律在ts+τ时刻解算出的升降舵舵偏指令δzcq等于δzch。
以下按照由内到外的顺序分别计算积分器初值x2、x1、x0:
1)姿控回路
在控制律切换时刻,使得以下等式成立:
则校正网络2的输入信号xin2和积分器的初值x2分别为
由于
其中为俯仰角,ωz为俯仰角速度,为俯仰角反馈增益系数,为俯仰角速度反馈增益系数,
则姿控回路的输入信号为
2)制导回路
在控制律切换时刻,使得以下等式成立:
则校正网络1的输入信号xin1和积分器的初值x1分别为
由于
xin1=kph·(hg-h)+x0+kdh·(vyg-vy),
其中hg为高度给定指令,h为高度,vyg为升降速度给定指令,vy为升降速度,kph为高度反馈增益系数,kdh为升降速度反馈增益系数,
则制导回路PID控制器中积分器的初值x0为
所以积分器初值x0、x1、x2的表达式分别为
其中
步骤四:将步骤三中计算的积分器初值x0、x1、x2代入五边飞行后段控制律中的积分器,最终实现两组纵向控制律之间的平稳切换。
图8~图13给出了分别在不采用平滑切换方法、采用单模态方法、采用双模态方法以及采用本发明提供的方法的情况下无人机相关参数的对比图,包括高度、俯仰角、俯仰角速度、迎角、过载、升降舵等参数。
由仿真结果可看出,控制律切换时刻,在采用本发明提出的方法时,无人机高度和俯仰角均具有较小的波动,其中高度差仅为采用其他舵面抑制措施时的四分之一,最小俯仰角仅为采用双模态抑制法时的三分之一。并且采用本发明提出的方法后无人机迎角、过载及升降舵均具有较小的波动,可很快达到新的平衡状态。本发明提出的控制律平滑切换方法效果明显,简单有效,不需要通过仿真调整淡入淡出参数,不依赖工程经验。