MATLAB实验三,四_原创文档.pdfVIP

%bode图

figure(1);n=[1];d=[1];sys=tf(n,d);subplot(2,3,1);bode(sys)

n1=[1];d1=[1,1];sys1=tf(n1,d1);subplot(2,3,2);bode(sys1)

n2=[1];d2=[1,0];sys2=tf(n2,d2);subplot(2,3,3);bode(sys2)

n3=[1,0];d3=[1];sys3=tf(n3,d3);subplot(2,3,4);bode(sys3)

n4=[10];d4=[1,2,10];sys4=tf(n4,d4);subplot(2,3,5);bode(sys4)

%Nyquist图

figure(2);n=[1];d=[1];sys=tf(n,d);subplot(2,3,1);nyquist(sys)

n1=[1];d1=[1,1];sys1=tf(n1,d1);subplot(2,3,2);nyquist(sys1)

n2=[1];d2=[1,0];sys2=tf(n2,d2);subplot(2,3,3);nyquist(sys2)

n3=[1,0];d3=[1];sys3=tf(n3,d3);subplot(2,3,4);nyquist(sys3)

n4=[10];d4=[1,2,10];sys4=tf(n4,d4);subplot(2,3,5);nyquist(sys4)

%幅值裕度和相角裕度

n5=20;d5=conv([1,0],[0.5,1]);sys5=tf(n5,d5);figure(3);margin(sys5)

figure(4);nyquist(sys5);v=[-10,10,-20,20];axis(v)

%自定义一个开环传递函数（高阶），绘制其闭环根轨迹。

figure(1);z=[];p=[0,-1,-2];k=1;sys=zpk(z,p,k);rlocus(sys);

z1=[];p1=[-1,-1,-1];k1=1;sys1=zpk(z1,p1,k1);holdon;rlocus(sys1)

%分析闭环系统稳定性，时域响应情况，确定主导极点，计算主要性能指标。

figure(2);G1=1;G=feedback(sys,G1,-1);step(G);

figure(3);z=[];p=[0,-1,-2];k=1;sys=zpk(z,p,k);rlocus(sys);[k,poles]=rlo

cfind(sys)

figure(4);z2=[];p2=[-2.8298];k2=1;sys2=zpk(z2,p2,k2);rlocus(sys2);

figure(5);z3=[];p3=[-0.0851+1.2293i,-0.0851-

1.2293i];k3=1;sys3=zpk(z3,p3,k3);rlocus(sys3);

figure(6);subplot(2,1,1);G1=1;G2=feedback(sys2,G1,-1);step(G2);subplot(

2,1,2);

G1=1;G3=feedback(sys3,G1,-1);step(G3)

%给定开环传递函数，绘制其闭环根轨迹，分析闭环系统稳定性，时域响应情况，

%计算具有最小阻尼比时对应的闭环极点与K1值

figure(7);z4=[-7];p4=[0,-4];k4=1;sys4=zpk(z4,p4,k4);rlocus(sys4);[k,pol

es]=rlocfind(sys)

figure(8);G1=1;G4=feedback(sys4,G1,-1);step(G4)

figure(9);z5=[-7];p5=[-2.1507];k5=1;sys5=zpk(z5,p5,k5);rlocus(sys5);

figure(10);z6=[-7];p6=[-0.5081];k6=1;sys6=zpk(z6,p6,k6);rlocus(sys6);

figure(11);z7=[-7];p7=[-0.3412];k7=1;sys7=zpk(z7,p7,k7);rlocus(sys7);

figure(12);G1=1;G5=feedback(sys5,G1,-1);step(G5);

holdon;G1=1;G6=feedback(sys6,G1,-1);step(G6);

holdon;G1=1;G7=feedback(sys7,G1,-1);step(G7);