如何用小波重构信号

realhappy

load leleccum
s=leleccum(1:3920);
ls=length(s);
[C,L]=wavedec(s,3,'db1');%多尺度一维分解
cA3=appcoef(C,L,'db1',3);%提取第三层的低频系数
cD3=detcoef(C,L,'db1',3);%提取第三层的高频系数
cD2=detcoef(C,L,'db1',2);%提取第二层的高频系数
cD1=detcoef(C,L,'db1',1);%提取第一层的高频系数
A3=wrcoef('a',C,L,'db1',3);%重构第三层的低频信号
D3=wrcoef('d',C,L,'db1',3);%重构第三层的高频信号
D2=wrcoef('d',C,L,'db1',2);%重构第二层的高频信号
D1=wrcoef('d',C,L,'db1',1);%重构第一层的高频信号
如何用A3,D1,D2,D3,重构信号。请版主和热心人帮忙。谢谢！

[ 本帖最后由 eight 于 2008-3-19 13:23 编辑 ]

雪薇梦寒

s=a3+d1+d2+d3 即可重构

realhappy

多谢楼上的帮忙，问题解决了。

qiuzhuwenti

:victory: 赞楼主的这种做法： 就是要把问题最后咋样了告诉大家一声

ltmx01

原帖由 realhappy 于 2006-8-9 10:00 发表

                               
登录/注册后可看大图

多谢楼上的帮忙，问题解决了。

楼主，怎么解决的说一下吧，能把程序发到我邮箱吗？ltmx01@sina.com

fushibao805825

原帖由 ltmx01 于 2008-3-19 11:27 发表

                               
登录/注册后可看大图

楼主，怎么解决的说一下吧，能把程序发到我邮箱吗？ltmx01@sina.com

  楼主顺便把程序发给我吧！  万分感谢！   lindefu421@163.com

yehe2000

帖子太老，怕是收不到了的:@P

cookie2525

原帖由 realhappy 于 2006-8-9 10:00 发表

                               
登录/注册后可看大图

多谢楼上的帮忙，问题解决了。

万分感谢 我也很想要 email:   rayzhong25@hotmail.com

粤语残片

麻烦楼主帮忙啊 email：pengcheng391514469@126.com

happy

大家可以参考下面的程序

1. %装载采集的信号leleccum.mat
   
2. load leleccum;
   
3. %=============================
   
4. %将信号中第2000到第3450个采样点赋给s
   
5. indx=2000:3450;
   
6. s=leleccum(indx);
   
7. %=============================
   
8. %画出原始信号
   
9. subplot(2,2,1);
   
10. plot(s);
    
11. title('原始信号');
    
12. %=============================
    
13. %用db1小波对原始信号进行3层分解并提取系数
    
14. [c,l]=wavedec(s,3,'db1');
    
15. a3=appcoef(c,l,'db1',3);
    
16. d3=detcoef(c,l,3);
    
17. d2=detcoef(c,l,2);
    
18. d1=detcoef(c,l,1);
    
19. %=============================
    
20. %对信号进行强制性消噪处理并图示结果
    
21. dd3=zeros(1,length(d3));
    
22. dd2=zeros(1,length(d2));
    
23. dd1=zeros(1,length(d1));
    
24. c1=[a3 dd3 dd2 dd1];
    
25. s1=waverec(c1,l,'db1');
    
26. subplot(2,2,2);
    
27. plot(s1);grid;
    
28. title('强制消噪后的信号');
    
29. %=============================
    
30. %用默认阈值对信号进行消噪处理并图示结果
    
31. %用ddencmp函数获得信号的默认阈值
    
32. [thr,sorh,keepapp]=ddencmp('den','wv',s);
    
33. s2=wdencmp('gbl',c,l,'db1',3,thr,sorh,keepapp);
    
34. subplot(2,2,3);
    
35. plot(s2);grid;
    
36. title('默认阈值消噪后的信号');
    
37. %=============================
    
38. %用给定的软阈值进行消噪处理
    
39. softd1=wthresh(d1,'s',1.465);
    
40. softd2=wthresh(d2,'s',1.823);
    
41. softd3=wthresh(d3,'s',2.768);
    
42. c2=[a3 softd3 softd2 softd1];
    
43. s3=waverec(c2,l,'db1');
    
44. subplot(2,2,4);
    
45. plot(s3);grid;
    
46. title('给定软阈值消噪后的信号');

复制代码