Řešení vlastního přechodového děje
Zadání obvodu pro přechodovou analýzy včetně počáteční podmínky určené výše.
> | obvod_tran:=" |
> | V1 0 1 SIN 0 5 50/Pi -Pi/4/100 |
> | *B1 0 1 V=5*sin(100*t+Pi/4) |
> | R1 1 2 1000 |
> | C 2 0 5e-6 ic=uC0 |
> | .end": |
Lze vůbec takto?? Asi NE!! Rozhodně to je blbě.
> | reseni_tran:=PraCAn(obvod_tran,tran,SAVE={v("V1"),v("C"),i("C")}); |
Přiřazení výsledků do proměnných.
> | uCt:=evalf(subs(reseni_tran,v("C"))); |
> | iCt:=evalf(subs(reseni_tran,i("C"))); |
Určení celkových funkcí pro obvodové veličiny kapacitoru a jejích vykreslení.
> | uCt_celk:=piecewise(t<0,uCt_pred,uCt); |
> | iCt_celk:=piecewise(t<0,iCt_pred,iCt); |
> | plot([uCt_celk,iCt_celk*1000],t=-0.1..0.2,color=[red,blue],thickness=3,legend=["uC [V]","iC [mA]"],title="Výsledný průběh napětí a proudu kapacitorem."); |