自动控制原理线性系统的频域分析实验报告

实验四控制系统的频域分析一实验目的1. 利用计算机作出开环系统的波特图2. 观察记录控制系统的开环频率特性3. 控制系统的开环频率特性分析二预习要点1.预习Bode图和Nyquist图的画法；2.映射定理的内容；3.Nyquist稳定性判据内容。

三实验方法1、奈奎斯特图（幅相频率特性图）❑对于频率特性函数G(jw)，给出w从负无穷到正无穷的一系列数值，分别求出Im(G(jw))和Re(G(jw))。

以Re(G(jw)) 为横坐标， Im(G(jw)) 为纵坐标绘制成为极坐标频率特性图。

MATLAB提供了函数nyquist()来绘制系统的极坐标图，其用法如下：❑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从负无穷到零变化的部分。

2、对数频率特性图（波特图）对数频率特性图包括了对数幅频特性图和对数相频特性图。

横坐标为频率w ，采用对数分度，单位为弧度/秒；纵坐标均匀分度，分别为幅值函数20lgA(w)，以dB 表示；相角，以度表示。

MATLAB 提供了函数bode()来绘制系统的波特图，其用法如下：❑ bode(a,b,c,d,iu)：可得到从系统第iu 个输入到所有输出的波特图。

实验三、线性系统的频域分析法一，实验目的1，掌握matlab绘制波特图以及奈奎斯特图的方法。

2，学会从波特图以及奈奎斯特图判定系统的稳定性。

3，学会从波特图上求系统的稳定裕度。

4，了解k值变化时对波特图幅频和相频曲线的影响。

5，掌握matalab绘制系统零极点分布图的方法。

6，学会从系统的零极点分布图判断系统的稳定性。

二，实验原理1，从奈奎斯特图判定系统是否稳定的原理奈式稳定判据：反馈控制系统稳定的充分必要条件是半闭合曲线ΓGH不穿过（-1，0j）点，且逆时针包围临界点（-1，0j）点的圈数R 等于开环传递函数正实部极点数P具体方法是，先观察系统传递函数得出系统是否在s平面的右半开平面由极点，得出P的值，在观察曲线从（-1，0j）点右侧穿越的次数，其中自上而下为正穿越，自下而上为负穿越，完整的一次穿越记为N 半次穿越记为0.5N,R=2N=2(N+ -N-) 而Z=P-R，观察Z是否为零，Z 为零则系统是稳定的，Z不为零时则系统是不稳定的。

2，从波特图判定系统是否稳定的原理。

从奈奎斯特稳定判定我们可以知道，要判定系统是否稳定就要观察曲线穿越（-1，0j）点次数，对应在波特图中，当取w=wc时，要满足A（wc）=|G（jwc）H(jwc)|=1 L（wc）=20logA（wc）=0因此wc为分界点，对应到相频曲线上，观察在w<wc时曲线穿越-180度的次数。

然后计算方法和上面相同，既可以判定系统的稳定性。

3，根据系统的零极点分布判断系统稳定性的原理三，实验内容A、设单位负反馈系统的开环传递函数为K（S+1）/S(S+2)(S^2+17S+4000) 其中K=1000（1）绘制波特图。

（2）观察绘制出的bode 图，分析系统的稳定性，并在图上求稳定裕度；（3）绘制K=2000 时系统的bode 图，分析曲线的改变情况，并分析K 值变化时，对系统幅频响应和相频响应的影响。

分析：1，绘制波特图matlab 文本命令为：s=tf(‘s’);G=1000*(s+1)/(s*(s+2)*(s^2+17*s+4000))Bode(G)Grid onMargin(G),2，绘制出的波形为2,由于传递函数中可知v=1所以要在相频中增补从-90度到0度的相频曲线，由波特图可以看出当L（w）=0dB时对应的频率值为wc，在w<wc 时，在相频曲线中没有穿越-180度，所以可以知道R=0，又由传递函数可以知道P=0，所以Z=0，从而我们知道系统此时是稳定的，由裕度函数我们可以在图中求出幅值裕度Gm=36.7dB，相角裕度Pm=93.5度.剪切频率wc=0.126rad/s.3，改变系统的k值，令k=2000绘制此时的波特图，matlab文本命令为；s=tf(‘s’);G=2000*(s+1)/(s*(s+2)*(s^2+17*s+4000))Bode(G)margin（G）grid on得到系统的波特图为：由波特图可以看出，当k值变大后，对相频曲线没有影响，因为k环节不提供相角，而对于幅频曲线来说当k值变为2000后相当于整个曲线向上平移了20lg2，从而使得幅值裕度和相角裕度改变了，幅值裕度为Gm=30.7dB，相角裕度为Pm=97度，剪切频率wc=0.256rad/s.B,设单位负反馈的开环传递函数为G（s）=10/(s+5)/(s-1)（1）绘制系统的Nyquist 曲线（2）分析系统的稳定性（3）根据系统的闭环零极点的分布图来分析系统的稳定性，和（2）得到的结果比较；1,绘制Nyquist 曲线的matlab文本命令为：num=10;den=conv([1 5],[1 -1]);nyquist(num,den)绘制出的图形为：2，分析系统的稳定性，当w趋于零时G(Jw)等于-2所以曲线的起点在（-2，0j），由曲线我们可以看出，曲线在（-1，0j）左边有半次自上而下的正穿越所以N+=0.5，N=2(N+-N-)=1，所以R=1，由系统的传递函数可以知道P=1，所以Z=P-R=0，从而得出系统是稳定的。

