控制系统的频率特性分析.docx

文档介绍：实验 控制系统的频特性分析
1-已知系统传递函数为：G(s)= ，要求：
+1
(1) 使用 simulink 进仿真， 改变正弦输入信号的频， 用示波器观察输 出信号，记录同频下输出信号与输入信号的幅值比和相位差，即
£
1:.
o
lllp「£dE-=]
File Edit View Insert Tools Desktop Window 旦Elp 竺
□ = 口寻「氐[題返「口 bT® 口
. o .
lllp「£_dE-=]
2)使用Mat lab函数bode ()绘制系统的^寸数频特性曲线，记录系统校正前 后的幅值裕和相角裕。
k=1/;
G0=zpk([],[0 -1],k);
[h0,r,wx,wc]=margin(G0);
wm=;
L=bode(G0,wm);
Lwc=20*log10(L);
a=10'(- *Lwc);
T=1/(wm*sqrt(a));
phi=asin((a-1)/(a+1));
Gc=(1/a)*tf([a*T 1],[T 1]);
Gc=a*Gc;
G=Gc*G0;
bode(G,'r',G0,'b--');grid;
[h,r,wx,wc]=margin(G)
2．联滞后校正
系统设计要求见课本题 6-4，要求按题目要求设计合的滞后校正环节，并完
成以下内容
1)用 matlab 画出系统校正前后的阶跃相应，并记录系统校正前后的超调及调 节时间
2)使用Mat lab函数bode ()绘制系统的^寸数频特性曲线，记录系统校正前 后的幅值裕和相角裕。
num=30; num=30;
den=conv([1 0],conv([ 1],[ 1])); den=conv([1 0],conv([
1],[ 1]));
Gc=tf(num,den);
G1=tf(num,den);
G1=feedback(Gc,1);
Gd=tf([,1],[41,1]);
Gd=tf([,1],[41,1]);
Ge=tf(num,den);
Ge=tf(num,den);
Gs=series(Gd,Ge);
G2=series(Gd,Ge);
G2=feedback(Gs,1);
subplot(2,1,1);step(G1);grid;
subplot(1,2,1);bode(G1);grid;
subplot(2,1,2);step(G2);grid;
subplot(1,2,2);bode(G2);grid;
File 丨 Edit View Insert Tools Desktop Window Help
File 1 Edit View Insert Tools Desktop Window Help
D已凶宀R现钦紗⑥陽謠•復□ E ■口
®IP< ii
已日©|氐|现弦巴⑥型凰| 口匡1|■四
mp®PNU6ew
Time (sec)
Frequency (rad/sec)
o o 5 o 5 o -15-9-135-18-22-27
(6ep) ®seud
Frequency (rad/sec)
3．联超前—滞后校正 系统设计要求见课本 容
1)用 matlab 画出系统校正前 节时间 2)使用
要求设计合的超前—滞后校正环节，并完成以下内
后的阶跃相应，并记录系统校正前后的超调及调
Matlab函数bode ()绘制系统的^寸数频特性曲线，记录系统校正前
num=180;
1],[
1])); den=conv([1 0],conv([1/6
G1=tf(num,den);
num1=conv([ 1],[ 1]);
den1=conv([64 1],[ 1]);
后的幅值裕和相角裕。 num=180;
den=conv([1 0],conv([1/6 1],[ 1]));
Gc=tf(num,den); G1=feedback(Gc,1);
num1=conv([ 1],[ 1]); den1=conv([64 1],[ 1]