-.z.實驗報告課程線性系統理論根底實驗日期2016年6月6日專業班級**同組人實驗名稱全維狀態觀測器的設計評分批閱教師簽字一、實驗目的1.學習用狀態觀測器獲取系統狀態估計值的方法，了解全維狀態觀測器的極點對狀態的估計誤差的影響；2.掌握全維狀態觀測器的設計方法；3.掌握帶有狀態觀測器的狀態反應系統設計方法。二、實驗容開環系統，其中用狀態反應配置系統的閉環極點：；設計全維狀態觀測器，觀測器的極點為：；研究觀測器極點位置對估計狀態逼近被估計值的影響；求系統的傳遞函數〔帶觀測器及不帶觀測器時〕；繪制系統的輸出階躍響應曲線。三、實驗環境MATLAB6.5四、實驗原理〔或程序框圖〕及步驟利用狀態反應可以使閉環系統的極點配置在所希望的位置上，其條件是必須對全部狀態變量都能進展測量，但在實際系統中，并不是所有狀態變量都能測量的，這就給狀態反應的實現造成了困難。因此要設法利用的信息(輸出量y和輸入量*)，通過一個模型重新構造系統狀態以對狀態變量進展估計。該模型就稱為狀態觀測器。假設狀態觀測器的階次與系統的階次是一樣的，這樣的狀態觀測器就稱為全維狀態觀測器或全階觀測器。設系統完全可觀，則可構造如圖4-1所示的狀態觀測器圖4-1全維狀態觀測器為求出狀態觀測器的反應ke增益，與極點配置類似，也可有兩種方法：方法一：構造變換矩陣Q，使系統變成標準能觀型，然后根據特征方程求出ke；方法二：是可采用Ackermann公式：，其中為可觀性矩陣。利用對偶原理，可使設計問題大為簡化。首先構造對偶系統然后可由變換法或Ackermann公式求出極點配置的反應k增益，這也可由MATLAB的place和acker函數得到；最后求出狀態觀測器的反應增益。五、程序源代碼、實驗數據、結果分析〔a〕源程序：A=[010;001;-6-116];B=[0;0;1];C=[100];D=0;P1=[-2+2*sqrt(3)*i;-2-2*sqrt(3)*i;-5];K1=place(A,B,P1)sysnew=ss(A-B*K1,B,C,D)運行結果：K1=74.000025.000015.0000a=*1*2*3*1010*2001*3-80-36-9b=u1*10*20*31c=*1*2*3y110d=u1y10〔b〕源程序：A=[010;001;-6-116];B=[0;0;1];C=[100];D=0;P2=[-5+2*sqrt(3)*i;-5-2*sqrt(3)*i;-10];K2=acker(A',C',P2);L=K2'Anew=A-L*C運行結果：L=262821770Anew=-2610-28201-1776-116〔c〕研究觀測器極點位置對估計狀態逼近被估計值的影響：觀測器極點距離虛軸越近，估計狀態逼近被估計值得速度越快。〔d〕不帶觀測器：源程序：A=[010;001;-6-116];B=[0;0;1];C=[100];D=0;P1=[-2+2*sqrt(3)*i;-2-2*sqrt(3)*i;-5];K1=place(A,B,P1)sysnew=ss(A-B*K1,B,C,D);[num,den]=ss2tf(A-B*K1,B,C,D);Gb=tf(num,den)step(Gb)gridon;title('不帶觀測器的系統的階躍響應曲線');運行結果：K1=74.000025.000015.0000Transferfunction:7.105e-015s^2+1.208e-013s+1--------------------------------------------s^3+9s^2+36s+80帶觀測器：源程序：A=[010;001;-6-116];B=[0;0;1];C=[100];D=0;P1=[-2+2*sqrt(3)*i;-2-2*sqrt(3)*i;-5];K1=place(A,B,P1);sysnew=ss(A-B*K1,B,C,D);P2=[-5+2*sqrt(3)*i;-5-2*sqrt(3)*i;-10];K2=acker(A',C',P2);L=K2';An=[A-B*K1;L*CA-B*K1-L*C]Bn=[B;B]=[C000]Dn=0;[num,den]=ss2tf(An,Bn,,Dn);Go=tf(num,den)step(Go)gridon;title('帶觀測器的系統的階躍響應曲線');運行結果：An=1.0e+003*00.00100000000.0010000-0.0060-0.01100.0060-0.0740-0.0250-0.01500.026000-0.02600.001000.282000-0.282000.00101.770000-1.8500-0.0360-0.0090Bn=001001=100000Transferfunction:-1.137e-013s^4+s^3+20s^2+137s+370-------------------------------------------------------------------------------s^6+