📚 Fundamentación de Transformadas de Laplace en Sistemas de Control de Microgrids
Introducción
El control de frecuencia en microgrids con alta penetración de fuentes de energía renovable (RES) requiere herramientas matemáticas avanzadas para modelar, analizar y diseñar estrategias de control robustas. La transformada de Laplace se constituye como la herramienta analítica fundamental que permite transformar las ecuaciones diferenciales del dominio del tiempo al dominio de la frecuencia compleja (plano s), facilitando el análisis de estabilidad, el diseño de controladores y la evaluación del desempeño del sistema.
1. Problemática del Sistema y Solución con Laplace
❌ Problemática
- Ecuaciones diferenciales acopladas que describen la dinámica del microgrid
- Múltiples variables de estado interconectadas
- Presencia de retardos de comunicación variables en el tiempo
- Perturbaciones estocásticas de RES y cargas
- Ataques de inyección de datos falsos (FDI)
- Condiciones iniciales no nulas
✅ Solución con Laplace
- Transforma EDOs en ecuaciones algebraicas en s
- Permite obtener funciones de transferencia
- Facilita el análisis de polos y ceros
- Ubicación regional de polos (RPP) en el plano s
- Especificación de desempeño H∞ en frecuencia
- Inclusión automática de condiciones iniciales
2. Modelado Matemático del Sistema
Ecuación (1) - Modelo en Espacio de Estado (Dominio del Tiempo)
y₁(t) = y(t) + v₁(t)
y₂(t) = y(t) + v₂(t)
donde y(t) = Cx(t)
Variables del Sistema:
Vector de estado: desviación de frecuencia, potencia mecánica, potencia BESS, integral de frecuencia, potencia
Señales de control secundario y auxiliar
Perturbaciones: generación eólica, solar y desviación de carga
Salidas medidas: desviación de frecuencia y su integral
Aplicando ℒ{·} a la ecuación (1) con condiciones iniciales x(0):
sX(s) - x(0) = AX(s) + BU(s) + EW(s)
Reorganizando: (sI - A)X(s) = x(0) + BU(s) + EW(s)
X(s) = (sI - A)⁻¹x(0) + (sI - A)⁻¹BU(s) + (sI - A)⁻¹EW(s)
La matriz (sI - A)⁻¹ es la matriz de transferencia del sistema y sus polos son las raíces de det(sI - A) = 0.
3. Matrices del Sistema en el Dominio de Laplace
Matrices A, B, E, C del Modelo
⎢ 0 -1/Tₜ 1/Tₜ 0 0 ⎥
⎢ -1/RTg 0 -1/Tg 0 0 ⎥
⎢ 0 0 0 -1/Tbess 0 ⎥
⎣ 0 1 0 0 0 ⎦
B = ⎡ 0 0 ⎤
⎢ 0 0 ⎥
⎢1/Tg 0 ⎥
⎢ 0 Kbess/Tbess ⎥
⎣ 0 0 ⎦
E = ⎡ 1/2H 1/2H -1/2H ⎤
⎢ 0 0 0 ⎥
⎢ 0 0 0 ⎥
⎢ 0 0 0 ⎥
⎣ 0 0 0 ⎦
C = [1 0 0 0 0]
[0 0 0 0 1]
• D = 0.015 pu MW/Hz (coeficiente de amortiguamiento)
• H = 0.083 pu MW s (constante de inercia)
• Tg = 0.1 s (constante de tiempo del gobernador)
• Tₜ = 0.4 s (constante de tiempo de la turbina)
• Tbess = 0.1 s (constante de tiempo del BESS)
• R = 2.4 Hz/pu MW (característica de caída de velocidad)
4. Diseño del Control Auxiliar con Ubicación Regional de Polos (RPP)
Ecuación (2) - Sistema en Lazo Cerrado con Retroalimentación de Estado
y(t) = Cx(t)
Transformada de Laplace:
X(s) = (sI - A - B₂K₂)⁻¹E W(s)
Función de transferencia:
G(s) = C(sI - A - B₂K₂)⁻¹E
Ecuaciones (4)-(7) - Restricciones RPP en el Plano s
AZ₂ + Z₂Aᵀ + B₂Q̄₂ + Q̄₂ᵀB₂ᵀ + 2αₚₚZ₂ < 0
(5) Frecuencia natural máxima:
⎡ -γₚₚZ₂ AZ₂ + B₂Q̄₂
⎢ ⎥ < 0
⎣ Z₂Aᵀ + Q̄₂ᵀB₂ᵀ -γₚₚZ₂ ⎦
(6) Amortiguamiento mínimo:
⎡ sinδₚₚ(AZ₂+Z₂Aᵀ+B₂Q̄₂+Q̄₂ᵀB₂ᵀ) cosδₚₚ(AZ₂-Z₂Aᵀ+B₂Q̄₂-Q̄₂ᵀB₂ᵀ) ⎤
⎢ ⎥ < 0
⎣ cosδₚₚ(Z₂Aᵀ+Q̄₂ᵀB₂ᵀ-AZ₂-B₂Q̄₂) sinδₚₚ(AZ₂+Z₂Aᵀ+B₂Q̄₂+Q̄₂B₂ᵀ)
(7) Región de polos ℛₘᵢ:
ā < -αₚₚ < 0
‖ā + jb̄‖ < γₚₚ
ā tanδₚₚ < -|b̄|
• αₚₚ = 3: Garantiza que todos los polos tengan parte real Re(s) < -3, asegurando decaimiento exponencial rápido
• γₚₚ = 5: Limita la frecuencia natural no amortiguada ωₙ < 5 rad/s
• δₚₚ = 0.8: Asegura factor de amortiguamiento ζ ≥ cos(0.8) ≈ 0.697
Polos obtenidos: -3.3737 ± 3.1838i, -3.5460 ± 1.4927i, -3.6912
Respuesta temporal: La transformada inversa de Laplace produce:
Δf(t) = Σ Aᵢe^(σᵢt)sin(ωᵢt + φᵢ) donde σᵢ < -3
5. Desempeño H∞ en el Dominio de la Frecuencia
Ecuación (3) - Condición LMI para Norma H∞
⎢ ⎥ ≤ 0
⎣ CZ₂ -γ²I ⎦
donde Q̄₂ = K₂Z₂
La norma H∞ de la función de transferencia se define como:
‖G(s)‖∞ = sup σ̄(G(jω))
donde σ̄ denota el valor singular máximo y s = jω (eje imaginario).
La condición garantiza que:
‖y‖₂ < γ‖w‖₂
es decir, la energía de la salida está acotada por γ veces la energía de la perturbación.
En frecuencia: Esto restringe el pico de resonancia de |G(jω)| para todas las frecuencias ω ∈ ℝ.
6. Estimador de Kalman y Observador
Ecuación (9) - Estimador de Kalman en el Dominio del Tiempo
Transformada de Laplace:
sX̂(s) = AX̂(s) + B₂U₂(s) + L₂[Y(s) - CX̂(s)]
(sI - A + L₂C)X̂(s) = B₂U₂(s) + L₂Y(s)
X̂(s) = (sI - A + L₂C)⁻¹B₂U₂(s) + (sI - A + L₂C)⁻¹L₂Y(s)
Ganancia óptima: L₂ = PCᵀV⁻¹
Ecuación de Riccati:
AP + PAᵀ + EWEᵀ - PCᵀV⁻¹CP = 0
7. Control Secundario Observer-based con Event-Triggering
Ecuación (12) - Controlador Observer-based (sin ataque)
u₁(t) = K₁x̂(sₖh), t ∈ [sₖh + τₛₖ, sₖ₊₁h + τₛₖ₊₁)
donde A_c = A + B₂K₂
El retardo τₛₖ introduce el término e^(-sτₛₖ) en el dominio de Laplace,
convirtiendo al sistema en dimensión infinita.
Por esta razón, se emplea el funcional de Lyapunov-Krasovskii en lugar de análisis directo en el dominio s, permitiendo derivar condiciones LMI que garantizan estabilidad a pesar de los retardos variables.
8. Modelado de Ataques FDI en el Dominio Discreto
Ecuaciones (22)-(23) - Sistema bajo Ataque FDI
u₁(t) = K₁x̂(sₖh)
yₐ(sₖh) = η(sₖh)aₛ(sₖh) + (1 - η(sₖh))y(sₖh)
Donde:
• aₛ(sₖh): señal de ataque acotada en ʟ₂[0, ∞)
• η(sₖh) ∈ {0, 1}: variable de Bernoulli
• E{η(sₖh)} = η₀ = 0.5 (probabilidad de ataque)
Si η = 1: yₐ = aₛ (ataque activo - datos falsos inyectados)
Si η = 0: yₐ = y (transmisión normal)
El ataque se modela como una señal estocástica en lugar de determinística, capturando la naturaleza impredecible de los ciberataques reales.
9. Dinámicas en Lazo Cerrado con Retardos
Ecuaciones (24)-(25) - Sistema Completo con Retardos y Ataques
ẋ(t) = (A_c + B₁K₁)x(t) + Ew(t) - B₁K₁e(t)
- B₁K₁∫ₜ₋χ₍ₜ₎ẋ(s)ds + B₁K₁∫ₜ₋χ₍ₜ₎ė(s)ds
Error de observación:
ė(t) = (A_c - L₁C)e(t) + Ew(t) - L₁η(sₖh)aₛ(sₖh)
+ L₁Cη(sₖh)x(t) + L₁C(1-η(sₖh))∫ₜ₋α₍ₜ₎ẋ(s)ds
+ L₁(1-η(sₖh))ζₖ(t) - L₁C∫ₜ₋χ₍ₜ₎ẋ(s)ds
+ L₁C∫ₜχ₍ₜ₎ė(s)ds
• χ(t) = t - sₖh: retardo debido al event-triggering
• α(t): retardo de comunicación variable
• ζₖ(t): error de event-triggering
Límites:
0 ≤ χ(t) ≤ q₁ = l̄h + τ̄
τₛₖ ≤ α(t) ≤ h + τ̄
Estos retardos convierten al sistema en infinito-dimensional en el dominio de Laplace, requiriendo análisis en el dominio del tiempo mediante funcionales de Lyapunov.
10. Análisis de Estabilidad con Lyapunov-Krasovskii
Ecuación (30) - Funcional Lyapunov-Krasovskii
V₁(x(t)) = xᵀ(t)Z₁x(t) + ∫ₜ₋q₁ᵗₛᵗẋᵀ(v)Z₁ẋ(v)dvds
V₂(e(t)) = eᵀ(t)Z₁e(t) + ∫ₜq₁ᵗ∫ₛᵗė(v)Z₁ė(v)dvds
donde q₁ = l̄h + τ̄ = 4(0.001) + 0.01 = 0.014 s
• Términos cuadráticos: xᵀZ₁x y eᵀZ₁e miden la "energía" instantánea del sistema
• Términos integrales: Capturan la historia del sistema en el intervalo [t-q₁, t]
Relación con Laplace:
Aunque el análisis se realiza en el dominio del tiempo, las condiciones LMI resultantes garantizan que los polos del sistema permanezcan en el semiplano izquierdo, asegurando estabilidad asintótica equivalente a Re(s) < 0 en el dominio de Laplace.
Ecuación (29) - Condición LMI de Estabilidad
⎢ ⎥
⎢ * -Z₁/q₁ 0 0 ⎥
⎢ ⎥ < 0
⎢ * * -Z₁/q₁ 0 ⎥
⎢
* * * -Φ/σ⎦
Ganancias:
L₁ = Z₁⁻¹P̄₁
K₁ = (B₁ᵀB₁)⁻¹B₁ᵀZ₁⁻¹Q₁
11. Desempeño H∞ bajo Ataques y Perturbaciones
Ecuación (28) - Índice de Desempeño H∞
Donde:
‖y‖₂² = ∫₀^∞ yᵀ(t)y(t)dt
‖w‖₂² = ∫₀^ wᵀ(t)w(t)dt
‖aₛ‖₂² = ∫₀^∞ aₛᵀ(t)aₛ(t)dt
Interpretación en frecuencia:
La función de transferencia desde [w, aₛ] hasta y debe satisfacer:
‖G(s)‖∞ < γ para todo s = jω
El sistema atenúa tanto las perturbaciones de RES/carga (w) como los ataques FDI (aₛ) por un factor γ. Valores menores de γ indican mayor robustez.
Valores usados:
• γ = 50 para el control auxiliar
• γ = 200 para el control secundario (más conservador debido a retardos y ataques)
12. Resultados Numéricos y Validación
| Parámetro | Valor | Significado en el Plano s |
|---|---|---|
| h | 0.001 s | Período de muestreo (frecuencia de muestreo: 1000 Hz) |
| l̄ | 4 | Máximo número de eventos perdidos consecutivos |
| τ̄ | 0.01 s | Retardo máximo de comunicación |
| η₀ | 0.5 | Probabilidad de ataque FDI (50%) |
| σ | 0.01 | Parámetro event-triggered (menor σ → más transmisiones) |
| αₚₚ | 3 | Parte real mínima de polos: Re(s) < -3 |
| γₚₚ | 5 | Magnitud máxima de polos: |s| < 5 |
| δₚₚ | 0.8 rad | Ángulo de amortiguamiento: ζ ≥ cos(0.8) ≈ 0.697 |
Ganancias Obtenidas
K₂ = [-5.9923 -2.0885 -8.1073 3.8922 -6.0050]
Observador Auxiliar (L₂):
L₂ = [2249.705 -0.00408 -4.14745 0.000 0.99879
0.99879 -0.27056 -0.37878 0.000 0.99999]
Polos en Lazo Cerrado (Auxiliar):
λ₁,₂ = -3.3737 ± 3.1838i
λ₃,₄ = -3.5460 ± 1.4927i
λ₅ = -3.6912
✓ Parte real: -3.37 < -αₚₚ = -3
✓ Magnitud: √(3.37² + 3.18²) ≈ 4.64 < γₚₚ = 5
✓ Amortiguamiento: ζ = 3.37/4.64 ≈ 0.726 > cos(0.8) ≈ 0.697
Respuesta temporal (transformada inversa de Laplace):
Δf(t) = A₁e^(-3.37t)sin(3.18t + φ₁) + A₂e^(-3.55t)sin(1.49t + φ₂) + A₃e^(-3.69t)
Esto produce un decaimiento rápido con oscilaciones moderadamente amortiguadas.
| Esquema de Control | Desviación Máxima | ITAE | Transmisión |
|---|---|---|---|
| VI RPP SF (auxiliar) | ±0.3178 Hz | 6449.191 | - |
| TTC (σ=0) | ±0.22390 Hz | 3882.627 | 100% (900,000 muestras) |
| OETC (σ=0.01, sin ataque) | ±0.22667 Hz | 3776.747 | 2.85% (25,661 muestras) |
| OETC bajo FDI | ±2.39652 Hz ❌ | 38744.06 | 1.32% (11,948 muestras) |
| FDI AROETC (propuesto) | ±0.22636 Hz ✓ | 3769.064 | 3.2% (28,800 muestras) |
13. Conexión Fundamental con Transformadas de Laplace
🎯 Síntesis: El Papel Central de Laplace
Aunque el documento emplea métodos modernos en el dominio del tiempo (LMIs, Lyapunov-Krasovskii, event-triggering), toda la fundamentación teórica se basa en la transformada de Laplace:
1. Del Modelo Físico al Dominio s
Las ecuaciones diferenciales del microgrid (Ec. 1) se transforman mediante ℒ{·} a:
X(s) = (sI - A)⁻¹[x(0) + BU(s) + EW(s)]
La matriz (sI - A)⁻¹ contiene toda la información dinámica del sistema.
2. Polos Determinan Estabilidad y Respuesta
Los autovalores de A (raíces de det(sI - A) = 0) son los polos del sistema. El diseño RPP posiciona explícitamente estos polos en regiones específicas del plano s mediante las restricciones LMI (Ecs. 4-7), garantizando:
- Estabilidad: Re(s) < -αₚₚ < 0
- Rapidez: Controlada por αₚₚ (tasa de decaimiento)
- Amortiguamiento: ζ ≥ cos(δₚₚ)
- Frecuencia: ωₙ < γₚₚ
3. Función de Transferencia y Desempeño en Frecuencia
La norma H∞ evalúa la función de transferencia G(s) = C(sI - A - B₂K₂)⁻¹E sobre el
eje imaginario s = jω, garantizando que:
sup_ω σ̄(G(jω)) < γ
Esto asegura atenuación de perturbaciones en todas las frecuencias.
4. Retardos y la Trascendencia del Método
Los retardos introducen términos e^(-sτ) en Laplace, haciendo el sistema infinito-dimensional. Por esto se recurre al funcional Lyapunov-Krasovskii (Ec. 30), que trabaja en el dominio del tiempo pero preserva las garantías de estabilidad del plano s.
5. Validación Numérica
Los polos obtenidos (-3.37 ± 3.18i, etc.) definen completamente la respuesta temporal mediante la
transformada inversa:
Δf(t) = ℒ⁻¹{ΔF(s)} = Σ Aᵢe^(σᵢt)sin(ωᵢt + φᵢ)
Los resultados de simulación (desviación ±0.3178 Hz, ITAE = 6449.191) son la manifestación práctica
de esta teoría.
✅ Conclusión
La transformada de Laplace es el puente matemático que conecta:
1. El modelo físico del microgrid ↔ El análisis de estabilidad
2. Las especificaciones de diseño ↔ La ubicación de polos en el plano s
3. El desempeño en frecuencia ↔ La norma H∞
4. La respuesta temporal ↔ La transformada inversa ℒ⁻¹
Los métodos modernos (LMIs, Lyapunov, event-triggering) son extensiones sofisticadas
que preservan los principios fundamentales establecidos por Laplace, permitiendo manejar
complejidades adicionales (retardos, ataques, comunicación eficiente) sin perder la conexión
con la teoría clásica de control.
• Sección 2.1: Modelo en espacio de estado (Ec. 1)
• Sección 2.2: Diseño RPP-based state feedback (Ecs. 2-7)
• Sección 3: Esquema event-triggered (Ecs. 10-11)
• Sección 4: Control FDI AROETC (Ecs. 12, 22-25, 28-30)
• Sección 5: Resultados de simulación y validación
• Figuras 1-7: Esquemas del sistema y respuestas de frecuencia
• Tablas 1-4: Parámetros del sistema y resultados comparativos