Výpočet nespojité časové odezvy (při prvním čtení lze přeskočit)
Výpočet nespojité časové odezvy v počátcích každé fáze.
> | res0:=PraCAn(sc3,TRAN,P=(2,10e3),O=0,SCIDEAL,SAVE=v("2")); |
Výpočet nespojité časové odezvy na konci každé fáze.
> | res1:=PraCAn(sc3,TRAN,P=(2,10e3),O=1,SCIDEAL,SAVE=v("2")); |
> | v210:=[seq([n/10e3,evalf(eval(v("2")[1],res0))],n=0..10)]; |
> | v220:=[seq([(n+1/2)/10e3,evalf(eval(v("2")[2],res0))],n=0..10)]; |
> | v211:=[seq([(n+1/2)/10e3,evalf(eval(v("2")[1],res1))],n=0..10)]; |
> | v221:=[seq([(n+1)/10e3,evalf(eval(v("2")[2],res1))],n=0..10)]; |
Vykreslení časové odezvy včetně výše vypočítaných bodů odezvy. Je evidentní, že výpočet byl proveden korektně.
> | plot(eval([v("2"),v210,v220,v211,v221],res),t=0..1e-3,style=[line,point$4],symbol=[circle$3,cross$2],symbolsize=[20$5],thickness=1,color=[red,green,blue,green,blue]); |