eh, si, è vero, oggi ho battuto un po' la fiacca, però qualcosina l'ho fatta, e sono arrivato ad un punto in cui credo di aver bisogno di un'altra spintarella per poter procedere..
Allora, ricapitolando, dopo aver stabilito il valore della corrente dell'induttore e della tensione del condensatore all'istante iniziale:
% Calcolo della tensione del condensatore e della corrente dell'induttore
% in t=0-
il_0=subs(il_t_hibrid_tneg,t,0);
vc_0=subs(vC_t_hibrid_tneg,t,0);
Ho determinato il polinomio caratteristico, che mi servirà successivamente per il calcolo del transitorio.. E qui, se non è troppo, chiedo un confronto sulla correttezza dello svolgimento:
%Determinazione del polinomio caratteristico
syms s
ec_s=vc-ic/(s*C)
el_s=vl-s*L*il
sol_reg_tpos_is=solve(lkt1,lkc1,lkc2,db1,db2,ec,el, 'V2, I1, I2, I3, Vl, Ic, Il');
FiL_tpos_is=subs(sol_sin_is_tneg.Il,[Is Vs],[Fis 0]);
HiL_tpos_is=simplify(FiL_tpos_is/Fis);
HiL_s=sym2tf(HiL_tpos_is);
poli=pole(HiL_s);
Fatto ciò, sono finalmente passato alla soluzione del circuito per t>0, analizzandolo sempre con il generatore di tensione stazionario e quello di corrente sinusoidale, e sono arrivato a queste soluzioni:
%c) Calcolo della soluzione in regime ibrido per t>0
%c.1) Soluzione rispetto al generatore di tensione stazionario
sol_hibrid_vs_tpos=solve(lkt1,lkc1,lkc2,db1,db2, 'V2, I1, I2, I3, Il');
il_hibrid_vs_tpos=subs(sol_hibrid_vs_tpos.Il, [Vs Vl Ic Is], [Vs 0 0 0]);
vc_hibrid_vs_tpos=subs(sol_hibrid_vs_tpos.V2, [Vs Vl Ic Is], [Vs 0 0 0]);
vr_hibrid_vs_temp_tpos=subs(sol_hibrid_vs_tpos.I1, [Vs Vl Ic Is], [Vs 0 0 0]);
vr_hibrid_vs_tpos=vr_hibrid_vs_temp_tpos*R;
%c.2) Soluzione rispetto al generatore di corrente sinusoidale
sol_sin_is_tpos=solve(lkt1,lkc1,lkc2,db1,db2,ec,el, 'V2, I1, I2, I3, Vl, Ic, Il');
Fil_is_tpos=subs(sol_sin_is_tpos.Il,[Is Vs],[Fis 0]);
Fvc_is_tpos=subs(sol_sin_is_tpos.V2,[Is Vs],[Fis 0]);
Fvr_is_tpos=subs(sol_sin_is_tpos.Vl,[Is Vs],[Fis 0]);
%c.3) Calcolo delle funzioni sinusoidali
il_t_is_tpos=abs(Fil_is_tpos)*cos(w*t+angle(Fil_is_pos));
vc_t_is_tpos=abs(Fvc_is_tpos)*cos(w*t+angle(Fvc_is_tpos));
vr_t_is_tpos=abs(Fvr_is_tpos)*cos(w*t+angle(Fvr_is_tpos));
%c.4) Soluzioni complete per il regime ibrido per t>0
il_t_hibrid_tpos=il_hibrid_vs_tpos+il_t_is_tpos;
vc_t_hibrid_tpos=vc_hibrid_vs_tpos+vc_t_is_tpos;
vr_t_hibrid_tpos=vr_hibrid_vs_tpos+vr_t_is_tpos;
Anche qui, se avete tempo, voglia e pazienza, gradirei un confronto sulla validità di ciò che ho scritto!
A questo punto, sempre dando un'occhio al codice scritto dal prof relativamente ad un progetto simile, si giunge al calcolo delle condizioni iniziali.. Immagino di non poter copiare spudoratamente, anche perchè la topologia del circuito è abbastanza diversa, quindi dovrò lambiccarmi il cervello per scrivere le mie condizioni iniziali.. e la vedo nera... quindi, anche qui, sarei felicissimo di un confronto... ma credo che questo dolce lavoro lo farò più tardi, magari dopo cena, perchè al momento altri impegni mi chiamano
p.s. Folgore, grazie mille dei feedback che ogni volta, con pazienza e gentilezza mi lasci!
Peppe, mi fa sempre molto piacere dare feedback su una materia che sento, permettimelo di dirlo, particolarmente mia.....
Ti dico la verità, ho guardato il codice "a campione", ovvero ho guardato solo alcune parti. Mi pare tutto ok, sulla programmazione non c'è che dire.....Bisogna vedere se nasce da un ragionamento sensato. Ma credo di si.
Giusto un'osservazione piccolissima: ho capito il calcolo del polinomio caratteristico, ma meglio usare questa porzione di codice dopo.....
Se il prof legge il progetto, vede una logica che improvvisamente si blocca perché il calcolo del polinomio riguarda una parte successiva della tua progettazione.
Credo che, quando depurerai il codice, sarebbe logico spostare quella parte....conoscendo Femia, non mi stupirebbe se ti chiedesse perché hai messo quella parte di codice lì e non altrove.
Ovviamente, è sempre un consiglio.
Ciao.