Varios tipos estándar de sistemas de control automático (con ejemplos)

1. Allanando el camino para el conocimiento relacionado con el espacio de estados

Antes de leer este artículo, el lector predeterminado ya conoce los conocimientos relevantes del espacio de estados, que incluyen:

  • Existe más de una forma de expresión del espacio de estados, la cual se puede expresar en x ~ = T x \tilde x = TxX~=Transformar sobre la base de T x ;
  • El espacio de estados es un medio para estudiar las variables internas del proceso dinámico del sistema y las variables de estado seleccionadas xi x_iXyoen más de una forma;
  • Todas las variables de estado xi x_i en el sistemaXyoNo tiene por qué tener un significado físico, a veces sólo se establece por conveniencia de la representación matemática; además, no todas las variables de estado serán observadas o medidas.

Todas las expresiones del espacio de estados en este artículo se establecen con base en la siguiente ecuación diferencial:
y ( n ) + an − 1 y ( n − 1 ) + an − 2 y ( n − 2 ) + ⋯ + a 2 y ¨ + a 1 y ˙ + a 0 y = bnu ( n ) + bn − 1 u ( n − 1 ) + ⋯ + b 1 u ˙ + b 0 u (1) y^{(n)} + a_{n-1} y ^{(n-1)} + a_{n-2} y^{(n-2)} + \cdots + a_2 \ddot y + a_1 \dot y + a_0 y = \\ b_n u^{(n ) } + b_{n-1} u^{(n-1)} + \cdots + b_1 \dot u + b_0 u \tag{1}y( n )+anorte 1y( norte - 1 )+anorte 2y( norte - 2 )++a2y¨+a1y˙+a0y=bnortetu( n )+bnorte 1tu( norte - 1 )++b1tu˙+b0tu( 1 ) Tenga en cuenta queen la fórmula anterior,y ( n ) y^{(n)}yCoeficiente an = 1antes ( n ) a_n = 1anorte=1 .

En cambio, el objetivo es escribirlo en la siguiente forma de espacio de estados:
X ˙ = AX + BUY = CX + DU \dot {\bm X} = {\bm A} {\bm X} + {\bm B} {\bm U} \\ {\bm Y} = {\bm C} {\bm X} + {\bm D} {\bm U}X˙=una x+BU _Y=CX _+DU_ _

2. La entrada del sistema [no contiene] la partida derivada

Es decir, el lado derecho de la fórmula (1) solo contiene uuu,式(1)变为:
y ( n ) + a n − 1 y ( n − 1 ) + a n − 2 y ( n − 2 ) + ⋯ + a 2 y ¨ + a 1 y ˙ + a 0 y = b 0 u (2) y^{(n)} + a_{n-1} y^{(n-1)} + a_{n-2} y^{(n-2)} + \cdots + a_2 \ddot y + a_1 \dot y + a_0 y = b_0 u \tag{2} y( n )+anorte 1y( norte - 1 )+anorte 2y( norte - 2 )++a2y¨+a1y˙+a0y=b0tu( 2 ) En este momento, la cantidad de estado se selecciona como:
{ x ˙ 1 = x 2 x ˙ 2 = x 3 ⋮ x ˙ n − 1 = xnx ˙ n = − a 0 x 1 − a 1 x 2 − ⋯ − an − 1 xn + b 0 uy = x 1 \begin{cases} \dot x_1 = x_2 \\ \dot x_2 = x_3 \\ \vdots \\ \dot x_{n-1} = x_n \\ \dot x_n = -a_0 x_1 - a_1 x_2 - \cdots - a_{n-1} x_n + b_0 u \\ y = x_1 \end{cases} X˙1=X2X˙2=X3X˙norte 1=XnorteX˙norte=un0X1a1X2anorte 1Xnorte+b0tuy=X1Dado el infinitesimal:
A = [ 0 1 0 ⋯ 0 0 0 1 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 ⋯ 1 − a 0 − a 1 − a 2 ⋯ − an − 1 ] {\bm A} = \left [ \begin{matrix} 0 & 1 & 0 & \cdots & 0 \\ 0 & 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 0 & 0&0&\ cdots&1\\-a_0&-a_1&-a_2&\cdots&-a_{n-1}\end{matrix}\right]A= 000un0100un1010un2001unnorte 1 B = [ 0 0 ⋮ 0 b 0 ] , C = [ 1 0 ⋯ 0 ] , D = 0 (3) {\bm B} = \left[ \begin{matrix} 0 \\ 0 \\ \vdots \\ 0 \\ b_0 \end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} 1 & 0 & \cdots & 0 \end{matrix} \right], \quad {\ bm D} = 0 \etiqueta{3}B= 000b0 ,C=[100],D=0( 3 )

3. [Contiene] el elemento derivado en la entrada del sistema.

