Exemple #1
0
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
Exemple #2
0
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