1

a0=[0 0 0 50];
b1=[1 1];
b2=[1 5];
b3=[1 -2];
b0=conv(b1,conv(b2,b3));
H1=tf(a0,b0);
subplot(2,2,1);
bode(H1);title('Q1-Bode diagram of open loop');
xlabel('w');
grid on;
aa=a0;
bb=b0+a0;
H=tf(aa,bb);
subplot(2,2,3);
impulse(H,20);title('Q1-Unit impulse response of closed loop');
xlabel('t');ylabel('c(t)');
grid on;
c0=[0 0 16.7 0];
d0=conv(conv([0.85,1],[0.25,1]),[0.0625,1]);
d1=c0+d0;
G1=tf(c0,d0);
Gs=tf(10*c0,d1);
cc=10*c0;
dd=10*c0+d1;
G=tf(cc,dd);
subplot(2,2,[2 4]);
nyquist(Gs);title('Q2-Nyquist diagram');