实验六Matlab频域分析

1、实验六 MATLAB频域分析5.1频率特性的概念系统的频率响应是在正弦信号作用下系统的稳态输出响应。对于线性定常系 统，在正弦信号作用下，稳态输出是与输入同频率的正弦信号， 仅是幅值和相位 不同。设系统传递函数为G(s)，其频率特性为G(j ) G(s) |s j例5-1对系统G(s)飞，在输入信号r(t) sint和r(t) sin3t下可由s 2s 3Matlab求系统的输出信号，其程序如下：num=2;den=1 2 3;G=tf(num,den);t=0:0.1:6*pi;u=sin(t);/ u=sin(3*t);y=lsim(G,u,t);plot(t,u,t,y)运行程序显示系统

2、响应如图5-1所示。a) sin t的响应b)sin (3t)的响应图5-1正弦信号输入系统的稳态响应5.2用nyquist (sys)绘制极坐标图频率特性中的奈奎斯特图是奈奎斯特(Nyquist)稳定性判据的基础。反馈 控制系统稳定的充分必要条件为：奈奎斯特曲线逆时针包围(1,j0)点的次数等于系统开环右极点个数。调用Matlab中nyquist()函数可绘出奈奎斯特图，其调用格式为：re,im, 3 = nyquist(num,den, w)或 sys = tf(nu m,de n); nyq uist(sys)式中，G(s) num/den ；为用户提供的频率范围；re为极坐标的实部；i

3、m为极坐标的虚部。若不指定频率范围，则为nyquist (n um,de n)。在输入指令中，如果缺省了左边的参数说明，奈奎斯特函数将直接生成奈奎斯特图；当命令包含左端变量时，即re,im, w = nyquist(num,den)时，则奈奎斯特函数将只计算频率 响应的实部和虚部，并将计算结果放在数据向量re和im中。在此情况下，只有调用plot函数和向量re、im，才能生成奈奎斯特图。例5-2设系统的传递函数为G(s)丄，绘制其奈奎斯特图图5-3奈奎斯特局部图-5 -解程序如下:num=1;den=1,2,2;nyquist(num,den)运行程序，显示奈奎斯特曲线如图5-2所示。plJI

4、tshAl-I01-i值得注意的是，由于nyquist ()函数自动生成的坐标尺度固定不变，nyquist ()函数可能会生成异常的奈奎斯特图，也可能会丢失一些重要的信息。在这种情 况下，为了重点关注奈奎斯特图在点(-1, jO)附近的形状，着重分析系统的稳 定性，需要首先调用轴函数axis(),自行定义坐标轴的显示尺度，以提高图形的分 辨率；或用放大镜工具放大，以便进行稳定性分析。例5-3设某系统的传递函数为G(s)飞 2，则绘制其奈奎斯特s 8s 17s 10图的程序如下：num=1000;den=1,8,17,10;nyquist (num,den);grid或num=1000;de n

5、=1,8,17,10;sys=tf( nu m,de n);nyq uist (sys);grid运行程序，显示奈奎斯特曲线如图 5-3 a)所示。可以看出在点(-1, j0)附 近,奈奎斯特图很不清楚，可利用放大镜对得出的奈奎斯特图进行局部放大，或利用如下Matlab命令v=-10,0,-1.5,1.5;axis(v)b)例5-4设某系统的开环传递函数为G(s)10(s 2)22(s 1)(s 2s 9)则绘制其奈奎斯特图的程序如下： num=10*1,4, 4;den=conv(1, 1,1,-2,9); nyquist (num,den); grid或num=10*1,4,4; den=

6、conv(1,1,1,-2,9);sys=tf(nu m,de n);nyq uist (sys);grid运行程序，显示得图5-4 a)。若规定实轴、虚轴范围(10, 10), (-10, 10),则绘制其奈奎斯特图的程序如下: num=10*1 4 4; den=conv(1 1,1 -2 9); nyquist (num,den); axis(-10,10,-10,10)运行程序，显示得图5-4 b)。0&*-7 0-2 KQ Aluursl匕-NcfUF3ilUD-9Ti194 E i -1 D 24( Qi 19彌3a)5.3用bode(sys)画对数坐标图伯德图由对数幅频和对数相频

