实验五Matlab频域分析2013

1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流实验五Matlab频域分析2013.精品文档.实验五 MATLAB频域特性分析5.1 频率特性的概念系统的频率响应是在正弦信号作用下系统的稳态输出响应。对于线性定常系统，在正弦信号作用下，稳态输出是与输入同频率的正弦信号，仅是幅值和相位不同。设系统传递函数为，其频率特性为例5-1 对系统，在输入信号和下可由Matlab求系统的输出信号，其程序如下：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）稳定性判据的基础。反馈控制系统稳定的充分必要条件为：奈奎斯特曲线逆时针包围点的次数等于系统开环右极点个数。调用Matlab中nyquist() 函数可绘出奈奎斯特图，其调用格式为：或式中，；为用户提供的频率范围；re为极坐标的实部；im为极坐标的虚部。若不指定频率范围，则为。在输入指令中，如果缺省了左边的参数说明，奈奎斯特函数将直接生成奈奎斯特图；当命令包含左端变量时，即时，则奈奎斯特函数将只计算频率响应的实部

3、和虚部，并将计算结果放在数据向量和中。在此情况下，只有调用plot函数和向量、，才能生成奈奎斯特图。例5-2 设系统的传递函数为，绘制其奈奎斯特图。解 程序如下：num=1;den=1,2,2;nyquist(num,den)运行程序，显示奈奎斯特曲线如图5-2所示。图5-2 奈奎斯特曲线值得注意的是，由于nyquist（）函数自动生成的坐标尺度固定不变，nyquist ()函数可能会生成异常的奈奎斯特图，也可能会丢失一些重要的信息。在这种情况下，为了重点关注奈奎斯特图在点（-1，j0）附近的形状，着重分析系统的稳定性，需要首先调用轴函数axis(),自行定义坐标轴的显示尺度，以提高图形的分辨

4、率；或用放大镜工具放大，以便进行稳定性分析。例5-3 设某系统的传递函数为，则绘制其奈奎斯特图的程序如下：num=1000;den=1,8,17,10;nyquist (num,den);grid或num=1000;den=1,8,17,10;sys=tf(num,den);nyquist (sys);grid运行程序，显示奈奎斯特曲线如图5-3 a)所示。可以看出在点（-1，j0）附近，奈奎斯特图很不清楚，可利用放大镜对得出的奈奎斯特图进行局部放大，或利用如下Matlab命令v=-10,0,-1.5,1.5;axis(v) a) b)图5-3 奈奎斯特局部图例5-4 设某系统的开环传递函数为

5、则绘制其奈奎斯特图的程序如下：num=10*1, 4, 4;den=conv(1, 1,1, -2,9);nyquist (num,den);grid或num=10*1,4,4;den=conv(1,1,1,-2,9);sys=tf(num,den);nyquist (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)。 a） b） 图 5

6、-4 奈奎斯特曲线5.3 用bode(sys)画对数坐标图伯德图由对数幅频和对数相频两张图构成，轴采用对数分度，幅值为对数增益即分贝()，相位为线性分度。Matlab中绘制伯德图的函数为bode( )，其调用格式为mag,phase,w=bode(num,den,w) 或sys=tf(num,den);bode(sys)式中，G(s)=num/den，频率自动选择范围从，若自行选择频率范围，可应用logspace()函数，其格式为 式中，a表示最小频率10a，b表示最大频率10b，n表示10a10b之间频率点数。例5-5 设某系统的传递函数为，则绘制其伯德图的Matlab程序如下:>&g

7、t; num=5*0.1 1;>> den=conv(1 0,conv(0.5 1,1/2500 0.6/50 1);>> bode(num,den)函数 的应用。程序如下：>> w=logspace(-1,4,300); %确定频率范围及点数>> mag,phase,w=bode(num,den,w);>> semilogx(w,20*log(mag);grid %绘图坐标及大小>> xlabel('Frequencyrad/s'),ylabel('20*log(mag)')图5-6 伯德

8、图5.4 用margin（sys）计算增益裕量和相位裕量Matlab中采用裕量函数 margin ()来确定相对稳定性，其调用格式为Gm,Pm, Wcg, Wcp=margin(sys) 或margin(sys)式中，Gm为增益裕量；Pm为相为裕量；Wcg为相角穿越频率；Wcp为幅值穿越频率。在输入指令中，如果缺省了左边的参数说明，margin ()函数将在伯德图上自动标注系统的增益裕量和相位裕量.例5-7 设某系统的传递函数为 ，则计算其增益裕量和相位裕量的程序如下：>> num=0.8;den=1 2 1 0.5;>> sys=tf(num,den);>>

9、; margin(sys)执行程序，显示该系统的伯德图及相对稳定裕度如图5-7所示。 图5-7 例5-7中系统的伯德图及相对稳定裕度若执行Gm,Pm, Wcg, Wcp=margin(sys)则可得Gm = 1.8772,Pm = 21.9176，Wcg = 1.0004，Wcp = 0.7881纯滞后环节的表示Sys.iodelay=0.1;Margin(sys);闭环系统带宽、谐振峰值、谐振频率>> BW=bandwidth(sys)作业：MP8.1 用MATLAB绘制的Bode 图，并验证其谐振频率为5，谐振峰值为14dB。>> num=25;den=1 1 25

10、; sys=tf(num,den); bode(sys)>> BW=bandwidth(sys)，运行结果 BW = 7.7112MP8.2 先手工绘制下列传递函数的Bode图，然后用MATLAB加以验证。MP8.3某单位负反馈系统的开环传递函数为:用MATLAB绘制闭环系统的Bode图，根据Bode图估计闭环带宽，并在图上标注所得结果。MP8.4某2阶系统框图如图MP8.4所示，图MP8.4 2阶反馈控制系统(a) 在到之间，用logspace函数生成系统闭环Bode图。根据该Bode图，估计系统的谐振峰值、谐振频率、和带宽；(b) 用图8.11估算系统的阻尼系数和固有频率；(c) 根据闭环传递函数计算和的精确值，并与(b)的结果作比较。MP8. 5考虑图MP8.5给出的闭环反馈系统，用MATLAB绘制系统的开环和闭环Bode图。图MP8.5 闭环反馈系统MP9.1某单位负反馈系统的开环传递