Salta al contenuto
Note per Studenti Esercizio 1 · risposta di regime e transitoria di un modello del secondo ordine con derivata dell'ingresso (tema d'esame settembre 2025)

Esercizio 1risposta di regime e transitoria di un modello del secondo ordine con derivata dell'ingresso (tema d'esame settembre 2025)

Esame
In questa pagina 4

Testo (tema d'esame settembre 2025, quarto appello, esercizio 3 punto ii). Dato il modello ingresso/uscita lineare e tempo-invariante d2ydt2+4dydt+3y=dudt+2u,\frac{d^2y}{dt^2}+4\frac{dy}{dt}+3y=\frac{du}{dt}+2u, determinare, se esistono, la risposta di regime permanente e la risposta transitoria al segnale di ingresso u(t)=sin⁡(t) δ−1(t)u(t)=\sin(t)\,\delta_{-1}(t), con condizioni iniziali nulle.


Teoria usata: Risposta di regime permanente e risposta transitoriaCon ingresso $u=e^{j\tilde\omega t}\delta_{-1}(t)$ e $\tilde\omega$ non polo di $W$, l'uscita è $y=y_{rp}+y_{tr}$ con $y_{rp}(t)=W(j\tilde\omega)e^{j\tilde\omega t}$ (regime permanente) e $y_{tr}$ combinazione dei modi di $a(s)$ (transitorio, a condizioni iniziali nulle solo dei poli di $W$). Per $u=A\cos(\tilde\omega t+\varphi)$: $y_{rp}=A|W(j\tilde\omega)|\cos(\tilde\omega t+\varphi+\arg W(j\tilde\omega))$. Se il sistema è stabile il transitorio si spegne.Risposta di regime permanente e risposta transitoria →, Laplace per i modelli I-O - fratti semplici e risposta del sistemaPer calcolare l'uscita di un modello I/O con Laplace si scrive $Y(s)=\frac{b(s)U(s)+c(s)}{a(s)}$ (forzata + libera), si decompone in fratti semplici e si antitrasforma con la tabella dei modi: $\frac1{(s-\lambda)^{k+1}}\leftrightarrow\frac{t^k}{k!}e^{\lambda t}$, coppie complesse $\to e^{\sigma t}\cos/\sin\omega t$. I poli di $Y$ determinano i modi dell'uscita; il teorema del valore finale dà il limite a regime se $sY(s)$ ha poli solo in $\mathrm{Re}<0$.Laplace per i modelli I-O - fratti semplici e risposta del sistema →, Stabilità asintotica e BIBO dei modelli I-OIl modello è asintoticamente stabile se la libera tende a zero per ogni stato iniziale: tutte le radici di $a(s)$ hanno $\mathrm{Re}<0$ ($a(s)$ di Hurwitz). È BIBO stabile se a ogni ingresso limitato corrisponde un'uscita forzata limitata: equivale a $\int_0^\infty|w|<\infty$, cioè a poli di $W(s)$ (ridotta) tutti in $\mathrm{Re}<0$. Asintotica $\Rightarrow$ BIBO, non viceversa (cancellazioni polo-zero instabili). Poli immaginari semplici: non BIBO (risonanza).Stabilità asintotica e BIBO dei modelli I-O →.

(1) Funzione di trasferimento e stabilità

I coefficienti dell'equazione danno (Sistemi dinamici e modelli ingresso-uscita SISOUn modello ingresso/uscita SISO LTI è un'equazione differenziale $\sum_{i=0}^na_iy^{(i)}=\sum_{i=0}^mb_iu^{(i)}$ con $n\ge m$ (modello proprio). L'uscita è somma di evoluzione libera (dallo stato iniziale, ingresso nullo) e forzata (da stato nullo, $y_f=w*u$, $Y_f=W(s)U(s)$ con $W=b(s)/a(s)$). Le soluzioni libere sono combinazioni dei modi $t^ke^{\lambda t}$ delle radici di $a(s)$.Sistemi dinamici e modelli ingresso-uscita SISO →) W(s)=s+2s2+4s+3=s+2(s+1)(s+3).W(s)=\frac{s+2}{s^2+4s+3}=\frac{s+2}{(s+1)(s+3)}. Poli −1-1 e −3-3, zero in −2-2: nessuna cancellazione e tutti i poli hanno parte reale negativa. Il modello è BIBO stabile (e anche asintoticamente stabile): la risposta di regime permanente esiste.

(2) Risposta di regime permanente

Con u=sin⁡(ω~t)δ−1u=\sin(\tilde\omega t)\delta_{-1}, ω~=1\tilde\omega=1, si valuta WW in s=js=j: W(j)=j+2(j)2+4j+3=2+j2+4j=(2+j)(2−4j)20=8−6j20=0,4−0,3j.W(j)=\frac{j+2}{(j)^2+4j+3}=\frac{2+j}{2+4j}=\frac{(2+j)(2-4j)}{20}=\frac{8-6j}{20}=0{,}4-0{,}3j. Modulo ∣W(j)∣=0,16+0,09=0,5|W(j)|=\sqrt{0{,}16+0{,}09}=0{,}5, fase arg⁡W(j)=−arctan⁡0,30,4=−arctan⁡0,75=−36,87∘=−0,6435\arg W(j)=-\arctan\frac{0{,}3}{0{,}4}=-\arctan0{,}75=-36{,}87^\circ=-0{,}6435 rad. Quindi yrp(t)=0,5 sin⁡(t−0,6435)=0,4sin⁡t−0,3cos⁡t,y_{rp}(t)=0{,}5\,\sin(t-0{,}6435)=0{,}4\sin t-0{,}3\cos t, dove l'ultima forma si ottiene da 0,5(sin⁡tcos⁡φ−cos⁡tsin⁡φ)0{,}5(\sin t\cos\varphi-\cos t\sin\varphi) con cos⁡φ=0,8\cos\varphi=0{,}8, sin⁡φ=0,6\sin\varphi=0{,}6.

(3) Risposta transitoria

Con condizioni iniziali nulle Y(s)=W(s)1s2+1=s+2(s+1)(s+3)(s2+1)Y(s)=W(s)\frac1{s^2+1}=\frac{s+2}{(s+1)(s+3)(s^2+1)}. Fratti semplici:

  • in s=−1s=-1: (−1)+2(−1+3)((−1)2+1)=12⋅2=14\frac{(-1)+2}{(-1+3)((-1)^2+1)}=\frac1{2\cdot2}=\frac14;
  • in s=−3s=-3: −3+2(−3+1)(9+1)=−1(−2)(10)=120\frac{-3+2}{(-3+1)(9+1)}=\frac{-1}{(-2)(10)}=\frac1{20};
  • la parte con i poli ±j\pm j del seno è 4−3s10(s2+1)\frac{4-3s}{10(s^2+1)} (si trova imponendo l'uguaglianza dei numeratori), che antitrasformata dà 0,4sin⁡t−0,3cos⁡t0{,}4\sin t-0{,}3\cos t: coincide con yrpy_{rp} trovata al punto (2).

Quindi y(t)=0,4sin⁡t−0,3cos⁡t⏟yrp(t)+14e−t+120e−3t⏟ytr(t)(t≥0).y(t)=\underbrace{0{,}4\sin t-0{,}3\cos t}_{y_{rp}(t)}+\underbrace{\tfrac14e^{-t}+\tfrac1{20}e^{-3t}}_{y_{tr}(t)}\qquad(t\ge0). Il transitorio contiene solo i poli di WW (−1-1 e −3-3), come previsto per condizioni iniziali nulle. Controlli: y(0)=−0,3+14+120=0y(0)=-0{,}3+\frac14+\frac1{20}=0 ✓; y′(0)y'(0): yrp′(0)=0,4y_{rp}'(0)=0{,}4 e ytr′(0)=−14−320=−0,4y_{tr}'(0)=-\frac14-\frac3{20}=-0{,}4, somma 00 ✓ (l'ingresso ha u(0)=0u(0)=0, quindi u˙\dot u non contiene impulsi in t=0t=0). Confronto con la simulazione numerica dell'equazione con u=sin⁡tu=\sin t: scarto massimo 4⋅10−84\cdot10^{-8} su t∈[0,15]t\in[0,15].

Grafico interattivo: Uscita y(t) = 0,4 sin t − 0,3 cos t + e^(−t)/4 + e^(−3t)/20: il transitorio si spegne in pochi secondi e resta la sinusoide di regime di ampiezza 0,5

Errori comuni

  • Calcolare W(j)W(j) con il numeratore e il denominatore sbagliati: (j)2=−1(j)^2=-1, quindi il denominatore è −1+4j+3=2+4j-1+4j+3=2+4j.
  • Dire che il regime è 0,5sin⁡(t+36,87∘)0{,}5\sin(t+36{,}87^\circ): la fase di W(j)W(j) è negativa (ritardo).
  • Includere nel transitorio anche termini di yrpy_{rp}, o dimenticare i due termini esponenziali.
  • Trascurare il termine u˙\dot u nell'equazione: contribuisce al numeratore (s+2s+2, zero in −2-2) e quindi cambia modulo e fase di W(j)W(j).

Versione ripasso

Testo. y′′+4y′+3y=u′+2uy''+4y'+3y=u'+2u, u=sin⁡t δ−1u=\sin t\,\delta_{-1}, condizioni nulle: regime e transitorio (settembre 2025).

Teoria collegata