test_reed_solomon.m

source Matlab traites the UWB

%
%
clear all
%
%

% GF(2^6)
m = 6;
% Number of bits to send. Note: according to 15.4a, the maximum number of bit
% is 330;
DATALENGTH = 876; % Corresponds to 1020 encoded bit
fullblock = floor(DATALENGTH/330);
PACKETLENGTH = fullblock*378 + (DATALENGTH-fullblock*330) + 48;
fprintf('[RS test] DATALENGTH   = %d\n',DATALENGTH);
fprintf('[RS test] PACKETLENGTH = %d\n',PACKETLENGTH);
data_bit = round(rand(1,DATALENGTH));

% Seed for the interleaver
interleaver_seed = round(2^20*rand);

% Encode
tx_bit = rs_encoder_154a(data_bit,DATALENGTH,interleaver_seed);
fprintf('[RS test] length of tx_data = %d\n',length(tx_bit));

% Add noise, do whatever necessary
rx_data = tx_bit;

% Decode
dec_data = rs_decoder_154a(rx_data,interleaver_seed);

% Compute the bit-error rate
num_err = sum(abs(data_bit-dec_data));
ber_sim = num_err/DATALENGTH;

fprintf('BER = %.12f (%d errors)\n',ber_sim,num_err);