H无穷-鲁棒性LMI学习-matlab小作业5.3

本次作业参考课程里边教材的题目null控制系统设计的线性矩阵不等式方法共计39条视频，包括：控制系统设计的线性矩阵不等式方法引论（上）、控制系统设计的线性矩阵不等式方法引论（下）、2-线性矩阵不等式与LMI相关例子（上）等，UP主更多精彩视频，请关注UP账号。https://www.bilibili.com/video/BV1jt411U7xj/

 该教材给出了H无穷一般情况下LMI表达式

基于此，我们选择2情况下的LMI矩阵，完成题目5.3

A = [-5 1 2; 1 -9 1; -1 -10 -3];
B = [0; 1; 0];
C = [1 0 0; 0 0 1];
D = [1; 1];

setlmis([]);  % Initialize the LMI system

% Define the variables
gamma = lmivar(1, [1 0]);  % Gamma (scalar)
X = lmivar(1, [3 1]);  % State-feedback matrix (symmetric matrix)

% LMI for the H-infinity norm
LMI=newlmi;  % 指定下三角元素
lmiterm([LMI 1 1 X], 1, A','s');  % 第一行第一列，XA^T,右乘A^T, 对称元素，AX
lmiterm([LMI 2 1 0], B');  % 第二行第一列，无变量,B^T
lmiterm([LMI 2 2 gamma], -1, 1);  % 第二行第二列，变量gamma，左乘-1，右乘1
lmiterm([LMI 3 1 X], C, 1);     % 第三行第一列，CX，左乘C，右乘1
lmiterm([LMI 3 2 0], D);    % 第三行第二列，无变量，D
lmiterm([LMI 3 3 gamma], -1, 1);  % 第三行第三列，变量gamma，左乘-1，右乘1

LMISYS = getlmis;

% Solve the LMI
[tmin, xfeas] = feasp(LMISYS);

% Display the results
gamma_min = dec2mat(LMISYS, xfeas, gamma);
disp(['The minimal gamma is: ', num2str(gamma_min)]);

得到结果

附上H2情况下LMI简单系统计算gamma的表达式