realhappy 发表于 2006-8-8 22:00

如何用小波重构信号

load leleccum
s=leleccum(1:3920);
ls=length(s);
=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 编辑 ]

雪薇梦寒 发表于 2006-8-9 09:09

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

realhappy 发表于 2006-8-9 10:00

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

qiuzhuwenti 发表于 2007-5-14 11:08

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

ltmx01 发表于 2008-3-19 11:27

原帖由 realhappy 于 2006-8-9 10:00 发表 http://www.chinavib.com/forum/images/common/back.gif
多谢楼上的帮忙，问题解决了。


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

fushibao805825 发表于 2008-4-23 13:42

顺便也发给我！ 谢谢

原帖由 ltmx01 于 2008-3-19 11:27 发表 http://www.chinavib.com/forum/images/common/back.gif



楼主，怎么解决的说一下吧，能把程序发到我邮箱吗？ltmx01@sina.com
楼主顺便把程序发给我吧！万分感谢！   lindefu421@163.com

yehe2000 发表于 2008-4-23 19:12

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

cookie2525 发表于 2008-4-23 23:20

原帖由 realhappy 于 2006-8-9 10:00 发表 http://www.chinavib.com/forum/images/common/back.gif
多谢楼上的帮忙，问题解决了。

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

粤语残片 发表于 2012-11-18 21:53

麻烦楼主帮忙啊 email：pengcheng391514469@126.com{:{05}:}

happy 发表于 2012-11-28 15:35

大家可以参考下面的程序%装载采集的信号leleccum.mat
load leleccum;
%=============================
%将信号中第2000到第3450个采样点赋给s
indx=2000:3450;
s=leleccum(indx);
%=============================
%画出原始信号
subplot(2,2,1);
plot(s);
title('原始信号');
%=============================
%用db1小波对原始信号进行3层分解并提取系数
=wavedec(s,3,'db1');
a3=appcoef(c,l,'db1',3);
d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);
%=============================
%对信号进行强制性消噪处理并图示结果
dd3=zeros(1,length(d3));
dd2=zeros(1,length(d2));
dd1=zeros(1,length(d1));
c1=;
s1=waverec(c1,l,'db1');
subplot(2,2,2);
plot(s1);grid;
title('强制消噪后的信号');
%=============================
%用默认阈值对信号进行消噪处理并图示结果
%用ddencmp函数获得信号的默认阈值
=ddencmp('den','wv',s);
s2=wdencmp('gbl',c,l,'db1',3,thr,sorh,keepapp);
subplot(2,2,3);
plot(s2);grid;
title('默认阈值消噪后的信号');
%=============================
%用给定的软阈值进行消噪处理
softd1=wthresh(d1,'s',1.465);
softd2=wthresh(d2,'s',1.823);
softd3=wthresh(d3,'s',2.768);
c2=;
s3=waverec(c2,l,'db1');
subplot(2,2,4);
plot(s3);grid;
title('给定软阈值消噪后的信号');
页: [1]
查看完整版本: 如何用小波重构信号