02.10.2019 Views

UploadFile_6417

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Properties of the Discrete Fourier Transform 169<br />

real number. In addition, if N is even, then N/2 isalso an integer.<br />

Then from (5.32)<br />

X (N/2) = X ∗ ((−N/2)) N<br />

= X ∗ (N/2)<br />

which means that even the k = N/2 component is also real-valued.<br />

This component is called the Nyquist component since k = N/2 means<br />

that the frequency ω N/2 = (N/2)(2π/N) = π, which is the digital<br />

Nyquist frequency.<br />

The real-valued signals can also be decomposed into their even and odd<br />

components, x e (n) and x o (n), respectively, as discussed in Chapter 2.<br />

However, these components are not N-point sequences and therefore we<br />

cannot take their N-point DFTs. Hence we define a new set of components<br />

using the circular folding discussed above. These are called circular-even<br />

and circular-odd components defined by<br />

{<br />

x ec (n) = △ x(0), n =0<br />

1<br />

2 [x(n)+x ((−n)) N ]= [x (n)+x (N − n)] , 1 ≤ n ≤ N − 1<br />

1<br />

2<br />

{<br />

x oc (n) = △ 0, n =0<br />

1<br />

2 [x(n) − x ((−n)) N ]= 1<br />

2<br />

[x (n) − x (N − n)] , 1 ≤ n ≤ N − 1<br />

(5.33)<br />

Then<br />

DFT [x ec (n)] = Re [X(k)] = Re [X ((−k)) N<br />

]<br />

DFT [x oc (n)] = Im [X(k)] = Im [X ((−k)) N<br />

]<br />

(5.34)<br />

Implication: If x(n) isreal and circular-even, then its DFT is also real<br />

and circular-even. Hence only the first 0 ≤ n ≤ N/2 coefficients are<br />

necessary for complete representation.<br />

Using (5.33), it is easy to develop a function to decompose an N-point<br />

sequence into its circular-even and circular-odd components. The following<br />

circevod function uses the mod function given earlier to implement<br />

the modulo-N operation.<br />

function [xec, xoc] = circevod(x)<br />

% signal decomposition into circular-even and circular-odd parts<br />

% --------------------------------------------------------------<br />

% [xec, xoc] = circevod(x)<br />

%<br />

if any(imag(x) ~= 0)<br />

error(’x is not a real sequence’)<br />

Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).<br />

Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.

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

Saved successfully!

Ooh no, something went wrong!