canon

(No se recomienda) Realización canónica en espacio de estados

canon no se recomienda. En su lugar, utilice modalreal y compreal. (desde R2023b) Para más información sobre cómo actualizar el código, consulte Historial de versiones.

Sintaxis

csys = canon(sys,type)

csys = canon(sys,'modal',condt)

[csys,T]= canon(___)

Descripción

csys = canon(sys,type) transforma el modelo lineal sys en la realización canónica csys. type especifica si csys está en forma modal o compañera.

Para obtener información sobre formas canónicas controlables observables, consulte Realizaciones en espacio de estados.

csys = canon(sys,'modal',condt) especifica un límite superior condt en el número de condición de la transformación diagonal de bloques. Utilice condt si tiene valores propios cercanos en csys.

[csys,T]= canon(___) también devuelve la matriz de transformación de coordenadas de estado T que relaciona los estados del modelo de espacio de estados sys con los estados de csys.

Ejemplos

Convertir un modelo de espacio de estados en una forma canónica compañera

Supongamos que tiene un modelo de espacio de estados de una aeronave donde la entrada es el ángulo de deflexión del elevador $\delta$ y la salida es el ángulo de cabeceo de la aeronave $\theta$.

$$\begin{array}{l}\left[\begin{array}{c}\underset{}{\overset{\dot{}}{\alpha }}\\ \underset{}{\overset{\dot{}}{q}}\\ \underset{}{\overset{\dot{}}{\theta }}\end{array}\right]=\left[\begin{array}{ccc}-0.313& 56.7& 0\\ -0.0139& -0.426& 0\\ 0& 56.7& 0\end{array}\right]\left[\begin{array}{c}\alpha \\ q\\ \theta \end{array}\right]+\left[\begin{array}{c}0.232\\ 0.0203\\ 0\end{array}\right]\left[\delta \right]\\ y=\left[\begin{array}{ccc}0& 0& 1\end{array}\right]\left[\begin{array}{c}\alpha \\ q\\ \theta \end{array}\right]+\left[0\right]\left[\delta \right]\end{array}$$

Cargue los datos del modelo en el área de trabajo y cree el modelo de espacio de estados sys.

load('aircraftPitchSSModel.mat');
sys = ss(A,B,C,D)

sys =
 
  A = 
            x1       x2       x3
   x1   -0.313     56.7        0
   x2  -0.0139   -0.426        0
   x3        0     56.7        0
 
  B = 
           u1
   x1   0.232
   x2  0.0203
   x3       0
 
  C = 
       x1  x2  x3
   y1   0   0   1
 
  D = 
       u1
   y1   0
 
Continuous-time state-space model.

Convierta el modelo de espacio de estados sys resultante en la forma compañera controlable.

csys = canon(sys,'companion')

csys =
 
  A = 
              x1         x2         x3
   x1          0          0  1.914e-15
   x2          1          0    -0.9215
   x3          0          1     -0.739
 
  B = 
       u1
   x1   1
   x2   0
   x3   0
 
  C = 
            x1       x2       x3
   y1        0    1.151  -0.6732
 
  D = 
       u1
   y1   0
 
Continuous-time state-space model.

csys es la forma compañera controlable de sys.

Convertir un modelo de espacio de estados en una forma modal

Supongamos que tiene un modelo de espacio de estados para un péndulo invertido en un carro donde las salidas son el desplazamiento del carro $\mathit{x}$ y el ángulo del péndulo $\theta$. La entrada de control u es la fuerza horizontal sobre el carro.

$$\begin{array}{l}\left[\begin{array}{c}\underset{}{\overset{\dot{}}{x}}\\ \underset{}{\overset{¨}{x}}\\ \underset{}{\overset{\dot{}}{\theta }}\\ \underset{}{\overset{¨}{\theta }}\end{array}\right]=\left[\begin{array}{cccc}0& 1& 0& 0\\ 0& -0.1& 3& 0\\ 0& 0& 0& 1\\ 0& -0.5& 30& 0\end{array}\right]\left[\begin{array}{c}x\\ \underset{}{\overset{\dot{}}{x}}\\ \theta \\ \underset{}{\overset{\dot{}}{\theta }}\end{array}\right]+\left[\begin{array}{c}0\\ 2\\ 0\\ 5\end{array}\right]u\\ y=\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& 0& 1& 0\end{array}\right]\left[\begin{array}{c}x\\ \underset{}{\overset{\dot{}}{x}}\\ \theta \\ \underset{}{\overset{\dot{}}{\theta }}\end{array}\right]+\left[\begin{array}{c}0\\ 0\end{array}\right]u\end{array}$$

