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

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

Info

Publication number
CN109302264B
CN109302264B CN201710608817.0A CN201710608817A CN109302264B CN 109302264 B CN109302264 B CN 109302264B CN 201710608817 A CN201710608817 A CN 201710608817A CN 109302264 B CN109302264 B CN 109302264B
Authority
CN
China
Prior art keywords
sub
region
constellation
max
subregion
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
CN201710608817.0A
Other languages
English (en)
Other versions
CN109302264A (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

Images

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个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
可选地，所述判决得到与所述发送符号对应的判决星座点，包括以下任意一种：采用迫零准则判决得到与所述发送符号对应的判决星座点；采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
可选地，所述预设大小的子区域为长方形。
可选地，采用如下方式选取对应的预设大小子区域：
当满足
Figure BDA0001358733890000021
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足
Figure BDA0001358733890000022
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ,subregionxmax＝Cmax；
当满足
Figure BDA0001358733890000023
且
Figure BDA0001358733890000024
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
Figure BDA0001358733890000025
当满足
Figure BDA0001358733890000026
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足
Figure BDA0001358733890000027
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足
Figure BDA0001358733890000028
且
Figure BDA0001358733890000029
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
Figure BDA0001358733890000031
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，
Figure BDA0001358733890000032
表示下取整运算，
Figure BDA0001358733890000033
表示上取整运算。
本发明实施例提供了一种最大似然检测的装置，所述装置包括：分解单元，适于根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；计算单元，适于对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
可选地，所述计算单元，适于或者采用迫零准则判决得到与所述发送符号对应的判决星座点；或者采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
可选地，所述预设大小的子区域为长方形。
可选地，所述计算单元，适于采用如下方式选取对应的预设大小子区域：
当满足
Figure BDA0001358733890000041
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足
Figure BDA0001358733890000042
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ,subregionxmax＝Cmax；
当满足
Figure BDA0001358733890000043
且
Figure BDA0001358733890000044
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
Figure BDA0001358733890000045
当满足
Figure BDA0001358733890000046
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足
Figure BDA0001358733890000047
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足
Figure BDA0001358733890000048
且
Figure BDA0001358733890000049
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
Figure BDA0001358733890000051
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，
Figure BDA0001358733890000052
表示下取整运算，
Figure BDA0001358733890000053
表示上取整运算。
本发明实施例提供了一种计算机可读存储介质，其上存储有计算机指令，所述计算机指令运行时执行以上任一种所述的最大似然检测的方法的步骤。
本发明实施例提供了一种最大似然检测的装置，包括存储器和处理器，所述存储器上存储有能够在所述处理器上运行的计算机指令，所述处理器运行所述计算机指令时执行以上任一种所述的最大似然检测的方法的步骤。
与现有技术相比，本发明实施例的技术方案具有以下优点：
上述的方案，通过判决得到与所述发送符号对应的判决星座点，进而围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。
进一步，采用迫零准则判决得到与所述发送符号对应的判决星座点，可以以较低的计算复杂度获得相对较高的最大似然检测准确度。
进一步，通过设置子区域的形状为长方形，可以使得子区域内涵盖的星座点数最少，故可以降低最大似然检测的复杂度。
附图说明
图1是本发明实施例中的一种最大似然检测的方法的流程图；
图2是本发明实施例中的另一种最大似然检测的方法的流程图；
图3是本发明实施例中的一种星座图的结构示意图；
图4是本发明实施例中的一种最大似然检测的装置的结构示意图。
具体实施方式
目前，在计算对数似然比时，先对信道矩阵和接收信号进行QR分解，检测发送符号，然后分别计算发送符号对应的所有星座点的欧式距离，进而根据欧式距离计算得到每个发送符号的对数似然比。比如，以两发送天线，两接收天线的信道模型为例，经过QR分解后的系统模型为y＝Rs+n。
其中：
Figure BDA0001358733890000061
是接收信号，
Figure BDA0001358733890000062
是信道矩阵，
Figure BDA0001358733890000063
是功率为1的白化噪声，
Figure BDA0001358733890000064
是发送符号。
则得到的发送符号中存在一个符号S1，则符号S1的第i个比特的对数似然比LLR1,i的计算公式(1)为：
Figure BDA0001358733890000065
其中：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满足
Figure BDA0001358733890000081
时，所述子区域在星座图上的最小横坐标subregionxmin可以满足如下公式(2)：
subregionxmin＝-Cmax (2)
同时，最大横坐标subregionxmax满足如下公式(3)：
subregionxmax＝-Cmax+(N-1)·Δ (3)
其中：
Figure BDA0001358733890000082
表示下取整运算，
Figure BDA0001358733890000083
表示上取整运算
当满足
Figure BDA0001358733890000084
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(4)及(5)：
subregionxmin＝Cmax-(N-1)·Δ (4)
subregionxmax＝Cmax (5)
当满足
Figure BDA0001358733890000085
且
Figure BDA0001358733890000086
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(6)及(7)：
Figure BDA0001358733890000091
Figure BDA0001358733890000092
当满足
Figure BDA0001358733890000093
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(8)及(9)：
subregionymin＝-Rmax (8)
subregionymax＝-Rmax+(M-1)·Δ (9)
当满足
Figure BDA0001358733890000094
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(10)及(11)：
subregionymin＝Rmax-(M-1)·Δ (10)
subregionymax＝Rmax (11)
当满足
Figure BDA0001358733890000095
且
Figure BDA0001358733890000096
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(12)及(13)：
Figure BDA0001358733890000097
Figure BDA0001358733890000098
在本发明一实施例中，可以设置Δ为2。本领域技术人员根据实际需要，也可以设置Δ为其他数值，但无论Δ大小为多少，均在本发明的保护范围内。
目前，需要遍历一个符号对应的所有星座点，以计算得到符号对应的最大似然比，该方案存在计算复杂度高，且硬件实现困难的问题。
而本发明实施例中通过判决得到与所述发送符号对应的判决星座点，进而围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。
为使得本领域技术人员更好地理解和实现本发明，下面以2个发送符号为例，对本发明实施例中的最大似然检测的方法进行具体说明。图2示出了本发明实施例中的另一种最大似然检测的方法的流程图，下面参考图2对所述方法进行详细介绍：
步骤S201：接收信号。
步骤S202：QR分解。
以两个发送天线，两接收天线的信道模型为例，经过分解后的系统模型可以如下公式(14)所示：
Figure BDA0001358733890000101
其中
Figure BDA0001358733890000102
是接收信号，
Figure BDA0001358733890000103
是信道矩阵，
Figure BDA0001358733890000104
是功率为1的白化噪声，
Figure BDA0001358733890000105
是发送符号。
因此有公式(15)：
y1＝r11s1+n1 (15)
步骤S203：判决发送符号S1的星座点。
在具体实施中，可以根迫零准则按照如下公式(16)计算发送符号S1的迫零星座点：
Figure BDA0001358733890000106
其中：Q()表示将一个值映射到最近的星座点。
步骤S204：判决发送符号S0的星座点。
需要说明的是，步骤S204与步骤S203的实施方法相同，在此不再赘述。并且步骤S203及步骤S204的执行顺序并没有限定，本领域技术人员根据实际需要，可以先执行步骤S203，可以先执行步骤S204，也可以同时执行步骤S203及步骤S204。
步骤S205：选取与发送符号S1对应的子区域。
为便于说明和理解，图3示出了本发明实施例中的一种星座图，在图3中，调制方法为以256-QAM，子区域大小为81为例，子区域为
Figure BDA0001358733890000111
迫零星座点为
Figure BDA0001358733890000112
基本思想是选取围绕迫零星座点的长方形区域，256-QAM的星座点关于x轴和y轴都是对称的，星座点在x轴的位置为-15：15，间隔为2，星座点在y轴的位置为-15：15。相邻星座点的横坐标及纵坐标的最小间隔均为2。假设迫零星座点
Figure BDA0001358733890000113
在x轴的位置是a，在y轴的位置是b。用(subregionxmin,subregionxmax)表示子区域
Figure BDA0001358733890000114
在x轴的区间，用(subregionymin,subregionymax)表示子区域
Figure BDA0001358733890000115
在y轴的区间。由于子区域
Figure BDA0001358733890000116
是长方形，确定了这两个区间，就确定了子区域
Figure BDA0001358733890000117
的范围。
详细的说，子区域
Figure BDA0001358733890000118
选取过程如下：
如果a-8≥-15并且a+8≤15，则所述子区域
Figure BDA00013587338900001110
在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(17)及(18)：
subregionxmin＝a-8 (17)
subregionxmax＝a+8 (18)
如果a-8＜-15，则所述子区域
Figure BDA00013587338900001111
在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(19)及(20)：
subregionxmin＝-15 (19)
subregionxmax＝1 (20)
如果a+8＞15，则所述子区域
Figure BDA0001358733890000119
在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式(21)及(22)：
subregionxmin＝-1 (21)
subregionxmax＝15 (22)
如果b-8≥-15并且b+8≤15，则所述子区域
Figure BDA0001358733890000121
在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(23)及(24)：
subregionymin＝b-8 (23)
subregionymax＝b+8 (24)
如果b-8＜-15，则所述子区域
Figure BDA0001358733890000122
在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(25)及(26)：
subregionymin＝-15 (25)
subregionymax＝1 (26)
如果b+8＞15，则所述子区域
Figure BDA0001358733890000123
在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式(27)及(28)：
subregionymin＝-1 (27)
subregionymax＝15 (28)
因此，按照如上方法选取子区域
Figure BDA0001358733890000124
可以降低最大似然检测过程的复杂度，且兼顾提高检测性能。当发送符号S1的调制方式为别的高阶调制方式或者子区域大小不为81个点时，S1的子区域
Figure BDA0001358733890000125
选择方法和上述步骤类似，在此不再赘述。
步骤S206：选取与发送符号S0对应的子区域。
需要说明的是，步骤S206与步骤S205的实施方法相同，在此不再赘述。并且步骤S205及步骤S206的执行顺序并没有限定，本领域技术人员根据实际需要，可以先执行步骤S205，可以先执行步骤S206，也可以同时执行步骤S206及步骤S205。
步骤S207：对发送符号S0对应的子区域内的星座点计算欧式距离。
步骤S208：对发送符号S1对应的子区域内的星座点计算欧式距离。
需要说明的是，步骤S207及步骤S208的执行顺序并没有限定，本领域技术人员根据实际需要，可以先执行步骤S207，可以先执行步骤S208，也可以同时执行步骤S207及步骤S208。
步骤S209：根据子区域内的星座点计算得到的欧式距离，计算发送符号S1对应的对数似然比。
在具体实施中，计算得到迫零星座点
Figure BDA0001358733890000131
和选取子区域
Figure BDA0001358733890000132
之后，可以按照如下公式(29)计算发送符号S1的第i个比特的对数似然比为：
Figure BDA0001358733890000133
步骤S210：根据子区域内的星座点计算得到的欧式距离，计算发送符号S0对应的对数似然比。
需要说明的是，步骤S210与步骤S209的实施方法相同，在此不再赘述。并且，步骤S209及步骤S210的执行顺序并没有限定，本领域技术人员根据实际需要，可以先执行步骤S209，可以先执行步骤S210，也可以同时执行步骤S209及步骤S210。
综上，本发明实施例通过判决得到与所述发送符号对应的判决星座点，进而围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。
为使得本领域技术人员更好地理解和实现本发明，图4示出了本发明实施例中的一种最大似然检测的装置的结构示意图，如图4所示，所述装置可以包括：分解单元41及计算单元42，其中：
所述分解单元41，适于根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；
所述计算单元42，适于对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比。
综上，本发明实施例中通过最大似然检测装置的计算单元判决得到与所述发送符号对应的判决星座点，进而计算单元围绕所述判决星座点选取预设大小的子区域，然后根据所述预设大小的子区域内的星座点即可计算欧式距离，而无需遍历发送符号对应的所有的星座点，因此可以降低对数似然比计算的复杂度，进而降低最大似然检测的复杂度，并且便于硬件实现。。
在具体实施中，所述计算单元42，适于或者采用迫零准则判决得到与所述发送符号对应的判决星座点；或者采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
在具体实施中，所述预设大小的子区域为长方形。
在具体实施中，所述计算单元42，适于采用如下方式选取对应的预设大小子区域：
当满足
Figure BDA0001358733890000141
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足
Figure BDA0001358733890000142
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ,subregionxmax＝Cmax；
当满足
Figure BDA0001358733890000143
且
Figure BDA0001358733890000144
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
Figure BDA0001358733890000151
当满足
Figure BDA0001358733890000152
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足
Figure BDA0001358733890000153
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足
Figure BDA0001358733890000154
且
Figure BDA0001358733890000155
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
Figure BDA0001358733890000156
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin,subregionymax)表示所述子区域星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，
Figure BDA0001358733890000157
表示下取整运算，
Figure BDA0001358733890000158
表示上取整运算。
为便于理解，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机指令，所述计算机指令运行时执行上述本发明实施例中提供的任一种最大似然检测的方法的步骤。
本发明实施例中还提供了一种最大似然检测的装置，包括存储器和处理器，所述存储器上存储有能够在所述处理器上运行的计算机指令，所述处理器运行所述计算机指令时执行上述本发明实施例提供的任一种最大似然检测的方法的步骤。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于以计算机可读存储介质中，存储介质可以包括：ROM、RAM、磁盘或光盘等。
虽然本发明披露如上，但本发明并非限定于此。任何本领域技术人员，在不脱离本发明的精神和范围内，均可作各种更动与修改，因此本发明的保护范围应当以权利要求所限定的范围为准。