En este caso, la ecuación diferencial del sistema tiene la forma de fórmula (1):
y ( n ) + an − 1 y ( n − 1 ) + an − 2 y ( n − 2 ) + ⋯ + a 2 y ¨ + a 1 y ˙ + a 0 y = bnu ( n ) + bn − 1 u ( n − 1 ) + ⋯ + b 1 u ˙ + b 0 u (1) y^{(n)} + a_{n- 1} y^ {(n-1)} + a_{n-2} y^{(n-2)} + \cdots + a_2 \ddot y + a_1 \dot y + a_0 y = \\ b_n u^{ (n)} + b_{n-1} u^{(n-1)} + \cdots + b_1 \dot u + b_0 u \tag{1}y( n )+anorte 1y( norte - 1 )+anorte 2y( norte - 2 )++a2y¨+a1y˙+a0y=bnortetu( n )+bnorte 1tu( norte - 1 )++b1tu˙+b0tu( 1 )此时状态量选择为:
{ x ˙ 1 = x 2 + h 1 ux ˙ 2 = x 3 + h 2 u ⋮ x ˙ n − 1 = xn + hn − 1 ux ˙ n = − a 0 x 1 − a 1 x 2 − ⋯ − an − 1 xn + hnuy = x 1 + h 0 u \begin{cases} \dot x_1 = x_2 + h_1 u \\ \dot x_2 = x_3 + h_2 u \\ \vdots \ \ \dot x_{n-1} = x_n + h_{n-1} u \\ \dot x_n = -a_0 x_1 - a_1 x_2 - \cdots - a_{n-1} x_n + h_n u \\ y = x_1 + h_0 u \end{casos} X˙1=X2+h1tuX˙2=X3+h2tuX˙norte 1=Xnorte+hnorte 1tuX˙norte=un0X1a1X2anorte 1Xnorte+hnortetuy=X1+h0tudonde el nuevo parámetro hola h_ihyo为:
h 0 = bnhi = bn − i − ∑ j = 1 ian − jhi − j (4) h_0 = b_n \\ h_i = b_{ni} - \sum _{j=1} ^i a_{nj} h_ {ij} \etiqueta{4}h0=bnortehyo=bnorte yoj = 1yoanorte - jhyo j( 4 ) Especifica la ecuación:
A = [ 0 1 0 ⋯ 0 0 0 1 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 ⋯ 1 − a 0 − a 1 − a 2 ⋯ − an − 1 ] {\bm A} = \left[ \begin{matrix} 0 & 1 & 0 & \cdots & 0 \\ 0 & 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots\ \0&0&0& \cdots&1\\-a_0&-a_1&-a_2&\cdots&-a_{n-1}\end{matrix}\right]A= 000un0100un1010un2001unnorte 1 B = [ h 1 h 2 ⋮ hn − 1 hn ] , C = [ 1 0 ⋯ 0 ] , D = h 0 (5) {\bm B} = \left[ \begin{matrix} h_1 \\ h_2 \\ \vdots \\ h_{n-1} \\ h_n \end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} 1 y 0 & \cdots & 0 \end{matrix } \right], \quad {\bm D} = h_0 \tag{5}B= h1h2hnorte 1hnorte ,C=[100],D=h0( 5 ) Matriz visibleA {\bm A}A es el mismo que en la fórmula (3).

4. La entrada del sistema [contiene] término derivado, pero bn = 0 {\bm b}_{\bm n} {\bm =} {\bm 0}bnorte= 0

En este momento, el lado derecho de la ecuación diferencial (1) del sistema comienza desde bn − 1 u ( n − 1 ) b_{n-1} u^{(n-1)}bnorte 1tu( n 1 )开始:
y ( n ) + an − 1 y ( n − 1 ) + an − 2 y ( n − 2 ) + ⋯ + a 2 y ¨ + a 1 y ˙ + a 0 y = bn − 1 u ( n − 1 ) + ⋯ + b 1 u ˙ + b 0 u (1) y^{(n)} + a_{n-1} y^{(n-1)} + a_{n-2 } y^{(n-2)} + \cdots + a_2 \ddot y + a_1 \dot y + a_0 y = \\ b_{n-1} u^{(n-1)} + \cdots + b_1 \ punto u + b_0 u \tag{1}y( n )+anorte 1y( norte - 1 )+anorte 2y( norte - 2 )++a2y¨+a1y˙+a0y=bnorte 1tu( norte - 1 )++b1tu˙+b0tu( 1 ) Sustituir la fórmula (4) es fácil de saber,h 0 = 0 h_0 = 0h0=0 , entonces se puede seleccionar un nuevo conjunto de variables de estado de acuerdo con las siguientes reglas:
{ x ˙ 1 = − a 0 xn + b 0 ux ˙ 2 = x 1 − a 1 xn + b 1 u ⋮ x ˙ n − 1 = xn − 2 − an − 2 xn + bn − 2 ux ˙ n = xn − 1 − an − 1 xn + bn − 1 uy = xn \begin{cases} \dot x_1 = -a_0 x_n + b_0 u \\ \ punto x_2 = x_1 - a_1 x_n + b_1 u \\ \vdots \\ \dot x_{n-1} = x_{n-2} - a_{n-2} x_n + b_{n-2} u \\ \ punto x_n = x_{n-1} - a_{n-1} x_n + b_{n-1} u \\ y = x_n \end{cases} X˙1=un0Xnorte+b0tuX˙2=X1a1Xnorte+b1tuX˙norte 1=Xnorte 2anorte 2Xnorte+bnorte 2tuX˙norte=Xnorte 1anorte 1Xnorte+bnorte 1tuy=XnorteDa la función infinitesimal:
A = [ 0 0 ⋯ 0 − a 0 1 0 ⋯ 0 − a 1 0 1 ⋯ 0 − a 2 ⋮ ⋮ ⋱ ⋮ ⋮ 0 0 ⋯ 1 − an − 1 ] {\bm A} = \ izquierda[\begin{matrix} 0&0&\cdots&0&-a_0\\1&0&\cdots&0&-a_1\\0&1&\cdots&0&-a_2\\\\vdots& \vdots&\ddots&\vdots&\vdots\\0&0&\cdots&1&-a_{n-1 }\end{matriz}\right]A= 010000100001un0un1un2unnorte 1 B = [ b 0 b 1 ⋮ bn − 1 ] , C = [ 0 0 ⋯ 1 ] , D = 0 (6) {\bm B} = \left[ \begin{matrix} b_0 \\ b_1 \\ \vdots \\ b_{n-1} \end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} 0 & 0 & \cdots & 1 \end{matrix} \right], \quad {\bm D} = 0 \tag{6}B= b0b1bnorte 1 ,C=[001],D=0( 6 )