北京XX大学实验报告课程（项目）名称：线性控制系统的频域分析学院：专业：姓名：学号：指导教师：成绩：2013年12 月12 日实验三 线性控制系统的频域分析3. 1 频率特性测试一．实验目的1．了解线性系统频率特性的基本概念。

2．了解和掌握对数幅频曲线和相频曲线（波德图）的构造及绘制方法。

二．实验内容及步骤被测系统是一阶惯性的模拟电路图见图3-1，观测被测系统的幅频特性和相频特性，填入实验报告，並在对数座标纸上画出幅频特性和相频特性曲线。

本实验将正弦波发生器（B5）单元的正弦波加于被测系统的输入端，用虚拟示波器观测被测系统的幅频特性和相频特性，了解各种正弦波输入频率的被测系统的幅频特性和相频特性。

图3-1 被测系统的模拟电路图实验步骤：（1）将函数发生器（B5）单元的正弦波输出作为系统输入。

① 在显示与功能选择（D1）单元中，通过波形选择按键选中‘正弦波’（正弦波指示灯亮）。

② 量程选择开关S2置下档，调节“设定电位器2”，使之正弦波频率为8Hz （D1单元右显示）。

③ 调节B5单元的“正弦波调幅”电位器，使之正弦波振幅值输出为2V 左右（D1单元左显示）。

（2）构造模拟电路：按图3-1安置短路套及测孔联线，表如下。

（a ）安置短路套 （b ）测孔联线（3）运行、观察、记录：① 运行LABACT 程序，在界面的自动控制菜单下的线性控制系统的频率响应分析实验项目，选择时域分析，就会弹出虚拟示波器的界面，点击开始，用示波器观察波形，应避免系统进入非线性状态。

②点击停止键后，可拖动时间量程（在运行过程中，时间量程无法改变），以满足观察要求。

示波器的截图详见虚拟示波器的使用。

三．实验报告要求：按下表改变实验被测系统正弦波输入频率：（输入振幅为2V）。

实验截图：频率为1Hz时：频率为1.6Hz时：频率为3.2Hz时：频率为4.5Hz时：频率为8Hz时：3.2 一阶惯性环节的频率特性曲线一．实验目的1．了解和掌握一阶惯性环节的对数幅频特性)(ωL 和相频特性)(ωϕ，实频特性)Re(ω和虚频特性)Im(ω的计算。

线性系统的频域分析实验心得
1·熟练掌握用 MATLA语句绘制频域曲线。

2·掌握控制系统频域范围内的分析校正方法。

3掌握用频率特性法进行串联校正设计的思路和步骤
某单位负反馈控制系统的开环传递函数4为，试设计一超前校正装置,G(s)1、' s(s 1)K. 20s 150使校正后系统的静态速度误差系数，相位裕量，增益裕量20lgK10dB
绘制伯德图程序，以及计算穿越频率，相位裕量ans =相位 Inf 9.0406频率Inf 3.1425>e=5; r=50; rO=9; >>[gm1,pm1,wcg1,wcp1]=marg in(num 0,de nO);phic=(r-rO+e)*pi/180;
[gm1,pm1,wcg1,wcp1]=margi n(num 0,de nO);>>alpha=(1+s in (phic))/(1-si n(phic))[gm1,pm1,wcg1,wcp1]=margin(num 0,de n0); alpha =6.1261 [gm1,pm1,wcg1,wcp1]=marg in(num 0,de n0);lgm1,pm1,wcg1,wcp1]
通过MATLAB寸系统进行校正，可以清晰明了的显示矫正过程，以及矫正结果，方便快捷。

