def test_002_pcfich(self):
print "test_002_pcfich"
# some constants
cell_id = 124
N_ant = 2
style = "tx_diversity"
vlen = 16
ns = 0
# new top_block because even the interface changes
self.tb2 = gr.top_block()
# generate test data together with the expected output
data = []
exp_res = []
for cfi in range(4):
cfi_seq = lte_test.get_cfi_sequence(cfi + 1)
scr_cfi_seq = lte_test.scramble_cfi_sequence(cfi_seq, cell_id, ns)
mod_cfi_seq = lte_test.qpsk_modulation(scr_cfi_seq)
lay_cfi_seq = lte_test.layer_mapping(mod_cfi_seq, N_ant, style)
lay_cfi_prep = lte_test.prepare_for_demapper_block(
lay_cfi_seq, N_ant, style)
exp_res.extend(lay_cfi_prep)
pc_cfi_seq = lte_test.pre_coding(lay_cfi_seq, N_ant, style)
pc_cfi_seq = [
pc_cfi_seq[0][i] + pc_cfi_seq[1][i]
for i in range(len(pc_cfi_seq[0]))
]
data.extend(pc_cfi_seq)
# dummy channel estimates
intu2 = [complex(1, 0)] * len(data)
intu3 = [complex(1, 0)] * len(data)
# get blocks
self.src1 = blocks.vector_source_c(data, False, vlen)
self.src2 = blocks.vector_source_c(intu2, False, vlen)
self.src3 = blocks.vector_source_c(intu3, False, vlen)
self.pd = lte.pre_decoder_vcvc(1, 1, vlen, style)
self.snk = blocks.vector_sink_c(vlen)
# connect all blocks
self.tb2.connect(self.src1, (self.pd, 0))
self.tb2.connect(self.src2, (self.pd, 1))
self.tb2.connect(self.src3, (self.pd, 2))
self.tb2.connect(self.pd, self.snk)
self.pd.set_N_ant(N_ant)
# run flowgraph
self.tb2.run()
# compare result with expected result
res = self.snk.data()
self.assertComplexTuplesAlmostEqual(res, exp_res)
def test_002_pcfich(self):
print "test_002_pcfich"
# some constants
cell_id = 124
N_ant = 2
style = "tx_diversity"
vlen = 16
ns = 0
# new top_block because even the interface changes
self.tb2 = gr.top_block()
# generate test data together with the expected output
data = []
exp_res = []
for cfi in range(4):
cfi_seq = lte_test.get_cfi_sequence(cfi + 1)
scr_cfi_seq = lte_test.scramble_cfi_sequence(cfi_seq, cell_id, ns)
mod_cfi_seq = lte_test.qpsk_modulation(scr_cfi_seq)
lay_cfi_seq = lte_test.layer_mapping(mod_cfi_seq, N_ant, style)
lay_cfi_prep = lte_test.prepare_for_demapper_block(lay_cfi_seq, N_ant, style)
exp_res.extend(lay_cfi_prep)
pc_cfi_seq = lte_test.pre_coding(lay_cfi_seq, N_ant, style)
pc_cfi_seq = [pc_cfi_seq[0][i] + pc_cfi_seq[1][i] for i in range(len(pc_cfi_seq[0]))]
data.extend(pc_cfi_seq)
# dummy channel estimates
intu2 = [complex(1, 0)] * len(data)
intu3 = [complex(1, 0)] * len(data)
# get blocks
self.src1 = blocks.vector_source_c(data, False, vlen)
self.src2 = blocks.vector_source_c(intu2, False, vlen)
self.src3 = blocks.vector_source_c(intu3, False, vlen)
self.pd = lte.pre_decoder_vcvc(1, vlen, style)
self.snk = blocks.vector_sink_c(vlen)
# connect all blocks
self.tb2.connect(self.src1, (self.pd, 0))
self.tb2.connect(self.src2, (self.pd, 1))
self.tb2.connect(self.src3, (self.pd, 2))
self.tb2.connect(self.pd, self.snk)
self.pd.set_N_ant(N_ant)
# run flowgraph
self.tb2.run()
# compare result with expected result
res = self.snk.data()
self.assertComplexTuplesAlmostEqual(res, exp_res)
def setUp (self):
self.tb = gr.top_block ()
# Input 1, PBCH frame
mod=scipy.io.loadmat('/home/demel/exchange/matlab_pbch_frame.mat')
#mat_u1=range(240)
mat_u1=tuple(mod['pbch_frame1'].flatten())
mat_d=range(240)
for idx, val in enumerate(mat_u1):
mat_d[idx]=val
intu1=tuple(mat_d)
# Input 2, CE values for antenna port 1
mod=scipy.io.loadmat('/home/demel/exchange/matlab_pbch_ce1.mat')
#mat_u2=range(240)
mat_u2=tuple(mod['pbch_ce1'].flatten())
mat_d=range(240)
for idx, val in enumerate(mat_u2):
mat_d[idx]=val
intu2=tuple(mat_d)
# Input 2, CE values for antenna port 1
mod=scipy.io.loadmat('/home/demel/exchange/matlab_pbch_ce2.mat')
#mat_u2=range(240)
mat_u2=tuple(mod['pbch_ce1'].flatten())
mat_d=range(240)
for idx, val in enumerate(mat_u2):
mat_d[idx]=val
intu3=tuple(mat_d)
print str(intu1[0]) + "\t" + str(intu2[0]) + "\t" + str(intu3[0])
self.src1 = gr.vector_source_c( intu1, False, 240)
self.src2 = gr.vector_source_c( intu2, False, 240)
#self.src3 = gr.vector_source_c( intu3, False, 240)
self.pd = lte_swig.pre_decoder_vcvc(1,'tx_diversity')
self.snk = gr.vector_sink_c(240)
self.tb.connect(self.src1,(self.pd,0) )
self.tb.connect(self.src2,(self.pd,1) )
#self.tb.connect(self.src3,(self.pd,2) )
self.tb.connect(self.pd,self.snk)
def setUp(self):
self.tb = gr.top_block()
N_ant = 2
vlen = 240
style = 'tx_diversity'
intu1 = [0] * vlen
intu2 = [0] * vlen
intu3 = [0] * vlen
self.src1 = blocks.vector_source_c(intu1, False, vlen)
self.src2 = blocks.vector_source_c(intu2, False, vlen)
self.src3 = blocks.vector_source_c(intu3, False, vlen)
self.pd = lte.pre_decoder_vcvc(1, N_ant, vlen, style)
self.snk = blocks.vector_sink_c(vlen)
self.tb.connect(self.src1, (self.pd, 0))
self.tb.connect(self.src2, (self.pd, 1))
self.tb.connect(self.src3, (self.pd, 2))
self.tb.connect(self.pd, self.snk)
def setUp(self):
self.tb = gr.top_block()
N_ant = 2
vlen = 240
style = 'tx_diversity'
intu1 = [0] * vlen
intu2 = [0] * vlen
intu3 = [0] * vlen
self.src1 = blocks.vector_source_c(intu1, False, vlen)
self.src2 = blocks.vector_source_c(intu2, False, vlen)
self.src3 = blocks.vector_source_c(intu3, False, vlen)
self.pd = lte.pre_decoder_vcvc(N_ant, vlen, style)
self.snk = blocks.vector_sink_c(vlen)
self.tb.connect(self.src1, (self.pd, 0))
self.tb.connect(self.src2, (self.pd, 1))
self.tb.connect(self.src3, (self.pd, 2))
self.tb.connect(self.pd, self.snk)
def setUp (self):
#print os.getpid()
#raw_input("Press the ANY key to continue")
setuptime = time.clock()
self.tb = gr.top_block ()
print "\nqa_multi_block_test START"
offset = 1133
fftl = 512
cell_id = 124
N_rb_dl = 6
I = 768
D = 1000
# This first part is for syncing!
#self.src = gr.vector_source_c(intu,False,1)
infile1 = '/home/demel/Dokumente/Messungen/Messung_LTE_2012-05-23_12:48:19.dat' # successful!
infile2 = '/home/demel/Dokumente/Messungen/Messung_LTE_2012-05-23_12:48:07.dat' # failed!
infile3 = '/home/demel/Dokumente/Messungen/Messung_LTE_2012-05-23_12:46:30.dat' # successful!
infile4 = '/home/johannes/lte/data/Messung_LTE_2012-05-23_12:47:32.dat' # successful!
infile5 = '/home/demel/Dokumente/Messungen/Messung_LTE_2012-05-23_12:49:57.dat' # successful!
samps = (fftl*7*2*10+10*fftl)*300 +offset
input_file = infile3
print input_file
self.src = gr.file_source(gr.sizeof_gr_complex,input_file, False)
self.resample = blks2.rational_resampler_ccc(I, D)
self.head1 = gr.head(gr.sizeof_gr_complex, samps)
#self.tag = lte_swig.tag_symbol_cc(offset,fftl)
self.sync = lte_swig.cp_time_freq_sync_cc(fftl)
self.pss = lte.hier_pss_sync_cc(fftl)
self.est = lte.hier_freq_estimate_cc(fftl)
# This is the actual place to initialize self.sss
# Nevertheless it it the last block to be initialized because it needs pointers to some others.
# self.sss = lte.hier_sss_sync_cc(fftl)
#This part does still process all data (more ore less)
self.rcp = lte_swig.remove_cp_cvc(fftl)
self.fft = gr.fft_vcc(fftl,True,window.rectangular(fftl),False,1)
self.ext = lte_swig.extract_occupied_tones_vcvc(N_rb_dl,fftl)
self.eq = lte.linear_OFDM_estimator_hier_vcvc(N_rb_dl)#cell_id,
#self.eq.set_cell_id(cell_id)
#self.sh1 = gr.skiphead(gr.sizeof_gr_complex*12*N_rb_dl,7)
#self.sh2 = gr.skiphead(gr.sizeof_gr_complex*12*N_rb_dl,7)
#self.sh3 = gr.skiphead(gr.sizeof_gr_complex*12*N_rb_dl,7)
# After the next block only the PBCH is processed!
self.demux= lte_swig.pbch_demux_vcvc(N_rb_dl)#cell_id,
#self.demux.set_cell_id(cell_id)
self.pd1 = lte_swig.pre_decoder_vcvc(1,'tx_diversity')
self.pd2 = lte_swig.pre_decoder_vcvc(2,'tx_diversity')
self.ldm1 = lte_swig.layer_demapper_vcvc(1,'tx_diversity')
self.ldm2 = lte_swig.layer_demapper_vcvc(2,'tx_diversity')
self.int = gr.interleave(240*8)
self.dmd = lte_swig.qpsk_soft_demod_vcvf()
self.ld = lte_swig.descrambling_vfvf()
#self.ld.set_cell_id(cell_id)
self.lru = lte_swig.rate_unmatch_vff()
# After Viterbi decoding, only bits are processed!
self.vit = lte.viterbi_vfvb()
self.crc = lte_swig.crc_calculator_vbvb()
self.mib = lte_swig.mib_unpack_vb()
#############################################################
# This last block is just a daemon to propagate the calculated cell id
#############################################################
self.daemon = lte.cell_id_daemon(self.eq.get_eq(),self.demux, self.ld)
self.sss = lte.hier_sss_sync_cc(self.daemon, fftl)
# declaration of additional sink blocks for tests
#self.snk1 = gr.vector_sink_f(120)
#self.snk2 = gr.vector_sink_b(40)
# Connect all blocks together!
self.tb.connect(self.src,self.resample,self.head1,self.sync,self.pss,self.est,self.sss,self.rcp,self.fft,self.ext,self.eq)#
self.tb.connect( (self.eq,0), (self.demux,0) )
self.tb.connect( (self.eq,1), (self.demux,1) )
self.tb.connect( (self.eq,2), (self.demux,2) )
self.tb.connect( (self.demux,0) ,(self.pd1,0) )
self.tb.connect( (self.demux,1) ,(self.pd1,1) )
self.tb.connect( (self.demux,0) ,(self.pd2,0) )
self.tb.connect( (self.demux,1) ,(self.pd2,1) )
self.tb.connect( (self.demux,2) ,(self.pd2,2) )
self.tb.connect(self.pd1,self.ldm1,(self.int,0) )
self.tb.connect(self.pd2,self.ldm2,(self.int,1) )
self.tb.connect(self.int,self.dmd,self.ld,self.lru,self.vit,self.crc)
self.tb.connect( (self.crc,0),(self.mib,0) )
self.tb.connect( (self.crc,1),(self.mib,1) )
# additional sink blocks for testing
#self.tb.connect(self.lru,self.snk1)
#self.tb.connect(self.vit,self.snk2)
print "\n\nsetuptime = " + str(time.clock() - setuptime) + "\n\n"
def setUp(self):
#print os.getpid()
#raw_input("Press the ANY key to continue")
setuptime = time.clock()
self.tb = gr.top_block()
print "\nqa_multi_block_test START"
offset = 1133
fftl = 512
cell_id = 124
N_rb_dl = 6
I = 768
D = 1000
# This first part is for syncing!
#self.src = gr.vector_source_c(intu,False,1)
infile1 = '/home/demel/Dokumente/Messungen/Messung_LTE_2012-05-23_12:48:19.dat' # successful!
infile2 = '/home/demel/Dokumente/Messungen/Messung_LTE_2012-05-23_12:48:07.dat' # failed!
infile3 = '/home/demel/Dokumente/Messungen/Messung_LTE_2012-05-23_12:46:30.dat' # successful!
infile4 = '/home/johannes/lte/data/Messung_LTE_2012-05-23_12:47:32.dat' # successful!
infile5 = '/home/demel/Dokumente/Messungen/Messung_LTE_2012-05-23_12:49:57.dat' # successful!
samps = (fftl * 7 * 2 * 10 + 10 * fftl) * 300 + offset
input_file = infile3
print input_file
self.src = gr.file_source(gr.sizeof_gr_complex, input_file, False)
self.resample = blks2.rational_resampler_ccc(I, D)
self.head1 = gr.head(gr.sizeof_gr_complex, samps)
#self.tag = lte_swig.tag_symbol_cc(offset,fftl)
self.sync = lte_swig.cp_time_freq_sync_cc(fftl)
self.pss = lte.hier_pss_sync_cc(fftl)
self.est = lte.hier_freq_estimate_cc(fftl)
# This is the actual place to initialize self.sss
# Nevertheless it it the last block to be initialized because it needs pointers to some others.
# self.sss = lte.hier_sss_sync_cc(fftl)
#This part does still process all data (more ore less)
self.rcp = lte_swig.remove_cp_cvc(fftl)
self.fft = gr.fft_vcc(fftl, True, window.rectangular(fftl), False, 1)
self.ext = lte_swig.extract_occupied_tones_vcvc(N_rb_dl, fftl)
self.eq = lte.linear_OFDM_estimator_hier_vcvc(N_rb_dl) #cell_id,
#self.eq.set_cell_id(cell_id)
#self.sh1 = gr.skiphead(gr.sizeof_gr_complex*12*N_rb_dl,7)
#self.sh2 = gr.skiphead(gr.sizeof_gr_complex*12*N_rb_dl,7)
#self.sh3 = gr.skiphead(gr.sizeof_gr_complex*12*N_rb_dl,7)
# After the next block only the PBCH is processed!
self.demux = lte_swig.pbch_demux_vcvc(N_rb_dl) #cell_id,
#self.demux.set_cell_id(cell_id)
self.pd1 = lte_swig.pre_decoder_vcvc(1, 'tx_diversity')
self.pd2 = lte_swig.pre_decoder_vcvc(2, 'tx_diversity')
self.ldm1 = lte_swig.layer_demapper_vcvc(1, 'tx_diversity')
self.ldm2 = lte_swig.layer_demapper_vcvc(2, 'tx_diversity')
self.int = gr.interleave(240 * 8)
self.dmd = lte_swig.qpsk_soft_demod_vcvf()
self.ld = lte_swig.descrambling_vfvf()
#self.ld.set_cell_id(cell_id)
self.lru = lte_swig.rate_unmatch_vff()
# After Viterbi decoding, only bits are processed!
self.vit = lte.viterbi_vfvb()
self.crc = lte_swig.crc_calculator_vbvb()
self.mib = lte_swig.mib_unpack_vb()
#############################################################
# This last block is just a daemon to propagate the calculated cell id
#############################################################
self.daemon = lte.cell_id_daemon(self.eq.get_eq(), self.demux, self.ld)
self.sss = lte.hier_sss_sync_cc(self.daemon, fftl)
# declaration of additional sink blocks for tests
#self.snk1 = gr.vector_sink_f(120)
#self.snk2 = gr.vector_sink_b(40)
# Connect all blocks together!
self.tb.connect(self.src, self.resample, self.head1, self.sync,
self.pss, self.est, self.sss, self.rcp, self.fft,
self.ext, self.eq) #
self.tb.connect((self.eq, 0), (self.demux, 0))
self.tb.connect((self.eq, 1), (self.demux, 1))
self.tb.connect((self.eq, 2), (self.demux, 2))
self.tb.connect((self.demux, 0), (self.pd1, 0))
self.tb.connect((self.demux, 1), (self.pd1, 1))
self.tb.connect((self.demux, 0), (self.pd2, 0))
self.tb.connect((self.demux, 1), (self.pd2, 1))
self.tb.connect((self.demux, 2), (self.pd2, 2))
self.tb.connect(self.pd1, self.ldm1, (self.int, 0))
self.tb.connect(self.pd2, self.ldm2, (self.int, 1))
self.tb.connect(self.int, self.dmd, self.ld, self.lru, self.vit,
self.crc)
self.tb.connect((self.crc, 0), (self.mib, 0))
self.tb.connect((self.crc, 1), (self.mib, 1))
# additional sink blocks for testing
#self.tb.connect(self.lru,self.snk1)
#self.tb.connect(self.vit,self.snk2)
print "\n\nsetuptime = " + str(time.clock() - setuptime) + "\n\n"