新建文件夹 (2).rar_Jakes仿真_Jakes模型_jakes

function [h, t_state] = Jakes_Flat(fd, Ts, Ns, t_state, E0, phi_N) % INPUTS % fd : Doppler frequency % Ts : Sampling period % Ns : Sampling numbers % t0 : Initial time % E0 : Channel power % phi_N : Initial phase of the sin signal with Max Doppler frequency % OUTPUTS % h : Complex fading vector % t_state : current occasion if nargin < 6, phi_N = 0; end if nargin < 5, E0 = 1; end if nargin < 4, t_state = 0; end if nargin < 3 error('More arguments are needed for Jakes_Flat()'); end N0 = 8; N = 4*N0 + 2; wd = 2*pi*fd; % Max Doppler frequency[rad] t = t_state + [0:Ns-1] * Ts; % Time vector t_state = t(end) + Ts; % Final time coswt = [sqrt(2) * cos(wd*t); 2*cos(wd*cos(2*pi/N*[1:N0]')*t)]; h = E0/sqrt(2*N0 + 1) * exp(1i*[phi_N pi/(N0+1)*[1:N0]]) * coswt; help butter