用MATLAB行系统频率特性曲线绘制.docVIP

自动控制原理 课程验证性实验报告 实验名称 用MATLAB进行系统频率特性曲线绘制 实验时间 2014年 12月10日 学生姓名 实验地点 070312 同组人员 无 专业班级 新能源1201B 1、实验目的 1）熟练掌握使用MATLAB命令绘制控制系统奈氏图的方法； 2）熟练掌握使用MATLAB命令绘制控制系统伯德图的方法； 3）加深理解控制系统奈氏稳定判据的实际应用。 2、实验主要仪器设备和材料： 计算机一台 matlab软件2010a版本 3、实验内容和原理： 原理：（1）幅相频率特性曲线，以角频率为参变量，当从变化时，频率特性构成的向量在复平面上描绘的曲线称为幅相频率特性曲线。（2）奈氏判据，奈氏判据是利用系统开环频率特性来判断闭环系统稳定性的一个判据，便于研究当系统结构改变时对系统稳定性的影响。 其内容是：反馈控制系统稳定的充分必要条件是当从变到时，开环系统的奈氏判据不穿过点逆时针包围临界点的圈数等于开环传递函数的正实部极点数。（1）对于开环稳定的系统，闭环系统稳定的充分必要条件是：开环系统的奈氏曲线不包围点。反之，则闭环系统是不稳定的。 （2）对于开环不稳定的系统，有个开环极点位于右半平面，则闭环系统稳定的充分必要条件是：当从变到时，开环系统的奈氏判据逆时针包围点次。 内容：1）绘制控制系统奈氏图 格式一：nyquist(num,den) 格式二:nyquist(num,den,w)格式三：[re,im,w]=nyquist(num,den) a.开环开环传递函数，绘制其Nyquist图。 b.已知，绘制Nyquist图，判定系统的稳定性。 c.已知系统开环传递函数为，要求：分别作出和时的Nyquist图。比较两图的区别与特点。如果该系统变成Ⅱ型系统，即，情况又发生怎么样的变化？ d.若，令，分别绘制时系统的Nyquist图并保持，比较分析系统开环增益不同时，系统Nyquist图的差异，并得出结论。令，分别绘制时系统的Nyquist图并保持，比较分析不同时，系统Nyquist图的差异，并得出结论。 e.二阶系统传递函数为，试用MATLAB绘制出不同和的伯德图。 f.系统的开环传递函数为求系统的幅值裕度和相角裕度，并求其闭环阶跃响应。 g.系统的开环传递函数为，求系统的幅值裕度和相角裕度。 4、实验方法、步骤： a）num=10,den=[1 2 10];w=0:0.1:100;axis([-1,1.5,-2,2]);nyquist(num,den) b) num=0.5;den=[1 2 1 0.5];figure(1);nyquist(num,den) c）k=1,T1=3,T2=2,T1T2 num=4;den=[2 1 0];nyquist(den,num) k=1,T1=2,T2=3,T1T2 num=3;den=[3 1 0];nyquist(den,num) k=1,T1=3,T2=2,T1T2 num=[3 1];den=[2 1 0 0];nyquist(den,num)）k=1,T1=2,T2=3,T1T2 num=[2 1];den=[3 1 0 0];nyquist(den,num) d) (1)u=1,k=1 num=1;den=conv([1 0],conv([1 1],[1 2]));nyquist(den,num) u=1,k=2 num=2;den=conv([1 0],conv([1 1],[1 2]));nyquist(den,num) u=1,k=10 num=2;den=conv([1 0],conv([1 1],[1 2]));nyquist(den,num) (2)k=1,u=2 num=2;den=conv([1 0 0],conv([1 1],[1 2]));nyquist(den,num) k=1,u=3 num=2;den=conv([1 0 0 0],conv([1 1],[1 2]));nyquist(den,num) k=1,u=4 num=2;den=conv([1 0 0 0 0],conv([1 1],[1 2]));nyquist(den,num) e)(1) 为固定值，变化时，运行下面的程序， wn=1,zet=[0:0.1:1,2,3,4,5]; hold on for i=1:length(zet) num=wn^2;den=[1,2*zet(i)*wn,wn^2]; bode(num,den); end grid on hold off 2)为固定值，变化时，运行下面的程序 Wn=[0.1:0.1:1];zet=0.707; hold on for i=1:le