实验三线性系统的频域分析

1、-自动控制理论上机实验报告学院：机电工程学院 班级：13级电信一班：*：实验三 线性系统的频域分析一、实验目的1掌握用MATLAB语句绘制各种频域曲线。2掌握控制系统的频域分析方法。二、根底知识及MATLAB函数频域分析法是应用频域特性研究控制系统的一种经典方法。它是通过研究系统对正弦信号下的稳态和动态响应特性来分析系统的。采用这种方法可直观的表达出系统的频率特性，分析方法比拟简单，物理概念明确。1频率曲线主要包括三种:Nyquist图、Bode图和Nichols图。1Nyquist图的绘制与分析MATLAB中绘制系统Nyquist图的函数调用格式为:nyquist(num,den) 频率响应w的围由软件自动设定nyquist(num,den,w) 频率响应w的围由人工设定Re,Im= nyquist(num,den) 返回奈氏曲线的实部和虚部向量，不作图例4-1:系统的开环传递函数为，试绘制Nyquist图，并判断系统的稳定性。图4-1 开环极点的显示结果及Nyquist图num=2 6;den=1 2 5 2;z,p,k=tf2zp(num,den); pnyquist(num,d

2、en)极点的显示结果及绘制的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)2Bode图的绘制与分析系统的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 6 30;den=1161

3、00 0; w=logspace(-2,3,100);bode(num,den,w)grid绘制的Bode图如图4-2(a)所示，其频率围由人工选定，而伯德图的幅值围和相角围是自动确定的。当需要指定幅值围和相角围时，则需用下面的功能指令: mag,phase,w=bode(num,den,w)图4-2(a) 幅值和相角*围自动确定的Bode图图4-2(b) 指定幅值和相角*围的Bode图mag,phase是指系统频率响应的幅值和相角，由所选频率点的w值计算得出。其中，幅值的单位为dB，它的算式为magdB=20lg10(mag)。指定幅值围和相角围的MATLAB调用语句如下，图形如图4-2(b)所示。num=0 0 15 30;den=116100 0;w=logspace(-2,3,100);mag,phase,w=bode(num,den,w); %指定Bode图的幅值围和相角围subplot(2,1,1); %将图形窗口分为2*1个子图，在第1个子图处绘制图形semilog*(w,20*log10(mag); %使用半对数刻度绘图，*轴为log10刻度，Y轴为线性刻度grid o

4、n*label(w/s-1); ylabel(L(w)/dB);title(Bode Diagram of G(s)=30(1+0.2s)/s(s2+16s+100);subplot(2,1,2);%将图形窗口分为2*1个子图，在第2个子图处绘制图形semilog*(w,phase);grid on*label(w/s-1); ylabel(0);注意:半Bode图的绘制可用semilg*函数实现，其调用格式为semilog*(w,L)，其中L=20*log10(abs(mag)。3Nichols图的绘制在MATLAB中绘制Nichols图的函数调用格式为:mag,phase,w=nichols(num,den,w) Plot(phase,20*log10(mag)图4-3 Nichols图例4-3:单位负反应的开环传递函数为，绘制Nichols图。对应的MATLAB语句如下，所得图形如图4-3所示:num=10; den=1 3 9 0;w=logspace(-1,1,500);mag,phase=nichols(num,den,w);plot(phase,20*log10(mag)

5、ngrid %绘制nichols图线上的网格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.7000pm = 64.6998wcg = 3.0000wcp = 1.1936如果系统的频域响应数据，还可以由下面的格式调用函数:Gm,Pm,Wcg,Wcp=margin(mag,phase,w)其中mag,phase,w分别为频域响应的幅值

6、、相位与频率向量。三、实验容1典型二阶系统绘制出，0.3，0.5，0.8，2的bode图，记录并分析对系统bode图的影响。解：当wn=6，时 程序如下num=0 0 36;den=11.2 36; w=logspace(-2,3,100);bode(num,den,w)grid 结果如图当wn=6，阻尼系数为0.3时 程序如下num=36;den=13.6 36; w=logspace(-2,3,100);bode(num,den,w)grid结果如图当wn=6，阻尼系数为0.5时 程序如下num=36;den=16 36; w=logspace(-2,3,100);bode(num,den,w)Grid结果如图当wn=6，阻尼系数为0.8时 程序如下num=36;den=19.6 36; w=logspace(-2,3,100);bode(num,den,w)Grid结果如图当wn=6，阻尼系数为0.8时 程序如下num=36;den=124 36; w=logspace(-2,3,100);bode(num,den,w)Grid结果如图2系统的开环传递函数为绘制系统的Nyquis

7、t曲线、Bode图，说明系统的稳定性，并通过绘制阶跃响应曲线验证。解：开环传递函数为 程序如下num=10;den=5 24 -5 0 0;z,p,k=tf2zp(num,den); pnyquist(num,den)GridNyquist曲线如下num=0 0 0 0 36;den=5 24 -5 0 0; w=logspace(-2,3,100);bode(num,den,w)Grid波特图如下开环传递函数为程序如下num=8 8;den=5 21 100 150 0 0;z,p,k=tf2zp(num,den); pnyquist(num,den)GridNyquist曲线如下num=0 0 0 0 0 36;den=5 21 100 150 0 0; w=logspace(-2,3,100);bode(num,den,w)Grid波特图如下开环传递函数程序如下num=0 0 0 1.6 4;den=0.0001 0.08 0.17 1 0;z,p,k=tf2zp(num,den); pnyquist(num,den)GridNyquist曲线如下num=0 0 0 1.6 4;den=0.0001 0.08 0.17 1 0; w=logspace(-2,3,100);bode(num,den,w)grid波特图如下 3系统的开环传递函数为。求系统的开环截止频率、穿越频率、幅值裕度和相位裕度。应用频率稳定判据判定系统的稳定性。解：程序如下num=1 1; den=0.1 1 0 0;gm,pm,wcg,wcp=margin(num,den);gm,pm,wcg,wcp 所求结果为gm = Infpm = 44.4594wcg = Infwcp = 1.2647四、总结经过这次实验，我学习了用MATLAB进展线性系统的频域分析，其中最主要的是三种图的绘制：Nyquist图、Bode图和Nichols图。结合课堂上对图形的学习，可以方便的通过MATLAB绘制出来图形，然后进展相应的分析，也能和自己的草图相结合进展比照，得出结论。. z.

《实验三线性系统的频域分析》由会员新**分享，可在线阅读，更多相关《实验三线性系统的频域分析》请在金锄头文库上搜索。