diff --git a/tools/MATLAB_scripts/plot_Ascan.m b/tools/MATLAB_scripts/plot_Ascan.m index da332fb2..e4540be3 100644 --- a/tools/MATLAB_scripts/plot_Ascan.m +++ b/tools/MATLAB_scripts/plot_Ascan.m @@ -20,7 +20,7 @@ if filename ~= 0 header.nrx = h5readatt(fullfilename, '/', 'nrx'); % Time vector for plotting - time = linspace(0, (header.iterations)*(header.dt), header.iterations)'; + time = linspace(0, (header.iterations - 1) * header.dt, header.iterations)'; % Initialise structure for field arrays fields.ex = zeros(header.iterations, header.nrx); diff --git a/tools/MATLAB_scripts/plot_Bscan.m b/tools/MATLAB_scripts/plot_Bscan.m index b570d383..b863c875 100644 --- a/tools/MATLAB_scripts/plot_Bscan.m +++ b/tools/MATLAB_scripts/plot_Bscan.m @@ -17,8 +17,8 @@ if filename ~= 0 field = input(prompt,'s'); fieldpath = strcat('/rxs/rx1/', field); field = h5read(filename, fieldpath)'; - time = linspace(0, iterations*dt, iterations)'; - traces = 0:size(field,2); + time = linspace(0, (iterations - 1) * dt, iterations)'; + traces = 0:size(field, 2); fh1=figure('Name', filename); clims = [-max(max(abs(field))) max(max(abs(field)))]; @@ -50,4 +50,4 @@ if filename ~= 0 set(fh1,'PaperPosition', [xMargin yMargin xSize ySize]) set(fh1,'PaperOrientation', 'portrait') -end \ No newline at end of file +end diff --git a/tools/plot_Ascan.py b/tools/plot_Ascan.py index badad769..bee92a4d 100644 --- a/tools/plot_Ascan.py +++ b/tools/plot_Ascan.py @@ -47,8 +47,9 @@ def mpl_plot(filename, outputs=Rx.defaultoutputs, fft=False): nrx = f.attrs['nrx'] dt = f.attrs['dt'] iterations = f.attrs['Iterations'] - time = np.linspace(0, 1, iterations) - time *= (iterations * dt) + time = np.linspace(0, (iterations - 1) * dt, num=iterations) + print(time[0], time[1], time[-1]) + print(time.size) # Check there are any receivers if nrx == 0: diff --git a/tools/plot_antenna_params.py b/tools/plot_antenna_params.py index 3751d09c..3cfc84bb 100644 --- a/tools/plot_antenna_params.py +++ b/tools/plot_antenna_params.py @@ -49,8 +49,7 @@ def calculate_antenna_params(filename, tltxnumber=1, tlrxnumber=None, rxnumber=N iterations = f.attrs['Iterations'] # Calculate time array and frequency bin spacing - time = np.linspace(0, 1, iterations) - time *= (iterations * dt) + time = np.linspace(0, (iterations - 1) * dt, num=iterations) df = 1 / np.amax(time) print('Time window: {:g} s ({} iterations)'.format(np.amax(time), iterations)) diff --git a/tools/plot_source_wave.py b/tools/plot_source_wave.py index ad7eb0b9..918dc897 100644 --- a/tools/plot_source_wave.py +++ b/tools/plot_source_wave.py @@ -73,8 +73,7 @@ def mpl_plot(w, timewindow, dt, iterations, fft=False): plt (object): matplotlib plot object. """ - time = np.linspace(0, 1, iterations) - time *= (iterations * dt) + time = np.linspace(0, (iterations - 1) * dt, num=iterations) waveform = np.zeros(len(time)) timeiter = np.nditer(time, flags=['c_index'])