Transformée de Hilbert

Version imprimable

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:

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;

10/05/2011, 22h39
Invité

Bonsoir,

Citation:

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?
...

Code:

1
2
Error in ==> fft_hilb_ref at 12
Xnwh  = xn.* wh;

10/05/2011, 23h53
Invité
Le message d'erreur aurait était mieux en entier, mais bon je suppose que le reste était
```
??? Error using ==> times
Matrix dimensions must agree.
```
Un simple ' devrait régler l'affaire:

Code:

Xnwh = xn.* wh';

Une question cependant: la variable a1 obtenue par la ligne a1= wavread ('nesrine'); est-elle de dimension Nx2?
Si c'est le cas:
Code:

1 2 a1= wavread ('nesrine'); xn = a1(4e4:7e4,1); % ou 2 suivant le canal choisi
serait plus juste, sinon tu prends une valeur sur deux de chaque canal.