7、两张图构成，轴采用对数分度，幅值为对数增益即分贝(dB)，相位()为线性分度。Matlab中绘制伯德图的函数为bode(), 其调用格式为mag,phase,w=bode( nu m,de n,w)或sys=tf( nu m,de n);bode(sys)式中，G(s)=num/den,频率自动选择范围从0.1到1000rad /s，若自行选择频率范围，可应用logspace ()函数，其格式为w = logspace(a,b ,n)式中，a表示最小频率10a，b表示最大频率10b，n表示10a10b之间频率点数。例5-5设某系统的传递函数为G(s)攀心，则绘制其s(0.5s 1)( s s

8、1)50250伯德图的Matlab程序如下： num=5*0.1 1; den=conv(1 0,conv(0.5 1,1/2500 0.6/50 1); bode( nu m,de n)50-270Bode Diagram-150-90-135-180-225-1010110210Freque ncy (rad/sec)10m00di-310函数 log space(a,b,n)的应用。程序如下: w=logspace(-1,4,300);%确定频率范围及点数 mag,phase,w=bode( nu m,de n,w); semilogx(w,20*log(mag);grid %绘图坐标及

9、大小 xlabel(Freque ncyrad/s),ylabel(20*log(mag)图5-6伯德图-7 -5.4用margin (sys)计算增益裕量和相位裕量Matlab中采用裕量函数 margin ()来确定相对稳定性，其调用格式为Gm,Pm, Wcg, Wcp=margi n(sys) 或margi n(sys)式中，Gm为增益裕量；Pm为相为裕量；Wcg为相角穿越频率；Wcp为幅 值穿越频率。在输入指令中，如果缺省了左边的参数说明， margin ()函数将在伯 德图上自动标注系统的增益裕量和相位裕量.例5-7设某系统的传递函数为量和相位裕量的程序如下：G(s)s3 2s； s

10、0.5，则计算其增益裕-2-9 - num=0.8;de n=1 2 1 0.5; sys=tf( nu m,de n); margi n(sys)执行程序，显示该系统的伯德图及相对稳定裕度如图5-7所示。50edun aaMLyeace san pBode DiagramGm = 5.47 dB (at 1 rad/sec) , P m = 21.9 deg (at 0.788 rad/sec)0-50-100-1500-90-180-270-110010Frequency (rad/sec)11010210图5-7例5-7中系统的伯德图及相对稳定裕度飞的Bode图，并验证其谐振频s2 s

11、25Bode Diagram20T0-20-40-45-600System: sysFrequency (rad/sec): 7.72Magnitude (dB): -3.01-90LyeoTesanp-135-180-110010110210若执行Gm,Pm, Wcg, Wcp=margi n(sys)则可得 Gm = 1.8772,Pm = 21.9176 Wcg = 1.0004, Wcp = 0.7881纯滞后环节的表示Sys.iodelay=0.1;Margi n(sys);闭环系统带宽、谐振峰值、谐振频率 BW=ba ndwidth(sys)作业：MP8.1 用 MATLAB 绘制

12、 T s 率为5rad / s，谐振峰值为14dB。 num=25;de n=1 1 25; sys=tf( nu m,de n); bode(sys) BW=bandwidth(sys)，运行结果BW = 7.7112System: sysFrequency (rad/sec): 4.96Magnitude (dB): 14-11 -Frequency (rad/sec)MP8.2 先手工绘制下列传递函数的 Bode图，然后用MATLAB加以验证(a) G s1s 1 s 10(b) G ss 10s 1 s 20(c) G s1s2 2s 50(d) G ss 52s 1 s 12s 50

13、MP8.3某单位负反馈系统的开环传递函数为25s s 2用MATLAB绘制闭环系统的Bode图，根据Bode图估计闭环带宽，并在图 上标注所得结果。MP8.4某2阶系统框图如图MP8.4所示,阳 _石而g图MP8.42阶反馈控制系统(a)在 0.1到 1000rad / s之间，用logspace函数生成系统闭环Bode图。根据该Bode图，估计系统的谐振峰值 M PM、谐振频率r、和带B ；(b) 用图8.11估算系统的阻尼系数和固有频率n ；(c) 根据闭环传递函数计算 和n的精确值，并与(b)的结果作比较MP8. 5考虑图MP8.5给出的闭环反馈系统，用MATLAB绘制系统的开环和闭环Bode图。)02$ + 10石丰 i)Ci * 2)图MP8.5闭环反馈系统MP9.1某单位负反馈系统的开环传递函数为 G s2100/ s 4s 10 ,试用MATLAB程序验证