Sistema dinamico

Da Wikipedia, l'enciclopedia libera.
(Reindirizzamento da Sistemi dinamici)

In fisica, matematica e ingegneria, in particolare nella teoria dei sistemi dinamici, un sistema dinamico è un modello matematico che rappresenta un oggetto (sistema) con un numero finito di gradi di libertà che evolve nel tempo secondo una legge deterministica. Un sistema dinamico viene identificato da un vettore nello spazio delle fasi, lo spazio degli stati del sistema, dove "stato" è un termine che indica l'insieme delle grandezze fisiche, dette variabili di stato, che caratterizzano la dinamica del sistema.

Descrizione[modifica | modifica wikitesto]

Lo studio dei sistemi dinamici rappresenta uno dei più antichi e importanti settori della matematica e della fisica; si tratta di un modello matematico utilizzato per descrivere i sistemi meccanici nell'ambito della meccanica classica e nella sua riformulazione sviluppata dalla meccanica lagrangiana e dalla meccanica hamiltoniana, e che è presente in molti settori dell'ingegneria, come l'automatica e l'ingegneria dei sistemi. Le applicazioni sono molteplici, spaziando dai circuiti elettrici ai sistemi termodinamici.

Alla fine del diciannovesimo secolo, poi, Henri Poincaré osserva la possibilità di un comportamento fortemente irregolare di alcuni sistemi dinamici studiando il problema dei tre corpi: negli anni '50 del secolo successivo, in seguito agli esperimenti numerici del meteorologo Edward Lorenz, che studiando l'atmosfera terrestre rivelò la dipendenza sensibile dalle condizioni iniziali, i risultati di Poincaré vennero presi in grande considerazione dalla comunità scientifica e posero le basi alla teoria del caos. Il comportamento caotico dei sistemi dinamici, la cui controparte matematica può raggiungere gradi di complessità che rendono vincolante l'utilizzo del calcolatore, è stato riscontrato in molti e diversi ambiti dello studio della natura della civiltà umana, tra cui la biologia e l'economia. Si può definire sistema dinamico un sistema la cui modellizzazione matematica può essere espressa da un'equazione differenziale (ordinaria o alle derivata parziali). A partire da questo esistono diversi formalismi matematici utili alla sua descrizione e studio sia in ambito fisico che ingegneristico (ingegneria dei sistemi e automatica).

Si possono identificare due tipologie di sistema dinamico:

  • Se l'evoluzione avviene ad intervalli discreti di tempo il sistema viene chiamato sistema dinamico discreto ed è definito dall'iterazione di una funzione
  • Se l'evoluzione è continua e definita da un'equazione differenziale, il sistema viene chiamato sistema dinamico continuo.

Di particolare importanza sono i sistemi dinamici lineari, i più semplici da analizzare in quanto le equazioni non lineari non sono solitamente risolvibili in modo esatto. Tra i sistemi lineari, i sistemi lineari tempo-invarianti (sistemi LTI) vengono ampiamente utilizzati nella teoria dei segnali e nella teoria del controllo.

Una delle caratteristiche dei sistemi dinamici che viene studiata più spesso è la stabilità. Per esempio, è comune studiare la stabilità in termini di limitatezza delle uscite nei confronti di un ingresso limitato (stabilità esterna), oppure in termini di allontanamento da uno stato di equilibrio (stabilità interna).

Per analizzare matematicamente il comportamento di un sistema dinamico si utilizzano soprattutto due tipologie di descrizione, la rappresentazione in spazio di stato e il formalismo del dominio della frequenza (si veda la funzione di trasferimento nel caso di sistemi stazionari).

Esempi[modifica | modifica wikitesto]

Esempi di sistemi dinamici continui sono:

Esempi di sistemi dinamici discreti sono:

Definizione[modifica | modifica wikitesto]

Nello specifico, per ogni si può definire tale che:

dove:

