%--------------------------------------------------------------------------
%settings

length_source = 1000;
pulse_width = 1e-9;
sample_time = 1e-10;
duty_ratio = 0.25;
ds_code = [-1;-1;1;1;1;1;1;1;-1;1;1;1;-1;1;-1;]';
SNR = 10; %dB
delay_sample = 10;

chan = 1;
ch_idx = 1;

%--------------------------------------------------------------------------
%transmitter

bit_source = randsrc(1,length_source);% generate"1" and "-1"

pulse_sample = pulse(pulse_width,sample_time,duty_ratio);% pulse sample
length_pulse_sample = length(pulse_sample);

bit_sample = ds_modulation(ds_code,pulse_sample);% bit sample
length_bit_sample = length(bit_sample);

tx_sample = tx(bit_source,bit_sample);% tx_sample
length_tx_sample = length(tx_sample);

%--------------------------------------------------------------------------
%channel

channel = uwb_imr(chan, ch_idx,sample_time);% channel

rx_without_noise_sample = conv(tx_sample,channel);
bit_template_sample = rx_without_noise_sample(1:length_bit_sample)'; 

rx_sample = awgn(rx_without_noise_sample,SNR);
rx_sample = [zeros(delay_sample,1);rx_sample(1:length_tx_sample-delay_sample)];
 
%load data_ds_uwb;

%--------------------------------------------------------------------------
%receiver

rx_sample = reshape(rx_sample,length_bit_sample,length_source);
bit_out = (bit_template_sample*rx_sample)>0*2-1;
BER = sum(bit_out~=bit_source)/length_source;