CN109302264A - 最大似然检测的方法、装置及计算机可读介质 - Google Patents

最大似然检测的方法、装置及计算机可读介质 Download PDF

Info

Publication number
CN109302264A
CN109302264A CN201710608817.0A CN201710608817A CN109302264A CN 109302264 A CN109302264 A CN 109302264A CN 201710608817 A CN201710608817 A CN 201710608817A CN 109302264 A CN109302264 A CN 109302264A
Authority
CN
China
Prior art keywords
subregion
planisphere
max
constellation point
abscissa
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.)
Granted
Application number
CN201710608817.0A
Other languages
English (en)
Other versions
CN109302264B (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.)
Spreadtrum Communications Shanghai Co Ltd
Original Assignee
Spreadtrum Communications Shanghai Co Ltd
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 Spreadtrum Communications Shanghai Co Ltd filed Critical Spreadtrum Communications Shanghai Co Ltd
Priority to CN201710608817.0A priority Critical patent/CN109302264B/zh
Publication of CN109302264A publication Critical patent/CN109302264A/zh
Application granted granted Critical
Publication of CN109302264B publication Critical patent/CN109302264B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045—Arrangements at the receiver end
    • H04L1/0054—Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/38—Demodulator circuits; Receiver circuits

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

一种最大似然检测的方法、装置及计算机可读介质，所述方法包括：根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。采用上述方案可以降低最大似然检测的复杂度，便于硬件实现。

Description

