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

波变换

开发平台:

Matlab

  1. function Q = mpower(P,pow)
  2. %MPOWER Laurent polynomial exponentiation.
  3. %   MPOWER(P,POW) overloads Laurent polynomial P^POW.
  4. %   For a positive integer POW, Q = mpower(P,POW)  
  5. %   returns the Laurent polynomial Q = P^POW.
  6. %   For a negative integer POW and for a monomial P
  7. %   Q = mpower(P,POW) returns the monomial Q = P^POW.
  8. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 19-Jun-2003.
  9. %   Last Revision: 19-Jun-2003.
  10. %   Copyright 1995-2004 The MathWorks, Inc.
  11. %   $Revision: 1.1.6.2 $  $Date: 2004/04/13 00:38:56 $
  12. if pow~=fix(pow)
  13.     error('Invalid value for power argument.')
  14. end
  15. Q = 1;
  16. if pow>0
  17.     for k = 1:pow , Q = Q*P; end
  18. elseif pow<0
  19.     if ismonomial(P)
  20.         D = 1/P;
  21.         for k = 1:abs(pow) , Q = Q*D; end
  22.     else
  23.         error('For negative powers, P must be a monomial.')
  24.     end
  25. else
  26.     
  27. end