用matlab求系统的频响特性,matlab中控制系统的频域分析

求取系统对数频率特性图(波特图)：bode()

求取系统奈奎斯特图(幅相曲线图或极坐标图)：nyquist()

bode(a,b,c,d)：自动绘制出系统的一组Bode图，它们是针对连续状态空间系统[a,b,c,d]的每个输入的Bode图。其中频率范围由函数自动选取，而且在响应快速变化的位置会自动采用更多取样点。

bode(a,b,c,d,iu)：可得到从系统第iu个输入到所有输出的波特图。

bode(num,den)：可绘制出以连续时间多项式传递函数表示的系统的波特图。

bode(a,b,c,d,iu,w)或bode(num,den,w)：可利用指定的角频率矢量绘制出系统的波特图。

当带输出变量[mag,pha,w]或[mag,pha]引用函数时，可得到系统波特图相应的幅值mag、相角pha及角频率点w矢量或只是返回幅值与相角。相角以度为单位，幅值可转换为分贝单位：magdb=20×log10(mag)

nyquist(a,b,c,d)：绘制出系统的一组Nyquist曲线，每条曲线相应于连续状态空间系统[a,b,c,d]的输入/输出组合对。其中频率范围由函数自动选取，而且在响应快速变化的位置会自动采用更多取样点。

nyquist(a,b,c,d,iu)：可得到从系统第iu个输入到所有输出的极坐标图。

nyquist(num,den)：可绘制出以连续时间多项式传递函数表示的系统的极坐标图。

nyquist(a,b,c,d,iu,w)或nyquist(num,den,w)：可利用指定的角频率矢量绘制出系统的极坐标图。

当不带返回参数时，直接在屏幕上绘制出系统的极坐标图(图上用箭头表示w的变化方向，负无穷到正无穷)

。当带输出变量[re,im,w]引用函数时，可得到系统频率特性函数的实部re和虚部im及角频率点w矢量(为正的部分)。可以用plot(re,im)绘制出对应w从负无穷到零变化的部分。

margin：求幅值裕度和相角裕度及对应的转折频率

freqs：模拟滤波器特性

nichols：求连续系统的尼科尔斯频率响应曲线(即对数幅相曲线)

ngrid：尼科尔斯方格图

幅值裕度是在相角为-180度处使开环增益为1的增益量，如在-180度相频处的开环增益为g，则幅值裕度为1/g；若用分贝值表示幅值裕度，则等于：-20*log10(g)。类似地，相角裕度是当开环增益为1.0时，相应的相角与180度角的和。

margin(mag,phase,w)：由bode指令得到的幅值mag(不是以dB为单位)

、相角phase及角频率w矢量绘制出带有裕量及相应频率显示的bode图。

margin(num,den)

：可计算出连续系统传递函数表示的幅值裕度和相角裕度并绘制相应波特图。类似，margin(a,b,c,d)可以计算出连续状态空间系统表示的幅值裕度和相角裕度并绘制相应波特图。

[gm,pm,wcg,wcp]=margin(mag,phase,w)：由幅值mag(不是以dB为单位)

、相角phase及角频率w矢量计算出系统幅值裕度和相角裕度及相应的相角交界频率wcg、截止频率wcp，而不直接绘出Bode图曲线。

h=freqs(b,a,w)用于计算模拟滤波器的幅频响应，其中实矢量w用于指定频率值，返回值h为一个复数行向量，要得到幅值必须对它取绝对值，即求模。

[h,w]=freqs(b,a)自动设定200个频率点来计算频率响应，这200个频率值记录在w中。

[h,w]=freqs(b,a,n)设定n个频率点计算频率响应。

不带输出变量的freqs函数，将在当前图形窗口中绘制出幅频和相频曲线，其中幅相曲线对纵坐标与横坐标均为对数分度。