freqtemp.m
À propos du fichier
- 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 965 fois depuis l'envoi du fichier
- Page de téléchargement
-
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 freqtemp.m sur le Web et les réseaux sociaux:
Télécharger le fichier freqtemp.m