Matlab复杂LMI仿真实例

来源：互联网发布：linux培训多少钱编辑：程序博客网时间：2024/06/05 16:02

（1）仿真的论文题目：

A PD-Like Protocol With a Time Delay to Average Consensus Control for Multi-Agent Systems Under an Arbitrarily Fast Switching Topology

Dong Wang, Ning Zhang, Jianliang Wang, Senior Member, IEEE, and Wei Wang

（2）仿真内容：其中提出的Theorem 1，固定时滞下的一种Lyapunov方程

（3）这个式子中的第一个 LMI 比较复杂，一般的 LMI 描述方式难以描述。

分析：矩阵P的大小是3n*3n，将R,W,Q,Y,S都是n*n的，ra 和T矩阵都是4n*3n的。

1.将P矩阵转化为3*3分块矩阵，ra和T分割为4*3的分块矩阵。块矩阵都是n*n的:

%把P定义为9个n*n的子块矩阵[P11,~,SP11]=lmivar(1,[10 1]);[P22,~,SP22]=lmivar(1,[10 1]);[P33,~,SP33]=lmivar(1,[10 1]);[P12,~,SP12]=lmivar(2,[10 10]);[P13,~,SP13]=lmivar(2,[10 10]);[P23,~,SP23]=lmivar(2,[10 10]);[P21,~,SP21]=lmivar(3,SP12');[P31,~,SP31]=lmivar(3,SP13');[P32,~,SP32]=lmivar(3,SP23');%P是3n*3n的[P,~,SP]=lmivar(3,[SP11,SP12,SP13;SP21,SP22,SP23;SP31,SP32,SP33]);

RA=cell(4,3);T =cell(4,3);%把RA和T分成分块矩阵for i=1:4    for j=1:3        RA{i,j}=ra(10*i-9:10*i,10*j-9:10*j);%[1 2 4;1 1 1;2 3 4;]        T{i,j}=t(10*i-9:10*i,10*j-9:10*j);    endend

2.将P,RA,T的分块矩阵带入第一个矩阵不等式乘开：（下面是ra*P*T'的乘开的结果）

3.同理乘开T*P*ra'

4.之后就可以用一般的LMI方法输入第一个LMI的转化结果和其他的LMI，然后求解。

0 0