#!/usr/bin/python
import numpy as np
import scipy.signal
import sys
from struct import pack
N=800
 # number of pfb channels
P=8
 # number of filter taps per channel

if len(sys.argv) < 2:
    print "Please supply an output filename."
    sys.exit(1)
h0 = scipy.signal.firwin(P * N, 1. / N)
p = N * np.flipud(h0.reshape(P, N))
s = "".join([pack('f',x) for x in p.flatten()])
f = open(sys.argv[1],"w")
f.write(s)
f.close()
print "Channels: %i, Taps: %i, Written to: %s\n" % (N, P, sys.argv[1])