北京科技大学控制实验报告材料4

实验四线性系统的频域分析实验目的熟练掌握运用MATLAB命令绘制控制系统伯德图的方法;了解系统伯德图的一般规律及其频域指标的获取方法；熟练掌握运用伯德图分析控制系统稳定性的方法；设计超前校正环节并绘制Bode图；设计滞后校正环节并绘制Bode图。实验原理及容1、频率特性函数G(加)。设线性系统传递函数为：bsm+bsm-ihfbs+bG(S)=―01m—1masn+asn++as+a则频率特性函数为：b(加)m+b(加)m-1+•••+b(加)+bG(jw)=--01m—1ma(jw)n+a(jw)n+fa(jw)+a由下面的MATLAB语句可直接求出G(jw)。i=sqrt(-1)%求取-1的平方根GW=polyval(num，i*w)./polyval(den，i*w)其中(num,den)为系统的传递函数模型。而w为频率点构成的向量，点右除(./)运算符表示操作元素点对点的运算。从数值运算的角度来看，上述算法在系统的极点附近精度不会很理想，甚至出现无穷大值，运算结果是一系列复数返回到变量GW中。2、用MATLAB作伯德图控制系统工具箱里提供的bode()函数可以直接求取、绘制给定线性系统的伯德图。当命令不包含左端返回变量时，函数运行后会在屏幕上直接画出伯德图。如果命令表达式的左端含有返回变量，bode()函数计算出的幅值和相角将返回到相应的矩阵中，这时屏幕上不显示频率响应图。命令的调用格式为：[mag,phase,w]=bode(num,den)[mag,phase,w]=bode(num,den,w)或[mag,phase,w]=bode(G)[mag,phase,w]=bode(G,w)矩阵mag、phase包含系统频率响应的幅值和相角，这些幅值和相角是在用户指定的频率点上计算得到的。用户如果不指定频率w,MATLAB会自动产生w向量，并根据w向量上各点计算幅值和相角。这时的相角是以度来表示的，幅值为增益值，在画伯德图时要转换成分贝值，因为分贝是作幅频图时常用单位。可以由以下命令把幅值转变成分贝：magdb=20*log10(mag)绘图时的横坐标是以对数分度的。为了指定频率的围，可采用以下命令格式:logspace(d1,d2)或logspace（d1,d2,n）第一种格式是在指定频率围按对数距离分成50等分的，即在两个十进制数31=10d和①2=10^2之间产生一个由50个点组成的分量，向量中的点数50是一个默认值。例如要在31=0.1弧度/秒与32=100弧度/秒之间的频区画伯德图，则输入命令时，d=log（3）,d=log（3）在此频区自动按对数距离等分成50个频率点，返回到工作11012102空间中，即w=logspace（-1，2）要对计算点数进行人工设定，则采用第二种格式。例如，要在31=1与32=1000之间产生100个对数等分点，可输入以下命令：w=logspace（0，3，100）在画伯德图时利用以上各式产生的频率向量w可以很方便地画出希望频率的伯德图。由于伯德图是半对数坐标图且幅频图和相频图要同时在一个绘图窗口中绘制，因此，要用到半对数坐标绘图函数和子图命令。（1）对数坐标绘图函数利用工作空间中的向量x，y绘图，要调用plot函数，若要绘制对数或半对数坐标图，只需要用相应函数名取代plot即可，其余参数应用与plot完全一致。命令公式有：semilogx（x，y，s）;只对x轴进行对数变换，y轴仍为线性坐标。semilogy(x,y,s);y轴取对数变换的半对数坐标图。Loglog(x，y，s);全对数坐标图，即x轴和y轴均取对数变换。子图命令MATLAB允许将一个图形窗口分成多个子窗口，分别显示多个图形，这就要用到subplot()函数，其调用格式为：subplot(m,n,k)该函数将把一个图形窗口分割成mxn个子绘图区域，m为行数，n为列数，用户可以通过参数k调用各子绘图区域进行操作，子图区域编号为按行从左至右编号。对一个子图进行的图形设置不会影响到其它子图，而且允许各子图具有不同的坐标系。例如，subplot(4,3,6)则表示将窗口分割成4x3个部分。在第6部分上绘制图形。MATLAB最多允许9x9的分割。【例6-3】给定单位负反馈系统的开环传递函数为：G(s)=10(+,试画出伯德图。,>'||。s(s+7)利用以下MATLAB程序，可以直接在屏幕上绘出伯德图如图6-2。>>num=10*[1,1];den=[1,7,0];bode(num,den)gridtitle('BodeDiagramofG(s)=10*(s+1)/[s(s+7)]')

该程序绘图时的频率围是自动确定的从0.01弧度渺到30弧度渺且幅值取分贝值，轴取对数，图形分成2个子图，均是自动完成的。i_I_ij-i.||_ms省WLi沼巨(E心p〕心wraLldSystem:sysFrequency(rad/sec):7.25Magnitudei_I_ij-i.||_ms省WLi沼巨(E心p〕心wraLldtt~i~i~iri~rirr~ttrrini1III1I11I11I1i11Ii1III1I11I11I1i11I尸/irrri~iiiiiiiiii-■---4...11111>\|1System:sysFrequency(rad.|,sec):7.24Phase(deg):-53.9■|・|7TLIIIIIIIIIIIIIIIIIIIIIIII1l-i■十H111l/TIlliiijX1iiiiF—i-1-1-i-i-tiiiiiiiiiiiiiiiiii\1—-i-i-nmi1-r11\i11i11i1Mill-l-IT十f1111111I1111111111111111111111iiiiiiiiiiiiiiiiiiIII111IIIIIIIl-l———I——I—I—I—1—i4-14-—二T3-1111111I1111如果希望显示的频率围窄一点，则程序修改为：>>num=10*[1,1];den=[1,7,0];w=logspace(-1,2,50);%从0.1至100，取50个点。[mag,phase,w]=bode(num,den,w);magdB=20*log10(mag)%增益值转化为分贝值。%第一个图画伯德图幅频部分。subplot(2,1,1);semilogx(w,magdB,'-r')%用红线画gridtitle('BodeDiagramofG(s)=10*(s+1)/[s(s+7)],)xlabel(1Frequency(rad/s)1)ylabel(1Gain(dB)1)%第二个图画伯德图相频部分。subplot(2,1,2);semilogx(w,phase,W);gridxlabel(】Frequency(rad/s)i)ylabel('Phase(deg)')Freqiiencyfrad/s)_-'fF-t-m03_J_LL图6-3用户指定的频率点画出的伯德图Freqiiencyfrad/s)_-'fF-t-m03_J_LL修改程序后画出的伯德图如6-3所示。【例6-4】已知控制系统开环传递函数：g(s)H(s)=—1^，绘制其Bode图并确定谐s2+2s+10振峰值Mr和谐振频率wr。绘制Bode图(略)，插值函数spline(m,w,mr)求谐振频率>>[jtljPgw]="bode(iiunijden);>>mr=maa(m)mr二1.6667>>wr=spline(m.j虬mr)wr=2.82845、用MATLAB求取稳定裕量同前面介绍的求时域响应性能指标类似，由MATLAB里bode()函数绘制的伯德图也可以采用游动鼠标法求取系统的幅值裕量和相位裕量。此外，控制系统工具箱中提供了margin()函数来求取给定线性系统幅值裕量和相位裕量，该函数可以由下面格式来调用：margin(num,den);给定开环系统的数学模型，作Bode图，并在图上方标注幅值裕度Gm和对应频率3g,相位裕度Pm和对应的频率wco[Gm,Pm,Wcg,Wcp]=margin(G);可以看出，幅值裕量与相位裕量可以由LTI对象G求出，返回的变量对(Gm,Wcg)为幅值裕量的值与相应的相角穿越频率，而(Pm,Wcp)则为相位裕量的值与相应的幅值穿越频率。若得出的裕量为无穷大测其值为Inf,这时相应的频率值为NaN(表示非数值)，Inf和NaN均为MATLAB软件保留的常数。如果已知系统的频率响应数据，我们还可以由下面的格式调用此函数。[Gm,Pm,Wcg,Wcp]=margin(mag,phase,w);其中(mag,phase,w)分别为频率响应的幅值、相位与频率向量。【例6-6】考虑一个新的系统模型，开环传递函数为：G(s)=100((")2,由下面(s+1)(s2+s+9)，MATLAB语句可直接求出系统的幅值裕量和相位裕量>>G=tf(100*conv([1,5],[1,5]),conv([1,1],[1,1,9]));[Gm,Pm,Wcg,Wcp]=margin(G)结果显示Gm=Pm=Inf85.4365Wcg=Wcp=NaN100.3285再输入命令>>G_c=feedback(G,1);step(G_c)gridxlalel('Time(sec)')ylalel('Amplitude')可以看出，该系统有无穷大幅值裕量，且相角裕量高达85.4365o。所以系统的闭环响图6-5较理想的系统响应应是较理想的，闭环响应图如图6-5。【自我实践6-1】某单位负反馈系统的开环传递函数gG)=—k求(1)当k=4时，s(s+1)(s+2)计算系统的增益裕度，相位裕度，在Bode图上标注低频段斜率，高频段斜率及低频段、高频段的渐近相位角。(2)如果希望增益裕度为16dB，求出响应的k值，并验证(1)clcnum=[4];den=[1,3,2,0];G=tf(num,den)[Gm,Pm,Wcg,Wcp]=margin(G)bode(num,den)gridTransferfunction:4s人3+3s人2+2sGm=1.5000Pm=11.4304Wcg=1.4142Wcp=1.1431iw、«?ini$po<oapo9「iw¥po<oapousiiQ心-52」JA3uunl>B」Ll_E&IHaSPOE(Bep)asELiri(8P)apniiuBeuu(uCDp、Enu)pno二ocxrrnrlnu①P3Ln6.0TEnuupILn6.0H>l(")«寸。孚.-tinEd«钮归ffiu①po<oapoooLn.<—InEow钮湄觐。«OR1[Gm,Pm,Wcg,Wcp]=margin(G)bode(num,den)gridTransferfunction:0.951sA3+3sA2+2sGm=6.3091Pm=54.7839Wcg=1.4142Wcp=0.42766、系统对数频率稳定性分析【自我实践6-2】系统开环传递函数g(s)=k，试分析系统的稳定性。s(0.5s+1)(0.k+1)经计算，当k=12时系统的增益裕度、相角裕度为零。clcnum=[12];den=[0.05,0.6,1,0];G=tf(num,den)[Gm,Pm,Wcg,Wcp]=margin(G)bode(num,den)gridTransferfunction:120.05s人3+0.6s人2+sGm=1Pm=9.5374e-006Wcg=4.4721BodeDiagramFrequency(rad/ssc)BodeDiagramGm=0dIB(at4.47radteec),Pm=9.54&-006deg(at4.47radteec)50当k<12时clcnum=[1];den=[0.05,0.6,1,0];G=tf(num,den)[Gm,Pm,Wcg,Wcp]=margin(G)bode(num,den)gridTransferfunction:10.05s人3+0.6s人2+sGm=12.0000Pm=60.4231Wcg=4.4721Wcp=0.9070BodeDiagramO-oooO5O59--1-1-mssmLIMM1o~110s1o11广Frequency(rad/sec)BotieDiagramGm=21.6dB(at4.47rad/sec),Pm=60.4deg(at0.907rad/sec^50I■11…i11…|ii「i系统产生衰减震荡，增益裕度和相角裕度大于0,系统稳定。当K>12时clcnum=[50];den=[0.05,0.6,1,0];G=tf(num,den)[Gm,Pm,Wcg,Wcp]=margin(G)bode(num,den)gridTransferfunction:500.05s人3+0.6s人2+sGm=0.2400Pm=-27.5774Wcg=4.4721Wcp=8.5945Boiie-DiagramBoiie-Diagram1Q-110°101DO1广Frequency(rad/sec)BadeDiagramGm=-12.4dB(at4.47rad/sec),Pm=-27.6deg(atS.591DO10Frequancy(rad/sec)增益裕度相角裕度均小于0,系统不稳定。【自我实践6-3】某单位负反馈系统的开环传递函数g(s)=址求(1)绘制s(0.01s+1)(0.1s+1)Bode图，在幅频特性曲线上标出低频段斜率、高频段斜率、开环截止频率和中频段穿越频率；在幅频特性曲线标出：低频段渐近相位角、高频段渐近相位角和-180°线的穿越频率。(2)计算系统的相位裕度g和幅值裕度h，并确定系统的稳定性(1)低频段斜率-20dB/10dec高频段-60dB/10dec,幅值穿越频率Wcp=16.3053低频段渐进相位角-90度，高频段渐进相位角-270度；-180度线穿越频率为Wcg=31.6228clcnum=[31.6];den=[0.001,0.11,1,0];G=tf(num,den)[Gm,Pm,Wcg,Wcp]=margin(G)bode(num,den)gridTransferfunction:31.60.001s人3+0.11s人2+sGm=3.4810Pm=22.2599Wcg=31.6228Wcp=16.3053BorieDiaqramSystem:s.ysFrequency(rad/s.ec):163二；血gnituWE日)：-口,旌药：■：System:aysFrequency(rad/s-ec):0.127MagnitudefdB):47.9System:sysFregency(rad/sec)ude(d0):S1.2System:sysFrequencyCra^sec):Llignitude(dB):-7S.ZSystem:5-ysFrequency(rsd/sec):0.13Phase(deg}:-91日System:s-ysFFrequency(rad/sec):9.21e+0-03PPhass(<Jeg):-269System:aysFrequency(rad/sec):1.24ei-Qt}5-Mggnitude[dB):-121iFrequency[rad/sBG):^1.4'_i■Phas-e(deg):-1B0BodeDiagramGm=10.SdB(at31.6rari/sac),Pm=22.Sdeg[at16.Srad/sec)【自我实践6-4】某单位负反馈系统的开环传递函数g(s)=心1)，令k=1作bode图，2(0.1s+1)应用频域稳定判据确定系统的稳定性，并确定使系统获得最大相位裕度的增益k值。clcnum=[1,1];den=[0.1,1,0,0];G=tf(num,den)[Gm,Pm,Wcg,Wcp]=margin(G)bode(num,den)gridGm=0Pm=44.4594Wcg=0Wcp=1.2647>>

100BadeDiagramGm=-InfdB(atQrad/sec),Pm=44.5deg{.at1.26rad/sec)100BadeDiagramGm=-InfdB(atQrad/sec),Pm=44.5deg{.at1.26rad/sec)-2-1Cl12101010101010Frequency（rad/sec）（1）比例环节：K（K=10、K=30）（2）惯性环节：_^（K=1、K=10、T=0.1、1）Ts+1（3）积分环节：K（K=1、K=10）s（4）微分环节：K（K=1、K=10）（5）二阶惯性环节：s2+^^+.2（K=1、K=10、&=0.1、&=1&=5、«=1）

环节传递函数BODE图形参数变化影响的分析1K值变大，系统的幅频特性曲线向上平移，相频特性不变，即号直只影响系统幅频特性的起点，不改变其形状，且对相频特性无影响。增加积分环节，系统性能编号，但会使系统稳定性变差。比例系数增大使Bode图中幅频特性的幅值增益增加，同时高频段的斜率增大。相频特性而言，中频段的下降速度随比例系数的增大而增大。时间常数不同使图形拐点出现的地方不同，时间常数较大的地方较先出现拐点。2345(1)比例环节K=10K=3021505Q-.W9.21ms£m匚麟M505―1(顷OJP》■的的r-dBodeDiagramT1F+•-+-T土■1卜Frequency(rad/sec)1O131BodeDiagramsS5951Q-.3也21&3N2(顷OJP》心slgr-d-O…TTT-L.F++T王1010Frequency(rad/sec)(2)惯性环节K=1,T=0.11/(0.1s+1)1>S+1)Magnitude(dB)khleij32Phase(deg)F「EU催ncy■Tadfs-cpo)lcvo.ls+1)Magnitude(dB)i挡hPhase(deg)KH1QTHP1s+1Magnitude(dB)Phase(deg)KnlQTnlBMWDmwmF「EU催ncy■Tadfs-cpo)1I「1—-1T1□回BOBBodeDiagramFrequency(radysec)BodeDiagramFrequency(radysec)FileEditViewI_nsertToolsDesktopWindowHelp笛DdQ受鼠贮陷晏K-□s■口(3)积分环节K=11

