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
- 2. La entrada del sistema [no contiene] la partida derivada
- 3. [Contiene] el elemento derivado en la entrada del sistema.
- 4. La entrada del sistema [contiene] término derivado, pero bn = 0 {\bm b}_{\bm n} {\bm =} {\bm 0}bnorte= 0
- 5. Tipo estándar controlable
- 6. Estándar considerable
- 7. El tipo equivalente de la ecuación característica de la función de transferencia [cuando no hay raíz repetida]
- 8. El tipo equivalente de ecuación característica de la función de transferencia [cuando hay múltiples raíces]
- 9. 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=X3⋮X˙norte − 1=XnorteX˙norte=− un0X1−a1X2−⋯−anorte − 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=
00⋮0− un010⋮0− un101⋮0− un2⋯⋯⋱⋯⋯00⋮1− unnorte − 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=
00⋮0b0
,C=[10⋯0],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+h2tu⋮X˙norte − 1=Xnorte+hnorte − 1tuX˙norte=− un0X1−a1X2−⋯−anorte − 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 − yo−j = 1∑yoanorte - 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=
00⋮0− un010⋮0− un101⋮0− un2⋯⋯⋱⋯⋯00⋮1− unnorte − 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=
h1h2⋮hnorte − 1hnorte
,C=[10⋯0],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=X1−a1Xnorte+b1tu⋮X˙norte − 1=Xnorte − 2−anorte − 2Xnorte+bnorte − 2tuX˙norte=Xnorte − 1−anorte − 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=
010⋮0001⋮0⋯⋯⋯⋱⋯000⋮1− un0− un1− un2⋮− unnorte − 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=
b0b1⋮bnorte − 1
,C=[00⋯1],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=X3⋮X˙norte=− un0X1−a1X2−⋯−anorte − 1Xnorte+tuy=− segundo0X1−b1X2−⋯−bnorte − 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=b0−a0bnorteb1=b1−a1bnorte⋮byo=byo−ayobnorte⋮bnorte − 2=bnorte − 2−anorte − 2bnortebnorte − 1=bnorte − 1−anorte − 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=
00⋮0− un010⋮0− un101⋮0− un2⋯⋯⋱⋯⋯00⋮1− unnorte − 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=
00⋮01
,C=[b0b1⋯bnorte − 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=
010⋮0001⋮0⋯⋯⋯⋱⋯000⋮1− un0− un1− un2⋮− unnorte − 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=
b0b1b2⋮bnorte − 1
,C=[00⋯1],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 )=( s−yo1)( s−yo2)⋯( s−yonorte)此时传递函数即:
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 )=( s−yo1)( s−yo2)⋯( s−yonorte)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 = 1∑nortes−yoyoCyo( 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=( s−yoyo)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 = 1∑nortes−yoyoCyoU ( 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 )=s−yoyo1U ( 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+tu⋮X˙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=
yo1yo2⋱yonorte − 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=
11⋮11
,C=[C1C2⋯Cnorte],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 )=s−yoyoCyoU ( 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+C2tu⋮X˙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=
yo1yo2⋱yonorte − 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=
C1C2⋮Cnorte − 1Cnorte
,C=[11⋯1],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 )=( s−yo1)3( s−yo4)⋯( s−yonorte) 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)=( s−yo1)3C11+( s−yo1)2C12+( s−yo1)C13+yo = 4∑nortes−yoyoCyo( 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=
yo11yo101yo1yo40⋱yonorte
,B=
0011⋮1
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=[C11C12C13C4⋯Cnorte],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=
yo11yo110yo1yo40⋱yonorte
,B=
C11C12C13C4⋮Cnorte
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=[0011⋯1],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_1yo1在A{\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}A、B {\ bm B}segundo、C {\ 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}segundo、C {\ 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 − yo−j = 1∑yoanorte - jhyo − j⟹h0=bnorte=b3=0 ,h1=b3-1 _ _−j = 1∑1a3 − jh1 − j=b2−a2h0=2 ,h2=b3-2 _ _−j = 1∑2a3 − jh2 − j=b1−a2h1−a1h0=3−2×2=− 1h3=b3-3 _ _−j = 1∑3a3 − jh3 − j=b0−a2h2−a1h1−a0h0=1+2×1−5×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
=
2− 1− 7
,C=[10⋯0],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=
010001− 2− 5− 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=
b0b1⋮bnorte − 1
=
132
,C=[00⋯1],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=b0−a0bnorte=4−3×2=− 2 ,b1=b1−a1bnorte=3−3×2=− 3 ,b2=b2−a2bnorte=1−2×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=
010001− 3− 3− 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=
− 2− 3− 3
,C=[001],D=0( 8 )