plotdelay.m

来自「Jakes Model + Rayleigh Channel + Time Se」· M 代码 · 共 63 行

%%%%%%%%%%%%%    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);