def get_results(fade_in): excitation_windowed = excitation * windows.window_kaiser( N, fade_in, fade_out, fs, beta) excitation_windowed_zeropadded = generation.zero_padding( excitation_windowed, pad, fs) system_response = system(excitation_windowed_zeropadded) ir = calculation.deconv_process(excitation_windowed_zeropadded, system_response, fs) return ir
import numpy as np # Parameters of the measuring system fs = 44100 fstart = 1 fstop = 22050 duration = 1 # Generate excitation signal excitation = generation.log_sweep(fstart, fstop, duration, fs) N = len(excitation) excitation_windowed_beta2 = excitation * windows.window_kaiser( N, 500, 500, fs, beta=2) excitation_windowed_beta7 = excitation * windows.window_kaiser( N, 500, 500, fs, beta=7) # Windowed Sweep-Signal plotting.plot_time(excitation, fs) plt.savefig('log_excitation.png') plt.close() plotting.plot_time(excitation_windowed_beta7, fs) plt.savefig('log_excitation_windowed') plt.close() # AKF