Perturbation of monic matrix polynomials

1. Problema e Contesto

Il lavoro si concentra sull'analisi della stabilità spettrale dei polinomi matriciali monici sotto perturbazioni strutturate dei loro coefficienti.
Sia $P_u(\lambda) = A_d(u)\lambda^d + A_{d-1}(u)\lambda^{d-1} + \dots + A_0(u)$ un polinomio matriciale, dove $\lambda \in \mathbb{C}$ e i coefficienti $A_k(u) \in \mathbb{C}^{n \times n}$ dipendono da un vettore di parametri $u \in \mathbb{C}^p$.
L'ipotesi fondamentale è che il polinomio sia monic (cioè $A_d(u) \equiv I$, la matrice identità) e che i coefficienti siano funzioni semialgebriche continue rispetto al parametro $u$.

L'obiettivo principale è studiare la stabilità dei seguenti oggetti spettrali quando il parametro $u$ subisce piccole variazioni:

• Spettro $\sigma(u)$: l'insieme delle radici di $\det P_u(\lambda) = 0$.
• $\epsilon$-pseudospettro $\Lambda_\epsilon(u)$: l'insieme dei $\lambda$ che diventano autovalori per qualche perturbazione del parametro $v$ con $\|v-u\| \le \epsilon$.
• Campo di valori (Numerical Range) $W(u)$: l'insieme dei $\lambda$ tali che $x^* P_u(\lambda) x = 0$ per qualche vettore unitario $x$.
• Campo di valori congiunto (Joint Numerical Range) $JW(u)$: l'immagine dell'insieme delle forme quadratiche associate ai coefficienti del polinomio.
• Raggio spettrale $r(u)$.
• Coppie di Jordan: la struttura degli autovalori e delle catene di Jordan.

Il problema affronta come questi insiemi e strutture variano al variare di $u$, distinguendo tra perturbazioni arbitrarie (caso classico) e perturbazioni strutturate (indotte dalla variazione del parametro $u$).

2. Metodologia

L'approccio metodologico combina l'analisi spettrale classica con strumenti avanzati di geometria semialgebrica.

• Geometria Semialgebrica: Gli autori sfruttano il fatto che gli insiemi definiti da equazioni e disuguaglianze polinomiali (o semialgebriche) godono di proprietà topologiche e geometriche robuste. In particolare, utilizzano il Teorema di Tarski-Seidenberg per garantire che le proiezioni di insiemi semialgebrici (come gli spettri o i campi di valori definiti tramite proiezioni) rimangano insiemi semialgebrici.
• Disuguaglianza di Łojasiewicz: Per stabilire la continuità, si fa ricorso alla disuguaglianza di Łojasiewicz, che fornisce stime quantitative sulla velocità di variazione di funzioni semialgebriche continue. Questo permette di dimostrare la continuità di Hölder.
• Decomposizione in Varietà: Si utilizza il Teorema di Decomposizione delle Celle per suddividere lo spazio dei parametri $\mathbb{C}^p$ in un'unione finita di sottovarietà analitiche semialgebriche. Su ciascuna di queste varietà, le funzioni spettrali diventano analitiche.
• Linearizzazione: Per lo studio delle coppie di Jordan, il problema viene ridotto allo studio della matrice compagna (companion matrix) associata al polinomio, trasformando il problema non lineare in uno lineare.

3. Risultati Chiave

A. Continuità di Hölder delle Mappe a Insiemi

Il risultato principale è la dimostrazione che le mappe a valori insiemistici (spettro, pseudospettro, campo di valori) sono continue di Hölder rispetto al parametro $u$.

• Teorema 3.5: La mappa dello spettro $\sigma: \mathbb{C}^p \Rightarrow \mathbb{C}$ è continua di Hölder di grado $\alpha$ in un intorno di ogni punto.
• Teorema 3.7: La mappa dell'$\epsilon$-pseudospettro $\Lambda_\epsilon$ è continua di Hölder sia rispetto a $u$ che rispetto al parametro $\epsilon$.
• Teorema 3.12: La mappa del campo di valori $W$ è continua di Hölder.
• Teorema 3.15: La mappa del campo di valori congiunto $JW$ è continua di Hölder.

B. Stima Quantitativa e Genericità

Il Teorema 3.17 rafforza i risultati precedenti mostrando l'esistenza di un sottoinsieme semialgebrico aperto e denso $U \subset \mathbb{C}^p$ su cui:

• Le mappe $\sigma, \Lambda_\epsilon, W, \partial W$ sono continue di Hölder con gradi espliciti (es. $1/(nd)$ per lo spettro, $1/d$ per il campo di valori).
• La mappa $JW$ è Lipschitziana (grado 1).
• Il raggio spettrale $r(u)$ è una funzione analitica.

C. Proprietà Geometriche degli Insiemi

• Bordi Algebrici: È dimostrato che i bordi dell'$\epsilon$-pseudospettro e del campo di valori giacciono su curve algebriche. Di conseguenza, questi bordi hanno un numero finito di singolarità (punti di cuspide o auto-intersezioni) e punti non differenziabili.
• Componenti Connesse: L'$\epsilon$-pseudospettro ha al più $dn$ componenti connesse, e ogni componente contiene almeno un autovalore per ogni perturbazione ammissibile.

D. Stabilità delle Coppie di Jordan

Il Teorema 4.3 stabilisce che lo spazio dei parametri $\mathbb{C}^p$ può essere partizionato in un'unione finita di sottovarietà analitiche semialgebriche $M_1, \dots, M_N$ tali che, su ciascuna $M_t$:

1. Il numero e le dimensioni dei blocchi di Jordan sono costanti.
2. Gli autovalori $\lambda_k(u)$ sono funzioni analitiche di $u$.
3. Le matrici che formano le coppie di Jordan ($X_u, J_u$) dipendono analiticamente da $u$.

4. Significato e Contributi

Questo lavoro offre contributi significativi in diversi ambiti:

1. Perturbazioni Strutturate: A differenza della teoria classica che considera perturbazioni arbitrarie di norma limitata, questo studio si focalizza su perturbazioni strutturate indotte da parametri fisici o di progetto. Questo è cruciale per applicazioni reali (controllo, ingegneria) dove i coefficienti non variano indipendentemente ma seguono vincoli fisici.
2. Robustezza Quantitativa: Fornisce stime quantitative (esponenti di Hölder) sulla stabilità, permettendo di quantificare quanto un errore nel parametro $u$ si traduca in un errore nello spettro o nel campo di valori.
3. Struttura Geometrica: L'uso della geometria semialgebrica rivela che, nonostante la complessità dei polinomi matriciali, gli oggetti spettrali hanno una struttura geometrica "ben comportata" (bordi algebrici, numero finito di singolarità) e che la loro regolarità è garantita quasi ovunque (su un insieme denso).
4. Generalizzazione: Sebbene il paper assuma la monicità per semplificare l'esposizione, gli autori notano che la maggior parte dei risultati vale anche se il coefficiente principale $A_d(u)$ è solo localmente non singolare, rendendo i risultati ampiamente applicabili.

In sintesi, il paper fornisce un quadro teorico rigoroso per la stabilità dei polinomi matriciali in contesti parametrici, collegando l'analisi numerica alla geometria algebrica reale e complessa.