Calcolatore Funzione di Trasferimento
Guida Completa al Calcolo della Funzione di Trasferimento
La funzione di trasferimento è un concetto fondamentale nell’analisi dei sistemi dinamici lineari tempo-invarianti (LTI). Rappresenta il rapporto tra l’uscita di un sistema e il suo ingresso nel dominio di Laplace, fornendo una descrizione completa del comportamento del sistema.
1. Definizione Matematica
Una funzione di trasferimento G(s) è definita come:
G(s) = Y(s) / U(s) = N(s) / D(s)
Dove:
- Y(s) è la trasformata di Laplace dell’uscita
- U(s) è la trasformata di Laplace dell’ingresso
- N(s) è il polinomio al numeratore
- D(s) è il polinomio al denominatore
2. Proprietà Fondamentali
| Proprietà | Descrizione | Formula |
|---|---|---|
| Guadagno Statico | Valore della funzione di trasferimento a regime (s=0) | G(0) = N(0)/D(0) |
| Poli | Radici del denominatore che determinano la stabilità | D(s) = 0 |
| Zeri | Radici del numeratore che influenzano la risposta | N(s) = 0 |
| Ordine del Sistema | Grado del polinomio al denominatore | deg(D(s)) |
3. Analisi della Stabilità
La stabilità di un sistema LTI può essere determinata analizzando la posizione dei poli della funzione di trasferimento nel piano complesso:
- Sistema Stabile: Tutti i poli hanno parte reale negativa (si trovano nel semipiano sinistro)
- Sistema Instabile: Almeno un polo ha parte reale positiva
- Sistema Marginalmente Stabile: Polo sull’asse immaginario (parte reale zero)
Il criterio di Routh-Hurwitz fornisce un metodo sistematico per determinare la stabilità senza dover calcolare esplicitamente i poli:
- Costruire la tabella di Routh dal denominatore D(s)
- Contare il numero di cambi di segno nella prima colonna
- Il numero di cambi di segno corrisponde al numero di poli con parte reale positiva
4. Diagrammi di Bode
I diagrammi di Bode rappresentano graficamente:
- Diagramma del Modulo: 20·log|G(jω)| in dB vs log(ω)
- Diagramma della Fase: ∠G(jω) in gradi vs log(ω)
| Elemento | Funzione di Trasferimento | Diagramma del Modulo | Diagramma della Fase |
|---|---|---|---|
| Guadagno K | K | 20·log|K| (linea orizzontale) | 0° (se K>0) o -180° (se K<0) |
| Polo semplice | 1/(τs+1) | -20·log√(1+ω²τ²) | -arctan(ωτ) |
| Zero semplice | (τs+1) | 20·log√(1+ω²τ²) | arctan(ωτ) |
| Integratore | 1/s | -20·log|ω| | -90° |
5. Applicazioni Pratiche
Le funzioni di trasferimento trovano applicazione in numerosi campi:
- Controllo Automatico: Progetto di regolatori PID, analisi della risposta temporale
- Elettronica: Analisi di filtri attivi e passivi, amplificatori
- Meccanica: Modellazione di sistemi massa-molla-smorzatore
- Economia: Modelli dinamici di sistemi economici
- Biologia: Modelli farmacocinetici e fisiologici
6. Limiti e Considerazioni
Sebbene estremamente utile, la funzione di trasferimento presenta alcune limitazioni:
- Si applica solo a sistemi lineari tempo-invarianti
- Non fornisce informazioni sulle condizioni iniziali
- Non rappresenta direttamente le variabili di stato interne
- Può diventare complessa per sistemi di ordine elevato
Per sistemi non lineari, si ricorre spesso alla linearizzazione attorno a un punto di equilibrio. Per sistemi varianti nel tempo, sono necessari altri approcci come le equazioni di stato.
Metodologie di Calcolo
1. Da Equazioni Differenziali
Partendo dalle equazioni differenziali che descrivono il sistema:
- Applicare la trasformata di Laplace assumendo condizioni iniziali nulle
- Risolvere per Y(s) in funzione di U(s)
- Calcolare il rapporto G(s) = Y(s)/U(s)
Esempio: Consideriamo il sistema descritto da:
2ẍ + 3ẋ + x = 4u(t)
Applicando la trasformata di Laplace con condizioni iniziali nulle:
(2s² + 3s + 1)X(s) = 4U(s)
La funzione di trasferimento risulta:
G(s) = X(s)/U(s) = 4 / (2s² + 3s + 1)
2. Da Schemi a Blocchi
Per sistemi rappresentati mediante schemi a blocchi:
- Identificare le funzioni di trasferimento dei singoli blocchi
- Applicare le regole di algebra dei blocchi:
- Serie: G(s) = G₁(s)·G₂(s)
- Parallelo: G(s) = G₁(s) ± G₂(s)
- Retroazione negativa: G(s) = G₁(s)/(1 ± G₁(s)H(s))
- Semplificare l’espressione risultante
3. Da Dati Sperimentali
Quando non si dispone di un modello matematico:
- Acquisire la risposta all’impulso o al gradino del sistema
- Applicare tecniche di identificazione sistemistica:
- Metodo della risposta al gradino
- Analisi in frequenza
- Algoritmi di ottimizzazione (es. minimi quadrati)
- Validare il modello ottenuto
Strumenti Software per il Calcolo
Numerosi strumenti software facilitano il calcolo e l’analisi delle funzioni di trasferimento:
| Strumento | Caratteristiche | Vantaggi | Limitazioni |
|---|---|---|---|
| MATLAB/Simulink | Ambiente completo per l’analisi dei sistemi dinamici |
|
Costo elevato della licenza |
| Python (SciPy, Control) | Librerie open-source per il controllo automatico |
|
Curva di apprendimento più ripida |
| Scilab | Alternativa open-source a MATLAB |
|
Interfaccia meno raffinata |
| Octave | Compatibile con MATLAB |
|
Prestazioni inferiori per sistemi complessi |
Errori Comuni e Come Evitarli
Nel calcolo delle funzioni di trasferimento è facile incorrere in errori. Ecco i più frequenti:
- Dimenticare le condizioni iniziali:
Quando si applica la trasformata di Laplace, è fondamentale assumere condizioni iniziali nulle per ottenere la funzione di trasferimento. In caso contrario, il risultato conterrà termini aggiuntivi dipendenti dalle condizioni iniziali.
- Errata semplificazione dei polinomi:
Nella semplificazione di N(s)/D(s), assicurarsi che numeratore e denominatore non abbiano fattori comuni. Una semplificazione errata può portare a conclusioni sbagliate sulla stabilità del sistema.
- Confondere poli e zeri:
I poli sono le radici del denominatore, mentre gli zeri sono le radici del numeratore. Invertire questi concetti porta a errori nell’analisi della risposta del sistema.
- Trascurare l’unità di misura:
La funzione di trasferimento è adimensionale solo se ingresso e uscita hanno le stesse unità di misura. In caso contrario, il guadagno statico avrà delle unità di misura che devono essere considerate.
- Ignorare i ritardi puri:
I sistemi con ritardi puri (e-sT) richiedono un trattamento speciale. Un ritardo puro non può essere approssimato semplicemente con un polo o uno zero.
Casi Studio Reali
1. Sistema Massa-Molla-Smorzatore
Consideriamo un sistema meccanico composto da:
- Massa m = 2 kg
- Costante elastica k = 100 N/m
- Coefficiente di smorzamento b = 5 N·s/m
L’equazione differenziale che descrive il sistema è:
mẍ + bẋ + kx = F(t)
Applicando la trasformata di Laplace con condizioni iniziali nulle:
(2s² + 5s + 100)X(s) = F(s)
La funzione di trasferimento risulta:
G(s) = X(s)/F(s) = 1 / (2s² + 5s + 100)
Analisi:
- Poli: s = [-5 ± √(25 – 800)]/4 = -1.25 ± 13.98i
- Sistema stabile (poli nel semipiano sinistro)
- Frequenza naturale: ωn = √(100/2) = 7.07 rad/s
- Fattore di smorzamento: ζ = 5/(2√(2·100)) = 0.177
2. Filtro RC Passa-Basso
Un semplice filtro RC con:
- Resistenza R = 10 kΩ
- Condensatore C = 1 μF
La funzione di trasferimento è:
G(s) = Vout(s)/Vin(s) = 1 / (RCs + 1) = 1 / (0.01s + 1)
Analisi:
- Polo: s = -100 rad/s
- Frequenza di taglio: ω = 1/RC = 100 rad/s ≈ 15.9 Hz
- Risposta in frequenza: attenuazione di -20 dB/decade oltre la frequenza di taglio