morlet小波变换的M函数

lemon_1984

这是个morlet小波变换的m函数希望大家帮我解释一下，我看不太懂，他构造的morlet小波好像跟教材上的函数不一样，有看懂的帮忙指点哈，写出这个函数。

1. %某个尺度的连续小波变换的M函数
   
2. % delta 小波变换的尺度
   
3. % N     小波函数的长度
   
4. % s     原始信号
   
5. % g     原始信号某个尺度下的小波变换系数
   
6. % w0   中心频率
   
7. function g=Singularity_Detection(delta,N,s);
   
8. %原始信号长度
   
9. n=length(s);
   
10. w0=10;
    
11. %构造Morlet小波函数
    
12. 
    
13. for index_x=1:N;
    
14.     x=index_x-(N+1)/2;
    
15.    
    
16.     phi_x(index_x)=exp(-(x.*x)/(2*delta^2))*exp(j*w0*x/delta);
    
17. end;
    
18. phi_x=phi_x/norm(phi_x); %能量归一化
    
19. 
    
20. %对信号做卷积
    
21. g=conv(s,phi_x);  %卷积
    
22. g=wkeep(g,n);     %保持信号长度

复制代码

fudy10

这是自定义morlet小波函数，可能用于优化。

lemon_1984

谢谢回答，最近看小波纠结了好几天了，几个程序，打算找老师取经去。

韵天之色

嗯嗯