path="../csv" t1,s1=lib.step_response([0.75],[1,0.5,1]) t2,s2=lib.step_response([0.75],[1,-2,1]) np.savetxt("%s/fig1.csv" % path,np.matrix([t1,s1]).T,header="t,s", delimiter=",") np.savetxt("%s/fig2.csv" % path,np.matrix([t2,s2]).T,header="t,s", delimiter=",") #system1 t3,s3=lib.step_response([1,3,2],[1,10,33,46,30]) np.savetxt("%s/fig3.csv" % path,np.matrix([t3,s3]).T,header="t,s", delimiter=",") #system2 t4,s4=lib.step_response([1,3],[1,6,10,-2,-15]) np.savetxt("%s/fig4.csv" % path,np.matrix([t4,s4]).T,header="t,s", delimiter=",") #critere du revers w,mag,phase=lib.bode([2],[1,3,3,1]) np.savetxt("%s/fig5_bode.csv" % path,np.matrix([w,mag,phase]).T,header="w,mag,phase", delimiter=",") t5,s5=lib.step_response_BF([2],[1,3,3,1]) np.savetxt("%s/fig5_BF.csv" % path,np.matrix([t5,s5]).T,header="t,s", delimiter=",") w,mag,phase=lib.bode([20],[1,3,3,1]) np.savetxt("%s/fig6_bode.csv" % path,np.matrix([w,mag,phase]).T,header="w,mag,phase", delimiter=",") t6,s6=lib.step_response_BF([20],[1,3,3,1]) np.savetxt("%s/fig6_BF.csv" % path,np.matrix([t6,s6]).T,header="t,s", delimiter=",")
from matplotlib import pyplot as plt sys.path.append('../../control_engineering_lib') import control_engineering as lib path="../csv" num=1; den=[1,3,3,1] t=np.linspace(0,30,101) w,mag,phase=lib.bode(num,den) np.savetxt("%s/bode_sys.csv" % path,np.matrix([w,mag,phase]).T,header="w,mag,phase", delimiter=",") t,s=lib.step_response_BF(num,den,1,t) np.savetxt("%s/step_sys.csv" % path,np.matrix([t,s]).T,header="t,s", delimiter=",") #correcteur P Cnum=2.9 Cden=1 w,mag,phase=lib.bode(Cnum,Cden) np.savetxt("%s/cor1_bode.csv" % path,np.matrix([w,mag,phase]).T,header="w,mag,phase", delimiter=",") w,mag,phase=lib.bode_cor(num,den,Cnum,Cden) np.savetxt("%s/syscor1_bode.csv" % path,np.matrix([w,mag,phase]).T,header="w,mag,phase", delimiter=",") t,s=lib.step_response_BF_cor(Cnum,Cden,num,den,1,t) np.savetxt("%s/cor1_step.csv" % path,np.matrix([t,s]).T,header="t,s", delimiter=",")