自动调节系统频域分析
班级 11081801
学号 1108180135
姓名 王佳炜
日期 2014.1.5
线性系统的频域分析
一、实验目的
1．掌握用MATLAB语句绘制各种频域曲线。
2．掌握控制系统的频域分析方法。
二、实验内容
1．典型二阶系统
绘制出
，
，0.3，0.5，0.8，2的bode图，记录并分析
对系统bode图的影响。
解:
程序如下:
num=[0 0 36];den1=[1 1.2 36];den2=[1 3.6 36];
den3=[1 6 36];den4=[1 9.6 36];den5=[1 24 36];
w=logspace(-2,3,100);
bode(num,den1,w)
grid
hold
bode(num,den2,w)
bode(num,den3,w)
bode(num,den4,w)
bode(num,den5,w)
分析：随着
的增大,伯德图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大.
2．系统的开环传递函数为
绘制系统的Nyquist曲线、Bode图，说明系统的稳定性，并通过绘制阶跃响应曲线验证。
解: 程序如下
奈氏曲线:
(1) num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5])));
w=logspace(-1,1,100);
nyquist(num1,den1,w)
(2) num2=[8,8];den2=conv([1,0],conv([1,0],conv([1,15],[1,6,10])));
w=logspace(-1,1,100);
nyquist(num2,den2)
(3) num3=[4/3,4];den3=conv([1,0],conv([0.02,1],conv([0.05,1],[0.1,1])));
w=logspace(-1,1,100);
nyquist(num3,den3)
分析：系统1，2 不稳定，系统3稳定。
伯德图:
num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5])));
num2=[8,8];den2=conv([1,0],conv([1,0],conv([1,15],[1,6,10])));
num3=[4/3,4];den3=conv([1,0],conv([0.02,1],conv([0.05,1],[0.1,1])));
bode(num1,den1)
grid
hold
bode(num2,den2)
bode(num3,den3)
分析：系统1，2 不稳定，系统3稳定。
。
阶跃响应曲线
（1）num1=[0,0,10];den1=conv([1,0],conv([1,0],conv([5,-1],[1,5])));
step(num1,den1)
grid
（2） num2=[8,8];den2=conv([1,0],conv([1,0],conv([1,15],[1,6,10])));
step(num2,den2)
grid
（3） num3=[4/3,4];den3=conv([1,0],conv([0.02,1],conv([0.05,1],[0.1,1])));
step(num3,den3)
grid
3．已知系统的开环传递函数为
。求系统的开环截止频率、穿越频率、幅值裕度和相位裕度。应用频率稳定判据判定系统的稳定性。
解：绘出系统伯德图，程序如下
num=[0 0 1 1];
den=[0.1 1 0 0];
w=logspace(-2,3,100);
bode(num,den,w)
[gm,pm,wcg,wcp]=margin(num,den);
gm,pm,wcg,wcp
grid
gm =0
pm =44.4594
wcg =0
wcp =1.2647
分析： 系统截止频率Wc=1.2647，相角裕度r=44.4594,幅值裕度hg=0，穿越频率Wg=0
因此系统稳定。
三．实验心得与体会
总结：通过这次实验，我掌握了各种图形的matlab绘制方法，加深了对课本上各种稳定性判别方法的理解，学会了用软件作图判定系统稳定性，进一步了解了各种系统参数对系统性能的影响。
第二篇：自动控制原理实验四 线性系统的频域分析
实验四线性系统的频域分析
一、实验目的
1．掌握用MATLAB语句绘制各种频域曲线。
2．掌握控制系统的频域分析方法。
二、基础知识及MATLAB函数
频域分析法是应用频域特性研究控制系统的一种经典方法。它是通过研究系统对正弦信号下的稳态和动态响应特性来分析系统的。采用这种方法可直观的表达出系统的频率特性，分析方法比较简单，物理概念明确。
1．频率曲线主要包括：Nyquist图、Bode图。
1）Nyquist图的绘制与分析
MATLAB中绘制系统Nyquist图的函数调用格式为：
nyquist(num,den) 频率响应w的范围由软件自动设定
nyquist(num,den,w) 频率响应w的范围由人工设定
[Re,Im]= nyquist(num,den) 返回奈氏曲线的实部和虚部向量，不作图
例4-1：已知系统的开环传递函数为
，试绘制Nyquist图，并判断系统的稳定性。
num=[2 6];
den=[1 2 5 2];
[z,p,k]=tf2zp(num,den);
p
nyquist(num,den)
注意： [Z,P,K] = TF2ZP(NUM,DEN)：求出零点、极点和增益！
极点的显示结果及绘制的Nyquist图如图4-1所示。由于系统的开环右根数P=0，系统的Nyquist曲线没有逆时针包围（-1，j0）点，所以闭环系统稳定。
p =
-0.7666 + 1.9227i
-0.7666 - 1.9227i
-0.4668
若上例要求绘制
间的Nyquist图，则对应的MATLAB语句为：
num=[2 6];
den=[1 2 5 2];
w=logspace(-1,1,100); 即在10-1和101之间，产生100个等距离的点
nyquist(num,den,w)
2）Bode图的绘制与分析
系统的Bode图又称为系统频率特性的对数坐标图。Bode图有两张图，分别绘制开环频率特性的幅值和相位与角频率
的关系曲线，称为对数幅频特性曲线和对数相频特性曲线。
MATLAB中绘制系统Bode图的函数调用格式为：
bode(num,den) 频率响应w的范围由软件自动设定
bode(num,den,w) 频率响应w的范围由人工设定
[mag,phase,w]=bode(num,den,w) 指定幅值范围和相角范围的伯德图
例4-2：已知开环传递函数为
，试绘制系统的伯德图。
num=[0 0 15 30];
den=[1 16 100 0];
w=logspace(-2,3,100);
bode(num,den,w)
grid
绘制的Bode图如图4-2(a)所示，其频率范围由人工选定，而伯德图的幅值范围和相角范围是自动确定的。当需要指定幅值范围和相角范围时，则需用下面的功能指令：
[mag,phase,w]=bode(num,den,w)
mag,phase是指系统频率响应的幅值和相角，由所选频率点的w值计算得出。其中，幅值的单位为dB，它的算式为magdB=20lg10(mag)。
指定幅值范围和相角范围的MATLAB调用语句如下，图形如图4-2(b)所示。
num=[0 0 15 30];
den=[1 16 100 0];
w=logspace(-2,3,100);
[mag,phase,w]=bode(num,den,w); %指定Bode图的幅值范围和相角范围
subplot(2,1,1); %将图形窗口分为2*1个子图，在第1个子图处绘制图形
semilogx(w,20*log10(mag)); %使用半对数刻度绘图，X轴为log10刻度，Y轴为线性刻度
grid on
xlabel(‘w/s^-1’); ylabel(‘L(w)/dB’);
title(‘Bode Diagram of G(s)=30(1+0.5s)/[s(s^2+16s+100)]’);
subplot(2,1,2);%将图形窗口分为2*1个子图，在第2个子图处绘制图形
semilogx(w,phase);
grid on
xlabel(‘w/s^-1’); ylabel(‘
(0)’);
注意：半Bode图的绘制可用semilogx函数实现，其调用格式为semilogx(w,L)，其中L=20*log10(abs(mag))。
2．幅值裕量和相位裕量
幅值裕量和相位裕量是衡量控制系统相对稳定性的重要指标，需要经过复杂的运算求取。应用MATLAB功能指令可以方便地求解幅值裕量和相位裕量。
其MATLAB调用格式为：
[Gm,Pm,Wcg,Wcp]=margin(num,den)
其中，Gm,Pm分别为系统的幅值裕量和相位裕量，而Wcg,Wcp分别为幅值裕量和相位裕量处相应的频率值。
另外，还可以先作bode图，再在图上标注幅值裕量Gm和对应的频率Wcg，相位裕量Pm和对应的频率Wcp。其函数调用格式为：
margin(num,den)
例4-4：对于例4-3中的系统，求其稳定裕度，对应的MATLAB语句如下：
num=10; den=[1 3 9 0];
[gm,pm,wcg,wcp]=margin(num,den);
gm,pm,wcg,wcp
gm = 2.7000
pm = 64.6998
wcg = 3.0000
wcp = 1.1936
如果已知系统的频域响应数据，还可以由下面的格式调用函数：
[Gm,Pm,Wcg,Wcp]=margin(mag,phase,w)
其中（mag,phase,w）分别为频域响应的幅值、相位与频率向量。
三、实验内容
1．典型二阶系统
绘制出
，
，0.3，0.5，0.8，2的bode图，记录并分析
对系统bode图的影响。
2．系统的开环传递函数为
绘制系统的Nyquist曲线、Bode图，说明系统的稳定性，并通过绘制阶跃响应曲线验证。
3．已知系统的开环传递函数为
。求系统的开环截止频率、穿越频率、幅值裕度和相位裕度。应用频率稳定判据判定系统的稳定性。
四、实验报告
1．根据内容要求，写出调试好的MATLAB语言程序，及对应的结果。
2. 记录显示的图形，根据实验结果与各典型环节的频率曲线对比分析。
3. 记录并分析对二阶系统bode图的影响。
4．根据频域分析方法分析系统，说明频域法分析系统的优点。
5．写出实验的心得与体会。
五、预习要求
1. 预习实验中的基础知识，运行编制好的MATLAB语句，熟悉绘制频率曲线的图形函数nyquist（）、bode（）。
2. 掌握控制系统的频域分析方法，理解系统绝对稳定性和相对稳定性的判断方法。