求解小波母函数Morlet的时域中心,时域宽度 代码如下：

求解小波母函数Morlet的时域中心,时域宽度 代码如下：
代码如下：
%% Morlet小波母函数；
% Set effective support and grid parameters.
lb = -5; ub = 5; n = 1000;
% Compute and plot Morlet wavelet.
[psi,x] = Morlet(lb,ub,n,4);
figure;
plot(x,psi);xlabel('x');ylabel(' Morlet'); title('Morlet小波母函数') ;
[t0,T]=loctime(psi')
loctime函数代码：
function [tm,T]=loctime(sig);
%LOCTIME Time localization caracteristics.
% [TM,T]=LOCTIME(SIG) computes the time localization
% caracteristics of signal SIG.
%
% SIG is the signal.
% TM is the averaged time center.
% T is the time spreading.
[sigr,sigc]=size(sig);
if (sigc~=1),
error('The signal must have 1 column');
else
sig2=abs(sig).^2; sig2=sig2/mean(sig2);
t=(1:sigr)';
tm=mean(t.*sig2);
T=2*sqrt(pi*mean((t-tm).^2 .* sig2));
end;
为什么我调用loctime以后求解的时域中心,时宽不对为：t0 = 500.5000 T = 382.5109
小波时域波形