import numpy as np
t = np.linspace(0, 1.0, num=int(1.0*44100), endpoint=False)
freqs = list(range(1000, 21000, 1000))
ampls = [1.0/len(freqs)]*len(freqs)
wave = np.zeros(len(t))
for f, a in zip(freqs, ampls):
wave = wave + a * np.sin(f*2*np.pi*t)
wave= np.reshape(wave,(-1,1))
wave = np.concatenate((wave, wave), axis=1)
wave64 = wave.astype('float64')
wave32 = wave.astype('float32')
wave64.tofile("multi_44.1_f64_2ch_1.0s.raw")
wave32.tofile("multi_44.1_f32_2ch_1.0s.raw")