Transformée de Fourier

Transformée de Fourier - Python - Programmation

Marsh Posté le 16-04-2015 à 20:00:03    

Bonjour, j'ai vraiment besoin d'aire pour un programme dont je ne comprends pas l'erreur:  
 

Code :
  1. from numpy import linspace, mean
  2. from scipy.fftpack import fft, fftfreq
  3. import matplotlib.pyplot as plt
  4. from scipy.io.wavfile import read
  6. # nom du fichier à étudier
  7. nomfile = open('7991.wav','r')
  9. # lecture du fichier son 8 Khz mono
  10. rate,signal = read(nomfile)
  12. # définition du vecteur temps
  13. dt = 1./rate
  14. FFT_size = 2**18
  15. NbEch = signal.shape[0]
  17. t = linspace(0,(NbEch-1)*dt,NbEch)
  18. t = t[0:FFT_size]
  20. signal = signal[0:FFT_size]# soustraction de la valeur moyenne du signal
  21. signal = signal - mean(signal)# la fréquence nulle ne nous intéresse pas
  22. signal_FFT = abs(fft(signal))  # on ne récupère que les composantes réelles
  26. signal_freq = fftfreq(signal.size, dt)
  27. signal_FFT = signal_FFT[0:len(signal_FFT)//2]
  28. signal_freq = signal_freq[0:len(signal_freq)//2]
  31. #affichage du signal
  32. plt.subplot(211)
  33. plt.title('Signal reel et son spectre')
  34. plt.plot(t,signal)
  35. plt.xlabel('Temps (s)'); plt.ylabel('Amplitude')
  37. #affichage FFT
  38. plt.subplot(212)
  39. plt.plot(signal_freq,signal_FFT)
  40. plt.xlabel('Frequence (Hz)'); plt.ylabel('Amplitude')
  41. plt.show()


 
J'obtiens ceci comme erreur:
 
File "/media/A620-A5CF/Raspberry Pi/Programmation/ProjetISN.py", line 50, in <module>
plt.plot(signal_freq,signal_FF T)
File "/usr/lib/pymodules/python2.7/matplotlib/pyplot.py", line 2467, in plot
ret = ax.plot(*args, **kwargs)
File "/usr/lib/pymodules/python2.7/matplotlib/axes.py", line 3893, in plot
for line in self._get_lines(*args, **kwargs):
File "/usr/lib/pymodules/python2.7/matplotlib/axes.py", line 322, in _grab_next_args
for seg in self._plot_args(remaining, kwargs):
File "/usr/lib/pymodules/python2.7/matplotlib/axes.py", line 300, in _plot_args
x, y = self._xy_from_xy(x, y)
File "/usr/lib/pymodules/python2.7/matplotlib/axes.py", line 240, in _xy_from_xy
raise ValueError("x and y must have same first dimension" )
ValueError: x and y must have same first dimension
 
Merci d'avance, bonne soirée.

Reply

Marsh Posté le 16-04-2015 à 20:00:03   

Reply

Marsh Posté le 22-04-2015 à 12:41:51    

le code que tu nous as copié vient du fichier ProjetISN.py ?
D'après la trace, l'erreur est ligne 50...
 
D'après l'exception soulevée, c'est un problème de taille de deux tableaux qui ne matchent pas.
 
Que te donne l'affichage suivant ?
 

Code :
  1. print len( signal_freq )
  2. print len( signal_FFT)

Reply

Sujets relatifs:

Leave a Replay

Make sure you enter the(*)required information where indicate.HTML code is not allowed