Différences

Ci-dessous, les différences entre deux révisions de la page.

--- francois:fft_sur_un_fichier_pcm [2013/04/06 18:54] – francois
+++ francois:fft_sur_un_fichier_pcm [2013/10/23 09:21] (Version actuelle) – francois
@@ Ligne 6: / Ligne 6: @@
 http://www.audiocheck.net/audiofrequencysignalgenerator_sinetone.php
+<file scilab fft_pcm.sce>
+//Composantes fréquentielles d'un signal
+//----------------------------------
+// Construction d'un signal bruité échantilloné à 1000hz
+//        contenant deux fréquences pures à  50 et 70 Hz.
+sample_rate=48000;
+file_path = "/home/francois/develop/calcul/scilab/audio"
+file_name = "sin_1000Hz_-3dBFS_3s.wav"
+afd = mopen(file_path + "/" + file_name, "rb")
+data = []
+// Positionning of the seek pointer to the start of wave
+mseek(44, afd)
+while ~meof(afd)
+  data = [data, mget(10,'s',afd)]
+end
+mclose(afd)
+hsize=max(data)
+data = data / (2^15)
+N=size(data,'*'); //nombre d'échantillons
+y=fft(data, -1);
+DSP=10*log10((1/N*abs(y))^2);
+//y est symétrique, on ne garde que  N/2 points
+f=sample_rate*(0:(N/2))/N; //vecteur de fréquences associé
+n=size(f,'*')
+clf()
+plot(f,DSP(1:n))
+</file>