CN104869088A - 用于快速变化信道和窄带宽信道的低复杂性gmsk接收机 - Google Patents

Abstract

用于快速变化信道和窄带宽信道的低复杂性GMSK接收机。一种用于对遭受符号间干扰的接收信号中的符号进行解码的方法，该接收信号包括以符号间隔隔开的多个符号，该方法包括以下步骤：对所述接收信号进行滤波以形成表示所述符号的样本，该样本包括所述符号和在所述接收信号中紧先于该符号以及接续该符号的和；将所述样本乘以在所述符号间隔的整数倍处交替地具有以下各项的复正弦函数：(i)零实分量和非零虚分量以及(ii)非零实分量和零虚分量，以生成复数值；以及依赖于所述复数值对所述符号进行解码。

Description

用于快速变化信道和窄带宽信道的低复杂性GMSK接收机

技术领域

本发明涉及用于对遭受符号间干扰的信号进行解码的方法和装置。

背景技术

高斯最小移位键控(GMSK)是载波信号的相位根据要发送的信息而变化的连续相位调制方案。GMSK使用具有窄带宽和锐截止的预调制高斯滤波器。高斯滤波器抑制高频分量并且使输出功率谱变得更紧凑。另外，因为GMSK仅使用相位调制，所以它是恒定包络调制方案。因此GMSK是能实现高效低功率发送机设计并且因此对于必须实现非常长电池寿命的电池供电的通信设备而言有吸引力的调制方案。

对通过连续相位调制方案所生成的信号进行解调因各个符号的初始相位是根据预先发送的符号的累积相位来确定的事实而变复杂。因此，接收机在不考虑发送符号的整个序列的情况下不能够对一个符号做出判定。GMSK还特别易于发生符号间干扰(ISI)，这使接收机设计复杂化。均衡化通常是需要的，并且通常由诸如维特比解码器的最大似然均衡器来执行。然而，这些均衡器不很适于变化信道条件或由发送机和接收机引入的相位噪声，并且由于相位噪声而导致的小的多普勒漂移或相位改变能够导致显著的性能劣化。这在GMSK符号率低使得信号能够在窄带宽信道中适应时特别成问题，因为给定的多普勒漂移或相位噪声性能将在各个符号周期期间导致更大的相位改变。

因此，存在对于用于对遭受符号间干扰的信号进行解码的改进的方法的需要。

发明内容

根据一个实施方式，提供了一种用于对遭受符号间干扰的接收信号中的符号进行解码的方法，该接收信号包括由符号间隔隔开的多个符号，该方法包括以下步骤：对所述接收信号进行滤波以形成表示所述符号的样本，该样本包括所述符号和在所述接收信号中紧先于该符号以及接续该符号的符号的加权和；将所述样本乘以在所述符号间隔的整数倍处交替地具有以下各项的复正弦函数，以生成复数值：(i)零实分量和非零虚分量以及(ii)非零实分量和零虚分量；以及依赖于所述复数值对所述符号进行解码。

所述接收信号可以是这样的，即所述样本包括为显著地实的或显著地虚的符号以及在所述接收信号中紧先于该符号以及接续该符号的显著地为实的或虚的中的另一方的符号的加权和。

所述方法可以包括以下步骤：如果所述符号的加权分量是显著地实的，并且所述复正弦函数在对应于所述样本的所述符号间隔的整数倍处具有零实分量和非零虚分量，则依赖于所述复数值的虚部对所述符号进行解码。

所述方法可以包括以下步骤：如果所述符号的加权分量是显著地实的，并且所述复正弦函数在对应于所述样本的所述符号间隔的整数倍处具有非零实分量和零虚分量，则依赖于所述复数值的实部对所述符号进行解码。

所述方法可以包括以下步骤：如果所述符号的加权分量是显著地虚的，并且所述复正弦函数在对应于所述样本的所述符号间隔的整数倍处具有零实分量和非零虚分量，则依赖于所述复数值的实部对所述符号进行解码。

所述方法可以包括以下步骤：如果所述符号的加权分量是显著地虚的，并且所述复正弦函数在对应于所述样本的所述符号间隔的整数倍处具有非零实分量和非零虚分量，则依赖于所述复数值的虚部对所述符号进行解码。

