Una Representación Exacta y de Dimensión Finita para Multiplicadores de Criterio de Círculo de Bloque Completo

Información Básica

• ID del Artículo: 2511.20995
• Título: Una Representación Exacta y de Dimensión Finita para Multiplicadores de Criterio de Círculo de Bloque Completo
• Autores: Felix Biertümpfel (Universidad de Michigan & TU Dresden), Bin Hu (Universidad de Illinois en Urbana-Champaign), Geir Dullerud (Universidad de Minnesota), Peter Seiler (Universidad de Michigan)
• Clasificación: eess.SY (Sistemas y Control), cs.SY
• Fecha de Publicación: 26 de noviembre de 2025 (Presentación en arXiv)
• Enlace del Artículo: https://arxiv.org/abs/2511.20995

Resumen

Este artículo proporciona por primera vez una caracterización de dimensión finita para el conjunto completo de multiplicadores de criterio de círculo de bloque completo (multiplicadores de bloque completo) para no linealidades acotadas en sector no repetidas. El estudio considera sistemas lineales invariantes en el tiempo de tiempo discreto interconectados con retroalimentación no lineal acotada en sector no repetida. Las condiciones suficientes para estabilidad y desempeño pueden derivarse mediante: (i) desigualdades de disipatividad y (ii) restricciones cuadráticas (QC) sobre pares de entrada/salida no lineales acotados. Las clases más grandes de QC (o multiplicadores) pueden reducir el conservadurismo de las condiciones. Los multiplicadores de criterio de círculo de bloque completo definen el conjunto completo de todas las QC posibles para no linealidades acotadas en sector no repetidas, proporcionando la condición menos conservadora. Sin embargo, los multiplicadores de bloque completo se definen mediante un número infinito no numerable de restricciones, por lo que la forma original no es computacionalmente viable. La perspectiva teórica clave de este artículo es: el conjunto de todos los pares de entrada/salida para no linealidades acotadas en sector no repetidas es igual al conjunto de todos los pares incrementales de una función lineal por tramos apropiadamente construida. La nueva descripción requiere solo un número finito de restricciones de copositividad matricial, con una implementación computacional exacta viable para problemas con dimensión de entrada/salida no lineal ≤4.

Antecedentes y Motivación de la Investigación

Definición del Problema

Este artículo estudia sistemas formados por la interconexión de un sistema lineal invariante en el tiempo (LTI) de tiempo discreto conocido con retroalimentación no lineal estática sin memoria. El objetivo del análisis es evaluar la estabilidad y la ganancia inducida del sistema interconectado.

Importancia del Problema

1. Problema Fundamental en Control Robusto: Las no linealidades acotadas en sector están ampliamente presentes en sistemas de control prácticos (como saturación, zona muerta, funciones de activación ReLU, etc.), y su análisis de estabilidad es un problema central en teoría de control
2. Verificación de Redes Neuronales: Para redes neuronales recurrentes pequeñas (como aquellas con funciones de activación ReLU), el método de este artículo puede usarse para análisis de estabilidad incremental
3. Completitud Teórica: Proporcionar herramientas de análisis no conservadoras es crucial para comprender los límites de desempeño del sistema

Limitaciones de Métodos Existentes

1. Multiplicadores Diagonales (M_d): Computacionalmente simples pero excesivamente conservadores, incapaces de aprovechar adecuadamente las características estructurales de la no linealidad
2. Relajaciones Convexas (M_c): Mejoran el conservadurismo mediante la imposición de restricciones en los vértices del hipercubo del sector, pero siguen siendo incompletas
3. Multiplicadores de Bloque Completo (M_fb): Teóricamente proporcionan el conjunto completo de QC, pero se definen mediante un número infinito no numerable de restricciones (se requiere una restricción para cada punto en α,β^m), siendo completamente inviable computacionalmente
4. Métodos de Aproximación Existentes: Basados en generalizaciones del teorema de Polya o aproximaciones de cuadrícula finita, pero sin garantía de exactitud

Motivación de la Investigación

Existe una brecha enorme entre la teoría y la práctica computacional: el multiplicador de bloque completo, teóricamente óptimo, es computacionalmente inviable. Este artículo tiene como objetivo encontrar una caracterización equivalente pero de dimensión finita, haciendo que la solución teóricamente óptima sea computable en la práctica.

Contribuciones Principales

