该系列博客主要讲述Matlab软件在自动控制方面的应用,如无自动控制理论基础,请先学习自动控制系列博文,该系列博客不再详细讲解自动控制理论知识。
自动控制理论基础相关链接:https://blog.csdn.net/qq_39032096/category_10287468.html?spm=1001.2014.3001.5482
博客参考书籍:《MATLAB/Simulink与控制系统仿真》。
3.描述函数法
3.1 描述函数的基本概念
3.1.1 描述函数的定义
设非线性环节输入输出描述为:
y = f ( x ) y=f(x) y=f(x)
当非线性环节的输入信号为正弦信号:
x ( t ) = A sin ω t x(t)=A\sin\omega{t} x(t)=Asinωt
对非线性环节的稳态输出 y ( t ) y(t) y(t)进行谐波分析;
一般情况下, y ( t ) y(t) y(t)为非正弦的周期信号,可以展开成傅里叶级数:
y ( t ) = A 0 + ∑ n = 1 ∞ ( A n cos n ω t + B n sin n ω t ) = A 0 + ∑ n = 1 ∞ Y n sin ( n ω t + φ n ) y(t)=A_0+\sum_{n=1}^{\infty}(A_n\cos{n}\omega{t}+B_n\sin{n}\omega{t})=A_0+\sum_{n=1}^{\infty}Y_n\sin(n\omega{t}+\varphi_n) y(t)=A0+n=1∑∞(Ancosnωt+Bnsinnωt)=A0+n=1∑∞Ynsin(nωt+φn)
其中: A 0 A_0 A0为直流分量; Y n sin ( n ω t + φ n ) Y_n\sin(n\omega{t}+\varphi_n) Ynsin(nωt+φn)为第 n n n次谐波分量,且有:
Y n = A n 2 + B n 2 , φ n = arctan A n B n Y_n=\sqrt{A_n^2+B_n^2},\varphi_n=\arctan\frac{A_n}{B_n} Yn=An2+Bn2,φn=arctanBnAn
其中: A n 、 B n A_n、B_n An、Bn为傅里叶系数;
A n = 1 π ∫ 0 2 π y ( t ) cos n ω t d ( ω t ) , B n = 1 π ∫ 0 2 π y ( t ) sin n ω t d ( ω t ) ( n = 1 , 2 , … , ) A_n=\frac{1}{\pi}\int_{0}^{2\pi}y(t)\cos{n}\omega{t}{\rm d}(\omega{t}),B_n=\frac{1}{\pi}\int_{0}^{2\pi}y(t)\sin{n}\omega{t}{\rm d}(\omega{t})(n=1,2,\dots,) An=π1∫02πy(t)cosnωtd(ωt),Bn=π1∫02πy(t)sinnωtd(ωt)(n=1,2,…,)
直流分量:
A 0 = 1 2 π ∫ 0 2 π y ( t ) d ( ω t ) A_0=\frac{1}{2\pi}\int_{0}^{2\pi}y(t){\rm d}(\omega{t}) A0=2π1∫02πy(t)d(ωt)
若 A 0 = 0 A_0=0 A0=0且当 n > 1 n>1 n>1时, Y n Y_n Yn均很小,则可近似认为非线性环节的正弦响应应仅有基波分量:
y ( t ) ≈ A 1 cos ω t + B 1 sin ω t = Y 1 sin ( ω t + φ 1 ) y(t)≈A_1\cos\omega{t}+B_1\sin\omega{t}=Y_1\sin(\omega{t}+\varphi_1) y(t)≈A1cosωt+B1sinωt=Y1sin(ωt+φ1)
定义正弦输入信号作用下,非线性环节的稳态输出中基波分量和输入信号的复数比为非线性环节的描述函数,用 N ( A ) N(A) N(A)表示,即:
N ( A ) = ∣ N ( A ) ∣ e j ∠ N ( A ) = Y 1 A e j φ 1 = B 1 + j A 1 A N(A)=|N(A)|{\rm e}^{j\angle{N(A)}}=\frac{Y_1}{A}{\rm e}^{j\varphi_1}=\frac{B_1+{\rm j}A_1}{A} N(A)=∣N(A)∣ej∠N(A)=AY1ejφ1=AB1+jA1
3.1.2 非线性系统描述函数法分析的应用条件
-
非线性系统应简化成一个非线性环节和一个线性部分闭环连接的典型结构图,如下图所示:
-
非线性环节的输入输出特性 y ( x ) y(x) y(x)应是 x x x的奇函数,即 f ( x ) = − f ( − x ) f(x)=-f(-x) f(x)=−f(−x),或正弦输入下的输出为 t t t的奇对称函数,即 y ( t + π ω ) = − y ( t ) y\left(t+\displaystyle\frac{\pi}{\omega}\right)=-y(t) y(t+ωπ)=−y(t),以保证非线性环节的正弦响应不含有常值分量,即 A 0 = 0 A_0=0 A0=0;
-
系统的线性部分应具有较好的低通滤波性能;线性部分的阶次越高,低通滤波性能越好;欲具有低通滤波性能,线性部分的极点应位于复平面的左半平面;
3.2 典型非线性特性的描述函数
3.2.1 死区饱和非线性环节
输出 y ( t ) y(t) y(t)的数学表达式为:
y ( t ) = { 0 , 0 ≤ ω t ≤ Ψ 1 K ( A sin ω t − Δ ) , Ψ 1 < ω t ≤ Ψ 2 K ( a − Δ ) , Ψ 2 < ω t ≤ π 2 y(t)= \begin{cases} 0,&0≤\omega{t}≤\Psi_1\\\\ K(A\sin\omega{t}-\Delta),&\Psi_1<\omega{t}≤\Psi_2\\\\ K(a-\Delta),&\Psi_2<\omega{t}≤\displaystyle\frac{\pi}{2} \end{cases} y(t)=⎩
⎨
⎧0,K(Asinωt−Δ),K(a−Δ),0≤ωt≤Ψ1Ψ1<ωt≤Ψ2Ψ2<ωt≤2π
确定 y ( t ) y(t) y(t)产生不同线性变换的区间端点为:
Ψ 1 = arcsin Δ A , Ψ 2 = arcsin a A \Psi_1=\arcsin\frac{\Delta}{A}, \Psi_2=\arcsin\frac{a}{A} Ψ1=arcsinAΔ,Ψ2=arcsinAa
由于 y ( t ) y(t) y(t)为奇函数,因此: A 0 = 0 , A 1 = 0 A_0=0,A_1=0 A0=0,A1=0,又因为 y ( t ) y(t) y(t)为半周期内对称,因此:
B 1 = 1 π ∫ 0 2 π y ( t ) sin ω t d ( ω t ) = 4 π ∫ 0 π 2 y ( t ) sin ω t d ( ω t ) = 4 π [ ∫ Ψ 1 Ψ 2 [ K ( A sin ω t − Δ ) ] sin ω t d ( ω t ) + ∫ Ψ 2 π 2 K ( a − Δ ) sin ω t d ( ω t ) ] = 2 K A π [ arcsin a A − arcsin Δ A + a A 1 − ( a A ) 2 − Δ A 1 − ( Δ A ) 2 ] \begin{aligned} B_1&=\frac{1}{\pi}\int_{0}^{2\pi}y(t)\sin\omega{t}{\rm d}(\omega{t})=\frac{4}{\pi}\int_0^{\frac{\pi}{2}}y(t)\sin\omega{t}{\rm d}(\omega{t})\\\\ &=\frac{4}{\pi}\left[\int_{\Psi_1}^{\Psi_2}[K(A\sin\omega{t}-\Delta)]\sin\omega{t}{\rm d}(\omega{t})+\int_{\Psi_2}^{\frac{\pi}{2}}K(a-\Delta)\sin\omega{t}{\rm d}(\omega{t})\right]\\\\ &=\frac{2KA}{\pi}\left[\arcsin\frac{a}{A}-\arcsin\frac{\Delta}{A}+\frac{a}{A}\sqrt{1-\left(\frac{a}{A}\right)^2}-\frac{\Delta}{A}\sqrt{1-\left(\frac{\Delta}{A}\right)^2}\right] \end{aligned} B1=π1∫02πy(t)sinωtd(ωt)=π4∫02πy(t)sinωtd(ωt)=π4[∫Ψ1Ψ2[K(Asinωt−Δ)]sinωtd(ωt)+∫Ψ22πK(a−Δ)sinωtd(ωt)]=π2KA
arcsinAa−arcsinAΔ+Aa1−(Aa)2−AΔ1−(AΔ)2
死区饱和特性的描述函数为:
N ( A ) = 2 K π [ arcsin a A − arcsin Δ A + a A 1 − ( a A ) 2 − Δ A 1 − ( Δ A ) 2 ] , A ≥ a N(A)=\frac{2K}{\pi}\left[\arcsin\frac{a}{A}-\arcsin\frac{\Delta}{A}+\frac{a}{A}\sqrt{1-\left(\frac{a}{A}\right)^2}-\frac{\Delta}{A}\sqrt{1-\left(\frac{\Delta}{A}\right)^2}\right],A≥a N(A)=π2K
arcsinAa−arcsinAΔ+Aa1−(Aa)2−AΔ1−(AΔ)2
,A≥a
取 Δ = 0 \Delta=0 Δ=0时,饱和特性的描述函数为:
N ( A ) = 2 K π [ arcsin a A + a A 1 − ( a A ) 2 ] , A ≥ a N(A)=\frac{2K}{\pi}\left[\arcsin\frac{a}{A}+\frac{a}{A}\sqrt{1-\left(\frac{a}{A}\right)^2}\right],A≥a N(A)=π2K[arcsinAa+Aa1−(Aa)2],A≥a
死区特性描述函数:
N ( A ) = 2 K π [ π 2 − arcsin Δ A − Δ A 1 − ( Δ A ) 2 ] , A ≥ Δ N(A)=\frac{2K}{\pi}\left[\frac{\pi}{2}-\arcsin\frac{\Delta}{A}-\frac{\Delta}{A}\sqrt{1-\left(\frac{\Delta}{A}\right)^2}\right],A≥\Delta N(A)=π2K
2π−arcsinAΔ−AΔ1−(AΔ)2
,A≥Δ
3.2.2 死区与滞环继电特性非线性环节
输出 y ( t ) y(t) y(t)的数学表达式:
y ( t ) = { 0 , 0 ≤ ω t < Ψ 1 M , Ψ 1 ≤ ω t ≤ Ψ 2 0 , Ψ 2 < ω t ≤ π y(t)= \begin{cases} 0,&0≤\omega{t}<\Psi_1\\\\ M,&\Psi_1≤\omega{t}≤\Psi_2\\\\ 0,&\Psi_2<\omega{t}≤\pi \end{cases} y(t)=⎩
⎨
⎧0,M,0,0≤ωt<Ψ1Ψ1≤ωt≤Ψ2Ψ2<ωt≤π
区间端点:
Ψ 1 = arcsin h A , Ψ 2 = π − arcsin m h A \Psi_1=\arcsin\frac{h}{A},\Psi_2=\pi-\arcsin\frac{mh}{A} Ψ1=arcsinAh,Ψ2=π−arcsinAmh
y ( t ) y(t) y(t)为奇对称函数,而非奇函数,则有:
A 1 = 2 π ∫ 0 π y ( t ) cos ω t d ( ω t ) = 2 π ∫ Ψ 1 Ψ 2 M cos ω t d ( ω t ) = 2 M h π A ( m − 1 ) B 1 = 2 π ∫ 0 π y ( t ) sin ω t d ( ω t ) = 2 π ∫ Ψ 1 Ψ 2 M sin ω t d ( ω t ) = 2 M π [ 1 − ( m h A ) 2 + 1 − ( h A ) 2 ] \begin{aligned} &A_1=\frac{2}{\pi}\int_0^{\pi}y(t)\cos\omega{t}{\rm d}(\omega{t})=\frac{2}{\pi}\int_{\Psi_1}^{\Psi_2}M\cos\omega{t}{\rm d}(\omega{t})=\frac{2Mh}{\pi{A}}(m-1)\\\\ &B_1=\frac{2}{\pi}\int_0^{\pi}y(t)\sin\omega{t}{\rm d}(\omega{t})=\frac{2}{\pi}\int_{\Psi_1}^{\Psi_2}M\sin\omega{t}{\rm d}(\omega{t})=\frac{2M}{\pi}\left[\sqrt{1-\left(\frac{mh}{A}\right)^2}+\sqrt{1-\left(\frac{h}{A}\right)^2}\right] \end{aligned} A1=π2∫0πy(t)cosωtd(ωt)=π2∫Ψ1Ψ2Mcosωtd(ωt)=πA2Mh(m−1)B1=π2∫0πy(t)sinωtd(ωt)=π2∫Ψ1Ψ2Msinωtd(ωt)=π2M
1−(Amh)2+1−(Ah)2
因此,死区滞环继电特性描述函数为:
N ( A ) = 2 M π A [ 1 − ( m h A ) 2 + 1 − ( h A ) 2 ] + j 2 M h π A 2 ( m − 1 ) , A ≥ h N(A)=\frac{2M}{\pi{A}}\left[\sqrt{1-\left(\frac{mh}{A}\right)^2}+\sqrt{1-\left(\frac{h}{A}\right)^2}\right]+{\rm j}\frac{2Mh}{\pi{A^2}}(m-1),A≥h N(A)=πA2M
1−(Amh)2+1−(Ah)2
+jπA22Mh(m−1),A≥h
取 h = 0 h=0 h=0,理想继电特性描述函数:
N ( A ) = 4 M π A N(A)=\frac{4M}{\pi{A}} N(A)=πA4M
取 m = 1 m=1 m=1,死区继电特性描述函数:
N ( A ) = 4 M π A 1 − ( h A ) 2 , A ≥ h N(A)=\frac{4M}{\pi{A}}\sqrt{1-\left(\frac{h}{A}\right)^2},A≥h N(A)=πA4M1−(Ah)2,A≥h
取 m = − 1 m=-1 m=−1,滞环继电特性描述函数:
N ( A ) = 4 M π A 1 − ( h A ) 2 − j 4 M h π A 2 , A ≥ h N(A)=\frac{4M}{\pi{A}}\sqrt{1-\left(\frac{h}{A}\right)^2}-{\rm j}\frac{4Mh}{\pi{A^2}},A≥h N(A)=πA4M1−(Ah)2−jπA24Mh,A≥h
3.2.3 典型非线性特性的描述函数
-
理想继电特性(库仑摩擦)和有死区的继电特性
-
理想继电特性:
N ( A ) = 4 M π A N(A)=\frac{4M}{\pi{A}} N(A)=πA4M -
有死区的继电特性:
N ( A ) = 4 M π A 1 − ( h A ) 2 , A ≥ h N(A)=\frac{4M}{\pi{A}}\sqrt{1-\left(\frac{h}{A}\right)^2},A≥h N(A)=πA4M1−(Ah)2,A≥h
-
-
有滞环的继电特性和有死区与滞环的继电特性
-
有滞环的继电特性:
N ( A ) = 4 M π A 1 − ( h A ) 2 − j 4 M h π A 2 N(A)=\frac{4M}{\pi{A}}\sqrt{1-\left(\frac{h}{A}\right)^2}-{\rm j}\frac{4Mh}{\pi{A^2}} N(A)=πA4M1−(Ah)2−jπA24Mh -
有死区与滞环的继电特性:
N ( A ) = 2 M π A [ 1 − ( m h A ) 2 + 1 − ( h A ) 2 ] + j 2 M h π A 2 ( m − 1 ) , A ≥ h N(A)=\frac{2M}{\pi{A}}\left[\sqrt{1-\left(\frac{mh}{A}\right)^2}+\sqrt{1-\left(\frac{h}{A}\right)^2}\right]+{\rm j}\frac{2Mh}{\pi{A^2}}(m-1),A≥h N(A)=πA2M 1−(Amh)2+1−(Ah)2 +jπA22Mh(m−1),A≥h
-
-
饱和特性(幅值限制)和有死区的饱和特性
-
饱和特性:
N ( A ) = 2 K π [ arcsin a A + a A 1 − ( a A ) 2 ] , A ≥ a N(A)=\frac{2K}{\pi}\left[\arcsin\frac{a}{A}+\frac{a}{A}\sqrt{1-\left(\frac{a}{A}\right)^2}\right],A≥a N(A)=π2K[arcsinAa+Aa1−(Aa)2],A≥a -
有死区的饱和特性:
N ( A ) = 2 K π [ arcsin a A − arcsin Δ A + a A 1 − ( a A ) 2 − Δ A 1 − ( Δ A ) 2 ] , A ≥ a N(A)=\frac{2K}{\pi}\left[\arcsin\frac{a}{A}-\arcsin\frac{\Delta}{A}+\frac{a}{A}\sqrt{1-\left(\frac{a}{A}\right)^2}-\frac{\Delta}{A}\sqrt{1-\left(\frac{\Delta}{A}\right)^2}\right],A≥a N(A)=π2K arcsinAa−arcsinAΔ+Aa1−(Aa)2−AΔ1−(AΔ)2 ,A≥a
-
-
死区特性和间隙特性
-
死区特性:
N ( A ) = 2 K π [ π 2 − arcsin Δ A − Δ A 1 − ( Δ A ) 2 ] , A ≥ Δ N(A)=\frac{2K}{\pi}\left[\frac{\pi}{2}-\arcsin\frac{\Delta}{A}-\frac{\Delta}{A}\sqrt{1-\left(\frac{\Delta}{A}\right)^2}\right],A≥\Delta N(A)=π2K 2π−arcsinAΔ−AΔ1−(AΔ)2 ,A≥Δ -
间隙特性:
N ( A ) = K π [ π 2 + arcsin ( 1 − 2 b A ) + 2 ( 1 − 2 b A ) b A ( 1 − b A ) ] + j 4 K b π A ( b A − 1 ) , A ≥ b N(A)=\frac{K}{\pi}\left[\frac{\pi}{2}+\arcsin\left(1-\frac{2b}{A}\right)+2\left(1-\frac{2b}{A}\right)\sqrt{\frac{b}{A}\left(1-\frac{b}{A}\right)}\right]+{\rm j}\frac{4Kb}{\pi{A}}\left(\frac{b}{A}-1\right),A≥b N(A)=πK[2π+arcsin(1−A2b)+2(1−A2b)Ab(1−Ab)]+jπA4Kb(Ab−1),A≥b
-
-
变增益特性、有死区的线性特性、库仑摩擦加黏性摩擦
-
变增益特性:
N ( A ) = K 2 + 2 ( K 1 − K 2 ) π [ arcsin s A + s A 1 − ( s A ) 2 ] , A ≥ s N(A)=K_2+\frac{2(K_1-K_2)}{\pi}\left[\arcsin\frac{s}{A}+\frac{s}{A}\sqrt{1-\left(\frac{s}{A}\right)^2}\right],A≥s N(A)=K2+π2(K1−K2)[arcsinAs+As1−(As)2],A≥s -
有死区的线性特性:
N ( A ) = K − 2 K π arcsin Δ A + 4 M − 2 K Δ π A 1 − ( Δ A ) 2 , A ≥ Δ N(A)=K-\frac{2K}{\pi}\arcsin\frac{\Delta}{A}+\frac{4M-2K\Delta}{\pi{A}}\sqrt{1-\left(\frac{\Delta}{A}\right)^2},A≥\Delta N(A)=K−π2KarcsinAΔ+πA4M−2KΔ1−(AΔ)2,A≥Δ -
库仑摩擦加黏性摩擦:
N ( A ) = K + 4 M π A N(A)=K+\frac{4M}{\pi{A}} N(A)=K+πA4M
-
3.3 非线性系统简化
3.3.1 非线性特性的并联
若两个非线性特性输入相同,输出相加、减,则等效非线性为两个非线性特性的叠加,并联等效非线性特性的描述函数为各非线性描述函数的代数和;
3.3.2 非线性特性的串联
当两个非线性环节串联时,采用图解法简化。
参数确定:
Δ 2 = K 1 ( Δ − Δ 1 ) ⇒ Δ = Δ 1 + Δ 2 K 1 \Delta_2=K_1(\Delta-\Delta_1)\Rightarrow\Delta=\Delta_1+\frac{\Delta_2}{K_1} Δ2=K1(Δ−Δ1)⇒Δ=Δ1+K1Δ2
a 2 = K 1 ( a − Δ 1 ) ⇒ a = a 2 K 1 + Δ 1 a_2=K_1(a-\Delta_1)\Rightarrow{a}=\frac{a_2}{K_1}+\Delta_1 a2=K1(a−Δ1)⇒a=K1a2+Δ1
当 ∣ x ∣ ≤ Δ |x|≤\Delta ∣x∣≤Δ时, y ( x 1 ) y(x_1) y(x1)特性可知, y ( x ) = 0 y(x)=0 y(x)=0;
当 ∣ x ∣ ≥ a |x|≥a ∣x∣≥a时, y ( x 1 ) y(x_1) y(x1)特性可知, y ( x ) = K 2 ( a 2 − Δ 2 ) y(x)=K_2(a_2-\Delta_2) y(x)=K2(a2−Δ2);
当 Δ < ∣ x ∣ < a \Delta<|x|<a Δ<∣x∣<a时, y ( x 1 ) y(x_1) y(x1)位于线性区, y ( x ) y(x) y(x)呈线性,设斜率为 K K K,有:
y ( x ) = K ( x − Δ ) = K 2 ( x 1 − Δ 2 ) y(x)=K(x-\Delta)=K_2(x_1-\Delta_2) y(x)=K(x−Δ)=K2(x1−Δ2)
当 x = a x=a x=a时, x 1 = a 2 x_1=a_2 x1=a2,由于: x 1 = Δ 2 + K 1 ( a − Δ ) x_1=\Delta_2+K_1(a-\Delta) x1=Δ2+K1(a−Δ),因此: a − Δ = a 2 − Δ 2 K 1 a-\Delta=\displaystyle\frac{a_2-\Delta_2}{K_1} a−Δ=K1a2−Δ2,因此有: K = K 1 K 2 K=K_1K_2 K=K1K2;
两个非线性环节串联,等效特性还取决于前后次序;
3.3.3 线性部分的等效变换
3.4 非线性系统稳定性分析的描述函数法
3.4.1 变增益线性系统的稳定性分析
设 G ( s ) G(s) G(s)的极点均位于 s s s的左半平面,即 P = 0 , G ( j ω ) P=0,G(j\omega) P=0,G(jω)的奈奎斯特曲线 Γ G \Gamma_G ΓG如上图所示,闭环系统特征方程为:
1 + K G ( j ω ) = 0 或 G ( j ω ) = − 1 K + j 0 1+KG({\rm j}\omega)=0或G({\rm j}\omega)=-\frac{1}{K}+{\rm j}0 1+KG(jω)=0或G(jω)=−K1+j0
当 Γ G \Gamma_G ΓG曲线不包围点 ( − 1 K , j 0 ) (-\displaystyle\frac{1}{K},{\rm j}0) (−K1,j0)时,即: Z = P − 2 N = − 2 N = 0 Z=P-2N=-2N=0 Z=P−2N=−2N=0,系统闭环稳定;
当 Γ G \Gamma_G ΓG曲线包围点 ( − 1 K , j 0 ) (-\displaystyle\frac{1}{K},{\rm j}0) (−K1,j0)时,系统不稳定;
当 Γ G \Gamma_G ΓG曲线穿越点 ( − 1 K , j 0 ) (-\displaystyle\frac{1}{K},{\rm j}0) (−K1,j0)时,系统临界稳定,将产生等幅振荡;
若设 K K K在一定范围内可变,即有: K 1 ≤ K ≤ K 2 K_1≤K≤K_2 K1≤K≤K2,则 ( − 1 K , j 0 ) (-\displaystyle\frac{1}{K},{\rm j}0) (−K1,j0)为复平面实轴上的一段直线,若 Γ G \Gamma_G ΓG曲线不包围该直线,则系统闭环稳定,当 Γ G \Gamma_G ΓG包围该直线时,则系统闭环不稳定;
3.4.2 应用描述函数分析非线性系统的稳定性
当非线性特性采用描述函数近似等效时,闭环系统的特征方程为:
1 + N ( A ) G ( j ω ) = 0 ⇒ G ( j ω ) = − 1 N ( A ) 1+N(A)G({\rm j}\omega)=0\Rightarrow{G({\rm j}\omega)}=-\frac{1}{N(A)} 1+N(A)G(jω)=0⇒G(jω)=−N(A)1
称 − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1为非线性环节的负倒描述函数;
在复平面上绘制 Γ G \Gamma_G ΓG曲线和 − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1曲线时, − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1曲线上箭头表示随 A A A增大, − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1的变化方向;若 Γ G \Gamma_G ΓG曲线和 − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1曲线无交点,表明上式无 ω \omega ω的正实数解;
非线性系统的稳定性判据:若 Γ G \Gamma_G ΓG曲线不包围 − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1曲线,则非线性系统稳定;若 Γ G \Gamma_G ΓG曲线包围 − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1曲线,则非线性系统不稳定;
3.4.3 非线性系统存在周期运动时的稳定性分析
当 Γ G \Gamma_G ΓG曲线和 − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1曲线存在交点时,下式成立:
G ( j ω ) = − 1 N ( A ) G({\rm j}\omega)=-\frac{1}{N(A)} G(jω)=−N(A)1
可得:
∣ G ( j ω ) ∣ = ∣ 1 N ( A ) ∣ , ∠ [ G ( j ω ) ] = − π − ∠ [ N ( A ) ] |G({\rm j}\omega)|=\left|\frac{1}{N(A)}\right|,\angle[G({\rm j}\omega)]=-\pi-\angle[N(A)] ∣G(jω)∣=
N(A)1
,∠[G(jω)]=−π−∠[N(A)]
或
R e [ G ( j ω ) N ( A ) ] = − 1 , I m [ G ( j ω ) N ( A ) ] = 0 {\rm Re}[G({\rm j}\omega)N(A)]=-1,{\rm Im}[G({\rm j}\omega)N(A)]=0 Re[G(jω)N(A)]=−1,Im[G(jω)N(A)]=0
上两式可解得交点处的频率 ω \omega ω和幅值 A A A;
系统处于周期运动时,非线性环节的输入近似为等幅振荡:
x ( t ) = A sin ω t x(t)=A\sin\omega{t} x(t)=Asinωt
即每一个交点对应着一个周期运动,如果该周期运动能够维持,即在外界小扰动作用下使系统偏离该周期运动,当该扰动消失后,系统的运动仍能恢复原周期运动,则称为稳定的周期运动;
下图给出非线性系统存在周期运动的四种形式,图中 Γ G \Gamma_G ΓG曲线和 − 1 N ( A ) -\displaystyle\frac{1}{N(A)} −N(A)1的交点为 N 0 = − 1 N ( A 0 ) N_0=-\displaystyle\frac{1}{N(A_0)} N0=−N(A0)1,负倒描述函数上的一点 N i N_i Ni对应的幅值为 A i A_i Ai;
- 图 ( a ) {\rm (a)} (a):设系统周期运动的幅值为 A 0 A_0 A0,当外界扰动使非线性环节输入振幅减小为 A 1 A_1 A1时,由于 Γ G \Gamma_G ΓG曲线包围 ( − 1 / N ( A 1 ) , j 0 ) (-1/N(A_1),{\rm j}0) (−1/N(A1),j0)点,系统不稳定,振幅将增大,最终回到 N 0 N_0 N0点;当外界扰动使输入振幅增大为 A 2 A_2 A2,由于 Γ G \Gamma_G ΓG曲线不包围 ( − 1 / N ( A 2 ) , j 0 ) (-1/N(A_2),{\rm j}0) (−1/N(A2),j0)点,系统稳定,振幅将衰减,最终也回到 N 0 N_0 N0点;即 N 0 N_0 N0点对应的周期运动是稳定的;
- 图 ( b ) {\rm (b)} (b):当外扰动使系统偏离周期运动至 N 2 N_2 N2点,即使其幅值由 A 0 A_0 A0增大为 A 2 A_2 A2时,系统不稳定,振幅将进一步增大,最终发散至无穷;当外扰动使系统偏离周期运动至 N 1 N_1 N1点,即使其幅值由 A 0 A_0 A0减小为 A 1 A_1 A1时,系统稳定,振幅将进一步减小,最终衰减为零;即 N 0 N_0 N0点对应的周期运动是不稳定的;
- 图 ( c ) {\rm (c)} (c): Γ G \Gamma_G ΓG曲线和 − 1 / N ( A ) -1/N(A) −1/N(A)曲线有两个交点 N 10 N_{10} N10和 N 20 N_{20} N20,系统中存在两个周期运动,幅值分别为: A 10 A_{10} A10和 A 20 A_{20} A20,在 N 20 N_{20} N20点,外界小扰动使系统运动偏离该周期运动后,系统运动仍然能恢复该周期运动;在 N 10 N_{10} N10点,只要有外界扰动使系统运动偏离该周期运动,则系统运动或收敛至零,或趋向于 N 20 N_{20} N20点对应的周期运动; N 10 N_{10} N10点对应的周期运动是不稳定的, N 20 N_{20} N20点对应的周期运动是稳定的;
- 图 ( d ) {\rm (d)} (d): N 10 N_{10} N10点对应的周期运动是稳定的, N 20 N_{20} N20点对应的周期运动是不稳定的,外界小扰动或使系统运动发散至无穷,或趋向于幅值 N 10 N_{10} N10点对应的周期运动;
综上:在复平面上将 Γ G \Gamma_G ΓG曲线包围的区域视为不稳定区域, Γ G \Gamma_G ΓG曲线不包围的区域视为稳定区域,则周期运动稳定性判据:在 Γ G \Gamma_G ΓG曲线和 − 1 / N ( A ) -1/N(A) −1/N(A)曲线的交点处,若 − 1 / N ( A ) -1/N(A) −1/N(A)曲线沿着振幅 A A A增加的方向由不稳定区域进入稳定区域时,该交点对应的周期运动是稳定的;反之,若 − 1 / N ( A ) -1/N(A) −1/N(A)曲线沿着振幅 A A A增加的方向在交点处由稳定区域进入不稳定区域时,该交点对应的周期运动是不稳定的;
3.5 实战
实验要求:设非线性控制系统如下图所示:
其中: G 1 ( s ) = 5 , G 2 ( s ) = 1 0.4 s + 1 , G 3 ( s ) = 1 s G_1(s)=5,G_2(s)=\displaystyle\frac{1}{0.4s+1},G_3(s)=\displaystyle\frac{1}{s} G1(s)=5,G2(s)=0.4s+11,G3(s)=s1,非线性环节 N N N为死区非线性,其表达式为:
y = { x + 2 , x < − 2 0 , ∣ x ∣ ≤ 2 x − 2 , x > 2 y=\begin{cases} &x+2,&x<-2\\ &0,&|x|≤2\\ &x-2,&x>2 \end{cases} y=⎩
⎨
⎧x+2,0,x−2,x<−2∣x∣≤2x>2
使用 S I M U L I N K {\rm SIMULINK} SIMULINK分析系统单位阶跃响应,并绘制响应曲线。
解:
【 S I M U L I N K {\rm SIMULINK} SIMULINK模型】
# SIMULINK设置:
# 所需模块:
Continuous模块组下的Transfer Fcn模块;
Discontinuties模块组下的Dead Zone模块;
Sources模块组下的Step模块;
Sources模块组下的Clocks模块;
Math Operations模块组下的Add模块;
Sinks模块组下的Scope模块;
Sinks模块组下的To Workspace模块;
# 参数设置
Dead Zone模块的死区范围:[-2,2];
To workspace模块的save format选择"Array";
【阶跃响应曲线】