Commit 1be69de2 authored by luis_pereira87's avatar luis_pereira87

Code cleanup

parent e1d88d6f
...@@ -141,7 +141,7 @@ int pss_search_time_nr(int **rxdata, ///rx data in time domain ...@@ -141,7 +141,7 @@ int pss_search_time_nr(int **rxdata, ///rx data in time domain
int *eNB_id, int *eNB_id,
int *f_off); int *f_off);
int nr_adjust_pss_synch(PHY_VARS_NR_UE *ue, int *f_off); int nr_adjust_pss_synch(PHY_VARS_NR_UE *ue);
#endif #endif
#undef EXTERN #undef EXTERN
......
...@@ -976,7 +976,7 @@ int pss_search_time_nr(int **rxdata, ///rx data in time domain ...@@ -976,7 +976,7 @@ int pss_search_time_nr(int **rxdata, ///rx data in time domain
return(peak_position); return(peak_position);
} }
int nr_adjust_pss_synch(PHY_VARS_NR_UE *ue, int *f_off) { int nr_adjust_pss_synch(PHY_VARS_NR_UE *ue) {
int **rxdata = ue->common_vars.rxdata; int **rxdata = ue->common_vars.rxdata;
NR_DL_FRAME_PARMS *frame_parms = &ue->frame_parms; NR_DL_FRAME_PARMS *frame_parms = &ue->frame_parms;
...@@ -1030,36 +1030,8 @@ int nr_adjust_pss_synch(PHY_VARS_NR_UE *ue, int *f_off) { ...@@ -1030,36 +1030,8 @@ int nr_adjust_pss_synch(PHY_VARS_NR_UE *ue, int *f_off) {
} }
} }
// Fractional frequency offset computation according to Cross-correlation Synchronization Algorithm Using PSS LOG_D(NR_PHY,"nr_adjust_pss_synch: peak_position %d, peak_position_diff %d, peak_value %d dB, avg %d dB\n",
// Shoujun Huang, Yongtao Su, Ying He and Shan Tang, "Joint time and frequency offset estimation in LTE downlink," peak_position, peak_position_diff, dB_fixed64(peak_value),dB_fixed64(avg));
// 7th International Conference on Communications and Networking in China, 2012.
// Computing cross-correlation at peak on half the symbol size for first half of data
int64_t result1 = dot_product64((short*)primary_synchro_time_nr[Nid2],
(short*) &(rxdata[0][peak_position]),
frame_parms->ofdm_symbol_size>>1,
shift);
// Computing cross-correlation at peak on half the symbol size for data shifted by half symbol size
// as it is real and complex it is necessary to shift by a value equal to symbol size to obtain such shift
int64_t result2 = dot_product64((short*)primary_synchro_time_nr[Nid2]+(frame_parms->ofdm_symbol_size),
(short*) &(rxdata[0][peak_position])+frame_parms->ofdm_symbol_size,
frame_parms->ofdm_symbol_size>>1,
shift);
int64_t re1 = ((int*) &result1)[0];
int64_t re2 = ((int*) &result2)[0];
int64_t im1 = ((int*) &result1)[1];
int64_t im2 = ((int*) &result2)[1];
// Estimation of fractional frequency offset: angle[(result1)'*(result2)]/pi
double ffo_est = atan2(re1*im2-re2*im1,re1*re2+im1*im2)/M_PI;
// Computing absolute value of frequency offset
*f_off = ffo_est*frame_parms->subcarrier_spacing;
LOG_D(NR_PHY,"nr_adjust_pss_synch: peak_position %d, peak_position_diff %d, peak_value %d dB, avg %d dB, ffo %lf\n",
peak_position, peak_position_diff, dB_fixed64(peak_value),dB_fixed64(avg),ffo_est);
return peak_position_diff; return peak_position_diff;
} }
...@@ -1648,7 +1648,7 @@ int phy_procedures_nrUE_RX(PHY_VARS_NR_UE *ue, ...@@ -1648,7 +1648,7 @@ int phy_procedures_nrUE_RX(PHY_VARS_NR_UE *ue,
if ((ue->decode_MIB == 1) && slot_pbch) { if ((ue->decode_MIB == 1) && slot_pbch) {
nr_adjust_pss_synch(ue, (int *)&ue->common_vars.freq_offset); nr_adjust_pss_synch(ue);
LOG_D(PHY," ------ Decode MIB: frame.slot %d.%d ------ \n", frame_rx%1024, nr_slot_rx); LOG_D(PHY," ------ Decode MIB: frame.slot %d.%d ------ \n", frame_rx%1024, nr_slot_rx);
nr_ue_pbch_procedures(gNB_id, ue, proc, 0); nr_ue_pbch_procedures(gNB_id, ue, proc, 0);
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment