spectre m(t).py
À propos
- Type de fichier
- Fichier PY de 1 Ko (text/x-python)
- Confidentialité
- Fichier public, envoyé le 31 octobre 2015 à 17:55, depuis l'adresse IP 90.29.x.x (France)
- Sécurité
- Ne contient aucun Virus ou Malware connus - Dernière vérification: 28 heures
- Statistiques
- La présente page de téléchargement a été vue 309 fois depuis l'envoi du fichier
- Page de téléchargement
-
Aperçu du fichier
# importation des librairies
from numpy import pi, cos, linspace, log10
from scipy.fftpack import fft, fftfreq
import matplotlib.pyplot as plt
#constantes du signal
K = 2*pi
f0 = 96
f1 = 238
f2 = 0
f3 = 334
f4 = 37
f5 = 367
A0 = 1.65
A1 = 0.15
A2 = 0.6
t0 = 0
t1 = 10
FreqEch = 1024
PerEch = 1./FreqEch
N = FreqEch*(t1 - t0)
t = linspace(t0, t1, N)
signal = A0*cos(f0*K*t) + A0*cos(f1*K*t) + A1*cos(f2*K*t) + A1*cos(f3*K*t) + A2*cos(f4*K*t) + A2*cos(f5*K*t)
FenAcq = signal.size
signal_FFT = abs(fft(signal)/N)
signal_freq = fftfreq(FenAcq,PerEch)
signal_FFT = signal_FFT[0:len(signal_FFT)//2]
signal_freq = signal_freq[0:len(signal_freq)//2]
plt.subplot(211)
plt.title('Signal et son spectre')
plt.ylim(-(A1+5), A1+5)
plt.plot(t, signal)
plt.xlabel('Temps (s)'); plt.ylabel('Amplitude')
plt.subplot(212)
plt.xlim(-5,f5+5)
plt.plot(signal_freq,signal_FFT)
plt.xlabel('Frequence (Hz)'); plt.ylabel('Amplitude')
plt.title('Signal et son spectre')
plt.show()
Partager le fichier
Télécharger spectre m(t).py
Télécharger le fichier (1 Ko)
