LR Rabiner and RW Schafer, June 3
LR Rabiner and RW Schafer, June 3
LR Rabiner and RW Schafer, June 3
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
DRAFT: L. R. <strong>Rabiner</strong> <strong>and</strong> R. W. <strong>Schafer</strong>, <strong>June</strong> 3, 2009<br />
8.2. HOMOMORPHIC SYSTEMS FOR CONVOLUTION 431<br />
To represent signals as sequences, rather than in the frequency domain as in<br />
Figure 8.4, the characteristic system can be represented as depicted in Figure 8.5<br />
where the log function is surrounded by the DTFT operator <strong>and</strong> its inverse,<br />
where the three operations are defined by the equations<br />
D∗{<br />
}<br />
* • • + + +<br />
F {} log{<br />
}<br />
-1<br />
F { }<br />
( )<br />
ω j<br />
x [n]<br />
xˆ<br />
[ n]<br />
X e<br />
ˆ jω<br />
X ( e )<br />
Figure 8.5: Representation of the characteristic system for homomorphic deconvolution<br />
in terms of discrete-time Fourier transform operators (denoted F{·}<br />
<strong>and</strong> F −1 {·}).<br />
X(e jω ) =<br />
∞<br />
x[n]e jωn<br />
(8.10a)<br />
ˆX(e<br />
n=−∞<br />
jω ) = log{X(e jω )} (8.10b)<br />
ˆx[n] = 1<br />
2π<br />
π<br />
−π<br />
ˆX(e jω )e jωn dω. (8.10c)<br />
In this representation, the characteristic system for convolution is represented<br />
as a cascade of three homomorphic systems: the F operator (DTFT) maps<br />
convolution into multiplication, the complex logarithm maps multiplication into<br />
addition, <strong>and</strong> the F −1 operator (IDTFT) maps addition into addition.<br />
Similarly, Figure 8.6 depicts the inverse characteristic system for convolution<br />
in terms of discrete-time Fourier transform operators F <strong>and</strong> F −1 <strong>and</strong> the<br />
complex exponential.<br />
+<br />
yn ˆ[ ]<br />
F {}<br />
+<br />
−1<br />
D∗<br />
{<br />
+<br />
exp{<br />
}<br />
}<br />
• •<br />
-1<br />
F { }<br />
*<br />
yn [ ]<br />
ˆ jω<br />
Y ( e ) ( )<br />
ω j<br />
Y e<br />
Figure 8.6: Representation of the inverse characteristic system for homomorphic<br />
deconvolution in terms of discrete-time Fourier transform operators.<br />
The mathematical difficulties associated with the complex logarithm revolve<br />
around problems of uniqueness, which are discussed in some detail in [20]. For<br />
our purposes here, it is sufficient to state that an appropriate definition of the<br />
complex logarithm is