Método de Cálculo de Discretización para Sistema Continuo Estacionario Lineal bajo Transformada Z Estándar
*Este artículo solo presenta el método de cálculo, no el principio de cálculo.
1. Forma de sistema continuo
Para un sistema continuo lineal estable, su forma general es
W ( s ) = bmsm + bm − 1 sm − 1 + ⋯ + b 1 s + b 0 ansn + an − 1 sn − 1 + ⋯ + a 1 s + a 0 ( norte , metro ∈ Z , norte ≥ metro ) W(s) = \frac{b_m s^m + b_{m-1} s^{m-1} + \cdots + b_1 s + b_0}{a_n s^ n + a_{n-1} s^{n-1} + \cdots + a_1 s + a_0} \quad \left( n, m \in Z, n \geq m \right)W ( s )=ansnorte+an − 1snorte - 1+⋯+a1s+a0bmsmetro+bmetro − 1smetro - 1+⋯+b1s+b0( n ,metro∈Z ,norte≥m ) Sea su polos = − di ( i = 1 , ⋯ , n ) s = -d_i \left( i = 1, \cdots, n \right)s=- reyo( yo=1 ,⋯,n ) , entonces la fórmula anterior se puede reducir a
W ( s ) = ∑ incis + di (1) W(s) = \sum _i ^ n \frac{c_i}{s+d_i} \tag{1}W ( s )=i∑ns+dyoCyo( 1 ) Entre ellos, el coeficienteci c_iCyoSe puede calcular mediante el teorema del residuo:
ci = W ( s ) ( s + di ) ∣ s = − di c_i = W(s) \left( s+ d_i \right) \Big\rvert _{s = -d_i}Cyo=W ( s )( s+dyo)
s = - reyo
2. Sustitución de polos
sistema continuo es ssdominio s , el sistema discreto eszzdominio z , existe sustitución de variables entre los dos. La sustitución de variables tiene la siguiente forma:
1 s + di → 1 1 − z − 1 e − di T (2) \frac{1}{s+d_i} \rightarrow \frac{1}{1-z^{-1 } e^{-d_iT}} \tag{2}s+dyo1→1−z- 1 mi- reyoT1( 2 ) dondeTTT es la unidad de tiempo de discretización (como 0,01 s).
Sustituyendo cada elemento de la fórmula (1) en la forma de la fórmula (2), la función de transferencia del sistema continuoW ( s ) W(s)W ( s ) se transforma en su forma de sustitución de polosW d ( z ) W_d (z)Wre( z )。
3. Sistema discreto equivalente
Después de obtener la forma de sustitución de polos W d ( s ) W_d (s)Wre( s ) , el sistema continuo original W ( s ) W(s)se obtiene mediante la siguiente fórmulaEl sistema discreto equivalenteW ( z ) W(z) de W ( s )W ( z ):
W ( z ) = ∣ W ( j ω ) ∣ ∣ W re ( ej ω T ) ∣ W re ( z ) (3) W(z) = \frac{ \big\lvert W \left( j \omega \right) \big\rvert}{ \big\lvert W_d \left( e^{j \omega T} \right)\big\rvert} W_d (z) \tag{3}W ( z )=
Wre( mijω T )
W( jω )
Wre( z )( 3 )其中∣ W ( ⋅ ) ∣ \big\lvert W(\cdot) \big\rvert
W ( ⋅ )
指W ( ⋅ ) W(\cdot)Módulo de W ( ⋅ ) ; TTT es la unidad de tiempo discretizada, y se toma el valor de acuerdo al problema real; ω \omegaω 1.
4. Ejemplo de cálculo
Tomando como ejemplo la siguiente función de transferencia, calcule la función de transferencia de su sistema discreto equivalente. Intervalo de tiempo de discretización T = 0,01 s T = 0,01 sT=0.01 s。
W ( s ) = s + 1 s 2 + 12 s + 32 = s + 1 ( s + 4 ) ( s + 8 ) W(s) = \frac{s+1}{s^2 + 12s + 32} = \frac{s+1}{(s+4)(s+8)}W ( s )=s2+12 segundos+32s+1=( s+4 ) ( s+8 )s+1(1) Primero use el teorema del residuo para escribirlo en forma de fórmula (1).
c 1 = W ( s ) ( s + 4 ) ∣ s = − 4 = s + 1 s + 8 ∣ s = − 4 = − 0,75 c 2 = W ( s ) ( s + 8 ) ∣ s = − 8 = s + 1 s + 4 ∣ s = − 8 = 1,75 c_1 = W(s) (s+4) \Big\rvert_{s=-4} = \frac{s+1}{s+8} \Bigg\ rvert_{s=-4} = -0.75 \\ c_2 = W(s) (s+8) \Big\rvert_{s=-8} = \frac{s+1}{s+4} \Bigg\rvert_ {s=-8} = 1,75C1=W ( s ) ( s+4 )
s = − 4=s+8s+1
s = − 4=− 0,75C2=W ( s ) ( s+8 )
s = − 8=s+4s+1
s = − 8=1.75于是W ( s ) = c 1 s + d 1 + c 2 s + d 2 = − 0.75 s + 4 + 1.75 s + 8 W(s) = \frac{c_1}{s+d_1} + \frac{ c_2}{s+d_2} = \frac{-0,75}{s+4} + \frac{1,75}{s+8}W ( s )=s+d1C1+s+d2C2=s+4− 0,75+s+81.75
(2) Sustitución de polos. De acuerdo con la fórmula (2), la función de transferencia de la forma de sustitución de polos se puede obtener:
W d ( z ) = c 1 1 − z − 1 e − d 1 T + c 2 1 − z − 1 e − d 2 T = − 0,75 1 − z − 1 mi − 4 T + 1,75 1 − z − 1 mi − 8 T ∣ T = 0,01 = − 0,75 1 − z − 1 mi − 0,04 + 1,75 1 − z − 1 mi − 0,08 = − 0,75 1 − 0,96079 z − 1 + 1,75 1 − 0,92312 z − 1 = 1 − 0,98904 z − 1 1 − 1,88391 z − 1 + 0,88692 z − 2 = z 2 − 0,98904 zz 2 − 1,88391 z + 0 .88692 \begin{alineado} W_d (z ) &= \frac{c_1}{1-z^{-1} e^{-d_1T}} + \frac{c_2}{1-z^{-1} e^{-d_2T}} \ \ &= \frac{-0.75}{1-z^{-1} e^{-4T}} + \frac{1.75}{1-z^{-1} e^{-8T}} \Bigg\ rvert_{T =0.01} \\ &= \frac{-0.75}{1-z^{-1} e^{-0.04}} + \frac{1.75}{1-z^{-1} e^{ -0.08} } \\ &= \frac{-0.75}{1-0.96079z^{-1}} + \frac{1.75}{1-0.92312z^{-1}} \\ &= \frac{1 -0.98904z ^{-1}}{1-1.88391z^{-1}+0.88692z^{-2}} \\ &= \frac{z^2-0.98904z}{z^2-1.88391z+ 0.88692} \ final {alineado}Wre( z )=1−z- 1 mi- re1TC1+1−z- 1 mi- re2TC2=1−z- 1 mi− 4T _− 0,75+1−z- 1 mi- 8 T1.75
t = 0,01=1−z- 1 mi− 0,04− 0,75+1−z- 1 mi− 0,081.75=1−0.96079 z− 1− 0,75+1−0,92312 z− 11.75=1−1.88391z _− 1+0.88692z _− 21−0.98904 z− 1=z2−1.88391z _+0.88692z2−0.98904 z
(3) Sistema discreto equivalente. Según la fórmula (3), es necesario calcular dos valores de módulo ∣ W ( j ω ) ∣ \big\lvert W \left( j \omega \right) \big\rvert
W( jω )
和∣ W re ( ej ω T ) ∣ \big\lvert W_d \left( e^{j \omega T} \right)\big\rvert
Wre( mijω T )
。对于前者:
∣ W ( j ω ) ∣ = ∣ s + 1 ( s + 4 ) ( s + 8 ) ∣ s = j ω = ω 2 + 1 ω 2 + 4 2 ω 2 + 8 2 ∣ ω = 1 = 2 17 65 = 0.04254 \begin{alineado} \Big\lvert W \left( j \omega \right) \Big\rvert &= \Bigg\lvert \frac{s+1}{(s+4)(s +8)} \Bigg\rvert_{s=j \omega} \\ &= \frac{ \sqrt{\omega^2 + 1} }{ \sqrt{\omega^2 + 4^2} \sqrt{\ omega^2 + 8^2} } \Bigg\rvert_{\omega = 1} = \frac{ \sqrt{2} }{ \sqrt{17} \sqrt{65} } = 0,04254 \end{alineado}
W( jω )
=
( s+4 ) ( s+8 )s+1
s = jω=Vaya2+42Vaya2+82Vaya2+1
ω = 1=17sesenta y cinco2=0.04254Para despues:
∣ w d (eJ Ω t) ∣ = ∣ z 2 - 0.98904 zz 2 - 1.88391 z + 0.88692 ∣ z = eJ ω t = ∣ e 2 j ω t - 0.98904 eJ Ω t e 2 j ω t - 1.88391 eJ Ω + 0.88692 ∣ T = 0.01 , ω = 1 = ∣ e 0.02 j − 0.98904 e 0.01 je 0.02 j − 1.88391 e 0.01 j + 0.88692 ∣ = ∣ ( cos 0.02 + j sin 0.02 ) − 0.98904 ( porque 0.01 + j sin 0.01 ) ( cos 0.02 + j sin 0.02 ) − 1.88391 ( cos 0.01 + j sin 0.01 ) + 0.88692 ∣ = ∣ 0.01081 + 0.01011 j 0.002904 + 0. 0011599 j ∣ = ∣ 4,40952 + 1,72018 j ∣ = 4,73317 \ begin{alineado} \Big\lvert W_d \left( e^{j \omega T} \right) \Big\rvert &= \Bigg\lvert \frac{z^2-0.98904z}{z^2-1.88391z +0.88692} \Bigg\rvert_{z = e^{j \omega T} } \\ &= \Bigg\lvert \frac{e^{2 j \omega T}-0.98904e^{j \omega T}} {e^{2 j \omega T}-1.88391e^{j \omega T}+0.88692} \Bigg\rvert_{T=0.01, \\\ \omega = 1} \\ &= \Bigg\lvert \frac {e^{0.02j}-0.98904e^{0.01j}}{e^{0.02j}-1.88391e^{0.01j}+0.88692} \Bigg\rvert \\ &= \Bigg\lvert \frac{ \left( \cos 0.02+ j \sin 0.02 \right) -0.98904 \left( \cos 0.01+ j \sin 0.01 \right)}{\left( \cos 0.02+ j \sin 0.02 \right)-1.88391\left( \cos 0.01+ j \sin 0.01 \right)+0.88692 } \Bigg\rvert \\ &= \Bigg\lvert \frac{0.01081+0.01011j}{0.002904+0.0011599j} \Bigg\rvert \\ &= \Big\lvert 4.40952+1.72018j \Big\rvert = 4.73317 \ fin {alineado}
Wre( mijω T )
=
z2−1.88391z _+0.88692z2−0.98904 z
z = mijωT _=
mi2 jo T−1.88391 ejωT _+0.88692mi2 jo T−0.98904 mijωT _
t = 0,01 , ω = 1=
mi0,02 unidades−1.88391 e0.01j _+0.88692mi0,02 unidades−0.98904 mi0,01 j
=
( porque0.02+jpecado0.02 )−1.88391( porque0.01+jpecado0.01 )+0.88692( porque0.02+jpecado0.02 )−0.98904( porque0.01+jpecado0.01 )
=
0.002904+0.0011599j _0.01081+0,01011 días
=
4.40952+1.72018j _
=4.73317代入(3),即可得到离散系统的传函:
W ( z ) = ∣ W ( j ω ) ∣ ∣ W d ( ej ω T ) ∣ W d ( z ) = 0.04254 4.73317 W d ( z ) = 0.008988 z 2 − 0.98904 zz 2 − 1.88391 z + 0.88692 = 0.008988 z 2 − 0.008889 zz 2 − 1.88391 z + 0.88692 \begin{aligned} W(z) &= \frac{ \big\lvert W \left( j \omega \ right) \big\rvert}{ \big\lvert W_d \left( e^{j \omega T} \right)\big\rvert} W_d (z) = \frac{0.04254}{4.73317} W_d (z) \ \ &= 0.008988 \frac{z^2-0.98904z}{z^2-1.88391z+0.88692} \\ &= \frac{0.008988z^2-0.008889z}{z^2-1.88391z+0.88692} \ fin {alineado}W ( z )=
Wre( mijω T )
W( jω )
Wre( z )=4.733170.04254Wre( z )=0.008988z2−1.88391z _+0.88692z2−0.98904 z=z2−1.88391z _+0.886920.008988 z2−0.008889 z