tp2.py

def spectre(F, Deltat):
    Ne = len(Deltat)
    coeffComplex = np.fft.fft(Deltat)
    coeffReel = np.absolute(coeffComplex)
    freq = []
    for i in range(Ne):
        if Deltat[i] == 0:
            freq.append(0)
        else:
            freq.append(i / Deltat[i])
    plt.plot(freq, coeffReel, "b")
    plt.xlabel("frequence (Hz)")
    plt.ylabel("spectre du signal")
    plt.show()