praacond.m
上传用户:haiyisale
上传日期:2013-01-09
资源大小:3246k
文件大小:1k
源码类别:

波变换

开发平台:

Matlab

  1. function [PRCond,AACond] = praacond(Hs,Gs,Ha,Ga)
  2. %PRAACOND Perfect reconstruction and anti-aliasing conditions.
  3. %   If (Hs,Gs) and (Ha,Ga) are two pairs of Laurent polynomials,
  4. %   [PRCond,AACond] = PRAACOND(Hs,Gs,Ha,Ga) returns the values of
  5. %   PRCond(z) and AACond(z) which are the Perfect Reconstruction 
  6. %   and the Anti-Aliasing "values":
  7. %      PRCond(z) = Hs(z) * Ha(1/z)  + Gs(z) * Ga(1/z)
  8. %      AACond(z) = Hs(z) * Ha(-1/z) + Gs(z) * Ga(-1/z)
  9. %
  10. %   The pairs (Hs,Gs) and (Ha,Ga) are associated to perfect
  11. %   reconstruction filters if and only if: 
  12. %      PRCond(z) = 2 and AACond(z) = 0.
  13. %
  14. %   If PRCond(z) = 2 * z^d, a delay is introduced in the 
  15. %   reconstruction process.
  16. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 02-Jul-2003.
  17. %   Last Revision: 06-Jul-2003.
  18. %   Copyright 1995-2004 The MathWorks, Inc.
  19. %   $Revision: 1.1.6.1 $
  20. PRCond = Hs * reflect(Ha) + Gs * reflect(Ga);
  21. AACond = Hs * modulate(reflect(Ha)) + Gs * modulate(reflect(Ga));
  22. % In an equivalent way:
  23. %   PRCond = -z*(Hs(z)*Gs(-z)-Hs(-z)*Gs(z))