1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37
| fclose('all');
[sn,fe,nbit]=wavread('C:\Users\sofiene\Desktop\restau');
L = 512;
p=0.0009801;
sn=sn(:,1);
l=length(sn);
Nbre_tra = fix(l/L);
n=(0:l-1)/fe;
freq = (0:L-1)*fe/L;
for i=2:Nbre_tra
t = sn( (i-1)*L+1-(L/4) :i*L -(L/4));
S = fftshift(fft(t,L));
S_d=(abs(S)).^2;
S_db = 10*log10(S_d/p);
figure(2),plot(freq, S_db); title('DSP du signal enregistré'); xlabel('Hz'); ylabel('amplitude du signal')
freq1=[125 250 500 1000 1500 2000 3000 4000 6000 8000];
freq1=[freq1 9000:1000:16000];
h=seuil(freq);
[p,f]=perte('homme',60);
f=[f zeros(1,8)];
q=interp1(freq1,f,freq2,'linear');
t=q+h;
v=interp1(freq2,t,freq,'linear');
if (S_db') < v
m = 0 ;
else m = v;
end
figure(3),plot(freq,m);
pause
end |
Partager