5. Tipo estándar controlable

En este caso, la ecuación diferencial del sistema tiene la forma de fórmula (1):
y ( n ) + an − 1 y ( n − 1 ) + an − 2 y ( n − 2 ) + ⋯ + a 2 y ¨ + a 1 y ˙ + a 0 y = bnu ( n ) + bn − 1 u ( n − 1 ) + ⋯ + b 1 u ˙ + b 0 u (1) y^{(n)} + a_{n- 1} y^ {(n-1)} + a_{n-2} y^{(n-2)} + \cdots + a_2 \ddot y + a_1 \dot y + a_0 y = \\ b_n u^{ (n)} + b_{n-1} u^{(n-1)} + \cdots + b_1 \dot u + b_0 u \tag{1}y( n )+anorte 1y( norte - 1 )+anorte 2y( norte - 2 )++a2y¨+a1y˙+a0y=bnortetu( n )+bnorte 1tu( norte - 1 )++b1tu˙+b0tu( 1 )此时状态量选择为:
{ x ˙ 1 = x 2 x ˙ 2 = x 3 ⋮ x ˙ n = − a 0 x 1 − a 1 x 2 − ⋯ − an − 1 xn + uy = − β 0 x 1 − β 1 x 2 − ⋯ − β n − 1 xn \begin{cases} \dot x_1 = x_2 \\ \dot x_2 = x_3 \\ \vdots \\ \dot x_n = -a_0 x_1 - a_1 x_2 - \cdots - a_{n-1} x_n + u \\ y = -\beta_0 x_1 - \beta_1 x_2 - \cdots - \beta_{n-1} x_n \end{cases} X˙1=X2X˙2=X3X˙norte=un0X1a1X2anorte 1Xnorte+tuy=segundo0X1b1X2bnorte 1Xnortedonde el nuevo parámetro β i \beta_ibyo为:
{ β 0 = b 0 − a 0 bn β 1 = b 1 − a 1 bn ⋮ β i = bi − aibn ⋮ β n − 2 = bn − 2 − an − 2 bn β n − 1 = bn − 1 − an − 1 bn (7) \begin{casos} \beta_0 = b_0 - a_0 b_n \\ \beta_1 = b_1 - a_1 b_n \\ \vdots \\ \beta_i = b_i - a_i b_n \\ \vdots \\ \beta_ {n-2} = b_{n-2} - a_{n-2} b_n \\ \beta_{n-1} = b_{n-1} - a_{n-1} b_n \end{cases} \ etiqueta {7} b0=b0a0bnorteb1=b1a1bnortebyo=byoayobnortebnorte 2=bnorte 2anorte 2bnortebnorte 1=bnorte 1anorte 1bnorte( 7 ) Observe el infinitesimal:
A = [ 0 1 0 ⋯ 0 0 0 1 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 ⋯ 1 − a 0 − a 1 − a 2 ⋯ − an − 1 ] {\bm A} = \left[ \begin{matrix} 0 & 1 & 0 & \cdots & 0 \\ 0 & 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots\ \0&0&0& \cdots&1\\-a_0&-a_1&-a_2&\cdots&-a_{n-1}\end{matrix}\right]A= 000un0100un1010un2001unnorte 1 B = [ 0 0 ⋮ 0 1 ] , C = [ β 0 β 1 ⋯ β n − 1 ] , D = 0 (8) {\bm B} = \left[ \begin{matrix} 0 \\ 0 \\ \vdots \\ 0 \\ 1 \end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} \beta_0 & \beta_1 & \cdots & \beta_{n-1} \ fin{matriz} \right], \quad {\bm D} = 0 \tag{8}B= 0001 ,C=[b0b1bnorte 1],D=0( 8 )

6. Estándar considerable

La forma canónica observable es la transpuesta de la forma canónica controlable:
A c = A o T , B c = C o T , C c = B o T {\bm A_c} = {\bm A_o}^T, {\bm B_c} = {\bm C_o}^T, {\bm C_c} = {\bm B_o}^TAc=AoT ,Bc=CoT ,Cc=BoT donde los subíndicesc, oc, oc ,o denotan matrices controlables y observables, respectivamente. Esto da:
A = [ 0 0 ⋯ 0 − a 0 1 0 ⋯ 0 − a 1 0 1 ⋯ 0 − a 2 ⋮ ⋮ ⋱ ⋮ ⋮ 0 0 ⋯ 1 − an − 1 ] {\bm A} = \left [ \begin{matrix} 0 & 0 & \cdots & 0 & -a_0 \\ 1 & 0 & \cdots & 0 & -a_1 \\ 0 & 1 & \cdots & 0 & -a_2 \\ \vdots & \vdots & \ddots & \vdots & \vdots \\ 0 & 0 & \cdots & 1 & -a_{n-1} \end{matrix} \right]A= 010000100001un0un1un2unnorte 1 B = [ β 0 β 1 β 2 ⋮ β n − 1 ] , C = [ 0 0 ⋯ 1 ] , D = 0 (9) {\bm B} = \left[ \begin{matrix} \beta_0 \\ \ beta_1 \\ \beta_2 \\ \vdots \\ \beta_{n-1} \end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} 0 & 0 & \cdots & 1 \end{matrix} \right], \quad {\bm D} = 0 \tag{9}B= b0b1b2bnorte 1 ,C=[001],D=0( 9 )

