main.m

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

%%%%%%%%%%%%%    Main      %%%%%%%%%%%%%%%
%%% CarrierFrq :  Carrier Frequence  2GHz
%%% VehicSpd   :  Vehicular Speed  3 / 120(Km/h)  =  0.8333  / 33.3333(m/s)
%%% C          :  Speed of the light  3e8 m/s
%%% fm         :  Max Doppler shift
%%% Tc         :  Correlation Time  1/fm
%%% T          :  Sampling interval
%%% No         :  Sampling Point
%%% P          :  Factor of Correlation Time/Sampling Interval(if want to observe the time selective fading, observe time must longer than Correlation Time, so, total number of sampling points must larger than P)
%%% Model      :  Jakes Model outcome
%%% Envelop    :  Jakes Model Envelop
%%% Phase      :  Jakes Model Phase
%%% PDF        :  Probability Density Functions of the model
%%% Time       :  Simulation time

clear;
clc;
%%%%%%%%%%%%%%%%%%%%   Define Parameters   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
CarrierFrq=2e9;
VehicSpd=33.333333333333333;
C=3e8;
fm=VehicSpd*CarrierFrq/C;
Tc=1/fm;              %Correlation Time
P=10;           
T=Tc/P;               %     Sampling Interval Time << Correlation Time
%T=4.5000e-006;
discard=P*20;         %%First 10 Correlation Time can not be used;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[Model,Time]=JakesModel(P,fm,T);
No=length(Model);
Envelop=abs(Model);
Phase=angle(Model)/pi;
[PDF_E,x_axis_E]=hist(Envelop,20);
PDF_E=PDF_E/No;
[PDF_P,x_axis_P]=hist(Phase,100);
PDF_P=PDF_P/No;

%%%%%%%%%%%%%%%%%%%%%%%%%          PLOT             %%%%%%%%%%%%%%
figure(1);
PlotEnvelop(Envelop,Time);
figure(2);
PlotPdf(x_axis_E,PDF_E);
figure(3);
PlotPhase(Phase,Time);
figure(4);
PlotPdf(x_axis_P,PDF_P);

%%%%%%%%%%%%%%   Autocorr   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Auto_corr=Autocorr(Model);
figure(5);
PlotAutocorr(Auto_corr,Time);

%%%%%%%%%%%   Plot U_Ship Spectrum   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
PlotSpectrum(real(Model),Time(length(Time)));

%%%%%%%%%%%   Plot  Delay   %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
PlotDelay;