十年网站开发经验 + 多家企业客户 + 靠谱的建站团队
量身定制 + 运营维护+专业推广+无忧售后,网站问题一站解决
小波变换去噪的基本思路可以概括为:利用小波变换把含噪信号分解到多尺度中,小波变换多采用二进型,然后在每一尺度下把属于噪声的小波系数去除,保留并增强属于信号的小波系数,最后重构出小波消噪后的信号。其中关键是用什么准则来去除属于噪声的小波系数,增强属于信号的部分。
成都创新互联坚持“要么做到,要么别承诺”的工作理念,服务领域包括:成都网站制作、成都网站设计、企业官网、英文网站、手机端网站、网站推广等服务,满足客户于互联网时代的衢州网站设计、移动媒体设计的需求,帮助企业找到有效的互联网解决方案。努力成为您成熟可靠的网络建设合作伙伴!
小波有两个显著特点:一是在时域中都具有紧支集或近似紧支集;二是正负交替的波动性。小波分析是将信号分解成一系列小波函数的叠加,而这些小波函数都是由一个母小波通过平移和尺度伸缩得来的。
小波分析理论的一个重要特色是可以进行多分辨率分析。信号可通过多层分解为反映高频信息的细节部分和反映低频信息的概貌部分,通过这种多分辨率分解,信号和噪声通常会有不同的表现,从而达到信嗓分离的目的。金融时间序列去噪处理采用更广泛的方法:非线性阈值处理方法。
非线性阈值处理方法又称小波收缩法,该方法的基本原理是基于小波变换的集中能力。即通过小波变换后有用信号的能量集中于少数小波系数上,而白噪声在小波变换域上仍然分散在大量小波系数之上。因而相对来说,有用信号的小波系数值必然大于那些能量分散且幅值较小的噪声的小波系数值。因此,从谱的幅度上(不是谱的位置)看,有用信号和噪声可以实现分离。
该方法可分为以下3个步骤:
(1)选择合适的正交小波基和分解层数J,对含噪信号进行小波变换分解到J层。
(2)对分解得到的小波系数进行闭值处理,可以使用两种处理方法:硬阈值和软阈值法。硬阈值法保留较大的小波系数并将较小的小波系数置零,即:
(3)软阈值法将较小的小波系数置零,而对较大的小波系数向零进行收缩,即:
学者证明了用软阈值法能使估计信号实现最大均方误差最小化,即去噪后的估计信号是原始信号的近似最优估计。该方法具有广泛的适用性,是应用最为广泛的一种小波去噪方法,其计算速度也很快。
s=zeros(1,100);
s(50:100)=1;
subplot(2,2,1);plot(s);
title('原始信号');
[c,l]=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=[a3
dd3
dd2
dd1];
s1=waverec(c1,l,'db1');
subplot(2,2,2);
plot(s1);grid;
title('强制消噪后的信号');
[thr,sorh,keepapp]=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=[a3
softd3
softd2
softd1];
s3=waverec(c2,l,'db1');
subplot(2,2,4);
plot(s3);grid;
title('给定软阈值消噪后的信号');