# process the signal
output = mics.process()

# save to output file
inp = pra.normalize(pra.highpass(mics.signals[mics.M / 2], Fs))
out = pra.normalize(pra.highpass(output, Fs))

wavfile.write('output_samples/input.wav', Fs, inp)
wavfile.write('output_samples/output.wav', Fs, out)
'''
Plot Stuff
'''
f_size = (3.93, 1.57)

# plot the room and beamformer
room1.plot(img_order=np.minimum(room1.max_order, 1), freq=freq)

# plot the beamforming weights
plt.figure()
mics.plot(FD=False)

# plot before/after processing
plt.figure()
pra.comparePlot(inp, out, Fs)

# plot angle/frequency plot
plt.figure()
mics.plot_beam_response()

plt.show()
# process the signal
output = mics.process()

# save to output file
inp = pra.normalize(pra.highpass(mics.signals[mics.M/2], Fs))
out = pra.normalize(pra.highpass(output, Fs))

wavfile.write('output_samples/input.wav', Fs, inp)
wavfile.write('output_samples/output.wav', Fs, out)

'''
Plot Stuff
'''
# plot the room and beamformer
room1.plot(img_order=np.minimum(room1.max_order, 1), 
        freq=freq)

# plot the beamforming weights
plt.figure()
mics.plot(FD=True)

# plot before/after processing
plt.figure()
pra.comparePlot(inp, out, Fs)

# plot angle/frequency plot
plt.figure()
mics.plot_beam_response()

plt.show()