Discussion :

[noussa006]

j'ai voulu faire la transformée de Hilbert d'un signal audio. Voici le code que j'ai inséré dans matlab.je ne sais pas ou est le problème exactement

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
38
39
40
41
42
43
44
tacq  = 8;
TEC=1e-4;
Te    = TEC;
tn    = 0:Te:tacq-Te;
a1= wavread ('nesrine');
a1=a1';
xn    = (a1(4e4:7e4))'
Nfft  = 2^ceil(log2(length(xn)))
% fenetres rectangulaire
wh    = (hann(length(xn)))';
Xnwh  = xn.* wh;
% ajout de zeros
xnwh  = [Xnwh, zeros(1,Nfft-length(Xnwh))];
Xjfh_fft  = fft(xnwh)/length(xn);
Xjfh_hilbert  = hilbertbig(fft(Xnwh)/length(xn));
% domaine spectral
fmax  = 1/TEC;
df    = fmax/Nfft;
ff1    = 0:df:fmax-df;
% zoom spectral
fz1   = 0;
fz2   = 1500;
% domaine interessant
dbmax = 400;
 
figure(1)
%fft
 
plot(ff1, 20*log10(abs(Xjfh_fft)),'b');
axis([fz1,fz2,-dbmax,40]); grid;
legend('FFT fenetrage de hanning')
%Hilbert
 
figure(2)
 
plot(ff1, 20*log10(abs(Xjfh_hilbert)),'r');
axis([fz1,fz2,-dbmax,40]); grid;
legend('Hilbert_fenetrage de hanning');
figure(3)
plot(ff1, 20*log10(abs(Xjfh_fft)),'b');
hold on;
plot(ff1, 20*log10(abs(Xjfh_hilbert)),'r');
axis([fz1,fz2,-dbmax,100]); grid;
hold off;

[Invité]

Bonsoir,

je ne sais pas ou est le problème exactement

Et bien nous non plus...
Tu obtiens des erreurs? Si oui lesquelles?
Tu n'obtiens pas ce que tu veux?
...