%《MATLAB在電子信息課程中的應用》第六章例6.22程序q622
% 求離散系統狀態轉移矩陣的解
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear
A = input('離散狀態方程系數A=(N*N方陣） 書上取[0.5,0;0.25,0.25] ');
B = input('離散系統狀態方程系數B=(N*1列陣） 書上取[1;0] ');
x0 = input('初始狀態變量x0=(N*1數組） 書上取[-1;0.5] ');
n = input('要計算的步數n= 書上取10] ');
u = input('輸入信號u=（長度為n的數組） 書上取[0,0.5*ones(1,n-1)] ');
x(:,1)= x0;		% 因MATLAB下標不能取0，x(1)對應于k=0處的x
for i=1:n					% 遞推計算n次
   x(:,i+1)= A*x(:,i)+B*u(i);
end
% 畫出前兩個狀態變量的曲線，注意如何使下標還原至k=0為起點
subplot(2,1,1),stem([0:n],x(1,:)),
line([0,10],[0,0])
subplot(2,1,2),stem([0:n],x(2,:)),
line([0,10],[0,0])