def test02(self): # Test QPSK sync M = 4 theta = 0 loop_bw = cmath.pi/100.0 fmin = -0.5 fmax = 0.5 mu = 0.5 gain_mu = 0.01 omega = 2 gain_omega = 0.001 omega_rel = 0.001 self.test = digital.mpsk_receiver_cc(M, theta, loop_bw, fmin, fmax, mu, gain_mu, omega, gain_omega, omega_rel) data = 10000*[complex( 0.707, 0.707), complex(-0.707, 0.707), complex(-0.707, -0.707), complex( 0.707, -0.707)] data = [0.5*d for d in data] self.src = gr.vector_source_c(data, False) self.snk = gr.vector_sink_c() # pulse shaping interpolation filter nfilts = 32 excess_bw = 0.35 ntaps = 11 * int(omega*nfilts) rrc_taps0 = filter.firdes.root_raised_cosine( nfilts, nfilts, 1.0, excess_bw, ntaps) rrc_taps1 = filter.firdes.root_raised_cosine( 1, omega, 1.0, excess_bw, 11*omega) self.rrc0 = filter.pfb_arb_resampler_ccf(omega, rrc_taps0) self.rrc1 = filter.fir_filter_ccf(1, rrc_taps1) self.tb.connect(self.src, self.rrc0, self.rrc1, self.test, self.snk) self.tb.run() expected_result = 10000*[complex(-0.5, +0.0), complex(+0.0, -0.5), complex(+0.5, +0.0), complex(+0.0, +0.5)] # get data after a settling period dst_data = self.snk.data()[200:] # Only compare last Ncmp samples Ncmp = 1000 len_e = len(expected_result) len_d = len(dst_data) expected_result = expected_result[len_e - Ncmp - 1:-1] dst_data = dst_data[len_d - Ncmp:] #for e,d in zip(expected_result, dst_data): # print "{0:+.02f} {1:+.02f}".format(e, d) self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 1)
def test01(self): # Test BPSK sync M = 2 theta = 0 loop_bw = cmath.pi / 100.0 fmin = -0.5 fmax = 0.5 mu = 0.5 gain_mu = 0.01 omega = 2 gain_omega = 0.001 omega_rel = 0.001 self.test = digital.mpsk_receiver_cc(M, theta, loop_bw, fmin, fmax, mu, gain_mu, omega, gain_omega, omega_rel) data = 10000 * [complex(1, 0), complex(-1, 0)] #data = [2*random.randint(0,1)-1 for x in xrange(10000)] self.src = blocks.vector_source_c(data, False) self.snk = blocks.vector_sink_c() # pulse shaping interpolation filter nfilts = 32 excess_bw = 0.35 ntaps = 11 * int(omega * nfilts) rrc_taps0 = filter.firdes.root_raised_cosine(nfilts, nfilts, 1.0, excess_bw, ntaps) rrc_taps1 = filter.firdes.root_raised_cosine(1, omega, 1.0, excess_bw, 11 * omega) self.rrc0 = filter.pfb_arb_resampler_ccf(omega, rrc_taps0) self.rrc1 = filter.fir_filter_ccf(1, rrc_taps1) self.tb.connect(self.src, self.rrc0, self.rrc1, self.test, self.snk) self.tb.run() expected_result = [0.5 * d for d in data] dst_data = self.snk.data() # Only compare last Ncmp samples Ncmp = 1000 len_e = len(expected_result) len_d = len(dst_data) expected_result = expected_result[len_e - Ncmp - 1:-1] dst_data = dst_data[len_d - Ncmp:] #for e,d in zip(expected_result, dst_data): # print "{0:+.02f} {1:+.02f}".format(e, d) self.assertComplexTuplesAlmostEqual(expected_result, dst_data, 1)
def test_fir_filter_ccf_002(self): src_data = 40*[1+1j, 2+2j, 3+3j, 4+4j] expected_data = ((0.5+0.5j), (5.5+5.5j), (10.5+10.5j), (15.5+15.5j), (20.5+20.5j), (25.5+25.5j), (30.5+30.5j), (35.5+35.5j), (40.5+40.5j), (45.5+45.5j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j)) src = gr.vector_source_c(src_data) op = filter.fir_filter_ccf(4, 20*[0.5, 0.5]) dst = gr.vector_sink_c() self.tb.connect(src, op, dst) self.tb.run() result_data = dst.data() self.assertComplexTuplesAlmostEqual(expected_data, result_data, 5)
def test_fir_filter_ccf_002(self): src_data = 40*[1+1j, 2+2j, 3+3j, 4+4j] expected_data = ((0.5+0.5j), (5.5+5.5j), (10.5+10.5j), (15.5+15.5j), (20.5+20.5j), (25.5+25.5j), (30.5+30.5j), (35.5+35.5j), (40.5+40.5j), (45.5+45.5j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j)) src = blocks.vector_source_c(src_data) op = filter.fir_filter_ccf(4, 20*[0.5, 0.5]) dst = blocks.vector_sink_c() self.tb.connect(src, op, dst) self.tb.run() result_data = dst.data() self.assertComplexTuplesAlmostEqual(expected_data, result_data, 5)
def test_fir_filter_ccf_001(self): src_data = 40*[1+1j, 2+2j, 3+3j, 4+4j] expected_data = ((0.5+0.5j), (1.5+1.5j), (3+3j), (5+5j), (5.5+5.5j), (6.5+6.5j), (8+8j), (10+10j), (10.5+10.5j), (11.5+11.5j), (13+13j), (15+15j), (15.5+15.5j), (16.5+16.5j), (18+18j), (20+20j), (20.5+20.5j), (21.5+21.5j), (23+23j), (25+25j), (25.5+25.5j), (26.5+26.5j), (28+28j), (30+30j), (30.5+30.5j), (31.5+31.5j), (33+33j), (35+35j), (35.5+35.5j), (36.5+36.5j), (38+38j), (40+40j), (40.5+40.5j), (41.5+41.5j), (43+43j), (45+45j), (45.5+45.5j), (46.5+46.5j), (48+48j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j)) src = gr.vector_source_c(src_data) op = filter.fir_filter_ccf(1, 20*[0.5, 0.5]) dst = gr.vector_sink_c() self.tb.connect(src, op, dst) self.tb.run() result_data = dst.data() self.assertComplexTuplesAlmostEqual(expected_data, result_data, 5)
def test_fir_filter_ccf_001(self): src_data = 40*[1+1j, 2+2j, 3+3j, 4+4j] expected_data = ((0.5+0.5j), (1.5+1.5j), (3+3j), (5+5j), (5.5+5.5j), (6.5+6.5j), (8+8j), (10+10j), (10.5+10.5j), (11.5+11.5j), (13+13j), (15+15j), (15.5+15.5j), (16.5+16.5j), (18+18j), (20+20j), (20.5+20.5j), (21.5+21.5j), (23+23j), (25+25j), (25.5+25.5j), (26.5+26.5j), (28+28j), (30+30j), (30.5+30.5j), (31.5+31.5j), (33+33j), (35+35j), (35.5+35.5j), (36.5+36.5j), (38+38j), (40+40j), (40.5+40.5j), (41.5+41.5j), (43+43j), (45+45j), (45.5+45.5j), (46.5+46.5j), (48+48j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j), (50+50j)) src = blocks.vector_source_c(src_data) op = filter.fir_filter_ccf(1, 20*[0.5, 0.5]) dst = blocks.vector_sink_c() self.tb.connect(src, op, dst) self.tb.run() result_data = dst.data() self.assertComplexTuplesAlmostEqual(expected_data, result_data, 5)