<< Other Other contwtgnmir >>

stftb >> stftb > Other > contwtgn

contwtgn

computes a continuous wavelet transform

Calling Sequence

[scalo,f,T,a,wt,wavescaled] = contwtgn(x, fmin, fmax, N, wave)

Parameters

x :

a real vector of size Nx: the signal (in time) to be analyzed

fmin:

a positive scalar in ]0 0.5], the normalized lower frequency bound in (in Hz) of the analyzed signal. When unspecified, you have to enter it at the command line from the plot of the spectrum.

fmax :

a positive scalar in ]0 0.5], the normalized upper frequency bound (in Hz) of the analyzed signal. When unspecified, you have to enter it at the command line from the plot of the spectrum.

N :

a positive integer: number of analyzed voices. When unspecified, you have to enter it at the command line.

wave :

a positive integer: the analyzing wavelet order "wave" derivative of the Gaussian is chosen

scalo :

a real N by Nx matrix: the scalogram (squared magnitude of WT)

f :

a real column vector of size N: the frequency samples (geometrically sampled between fmax and fmin).

T :

a real row vector of size Nx: the time samples (sampling period =1)

a :

a real column vector of size N: the scale vector (geometrically sampled between 1 and fmax/fmin)

wt :

a real N by Nx matrix: the coefficients of the wavelet transform. X-axis corresponds to time (uniformly sampled), Y-axis corresponds to frequency (or scale) samples (geometrically sampled between fmin (resp. fmax/fmin and fmax (resp. 1) First row of wt corresponds to the lowest analyzed frequency.

wavescaled :

when the analyzing wavelet is Morlet or Mexican hat, wavescaled = wave.

Description

contwtgn computes a continuous wavelet transform.

Examples

S = altes(256,0.1,0.45,10000) ;
[scalo,f,T] = contwtgn(S,0.01,0.5,128,8) ;

clf; gcf().color_map= jetcolormap(128);
grayplot(T,f,scalo')

Authors


Report an issue
<< Other Other contwtgnmir >>