plotdelay.m

using Jakes model to simulate rayleigh channel fading

%%%%%%%%%%%%%    Delay Profile     %%%%%%%%%%%%%%% 
function Plot_Delay=PlotDelay  %%%%   
 
%%%%%%%%%%%%%%%%%%    Templet    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
Templet_ITUR_M1225_PedA=[0,110e-9,190e-9,410e-9;0,-9.7,-19.2,-22.8];   %% Time(ns); Power(dB) 
Templet_ITUR_M1225_VehA=[0,310e-9,710e-9,1090e-9,1730e-9,2510e-9;0,-1.0,-9.0,-10.0,-15,-20];   %% Time(ns); Power(dB) 
 
%%%%%%%%%%%%%%%%%%    Time Domain    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
figure(7); 
stem(Templet_ITUR_M1225_PedA(1,:),Templet_ITUR_M1225_PedA(2,:)); 
Title=title('Impulse Response'); 
set(Title,'FontSize',14); 
X=xlabel('Time (s) '); 
set(X,'FontSize',14); 
Y=ylabel('Normalized Average Power (dB)'); 
set(Y,'FontSize',14); 
grid on; 
 
figure(8); 
stem(Templet_ITUR_M1225_VehA(1,:),Templet_ITUR_M1225_VehA(2,:)); 
Title=title('Impulse Response'); 
set(Title,'FontSize',14); 
X=xlabel('Time (s) '); 
set(X,'FontSize',14); 
Y=ylabel('Normalized Average Power (dB)'); 
set(Y,'FontSize',14); 
grid on; 
 
%%%%%%%%%%%%%%%%%%    Freq Domain    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
dt=10e-9; 
fs=1/dt; 
 
figure(9); 
width=max(max(Templet_ITUR_M1225_PedA)); 
samples=zeros(1,max(max(Templet_ITUR_M1225_PedA))/dt); 
for i=1:length(samples) 
    for j=1:length(Templet_ITUR_M1225_PedA(1,:)) 
       if i*dt==Templet_ITUR_M1225_PedA(1,j) 
        samples(i)=10^(Templet_ITUR_M1225_PedA(2,j)/20); 
       end; 
    end; 
end; 
PlotSpectrum2(samples,width); 
 
figure(10); 
width=max(max(Templet_ITUR_M1225_VehA)); 
samples=zeros(1,max(max(Templet_ITUR_M1225_VehA))/dt); 
for i=1:length(samples) 
    for j=1:length(Templet_ITUR_M1225_VehA(1,:)) 
        if i*dt==Templet_ITUR_M1225_VehA(1,j) 
          samples(i)=10^(Templet_ITUR_M1225_VehA(2,j)/20); 
        end; 
    end; 
end; 
PlotSpectrum2(samples,width);