dyadlength.sci

小波分解源代码

function [n,J] = dyadlength(x)
// dyadlength -- Find length and dyadic length of array
//  Usage
//    [n,J] = dyadlength(x)
//  Inputs
//    x    array of length n = 2^J (hopefully)
//  Outputs
//    n    length(x)
//    J    least power of two greater than n
//
//  Side Effects
//    A warning is issued if n is not a power of 2.
//
//  See Also
//    quadlength, dyad, dyad2ix
//
//  Copyright Aldo I Maalouf

  n = length2(x) ;
  J = ceil(log(n)/log(2));
  if 2^J ~= n ,
      error('Warning in dyadlength: n != 2^J')
  end
  endfunction