7. El tipo equivalente de la ecuación característica de la función de transferencia [cuando no hay raíz repetida]

Nota: ¡Este tipo sólo se puede utilizar cuando la ecuación característica de la función de transferencia no tiene raíces repetidas!

Según los diferentes métodos de selección de las variables de estado, generalmente existen dos expresiones del tipo equivalente.

A partir de la fórmula (1), la relación de la función de transferencia en el plano complejo se puede escribir:
W ( s ) = bnsn + bn − 1 sn − 1 + ⋯ + b 1 s + b 0 sn + an − 1 sn − 1 + ⋯ + a 1 s + a 0 = N ( s ) D ( s ) W(s) = \frac{ b_n s^n + b_{n-1} s^{n-1} + \cdots + b_1 s + b_0 } { s^n + a_{n-1} s^{n-1} + \cdots + a_1 s + a_0 } = \frac{ N(s) }{ D(s) }W ( s )=snorte+anorte 1snorte - 1++a1s+a0bnortesnorte+bnorte 1snorte - 1++b1s+b0=D ( s )norte ( s )donde D ( s ) = 0 D(s)=0D ( s )=0 es la ecuación característica de la función de transferencia y su raíz característica esλ i \lambda_iyoyo. Entonces el denominador se puede escribir como:
D ( s ) = ( s − λ 1 ) ( s − λ 2 ) ⋯ ( s − λ n ) D(s) = \left(s - \lambda_1 \right) \left( s - \lambda_2 \right) \cdots \left(s - \lambda_n \right)D ( s )=( syo1)( syo2)( syonorte)此时传递函数即:
W ( s ) = N ( s ) D ( s ) = bnsn + bn − 1 sn − 1 + ⋯ + b 1 s + b 0 ( s − λ 1 ) ( s − λ 2 ) ⋯ ( s − λ n ) W(s) = \frac{ N(s) }{ D(s) } = \frac{ b_n s^n + b_{n-1} s^{n-1} + \ cdots + b_1 s + b_0 }{ \left(s - \lambda_1 \right) \left(s - \lambda_2 \right) \cdots \left(s - \lambda_n \right) }W ( s )=D ( s )norte ( s )=( syo1)( syo2)( syonorte)bnortesnorte+bnorte 1snorte - 1++b1s+b0展开即得:
W ( s ) = N ( s ) D ( s ) = Y ( s ) U ( s ) = ∑ i = 1 ncis − λ i (10) W(s) = \frac{ N(s) }{ D(s) } = \frac{ Y(s) }{ U(s) } = \sum _{i=1} ^n \frac{ c_i }{ s - \lambda_i } \tag{10}W ( s )=D ( s )norte ( s )=U ( s )Y(s)=yo = 1nortesyoyoCyo( 10 ) cic_iCyoes la función de transferencia W ( s ) W(s)W ( s ) en el poloλ i \lambda_iyoyo处的留数
ci = ( s − λ i ) W ( s ) ∣ s = λ i (11) c_i = \left(s - \lambda_i \right) W(s) \Big\rvert _{s = \ lambda_i} \etiqueta{11}Cyo=( syoyo)W ( s ) s = λyo( 11 )从而
Y ( s ) = ∑ i = 1 ncis − λ i U ( s ) Y(s) = \sum _{i=1} ^n \frac{ c_i }{ s - \lambda_i } U(s )Y(s)=yo = 1nortesyoyoCyoU ( s )

