P1=10; P2=20; Tr=30; Ts=1; beta=exp(-Ts/Tr); numm=[1]; denm=[3 8 1]; num=[2]; den=[16 9 1]; [A,B,C,D]=tf2ss(num,den); [Gd,Hd,Cd,Dd]=c2dm(A,B,C,D,Ts); [Am,Bm,Cm,Dm]=tf2ss(numm,denm); [Gmd,Hmd,Cmd,Dmd]=c2dm(Am,Bm,Cm,Dm,Ts);