求助：小波变换模极大值的求取

lsq001

小波变换模极大值的求取
如何实现呢？按照对小波变换，以及模极大值的理论，对应于信号的奇异点，是有相对应的模极大值的。
如此的话，如何求取模极大值呢。
比如对一个信号s进行4尺度分解，用DB小波。之后在D1尺度上，针对这一尺度小波分解的系数，找到对应的模极大值。
之前也在一些论坛上找到过一些，可惜不是很明白。主要是对编写程序比较头疼吧。
下面是找来的一个小程序
他的思想我能明白
但是是需要一定修改之后才能为所用吧
模极大值是要设置阀值，将信号的每点计算值与相邻点进行比较
可以参考我如下的程序实现：


function md=MaxMode(W,L)

Wmax=0;%求取第一个小波变换的模极大值
tfirst=0;%求第一个小波变换模极大值点的时刻

for i=1:5001
md(i)=0;
end

for(i=1:1:4009)
if(abs(W(i))>abs(Wmax))
Wmax=W(i);
tfirst=i;
end
end
%求取小波变换模极大值点。
Wflag=abs(0.02*Wmax);%设定小波变换模极大值的门坎值
tflag=tfirst+2*L/0.28;
for(i=1:1:tflag)
if((W(i+1)>Wflag)&(W(i+1)-W(i)>Wflag))
md(i+1)=W(i+1);
elseif((W(i+1)<-Wflag)&(W(i+1)-W(i)<-Wflag))
md(i+1)=W(i+1);
else
md(i+1)=0;
end

end

for(i=1:1:tflag)
if((md(i)>0)&(md(i+1)>md(i)))
md(i)=0;
elseif((md(i)<0)&(md(i+1)<md(i)))
md(i)=0;
end
end

W代表信号，L代表线路长度，Wflag=abs(0.02*Wmax)是阀值可以自己定义；
tflag=tfirst+2*L/0.28是选取信号分析区间，这也是根据需求而定

要是有人能明白全部，并且能就此给出仔细诠释
让我好好理解，接受就太好了。在此先谢过了

[ 本帖最后由 suffer 于 2006-10-20 14:55 编辑 ]

lsq001

要是有其它例子程序
或者有做好了的，能借鉴一下就好了
呵呵

lsq001

不是吧
难道这里没有人
做过关于
小波变换模极大值的东西么
还是这个问题容易了
我求助反而显得有点不合时宜了
还往有人给些指点建议

lsq001

我顶起
怎么没人知道啊
还是比较麻烦的说
还望有牛人劳力指点啊

lsq001

郁闷了
好不容易找到一个网站提供程序源代码
交流
可还要先上传5——10个
唉
自己在想办法吧

jimin

这个我不懂，支持一下你，别急，慢慢来
加油！
有什么心得和成果拿出来和大家共享哦

happy

参考下面的函数吧

1. function [signal,swa,swd,ddw,wpeak]=wave_peak(points,level,Lo_D,Hi_D,Lo_R,Hi_R,offset)
   
2. % 该函数用于读取ecg信号，找到小波变换模极大序列
   
3. 
   
4. warning off;
   
5. ecgdata=load('ecg.txt');
   
6. plot(ecgdata(1:points)),grid on,axis tight,axis([1,points,-50,300]);
   
7. signal=ecgdata(1:points)'+offset;
   
8. 
   
9. %  信号的小波变换，按级给出概貌和细节的波形
   
10. [swa,swd] = swt(signal,level,Lo_D,Hi_D);
    
11. figure;
    
12. subplot(level,1,1); plot(real(signal)); grid on;axis tight;
    
13. for i=1:level
    
14.     subplot(level+1,2,2*(i)+1);
    
15.     plot(swa(i,:)); axis tight;grid on;xlabel('time');
    
16.     ylabel(strcat('a   ',num2str(i)));
    
17.     subplot(level+1,2,2*(i)+2);
    
18.     plot(swd(i,:)); axis tight;grid on;
    
19. ylabel(strcat('d   ',num2str(i)));
    
20. end
    
21. 
    
22. %求小波变换的模极大值及其位置
    
23. ddw=zeros(size(swd));
    
24. pddw=ddw;
    
25. nddw=ddw;
    
26. posw=swd.*(swd>0);
    
27. pdw=((posw(:,1:points-1)-posw(:,2:points))<0);
    
28. pddw(:,2:points-1)=((pdw(:,1:points-2)-pdw(:,2:points-1))>0);
    
29. negw=swd.*(swd<0);
    
30. ndw=((negw(:,1:points-1)-negw(:,2:points))>0);
    
31. nddw(:,2:points-1)=((ndw(:,1:points-2)-ndw(:,2:points-1))>0);
    
32. ddw=pddw|nddw;
    
33. ddw(:,1)=1;
    
34. ddw(:,points)=1;
    
35. wpeak=ddw.*swd;
    
36. wpeak(:,1)=wpeak(:,1)+1e-10;
    
37. wpeak(:,points)=wpeak(:,points)+1e-10;
    
38. 
    
39. %按级给出小波变换模极大的波形
    
40. figure;
    
41. for i=1:level
    
42.     subplot(level,1,i);
    
43.     plot(wpeak(i,:)); axis tight;grid on;
    
44. ylabel(strcat('j=   ',num2str(i)));
    
45. end

复制代码

lsq001

嗯，多谢happy
我也最近找到这个代码了
费了不少功夫
现在还在钻研
希望很快就有个结果吧
附带说一句
这个程序是用来处理心电信号的
程序还没有完全读懂

visstic

很有道理！

hyl2323

搂住的心态可以理解，我也是这样的，遇到问题总是迫切需要解决，非常期望有人指点一二。
我也做过一阵子小波，不过应用在模态识别方面，你所说的只是大概了解点，呵呵，多交流。

lsq001

嗯，多谢了
我会抓紧时间的
争取早出结果

appier520

我刚才还发了个类似的帖子提问这个，没想到转眼就看到这个程序，多谢多谢！

lsq001

只是一个大概的思想吧
可以好好借鉴的
呵呵

happy

nod，很多论坛朋友给出来的东西大多只是提供一个参考，是否正确，或者和你需要的有多大差距完全可以探讨