7.1 Tipo I

El cambio de diseño es
X i ( s ) = 1 s − λ i U ( s ) X_i(s) = \frac{1}{s - \lambda_i} U(s)Xyo( s )=syoyo1U ( s )则有
{ x ˙ 1 = λ 1 x 1 + ux ˙ 2 = λ 2 x 2 + u ⋮ x ˙ n = λ nxn + uy = c 1 x 1 + c 2 x 2 + ⋯ + cnxn \ comenzar{casos} \dot x_1 = \lambda_1 x_1 + u \\ \dot x_2 = \lambda_2 x_2 + u \\ \vdots \\ \dot x_n = \lambda_n x_n + u \\ y = c_1 x_1 + c_2 x_2 + \ cdots + c_n x_n \end{casos} X˙1=yo1X1+tuX˙2=yo2X2+tuX˙norte=yonorteXnorte+tuy=C1X1+C2X2++CnorteXnorteEspecifique el valor límite:
A = [ λ 1 λ 2 ⋱ λ n − 1 λ n ] {\bm A} = \left[ \begin{matrix} \lambda_1 & & & & \\ & \lambda_2 & & & \ \ &&\ddots&&\\&&&\lambda_{n-1}&\\&&&&\lambda_n\end{matrix}\right]A= yo1yo2yonorte 1yonorte B = [ 1 1 ⋮ 1 1 ] , C = [ c 1 c 2 ⋯ cn ] , D = 0 (12) {\bm B} = \left[ \begin{matrix} 1 \\ 1 \\ \vdots \ \ 1 \\ 1 \end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} c_1 & c_2 & \cdots & c_n \end{matrix} \right], \quad { \bm D} = 0 \etiqueta{12}B= 1111 ,C=[C1C2Cnorte],D=0( 12 ) dondeA {\bm A}A esuna matriz diagonal y los elementos de la diagonal principal son todos raíces característicasλ i \lambda_iyoyo, mientras que el resto de elementos son todos 0 .

7.2 Tipo II

El cambio de diseño es
X i ( s ) = cis − λ i U ( s ) X_i(s) = \frac{c_i}{s - \lambda_i} U(s)Xyo( s )=syoyoCyoU ( s ) en comparación con el primer tipo,X i ( s ) X_i(s)XyoEl numerador de ( s ) ya no es 1 sinoci c_iCyo。则有
{ x ˙ 1 = λ 1 x 1 + c 1 ux ˙ 2 = λ 2 x 2 + c 2 u ⋮ x ˙ n = λ nxn + cnuy = x 1 + x 2 + ⋯ + xn \begin{casos } \dot x_1 = \lambda_1 x_1 + c_1 u \\ \dot x_2 = \lambda_2 x_2 + c_2 u \\ \vdots \\ \dot x_n = \lambda_n x_n + c_n u \\ y = x_1 + x_2 + \cdots + x_n \end{casos} X˙1=yo1X1+C1tuX˙2=yo2X2+C2tuX˙norte=yonorteXnorte+Cnortetuy=X1+X2++XnorteEspecifique el valor límite:
A = [ λ 1 λ 2 ⋱ λ n − 1 λ n ] {\bm A} = \left[ \begin{matrix} \lambda_1 & & & & \\ & \lambda_2 & & & \ \ &&\ddots&&\\&&&\lambda_{n-1}&\\&&&&\lambda_n\end{matrix}\right]A= yo1yo2yonorte 1yonorte B = [ c 1 c 2 ⋮ cn − 1 cn ] , C = [ 1 1 ⋯ 1 ] , D = 0 (13) {\bm B} = \left[ \begin{matrix} c_1 \\ c_2 \\ \ vdots \\ c_{n-1} \\ c_n \end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} 1 & 1 & \cdots & 1 \end{matrix} \right], \quad {\bm D} = 0 \tag{13}B= C1C2Cnorte 1Cnorte ,C=[111],D=0( 13 ) dondeA {\bm A}A es una matriz diagonal, igual que la del primer tipo.

8. El tipo equivalente de ecuación característica de la función de transferencia [cuando hay múltiples raíces]

