定义飞机模型

以下线性时不变模型是由一架飞机在3000英尺高度、速度为0.6马赫([1])时纵向动力学的线性化得到的。开环模型有以下状态空间矩阵:

A = [-0.0151 -60.5651 0 -32.174;-1.3411 0.9929 0;0.00018 43.2541 -0.86939 0;0 0 1 0];B = [-2.516 -13.136;-0.1689 - -0.2514;-17.251 - -1.5766;0 0];C = [0 1 0 0;0 0 0 1]; D = [0 0; 0 0];

线性模型的输入、状态和输出代表了它们在非线性模型工作点的标称值的偏差。

这里，状态变量是:

操纵的变量是升力角和襟副翼角，以度为单位。攻角和俯仰角是测量输出，需要调整。

创建植物，并将初始状态指定为零。

plant = ss(A,B,C,D);X0 = 0 (4,1);

开环系统不稳定。

潮湿(植物)

极阻尼频率时间常数(rad/秒)(秒)-7.50e-03 + 5.56e-02 1.33e+02 -7.50e-03 - 5.56e-02 1.34e-01 5.61e-02 1.33e+02 5.45e+ 02 5.45e+ 02 5.45e+00 -1.83e-01 -7.66e+00 1.00e+00 7.66e+00 1.30e-01

指定控制器约束

两个被操纵的变量都被限制在+/- 25度之间。使用比例因子来促进MPC调优。刻度因子的典型选择是工作范围的上/下限。

MV = struct(“最小值”{-25、-25},“马克斯”25},{25日,“ScaleFactor”, {50,50});

指定工厂输出的比例因子。

OV = struct(“ScaleFactor”, {60 60});

指定控制器调优权重

控制任务是在避免输入饱和导致的不稳定性的同时，使分段常数参考的偏移量为零。因为MV和OV变量都已经在MPC控制器中进行了缩放，所以MPC权重是无量纲的，并应用于缩放后的MV和OV值。对于本例，通过指定比用于俯仰角的权重更大的权重来强调攻击角的跟踪。

权重= struct(“MV”[0 0],“MVRate”(0.1 - 0.1),机汇的, 200年[10]);

创建MPC控制器

创建一个MPC控制器，使用指定的工厂模型，采样时间为0.05秒(20赫兹)，预测水平为10步，控制水平为2步，以及前面指定的权重、约束和比例因子。

mpcobj = mpc(plant,0.05,10,2,Weights,MV,OV);

计算闭环直流增益矩阵

计算闭环稳态输出灵敏度。零值意味着被测设备的输出可以跟踪所需的输出参考设定值。

cloffset (mpcobj)

——>转换模型到离散时间。假设在测量的输出通道1中添加的输出干扰是集成白噪声。假设输出干扰添加到测量输出通道2是集成白噪声。——>”模式。“噪音”属性为空。假设每个测量输出都有白噪声。Ans = 1.0e-11 * 0.1085 -0.0058 0.1224 -0.0355

使用Simulink®进行模拟万博1manbetx

利用Si万博1manbetxmulink分别对攻击角和俯仰角参考信号进行0.1度和2度步长闭环响应仿真。

打开Simulin万博1manbetxk模型并设置MPC控制器财产mpcobj．对于本例，已经设置了属性。

mdl =“mpc_aircraft”；open_system (mdl)

使用Simulink从命令行模拟系统万博1manbetxsim卡命令。

sim (mdl)

打开显示操作变量和飞机输出响应的范围

open_system (“mpc_aircraft / MV”) open_system ('mpc_aircraft/攻角响应') open_system ('mpc_aircraft/俯仰角响应'）

如预期的那样，闭环响应对两个通道都显示出良好的设定点跟踪性能。

参考文献

P. Kapasouris, M. Athans, G. Stein，“带饱和执行器的不稳定装置反馈控制系统设计”，IFAC Symp。非线性控制系统设计佩加蒙出版社，第302—307页，1990年

A. Bemporad, A. Casavola, E. Mosca，“基于预测参考管理的约束线性系统的非线性控制”，IEEE®反式。自动控制，卷AC-42, no。3，第340-349页，1997。

bdclose (mdl)