Primero, cargue el modelo de espacio de estados sys en el área de trabajo.

load('pendulumCartSSModel.mat','sys');

Convierta sys a una forma modal y extraiga los tamaños de los bloques.

[csys,T] = canon(sys,'modal')

msys =
 
  A = 
           x1      x2      x3      x4
   x1       0       0       0       0
   x2       0   -0.05       0       0
   x3       0       0  -5.503       0
   x4       0       0       0   5.453
 
  B = 
          u1
   x1  1.875
   x2  6.298
   x3   12.8
   x4  12.05
 
  C = 
              x1         x2         x3         x4
   y1         16     -4.763  -0.003696   0.003652
   y2          0   0.003969   -0.03663    0.03685
 
  D = 
       u1
   y1   0
   y2   0
 
Continuous-time state-space model.

T = 4×4

    0.0625    1.2500   -0.0000   -0.1250
         0    4.1986    0.0210   -0.4199
         0    0.2285  -13.5873    2.4693
         0   -0.2251   13.6287    2.4995

csys es la forma canónica modal de sys, mientras que T representa la transformación entre los vectores de estado de sys y csys.

Argumentos de entrada

sys — Sistema dinámico

modelo de sistema dinámico

Sistema dinámico, especificado como un modelo de sistema dinámico SISO o MIMO. Se admiten los siguientes tipos de sistemas dinámicos:

• Modelos LTI numéricos en tiempo continuo o en tiempo discreto, como modelos tf, zpk, ss o pid.

• Modelos LTI generalizados o con incertidumbre, como modelos genss o uss (Robust Control Toolbox). El uso de modelos con incertidumbre requiere Robust Control Toolbox™.

  El modelo de espacio de estados canónico resultante asume

  • valores actuales de los componentes ajustables para los bloques de diseño de control ajustables.

  • valores nominales del modelo para los bloques de diseño de control con incertidumbre.

• Modelos LTI identificados, como modelos idtf (System Identification Toolbox), idss (System Identification Toolbox), idproc (System Identification Toolbox), idpoly (System Identification Toolbox) y idgrey (System Identification Toolbox). El uso de modelos identificados requiere System Identification Toolbox™.

No se pueden utilizar modelos de datos de respuesta en frecuencia, como modelos frd.

type — Tipo de transformación

'modal' (predeterminado) | 'companion'

Tipo de transformación, especificado como 'modal' o 'companion'. Si no se especifica type, canon convierte el modelo de sistema dinámico especificado en la forma canónica modal de manera predeterminada.

La forma canónica compañera es la misma que la forma canónica observable. Para obtener información sobre formas canónicas controlables observables, consulte Realizaciones en espacio de estados.

• Forma modal

  En la forma modal, A es una matriz de diagonales de bloques. El tamaño del bloque normalmente es de 1 por 1 para los valores propios reales y de 2 por 2 para los valores propios complejos. Sin embargo, si hay valores propios repetidos o clusters de valores propios cercanos, el tamaño del bloque puede ser mayor.

  Por ejemplo, para un sistema con valores propios $$({\lambda }_1,\sigma \pm j\omega ,{\lambda }_2)$$, la matriz modal A presenta el siguiente formato:

  $$A_m=\left[\begin{array}{cccc}{\lambda }_1& 0& 0& 0\\ 0& \sigma & \omega & 0\\ 0& -\omega & \sigma & 0\\ 0& 0& 0& {\lambda }_2\end{array}\right].$$