Claims (8)

1.一种最大似然检测的方法，其特征在于，包括：
根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；
对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；
围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；
基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比；
采用如下方式选取对应的预设大小子区域：
当满足
Figure FDA0002927116990000011
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足
Figure FDA0002927116990000012
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝Cmax-(N-1)·Δ，subregionxmax＝Cmax；
当满足
Figure FDA0002927116990000013
且
Figure FDA0002927116990000014
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
Figure FDA0002927116990000015
当满足
Figure FDA0002927116990000021
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregiony max分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足
Figure FDA0002927116990000022
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregiony max分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足
Figure FDA0002927116990000023
且
Figure FDA0002927116990000024
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
Figure FDA0002927116990000025
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subregionxmin，subregionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin，subregionymax)表示所述子区域在星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，
Figure FDA0002927116990000026
表示下取整运算，
Figure FDA0002927116990000027
表示上取整运算。
2.如权利要求1所述的最大似然检测的方法，其特征在于，所述判决得到与所述发送符号对应的判决星座点，包括以下任意一种：
采用迫零准则判决得到与所述发送符号对应的判决星座点；
采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
3.如权利要求1或2任一项所述的最大似然检测的方法，其特征在于，所述预设大小的子区域为长方形。
4.一种最大似然检测的装置，其特征在于，包括：
分解单元，适于根据噪声白化后的接收信号及信道矩阵，进行QR分解检测p个发送符号，p为不小于1的正整数；
计算单元，适于对p个发送符号中的任意一个发送符号，执行如下操作得到对数似然比：判决得到与所述发送符号对应的判决星座点；围绕所述判决星座点选取预设大小的子区域；其中：所述子区域内包括的星座点数目小于所述发送符号对应的所有星座点数目；基于所述预设大小的子区域内的星座点计算欧式距离，且根据所述欧式距离计算所述发送符号的对数似然比；
所述计算单元，适于采用如下方式选取对应的预设大小子区域：
当满足
Figure FDA0002927116990000031
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin＝-Cmax，subregionxmax＝-Cmax+(N-1)·Δ；
当满足
Figure FDA0002927116990000032
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
subregionxmin-Cmax-(N-1)·Δ，subregionxmax-Cmax；
当满足
Figure FDA0002927116990000033
且
Figure FDA0002927116990000034
时，所述子区域在星座图上的最小横坐标subregionxmin及最大横坐标subregionxmax分别满足如下公式：
Figure FDA0002927116990000041
当满足
Figure FDA0002927116990000042
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝-Rmax，subregionymax＝-Rmax+(M-1)·Δ；
当满足
Figure FDA0002927116990000043
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
subregionymin＝Rmax-(M-1)·Δ，subregionymax＝Rmax；
当满足
Figure FDA0002927116990000044
且
Figure FDA0002927116990000045
时，所述子区域在星座图上的最小纵坐标subregionymin及最大纵坐标subregionymax分别满足如下公式：
Figure FDA0002927116990000046
其中：Δ表示相邻星座点的横坐标以及纵坐标的最小间隔，a表示所述判决星座点在星座图上的横坐标，b表示所述判决星座点在星座图上的纵坐标，星座图的横坐标范围为(-Cmax，Cmax)，星座图的纵坐标范围为(-Rmax，Rmax)，(subrcgionxmin，subrcgionxmax)表示所述子区域在星座图上的横坐标区间，(subregionymin，subregionymax)表示所述子区域在星座图上的纵坐标区间，N为所述子区域在星座图的横坐标上投影范围所包括的星座点数，M为所述子区域在星座图的纵坐标上投影范围所包括的星座点数，
Figure FDA0002927116990000047
表示下取整运算，
Figure FDA0002927116990000048
表示上取整运算。
5.如权利要求4所述的最大似然检测的装置，其特征在于，所述计算单元，适于或者采用迫零准则判决得到与所述发送符号对应的判决星座点；或者采用最小均方误差准则判决得到与所述发送符号对应的判决星座点。
6.如权利要求4或5任一项所述的最大似然检测的装置，其特征在于，所述预设大小的子区域为长方形。
7.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，所述计算机指令被处理器运行时执行权利要求1至3任一项所述的最大似然检测的方法的步骤。
8.一种最大似然检测的装置，其特征在于，包括存储器和处理器，所述存储器上存储有能够在所述处理器上运行的计算机指令，所述处理器运行所述计算机指令时执行权利要求1至3任一项所述的最大似然检测的方法的步骤。
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 CN109302264A (zh) 2019-02-01
CN109302264B true 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
CN109302264A (zh) 2019-02-01

