关于复小波cmor时频图的问题

pigzhu

用工具箱中tfrscalo函数做时频分析，主要是基于morlet或Mexcian Hat函数。请教各位大侠，怎样画出基于复morlet小波函数cmor系列的时频图呢？
看了一下tfrscalo函数的源程序，关键的这几句没看懂，请知情人指点，多谢！
源程序第135行至141行，也就是求tfr的这几行。
nha = wave*a(ptr);
  tha = -round(nha) : round(nha);
  ha  = exp(-(2*log(10)/nha^2)*tha.^2).*exp(i*2*pi*f(ptr)*tha);
  detail = conv(z,ha)./sqrt(a(ptr));
  detail = detail(round(nha)+1:length(detail)-round(nha)) ;
  wt(ptr,:)  = detail(time) ;
  tfr(ptr,:) = detail(time).*conj(detail(time)) ;

这里ha是什么？morlet的小波函数表达式为exp(-t^2/2)*exp(iwt)或者 exp(-t^2/2)*cos(5t),与ha的表达式相差较远。
如果要改成是cmor小波又该如何做呢？cmor小波函数为1/sqrt(pi*Fb)*exp(2*i*Fc*x)*exp(x^2/Fb)，其中Fb Fc分别为带宽参数和小波中心频率。

谢谢各位的帮忙！