work4.m

matlab例程

s=tf('s');
sys=1/s/(s*s+2*s+16);
sys1=feedback(sys,1);
nyquist(sys)