freqtemp.m

À propos

Type de fichier

Fichier M de < 1 Ko (text/plain)

Confidentialité

Fichier public, envoyé le 11 mai 2013 à 00:47, depuis l'adresse IP 80.236.x.x (Belgique)

Sécurité

Ne contient aucun Virus ou Malware connus - Dernière vérification: hier

Statistiques

La présente page de téléchargement a été vue 653 fois depuis l'envoi du fichier

Page de téléchargement

https://www.petit-fichier.fr/2013/05/11/freqtemp/   Copier

Aperçu du fichier

function [freqtable]= freqtemp(name,Tsub,nbfrequency)
if nargin < 3
    nbfrequency =16;
end,
if nargin < 2
    Tsub =0.5;
end,
[y,Fs,NBITS]=WAVREAD(name);

NBsample = Tsub*Fs;
i=1;
freqtable(round(size(y,1)/NBsample),nbfrequency)=0;
while i*NBsample < size(y,1)+1
    y1=y((i-1)*NBsample+1:i*NBsample);
    NFFT = 2^nextpow2(NBsample); 
    Y = fft(y1,NFFT)/NBsample;
    factor=0.3;
    while true
        freq=find(abs(Y(1:NFFT/2+1))>abs(factor*max(Y(1:NFFT/2+1))))*Fs/(NFFT);
       
        factor =factor+ 0.1;
        if size(freq,1)<nbfrequency || factor  >1
            break,            
        end,
    end,
    if size(freq,1)< nbfrequency
        freq(nbfrequency)=0;
    end,
    freqtable(i,:)=freq(1:nbfrequency);
    i=1+i;
end,
t=Tsub:Tsub:(size(y,1)/Fs)+Tsub;
t=t(1:size(freqtable,1));
bar(t,freqtable);
xlabel('Temps(s)')
ylabel('Frequency (Hz)')

Partager le fichier

Télécharger freqtemp.m