lcr_simch.m

来自「无线移动信道的matlab源码 用于对无线通信信道的编程实现」· M 代码 · 共 40 行

%--------------------------------------------------------------------
% lcr_sim.m ---------------------------------------------------------
%
% Program for the computation of the level-crossing rate N(r).
%
%--------------------------------------------------------------------
% N_r=lcr_sim(xi_t,r,T_sim,PLOT)
%--------------------------------------------------------------------
% Explanation of the input parameters:
%
% xi_t: deterministic process or time-domain signal to be analysed 
%       with respect to the level-crossing rate N(r)
% r: level vector
% T_sim: duration of the simulation
% PLOT: plot of the resulting level-crossing rate N(r), if PLOT==1

function N_r=lcr_simch(sigma_1_2,BEITA1,r,f_max,PLOT)

if nargin==3,
   PLOT=0;
end

N_r=zeros(size(r));

for k=1:length(r),
    a=sqrt(BEITA1/(2*pi));
    rk=10^(r(k)/20);
    b=rk/sigma_1_2;
    c=-rk*rk/(2*sigma_1_2);
    N_r(k)=a*b*exp(c);
end

N_r=N_r/f_max;

if PLOT==1,
   plot(r,N_r)
   xlabel('r')
   ylabel('N(r)')
end