calccwtscale.sci

小波分解源代码

function [scales,xtick,ytick] = CalcCWTScale(sz)
// CalcCWTScale -- Calculate Scales and TickMarks for CWT Display
//  Usage
//    [scales,xtick,ytick] = CalcCWTScale(sz);
//  Inputs
//    sz      size of CWT array
//  Outputs
//    scales  vector of scales in CWT
//    xtick   vector of positions
//    ytick   vector of log2(scales)
//  
//  Copyright Aldo I Maalouf

	n = sz(1); nscale = sz(2);
	noctave = floor(log2(n))-5;
	nvoice  = nscale ./ noctave;
	scales  = 2.^(2 + (1:nscale)./nvoice);
	xtick   = (.5:(n-.5))./n;
	ytick   = (2 + (1:nscale)./nvoice);
	endfunction