2025-11-19T00:28:13.831044

Stochastic Stability of Monotone Dynamical Systems. I. The Irreducible Cooperative Systems

Jiang, Sheng, Wang
The current series of papers is concerned with stochastic stability of monotone dynamical systems by identifying the basic dynamical units that can survive in the presence of noise interference. In the first of the series, for the cooperative and irreducible systems, we will establish the stochastic stability of a dynamical order, that is, the zero-noise limit of stochastic perturbations will be concentrated on a simply ordered set consisting of Lyapunov stable equilibria. In particular, we utilize the Freidlin--Wentzell large deviation theory to gauge the rare probability in the vicinity of unordered chain-transitive invariant set on a nonmonotone manifold. We further apply our theoretic results to the stochastic stability of classical positive feedback systems by showing that the zero-noise limit is a convex combination of the Dirac measures on a finite number of asymptotically stable equilibria although such system may possess nontrivial periodic orbits.
academic

Stabilità Stocastica dei Sistemi Dinamici Monotoni. I. I Sistemi Cooperativi Irriducibili

Informazioni Fondamentali

  • ID Articolo: 2412.19977
  • Titolo: Stochastic Stability of Monotone Dynamical Systems. I. The Irreducible Cooperative Systems
  • Autori: Jifa Jiang (Università Normale di Henan), Xi Sheng, Yi Wang (Università di Scienza e Tecnologia della Cina)
  • Classificazione: math.DS (Sistemi Dinamici)
  • Data di Pubblicazione: Dicembre 2024 (arXiv v2: 15 novembre 2025)
  • Link Articolo: https://arxiv.org/abs/2412.19977v2

Riassunto

Questa serie di articoli studia la stabilità stocastica dei sistemi dinamici monotoni, con l'obiettivo di identificare le unità dinamiche fondamentali che possono "sopravvivere" sotto perturbazioni rumorose. Nel primo articolo della serie, gli autori stabiliscono la stabilità stocastica dell'ordine dinamico per sistemi cooperativi irriducibili, dimostrando che il limite di rumore zero delle perturbazioni stocastiche si concentra su un insieme ordinato semplice composto da punti di equilibrio Lyapunov stabili. Gli autori utilizzano in particolare la teoria delle grandi deviazioni di Freidlin-Wentzell per stimare le probabilità rare vicino agli insiemi invarianti di catena transitiva disordinati su varietà non monotone. I risultati teorici vengono inoltre applicati alla stabilità stocastica dei sistemi classici di retroazione positiva, dimostrando che nonostante il sistema possa possedere orbite periodiche non banali, il limite di rumore zero rimane una combinazione convessa di misure di Dirac su un numero finito di punti di equilibrio asintoticamente stabili.

Contesto di Ricerca e Motivazione

Problemi di Ricerca

I sistemi dinamici monotoni derivano dal lavoro pioneristico di M. W. Hirsch, caratterizzati dall'esistenza di principi di confronto coerenti con relazioni di ordine parziale chiuse nello spazio degli stati (indotte da coni convessi). Sebbene le traiettorie in avanti dei sistemi monotoni deterministici convergano quasi ovunque (sia nel senso topologico che in quello della teoria della misura), il comportamento di questi sistemi sotto perturbazioni rumorose rimane poco chiaro.

Importanza

  1. Esigenze di Applicazione Pratica: I modelli matematici derivanti da equazioni differenziali e alle differenze sono frequentemente soggetti a perturbazioni rumorose. Lo studio di come il rumore influenzi la dinamica dei sistemi monotoni è cruciale sia per la modellazione teorica che per l'indagine del comportamento dei sistemi reali.
  2. Sfide Teoriche: La ricerca esistente adotta principalmente approcci basati su traiettorie (trajectory-based approach) nel quadro dei sistemi monotoni stocastici. Tuttavia, questo metodo può portare a un'eccessiva semplificazione della dinamica (come il fenomeno di sincronizzazione del rumore), senza rivelare pienamente gli effetti delle perturbazioni rumorose.
  3. Innovazione Metodologica: L'approccio basato sulla distribuzione (distribution-based approach) è considerato una strategia più efficace per studiare perturbazioni rumorose generali, con il problema centrale di identificare le unità dinamiche statistiche fondamentali che possono "sopravvivere" sotto interferenze rumorose.

