home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l455 / 8.ddi / SIGNAL.DI$ / CCEPS.M < prev    next >
Encoding:
Text File  |  1993-03-11  |  533 b   |  18 lines
  1. function xhat = cceps(x)
  2. %CCEPS    Complex cepstrum.
  3. %    CCEPS(x) returns the complex cepstrum of the sequence x.
  4. %    See also RCEPS, HILBERT, and FFT.
  5.  
  6. %    L. Shure 6-9-88
  7. %    Copyright (c) 1988 by the MathWorks, Inc.
  8.  
  9. %    References: 
  10. %      [1] A.V. Oppenheim and R.W. Schafer, Digital Signal 
  11. %          Processing, Prentice-Hall, 1975.
  12. %      [2] Programs for Digital Signal Processing, IEEE Press,
  13. %          John Wiley & Sons, 1979, algorithm 7.1.
  14.  
  15. h = fft(x);
  16. logh = log(abs(h))+sqrt(-1)*rcunwrap(angle(h));
  17. xhat = real(ifft(logh));
  18.