clear all; close all; s = tf('s'); for c=0.8:0.05:1.2 Ptilde = 1/((s - 1)*(c*s + 1)); figure(1); bode(Ptilde, 'b'); hold on; figure(2); pzmap(Ptilde); hold on; end P0 = 1/((s - 1)*(s + 1)); %P0 = 1/((s - 1)*(1.2*s + 1)); figure(1); bode(P0, 'r');