import sys import numpy as np from matplotlib import pyplot as plt sys.path.append('../../control_engineering_lib') import control_engineering as lib 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])
import sys import numpy as np from matplotlib import pyplot as plt sys.path.append("../../control_engineering_lib") import control_engineering as lib path = "../csv" t, s = lib.step_response([2], [1, 4, 3], amp=1) plt.plot(t, s) # normalized plt.xlabel("t") plt.title("Step response") plt.show() # save file np.savetxt("%s/fig1.csv" % path, np.matrix([t, s]).T, header="t,s", delimiter=",")
import sys import numpy as np from matplotlib import pyplot as plt sys.path.append('../../control_engineering_lib') import control_engineering as lib show_cor=1 #Cnum=12.66 num=3/2.0; den=[3,1] t=np.linspace(0,18,401) t,s=lib.step_response(num,den,2,t) plt.plot(t,s,label="Système orginal") plt.xlabel('Temps (s)') plt.ylabel('Sortie') if show_cor==1: #correcteur P Cnum = float(input('Valeur du gain: ')) Cden=1 t,s=lib.step_response_BF_cor(Cnum,Cden,num,den,1,t) plt.plot(t,2*s,label="Système corrigé") plt.axhline(y=2,color='k',ls='dashed',label="Entrée") plt.legend() plt.show()