?clear
?I=[12.490.670.06;0.6713.850.06;0.060.0615.75];???????%%可以拉偏的轉動慣量7890
?period_torqued=[111]*1e-4;
?const_torqued=[111]*1e-4;
?wn=1;
?ep=1;
?Kp=diag（[0.5^2*12.490.5^2*13.850.5^2*15.75]）;
?Kd=diag（[2*ep*0.5*12.492*ep*0.5*13.852*ep*0.5*15.75]）;
w0=0.001;
A0=0.000015;
Theat=0.7853;
Jsx=0.000637;Jsy=0.000637;Jsz=0.000637;Js=0.000637;
J=diag（[Jsx;Jsy;Jsz]）;%飛輪的轉動慣量
%I=[12.49?0.67?0.06;0.67?13.85?0.06;0.06?0.06?15.75];%%星體的轉動慣量
Altitude=500;	
W0=-631.3481*（6378+Altitude）^（-1.5）	;
i1=12.49;i2=13.85;i3=15.75;
di1=1.49;di2=1.85;di3=1.75;%%星體主軸轉動慣量
%?i1=12.49;i2=i1;i3=i1;
I=[i1?0?0;0?i2?0;0?0?i3];
%?dI=[di1?0?0;0?di2?0;0?0?di3];
A0=0.000015;
%%%%系統矩陣%%%%%%%%%%%%%
sigmaa=0*diag（[1?1?1]）;
sigmas=0*eye（3）;%只考慮飛輪加性故障，因為采用的是狀態反饋，量測故障暫時不考慮
A2=I;
A1=W0*[0?0?-（i

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

????.......??????4648??2014-12-04?14:58??Hoo\H2Hoo.m

????.......???????280??2014-12-04?14:58??Hoo\initial.m

?????文件?????127341??2014-12-04?17:00??Hoo\Rubust_control.mdl

????.......????187522??2014-12-04?14:58??Hoo\Rubust_control.mdl.r2009a

?????文件?????111893??2014-12-04?14:58??Hoo\Rubust_control.mdl.r2010b

????.......??????1343??2014-12-04?14:58??Hoo\with?failure.m

????.......??????1343??2014-12-04?14:58??Hoo\without?failure.m

?????目錄??????????0??2014-12-04?17:00??Hoo

-----------?---------??----------?-----??----

???????????????434370????????????????????8