18.07.2013 Views

LR Rabiner and RW Schafer, June 3

LR Rabiner and RW Schafer, June 3

LR Rabiner and RW Schafer, June 3

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!