最大似然检测的方法、装置及计算机可读介质
技术领域
本发明实施例涉及通信领域，尤其涉及一种最大似然检测的方法、装置及计算机可读介质。
背景技术
在对接收信号进行处理的过程中，会涉及到最大似然检测，需要计算每个发送符号的对数似然比(Log Likelihood Ratio，LLR)。
目前，在计算对数似然比时，先对信道矩阵和接收信号进行QR分解，检测发送符号，然后分别计算发送符号对应的所有星座点的欧式距离，进而根据欧式距离计算得到每个发送符号的对数似然比。
但是，上述的最大似然检测的方法，计算复杂度高，且硬件实现困难。
发明内容
本发明实施例解决的问题是如何降低最大似然检测的复杂度，及便于硬件实现。
为解决上述问题，本发明实施例提供了一种最大似然检测的方法，所述方法包括：根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
可选地，所述判决得到与所述发送符号对应的判决星座点，包括以下任意一种：采用迫零准则判决得到与所述发送符号对应的判决星座点；采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
可选地，所述预设大小的子区域为长方形。
可选地，采用如下方式选取对应的预设大小子区域：
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ,subregionxmax＝Cmax；
当满足且时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足且时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，表示下取整运算，表示上取整运算。
本发明实施例提供了一种最大似然检测的装置，所述装置包括：分解单元，适于根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；计算单元，适于对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
可选地，所述计算单元，适于或者采用迫零准则判决得到与所述发送符号对应的判决星座点；或者采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
可选地，所述预设大小的子区域为长方形。
可选地，所述计算单元，适于采用如下方式选取对应的预设大小子区域：
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ,subregionxmax＝Cmax；
当满足且时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足且时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，表示下取整运算，表示上取整运算。
本发明实施例提供了一种计算机可读存储介质，其上存储有计算机指令，所述计算机指令运行时执行以上任一种所述的最大似然检测的方法的步骤。
本发明实施例提供了一种最大似然检测的装置，包括存储器和处理器，所述存储器上存储有能够在所述处理器上运行的计算机指令，所述处理器运行所述计算机指令时执行以上任一种所述的最大似然检测的方法的步骤。
与现有技术相比，本发明实施例的技术方案具有以下优点：
上述的方案，通过判决得到与所述发送符号对应的判决星座点，进而围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。
进一步，采用迫零准则判决得到与所述发送符号对应的判决星座点，可以以较低的计算复杂度获得相对较高的最大似然检测准确度。
进一步，通过设置子区域的形状为长方形，可以使得子区域内涵盖的星座点数最少，故可以降低最大似然检测的复杂度。
附图说明
图1是本发明实施例中的一种最大似然检测的方法的流程图；
图2是本发明实施例中的另一种最大似然检测的方法的流程图；
图3是本发明实施例中的一种星座图的结构示意图；
图4是本发明实施例中的一种最大似然检测的装置的结构示意图。
具体实施方式
目前，在计算对数似然比时，先对信道矩阵和接收信号进行QR分解，检测发送符号，然后分别计算发送符号对应的所有星座点的欧式距离，进而根据欧式距离计算得到每个发送符号的对数似然比。比如，以两发送天线，两接收天线的信道模型为例，经过QR分解后的系统模型为y＝Rs+n。
其中：是接收信号，是信道矩阵，是功率为1的白化噪声，是发送符号。
则得到的发送符号中存在一个符号S1，则符号S1的第i个比特的对数似然比LLR1,i的计算公式(1)为：
其中：S1i表示发送符号S1的第i个比特，Rs表示符号S1所对应的信道矩阵，y表示接收信号。从公式(1)可以看出，计算符号S1的对数似然比，需要遍历符号S1的所有星座点。假设符S1的调制阶数为L，则需要计算L个欧式距离。而经过QR分解检测的其它符号的对数似然比的计算过程也类似于符号S1，此处不再赘述。由此可见，现在的最大似然检测的方法存在计算复杂度高且硬件实现困难的问题。
为解决上述问题，本发明实施例通过判决得到与所述发送符号对应的判决星座点，进而围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。
为使本发明实施例的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明实施例的具体实施例做详细的说明。
图1示出了本发明实施例中的一种最大似然检测的方法流程图，下面参考图1，对所述方法进行分步骤详细介绍，所述方法可以按照如下步骤实施：
步骤S11：根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号。
在具体实施中，由于发送符号经过处理以及信道的传输，最终被接收端接收得到接收信号，故根据噪声白化后的接收信号及信道矩阵，进行QR分解可以检测p个发送符号，并且p可以为不小于1的正整数。
步骤S12：对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
在具体实施中，所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目。
在具体实施中，可以有多种方法判决得到与所述发送符号对应的判决星座点。例如，为了提高最大似然检测的准确度，可以采用迫零准则判决得到与所述发送符号对应的判决星座点。又如，也可以采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
在具体实施中，所述围绕所述判决星座点选取预设大小的子区域，换言之，可以为以所述判决星座点为中心，进而选取预设大小的子区域。
在具体实施中，所述预设大小的子区域为长方形。因此可以使得子区域内所涵盖的星座点数最少，进而可以降低最大似然检测的复杂度。
在具体实施中，可以根据星座图的坐标信息、大小信息及判决星座点的坐标，来相应地选取对应的预设大小子区域。
详细地说，一般相邻星座点的横坐标的最小间隔与纵坐标标的最小间隔相同，故可以以Δ表示相邻星座点的横坐标以及纵坐标的最小间隔。a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M表示所述子区域在星座图的纵坐标上投影范围所包括的星座点数。
则当a满足时，所述子区域在星座图上的最小横坐标subregionxmin可以满足如下公式(2)：
subregionxmin＝-Cmax (2)
同时，最大横坐标subregionxmax满足如下公式(3)：
subregionxmax＝-Cmax+(N-1)·Δ (3)
其中：表示下取整运算，表示上取整运算
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(4)及(5)：
subregionxmin＝Cmax-(N-1)·Δ (4)
subregionxmax＝Cmax (5)
当满足且时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(6)及(7)：
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(8)及(9)：
subregionymin＝-Rmax (8)
subregionymax＝-Rmax+(M-1)·Δ (9)
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(10)及(11)：
subregionymin＝Rmax-(M-1)·Δ (10)
subregionymax＝Rmax (11)
当满足且时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(12)及(13)：
在本发明一实施例中，可以设置Δ为2。本领域技术人员根据实际需要，也可以设置Δ为其他数值，但无论Δ大小为多少，均在本发明的保护范围内。
目前，需要遍历一个符号对应的所有星座点，以计算得到符号对应的最大似然比，该方案存在计算复杂度高，且硬件实现困难的问题。
而本发明实施例中通过判决得到与所述发送符号对应的判决星座点，进而围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。
为使得本领域技术人员更好地理解和实现本发明，下面以2个发送符号为例，对本发明实施例中的最大似然检测的方法进行具体说明。图2示出了本发明实施例中的另一种最大似然检测的方法的流程图，下面参考图2对所述方法进行详细介绍：
步骤S201：接收信号。
步骤S202：QR分解。
以两个发送天线，两接收天线的信道模型为例，经过分解后的系统模型可以如下公式(14)所示：
其中是接收信号，是信道矩阵，是功率为1的白化噪声，是发送符号。
因此有公式(15)：
y1＝r11s1+n1 (15)
步骤S203：判决发送符号S1的星座点。
在具体实施中，可以根迫零准则按照如下公式(16)计算发送符号S1的迫零星座点：
其中：Q()表示将一个值映射到最近的星座点。
步骤S204：判决发送符号S0的星座点。
需要说明的是，步骤S204与步骤S203的实施方法相同，在此不再赘述。并且步骤S203及步骤S204的执行顺序并没有限定，本领域技术人员根据实际需要，可以先执行步骤S203，可以先执行步骤S204，也可以同时执行步骤S203及步骤S204。
步骤S205：选取与发送符号S1对应的子区域。
为便于说明和理解，图3示出了本发明实施例中的一种星座图，在图3中，调制方法为以256-QAM，子区域大小为81为例，子区域为迫零星座点为基本思想是选取围绕迫零星座点的长方形区域，256-QAM的星座点关于x轴和y轴都是对称的，星座点在x轴的位置为-15：15，间隔为2，星座点在y轴的位置为-15：15。相邻星座点的横坐标及纵坐标的最小间隔均为2。假设迫零星座点在x轴的位置是a，在y轴的位置是b。用(subregionxmin,subregionxmax)表示子区域在x轴的区间，用(subregionymin,subregionymax)表示子区域在y轴的区间。由于子区域是长方形，确定了这两个区间，就确定了子区域的范围。
详细的说，子区域选取过程如下：
如果a-8≥-15并且a+8≤15，则所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(17)及(18)：
subregionxmin＝a-8 (17)
subregionxmax＝a+8 (18)
如果a-8＜-15，则所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(19)及(20)：
subregionxmin＝-15 (19)
subregionxmax＝1 (20)
如果a+8＞15，则所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(21)及(22)：
subregionxmin＝-1 (21)
subregionxmax＝15 (22)
如果b-8≥-15并且b+8≤15，则所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(23)及(24)：
subregionymin＝b-8 (23)
subregionymax＝b+8 (24)
如果b-8＜-15，则所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(25)及(26)：
subregionymin＝-15 (25)
subregionymax＝1 (26)
如果b+8＞15，则所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(27)及(28)：
subregionymin＝-1 (27)
subregionymax＝15 (28)
因此，按照如上方法选取子区域可以降低最大似然检测过程的复杂度，且兼顾提高检测性能。当发送符号S1的调制方式为别的高阶调制方式或者子区域大小不为81个点时，S1的子区域选择方法和上述步骤类似，在此不再赘述。
步骤S206：选取与发送符号S0对应的子区域。
需要说明的是，步骤S206与步骤S205的实施方法相同，在此不再赘述。并且步骤S205及步骤S206的执行顺序并没有限定，本领域技术人员根据实际需要，可以先执行步骤S205，可以先执行步骤S206，也可以同时执行步骤S206及步骤S205。
步骤S207：对发送符号S0对应的子区域内的星座点计算欧式距离。
步骤S208：对发送符号S1对应的子区域内的星座点计算欧式距离。
需要说明的是，步骤S207及步骤S208的执行顺序并没有限定，本领域技术人员根据实际需要，可以先执行步骤S207，可以先执行步骤S208，也可以同时执行步骤S207及步骤S208。
步骤S209：根据子区域内的星座点计算得到的欧式距离，计算发送符号S1对应的对数似然比。
在具体实施中，计算得到迫零星座点和选取子区域之后，可以按照如下公式(29)计算发送符号S1的第i个比特的对数似然比为：
步骤S210：根据子区域内的星座点计算得到的欧式距离，计算发送符号S0对应的对数似然比。
需要说明的是，步骤S210与步骤S209的实施方法相同，在此不再赘述。并且，步骤S209及步骤S210的执行顺序并没有限定，本领域技术人员根据实际需要，可以先执行步骤S209，可以先执行步骤S210，也可以同时执行步骤S209及步骤S210。
综上，本发明实施例通过判决得到与所述发送符号对应的判决星座点，进而围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。
为使得本领域技术人员更好地理解和实现本发明，图4示出了本发明实施例中的一种最大似然检测的装置的结构示意图，如图4所示，所述装置可以包括：分解单元41及计算单元42，其中：
所述分解单元41，适于根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；
所述计算单元42，适于对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
综上，本发明实施例中通过最大似然检测装置的计算单元判决得到与所述发送符号对应的判决星座点，进而计算单元围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。。
在具体实施中，所述计算单元42，适于或者采用迫零准则判决得到与所述发送符号对应的判决星座点；或者采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
在具体实施中，所述预设大小的子区域为长方形。
在具体实施中，所述计算单元42，适于采用如下方式选取对应的预设大小子区域：
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ,subregionxmax＝Cmax；
当满足且时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足且时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，表示下取整运算，表示上取整运算。
为便于理解，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机指令，所述计算机指令运行时执行上述本发明实施例中提供的任一种最大似然检测的方法的步骤。
本发明实施例中还提供了一种最大似然检测的装置，包括存储器和处理器，所述存储器上存储有能够在所述处理器上运行的计算机指令，所述处理器运行所述计算机指令时执行上述本发明实施例提供的任一种最大似然检测的方法的步骤。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于以计算机可读存储介质中，存储介质可以包括：ROM、RAM、磁盘或光盘等。
虽然本发明披露如上，但本发明并非限定于此。任何本领域技术人员，在不脱离本发明的精神和范围内，均可作各种更动与修改，因此本发明的保护范围应当以权利要求所限定的范围为准。

