CN109672501B - 正交调制的软判决度量方法和装置 - Google Patents

正交调制的软判决度量方法和装置 Download PDF

Info

Publication number
CN109672501B
CN109672501B CN201811579456.2A CN201811579456A CN109672501B CN 109672501 B CN109672501 B CN 109672501B CN 201811579456 A CN201811579456 A CN 201811579456A CN 109672501 B CN109672501 B CN 109672501B
Authority
CN
China
Prior art keywords
matrix
value
signal sequence
row
bit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811579456.2A
Other languages
English (en)
Other versions
CN109672501A (zh
Inventor
王纪东
丁昊成
游行远
徐彬彬
付睿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
722th Research Institute of CSIC
Original Assignee
722th Research Institute of CSIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 722th Research Institute of CSIC filed Critical 722th Research Institute of CSIC
Priority to CN201811579456.2A priority Critical patent/CN109672501B/zh
Publication of CN109672501A publication Critical patent/CN109672501A/zh
Application granted granted Critical
Publication of CN109672501B publication Critical patent/CN109672501B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/004Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045Arrangements at the receiver end
    • H04L1/0054Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03312Arrangements specific to the provision of output signals
    • H04L25/03318Provision of soft decisions

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Power Engineering (AREA)
  • Error Detection And Correction (AREA)

Abstract

本发明公开了一种正交调制的软判决度量方法和装置,属于通信领域。所述方法包括:获取待判决信号序列,待判决信号序列为完成前端处理的信号序列,待判决信号序列的长度为L;将待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值序列,L阶标准矩阵是根据正交调制中约定的映射方式由调制过程中的正交调制序列矩阵映射得到的;采用WZ矩阵和WO矩阵中各行的值作为坐标将相关值映射为2个矩阵,得到Z矩阵和O矩阵,WZ矩阵和WO矩阵分别用于记录正交调制序列矩阵中值为0和1的比特的序号;根据Z矩阵和O矩阵计算待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值。

Description

正交调制的软判决度量方法和装置
技术领域
本发明涉及通信领域,特别涉及一种正交调制的软判决度量方法和装置。
背景技术
随着现代通信技术的发展,在无线通信系统中对频谱利用率的要求越来越高,因此通信过程中采用多址技术以满足需求。其中,码分多址指的是利用各不相同的、相互(准)正交的地址码调制所发送的信号的技术。在接收端利用码型的(准)正交性,通过地址识别(相关检测)从混合信号中选出相应的信号。
常见的码分多址译码包括硬判决译码和软判决译码。在软判决译码中,需要先根据软判决度量生成软判决译码的输入值,然后通过判决算法进行译码得到编码前的数据。传统的软判决度量生成算法多是基于最小欧氏距离进行计算,是一种基于最大似然准则的解调方式,其计算复杂度与调制阶数有关。
但在计算最小欧氏距离时,需要求复数模的平方,运算过程具有较高的复杂度,硬件开销较大。
发明内容
本发明实施例提供了一种正交调制的软判决度量方法和装置,降低软判决度量运算过程的复杂度,减小硬件开销。所述技术方案如下:
一方面,本发明实施例提供了一种正交调制的软判决度量方法,所述方法包括:
获取待判决信号序列,所述待判决信号序列为完成前端处理的信号序列,所述待判决信号序列的长度为L;
将所述待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值序列,所述L阶标准矩阵是根据正交调制中约定的映射方式由调制过程中的正交调制序列矩阵映射得到的;
采用WZ矩阵和WO矩阵中各行的值作为坐标将所述相关值映射为2个矩阵,得到Z矩阵和O矩阵,所述WZ矩阵和WO矩阵分别用于记录所述正交调制序列矩阵中值为0和1的比特的序号;
根据所述Z矩阵和O矩阵计算所述待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值。
在本发明实施例的一种实现方式中,所述将所述待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值,包括:
按照如下公式将所述待判决信号序列与L阶标准矩阵的第l行做相关,得到第l个取值:
Figure BDA0001917497730000021
其中,xM表示所述待判决信号序列,xM(i)表示xM的第i个元素,0≤i≤L-1;W为所述L阶标准矩阵,Wl表示W矩阵的第l行,Wl(i)为W矩阵的第l行的第i个元素,Wl(i)*表示Wl(i)的共轭值,0≤l≤L-1;
将所述第l个取值作为所述相关值序列的第l位。
在本发明实施例的一种实现方式中,所述采用WZ矩阵和WO矩阵中各行的值作为坐标将所述相关值映射为2个矩阵,得到Z矩阵和O矩阵,包括:
对于所述待判决信号序列的每个比特均按照如下方式进行处理:
以所述待判决信号序列的第n个比特的位置为行坐标,取出所述WZ矩阵和WO矩阵中第n行元素的值,0≤n≤L0-1;
分别以所述WZ矩阵和WO矩阵中第n行元素的值作为坐标取出所述相关值序列中的值,分别放入所述Z矩阵和所述O矩阵中以所述第n个比特的位置为行坐标的行。
在本发明实施例的一种实现方式中,所述根据所述Z矩阵和O矩阵计算所述待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值,包括:
按照如下公式计算所述待判决信号序列在编码前的每个比特的对数似然比:
LLR(n)=max[O(n)]-max[Z(n)];
其中,LLR(n)为所述待判决信号序列在编码前的第n比特的对数似然比,0≤n≤L0-1;max[O(n)]为所述O矩阵的第n行的最大值,max[Z(n)]为所述Z矩阵的第n行的最大值;
按照如下公式计算所述待判决信号序列在编码前的每个比特的软判决度量的度量值:
Figure BDA0001917497730000031
其中,scale(n)为所述待判决信号序列在编码前的第n比特的软判决度量的度量值。
在本发明实施例的一种实现方式中,所述方法还包括:
按照如下方式得到所述WZ矩阵和WO矩阵:
记录所述正交调制序列矩阵的第k列中为值为0的比特的索引值,写入所述WZ矩阵的第k行,0≤k≤L0-1;
记录所述正交调制序列矩阵的第k列中为值为1的比特的索引值,写入所述WO矩阵的第k行。
另一方面,本发明实施例提供了一种正交调制的软判决度量装置,所述装置包括:
获取模块,用于获取待判决信号序列,所述待判决信号序列为完成前端处理的信号序列,所述待判决信号序列的长度为L;
相关模块,用于将所述待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值序列,所述L阶标准矩阵是根据正交调制中约定的映射方式由调制过程中的正交调制序列矩阵映射得到的;
映射模块,用于采用WZ矩阵和WO矩阵中各行的值作为坐标将所述相关值映射为2个矩阵,得到Z矩阵和O矩阵,所述WZ矩阵和WO矩阵分别用于记录所述正交调制序列矩阵中值为0和1的比特的序号;
处理模块,用于根据所述Z矩阵和O矩阵计算所述待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值。
在本发明实施例的一种实现方式中,所述相关模块,用于按照如下公式将所述待判决信号序列与L阶标准矩阵的第l行做相关,得到第l个取值:
Figure BDA0001917497730000041
其中,xM表示所述待判决信号序列,xM(i)表示xM的第i个元素,0≤i≤L-1;W为所述L阶标准矩阵,Wl表示W矩阵的第l行,Wl(i)为W矩阵的第l行的第i个元素,Wl(i)*表示Wl(i)的共轭值,0≤l≤L-1;
将所述第l个取值作为所述相关值序列的第l位。
在本发明实施例的一种实现方式中,所述映射模块,用于对于所述待判决信号序列的每个比特均按照如下方式进行处理:
以所述待判决信号序列的第n个比特的位置为行坐标,取出所述WZ矩阵和WO矩阵中第n行元素的值,0≤n≤L0-1;
分别以所述WZ矩阵和WO矩阵中第n行元素的值作为坐标取出所述相关值序列中的值,分别放入所述Z矩阵和所述O矩阵中以所述第n个比特的位置为行坐标的行。
在本发明实施例的一种实现方式中,所述处理模块,用于按照如下公式计算所述待判决信号序列在编码前的每个比特的对数似然比:
LLR(n)=max[O(n)]-max[Z(n)];
其中,LLR(n)为所述待判决信号序列在编码前的第n比特的对数似然比,0≤n≤L0-1;max[O(n)]为所述O矩阵的第n行的最大值,max[Z(n)]为所述Z矩阵的第n行的最大值;
所述处理模块,还用于按照如下公式计算所述待判决信号序列在编码前的每个比特的软判决度量的度量值:
Figure BDA0001917497730000042
其中,scale(n)为所述待判决信号序列在编码前的第n比特的软判决度量的度量值。
在本发明实施例的一种实现方式中,所述映射模块,还用于按照如下方式得到所述WZ矩阵和WO矩阵:
记录所述正交调制序列矩阵的第k列中为值为0的比特的索引值,写入所述WZ矩阵的第k行,0≤k≤L0-1;
记录所述正交调制序列矩阵的第k列中为值为1的比特的索引值,写入所述WO矩阵的第k行。
本发明实施例提供的技术方案带来的有益效果是:
通过在求出待判决信号序列与L阶标准矩阵的相关值的基础上,计算待判决信号序列在编码前的每个比特的对数似然比,采用归一化后的对数似然比作为软判决度量的度量值,这种方案能够简化软判决度量的计算及输出,大幅度减少了对复数序列取模值和开平方计算的次数,从而降低了运算的复杂度,能够较大程度上减少软判决过程中硬件的开销,且不会使性能有明显下降。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种正交调制的软判决度量方法的流程图;
图2是本发明实施例提供的另一种正交调制的软判决度量方法的流程图;
图3是本发明实施例提供的一种WZ矩阵和WO矩阵形成示意图;
图4是本发明实施例提供的信噪比为0-10dB环境下的误码率曲线图;
图5是本发明实施例提供的一种正交调制的软判决度量装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种正交调制的软判决度量方法的流程图。参见图1,该方法包括:
步骤101:获取待判决信号序列。
其中,待判决信号序列为完成前端处理的信号序列,待判决信号序列的长度为L。
这里,前端处理包括信号接收端在进行判决译码前的处理过程,包括模数转换、子载波解调、同步等处理。本发明实施例提供的软判决度量方法不需要特定的硬件设备,采用数字信号处理器(Digital Signal Processor,DSP)等器件即可执行。
步骤102:将待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值序列。
其中,L阶标准矩阵是根据正交调制中约定的映射方式由调制过程中的正交调制序列矩阵映射得到的。
例如,在编码端每组待编码数据比特的长度为L0,记为
Figure BDA0001917497730000061
每组数据比特将映射成为一组长度为L的正交调制符号,也就是前述正交调制序列矩阵,正交调制序列矩阵的尺寸为L行L0列,
Figure BDA0001917497730000062
并且该正交调制序列矩阵是在确定编码方式后收发双方已知的,并且在通信过程中是固定不变的。根据正交调制的标准编码和编码过程中约定的映射方式形成用于后续度量的L阶标准矩阵WL*L[·]。比如,以8移相键控(8Phase ShiftKeying,8PSK)为例,在调制过程中选择将[0]值映射为常数a0,[4]值映射为常数a1;则按此映射方式,将正交调制序列矩阵映射为L阶的标准矩阵W时,采用[0]值代替正交调制序列矩阵中的常数a0,采用[4]值代替正交调制序列矩阵中的常数a1,得到L阶的标准矩阵W。
步骤103:采用WZ矩阵和WO矩阵中各行的值作为坐标将相关值映射为2个矩阵,得到Z矩阵和O矩阵。
其中,WZ矩阵和WO矩阵分别用于记录正交调制序列矩阵中值为0和1的比特的序号。
具体地,WZ矩阵和WO矩阵均为L0行L/2列的矩阵,WZ矩阵的各行依次记录正交调制序列矩阵中各列里值为0的比特的序号,WO矩阵的各行依次记录正交调制序列矩阵中各列里值为0的比特的序号。
步骤104:根据Z矩阵和O矩阵计算待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值。
这里,每个比特的对数似然比用来表示该比特在判决中判断为0或1的概率。由于软判决度量的度量值表示的也是判决中将一位比特判断为0或1的概率,或者称为比特对应的数据点在坐标空间中距离0点和1点的距离,因此,采用每个比特的对数似然比可以表示该比特的软判决度量的度量值。
在该步骤中,归一化是指将对数似然比化为(0,1)之间的小数。因为在判决译码时,维特比(viterbi)译码器的输入要求是(-1,1)之间的值,通过归一化使得本方案提供的软判决度量的度量值可以满足后续译码要求。
本发明提供了一种适用于正交调制的软判决度量方法,能够在求出待判决信号序列与L阶标准矩阵的相关值的基础上,计算待判决信号序列在编码前的每个比特的对数似然比,采用归一化后的对数似然比作为软判决度量的度量值,这种方案能够简化软判决度量的计算及输出,大幅度减少了对复数序列取模值和开平方计算的次数,从而降低了运算的复杂度,能够较大程度上减少软判决过程中硬件的开销,且不会使性能有明显下降。
图2是本发明实施例提供的另一种正交调制的软判决度量方法的流程图。
参见图2,该方法包括:
步骤200:根据正交调制中约定的映射方式,将调制过程中的正交调制序列矩阵映射得到L阶标准矩阵。
例如,在编码端每组待编码数据比特的长度为L0,记为
Figure BDA0001917497730000071
每组数据比特将映射成为一组长度为L的正交调制符号,也就是前述正交调制序列矩阵,正交调制序列矩阵的尺寸为L行L0列,
Figure BDA0001917497730000072
并且该正交调制序列矩阵是在确定编码方式后收发双方已知的,并且在通信过程中是固定不变的。根据正交调制的标准编码和编码过程中约定的映射方式形成用于后续度量的L阶标准矩阵WL*L[·]。比如,以8PSK为例,在调制过程中选择将[0]值映射为常数a0,[4]值映射为常数a1;则按此映射方式,将正交调制序列矩阵映射为L阶的标准矩阵W时,采用[0]值代替正交调制序列矩阵中的常数a0,采用[4]值代替正交调制序列矩阵中的常数a1,得到L阶的标准矩阵W。
步骤201:根据正交调制序列矩阵得到WZ矩阵和WO矩阵。
其中,WZ矩阵和WO矩阵分别用于记录正交调制序列矩阵中值为0和1的比特的序号。
具体地,该步骤可以包括:
按照如下方式得到WZ矩阵和WO矩阵:
记录正交调制序列矩阵的第k列中为值为0的比特的索引值(位置),写入WZ矩阵的第k行,0≤k≤L0-1;
记录正交调制序列矩阵的第k列中为值为1的比特的索引值,写入WO矩阵的第k行。
按照上述方式,得到两个L0行L/2列的WZ矩阵和WO矩阵。
图3是本发明实施例提供的一种WZ矩阵和WO矩阵形成示意图。参见图3,待编码数据比特的每组长度为4,每组比特为b3b2b1b0
按照前述方式,在第0列(b3列)中,第0~7个比特的取值为0,第8~15个比特的取值为1,则WZ矩阵的第0行(图中WZ[0][·])为0~7,WO矩阵的第0行(图中WO[0][·])为8~15。后续各列,也按照此方式得到,得到的矩阵为4行8列(L0行L/2列)。
步骤202:获取待判决信号序列。
其中,待判决信号序列为完成前端处理的信号序列,待判决信号序列的长度为L。
这里,前端处理包括信号接收端在进行判决译码前的处理过程,包括模数转换、子载波解调、同步等处理。本发明实施例提供的软判决度量方法不需要特定的硬件设备,采用DSP等器件即可执行。
步骤203:将待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值序列。
在本发明实施例中,可以按照如下公式将待判决信号序列与L阶标准矩阵的第l行做相关,得到第l个取值:
Figure BDA0001917497730000081
其中,xM表示待判决信号序列,xM(i)表示xM的第i个元素,0≤i≤L-1;W为L阶标准矩阵,Wl表示W矩阵的第l行,Wl(i)为W矩阵的第l行的第i个元素,Wl(i)*表示Wl(i)的共轭值,0≤l≤L-1;
将第l个取值作为相关值序列的第l位。
也即,可以按照将待判决信号序列依次与L阶标准矩阵的各行依次相关,得到L个取值,将各个取值分别放入相关值序列的对应位置,得到长度为L的相关值序列。
步骤204:采用WZ矩阵和WO矩阵中各行的值作为坐标将相关值映射为2个矩阵,得到Z矩阵和O矩阵。
在本发明实施例中,对于待判决信号序列的每个比特均按照如下方式进行处理:
以待判决信号序列的第n个比特的位置为行坐标,取出WZ矩阵和WO矩阵中第n行元素的值,0≤n≤L0-1;其中,第0个比特的位置为0,则对应的行坐标为0,也即WZ矩阵和WO矩阵中的第0行,第1个比特的位置为1,则对应的行坐标为1,也即WZ矩阵和WO矩阵中的第1行,依次类推。
分别以WZ矩阵和WO矩阵中第n行元素的值作为坐标取出相关值序列中的值,分别放入Z矩阵和O矩阵中以第n个比特的位置为行坐标的行。也即,将WZ矩阵第n行的各个元素的值作为坐标取出相关值序列中的值,分别放入Z矩阵第n行的各个位置,将WO矩阵第n行的各个元素的值作为坐标取出相关值序列中的值,分别放入O矩阵第n行的各个位置。例如,WZ矩阵第0行的第0元素(第0列)的值为0,则取出相关值序列中的第0个元素的值,放入Z矩阵第0行的第0列。
在本发明实施例中,Z矩阵和O矩阵的尺寸与均为L0行L/2列。
上述过程可以采用如下公式表示:
Figure BDA0001917497730000091
其中,rM为相关值序列,WZ(n,i)为WZ矩阵第n行第i个元素,Z(n,i)为Z矩阵第n行第i个元素,WO(n,i)为WO矩阵第n行第i个元素,O(n,i)为O矩阵第n行第i个元素。
步骤205:根据Z矩阵和O矩阵计算待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值。
这里,每个比特的对数似然比用来表示该比特在判决中判断为0或1的概率。由于软判决度量的度量值表示的也是判决中将一位比特判断为0或1的概率,或者称为比特对应的数据点在坐标空间中距离0点和1点的距离,因此,采用每个比特的对数似然比可以表示该比特的软判决度量的度量值。
在该步骤中,归一化是指将对数似然比化为(0,1)之间的小数。因为在判决译码时,viterbi译码器的输入要求是(-1,1)之间的值,通过归一化使得本方案提供的软判决度量的度量值可以满足后续译码要求。
在本发明实施例中,可以按照如下公式计算待判决信号序列在编码前的每个比特的对数似然比:
Figure BDA0001917497730000101
其中,LLR(n)为待判决信号序列在编码前的第n比特的对数似然比,编码前的比特一共L0个,也即0≤n≤L0-1;exp(Z(n))表示e为底Z矩阵的第n行的各个取值为指数的指数函数,
Figure BDA0001917497730000102
表示该指数函数在各个取值时求和,0≤i≤L/2-1。exp(O(n))表示e为底O矩阵的第n行的各个取值为指数的指数函数,
Figure BDA0001917497730000103
表示该指数函数在各个取值时求和。
exp(Z(n))、exp(O(n))物理意义上表示一种度量,意味着相关值越大,这个比特的点距离0点越近,判决为0的概率越大。
exp是指数运算,对于指数的大小变化比较敏感并且是单调增函数,当Z(n)和O(n)的值较小时其指数函数趋近于0,因此实际上这里的求和值可以认为约等于max[exp(Z(n))]、max[exp(O(n))],即:
Figure BDA0001917497730000104
其中,max[O(n)]为O矩阵的第n行的最大值,max[Z(n)]为Z矩阵的第n行的最大值。
按照如下公式计算待判决信号序列在编码前的每个比特的软判决度量的度量值:
Figure BDA0001917497730000105
其中,scale(n)为待判决信号序列在编码前的第n比特的软判决度量的度量值。
在得到度量值scale(i)可直接送入解交织器或viterbi译码器进行后续处理,后续处理指的是在送入译码器后按viterbi法进行解码,相当于进行逆运算得到编码前的信息。
图4为采用本发明实施例提供的正交调制的软判决度量方法后,信噪比为0-10dB环境下的误码率曲线图。在图4中,纵坐标BER,表示误码率;横坐标Eb/n0,表示归一化信噪比。参见图4,可以看出在不同的输入归一化信噪比下解码时的误码率变化,显然在信噪比越高的时候译码过程中的错误越少,误码率越低,符合理论规律,因此本发明提供的软判决度量方法有效。
图5是本发明实施例提供的一种正交调制的软判决度量装置的结构示意图。
参见图5,该装置包括:
获取模块301,用于获取待判决信号序列,待判决信号序列为完成前端处理的信号序列,待判决信号序列的长度为L;
相关模块302,用于将待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值序列,L阶标准矩阵是根据正交调制中约定的映射方式由调制过程中的正交调制序列矩阵映射得到的;
映射模块303,用于采用WZ矩阵和WO矩阵中各行的值作为坐标将相关值映射为2个矩阵,得到Z矩阵和O矩阵,WZ矩阵和WO矩阵分别用于记录正交调制序列矩阵中值为0和1的比特的序号;
处理模块304,用于根据Z矩阵和O矩阵计算待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值。
在本发明实施例的一种实现方式中,相关模块302,用于按照如下公式将待判决信号序列与L阶标准矩阵的第l行做相关,得到第l个取值:
Figure BDA0001917497730000111
其中,xM表示待判决信号序列,xM(i)表示xM的第i个元素,0≤i≤L-1;W为L阶标准矩阵,Wl表示W矩阵的第l行,Wl(i)为W矩阵的第l行的第i个元素,Wl(i)*表示Wl(i)的共轭值,0≤l≤L-1;
将第l个取值作为相关值序列的第l位。
在本发明实施例的一种实现方式中,映射模块303,用于对于待判决信号序列的每个比特均按照如下方式进行处理:
以待判决信号序列的第n个比特的位置为行坐标,取出WZ矩阵和WO矩阵中第n行元素的值,0≤n≤L0-1;
分别以WZ矩阵和WO矩阵中第n行元素的值作为坐标取出相关值序列中的值,分别放入Z矩阵和O矩阵中以第n个比特的位置为行坐标的行。
在本发明实施例的一种实现方式中,处理模块304,用于按照如下公式计算待判决信号序列在编码前的每个比特的对数似然比:
LLR(n)=max[O(n)]-max[Z(n)];
其中,LLR(n)为待判决信号序列在编码前的第n比特的对数似然比,0≤n≤L0-1;max[O(n)]为O矩阵的第n行的最大值,max[Z(n)]为Z矩阵的第n行的最大值;
处理模块304,还用于按照如下公式计算待判决信号序列在编码前的每个比特的软判决度量的度量值:
Figure BDA0001917497730000121
其中,scale(n)为待判决信号序列在编码前的第n比特的软判决度量的度量值。
在本发明实施例的一种实现方式中,映射模块303,还用于按照如下方式得到WZ矩阵和WO矩阵:
记录正交调制序列矩阵的第k列中为值为0的比特的索引值,写入WZ矩阵的第k行,0≤k≤L0-1;
记录正交调制序列矩阵的第k列中为值为1的比特的索引值,写入WO矩阵的第k行。
需要说明的是:上述实施例提供的正交调制的软判决度量装置在计算正交调制的软判决度量时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的正交调制的软判决度量装置与正交调制的软判决度量方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种正交调制的软判决度量方法,其特征在于,所述方法包括:
获取待判决信号序列,所述待判决信号序列为完成前端处理的信号序列,所述待判决信号序列的长度为L;
将所述待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值序列,所述L阶标准矩阵是根据正交调制中约定的映射方式由调制过程中的正交调制序列矩阵映射得到的;
采用WZ矩阵和WO矩阵中各行的值作为坐标将所述相关值映射为2个矩阵,得到Z矩阵和O矩阵,所述WZ矩阵和WO矩阵分别用于记录所述正交调制序列矩阵中值为0和1的比特的序号;
根据所述Z矩阵和O矩阵计算所述待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值;
其中,所述采用WZ矩阵和WO矩阵中各行的值作为坐标将所述相关值映射为2个矩阵,得到Z矩阵和O矩阵,包括:
对于所述待判决信号序列的每个比特均按照如下方式进行处理:
以所述待判决信号序列的第n个比特的位置为行坐标,取出所述WZ矩阵和WO矩阵中第n行元素的值,0≤n≤L0-1;
分别以所述WZ矩阵和WO矩阵中第n行元素的值作为坐标取出所述相关值序列中的值,分别放入所述Z矩阵和所述O矩阵中以所述第n个比特的位置为行坐标的行。
2.根据权利要求1所述的方法,其特征在于,所述将所述待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值,包括:
按照如下公式将所述待判决信号序列与L阶标准矩阵的第l行做相关,得到第l个取值:
Figure FDA0002930360510000011
其中,xM表示所述待判决信号序列,xM(i)表示xM的第i个元素,0≤i≤L-1;W为所述L阶标准矩阵,Wl表示W矩阵的第l行,Wl(i)为W矩阵的第l行的第i个元素,Wl(i)*表示Wl(i)的共轭值,0≤l≤L-1;
将所述第l个取值作为所述相关值序列的第l位。
3.根据权利要求1所述的方法,其特征在于,所述根据所述Z矩阵和O矩阵计算所述待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值,包括:
按照如下公式计算所述待判决信号序列在编码前的每个比特的对数似然比:
LLR(n)=max[O(n)]-max[Z(n)];
其中,LLR(n)为所述待判决信号序列在编码前的第n比特的对数似然比,0≤n≤L0-1;max[O(n)]为所述O矩阵的第n行的最大值,max[Z(n)]为所述Z矩阵的第n行的最大值;
按照如下公式计算所述待判决信号序列在编码前的每个比特的软判决度量的度量值:
Figure FDA0002930360510000021
其中,scale(n)为所述待判决信号序列在编码前的第n比特的软判决度量的度量值。
4.根据权利要求1-3任一项所述的方法,其特征在于,所述方法还包括:
按照如下方式得到所述WZ矩阵和WO矩阵:
记录所述正交调制序列矩阵的第k列中为值为0的比特的索引值,写入所述WZ矩阵的第k行,0≤k≤L0-1;
记录所述正交调制序列矩阵的第k列中为值为1的比特的索引值,写入所述WO矩阵的第k行。
5.一种正交调制的软判决度量装置,其特征在于,所述装置包括:
获取模块,用于获取待判决信号序列,所述待判决信号序列为完成前端处理的信号序列,所述待判决信号序列的长度为L;
相关模块,用于将所述待判决信号序列与L阶标准矩阵的每一行依次做相关,得到一列长度为L的相关值序列,所述L阶标准矩阵是根据正交调制中约定的映射方式由调制过程中的正交调制序列矩阵映射得到的;
映射模块,用于采用WZ矩阵和WO矩阵中各行的值作为坐标将所述相关值映射为2个矩阵,得到Z矩阵和O矩阵,所述WZ矩阵和WO矩阵分别用于记录所述正交调制序列矩阵中值为0和1的比特的序号;
处理模块,用于根据所述Z矩阵和O矩阵计算所述待判决信号序列在编码前的每个比特的对数似然比,并进行归一化作为软判决度量的度量值;
所述映射模块,用于对于所述待判决信号序列的每个比特均按照如下方式进行处理:
以所述待判决信号序列的第n个比特的位置为行坐标,取出所述WZ矩阵和WO矩阵中第n行元素的值,0≤n≤L0-1;
分别以所述WZ矩阵和WO矩阵中第n行元素的值作为坐标取出所述相关值序列中的值,分别放入所述Z矩阵和所述O矩阵中以所述第n个比特的位置为行坐标的行。
6.根据权利要求5所述的装置,其特征在于,所述相关模块,用于按照如下公式将所述待判决信号序列与L阶标准矩阵的第l行做相关,得到第l个取值:
Figure FDA0002930360510000031
其中,xM表示所述待判决信号序列,xM(i)表示xM的第i个元素,0≤i≤L-1;W为所述L阶标准矩阵,Wl表示W矩阵的第l行,Wl(i)为W矩阵的第l行的第i个元素,Wl(i)*表示Wl(i)的共轭值,0≤l≤L-1;
将所述第l个取值作为所述相关值序列的第l位。
7.根据权利要求5所述的装置,其特征在于,所述处理模块,用于按照如下公式计算所述待判决信号序列在编码前的每个比特的对数似然比:
LLR(n)=max[O(n)]-max[Z(n)];
其中,LLR(n)为所述待判决信号序列在编码前的第n比特的对数似然比,0≤n≤L0-1;max[O(n)]为所述O矩阵的第n行的最大值,max[Z(n)]为所述Z矩阵的第n行的最大值;
所述处理模块,还用于按照如下公式计算所述待判决信号序列在编码前的每个比特的软判决度量的度量值:
Figure FDA0002930360510000032
其中,scale(n)为所述待判决信号序列在编码前的第n比特的软判决度量的度量值。
8.根据权利要求5-7任一项所述的装置,其特征在于,所述映射模块,还用于按照如下方式得到所述WZ矩阵和WO矩阵:
记录所述正交调制序列矩阵的第k列中为值为0的比特的索引值,写入所述WZ矩阵的第k行,0≤k≤L0-1;
记录所述正交调制序列矩阵的第k列中为值为1的比特的索引值,写入所述WO矩阵的第k行。
CN201811579456.2A 2018-12-24 2018-12-24 正交调制的软判决度量方法和装置 Active CN109672501B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811579456.2A CN109672501B (zh) 2018-12-24 2018-12-24 正交调制的软判决度量方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811579456.2A CN109672501B (zh) 2018-12-24 2018-12-24 正交调制的软判决度量方法和装置

Publications (2)

Publication Number Publication Date
CN109672501A CN109672501A (zh) 2019-04-23
CN109672501B true CN109672501B (zh) 2021-07-20

Family

ID=66145988

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811579456.2A Active CN109672501B (zh) 2018-12-24 2018-12-24 正交调制的软判决度量方法和装置

Country Status (1)

Country Link
CN (1) CN109672501B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1429427A (zh) * 2000-03-08 2003-07-09 高通股份有限公司 用于计算Turbo解码器的软判决输入度量的方法和装置
CN101217355A (zh) * 2008-01-11 2008-07-09 清华大学 基于低密度奇偶校验译码软判决信息的迭代定时同步方法
CN101232349A (zh) * 2007-12-27 2008-07-30 复旦大学 快速生成qam比特置信度软判决度量的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5326976B2 (ja) * 2009-09-30 2013-10-30 富士通株式会社 無線通信装置、誤り訂正方法及び誤り訂正プログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1429427A (zh) * 2000-03-08 2003-07-09 高通股份有限公司 用于计算Turbo解码器的软判决输入度量的方法和装置
CN101098215A (zh) * 2000-03-08 2008-01-02 高通股份有限公司 用于计算Turbo解码器的软判决输入度量的方法和装置
CN101232349A (zh) * 2007-12-27 2008-07-30 复旦大学 快速生成qam比特置信度软判决度量的方法
CN101217355A (zh) * 2008-01-11 2008-07-09 清华大学 基于低密度奇偶校验译码软判决信息的迭代定时同步方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于软判决反馈的猝发通信同步算法;游行远 等;《计算机工程与应用》;20180531;全文 *
快速生成QAM比特置信度软判决度量的算法;赵晋 等;《电子与信息学报》;20090430;全文 *

Also Published As

Publication number Publication date
CN109672501A (zh) 2019-04-23

Similar Documents

Publication Publication Date Title
EP3614637B1 (en) Systems and methods for adjusting the sample timing of a gfsk modulated signal
US8166379B1 (en) Calculating soft information from a multi-level modulation signal
CN103124250B (zh) 接收设备、信号处理设备以及信号处理方法
EP1453262A1 (en) Iterative MMSE detection
JP2008016959A (ja) 復号装置、方法およびプログラム
CN110336639B (zh) 一种基于容量分布的极化码多层编码调制方法及应用
CN105144650A (zh) 用于qam数据信号的增强式解码和解映射方法和装置
KR100943519B1 (ko) 로그우도비 산출 방법 및 송신 신호 검출 방법
CN109672501B (zh) 正交调制的软判决度量方法和装置
CN111585698A (zh) 一种多维编码调制方法及译码解调方法
KR100926563B1 (ko) 로그우도비 산출 방법, 송신 신호 검출 방법 및 수신기
CN107800511A (zh) 一种多级编码调制系统映射方法及装置
KR101072559B1 (ko) 대수우도비를 산출하는 방법 및 복조 장치
US20230344685A1 (en) A communication unit for soft-decision demodulation and method therefor
CN111865867B (zh) 一种256apsk调制信号的解调方法、系统及装置
CN108734188A (zh) 一种聚类方法、设备及存储介质
CN107070602A (zh) 一种基于k均值聚类算法的空间调制系统盲检测方法
CN109639618B (zh) 一种适用于高阶正交幅度调制信号的低复杂度软输出解调方法
CN114826853B (zh) 基于fpga的对数似然比方法及系统
KR100926566B1 (ko) 연판정값 산출 방법 및 송신 신호 검출 방법
JP2008092519A (ja) 受信装置
CN1175581C (zh) 相加-比较选择电路
CN111277527B (zh) N阶维特比思想应用在m-gfsk相位域的解调方法
WO2024091727A1 (en) Decoding correlated received multi-symbol data by guessing noise effect
CN115941419A (zh) 一种广义星座调制映射方法、装置和通信系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant