berrayleighdoppler.m

来自「这是一个关于完成QAM调制的Matlab示例程序」· M 代码 · 共 20 行

function [ber] = BerRayleighDoppler(gammaB,fDT)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% This is just a utility function to compute the expected BER
% of a signal in Rayleigh fading with doppler error in the receiver
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

MIN_FDT = 1/10000;  % Not sure if this is a reasonable value.  Using for now.

if (fDT < MIN_FDT)
    % Doppler is essentially zero.  Just use standard Rayleigh fading equation
    ber = 0.5*(1-sqrt(gammaB/(1+gammaB)));
else
    ber = 0.5 + (1/(8*pi*fDT))*(asin(sqrt(2*gammaB)*cos((pi/4)+2*pi*fDT)/sqrt(1+2*gammaB)) - ...
        asin(sqrt(2*gammaB)*cos((pi/4)-2*pi*fDT)/sqrt(1+2*gammaB)));
end

return