Ciò rispecchia il fatto che la legge di evoluzione del sistema non cambia essa stessa nel tempo. Le funzioni parametrizzate da , con la legge di composizione , formano un gruppo commutativo ad un parametro. Frequentemente nel caso discreto coincide con , mentre nel caso continuo coincide con .[1]

Il grafico di è la traiettoria del sistema nel tempo e l'insieme:

è l'orbita passante per (ovvero l'immagine del flusso in ).

Un sottoinsieme è detto -invariante se:

In particolare, affinché sia invariante si deve verificare per tutti gli , ovvero il flusso lungo deve essere definito per tutti i punti di ad ogni tempo.

Sistemi fisici[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Meccanica lagrangiana e Meccanica hamiltoniana.

La dinamica dei sistemi fisici può essere caratterizzata dal fatto che il loro moto tra due punti di coordinate generalizzate e segue un cammino che rende stazionario, ovvero a variazione nulla, il funzionale azione:[2]

in accordo con il principio di minima azione (principio variazionale di Hamilton). L'azione è l'integrale nel tempo della lagrangiana :[3]

dove . Si dimostra che così definita soddisfa le equazioni di Eulero-Lagrange:

Rendere stazionaria l'azione corrisponde a minimizzare l'energia del sistema considerato, e solitamente si fa corrispondere all'energia totale del sistema una funzione , detta hamiltoniana e introdotta nel 1835 da William Rowan Hamilton, che dipende dalle coordinate generalizzate e dai rispettivi momenti coniugati:

L'hamiltoniana è data dalla somma dell'energia cinetica e dell'energia potenziale del sistema, ed è la trasformata di Legendre della lagrangiana :[4][5]

La formalizzazione di un problema dinamico tramite il principio di minima azione (valido per sistemi olonomi e monogenici) è alla base della riformulazione della meccanica classica sviluppata dalla meccanica hamiltoniana e lagrangiana.

In particolare le equazioni di Hamilton:

sono equivalenti alle equazioni del moto di Eulero-Lagrange, a loro volta equivalenti alla legge di Newton.[6]

Il principio di conservazione dell'energia viene poi espresso, in tale contesto, dicendo che è un integrale primo delle equazioni di Hamilton, oppure con il fatto che la lagrangiana non dipende esplicitamente dal tempo:

Più in generale, per il teorema di Noether ad ogni simmetria della lagrangiana, ovvero ad ogni trasformazione infinitesima continua delle coordinate che lascia inalterata , corrisponde una quantità conservata.

Esempio[modifica | modifica wikitesto]

In meccanica classica un esempio elementare di sistema dinamico è fornito da un punto che si muove nello spazio. Il punto viene completamente caratterizzato dalla sua posizione (un vettore dipendente da ) e dalla sua velocità . Lo stato di tale sistema è il vettore , dove è lo spazio delle fasi utilizzato. Lo spazio delle fasi viene anche detto spazio delle configurazioni per il fatto che i suoi elementi rappresentano tutti gli stati possibili che il sistema può assumere. L'evoluzione temporale del punto è quindi data dalle due derivate:

dove è l'accelerazione del punto (che dipende dalla somma delle forze a cui è soggetto). Definendo:

il moto del punto può essere scritto con l'equazione ordinaria autonoma:

Scegliendo un punto e una velocità iniziali , ovvero ponendo , si ottiene l'evoluzione del sistema a partire da (problema di cauchy per l'equazione differenziale).

Tutti i sistemi dinamici a tempo continuo vengono scritti in modo analogo, eventualmente con che dipende esplicitamente dal tempo:

dove è una funzione almeno differenziabile. Tale sistema può essere ricondotto a quello autonomo () con un cambio di variabili.

La soluzione al variare di è la traiettoria (orbita) seguita dal sistema nello spazio delle fasi a partire da . Nell'impostare formalmente lo studio di un sistema dinamico si fa in modo che la funzione sia sufficientemente regolare da fornire una soluzione unica (teorema di esistenza e unicità), in accordo con il fatto che l'evoluzione del sistema a partire da un punto dato è unica. In generale, un sistema dinamico è definito da un gruppo (o un semigruppo) , che è l'insieme dei valori del parametro tempo , e un insieme , detto lo spazio delle fasi o spazio degli stati. La funzione di evoluzione temporale (flusso) determina l'azione di su . Nella teoria ergodica è uno spazio misurabile con misura di probabilità e è una funzione misurabile che preserva , mentre nella cosiddetta topologia dinamica è uno spazio topologico completo e è una funzione continua (spesso anche invertibile).[7]