Claims (10)

1.一种最大似然检测的方法，其特征在于，包括：
根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，
p为不小于1的正整数；
对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：
判决得到与所述发送符号对应的判决星座点；
围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；
基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
2.如权利要求1所述的最大似然检测的方法，其特征在于，所述判决得到与所述发送符号对应的判决星座点，包括以下任意一种：
采用迫零准则判决得到与所述发送符号对应的判决星座点；
采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
3.如权利要求1或2任一项所述的最大似然检测的方法，其特征在于，所述预设大小的子区域为长方形。
4.如权利要求3所述的最大似然检测的方法，其特征在于，采用如下方式选取对应的预设大小子区域：
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ,subregionxmax＝Cmax；
当满足且时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足且时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，表示下取整运算，表示上取整运算。
5.一种最大似然检测的装置，其特征在于，包括：
分解单元，适于根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；
计算单元，适于对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
6.如权利要求5所述的最大似然检测的装置，其特征在于，所述计算单元，适于或者采用迫零准则判决得到与所述发送符号对应的判决星座点；或者采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
7.如权利要求5或6任一项所述的最大似然检测的装置，其特征在于，所述预设大小的子区域为长方形。
8.如权利要求7所述的最大似然检测的装置，其特征在于，所述计算单元，适于采用如下方式选取对应的预设大小子区域：
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ,subregionxmax＝Cmax；
当满足且时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足且时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，表示下取整运算，表示上取整运算。
9.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，所述计算机指令运行时执行权利要求1至4任一项所述的最大似然检测的方法的步骤。
10.一种最大似然检测的装置，其特征在于，包括存储器和处理器，所述存储器上存储有能够在所述处理器上运行的计算机指令，所述处理器运行所述计算机指令时执行权利要求1至4任一项所述的最大似然检测的方法的步骤。
CN201710608817.0A 2017-07-24 2017-07-24 最大似然检测的方法、装置及计算机可读介质 Active CN109302264B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710608817.0A CN109302264B (zh) 2017-07-24 2017-07-24 最大似然检测的方法、装置及计算机可读介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710608817.0A CN109302264B (zh) 2017-07-24 2017-07-24 最大似然检测的方法、装置及计算机可读介质

