自动控制原理5.2---典型环节与开环系统的频率特性
参考书籍:《自动控制原理》(第七版).胡寿松主编.
《自动控制原理PDF版下载》
2.典型环节与开环系统的频率特性
2.1 典型环节
典型环节分为:最小相位环节、非最小相位环节;
最小相位环节:
- 比例环节:K,(K>0)K,(K>0)K,(K>0);
- 惯性环节:1/(Ts+1),(T>0)1/(Ts+1),(T>0)1/(Ts+1),(T>0);
- 一阶微分环节:Ts+1,(T>0)Ts+1,(T>0)Ts+1,(T>0);
- 振荡环节:1/(s2/ωn2+2ζs/ωn+1),(ωn>0,0≤ζ<1)1/(s^2/\omega_n^2+2\zeta{s}/\omega_n+1),(\omega_n>0,0≤\zeta<1)1/(s2/ωn2+2ζs/ωn+1),(ωn>0,0≤ζ<1);
- 二阶微分环节:s2/ωn2+2ζs/ωn+1,(ωn>0,0≤ζ<1)s^2/\omega_n^2+2\zeta{s}/\omega_n+1,(\omega_n>0,0≤\zeta<1)s2/ωn2+2ζs/ωn+1,(ωn>0,0≤ζ<1);
- 积分环节:1/s1/s1/s;
- 微分环节:sss;
非最小相位环节:
- 比例环节:K,(K<0)K,(K<0)K,(K<0);
- 惯性环节:1/(−Ts+1),(T>0)1/(-Ts+1),(T>0)1/(−Ts+1),(T>0);
- 一阶微分环节:−Ts+1,(T>0)-Ts+1,(T>0)−Ts+1,(T>0);
- 振荡环节:1/(s2/ωn2−2ζs/ωn+1),(ωn>0,0<ζ<1)1/(s^2/\omega_n^2-2\zeta{s}/\omega_n+1),(\omega_n>0,0<\zeta<1)1/(s2/ωn2−2ζs/ωn+1),(ωn>0,0<ζ<1);
- 二阶微分环节:s2/ωn2−2ζs/ωn+1,(ωn>0,0<ζ<1)s^2/\omega_n^2-2\zeta{s}/\omega_n+1,(\omega_n>0,0<\zeta<1)s2/ωn2−2ζs/ωn+1,(ωn>0,0<ζ<1);
开环传递函数的典型环节分解可将开环系统表示为若干个典型环节的串联形式:
G(s)H(s)=∏i=1NGi(s)(1)G(s)H(s)=\prod_{i=1}^NG_i(s)\tag{1} G(s)H(s)=i=1∏NGi(s)(1)
设典型环节的频率特性为:
Gi(jω)=Ai(ω)ejφi(ω)(2)G_i(j\omega)=A_i(\omega)e^{j\varphi_i(\omega)}\tag{2} Gi(jω)=Ai(ω)ejφi(ω)(2)
则系统开环频率特性为:
G(jω)H(jω)=[∏i=1NAi(ω)]ej[∑i=1Nφi(ω)](3)G(j\omega)H(j\omega)=\left[\prod_{i=1}^NA_i(\omega)\right]e^{j\left[\sum_{i=1}^N\varphi_i(\omega)\right]}\tag{3} G(jω)H(jω)=[i=1∏NAi(ω)]ej[∑i=1Nφi(ω)](3)
系统开环幅频特性和开环相频特性为:
A(ω)=∏i=1NAi(ω),φ(ω)=∑i=1Nφi(ω)(4)A(\omega)=\prod_{i=1}^NA_i(\omega),\varphi(\omega)=\sum_{i=1}^N\varphi_i(\omega)\tag{4} A(ω)=i=1∏NAi(ω),φ(ω)=i=1∑Nφi(ω)(4)
系统开环对数幅频特性为:
L(ω)=20lgA(ω)=∑i=1N20lgAi(ω)=∑i=1NLi(ω)(5)L(\omega)=20\lg{A(\omega)}=\sum_{i=1}^N20\lg{A_i(\omega)}=\sum_{i=1}^NL_i(\omega)\tag{5} L(ω)=20lgA(ω)=i=1∑N20lgAi(ω)=i=1∑NLi(ω)(5)
- 系统开环频率特性表现为组成开环系统的诸典型环节频率特性的合成;
- 系统开环对数频率特性,表现为诸典型环节对数频率特性叠加;
2.2 典型环节的频率特性
典型环节频率特性曲线的若干特点:
非最小相位与对应的最小相位环节
最小相位的比例环节G(s)=K(K>0)G(s)=K(K>0)G(s)=K(K>0),简称比例环节,幅频和相频特性如下:
A(ω)=K,φ(ω)=0°(6)A(\omega)=K,\varphi(\omega)=0°\tag{6} A(ω)=K,φ(ω)=0°(6)
非最小相位的比例环节G(s)=−K(K>0)G(s)=-K(K>0)G(s)=−K(K>0),幅频和相频特性如下:
A(ω)=K,φ(ω)=−180°(7)A(\omega)=K,\varphi(\omega)=-180°\tag{7} A(ω)=K,φ(ω)=−180°(7)
最小相位的惯性环节G(s)=1Ts+1,(T>0)G(s)=\frac{1}{Ts+1},(T>0)G(s)=Ts+11,(T>0),幅频和相频特性为:
A(ω)=11+T2ω2,φ(ω)=−arctanTω(8)A(\omega)=\frac{1}{\sqrt{1+T^2\omega^2}},\varphi(\omega)=-\arctan{T\omega}\tag{8} A(ω)=1+T2ω21,φ(ω)=−arctanTω(8)
非最小相位的惯性环节G(s)=1−Ts+1,(T>0)G(s)=\frac{1}{-Ts+1},(T>0)G(s)=−Ts+11,(T>0),幅频和相频特性为:
A(ω)=11+T2ω2,φ(ω)=arctanTω(9)A(\omega)=\frac{1}{\sqrt{1+T^2\omega^2}},\varphi(\omega)=\arctan{T\omega}\tag{9} A(ω)=1+T2ω21,φ(ω)=arctanTω(9)- 最小相位惯性环节和非最小相位的惯性环节,幅频特性相同,相频特性符号相反,幅相曲线关于实轴对称;
- 最小相位惯性环节和非最小相位惯性环节,对数幅频曲线相同,对数相频曲线关于0°线对称;
- 上述两点特点,对于振荡环节和非最小相位振荡环节、一阶微分环节和非最小相位一阶微分环节、二阶微分环节和非最小相位二阶微分环节均适用;
传递函数互为倒数的典型环节
最小相位典型环节中,积分环节和微分环节、惯性环节和一阶微分环节、振荡环节和二阶微分环节的传递函数互为倒数,有如下关系:
G1(s)=1G2(s)(10)G_1(s)=\frac{1}{G_2(s)}\tag{10} G1(s)=G2(s)1(10)
设G1(jω)=A1(ω)ejφ1(ω)G_1(j\omega)=A_1(\omega)e^{j\varphi_1(\omega)}G1(jω)=A1(ω)ejφ1(ω),则有:
{φ2(ω)=−φ1(ω)L2(ω)=20lgA2(ω)=20lg1A1(ω)=−L1(ω)(11)\begin{cases} & \varphi_2(\omega)=-\varphi_1(\omega) \\ & L_2(\omega)=20\lg{A_2(\omega)}=20\lg\frac{1}{A_1(\omega)}=-L_1(\omega) \end{cases}\tag{11} {φ2(ω)=−φ1(ω)L2(ω)=20lgA2(ω)=20lgA1(ω)1=−L1(ω)(11)
传递函数互为倒数的典型环节,对数幅频曲线关于0dB线对称,对数相频曲线关于0°线对称;此结论在非最小相位环节中亦适用;振荡环节和二阶微分环节
振荡环节的传递函数为:
G(s)=1(s/ωn)2+2ζ(s/ωn)+1;ωn>0,0<ζ<1(12)G(s)=\frac{1}{(s/\omega_n)^2+2\zeta(s/\omega_n)+1};\omega_n>0,0<\zeta<1\tag{12} G(s)=(s/ωn)2+2ζ(s/ωn)+11;ωn>0,0<ζ<1(12)
振荡环节的频率特性:
A(ω)=1(1−ω2ωn2)2+4ζ2ω2ωn2(13)A(\omega)=\frac{1}{\sqrt{\begin{pmatrix}1-\frac{\omega^2}{\omega_n^2}\end{pmatrix}^2+4\zeta^2\frac{\omega^2}{\omega_n^2}}}\tag{13} A(ω)=(1−ωn2ω2)2+4ζ2ωn2ω21(13)φ(ω)=−arctan2ζωωn1−ω2ωn2={−arctan2ζωωn1−ω2ωn2,ω≤ωn−(180−arctan2ζωωnω2ωn2−1),ω>ωn(14)\varphi(\omega)=-\arctan\frac{2\zeta\frac{\omega}{\omega_n}}{1-\frac{\omega^2}{\omega_n^2}}= \begin{cases} &-\arctan\frac{2\zeta\frac{\omega}{\omega_n}}{1-\frac{\omega^2}{\omega_n^2}},\omega≤\omega_n \\ &-\begin{pmatrix} 180-\arctan\frac{2\zeta\frac{\omega}{\omega_n}}{\frac{\omega^2}{\omega_n^2}-1} \end{pmatrix},\omega>\omega_n \end{cases}\tag{14} φ(ω)=−arctan1−ωn2ω22ζωnω=⎩⎨⎧−arctan1−ωn2ω22ζωnω,ω≤ωn−(180−arctanωn2ω2−12ζωnω),ω>ωn(14)
φ(0)=0°,φ(∞)=−180°\varphi(0)=0°,\varphi(\infty)=-180°φ(0)=0°,φ(∞)=−180°,相频特性曲线从0°单调减至-180°;当ω=ωn\omega=\omega_nω=ωn时,φ(ωn)=−90°\varphi(\omega_n)=-90°φ(ωn)=−90°,A(ωn)=12ζ,A(\omega_n)=\frac{1}{2\zeta},A(ωn)=2ζ1,振荡环节与虚轴的交点为−j12ζ-j\frac{1}{2\zeta}−j2ζ1;
A(0)=1,A(∞)=0A(0)=1,A(\infty)=0A(0)=1,A(∞)=0,求A(ω)A(\omega)A(ω)的极值,
dA(ω)dω=−[−2ωωn2(1−ω2ωn2)+4ζ2ωωn2][(1−ω2ωn2)2+4ζ2ω2ωn2]32=0(15)\frac{dA(\omega)}{d\omega}=\frac{-\left[-\frac{2\omega}{\omega_n^2}\left(1-\frac{\omega^2}{\omega_n^2}\right)+4\zeta^2\frac{\omega}{\omega_n^2}\right]}{\left[\left(1-\frac{\omega^2}{\omega_n^2}\right)^2+4\zeta^2\frac{\omega^2}{\omega_n^2}\right]^{\frac{3}{2}}}=0\tag{15} dωdA(ω)=[(1−ωn2ω2)2+4ζ2ωn2ω2]23−[−ωn22ω(1−ωn2ω2)+4ζ2ωn2ω]=0(15)
得谐振频率:
ωr=ωn1−2ζ2,0<ζ≤2/2(16)\omega_r=\omega_n\sqrt{1-2\zeta^2},0<\zeta≤\sqrt{2}/2\tag{16} ωr=ωn1−2ζ2,0<ζ≤2/2(16)
谐振峰值:
Mr=A(ωr)=12ζ1−ζ2,0<ζ≤2/2(17)M_r=A(\omega_r)=\frac{1}{2\zeta\sqrt{1-\zeta^2}},0<\zeta≤\sqrt{2}/2\tag{17} Mr=A(ωr)=2ζ1−ζ21,0<ζ≤2/2(17)
当0<ζ<220<\zeta<\frac{\sqrt{2}}{2}0<ζ<22时,
dMrdζ=−(1−2ζ2)ζ2(1−ζ2)32<0(18)\frac{dM_r}{d\zeta}=\frac{-(1-2\zeta^2)}{\zeta^2(1-\zeta^2)^{\frac{3}{2}}}<0\tag{18} dζdMr=ζ2(1−ζ2)23−(1−2ζ2)<0(18)
ωr,Mr\omega_r,M_rωr,Mr均为阻尼比ζ\zetaζ的减函数(0<ζ≤22)(0<\zeta≤\frac{\sqrt{2}}{2})(0<ζ≤22);当0<ζ<220<\zeta<\frac{\sqrt{2}}{2}0<ζ<22时,且ω∈(0,ωr)\omega\in(0,\omega_r)ω∈(0,ωr)时,A(ω)A(\omega)A(ω)单调增;ω∈(ωr,∞)\omega\in(\omega_r,\infty)ω∈(ωr,∞)时,A(ω)A(\omega)A(ω)单调减;当22<ζ<1\frac{\sqrt{2}}{2}<\zeta<122<ζ<1时,A(ω)A(\omega)A(ω)单调减;二阶微分环节的传递函数为振荡环节传递函数的倒数,按对称性可得二阶微分环节的对数频率曲线,有:
{A(0)=1φ(0)=0°,{A(ωn)=2ζφ(ωn)=90°,{A(∞)=∞φ(∞)=180°(19)\begin{cases} &A(0)=1\\ &\varphi(0)=0° \end{cases}, \begin{cases} &A(\omega_n)=2\zeta\\ &\varphi(\omega_n)=90° \end{cases}, \begin{cases} &A(\infty)=\infty\\ &\varphi(\infty)=180° \end{cases}\tag{19} {A(0)=1φ(0)=0°,{A(ωn)=2ζφ(ωn)=90°,{A(∞)=∞φ(∞)=180°(19)
当阻尼比22<ζ<1\frac{\sqrt{2}}{2}<\zeta<122<ζ<1时,A(ω)A(\omega)A(ω)从1单调增至∞\infty∞;当阻尼比0<ζ<220<\zeta<\frac{\sqrt{2}}{2}0<ζ<22,且ω∈(0,ωr)\omega\in(0,\omega_r)ω∈(0,ωr)时,A(ω)A(\omega)A(ω)从1单调减至
{A(ωr)=2ζ1−ζ2<1ωr=ω1−2ζ2(20)\begin{cases} &A(\omega_r)=2\zeta\sqrt{1-\zeta^2}<1\\ &\omega_r=\omega\sqrt{1-2\zeta^2} \end{cases}\tag{20} {A(ωr)=2ζ1−ζ2<1ωr=ω1−2ζ2(20)
在ω∈(ωr,∞)\omega\in(\omega_r,\infty)ω∈(ωr,∞)时,A(ω)A(\omega)A(ω)单调增;对数幅频渐近特性曲线
为简化惯性环节、一阶微分环节、振荡环节和二阶微分环节的对数幅频曲线的作图,常用低频和高频渐近线近似表示对数幅频曲线,称为对数幅频渐近特性曲线;
惯性环节的对数幅频渐近特性为:
La(ω)={0,ω<1T−20lgωT,ω>1T(21)L_a(\omega)= \begin{cases} 0,&\omega<\frac{1}{T}\\ -20\lg\omega{T},&\omega>\frac{1}{T} \end{cases}\tag{21} La(ω)={0,−20lgωT,ω<T1ω>T1(21)
低频部分是零分贝线,高频部分是斜率为−20dB/dec-20dB/dec−20dB/dec的直线,两条直线交于ω=1T\omega=\frac{1}{T}ω=T1处,称频率1T\frac{1}{T}T1为惯性环节的交接频率;用渐近特性近似表示对数幅频特性存在误差:ΔL(ω)=L(ω)−La(ω)\Delta{L(\omega)=L(\omega)-L_a(\omega)}ΔL(ω)=L(ω)−La(ω),在交接频率处误差最大,约为−3dB-3dB−3dB;一阶微分环节和非最小相位一阶微分环节与惯性环节的对数幅频渐近特性曲线以0dB0dB0dB线互为镜像;振荡环节的对数幅频特性为:
L(ω)=−20lg(1−ω2ωn2)2+4ζ2ω2ωn2(22)L(\omega)=-20\lg\sqrt{\left(1-\frac{\omega^2}{\omega_n^2}\right)^2+4\zeta^2\frac{\omega^2}{\omega_n^2}}\tag{22} L(ω)=−20lg(1−ωn2ω2)2+4ζ2ωn2ω2(22)
当ω<<ωn\omega<<\omega_nω<<ωn时,L(ω)≈0L(\omega)≈0L(ω)≈0,低频渐近线为0dB线;当ω>>ωn\omega>>\omega_nω>>ωn时,L(ω)=−40lgωωnL(\omega)=-40\lg\frac{\omega}{\omega_n}L(ω)=−40lgωnω,高频渐近线为过(ωn,0)(\omega_n,0)(ωn,0)点,斜率为-40dB/dec的直线,振荡环节的交接频率为ωn\omega_nωn,对数幅频渐近特性为:
La(ω)={0,ω<ωn−40lgωωn,ω>ωn(23)L_a(\omega)= \begin{cases} 0,&\omega<\omega_n\\ -40\lg\frac{\omega}{\omega_n},&\omega>\omega_n \end{cases}\tag{23} La(ω)={0,−40lgωnω,ω<ωnω>ωn(23)
非最小相位振荡环节与振荡环节对数幅频渐近特性曲线相同,二阶微分环节和非最小相位二阶微分环节与振荡环节的对数幅频渐近特性曲线关于0dB线对称;半对数坐标系中的直线方程为:
k=La(ω2)−La(ω1)lgω2−lgω1(24)k=\frac{L_a(\omega_2)-L_a(\omega_1)}{\lg\omega_2-\lg\omega_1}\tag{24} k=lgω2−lgω1La(ω2)−La(ω1)(24)
其中:[ω1,lg(ω1)]、[ω2,lg(ω2)][\omega_1,\lg(\omega_1)]、[\omega_2,\lg(\omega_2)][ω1,lg(ω1)]、[ω2,lg(ω2)]为直线上的两点,k(dB/dec)k(dB/dec)k(dB/dec)为直线斜率;
2.3 开环幅相特性曲线
绘制概略开环幅相特性曲线的方法:
开环幅相特性曲线的起点(ω=0+)(\omega=0_+)(ω=0+)和终点(ω=∞)(\omega=\infty)(ω=∞);
开环幅相特性曲线与实轴的交点。设ω=ωx\omega=\omega_xω=ωx时,G(jωx)H(jωx)G(j\omega_x)H(j\omega_x)G(jωx)H(jωx)的虚部为:
Im[G(jωx)H(jωx)]=0或φ(ωx)=∠[G(jωx)H(jωx)]=kπ,k=0,±1,±2,…(25)Im[G(j\omega_x)H(j\omega_x)]=0或\varphi(\omega_x)=\angle[G(j\omega_x)H(j\omega_x)]=k\pi,k=0,±1,±2,\dots\tag{25} Im[G(jωx)H(jωx)]=0或φ(ωx)=∠[G(jωx)H(jωx)]=kπ,k=0,±1,±2,…(25)
其中:ωx\omega_xωx称为穿越频率;开环频率特性曲线与实轴交点坐标值为:
Re[G(jωx)H(jωx)]=G(jωx)H(jωx)(26)Re[G(j\omega_x)H(j\omega_x)]=G(j\omega_x)H(j\omega_x)\tag{26} Re[G(jωx)H(jωx)]=G(jωx)H(jωx)(26)开环幅相特性曲线的变化范围(象限、单调性);
实例分析:
Example1: 设系统开环传递函数为:
G(s)H(s)=Ks(T1s+1)(T2s+1);K,T1,T2>0G(s)H(s)=\frac{K}{s(T_1s+1)(T_2s+1)};K,T_1,T_2>0 G(s)H(s)=s(T1s+1)(T2s+1)K;K,T1,T2>0
绘制系统概略开环幅相特性曲线。
解:
系统开环频率特性:
G(jω)H(jω)=K(1−jT1ω)(1−jT2ω)(−j)ω(1+T12ω2)(1+T22ω2)=K[−(T1+T2)ω+j(−1+T1T2ω2)]ω(1+T12ω2)(1+T22ω2)G(j\omega)H(j\omega)=\frac{K(1-jT_1\omega)(1-jT_2\omega)(-j)}{\omega(1+T_1^2\omega^2)(1+T_2^2\omega^2)}=\frac{K\left[-(T_1+T_2)\omega+j(-1+T_1T_2\omega^2)\right]}{\omega(1+T_1^2\omega^2)(1+T_2^2\omega^2)} G(jω)H(jω)=ω(1+T12ω2)(1+T22ω2)K(1−jT1ω)(1−jT2ω)(−j)=ω(1+T12ω2)(1+T22ω2)K[−(T1+T2)ω+j(−1+T1T2ω2)]
幅值变化:A(0+)=∞,A(∞)=0A(0_+)=\infty,A(\infty)=0A(0+)=∞,A(∞)=0;
相角变化:
∠(1jω): −90°~−90°∠(11+jT1ω):0°~−90°∠(11+jT2ω):0°~−90°∠K: 0°~0°φ(ω): −90°~−270°\begin{aligned} &\angle\left(\frac{1}{j\omega}\right):\space\space\space\space\space\space-90°~-90°\\ &\angle\left(\frac{1}{1+jT_1\omega}\right):0°~-90°\\ &\angle\left(\frac{1}{1+jT_2\omega}\right):0°~-90°\\ &\angle{K}:{\space\space\space\space\space\space\space\space\space\space\space\space\space\space\space\space\space\space\space\space}0°~0°\\ &\varphi(\omega):\space\space\space\space\space\space\space\space\space\space\space\space\space\space-90°~-270° \end{aligned} ∠(jω1): −90°~−90°∠(1+jT1ω1):0°~−90°∠(1+jT2ω1):0°~−90°∠K: 0°~0°φ(ω): −90°~−270°
起点处:
Re[G(j0+)H(j0+)]=−K(T1+T2),Im[G(j0+)H(j0+)]=−∞Re[G(j0_+)H(j0_+)]=-K(T_1+T_2),Im[G(j0_+)H(j0_+)]=-\infty Re[G(j0+)H(j0+)]=−K(T1+T2),Im[G(j0+)H(j0+)]=−∞
与实轴交点:
Im[G(jω)H(jω)]=0,⇒ωx=1T1T2Im[G(j\omega)H(j\omega)]=0,\Rightarrow\omega_x=\frac{1}{T_1T_2} Im[G(jω)H(jω)]=0,⇒ωx=T1T21
G(jωx)H(jωx)=Re[G(jωx)H(jωx)]=−KT1T2T1+T2G(j\omega_x)H(j\omega_x)=Re[G(j\omega_x)H(j\omega_x)]=-\frac{KT_1T_2}{T_1+T_2} G(jωx)H(jωx)=Re[G(jωx)H(jωx)]=−T1+T2KT1T2
绘制概略开环幅相特性曲线规律小结:
开环幅相特性曲线的起点,取决于比例环节KKK和系统积分或微分环节的个数ν\nuν(系统型别);
- ν<0\nu<0ν<0,起点为原点;
- ν=0\nu=0ν=0,起点为实轴上的点KKK处,KKK为系统开环增益,KKK有正负之分;
- ν>0\nu>0ν>0,设ν=4k+i(k=0,1,2,…;i=1,2,3,4),\nu=4k+i(k=0,1,2,\dots;i=1,2,3,4),ν=4k+i(k=0,1,2,…;i=1,2,3,4),则K>0K>0K>0时为i×(−90°)i\times(-90°)i×(−90°)的无穷远处,K<0K<0K<0时为i×(−90°)−180°i\times(-90°)-180°i×(−90°)−180°的无穷远处;
开环幅相特性曲线的终点,取决于开环传递函数分子、分母多项式中最小相位环节和非最小相位环节的阶次和;
设系统开环传递函数的分子、分母多项式的阶次分别为mmm和nnn,记除KKK外,分子多项式中最小相位环节的阶次和为m1m_1m1,非最小相位环节的阶次和为m2m_2m2,分母多项式中最小相位环节的阶次和为n1n_1n1,非最小相位环节的阶次和为n2n_2n2,则有:
m=m1+m2,n=n1+n2m=m_1+m_2,n=n_1+n_2 m=m1+m2,n=n1+n2φ(∞)={[(m1−m2)−(n1−n2)]×90°,K>0[(m1−m2)−(n1−n2)]×90°−180°,K<0(27)\varphi(\infty)= \begin{cases} [(m_1-m_2)-(n_1-n_2)]\times90°,&K>0\\ [(m_1-m_2)-(n_1-n_2)]\times90°-180°,&K<0 \end{cases}\tag{27} φ(∞)={[(m1−m2)−(n1−n2)]×90°,[(m1−m2)−(n1−n2)]×90°−180°,K>0K<0(27)
当开环系统为最小相位系统时,若:
n=m,G(j∞)H(j∞)=K∗n>m,Gj∞H(j∞)=0∠[(n−m)×(−90°)](28)\begin{aligned} n=m, {\space\space\space\space\space\space\space\space\space\space}& G(j\infty)H(j\infty)=K^*\\ n>m,{\space\space\space\space\space\space\space\space\space\space}& G{j\infty}H(j\infty)=0\angle[(n-m)\times(-90°)] \end{aligned}\tag{28} n=m, n>m, G(j∞)H(j∞)=K∗Gj∞H(j∞)=0∠[(n−m)×(−90°)](28)
其中:K∗K^*K∗为系统开环根轨迹增益;若开环系统存在等幅振荡环节,重数lll为正整数,即开环传递函数具有如下形式:
G(s)H(s)=1(s2ωn2+1)lG1(s)H1(s)(29)G(s)H(s)=\frac{1}{(\frac{s^2}{\omega_n^2}+1)^l}G_1(s)H_1(s)\tag{29} G(s)H(s)=(ωn2s2+1)l1G1(s)H1(s)(29)
G1(s)H1(s)G_1(s)H_1(s)G1(s)H1(s)不含±jωn±j\omega_n±jωn的极点,则当ω\omegaω趋于ωn\omega_nωn时,A(ω)A(\omega)A(ω)趋于无穷,而:
φ(ωn−)≈φ1(ωn)=∠[G1(jωn)H1(jωn)]φ(ωn+)≈φ1(ωn)−l×180°(30)\begin{aligned} &\varphi(\omega_{n^-})≈\varphi_1(\omega_n)=\angle[G_1(j\omega_n)H_1(j\omega_n)]\\ &\varphi(\omega_{n^+})≈\varphi_1(\omega_n)-l\times180° \end{aligned}\tag{30} φ(ωn−)≈φ1(ωn)=∠[G1(jωn)H1(jωn)]φ(ωn+)≈φ1(ωn)−l×180°(30)
即φ(ω)\varphi(\omega)φ(ω)在ω=ωn\omega=\omega_nω=ωn附近,相角突变−l×180°-l\times180°−l×180°;
2.4 开环对数频率特性曲线
系统开环对数幅频渐近特性:
La(ω)=∑i=1NLai(ω)(31)L_a(\omega)=\sum_{i=1}^NL_{a_i}(\omega)\tag{31} La(ω)=i=1∑NLai(ω)(31)
对于任意的开环传递函数,按典型环节分解,将组成系统的典型环节分为三部分:
- Ksν\displaystyle\frac{K}{s^{\nu}}sνK或−Ksν(K>0)\displaystyle\frac{-K}{s^{\nu}}(K>0)sν−K(K>0);
- 一阶环节,包括惯性环节、一阶微分环节及对应的非最小相位环节,交接频率为1T\displaystyle\frac{1}{T}T1;
- 二阶环节,包括振荡环节、二阶微分环节及对应的非最小相位环节,交接频率为ωn\omega_nωn;
记ωmin\omega_{min}ωmin为最小交接频率,称ω<ωmin\omega<\omega_{min}ω<ωmin的频率范围为低频段;
开环对数幅频渐近特性曲线绘制步骤:
开环传递函数典型环节分解;
确定一阶环节、二阶环节的交接频率,将各交接频率标注在半对数坐标图的ω\omegaω轴上;
绘制低频段渐近特性线:在ω<ωmin\omega<\omega_{min}ω<ωmin频段内,开环系统幅频渐近特性的斜率取决于Kων\frac{K}{\omega^{\nu}}ωνK,因而直线斜率为-20ν\nuνdB/dec;获得低频渐近线,需要确定直线上的一点,方法如下:
- 方法1:在ω<ωmin\omega<\omega_{min}ω<ωmin范围内,任选一点ω0\omega_0ω0,计算La(ω0)=20lgK−20νlgω0L_a(\omega_0)=20\lg{K}-20\nu\lg\omega_0La(ω0)=20lgK−20νlgω0;
- 方法2:取频率为特定值ω0=1\omega_0=1ω0=1,则La(1)=20lgKL_a(1)=20\lg{K}La(1)=20lgK;
- 方法3:取La(ω0)L_a(\omega_0)La(ω0)为特殊值0,有Kω0ν=1\frac{K}{\omega_0^{\nu}}=1ω0νK=1,则ω0=K1ν\omega_0=K^{\frac{1}{\nu}}ω0=Kν1;
过点(ω0,La(ω0))(\omega_0,L_a(\omega_0))(ω0,La(ω0))在ω<ωmin\omega<\omega_{min}ω<ωmin范围内可作斜率为-20ν\nuνdB/dec的直线;
作ω≥ωmin\omega≥\omega_{min}ω≥ωmin频段渐近特性线;在每个交接频率点处,斜率发生变化,变化规律取决于该交接频率对应的典型环节的种类,当系统的多个环节具有相同交接频率时,该交接频率点处斜率的变化应为各个环节对应的斜率变化值的代数和;
交接频率点处斜率的变化表:
实例分析:
Example2: 已知系统开环传递函数为:
G(s)H(s)=2000s−4000s2(s+1)(s2+10s+400)G(s)H(s)=\frac{2000s-4000}{s^2(s+1)(s^2+10s+400)} G(s)H(s)=s2(s+1)(s2+10s+400)2000s−4000
绘制系统开环对数幅频渐近特性曲线。
解:
开环传递函数的典型环节分解形式为:
G(s)H(s)=−10(1−s2)s2(s+1)(s2202+12s20+1)G(s)H(s)=\frac{-10(1-\frac{s}{2})}{s^2(s+1)(\frac{s^2}{20^2}+\frac{1}{2}\frac{s}{20}+1)} G(s)H(s)=s2(s+1)(202s2+2120s+1)−10(1−2s)
开环系统由6个典型环节串联而成:非最小相位比例环节、两个积分环节、非最小相位一阶微分环节、惯性环节和振荡环节;
STEP1:确定各交接频率ωi,i=1,2,3\omega_i,i=1,2,3ωi,i=1,2,3及斜率变化值;
非最小相位一阶微分环节:ω2=2\omega_2=2ω2=2,斜率增加20dB/dec;
惯性环节:ω1=1\omega_1=1ω1=1,斜率减小20dB/dec;
振荡环节:ω3=20\omega_3=20ω3=20,斜率减小40dB/dec;
STEP2:绘制低频段渐近特性曲线;
因为ν=2\nu=2ν=2,低频渐近线斜率k=−40k=-40k=−40dB/dec,直线上一点(ω0,La(ω0))=(1,20dB)(\omega_0,L_a(\omega_0))=(1,20dB)(ω0,La(ω0))=(1,20dB);
STEP3:绘制中高频段渐近特性曲线;
ωmin≤ω<ω2,k=−60dB/decω2≤ω<ω3,k=−40dB/decω≥ω3,k=−80dB/dec\begin{aligned} &\omega_{min}≤\omega<\omega_2,k=-60dB/dec\\ &\omega_2≤\omega<\omega_3,k=-40dB/dec\\ &\omega≥\omega_3,k=-80dB/dec \end{aligned} ωmin≤ω<ω2,k=−60dB/decω2≤ω<ω3,k=−40dB/decω≥ω3,k=−80dB/dec
2.5 延迟环节和延迟系统
输出量经恒定延时后不失真地复现输入量变化的环节称为延迟环节;含有延迟环节的系统称为延迟系统;
延迟环节的输入输出时域表达式为:
c(t)=1(t−τ)r(t−τ)(32)c(t)=1(t-\tau)r(t-\tau)\tag{32} c(t)=1(t−τ)r(t−τ)(32)
其中:τ\tauτ为延迟时间;
延迟环节的传递函数为:
G(s)=C(s)R(s)=e−τs(33)G(s)=\frac{C(s)}{R(s)}=e^{-\tau{s}}\tag{33} G(s)=R(s)C(s)=e−τs(33)
延迟环节的频率特性为:
G(jω)=e−jτω=1⋅∠(−57.3τω)(34)G(j\omega)=e^{-j\tau\omega}=1·\angle(-57.3\tau\omega)\tag{34} G(jω)=e−jτω=1⋅∠(−57.3τω)(34)
2.6 传递函数的频域实验确定
频率响应实验
- 先选择信号源输出的正弦信号的幅值,以使系统处于非饱和状态;
- 在一定频率范围内,改变输入正弦信号的频率,记录各频率点处系统输出信号的波形;
- 由稳态段的输入输出信号的幅值比和相位差绘制对数频率特性曲线;
传递函数确定
从低频段起,将实验所得的对数幅频特性曲线用斜率为0dB/dec,±20dB/dec,±40dB/dec,…直线分段近似,获得对数幅频渐近特性曲线;
由对数幅频渐近特性曲线确定最小相位条件下系统的传递函数,是对数幅频渐近特性曲线绘制的逆问题;
实例分析:
Example3: 下图由频率响应实验获得的某最小相位系统的开环对数幅频特性曲线和对数幅频渐近特性曲线,确定系统的传递函数。
解:
STEP1:确定系统积分或微分环节的个数。
由图可知,对数幅频渐近特性曲线低频渐近线的斜率为+20dB/dec,因此ν=−1\nu=-1ν=−1,即系统含有一个微分环节;
STEP2:确定系统传递函数结构形式。
本题中,有两个交接频率:
ω=ω1\omega=\omega_1ω=ω1处,斜率变化-20dB/dec,对应惯性环节;
ω=ω2\omega=\omega_2ω=ω2处,斜率变化-40dB/dec,可以对应振荡环节,也可以是重惯性环节,由图可知,在ω2\omega_2ω2附近存在谐振现象,因此为振荡环节;
系统传递函数形式为:
G(s)=Ks(1+sω1)(s2ω22+2ζsω2+1)G(s)=\frac{Ks}{\left(1+\displaystyle\frac{s}{\omega_1}\right)\left(\displaystyle\frac{s^2}{\omega_2^2}+2\zeta\frac{s}{\omega_2}+1\right)} G(s)=(1+ω1s)(ω22s2+2ζω2s+1)Ks
其中参数ω1、ω2、ζ、K\omega_1、\omega_2、\zeta、Kω1、ω2、ζ、K待定;
STEP3:由给定条件确定传递函数参数。
低频渐近线方程为:
La(ω)=20lgKων=20lgK−20νlgωL_a(\omega)=20\lg\frac{K}{\omega^{\nu}}=20\lg{K}-20\nu\lg\omega La(ω)=20lgωνK=20lgK−20νlgω
由给定点(ω,La(ω))=(1,0)(\omega,L_a(\omega))=(1,0)(ω,La(ω))=(1,0)及ν=−1\nu=-1ν=−1可得:K=1K=1K=1
根据直线方程式:
La(ωa)−La(ωb)=k(lgωa−lgωb)L_a(\omega_a)-L_a(\omega_b)=k(\lg\omega_a-\lg\omega_b) La(ωa)−La(ωb)=k(lgωa−lgωb)
由给定点:
ωa=1,La(ωa)=0,ωb=ω1,La(ωb)=12,k=20\omega_a=1,L_a(\omega_a)=0,\omega_b=\omega_1,L_a(\omega_b)=12,k=20 ωa=1,La(ωa)=0,ωb=ω1,La(ωb)=12,k=20
可得:
ω1=101220=3.98\omega_1=10^{\frac{12}{20}}=3.98 ω1=102012=3.98
由给顶点:
ωa=100,La(ωa)=0,ωb=ω2,La(ωb)=12,k=−40\omega_a=100,L_a(\omega_a)=0,\omega_b=\omega_2,L_a(\omega_b)=12,k=-40 ωa=100,La(ωa)=0,ωb=ω2,La(ωb)=12,k=−40
可得:
ω2=10(−1240+lg100)=50.1\omega_2=10^{\left(-\frac{12}{40}+\lg100\right)}=50.1 ω2=10(−4012+lg100)=50.1
在谐振频率ωr\omega_rωr处,振荡环节的谐振峰值为:
20lgMr=20lg12ζ1−ζ220\lg{M_r}=20\lg\frac{1}{2\zeta\sqrt{1-\zeta^2}} 20lgMr=20lg2ζ1−ζ21
根据:20lgMr=20−12=8(dB)20\lg{M_r}=20-12=8(dB)20lgMr=20−12=8(dB),因此:Mr=2.512M_r=2.512Mr=2.512,则有:
ζ4−ζ2+0.04=0⇒ζ1=0.204,ζ2=0.979\zeta^4-\zeta^2+0.04=0\Rightarrow\zeta_1=0.204,\zeta_2=0.979 ζ4−ζ2+0.04=0⇒ζ1=0.204,ζ2=0.979
因为0<ζ<0.7070<\zeta<0.7070<ζ<0.707才存在谐振峰值,因此ζ=0.204\zeta=0.204ζ=0.204;
综上,系统的传递函数为:
G(s)=s(s3.98+1)(s250.12+0.408s50.1)+1G(s)=\frac{s}{\left(\displaystyle\frac{s}{3.98}+1\right)\left(\displaystyle\frac{s^2}{50.1^2}+0.408\frac{s}{50.1}\right)+1} G(s)=(3.98s+1)(50.12s2+0.40850.1s)+1s
自动控制原理5.2---典型环节与开环系统的频率特性相关推荐
- matlab振荡环节的频率特性,5-2典型环节与开环系统的频率特性ppt2010.ppt
<5-2典型环节与开环系统的频率特性ppt2010.ppt>由会员分享,可在线阅读,更多相关<5-2典型环节与开环系统的频率特性ppt2010.ppt(43页珍藏版)>请在人人 ...
- 自动控制原理实验一 典型环节及其阶跃响应
实验内容 构成下述典型一阶系统的模拟电路,并测量其阶跃响应: 比例环节的模拟电路及其传递函数如图1-1. G(S)= -R2/R1 2.惯性环节的模拟电路及其传递函数如图1-2. G(S)= - ...
- 自动控制原理上课笔记(不定期更新)
文章目录 前言 自动控制的一般概念 经典控制理论和现代控制理论的区别 自动控制的基本原理与方式 自动控制系统基本控制方式 对自动控制系统的基本要求 控制系统的数学模型 数学模型 微分方程 小偏差法 传 ...
- MATLAB rolcus函数,自动控制原理实验报告 .doc
<自动控制原理实验报告 .doc>由会员分享,可在线阅读,更多相关<自动控制原理实验报告 .doc(32页珍藏版)>请在装配图网上搜索. 1. 实验报告课程名称 自动控制原理 ...
- 【控制】《自动控制原理》胡寿松老师-第5章-线性系统的频域分析法
第4章 回到目录 第6章 第5章-线性系统的频域分析法 5.1 频率特性 5.2 典型环节与开环系统的频率特性 5.3 频率域稳定判据 5.3.1 奈氏判据的数学基础 5.3.2 奈奎斯特稳定判据(奈 ...
- matlab振荡环节的频率特性,自动控制原理_典型环节的频率特性教材教学课件.ppt...
演示文稿演讲PPT学习教学课件医学文件教学培训课件 5.3 控制系统开环频率特性曲线的绘制 控制系统开环频率特性的典型环节分解 开环幅相特性曲线的绘制(Nyquist图) 开环对数频率特性曲线的绘制( ...
- 自动控制原理(4)——传递函数、典型环节的传递函数
自动控制原理(4)--传递函数.典型环节的传递函数 微分方程模型 优点:是时间域的数学模型,比较直观,它用时间域的方式,描述系统输入和输出变量之间的关系 在给定初始条件和输入信号后,借助计算机可以 ...
- 典型环节的matlab仿真实验报告,自动控制原理实验报告
<自动控制原理实验报告>由会员分享,可在线阅读,更多相关<自动控制原理实验报告(39页珍藏版)>请在人人文库网上搜索. 1.自动控制原理实验报告学院:电子信息工程学院班级:自动 ...
- 计算机控制实验比例环节,自动控制原理实验教学课件.ppt
<自动控制原理实验教学课件.ppt>由会员分享,可在线阅读,更多相关<自动控制原理实验教学课件.ppt(45页珍藏版)>请在装配图网上搜索. 1.自动控制原理实验教学课件 信息 ...
最新文章
- 中文版-He Knows My Name(他知道我的名字)-祂认得我-陈熙(音乐河4)
- Framebuffer的配置及应用——先转载留着,以后一定要弄懂
- FlowDroid工具的构建与运行
- python3 for sum_Python for循环和“sum13”方法
- js 字符串截取 获取固定标识字段
- 独立博客网站FansUnion.cn运营2年的经验和教训以及未来规划
- [詹兴致矩阵论习题参考解答]习题1.10
- Finding distance between two curves
- 抖音python上的代码_抖音代码舞「图文推荐」,python实例代码
- 关于新加坡的身份证与电话号码验证
- Mac访问NTFS文件系统的移动硬盘
- chapter02作业
- H3C防火墙基础配置1-登录配置、安全域配置
- Visual C#程序设计基础pdf
- android设计模式面试,从外包公司到今日头条offer,吐血整理
- 体检报告录入有误,到底是谁的错?
- php去掉二维数组中某key的值
- (C++实例)实现people类、student类,teacher类、graduate类、助教类继承和派生并测试
- mysql inet aton ipv6_Linux网络编程IPv4和IPv6的inet_addr、inet_aton、inet_pton等函数小结
- 2022年了,密码该如何保存都不会?
热门文章
- 触景无限CEO肖洪波:冗杂的AI安防行业,需追寻“精益之道”
- js获取当前日期凌晨开始到当前日期24点的时间,转换时间戳
- 安德普泰靠什么跻身“毕马威新国货50”榜单?
- 数字孪生钢铁行业研究案例
- 丝印代码html,丝印代码反查型号:ic上的丝印为HH1YA,封装是SOT23-6的,请问是什么型号呢?...
- 机器人战争显示服务器连接,飞机机器人战争城市救援
- 【macos vs2022】记录一次.NET6 webapi 证书问题
- Java--关于ROUND_HALF_DOEN五舍六入
- ctfshow phps源码泄露
- 作为投资人,我的两次投资失败经历