|
是不是把[c,l]=wavedec(f',lev,wname);这一句的f 后的' 去掉呢?语法上没f'
我不知道你的f 是什么信号,请看看我调试的,如下:
load leleccum;
indx=1:1024;
x=leleccum(indx);
init=2055615866;
randn('seed',init);
f=x+18*randn(size(x));
wname='sym6';
lev=5;
[c,l]=wavedec(f,lev,wname);
sigma=wnoisest(c,l,1)
sigma=median(abs(sigma))/0.6745;
alpha=2;
thr=wbmpen(c,l,sigma,alpha);
keepapp=1;
xd=wdencmp('gbl',c,l,wname,lev,thr,'s',keepapp);
subplot(2,2,1);
plot(x);
title('yuanxinhao');
subplot(2,2,2);
plot(f);
title('hanzaoxinhao');
subplot(2,2,3);
plot(xd);
title('xiaozaohou');
我也刚学这个,多多指教!:@) |
|