multi_channele.m

B3g_phase2_C语言_Matlab程序及说明

function CH_Data=Multi_CHannelE(ufd,ut,uo,pathnum,Frame_begin,Fram_length,banwid,Fc,V,delaymax)
% delaylength=40;
% Fram_length=1000;
% T=0.05*10e-6;
delaylength=delaymax*banwid;
used=delaylength/pathnum;
T=1/banwid*10e-6;
fdmax=Fc*V/1.08;
tmax=delaymax*10e-6;
qt2=(0.3279*T)^2;
% 10=10;
% ufd=rand(10,1);
% ut=rand(10,1);
% uo=rand(10,1)*2*pi;

fd=fdmax*cos(pi*ufd)/6;
t=-tmax/3*log10(1-ut);

delay=(0:pathnum-1)*used;

for p=1:pathnum
%     CH_Data(p,:)=1.313/sqrt(40)*reshape(sum(exp(j*permute(repmat((repmat(uo,1,Fram_length)+fd*(Frame_begin:Frame_begin+Fram_length-1)*2*pi*T),[1,1,pathnum]),[1 3 2])-repmat((repmat((0:pathnum-1)*T*uses,40,1)-repmat(t,1,pathnum)),[1,1,Fram_length]).^2/2/qt2)),pathnum,Fram_length);
    
    CH_Data(p,:)=1.313/sqrt(40)*sum(diag(exp(j*uo-(delay(p)*T-t).^2/2/qt2))*exp(j*2*pi*fd*(Frame_begin:Frame_begin+Fram_length-1)*T));
    
  plot(real(CH_Data(p,:)))
    pause
   
end
norm_factor=1./sqrt(sum(abs(CH_Data.^2)));
CH_Data=CH_Data.*(ones(pathnum,1)*norm_factor);