实验四用MATLAB实现线性系统的频域分析(已完成)8901.pdfVIP

实验四 用MATLAB 实现线性系统的频域分析 [实验目的] 1．掌握MATLAB 平台下绘制典型环节及系统开环传递函数的Bode 图和Nyquist 图（极 坐标图）绘制方法； 2．掌握利用Bode 图和Nyquist 图对系统性能进行分析的理论和方法。 [实验原理] 一、绘制Bode 图和Nyquis 图 1．Bode 图绘制 采用bode() 函数，调用格式： ①bode(sys) ；bode(num,den) ； 系统自动地选择一个合适的频率范围。 ②bode(sys ，w) ； 其中w( 即ω)是需要人工给出频率范围，一般由语句w=logspace(a,b,n)给出。 a b logspace(a,b,n) ：表示在10 到10 之间的 n 个点,得到对数等分的w 值。 ③bode(sys,{wmin,wmax}) ； 其中{wmin,wmax}是在命令中直接给定的频率w 的区间。 以上这两种格式可直接画出规范化的图形。 ④[mag,phase,ω]=bode(sys)或[m,p]=bode(sys) 这种格式只计算Bode 图的幅值向量和相位向量，不画出图形。 m 为频率特性G(jω ) 的幅值向量; p 为频率特性G(jω )的幅角向量，单位为角度 （° ）。 w 为频率向量，单位为[弧度]/秒 。 在此基础上再画图，可用： subplot(211);semilogx(w,20*log10(m) % 对数幅频曲线 subplot(212);semilogx(w,p) %对数相 频曲线 ⑤bode(sys1,sys2,…,sysN) ; ⑥bode((sys1,sys2,…,sysN，w); 这两种格式可在一个图形窗口同时绘多个系统的bode 图。 2. Nyquist 曲线的绘制 采用nyquist() 函数调用格式： ① nyquist(sys) ; ② nyquist(sys,w) ; 其中频率范围w 由语句w=w1:Δw:w2确定。 ③nyquist(sys1,sys2,…,sysN) ; ④nyquist(sys1,sys2,…,sysN,w); ⑤ [re,im,w]=nyquist(sys) ; re—频率响应实部 im—频率响应虚部 使用命令axis()改变坐标显示范围，例如axis([-1,1.5,-2,2])。 ⑥当传递函数串有积分环节时 ω=0 处会出现幅频特性为无穷大的情况，可用命令 axis()，自定义图形显示范围，避开无穷大点。 [实验内容] 1．作各典型环节的Bode 图和Nyquist 图。 1 G(s )  s  1 程序： n=[1];d=[1 1];sys=tf(n,d);subplot(211);bode(sys);subplot(212);nyquist(sys) 运行结果: Bode Diagram 0 ) B d ( e d -20 u t i n g a M -40 0