<p>pespin <strong>merged</strong> this change.</p><p><a href="https://gerrit.osmocom.org/c/osmo-trx/+/14674">View Change</a></p><div style="white-space:pre-wrap">Approvals:
laforge: Looks good to me, approved
fixeria: Looks good to me, but someone else must approve
Jenkins Builder: Verified
</div><pre style="font-family: monospace,monospace; white-space: pre-wrap;">sigProcLib: Add C/I (Carrier-to-Interference ratio) computation<br><br>Related: OS#4006<br>Change-Id: Ib4ceec553f2e5f77bf3f6777724968456a180f5e<br>---<br>M Transceiver52M/sigProcLib.cpp<br>M Transceiver52M/sigProcLib.h<br>2 files changed, 51 insertions(+), 13 deletions(-)<br><br></pre><pre style="font-family: monospace,monospace; white-space: pre-wrap;"><span>diff --git a/Transceiver52M/sigProcLib.cpp b/Transceiver52M/sigProcLib.cpp</span><br><span>index c22009a..860bbe4 100644</span><br><span>--- a/Transceiver52M/sigProcLib.cpp</span><br><span>+++ b/Transceiver52M/sigProcLib.cpp</span><br><span>@@ -1454,6 +1454,33 @@</span><br><span> };</span><br><span> </span><br><span> /*</span><br><span style="color: hsl(120, 100%, 40%);">+ * Computes C/I (Carrier-to-Interference ratio) in dB (deciBels).</span><br><span style="color: hsl(120, 100%, 40%);">+ * It is computed from the training sequence of each received burst,</span><br><span style="color: hsl(120, 100%, 40%);">+ * by comparing the "ideal" training sequence with the actual one.</span><br><span style="color: hsl(120, 100%, 40%);">+ */</span><br><span style="color: hsl(120, 100%, 40%);">+static float computeCI(const signalVector *burst, CorrelationSequence *sync,</span><br><span style="color: hsl(120, 100%, 40%);">+ float toa, int start, complex xcorr)</span><br><span style="color: hsl(120, 100%, 40%);">+{</span><br><span style="color: hsl(120, 100%, 40%);">+ float S, C;</span><br><span style="color: hsl(120, 100%, 40%);">+ int ps;</span><br><span style="color: hsl(120, 100%, 40%);">+</span><br><span style="color: hsl(120, 100%, 40%);">+ /* Integer position where the sequence starts */</span><br><span style="color: hsl(120, 100%, 40%);">+ ps = start + 1 - sync->sequence->size() + (int)roundf(toa);</span><br><span style="color: hsl(120, 100%, 40%);">+</span><br><span style="color: hsl(120, 100%, 40%);">+ /* Estimate Signal power */</span><br><span style="color: hsl(120, 100%, 40%);">+ S = 0.0f;</span><br><span style="color: hsl(120, 100%, 40%);">+ for (int i=0, j=ps; i<(int)sync->sequence->size(); i++,j++)</span><br><span style="color: hsl(120, 100%, 40%);">+ S += (*burst)[j].norm2();</span><br><span style="color: hsl(120, 100%, 40%);">+ S /= sync->sequence->size();</span><br><span style="color: hsl(120, 100%, 40%);">+</span><br><span style="color: hsl(120, 100%, 40%);">+ /* Esimate Carrier power */</span><br><span style="color: hsl(120, 100%, 40%);">+ C = xcorr.norm2() / ((sync->sequence->size() - 1) * sync->gain.abs());</span><br><span style="color: hsl(120, 100%, 40%);">+</span><br><span style="color: hsl(120, 100%, 40%);">+ /* Interference = Signal - Carrier, so C/I = C / (S - C) */</span><br><span style="color: hsl(120, 100%, 40%);">+ return 3.0103f * log2f(C / (S - C));</span><br><span style="color: hsl(120, 100%, 40%);">+}</span><br><span style="color: hsl(120, 100%, 40%);">+</span><br><span style="color: hsl(120, 100%, 40%);">+/*</span><br><span> * Detect a burst based on correlation and peak-to-average ratio</span><br><span> *</span><br><span> * For one sampler-per-symbol, perform fast peak detection (no interpolation)</span><br><span>@@ -1468,6 +1495,8 @@</span><br><span> {</span><br><span> const signalVector *corr_in;</span><br><span> signalVector *dec = NULL;</span><br><span style="color: hsl(120, 100%, 40%);">+ complex xcorr;</span><br><span style="color: hsl(120, 100%, 40%);">+ int rc = 1;</span><br><span> </span><br><span> if (sps == 4) {</span><br><span> dec = downsampleBurst(burst);</span><br><span>@@ -1480,35 +1509,42 @@</span><br><span> /* Correlate */</span><br><span> if (!convolve(corr_in, sync->sequence, &corr,</span><br><span> CUSTOM, start, len)) {</span><br><span style="color: hsl(0, 100%, 40%);">- delete dec;</span><br><span style="color: hsl(0, 100%, 40%);">- return -1;</span><br><span style="color: hsl(120, 100%, 40%);">+ rc = -1;</span><br><span style="color: hsl(120, 100%, 40%);">+ goto del_ret;</span><br><span> }</span><br><span> </span><br><span style="color: hsl(0, 100%, 40%);">- delete dec;</span><br><span style="color: hsl(0, 100%, 40%);">-</span><br><span> /* Running at the downsampled rate at this point */</span><br><span> sps = 1;</span><br><span> </span><br><span> /* Peak detection - place restrictions at correlation edges */</span><br><span> ebp->amp = fastPeakDetect(corr, &ebp->toa);</span><br><span> </span><br><span style="color: hsl(0, 100%, 40%);">- if ((ebp->toa < 3 * sps) || (ebp->toa > len - 3 * sps))</span><br><span style="color: hsl(0, 100%, 40%);">- return 0;</span><br><span style="color: hsl(120, 100%, 40%);">+ if ((ebp->toa < 3 * sps) || (ebp->toa > len - 3 * sps)) {</span><br><span style="color: hsl(120, 100%, 40%);">+ rc = 0;</span><br><span style="color: hsl(120, 100%, 40%);">+ goto del_ret;</span><br><span style="color: hsl(120, 100%, 40%);">+ }</span><br><span> </span><br><span style="color: hsl(0, 100%, 40%);">- /* Peak -to-average ratio */</span><br><span style="color: hsl(0, 100%, 40%);">- if (computePeakRatio(&corr, sps, ebp->toa, ebp->amp) < thresh)</span><br><span style="color: hsl(0, 100%, 40%);">- return 0;</span><br><span style="color: hsl(120, 100%, 40%);">+ /* Peak-to-average ratio */</span><br><span style="color: hsl(120, 100%, 40%);">+ if (computePeakRatio(&corr, sps, ebp->toa, ebp->amp) < thresh) {</span><br><span style="color: hsl(120, 100%, 40%);">+ rc = 0;</span><br><span style="color: hsl(120, 100%, 40%);">+ goto del_ret;</span><br><span style="color: hsl(120, 100%, 40%);">+ }</span><br><span> </span><br><span style="color: hsl(0, 100%, 40%);">- /* Compute peak-to-average ratio. Reject if we don't have enough values */</span><br><span style="color: hsl(0, 100%, 40%);">- ebp->amp = peakDetect(corr, &ebp->toa, NULL);</span><br><span style="color: hsl(120, 100%, 40%);">+ /* Refine TOA and correlation value */</span><br><span style="color: hsl(120, 100%, 40%);">+ xcorr = peakDetect(corr, &ebp->toa, NULL);</span><br><span style="color: hsl(120, 100%, 40%);">+</span><br><span style="color: hsl(120, 100%, 40%);">+ /* Compute C/I */</span><br><span style="color: hsl(120, 100%, 40%);">+ ebp->ci = computeCI(corr_in, sync, ebp->toa, start, xcorr);</span><br><span> </span><br><span> /* Normalize our channel gain */</span><br><span style="color: hsl(0, 100%, 40%);">- ebp->amp = ebp->amp / sync->gain;</span><br><span style="color: hsl(120, 100%, 40%);">+ ebp->amp = xcorr / sync->gain;</span><br><span> </span><br><span> /* Compensate for residuate time lag */</span><br><span> ebp->toa = ebp->toa - sync->toa;</span><br><span> </span><br><span style="color: hsl(0, 100%, 40%);">- return 1;</span><br><span style="color: hsl(120, 100%, 40%);">+del_ret:</span><br><span style="color: hsl(120, 100%, 40%);">+ delete dec;</span><br><span style="color: hsl(120, 100%, 40%);">+ return rc;</span><br><span> }</span><br><span> </span><br><span> static float maxAmplitude(const signalVector &burst)</span><br><span>@@ -1563,6 +1599,7 @@</span><br><span> } else if (!rc) {</span><br><span> ebp->amp = 0.0f;</span><br><span> ebp->toa = 0.0f;</span><br><span style="color: hsl(120, 100%, 40%);">+ ebp->ci = 0.0f;</span><br><span> return clipping?-SIGERR_CLIP:SIGERR_NONE;</span><br><span> }</span><br><span> </span><br><span>diff --git a/Transceiver52M/sigProcLib.h b/Transceiver52M/sigProcLib.h</span><br><span>index 6d4ca9b..fd9a5f0 100644</span><br><span>--- a/Transceiver52M/sigProcLib.h</span><br><span>+++ b/Transceiver52M/sigProcLib.h</span><br><span>@@ -111,6 +111,7 @@</span><br><span> complex amp;</span><br><span> float toa;</span><br><span> uint8_t tsc;</span><br><span style="color: hsl(120, 100%, 40%);">+ float ci;</span><br><span> };</span><br><span> /**</span><br><span> 8-PSK/GMSK/RACH burst detector</span><br><span></span><br></pre><p>To view, visit <a href="https://gerrit.osmocom.org/c/osmo-trx/+/14674">change 14674</a>. To unsubscribe, or for help writing mail filters, visit <a href="https://gerrit.osmocom.org/settings">settings</a>.</p><div itemscope itemtype="http://schema.org/EmailMessage"><div itemscope itemprop="action" itemtype="http://schema.org/ViewAction"><link itemprop="url" href="https://gerrit.osmocom.org/c/osmo-trx/+/14674"/><meta itemprop="name" content="View Change"/></div></div>
<div style="display:none"> Gerrit-Project: osmo-trx </div>
<div style="display:none"> Gerrit-Branch: master </div>
<div style="display:none"> Gerrit-Change-Id: Ib4ceec553f2e5f77bf3f6777724968456a180f5e </div>
<div style="display:none"> Gerrit-Change-Number: 14674 </div>
<div style="display:none"> Gerrit-PatchSet: 6 </div>
<div style="display:none"> Gerrit-Owner: pespin <pespin@sysmocom.de> </div>
<div style="display:none"> Gerrit-Reviewer: Jenkins Builder </div>
<div style="display:none"> Gerrit-Reviewer: fixeria <axilirator@gmail.com> </div>
<div style="display:none"> Gerrit-Reviewer: laforge <laforge@gnumonks.org> </div>
<div style="display:none"> Gerrit-Reviewer: pespin <pespin@sysmocom.de> </div>
<div style="display:none"> Gerrit-Reviewer: tnt <tnt@246tNt.com> </div>
<div style="display:none"> Gerrit-MessageType: merged </div>