所述方法可以包括以下步骤：依赖于所述复数值的实分量和虚分量中的一方对所述符号进行解码并且依赖于所述复数值的实分量和虚分量中的另一方来估计所述符号间干扰。

所述复正弦函数的一个周期可以等于所述符号间隔的四倍的整数倍。

所述调制方法可以是GMSK。

所述GMSK调制的时间带宽乘积可以大于或等于0.3。

所述方法可以包括以下步骤：依赖于所述复正弦函数来生成期望符号的星座；监测所述多个符号与期望符号的星座之间的偏差；以及依赖于所监测到的偏差来跟踪通过其接收到所述信号的信道的改变。

所述方法可以包括：生成期望符号，该期望符号表示在假设所接收到的符号和经解码的符号相同的情况下将已生成的复数值。

所述方法可以包括通过将经解码的符号与从在所述接收信号中紧先于所接收到的符号以及接续所接收到的符号的所述符号生成的解码的符号合计来生成所述期望符号。

所述方法可以包括以下步骤：将所述期望符号除以所接收到的符号；将信道估计减去除法的输出；将已减信号提供给自适应滤波器作为所述信道估计的误差的表示，该自适应滤波器利用该表示来处理所述接收信号；以及依赖于已减信号来更新所述信道估计。

根据第二实施方式，提供了一种用于被构造为接收遭受符号间干扰的信号的接收机的解码器，所述接收信号包括以规则间隔隔开的多个符号，该解码器包括：滤波器，该滤波器被构造为对所述接收信号进行滤波以形成表示所述符号的样本，该样本包括所述符号和在所述接收信号中紧先于该符号以及接续该符号的符号的加权和；乘法器，该乘法器被构造为将所述样本乘以在符号间隔的整数倍处交替地具有以下各项的复正弦函数，以生成复数值：(i)零实分量和非零虚分量以及(ii)非零实分量和零虚分量；以及判定单元，该判定单元被构造为依赖于所述复数值对所述符号进行解码。

附图说明

现在将参照附图通过示例描述本发明。附图中：

图1示出了用于对接收信号中的符号进行解码的方法；

图2a示出了包括数据符号和符号间干扰的完整星座的示例；

图2b示出了在存在噪声的情况下包括数据符号和符号间干扰的完整星座的示例；

图3示出了用于跟踪信道改变的方法的示例；以及

图4示出了接收机结构的示例。

具体实施方式

在图1中示出了用于对遭受符号间干扰的接收信号中的符号进行解码的方法的示例。该方法包括接收表示多个符号的信号，在各个符号之间具有给定间隔(步骤101)。该符号间隔可以被表示为T_s。在这个上下文中措辞“符号”被用来意指数据的任何单元，所以各个“符号”可以表示任何数量的一个或更多个比特。接收信号被滤波以产生样本的序列(步骤102)。序列中的各个样本优选地表示这些符号中的一个以及至少在接收信号中紧先于该符号以及接续该符号的符号的和。样本可以由符号的加权和形成。组合邻近符号以形成各个符号反映了符号间干扰：各个符号受在它任一侧的那些符号影响。

样本的序列然后乘以复正弦函数(步骤103)。优选地复正弦函数在符号间隔的整数倍处具有零实分量和非零虚分量，或非零实分量和零虚分量。这个乘法的结果是与样本对应的复数值的序列。该方法然后依赖于序列中与其对应的复数值对各个符号进行解码(步骤104)。

该方法导致特别适合于遭受平坦衰落的信道的低复杂性接收机结构。该方法对变化信道非常鲁棒并且对由发送机和接收机引入的相位噪声鲁棒。该方法还使得针对线性调制方案(诸如BPSK、QPSK等)而开发的标准信号处理算法能够被容易地适用于GMSK或其它连续相位调制方案。

现在将具体参照使用了GMSK对接收信号进行调制的实施方式来描述方法和接收机装置。然而，这不旨在为限制性的，因为本文所描述的方法和装置可以同样地适用于其它调制方法，特别是连续相位调制方法。

首先，将对GMSK调制进行描述。

要发送的数据被适合地形成为非归零(NRZ)序列α_k∈{-1，+1}。这是调制数据，并且能够对它区别地进行编码。

调制数据被适合地具体实现为矩形脉冲：

然后通过高斯滤波器对调制数据进行滤波以获得频率脉冲高斯脉冲由下式给出：

