Gs11=tf (32,[1 6 34]);Gs12=tf ([5 0],[1 6 34]);Gs21=tf ([-4 0],[1 6 34]);Gs22=tf (34,[1 6 34]);t=0:0.001:2.5;y1u1=step(10*Gs11,t);y1u2=step(8*Gs12,t);y2u1=step(10*Gs21,t);y2u2=step(8*Gs22,t);figure(1)subplot(2,2,1)plot(t,y1u1)subplot(2,2,2)plot(t,y1u2)subplot(2,2,3)plot(t,y2u1)subplot(2,2,4)plot(t,y2u2)figure(2)subplot(2,1,1)plot(t,y1u1+y1u2)subplot(2,1,2)plot(t,y2u1+y2u2)G11=[10*Gs11 8*Gs12;10*Gs21 8*Gs22]figure(3)step(G11,t)Y=step(G11,t);figure(4)subplot(2,1,1)plot(t,Y(:,1,1))subplot(2,1,2)plot(t,Y(:,2,2))Gs=[Gs11 Gs12;Gs21 Gs22];u1=10*ones(size(t));u2=8*ones(size(t));U=[u1;u2];figure(4)lsim(Gs,U,t)uu1=5*sign(sin(2/3*pi*t));% varia la frecuenciauu2=4*sin(2/3*pi*t);% varia la frecuenciaUU=[uu1;uu2];figure(5)lsim(Gs,UU,t)