24小时热门版块排行榜

返回列表

max3008

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 398
帖子: 86
在线: 55.9小时
虫号: 1630788
注册: 2012-02-21
性别: GG
专业: 计算机网络

[求助] 反变换出了问题，帮我改改了

load leleccum;
indx=1:2048;
x=leleccum(indx);% 将信号中的第1到第2048个采样点赋给x
init=2055615866;
randn('seed',init);% 产生含噪信号
sn=x+18*randn(size(x));
len=length(sn);
L=4;
lsbior=liftwave('db4');
els={'p',[-0.125 0.125],0};
lsnew=addlift(lsbior,els);
xDec=lwt(sn,lsnew,L);
ca=zeros(L,len);
cd=zeros(L,len);

for i=1:L
ca(i,

=lwtcoef('a',xDec,lsnew,L,i); %提取各层近似系数与细节系数
cd(i,

=lwtcoef('d',xDec,lsnew,L,i);
end

for j=1:L %硬阈值处理细节系数
len=length(cd(j,

);
   sigma=median(abs(cd(j)))/0.6745;
r=sigma*sqrt(2*log(len));
   for k=1:len
      if abs(cd(j,k))>r
         cd(j,k)=sign(cd(j,k))*(abs(cd(j,k))-r);
      else
         cd(j,k)=0;
      end
   end
end

ca3=ilwt(ca(4,

,cd(4,

,lsnew,4); %后面是有问题的？？
ca2=ilwt(ca3,cd(3,

,lsnew,3);
ca1=ilwt(ca2,cd(2,

,lsnew,2);
S=ilwt(ca1,cd(1,

,lsnew,1);
??? In an assignment A(I) = B, the number of elements in B and
I must be the same.

Error in ==> ilwt at 57
x(firstIdxAPP:2:end) = a;

回复此楼