卫星对地定向动力学模型控制器建模仿真MATLAB SIMUlink

卫星对地定向动力学模型控制器建模仿真MATLAB SIMULINK

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  卫星对地定向动力学模型控制器建模仿真\