Rx=2;
Tx=1;
L=2;%實際多徑數為2
Nc=6;
Cp=3;
j=sqrt(-1);
a=rand(Rx,Tx,(L+1));
e=rand(Tx,Nc+Cp);
E=zeros(Tx,Nc);
for i=1:Tx
    E(i,:)=fft(e(i,1:Nc));
end
E
%pause
f=zeros(Tx,(Nc+L));
f=e(:,(Cp-(L-1)):(Nc+Cp));
%pause
%b=zeros(Rx,Tx*L)
b=[];
%a(:,:,L)
for path=1:L+1
    b=[b,a(:,:,(L-path+2))];
end
%b%轉置后的多徑信道

tem1=zeros(Rx,(Nc+L)*Tx);
tem=zeros(Rx,1);
re=zeros(Rx*Nc,1);
ff=zeros(Rx,(Nc+L)*Tx,Nc);
for i=1:Nc  
    ff(:,:,i)=[zeros(Rx,(i-1)*Tx),b,zeros(Rx,(Nc-i)*Tx)];%擴展后的信道卷積矩陣
    re((i-1)*Rx+1:i*Rx)= ff(:,:,i)*(f.');
    %re
    %pause
end
%re
%Re=fft(re) 
%re
re=reshape(re,Rx,Nc)
%pause
temp=zeros(Rx,Tx,L+1);
H=zeros(Rx,Tx,Nc);
H_pinv=zeros(Tx,Rx,Nc);
%for t=1:N
    for path=1:L+1       
        temp(:,:,path)=a(:,:,path);
    end
    for i=1:Nc  
        for path=1:L+1
            H(:,:,i)=H(:,:,i)+temp(:,:,path)*exp(-j*2*pi*(path-1)*i/Nc);%temp2*exp(-j*2*pi*1*i/carrier_number)+temp3*exp(-j*2*pi*2*i/carrier_number);
        end
       % H
       % pause
            H_pinv(:,:,i)=pinv(H(:,:,i));   
   %         H(:,:,i,t)*H_pinv(:,:,i,t)
    end
    %H
    %end
   % H_pinv
    %ee=zeros(Nc*Rx,1)
    
  
    for jj=1:Nc
        %H_pinv(:,:,jj)
        S_de(:,jj)=H_pinv(:,:,jj)*re(:,jj);
        %pause 
    end
    S_de