这种基于MATLAB的方法对于系统的设计非常实用。

值得以后再学习过程中认真领悟学习!! ! ! !。

实验1 控制系统典型环节的模拟实验（一）实验目的：1．掌握控制系统中各典型环节的电路模拟及其参数的测定方法。

2．测量典型环节的阶跃响应曲线，了解参数变化对环节输出性能的影响。

实验原理：控制系统模拟实验采用复合网络法来模拟各种典型环节，即利用运算放大器不同的输入网络和反馈网络模拟各种典型环节，然后按照给定系统的结构图将这些模拟环节连接起来，便得到了相应的模拟系统。

再将输入信号加到模拟系统的输入端，并利用计算机等测量仪器，测量系统的输出，便可得到系统的动态响应曲线及性能指标。

实验内容及步骤实验内容：观测比例、惯性和积分环节的阶跃响应曲线。

实验步骤：分别按比例，惯性和积分实验电路原理图连线，完成相关参数设置，运行。

①按各典型环节的模拟电路图将线接好(先接比例)。

(PID先不接)②将模拟电路输入端(U i)与阶跃信号的输出端Y相连接；模拟电路的输出端(Uo)接至示波器。

③按下按钮(或松开按钮)SP时，用示波器观测输出端的实际响应曲线Uo(t)，且将结果记下。

改变比例参数，重新观测结果。

④同理得积分和惯性环节的实际响应曲线，它们的理想曲线和实际响应曲线。

实验数据实验二控制系统典型环节的模拟实验（二）实验目的1．掌握控制系统中各典型环节的电路模拟及其参数的测定方法。

2．测量典型环节的阶跃响应曲线，了解参数变化对环节输出性能的影响。

实验仪器1．自动控制系统实验箱一台2．计算机一台实验原理控制系统模拟实验采用复合网络法来模拟各种典型环节，即利用运算放大器不同的输入网络和反馈网络模拟各种典型环节，然后按照给定系统的结构图将这些模拟环节连接起来，便得到了相应的模拟系统。

再将输入信号加到模拟系统的输入端，并利用计算机等测量仪器，测量系统的输出，便可得到系统的动态响应曲线及性能指标。

实验内容及步骤内容：观测PI，PD和PID环节的阶跃响应曲线。

步骤：分别按PI，PD和PID实验电路原理图连线，完成相关参数设置，运行①按各典型环节的模拟电路图将线接好。

1γ=50 20-=sK0原系统的伯德图：num/den =1.2347 s + 1 ------------- 0.20154 s + 1校正之后的系统开环传递函数为:num/den =6.1734 s + 5 ------------------------------------------- 0.20154 s^4 + 1.6046 s^3 + 3.4031 s^2 + 2 sP h a s e (d e g )Bode DiagramGm = Inf dB (at Inf rad/sec) , P m = 9.04 deg (at 3.14 rad/sec)-20020406080M a g n i t u d e (d B )alpha =6.1261;[il,ii]=min(abs(mag1-1/sqrt(alpha))); wc=w( ii); T=1/(wc*sqrt(alpha)); numc=[alpha*T,1]; denc=[T,1];[num,den]=series(num0,den0,numc,denc);[gm,pm,wcg,wcp]=margin(num,den); printsys(numc,denc)disp('УÕýÖ®ºóµÄϵͳ¿ª»·´«µÝº¯ÊýΪ:');printsys(num,den) [mag2,phase2]=bode(numc,denc,w); [mag,phase]=bode(num,den,w); subplot(2,1,1);semilogx(w,20*log10(mag),w,20*log10(mag1),'--',w,20*log10(mag2),'-.');grid; ylabel('·ùÖµ(db)'); title('--Go,-Gc,GoGc'); subplot(2,1,2);semilogx(w,phase,w,phase1,'--',w,phase2,'-',w,(w-180-w),':'); grid; ylabel('Ïàλ(0)'); xlabel('ƵÂÊ(rad/sec)');title(['УÕýÇ°£º·ùÖµÔ£Á¿=',num2str(20*log10(gm1)),'db','ÏàλԣÁ¿=',num2str(pm1),'0';'УÕýºó£º·ùÖµÔ£Á¿=',num2str(20*log10(gm)),'db','ÏàλԣÁ¿=',num2s tr(pm),'0']);10-110101102-60-40-2002040幅值(d b )--Go,-Gc,GoGc10-110101102-300-200-1000100相位(0)频率(rad/sec)矫正后系统的伯德图矫正之前系统单位阶跃响应矫正之后系统的单位阶跃响应：比较矫正前后系统的响应情况：可以看出超前矫正使系统的调节时间变短，响应更加迅速，但是超调量偏大，对改善系统的动态性能起到了巨大的作用。

