awgn_test.m

来自「convolution coding program using viterab」· M 代码 · 共 35 行

% setup
e = 0;
l = 50000;
snr = 0;
% Data generator
temp=rand(1,l);
for a= 1:l
if(temp(a)<0.5)
data(a)=0;
else
data(a)=1;
end
end
% bpsk modulation
for b = 1 : l
if data(b) == 0
temp1(b) = -1;
else temp1(b) = 1;
end
end
% AWGN channel
temp2 = awgn(temp1, snr);
for c= 1 : l
if temp2(c) <= 0
data_dec(c) = 0 ;
else data_dec(c) = 1;
end
end
% BER caculation
for m= 1:l
if (data(m) ~=data_dec(m)),
e = e + 1;
end
end
BER = e/l;