-
Notifications
You must be signed in to change notification settings - Fork 12
Expand file tree
/
Copy pathimport numpy as np.py
More file actions
127 lines (76 loc) · 2.47 KB
/
import numpy as np.py
File metadata and controls
127 lines (76 loc) · 2.47 KB
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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
import numpy as np
from scipy.io import wavfile
import random
import sys
fs = 43750
# t = 1.0 / (520 + (5/6))
# f = 2083.33333
samples = np.zeros(0)
def markBit():
f = 2083.33333
t = 1.0 / (520 + (5/6))
samples = np.arange(t * fs) / fs
roffle = np.sin(2 * np.pi * f * samples)
return roffle * 0.8
def spaceBit():
f = 1562.5
t = 1.0 / (520 + (5/6))
samples = np.arange(t * fs) / fs
return np.sin(2 * np.pi * f * samples)
signal = np.zeros(20000)
def byte(the_byte):
sys.stdout.write(the_byte)
sys.stdout.write(" ")
byte_data = np.zeros(0)
for i in range(0, 8):
if ord(the_byte) >> i & 1:
sys.stdout.write("1")
byte_data = np.append(byte_data, markBit())
else:
sys.stdout.write("0")
byte_data = np.append(byte_data, spaceBit())
sys.stdout.write("\n")
sys.stdout.flush()
return byte_data
def extramarks(numberOfMarks):
"""SAGE encoders seem to add a few mark bits at the beginning and end"""
byte_data = np.zeros(0)
for i in range(0, numberOfMarks):
byte_data = np.append(byte_data, markBit())
return byte_data
def preamble():
byte_data = np.zeros(0)
for i in range(0, 16):
byte_data = np.append(byte_data, markBit())
byte_data = np.append(byte_data, markBit())
byte_data = np.append(byte_data, spaceBit())
byte_data = np.append(byte_data, markBit())
byte_data = np.append(byte_data, spaceBit())
byte_data = np.append(byte_data, markBit())
byte_data = np.append(byte_data, spaceBit())
byte_data = np.append(byte_data, markBit())
return byte_data
# ZCZC-WXR-RWT-020103-020209-020091-020121-029047-029165-029095-029037+0030-1051700-KEAX/NWS
# code = "ZCZC-EAS-RMT-011000+0100-2141800-SCIENCE-"
code = "ZCZC-WXR-HUW-024031+0030-2142201-SCIENCE -"
# code = "ZCZC-PEP-EAT-000000+0400-2142300-SCIENCE -"
# code = "SUCK MY FUCKING BALLS YOU FUCKING COCKSUCKERS"
# control string
# code = "ZCZC-EAS-RMT-011000+0100-2142200-KMMS FM -"
for i in range(0, 3):
signal = np.append(signal, extramarks(10))
signal = np.append(signal, preamble())
for char in code:
signal = np.append(signal, byte(char))
signal = np.append(signal, extramarks(6))
signal = np.append(signal, np.zeros(43750))
for i in range(0, 3):
signal = np.append(signal, extramarks(10))
signal = np.append(signal, preamble())
for char in "NNNN":
signal = np.append(signal, byte(char))
signal = np.append(signal, extramarks(6))
signal = np.append(signal, np.zeros(43750))
signal *= -32767
signal = np.int16(signal)
wavfile.write(str("same.wav"), fs, signal)