MFigureL_________I1=1回・|£3^1FileEditViewI_nsertToolsDesktopWindowHelpQSidaBS■口10°1O1Frequencytrsd/sec)LK=10Diagram10Diagram(33邳ps」止)听nE」jujbjBsio事住口日MCJgnrtudEas(□■IE3口I冒卜队替⑥总在惟国I胃Fl日口di9|-|A^opui；\\dopisaQsiooj_^issujmwi/\;ip39ii-|虱⑥曼贰-a□s1010Frequencyfrsd/sec)K=10FileEditViewInsertToolsDesktopWindowHelp笛OSfeRR凹倒土丈'息□E■0Diagram1oa1o1Frequency(rad^sec)(5)二阶惯性环节K=1、&=0.1、®=1Diagram1sA2+0.2s+1BodeDiagram2O1OQQooo-1-2-2-4(CQS心口n菖raM53-1Bo-1(50JP)心切国L-_LLK=1、&=1、«=11■■Fjgure1□□回7-EditViev弓98/I_nsertToolsDesktopWindow0Help□■■口BodeDiagramDoo

-2-4-S

([ns理muDoo

-2-4-S

([ns理muojra巨Frequency(rad/sec)K=1、&=5、®=1111L11JIIIIIIIIIIIiiiiII111111111111111111IIIIIIIIIIIIIIIIIIii■iii11|-J_111■1iiiiii11iiiiiii■|aiiiiia1111■iiii■i■1I1111II1111111II1iiiiia1111■iiii■i■lllll11III1I1111111IIIII1I11」_L1LUU1-1_I_LLI1U___i__LJJJPtekL_-L1J1LUI」_」_LLU」UL_1_l_l_ULLiiii■ii11i■iiiiiii_iiiii■11ii■iiii■i■iiiiii■IIIIIII11IIIIIIIIIiiiiiiiii\r-._iiii1111iiii111iiiIIIIIIIiiii■ii11i■iiiiiiiiitiii111r-—j_iii■11ii■i111■i■iiiiii■IIIIIII11IIIIIIIIIiiiiiiiii11TXJ■111111111II1IIIIIII—一_-i一r-rtrmti~~-t-t-i—r~mtn~1I--!-I~ITnT-----i-titi--i-b-rmn--r-T-Li-i-1n-1111111II111111111111111II1IIIIIIii-'..IIIIIIIiiI11111IIiiii■ii11i■iiiiiii1111I111111ii■III1I11111111111111111111111111II1IIIIIIiiIIIIII11I11111IIiiii■ii11i■iiiiiii1111I1111iiiii■11iF■iiii11iIII1I11111111111111111111111111II1iTiT11iii11111IIiiii■ii11i■iiiiiii1111I1111iiiii■11ii■'■■j11111iIII1I11111111111111111111111111II1111111II111I'H1111111111IIiiii■ii11i■iiiiiii11I111111iiiii■iiii■iii■IIIIII■---■<-■卓*—~-4-*-1-hM+M-M+■■■--1-+d+bldfe---►-4-l-l-IH+iiii■ii11i■iiiiiii11I11111111ii■IIIIII■111111111111111111111111II1II1111111IIiiii■ii11i■iiiiiii11I111111iiiii■iiii■iiii■iTx—1111111II111111111111111II1111111II1111111II11Mill1111■11111■11111111111I1111IIIIII1111■IIII11II111111111111111111111111II11I1111II1111111II11V-ii1II1111■11111■11111111111I1111IIIIII1111■IIII11IiiZ、i1111111111111111111I111111II11I1111II1111111II11111MJIiiiiii11iiiii■ii1111111111liiiii11ii■iiii■1■IiiiiiT-.K=10、&=0.1、®=1sA2+0.2s+1FileEditView[nssErtToolsDesktopWindowHelp也BS■口BodeDiagram4O3Q20WD-10-2Q(CQS心口n菖raM53-■1Bo-1(50JP)心切国L-CLK=10、&=1BodeDiagram4O3Q20WD-10-2Q(CQS心口n菖raM53-■1Bo-1(50JP)心切国L-CL10ri-rttr|iiiii~T~-TT~1~II~ITiiri-rttr|iiiii~T~-TT~1~II~ITiiiiiii1111111■iiiiiiIIIIIIIiiiiiii1111111■i■iiiiIIIIIIIiiiiiii1111111■-~i~rttri一TLT~i~~rr*-.iiiiii11111ii---.L■iiiiIIIII■i■r■..Iiiii11111ii■ir-.jiiiI11IIiiiiinili11111ii■iiii'r.j■11IIiiiiiitir-._11111ii■iI11IIiiiiii111s■1111ii■i■iiiiJIIII■i■iiiiii1111ii■i■iiiii气II■i■_j_L11Li-1--X-j_1J_-i-L■i■iiiiII1ii■iIlllI11IX11iI■iiiiiiIi1X■iiiilIIII1iiisIi■iiiiiiI•i1>1■J1iilIIII1ii11T-.J■iiiiiiI■i11ih06n吨氏凹倒V.北』也口启Frequency(rad/sec)K=10、&=5、®=110111L11JIIIIIIIIIIIiiiiII111111111111111111IIIIIIIIIIIIIIIIIIii■iii11|-J_111■1iiiiii11iiiiiii■|aiiiiia1111■iiii■i■1I1111II1111111II1iiiiia1111■iiii■i■lllll11III1I1111111IIIII1I11」_L1LUU1-1_I_LLI1U___i__LJJJPtekL_-L1J1LUI」_」_LLU」UL_1_l_l_ULLiiii■ii11i■iiiiiii_iiiii■11ii■iiii■i■iiiiii■IIIIIII11IIIIIIIIIiiiiiiiii\r-._iiii1111iiii111iiiIIIIIIIiiii■ii11i■iiiiiiiiitiii111r-—j_iii■11ii■i111■i■iiiiii■IIIIIII11IIIIIIIIIiiiiiiiii11TXJ■111111111II1IIIIIII—一_-i一r-rtrmti~~-t-t-i—r~mtn~1I--!-I~ITnT-----i-titi--i-b-rmn--r-T-Li-i-1n-1111111II111111111111111II1IIIIIIii-'..IIIIIIIiiI11111IIiiii■ii11i■iiiiiii1111I111111ii■III1I11111111111111111111111111II1IIIIIIiiIIIIII11I11111IIiiii■ii11i■iiiiiii1111I1111iiiii■11iF■iiii11iIII1I11111111111111111111111111II1iTiT11iii11111IIiiii■ii11i■iiiiiii1111I1111iiiii■11ii■'■■j11111iIII1I11111111111111111111111111II1111111II111I'H1111111111IIiiii■ii11i■iiiiiii11I111111iiiii■iiii■iii■IIIIII■---■<-■卓*—~-4-*-1-hM+M-M+■■■--1-+d+bldfe---►-4-l-l-IH+iiii■ii11i■iiiiiii11I11111111ii■IIIIII■111111111111111111111111II1II1111111IIiiii■ii11i■iiiiiii11I111111iiiii■iiii■iiii■iTx—1111111II111111111111111II1111111II1111111II11Mill1111■11111■11111111111I1111IIIIII1111■IIII11II111111111111111111111111II11I1111II1111111II11V-ii1II1111■11111■11111111111I1111IIIIII1111■IIII11IiiZ、i1111111111111111111I111111II11I1111II1111111II11111MJIiiiiii11iiiii■ii1111111111liiiii11ii■iiii■1■IiiiiiT-.【综合实践】已知系统开环传递函数:1)0型系统:K(Ts+1)(T2s+1),(K=10、七二1、t2=12)Transferfunction:1012sA2+13s+1Gm=InfPm=62.6208Wcg=InfWcp=0.6830E1™碧n拦冒w(5心p)口！sfsLhLL__\ME1™碧n拦冒w(5心p)口！sfsLhLLEditViewI_nsertToolsDesktopWindowHelp号晏4a□匡］■口2)1型系统：一K(K=10、T=1、T=15)、5(Ts+1)s(Ts+1)(Ts+1)1、2)、112s(s2+2&①s+①2)Transferfunction:10(K=2、&=0.4、®=1),sA2+sGm=InfPm=17.9642Wcg=InfWcp=3.08421015sA3+16sA2+sGm=0.1067Pm=-30.9978Wcg=0.2582

i=iiiiiiii11iii=iiiiiiii11iiiiiii■iiiiii111iiiiiiiiiiiiiiiBodeDiagram__I_J_LLLU1_L_l_L1LLI1iiiiiiiiii■■iiiiii1-"r「「nTt~r~i~rtrmEditViev弓9焰/I_nsertToolsDesktopWin虱鼠急⑥曼贰-dowsHelp□■■口-3-3-1fl1210101010IQ10Frequency(radysec)」」2sA3+0.8sA2+sGm=0.4000Pm=-39.7436Wcg=1r-r-Figure1BodeDiagram11BodeDiagram111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111■i■iiii■|ii■iiii■iiiiiiii11ii■iiii■iiiiiiii11iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii-I—-1—J_UULI_I__-—-