$g\left(t\right)=\frac{1}{\sqrt{2\mathrm{\π}{\mathrm{\δ}}^2}}\exp \left(\frac{{-t}^2}{{2\mathrm{\δ}}^2}\right)---\left(2\right)$

其中：

$\mathrm{\δ}=\frac{\sqrt{2}}{2\mathrm{\π\β}}$

β是带宽-时间乘积。它适合地为0.3。

然后根据下式对信号进行调制：

其中h是0.5的调制指数并且T_s是符号间隔。

利用劳伦特分解，能够将这个信号分解成脉冲振幅调制信号的和。这使得GMSK的非线性调制能够以可辨识的线性方式表达。GMSK的分解能够被写为：

$s\left(t\right)={\mathrm{\Σ}}_{k=-\∞}^{\∞}{\mathrm{\Σ}}_{q=0}^Q{c}_q\left(k\right)h_q(t-{\mathrm{kT}}_s)---\left(4\right)$

其中，c_q(t)和hq₍t)分别是针对第q个脉冲振幅调制信号的信号和整形滤波器。

对于时间带宽乘积0.3，q＝0的信号包含能量的大多数。所发送的信号因此能够近似为：

$s\left(t\right)\≈{\mathrm{\Σ}}_{k=-\∞}^{\∞}{c}_0\left(k\right)h_0(t-{\mathrm{kT}}_s)---\left(5\right)$

其中

$c_0={\mathrm{\Σ}}_{k=-\∞}^{\∞}{j}^{k-1}{\mathrm{\α}}_k{\mathrm{\α}}_{k-1}---\left(6\right)$

还能够区别地对用于传输的数据进行编码(即α_k＝d_kα_k-1，其中d_k∈{-1，+1}是信息序列)。这暗示d_k＝α_kα_k-1。下式的序列假定已应用了差分编码。然而，差分编码不是必需的，并且可以同样地在没有差分编码的情况下实现本文所描述的方法和装置。

接收信号通常由和h₀(t)匹配的滤波器滤波以将以下输入提供给解调器：

$r\left(t\right)={\mathrm{\Σ}}_{k=-\∞}^{\∞}{j}^{k-1}{d}_k{h}_M(t-{\mathrm{kT}}_s)---\left(7\right)$

式7表示理想接收信号的近似。可以从式7看到，如果以符号间隔对理想信号的这个近似进行采样，则那些符号将在实的与虚的之间交替。在实系统中，不仅近似很可能不是完全准确的而且信号将已在传输期间经受噪声和其它劣化。因此，对信号进行采样不太可能给出全实样本或全虚样本，而是替代地为显著地实的或显著地虚的样本。

与诸如BPSK的一些其它调制方案不同，匹配滤波器响应h_M(t)不是奈奎斯特的。因此，存在符号间干扰，这意味着需要某种形式的均衡化以获得输入数据d_k的最大似然估计。适合的均衡器的示例是维特比均衡器。

现在将描述解码方法的具体示例。

接收信号r(t)表示以规则间隔T_s隔开的多个数据符号d_k。该方法通过使用匹配滤波器对信号进行滤波而适当地开始。能够通过离散的三抽头滤波器近似匹配滤波器h_M(t)。抽头优选地在符号间距T_s处。对于时间带宽乘积为0.3的示例，抽头可以被加权为1/2、1、1/2。对于较高的时间带宽乘积，能够做出准确的三抽头近似，并且如果权重是对称的则能够应用以上所述的星座算式。然而，在结果得到的ISI大小方面将存在差异，因为ISI抽头对于较高的时间带宽乘积优选地小于1/2。能够对较低的时间带宽乘积应用近似。然而由于三抽头近似对于较低的时间带宽乘积不太准确(即，ISI分量在三个抽头外面较大)，可能观测到性能的损失。

基于理想接收信号(所以没有加入的噪声或其它劣化)的滤波器的输出能够被写为：

$\begin{array}{c}r\left(t\right)=\\ \frac{1}{2}{\mathrm{\Σ}}_{k=-\∞}^{\∞}{j}^{k-2}{d}_{k-1}\mathrm{\δ}(t-(k-1)T_s)+{\mathrm{\Σ}}_{k=-\∞}^{\∞}{j}^{k-1}{d}_k\mathrm{\δ}(t-{\mathrm{kT}}_s)+\frac{1}{2}{\mathrm{\Σ}}_{k=-\∞}^{\∞}{j}^k{d}_{k-2}\mathrm{\δ}(t-\\ (k+1)T_s)\end{array}---\left(8\right)$

