Complex frequency B-spline wavelet
[PSI,X] = fbspwavf(LB,UB,N,M,FB,FC)
[PSI,X] = fbspwavf(LB,UB,N,M,FB,FC) returns values of the complex frequency B-Spline wavelet defined by the order parameter M (M is an integer such that 1 ≤ M), a bandwidth parameter FB, and a wavelet center frequency FC.
The function PSI is computed using the explicit expression
PSI(X) = (FB^0.5)*((sinc(FB*X/M).^M).*exp(2*i*pi*FC*X))
on an N point regular grid in the interval [LB,UB].
FB and FC must be such that FC > 0 and > FB > 0.
Output arguments are the wavelet function PSI computed on the grid X.
% Set order, bandwidth and center frequency parameters. m = 2; fb = 0.5; fc = 1; % Set effective support and grid parameters. lb = -20; ub = 20; n = 1000; % Compute complex Frequency B-Spline wavelet fbsp2-0.5-1. [psi,x] = fbspwavf(lb,ub,n,m,fb,fc); % Plot complex Frequency B-Spline wavelet. subplot(211) plot(x,real(psi)) title('Complex Frequency B-Spline wavelet fbsp2-0.5-1') xlabel('Real part'), grid subplot(212) plot(x,imag(psi)) xlabel('Imaginary part'), grid