消噪加提取基波

kevin19821

小波用的少，在做一个课题时候，需和小波作比较，不知道下面这个程序错在那里？提取的基波怎么和别来的差别那么大
先进行默认阈值消噪，再提取基波，还有小波分解结果的好坏与频率分辨率有何关系
clc
clear all
close all
fs=10000;
n=20480;
dt=1/fs;
t=0:dt:(n-1)*dt;
randn('state',0);
tt=randn(1,n);
s=sin(2*pi*50*t+pi/2-0.0087)+0.02*sin(6*pi*50*t+pi/6)+0.01*sin(10*pi*50*t+pi/8)+0.02*tt;
  s1=sin(2*pi*50*t+pi/2-0.0087);
[c,l]=wavedec(s,7,'db8');
  %%%%%%%===========使用小波默认阈值降噪
  [thr,sorh,keepapp]=ddencmp('den','wv',s)
  s2=wdencmp('gbl',c,l,'db8',7,thr,'s',0)；
[c2 l2]=wavedec(s2,7,'db8');
  si=wrcoef('a',c2,l2,'db8',7);
figure(2)
  plot(si)    ％％％提取的基波信号
  hold on
  plot(s1,'r')％％％原基波信号

[ 本帖最后由 eight 于 2008-3-24 16:28 编辑 ]

kevin19821

每次小波分解的时候每次都是点数减半分解，那如果进行很多层分解的话，那岂不是要求点数很多才行？

kevin19821

晕 怎么没有人帮我看看这个结果为什么不对呀

lanhucx

如果小波基选择得当的话，分解层数为4～5层即可！！！！

byrsky

怎么这么多人做小波分析啊

jinnian

小波分析还是几年前看过一点，全忘了！