1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
# Make a signal for testing purposes
import numpy as np
t = np.linspace(0, 1.0, num=int(1.0*44100), endpoint=False)
#freqs = [1000, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000]
#ampls = [500000/(f*f) for f in freqs]
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')
#print(wave64)
wave64.tofile("multi_44.1_f64_2ch_1.0s.raw")
wave32.tofile("multi_44.1_f32_2ch_1.0s.raw")
|
