%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?All?rights?reserved?by?Krishna?Pillai?http://www.dsplog.com
%?The?file?may?not?be?re-distributed?without?explicit?authorization
%?from?Krishna?Pillai.
%?Checked?for?proper?operation?with?Octave?Version?3.0.0
%?Author????????:?Krishna?Pillai
%?Email?????????:?krishna@dsplog.com
%?Version???????:?1.0
%?Date??????????:?8?August?2008
%?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%?script?for?computing?the?BER?for?BPSK?modulation?in?a
%?Rayleigh?fading?channel

clear
N?=?10^6?%?number?of?bits?or?symbols

%?Transmitter
ip?=?rand（1N）>0.5;?%?generating?01?with?equal?probability
s?=?2*ip-1;?%?BPSK?modulation?0?->?-1;?1?->?0?

?
Eb_N0_dB?=?[-3:35];?%?multiple?Eb/N0?values

for?ii?=?1:length（Eb_N0_dB）
???
???n?=?1/sqrt（2）*[randn（1N）?+?j*randn（1N）];?%?white?gaussian?noise?0dB?variance?
???h?=?1/sqrt（2）*[randn（1N）?+?j*randn（1N）];?%?Rayleigh?channel
???
???%?Channel?and?noise?Noise?addition
???y?=?h.*s?+?10^（-Eb_N0_dB（ii）/20）*n;?

???%?equalization
???yHat?=?y./h;

???%?receiver?-?hard?decision?decoding
???ipHat?=?real（yHat）>0;

???%?counting?the?errors
???nErr（ii）?=?size（find（[ip-?ipHat]）2）;

end

simBer?=?nErr/N;?%?simulated?ber
theoryBerAWGN?=?0.5*erfc（sqrt（10.^（Eb_N0_dB/10）））;?%?theoretical?ber
EbN0Lin?=?10.^（Eb_N0_dB/10）;
theoryBer?=?0.5.*（1-sqrt（EbN0Lin./（EbN0Lin+1）））;

%?plot
close?all
figure
semilogy（Eb_N0_dBtheoryBerAWGN‘cd-‘‘LineWidth‘2）;
hold?on
semilogy（Eb_N0_dBtheoryBer‘bp-‘‘LineWidth‘2）;
semilogy（Eb_N0_dBsimBer‘mx-‘‘LineWidth‘2）;
axis（[-3?35?10^-5?0.5]）
grid?on
legend（‘AWGN-Theory‘‘Rayleigh-Theory‘?‘Rayleigh-Simulation‘）;
xlabel（‘Eb/No?dB‘）;
ylabel（‘Bit?Error?Rate‘）;
title（‘BER?for?BPSK?modulation?in?Rayleigh?channel‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1813??2010-04-16?03:23??BER?for?BPSK?in?Rayleigh?channel.m

-----------?---------??----------?-----??----

?????????????????1813????????????????????1