Aquí tomamos como ejemplo una ecuación característica con una raíz triple. Sea la ecuación característica en este momento
D ( s ) = ( s − λ 1 ) 3 ( s − λ 4 ) ⋯ ( s − λ n ) D(s) = \left(s - \lambda_1 \right) ^3 \ izquierda(s - \lambda_4 \right) \cdots \left(s - \lambda_n \right)D ( s )=( syo1)3( syo4)( syonorte) dondeλ 1 \lambda_1yo1es un polo real triple, y el restante λ i ( i ≠ 1 ) \lambda_i ( i \neq 1)yoyo( yo=1 ) es un polo real simple. Entonces la función de transferencia se puede descomponer en:
W ( s ) = N ( s ) D ( s ) = Y ( s ) U ( s ) = c 11 ( s − λ 1 ) 3 + c 12 ( s − λ 1 ) 2 + c 13 ( s − λ 1 ) + ∑ i = 4 ncis − λ i (14) W(s) = \frac{ N(s) }{ D(s) } = \frac{ Y(s) } { U(s) } = \frac{ c_{11} }{ \left(s - \lambda_1 \right)^3 } + \frac{ c_{12} }{ \left(s - \lambda_1 \right)^ 2 } + \frac{ c_{13} }{ \left(s - \lambda_1 \right) } + \sum_{i=4} ^n \frac{ c_i} {s - \lambda_i } \tag{14}W ( s )=D ( s )norte ( s )=U ( s )Y(s)=( syo1)3C11+( syo1)2C12+( syo1)C13+yo = 4nortesyoyoCyo( 14 ) El método de selección de la variable de estado es el mismo que el del primer tipo y el del segundo tipo. La matriz del espacio de estados en este momento es:
A = [ λ 1 1 λ 1 1 0 λ 1 λ 4 0 ⋱ λ n ] , B = [ 0 0 1 1 ⋮ 1 ] {\bm A} = \left[ \begin { matriz}{ ccc | ccc } \lambda_1 & 1 & {} & {} & {} & {} \\ {} & \lambda_1 & 1 & {} & {\bm 0} & {} \\ {} & { } & \lambda_1 & {} & {} & {} \\ \hline {} & {} & {} & \lambda_4 & {} & {} \\ {} & {\bm 0} & {} & { } & \ddots & {} \\ {} & {} & {} & {} & {} & \lambda_n \\ \end{array} \right], \quad {\bm B} = \left[ \begin { matriz}{ c } 0 \\ 0 \\ 1 \\ \hline 1 \\ \vdots \\ 1 \end{array} \right]A= yo11yo101yo1yo40yonorte ,B= 00111 C = [ c 11 c 12 c 13 c 4 ⋯ cn ] , D = 0 (15) {\bm C} = \left[ \begin{array}{ ccc | ccc } c_{11} & c_{12} & c_{13} & c_4 & \cdots & c_n \end{array} \right], \quad {\bm D} = 0 \tag{15}C=[C11C12C13C4Cnorte],D=0( 15 ) Especifique las ecuaciones:
A = [ λ 1 1 λ 1 0 1 λ 1 λ 4 0 ⋱ λ n ] , B = [ c 11 c 12 c 13 c 4 ⋮ cn ] {\bm A} = \ . left [ \begin{array}{ccc | ccc}\lambda_1&{}&{}&{}&{}&{}\\1&\lambda_1&{}&{}&{\bm0}&{}\\{}&1&\lambda_1&{} {}&{} &{}\\\hline{}&{}&{}&\lambda_4&{}&{}\\{}&{\bm0}&{}&{}&\ddots&{ } \\{}&{} &{}&{}&{}&\lambda_n\\end{array}\right],\quad{\bmB}=\left[\begin{array}{c} c_{11}\\c_{12} \\c_{13}\\\hline c_4\\\vdots\\c_n\end{array}\right]A= yo11yo110yo1yo40yonorte ,B= C11C12C13C4Cnorte C = [ 0 0 1 1 ⋯ 1 ] , D = 0 (16) {\bm C} = \left[ \begin{array}{ ccc | ccc } 0 & 0 & 1 & 1 & \cdots & 1 \end{array} \right], \quad {\bm D} = 0 \tag{16}C=[00111],D=0( 16 ) Conviene explicar que, en las fórmulas (15) y (16),A {\bm A}A ya no es una matriz diagonal completa, debido a la raíz repetidaλ 1 \lambda_1yo1La existencia de , λ 1 \lambda_1yo1A{\bmA}Hay unsubbloque correspondiente en A , y su interior es una matriz con 1 en el triángulo superior o en el triángulo inferior (comoA {\bm A}subbloque en la esquina superior izquierda de A ). Además, en B o C, los subbloques correspondientes también tienen sus restos. A {\bm A}AB {\ bm B}segundoC {\ bm C}El orden de los subbloques en C es el mismo que λ 1 \lambda_1yo1tienen la misma multiplicidad, como en este ejemplo λ 1 \lambda_1yo1es una raíz triple, entonces A {\bm A}El subbloque en la esquina superior izquierda de A es de tercer orden, B {\bm B}segundoC {\ bm C}Los subbloques en C también son de tercer orden. En este artículo, los subbloques están separados por líneas finas para una fácil identificación.

9. Ejemplos

9.1 Partida derivada [No incluida] en la entrada del sistema

