阶变系统的开环传递函数

阶变系统的开环传递函数

clear all;

Ap=1.68e-2;

In=0.03;

ps=4e6;

pL=2*ps/3;

Ki=188.6;

Vt=2.873e-3;

Kf=1;

bate=6900e5;

m=35000;

Wh=sqrt(4*bate*Ap^2/(m*Vt))

zuni1=0.3;

sys1=tf(1/Ap,[1/Wh^2 2*zuni1/Wh 1 0])

Wsv=157;

zuni2=0.7;

Ksv=1.96e-3;

sys2=tf(Ksv,[1/Wsv^2 2*zuni1/Wsv 1])

%系统的开环传递函数

sys_open=Ki*sys1*sys2

sysclose=feedback(sys_open,1);

figure; %绘制nyquist曲线

subplot(121);pzmap(sys_open);

grid on;

xlabel('实轴');ylabel('虚轴');title('零极点图'); subplot(122);

nyquist(sys_open);

grid on;

xlabel('实轴');ylabel('虚轴');title('Nyquist图'); figure; %时域分析

subplot(121);step(sysclose);

grid on;

xlabel('时间');ylabel('振幅');title('阶跃响应'); subplot(122);impulse(sysclose);

grid on;

xlabel('时间');ylabel('振幅');title('脉冲图响应'); figure; %绘制Bode图及其参数求解

w=logspace(-1,2);

grid on;

margin(sys_open);

xlabel('频率');title('Bode图');

[mag,phase,W]=bode(sysclose,w);

[l,c]=size(mag);

mag1=zeros(c,l);

for i=1:c

mag1(i)=10*log10(mag(1,1,i));

end

disp('-3dB频率：');

W_3dB=abs(interp1(mag1,W,-3,'spline')); grid on;

Kv=Ki*Ksv/Ap;

ef=0.02*0.003/Ki*Kf;

vm=2.2e-2;

er=vm/Kv;

e=ef+er;

[W_3dB e]