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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
#!/usr/bin/python
import scipy, sys
import numpy as np
ch = int(sys.argv[1])
out = sys.argv[2]
FS = 48000
PILOT_DB = -30
PILOT_FREQ = 43
TEST0 = (500, 1, -10)
TEST1 = (500, 10, -15)
TEST2 = (43, 1.5, -6)
TEST3 = (1000, 1.5, -6)
def db(x):
return 10 ** (x / 20)
def silence(t):
return np.zeros(int(FS * t))
def sine(f, t, v):
space = np.linspace(0, t, int(FS * t), endpoint=False)
return np.sin(2 * np.pi * f * space) * db(v)
signal = np.concatenate((
silence(3),
sine(*TEST0),
silence(2),
sine(*TEST1),
silence(2),
sine(*TEST2),
silence(2),
sine(*TEST3),
silence(5)
))
space = np.linspace(0, len(signal) / FS, len(signal), endpoint=False)
signal += np.sin(2 * np.pi * PILOT_FREQ * space) * db(PILOT_DB)
signal = np.concatenate((silence(60), signal))
signal = np.repeat(signal, ch).reshape((-1, ch))
scipy.io.wavfile.write(out, FS, signal.astype("float32"))
|