Sea la ecuación diferencial del sistema
y ( 3 ) + 2 y ¨ + 5 y ˙ + 2 y = 2 uy^{(3)} + 2 \ddot y + 5 \dot y + 2y= 2uy( 3 )+2y¨+5y˙+2 años=2 ves n= 3 n=3norte=3 , entonces sus coeficientes se pueden escribir inmediatamente:a 0 = 2, a 1 = 5, a 2 = 2, a 3 = 1, b 0 = 2 a_0 = 2, a_1 = 5, a_2 = 2, a_3 = 1; b_0 = 2a0=2 ,a1=5 ,a2=2 ,a3=1 ;b0=2。Ver (3):
A = [ 0 1 0 0 0 1 − 2 − 5 − 2 ] {\bm A} = \left[ \begin{matrix} 0 & 1 & 0 \\ 0 & 0 & 1\ \-2&-5&-2\end{matrix}\right]A= 00 210 501 2 B = [ 0 0 b 0 ] , C = [ 1 0 0 ] , D = 0 {\bm B} = \left[ \begin{matrix} 0 \\ 0 \\ b_0 \end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} 1 & 0 & 0 \end{matrix} \right], \quad {\bm D} = 0B= 00b0 ,C=[100],D=0

9.2 [contiene] elementos derivados en la entrada del sistema

