[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN113098571B - Massive MIMO system signal detection method, system, base station and storage medium - Google Patents

Massive MIMO system signal detection method, system, base station and storage medium Download PDF

Info

Publication number
CN113098571B
CN113098571B CN202110330579.8A CN202110330579A CN113098571B CN 113098571 B CN113098571 B CN 113098571B CN 202110330579 A CN202110330579 A CN 202110330579A CN 113098571 B CN113098571 B CN 113098571B
Authority
CN
China
Prior art keywords
vector
signal
initial
symbol vector
output
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
CN202110330579.8A
Other languages
Chinese (zh)
Other versions
CN113098571A (en
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.)
Anhui University
Original Assignee
Anhui University
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 Anhui University filed Critical Anhui University
Priority to CN202110330579.8A priority Critical patent/CN113098571B/en
Publication of CN113098571A publication Critical patent/CN113098571A/en
Application granted granted Critical
Publication of CN113098571B publication Critical patent/CN113098571B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明提供一种大规模MIMO系统信号检测方法、系统、基站及存储介质,所述检测方法包括对第一初始信号进行主动禁忌搜索检测,以获取初始估计向量;根据所述初始估计向量消除所述第一初始信号中的干扰信号,以获取第二初始信号;对所述第二初始信号进行消息传递检测以获取输出向量估计;根据所述输出向量估计对符号向量进行重构,以获取符号向量重构值;将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作,迭代结束后的最终的所述符号向量重构值作为检测结果输出。利用本发明,可以改善RTS算法在高阶调制下性能不佳的情况,且拥有更低的复杂度。

Figure 202110330579

The present invention provides a signal detection method, system, base station and storage medium for a massive MIMO system. The detection method includes performing active tabu search detection on a first initial signal to obtain an initial estimation vector; the interference signal in the first initial signal to obtain the second initial signal; the message passing detection is performed on the second initial signal to obtain the output vector estimate; the symbol vector is reconstructed according to the output vector estimate to obtain the symbol vector reconstruction value; the symbol vector reconstruction value is used as the input of the active tabu search to perform an iterative operation, and the final symbol vector reconstruction value after the iteration is output as the detection result. By using the present invention, the poor performance of the RTS algorithm under high-order modulation can be improved, and it has lower complexity.

Figure 202110330579

Description

大规模MIMO系统信号检测方法、系统、基站及存储介质Massive MIMO system signal detection method, system, base station and storage medium

技术领域technical field

本发明涉及多输入多输出信号检测技术领域,特别涉及大规模MIMO系统信号检测方法、系统、基站及存储介质。The present invention relates to the technical field of multiple-input multiple-output signal detection, in particular to a signal detection method, system, base station and storage medium of a massive MIMO system.

背景技术Background technique

近年来移动通信发生了巨大变化,智能终端设备的迅速普及、无线数据通信业务的快速加入以及新兴产业的物联网技术新服务等,都对数据通信发明了更高的要求。为解决频谱资源短缺的问题,大规模多输入多输出(Multiple Input Multiple Output,MIMO)技术随之出现。MIMO技术是通过在基站上配置上百根天线同时向多位用户提供服务。相比传统小规模MIMO系统,大规模MIMO系统显著提升了频谱效率,可以满足高数据传输速率、稳定连接和低延迟等方面的要求。具有极高能效、频谱利用率的大规模MIMO技术被广泛认为是下一代通信网络架构的关键技术。In recent years, great changes have taken place in mobile communication. The rapid popularization of intelligent terminal equipment, the rapid addition of wireless data communication services, and the new services of Internet of Things technology in emerging industries have created higher requirements for data communication. In order to solve the problem of shortage of spectrum resources, a large-scale multiple input multiple output (Multiple Input Multiple Output, MIMO) technology emerges. MIMO technology provides services to multiple users at the same time by configuring hundreds of antennas on the base station. Compared with traditional small-scale MIMO systems, massive MIMO systems significantly improve spectral efficiency and can meet the requirements of high data transmission rates, stable connections, and low latency. Massive MIMO technology with extremely high energy efficiency and spectrum utilization is widely regarded as the key technology of next-generation communication network architecture.

MIMO带来显著性能增益的同时,大规模天线数目、高维度的信道矩阵使得检测信号算法复杂度呈指数增长,给接收端的信号检测带来了巨大挑战。主动禁忌搜索算法是一种启发式算法,具有复杂度低,性能良好的特点,在Massive MIMO信号检测中有很大的适用性,是大规模MIMO检测的研究热点,但是在高阶调制下仍存在着性能不佳的问题。While MIMO brings significant performance gains, the large number of antennas and the high-dimensional channel matrix increase the complexity of the detection signal algorithm exponentially, which brings great challenges to the signal detection at the receiving end. The active tabu search algorithm is a heuristic algorithm with low complexity and good performance. It has great applicability in Massive MIMO signal detection and is a research hotspot in massive MIMO detection. There is a problem with poor performance.

发明内容SUMMARY OF THE INVENTION

鉴于以上所述现有技术的缺点,本发明的目的在于提供一种大规模MIMO系统信号检测方法、系统、基站及存储介质,用于解决现有技术中的大规模MIMO系统信号检测中RTS算法在高阶调制系统中性能不佳的技术问题。In view of the above-mentioned shortcomings of the prior art, the purpose of the present invention is to provide a massive MIMO system signal detection method, system, base station and storage medium for solving the RTS algorithm in the massive MIMO system signal detection in the prior art Technical issues with poor performance in higher order modulation systems.

为实现上述目的及其他相关目的,本发明提供一种大规模MIMO系统信号检测方法,包括:To achieve the above object and other related objects, the present invention provides a signal detection method for a massive MIMO system, including:

对第一初始信号进行主动禁忌搜索检测,以获取初始估计向量;Perform active tabu search detection on the first initial signal to obtain an initial estimated vector;

根据所述初始估计向量消除所述第一初始信号中的干扰信号,以获取第二初始信号;cancel the interference signal in the first initial signal according to the initial estimation vector to obtain a second initial signal;

对所述第二初始信号进行消息传递检测以获取输出向量估计;performing message passing detection on the second initial signal to obtain an output vector estimate;

根据所述输出向量估计对符号向量进行重构,以获取符号向量重构值;Reconstructing the symbol vector according to the output vector estimate to obtain a symbol vector reconstruction value;

将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作,迭代结束后的最终的所述符号向量重构值作为检测结果输出。The iterative operation is performed using the reconstructed value of the symbol vector as the input of the active tabu search, and the final reconstructed value of the symbol vector after the iteration is output as the detection result.

在一可选实施例中,所述根据所述初始估计向量消除所述第一初始信号中的干扰信号,以获取第二初始信号的步骤包括:In an optional embodiment, the step of canceling the interference signal in the first initial signal according to the initial estimation vector to obtain the second initial signal includes:

根据所述初始估计向量获取输出比特值;Obtain output bit values according to the initial estimation vector;

根据所述输出比特值获取信号间总干扰;Obtain the total interference between signals according to the output bit value;

用所述第一初始信号减去所述信号间总干扰,以获取所述第二初始信号。The total inter-signal interference is subtracted from the first initial signal to obtain the second initial signal.

在一可选实施例中,根据下式来获取所述输出比特值

Figure BDA0002988681230000021
In an optional embodiment, the output bit value is obtained according to the following formula
Figure BDA0002988681230000021

Figure BDA0002988681230000022
Figure BDA0002988681230000022

其中,

Figure BDA0002988681230000023
Q是星座图中点的个数,K为MIMO系统在上行链路中用户的个数,
Figure BDA0002988681230000024
为初始估计向量
Figure BDA0002988681230000025
的第i个分量。in,
Figure BDA0002988681230000023
Q is the number of points in the constellation diagram, K is the number of users in the uplink of the MIMO system,
Figure BDA0002988681230000024
is the initial estimate vector
Figure BDA0002988681230000025
the ith component of .

在一可选实施例中,根据所述输出比特值获取信号间总干扰的步骤中,所述信号间总干扰

Figure BDA0002988681230000026
的表达式如下:In an optional embodiment, in the step of obtaining the total inter-signal interference according to the output bit value, the total inter-signal interference
Figure BDA0002988681230000026
The expression is as follows:

Figure BDA0002988681230000027
Figure BDA0002988681230000027

其中,

Figure BDA0002988681230000028
H为信道增益矩阵。in,
Figure BDA0002988681230000028
H is the channel gain matrix.

在一可选实施例中,根据所述输出向量估计对符号向量进行重构,以获取符号向量重构值的步骤中,所述符号向量重构值

Figure BDA0002988681230000029
的表达式如下:In an optional embodiment, in the step of reconstructing the symbol vector according to the output vector estimate to obtain the reconstructed value of the symbol vector, the reconstructed value of the symbol vector is
Figure BDA0002988681230000029
The expression is as follows:

Figure BDA00029886812300000210
Figure BDA00029886812300000210

其中,

Figure BDA00029886812300000211
为输出向量估计。in,
Figure BDA00029886812300000211
Estimated for the output vector.

在一可选实施例中,将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作的步骤中,迭代次数大于等于三次。In an optional embodiment, in the step of performing the iterative operation using the reconstructed value of the symbol vector as the input of the active tabu search, the number of iterations is greater than or equal to three times.

在一可选实施例中,将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作的步骤中,迭代次数为三次。In an optional embodiment, in the step of performing the iterative operation using the reconstructed value of the symbol vector as the input of the active tabu search, the number of iterations is three times.

为实现上述目的及其他相关目的,本发明还提供一种大规模MIMO系统信号检测系统,包括:To achieve the above object and other related objects, the present invention also provides a massive MIMO system signal detection system, including:

主动禁忌搜索模块,用于对第一初始信号进行主动禁忌搜索检测,以获取初始估计向量;an active tabu search module, which is used to perform active tabu search detection on the first initial signal to obtain an initial estimated vector;

干扰信号消除模块,用于根据所述初始估计向量消除所述第一初始信号中的干扰信号,以获取第二初始信号;an interference signal cancellation module, configured to cancel the interference signal in the first initial signal according to the initial estimation vector to obtain a second initial signal;

向量估计获取模块,用于对所述第二初始信号进行消息传递检测以获取输出向量估计;a vector estimation obtaining module, configured to perform message passing detection on the second initial signal to obtain an output vector estimation;

符号向量重构模块,用于根据所述输出向量估计对符号向量进行重构,以获取符号向量重构值;a symbol vector reconstruction module, configured to reconstruct the symbol vector according to the output vector estimation to obtain a symbol vector reconstruction value;

迭代输出模块,用于将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作,迭代结束后的最终的所述符号向量重构值作为检测结果输出。The iterative output module is configured to perform an iterative operation using the reconstructed value of the symbol vector as an input of the active tabu search, and output the final reconstructed value of the symbol vector after the iteration is completed as a detection result.

为实现上述目的及其他相关目的,本发明还提供一种基站,所述基站包括:To achieve the above object and other related objects, the present invention also provides a base station, the base station includes:

基站本体,设置于基本本体上的若干接收天线和控制单元;The base station body, a plurality of receiving antennas and control units arranged on the base body;

其中,所述接收天线用于接收MIMO系统的用户发送的符号向量,并将所述符号向量作为第一初始信号传输给所述控制单元进行信号检测;所述控制单元包括相互耦合的处理器和存储器,所述存储器存储有程序指令,当所述存储器存储的程序指令被所述处理器执行时实现上述任意一项所述的大规模MIMO系统信号检测方法。Wherein, the receiving antenna is used to receive a symbol vector sent by a user of the MIMO system, and transmit the symbol vector as a first initial signal to the control unit for signal detection; the control unit includes a mutually coupled processor and A memory, where program instructions are stored in the memory, and when the program instructions stored in the memory are executed by the processor, any one of the above-mentioned methods for detecting a signal in a massive MIMO system is implemented.

为实现上述目的及其他相关目的,本发明还提供一种存储介质,包括程序,当所述程序在计算机上运行时,使得所述计算机执行所述的大规模MIMO系统信号检测方法。In order to achieve the above object and other related objects, the present invention also provides a storage medium including a program, when the program runs on a computer, the computer enables the computer to execute the signal detection method for a massive MIMO system.

本发明的大规模MIMO系统信号检测方法、系统、基站及存储介质,提出了一种主动禁忌搜索与消息传递检测联合算法用于高阶调制下大规模MIMO信号检测,将MPD算法结合到RTS算法中,使用MPD算法来纠正符号的最后一位的错误,可以改善RTS算法在高阶调制下性能不佳的情况,且拥有更低的复杂度。The present invention provides a massive MIMO system signal detection method, system, base station and storage medium, proposes a combined active tabu search and message passing detection algorithm for massive MIMO signal detection under high-order modulation, and combines MPD algorithm with RTS algorithm In , the MPD algorithm is used to correct the error of the last bit of the symbol, which can improve the poor performance of the RTS algorithm under high-order modulation, and has a lower complexity.

附图说明Description of drawings

图1显示为本发明的大规模MIMO系统模型图。FIG. 1 shows a model diagram of a massive MIMO system of the present invention.

图2显示为本发明的一实施例的大规模MIMO系统信号检测方法流程图。FIG. 2 is a flowchart of a signal detection method for a massive MIMO system according to an embodiment of the present invention.

图3显示为本发明的一实施例的大规模MIMO系统信号检测系统的结构框图。FIG. 3 is a structural block diagram of a signal detection system of a massive MIMO system according to an embodiment of the present invention.

图4显示为本发明的一实施例的控制单元的结构框图。FIG. 4 is a structural block diagram of a control unit according to an embodiment of the present invention.

图5显示为4-QAM调制下不同天线数目时RTS算法性能曲线图。Figure 5 shows the performance curve of the RTS algorithm with different numbers of antennas under 4-QAM modulation.

图6显示为16-QAM与64-QAM调制下RTS算法性能比较图。Figure 6 shows the performance comparison of the RTS algorithm under 16-QAM and 64-QAM modulation.

图7显示为本发明的RTS-MPD算法与RTS算法的检测性能对比图。FIG. 7 is a graph showing the comparison of detection performance between the RTS-MPD algorithm of the present invention and the RTS algorithm.

具体实施方式Detailed ways

以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。The embodiments of the present invention are described below through specific specific examples, and those skilled in the art can easily understand other advantages and effects of the present invention from the contents disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention.

请参阅图1-7。需要说明的是,本实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图示中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。See Figures 1-7. It should be noted that the diagrams provided in this embodiment are only to illustrate the basic concept of the present invention in a schematic way, so the diagrams only show the components related to the present invention rather than the number, shape and the number of components in the actual implementation. For dimension drawing, the type, quantity and proportion of each component can be changed at will in actual implementation, and the component layout may also be more complicated.

图1示出了本发明的大规模MIMO(multiple-input multiple-output,多输入多输出)系统模型图。考虑大规模MIMO系统在上行链路有K个用户1,每个用户有一根发射天线11,与一个装备N个接收天线21的基站2进行通信,N的取值范围在数十到数百之间。α=K/N是系统的加载因子(α≤1)。系统模型如图1所示。FIG. 1 shows a model diagram of a massive MIMO (multiple-input multiple-output, multiple-input multiple-output) system of the present invention. Consider that a massive MIMO system has K users 1 in the uplink, each user has a transmit antenna 11, and communicates with a base station 2 equipped with N receive antennas 21, and the value of N ranges from tens to hundreds between. α=K/N is the loading factor of the system (α≤1). The system model is shown in Figure 1.

Figure BDA0002988681230000041
表示在第t个信道的信道增益矩阵,
Figure BDA0002988681230000042
表示从第j个用户天线(发射天线11)到第i个基站天线(接收天线21)的复数信道增益,信道增益
Figure BDA0002988681230000043
是服从均值为0方差为
Figure BDA0002988681230000044
Figure BDA0002988681230000045
模型由于路径损耗等导致来自用户j的接收功率不平衡,
Figure BDA0002988681230000046
对应于完全功率控制情况。用
Figure BDA0002988681230000047
表示调制符号向量传输在第t个信道中,其中
Figure BDA0002988681230000048
的第j个元素表示通过第j个用户传输的调制符号。假设完全同步,在基站2第t个信道接收向量为
Figure BDA0002988681230000049
可以写为:use
Figure BDA0002988681230000041
represents the channel gain matrix at the t-th channel,
Figure BDA0002988681230000042
Represents the complex channel gain from the jth user antenna (transmitting antenna 11) to the ith base station antenna (receiving antenna 21), the channel gain
Figure BDA0002988681230000043
is subject to the mean of 0 and the variance is
Figure BDA0002988681230000044
Figure BDA0002988681230000045
The model causes the received power from user j to be unbalanced due to path loss, etc.,
Figure BDA0002988681230000046
Corresponds to the full power control case. use
Figure BDA0002988681230000047
represents that the modulation symbol vector is transmitted in the t-th channel, where
Figure BDA0002988681230000048
The jth element of is the modulation symbol transmitted by the jth user. Assuming complete synchronization, the received vector on the t-th channel of base station 2 is
Figure BDA0002988681230000049
can be written as:

Figure BDA00029886812300000410
Figure BDA00029886812300000410

其中

Figure BDA00029886812300000411
为噪声向量,为了方便删除信道索引,(1)实数系统模型写为:in
Figure BDA00029886812300000411
is the noise vector, in order to conveniently delete the channel index, (1) the real number system model is written as:

y=Hx+n (2)y=Hx+n (2)

其中

Figure BDA00029886812300000412
Figure BDA00029886812300000413
Figure BDA00029886812300000414
分别表示实数和复数部分,同时,
Figure BDA00029886812300000415
对于
Figure BDA00029886812300000416
Figure BDA00029886812300000417
其中
Figure BDA00029886812300000418
是脉冲幅度调制(PluseAmplitude Modulation,PAM)字母表,表示在复数星座图Θ中同向或正交部分的点。由于
Figure BDA00029886812300000419
传输符号的星座图大小从Q减小到
Figure BDA00029886812300000420
其中,Q为星座图中点的个数。现在一个复数N×K MIMO复数星座图Θ等价于一个实数2N×2K MIMO实数星座图
Figure BDA00029886812300000421
n是噪声向量服从独立同分布
Figure BDA00029886812300000422
每根接收天线21的平均接收信噪比
Figure BDA00029886812300000423
其中Es表示传输符号的平均能量。in
Figure BDA00029886812300000412
Figure BDA00029886812300000413
and
Figure BDA00029886812300000414
represent the real and complex parts, respectively, and, at the same time,
Figure BDA00029886812300000415
for
Figure BDA00029886812300000416
Figure BDA00029886812300000417
in
Figure BDA00029886812300000418
is the Pulse Amplitude Modulation (PAM) alphabet, representing points in the in-direction or quadrature portion of the complex constellation Θ. because
Figure BDA00029886812300000419
The constellation size of the transmitted symbols is reduced from Q to
Figure BDA00029886812300000420
Among them, Q is the number of points in the constellation diagram. Now a complex N×K MIMO complex constellation Θ is equivalent to a real 2N×2K MIMO real constellation
Figure BDA00029886812300000421
n is the noise vector that is independent and identically distributed
Figure BDA00029886812300000422
Average received signal-to-noise ratio of each receive antenna 21
Figure BDA00029886812300000423
where Es represents the average energy of the transmitted symbols.

下面将分别对主动禁忌搜索(Reactive Tabu Search,简称RTS)算法、消息传递检测(Message Passing Detection,简称MPD)算法及本发明的RTS-MPD联合算法的原理分别进行说明。The principles of the Reactive Tabu Search (RTS for short) algorithm, the Message Passing Detection (MPD for short) algorithm and the RTS-MPD joint algorithm of the present invention will be respectively described below.

A、RTS检测算法A. RTS detection algorithm

对于实数系统模型,最大似然(Maximum likelihood,简称ML)检测模型,可表示为For the real number system model, the maximum likelihood (ML) detection model can be expressed as

Figure BDA0002988681230000051
Figure BDA0002988681230000051

当传输比特的先验概率相同时,似然检测等价于最大后验检测,可表示为:When the prior probability of the transmitted bits is the same, the likelihood detection is equivalent to the maximum a posteriori detection, which can be expressed as:

Figure BDA0002988681230000052
Figure BDA0002988681230000052

其中式(3)和式(4)精确计算的复杂度是K的指数级。Among them, the exact calculation complexity of formula (3) and formula (4) is the exponential level of K.

由公式(3),最大似然模型可以写成如下形式:From formula (3), the maximum likelihood model can be written in the following form:

Figure BDA0002988681230000053
Figure BDA0002988681230000053

其中,

Figure BDA0002988681230000054
f(x)为传输信号向量的代价函数。in,
Figure BDA0002988681230000054
f(x) is the cost function of the transmitted signal vector.

假设g(m)代表m次迭代以后的最优的解向量(ML)代价函数值最小,Lrep代表两个重复的解向量之间的平均长度(即平均迭代次数),Nrep代表解向量重复的次数,P代表禁忌周期,lflag∈{0,1}是局部最小的标志,用来判断在本次迭代中是否到达局部最小点。Suppose g (m) represents the optimal solution vector (ML) after m iterations with the smallest cost function value, L rep represents the average length between two repeated solution vectors (ie, the average number of iterations), and N rep represents the solution vector The number of repetitions, P represents the taboo period, and lflag∈{0,1} is the sign of the local minimum, which is used to judge whether the local minimum is reached in this iteration.

RTS算法从一个初始解向量开始,记为x(0),其可以是已知检测器ZF(ZeroForcing,迫零)/MMSE(Minimum Mean Squared Error,最小均方误差)的输出,也可以是随机生成的。令g(0)=x(0),Lrep=0,Nrep=0,p=p0(p0是初始禁忌周期值,是一个正整数),且禁忌矩阵T=0,定义

Figure BDA0002988681230000055
并计算出yMF和R的值。下面给出了在每一次迭代时所需要的步骤,考虑RTS算法检测过程的第m(m>0)次迭代:The RTS algorithm starts with an initial solution vector, denoted as x (0) , which can be the output of a known detector ZF (ZeroForcing)/MMSE (Minimum Mean Squared Error), or a random Generated. Let g (0) = x (0) , L rep = 0, N rep = 0, p = p 0 (p 0 is the initial taboo period value, which is a positive integer), and the taboo matrix T=0, define
Figure BDA0002988681230000055
And calculate the value of y MF and R. The steps required at each iteration are given below, considering the mth (m>0) iteration of the RTS algorithm detection process:

第一步:初始化lflag=0。设

Figure BDA0002988681230000056
e=Z(m)(u,v)-x(m),则x(m)的KM(M为符号邻居的数目)个向量邻居为Z(m)(u,v),u=1,2,…,K,v=1,2,…,M的ML代价函数值f(z(m)(u,v))计算如下:Step 1: Initialize lflag=0. Assume
Figure BDA0002988681230000056
e=Z (m) (u,v)-x (m) , then the KM (M is the number of symbol neighbors) vector neighbors of x (m) is Z (m) (u,v), u=1, The ML cost function value f(z (m) (u,v)) for 2,…,K, v=1,2,…,M is calculated as follows:

Figure BDA0002988681230000057
Figure BDA0002988681230000057

其中,

Figure BDA0002988681230000058
in,
Figure BDA0002988681230000058

(u1,v1)=arg minu,vf(z(m)(u,v))=arg minu,v(u,v) (7)(u 1 ,v 1 )=arg min u,v f(z (m) (u,v))=arg min u,v (u,v) (7)

则,判断移动move(u1,v1)可以发生的条件如下:Then, the conditions for judging that the move move(u 1 , v 1 ) can occur are as follows:

f(z(m)(u1,v1))<f(g(m)) (8)f(z (m) (u 1 ,v 1 ))<f(g (m) ) (8)

T((u1-1)M+q,v1)=0 (9)T((u 1 -1)M+q,v 1 )=0 (9)

只要满足上面一个条件中的任意一个,移动move(u1,v1)就可以发生,即执行当前解向量x(m)The move move(u 1 , v 1 ) can occur as long as any one of the above conditions is satisfied, ie the current solution vector x (m) is executed.

向第(u1,v1)个向量邻居z(m)(u1,v1)的转移。上式(9)中q是根据

Figure BDA0002988681230000061
推算出来的。如果移动move(u1,v1)对应的向量邻居不满足式(8)中的条件,且禁忌矩阵T也不满足式(9)中的条件,则用(u2,v2)进行条件判断,其中Transition to the (u 1 , v 1 )-th vector neighbor z (m) (u 1 , v 1 ). In the above formula (9), q is based on
Figure BDA0002988681230000061
inferred. If the vector neighbor corresponding to move move(u 1 , v 1 ) does not satisfy the condition in equation (8), and the taboo matrix T does not satisfy the condition in equation (9), then use (u 2 , v 2 ) to condition judgment, which

Figure BDA0002988681230000062
Figure BDA0002988681230000062

如果move(u2,v2)也不能够发生,则继续使用(uk,vk),k=3,4,…,kM进行条件判断,直到x(m)能够向一个向量邻居方向转移为止。如果KM个待转移方向全部被禁忌,则找出禁忌矩阵T中的最小值,并使禁忌矩阵T中的所有元素均减去这个最小值,然后重新重复上述过程进行是否可以转移的条件判断。假设(u′,v′)所对应向量邻居z(m)(u′,v′)的代价函数值是当前最小,且移动move(u′,v′)可以发生,则有If move(u 2 , v 2 ) cannot occur, continue to use (u k , v k ), k=3,4,...,kM to perform conditional judgment until x (m) can be transferred to a vector neighbor direction until. If all the KM directions to be transferred are taboo, find the minimum value in the taboo matrix T, subtract this minimum value from all elements in the taboo matrix T, and then repeat the above process again to judge whether the transfer can be performed. Assuming that the cost function value of the vector neighbor z (m) (u', v') corresponding to (u', v') is the current minimum, and the move (u', v') can occur, then we have

x(m+1)=z(m)(u′,v′) (11)x (m+1) = z (m) (u′,v′) (11)

第二步:对第一步获得的解向量进行解的重复性检验。如果新获得的解向量和之前迭代得到的解向量出现重复,则Nrep=Nrep+1,同时计算并更新lrep的值。禁忌周期P的值调整为P+1。但如果对于一个固定的β,当前P值超过βlrep的值,则令P=max(1,P-1),即Step 2: Check the repeatability of the solution on the solution vector obtained in the first step. If the newly obtained solution vector and the solution vector obtained by the previous iteration are repeated, N rep =N rep +1, and the value of l rep is calculated and updated at the same time. The value of the taboo period P is adjusted to P+1. But if for a fixed β, the current value of P exceeds the value of βl rep , then let P=max(1,P-1), that is

Figure BDA0002988681230000063
Figure BDA0002988681230000063

如果f(x(m+1))<f(g(m))则按照(13)和(14)对禁忌矩阵的某些禁忌值和最优解向量进行更新:If f(x (m+1) )<f(g (m) ), update some tabu values of the tabu matrix and the optimal solution vector according to (13) and (14):

T((u′-1)M+q′,v′)=T((u′-1)M+q″,v″)=0 (13)T((u′-1)M+q′,v′)=T((u′-1)M+q″,v″)=0 (13)

lflag=0,g(m+1)=x(m+1) (14)lflag=0,g (m+1) =x (m+1) (14)

否则,则按照式(15)和(16)进行更新:Otherwise, update according to equations (15) and (16):

T((u′-1)M+q′,v′)=T((u′-1)M+q″,v″)=P+1 (15)T((u′-1)M+q′,v′)=T((u′-1)M+q″,v″)=P+1 (15)

lflag=1,g(m+1)=g(m+1) (16)lflag=1,g (m+1) =g (m+1) (16)

第三步:根据式(17)更新禁忌矩阵T的值Step 3: Update the value of the taboo matrix T according to equation (17)

T=max{T-1,0} (17)T=max{T-1,0} (17)

同时根据式(18)跟新f(m)的值,如下At the same time, according to formula (18) and the new value of f (m) , as follows

Figure BDA0002988681230000071
Figure BDA0002988681230000071

其中Ru′表示R的第u′列。此时,判断是否满足终止条件,如果满足则终止迭代过程,并返回检测结果作为最终解向量;如果不满足终止条件,则跳到第一步继续执行上述三个步骤。where R u' denotes the u'th column of R. At this point, it is judged whether the termination condition is met, if so, the iteration process is terminated, and the detection result is returned as the final solution vector; if the termination condition is not met, skip to the first step and continue to execute the above three steps.

B、MPD检测算法B. MPD detection algorithm

消息传递检测(MPD)算法是一种利用信道硬化现象低复杂度的消息传递方法。The Message Passing Detection (MPD) algorithm is a low-complexity message passing method utilizing the channel hardening phenomenon.

用y=Hx+n乘以HT在除以N,我们可以得到Multiplying H T by y=Hx+n and dividing by N, we get

Z=Jx+V (19)Z=Jx+V (19)

其中Z=HTy/N,J=HTH/N,V=HTn/N。Z的第i个元素是

Figure BDA0002988681230000072
Figure BDA0002988681230000073
其中
Figure BDA0002988681230000074
Figure BDA0002988681230000075
gi被假设是一个服从均值为μi和方差为
Figure BDA0002988681230000076
高斯分布Wherein Z= HT y/N, J= HT H/N, V= HT n/N. The ith element of Z is
Figure BDA0002988681230000072
Figure BDA0002988681230000073
in
Figure BDA0002988681230000074
and
Figure BDA0002988681230000075
g i is assumed to be a subject with mean μ i and variance of
Figure BDA0002988681230000076
Gaussian distribution

Figure BDA0002988681230000077
Figure BDA0002988681230000077

Figure BDA0002988681230000078
Figure BDA0002988681230000078

其中

Figure BDA0002988681230000079
计算E(gi)和Var(gi),需要计算pj(sk)的后验概率,即xj成为sk∈θ的概率。符号xj的对数似然比定义为
Figure BDA00029886812300000710
其中in
Figure BDA0002988681230000079
To calculate E( gi ) and Var( gi ), we need to calculate the posterior probability of p j (s k ), that is, the probability that x j becomes s k ∈ θ. The log-likelihood ratio of the symbol x j is defined as
Figure BDA00029886812300000710
in

Figure BDA00029886812300000711
Figure BDA00029886812300000711

Figure BDA00029886812300000712
此外
Figure BDA00029886812300000712
also

Figure BDA00029886812300000713
Figure BDA00029886812300000713

最后,符号xj被确定为sk的最大概率。Finally, the symbol xj is determined as the maximum probability of sk .

C、RTS-MPD联合算法的基本原理C. The basic principle of RTS-MPD joint algorithm

定义x为传输向量、

Figure BDA00029886812300000714
为RTS检测器的输出、
Figure BDA00029886812300000715
Figure BDA00029886812300000716
调制的字母表并且x的取值来自此集合。考虑符号与比特之间的映射,这里将
Figure BDA0002988681230000081
的每一项的值写成成比特的线性组合:Define x as the transmission vector,
Figure BDA00029886812300000714
is the output of the RTS detector,
Figure BDA00029886812300000715
for
Figure BDA00029886812300000716
The modulated alphabet and the value of x comes from this set. Consider the mapping between symbols and bits, here will be
Figure BDA0002988681230000081
The value of each term is written as a linear combination of bits:

Figure BDA0002988681230000082
Figure BDA0002988681230000082

其中

Figure BDA0002988681230000083
Q是星座图中点的个数,
Figure BDA0002988681230000084
Figure BDA0002988681230000085
的第i个分量,根据前面介绍的RTS检测算法可知,RTS检测算法的输出的解向量为局部最小值。所以,
Figure BDA0002988681230000086
满足以下不等式:in
Figure BDA0002988681230000083
Q is the number of points in the constellation diagram,
Figure BDA0002988681230000084
for
Figure BDA0002988681230000085
The i-th component of , according to the RTS detection algorithm introduced above, the solution vector output by the RTS detection algorithm is a local minimum value. so,
Figure BDA0002988681230000086
Satisfy the following inequalities:

Figure BDA0002988681230000087
Figure BDA0002988681230000087

其中

Figure BDA0002988681230000088
ei表示特征矩阵的第i列。定义hi表示H的第i列、
Figure BDA0002988681230000089
Figure BDA00029886812300000810
则式(25)可表示为in
Figure BDA0002988681230000088
e i represents the i-th column of the feature matrix. Define h i to represent the i-th column of H,
Figure BDA0002988681230000089
Figure BDA00029886812300000810
The formula (25) can be expressed as

Figure BDA00029886812300000811
Figure BDA00029886812300000811

上式(26)中,fij是F的第i行j列的元素,在高信噪比的情况下忽略噪声,可进一步化为In the above formula (26), f ij is the element of the i-th row and j-column of F. In the case of high signal-to-noise ratio, ignoring the noise, it can be further transformed into

Figure BDA00029886812300000812
Figure BDA00029886812300000812

其中,fi为F的第i行,在瑞丽衰落信道中,fii服从自由度为2K,均值为N的卡方分布。当i≠j时,由中心极限定理可知fii逼近于0均值,方差为K/4的正态分布。式(27)可变成Among them, f i is the ith row of F, and in the Rayleigh fading channel, f ii obeys a chi-square distribution with 2K degrees of freedom and a mean of N. When i≠j, it can be known from the central limit theorem that f ii is close to 0 mean, and the variance is a normal distribution of K/4. Equation (27) can be transformed into

Figure BDA00029886812300000813
Figure BDA00029886812300000813

式(28)中:

Figure BDA00029886812300000814
Figure BDA00029886812300000815
则:In formula (28):
Figure BDA00029886812300000814
like
Figure BDA00029886812300000815
but:

Figure BDA00029886812300000816
Figure BDA00029886812300000816

图2示出了本发明的大规模MIMO系统信号检测方法的流程示意图,下面将结合附图来详细阐述本发明的实施例的大规模MIMO系统信号检测方法的实现流程。FIG. 2 shows a schematic flowchart of the signal detection method for a massive MIMO system according to the present invention, and the implementation flow of the signal detection method for a massive MIMO system according to an embodiment of the present invention will be described in detail below with reference to the accompanying drawings.

请参阅图2,所述大规模MIMO系统信号检测方法包括:Please refer to FIG. 2, the signal detection method of the massive MIMO system includes:

步骤S10、对第一初始信号进行主动禁忌搜索检测,以获取初始估计向量

Figure BDA00029886812300000817
具体地,单天线用户1通过发射天线11发射信号x,发射信号x经过调制,在经过信道传到基站2,生成接收信号y,作为所述第一初始信号y,将接收信号y经过RTS检测器进行RTS检测,以获取初始估计向量
Figure BDA00029886812300000818
Step S10, perform active tabu search detection on the first initial signal to obtain an initial estimated vector
Figure BDA00029886812300000817
Specifically, the single-antenna user 1 transmits a signal x through the transmit antenna 11, the transmit signal x is modulated, and then transmitted to the base station 2 through a channel to generate a received signal y, which is used as the first initial signal y, and the received signal y is detected by RTS the RTS detection to obtain the initial estimate vector
Figure BDA00029886812300000818

步骤S20、根据所述初始估计向量

Figure BDA00029886812300000819
消除所述第一初始信号中的干扰信号,以获取第二初始信号
Figure BDA00029886812300000820
具体地,可通过式(24)得到输出的比特值
Figure BDA00029886812300000821
i=1,…,2K,
Figure BDA00029886812300000822
由比特值
Figure BDA00029886812300000823
得出总的信号间干扰(也即信号间总干扰):
Figure BDA00029886812300000824
其中
Figure BDA00029886812300000825
Figure BDA0002988681230000091
H为信道增益矩阵;用所述接收信号y减去所述信号间总干扰
Figure BDA0002988681230000092
以获取所述第二初始信号
Figure BDA0002988681230000093
也即
Figure BDA0002988681230000094
Step S20, according to the initial estimation vector
Figure BDA00029886812300000819
Eliminate the interference signal in the first initial signal to obtain the second initial signal
Figure BDA00029886812300000820
Specifically, the output bit value can be obtained by formula (24)
Figure BDA00029886812300000821
i=1,...,2K,
Figure BDA00029886812300000822
by bit value
Figure BDA00029886812300000823
The total inter-signal interference (that is, the total inter-signal interference) is obtained:
Figure BDA00029886812300000824
in
Figure BDA00029886812300000825
Figure BDA0002988681230000091
H is the channel gain matrix; subtract the total interference between the signals from the received signal y
Figure BDA0002988681230000092
to obtain the second initial signal
Figure BDA0002988681230000093
that is
Figure BDA0002988681230000094

步骤S30、对所述第二初始信号

Figure BDA0002988681230000095
进行消息传递检测以获取输出向量估计
Figure BDA0002988681230000096
具体地,将步骤S20中获取的第二初始信号
Figure BDA0002988681230000097
输入到MPD检测模块进行MPD检测可以得到得到输出向量估计为
Figure BDA0002988681230000098
Step S30, for the second initial signal
Figure BDA0002988681230000095
Do message passing detection to get output vector estimates
Figure BDA0002988681230000096
Specifically, the second initial signal obtained in step S20 is
Figure BDA0002988681230000097
Input to the MPD detection module for MPD detection, the output vector can be estimated as
Figure BDA0002988681230000098

步骤S40、根据所述输出向量估计

Figure BDA0002988681230000099
对符号向量x进行重构,以获取符号向量重构值
Figure BDA00029886812300000910
其中,
Figure BDA00029886812300000911
Step S40, estimate according to the output vector
Figure BDA0002988681230000099
Reconstruct the symbol vector x to obtain the symbol vector reconstruction value
Figure BDA00029886812300000910
in,
Figure BDA00029886812300000911

步骤S50、将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作,直到满足迭代条件,迭代结束后的最终的所述符号向量重构值

Figure BDA00029886812300000912
作为检测结果输出。具体地,将步骤S40获取的符号向量重构值
Figure BDA00029886812300000913
作为RTS检测器的输入,重复执行执行步骤S10-S40的步骤,直到迭代结束,其中,迭代结束的依据是迭代次数达到三次或者大于三的任一整数值。Step S50, performing an iterative operation using the reconstructed value of the symbol vector as the input of the active tabu search, until the iteration conditions are met, and the final reconstructed value of the symbol vector after the iteration ends
Figure BDA00029886812300000912
output as the detection result. Specifically, reconstruct the value of the symbol vector obtained in step S40
Figure BDA00029886812300000913
As the input of the RTS detector, the steps of performing steps S10-S40 are repeatedly performed until the iteration ends, where the iteration is terminated based on the number of iterations reaching three or any integer value greater than three.

如图3所示,本发明的实施例还公开一种大规模MIMO系统信号检测系统3,所述大规模MIMO系统信号检测系统3包括主动禁忌搜索模块31、干扰信号消除模块32、向量估计获取模块33、符号向量重构模块34及迭代输出模块35。主动禁忌搜索模块,用于对第一初始信号进行主动禁忌搜索检测,以获取初始估计向量;该干扰信号消除模块31用于根据所述初始估计向量消除所述第一初始信号中的干扰信号,以获取第二初始信号;该向量估计获取模块32用于对所述第二初始信号进行消息传递检测以获取输出向量估计;该符号向量重构模块33用于根据所述输出向量估计对符号向量进行重构,以获取符号向量重构值;该迭代输出模块34用于将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作,迭代结束后的最终的所述符号向量重构值作为检测结果输出。As shown in FIG. 3 , an embodiment of the present invention further discloses a massive MIMO system signal detection system 3 . The massive MIMO system signal detection system 3 includes an active taboo search module 31 , an interference signal elimination module 32 , and a vector estimation acquisition module 31 . module 33 , symbol vector reconstruction module 34 and iterative output module 35 . an active tabu search module, configured to perform active tabu search detection on the first initial signal to obtain an initial estimation vector; the interference signal elimination module 31 is used to eliminate interference signals in the first initial signal according to the initial estimation vector, to obtain the second initial signal; the vector estimation obtaining module 32 is used to perform message passing detection on the second initial signal to obtain an output vector estimate; the symbol vector reconstruction module 33 is used to estimate the symbol vector according to the output vector estimation Reconstruction is performed to obtain the reconstruction value of the symbol vector; the iterative output module 34 is used to perform an iterative operation using the reconstructed value of the symbol vector as the input of the active tabu search, and the final reconstructed value of the symbol vector after the iteration ends output as the detection result.

需要说明的是,本发明的大规模MIMO系统信号检测系统3是与上述大规模MIMO系统信号检测方法相对应的系统,大规模MIMO系统信号检测系统3中的功能模块分别对应大规模MIMO系统信号检测方法中的相应步骤。本发明的大规模MIMO系统信号检测系统3可与大规模MIMO系统信号检测方法相互相配合实施。本发明的大规模MIMO系统信号检测方法中提到的相关技术细节在大规模MIMO系统信号检测系统3中依然有效,为了减少重复,这里不再赘述。相应地,本发明的大规模MIMO系统信号检测系统3中提到的相关技术细节也可应用在上述大规模MIMO系统信号检测方法中。It should be noted that the massive MIMO system signal detection system 3 of the present invention is a system corresponding to the above-mentioned massive MIMO system signal detection method, and the functional modules in the massive MIMO system signal detection system 3 respectively correspond to massive MIMO system signals corresponding steps in the detection method. The massive MIMO system signal detection system 3 of the present invention can be implemented in cooperation with the massive MIMO system signal detection method. The related technical details mentioned in the massive MIMO system signal detection method of the present invention are still valid in the massive MIMO system signal detection system 3, and are not repeated here in order to reduce repetition. Correspondingly, the relevant technical details mentioned in the massive MIMO system signal detection system 3 of the present invention can also be applied to the above-mentioned massive MIMO system signal detection method.

需要说明的是,上述的各功能模块,实际实现时可以全部或部分集成到一个物理实体上,也可以物理上分开。且这些单元可以全部以软件通过处理元件调用的形式实现;也可以全部以硬件的形式实现;还可以部分单元通过处理元件调用软件的形式实现,部分单元通过硬件的形式实现。此外这些单元全部或部分可以集成在一起,也可以独立实现。这里所述的处理元件可以是一种集成电路,具有信号的处理能力。在实现过程中,上述方法的各步骤或以上各个模块可以通过处理器41元件中的硬件的集成逻辑电路或者软件形式的指令完成。It should be noted that the above-mentioned functional modules may be integrated into one physical entity in whole or in part during actual implementation, or may be physically separated. And these units can all be implemented in the form of software calling through processing elements; they can also all be implemented in hardware; some units can also be implemented in the form of calling software through processing elements, and some units can be implemented in hardware. In addition, all or part of these units can be integrated together, and can also be implemented independently. The processing element described here may be an integrated circuit with signal processing capability. In the implementation process, each step of the above-mentioned method or each of the above-mentioned modules can be completed by an integrated logic circuit of hardware in the element of the processor 41 or an instruction in the form of software.

需要说明的是,如图4所示,本发明的大规模MIMO系统信号检测方法还可以通过一设置于基站(当然也可以是设置于其他主体上)上的控制单元4实现,所述控制单元4包括相互连接的存储器43和处理器41,所述存储器43存储有程序指令,该程序指令被所述处理器41执行时实现上述的大规模MIMO系统信号检测方法。需要说明的是,当需要和外部进行通信时,所述控制单元4还包括通信器42,所述通信器42与所述处理器41连接。It should be noted that, as shown in FIG. 4 , the signal detection method of the massive MIMO system of the present invention can also be implemented by a control unit 4 arranged on the base station (of course, it can also be arranged on other main bodies). 4 includes a memory 43 and a processor 41 that are connected to each other, the memory 43 stores program instructions, and when the program instructions are executed by the processor 41, the above-mentioned massive MIMO system signal detection method is implemented. It should be noted that, when communication with the outside is required, the control unit 4 further includes a communicator 42 , and the communicator 42 is connected to the processor 41 .

上述的处理器41可以是通用处理器,包括中央处理器(CentralProcessingUnit,简称CPU)、网络处理器(NetworkProcessor,简称NP)等;还可以是数字信号处理器(DigitalSignal Processing,简称DSP)、专用集成电路(Application Specific IntegratedCircuit,简称ASIC)、现场可编程门阵列(Field-Programmable Gate Array,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件;上述的存储器43可能包含随机存取存储器(Random Access Memory,简称RAM),也可能还包括非易失性存储器(Non-volatile Memory),例如至少一个磁盘存储器。The above-mentioned processor 41 may be a general-purpose processor, including a central processing unit (Central Processing Unit, referred to as CPU), a network processor (Network Processor, referred to as NP), etc.; may also be a digital signal processor (Digital Signal Processing, referred to as DSP), dedicated integrated Circuit (Application Specific Integrated Circuit, referred to as ASIC), Field-Programmable Gate Array (Field-Programmable Gate Array, referred to as FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components; the above-mentioned memory 43 may contain random The access memory (Random Access Memory, RAM for short) may also include a non-volatile memory (Non-volatile Memory), such as at least one disk memory.

需要说明的是,上述存储器43中的控制单元4可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,电子设备,或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。It should be noted that the control unit 4 in the above-mentioned memory 43 may be implemented in the form of software functional units and may be stored in a computer-readable storage medium when sold or used as an independent product. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including several The instructions are used to cause a computer device (which may be a personal computer, an electronic device, or a network device, etc.) to execute all or part of the steps of the methods of various embodiments of the present invention.

本发明还可以提供一种存储介质,其存储有程序,该程序被处理器41执行时实现上述的大规模MIMO系统信号检测方法;所述存储介质包括所有形式的非易失性存储器、介质和存储器设备,包括例如:半导体存储器设备,例如EPROM、EEPROM和闪存设备;磁盘,例如内部硬盘或可移动盘;磁光盘;以及CD-ROM和DVD-ROM盘。The present invention can also provide a storage medium, which stores a program, and when the program is executed by the processor 41, implements the above-mentioned massive MIMO system signal detection method; the storage medium includes all forms of non-volatile memory, medium and Memory devices, including, for example: semiconductor memory devices, such as EPROM, EEPROM, and flash memory devices; magnetic disks, such as internal hard disks or removable disks; magneto-optical disks; and CD-ROM and DVD-ROM disks.

为了验证本发明的大规模MIMO系统信号检测方法,发明人进行了仿真实验。本实施例给出了Matlab蒙特卡洛仿真结果。仿真参数:调制方式选择QAM;天线规模为16×16、32×32、64×64。传输信道是瑞丽衰落信道,信道特性时变,噪声为加性高斯白噪声。本节主要研究不同调制阶数和不同的天线配置下RTS与RTS-MPD算法误码率的对比。In order to verify the signal detection method of the massive MIMO system of the present invention, the inventor has conducted a simulation experiment. This embodiment presents the Matlab Monte Carlo simulation results. Simulation parameters: QAM is selected as the modulation mode; the antenna scale is 16×16, 32×32, 64×64. The transmission channel is a Rayleigh fading channel, the channel characteristics are time-varying, and the noise is additive white Gaussian noise. This section mainly studies the comparison of bit error rates between RTS and RTS-MPD algorithms under different modulation orders and different antenna configurations.

图5绘制了在调制方式为4QAM的情况下,天线配置分别为16×16,32×32,64×64情况下RTS检测算法的性能。图中还加入了SISO AWGN性能曲线进行对比。从图中可以看出RTS在低阶调制情况下具有良好的性能,且随着天线数目的增加,其性能提升明显,接近SISO AWGN。所以RTS是适用于调制阶数低,天线数目多的大规模MIMO系统中的信号检测。Figure 5 plots the performance of the RTS detection algorithm when the modulation scheme is 4QAM and the antenna configurations are 16×16, 32×32, and 64×64, respectively. The SISO AWGN performance curve is also added to the figure for comparison. It can be seen from the figure that RTS has good performance in the case of low-order modulation, and with the increase of the number of antennas, its performance improves significantly, which is close to SISO AWGN. Therefore, RTS is suitable for signal detection in massive MIMO systems with a low modulation order and a large number of antennas.

图6绘制了在高阶调制下RTS算法性能曲线图,由图5可知在4-QAM(QuadratureAmplitude Modulation,正交振幅调制)调制情况下,RTS算法在信噪比大约为11db时,误码率就可达到10-3,但是从图6中可以看出,在16-QAM高阶调制下,RTS算法性能不佳,在信噪比高达30时才能达到10-3;在图6中,还可以看出当调制阶数变得更高时,RTS算法的检测性能进一步明显变差,RTS算法在高阶调制下性能表现的非常不好。Figure 6 shows the performance curve of the RTS algorithm under high-order modulation. It can be seen from Figure 5 that in the case of 4-QAM (Quadrature Amplitude Modulation, quadrature amplitude modulation) modulation, the bit error rate of the RTS algorithm is about 11db when the signal-to-noise ratio is about 11db. It can reach 10 -3 , but it can be seen from Figure 6 that under 16-QAM high-order modulation, the performance of the RTS algorithm is not good, and it can only reach 10 -3 when the signal-to-noise ratio is as high as 30; in Figure 6, also It can be seen that when the modulation order becomes higher, the detection performance of the RTS algorithm further deteriorates significantly, and the performance of the RTS algorithm under high-order modulation is very poor.

仿真试验中各试验参数见下为了简单直观验证所介绍的RTS-MPD算法检测的性能,给出了RTS算法与RTS-MPD算法分别在16QAM和64QAM调制情况下的实验仿真对比图,如图7所示。可以从图中看出,在16QAM调制的情况下,当BER=10-3时,RTS-MPD联合算法优于RTS算法性能,大约提升有1db的性能;在64QAM调制的情况,当BER=10-3时,RTS-MPD联合算法优于RTS算法性能,性能提升更加明显,达到2db左右。The test parameters in the simulation test are shown below. In order to simply and intuitively verify the detection performance of the RTS-MPD algorithm introduced, the experimental simulation comparison diagrams of the RTS algorithm and the RTS-MPD algorithm under 16QAM and 64QAM modulation are given, as shown in Figure 7. shown. It can be seen from the figure that in the case of 16QAM modulation, when BER=10 -3 , the RTS-MPD joint algorithm is better than the RTS algorithm, and the performance is improved by about 1db; in the case of 64QAM modulation, when BER=10 When -3 , the RTS-MPD joint algorithm is better than the RTS algorithm, and the performance improvement is more obvious, reaching about 2db.

综上所述,本发明的大规模MIMO系统信号检测方法、系统、基站及存储介质,提出了一种主动禁忌搜索与消息传递检测联合算法用于高阶调制下大规模MIMO信号检测,将MPD算法结合到RTS算法中,使用MPD算法来纠正符号的最后一位的错误,可以改善RTS算法在高阶调制下性能不佳的情况,且拥有更低的复杂度。To sum up, the method, system, base station and storage medium for massive MIMO system signal detection of the present invention propose a joint algorithm of active tabu search and message passing detection for massive MIMO signal detection under high-order modulation. The algorithm is integrated into the RTS algorithm, and the MPD algorithm is used to correct the error of the last bit of the symbol, which can improve the poor performance of the RTS algorithm under high-order modulation, and has a lower complexity.

在本文的描述中,提供了许多特定细节,诸如部件和/或方法的实例,以提供对本发明实施例的完全理解。然而,本领域技术人员将认识到可以在没有一项或多项具体细节的情况下或通过其他设备、系统、组件、方法、部件、材料、零件等等来实践本发明的实施例。在其他情况下,未具体示出或详细描述公知的结构、材料或操作,以避免使本发明实施例的方面变模糊。In the description herein, numerous specific details are provided, such as examples of components and/or methods, in order to provide a thorough understanding of embodiments of the present invention. However, one skilled in the art will recognize that embodiments of the invention may be practiced without one or more of the specific details or with other devices, systems, assemblies, methods, components, materials, parts, and the like. In other instances, well-known structures, materials, or operations are not specifically shown or described in detail to avoid obscuring aspects of the embodiments of the invention.

在整篇说明书中提到“一个实施例(one embodiment)”、“实施例(anembodiment)”或“具体实施例(a specific embodiment)”意指与结合实施例描述的特定特征、结构或特性包括在本发明的至少一个实施例中,并且不一定在所有实施例中。因而,在整篇说明书中不同地方的短语“在一个实施例中(in one embodiment)”、“在实施例中(inan embodiment)”或“在具体实施例中(in a specific embodiment)”的各个表象不一定是指相同的实施例。此外,本发明的任何具体实施例的特定特征、结构或特性可以按任何合适的方式与一个或多个其他实施例结合。应当理解本文所述和所示的发明实施例的其他变型和修改可能是根据本文教导的,并将被视作本发明精神和范围的一部分。Reference throughout this specification to "one embodiment," "anembodiment," or "a specific embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment includes In at least one embodiment of the invention, and not necessarily in all embodiments. Thus, the phrases "in one embodiment", "in an embodiment" or "in a specific embodiment" are used in various places throughout the specification. Appearances are not necessarily referring to the same embodiment. Furthermore, the particular features, structures, or characteristics of any particular embodiment of the present invention may be combined in any suitable manner with one or more other embodiments. It should be understood that other variations and modifications of the embodiments of the invention described and illustrated herein are possible in light of the teachings herein and are to be considered part of the spirit and scope of the invention.

还应当理解还可以以更分离或更整合的方式实施附图所示元件中的一个或多个,或者甚至因为在某些情况下不能操作而被移除或因为可以根据特定应用是有用的而被提供。It should also be understood that one or more of the elements shown in the figures may also be implemented in a more discrete or integrated manner, or even removed as inoperable in certain circumstances or as may be useful according to a particular application. Provided.

另外,除非另外明确指明,附图中的任何标志箭头应当仅被视为示例性的,而并非限制。此外,除非另外指明,本文所用的术语“或”一般意在表示“和/或”。在术语因提供分离或组合能力是不清楚的而被预见的情况下,部件或步骤的组合也将视为已被指明。Additionally, any identifying arrows in the accompanying drawings should be regarded as illustrative only and not restrictive unless expressly indicated otherwise. In addition, the term "or" as used herein is generally intended to mean "and/or" unless stated otherwise. Combinations of components or steps will also be considered to have been specified where the term is foreseen because the ability to provide separation or combination is unclear.

如在本文的描述和在下面整篇权利要求书中所用,除非另外指明,“一个(a)”、“一个(an)”和“该(the)”包括复数参考物。同样,如在本文的描述和在下面整篇权利要求书中所用,除非另外指明,“在…中(in)”的意思包括“在…中(in)”和“在…上(on)”。As used in the description herein and throughout the claims below, "a (a)," "an (an)," and "the (the)" include plural references unless otherwise indicated. Likewise, as used in the description herein and throughout the claims below, unless otherwise specified, the meaning of "in" includes "in" and "on" .

本发明所示实施例的上述描述(包括在说明书摘要中所述的内容)并非意在详尽列举或将本发明限制到本文所公开的精确形式。尽管在本文仅为说明的目的而描述了本发明的具体实施例和本发明的实例,但是正如本领域技术人员将认识和理解的,各种等效修改是可以在本发明的精神和范围内的。如所指出的,可以按照本发明所述实施例的上述描述来对本发明进行这些修改,并且这些修改将在本发明的精神和范围内。The above description of illustrated embodiments of the present invention, including what is described in the Abstract, is not intended to be exhaustive or to limit the invention to the precise form disclosed herein. While specific embodiments of, and examples for, the invention are described herein for illustrative purposes only, various equivalent modifications are possible within the spirit and scope of the invention, as those skilled in the art will recognize and appreciate of. As indicated, these modifications may be made to the present invention in light of the foregoing description of the described embodiments of the present invention and are intended to be within the spirit and scope of the present invention.

本文已经在总体上将系统和方法描述为有助于理解本发明的细节。此外,已经给出了各种具体细节以提供本发明实施例的总体理解。然而,相关领域的技术人员将会认识到,本发明的实施例可以在没有一个或多个具体细节的情况下进行实践,或者利用其它装置、系统、配件、方法、组件、材料、部分等进行实践。在其它情况下,并未特别示出或详细描述公知结构、材料和/或操作以避免对本发明实施例的各方面造成混淆。The systems and methods have generally been described herein with details that are helpful in understanding the invention. Furthermore, various specific details have been set forth in order to provide a general understanding of embodiments of the present invention. One skilled in the relevant art will recognize, however, that embodiments of the invention may be practiced without one or more of the specific details, or with other devices, systems, accessories, methods, components, materials, parts, etc. practice. In other instances, well-known structures, materials and/or operations have not been specifically shown or described in detail to avoid obscuring aspects of the embodiments of the invention.

因而,尽管本发明在本文已参照其具体实施例进行描述,但是修改自由、各种改变和替换亦在上述公开内,并且应当理解,在某些情况下,在未背离所提出发明的范围和精神的前提下,在没有对应使用其他特征的情况下将采用本发明的一些特征。因此,可以进行许多修改,以使特定环境或材料适应本发明的实质范围和精神。本发明并非意在限制到在下面权利要求书中使用的特定术语和/或作为设想用以执行本发明的最佳方式公开的具体实施例,但是本发明将包括落入所附权利要求书范围内的任何和所有实施例及等同物。因而,本发明的范围将只由所附的权利要求书进行确定。Thus, although the invention has been described herein with reference to specific embodiments thereof, freedom of modification, various changes and substitutions is within the above disclosure, and it should be understood that, in certain circumstances, without departing from the scope and scope of the proposed invention, Some features of the present invention will be employed without the corresponding use of other features in the spirit of the present invention. Therefore, many modifications may be made to adapt a particular environment or material to the essential scope and spirit of the invention. It is not intended that the invention be limited to the specific terms used in the following claims and/or the specific embodiments disclosed as the best modes contemplated for carrying out the invention, but the invention is to be included within the scope of the appended claims any and all embodiments and equivalents within. Accordingly, the scope of the present invention should be determined only by the appended claims.

Claims (6)

1.一种大规模MIMO系统信号检测方法,其特征在于,包括:1. a massive MIMO system signal detection method, is characterized in that, comprises: 对第一初始信号进行主动禁忌搜索检测,以获取初始估计向量;Perform active tabu search detection on the first initial signal to obtain an initial estimated vector; 根据所述初始估计向量获取输出比特值;Obtain output bit values according to the initial estimation vector; 根据所述输出比特值获取信号间总干扰;Obtain the total interference between signals according to the output bit value; 用所述第一初始信号减去所述信号间总干扰,以获取第二初始信号;subtracting the total inter-signal interference from the first initial signal to obtain a second initial signal; 对所述第二初始信号进行消息传递检测以获取输出向量估计;performing message passing detection on the second initial signal to obtain an output vector estimate; 根据所述输出向量估计对符号向量进行重构,以获取符号向量重构值;Reconstructing the symbol vector according to the output vector estimate to obtain a symbol vector reconstruction value; 将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作,迭代结束后的最终的所述符号向量重构值作为检测结果输出;Performing an iterative operation using the symbol vector reconstruction value as the input of the active tabu search, and outputting the final symbol vector reconstruction value after the iteration ends as a detection result; 其中,根据下式来获取所述输出比特值
Figure FDA0003426212110000011
Wherein, the output bit value is obtained according to the following formula
Figure FDA0003426212110000011
Figure FDA0003426212110000012
Figure FDA0003426212110000012
其中,
Figure FDA0003426212110000013
Q是星座图中点的个数,K为MIMO系统在上行链路中用户的个数,
Figure FDA0003426212110000014
为初始估计向量
Figure FDA0003426212110000015
的第i个分量;
in,
Figure FDA0003426212110000013
Q is the number of points in the constellation diagram, K is the number of users in the uplink of the MIMO system,
Figure FDA0003426212110000014
is the initial estimate vector
Figure FDA0003426212110000015
the ith component of ;
根据所述输出比特值获取信号间总干扰的步骤中,所述信号间总干扰
Figure FDA0003426212110000016
的表达式如下:
In the step of obtaining the total inter-signal interference according to the output bit value, the total inter-signal interference
Figure FDA0003426212110000016
The expression is as follows:
Figure FDA0003426212110000017
Figure FDA0003426212110000017
其中,
Figure FDA0003426212110000018
H为信道增益矩阵。
in,
Figure FDA0003426212110000018
H is the channel gain matrix.
2.根据权利要求1所述的大规模MIMO系统信号检测方法,其特征在于,根据所述输出向量估计对符号向量进行重构,以获取符号向量重构值的步骤中,所述符号向量重构值
Figure FDA00034262121100000111
的表达式如下:
2. The signal detection method for a massive MIMO system according to claim 1, wherein, in the step of reconstructing the symbol vector according to the output vector estimate to obtain a reconstructed value of the symbol vector, the symbol vector repeats construct value
Figure FDA00034262121100000111
The expression is as follows:
Figure FDA0003426212110000019
Figure FDA0003426212110000019
其中,
Figure FDA00034262121100000110
为输出向量估计。
in,
Figure FDA00034262121100000110
Estimated for the output vector.
3.根据权利要求1或2所述的大规模MIMO系统信号检测方法,其特征在于,将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作的步骤中,迭代次数大于等于三次。3 . The signal detection method for a massive MIMO system according to claim 1 or 2 , wherein, in the step of performing an iterative operation using the symbol vector reconstruction value as the input of the active tabu search, the number of iterations is greater than or equal to three times. 4 . 4.一种大规模MIMO系统信号检测系统,其特征在于,包括:4. A massive MIMO system signal detection system, characterized in that, comprising: 主动禁忌搜索模块,用于对第一初始信号进行主动禁忌搜索检测,以获取初始估计向量;an active tabu search module, which is used to perform active tabu search detection on the first initial signal to obtain an initial estimated vector; 干扰信号消除模块,用于根据所述初始估计向量获取输出比特值,根据所述输出比特值获取信号间总干扰,用所述第一初始信号减去所述信号间总干扰,以获取第二初始信号,其中,根据下式来获取所述输出比特值
Figure FDA0003426212110000021
an interference signal elimination module, configured to obtain an output bit value according to the initial estimation vector, obtain the total interference between signals according to the output bit value, and subtract the total interference between signals from the first initial signal to obtain a second initial signal, where the output bit value is obtained according to
Figure FDA0003426212110000021
Figure FDA0003426212110000022
Figure FDA0003426212110000022
其中,
Figure FDA0003426212110000023
Q是星座图中点的个数,K为MIMO系统在上行链路中用户的个数,
Figure FDA0003426212110000024
为初始估计向量
Figure FDA0003426212110000025
的第i个分量;
in,
Figure FDA0003426212110000023
Q is the number of points in the constellation diagram, K is the number of users in the uplink of the MIMO system,
Figure FDA0003426212110000024
is the initial estimate vector
Figure FDA0003426212110000025
the ith component of ;
根据所述输出比特值获取信号间总干扰的步骤中,所述信号间总干扰
Figure FDA0003426212110000026
的表达式如下:
In the step of obtaining the total inter-signal interference according to the output bit value, the total inter-signal interference
Figure FDA0003426212110000026
The expression is as follows:
Figure FDA0003426212110000027
Figure FDA0003426212110000027
其中,
Figure FDA0003426212110000028
H为信道增益矩阵;
in,
Figure FDA0003426212110000028
H is the channel gain matrix;
向量估计获取模块,用于对所述第二初始信号进行消息传递检测以获取输出向量估计;a vector estimation obtaining module, configured to perform message passing detection on the second initial signal to obtain an output vector estimation; 符号向量重构模块,用于根据所述输出向量估计对符号向量进行重构,以获取符号向量重构值;a symbol vector reconstruction module, configured to reconstruct the symbol vector according to the output vector estimation to obtain a symbol vector reconstruction value; 迭代输出模块,用于将所述符号向量重构值作为主动禁忌搜索的输入进行迭代操作,迭代结束后的最终的所述符号向量重构值作为检测结果输出。The iterative output module is configured to perform an iterative operation using the reconstructed value of the symbol vector as an input of the active tabu search, and output the final reconstructed value of the symbol vector after the iteration is completed as a detection result.
5.一种基站,其特征在于,所述基站包括:5. A base station, wherein the base station comprises: 基站本体,设置于基站本体上的若干接收天线和控制单元;The base station body, a plurality of receiving antennas and control units arranged on the base station body; 其中,所述接收天线用于接收MIMO系统的用户发送的符号向量,并将所述符号向量作为第一初始信号传输给所述控制单元进行信号检测;所述控制单元包括相互耦合的处理器和存储器,所述存储器存储有计算机程序,当所述存储器存储的计算机程序被所述处理器执行时实现权利要求1-3中任意一项所述的大规模MIMO系统信号检测方法。The receiving antenna is used to receive a symbol vector sent by a user of the MIMO system, and transmit the symbol vector as a first initial signal to the control unit for signal detection; the control unit includes a processor coupled to each other and a A memory, where a computer program is stored in the memory, and when the computer program stored in the memory is executed by the processor, the signal detection method for a massive MIMO system according to any one of claims 1-3 is implemented. 6.一种存储介质,其特征在于,包括程序,当所述程序在计算机上运行时,使得所述计算机执行如权利要求1-3中任意一项所述的大规模MIMO系统信号检测方法。6 . A storage medium, characterized by comprising a program, which, when the program runs on a computer, causes the computer to execute the signal detection method for a massive MIMO system according to any one of claims 1 to 3 .
CN202110330579.8A 2021-03-23 2021-03-23 Massive MIMO system signal detection method, system, base station and storage medium Active CN113098571B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110330579.8A CN113098571B (en) 2021-03-23 2021-03-23 Massive MIMO system signal detection method, system, base station and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110330579.8A CN113098571B (en) 2021-03-23 2021-03-23 Massive MIMO system signal detection method, system, base station and storage medium

Publications (2)

Publication Number Publication Date
CN113098571A CN113098571A (en) 2021-07-09
CN113098571B true CN113098571B (en) 2022-03-01

Family

ID=76670232

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110330579.8A Active CN113098571B (en) 2021-03-23 2021-03-23 Massive MIMO system signal detection method, system, base station and storage medium

Country Status (1)

Country Link
CN (1) CN113098571B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113726368A (en) * 2021-09-01 2021-11-30 南京信息工程大学 LTS signal detection method based on constellation constraint interference elimination
CN115001544B (en) * 2022-05-10 2024-06-28 天航长鹰(江苏)科技有限公司 Signal detection method based on improved Richardson algorithm
CN115051900B (en) * 2022-05-10 2024-01-16 四川创智联恒科技有限公司 Wireless MIMO receiver detection method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105406909A (en) * 2015-10-27 2016-03-16 电子科技大学 MIMO (Multiple Input Multiple Output) detection method and system based on internal iterative interference cancellation
CN109951214A (en) * 2019-02-22 2019-06-28 南京邮电大学 A Signal Detection Method Applicable to Massive MIMO System

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8503544B2 (en) * 2010-04-30 2013-08-06 Indian Institute Of Science Techniques for decoding transmitted signals using reactive taboo searches (RTS)
CN104661239B (en) * 2015-03-09 2018-06-19 北京邮电大学 A kind of RI-RTS detection method and device
CN112235216A (en) * 2020-10-16 2021-01-15 哈尔滨工程大学 Compressed sensing pilot frequency design method based on tabu search

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105406909A (en) * 2015-10-27 2016-03-16 电子科技大学 MIMO (Multiple Input Multiple Output) detection method and system based on internal iterative interference cancellation
CN109951214A (en) * 2019-02-22 2019-06-28 南京邮电大学 A Signal Detection Method Applicable to Massive MIMO System

Also Published As

Publication number Publication date
CN113098571A (en) 2021-07-09

Similar Documents

Publication Publication Date Title
Cammerer et al. Trainable communication systems: Concepts and prototype
CN113098571B (en) Massive MIMO system signal detection method, system, base station and storage medium
Li et al. Multiple output selection-LAS algorithm in large MIMO systems
Wang et al. Multiuser detection for uplink large-scale MIMO under one-bit quantization
Lee et al. Deep learning based transmit power control in underlaid device-to-device communication
CN107241167B (en) An Improved Method of BP Detection Based on Massive MIMO System
He et al. Robust BICM design for the LDPC coded DCO-OFDM: A deep learning approach
CN107276935B (en) Method and apparatus for sequential spherical decoding
Qiao et al. Massive access in media modulation based massive machine-type communications
Tian et al. A low complexity detector for downlink SCMA systems
Bian et al. Joint activity detection, channel estimation, and data decoding for grant-free massive random access
Narasimhan et al. Factor graph based joint detection/decoding for LDPC coded large-MIMO systems
Qi et al. When NOMA meets sparse signal processing: Asymptotic performance analysis and optimal sequence design
CN113067667B (en) User activity and multi-user joint detection method
Nguyen et al. QR-decomposition-aided tabu search detection for large MIMO systems
CN109995403A (en) Improved LAS detection algorithm based on simulated annealing idea in massive MIMO system
CN112383328B (en) Improved matched filtering message transmission detection method based on probability cutting in communication system
Azari et al. Automated deep learning-based wide-band receiver
Nguyen et al. Groupwise neighbor examination for tabu search detection in large MIMO systems
CN113691353B (en) Low-complexity signal detection method and system
Narasimhan et al. Channel hardening-exploiting message passing (CHEMP) receiver in large MIMO systems
Chen et al. Memory AMP for generalized MIMO: Coding principle and information-theoretic optimality
CN109067683B (en) Blind detection and modulation constellation optimization method and storage medium in wireless communication
CN108199752A (en) A kind of dual polarization satellite MIMO signal detection method of low complex degree
Eswaramoorthi et al. Deep graph neural network optimized with fertile field algorithm based detection model for uplink multiuser massive multiple‐input and multiple‐output system

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