经滤波的信号因此是符号及其近邻居的加权和。这是被输入到解调器中的信号。

以每符号间隔T_s对r(t)进行采样给出信号：

$r\left[m\right]=\frac{1}{2}{j}^{m-2}{d}_{m-1}+j^m{d}_m+\frac{1}{2}{j}^m{d}_{m+1}---\left(9\right)$

并且最终乘以具有等于-1/4T_s的频率的复正弦函数，以给出：

$\begin{array}{c}\stackrel{\‾}{r}\left[m\right]=j^{-(m-1)}r\left[m\right]\\ =\frac{1}{2}{j}^{-1}{d}_{m-1}+d_m+\frac{1}{2}{j}^1{d}_{m+1}\\ =d_m+\frac{1}{2}{j}^{-1}{d}_{m-1}+\frac{1}{2}{j}^1{d}_{m+1}\\ =d_m-\frac{j}{2}{d}_{m-1}+\frac{j}{2}{d}_{m+1}\\ =d_m-\frac{j}{2}{d}_{m-1}+\frac{j}{2}{d}_{m+1}\end{array}---\left(10\right)$

复正弦函数的意义是它在各个符号间隔具有为全实的或全虚的值。因为经滤波的信号的样本r[m]在实与虚之间交替(对于没有加入的噪声或其它信号劣化的理想情况)，将经采样滤波的信号乘以复正弦函数将信号及其符号间干扰在复数乘积的实部与虚部之间分割。复正弦函数可以是任何适合的频率，但是优选地它的周期是符号间隔的四倍的整数倍。

式10的实部表示数据符号。虚部表示符号间干扰。在实践中，实部和虚部这二者将以与BPSK信号相似的方式经受噪声。数据的软估计能够被记为并且这能够以与其它调制方案(诸如BPSK)中的数据估计相同的方式用在后续阶段中。符号间干扰能够被记为并且这能够被用来适配解码过程以去除符号间干扰。

应当理解在其它示例中，式10的虚部可以表示数据符号而实部可以表示符号间干扰：至于数据符号乘以j的奇数幂还是j的偶数幂(并且类似地对于符号间干扰)，它取决于特定复正弦函数与特定采样滤波的信号之间的对应关系。

另一重要点是符号间干扰仅能够取特定值。在式10的示例中，这些值是j、0、-j。因此完整信号星座是已知的。在图2a中示出了没有噪声的理想星座。在图2b中示出了有噪声的星座的示例。信号星座能够例如通过监测接收到的符号与它们在解码过程期间被映射到的星座点的偏差而被用于信道跟踪。在这个上下文中，信道跟踪包括跟踪由发送机和接收机实现方式所引入的多普勒漂移和相位噪声或频率偏移这二者。

维特比均衡化需要分支度量的计算然后向后追踪以确定最可能的状态。对于低带宽GMSK信号，信道能够相对于符号周期相当迅速地改变并且频率误差还可能漂移。这导致随着分支度量正被计算而潜在显著的信道改变。然后可能在误差中发生向后追踪操作。另外，判定定向跟踪因为由向后追踪导致的延迟而是不准确的。

利用所导出的算式，信号的实部包含能够被即时估计的信息。这允许信道估计的即时更新。还能够通过计算期望符号(包括复分量)并且将它与所接收到的值相比较来估计信道或频率值。

在使用本文所描述的方法和装置时的优点是它跟踪信道和频率改变的低复杂性能力。在图3中示出了一个可能的跟踪方案的示例，将参照图所示的接收机结构对图3进行描述。

接收机结构包括匹配滤波器401、乘法器402和均衡器403(例如，单抽头MMSE均衡器)。该接收机还包括用于适当地识别被乘信号的实部或虚部并且将其作为软数据符号输出的判定单元404。该接收机结构还包括反馈回路以用于将信号反馈给均衡器，从而使得它能够适于变化的信道条件和多径。反馈回路包括重建单元405、除法器406和跟踪回路(407)。在下面参照图3描述这三个单元的操作。