Publications (2)

Publication Number Publication Date
CN109302264A true CN109302264A (zh) 2019-02-01
CN109302264B CN109302264B (zh) 2021-04-02

Family

ID=65167243

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710608817.0A Active CN109302264B (zh) 2017-07-24 2017-07-24 最大似然检测的方法、装置及计算机可读介质

Country Status (1)

Country Link
CN (1) CN109302264B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1697433A (zh) * 2004-04-30 2005-11-16 索尼国际(欧洲)股份有限公司 采用最大似然检测的均衡结构和方法
CN1863184A (zh) * 2004-05-13 2006-11-15 株式会社Ntt都科摩 信号分解装置和信号分解方法
CN1968067A (zh) * 2005-11-16 2007-05-23 松下电器产业株式会社 简化最大似然检测方法及检测设备
US7453965B1 (en) * 1999-02-10 2008-11-18 At&T Corp. Differential transmitter diversity technique for wireless communications
CN101394392A (zh) * 2008-11-12 2009-03-25 北京邮电大学 一种ofdm系统的信号分集的方法
CN101741514A (zh) * 2009-11-24 2010-06-16 北京天碁科技有限公司 在信号检测过程中确定需要检测的星座点的方法和装置
CN103166903A (zh) * 2013-03-19 2013-06-19 清华大学 星座映射的软解预处理方法以及软解方法
CN103378898A (zh) * 2012-04-24 2013-10-30 马维尔国际有限公司 Mimo系统中的信号检测方法和设备
CN104283835A (zh) * 2013-07-02 2015-01-14 大唐联诚信息系统技术有限公司 一种正交幅度调制软比特解调方法及其装置
CN104601298A (zh) * 2015-02-17 2015-05-06 美商威睿电通公司 在多重输入输出通信系统中检测信号的方法及装置

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7453965B1 (en) * 1999-02-10 2008-11-18 At&T Corp. Differential transmitter diversity technique for wireless communications
CN1697433A (zh) * 2004-04-30 2005-11-16 索尼国际(欧洲)股份有限公司 采用最大似然检测的均衡结构和方法
CN1863184A (zh) * 2004-05-13 2006-11-15 株式会社Ntt都科摩 信号分解装置和信号分解方法
CN1968067A (zh) * 2005-11-16 2007-05-23 松下电器产业株式会社 简化最大似然检测方法及检测设备
CN101394392A (zh) * 2008-11-12 2009-03-25 北京邮电大学 一种ofdm系统的信号分集的方法
CN101741514A (zh) * 2009-11-24 2010-06-16 北京天碁科技有限公司 在信号检测过程中确定需要检测的星座点的方法和装置
CN103378898A (zh) * 2012-04-24 2013-10-30 马维尔国际有限公司 Mimo系统中的信号检测方法和设备
CN103166903A (zh) * 2013-03-19 2013-06-19 清华大学 星座映射的软解预处理方法以及软解方法
CN104283835A (zh) * 2013-07-02 2015-01-14 大唐联诚信息系统技术有限公司 一种正交幅度调制软比特解调方法及其装置
CN104601298A (zh) * 2015-02-17 2015-05-06 美商威睿电通公司 在多重输入输出通信系统中检测信号的方法及装置

