Método de Cálculo de Discretización para Sistema Continuo bajo Transformación Z Estándar

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 ,metroZ ,nortem ) 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 )=ins+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+dyo11z- 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( miT ) W( ) 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 )=1z- 1 mi- re1TC1+1z- 1 mi- re2TC2=1z- 1 mi4T _ 0,75+1z- 1 mi- 8 T1.75 t = 0,01=1z- 1 mi0,04 0,75+1z- 1 mi0,081.75=10.96079 z1 0,75+10,92312 z11.75=11.88391z _1+0.88692z _210.98904 z 1=z21.88391z _+0.88692z20.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( ) ∣ W re ( ej ω T ) ∣ \big\lvert W_d \left( e^{j \omega T} \right)\big\rvert Wre( miT ) 。对于前者:
∣ 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( ) = ( s+4 ) ( s+8 )s+1 s ==Vaya2+42 Vaya2+82 Vaya2+1 ω = 1=17 sesenta y cinco 2 =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( miT ) = z21.88391z _+0.88692z20.98904 z z = mijωT _= mi2 jo T1.88391 ejωT _+0.88692mi2 jo T0.98904 mijωT _ t = 0,01 , ω = 1= mi0,02 unidades1.88391 e0.01j _+0.88692mi0,02 unidades0.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( miT ) W( ) Wre( z )=4.733170.04254Wre( z )=0.008988z21.88391z _+0.88692z20.98904 z=z21.88391z _+0.886920.008988 z20.008889 z

Supongo que te gusta

Origin blog.csdn.net/weixin_58399148/article/details/131056017
Recomendado
Clasificación