跟踪方案由图4中的虚线表示的三个主要阶段(I至III)组成。在对符号进行解码(步骤301)之后，在阶段I中使用影响这个符号的三个判定来计算期望星座符号(步骤302)。这个步骤基本上取得经解码的符号并且估计假若所接收到的符号和经解码的符号相同则乘法器的输出将象什么。假定符号已被正确地解码，期望符号与实际接收到的符号之间的任何差是由于传输期间的失真而导致的。两个值的比较因此提供信道条件的指示。(注意图4示出了因果系统，所以期望星座符号比判定滞后了一。这仅用于示例的目的。将本文所描述的原理扩展到非因果系统是容易的。)

通过将期望符号划分成接收信号来计算信道响应(步骤303；阶段II)。计算其与当前跟踪的信道之间的差(步骤304)。该差然后被用作进入阶段III的自适应跟踪回路的误差信号，以将信道估计更新为其新值(步骤305)。也就是说，在时间n的信道抽头估计a_n被用来通过a_n+1＝a_n+μe_n来更新在时间n+1的信道估计，其中μ是控制跟踪速度的跟踪参数并且e_n是在时间n的误差。最终在均衡器中使用新的信道值(步骤306)，所述均衡器可以是例如单抽头MMSE均衡器(例如，h*/(|h|²+σ²)，其中σ²是噪声方差)或任何其它均衡器。

能够容易地通过引入新的星座算式来适配被用于其它星座类型的其它信道和频率跟踪方案(BPSK、QPSK等)。

图4示出了具体接收机结构。这仅用于例示性目的。图中所示的各个结构对应于可能由任何适合的功能单元、部件或这些部件的合集执行的功能。这样的功能可能用硬件或软件或这二者的组合加以实现。图4所示的结构不旨在定义在芯片上的硬件的不同部分之间或在软件中的不同程序、过程或函数之间的严格划分。在一些实施方式中，本文所描述的算法中的一些或全部可以全部或部分地用硬件加以执行。在许多实施方式中，算法中的至少一部分可以由在软件控制下行动的处理器(例如通信设备的CPU或DSP)来实现。任何这样的软件被优选地存储在非暂时计算机可读介质上，所述非暂时性计算机可读介质诸如存储器(RAM、高速缓存、硬盘等)或其它存储装置(USB棒、CD、磁盘等)。

在大多数实施方式中接收机结构将形成更大通信设备的一部分。示例包括M2M设备、移动电话、智能电话、线路连接电话、膝上型电脑、平板等。典型的通信设备包括天线、CPU、存储器、信号处理电路，诸如DSP和滤波器等。

本文所描述的方法可以应用于针对物联网(IoT)通信构造的通信网络。示例将包括被构造为根据Weightless^TM协议操作的网络(但是本文所描述的方法可以由被构造为根据不同的协议(诸如例如LTE、蓝牙、WiFi、VoIP)操作的网络容易地实现)。通常，网络将包括被各自构造为与大量地理上隔开的终端进行通信的许多通信设备(例如基站)。网络可以是蜂窝网络，其中各个通信设备负责与位于相应小区中的终端的空中通信。所描述的方法对于其中由GMSK信号使用的信道带宽相对低例如小于100kHz的通信系统可能是特别有利的。

在一个示例中，接收机结构可以被构造为根据Weightless^TM IoT规范操作。Weightless^TM使用蜂窝WAN架构，其中针对IoT系统的要求(低终端成本、低终端占空度以及因此低功耗和对非常低数据速率的可伸缩性)优化了协议。它原先被设计为在从470MHz到790MHz的TV空白频谱中操作，但是PHY被一般化为在变化带宽的许可频带、共享许可接入频带和免许可频带中操作。

本申请人因此孤立地公开本文所描述的各个单独的特征以及两个或更多个这样的特征的任何组合，在这个意义上这样的特征或组合能够总体上鉴于本领域技术人员的公共一般知识基于本说明书被执行，而不管这样的特征或这些特征的组合是否解决本文所公开的任何问题，并且不限于权利要求的范围。本申请人指示本发明的方面可以包括任何这样的单独特征或这些特征的组合。鉴于上述描述，对于本领域技术人员将显然的是，可以在本发明的范围内做出各种修改。

Claims (14)

