code so far

fs$=100$

dur $=13$

tt$13=$ np$.$arange$(0,$ dur, $1/$fs$)$

n$13=$ tt$13*$fs

b$13=$ np$.$array$([1,1,1,-1,-1,-1,1,-1,-1,1,-1,1,-1])$

rect$13=$np$.$ones$($len$($tt$13))$

# nreceptions $=$ number of receptions

# dur$\_$reception is the total duration $($seconds$)$ of the reception that the signal pulse is embedded in

# datfile $=$ name of data file to save receptions in

# First set sigpulse equal to whichever of the two pulses we want to use $($b$13$ or rect$13)$

# and define the output data filename accordingly

# Uncomment the correct lines to make either the Barker or rectangular pulses

#sigpulse$=$b$13$

#datfile $=$ 'proj$\_$part$2\_$data$\_$barker$13.$npz$\text{'}$

sigpulse$=$rect$13$

datfile $=$ 'proj$\_$part$2\_$data$\_$rect$13.$npz$\text{'}$

# The settings below make $10$ receptions all with an arrival time of $1\mathrm{.}5$ seconds

# Note that to increase or decrease the noise, change the noise$\_$std variable $($noise standard deviation$)$

nreceptions$=10$

starttime$=0$

arrivetime$=1\mathrm{.}5$

indarrive$=($np$.$rint$($arrivetime$*$fs$)).$astype$($int$)$

dur$\_$reception$=30$

noise$\_$std$=0\mathrm{.}05$

# Initialize receptions matrix with zeros

receptions$=$np$.$zeros$(($nreceptions$,$dur$\_$reception$*$fs$))$

# Create the signal piece $($sigpulse embedded in a zero vector$)$

sigrecept$=$np$.$zeros$($dur$\_$reception$*$fs$)$

sigrecept$[$indarrive:indarrive$+$len$($sigpulse$)]=$sigpulse

# For loop to generate receptions

for ind in range$(0,$nreceptions$)$:

# Create the Gaussian random noise to add to signal

noiserecept$=$noise$\_$std$*$np$.$random.standard$\_$normal$($len$($sigrecept$))$

# Store in reception in array

receptions$[$ind$,$:$]=$sigrecept$+$noiserecept

# Save the data

np$.$savez$($datfile$,$fs$=$fs$,$starttime$=$starttime,nreceptions$=$nreceptions,receptions$=$receptions$)$

# You can change the name of indatfile to reflect whatever data set you're using

indatfile$=$datfile

npz$\_$file $=$ np$.$load$($indatfile$)$

fs$=$npz$\_$file$[\text{'}$fs$\text{'}]$

starttime$=$npz$\_$file$[\text{'}$starttime$\text{'}]$

nreceptions$=$npz$\_$file$[\text{'}$nreceptions$\text{'}]$

receptions$=$npz$\_$file$[\text{'}$receptions$\text{'}]$

plt$.$plot$($receptions$[0,$:$])$

plt$.$plot$($receptions$[1,$:$])$

plt$.$show$()$