Classificazione e tipologie di sistemi[modifica | modifica wikitesto]

Sistemi continui[modifica | modifica wikitesto]

Data una varietà , sia un campo vettoriale differenziabile, cioè che associa ad ogni punto un vettore le cui coordinate sono legate alle coordinate di (definite in un suo intorno rispetto a qualche base) tramite una funzione differenziabile. Un sistema dinamico è definito dall'equazione autonoma (l'equazione del moto per sistemi meccanici):

Trattandosi di un'equazione differenziale ordinaria, il relativo teorema di esistenza e unicità della soluzione stabilisce che preso un punto iniziale esiste un intervallo , con , in cui il sistema dinamico ha una soluzione unica .

Se la soluzione (traiettoria) esiste per tutti i tempi e per qualsiasi scelta del punto iniziale si ha che il tempo può scorrere nel verso contrario, ovvero è possibile predire il passato conoscendo uno stato del sistema nel futuro. In particolare, si verifica che e l'insieme delle forma un gruppo continuo ad un parametro di diffeomorfismi su .

La struttura matematica che viene assegnata allo spazio delle fasi dipende comunque dal contesto; solitamente è uno spazio topologico, in cui ha senso parlare di continuità nell'evoluzione temporale dello stato. Uno spazio topologico in cui è possibile l'utilizzo di strumenti metrici e differenziali è ad esempio la varietà differenziabile, una delle strutture più utilizzate in quanto risulta particolarmente adatta per modellare i sistemi fisici. Per i sistemi nei quali allo stato viene associata una nozione di misura, ad esempio una probabilità, si utilizza uno spazio misurabile. Si richiede inoltre che il flusso sia compatibile con la struttura di : nel caso in cui sia rispettivamente uno spazio topologico, uno spazio misurabile, una varietà differenziabile o una varietà complessa, è un omeomorfismo, una funzione misurabile, un diffeomorfismo o una funzione olomorfa.

Sistemi discreti[modifica | modifica wikitesto]

I sistemi dinamici discreti sono definiti da un'iterazione del tipo:

di una funzione , con . Può essere vista come un'equazione alle differenze:

che definendo assume la stessa forma dell'equazione differenziale ordinaria del caso continuo.

Le orbite di un sistema discreto sono una successione di stati . Il gruppo di trasformazioni è quindi dato dall'insieme:

dove l'espressione indica la composizione di funzioni di con sé stessa iterata volte.

Classificazione in base a ingressi e uscite[modifica | modifica wikitesto]