Limitazioni dei Metodi Esistenti

  1. Semplificazione dell'Approccio Orbitale: Il comportamento a lungo termine dei sistemi monotoni stocastici è molto più semplice rispetto ai corrispondenti deterministici (come la sincronizzazione del rumore), senza fornire intuizioni approfondite.
  2. Dipendenza da Obiettivi Predefiniti: La ricerca precedente sulla stabilità stocastica si è principalmente concentrata su obiettivi dinamici predefiniti (come attrattori/repellenti, classi di equivalenza quasipotenziale), mancando di analisi della struttura su larga scala del sistema.
  3. Strutture Atipiche: I sistemi monotoni non sono né sistemi gradiente né sistemi assiomaticamente A, e il comportamento a lungo termine potrebbe non essere rappresentabile come unione di un numero finito di classi di equivalenza (a causa della possibile esistenza di infinite classi di equivalenza su varietà non monotone).

Motivazione della Ricerca

Questo articolo adotta un'idea strutturale per analizzare la stabilità stocastica dei sistemi dinamici monotoni, identificando dove si concentra il limite di rumore zero e stabilendo la stabilità stocastica dell'ordine dinamico per sistemi cooperativi irriducibili. Questo è il primo studio sulla stabilità stocastica di sistemi che potrebbero possedere potenzialmente infinite classi di equivalenza.

Contributi Fondamentali

  1. Stabilità Stocastica dell'Ordine Dinamico: Si dimostra che il limite di rumore zero dei sistemi cooperativi irriducibili si concentra su un insieme ordinato semplice (stationary p-arc) composto da punti di equilibrio Lyapunov stabili (Teorema 3.1).
  2. Proprietà Critiche del Quasipotenziale: Si dimostra che per qualsiasi insieme di catena transitiva disordinato K su una varietà non monotone, il quasipotenziale da K al suo attrattore duale superiore (inferiore) è sempre zero (Proposizioni 4.2 e 4.3), rappresentando una svolta tecnica cruciale.
  3. Caso Analitico/Finito di Punti di Equilibrio: Quando il sistema è analitico o possiede solo un numero finito di punti di equilibrio, il limite di rumore zero è una combinazione convessa di misure di Dirac su punti di equilibrio asintoticamente stabili (Corollario 3.2).
  4. Applicazione ai Sistemi di Retroazione Positiva: Si dimostra che il limite di rumore zero dei sistemi classici di retroazione positiva (inclusi i circuiti di controllo biochimico di tipo Griffith) si concentra su un numero finito di punti di equilibrio asintoticamente stabili, anche se il sistema potrebbe possedere orbite periodiche non banali (Proposizione 6.1).
  5. Innovazione Metodologica: Fusione delle idee di Benäım-Hirsch sulla teoria degli algoritmi di approssimazione stocastica e della teoria delle grandi deviazioni di Freidlin-Wentzell, utilizzando direttamente la teoria delle grandi deviazioni (piuttosto che metodi indiretti a tempo discreto) per stimare le probabilità rare vicino agli insiemi invarianti disordinati.

Spiegazione Dettagliata dei Metodi

Definizione del Compito

Studio delle perturbazioni stocastiche di sistemi di equazioni differenziali ordinarie cooperative:

  • Sistema Deterministico: dxdt=b(x)\frac{dx}{dt} = b(x), dove b:RrRrb: \mathbb{R}^r \to \mathbb{R}^r è continuamente differenziabile, soddisfa la cooperatività (bixj0,ij\frac{\partial b_i}{\partial x_j} \geq 0, i \neq j) e l'irriducibilità (matrice jacobiana irriducibile).
  • Perturbazione Stocastica: dXtε=b(Xtε)dt+εσ(Xtε)dWtdX^\varepsilon_t = b(X^\varepsilon_t)dt + \varepsilon\sigma(X^\varepsilon_t)dW_t, dove ε\varepsilon è un piccolo parametro di perturbazione e WtW_t è un processo di Wiener standard r-dimensionale.
  • Obiettivo: Determinare il supporto del limite di rumore zero μ\mu della misura invariante με\mu_\varepsilon (il limite debole di με\mu_\varepsilon quando ε0\varepsilon \to 0).

Quadro Teorico