• Forma compañera

  En la realización compañera, el polinomio característico del sistema aparece explícitamente en la columna situada más a la derecha de la matriz A. Para un sistema con polinomio característico

  $$P(s)=s^n+{\alpha }_{n-1}{s}^{n-1}+{\alpha }_{n-2}{s}^{n-2}+\dots +{\alpha }_{1}s+{\alpha }_0,$$

  la matriz compañera correspondiente A es

  $$A_{ccom}=\left[\begin{array}{ccc}\begin{array}{l}0\\ 1\\ 0\\ 0\\ ⋮\\ 0\end{array}& \begin{array}{l}0\\ 0\\ 1\\ 0\\ ⋮\\ 0\end{array}& \begin{array}{ccc}\begin{array}{cc}\begin{array}{l}0\\ 0\\ 0\\ 1\\ ⋮\\ 0\end{array}& \begin{array}{l}\dots \\ \dots \\ \dots \\ \dots \\ \ddots \\ \dots \end{array}\end{array}& \begin{array}{l}0\\ 0\\ 0\\ 0\\ ⋮\\ 1\end{array}& \begin{array}{l}-{\alpha }_0\\ -{\alpha }_1\\ -{\alpha }_2\\ -{\alpha }_3\\ ⋮\\ -{\alpha }_{n-1}\end{array}\end{array}\end{array}\right],B_{ccom}=\left[\begin{array}{c}1\\ 0\\ ⋮\\ 0\end{array}\right].$$

  La transformación compañera exige que el sistema sea controlable desde la primera entrada. La transformación en forma compañera se basa en la matriz de controlabilidad, que casi siempre es numéricamente singular para los órdenes medios. Por tanto, evite utilizarla siempre que sea posible.

  La realización compañera devuelta por canon a veces se conoce como forma canónica controlable. Para obtener más información sobre formas canónicas controlables y observables, consulte Realizaciones en espacio de estados.

condt — Límite superior del número de condición de la transformación diagonal de bloques

1e4 (predeterminado) | escalar positivo

Límite superior del número de condición de la transformación diagonal de bloques, especificado como un escalar positivo. Este argumento solo está disponible cuando type se establece en 'modal'.

Aumente condt para reducir el tamaño de los grupos de valores propios en la matriz A de csys. Si se establece condt = Inf, se diagonaliza la matriz A.

Argumentos de salida

csys — Forma canónica del espacio de estados del modelo dinámico

Objeto de modelo ss

Forma canónica del espacio de estados del modelo dinámico, devuelta como un objeto de modelo ss. csys es una realización en espacio de estados de sys en la forma canónica especificada por type.

T — Matriz de transformación

matriz

Matriz de transformación, devuelta como una matriz de n por n, donde n es el número de estados. T es la transformación entre el vector de estados x del modelo de espacio de estados sys y el vector de estados x_c de csys:

x_c = Tx

.

Este argumento solo está disponible cuando sys es un objeto de modelo ss.

Limitaciones

• No se pueden utilizar modelos de datos de respuesta en frecuencia para convertir a la forma canónica de espacio de estados.

• La forma compañera está mal condicionada para la mayoría de los cálculos del espacio de estados; es decir, la transformación a la forma compañera se basa en la matriz de controlabilidad, que casi siempre es numéricamente singular para los órdenes medios. Por tanto, evite utilizarla siempre que sea posible.

Algoritmos

El comando canon utiliza el comando bdschur para convertir sys en una forma modal y calcular la transformación T. Si sys no es un modelo de espacio de estados, canon primero lo convierte al espacio de estados usando ss.

La reducción a la forma compañera utiliza una transformación de similitud de estados basada en la matriz de controlabilidad [1].

Referencias

[1] Kailath, T. Linear Systems, Prentice-Hall, 1980.

Historial de versiones

Introducido antes de R2006a

R2023b: No recomendado

El comando canon no se recomienda. En su lugar, utilice la funcionalidad descrita en esta tabla.

Realización	Comando antiguo	Comando nuevo
Modal	sysT = canon(sys,'modal')	sysT = modalreal(sys)
Compañera	sysT = canon(sys,'companion')	sysT = compreal(sys)

Para obtener más información sobre los comandos nuevos, consulte modalreal y compreal.

R2022a: El valor predeterminado del argumento de entrada condt ha cambiado

El valor predeterminado del argumento de entrada condt ahora es 1e4. Anteriormente, el valor predeterminado era 1e8. Aumentar condt reduce el tamaño de los grupos de valores propios en el modelo transformado, pero también disminuye la precisión de la transformación. El nuevo valor predeterminado equilibra el condicionamiento y la precisión de la transformación. Si tiene código que se basa en el valor predeterminado de condt establecido como 1e8, actualice el código para establecer explícitamente este argumento de entrada.

Consulte también

ctrb | ctrbf | ss2ss | tf | zpk | ss | pid | genss | uss (Robust Control Toolbox) | idtf (System Identification Toolbox) | idss (System Identification Toolbox) | idproc (System Identification Toolbox) | idpoly (System Identification Toolbox) | idgrey (System Identification Toolbox)

Temas

• Realizaciones en espacio de estados