《线性系统频域分析实验报告》由会员分享,可在线阅读,更多相关《线性系统频域分析实验报告(13页珍藏版)》请在金锄头文库上搜索。
1、实验三、线性系统的频域分析法一,实验目的1,掌握 matlab 绘制波特图以及奈奎斯特图的方法。2,学会从波特图以及奈奎斯特图判定系统的稳定性。3,学会从波特图上求系统的稳定裕度。4,了解 k 值变化时对波特图幅频和相频曲线的影响。5,掌握 matalab 绘制系统零极点分布图的方法。6,学会从系统的零极点分布图判断系统的稳定性。二,实验原理1,从奈奎斯特图判定系统是否稳定的原理奈式稳定判据:反馈控制系统稳定的充分必要条件是半闭合曲线GH 不穿过(-1,0j )点,且逆时针包围临界点(-1,0j)点的圈数 R 等于开环传递函数正实部极点数 P具体方法是,先观察系统传递函数得出系统是否在 s 平
2、面的右半开平面由极点,得出 P 的值,在观察曲线从(-1,0j)点右侧穿越的次数,其中自上而下为正穿越,自下而上为负穿越,完整的一次穿越记为 N 半次穿越记为 0.5N,R=2N=2(N+ -N-) 而 Z=P-R,观察 Z 是否为零,Z 为零则系统是稳定的,Z 不为零时则系统是不稳定的。2,从波特图判定系统是否稳定的原理。从奈奎斯特稳定判定我们可以知道,要判定系统是否稳定就要观察曲线穿越(-1,0j)点次数,对应在波特图中,当取 w=wc 时,要满足 A(wc)=|G(jwc)H(jwc)|=1 L(wc)=20logA(wc)=0因此 wc 为分界点,对应到相频曲线上,观察在 wwc 时曲
3、线穿越-180 度的次数。然后计算方法和上面相同,既可以判定系统的稳定性。3,根据系统的零极点分布判断系统稳定性的原理三,实验内容A、 设单位负反馈系统的开环传递函数为K(S+1)/S(S+2)(S2+17S+4000) 其中K=1000(1)绘制波特图。(2) 观察绘制出的bode 图,分析系统的稳定性,并在图上求稳定裕度; (3) 绘制K=2000 时系统的bode 图,分析曲线的改变情况,并分析K 值变化时,对系统幅频响应和相频响应的影响。分析:1,绘制波特图matlab 文本命令为:s=tf(s); G=1000*(s+1)/(s*(s+2)*(s2+17*s+4000) Bode(G
4、) Grid onMargin(G),2,绘制出的波形为2,由于传递函数中可知v=1所以要在相频中增补从-90 度到0度的相频曲线,由波特图可以看出当L(w)=0dB时对应的频率值为 wc,在wwc时,在相频曲线中没有穿越-180度,所以可以知道R=0,又由传递函数可以知道P=0,所以Z=0 ,从而我们知道系统此时是稳定的,由裕度函数我们可以在图中求出幅值裕度Gm=36.7dB,相角裕度Pm=93.5度.剪切频率wc=0.126rad/s.3,改变系统的k值,令k=2000绘制此时的波特图,matlab 文本命令为;s=tf(s); G=2000*(s+1)/(s*(s+2)*(s2+17*s
5、+4000) Bode(G) margin(G) grid on得到系统的波特图为:由波特图可以看出,当k值变大后,对相频曲线没有影响,因为k 环节不提供相角,而对于幅频曲线来说当k值变为2000后相当于整个曲线向上平移了20lg2,从而使得幅值裕度和相角裕度改变了,幅值裕度为Gm=30.7dB,相角裕度为Pm=97度,剪切频率wc=0.256rad/s.B,设单位负反馈的开环传递函数为 G(s )=10/(s+5)/(s-1)(1) 绘制系统的Nyquist 曲线 (2) 分析系统的稳定性 (3) 根据系统的闭环零极点的分布图来分析系统的稳定性,和(2)得到的结果比较; 1,绘制Nyquis
6、t 曲线的matlab文本命令为:num=10; den=conv(1 5,1 -1); nyquist(num,den)绘制出的图形为:2,分析系统的稳定性,当w趋于零时G(Jw)等于-2所以曲线的起点在(-2,0j),由曲线我们可以看出,曲线在(-1,0j)左边有半次自上而下的正穿越所以N+=0.5,N=2(N+-N-)=1,所以R=1,由系统的传递函数可以知道P=1,所以Z=P-R=0,从而得出系统是稳定的。3,绘制系统的闭环零极点的分布图的matlab文本命令为:num=10; den=conv(1 5,1 -1);G=tf(num,den) pzmap(feedback(G,1)可以
7、得到系统的零极点分布图为:由图可知闭环传递在s平面的右开平面没有零极点分布,所以可以知道系统是稳定的。与2得到的结果是一致的。C,完成教材第五章练习5-11中各系统的Nyquist 图并分析系统的稳定性和闭环系统零极点分布图。1,G(s)=2/(2s+1)/(8s+1)2, G(s)=200/s2/(s+1)/(10s+1)3, G(s)=(80s+8)/s/(s2+s+1)/(0.5s+1)4,G(s)=(0.025s2+s+10)/s/(s+1)/(10s+1)绘制奈奎斯特图的matlab文本命令为;subplot(2,2,1);num=2;den=conv(2 1,8 1);nyquis
8、t(num,den);subplot(2,2,2);num=200;den=conv(1 1 0 0,10 1);nyquist(num,den);subplot(2,2,3);num=80 8;den=conv(1 1 1 0,0.5 1);nyquist(num,den);subplot(2,2,4);num=0.025 1 10;den=conv(1 1 0,10 1);nyquist(num,den);得到各传递函数的nyquist 图为:绘制系统的零极点分布图的matlab文本命令为:subplot(2,2,1);num=2;den=conv(2 1,8 1);G=tf(num,de
9、n);pzmap(feedback(G,1);subplot(2,2,2);num=200;den=conv(1 1 0 0,10 1);G=tf(num,den);pzmap(feedback(G,1);subplot(2,2,3);num=80 8;den=conv(1 1 1 0,0.5 1);G=tf(num,den);pzmap(feedback(G,1);subplot(2,2,4);num=0.025 1 10;den=conv(1 1 0,10 1);G=tf(num,den);pzmap(feedback(G,1);得到各系统的零极点分布图为:分析:1,在奈奎斯特图上可以看到
10、曲线没有从(-1,0j)的左侧穿越,所以R=0,由系统的传递函数可以知道P=0,所以Z=0,所以系统是稳定的,从闭环函数的零极点分布图可以知道,闭环系统在s平面的右半开平面没有零极点的分布同样可以得出系统是稳定的结论。2,系统是2型系统增补从0到-180度一段曲线,在奈奎斯特图上可以看到曲线有一次自下而上的负穿越,所以R=-2,由系统的传递函数可以知道P=0,所以Z=2,所以系统是不稳定的,从闭环函数的零极点分布图可以知道,闭环系统在s平面的右半开平面有两个极点的分布同样可以得出系统是不稳定的结论。3,在奈奎斯特图上可以看到曲线有从(-1,0j)的一次负穿越,所以R=2N=-2,由系统的传递函
11、数可以知道P=0,所以Z=p-R=2,所以系统是不稳定的,从闭环函数的零极点分布图可以知道,闭环系统在s平面的右半开平面有两个闭环极点,同样可以得出系统是不稳定的结论。4,因为系统是一型的所以要增补0到-90一段曲线,在奈奎斯特图上可以看到曲线有自上而下从(-1,0j)的左侧穿越,所以R=-2,由系统的传递函数可以知道P=0,所以Z=2,所以系统是不稳定的,从闭环函数的零极点分布图可以知道,闭环系统在s平面的右半开平面有两个极点,同样可以得出系统是不稳定的结论。D, 完成教材第五章练习5-11中各系统的波特图,利用波特图分析系统的稳定性与C的结果比较。绘制波特图的matlab文本命令为:sub
12、plot(2,2,1);num=2;den=conv(2 1,8 1);bode(num,den);grid onsubplot(2,2,2);num=200;den=conv(1 1 0 0,10 1);bode(num,den);grid onsubplot(2,2,3);num=80 8;den=conv(1 1 1 0,0.5 1);bode(num,den);grid onsubplot(2,2,4);num=0.025 1 1;den=conv(1 1 0,10 1);bode(num,den);grid on得到各系统的波特图为:分析:1,令L(w)=0得到wc在观察当wwc时,
13、相频曲线中穿越-180度的次数,从1图中可以看出穿越-180度的次数为0次所以R=0,由系统的开环传递函数可以知道P=0所以Z=0,系统是稳定的,与C实验中所得结果一致。2,令L(w)=0得到wc在观察当wwc时,相频曲线中穿越-180度的次数从2图中再增补相频曲线从-180度到0度,所以可以知道R=2,而P=0,所以Z=-2,系统是不稳定的,与C实验中所得结果一致。3,令L(w)=0得到wc在观察当wwc时,相频曲线中穿越-180度的次数从2图中再增补相频曲线从-90度到0度,又一次穿越-180度所以可以知道R=2,而P=0,所以Z=-2,系统是不稳定的,与C实验中所得结果一致。4,令L(w
14、)=0得到wc在观察当wwc时,相频曲线中穿越-180度的次数从2图中再增补相频曲线从-90度到0度,其中没有穿越=180度所以可以知道R=0,而P=0,所以Z=0,系统是稳定的,与C实验中所得结果一致。四,实验结论,及心得体会。1,在实验A中改变k的值后我们发现系统的幅频特性曲线改变,从而使得系统的稳定裕度也会改变,所以改变k值会影响系统的稳定性。2,从系统的波特图和奈奎斯特图我们很容易就能看出系统的稳定性,这给我们提供了一个很好判定系统稳定性的方法。3,通过这次实验,我掌握了各种图形的matlab绘制方法,加深了对课本上各种稳定性判别方法的理解,学会了用软件作图判定系统稳定性,进一步了解了各种系统参数对系统性能的影响。 4,由实验B我们可以知道,判断系统的稳定性还可以通过观察闭环传递函数的零极点分布,若闭环传递函数在虚轴右侧没有零极点,则系统是稳定的,若在虚轴右侧有极点分布,则系统是不稳定的。当极点在虚轴时系统是临界稳定的。