1. Fondamenti dei Sistemi Dinamici Monotoni

  • Relazione di Ordine Parziale: xyxiyi,ix \leq y \Leftrightarrow x_i \leq y_i, \forall i; xyxi<yi,ix \ll y \Leftrightarrow x_i < y_i, \forall i
  • Forte Monotonia: Il flusso di soluzione Φt\Phi_t dei sistemi cooperativi irriducibili soddisfa: x<yΦt(x)Φt(y),t>0x < y \Rightarrow \Phi_t(x) \ll \Phi_t(y), \forall t > 0
  • Dicotomia degli Insiemi di Catena Transitiva (Lemma 4.1): Un insieme di catena transitiva KK è o disordinato o un p-arc statico (arco di punti di equilibrio)

2. Teoria delle Grandi Deviazioni di Freidlin-Wentzell

Funzione di Energia di Velocità:

\int_0^T L(\phi(t), \dot{\phi}(t))dt, & \phi \in AC_T \\ \infty, & \text{altrimenti} \end{cases}$$ dove $L(u,\beta) = \frac{1}{2}(\beta - b(u))^T a^{-1}(u)(\beta - b(u))$, $a(x) = \sigma(x)\sigma^T(x)$. **Quasipotenziale**: Il quasipotenziale da $x$ a $y$ è definito come $$V(x,y) = \inf_{T>0} \inf_{\phi \in AC^x_T} \{S^x_T(\phi): \phi(0)=x, \phi(T)=y\}$$ **Principio delle Grandi Deviazioni**: Sotto le ipotesi (H2)-(H3), il sistema soddisfa il principio uniforme delle grandi deviazioni di Freidlin-Wentzell, includendo: - (C) Condizione di Compattezza - (L) Stima Inferiore: $P\{\rho_T(X^{\varepsilon,x}_\cdot, \phi) < \delta\} \geq \exp\{-\frac{S^x_T(\phi)+\gamma}{\varepsilon^2}\}$ - (U) Stima Superiore: $P\{\rho_T(X^{\varepsilon,x}_\cdot, F^x_T(s)) \geq \delta\} \leq \exp\{-\frac{s-\gamma}{\varepsilon^2}\}$ ### Punti di Innovazione Tecnica #### 1. Proprietà di Nullità del Quasipotenziale (Svolta Cruciale) **Proposizione 4.2**: Sia $K$ un insieme di catena transitiva disordinato, allora esiste un attrattore $A = A(K)$ tale che $K \cap A = \emptyset$ e $$V(y,A) = 0, \quad \forall y \in K$$ **Strategia di Dimostrazione**: 1. Costruzione dell'attrattore duale superiore: $A = \Lambda \cap \bigcap_{t>0} \Phi_t[p,+\infty]]$, dove $p = \sup K$ 2. Definizione di $H^+(K) = \{x: y \ll \Phi_s(x) \text{ per qualche } y \in K, s \geq 0\}$ 3. Dimostrazione che $H^+(K) \subset \text{Basin}(A)$ 4. Costruzione di funzioni di interpolazione lineare (LIF) per costruire percorsi con energia quasi nulla **Intuizione Chiave**: Gli insiemi di catena transitiva disordinati si trovano su varietà non monotone, e il trasferimento da tale insieme al suo attrattore duale richiede quasi nessun consumo di energia. #### 2. Stima della Probabilità di Eventi Rari (Lemma 5.1) Per insiemi di catena transitiva disordinati o p-arc statici con punti finali instabili $K$, esiste un attrattore $A$ e costanti $\kappa_2 > \kappa_1 > 0$ tali che: - (i) $\liminf_{\varepsilon \to 0} \varepsilon^2 \log P\{X^{\varepsilon,z}_T \in A_\eta\} \geq -\kappa_1$ - (ii) $\limsup_{\varepsilon \to 0} \varepsilon^2 \log P\{X^{\varepsilon,z}_T \notin A_\eta\} \leq -\kappa_2$ Questo indica che la probabilità di fuga da insiemi instabili verso un attrattore è molto maggiore della probabilità di fuga dall'attrattore. #### 3. Strategia di Dimostrazione del Teorema Principale Attraverso la catena di disuguaglianze: $$\mu_\varepsilon(B_\delta(y)) \leq I_1 \exp\{\frac{\kappa_1+\gamma}{\varepsilon^2}\}, \quad \mu_\varepsilon(A_\eta) \geq I_2 \exp\{\frac{\kappa_2-\gamma}{\varepsilon^2}\}$$ dove $I_1 \leq I_2$ (per proprietà della misura invariante), si ottiene $$\mu_\varepsilon(B_\delta(y)) \leq \exp\{-\frac{(\kappa_2-\kappa_1)-2\gamma}{\varepsilon^2}\} \to 0$$ pertanto $\mu(K) = 0$, escludendo insiemi disordinati e p-arc instabili. ## Configurazione Sperimentale ### Modello di Applicazione: Circuito di Controllo Biochimico Studio del sistema di retroazione positiva di tipo Griffith: $$\begin{cases} \dot{x}_1 = f(x_r) - \alpha_1 x_1 \\ \dot{x}_j = x_{j-1} - \alpha_j x_j, \quad 2 \leq j \leq r \end{cases}$$ dove $f(z) = \text{sgn}(z)\frac{|z|^m}{1+|z|^m}$ (parametro $m \geq 1$). ### Modello di Perturbazione Stocastica $$\begin{cases} dx_1 = (f(x_r) - \alpha_1 x_1)dt + \varepsilon\sigma_1(x_1)dW^1_t \\ dx_j = (x_{j-1} - \alpha_j x_j)dt + \varepsilon\sigma_j(x_j)dW^j_t \end{cases}$$ ### Verifica delle Condizioni Costruzione di una funzione di Lyapunov $V(x) = x^T Bx$ ($B$ matrice definita positiva), verificando: 1. Dissipazione del sistema (H1) 2. Cooperatività e irriducibilità (o forte monotonia) (H2) 3. Condizioni (H3) e (H4) ## Risultati Sperimentali ### Risultato Principale: Classificazione del Limite di Rumore Zero **Proposizione 6.1**: Sotto le ipotesi di tipo Griffith, il limite di rumore zero $\mu$ è una combinazione convessa di misure di Dirac su un numero finito di punti di equilibrio asintoticamente stabili. #### Caso 1: $m = 1$ (Retroazione Lineare) Definendo $\varphi = \prod_{i=1}^r \alpha_i$, $h(z) = \frac{1}{1+z}$: - Se $\varphi \geq 1$: $\mu = \delta_O(\cdot)$ (origine unico punto stabile) - Se $0 < \varphi < 1$: $\mu = \lambda_1 \delta_{h^{-1}(\varphi)V_0}(\cdot) + \lambda_2 \delta_{-h^{-1}(\varphi)V_0}(\cdot)$ #### Caso 2: $m > 1$ (Retroazione Non Lineare) Definendo $\varphi_m = \max_{z>0} h(z) = \frac{m-1}{m} \cdot m^{-\frac{1}{m}}$: - Se $\varphi \geq \varphi_m$: $\mu = \delta_O(\cdot)$ (origine globalmente attrattiva) - Se $0 < \varphi < \varphi_m$: esiste $z_2 > 0$ tale che $$\mu = \lambda_1 \delta_O(\cdot) + \lambda_2 \delta_{z_2 V_0}(\cdot) + \lambda_3 \delta_{-z_2 V_0}(\cdot)$$ ### Scoperte Chiave 1. **Esclusione delle Orbite Periodiche**: Per $m > 1$ e $0 < \varphi < \varphi_m$, il sistema potrebbe possedere orbite periodiche prodotte dalla biforcazione di Hopf (Selgrade, 1982), ma il limite di rumore zero rimane concentrato sui punti di equilibrio, non sulle orbite periodiche. 2. **Dipendenza dai Parametri**: La struttura del limite di rumore zero è completamente determinata dai parametri $\varphi$ e $m$, presentando una chiara struttura di biforcazione. 3. **Simmetria**: A causa della simmetria dispari di $f$, i punti di equilibrio non nulli appaiono sempre in coppia ($\pm z V_0$). ## Lavori Correlati ### Approccio Orbitale per Sistemi Monotoni Stocastici - **Arnold & Chueshov (1998)**: Dinamica a lungo termine di sistemi monotoni stocastici semplificata - **Flandoli et al. (2017)**: Fenomeno di sincronizzazione del rumore - **Chueshov & Scheutzow (2004)**: Struttura di attrattori e misure invarianti ### Approccio Basato sulla Distribuzione - **Freidlin & Wentzell (1970, 2012)**: Teoria delle grandi deviazioni e quasipotenziale, stabilità stocastica delle classi di equivalenza - **Huang et al. (2015-2018)**: Stati stazionari dell'equazione di Fokker-Planck, stabilità stocastica dei sistemi gradiente - **Young (1986, 2002)**: Stabilità stocastica degli attrattori iperbolici, misure SRB come limiti di rumore zero - **Kifer (1974, 1988)**: Perturbazioni stocastiche di sistemi assiomaticamente A ### Teoria dei Sistemi Dinamici Monotoni - **Hirsch (1984-1988)**: Teoria della convergenza per sistemi monotoni - **Smith (1995, 2017)**: Monografia sui sistemi cooperativi - **Mierczyński (1994)**: Teoria dei p-arcs - **Poláčik (1989-1992)**: Flussi monotoni di equazioni paraboliche ### Vantaggi di Questo Articolo 1. Primo trattamento di sistemi che potrebbero possedere infinite classi di equivalenza 2. Non dipende da obiettivi dinamici predefiniti, basato sulla struttura su larga scala del sistema 3. Utilizza direttamente la teoria delle grandi deviazioni a tempo continuo (non metodi di discretizzazione) 4. Fornisce una classificazione completa per sistemi di retroazione positiva ## Conclusioni e Discussione ### Conclusioni Principali 1. **Contributo Teorico**: Stabilisce che il limite di rumore zero dei sistemi cooperativi irriducibili si concentra su un insieme ordinato semplice composto da punti di equilibrio Lyapunov stabili, provando la stabilità stocastica dell'ordine dinamico. 2. **Svolta Tecnica**: Dimostra che il quasipotenziale da insiemi di catena transitiva disordinati su varietà non monotone al loro attrattore duale è zero, rappresentando la chiave per superare le difficoltà nel calcolo dei quasipotenziali ad alta dimensione. 3. **Valore Applicativo**: Classifica completamente il limite di rumore zero dei sistemi di retroazione positiva di tipo Griffith, provando che anche in presenza di orbite periodiche, il limite di rumore zero rimane concentrato su un numero finito di punti di equilibrio stabili. ### Limitazioni 1. **Condizioni di Ipotesi**: - Richiede dissipazione del sistema (H1) - Richiede che la funzione di Lyapunov soddisfi condizioni tecniche (H2)-(H3) - L'irriducibilità (o forte monotonia) è necessaria 2. **Restrizioni di Dimensionalità**: Sebbene la teoria sia applicabile a qualsiasi dimensione finita, la costruzione della funzione di Lyapunov potrebbe essere difficile in dimensioni elevate. 3. **Sistemi Non Cooperativi**: Il metodo non si applica direttamente a sistemi non monotoni o competitivi. 4. **Coefficienti di Ponderazione Sconosciuti**: Sebbene sia determinato il supporto del limite di rumore zero, non sono fornite formule esplicite per i coefficienti $\lambda_i$ della combinazione convessa. ### Direzioni Future 1. **Lavori Successivi della Serie**: Questo articolo è il primo della serie; i lavori successivi studieranno: - Sistemi monotoni più generali (caso non irriducibile) - Stabilità stocastica di sistemi competitivi - Sistemi monotoni dipendenti dal tempo (non autonomi) 2. **Analisi Quantitativa**: - Calcolo dei coefficienti di combinazione convessa $\lambda_i$ - Stima dei tassi di convergenza - Studio della struttura fine della funzione di velocità delle grandi deviazioni 3. **Estensione Applicativa**: - Reti biochimiche più complesse - Dinamica delle reti neurali - Modelli di sistemi ecologici 4. **Metodi Computazionali**: Sviluppo di algoritmi numerici per il calcolo del quasipotenziale e del limite di rumore zero. ## Valutazione Approfondita ### Punti di Forza 1. **Profondità Teorica**: - Primo studio sistematico della stabilità stocastica di sistemi che potrebbero possedere infinite classi di equivalenza - Fusione ingegnosa della teoria dei sistemi dinamici monotoni e della teoria delle grandi deviazioni - La dimostrazione della proprietà di nullità del quasipotenziale ha significato teorico importante 2. **Innovazione Metodologica**: - Utilizzo diretto della teoria di Freidlin-Wentzell a tempo continuo, evitando la complessità della discretizzazione - Non dipende da obiettivi dinamici predefiniti, basato su proprietà strutturali del sistema - La dimostrazione costruttiva fornisce un'immagine geometrica chiara 3. **Rigore Matematico**: - Condizioni di ipotesi esplicite e verificabili - Dimostrazione completa e logica chiara - Lemmi e proposizioni ben organizzati gerarchicamente 4. **Valore Applicativo**: - La classificazione completa dei sistemi di retroazione positiva ha significato pratico - Il risultato di esclusione delle orbite periodiche ha importante significato biologico - Fornisce fondamenti teorici per la modellazione stocastica di reti biochimiche ### Insufficienze 1. **Complessità Tecnica**: - Richiede verifica di molteplici ipotesi tecniche (H1)-(H3) - La costruzione della funzione di Lyapunov potrebbe essere difficile nel caso generale - Soglia di accesso elevata per lettori non specialisti 2. **Informazioni Quantitative Limitate**: - Non fornisce espressioni esplicite per i coefficienti di combinazione convessa - Mancano stime dei tassi di convergenza - Assenza di verifica sperimentale numerica 3. **Ambito di Applicabilità**: - Limitato a sistemi cooperativi irriducibili - L'ipotesi di dissipazione esclude alcuni modelli importanti - Il caso non autonomo non è affrontato 4. **Calcolabilità**: - Il calcolo pratico del quasipotenziale rimane difficile - Mancano algoritmi numerici e implementazioni ### Impatto 1. **Contributo Accademico**: - Inaugura lo studio sistematico della stabilità stocastica dei sistemi monotoni - Fornisce quadro e strumenti per ricerche successive - Connette i due grandi campi della teoria dei sistemi dinamici e dell'analisi stocastica 2. **Valore Pratico**: - Fornisce guida teorica per sistemi biochimici, reti neurali, ecc. - Predice il comportamento a lungo termine del sistema in ambiente rumoroso - Guida la progettazione della robustezza 3. **Riproducibilità**: - I risultati teorici sono verificabili - Gli esempi di sistemi di retroazione positiva sono riproducibili - Necessita di integrazione con codice di implementazione numerica ### Scenari Applicabili 1. **Sistemi Biologici**: - Reti di regolazione genica (circuiti di retroazione positiva/negativa) - Trasduzione del segnale cellulare - Dinamica delle popolazioni 2. **Sistemi Chimici**: - Reti di reazioni catalitiche - Reazioni oscillanti (come la reazione di Belousov-Zhabotinsky) 3. **Sistemi di Ingegneria**: - Progettazione di sistemi di controllo - Dinamica delle reti neurali - Modellazione di sistemi economici 4. **Ricerca Teorica**: - Perturbazioni stocastiche di sistemi monotoni - Applicazioni della teoria delle grandi deviazioni - Fisica statistica di non equilibrio ## Bibliografia (Riferimenti Chiave) 1. **Freidlin & Wentzell (2012)**: Random Perturbations of Dynamical Systems - Fondamenti della teoria delle grandi deviazioni 2. **Hirsch (1985, 1988)**: Lavori pioneristici sui sistemi dinamici monotoni 3. **Smith (1995)**: Monografia Monotone Dynamical Systems 4. **Benäım & Hirsch (1999)**: Algoritmi di approssimazione stocastica per sistemi cooperativi 5. **Huang et al. (2015-2018)**: Serie di lavori su equazione di Fokker-Planck e stabilità stocastica 6. **Mierczyński (1994)**: Teoria dei p-arcs 7. **Selgrade (1979-1982)**: Comportamento asintotico di sistemi di retroazione positiva e biforcazione di Hopf --- **Valutazione Complessiva**: Questo è un articolo di matematica teorica di alta qualità che raggiunge importanti progressi nello studio della stabilità stocastica dei sistemi dinamici monotoni. Attraverso la combinazione ingegnosa della teoria dei sistemi dinamici e della teoria delle grandi deviazioni, risolve difficoltà tecniche di lunga data. Sebbene sia di natura altamente tecnica, fornisce intuizioni profonde per comprendere il comportamento a lungo termine di sistemi complessi in ambienti rumorosi, possedendo importante valore sia teorico che applicativo. Si consiglia che i lavori successivi integrino esperimenti numerici e più casi di applicazione pratica.