jfilter.m

来自「JAKES 模型的具体实现,应用瑞利衰落信道的特性」· M 代码 · 共 14 行

function[hnw]=Jfilter(Fd)
% generate the time domain response according the Jakes Spectrum Model;
% fd is the maximum Doppler frequency;

NU=0.25;
Xt=20/(2*pi*Fd);
C=1.457;                                    %% normolizing factor;
t=-Xt:(2*Xt/63):Xt;
x=2*pi*Fd*abs(t);
ht=C*Fd*BESSELJ(NU,x)./sqrt(sqrt(x));       %% Bessel function;

n=0:63;                                     %% 64 taps Jakes' filter;
hn=ht;
hnw=hn.*hamming(64)';                       %% hamming windowed;