def test_006(self):
''' Test the complex AGC loop (attack and decay rate inputs) '''
tb = self.tb
sampling_freq = 100
N = int(5 * sampling_freq)
src1 = analog.sig_source_c(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100)
dst1 = blocks.vector_sink_c()
head = blocks.head(gr.sizeof_gr_complex, N)
ref = 1
agc = analog.agc3_cc(1e-2, 1e-3, ref)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
M = 100
result = map(lambda x: abs(x), dst_data[N - M:])
self.assertFloatTuplesAlmostEqual(result, M * [
ref,
], 4)
def test_000(self):
N = 1000 # number of samples to use
fs = 1000 # baseband sampling rate
freq = 100
signal = analog.sig_source_c(fs, analog.GR_SIN_WAVE, freq, 1)
head = blocks.head(gr.sizeof_gr_complex, N)
op = channels.channel_model(0.0, 0.0, 1.0, [
1,
], 0)
snk = blocks.vector_sink_c()
snk1 = blocks.vector_sink_c()
op.set_noise_voltage(0.0)
op.set_frequency_offset(0.0)
op.set_taps([
1,
])
op.set_timing_offset(1.0)
self.tb.connect(signal, head, op, snk)
self.tb.connect(op, snk1)
self.tb.run()
dst_data = snk.data()
exp_data = snk1.data()
self.assertComplexTuplesAlmostEqual(exp_data, dst_data, 5)
def test_001(self):
# Just running some data through null source/sink
src = blocks.null_source(gr.sizeof_float)
hed = blocks.head(gr.sizeof_float, 100)
dst = blocks.null_sink(gr.sizeof_float)
self.tb.connect(src, hed, dst)
self.tb.run()
def test_tags(self):
src = tag_source()
sink = tag_sink()
head = blocks.head(gr.sizeof_float, 50000) #should be enough items to get a tag through
tb = gr.top_block()
tb.connect(src, head, sink)
tb.run()
self.assertEqual(sink.key, "example_key")
def test_001(self):
# Just running some data through null source/sink
src = blocks.null_source(gr.sizeof_float)
hed = blocks.head(gr.sizeof_float, 100)
dst = blocks.null_sink(gr.sizeof_float)
self.tb.connect(src, hed, dst)
self.tb.run()
def test_tags(self):
src = tag_source()
sink = tag_sink()
head = blocks.head(gr.sizeof_float,
50000) #should be enough items to get a tag through
tb = gr.top_block()
tb.connect(src, head, sink)
tb.run()
self.assertEqual(sink.key, "example_key")
def test_head(self):
src_data = (1, 2, 3, 4, 5, 6, 7, 8, 9, 10)
expected_result = (1, 2, 3, 4)
src1 = blocks.vector_source_i(src_data)
op = blocks.head(gr.sizeof_int, 4)
dst1 = blocks.vector_sink_i()
self.tb.connect(src1, op)
self.tb.connect(op, dst1)
self.tb.run()
dst_data = dst1.data()
self.assertEqual(expected_result, dst_data)
def test_head(self):
src_data = (1, 2, 3, 4, 5, 6, 7, 8, 9, 10)
expected_result = (1, 2, 3, 4)
src1 = blocks.vector_source_i(src_data)
op = blocks.head(gr.sizeof_int, 4)
dst1 = blocks.vector_sink_i()
self.tb.connect(src1, op)
self.tb.connect(op, dst1)
self.tb.run()
dst_data = dst1.data()
self.assertEqual(expected_result, dst_data)
def test_sqr_f(self):
tb = self.tb
expected_result = (0, 0, 0, 0, 1, 1, 1, 1, 0)
src1 = analog.sig_source_f(8, analog.GR_SQR_WAVE, 1.0, 1.0)
op = blocks.head(gr.sizeof_float, 9)
dst1 = blocks.vector_sink_f()
tb.connect(src1, op)
tb.connect(op, dst1)
tb.run()
dst_data = dst1.data()
self.assertEqual(expected_result, dst_data)
def test_saw_f(self):
tb = self.tb
expected_result = (.5, .625, .75, .875, 0, .125, .25, .375, .5)
src1 = analog.sig_source_f(8, analog.GR_SAW_WAVE, 1.0, 1.0)
op = blocks.head(gr.sizeof_float, 9)
dst1 = blocks.vector_sink_f()
tb.connect(src1, op)
tb.connect(op, dst1)
tb.run()
dst_data = dst1.data()
self.assertFloatTuplesAlmostEqual(expected_result, dst_data, 5)
def test_const_i(self):
tb = self.tb
expected_result = (1, 1, 1, 1)
src1 = analog.sig_source_i(1e6, analog.GR_CONST_WAVE, 0, 1)
op = blocks.head(gr.sizeof_int, 4)
dst1 = blocks.vector_sink_i()
tb.connect(src1, op)
tb.connect(op, dst1)
tb.run()
dst_data = dst1.data()
self.assertEqual(expected_result, dst_data)
def test_sqr_c(self):
tb = self.tb #arg6 is a bit before -PI/2
expected_result = (1j, 1j, 0, 0, 1, 1, 1+0j, 1+1j, 1j)
src1 = analog.sig_source_c(8, analog.GR_SQR_WAVE, 1.0, 1.0)
op = blocks.head(gr.sizeof_gr_complex, 9)
dst1 = blocks.vector_sink_c()
tb.connect(src1, op)
tb.connect(op, dst1)
tb.run()
dst_data = dst1.data()
self.assertEqual(expected_result, dst_data)
def test_001_correlate(self):
degree = 10
length = 2**degree - 1
src = digital.glfsr_source_f(degree)
head = blocks.head(gr.sizeof_float, length * length)
f2c = blocks.float_to_complex()
corr = digital.pn_correlator_cc(degree)
dst = blocks.vector_sink_c()
self.tb.connect(src, head, f2c, corr, dst)
self.tb.run()
data = dst.data()
self.assertEqual(data[-1], (1.0 + 0j))
def test_001_correlate(self):
degree = 10
length = 2**degree-1
src = digital.glfsr_source_f(degree)
head = blocks.head(gr.sizeof_float, length*length)
f2c = blocks.float_to_complex()
corr = digital.pn_correlator_cc(degree)
dst = blocks.vector_sink_c()
self.tb.connect(src, head, f2c, corr, dst)
self.tb.run()
data = dst.data()
self.assertEqual(data[-1], (1.0+0j))
def test_cosine_f(self):
tb = self.tb
sqrt2 = math.sqrt(2) / 2
expected_result = (1, sqrt2, 0, -sqrt2, -1, -sqrt2, 0, sqrt2, 1)
src1 = analog.sig_source_f(8, analog.GR_COS_WAVE, 1.0, 1.0)
op = blocks.head(gr.sizeof_float, 9)
dst1 = blocks.vector_sink_f()
tb.connect(src1, op)
tb.connect(op, dst1)
tb.run()
dst_data = dst1.data()
self.assertFloatTuplesAlmostEqual(expected_result, dst_data, 5)
def test_tri_c(self):
tb = self.tb
expected_result = (1+.5j, .75+.75j, .5+1j, .25+.75j, 0+.5j,
.25+.25j, .5+0j, .75+.25j, 1+.5j)
src1 = analog.sig_source_c(8, analog.GR_TRI_WAVE, 1.0, 1.0)
op = blocks.head(gr.sizeof_gr_complex, 9)
dst1 = blocks.vector_sink_c()
tb.connect(src1, op)
tb.connect(op, dst1)
tb.run()
dst_data = dst1.data()
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 5)
def test_000(self):
N = 1000 # number of samples to use
M = 5 # Number of channels
fs = 1000 # baseband sampling rate
ifs = M * fs # input samp rate to decimator
channel = 0 # Extract channel 0
taps = filter.firdes.low_pass_2(
1,
ifs,
fs / 2,
fs / 10,
attenuation_dB=80,
window=filter.firdes.WIN_BLACKMAN_hARRIS)
signals = list()
add = blocks.add_cc()
freqs = [-200, -100, 0, 100, 200]
for i in xrange(len(freqs)):
f = freqs[i] + (M / 2 - M + i + 1) * fs
data = sig_source_c(ifs, f, 1, N)
signals.append(blocks.vector_source_c(data))
self.tb.connect(signals[i], (add, i))
head = blocks.head(gr.sizeof_gr_complex, N)
s2ss = blocks.stream_to_streams(gr.sizeof_gr_complex, M)
pfb = filter.pfb_decimator_ccf(M, taps, channel)
snk = blocks.vector_sink_c()
self.tb.connect(add, head, s2ss)
for i in xrange(M):
self.tb.connect((s2ss, i), (pfb, i))
self.tb.connect(pfb, snk)
self.tb.run()
Ntest = 50
L = len(snk.data())
t = map(lambda x: float(x) / fs, xrange(L))
# Create known data as complex sinusoids for the baseband freq
# of the extracted channel is due to decimator output order.
phase = 0
expected_data = map(lambda x: math.cos(2.*math.pi*freqs[2]*x+phase) + \
1j*math.sin(2.*math.pi*freqs[2]*x+phase), t)
dst_data = snk.data()
self.assertComplexTuplesAlmostEqual(expected_data[-Ntest:],
dst_data[-Ntest:], 4)
def est_fmdet_cf_002(self):
N = 100
src = analog.sig_source_c(1, analog.GR_SIN_WAVE, 0.2, 1)
head = blocks.head(gr.sizeof_gr_complex, N)
op = analog.fmdet_cf(1, 0.1, 0.3, 0.1)
dst = blocks.vector_sink_f()
self.tb.connect(src, head, op)
self.tb.connect(op, dst)
self.tb.run()
result_data = dst.data()[4:N]
expected_result = (100-4)*[-0.21755,]
self.assertFloatTuplesAlmostEqual(expected_result, result_data, 4)
def test_000(self):
N = 1000 # number of samples to use
fs = 1000 # baseband sampling rate
freq = 100
fDTs = 0.01
K = 4
signal = analog.sig_source_c(fs, analog.GR_SIN_WAVE, freq, 1)
head = blocks.head(gr.sizeof_gr_complex, N)
op = channels.fading_model(8, fDTs=fDTs, LOS=True, K=K, seed=0)
snk = blocks.vector_sink_c()
snk1 = blocks.vector_sink_c()
self.assertAlmostEqual(K, op.K(), 4)
self.assertAlmostEqual(fDTs, op.fDTs(), 4)
def test_000(self):
N = 1000 # number of samples to use
fs = 1000 # baseband sampling rate
freq = 100
fDTs = 0.01
K = 4
signal = analog.sig_source_c(fs, analog.GR_SIN_WAVE, freq, 1)
head = blocks.head(gr.sizeof_gr_complex, N)
op = channels.fading_model(8, fDTs=fDTs, LOS=True,
K=K, seed=0)
snk = blocks.vector_sink_c()
snk1 = blocks.vector_sink_c()
self.assertAlmostEqual(K, op.K(), 4)
self.assertAlmostEqual(fDTs, op.fDTs(), 4)
def test_000(self):
dst_data = [0, 0, 1, 2, 0, 2, 1, 0]
src0 = blocks.vector_source_f(200 * [0])
src1 = blocks.vector_source_f(200 * [1])
src2 = blocks.vector_source_f(200 * [2])
itg = blocks.patterned_interleaver(gr.sizeof_float, dst_data)
dst = blocks.vector_sink_f()
head = blocks.head(gr.sizeof_float, 8)
self.tb.connect(src0, (itg, 0))
self.tb.connect(src1, (itg, 1))
self.tb.connect(src2, (itg, 2))
self.tb.connect(itg, head, dst)
self.tb.run()
self.assertEqual(list(dst_data), list(dst.data()))
def test_000(self):
dst_data = [0, 0, 1, 2, 0, 2, 1, 0]
src0 = blocks.vector_source_f(200 * [0])
src1 = blocks.vector_source_f(200 * [1])
src2 = blocks.vector_source_f(200 * [2])
itg = blocks.patterned_interleaver(gr.sizeof_float, dst_data)
dst = blocks.vector_sink_f()
head = blocks.head(gr.sizeof_float, 8)
self.tb.connect(src0, (itg, 0))
self.tb.connect(src1, (itg, 1))
self.tb.connect(src2, (itg, 2))
self.tb.connect(itg, head, dst)
self.tb.run()
self.assertEqual(list(dst_data), list(dst.data()))
def test_ff(self):
amp = 2
src_data = sig_source_f(1, 0.01, amp, 200)
N = 750000
expected_data = amp / math.sqrt(2.0)
src = blocks.vector_source_f(src_data, True)
head = blocks.head(gr.sizeof_float, N)
op = blocks.rms_ff(0.0001)
dst = blocks.vector_sink_f()
self.tb.connect(src, head, op, dst)
self.tb.run()
dst_data = dst.data()
self.assertAlmostEqual(dst_data[-1], expected_data, 4)
def test_cf(self):
amp = 4
src_data = sig_source_c(1, 0.01, amp, 200)
N = 750000
expected_data = amp
src = blocks.vector_source_c(src_data, True)
head = blocks.head(gr.sizeof_gr_complex, N)
op = blocks.rms_cf(0.0001)
dst = blocks.vector_sink_f()
self.tb.connect(src, head, op, dst)
self.tb.run()
dst_data = dst.data()
self.assertAlmostEqual(dst_data[-1], expected_data, 4)
def test_cf(self):
amp = 4
src_data = sig_source_c(1, 0.01, amp, 200)
N = 750000
expected_data = amp
src = blocks.vector_source_c(src_data, True)
head = blocks.head(gr.sizeof_gr_complex, N)
op = blocks.rms_cf(0.0001)
dst = blocks.vector_sink_f()
self.tb.connect(src, head, op, dst)
self.tb.run()
dst_data = dst.data()
self.assertAlmostEqual(dst_data[-1], expected_data, 4)
def test_ff(self):
amp = 2
src_data = sig_source_f(1, 0.01, amp, 200)
N = 750000
expected_data = amp / math.sqrt(2.0)
src = blocks.vector_source_f(src_data, True)
head = blocks.head(gr.sizeof_float, N)
op = blocks.rms_ff(0.0001)
dst = blocks.vector_sink_f()
self.tb.connect(src, head, op, dst)
self.tb.run()
dst_data = dst.data()
self.assertAlmostEqual(dst_data[-1], expected_data, 4)
def __init__(self, constellation, f, N0=0.25, seed=-666L):
"""
constellation - a constellation object used for modulation.
f - a finite state machine specification used for coding.
N0 - noise level
seed - random seed
"""
super(trellis_tb, self).__init__()
# packet size in bits (make it multiple of 16 so it can be packed in a short)
packet_size = 1024 * 16
# bits per FSM input symbol
bitspersymbol = int(round(math.log(f.I()) /
math.log(2))) # bits per FSM input symbol
# packet size in trellis steps
K = packet_size / bitspersymbol
# TX
src = blocks.lfsr_32k_source_s()
# packet size in shorts
src_head = blocks.head(gr.sizeof_short, packet_size / 16)
# unpack shorts to symbols compatible with the FSM input cardinality
s2fsmi = blocks.packed_to_unpacked_ss(bitspersymbol, gr.GR_MSB_FIRST)
# initial FSM state = 0
enc = trellis.encoder_ss(f, 0)
mod = digital.chunks_to_symbols_sc(constellation.points(), 1)
# CHANNEL
add = blocks.add_cc()
noise = analog.noise_source_c(analog.GR_GAUSSIAN, math.sqrt(N0 / 2),
seed)
# RX
# data preprocessing to generate metrics for Viterbi
metrics = trellis.constellation_metrics_cf(constellation.base(),
digital.TRELLIS_EUCLIDEAN)
# Put -1 if the Initial/Final states are not set.
va = trellis.viterbi_s(f, K, 0, -1)
# pack FSM input symbols to shorts
fsmi2s = blocks.unpacked_to_packed_ss(bitspersymbol, gr.GR_MSB_FIRST)
# check the output
self.dst = blocks.check_lfsr_32k_s()
self.connect(src, src_head, s2fsmi, enc, mod)
self.connect(mod, (add, 0))
self.connect(noise, (add, 1))
self.connect(add, metrics, va, fsmi2s, self.dst)
def test_000(self):
N = 1000 # number of samples to use
M = 5 # Number of channels
fs = 1000 # baseband sampling rate
ifs = M*fs # input samp rate to decimator
channel = 0 # Extract channel 0
taps = filter.firdes.low_pass_2(1, ifs, fs/2, fs/10,
attenuation_dB=80,
window=filter.firdes.WIN_BLACKMAN_hARRIS)
signals = list()
add = blocks.add_cc()
freqs = [-200, -100, 0, 100, 200]
for i in xrange(len(freqs)):
f = freqs[i] + (M/2-M+i+1)*fs
data = sig_source_c(ifs, f, 1, N)
signals.append(blocks.vector_source_c(data))
self.tb.connect(signals[i], (add,i))
head = blocks.head(gr.sizeof_gr_complex, N)
s2ss = blocks.stream_to_streams(gr.sizeof_gr_complex, M)
pfb = filter.pfb_decimator_ccf(M, taps, channel)
snk = blocks.vector_sink_c()
self.tb.connect(add, head, s2ss)
for i in xrange(M):
self.tb.connect((s2ss,i), (pfb,i))
self.tb.connect(pfb, snk)
self.tb.run()
Ntest = 50
L = len(snk.data())
t = map(lambda x: float(x)/fs, xrange(L))
# Create known data as complex sinusoids for the baseband freq
# of the extracted channel is due to decimator output order.
phase = 0
expected_data = map(lambda x: math.cos(2.*math.pi*freqs[2]*x+phase) + \
1j*math.sin(2.*math.pi*freqs[2]*x+phase), t)
dst_data = snk.data()
self.assertComplexTuplesAlmostEqual(expected_data[-Ntest:], dst_data[-Ntest:], 4)
def __init__(self, constellation, f, N0=0.25, seed=-666L):
"""
constellation - a constellation object used for modulation.
f - a finite state machine specification used for coding.
N0 - noise level
seed - random seed
"""
super(trellis_tb, self).__init__()
# packet size in bits (make it multiple of 16 so it can be packed in a short)
packet_size = 1024 * 16
# bits per FSM input symbol
bitspersymbol = int(round(math.log(f.I()) / math.log(2))) # bits per FSM input symbol
# packet size in trellis steps
K = packet_size / bitspersymbol
# TX
src = blocks.lfsr_32k_source_s()
# packet size in shorts
src_head = blocks.head(gr.sizeof_short, packet_size / 16)
# unpack shorts to symbols compatible with the FSM input cardinality
s2fsmi = blocks.packed_to_unpacked_ss(bitspersymbol, gr.GR_MSB_FIRST)
# initial FSM state = 0
enc = trellis.encoder_ss(f, 0)
mod = digital.chunks_to_symbols_sc(constellation.points(), 1)
# CHANNEL
add = blocks.add_cc()
noise = analog.noise_source_c(analog.GR_GAUSSIAN, math.sqrt(N0 / 2), seed)
# RX
# data preprocessing to generate metrics for Viterbi
metrics = trellis.constellation_metrics_cf(constellation.base(), digital.TRELLIS_EUCLIDEAN)
# Put -1 if the Initial/Final states are not set.
va = trellis.viterbi_s(f, K, 0, -1)
# pack FSM input symbols to shorts
fsmi2s = blocks.unpacked_to_packed_ss(bitspersymbol, gr.GR_MSB_FIRST)
# check the output
self.dst = blocks.check_lfsr_32k_s()
self.connect(src, src_head, s2fsmi, enc, mod)
self.connect(mod, (add, 0))
self.connect(noise, (add, 1))
self.connect(add, metrics, va, fsmi2s, self.dst)
def test_001(self):
src_data = [float(x) for x in range(16)]
expected_result = tuple(src_data)
period = 9177
offset = 0
src = blocks.null_source(1)
head = blocks.head(1, 10000000)
ins = blocks.vector_insert_b([1], period, offset)
dst = blocks.vector_sink_b()
self.tb.connect(src, head, ins, dst)
self.tb.run()
result_data = dst.data()
for i in range(10000):
if (i % period == offset):
self.assertEqual(1, result_data[i])
else:
self.assertEqual(0, result_data[i])
def test_001(self):
src_data = [float(x) for x in range(16)]
expected_result = tuple(src_data)
period = 9177;
offset = 0;
src = blocks.null_source(1)
head = blocks.head(1, 10000000);
ins = blocks.vector_insert_b([1], period, offset);
dst = blocks.vector_sink_b()
self.tb.connect(src, head, ins, dst)
self.tb.run()
result_data = dst.data()
for i in range(10000):
if(i%period == offset):
self.assertEqual(1, result_data[i])
else:
self.assertEqual(0, result_data[i])
def test_002(self):
''' Test the floating point AGC loop (single rate input) '''
tb = self.tb
expected_result = (
7.2191943445432116e-07, 58.837181091308594, 89.700050354003906,
81.264183044433594, 45.506141662597656, 4.269894304798072e-07,
-42.948936462402344, -65.50335693359375, -59.368724822998047,
-33.261005401611328, -4.683740257860336e-07, 31.423542022705078,
47.950984954833984, 43.485683441162109, 24.378345489501953,
5.7254135299444897e-07, -23.062990188598633, -35.218441009521484,
-31.964075088500977, -17.934831619262695, -5.0591745548445033e-07,
16.998210906982422, 25.982204437255859, 23.606258392333984,
13.260685920715332, 4.9936483037527069e-07, -12.59880542755127,
-19.28221321105957, -17.54347038269043, -9.8700437545776367,
-4.188150626305287e-07, 9.4074573516845703, 14.422011375427246,
13.145503044128418, 7.41046142578125, 3.8512698097292741e-07,
-7.0924453735351562, -10.896408081054688, -9.9552040100097656,
-5.6262712478637695, -3.1982864356905338e-07, 5.4131259918212891,
8.3389215469360352, 7.6409502029418945, 4.3320145606994629,
2.882407841298118e-07, -4.194943904876709, -6.4837145805358887,
-5.9621825218200684, -3.3931560516357422)
sampling_freq = 100
src1 = analog.sig_source_f(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100.0)
dst1 = blocks.vector_sink_f()
head = blocks.head(gr.sizeof_float, int(5 * sampling_freq * 0.10))
agc = analog.agc_ff(1e-3, 1, 1)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
self.assertFloatTuplesAlmostEqual(expected_result, dst_data, 4)
def test_000(self):
N = 1000 # number of samples to use
fs = 1000 # baseband sampling rate
freq = 100
signal = analog.sig_source_c(fs, analog.GR_SIN_WAVE, freq, 1)
head = blocks.head(gr.sizeof_gr_complex, N)
op = channels.channel_model(0.0, 0.0, 1.0, [1,], 0)
snk = blocks.vector_sink_c()
snk1 = blocks.vector_sink_c()
op.set_noise_voltage(0.0)
op.set_frequency_offset(0.0)
op.set_taps([1,])
op.set_timing_offset(1.0)
self.tb.connect(signal, head, op, snk)
self.tb.connect(op, snk1)
self.tb.run()
dst_data = snk.data()
exp_data = snk1.data()
self.assertComplexTuplesAlmostEqual(exp_data, dst_data, 5)
def test_006(self):
''' Test the complex AGC loop (attack and decay rate inputs) '''
tb = self.tb
sampling_freq = 100
N = int(5*sampling_freq)
src1 = analog.sig_source_c(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100)
dst1 = blocks.vector_sink_c()
head = blocks.head(gr.sizeof_gr_complex, N)
ref = 1
agc = analog.agc3_cc(1e-2, 1e-3, ref)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
M = 100
result = map(lambda x: abs(x), dst_data[N-M:])
self.assertFloatTuplesAlmostEqual(result, M*[ref,], 4)
def test_004(self):
''' Test the floating point AGC loop (attack and decay rate inputs) '''
tb = self.tb
expected_result = \
(7.2191943445432116e-07,
58.837181091308594,
40.194305419921875,
2.9183335304260254,
0.67606079578399658,
8.6260438791896377e-09,
-1.4542514085769653,
-1.9210131168365479,
-1.0450780391693115,
-0.61939650774002075,
-1.2590258613442984e-08,
1.4308931827545166,
1.9054338932037354,
1.0443156957626343,
0.61937344074249268,
2.0983527804219193e-08,
-1.4308838844299316,
-1.9054274559020996,
-1.0443152189254761,
-0.61937344074249268,
-2.5180233009791664e-08,
1.4308837652206421,
1.9054274559020996,
1.0443154573440552,
0.61937344074249268,
3.3573645197293445e-08,
-1.4308838844299316,
-1.9054274559020996,
-1.0443152189254761,
-0.61937350034713745,
-3.7770352179222755e-08,
1.4308837652206421,
1.9054274559020996,
1.0443154573440552,
0.61937350034713745,
4.6163762590367696e-08,
-1.4308838844299316,
-1.9054274559020996,
-1.0443153381347656,
-0.61937344074249268,
-5.0360466019583328e-08,
1.4308837652206421,
1.9054274559020996,
1.0443155765533447,
0.61937344074249268,
5.8753879983441948e-08,
-1.4308837652206421,
-1.9054274559020996,
-1.0443153381347656,
-0.61937344074249268)
sampling_freq = 100
src1 = analog.sig_source_f(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100)
dst1 = blocks.vector_sink_f()
head = blocks.head(gr.sizeof_float, int(5 * sampling_freq * 0.10))
agc = analog.agc2_ff(1e-2, 1e-3, 1, 1)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
self.assertFloatTuplesAlmostEqual(expected_result, dst_data, 4)
def test_005(self):
''' Test the complex AGC loop (attack and decay rate inputs) '''
tb = self.tb
expected_result = \
((100.000244140625+7.2191943445432116e-07j),
(0.80881959199905396+0.58764183521270752j),
(0.30894950032234192+0.95084899663925171j),
(-0.30895623564720154+0.95086973905563354j),
(-0.80887287855148315+0.58768033981323242j),
(-0.99984413385391235+5.850709250410091e-09j),
(-0.80889981985092163-0.58770018815994263j),
(-0.30897706747055054-0.95093393325805664j),
(0.30898112058639526-0.95094609260559082j),
(0.80893135070800781-0.58772283792495728j),
(0.99990922212600708-8.7766354184282136e-09j),
(0.80894720554351807+0.58773452043533325j),
(0.30899339914321899+0.95098406076431274j),
(-0.30899572372436523+0.95099133253097534j),
(-0.80896598100662231+0.58774799108505249j),
(-0.99994778633117676+1.4628290578855285e-08j),
(-0.80897533893585205-0.58775502443313599j),
(-0.30900305509567261-0.95101380348205566j),
(0.30900448560714722-0.95101797580718994j),
(0.80898630619049072-0.58776277303695679j),
(0.99997037649154663-1.7554345532744264e-08j),
(0.80899184942245483+0.58776694536209106j),
(0.30900871753692627+0.95103120803833008j),
(-0.30900952219963074+0.95103377103805542j),
(-0.8089984655380249+0.58777159452438354j),
(-0.99998390674591064+2.3406109050938539e-08j),
(-0.809001624584198-0.58777409791946411j),
(-0.30901208519935608-0.95104163885116577j),
(0.30901262164115906-0.95104306936264038j),
(0.80900543928146362-0.587776780128479j),
(0.99999171495437622-2.6332081404234486e-08j),
(0.80900734663009644+0.58777821063995361j),
(0.30901408195495605+0.95104765892028809j),
(-0.30901429057121277+0.95104855298995972j),
(-0.80900967121124268+0.58777981996536255j),
(-0.99999648332595825+3.2183805842578295e-08j),
(-0.80901080369949341-0.58778077363967896j),
(-0.30901527404785156-0.95105135440826416j),
(0.30901545286178589-0.95105189085006714j),
(0.80901217460632324-0.58778166770935059j),
(0.99999916553497314-3.5109700036173308e-08j),
(0.809012770652771+0.58778214454650879j),
(0.30901595950126648+0.9510534405708313j),
(-0.30901598930358887+0.95105385780334473j),
(-0.80901366472244263+0.58778274059295654j),
(-1.0000008344650269+4.0961388947380328e-08j),
(-0.8090139627456665-0.58778303861618042j),
(-0.30901634693145752-0.95105475187301636j),
(0.30901640653610229-0.95105493068695068j),
(0.80901449918746948-0.5877833366394043j))
sampling_freq = 100
src1 = analog.sig_source_c(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100)
dst1 = blocks.vector_sink_c()
head = blocks.head(gr.sizeof_gr_complex, int(5 * sampling_freq * 0.10))
agc = analog.agc2_cc(1e-2, 1e-3, 1, 1)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 4)
def test_001(self):
''' Test the complex AGC loop (single rate input) '''
tb = self.tb
expected_result = (
(100.000244140625 + 7.2191943445432116e-07j),
(72.892257690429688 + 52.959323883056641j), (25.089065551757812 +
77.216217041015625j),
(-22.611061096191406 + 69.589706420898438j), (-53.357715606689453 +
38.766635894775391j),
(-59.458671569824219 + 3.4792964243024471e-07j),
(-43.373462677001953 - 31.512666702270508j), (-14.94139289855957 -
45.984889984130859j),
(13.478158950805664 - 41.48150634765625j), (31.838506698608398 -
23.132022857666016j),
(35.519271850585938 - 3.1176801940091536e-07j),
(25.942903518676758 + 18.848621368408203j), (8.9492912292480469 +
27.5430908203125j),
(-8.0852642059326172 + 24.883890151977539j), (-19.131628036499023 +
13.899936676025391j),
(-21.383295059204102 + 3.1281737733479531e-07j),
(-15.650330543518066 - 11.370632171630859j), (-5.4110145568847656 -
16.65339469909668j),
(4.9008159637451172 - 15.083160400390625j), (11.628337860107422 -
8.4484796524047852j),
(13.036135673522949 - 2.288476110834381e-07j),
(9.5726661682128906 + 6.954948902130127j), (3.3216962814331055 +
10.223132133483887j),
(-3.0204284191131592 + 9.2959251403808594j), (-7.1977195739746094 +
5.2294478416442871j),
(-8.1072216033935547 + 1.8976157889483147e-07j),
(-5.9838657379150391 - 4.3475332260131836j), (-2.0879747867584229 -
6.4261269569396973j),
(1.9100792407989502 - 5.8786196708679199j), (4.5814824104309082 -
3.3286411762237549j),
(5.1967458724975586 - 1.3684227440080576e-07j),
(3.8647139072418213 + 2.8078789710998535j), (1.3594740629196167 +
4.1840314865112305j),
(-1.2544282674789429 + 3.8607344627380371j), (-3.0366206169128418 +
2.2062335014343262j),
(-3.4781389236450195 + 1.1194014604143376e-07j),
(-2.6133756637573242 - 1.8987287282943726j), (-0.9293016791343689 -
2.8600969314575195j),
(0.86727333068847656 - 2.6691930294036865j), (2.1243946552276611 -
1.5434627532958984j),
(2.4633183479309082 - 8.6486437567145913e-08j),
(1.8744727373123169 + 1.3618841171264648j), (0.67528903484344482 +
2.0783262252807617j),
(-0.63866174221038818 + 1.965599536895752j), (-1.5857341289520264 +
1.152103066444397j),
(-1.8640764951705933 +
7.6355092915036948e-08j), (-1.4381576776504517 -
1.0448826551437378j),
(-0.52529704570770264 -
1.6166983842849731j), (0.50366902351379395 - 1.5501341819763184j),
(1.26766037940979 - 0.92100900411605835j))
sampling_freq = 100
src1 = analog.sig_source_c(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100.0)
dst1 = blocks.vector_sink_c()
head = blocks.head(gr.sizeof_gr_complex, int(5 * sampling_freq * 0.10))
agc = analog.agc_cc(1e-3, 1, 1)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 4)
def test_005(self):
''' Test the complex AGC loop (attack and decay rate inputs) '''
tb = self.tb
expected_result = \
((100.000244140625+7.2191943445432116e-07j),
(0.80881959199905396+0.58764183521270752j),
(0.30894950032234192+0.95084899663925171j),
(-0.30895623564720154+0.95086973905563354j),
(-0.80887287855148315+0.58768033981323242j),
(-0.99984413385391235+5.850709250410091e-09j),
(-0.80889981985092163-0.58770018815994263j),
(-0.30897706747055054-0.95093393325805664j),
(0.30898112058639526-0.95094609260559082j),
(0.80893135070800781-0.58772283792495728j),
(0.99990922212600708-8.7766354184282136e-09j),
(0.80894720554351807+0.58773452043533325j),
(0.30899339914321899+0.95098406076431274j),
(-0.30899572372436523+0.95099133253097534j),
(-0.80896598100662231+0.58774799108505249j),
(-0.99994778633117676+1.4628290578855285e-08j),
(-0.80897533893585205-0.58775502443313599j),
(-0.30900305509567261-0.95101380348205566j),
(0.30900448560714722-0.95101797580718994j),
(0.80898630619049072-0.58776277303695679j),
(0.99997037649154663-1.7554345532744264e-08j),
(0.80899184942245483+0.58776694536209106j),
(0.30900871753692627+0.95103120803833008j),
(-0.30900952219963074+0.95103377103805542j),
(-0.8089984655380249+0.58777159452438354j),
(-0.99998390674591064+2.3406109050938539e-08j),
(-0.809001624584198-0.58777409791946411j),
(-0.30901208519935608-0.95104163885116577j),
(0.30901262164115906-0.95104306936264038j),
(0.80900543928146362-0.587776780128479j),
(0.99999171495437622-2.6332081404234486e-08j),
(0.80900734663009644+0.58777821063995361j),
(0.30901408195495605+0.95104765892028809j),
(-0.30901429057121277+0.95104855298995972j),
(-0.80900967121124268+0.58777981996536255j),
(-0.99999648332595825+3.2183805842578295e-08j),
(-0.80901080369949341-0.58778077363967896j),
(-0.30901527404785156-0.95105135440826416j),
(0.30901545286178589-0.95105189085006714j),
(0.80901217460632324-0.58778166770935059j),
(0.99999916553497314-3.5109700036173308e-08j),
(0.809012770652771+0.58778214454650879j),
(0.30901595950126648+0.9510534405708313j),
(-0.30901598930358887+0.95105385780334473j),
(-0.80901366472244263+0.58778274059295654j),
(-1.0000008344650269+4.0961388947380328e-08j),
(-0.8090139627456665-0.58778303861618042j),
(-0.30901634693145752-0.95105475187301636j),
(0.30901640653610229-0.95105493068695068j),
(0.80901449918746948-0.5877833366394043j))
sampling_freq = 100
src1 = analog.sig_source_c(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100)
dst1 = blocks.vector_sink_c()
head = blocks.head(gr.sizeof_gr_complex, int(5*sampling_freq * 0.10))
agc = analog.agc2_cc(1e-2, 1e-3, 1, 1)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 4)
def test_pll_carriertracking(self):
expected_result = ((1.00000238419+7.21919457547e-09j),
(0.998025715351+0.062790453434j),
(0.992777824402+0.119947694242j),
(0.985192835331+0.171441286802j),
(0.976061582565+0.217501848936j),
(0.966034710407+0.258409559727j),
(0.95565611124+0.294477283955j),
(0.945357382298+0.326030552387j),
(0.935475051403+0.353395611048j),
(0.926258146763+0.376889169216j),
(0.917895197868+0.39681750536j),
(0.910515546799+0.413470208645j),
(0.904196679592+0.427117019892j),
(0.898972511292+0.438006043434j),
(0.894769787788+0.446523308754j),
(0.891652584076+0.452715367079j),
(0.8895829916+0.456773489714j),
(0.888502895832+0.458873122931j),
(0.888343691826+0.459175437689j),
(0.889035582542+0.457833081484j),
(0.890497922897+0.454985737801j),
(0.892645597458+0.450762689114j),
(0.895388305187+0.445282936096j),
(0.898648142815+0.438664674759j),
(0.902342617512+0.431016951799j),
(0.906392872334+0.422441422939j),
(0.910642921925+0.413191765547j),
(0.915039420128+0.403358519077j),
(0.919594764709+0.392864197493j),
(0.92425006628+0.381792247295j),
(0.928944349289+0.370217680931j),
(0.933634519577+0.358220815659j),
(0.938279032707+0.345874190331j),
(0.942840516567+0.333247303963j),
(0.947280526161+0.32040438056j),
(0.951574921608+0.307409763336j),
(0.955703914165+0.294323593378j),
(0.959648966789+0.281201630831j),
(0.963392794132+0.268095195293j),
(0.966880619526+0.255221515894j),
(0.970162451267+0.242447137833j),
(0.973235487938+0.229809194803j),
(0.97609680891+0.217341512442j),
(0.978744983673+0.20507311821j),
(0.981189727783+0.193033605814j),
(0.983436584473+0.181248426437j),
(0.985490739346+0.169738590717j),
(0.987353682518+0.158523857594j),
(0.989041447639+0.147622272372j),
(0.990563035011+0.137049794197j),
(0.991928339005+0.126818582416j),
(0.993117690086+0.117111675441j),
(0.994156062603+0.107930034399j),
(0.995076179504+0.0990980416536j),
(0.995887458324+0.0906178802252j),
(0.996591091156+0.0824909061193j),
(0.997202515602+0.0747182965279j),
(0.997730851173+0.0672992765903j),
(0.998185396194+0.0602316558361j),
(0.99856698513+0.0535135567188j),
(0.998885989189+0.0471420884132j),
(0.99915266037+0.0411129891872j),
(0.999372899532+0.0354214012623j),
(0.999548316002+0.0300626158714j),
(0.999680638313+0.0252036750317j),
(0.999784469604+0.020652115345j),
(0.999865531921+0.0163950324059j),
(0.999923825264+0.0124222636223j),
(0.999960243702+0.00872156023979j),
(0.999983668327+0.00528120994568j),
(0.999997138977+0.00209015607834j),
(1.00000119209-0.00086285173893j),
(0.999992132187-0.00358882546425j),
(0.999979138374-0.00609711557627j),
(0.999963641167-0.00839691981673j),
(0.999947249889-0.0104993218556j),
(0.999924004078-0.0122378543019j),
(0.999904811382-0.0136305987835j),
(0.999888062477-0.0148707330227j),
(0.9998742342-0.0159679055214j),
(0.999856114388-0.0169314742088j),
(0.999839782715-0.0177700817585j),
(0.999826967716-0.0184917747974j),
(0.999818325043-0.0191045701504j),
(0.999807476997-0.0196143388748j),
(0.999797284603-0.0200265944004j),
(0.999791204929-0.0203481912613j),
(0.99978852272-0.0205836892128j),
(0.99978530407-0.0207380950451j),
(0.999785065651-0.0206423997879j),
(0.999787807465-0.0204866230488j),
(0.999794304371-0.0202808082104j),
(0.999800384045-0.0200312435627j),
(0.999803245068-0.0197458267212j),
(0.9998087883-0.0194311738014j),
(0.999816894531-0.0190933048725j),
(0.999825954437-0.0187371373177j),
(0.999829888344-0.0183679759502j),
(0.999835848808-0.017987690866j),
(0.999844014645-0.0176006518304j))
sampling_freq = 10e3
freq = sampling_freq / 100
loop_bw = math.pi/100.0
maxf = 1
minf = -1
src = analog.sig_source_c(sampling_freq, analog.GR_COS_WAVE, freq, 1.0)
pll = analog.pll_carriertracking_cc(loop_bw, maxf, minf)
head = blocks.head(gr.sizeof_gr_complex, int (freq))
dst = blocks.vector_sink_c()
self.tb.connect(src, pll, head)
self.tb.connect(head, dst)
self.tb.run()
dst_data = dst.data()
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 5)
def test_pll_refout(self):
expected_result = (
(1 + 0j), (1 + 6.4087357643e-10j),
(0.999985277653 + 0.00542619498447j),
(0.999868750572 + 0.0162021834403j),
(0.99948567152 + 0.0320679470897j), (0.99860727787 +
0.0527590736747j),
(0.996953129768 + 0.0780025869608j), (0.994203746319 +
0.107512556016j),
(0.990011692047 + 0.140985429287j), (0.984013140202 +
0.178095817566j),
(0.975838363171 + 0.218493551016j),
(0.965121984482 + 0.261800557375j), (0.95151245594 +
0.307610183954j),
(0.934681296349 + 0.355486690998j), (0.914401650429 +
0.404808044434j),
(0.890356600285 + 0.455263823271j), (0.862329125404 +
0.506348133087j),
(0.830152392387 + 0.557536482811j), (0.793714106083 +
0.608290970325j),
(0.752960026264 + 0.658066213131j), (0.707896590233 +
0.706316053867j),
(0.658591926098 + 0.752500295639j), (0.605175673962 +
0.796091973782j),
(0.547837555408 + 0.836584687233j), (0.48682525754 +
0.873499393463j),
(0.42244040966 + 0.906390726566j), (0.355197101831 +
0.934791445732j),
(0.285494059324 + 0.958380460739j), (0.213591173291 +
0.976923108101j),
(0.139945343137 + 0.990159213543j), (0.065038472414 +
0.997882783413j),
(-0.0106285437942 + 0.999943494797j), (-0.0865436866879 +
0.996248066425j),
(-0.162189796567 + 0.986759603024j), (-0.23705175519 +
0.971496999264j),
(-0.310622543097 + 0.950533330441j), (-0.38240903616 +
0.923993110657j),
(-0.451937526464 + 0.89204955101j), (-0.518758952618 +
0.854920566082j),
(-0.582311093807 + 0.812966048717j), (-0.642372369766 +
0.76639264822j),
(-0.698591887951 + 0.715520322323j), (-0.750654160976 +
0.660695314407j),
(-0.798280358315 + 0.602286040783j), (-0.841228663921 +
0.540679454803j),
(-0.87929558754 + 0.476276367903j), (-0.912315964699 +
0.409486919641j),
(-0.940161883831 + 0.340728074312j), (-0.962742805481 +
0.270418733358j),
(-0.980004072189 + 0.198977485299j), (-0.991925954819 +
0.126818284392j),
(-0.99851256609 +
0.0545223206282j), (-0.999846458435 -
0.0175215266645j), (-0.996021270752 -
0.0891158208251j),
(-0.987133920193 -
0.159895718098j), (-0.973306238651 -
0.2295101583j), (-0.954683184624 -
0.297624111176j),
(-0.931430280209 -
0.363919824362j), (-0.903732538223 -
0.428097635508j), (-0.871792256832 -
0.489875763655j),
(-0.835827112198 -
0.548992812634j), (-0.796068251133 -
0.605206847191j), (-0.752758979797 -
0.658296227455j),
(-0.706152498722 -
0.70805978775j), (-0.656641483307 -
0.754202902317j), (-0.604367733002 -
0.79670548439j),
(-0.549597978592 -
0.835429251194j), (-0.492602348328 -
0.870254516602j), (-0.433654457331 -
0.901079237461j),
(-0.373029649258 -
0.927819430828j), (-0.31100410223 -
0.950408577919j), (-0.247853919864 -
0.968797445297j),
(-0.183855071664 -
0.982953369617j), (-0.119282215834 -
0.992860376835j), (-0.0544078871608 -
0.998518764973j),
(0.0104992967099 -
0.999944865704j), (0.0749994292855 -
0.997183561325j), (0.138844624162 -
0.990314185619j),
(0.201967850327 -
0.979392170906j), (0.264124274254 -
0.964488625526j), (0.325075358152 -
0.945688128471j),
(0.3845885396 - 0.92308807373j), (0.442438393831 - 0.89679890871j),
(0.498407125473 -
0.866943061352j), (0.552284479141 -
0.833655714989j), (0.603869199753 -
0.797083437443j),
(0.652970373631 -
0.757383465767j), (0.69940674305 -
0.714723825455j), (0.743007957935 -
0.66928255558j),
(0.78350687027 -
0.62138313055j), (0.820889055729 -
0.571087777615j), (0.855021059513 -
0.51859331131j),
(0.885780930519 -
0.46410369873j), (0.913058102131 -
0.407829582691j), (0.936754107475 -
0.349988251925j),
(0.956783294678 -
0.290801793337j), (0.973072886467 -
0.230497643352j), (0.985563337803 -
0.169307261705j),
(0.9942086339 - 0.1074674353j), (0.9989772439 - 0.0452152714133j))
sampling_freq = 10e3
freq = sampling_freq / 100
loop_bw = math.pi / 100.0
maxf = 1
minf = -1
src = analog.sig_source_c(sampling_freq, analog.GR_COS_WAVE, freq, 1.0)
pll = analog.pll_refout_cc(loop_bw, maxf, minf)
head = blocks.head(gr.sizeof_gr_complex, int(freq))
dst = blocks.vector_sink_c()
self.tb.connect(src, pll, head)
self.tb.connect(head, dst)
self.tb.run()
dst_data = dst.data()
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 4)
def test_pll_refout(self):
expected_result = ((1+0j),
(1+6.4087357643e-10j),
(0.999985277653+0.00542619498447j),
(0.999868750572+0.0162021834403j),
(0.99948567152+0.0320679470897j),
(0.99860727787+0.0527590736747j),
(0.996953129768+0.0780025869608j),
(0.994203746319+0.107512556016j),
(0.990011692047+0.140985429287j),
(0.984013140202+0.178095817566j),
(0.975838363171+0.218493551016j),
(0.965121984482+0.261800557375j),
(0.95151245594+0.307610183954j),
(0.934681296349+0.355486690998j),
(0.914401650429+0.404808044434j),
(0.890356600285+0.455263823271j),
(0.862329125404+0.506348133087j),
(0.830152392387+0.557536482811j),
(0.793714106083+0.608290970325j),
(0.752960026264+0.658066213131j),
(0.707896590233+0.706316053867j),
(0.658591926098+0.752500295639j),
(0.605175673962+0.796091973782j),
(0.547837555408+0.836584687233j),
(0.48682525754+0.873499393463j),
(0.42244040966+0.906390726566j),
(0.355197101831+0.934791445732j),
(0.285494059324+0.958380460739j),
(0.213591173291+0.976923108101j),
(0.139945343137+0.990159213543j),
(0.065038472414+0.997882783413j),
(-0.0106285437942+0.999943494797j),
(-0.0865436866879+0.996248066425j),
(-0.162189796567+0.986759603024j),
(-0.23705175519+0.971496999264j),
(-0.310622543097+0.950533330441j),
(-0.38240903616+0.923993110657j),
(-0.451937526464+0.89204955101j),
(-0.518758952618+0.854920566082j),
(-0.582311093807+0.812966048717j),
(-0.642372369766+0.76639264822j),
(-0.698591887951+0.715520322323j),
(-0.750654160976+0.660695314407j),
(-0.798280358315+0.602286040783j),
(-0.841228663921+0.540679454803j),
(-0.87929558754+0.476276367903j),
(-0.912315964699+0.409486919641j),
(-0.940161883831+0.340728074312j),
(-0.962742805481+0.270418733358j),
(-0.980004072189+0.198977485299j),
(-0.991925954819+0.126818284392j),
(-0.99851256609+0.0545223206282j),
(-0.999846458435-0.0175215266645j),
(-0.996021270752-0.0891158208251j),
(-0.987133920193-0.159895718098j),
(-0.973306238651-0.2295101583j),
(-0.954683184624-0.297624111176j),
(-0.931430280209-0.363919824362j),
(-0.903732538223-0.428097635508j),
(-0.871792256832-0.489875763655j),
(-0.835827112198-0.548992812634j),
(-0.796068251133-0.605206847191j),
(-0.752758979797-0.658296227455j),
(-0.706152498722-0.70805978775j),
(-0.656641483307-0.754202902317j),
(-0.604367733002-0.79670548439j),
(-0.549597978592-0.835429251194j),
(-0.492602348328-0.870254516602j),
(-0.433654457331-0.901079237461j),
(-0.373029649258-0.927819430828j),
(-0.31100410223-0.950408577919j),
(-0.247853919864-0.968797445297j),
(-0.183855071664-0.982953369617j),
(-0.119282215834-0.992860376835j),
(-0.0544078871608-0.998518764973j),
(0.0104992967099-0.999944865704j),
(0.0749994292855-0.997183561325j),
(0.138844624162-0.990314185619j),
(0.201967850327-0.979392170906j),
(0.264124274254-0.964488625526j),
(0.325075358152-0.945688128471j),
(0.3845885396-0.92308807373j),
(0.442438393831-0.89679890871j),
(0.498407125473-0.866943061352j),
(0.552284479141-0.833655714989j),
(0.603869199753-0.797083437443j),
(0.652970373631-0.757383465767j),
(0.69940674305-0.714723825455j),
(0.743007957935-0.66928255558j),
(0.78350687027-0.62138313055j),
(0.820889055729-0.571087777615j),
(0.855021059513-0.51859331131j),
(0.885780930519-0.46410369873j),
(0.913058102131-0.407829582691j),
(0.936754107475-0.349988251925j),
(0.956783294678-0.290801793337j),
(0.973072886467-0.230497643352j),
(0.985563337803-0.169307261705j),
(0.9942086339-0.1074674353j),
(0.9989772439-0.0452152714133j))
sampling_freq = 10e3
freq = sampling_freq / 100
loop_bw = math.pi/100.0
maxf = 1
minf = -1
src = analog.sig_source_c(sampling_freq, analog.GR_COS_WAVE, freq, 1.0)
pll = analog.pll_refout_cc(loop_bw, maxf, minf)
head = blocks.head(gr.sizeof_gr_complex, int (freq))
dst = blocks.vector_sink_c()
self.tb.connect(src, pll, head)
self.tb.connect(head, dst)
self.tb.run()
dst_data = dst.data()
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 4)
def test_002(self):
''' Test the floating point AGC loop (single rate input) '''
tb = self.tb
expected_result = (
7.2191943445432116e-07,
58.837181091308594,
89.700050354003906,
81.264183044433594,
45.506141662597656,
4.269894304798072e-07,
-42.948936462402344,
-65.50335693359375,
-59.368724822998047,
-33.261005401611328,
-4.683740257860336e-07,
31.423542022705078,
47.950984954833984,
43.485683441162109,
24.378345489501953,
5.7254135299444897e-07,
-23.062990188598633,
-35.218441009521484,
-31.964075088500977,
-17.934831619262695,
-5.0591745548445033e-07,
16.998210906982422,
25.982204437255859,
23.606258392333984,
13.260685920715332,
4.9936483037527069e-07,
-12.59880542755127,
-19.28221321105957,
-17.54347038269043,
-9.8700437545776367,
-4.188150626305287e-07,
9.4074573516845703,
14.422011375427246,
13.145503044128418,
7.41046142578125,
3.8512698097292741e-07,
-7.0924453735351562,
-10.896408081054688,
-9.9552040100097656,
-5.6262712478637695,
-3.1982864356905338e-07,
5.4131259918212891,
8.3389215469360352,
7.6409502029418945,
4.3320145606994629,
2.882407841298118e-07,
-4.194943904876709,
-6.4837145805358887,
-5.9621825218200684,
-3.3931560516357422)
sampling_freq = 100
src1 = analog.sig_source_f(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100.0)
dst1 = blocks.vector_sink_f ()
head = blocks.head (gr.sizeof_float, int (5*sampling_freq * 0.10))
agc = analog.agc_ff(1e-3, 1, 1)
tb.connect (src1, head)
tb.connect (head, agc)
tb.connect (agc, dst1)
tb.run ()
dst_data = dst1.data ()
self.assertFloatTuplesAlmostEqual (expected_result, dst_data, 4)
def test_pll_freqdet(self):
expected_result = (0.0,
4.33888922882e-08,
0.367369994515,
1.08135249597,
2.10983253908,
3.42221529438,
4.98940390402,
6.78379190842,
8.77923286024,
10.9510106794,
13.2758363182,
15.7317829127,
18.2982902299,
20.9561068599,
23.6755271122,
26.452952094,
29.2731265301,
32.1219053479,
34.9862418188,
37.8540971414,
40.7144315483,
43.5571390869,
46.3730179743,
49.1537231663,
51.8917218889,
54.58026103,
57.2015358514,
59.7513664199,
62.2380533124,
64.657612252,
67.006640002,
69.2822432184,
71.4820384499,
73.6041047056,
75.6469478817,
77.6094829742,
79.4909866472,
81.2911031615,
83.0097850853,
84.6355598352,
86.1820937186,
87.6504420946,
89.0418441206,
90.3577286819,
91.5996432431,
92.7692775646,
93.8684162704,
94.8989269904,
95.8627662892,
96.7619381633,
97.598505899,
98.362769679,
99.0579904444,
99.6992633875,
100.288805948,
100.828805921,
101.321421457,
101.76878699,
102.17300138,
102.536116055,
102.860158727,
103.147085962,
103.398830608,
103.617254366,
103.792467691,
103.939387906,
104.060030865,
104.15631756,
104.230085975,
104.283067372,
104.316933727,
104.333238432,
104.333440018,
104.318914008,
104.290941063,
104.250742554,
104.187634452,
104.103822339,
104.013227468,
103.916810336,
103.815448432,
103.709936239,
103.600997093,
103.489283183,
103.375351833,
103.259712936,
103.142828952,
103.025091195,
102.90686726,
102.776726069,
102.648078982,
102.521459607,
102.397294831,
102.275999684,
102.157882471,
102.043215927,
101.93218978,
101.824958181,
101.72159228,
101.622151366)
sampling_freq = 10e3
freq = sampling_freq / 100
loop_bw = math.pi/100.0
maxf = 1
minf = -1
src = analog.sig_source_c(sampling_freq, analog.GR_COS_WAVE, freq, 1.0)
pll = analog.pll_freqdet_cf(loop_bw, maxf, minf)
head = blocks.head(gr.sizeof_float, int (freq))
dst = blocks.vector_sink_f()
self.tb.connect(src, pll, head)
self.tb.connect(head, dst)
self.tb.run()
dst_data = dst.data()
# convert it from normalized frequency to absolute frequency (Hz)
dst_data = [i*(sampling_freq/(2*math.pi)) for i in dst_data]
self.assertFloatTuplesAlmostEqual(expected_result, dst_data, 3)
def test_001(self):
''' Test the complex AGC loop (single rate input) '''
tb = self.tb
expected_result = (
(100.000244140625+7.2191943445432116e-07j),
(72.892257690429688+52.959323883056641j),
(25.089065551757812+77.216217041015625j),
(-22.611061096191406+69.589706420898438j),
(-53.357715606689453+38.766635894775391j),
(-59.458671569824219+3.4792964243024471e-07j),
(-43.373462677001953-31.512666702270508j),
(-14.94139289855957-45.984889984130859j),
(13.478158950805664-41.48150634765625j),
(31.838506698608398-23.132022857666016j),
(35.519271850585938-3.1176801940091536e-07j),
(25.942903518676758+18.848621368408203j),
(8.9492912292480469+27.5430908203125j),
(-8.0852642059326172+24.883890151977539j),
(-19.131628036499023+13.899936676025391j),
(-21.383295059204102+3.1281737733479531e-07j),
(-15.650330543518066-11.370632171630859j),
(-5.4110145568847656-16.65339469909668j),
(4.9008159637451172-15.083160400390625j),
(11.628337860107422-8.4484796524047852j),
(13.036135673522949-2.288476110834381e-07j),
(9.5726661682128906+6.954948902130127j),
(3.3216962814331055+10.223132133483887j),
(-3.0204284191131592+9.2959251403808594j),
(-7.1977195739746094+5.2294478416442871j),
(-8.1072216033935547+1.8976157889483147e-07j),
(-5.9838657379150391-4.3475332260131836j),
(-2.0879747867584229-6.4261269569396973j),
(1.9100792407989502-5.8786196708679199j),
(4.5814824104309082-3.3286411762237549j),
(5.1967458724975586-1.3684227440080576e-07j),
(3.8647139072418213+2.8078789710998535j),
(1.3594740629196167+4.1840314865112305j),
(-1.2544282674789429+3.8607344627380371j),
(-3.0366206169128418+2.2062335014343262j),
(-3.4781389236450195+1.1194014604143376e-07j),
(-2.6133756637573242-1.8987287282943726j),
(-0.9293016791343689-2.8600969314575195j),
(0.86727333068847656-2.6691930294036865j),
(2.1243946552276611-1.5434627532958984j),
(2.4633183479309082-8.6486437567145913e-08j),
(1.8744727373123169+1.3618841171264648j),
(0.67528903484344482+2.0783262252807617j),
(-0.63866174221038818+1.965599536895752j),
(-1.5857341289520264+1.152103066444397j),
(-1.8640764951705933+7.6355092915036948e-08j),
(-1.4381576776504517-1.0448826551437378j),
(-0.52529704570770264-1.6166983842849731j),
(0.50366902351379395-1.5501341819763184j),
(1.26766037940979-0.92100900411605835j))
sampling_freq = 100
src1 = analog.sig_source_c(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100.0)
dst1 = blocks.vector_sink_c()
head = blocks.head(gr.sizeof_gr_complex, int (5*sampling_freq * 0.10))
agc = analog.agc_cc(1e-3, 1, 1)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 4)
def test_004(self):
''' Test the floating point AGC loop (attack and decay rate inputs) '''
tb = self.tb
expected_result = \
(7.2191943445432116e-07,
58.837181091308594,
40.194305419921875,
2.9183335304260254,
0.67606079578399658,
8.6260438791896377e-09,
-1.4542514085769653,
-1.9210131168365479,
-1.0450780391693115,
-0.61939650774002075,
-1.2590258613442984e-08,
1.4308931827545166,
1.9054338932037354,
1.0443156957626343,
0.61937344074249268,
2.0983527804219193e-08,
-1.4308838844299316,
-1.9054274559020996,
-1.0443152189254761,
-0.61937344074249268,
-2.5180233009791664e-08,
1.4308837652206421,
1.9054274559020996,
1.0443154573440552,
0.61937344074249268,
3.3573645197293445e-08,
-1.4308838844299316,
-1.9054274559020996,
-1.0443152189254761,
-0.61937350034713745,
-3.7770352179222755e-08,
1.4308837652206421,
1.9054274559020996,
1.0443154573440552,
0.61937350034713745,
4.6163762590367696e-08,
-1.4308838844299316,
-1.9054274559020996,
-1.0443153381347656,
-0.61937344074249268,
-5.0360466019583328e-08,
1.4308837652206421,
1.9054274559020996,
1.0443155765533447,
0.61937344074249268,
5.8753879983441948e-08,
-1.4308837652206421,
-1.9054274559020996,
-1.0443153381347656,
-0.61937344074249268)
sampling_freq = 100
src1 = analog.sig_source_f(sampling_freq, analog.GR_SIN_WAVE,
sampling_freq * 0.10, 100)
dst1 = blocks.vector_sink_f()
head = blocks.head(gr.sizeof_float, int(5*sampling_freq * 0.10))
agc = analog.agc2_ff(1e-2, 1e-3, 1, 1)
tb.connect(src1, head)
tb.connect(head, agc)
tb.connect(agc, dst1)
tb.run()
dst_data = dst1.data()
self.assertFloatTuplesAlmostEqual(expected_result, dst_data, 4)
