Ciao ragazzi, volevo postarvi un piccolo problema. Ho due sorgenti stazionarie per il calcolo della soluzione particolare e della sua derivata. Scrivo il codice, ovviamente nel dominio del tempo senza passare per i fasori, ma mi da un errore che davvero non capisco. Ho provato a visualizzare i valori del mio vettore X_p e semplicemente mi dice: FAIL. Vi allego il codice che ho utilizzato e l'errore che mi posta nella command window. Qualcuno sarebbe così gentile da chiarirmi questo problema? premetto che le matrici che vado ad utilizzare A_p e B_p sono entrambe 2*2, mentre invece il vettore Up è un vettore colonna.
%calcolo condizioni particolari
%quando t>0 abbiamo due contributi di tipo stazionario
Up=[IS;VS1];
Xpos=inv(A_p)*(-B_p)*Up;
VC_p=Xpos(1);
IL_p=Xpos(2);
X_p=[VC_p;IL_p];
Xp_t0=subs(X_p,t,'0');
dXp_t0=subs(diff(Xp_t0,t),t,'0');
"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""
Error in sym/diff (line 96)
R = mupadmex('symobj::diff', S.s, x.s, n.s);
Error in TracciaIbrida (line 114)
dXp_t0=subs(diff(Xp_t0,t),t,'0');
"""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""""
Grazie per la disponibilità!