gama=2;????%?魯棒性能參數
A=[0?1;?0?0];??%?等式（4-9）中變量賦值
B1=[0?1]‘;
B2=[0?1]‘;
C0=[1?0];
D0=[0];
C1=[1?0];
D1=[0];
?
setlmis（[]）;???%?開始建立LMI
X?=?lmivar（1[1?1]）;??%?為變量設置存儲空間
W?=?lmivar（2[1?2]）;??%?lmivar的用法可參見4.1.3小節介紹
Y?=?lmivar（1[2?1]）;
lmiterm（[-1?1?1?X]11）;??%?建立線性矩陣不等式對應項
lmiterm（[-1?1?2?Y]C01）;?%?lmiterm的用法可參見4.1.3小節介紹
lmiterm（[-1?1?2?W]D01）;
lmiterm（[-1?2?1?-Y]1C0‘）;
lmiterm（[-1?2?1?-W]1D0‘）;
lmiterm（[-1?2?2?Y]1010）;
lmiterm（[-2?1?1?Y]-A1）;
lmiterm（[-2?1?1?Y]-1A‘）;
lmiterm（[-2?1?1?W]-B21）;
lmiterm（[-2?1?1?-W]-1B2‘）;
lmiterm（[-2?1?1?0]-B1*B1‘）;
lmiterm（[-2?1?2?-Y]1C1‘）;
lmiterm（[-2?1?2?-W]1D1‘）;
lmiterm（[-2?2?1?Y]C11）;
lmiterm（[-2?2?1?W]D11）;
lmiterm（[-2?2?2?0]gama*gama*eye（1））;
LMIs?=?getlmis;???????????%?完成LMI建立
[tminxfeas]?=?feasp（LMIs）??%?求解LMI
X?=?dec2mat（LMIsxfeasX）;?%?取出變量值
W?=?dec2mat（LMIsxfeasW）;
Y?=?dec2mat（LMIsxfeasY）;
K?=?W*inv（Y）????????????%?確定控制器


Jx?=?2000;?Jy?=?1800;?Jz?=?250;
w0?=?2*pi/5400;
sim（‘satellite‘）???%?利用Matlab語言驅動Simulink文件仿真
subplot（311）?plot（simout.time?simout.signals.values（:1）*57.3）?ylabel（‘fai‘）
subplot（312）?plot（simout.time?simout.signals.values（:2）*57.3）?ylabel（‘zeta‘）
subplot（313）?plot（simout.time?simout.signals.values（:3）*57.3）?xlabel（‘t/sec‘）ylabel（‘bosai‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????41525??2013-06-02?21:27??衛星對地定向動力學模型控制器建模仿真\satellite.mdl
?????文件????????1408??2013-06-02?19:54??衛星對地定向動力學模型控制器建模仿真\satellite_duidi.m
?????文件?????????104??2013-06-02?21:29??衛星對地定向動力學模型控制器建模仿真\仿真說明.txt
?????目錄???????????0??2017-03-26?21:55??衛星對地定向動力學模型控制器建模仿真\