二阶系统v电气与电子测量技术.pdf

二阶系统 matlab 程序1： kosai=[0.1:0.1:1.0,2.0]; figure(1) hold on for i=kosai num=1; den=[1,2*i,1]; step(num,den) end title(The Step Response of Two Order System); hold off 仿真结果： The Step Response of Tw o Order System 1.8 1.6 1.4 1.2 e 1 d u t i l p m 0.8 A 0.6 0.4 0.2 0 0 10 20 30 40 50 60 Time (sec) matlab 程序2： fig_size = [232 84 774 624]; zeta = 0.5; % 原始系统参数 wn = 4; num0 = wn^2; den0 = [1 2*zeta*wn 0]; % 原始系统传递函数 [ncl0,dcl0] = feedback(num0,den0,1,1,-1); % t = linspace(0,4,1001); % 时间向量 ys0 = step(ncl0,dcl0,t); % 原始系统的阶跃响应 % figure(1),clf,plot(t,ys0),grid,xlabel(Time (s)),... ylabel(Amplitude),title(Step Response: zeta = 0.5, wn = 4 r/s),... text(1.25,0.7,P.O. = 16.3%, Ts = 2 s),set(gcf,Position,fig_size) % % 为系统增加一个开环零点 % z = [0.2 0.5 1 2 5 10]; % 增加了零点后的系统参数 % for i = 1:length(z) numz(i,:) = (wn^2 / z(i)) * [1 z(i)]; [nclz(i,:),dclz(i,:)] = feedback(numz(i,:),den0,1,1,-1); clpz(:,i) = roots(dclz(i,:)); ysz(:,i) = step(nclz(i,:),dclz(i,:),t); % 增加零点后的阶跃响应 end % figure(2),clf,plot(t,ysz),grid,xlabel(Time (s)),... ylabel(Amplitude),title(Step Response with One Zero in Open-Loop System),... text(0.1,0.98,z = 0.2),text(1.2,1.07,z = 10),set(gcf,Position,fig_size) % % 研究系统增加一个开环极点后的效果 % p = [1 2 5 10 20]; % 增加了一个极点的参数 t = 5*t; % 新的时间参数 % for i = 1:length(p) nump(i) = p(i) * num0; denp(i,:) = conv(den0,[1 p(i)]); [nclp(i,:),dclp(i,:)] = feedback(nump(i),denp(i,:),1,1,-1); clpp(:,i) = roots(dclp(i,:)); ysp(:,i) = step(nclp(i,:),dclp(i,:),t); % 增加极点后的阶跃响应 end % figure(3),clf,plot(t,ysp),grid,xlabel(Time (s)),... ylabel(Amplitude),title(Step Response with Extra Pole in Open-Loop System),... text(2.6,1.7,p = 1),set(gcf,Position,fig_size) % 仿真结果：