M1.3-Matlab频域分析资料课件.ppt
《M1.3-Matlab频域分析资料课件.ppt》由会员分享,可在线阅读,更多相关《M1.3-Matlab频域分析资料课件.ppt(25页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、机械机械机械机械学学学学院院院院机械控制工程机械控制工程控制系统的计算机辅助分析主要内容:主要内容:一、概述二、控制系统的数学描述与建模 三、模型的转换四、模型连接和化简五、控制系统的时域分析六、六、控制系统的频域分析控制系统的频域分析1机械机械机械机械学学学学院院院院机械控制工程机械控制工程控制系统的分析和设控制系统的分析和设计流程计流程物理系统 假 设数 学 模 型数学分析计算机计算机仿真仿真模型响应预测物理系统的期望响应改进系统结构修改系统参数2机械机械机械机械学学学学院院院院机械控制工程机械控制工程q频率响应是指系统对正弦输入信号的稳态响应,从频率响应中可以得出带宽、增益、转折频率、闭
2、环稳定性等系统特征。q频率特性是指系统在正弦信号作用下,稳态输出与输入之比对频率的关系特性。频率特性函数与传递函数有直接的关系,记为:一、频域分析的一般方法q求取系统对数频率特性图(波特图):bode()q求取系统奈奎斯特图(幅相曲线图或极坐标图):nyquist()q频域分析法是应用频率特性研究控制系统的一种典型方法。采用这种方法可直观地表达出系统的频率特性,分析方法比较简单,物理概念比较明确,对于诸如防止结构谐振、抑制噪声、改善系统稳定性和暂态性能等问题,都可以从系统的频率特性上明确地看出其物理实质和解决途经。通常将频率特性用曲线的形式进行表示,包括对数频率特性曲线和幅相频率特性曲线简称幅
3、相曲线,MATLAB提供了绘制这两种曲线的函数。3机械机械机械机械学学学学院院院院机械控制工程机械控制工程1、对数频率特性图(波特图)exp4_10.m exp4_10_.mq对数频率特性图包括了对数幅频特性图和对数相频特性图。横坐标为频率w,采用对数分度,单位为弧度/秒;纵坐标均匀分度,分别为幅值函数20lgA(w),以dB表示;相角,以度表示。MATLAB提供了函数bode()来绘制系统的波特图,其用法如下:qbode(a,b,c,d):自动绘制出系统的一组Bode图,它们是针对连续状态空间系统a,b,c,d的每个输入的Bode图。其中频率范围由函数自动选取,而且在响应快速变化的位置会自动
4、采用更多取样点。qbode(a,b,c,d,iu):可得到从系统第iu个输入到所有输出的波特图。qbode(num,den):可绘制出以连续时间多项式传递函数表示的系统的波特图。qbode(a,b,c,d,iu,w)或bode(num,den,w):可利用指定的角频率矢量绘制出系统的波特图。q当带输出变量mag,pha,w或mag,pha引用函数时,可得到系统波特图相应的幅值mag、相角pha及角频率点w矢量或只是返回幅值与相角。相角以度为单位,幅值可转换为分贝单位:magdb=20log10(mag)4机械机械机械机械学学学学院院院院机械控制工程机械控制工程2、奈奎斯特图(幅相频率特性图)e
5、xp4_11.m exp4_11_.mq对于频率特性函数G(jw),给出w从负无穷到正无穷的一系列数值,分别求出Im(G(jw)和Re(G(jw)。以Re(G(jw)为横坐标,Im(G(jw)为纵坐标绘制成为极坐标频率特性图。MATLAB提供了函数nyquist()来绘制系统的极坐标图,其用法如下:qnyquist(a,b,c,d):绘制出系统的一组Nyquist曲线,每条曲线相应于连续状态空间系统a,b,c,d的输入/输出组合对。其中频率范围由函数自动选取,而且在响应快速变化的位置会自动采用更多取样点。qnyquist(a,b,c,d,iu):可得到从系统第iu个输入到所有输出的极坐标图。q
6、nyquist(num,den):可绘制出以连续时间多项式传递函数表示的系统的极坐标图。qnyquist(a,b,c,d,iu,w)或nyquist(num,den,w):可利用指定的角频率矢量绘制出系统的极坐标图。q当不带返回参数时,直接在屏幕上绘制出系统的极坐标图(图上用箭头表示w的变化方向,负无穷到正无穷)。当带输出变量re,im,w引用函数时,可得到系统频率特性函数的实部re和虚部im及角频率点w矢量(为正的部分)。可以用plot(re,im)绘制出对应w从负无穷到零变化的部分。5机械机械机械机械学学学学院院院院机械控制工程机械控制工程二、常用频域分析函数MATLAB除了提供前面介绍的
7、基本频域分析函数外,还提供了大量在工程实际中广泛应用的库函数,由这些函数可以求得系统的各种频率响应曲线和 特征值。如:margin:求幅值裕度和相角裕度及对应的转折频率freqs:模拟滤波器特性nichols:求连续系统的尼科尔斯频率响应曲线(即对数幅相曲线)ngrid:尼科尔斯方格图6机械机械机械机械学学学学院院院院机械控制工程机械控制工程margin()函数 exp4_12.m exp4_12_.mqmargin函数可以从频率响应数据中计算出幅值裕度、相角裕度以及对应的频率。幅值裕度和相角裕度是针对开环SISO系统而言,它指示出系统闭环时的相对稳定性。当不带输出变量引用时,margin可在
8、当前图形窗口中绘制出带有裕量及相应频率显示的Bode图,其中幅值裕度以分贝为单位。q幅值裕度是在相角为-180度处使开环增益为1的增益量,如在-180度相频处的开环增益为g,则幅值裕度为1/g;若用分贝值表示幅值裕度,则等于:-20*log10(g)。类似地,相角裕度是当开环增益为1.0时,相应的相角与180度角的和。qmargin(mag,phase,w):由bode指令得到的幅值mag(不是以dB为单位)、相角phase及角频率w矢量绘制出带有裕量及相应频率显示的bode图。qmargin(num,den):可计算出连续系统传递函数表示的幅值裕度和相角裕度并绘制相应波特图。类似,margi
9、n(a,b,c,d)可以计算出连续状态空间系统表示的幅值裕度和相角裕度并绘制相应波特图。qgm,pm,wcg,wcp=margin(mag,phase,w):由幅值mag(不是以dB为单位)、相角phase及角频率w矢量计算出系统幅值裕度和相角裕度及相应的相角交界频率wcg、截止频率wcp,而不直接绘出Bode图曲线。7机械机械机械机械学学学学院院院院机械控制工程机械控制工程freqs()函数 exp4_13.mqfreqs用于计算由矢量a和b构成的模拟滤波器H(s)=B(s)/A(s)的幅频响应。qh=freqs(b,a,w)用于计算模拟滤波器的幅频响应,其中实矢量w用于指定频率值,返回值h
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- M1 Matlab 分析 资料 课件
限制150内