1.一种用于对遭受符号间干扰的接收信号中的符号进行解码的方法，该接收信号包括由符号间隔隔开的多个符号，该方法包括以下步骤：

对所述接收信号进行滤波以形成表示所述符号的样本，该样本包括所述符号和在所述接收信号中紧先于该符号以及接续该符号的符号的加权和；

将所述样本乘以在所述符号间隔的整数倍处交替地具有以下各项的复正弦函数，以生成复数值：(i)零实分量和非零虚分量以及(ii)非零实分量和零虚分量；以及

依赖于所述复数值对所述符号进行解码。

2.根据权利要求1中所述的方法，其中，所述接收信号是这样的，即所述样本包括为显著地实的或显著地虚的符号以及在所述接收信号中紧先于该符号以及接续该符号的显著地为实的或虚的中的另一方的符号的加权和，所述方法包括以下步骤：

如果所述符号的加权分量是显著地实的，并且所述复正弦函数在对应于所述样本的所述符号间隔的整数倍处具有零实分量和非零虚分量，则依赖于所述复数值的虚部对所述符号进行解码；

如果所述符号的加权分量是显著地实的，并且所述复正弦函数在对应于所述样本的所述符号间隔的整数倍处具有非零实分量和零虚分量，则依赖于所述复数值的实部对所述符号进行解码；

如果所述符号的加权分量是显著地虚的，并且所述复正弦函数在对应于所述样本的所述符号间隔的整数倍处具有零实分量和非零虚分量，则依赖于所述复数值的实部对所述符号进行解码；以及

如果所述符号的加权分量是显著地虚的，并且所述复正弦函数在对应于所述样本的所述符号间隔的整数倍处具有非零实分量和零虚分量，则依赖于所述复数值的虚部对所述符号进行解码。

3.根据权利要求1或2中所述的方法，该方法包括以下步骤：依赖于所述复数值的实分量和虚分量中的一方对所述符号进行解码并且依赖于所述复数值的实分量和虚分量中的另一方来估计所述符号间干扰。

4.根据前述权利要求中任一项所述的方法，其中，所述复正弦函数的一个周期等于所述符号间隔的四倍的整数倍。

5.根据前述权利要求中任一项所述的方法，其中，所述调制方法是GMSK。

6.根据前述权利要求中任一项所述的方法，其中，所述GMSK调制的时间带宽乘积大于或等于0.3。

7.根据前述权利要求中任一项所述的方法，该方法包括以下步骤：

依赖于所述复正弦函数来生成期望符号的星座；

监测所述多个符号与期望符号的星座之间的偏差；以及

依赖于所监测到的偏差来跟踪通过其接收到所述信号的信道的改变。

8.根据前述权利要求中任一项所述的方法，该方法包括：生成期望符号，该期望符号表示在假设所接收到的符号和经解码的符号相同的情况下将生成的复数值。

9.根据权利要求6中所述的方法，该方法包括：通过将经解码的符号与从在所述接收信号中紧先于所接收到的符号以及接续所接收到的符号的符号生成的解码的符号合计来生成所述期望符号。

10.根据权利要求7中所述的方法，该方法包括以下步骤：

将所述期望符号除以所接收到的符号；

将信道估计减去除法的输出；

将已减信号提供给自适应滤波器作为所述信道估计的误差的表示，该自适应滤波器利用该表示来处理所述接收信号；以及

依赖于已减信号来更新所述信道估计。

11.一种用于被构造为接收遭受符号间干扰的信号的接收机的解码器，所述接收信号包括以规则间隔隔开的多个符号，该解码器包括：

滤波器，该滤波器被构造为对所述接收信号进行滤波以形成表示所述符号的样本，该样本包括所述符号和在所述接收信号中紧先于该符号以及接续该符号的符号的加权和；

乘法器，该乘法器被构造为将所述样本乘以在符号间隔的整数倍处交替地具有以下各项的复正弦函数，以生成复数值：(i)零实分量和非零虚分量以及(ii)非零实分量和零虚分量；以及

判定单元，该判定单元被构造为依赖于所述复数值对所述符号进行解码。

12.一种被构造为操作根据权利要求1至7中的任一项所述的方法的装置。

13.一种基本上如参照附图本文所描述的方法。

14.一种基本上如参照附图本文所描述的解码器。