Sea la ecuación diferencial del sistema
y ( 3 ) + 2 y ¨ + 5 y ˙ + 2 y = 2 u ¨ + 3 u ˙ + uy^{(3)} + 2 \ddot y + 5 \dot y + 2y= 2 \ddot u +3 \punto u + uy( 3 )+2y¨+5y˙+2 años=2tu¨+3tu˙+puedes vern = 3 n=3norte=3 , entonces sus coeficientes se pueden escribir inmediatamente:a 0 = 2, a 1 = 5, a 2 = 2, a 3 = 1; b 0 = 1, b 1 = 3, b 2 = 2, b 3 = 0 a_0 = 2, a_1 = 5, a_2 = 2, a_3 = 1; b_0 = 1, b_1 = 3, b_2 = 2, b_3 = 0a0=2 ,a1=5 ,a2=2 ,a3=1 ;b0=1 ,b1=3 ,b2=2 ,b3=0 . Sustituyendo en la fórmula (4) primero calculehi h_ihyo:
hi = bn − i − ∑ j = 1 ian − jhi − j ⟹ h 0 = bn = b 3 = 0 , h 1 = b 3 − 1 − ∑ j = 1 1 a 3 − jh 1 − j = b 2 − una 2 h 0 = 2 , h 2 = segundo 3 − 2 − ∑ j = 1 2 una 3 − jh 2 − j = segundo 1 − una 2 h 1 − una 1 h 0 = 3 − 2 × 2 = − 1 h 3 = segundo 3 − 3 − ∑ j = 1 3 un 3 − jh 3 − j = segundo 0 − un 2 h 2 − un 1 h 1 − un 0 h 0 = 1 + 2 × 1 − 5 × 2 = − 7 h_i = b_{ni} - \sum _{j=1} ^i a_{nj} h_{ij} \Longrightarrow \\ h_0 = b_n = b_3 = 0, \\ h_1 = b_{3-1} - \ suma _{j=1} ^1 a_{3-j} h_{1-j} = b_2 - a_2 h_0 = 2, \\ h_2 = b_{3-2} - \sum _{j=1} ^2 a_{3-j} h_{2-j} = b_1 - a_2 h_1 - a_1 h_0 = 3 - 2 \times 2 = -1 \\ h_3 = b_{3-3} - \sum_{j=1} ^ 3 a_{3-j} h_{3-j} = b_0 - a_2 h_2 - a_1 h_1 - a_0 h_0 = 1 + 2 \times 1 - 5 \times 2 = -7hyo=bnorte yoj = 1yoanorte - jhyo jh0=bnorte=b3=0 ,h1=b3-1 _ _j = 11a3 jh1 j=b2a2h0=2 ,h2=b3-2 _ _j = 12a3 jh2 j=b1a2h1a1h0=32×2=1h3=b3-3 _ _j = 13a3 jh3 j=b0a2h2a1h1a0h0=1+2×15×2=7 Sustituyendo en la fórmula (5):
su matriz de espacio de estados es:
A = [ 0 1 0 0 0 1 − 2 − 5 − 2 ] {\bm A} = \left[ \begin{matrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -2 & -5 & -2 \end{matrix} \right]A= 00 210 501 2 B = [ h 1 h 2 h 3 ] = [ 2 − 1 − 7 ] , C = [ 1 0 ⋯ 0 ] , D = h 0 = 0 {\bm B} = \left[ \begin{matrix} h_1 \ \ h_2 \\ h_3 \end{matrix} \right] = \left[ \begin{matrix} 2 \\ -1 \\ -7 \end{matrix} \right], \quad {\bm C} = \left [ \begin{matrix} 1 & 0 & \cdots & 0 \end{matrix} \right], \quad {\bm D} = h_0 = 0B= h1h2h3 = 21 7 ,C=[100],D=h0=0

9.3 La entrada del sistema [contiene] término derivado, pero bn = 0 {\bm b}_{\bm n} {\bm =} {\bm 0}bnorte= 0

Aún tomando la ecuación diferencial de 9.3 como ejemplo, observe que bn = b 3 = 0 b_n = b_3 = 0bnorte=b3=0 . Luego se puede sustituir en la fórmula (6):
A = [ 0 0 − 2 1 0 − 5 0 1 − 2 ] {\bm A} = \left[ \begin{matrix} 0 & 0 & -2 \\ 1 & 0 & -5 \\ 0 & 1 & -2 \\ \end{matrix} \right]A= 01000125 2 B = [ b 0 b 1 ⋮ bn − 1 ] = [ 1 3 2 ] , C = [ 0 0 ⋯ 1 ] , D = 0 (6) {\bm B} = \left[ \begin{matrix} b_0 \ \ b_1 \\ \vdots \\ b_{n-1} \end{matrix} \right] = \left[ \begin{matrix} 1 \\ 3 \\ 2 \end{matrix} \right], \quad { \bm C} = \left[ \begin{matrix} 0 & 0 & \cdots & 1 \end{matrix} \right], \quad {\bm D} = 0 \tag{6}B= b0b1bnorte 1 = 132 ,C=[001],D=0( 6 )

9.4 Tipo estándar controlable

Sea la ecuación diferencial del sistema:
y ( 3 ) + 2 y ¨ + 3 y ˙ + 3 y = 2 u ( 3 ) + u ¨ + 3 u ˙ + 4 uy^{(3)} + 2 \ddot y + 3 \punto y + 3y = 2 u^{(3)} + \ddot u + 3 \punto u + 4uy( 3 )+2y¨+3y˙+3 años=2 y( 3 )+tu¨+3tu˙+4 puedes vern = 3 , a 0 = 3 , a 1 = 3 , a 2 = 2 , a 3 = 1 ; b 0 = 4 , b 1 = 3 , b 2 = 1 , b 3 = 2 n=3 , a_0 = 3, a_1 = 3, a_2 = 2, a_3 = 1; b_0 = 4, b_1 = 3, b_2 = 1, b_3 = 2norte=3 ,a0=3 ,a1=3 ,a2=2 ,a3=1 ;b0=4 ,b1=3 ,b2=1 ,b3=2 . Sustituyendo la fórmula (7) para calcular el parámetroβ i \beta_ibyo
{ β 0 = b 0 − a 0 bn = 4 − 3 × 2 = − 2 , β 1 = b 1 − a 1 bn = 3 − 3 × 2 = − 3 , β 2 = b 2 − a 2 bn = 1 − 2 × 2 = − 3 (7) \begin{casos} \beta_0 = b_0 - a_0 b_n = 4 - 3 \times 2 = -2, \\ \beta_1 = b_1 - a_1 b_n = 3 - 3 \times 2 = -3, \\ \beta_2 = b_2 - a_2 b_n = 1 - 2 \times 2 = -3 \end{cases} \tag{7} b0=b0a0bnorte=43×2=2 ,b1=b1a1bnorte=33×2=3 ,b2=b2a2bnorte=12×2= 3( 7 ) Sustituya en la fórmula (8) para obtener:
A = [ 0 1 0 0 0 1 − 3 − 3 − 2 ] {\bm A} = \left[ \begin{matrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -3 & -3 & -2 \end{matrix} \right]A= 00 310 301 2 B = [ 0 0 1 ] , C = [ − 2 − 3 − 3 ] , D = 0 (8) {\bm B} = \left[ \begin{matrix} 0 \\ 0 \\ 1 \end{matrix } \right], \quad {\bm C} = \left[ \begin{matrix} -2 & -3 & -3 \end{matrix} \right], \quad {\bm D} = 0 \tag{ 8}B= 001 ,C=[ 2 3 3],D=0( 8 )

9.5 Forma estándar considerable

Dado que la forma canónica observable es la transpuesta de la forma canónica controlable, podemos escribir inmediatamente:
A = [ 0 0 − 3 1 0 − 3 0 1 − 2 ] {\bm A} = \left[ \begin{matrix} 0 & 0 & -3 \\ 1 & 0 & -3 \\ 0 & 1 & -2 \end{matrix} \right]A= 01000133 2 B = [ − 2 − 3 − 3 ] , C = [ 0 0 1 ] , D = 0 (8) {\bm B} = \left[ \begin{matrix} -2 \\ -3 \\ -3 \ end{matrix} \right], \quad {\bm C} = \left[ \begin{matrix} 0 & 0 & 1 \end{matrix} \right], \quad {\bm D} = 0 \tag{ 8}B= 23 3 ,C=[001],D=0( 8 )

Supongo que te gusta

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