moduidwpt.asv

block_mimo_ ofdm_channelestimation编码 很好的程序

function   [t_data,lextop]=moduidwpt(T_data,dep)

 
T_data=reshape(T_data,size(T_data,1)*size(T_data,2)/(2^dep),2^dep);
chn=[2^dep-1:2^(dep+1)-2]; % dep=3:[7:14];dep=2:[3:6]
dwtmode('per');    
  lextop= size(T_data,1)*(2^dep); %小波变换需要的点数, 实际上是星座映射后的符号数
  t = wptree(2,dep,zeros(lextop,1),'db4');
 

         for wavii=1:2^dep
            t = write(t,'cfs',chn(wavii),T_data(:,wavii));  %分别返回节点chn（ii）的系数
 
  	   end
 	t_data=read(wpjoin(t),'data')'; %把小波系数的数据读出来
% a11=size(t_data,1)%10240
% a12=size(t_data,2)%1
  	t_data = upcoef('a',t_data,'db4',6);%filtering
% a21=size(t_data,1)%655738
% a22=size(t_data,2)%1