UNIVERSIDAD AUTÓNOMA DE NUEVO LEÓN
FACULTAD DE INGENIERÍA MECÁNICA ELÉCTRICA
Laboratorio de identificación de sistemas Práctica N°7:
Identificación por Subespacios.
Catedrático: Dr. Luis Humberto Rodríguez Alfaro
Matrícula Nombre Carrera
1992242 Mauricio Andrés Treviño López IEA
1802000 Christian Axel De La Cruz Puente IEA
1843266 André Marín Almaguer Guerra IEA
Hora: N1 Día: Martes Grupo: 237
Semestre: Agosto - diciembre 2023
Ciudad Universitaria, N.L, a 21 de noviembre del 2023
Objetivo: Aplicar los métodos de identificación por subespacios disponibles en el toolbox de
identificación de Sistemas de Scilab.
Desarrollo de la práctica.
Utilice las señales de entrada U y salida Y proporcionadas por el instructor, estas señales
están definidas por 1000 muestras agrupadas en un vector y cada una de estas señales esta
guardada en un archivo.
Para estas señales de entrada salida realice lo siguiente:
i) Procesamiento de las señales de entrada salida.
a) Gráfique las señales de entrada salida.
Y:
U:
b) Aplique un comando de procesamiento de señales que permita eliminar la
deficiencia que aparece en estas señales.
c) Grafique las señales resultantes de aplicar el comando de procesamientode
señales.
U procesada:
Y procesada:
d) Defina cual es la deficiencia de estas señales y las diferencias entre lasseñales
procesadas y no procesadas.
Como se puede apreciar en las gráficas entrada salida nuestras señales muestran niveles de
continua que en la salida (Y) es de magnitud 2 aproximadamente mientras que en la entrada
(U) es de 1, con el comando “detrend” eliminamos estosniveles de continúa haciendo que
nuestras señales procesadas estén en el nivel demagnitud 0.
ii) Obtención de modelos estimados en tiempo discreto utilizando las
metodologías MOESP y N4sid, para cada uno de los modelos obtenidos
reporte lo siguiente:
a) Las matrices del sistema.
MOESP:
N4sid:
b) La ecuación de estado y la ecuación de salida que forman al modelo
estimado.
𝑥̇ = 𝐴𝑥̇ + 𝐵𝑢
𝑦 = 𝐶𝑥̇ + 𝐷𝑢
MOESP:
𝑥̇1 −0.98.3339 −0.213479 −0.2878645 𝑥̇1 −0.5877213
[𝑥̇2 ] = [ 0.0296784 0.7076414 −0.4148243] [𝑥̇2 ] + [ 0.6626815 ] 𝑈
𝑥̇3 −0.0189488 0.1876472 −0.7135383 𝑥̇3 1.5187575
𝑥̇1
𝑦 = [−0.3544169 0.645801 0.5585025] [𝑥̇2] + [0.6082591]𝑈
𝑥̇3
N4sid:
𝑥̇1 0 0 0 𝑥̇1 −0.2864522
[𝑥̇2 ] = [0 0 0] [𝑥̇2 ] + [ 0.0460335 ] 𝑢
𝑥̇3 0 0 0 𝑥̇3 −0.1118977
𝑥̇1
𝑦 = [0 0 0] [𝑥̇2] + [0.1072546]𝑢
𝑥̇3
c) La respuesta en el tiempo que resulta utilizando el modelo estimado.
MOESP:
N4sid:
iii) Comparación. Para cada uno de los modelos obtenga:
a) La gráfica que compara la señal de salida real con la señal de salida
obtenida a partir del modelo estimado.
MOESP:
N4sid:
b) Defina las diferencias entre la señal de salida estimada y la real obtenidas en a).
Puede ver la diferencia en el rango de mediciones realizadas. Por ejemplo, para el espectro
N4sid, puede ver una línea azul que representa una salida real mucho más pequeña en
comparación con el gráfico obtenido con el comando MOESP. Esto es claramente similar.
Como una luz verde. Al medirlo verás que hay una gran diferencia en las medidas.
iv)Defina con cuál de los modelos se obtiene un mejor resultado y expliquelas
razones.
Como se ve en el gráfico, los resultados reales y estimados son muy similares para el
modelo MOESP, pero los resultados reales y estimados son muy diferentes de los
resultados estimados para el modelo N4sid.