Capítulo 4
Respuesta temporal

La respuesta temporal de un sistema físico ante unas ciertas entradas requiere resolver la ecuación diferencial que constituye el modelo matemático de dicho sistema.

Algunos tipos de ecuaciones diferenciales (en concreto, las lineales) pueden ser resueltas explícitamente (obteniendo una fórmula en función del tiempo para la solución) y otras no, y sólo pueden ser resueltas numéricamente (integración numérica, simulación, sección 2.1).

El siguiente vídeo ilustra el uso del comando dsolve, de la Symbolic toolbox de Matlab, para obtener la solución explícita de algunos tipos de ecuaciones diferenciales.

 Resolución de ecuaciones diferenciales con Matlab : [Fernando Giménez] [Enlace]

4.1 Transformada de Laplace

El primer método que se suele explicar en cursos de grado para el cálculo de la respuesta temporal se basa en la Transformada de Laplace.

4.1.1 Función de transferencia (SISO), Matriz de transferencia (MIMO)

[98: fdtLa función de transferencia controltheoryorg (canal YouTube) ** PIC 16:03

Videos de terceros relacionados:

Aquí se enlazan algunos materiales de terceros relacionados, utilizando la Transformada de Laplace:

 Cálculo de la antitransformada de Laplace por descomposición en fracciones simples: [Macarena Boix] [Enlace]

 Transformada de Laplace con Mathematica: [Elena Sánchez] [Enlace]

 Aplicación de la Transformada de Laplace a la resolución de problemas de valores iniciales: [P. Carmen Coll] [Enlace]

4.2 Respuesta temporal de sistemas sencillos

[99: stdinpEntradas estándar para el análisis de sistemas (escalón, rampa, senoidal): definición y utilidad controltheoryorg (canal YouTube) ** PIC 11:58

4.2.1 Caso de estudio masa-muelle-amortiguador

[100: mmamexpmasa resorte amortiguador real Señales y Sistemas (canal YouTube) * PIC 03:00

[101: masmusym1Resolución ecuaciones diferenciales ordinarias (EDO) con Matlab (dsolve): respuesta libre masa-muelle ** PIC 11:30 *Link to English version

[102: masmuAnimAnimación respuesta temporal masa-muelle en Matlab, efecto del amortiguamiento *** PIC 10:36 *Link to English version

[103: masmusym2Resolución ecuaciones diferenciales ordinarias (EDO) con Matlab (dsolve): respuesta libre masa-muelle representación interna ** PIC 07:53 *Link to English version

[104: masmusymForzResolución EDO con Matlab (dsolve): respuesta forzada masa-muelle (constante + senoidal) *** PIC 12:59 *Link to English version

[105: dsolvevsode45Resolución simbólica (dsolve) versus numérica (ode45): masa muelle amortiguador, ventajas e inconvenientes *** PIC 14:10

[106: sssmlkSistema masa-muelle-amortiguador: simulación en Simulink, extracción de representación interna y FdT Eduardo Giraldo *** PIC 11:46

[107: masmuLaplResolución EDO por transformada de Laplace: ejemplo masa-muelle ante escalón *** PIC 17:59

4.2.2 Sistemas genéricos

[108: smktfSimulink: simulación de sistemas dinámicos en función de transferencia Juan M. Herrero Durá (UPV) ** PIC 06:27

[109: stepord1Respuesta ante escalón unitario de sistemas de primer orden controltheoryorg (canal YouTube) ** PIC 14:16

[110: ejord1Ejemplos de sistemas de primer orden (experimentales) Antonio Barrientos (UPM) * PIC 06:09

[111: ilaplaceex1Respuesta temporal mediante fracciones simples y tablas de transformada inversa: sistemas 1er orden ** PIC 13:53

[112: ilaplaceex2Respuesta temporal mediante fracciones simples y tablas de transformada inversa: sistemas 2o orden polos reales ** PIC 15:58

[113: ilaplaceex3Respuesta temporal mediante fracciones simples y tablas de transformada inversa: sistemas 2o orden polos complejos *** PIC 14:57

[114: ilaplaceex4Respuesta temporal mediante fracciones simples y tablas de transformada inversa: sistema 4o orden ante rampa **** PIC 12:39

4.2.3 Otros ejemplos:

[115: rcstepsinCircuito RC serie: modelado y respuesta entrada escalón/senoidal por transformada de Laplace *** PIC 18:48

[116: fugeqlinDinámica fugoide aeronave: equilibrio, linealización, estabilidad *** PIC 19:30 *Link to English version

4.2.4 Sistemas con retardo en dinámica, o entradas “por tramos”

[117: retardSistemas con retardo: introducción y aproximación Antonio Barrientos (UPM) ** PIC 13:17

Sistema de primer orden con retardo

[118: dly1er1Respuesta temporal sistema primer orden con retardo de transporte: dsolve, discusión sobre resultados ** PIC 05:57

[119: dly1er2Respuesta temporal sistema primer orden con retardo: resolución por transformada de Laplace *** PIC 12:29

Circuito de 1er orden ante tren de pulsos senoidales

[120: sinpulRCRRespuesta temporal circuito RCR ante pulso senoidal y condiciones iniciales no nulas (1: Laplace, superposición) ** PIC 13:51 *Link to English version

[121: sinpulRCR2Respuesta temporal circuito RCR ante pulso senoidal y condiciones iniciales no nulas (2: por tramos) *** PIC 12:43 *Link to English version

[122: trenpulRCRRespuesta temporal de un circuito RCR ante tren de pulsos senoidales **** PIC 15:45 *Link to English version

Más ejemplos de formas de onda por tramos

[123: pwlap1Respuesta temporal con entrada a tramos: planteamiento general, ejemplo tren escalones (superposición) ** PIC 12:53

[124: pwlap2Respuesta temporal: entrada a tramos, ejemplo tren de rampas ** PIC 09:09

[125: pwlap3Respuesta temporal entrada a tramos, ejemplo rampas y escalones combinados *** PIC 08:35

Respuesta temporal ante entradas definidas por tramos (piecewise)

La transformada de Laplace de entradas definidas “por tramos” incorpora retardos en cada “cambio de fórmula” de u(t), como en el siguiente ejemplo:

[126: sinpulLTransformada de Laplace de un pulso senoidal (semiperíodo) ** PIC 11:33 *Link to English version

En ocasiones, la transformada de Laplace de una entrada piecewise nunca es “explícitamente” utilizada, porque la salida se puede obtener por superposición. El siguiente video ejemplifica el caso de un pulso rectangular de duración finita.

[127: dly1erpulRespuesta temporal primer orden inestable ante pulso rectangular (transformada de Laplace) *** PIC 19:38

La idea también puede aplicarse a otras formas de onda, como rampas:

[128: dlyrampRespuesta temporal sistema segundo orden ante rampa truncada, por transf. Laplace; comparación con escalón **** PIC 15:58

4.2.5 Matriz de transferencia (MIMO)

[129: mdtRepresentación en Matriz de Transferencia de sistemas multivariables ** PIC 10:27

4.3 Teorema del Valor Final

[130: tvf1Teorema valor final (Laplace) [1]: enunciado, uso y ejemplo sencillo *** PIC 08:40

[131: tvf2Teorema valor final (Laplace) [2]: enunciado y demostración *** PIC 07:28

[132: tvf3Teorema valor final (Laplace) [3]: relación con ganancia y ecuaciones en equilibrio *** PIC 14:45

[133: tvf4Teorema valor final (Laplace) [4]: ejemplos Matlab + generalización a entradas polinomiales *** PIC 15:30

4.4 Recapitulación: respuesta de sistemas sencillos ante entradas típicas

4.4.1 Sistemas de primer orden (+retardo)

[134: ord1teoRespuesta sistemas primer orden ante escalón: ganancia, constante de tiempo (teoría). ** PIC 14:49

[135: ord1trampRespuesta sistemas primer orden ante rampa en función de ganancia y constante de tiempo: teoría + ejemplo Matlab ** PIC 14:31

4.4.2 Sistemas de segundo orden

[136: ord2motiRespuesta escalón sistemas de segundo orden: motivación, caso polos reales, coeficiente de amortiguamiento *** PIC 13:48

[137: ord2stepRespuesta escalón sistemas segundo orden subamortiguados: sobreoscilación, frecuencia propia, tiempo pico, … *** PIC 19:01

[138: paramsor2Sistema segundo orden subamortiguado: efecto K, wn, wp, sigma, amortiguamiento (simulación/animación) ** PIC 19:56

4.4.3 Efecto de los ceros en la dinámica de los sistemas

[139: zeros1Efecto de ceros en respuesta de un sistema: inversión de dinámica, superposición con derivada, ceros en el origen *** PIC 15:25

[140: zerosgrCeros de una función de transferencia: grado relativo, efecto en respuesta temporal y en frecuencia **** PIC 15:32

 Efectos de la adición de un cero: [14:33] [Enlace]

[141: zerosord1Efecto de ceros en respuesta temporal: sistemas de primer orden con cero adicional (fase mínima y no mínima) ** PIC 12:47

[142: expPZRespuesta experimental circuito polo-cero (RC+C) Antonio Barrientos (UPM) ** PIC 03:40

[143: zerosord2Efecto de ceros en respuesta temporal: sistema segundo orden subamortiguado con cero adicional (fase mínima y no mínima) ** PIC 09:30

Sistemas de orden superior:

Los polos y ceros asociados a constantes de tiempo más rápidas pueden ser grosso modo despreciados para hacerse una idea cualitativa de la respuesta temporal.

[144: ordsupej1Análisis de respuesta escalón (aproximada) orden superior: ejemplos Matlab (1: polos reales) ** PIC 13:48

[145: ordsupej2Análisis de respuesta escalón (aproximada) orden superior: ejemplos Matlab (2: mezcla reales/complejos) *** PIC 13:50

[146: ordsupej3Análisis de respuesta escalón (aproximada) orden superior: ejemplos Matlab (3: ceros, casos sin dominancia clara) *** PIC 13:54

Este tipo de simplificaciones es la base de la reducción modal y de los conceptos de dominancia y cancelación discutidos en la Sección 9.1; también pueden aproximarse a un retraso (vídeos [ fopd(10:40)] y [ ordnafopd(10:25)]). Los detalles sobre todos esos enfoques a la reducción/simplificación de modelos de sistemas dinámicos no son ojetivo de este capítulo.

Régimen estacionario ante rampa de una FdT

[147: regestramp1Respuesta estacionaria ante rampa de sistemas genéricos en función de transferencia (1): teoría *** PIC 11:55

[148: regestramp2Respuesta estacionaria ante rampa de sistemas genéricos en función de transferencia (2): interpretación y ejemplos *** PIC 13:18

Caso de estudio linealidad + invarianza temporal + superposición

[149: linregla3Utilización intuitiva de dinámica lineal: tren de escalones con regla de tres, control bucle abierto * PIC 15:26 *Link to English version

[150: linregla3ord2Utilización intuitiva linealidad: tren de escalones regla de tres, proceso segundo orden (no va bien) *** PIC 08:16 *Link to English version

4.4.4 Identificación experimental basada en respuesta ante entradas típicas

Sistemas de primer orden + retardo

[151: fopdModelos primer orden + retardo: propiedades resp. escalón e identificación experimental (manual) *** PIC 10:40

[152: ord1exgIdentificación experimental de sistemas de 1er orden (+retardo): tres ejemplos que salen perfectos ** PIC 10:57

[153: ord1exmalIdentificación experimental de sistemas de 1er orden (+retardo): ejemplos no tan perfectos *** PIC 17:36

Sistemas de segundo orden

[154: ord2idgRespuesta escalón sistemas de segundo orden subamortiguado: identificación experimental (ejemplo 1) *** PIC 12:49

[155: ord2idbRespuesta escalón segundo orden subamortiguado: identificación experimental (2), comparación con procest *** PIC 16:32

Identificación en bucle cerrado (introducción)

[156: idbcind1Identificación en bucle cerrado: introducción, planteamiento del problema, identificación directa vs. indirecta *** PIC 11:52

[157: idbcind2Identificación indirecta en bucle cerrado ante ensayo escalón en referencia, ejemplo eje *** PIC 14:49

[158: idbcind3Identificación indirecta en bucle cerrado: ensayo escalón en acción control, ejemplo eje *** PIC 13:02

4.5 Respuesta temporal en representación interna

4.5.1 Respuesta libre

[159: rtss1Resolución dx=Ax: exponencial de matriz ** PIC 11:00

[160: rtssexplFórmulas explícitas de la exponencial de una matriz (via forma de Jordan) *** PIC 10:32

[161: rtssmlRespuesta temporal (exponencial de matriz), ejemplo Matlab ** PIC 10:38

4.5.2 Régimen Estacionario

[162: ressRégimen estacionario de sistemas en representación en variables de estado (repr. interna): constante/rampa/senoidal *** PIC 10:57

[163: ressmlRégimen estacionario constante/rampa/senoidal, repr. interna: ejemplo Matlab *** PIC 07:02

4.6 Ejemplos de modelado y cálculo de respuesta temporal

[164: rcrmlModelado de un circuito R-C-R con dos fuentes de tensión: matriz de transferencia, término de condiciones iniciales, resp. escalón *** PIC 08:43

[165: rcrssmlModelado/simulación de un circuito R-C-R con dos fuentes de tensión: representación interna *** PIC 07:49

[166: ep1Sistema engranajes y poleas (lineal). Modelado y Función de transferencia ** PIC 19:37

[167: ep2Sistema engranajes y poleas (lineal): representación en variables de estado ***** PIC 18:20

[168: moll3freeSistema 3 masas 4 muelles: análisis modal de la respuesta libre *** PIC 14:50 *Link to English version

4.6.1 Caso de estudio: motor de corriente continua con carga mecánica

[169: motccmodelModelado motor corriente continua haciendo girar una carga (representación en variables de estado) ** PIC 13:59

[170: motengModelado motor CC con reductora *** PIC 18:54

[171: motcctfstepMotor corriente continua haciendo girar carga: función de transferencia y respuesta escalón (step) *** PIC 14:58

[172: motccLaplMotor corriente continua que gira carga: respuesta ante escalón (método Laplace) y tren de escalones *** PIC 15:53

4.6.2 Caso de estudio: calentamiento pieza de dos capas

En elaboración...

[173: term2mod1Modelado de calentamiento pieza de 2 capas: representación interna (espacio de estado) ** PIC 13:35

[174: term2mod2Modelado de calentamiento pieza de 2 capas, opciones alternativas: 1er orden, EDP, elementos finitos *** PIC 10:58

[175: term2stepCalentamiento pieza de 2 capas: función (matriz) de transferencia, respuesta escalón control toolbox *** PIC 15:44

[176: term2rampCalentamiento pieza de 2 capas: tratamiento térmico en rampa para limitar diferencia de temperaturas **** PIC 17:51

4.6.3 Caso de estudio: transformador eléctrico

[177: tni1Modelado de un transformador eléctrico no ideal: ecuaciones físicas *** PIC 09:37

[178: tni2Representación transformada de Laplace (func. de transferencia) y respuesta temporal de transformador eléctrico no ideal *** PIC 06:16

[179: tni3Modelado de un transformador eléctrico no ideal: representación en variables de estado *** PIC 07:45

El estudio del régimen senoidal (corriente alterna) de este transformador se aborda en el vídeo [ tni4(09:21)] en la Sección 5.1.

4.6.4 Caso de estudio: calentador tubular

Una primera aproximación al modelado se vio en la sección 1.2.4.

[180: term1ePtfModelado dinámico de tanque de calentamiento como sistema de primer orden: linealización y Función de Transferencia (2 perfiles diferentes) *** PIC 12:58

[181: termedpsolEcs. derivadas parciales de calentador tubular de líquido: solución EDP con método de Laplace ***** PIC 19:57 *Link to English version

[182: termedpstepRespuesta escalón de solución (función transferencia) a la EDP de un calentador tubular (Matlab) *** PIC 09:58 *Link to English version

[183: term1evsedpComparación entre solución EDP y diferentes aproximaciones de orden 1 y 3 de un calentador tubular **** PIC 16:55 *Link to English version

[184: tubulFE1Calentador/intercambiador tubular: modelado elementos finitos (I) ***** PIC 14:47 *Link to English version

[185: tubulFEsim1Calentador tubular, elementos finitos: simulación numérica (1) **** PIC 14:39 *Link to English version

[186: tubulFEsim2Calentador tubular, elementos finitos: simulación numérica (2) ***** PIC 11:38 *Link to English version

[187: tubulconcCalentador tubular: recapitulación y conclusiones sobre complejidad de modelos *** PIC 14:55 *Link to English version

4.6.5 Sistemas con retardo

[188: aproxretAproximación del retardo: simulación y comparación de opciones (PDE, Padé, elementos finitos) *** PIC 15:59 *Link to English version