控制系统的频域分析实验报告
摘要：
本实验旨在通过频域分析的方法来研究和评估控制系统的特性和性能。

在实验中，我们采用了频域分析的基本工具——Bode图和Nyquist图，通过对控制系统的幅频特性和相频特性进行分析，得出了系统的稳定性、干扰抑制能力和稳态性精度等方面的结论。

实验结果表明，频域分析是评估和优化控制系统的一种有效方法。

一、引言
频域分析是控制系统分析中常用的一种方法，通过对系统的频率响应进行研究，可以揭示系统的动态特性和性能，为控制系统的设计和优化提供指导。

在本实验中，我们将利用频域分析方法对一个具体的控制系统进行分析，通过实验验证频域分析的有效性。

二、实验装置和方法
实验所用控制系统包括一个控制对象（如电动机或水流系统）和一个控制器（如PID控制器）。

在实验中，我们将通过改变输入信号的频率来研究系统的频率响应。

实验步骤如下：
1. 连接实验装置，确保控制系统可正常工作。

2. 设计和设置适当的输入信号，包括常值信号、正弦信号和随
机信号等。

3. 改变输入信号的频率，记录系统的输出信号。

4. 利用实验记录的数据，绘制系统的幅频特性曲线和相频特性
曲线。

三、实验结果与讨论
根据实验记录的数据，我们绘制了控制系统的幅频特性曲线和
相频特性曲线，并对实验结果进行了分析和讨论。

1. 幅频特性分析
幅频特性曲线描述了控制系统对不同频率输入信号的增益特性。

在幅频特性曲线中，频率越高，输出信号的幅值越低，说明系统对
高频信号具有抑制作用。

本科实验报告课程名称：自动控制原理实验项目：控制系统的根轨迹和频域特性分析实验地点：多学科楼机房专业班级：学号：学生姓名：指导教师：2012 年5 月15 日一、实验目的和要求：1．学会利用MATLAB 绘制系统的根轨迹，并对系统进行分析； 2．学会利用MATLAB 对系统进行频域特性分析。

二、实验内容和原理：1．基于MATLAB 的控制系统根轨迹分析 1）利用MATLAB 绘制系统的根轨迹利用rlocus( )函数可绘制出当根轨迹增益k 由0至+∝变化时，闭环系统的特征根在s 平面变化的轨迹，该函数的调用格式为[r,k]=rlocus(num,den) 或 [r,k]=rlocus(num,den,k)其中，返回值r 为系统的闭环极点，k 为相应的增益。

rlocus( )函数既适用于连续系统，也适用于离散系统。

rlocus(num,den)绘制系统根轨迹时，增益k 是自动选取的，rlocus(num,den, k)可利用指定的增益k 来绘制系统的根轨迹。

在不带输出变量引用函数时，rolcus( )可在当前图形窗口中绘制出系统的根轨迹图。

当带有输出变量引用函数时，可得到根轨迹的位置列向量r 及相应的增益k 列向量，再利用plot(r,‘x’)可绘制出根轨迹。

2）利用MATLAB 获得系统的根轨迹增益 在系统分析中，常常希望确定根轨迹上某一点处的增益值k ，这时可利用MATLAB 中的rlocfind( )函数，在使用此函数前要首先得到系统的根轨迹，然后再执行如下命令[k,poles]=rlocfind(num,den) 或 [k,poles]=rlocfind(num,den,p)其中，num 和den 分别为系统开环传递函数的分子和分母多项式的系数按降幂排列构成的系数向量；poles 为所求系统的闭环极点；k 为相应的根轨迹增益；p 为系统给定的闭环极点。

