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