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))