Also Published As

Publication number Publication date
CN109302264B (zh) 2021-04-02

Similar Documents

Publication Publication Date Title
EP3171564B1 (en) Digital modulation method and device
CN108551371B (zh) 一种qam信号盲载波相位估计方法及系统
CN103166903A (zh) 星座映射的软解预处理方法以及软解方法
CN114677565A (zh) 特征提取网络的训练方法和图像处理方法、装置
Abdel‐Moneim et al. An efficient modulation classification method using signal constellation diagrams with convolutional neural networks, Gabor filtering, and thresholding
KR102664025B1 (ko) 초기 후보 감축에 기반한 심볼 검출을 위한 중첩 룩업 테이블
CN109302264A (zh) 最大似然检测的方法、装置及计算机可读介质
CN107800491B (zh) 一种基于信道状态信息的无线感知抗干扰方法及装置
CN109461171A (zh) 基于多通道改进DoG滤波的红外弱小目标检测算法
CN114389757A (zh) 球形译码检测方法和装置、电子设备、存储介质
CN107329473A (zh) 基于均值和声搜索的机器人巡检路径规划方法
KR101586000B1 (ko) 대수우도비 산출 장치 및 방법
Cheon et al. Fast descriptor extraction method for a SURF‐based interest point
CN106683107A (zh) 基于roewa改进的sar图像边缘检测方法
CN110571825A (zh) 一种静止同步补偿器模型参数辨识方法及系统
CN116016225A (zh) 信息处理方法、信号定位方法、装置、设备和介质
CN102769600B (zh) 一种ofdm系统接收机中缩放因子解调制的方法及装置
CN110954879B (zh) 一种移动门限的数字检波方法及系统
CN104765776A (zh) 一种数据样本的聚类方法和装置
CN108896967A (zh) 基于杂波协方差矩阵估计的距离扩展目标检测方法及装置
WO2014031412A1 (en) Vectorization approach to isolating local maxima in an n-dimensional dataset
Hong Medical image segmentation based on accelerated Dijkstra algorithm
Sembiring et al. Performance analysis of discrete hartley transform based ofdm modulator and demodulator
Feng et al. An accurate modulation recognition method of QPSK signal
Shi et al. STTMC: A Few-shot Spatial Temporal Transductive Modulation Classifier

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