- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验5伯德图串联校正设计
一实验要求
通过绘制系统校正前后的伯德图 ,明确幅值稳定裕量、相角稳定裕量的定义 ,观察相位超
前、滞后环节对系统伯德图的影响,掌握用伯德图串联校正系统。
二实验步骤
1串联超前校正
掌握计算出幅值稳定裕量、相角稳定裕量 以及对应的频率的函数margin ()及其参数的
使用方法。
在Matlab中输入下面例子的程序,观察并记录结果,分析校正结果。
例:设单位负反馈系统的开环传递函数为
G°
G°(s)
K
s(0.1s 1)
要求系统的静态速度误差系数 Kv 100s 1,相角稳定裕量 550 ,增益稳定裕量
kg 10dB,试确定串联校正装置。
Matlab命令窗口输入:
n um=100;%确定静态速度误差系数 Kv=100
den=[0、1 1 0];
g0=tf(num,den);% 求 G(s)
dpm=55+7、5;%满足相角稳定裕量
[mag,phase,w]=bode(g0);
Mag=20*log10(mag);% 求 L( w ),单位 dB
[Gm,Pm,Wcg,Wcp]=margin(g0);% 求满足Kv时系统的相角稳定裕量 Gm与增益稳定裕
量Pm
phi=(dpm-Pm)*pi/180;% 求最大相位超前角 m,单位弧度
alpha=(1+si n( phi))/(1-si n( phi));% 求校正强度
1
Mn=-10*log10(alpha);% 求 —L( ) 10lg
2 m
Wcgn=spli ne(Mag,w,M n);% 求 m
T=1/Wcgn/sqrt(alpha);% 求 T
Tz=alpha*T;Gc=tf([Tz,1],[T 1]);bode(Gc);hold on% 求 Gc(s)
bode(g0);hold on% 绘制校正前系统的 bode图
bode(g0*Gc);grid on ;margi n(g0*Gc)
Gc%合出校正装置的传递函数
2串联滞后校正
(1)掌握系统对数频率特性曲线 (Bode)图绘制的函数 bode()及其参数的使用方法。(可通过
help方法)
Ks(0.04s 1):K 100, 450。单位dB(2)在Matlab中输入下面例子的程序,
K
s(0.04s 1)
:K 100, 450。
单位dB
G°(s)
设计串联校正装置,使系统满足下列指标
Matlab命令窗口输入:
num=100;
den=[0、04 1 0];
gO=tf( nu m,de n);
[mag,phase,w]=bode(gO);
magdb=20*log10(mag);% 求 L( w ),
wcg=21;%初选截止频率w c
gr=spline(w,magdb,wcg); % 计算 L( w c)
beta=10A(gr/20);% 计算B
T=10/wcg;gc=tf([T 1],[beta*T,1]);
FO=feedback(gO,1);
F=feedback(g0*gc,1);
figure(1);
margin( g0);grid on;
figure(2);
bode(gO,gO*gc);hold on; marg in( g0*gc);grid on;
figure(3);
step(F0,F);
gc
三思考题
(1)编程实现课本上第六章例 6-1与例6-2,要求给出校正装置的传递函数 ,校正后的相
角稳定裕量与增益稳定裕量,并绘制系统校正前与校正后的 Bode图,分析校正结果。
1、 程序如下:
num=100;% 确定静态速度误差系数 Kv=100
den=[0 、001 0 、11 1 0];
g0=tf( num,de n);% 求 G0(s)
[mag,phase,w]=bode(g0);
Mag=20*log10(mag);% 求 L( w),单位 dB
wcg=50;
pag=spli ne(w,phase,wcg);
dpm=30-180-pag+5;
phi=dpm*pi/180;
alpha=(1+si n( phi))/(1-si n( phi));% 求校正强度
T=1/sqrt(alpha)/wcg;% 求 T
Tz=alpha*T;Gc=tf([Tz,1],[T 1]);
bode(g0*Gc);grid on ;hold on
margi n(gO*Gc)
Gc%给岀校正装置的传递函数
Tran sfer function:
0、05533 s + 1
0
007229S+1
Bode DiagramGm = 11.3 dB (at 109 rad/sec) , P m = 35.7 deg (at 48.9
原创力文档


文档评论(0)