基于Matlab控制系统频率特性分析法
实验三 系统频率特性曲线的绘制及系统分析
《自动控制原理》实践报告实验三系统频率特性曲线的绘制及系统分析熟悉利用计算机绘制系统伯德图、乃奎斯特曲线的方法，并利用所绘制图形分析系统性能。
一、实验目的1．熟练掌握使用MATLAB软件绘制Bode图及Nyquist曲线的方法；2．进一步加深对Bode图及Nyquist曲线的了解；3．利用所绘制Bode图及Nyquist曲线分析系统性能。
二、主要实验设备及仪器实验设备：每人一台计算机奔腾系列以上计算机，配置硬盘≥2G，内存≥64M。
实验软件：WINDOWS操作系统（WINDOWS XP 或WINDOWS 2000），并安装MATLAB 语言编程环境。
三、实验内容已知系统开环传递函数分别为如下形式， （1）)2)(5(50)(++=s s s G （2）)15)(5(250)(++=s s s s G（3）210()(21)s G s s s s +=++ （4）)12.0)(12(8)(++=s s s s G （5）23221()0.21s s G s s s s ++=+++ （6）)]105.0)125.0)[(12()15.0(4)(2++++=s s s s s s G 1．绘制其Nyquist 曲线和Bode 图，记录或拷贝所绘制系统的各种图形； 1、 程序代码： num=[50];den=conv([1 5],[1 2]); bode(num,den)num=[50];den=conv([1 5],[1 2]); nyquist(num,den)-80-60-40-20020M a g n i t u d e (d B)10-210-110101102103-180-135-90-450P h a s e (d e g )Bode DiagramFrequency (rad/sec)-1012345-4-3-2-11234Nyquist DiagramReal AxisI m a g i n a r y A x i s2、 程序代码： num=[250];den=conv(conv([1 0],[1 5]),[1 15]); bode(num,den)num=[250];den=conv(conv([1 0],[1 5]),[1 15]);-150-100-5050M a g n i t u d e (d B )10-110101102103-270-225-180-135-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)nyquist(num,den)3、 程序代码： num=[1 10];den=conv([1 0],[2 1 1]); bode(num,den)-150-100-50050100M a g n i t u d e (d B)10-210-110101102103-270-225-180-135-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)-1-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.10-15-10-551015System: sys Real: -0.132Imag: -0.0124Frequency (rad/sec): -10.3Nyquist DiagramReal AxisI m a g i n a r y A x i snum=[1 10];den=conv([1 0],[2 1 1]); nyquist(num,den)-25-20-15-10-5-200-150-100-5050100150200Nyquist DiagramReal AxisI m a g i n a r y A x i s-100-5050100M a g n i t u d e (d B )10-210-110101102-270-225-180-135-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)4、 程序代码： num=[8];den=conv(conv([1 0],[2 1]),[0.2 1]); bode(num,den)-18-16-14-12-10-8-6-4-20-250-200-150-100-50050100150200250Nyquist DiagramReal AxisI m a g i n a r y A x i snum=[8];den=conv(conv([1 0],[2 1]),[0.2 1]); nyquist(num,den)5、 程序代码： num=[1 2 1]; den=[1 0.2 1 1]; bode(num,den)num=[1 2 1];den=[1 0.2 1 1]; nyquist(num,den)-40-30-20-10010M a g n i t u d e (d B )10-210-110101102-360-270-180-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)-2.5-2-1.5-1-0.500.51 1.5-3-2-1123Nyquist DiagramReal AxisI m a g i n a r y A x i s-100-5050100M a g n i t u d e (d B )10-210-110101102-270-225-180-135-90P h a s e (d e g )Bode DiagramFrequency (rad/sec)6、 num=[2 4];den=conv(conv([1 0],[2 1]),[0.015625 0.05 1]); bode(num,den)num=[2 4];den=conv(conv([1 0],[2 1]),[0.015625 0.05 1]); nyquist(num,den)2．利用所绘制出的Nyquist 曲线及Bode 图对系统的性能进行分析：（1）利用以上任意一种方法绘制的图形判断系统的稳定性； 由Nyquist 曲线判断系统的稳定性，Z=P-2N 。
用MATLAB分析闭环系统的频率特性
用MATLAB分析闭环系统的频率特性闭环系统的频率特性指的是系统在不同频率下的响应特性。
在MATLAB中，可以通过不同的函数和工具箱来分析闭环系统的频率特性。
下面将介绍一些常用的方法。
1. 传递函数分析法（Transfer Function Analysis Method）：传递函数描述了系统的输入和输出之间的关系。
在MATLAB中，可以使用tf函数创建传递函数对象，并利用bode函数绘制系统的频率响应曲线。
例如，假设有一个传递函数G(s) = 1/(s^2 + s + 1)，可以用以下代码创建传递函数对象并绘制其频率响应曲线：```matlabG = tf([1], [1, 1, 1]);bode(G);```运行上述代码，将会显示出频率响应曲线，并且可以通过该函数的增益曲线和相位曲线来分析系统在不同频率下的响应特性。
2. 状态空间分析法（State-Space Analysis Method）：状态空间模型描述了系统的状态变量之间的关系。
在MATLAB中，可以使用ss函数创建状态空间模型，并利用bode函数绘制系统的频率响应曲线。
例如，假设有一个状态空间模型A、B、C和D分别为：```matlabA=[01;-1-1];B=[0;1];C=[10];D=0;sys = ss(A, B, C, D);bode(sys);```运行上述代码，将会显示出频率响应曲线，并且可以通过该函数的增益曲线和相位曲线来分析系统在不同频率下的响应特性。
3. 伯德图法（Bode Plot Method）：Bode图可以直观地表示系统的频率响应曲线。
在MATLAB中，可以使用bode函数绘制系统的Bode图。
例如，假设有一个传递函数G(s) =1/(s^2 + s + 1)，可以用以下代码绘制其Bode图：```matlabG = tf([1], [1, 1, 1]);bode(G);```运行上述代码，将会显示出Bode图，并且可以通过该图来分析系统在不同频率下的增益和相位特性。
实验2利用MATLAB分析信号频谱及系统的频率特性
实验2利用MATLAB分析信号频谱及系统的频率特性引言：在信号处理和通信领域中，频谱分析是一项非常重要的技术。
频谱分析可以帮助我们了解信号的频率特性，包括频率成分和幅度。
MATLAB是一款功能强大的数学软件，提供了多种工具和函数用于信号处理和频谱分析。
本实验旨在通过MATLAB分析信号频谱及系统的频率特性，深入理解信号处理和频域分析的原理和应用。
实验步骤：1.生成一个信号并绘制其时域波形。
首先，我们可以使用MATLAB提供的函数生成一个信号。
例如，我们可以生成一个用正弦函数表示的周期信号。
```matlabt=0:0.001:1;%时间范围为0到1秒，采样率为1000Hzf=10;%信号频率为10Hzx = sin(2*pi*f*t); % 生成正弦信号plot(t,x) % 绘制信号的时域波形图title('Time domain waveform') % 添加标题```2.计算信号的频谱并绘制频谱图。
使用MATLAB中的FFT函数可以计算信号的频谱。
FFT函数将信号从时域转换为频域。
```matlabFs=1000;%采样率为1000HzL = length(x); % 信号长度NFFT = 2^nextpow2(L); % FFT长度X = fft(x,NFFT)/L; % 计算X(k)f = Fs/2*linspace(0,1,NFFT/2+1); % 计算频率轴plot(f,2*abs(X(1:NFFT/2+1))) % 绘制频谱图title('Frequency spectrum') % 添加标题```3.使用MATLAB分析系统的频率特性。
MATLAB提供了Signal Processing Toolbox，其中包含了分析系统频率特性的函数和工具。
```matlabHd = designfilt('lowpassfir', 'FilterOrder', 6,'CutoffFrequency', 0.3, 'SampleRate', Fs); % 设计一个低通滤波器fvtool(Hd) % 显示滤波器的频率响应``````matlab[W,F] = freqz(Hd); % 计算滤波器的频率响应plot(F,abs(W)) % 绘制滤波器的振幅响应title('Frequency response of lowpass filter') % 添加标题```实验结果：运行上述代码后，我们可以得到如下结果：1.时域波形图2.频谱图3.滤波器频率响应讨论与结论：本实验通过MATLAB分析信号频谱及系统的频率特性，深入理解了信号处理和频域分析的原理和应用。
Matlab技术控制系统性能分析指南
MatIab技术控制系统性能分析指南概论当今社会，控制系统已成为各种领域中重要的技术和应用之一。
它们被广泛用于工业自动化、机电设备、航天航空等众多领域中。
控制系统的性能分析是确保系统正常运行和提高系统性能的必要步骤。
Mat1ab作为一种功能强大的工具，为控制系统性能分析提供了多种方法和技术。
本文将介绍一些基本的MaIIab技术，帮助读者进行控制系统性能分析。
一、系统建模在进行控制系统性能分析之前，首先需要进行系统建模。
系统建模是将实际物理系统抽象为数学模型的过程。
掌握系统建模方法对于准确分析系统性能至关重要。
Mat1ab提供了一系列工具和函数，可以用于快速建立系统模型。
有两种常用的系统建模方法：时域建模和频域建模。
1.时域建模时域建模基于系统的时间响应特性。
通过测量系统的输入和输出信号，并对其进行采样和离散化，可以得到系统的差分方程。
MaUab中的State-space函数是进行时域建模的常用工具。
它可以根据系统的状态方程和输出方程生成系统模型。
可以使用如下代码进行建模：A=∏2;341;B=[1;1];C=[10];D=O;sys=ss(A,B,C,D);其中，A、B、C和D分别表示状态空间方程的系数矩阵。
利用该函数建立的系统模型可以方便地进行时域性能分析。
2.频域建模频域建模基于系统的频率响应特性。
通过测量系统的输入和输出信号的频谱，并进行信号处理，可以得到系统的传递函数。
Mat1ab中的tf函数是进行频域建模的常用工具。
它可以根据系统的传递函数生成系统模型。
可以使用如下代码进行建模：num=[1];den=[11];sys=tf(num,den);其中，num和den分别表示传递函数的分子和分母系数。
利用该函数建立的系统模型可以方便地进行频域性能分析。
二、系统性能评估建立了系统模型之后，就可以进行系统性能的评估了。
针对不同的性能指标，可以使用不同的分析方法。
1稳态误差分析稳态误差衡量了系统在输入信号为稳态信号时的输出误差。
用MATLAB进行控制系统的动态性能的分析
用MATLAB进行控制系统的动态性能的分析MATLAB是一款功能强大的工具，可用于控制系统的动态性能分析。
本文将介绍使用MATLAB进行动态性能分析的常用方法和技巧，并提供实例来说明如何使用MATLAB来评估和改进控制系统的性能。
控制系统的动态性能是指系统对输入信号的响应速度、稳定性和精度。
评估控制系统的动态性能往往需要分析系统的阶跃响应、频率响应和稳态误差等指标。
一、阶跃响应分析在MATLAB中，可以使用step函数来绘制控制系统的阶跃响应曲线。
假设我们有一个系统的传递函数为：G(s)=(s+1)/(s^2+s+1)要绘制阶跃响应曲线，可以按照以下步骤操作：1.自动生成传递函数：num = [1 1];den = [1 1 1];G = tf(num,den);2.绘制阶跃响应曲线：step(G);二、频率响应分析频率响应分析用于研究控制系统对不同频率输入信号的响应特性。
在MATLAB中，可以使用bode函数来绘制控制系统的频率响应曲线。
假设我们有一个传递函数为：G(s)=1/(s+1)要绘制频率响应曲线，可以按照以下步骤操作：1.自动生成传递函数：num = [1];den = [1 1];G = tf(num,den);2.绘制频率响应曲线：bode(G);运行以上代码，MATLAB将生成一个包含系统幅频特性和相频特性的图形窗口。
通过观察频率响应曲线，可以评估系统的增益裕度（gain margin）和相位裕度（phase margin）等指标。
三、稳态误差分析稳态误差分析用于研究控制系统在稳态下对输入信号的误差。
在MATLAB中，可以使用step函数结合stepinfo函数来计算控制系统的稳态误差。
假设我们有一个传递函数为：G(s)=1/s要计算稳态误差，可以按照以下步骤操作：1.自动生成传递函数：num = [1];den = [1 0];G = tf(num,den);2.计算稳态误差：step(G);info = stepinfo(G);运行以上代码，MATLAB将生成一个阶跃响应曲线的图形窗口，并输出稳态误差等信息。
自动控制原理 matlab实验报告
自动控制原理实验（二）一、实验名称：基于MATLAB的控制系统频域及根轨迹分析二、实验目的：（1）、了解频率特性的测试原理及方法；（2）、理解如何用MATLAB对根轨迹和频率特性进行仿真和分析；（3）、掌握控制系统的根轨迹和频率特性两大分析和设计方法。
三、实验要求：（1）、观察给定传递函数的根轨迹图和频率特性曲线；（2）、分析同一传递函数形式，当K值不同时，系统闭环极点和单位阶跃响应的变化情况；（3）、K值的大小对系统的稳定性和稳态误差的影响；（4）、分析增加系统开环零点或极点对系统的根轨迹和性能的影响。
四、实验内容及步骤（1）、实验指导书：实验四（1）、“rlocus”命令来计算及绘制根轨迹。
会出根轨迹后，可以交互地使用“rlocfind”命令来确定点击鼠标所选择的根轨迹上任意点所对应的K值，K值所对应的所有闭环极点值也可以使用形如“[K, PCL] = rlocfind(G1)”命令来显示。
（2）、波特图：bode(G1, omga)另外，bode图还可以通过下列指令得出相位和裕角：[mag,phase,w] = bode(sys)（3）、奈奎斯特图：nuquist(G, omega)（2）课本：例4-1、4-2、4-7五实验报告要求（1）、实验指导书：实验四思考题请绘制下述传递函数的bode图和nyquist图。
1. 根据实验所测数据分别作出相应的幅频和相频特性曲线；2. 将思考题的解题过程(含源程序)写在实验报告中。
幅频特性曲线相频特性曲线Gs = zpk([10], [-5; -16; 9], 200)subplot(1, 2, 1)bode(Gs)gridsubplot(1, 2, 2)nyquist(Gs)grid（2）课本：例4-1、4-2、4-7图像结果：程序：Gs = zpk([-1], [0; -2; -3],1) rlocus(Gs)图像结果：程序：Gs = zpk([-2], [-1-j; -1+j],1) rlocus(Gs)程序：K=[0.5 1 2]for i=1:1:3num=[1,1,0,0]; den=[1,1,K(i)]; sys=tf(num,den); rlocus(sys); hold ongrid onend图像结果：目标：改变增益K和转折频率依次调节源程序：k1=[4.44,10,20];num=[1,2];den=conv([1,1],[1,2,4]);%一阶转折频率 1/T（wn1=2,wn2=1）二阶转折频率 wn3=wn'=2，伊布西塔=1/2 num1=[1,1];den1=conv([1,2],[1,2,4]);%一阶转折频率 1/T（wn1=1,wn2=2）二阶转折频率 wn3=wn'=2，伊布西塔=1/2 t=[0:0.1:7]; %for i=1:3g0=tf(k1(i)*num,den);g=feedback(g0,1);[y,x]=step(g,t);c(:,i)=y;g1=tf(k1(i)*num1,den1);g(1)=feedback(g1,1);[y1,x]=step(g(1),t);c1(:,i)=y1;endplot(t,c(:,1),'-',t,c(:,2),'-',t,c(:,3),'-',t,c1(:,1),'-',t,c1(:,2), '-',t,c1(:,3),'-');gridxlabel('Time/sec'),ylabel('out')结果分析：在本题中（1）改变k值：k值越大，超调量越大，调节时间越长，峰值时间越短，稳态误差越小（2）改变转折频率：超调量，调节时间，峰值时间，稳态误差同样有相应的变化。
基于Matlab的DFT及FFT频谱分析
基于Matlab的DFT及FFT频谱分析基于Matlab的DFT及FFT频谱分析一、引言频谱分析是信号处理中的重要任务之一，它可以揭示信号的频率特性和能量分布。
离散傅里叶变换（DFT）及快速傅里叶变换（FFT）是常用的频谱分析工具，广泛应用于许多领域。
本文将介绍通过Matlab进行DFT及FFT频谱分析的方法和步骤，并以实例详细说明。
二、DFT及FFT原理DFT是一种将时域信号转换为频域信号的离散变换方法。
它将信号分解成若干个正弦和余弦函数的叠加，得到频率和幅度信息。
FFT是一种高效的计算DFT的算法，它利用信号的对称性和周期性，将计算复杂度从O(N^2)降低到O(NlogN)。
FFT通过将信号分解成不同长度的子序列，递归地进行计算，最终得到频谱信息。
三、Matlab中的DFT及FFT函数在Matlab中，DFT及FFT可以通过内置函数进行计算。
其中，DFT使用函数fft，FFT使用函数fftshift。
fft函数可直接计算信号的频谱，fftshift函数对频谱进行频移操作，将低频移到频谱中心。
四、Matlab中DFT及FFT频谱分析步骤1. 读取信号数据首先，将待分析的信号数据读入到Matlab中。
可以使用内置函数load读取文本文件中的数据，或通过自定义函数生成模拟信号数据。
2. 时域分析通过plot函数将信号数据在时域进行绘制，以观察信号的波形。
可以设置合适的坐标轴范围和标签，使图像更加清晰。
3. 信号预处理针对不同的信号特点，可以进行预处理操作，例如去除直流分量、滤波等。
这些操作可提高信号的频谱分析效果。
4. 计算DFT/FFT使用fft函数计算信号数据的DFT/FFT，并得到频谱。
将信号数据作为输入参数，设置采样频率和点数，计算得到频谱数据。
5. 频域分析通过plot函数将频谱数据在频域进行绘制，观察信号的频率特性。
可以设置合适的坐标轴范围和标签，使图像更加清晰。
6. 结果解读根据频谱图像，分析信号的频率成分、幅度分布和峰值位置。
应用MATLAB绘制系统频率特性曲线
G(s)
稳定性。
s3
5系s22.7统4s的 2奈奎斯特曲线，并利用曲线来判别闭环系统的
解 MATLAB仿真程序代码如下：
num1=[ 2.7]；
den1=[1 5 4 2]；
sys1=tf(num1,den1)；
nyquist(sys1)
title('Nyquist图');
运行后，获得如图1-46所示曲线。
自动控制原理
应用MATLAB绘制系统频率特性曲线
1.1用MATLAB绘制系统开环对数频率特性 对于连续系统，用MATLAB函数绘制系统开环对 数频率特性的函数命令调用格式有 Bode（sys） Bode（sys，w） Bode（sys1，sys2,…，sysN） Bode（sys1，sys2,…，sysN，w） [mag，phase，w]=Bode(sys)
例1-15 绘制一阶惯性环节
的G(s奈) 奎3斯特图。
5s 1
解 MATLAB仿真程序代码如下：
G=tf(3，[5 1])；
nyquist(G)；Fra bibliotekhold on；
title('Nyquist图')；
运行后，获得如图1-45所示曲线。
图1-45 例1-15系统极坐标曲线图
例1-16 用函数nyquist(sys)绘制开环传递函数为
bode(num,den);hold on；
end
grid
获得振荡环节伯德图如图1-43所示，
图1-43 例1-13开环系统伯德图
如果希望求取控制系统的增益裕量 、相位GM裕量 界频率（也称PM交叉频
率） 、穿越g频率 ，可以使c 用margin函数 计算控制系统的相关稳定裕度值。
浅议用Matlab绘制Bode图及频率特性分析
频率特性的描述： 极坐标形式：
这两种描述之间的关系：
系统的频率特性图主要有 Ｂｏｄｅ（波德）图、Ｎｙｑｕｉｓｔ（奈 奎斯特）曲线、Ｎｉｃｈｏｌｓ（尼柯尔斯）图。
限于篇幅，本文以 Ｂｏｄｅ 图为例来说明。 ２ 用 Ｍａｔｌａｂ 来绘制 Ｂｏｄｅ 图
Ｂｏｄｅ图是指由对数幅频特性和相频特性组成的对数频率 特性图。对数幅频特性和相频特性曲线的横轴都是频率，且 均是按对数分度的。对数幅频特性的纵轴为：
通过设计，构思出程序流程框图，以最典型的二阶系统 为例来说明。系统的传பைடு நூலகம்函数为
绘制出 ｗ 取一系列不同值时的 Ｂｏｄｅ 图，如图一所示。
图一 Ｂｏｄｅ 图
２５１
科技广场 ２００７．１１
浅议用 Ｍａｔｌａｂ 绘制 Ｂｏｄｅ 图及频率特性分析
Ｄｒａｗｉｎｇ Ｂｏｄｅ Ｇｒａｐｈｉｃ ｗｉｔｈ Ｍａｔｌａｂ ａｎｄ Ａｎａｌｙｓｉｓ ｏｆ Ｆｒｅｑｕｅｎｃｙ
黄 涛 郭玉婷 Ｈｕａｎｇ Ｔａｏ Ｇｕｏ Ｙｕｔｉｎｇ （１．吉林师范大学信息技术学院， 吉林 四平 １３６０００； ２．四平市第三高级中学， 吉林 四平 １３０００１） （１．Ｃｏｌｌｅｇｅ ｏｆ Ｉｎｆｏｒｍａｔｉｏｎ Ｔｅｃｈｎｏｌｏｇｙ， Ｊｉｌｉｎ Ｎｏｒｍａｌ Ｕｎｉｖｅｒｓｉｔｙ， Ｊｉｌｉｎ Ｓｈｐｉｎｇ １３６０００； ２．Ｔｈｅ Ｔｈｉｒｄ Ｈｉｇｈ Ｓｃｈｏｏｌ， Ｊｉｌｉｎ Ｓｉｐｉｎｇ １３６００１）
分
版社．
析
［３］ＭＡＴＬＡＢ电子仿真与应用［Ｍ］．北京：国防工业出版社．
［４］电子系统仿真与ＭＡＴＬＡＢ［Ｍ］．北京：北京广播学院出版
社．
作者简介 黄涛，男，吉林师范大学信息技术学院教师，讲师，吉
基于MATLAB的控制系统稳定性分析报告
四川师范大学本科毕业设计基于MATLAB的控制系统稳定性分析学生姓名宋宇院系名称工学院专业名称电气工程及其自动化班级 2010 级 1 班学号**********指导教师杨楠完成时间2014年 5月 12日基于MATLAB的控制系统稳定性分析电气工程及其自动化本科生宋宇指导老师杨楠摘要系统是指具有某些特定功能，相互联系、相互作用的元素的集合。
一般来说，稳定性是系统的重要性能，也是系统能够正常运行的首要条件。
如果系统是不稳定，它可以使电机不工作，汽车失去控制等等。
因此，只有稳定的系统，才有价值分析与研究系统的自动控制的其它问题。
为了加深对稳定性方面的研究，本设计运用了MATLAB软件采用时域、频域与根轨迹的方法对系统稳定性的判定和分析。
关键词：系统稳定性 MATLAB MATLAB稳定性分析ABSTRACT System is to point to have certain function, connect with each other, a collection of interacting elements. Generally speaking, the stability is an important performance of system, also is the first condition of system can run normally. If the system is not stable, it could lead to motor cannot work normally, the car run out of control, and so on. Only the stability of the system, therefore, have a value analysis and the research system of the automatic control of other problems. In order to deepen the study of stability, this design USES the MATLAB software using the time domain, frequency domain and the root locus method determination and analysis of the system stability.Keywords: system stability MATLAB MATLAB stability analysis目录摘要 (I)ABSTRACT .......................................................... I I 目录1.绪论 (1)1.1自动控制理论发展概述 (1)1.1.1经典控制理论的发展及其基本内容 (1)1.1.2现代控制理论的发展及其基本内容 (1)1.1.3智能控制理论的发展及其主要内容 (2)1.2本文的章节安排 (2)2控制系统的理论基础 (3)2.1控制系统的基本形式 (3)2.1.1闭环控制系统 (3)2.1.2开环控制系统 (4)2.1.3小结 (4)2.2控制系统的分类 (4)2.3控制系统的稳定性 (5)3 MATLAB基础介绍 (6)3.1MALTAB概述 (6)3.2MATLAB的特点 (6)4稳定性分析的方法介绍 (7)4.1时域分析法 (7)4.1.1时域分析法的概念 (7)4.1.2控制系统的性能指标 (7)4.1.3典型的输入信号 (7)4.1.4系统时域分析函数-Step函数 (8)4.1.5控制系统的时域分析-impulse函数 (10)5根轨迹分析法 (12)5.1根轨迹分析法的概念 (12)5.1.1一般控制系统 (12)5.2绘制控制系统的根轨迹图的一般规则 (12)5.3pzmap函数 (13)5.4rlocus函数 (14)6频域法分析 (16)6.2奈氏图（Nyquist） (16)6.3波德图（Bode） (18)7总结 (22)参考文献 (23)致谢 (24)基于MATLAB的控制系统稳定性分析1.绪论这章讲述了自动控制理论与控制技术概述，主要介绍了几种自动控制理论的发展概况以及基本的内容。
实验七 用MATLAB进行系统频率特性曲线绘制
自动控制原理课程验证性实验报告
实验名称
用MATLAB进行系统频率特性曲线绘制
实验时间
2013年05月30日
学生姓名
实验地点
实验数据分析时可围绕要达成的实验目标，借助公式、模型或必要软件，结合作图、计算等方法对实验数据进行数据处理后，得出实验结果。实验数据分析过程应条理清晰，层次清楚。
必要时，可对实验现象和实验结果综合分析，相互印证。
7、实验结论：
通过对实验现象和实验数据分析后得_close),grid on
g.G=tf(100*conv([1,5],[1,5]),conv([1,1],[1,1,9]));
[Gm,Pm,Wcg,Wcp]=margin(G) G_close=feedback(G,1);step(G_close),grid on
5、实验现象、实验数据记录：
原理：（1）幅相频率特性曲线，以角频率 为参变量，当 从 变化时，频率特性构成的向量在复平面上描绘的曲线称为幅相频率特性曲线。（2）奈氏判据，奈氏判据是利用系统开环频率特性来判断闭环系统稳定性的一个判据，便于研究当系统结构改变时对系统稳定性的影响。其内容是：反馈控制系统稳定的充分必要条件是当 从 变到 时，开环系统的奈氏判据 不穿过 点逆时针包围临界点 的圈数 等于开环传递函数的正实部极点数 。1）对于开环稳定的系统，闭环系统稳定的充分必要条件是：开环系统的奈氏曲线 不包围 点。反之，则闭环系统是不稳定的。2）对于开环不稳定的系统，有 个开环极点位于右半 平面，则闭环系统稳定的充分必要条件是：当 从 变到 时，开环系统的奈氏判据 逆时针包围 点 次。
《自动控制原理》Matlab求解控制系统频域分析实验
《自动控制原理》Matlab求解控制系统频域分析实验
一、实验目的
1、加深了解系统频率特性的概念。
2、学习使用Matlab软件绘制Nyquist图、Bode图的基本方法。
3、掌握典型环节的频率特性。
二、实验仪器
Matlab2014b版
三、实验原理
1.奈奎斯特图（幅相频率特性图）
MATLAB为用户提供了专门用于绘制奈奎斯特图的函数nyquist
axis([-2,0.4,-1.5,1.5]);
k=500;
num=[1,10];
den=conv([1,0],conv([1,1],conv([1,20],[h,50])))；
w=logspace(-1,3,200)
bode(k*num,den,w);
grid;
五、实验原始数据记录与数据处理
六、实验结果与分析讨论
范围是自动确定的。当需要指定幅值范围和相角范围时，则需用下面的功能指令：
[mag,phase,w]=bode(num,den,w)
四、实验内容及步骤
z=[]:
p=[0,-1,-2]:
k=5;
g=zpk(z,p,k):
nyquist(g);
w=0.5：0.1:10:
figure(2):
nyquist(g:w);
基于MATLAB自动控制系统时域频域分析与仿真
基于MATLAB自动控制系统时域频域分析与仿真MATLAB是一款强大的数学软件，也是自动控制系统设计的常用工具。
它不仅可以进行时域分析和频域分析，还可以进行相关仿真实验。
本文将详细介绍MATLAB如何进行自动控制系统的时域和频域分析，以及如何进行仿真实验。
一、时域分析时域分析是指对系统的输入信号和输出信号进行时域上的观察和分析，以了解系统的动态特性和稳定性。
MATLAB提供了一系列的时域分析工具，如时域响应分析、稳态分析和步骤响应分析等。
1.时域响应分析通过时域响应分析，可以观察系统对于不同的输入信号的响应情况。
在MATLAB中，可以使用`lsim`函数进行系统的时域仿真。
具体步骤如下：- 利用`tf`函数或`ss`函数创建系统模型。
-定义输入信号。
- 使用`lsim`函数进行时域仿真，并绘制系统输出信号。
例如，假设我们有一个二阶传递函数模型，并且输入信号为一个单位阶跃函数，可以通过以下代码进行时域仿真：```num = [1];den = [1, 1, 1];sys = tf(num, den);t=0:0.1:10;u = ones(size(t));[y, t, x] = lsim(sys, u, t);plot(t, y)```上述代码中，`num`和`den`分别表示系统的分子和分母多项式系数，`sys`表示系统模型，`t`表示时间序列，`u`表示输入信号，`y`表示输出信号。
通过绘制输出信号与时间的关系，可以观察到系统的响应情况。
2.稳态分析稳态分析用于研究系统在稳态下的性能指标，如稳态误差和稳态标准差。
在MATLAB中，可以使用`step`函数进行稳态分析。
具体步骤如下：- 利用`tf`函数或`ss`函数创建系统模型。
- 使用`step`函数进行稳态分析，并绘制系统的阶跃响应曲线。
例如，假设我们有一个一阶传递函数模型，可以通过以下代码进行稳态分析：```num = [1];den = [1, 1];sys = tf(num, den);step(sys)```通过绘制系统的阶跃响应曲线，我们可以观察到系统的稳态特性。
MATLAB与控制系统的频域分析
、
对数 频率特 性 图 （ 德 图 ） 伯
Ｂ ｄ （ｕ ｄ ｎｗ） ｏ ｅｎ ｍ，ｅ， 函数 可利用 指 定 的角 频率 矢量绘 制 出系 统 的伯德 图。一 典型 二阶系 统 开环传 递 函数
Ｇ 南 （
要改参 制统雌 并析 对绕 影。ＴＢ 求变数 系 分参 ／ 响 Ａ 系 ＭＬ Ａ
应用。
关键 词：ＭＡ Ｌ Ｔ ＡＢ；伯 德 图；奈 奎斯特 图；频率 特性
中图分类号 ：Ｔ ２ ４ Ｐ ７ 文献标 识码 ：Ａ 文章 编号 ： １ ０－ ６ ２ （ ０ ）Ｏ ０ ５ ． ４ ０ ８－ １ ９ ２ １ １ ｌ ８－ ０
工程 上 ，控 制系 统性 能分 析 目前 广泛应 用 的有 两种 方法 ：一种 是频 率特 性法 ：另一种 是 根轨 迹法 。 这些方 法 的使 用 需要绘 制 时间 响应 曲线 、 轨迹 图及频 率特 性 中 的 Ｂ ｄ 、 ｑ ｉ 图和 Ｎｉ ｏｓ曲线 等 。 根 ｏ ｅ图 Ｎｙ ｕｓ ｔ ｃ ｌ ｈ 这些 曲线 的绘制 过去 通 常是 由手 工完 成 的 ，其 误差较 大 、效率 低 ，对 于 某 些较 复杂 的高 阶系 统 ，很难 用
特 性 的幅值和 角度值 ，当该函 数不 带返 回参数 时 ，将直 接在 屏幕 上绘 制 出系统 的伯 德 图 。取 ｗｎ１，，取 ＝０ 亏 ［．：０６ ］ Ｏ ２ ．：１时的二 阶系 统 的伯 德 图可直 接 由 ｂ ｄ ０＂ 得 ￣ ． 行 下面 的 Ｍ 文件 ，运行 后得 到 的伯德 ｏｅ ｉ ￣数 Ｊ ｌ ｏ执
ｆ ｒｉ＝ｋｏ ａ ｏ ｉ ｓｉ
ｓｂ ｌｔ ，，） ｕ ｐｏ（ １ ２ １
ｇｉ ｎ ｒｄｏ
ｔｌ ‘ ｏ ｅ ａｒｍ’ ｉｅ Ｂ ｄ ｇａ ） ｔ（ Ｄｉ ； ｘａ ｌ ‘ｒｑ ｅ ｃ （ａ ／ｅ ） ； ｌ ｅ Ｆｅ ｕ ｎｙ ｒｄ ｓｃ’ ｂ（ ） ｙａ ｌ ‘ ｇ ｉ ｄ （Ｂ ’ ｌ ｅ Ｍａ ｎｔ ｅ ） ； ｂ（ ｕ ｄ ）
MATLAB进行控制系统频域分析
一、基于MATLAB 的线性系统的频域分析基本知识（１）频率特性函数)(ωj G 。
设线性系统传递函数为：nn n n m m m m a s a s a s a b s b s b s b s G ++⋅⋅⋅++++⋅⋅⋅++=---1101110)( 则频率特性函数为：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 ）其中（num ，den ）为系统的传递函数模型。
而w 为频率点构成的向量，点右除（./)运算符表示操作元素点对点的运算.从数值运算的角度来看，上述算法在系统的极点附近精度不会很理想，甚至出现无穷大值，运算结果是一系列复数返回到变量GW 中。
（２）用MATLAB 作奈魁斯特图。
控制系统工具箱中提供了一个MATLAB 函数nyquist( ），该函数可以用来直接求解Nyquist 阵列或绘制奈氏图。
当命令中不包含左端返回变量时，nyquist （）函数仅在屏幕上产生奈氏图，命令调用格式为:nyquist(num ，den) nyquist （num,den ，w) 或者nyquist(G) nyquist(G,w ） 该命令将画出下列开环系统传递函数的奈氏曲线： )()()(s den s num s G = 如果用户给出频率向量w ，则w 包含了要分析的以弧度/秒表示的诸频率点。
在这些频率点上,将对系统的频率响应进行计算，若没有指定的w 向量，则该函数自动选择频率向量进行计算。
w 包含了用户要分析的以弧度/秒表示的诸频率点，MATLAB 会自动计算这些点的频率响应。
当命令中包含了左端的返回变量时，即：[re,im ，w]=nyquist （G ）或[re ，im,w ］=nyquist （G ，w ） 函数运行后不在屏幕上产生图形，而是将计算结果返回到矩阵re 、im 和w 中。
控制系统的频率特性分析
【实验名称】控制系统的频率特性分析【实验目的】1) 掌握运用MATLAB 软件绘制控制系统波特图的方法； 2) 掌握MATLAB 软件绘制奈奎斯特图的方法； 3) 利用波特图和奈奎斯特图对控制系统性能进行分析。
【实验仪器】1) PC 机一台 2) MATLAB 软件【实验原理】1. 奈奎斯特稳定判据及稳定裕量（1）奈氏（Nyquist ）判据：反馈控制系统稳定的充要条件是奈氏曲线逆时针包围临界点的圈数R 等于开环传递函数右半s 平面的极点数P , 即R=P ;否则闭环系统不稳定, 闭环正实部特征根个数Z 可按下式确定Z=P-R=P-2N （2）稳定裕量利用)()(ωωj H j G 轨迹上两个特殊点的位置来度量相角裕度和增益裕度。
其中)()(ωωj H j G 与单位圆的交点处的频率为c ω（截止频率）；)()(ωωj H j G 与负实轴的交点频率为x ω（穿越频率）。
则相角裕度：)(180)()(180c c c j H j G ωϕωωγ+=∠+= 增益裕度：)(1)()(1x x x A j H j G h ωωω==（对数形式：)(lg 20)()(lg 20x x x A j H j G h ωωω-=-= 2. 对数频率稳定判据将系统开环频率特性曲线分为幅频特性和相频特性，分别画在两个坐标上，横轴都用频率ω，纵轴一个用对数幅值和相角，这两条曲线画成的图就是Bode 图，即对数频率特性图。
因为Bode 图与奈氏图有一一对应关系，因此，奈氏稳定判据就可描述为基于Bode 图的对数频率稳定判据：（1）开环系统稳定，即开环系统没有极点在正右半根平面，如果其对数幅频曲线大于0dB 的区域内，相频曲线对180-线正负穿越次数相等，那么闭环系统就是稳定的，否则是不稳定的。
（2）开环系统不稳定，有P 个极点在正右半平面，如果其对数幅频曲线大于0dB 的区域内，相频曲线对180-线正穿越次数大于负穿越次数P/2，闭环系统就是稳定的，否则是不稳定的。
自动控制原理的MATLAB仿真与实践第5章 线性系统的频域分析
7
【例5-1】 试绘制惯性环节G(jω)=1/(2s+1)的Nyquist曲线 和Bode图。
解：程序如下：
>>clear
G=tf(1,[2,1]); %建立模型
nyquist (G); %绘制Nyquist图
figure(2); bode (G); %绘制Bode图
4
ngrid；ngrid(‘new’)：绘制尼科尔斯坐标网格即等 20lgM圆和等角曲线组成的网格。‘new’代表清除以前 的图形，与后一个nichols()一起绘制网格。
semilogx(w,20*log10(mag))：绘制半对数坐标下的幅 频特性曲线。
semilogx(w,phase*180/pi)：绘制半对数坐标下的相频 特性曲线。
MATLAB提供了许多用于线性系统频率分析 的函数命令，可用于系统频域的响应曲线、参数分析 和系统设计等。常用的频率特性函数命令格式及其功 能见表5-1。 bode (G)：绘制传递函数的伯德图。其中：G为传递
函数模型，如：tf(), zpk(), ss()。 bode(num,den)：num,den分别为传递函数的分子与
本节分别介绍利用MATLAB进行频域绘图和频 率分析的基本方法。
6
5.2.1 Nyquist曲线和Bode图
MATLAB频率特性包括幅频特性和相频特性。 当用极坐标图描述系统的幅相频特性时，通常称为 奈奎斯特（Nyquist）曲线；用半对数坐标描述系 统的幅频特性和相频特性时，称为伯德（Bode） 图；在对数幅值-相角坐标系上绘制等闭环参数（ M和N）轨迹图，称为尼克尔斯（Nichols）图。
运行结果如图5-2所示。
实验二利用MATLAB进行系统动态特性分析(任务)
实验二利用MATLAB进行系统动态特性分析（任务）引言：系统动态特性分析是指通过研究系统的动态响应，来了解系统的性能和稳定性。
在工程领域中，对不同系统进行动态特性分析是非常重要的，可以帮助我们了解系统的稳定性、响应特性以及对外部输入的敏感度等，并且可以为系统设计和控制提供重要的依据。
实验目的：通过数据采集的方法，运用MATLAB工具对动态系统进行特性分析，掌握系统的稳态特性和暂态特性，并对系统性能进行评估。
实验器材和原理：实验器材：电脑、MATLAB软件实验步骤：1. 导入数据：将实验得到的数据导入MATLAB中，可以通过Excel等工具将数据保存为文本格式，然后使用MATLAB的读取函数导入数据。
2.绘制时域响应曲线：根据导入的数据，使用MATLAB中的绘图函数绘制出时域响应曲线。
根据实验需要，选择绘制的曲线类型，如步跃响应曲线、阶跃响应曲线等。
3.基本特性分析：-稳态误差：通过分析曲线的极限值和最终值，计算出系统的稳态误差。
-加载响应：通过观察曲线的上升时间、峰值时间、峰值以及超调量等指标，来评估系统的负载能力。
-过渡过程：观察曲线的上升时间、峰值时间以及超调量等指标，来评估系统的动态响应特性。
4.绘制频域响应曲线：通过数据采集得到的数据，使用MATLAB中的频域分析工具绘制频域响应曲线，观察系统的频域特性。
5.使用MATLAB进行数据处理和分析：根据实验需要，对导入的数据进行处理和分析，如计算系统的传递函数、计算系统的频域性能等。
6.实验结果分析：根据绘制的曲线和计算的数据，分析系统的稳态特性和暂态特性，并对系统的性能进行评估。
可以根据实验结果，进行系统设计改进或控制参数调整。
实验注意事项：1.数据采集过程中要注意信号的采样频率和采样精度，以保证数据的准确性。
2.在绘制曲线时要选择合适的曲线类型和参数，使得曲线能够准确表达系统的动态特性。
3.在数据处理和分析过程中要注意使用合适的算法和公式，确保结果的准确性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
基于Matlab控制系统频率特性分析法本文主要介绍了基于Matlab控制系统的频率特性分析方法、频域稳定性判据以及开环频域性能分析，并获得频率响应曲线等。
通过本章的学习，可以利用MATLAB对各种复杂控制系统进行频率分析，以此获得系统稳定性及其它性能指标。
一、频率特性基本概念如果将控制系统中的各个变量看成是一些信号，而这些信号又是由许多不同频率的正弦信号合成的，则各个变量的运动就是系统对各个不同频率信号响应的总和。
系统对正弦输入的稳态响应称频率响应。
利用这种思想研究控制系统稳定性和动态特性的方法即为频率响应法。
频率响应法的优点为：⑴物理意义明确；⑵可利用试验方法求出系统的数学模型，易于研究机理复杂或不明的系统，也适用于某些非线性系统；⑶采用作图方法，非常直观。
1. 频率特性函数的定义对于稳定的线性系统或者环节，在正弦输入的作用下，其输出的稳态分量是与输入信号相同频率的正弦函数。
输出稳态分量与输入正弦信号的复数比，称为该系统或环节的频率特性函数，简称为频率特性，记作G（jω）=Y(jω)/R(jω)对于不稳定系统，上述定义可以作如下推广。
在正弦输入信号的作用下，系统输出响应中与输入信号同频率的正弦函数分量和输入正弦信号的复数比，称为该系统或环节的频率特性函数。
当输入信号和输出信号为非周期函数时，则有如下定义。
系统或者环节的频率特性函数，是其输出信号的傅里叶变换像函数与输入信号的傅里叶变换像函数之比。
2. 频率特性函数的表示方法系统的频率特性函数可以由微分方程的傅里叶变换求得，也可以由传递函数求得。
这三种形式都是系统数学模型的输入输出模式。
当传递函数G(s)的复数自变量s沿复平面的虚轴变化时，就得到频率特性函数G(jω)=G(s)|s=jω。
所以频率特性是传递函数的特殊形式。
代数式：G(jω)=R(w)+jI(ω)R(w)和I(w)称为频率特性函数G(jw)的实频特性和虚频特性。
指数式：G(jω)=A(w)eΦ(ω)式中A(ω)=| G(jω)|是频率特性函数G(jw)的模，称为幅频特性函数。
Φ(w)=arg G(jω)是频率特性函数G(jω)的幅角，称为相频特性函数。
3．频率响应曲线系统的频率响应可以用复数形式表示为G(jω)，常用的频率响应表示方法是图形表示法。
根据系统频率响应幅值、相位和频率之间的不同显示形式，有伯德（Bode）图、奈魁斯特（Nyquist）图和尼柯尔斯（Nichols）图。
3.1 伯德图伯德（Bode）图又称对数频率特性图，由对数幅频特性图和相频特性图组成。
伯德图的横坐标为角频率ω，按常对数lgω分度。
对数幅频特性的纵坐标是对数幅值。
L(ω)=20lg A(ω)，单位为分贝（dB），线性分度。
对数相频特性的纵坐标为φ(ω)，单位为度，线性分度。
一般情况下，控制系统开环对数频率特性图的绘制步骤如下：⑴将开环频率特性按典型环节分解，并写成时间常数形式；⑵求出各转角频率（交接频率），将其从小到大排列为ω1，ω2，ω3，……，并标注在ω轴上；⑶绘制低频渐近线（ω1左边的部分），这是一条斜率为-20νdB/decade（ν为系统开环频率特性所含1/jw因子的个数）的直线，它或者它的延长线应通过点（1,20lgK）；⑷各转角频率间的渐近线都是直线，但自最小的转角频率ω1起，渐近线斜率发生变化，斜率变化取决于各转角频率对应的典型环节的频率特性函数。
例1 绘制一阶惯性环节G(s)=1/(4s+1)的伯德图。
程序代码如下：num=1;den=[4 1];G=tf(num,den);bode(G,'r')3.2奈魁斯特图奈魁斯特图又称为极坐标图或者幅相频率特性图。
频率特性函数G(jω)的奈魁斯特图是角频率ω由0变化到∞时，频率特性函数在复平面上的图像。
它以ω为参变量，以复平面上的向量表示G(jω)的一种方法。
G(jω)曲线的每一点都表示与特定ω值相应的向量端点，向量的幅值为|G(jω)|，相角为argG(jω)；向量在实轴和虚轴上的投影分别为实频特性R(ω)和虚频特性I(ω)。
一般情况下，系统开环频率特性函数奈魁斯特图的绘制步骤如下：⑴将系统的开环频率特性函数G0(jω)写成G(jω)=A(w)eΦ(ω)；⑵确定奈魁斯特图的起点（ω＝0+）和（ω→+∞）。
起点与系统所包含的积分环节个数（ ）有关，终点的A(ω)与系统开环传递函数分母和分子多项式阶次的差有关；⑶确定奈魁斯特图与坐标轴的交点；⑷根据以上的分析并且结合开环频率特性的变化趋势绘制奈魁斯特图。
例5-3 绘制一阶惯性环节G（s)=3/(5s+1)的奈魁斯特图。
程序代码如下：G=tf(3,[5 1]);nyquist(G);hold on;set(G,'inputdelay',5);nyquist(G);hold on;set(G,'inputdelay',10);nyquist(G);hold on;title('Nyquist图');3.3 尼柯尔斯图尼科尔斯图又称为对数幅频率特性图，它以开环频率特性函数的对数幅值为纵坐标，以相角值为横坐标，以角频率为参变量绘制的频率特性图。
采用直角坐标。
纵坐标表示20lg| G(jω)|，单位是dB，线性刻度。
横坐标表示∠G(jω)，单位是度，线性分度。
在曲线上一般标注角频率ω的值作为参变量。
通常是先画出Bode图，再根据Bode图绘制尼科尔斯图。
4频率响应分析时域分析中的性能指标直观反映控制系统动态相应的特征，属于直接性能指标，而系统频率特性函数的某些特征可以用作间接性能指标。
⑴开环频率特性的性能分析基于开环频率特性函数的性能分析指标有如下两个：一是相角裕量γ，反映系统的相对稳定性；另一个是截止频率c ω，反映系统的快速性。
c ω是A(c ω)=1所对应的角频率，或对数幅频特性图上L(ω)穿越0分贝线的斜率，在采用渐近线作图时，两者略有不同。
⑵闭环频率特性的性能分析基于闭环频率特性函数的常用指标有两个：一是谐振峰值Mr ，反映系统的相对稳定性；另一个是频带宽度或者带宽频率B ω，定义为闭环幅频特性幅值M(ω)下降到0.707M(0)时对应的角频率，它反映了系统的快速性。
例2 用直接计算法，确定系统的谐振振幅和谐振频率。
已知一控制系统开环传递函数Go (s)=5.5/(s2+3s+5)，试求此系统的谐振振幅Mr 和谐振频率ωr 。
程序代码如下：Go=tf(5.5,[1 3 5]);[Mr,Pr,Wr]=mwr(Go)mwr 函数程序如下：function [Mr,Pr,Wr]=mwr(G)[mag,pha,w]=bode(G);magn=mag(1,:);phase=pha(1,:);[M,i]=max(magn);Mr=20*log10(M);Pr=phase(1,i);Wr=w(i,1);运行结果：Wr = 0.6915Mr = 0.8714Pr = -24.6446结果中的单位分别是分贝（dB ）、rad/s 和度。
例3 利用LTIView 工具，获得系统频率响应的谐振振幅和谐振频率。
以例2中的传递函数为例，确定系统的谐振振幅Mr 和谐振频率ωr 。
在命令窗口界面键入如下命令：>> Go=tf(5.5,[1 3 5])>> ltiview进入LTIView工具箱界面，对此系统进行分析。
Matlab控制仿真的简单例子sys=tf(1,[3 1 0],'inputdelay',6/100); %10/s(s+2) tf-----建立传递函数模型Ts=0.25;sysd=c2d(sys,Ts,'z'); % c2d ---连续的时间系统转化为离散的时间系统,即离散化.'t'表示双线性变换法% Find critical valuecrs Kcr,Wcr % 'z'表示带零阶保持器的Z变换法,可缺省% =============================figure(1);bode(sysd);grid on % 画伯德图[Gm,Pm,Wc,Wg]=margin(sysd) % margin---计算增益和相位裕度函数% Gm--增益裕度(即相位-180度处所对应的幅频特性绝对值的倒数)% Wg--增益裕度对应的频率(即相位-180度处的频率)% Wc--剪切频率(即幅频特性增益为0时的频率)% Pm--相角裕度(即剪切频率处,使系统达到临界稳定状态时所需的附加相移)Kcr=Gm;Wcr=Wg;Tcr=2*pi/Wcr; % 求临界状态时对应的周期pause % 暂停% Design PID Controller% =====================disp('parameters of PID')Kp=0.6*Kcr;Ti=0.5*Tcr;Td=0.125*Tcr;Ki=Kp*Ts/Ti;Kd=Kp*Td/Ts;disp('Model of PID Controller:')sysc=Kp+tf(Ki*[1 0],[1 -1],Ts)+tf(Kd*[1 -1],[1 0],Ts)% sysc=Kp+Ki*z/(z-1)+Kd*(z-1)/z% Check the Controller% ====================sysopen=sysc*sysd;sysclose2=feedback(sysopen,1); % feedback-----求闭环反馈figure(2);step(sysclose2,'b') % step----------作单位阶跃响应函数曲线grid onhold on。