morlet.sci

小波分解源代码

function [psix,x]=Morlet(varargin)
[lhs,rhs]=argn();

if (rhs>3)|(rhs<2), error('*');end
if lhs>2, error('*');end

x = linspace(varargin(1:3));      
psix = exp(-(x.^2)/2) .* cos(5*x);
endfunction