1. Avance Teórico: Se demuestra que el conjunto de todos los pares de entrada/salida para no linealidades acotadas en sector no repetidas es igual al conjunto de pares incrementales de una función lineal por tramos específica (Lema 5)
2. Caracterización de Dimensión Finita: Se proporciona una nueva caracterización de dimensión finita para multiplicadores de criterio de círculo de bloque completo (Corolario 1), transformando el problema de restricciones infinitas en 4^m restricciones de copositividad
3. Método de Cálculo Exacto: Para problemas de dimensión m≤4, utilizando relajaciones exactas de copositividad conocidas, se logra una implementación computacional exacta del conjunto completo de QC
4. Verificación de Practicidad: Mediante ejemplos numéricos se demuestra la ventaja del nuevo método sobre métodos existentes, verificando márgenes de estabilidad más grandes e intervalos de desempeño más ajustados

Explicación Detallada del Método

Definición de la Tarea

Configuración del Sistema: Considere la interconexión de retroalimentación Fu(G,ΔΦ) mostrada en la Figura 1:

• Sistema Nominal G: Sistema LTI de tiempo discreto

  x(k+1) = Ax(k) + B₁w(k) + B₂u(k)
  v(k) = C₁x(k) + D₁₁w(k) + D₁₂u(k)
  y(k) = C₂x(k) + D₂₁w(k) + D₂₂u(k)
  

• Perturbación ΔΦ: Definida por no linealidades escalares no repetidas, w(k) = Φ(v(k)), donde Φ(v) = φ₁(v₁),...,φₘ(vₘ)ᵀ

Definición de Acotamiento en Sector: Una función escalar φ está acotada en el sector α,β si satisface:

(φ(x) - αx)(βx - φ(x)) ≥ 0, ∀x ∈ ℝ

Objetivos del Análisis:

1. Bien-planteamiento (well-posedness)
2. Estabilidad interna
3. Ganancia inducida ℓ₂ finita: ‖Fu(G,ΔΦ)‖₂→₂ < γ

Marco Teórico Principal

1. Fundamentos de Restricciones Cuadráticas (QC)

Definición: Un conjunto G ⊂ ℝ^(2m) satisface una QC definida por M ∈ S^(2m) si:

zᵀMz ≥ 0, ∀z ∈ G

Aplicación a No linealidades: Para la gráfica de Φ, G(Φ) = {v; w : v ∈ ℝᵐ, w = Φ(v)}, si:

[v; w]ᵀM[v; w] ≥ 0, ∀v ∈ ℝᵐ, w = Φ(v)

entonces M define una QC válida para Φ.

Conjunto Completo:

G(sec[α,β]ᵐ) := ⋃_{Φ∈sec[α,β]ᵐ} G(Φ)

2. Condición de Estabilidad (Teorema 1)

Dadas las suposiciones:

1. D₁₁ = 0 (garantiza bien-planteamiento)
2. G(secα,βᵐ) ⊂ QC(M) (M es un multiplicador válido)
3. ∃P ≽ 0, γ > 0 tal que L(P,M,γ²) ≺ 0

donde L es una función afín que contiene matrices del sistema, matriz de Lyapunov P y multiplicador M.

Conclusión: El sistema Fu(G,ΔΦ) está bien-planteado, es internamente estable y ‖Fu(G,ΔΦ)‖₂→₂ < γ.

Innovación Teórica Clave

Construcción de Función Lineal por Tramos

Defina la función lineal por tramos especial:

