s = tf('s'); G11 = -0.1/(18*s+1)*exp(-10*s); G12 = 0.02/(35*s^2+33*s+1)*exp(-18*s); G13 = 0.3/(34*s^2+29*s+1)*exp(-4*s); G21 = 0.8/(30*s+1)*exp(-15*s); G22 = (0.02-0.12*s)/(280*s^2+33*s+1); G23 = -0.25*s/(20*s+1)*exp(-12*s); Kolonne = [G11 G12 G13; G21 G22 G23]; % Matrix der Gij(s) prb_signal = idinput([5000 3],'PRBS',[0 0.03],[-1 1]); figure(1) subplot(3,1,1) plot(prb_signal(:,1)) subplot(3,1,2) plot(prb_signal(:,2)) subplot(3,1,3) plot(prb_signal(:,3)) zeit = 0:1:4999'; %testsignal = [zeit prb_signal]; %testsignal = [zeit prb_signal]; %testsignal = [zeit prb_signal]; y = lsim(Kolonne,prb_signal,zeit); figure(2) subplot(3,1,1) plot(y(:,1)) subplot(3,1,2) plot(y(:,2)) subplot(3,1,3) plot(y(:,3))