Similar Documents

Publication Publication Date Title
EP4047511A1 (en) Method and device for detecting vehicle
CN109493367B (zh) 一种对目标对象进行跟踪的方法和设备
CN110139359B (zh) 一种干扰源定位处理方法及装置
CN112560684B (zh) 车道线检测方法、装置、电子设备、存储介质以及车辆
CN109145251B (zh) 一种改进型同步扰动随机逼近算法的大气参数求解方法
CN109785247B (zh) 激光雷达异常点云数据的修正方法、装置及存储介质
EP3872762A1 (en) Edge-based three-dimensional tracking and registration method and apparatus for augmented reality, and storage medium
CN108627875A (zh) 雷达偏移速度确定和偏移处理方法、存储介质、电子设备
CN111738396A (zh) 一种应用于潜艇路径规划的自适应栅格颗粒度蚁群算法
CN109284062B (zh) 一种触摸数据处理方法、装置、终端和介质
CN116032823B (zh) 基于数字高程地图的散射通信链路快速计算方法及系统
CN114567396A (zh) 无线通信方法、非线性函数的拟合方法、终端及设备
US9569866B2 (en) Flexible video object boundary tracking
CN109302264B (zh) 最大似然检测的方法、装置及计算机可读介质
CN113781569A (zh) 一种回环检测方法和装置
CN112087235A (zh) 基于伪逆感知字典的稀疏度自适应doa估计方法及系统
JP5078669B2 (ja) 目標検出装置、目標検出方法及び目標検出プログラム
JP5266976B2 (ja) 端末位置決定方法、端末位置決定プログラム、端末位置決定装置
CN105631856B (zh) 自适应调整的红外船舶图像活动轮廓分割方法
CN111224908B (zh) 信号序列的检测方法及装置、存储介质、终端
CN110705334A (zh) 目标追踪方法、装置、设备和介质
US11689693B2 (en) Video frame interpolation method and device, computer readable storage medium
CN110602635B (zh) 一种室内地图匹配增强定位方法、设备及存储设备
CN111757242B (zh) 一种基于ais信息计算的海域通信波束控制方法及装置
CN114417946A (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