f_{αβ}(x) = {αx,  x ≤ 0
            {βx,  x > 0

Función correspondiente repetida: Fαβ(v) = f_{αβ}(v₁),...,f_{αβ}(vₘ)ᵀ

Gráfica Incremental:

I(Fαβ) := {[v̄-v̂; w̄-ŵ] : ∃v̄,v̂ ∈ ℝᵐ s.t. w̄=Fαβ(v̄), ŵ=Fαβ(v̂)}

Equivalencia Principal (Lema 5)

Conclusión Clave: I(Fαβ) = G(secα,βᵐ)

Esquema de Prueba:

1. I(Fαβ) ⊆ G(secα,βᵐ): Para cualquier par incremental (dv,dw), puede construirse una matriz diagonal Γ tal que dw_i = γ_i·dv_i, donde γ_i ∈ α,β
2. G(secα,βᵐ) ⊆ I(Fαβ): Para cualquier par entrada/salida (v,w), defina:

   v̄ᵢ = (γᵢ-α)/(β-α)·vᵢ, v̂ᵢ = (γᵢ-β)/(β-α)·vᵢ
   

   donde γᵢ = wᵢ/vᵢ (o (α+β)/2 si vᵢ=0), puede verificarse que v̄-v̂=v y w̄-ŵ=w

Intuición Geométrica (Figura 2): Los incrementos entre dos puntos cualesquiera de la función lineal por tramos pueden realizar cualquier pendiente dentro de α,β, cubriendo así todos los pares entrada/salida de no linealidades acotadas en sector.

Nueva Caracterización de Dimensión Finita

Conjunto de QC Incrementales (Teorema 2)

Defina la función:

g_M(Γ̄,Γ̂) = [Γ̄  -Γ̂      ]ᵀ   [Γ̄  -Γ̂      ]
             [rΓ̄+cI -rΓ̂-cI]  M [rΓ̄+cI -rΓ̂-cI]

donde c = (α+β)/2, r = (β-α)/2

Nuevo Conjunto de Multiplicadores:

M_inc := {M ∈ S^(2m) : g_M(Γ̄,Γ̂) ∈ COP^(2m), 
          ∀Γ̄,Γ̂ ∈ diag({-1,1}ᵐ)}

Teorema: I(Fαβ) ⊂ QC(M) ⟺ M ∈ M_inc

Puntos Clave de la Prueba:

• Utiliza la descomposición de |v̄|,|v̂| ∈ ℝ^m_≥0 y matrices de signos Γ̄,Γ̂
• Transforma la condición de QC incremental en condiciones de copositividad
• La copositividad garantiza no negatividad para todos los vectores no negativos

Resultado de Equivalencia (Corolario 1)

Conclusión Principal: M_inc = M_fb

Cadena de Prueba:

1. Por Teorema 2: M ∈ M_inc ⟺ I(Fαβ) ⊂ QC(M)
2. Por Lema 5: I(Fαβ) = G(secα,βᵐ)
3. Por Lema 2: G(secα,βᵐ) ⊂ QC(M) ⟺ M ∈ M_fb

Características Computacionales

Número de Restricciones: M_inc requiere 4^m = 2^(2m) restricciones de copositividad (una para cada par Γ̄,Γ̂ ∈ diag({-1,1}ᵐ))

Complejidad Computacional:

• La verificación de copositividad es un problema co-NP-completo
• Relajación Exacta: Para m≤4, cualquier matriz copositiviva es ordinariamente copositiviva (suma de matriz semidefinida positiva y matriz no negativa), la relajación es exacta
• Para m≥5, existen matrices copositivivas no ordinarias (como matrices de Horn), requiriendo aproximaciones conservadoras

Relación de Inclusión de Conjuntos:

M_d ⊂ M_c ⊂ M_fb = M_inc

Configuración Experimental

Configuración del Sistema

• Dimensión del Sistema: Sistema LTI de tiempo discreto de orden 3
• Dimensión de No linealidad: m = 3 (tres no linealidades acotadas en sector independientes)
• Límites de Sector: 0,β, con β variando en el rango 0,1.3
• Contexto de Aplicación: Similar al análisis de estabilidad incremental de funciones de activación ReLU en pequeñas redes neuronales recurrentes

Matrices de Espacio de Estados

Las matrices específicas del sistema G (A,B₁,B₂,C₁,C₂,D₁₁,D₁₂,D₂₁,D₂₂) se proporcionan en el artículo, nótese que D₁₁≠0 (sistema asumido bien-planteado).

Métodos de Comparación

1. M_d: Multiplicador diagonal (más conservador pero computacionalmente más rápido)
2. M_c: Relajación convexa (conservadurismo medio)
3. M_inc: Conjunto completo propuesto (menos conservador)

Detalles de Implementación

• Solucionador: MATLAB + CVX + MOSEK
• Hardware: Procesador de 4.51GHz, 24GB de memoria
• Implementación de Copositividad: Para m=3 se utiliza relajación exacta (descomposición semidefinida positiva + no negativa)

Métricas de Evaluación

• Límite Superior de Ganancia: γ tal que ‖Fu(G,ΔΦ)‖₂→₂ < γ
• Margen de Estabilidad: Máximo límite de sector β_max verificable
• Tiempo de Cálculo: Tiempo requerido para resolver SDP

Resultados Experimentales

Resultados Principales (caso β=1)

Hallazgos Clave:

1. Conservadurismo: M_inc reduce el conservadurismo en 12.87% comparado con M_c, y 47.3% comparado con M_d
2. Costo Computacional: El tiempo de cálculo de M_inc es aproximadamente 2 veces el de M_c, pero sigue siendo aceptable (<0.5s)
3. Exactitud: M_inc proporciona el límite teóricamente óptimo (para QC estáticas)

Análisis de Variación del Límite de Sector (Figura 3)

Prueba de 15 puntos equidistantes en β ∈ 0,1.3:

Comportamiento de Convergencia (β→0):

• Todos los métodos convergen a la ganancia nominal γ_nom = 1.396
• Verifica la corrección de los métodos

Evolución de la Brecha de Conservadurismo:

• β ∈ 0,0.6: La brecha entre los tres métodos es relativamente pequeña
• β > 0.6: La brecha entre M_d y multiplicadores de bloque completo aumenta significativamente
• β > 0.9: La brecha entre M_c y M_inc se abre claramente

Margen de Estabilidad:

Significancia Estadística: M_inc proporciona el límite más ajustado en todos los puntos de prueba, con la diferencia siendo especialmente significativa en límites de sector grandes.

Hallazgos Experimentales

1. Sensibilidad a Dimensión: Para problemas de m=3, la relajación exacta de copositividad funciona bien
2. Viabilidad Computacional: Incluso para el conjunto completo M_inc, el tiempo de cálculo sigue siendo práctico
3. Mejora de Conservadurismo: La mejora es más significativa en regiones críticas (cerca del límite de estabilidad)
4. Verificación Teórica: Los resultados experimentales son consistentes con predicciones teóricas (M_inc = M_fb proporciona el límite óptimo)

Trabajo Relacionado

Teoría de Restricciones Cuadráticas y Multiplicadores

1. Marco IQC: Megretski & Rantzer (2002) establecen la teoría fundamental de restricciones cuadráticas integrales
2. Multiplicadores de Bloque Completo: Veenman & Scherer (2014-2017) estudian sistemáticamente las propiedades y aplicaciones de multiplicadores de bloque completo
3. Sistemas de Tiempo Discreto: Fetzer & Scherer (2017) analizan específicamente la estabilidad absoluta de interconexiones de retroalimentación de tiempo discreto

Métodos Computacionales

1. Generalización del Teorema de Polya: Fetzer & Scherer (2017) utilizan el teorema de Polya para proporcionar aproximaciones de M_fb
2. Cuadrícula Finita: Aproximación de M_fb mediante imposición de restricciones en puntos de cuadrícula finita en α,β^m
3. Relajación Convexa: M_c impone restricciones en vértices del hipercubo (línea base de comparación de este artículo)

Ventajas de Este Artículo

1. Completitud Teórica: Primera caracterización exacta de dimensión finita de M_fb
2. Viabilidad Computacional: Implementación exacta computable para m≤4 (métodos anteriores son inexactos o no computables)
3. Marco Unificado: Establece conexión profunda entre no linealidades de sector y funciones lineales por tramos mediante análisis incremental

Direcciones de Investigación Relacionadas pero Distintas

1. No linealidades Repetidas: Lema 4 demuestra que no linealidades repetidas y no repetidas tienen el mismo conjunto completo de QC
2. Análisis de Redes ReLU: Noori et al. (2024) estudian el conjunto completo de QC para ReLU repetidas
3. Análisis Regional: Valmorbida et al. (2018) estudian análisis regional de sistemas Lurie con pendiente limitada

Conclusiones y Discusión

Conclusiones Principales

1. Contribución Teórica: Se demuestra que el conjunto de pares entrada/salida para no linealidades acotadas en sector no repetidas es igual al conjunto de pares incrementales de una función lineal por tramos específica
2. Contribución Computacional: Transformación de la caracterización infinita-dimensional de multiplicadores de bloque completo en una caracterización finita-dimensional con 4^m restricciones de copositividad
3. Contribución Práctica: Proporciona implementación exacta computable para problemas con m≤4
4. Conclusión de Verificación: Experimentos numéricos confirman ventajas significativas del nuevo método en conservadurismo y margen de estabilidad

Limitaciones

1. Restricción de Dimensión:
   • La relajación exacta es válida solo para m≤4
   • Para m≥5 se requieren aproximaciones conservadoras de copositividad (como descomposición semidefinida positiva + no negativa)
   • El número de restricciones crece exponencialmente (4^m)
2. Complejidad Computacional:
   • La verificación de copositividad es esencialmente co-NP-completa
   • Para problemas grandes (m>4), la carga computacional puede ser excesiva
3. Condiciones de Suposición:
   • Requiere D₁₁=0 para garantizar bien-planteamiento (puede relajarse bajo suposiciones adicionales)
   • Solo aplicable a no linealidades estáticas sin memoria
   • Limitado a no linealidades del tipo acotadas en sector
4. Fuentes de Conservadurismo:
   • Utiliza QC estáticas (multiplicadores dinámicos podrían ser menos conservadores)
   • Función de Lyapunov de tipo cuadrático (funciones de Lyapunov no cuadráticas podrían ser mejores)

Direcciones Futuras

1. Extensión a Dimensiones Altas:
   • Investigar aproximaciones de copositividad más ajustadas para m>4
   • Explorar condiciones de copositividad estructuradas para reducir carga computacional
   • Desarrollar algoritmos numéricos escalables
2. Multiplicadores Dinámicos:
   • Extender caracterización de dimensión finita a multiplicadores dinámicos
   • Investigar utilización de características en dominio de frecuencia
3. Otros Tipos de No linealidades:
   • Extender a pendiente limitada, monótona, simétricamente impar y otros tipos
   • Investigar no linealidades con restricciones mixtas
4. Expansión de Aplicaciones:
   • Verificación y certificación de redes neuronales
   • Análisis de sistemas de control adaptativo
   • Sistemas de control en red
5. Profundización Teórica:
   • Investigar diferencias esenciales entre no linealidades repetidas y no repetidas
   • Explorar aplicabilidad del método de gráfica incremental en otros problemas

Evaluación Profunda

Fortalezas

1. Innovación Teórica (★★★★★):
   • Primera solución al problema de larga data de caracterización de dimensión finita de multiplicadores de bloque completo
   • La equivalencia de gráfica incremental (Lema 5) proporciona perspectiva teórica profunda
   • Prueba rigurosa y completa, lógica clara
2. Valor Práctico (★★★★☆):
   • Proporciona método exacto computable para m≤4, cubriendo muchas aplicaciones prácticas
   • Tiempo de cálculo en rango aceptable (<0.5s)
   • Buena estabilidad numérica (basada en solucionadores SDP maduros)
3. Elegancia del Método (★★★★★):
   • Transformación ingeniosa del problema mediante construcción de función lineal por tramos especial
   • Simplificación de prueba mediante descomposición de matriz de signos y valores absolutos
   • Unificación de perspectivas de no linealidades de sector y análisis incremental
4. Suficiencia Experimental (★★★☆☆):
   • Demuestra claramente ventajas sobre métodos existentes
   • Cubre rango de parámetros críticos
   • Pero solo un ejemplo numérico, falta diversidad
5. Calidad de Escritura (★★★★★):
   • Estructura clara, lógica rigurosa
   • Definición de símbolos normalizada y uniforme
   • Figura 2 proporciona excelente intuición geométrica

Insuficiencias

1. Escalabilidad de Dimensión:
   • Imposibilidad de implementación exacta para m>4 es limitación seria
   • Aplicación limitada para redes neuronales modernas (típicamente m≫4)
   • Número de restricciones exponencial (4^m) es cuello de botella fundamental
2. Limitaciones Experimentales:
   • Solo un ejemplo de 3 dimensiones, falta experimentos de 4 dimensiones o superiores
   • No se prueba conservadurismo de relajación de copositividad en m=5
   • Falta comparación detallada con métodos basados en teorema de Polya
3. Detalles Computacionales Insuficientes:
   • No se discute estabilidad numérica de relajación de copositividad
   • No se analiza convergencia y precisión de solucionador SDP
   • Falta análisis de complejidad computacional para problemas grandes
4. Rango de Aplicación:
   • Limitado a sistemas de tiempo discreto (tiempo continuo requiere tratamiento separado)
   • Suposición de no linealidad estática excluye incertidumbre dinámica
   • Suposición de acotamiento en sector puede no satisfacerse en algunas aplicaciones
5. Profundidad Teórica:
   • No explora por qué m≤4 es dimensión crítica para relajación exacta
   • No discute aplicabilidad más amplia del método de gráfica incremental
   • Relación con multiplicadores dinámicos no analizada profundamente

Evaluación de Impacto

Impacto a Corto Plazo (1-2 años):

• Proporciona herramienta de análisis óptima para sistemas de control pequeños
• Impulsa adopción de métodos más precisos en verificación de redes neuronales
• Estimula aplicación adicional de teoría de copositividad en control

Impacto a Mediano Plazo (3-5 años):

• Posible contenido estándar en libros de texto de control robusto
• Inspira nuevos métodos de aproximación para problemas de alta dimensión
• Promueve aplicación del paradigma de análisis incremental en otros problemas

Impacto a Largo Plazo (5+ años):

• Si extensión a alta dimensión tiene éxito, podría revolucionar análisis de sistemas no lineales
• Marco teórico podría extenderse a tipos de restricciones más amplios
• Se convierte en trabajo clásico en intersección de teoría de control y optimización

Escenarios de Cita Potencial:

1. Investigación de conjuntos de multiplicadores completos para otros tipos de no linealidades
2. Desarrollo de nuevos algoritmos de aproximación de copositividad
3. Certificación de robustez de redes neuronales
4. Análisis de estabilidad de control adaptativo y de aprendizaje

Escenarios de Aplicabilidad

Aplicaciones Más Adecuadas:

1. Sistemas de Control de Baja Dimensión (m≤4):
   • Control de actitud en aeroespacial (típicamente 3D)
   • Control de articulaciones robóticas (articulación individual o robot pequeño)
   • Convertidores de electrónica de potencia (pocos elementos de conmutación)
2. Verificación de Redes Neuronales Pequeñas:
   • RNN pequeñas en controladores embebidos
   • Redes neuronales simples en sistemas críticos de seguridad
   • Certificación de estabilidad incremental y robustez
3. Investigación Teórica:
   • Prueba de referencia para probar conservadurismo de nuevos algoritmos
   • Investigación de problemas fundamentales en teoría de multiplicadores
   • Enseñanza y demostración de método de análisis óptimo

Escenarios Menos Adecuados:

1. Sistemas de Alta Dimensión (m>4): Requiere métodos de aproximación
2. Redes Neuronales Grandes: Computacionalmente inviable
3. Aplicaciones en Tiempo Real: Tiempo de cálculo puede ser excesivo
4. Incertidumbre Dinámica: Requiere extensión teórica

Evaluación de Reproducibilidad

Disponibilidad de Código: El artículo no proporciona código, pero la implementación es relativamente directa

Dificultad de Reproducción (★★☆☆☆):

• Marco SDP estándar y maduro (CVX+MOSEK)
• Teoría de relajación de copositividad disponible (libro de Berman & Shaked-Monderer)
• Desafío principal es implementación correcta de generación de 4^m restricciones

Recomendaciones:

• Los autores deberían proporcionar implementación de código abierto
• Se necesita descripción detallada del algoritmo de descomposición de copositividad
• Deberían proporcionarse más casos de prueba

Referencias Bibliográficas (Literatura Clave Seleccionada)

1. Megretski & Rantzer (2002): "System analysis via integral quadratic constraints" - Fundamentos de teoría IQC
2. Veenman et al. (2016): "Robust stability and performance analysis based on integral quadratic constraints" - Revisión exhaustiva
3. Fetzer & Scherer (2017): "Full-block multipliers for repeated, slope-restricted scalar nonlinearities" - Teoría de multiplicadores de bloque completo
4. Berman & Shaked-Monderer (2003): "Completely Positive Matrices" - Obra clásica en teoría de copositividad
5. Boyd et al. (1994): "Linear Matrix Inequalities in System and Control Theory" - Fundamentos de métodos LMI

Puntuación General: ★★★★☆ (4.5/5)

Este es un artículo de alta calidad en teoría de control que resuelve un problema abierto importante. La innovación teórica es significativa, el método es elegante y la prueba es rigurosa. Las limitaciones principales son escalabilidad de dimensión y amplitud de experimentos. Para problemas pequeños, este es actualmente el mejor método; para problemas grandes, proporciona un importante punto de referencia teórico y nuevas direcciones de investigación. Se recomienda publicación en revista de control de primer nivel (como Automatica o IEEE TAC).