In ambito ingegneristico i sistemi dinamici vengono classificati in base al numero di variabili d'ingresso e d'uscita, si hanno infatti:

  • sistemi a singolo ingresso e singola uscita (SISO, dall'inglese single input-single output);
  • sistemi a ingresso multiplo e uscita multipla (MIMO, dall'inglese multiple input-multiple output);

e meno frequentemente:

  • sistemi a singolo ingresso e uscita multipla (SIMO, dall'inglese single input-multiple output);
  • sistemi a ingresso multiplo e singola uscita (MISO, dall'inglese multiple input-single output).

Sistemi lineari[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Sistema dinamico lineare.

Una classe molto importante di sistemi dinamici è quella dei sistemi lineari, in cui il legame tra variabili di ingresso e l'uscita è lineare. Sono utilizzati ad esempio nella teoria dei segnali o nella teoria dei circuiti, e spesso sono analizzati in frequenza tramite l'utilizzo di trasformate integrali, come la trasformata di Fourier o la trasformata di Laplace.

Un sistema lineare di stati , input e uscite viene descritto da un'equazione del tipo:[8]

dove , , e sono matrici (che nel caso stazionario non dipendono dal tempo).

Sistemi lineari e stazionari[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Sistema dinamico lineare stazionario e Sistema dinamico lineare stazionario discreto.

Un sistema dinamico lineare e stazionario è anche detto lineare tempo-invariante, abbreviato spesso con la sigla LTI (dall'inglese Linear Time-Invariant). Nel caso di un sistema continuo, è caratterizzato dal fatto che l'uscita per un segnale in ingresso è descritta dalla convoluzione:

dove è la risposta impulsiva, ovvero la risposta del sistema quando l'ingresso è una funzione a delta di Dirac. Se la funzione è nulla quando allora dipende soltanto dai valori assunti da precedentemente al tempo , ed il sistema è detto causale.

Un sistema a tempo discreto trasforma la successione in ingresso in un'altra successione , data dalla convoluzione discreta con la risposta alla delta di Kronecker:

Gli elementi di possono dipendere da ogni elemento di . Solitamente dipende maggiormente dagli elementi in prossimità del tempo .

I sistemi lineari stazionari sono spesso descritti nel dominio della frequenza (risposta in frequenza) attraverso la funzione di trasferimento, definita come la trasformata di Laplace della risposta all'impulso a Delta.

Sistemi strettamente propri[modifica | modifica wikitesto]

Un ulteriore classificazione per i sistemi lineari li divide in strettamente propri (o puramente dinamici) quando l'uscita dipende esclusivamente dagli stati del sistema, e in tal caso nella rappresentazione matriciale ciò corrisponde a una matrice nulla, mentre si parla di sistema proprio in tutti gli altri casi. Un caso particolare di sistema proprio si ha quando è la matrice ad azzerarsi, in tal caso il sistema è detto non dinamico e non è necessario ricorrere a variabili di stato per rappresentarlo, poiché il legame fra ingresso e uscita è istantaneo.[9] È possibile dimostrare che un sistema puramente dinamico ha funzione di trasferimento con grado del numeratore minore a quello del denominatore mentre un sistema non dinamico ha, ovviamente, funzione di trasferimento con grado zero.

Sistemi non lineari[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Sistema dinamico non lineare.

Sistemi complessi[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Teoria della complessità e Sistema complesso.

Analisi[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Analisi dei sistemi dinamici.

Strumenti matematici e rappresentazioni[modifica | modifica wikitesto]

Dominio del tempo (spazio di stato o delle fasi)[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Spazio di stato.
Traiettorie di stato[modifica | modifica wikitesto]

Supponendo di perturbare un sistema ed osservando la traiettoria di una grandezza di interesse, si verificano casi di particolare interesse quando l'evoluzione tenderà a stabilizzarsi in una posizione di equilibrio, ovvero un punto fisso dell'evoluzione del sistema.

Gli equilibri di un sistema cambiano al variare di ingressi e disturbi (supposti costanti), ad esempio modificando la tensione ai capi di un motore varia la velocità raggiunta a regime. Lo studio degli equilibri di un sistema dinamico è di estremo interesse, tipicamente i problemi di controllo possono essere interpretati come una modifica del punto di equilibrio di un dato sistema. Un esempio semplice è dato dall'equilibrio termico di un appartamento, la cui temperatura interna è l'equilibrio imposto dalle condizioni ambientali ed interne. L'utilizzo di un condizionatore d'aria (sistema di controllo) modificando la temperatura interna alla stanza non fa altro che modificare il punto di equilibrio del sistema.

Dominio della frequenza[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Dominio della frequenza e Rappresentazione spettrale dei segnali.

Rappresentazione grafica[modifica | modifica wikitesto]

Nell'ingegneria dei sistemi un sistema può essere modellizzato graficamente tramite una scomposizione in un insieme di sottosistemi collegati tra loro in vario modo (serie, parallelo, retroazione ecc...), ciascuno dei quali è identificato da uno scatolotto il cui funzionamento o comportamento è descritto da una funzione di sottoprocesso che esso svolge all'interno del sistema generale. Lo schema risultante si darà schema a blocchi del sistema (si veda Modello black box).

Stabilità e punti di equilibrio[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Teoria della stabilità.

Si possono definire diversi tipi di stabilità per un sistema dinamico, ad esempio la stabilità esterna, anche detta stabilità BIBO (da Bounded Input, Bounded Output), ovvero la proprietà di avere un'uscita limitata se l'ingresso è limitato, oppure la stabilità interna, che si riferisce alla capacità di tornare in una configurazione di equilibrio dopo una perturbazione dello stato di equilibrio stesso. La stabilità esterna viene generalmente utilizzata per analizzare il comportamento di sistemi lineari stazionari (per i quali si valutano i poli della funzione di trasferimento), mentre la stabilità interna sfrutta la rappresentazione in spazio di stato del sistema ed è stata studiata in particolare da Aleksandr Michajlovič Ljapunov.

L'analisi della stabilità di un sistema meccanico è collegata con il fatto che il sistema, se lasciato libero di evolvere, tende spontaneamente a portarsi in una configurazione dove la sua energia potenziale è minima: tale configurazione che corrisponde ad uno stato di equilibrio stabile (si veda il teorema di Lagrange-Dirichlet).

Stabilità interna[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Stabilità interna.

Stabilità esterna[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Stabilità esterna.

Un sistema è stabile esternamente (BIBO stabile) se ad un ingresso limitato corrisponde una uscita limitata. La limitatezza di una funzione scalare è generalmente definita in tale contesto dal fatto che esiste un tale che:

Nel caso di sistemi dinamici lineari, un sistema lineare è BIBO stabile se e solo se la risposta impulsiva è assolutamente integrabile, cioè esiste un tale che:[10]

Stabilità strutturale[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Stabilità strutturale.

Controllabilità e osservabilità[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Controllabilità e Osservabilità.

I concetti di controllabilità e osservabilità di un sistema dinamico sono stati introdotti da Kalman nel 1960 e sono alla base della teoria del controllo. Informalmente, un sistema è controllabile se è possibile portarlo in qualsiasi configurazione finale agendo opportunamente sull'ingresso in un tempo finito; viceversa, è osservabile se dall'uscita è possibile risalire allo stato del sistema. Nei sistemi lineari controllabilità e osservabilità sono due proprietà duali.

Sistemi lineari[modifica | modifica wikitesto]

Dato un sistema dinamico lineare:

dove è un vettore costante, si consideri la matrice:

Il sistema è completamente osservabile se il rango di è massimo.

Considerando invece la matrice:

il sistema è completamente controllabile se la matrice ha rango massimo.

Definendo il sistema duale:[11]

si dimostra che il sistema di partenza è completamente osservabile se e solo se il sistema duale è completamente controllabile, ed è completamente controllabile se e solo se il sistema duale è completamente osservabile.

Sistemi non lineari[modifica | modifica wikitesto]

Dato un sistema dinamico definito su una varietà di dimensione :

con l'ingresso, l'uscita e , i problemi di controllabilità si traducono nel verificare se lo spazio delle fasi è sufficientemente grande da contenere tutti gli stati possibili (altrimenti il sistema non è osservabile) o se, al contrario, contiene stati che il sistema non può raggiungere (il sistema non è controllabile).

Una descrizione matematica comunemente utilizzata considera l'algebra di Lie di campi vettoriali sullo spazio delle fasi generata dal campo vettoriale , con un controllo costante: se la dimensione dell'algebra è costante esiste un'unica sotto-varietà tangente lo stato iniziale contenente tutte le orbite raggiungibili dal sistema (andando avanti o all'indietro nel tempo) passanti per . Se la dimensione di è allora e il sistema è in qualche modo controllabile; in caso contrario, se la dimensione è minore di si considera solo l'insieme in cui il sistema è controllabile.[12]

Sistemi ergodici[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Teoria ergodica.

Teoria delle biforcazioni[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Teoria delle biforcazioni.

Un punto fisso (in generale un punto periodico) è un punto nello spazio delle fasi (cioè uno stato) che rimane invariato durante l'evoluzione del sistema. Un insieme invariante è un insieme di stati che viene mandato in sé stesso dall'evoluzione del sistema, eventualmente spostando i singoli stati all'interno dell'insieme, e un attrattore è un insieme invariante a cui le orbite si avvicinano per tempi che tendono all'infinito.

Caos e attrattori[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Teoria del caos e Attrattore.

Teoria del controllo[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Controllo automatico.

Esempio[modifica | modifica wikitesto]

Per introdurre l'analisi di un sistema dinamico possiamo fare riferimento al modello costituito da un serbatoio d'acqua forato. In tale modello fissiamo le variabili e le costanti del sistema che si è creato. Abbiamo:

  • La sezione del serbatoio (S) che rimane costante nel tempo
  • Una costante generale K del liquido considerato che riassume: Densità liquido, dimensione foro ecc
  • Il livello di acqua nel serbatoio x(t) che definiamo come variabile di stato del sistema
  • La portata d'acqua entrante che definiamo come ingresso del sistema u(t)
  • La portata uscente dell'acqua che definiamo uscita del sistema y(t) che è proporzionale alla quantità di liquido sovrastante (ossia livello d'acqua per la sezione del serbatoio) e alla costante del sistema. Infatti y(t)=K*x(t)

Sappiamo che, essendo un serbatoio un sistema dinamico, il suo stato al tempo t è definito sia dalla variabile di ingresso, sia dalla variabile di uscita, sia dallo stato precedente del sistema x(t-∆t). Possiamo quindi definire la formula generale dei sistemi dinamici (del primo ordine: ossia quelli definiti da una sola variabile di uscita) per i quali:
∆x/∆t= A*x(t)+Bu(t)

Se voglio sapere il livello di acqua nel serbatoio all'istante t posso ragionare sulle variabili del sistema:

  1. So che u(t)-Kx(t) corrisponde alla quantità di liquido del serbatoio (quantità entrante meno quantità uscente)
  2. So che tale valore è uguale a S*∆x/∆t (in quanto tale valore corrisponde anch'esso alla variazione di livello di liquido all'interno del serbatoio nell'unità di tempo), quindi
  1. u(t)-Kx(t)= S*∆x/∆t
  2. Ricavo il rapporto ∆x/∆t e ottengo
  1. ∆x/∆t=kx(t)/S + u(t)/S che si ritrova perfettamente con la formula generale dei sistemi di primo ordine.

Se volessimo analizzare graficamente l'andamento dello Stato del sistema potremmo, tramite foglio di calcolo, determinare l'avanzare del sistema in funzione di un intervallo di tempo ∆t che viene scelto "empiricamente" tramite la formula ∆t=0,1/ASS(A) ossia 0,1 fratto il valore assoluto del coefficiente moltiplicante lo stato del sistema nella formula generale dei sistemi.
Graficamente otterrei un iniziale andamento esponenziale del sistema seguito da un EQUILIBRIO dello stato del sistema. Tendenza dei sistemi dinamici è infatti il raggiungimento di uno stato equilibrato che si conservi nel tempo.

Soluzioni numeriche[modifica | modifica wikitesto]

Magnifying glass icon mgx2.svgLo stesso argomento in dettaglio: Metodi di soluzione numerica per equazioni differenziali ordinarie.

Le soluzioni numeriche sono degli algoritmi che permettono di approssimare la soluzione del sistema di equazioni differenziali che costituiscono il modello matematico del sistema. Questi algoritmi sono alla base dei software di simulazione come MATLAB/Simulink ed in linea generale possono risolvere anche problemi che non ammettono soluzioni in forma chiusa.

Note[modifica | modifica wikitesto]

Bibliografia[modifica | modifica wikitesto]

Voci correlate[modifica | modifica wikitesto]

Altri progetti[modifica | modifica wikitesto]

Collegamenti esterni[modifica | modifica wikitesto]