Digital Signal Processing Reference
In-Depth Information
|
H
( e
j2
π
f
) |
1
f
, Hz
0
−
f
s
/ 2 = −500
f
s
/ 2 = 500
(a)
x
(
t
) = sin(
ω
o
t
)
t
0
0.2
x
(
t
) = cos(
ω
o
t
)
(b)
Fig. 3.6 FIR implementation of the digital HT. a The HT magnitude response. b HT applied to
the sinusoid x
ð
t
Þ¼
cos
ð
x
o
t
Þ;
which results in a 90 phase shift. Note that the first few samples of
x
ð
t
Þ
are incorrect, hence ignored in the plot
h=hd.*w_han;
num=h;
den=[1 zeros(1,M-1)];
H=freqz(num,den,f,fs);
Figure
3.6
shows the results of this implementation and its application to a
sinusoid.
3.3.2 The Analytic Signal
Most practical signals are real and have a positive frequency spectrum as well as a
mirror image negative frequency spectrum. It is possible to synthesize a so-called
analytic signal which has a no negative frequency content. The analytic signal
z(t) associated with a real signal x(t) is defined as:
z
ð
t
Þ¼
x
ð
t
Þþ
HT
f
x
ð
t
Þg ¼
x
ð
t
Þþ
j
x
ð
t
Þ
ð
3
:
9
Þ
) Z
ð
f
Þ¼
X
ð
f
Þþ
j
½
jsgn
ð
f
Þ
X
ð
f
Þ¼
X
ð
f
Þ½
1
þ
sgn
ð
f
Þ
) Z
ð
f
Þ¼
2X
ð
f
Þ;
f
0
ð
3
:
10
Þ
0
;
f \0
:
The analytic signal z(t) associated with the original signal z(t) can be generated
in MATLAB using the command:
z=x+j*hilbert(x)
Search WWH ::
Custom Search