%《MATLAB在電子信息課程中的應用》第八章例8.20程序q820
% 離散系統狀態觀測器設計
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear,
Ad=[0,0,-0.5;1,0,0.2;0,1,0.8]; 	    % 設定系統原有參數
Bd=[0;0;1];Cd=[1,0,0]; 			
pbd=[0.5+0.7j,0.5-0.7j,-0.1];		% 要求實現的狀態觀測器極點
Lp=place(Ad',Cd',pbd)				% 用兩種函數來求
La=acker(Ad',Cd',pbd)
pd=eig(Ad);							% 求系統原有極點
figure(1),zplane([ ],pd)		% 畫z平面上原有零極點分布圖
figure(2),zplane([ ],pbd')		% 把極點pbd化為列向量畫圖