Just for fun, ti presento una soluzione basata sulla famosa formula di Lagrange della teoria dei sistemi, molto diffusa in automatica.
Prima di tutto faccio l'equivalente di Thevenin alla maglia di destra,
.
Poi, traduco il circuito in equazioni di stato,
[unparseable or potentially dangerous latex formula],
[unparseable or potentially dangerous latex formula].
Sostituisco i valori numerici e riscrivo il sistema in forma matriciale,
[unparseable or potentially dangerous latex formula].
Ricordo che le condizioni iniziali sono:
[unparseable or potentially dangerous latex formula].
La formula di Lagrange per un modello in variabili di stato, scritto nella forma,
[unparseable or potentially dangerous latex formula],
è
[unparseable or potentially dangerous latex formula].
Dove, il termine

è la matrice di transizione dello stato, che ho calcolato usando matlab con il seguente script;
- Code: Seleziona tutto
A=[-33.33, 12.5; 11.11, -27.78];
syms t;
digits(4);
expm(A*t);
vpa(ans)
che restituisce
- Code: Seleziona tutto
>> transA
ans =
[ 0.3854*exp(-18.45*t) + 0.6146*exp(-42.66*t), 0.5162*exp(-18.45*t) - 0.5162*exp(-42.66*t)]
[ 0.4588*exp(-18.45*t) - 0.4588*exp(-42.66*t), 0.6146*exp(-18.45*t) + 0.3854*exp(-42.66*t)]
>>
A questo punto ho calcolato l'integrale di convoluzione, per determinare la risposta forzata. Vale la pena notare che non è necessario calcolare tutto quanto, infatti noi siamo interessati alla

.

.
Quindi sostanzialmente a noi interessa calcolare l'integrale,
[unparseable or potentially dangerous latex formula].
Mettendo assieme i tasselli, si ottiene:
[unparseable or potentially dangerous latex formula].