final.m

state space discreet

%program untuk mencari y(k) dan mengambar y(k) terhadap k
%x0=inisialisasi awal x
%T_x=(zI-G)^-1
%E_f=F(z)
%U_xz=X(z)
%H_xk=x(k)
%B_yk=y(k)
%B_g=y(k)
syms z
G=[0.527 0.159;-0.477 -0.109];
H=[0.158;0.159];
C=[1 0];
u=z/(z-1);
x0=[1;-1];
a1=-trace(G);
H1=G+a1*eye(2);
a2=-0.5*trace(G*H1);
adj=eye(2)*z+H1;
det=z^2+a1*z+a2;
T_x=adj/det; 
E_f=C*T_x*H; 
G_m=(z*x0)+(H*u);
U_xz=T_x*Gm; 
H_xk=iztrans(xz); 
B_yk=C*xk 
O_n=0:0.1:20;
B_g=200503/600400+317/1900*((1/20).^n)+631/1264*((46/125).^n)
plot(n,B_g,'o')
title('y(k) terhadap k')
xlabel('k')
ylabel('y(k)')
%created by: Teguh "Bob" Budi Wibowo