%求解正定矩陣P
E=[1?-1?0?0?;1?0?-1?0?;1?0?0?-1?];
L=[1?0?-1?0;-1?1?0?0;-1?-1?2?0;0?-1?-1?2];
F=[0?0?0;?-1?0?0;0?-1?0;0?0?-1];
H=E*L;
M=H*F;
Q=-M;
setlmis（[]）;??????????????????????????
P=lmivar（1[3?1]）;?????????????????????????????
lmiterm（[1?1?1?P]1Q‘s‘）;
lmiterm（[-2?1?1?P]11）;
lmis=getlmis;????????????????????
[tminxfeas]=feasp（lmis）????????
PP?=?dec2mat（lmisxfeasP）
eig（PP）
%定義系統并賦予初始值
A=-L;B=[0;0;0;0];?D=B;?C=eye（4）;?????????????????????????????????
x0=[1.1;-1.2;-0.5;2];?????????????????????????
[yxt]=initial（ABCDx010）;
%每個智能體的狀態軌跡圖
figure（1）;
plot（tx（:1）‘b‘tx（:2）‘-.‘tx（:3）‘--‘tx（:4））;?
hleg?=?legend（‘x1‘‘x2‘‘x3‘‘x4‘）;
grid;?
title（‘Response?to?Initial?Condition‘）;???
xlabel（‘Time?（sec）‘）;??????????????????
ylabel（‘state‘）;
%智能體間的狀態差軌跡圖
figure（2）;
plot（tx（:2）-x（:1）tx（:3）-x（:1）tx（:4）-x（:1））;
hleg?=?legend（‘e1‘‘e2‘‘e3‘）;
grid;?????????????????????????????????????
title（‘Response?to?Initial?Condition‘）;???
xlabel（‘Time?（sec）‘）;??????????????????
ylabel（‘error‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????485??2018-11-14?08:11??ExperFigurePlot.m
?????文件?????????354??2018-06-15?21:36??feasptest.m
?????文件????????1090??2018-06-15?22:19??AgentTest1.m