例3-1 已知某反馈系统的开环传递函数为)2)(1()()(++=s s s ks H s G试绘制该系统根轨迹，并利用根轨迹分析系统稳定的k 值范围。

自动控制原理实验报告实验题目：实验三线性系统的频域分析班级：学号：姓名：指导老师：实验时间：实验三 线性系统的频域分析一、实验目的1．掌握用MATLAB 语句绘制各种频域曲线。

2．掌握控制系统的频域分析方法。

二、实验内容1．典型二阶系统2222)(nn ns s s G ωζωω++= 绘制出6=n ω，1.0=ζ，0.3，0.5，0.8，2的bode 图，记录并分析ζ对系统bode 图的影响。

程序：6=n ω，1.0=ζnum=[0 0 36];den=[1 1.2 36] w=logspace(-2,3,100); bode(num,den,w) grid3.0=ζnum=[0 0 36]; den=[1 3.6 36] w=logspace(-2,3,100); bode(num,den,w) grid5.0=num=[0 0 36];den=[1 6 36]w=logspace(-2,3,100); bode(num,den,w)grid8.0=num=[0 0 36];den=[1 9.6 36]w=logspace(-2,3,100);bode(num,den,w)grid=2num=[0 0 36];den=[1 24 36]w=logspace(-2,3,100);bode(num,den,w)grid分析ζ对系统bode图的影响：比较上述bode图知：根据其传递函数可知为振荡环节，其谐振频率，谐振峰值，系统的幅值和相位都与ζ有关。

系统的自然频率不变，当ζ逐渐增大时，其谐振频率逐渐减小。

根据bode图中的对数幅频特性图可知，在一定的范围内，随着ζ的增大，其谐振峰值逐渐减小。

ζ也影响其相频值，在一定的范围内，随着ζ的增大，其对数相频特性图曲线越平缓。

2．系统的开环传递函数为)5)(15(10)(2+-=s s s s G)106)(15()1(8)(22++++=s s s s s s G )11.0)(105.0)(102.0()13/(4)(++++=s s s s s s G绘制系统的Nyquist 曲线、Bode 图，说明系统的稳定性。

自动控制原理实验报告实验名称：线性系统的时域分析线性系统的频域分析线性系统的校正与状态反馈班级：学号：姓名：指导老师：2013 年12 月15日典型环节的模拟研究一. 实验目的1．了解和掌握各典型环节模拟电路的构成方法、传递函数表达式及输出时域函数表达式2．观察和分析各典型环节的阶跃响应曲线，了解各项电路参数对典型环节动态特性的影响二．实验内容及步骤观察和分析各典型环节的阶跃响应曲线，了解各项电路参数对典型环节动态特性的影响.。

改变被测环节的各项电路参数，画出模拟电路图，阶跃响应曲线，观测结果，填入实验报告运行LABACT 程序，选择自动控制菜单下的线性系统的时域分析下的典型环节的模拟研究中的相应实验项目，就会弹出虚拟示波器的界面，点击开始即可使用本实验机配套的虚拟示波器（B3）单元的CH1测孔测量波形。

具体用法参见用户手册中的示波器部分1)．观察比例环节的阶跃响应曲线典型比例环节模拟电路如图3-1-1所示。

图3-1-1 典型比例环节模拟电路传递函数：01(S)(S)(S)R R K KU U G i O === ； 单位阶跃响应： K )t (U = 实验步骤：注：‘S ST ’用短路套短接！（1）将函数发生器（B5）所产生的周期性矩形波信号（OUT ），作为系统的信号输入（Ui ）；该信号为零输出时，将自动对模拟电路锁零。

① 在显示与功能选择（D1）单元中，通过波形选择按键选中矩形波’（矩形波指示灯亮）。

② 量程选择开关S2置下档，调节“设定电位器1”，使之矩形波宽度＞1秒（D1单元左显示）。

③ 调节B5单元的“矩形波调幅”电位器使矩形波输出电压= 4V （D1单元‘右显示）。

（2）构造模拟电路：按图3-1-1安置短路套及测孔联线，表如下。

（a ）安置短路套 （b ）测孔联线（3）运行、观察、记录：打开虚拟示波器的界面，点击开始，按下信号发生器（B1）阶跃信号按钮（0→+4V 阶跃），观测A5B 输出端（Uo ）的实际响应曲线。

第五章 线性系统的频域分析例5-1 已知一控制系统结构图如图5-1所示，当输入r (t ) = 2sin t 时，测得输出c (t )=4sin(t -45︒)，试确定系统的参数ξ ，ωn 。

图5-1 系统结构图解：系统闭环传递函数为222()2nn ns s s ωΦξωω=++ 系统幅频特性为2()j Φω相频特性为222arctan)(ωωωξωωϕ--=nn由题设条件知c (t ) = 4sin( t -45︒) =2 A (1) sin(t + ϕ(1)) 即122222224)()1(=+-=ωωωξωωωnnnA 24)1(22222=+-=nnnωξωω1222arctan)1(=--=ωωωωξωϕn n ︒-=--=4512arctan2n nωξω整理得]4)1[(422224n n n ωξωω+-= 122-=n n ωξω解得 1.244n ω=，0.22ζ=例5-2 已知系统传递函数为)5.0)(2()52(10)(2-++-=s s s s s G ，试绘制系统的概略幅相特性曲线。

解：（1）传递函数按典型环节分解)15.0)(12()1)5(51251(50)(2+-++--=s s s s s G（2）计算起点和终点50)(lim 0-=→ωωj G ，10)(lim =∞→ωωj G相角变化范围：不稳定比例环节-50：-180︒ ~ -180︒；惯性环节1/(0.2s +1)：0︒~ -90︒；不稳定惯性环节1/(-2s +1)：0︒~ +90︒；不稳定二阶微分环节0.2s 2-0.4s +1：0︒~ -180︒（3）计算与实轴的交点22222)5.1()1()5.11)(25(10)(ωωωωωωω++-----=j j j G2222222)5.1()1()]5.35.5(3)1)(5([10ωωωωωωω+++-+++--=j（4） 确定变化趋势 根据G (j ω)的表达式，当ω <ωx 时，I m [G (j ω)] < 0；当ω >ωx 时，I m [G (j ω)] > 0。

实验二 频域分析一、实验目的1、熟练绘制系统的Nyquist 和Bode 图。

2、熟练掌握频域法分析系统。

二、 实验设备PC 机及MATLAB 平台三、实验原理及方法1、绘制Bode 图绘制bode 图的指令为bode(num,den)该指令表示在同一幅图中，分上、下两部分生成幅频特性和相频特性曲线。

虽未明确给出频率的取值范围，但Matlab 在频率响应范围内能自动选取频率值绘图。

若要具体给出频率的范围，可调用指令logspace(a,b,n)和bode(num,den,w)来绘制bode 图。

其中logspace(a,b,n)是产生频率响应自变量w 的采样点，即在十进制数10a 和10b 之间产生n 个十进制对数分度的等距离点，采样点n 的具体值由用户确定。

例：已知一单位反馈控制系统的开环传递函数为)15.0()11.0(10)(++=s s s s G ，试绘制该传递函数对应的bode 图。

解：先将传递函数改写为如下形式：s s s s den s num s G ++==25.010)()()(，然后使用bode(num,den)指令。

程序如下：num=[1 10];den=[0.5 1 0];w=logspace(-2,3,100);bode(num,den,w);grid on;title(‘bode diagram of G(s)=10(1+0.1s)/s(1+0.5s)’);若需要指定幅值和相角的取值范围，则需要调用如下的指令：[mag,phase,w]=bode(num,den,w)该指令等号左方的变量mag 和phase 试表示频率响应的幅值和相角，这些幅值和相角均由所选频率点的频率值计算后得出。

由于幅值的单位不是dB ，需增加一条指令：mag dB=20*lg10(mag)上述两条指令在应用时，还需加上如下两条指令，才能在屏幕上显示完整的bode 图。

Subplot(2 1 1) semilgx(w,20*lg10(mag));Subplot( 2 1 2) semilgx(w,phase); 例：已知系统的开环传递函数为)10016()12.0(30)(2+++=s s s s s G ，要求ω在3210~10-间作出该系统的bode 图。

系统频域分析实验报告1. 引言系统频域分析是一种用于研究线性时不变系统的方法，通过对系统的输入和输出信号在频域上的分析，可以得到系统的频率响应特性。

本实验旨在通过实际测量和分析，了解系统频域分析的基本原理和方法。

2. 实验设备和原理2.1 实验设备本实验所用设备包括： - 函数发生器 - 数字示波器 - 电阻、电容和电感等被测元件 - 电缆和连接线等连接配件2.2 实验原理系统频域分析是基于傅里叶变换的原理，通过将时域上的信号转换到频域上进行分析。

在本实验中，我们将使用函数发生器产生不同频率和幅度的正弦信号作为输入信号，通过被测系统输出的信号，使用数字示波器进行采集和分析。

3. 实验步骤3.1 连接实验设备将函数发生器的输出端与被测系统的输入端相连，将被测系统的输出端与数字示波器的输入端相连，确保连接正确可靠。

3.2 设置函数发生器调整函数发生器的频率、幅度和波形等参数，以产生不同频率和幅度的正弦信号作为输入信号。

3.3 采集数据使用数字示波器对被测系统的输出信号进行采集和记录。

可以选择适当的采样频率和采样时间，确保得到足够的数据点。

3.4 数据分析使用计算机软件或编程语言，对采集到的数据进行频域分析。

可以使用离散傅里叶变换（DFT）等方法，将时域上的信号转换到频域上，得到信号的频谱图。

3.5 分析结果根据得到的频谱图，可以分析出被测系统的频率响应特性。

可以通过找到频率响应曲线的极值点、截止频率等特征，来判断系统的性能和特点。

4. 实验结果和讨论4.1 频谱图展示根据采集到的数据和进行频域分析的结果，绘制出被测系统的频谱图。

4.2 频率响应特性分析根据频谱图的分析结果，可以得到被测系统的频率响应特性。

比如，可以观察到系统在不同频率下的增益特性、相位特性等。

4.3 讨论实验误差在实际实验中，可能存在各种误差的影响。

可以对实验误差进行分析和讨论，比如测量误差、系统本身的非线性特性等。

5. 结论通过本实验，我们了解了系统频域分析的基本原理和方法。

.专业整理 .实验三·线性系统的频域分析一、实验目的1．掌握用 MATLAB 语句绘制各种频域曲线。

2．掌握控制系统的频域分析方法。

二、实验内容1.典型二阶系统2G(s)ns2 2 n s2n绘制出n 6 ，0.1 ，0.3，0.5，0.8，2的bode图，记录并分析对系统bode 图的影响。

2.系统的开环传递函数为G( s)10s2(5s 1)(s5)G( s)8( s1)s2(s 15)( s26s10)G( s)4( s /31)s(0.02s 1)(0.05s 1)(0.1s1)绘制系统的 Nyquist 曲线、Bode 图和 Nichols 图，说明系统的稳定性，并通过绘制阶跃响应曲线验证。

s13.已知系统的开环传递函数为G(s)s2(0.1s 1)。

求系统的开环截止频率穿越频率、幅值裕度和相位裕度。

应用频率稳定判据判定系统的稳定性。

三、实验内容及分析.专业整理 .1. 系统 1： G (s)26，（1）0.1时n2中n 2s 2 n snMatlab 文本如下 ：num=[36 0 0];den=[11.236];w=logspace(-2,3,100); bode(num,den,w) Grid得到图像 ：同理 ，得到其他值情况下的波特图：ξ=0.3 时ξ=0.5 时ξ=0.8 时ξ=2 时从上面的图像中可以看出：随着的不断增大，波特图中震荡的部分变得越来越平滑。

而且，对幅频特性曲线来说，其上升的斜率越来越慢；对相频特性曲线来说，下降的幅度也在变缓。

2. 开环传递函数 1：G(s)10s2 (5s 1)(s 5)奈奎斯特图函数及图像如下：num=[0 10];den=[conv([5,-1],[1,5]),0,0];[z,p,k]=tf2zp(num,den); pnyquist(num,den)结果：p =0-5.00000.2000从上面的结果可知：在右半平面根的个数P=1 。

实验六 线性系统的频域分析一. 实验目的（1）熟练掌握使用MA TLAB 命令绘制控制系统Nyquist 图的方法； （2）能够分析控制系统Nyquist 图的基本规律；（3）加深理解控制系统乃奎斯特稳定性判据的实际应用； （4）学会利用奈氏图设计控制系统；（5）熟练掌握运用MA TLAB 命令绘制控制系统伯德图的方法； （6）了解系统伯德图的一般规律及其频域指标的获取方法； （7）熟练掌握运用伯德图分析控制系统稳定性的方法； （8）设计超前校正环节并绘制Bode 图； （9）设计滞后校正环节并绘制Bode 图。

二. 实验原理及内容 1、频率特性函数)(ωj G 。

频率特性函数为：nn n n m m m m a j a j a j a b j b j b j b jw G ++⋅⋅⋅++++⋅⋅⋅++=---)()()()()()()(1101110ωωωωωω由下面的MATLAB 语句可直接求出G(jw)。

i=sqrt(-1) % 求取-1的平方根GW=polyval(num ，i*w)./polyval(den ，i*w) 2、用MATLAB 作奈魁斯特图。

控制系统工具箱中提供了一个MATLAB 函数nyquist( )，该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。

当命令中不包含左端返回变量时，nyquist （）函数仅在屏幕上产生奈氏图，命令调用格式为：nyquist(num,den) ;作Nyquist 图，nyquist(num,den,w);作开环系统的奈氏曲线， 3、奈奎斯特稳定性判据（又称奈氏判据）反馈控制系统稳定的充分必要条件是当ω从-∞变到∞时，开环系统的奈氏曲线不穿过点（-1，j0）且逆时针包围临界点（-1，j0）点的圈数R 等于开环传递函数的正实部极点数。

4、用MATLAB 作伯德图控制系统工具箱里提供的bode()函数可以直接求取、绘制给定线性系统的伯德图。

实验三 线性系统的根轨迹一、实验目的1. 熟悉MATLAB 用于控制系统中的一些基本编程语句和格式。

2. 利用MATLAB 语句绘制系统的根轨迹。

3. 掌握用根轨迹分析系统性能的图解方法。

4. 掌握系统参数变化对特征根位置的影响。

二、基础知识及MATLAB 函数根轨迹是指系统的某一参数从零变到无穷大时，特征方程的根在s 平面上的变化轨迹。

这个参数一般选为开环系统的增益K 。

课本中介绍的手工绘制根轨迹的方法，只能绘制根轨迹草图。

而用MATLAB 可以方便地绘制精确的根轨迹图，并可观测参数变化对特征根位置的影响。

假设系统的对象模型可以表示为11210111()()m m m m n n n nb s b s b s b G s KG s K s a s b s a -+--++++==++++系统的闭环特征方程可以写成01()0KG s +=对每一个K 的取值，我们可以得到一组系统的闭环极点。

如果我们改变K 的数值，则可以得到一系列这样的极点集合。

若将这些K 的取值下得出的极点位置按照各个分支连接起来，则可以得到一些描述系统闭环位置的曲线，这些曲线又称为系统的根轨迹。

绘制系统的根轨迹rlocus （）MATLAB 中绘制根轨迹的函数调用格式为：rlocus(num,den) 开环增益k 的范围自动设定。

rlocus(num,d en,k) 开环增益k 的范围人工设定。

rlocus(p,z) 依据开环零极点绘制根轨迹。

r=rlocus(num,den) 不作图，返回闭环根矩阵。

[r,k]=rlocus(num,den) 不作图，返回闭环根矩阵r 和对应的开环增益向量k 。

其中，num,den 分别为系统开环传递函数的分子、分母多项式系数，按s 的降幂排列。

K 为根轨迹增益，可设定增益范围。

例3-1：已知系统的开环传递函数32(1)()429s G s K s s s *+=+++，绘制系统的根轨迹的matlab 的调用语句如下：num=[1 1]; %定义分子多项式 den=[1 4 2 9]; %定义分母多项式 rlocus (num,den) %绘制系统的根轨迹 grid %画网格标度线xlabel(‘Real Axis ’),ylabel(‘Imaginary Axis ’) %给坐标轴加上说明 title(‘Root Locus ’) %给图形加上标题名 若上例要绘制K 在（1，10）的根轨迹图，则此时的matlab 的调用格式如下。

线性系统的频率响应分析一、实验原理及内容频率特性当输入正弦信号时，线性系统的稳态响应具有随频率（ω由0变至∞）而变化的特性。

频率响应法的基本思想是：尽管控制系统的输入信号不是正弦函数，而是其它形式的周期函数或非周期函数，但是，实际上的周期信号，都能满足狄利克莱条件，可以用富氏级数展开为各种谐波分量；而非周期信号也可以使用富氏积分表示为连续的频谱函数。

因此，根据控制系统对正弦输入信号的响应，可推算出系统在任意周期信号或非周期信号作用下的运动情况。

线性系统的频率特性系统的正弦稳态响应具有和正弦输入信号的幅值比()Φ和相jw位差jwΦ随角频率（ω由0 变到∞）变化的特性。

而幅值比()Φjw 和相位差∠jwΦ的模和幅角。

所以只要把系统的传Φ恰好是函数jw递函数Φ（s）,令s =jw，即可得到Φ（jw）。

我们把Φ（jw）称为系统的频率特性或频率传递函数。

当w由0 到∞变化时，Φ（jw）随频率ω 的变化特性成为幅频特性，Φ（jw）随频率w的变化特性称为相频特性。

幅频特性和相频特性结合在一起时称为频率特性。

频率特性的表达式对数频率特性：又称波特图，它包括对数幅频和对数相频两条曲线，是频率响应法中广泛使用的一组曲线。

这两组曲线连同它们的坐标组成了对数坐标图。

对数频率特性图的优点：1.它把各串联环节幅值的乘除化为加减运算，从而简化了开环频率特性的计算与作图。

2.利用渐近直线来绘制近似的对数幅频特性曲线，而且对数相频特性曲线具有奇对称于转折频率点的性质，这些可使作图大为简化。

3.通过对数的表达式，可以在一张图上既能绘制出频率特性的中、高频率特性，又能清晰地画出其低频特性。

极坐标图（或称为奈奎斯特图）对数幅相图（或称为尼柯尔斯图）本次实验中，采用对数频率特性图来进行频域响应的分析研究。

实验中提供了两种实验测试方法：直接测量和间接测量。

直接频率特性的测量用来直接测量对象的输出频率特性，适用于时域响应曲线收敛的对象（如：惯性环节）。