WO2019134627A1

WO2019134627A1 - Receiving-end signal processing method and apparatus

Info

Publication number: WO2019134627A1
Application number: PCT/CN2018/125559
Authority: WO
Inventors: 刘若鹏; 季春霖; 尤琳
Original assignee: 深圳超级数据链技术有限公司
Priority date: 2018-01-03
Filing date: 2018-12-29
Publication date: 2019-07-11

Abstract

Disclosed are a receiving-end signal processing method and apparatus. The method comprises: performing QR decomposition on a multiplexed waveform matrix H; calculating a signal sequence y according to a result of the QR decomposition and a pre-processed signal sequence r, wherein the pre-processed signal sequence r is obtained by pre-processing a received signal sequence; and calculating, according to the signal sequence y, a sequence X to be transmitted. By means of the present invention, the effect of reducing processing complexity of a receiving-end signal is achieved.

Description

接收端信号处理方法和装置Receiver signal processing method and device

技术领域Technical field

本发明涉及通信领域，具体而言，涉及一种接收端信号处理方法和装置。The present invention relates to the field of communications, and in particular to a method and apparatus for processing signal at a receiving end.

背景技术Background technique

OvXDM(X代表时间T，频率F，码分C，空间S或混合H等)***常见的译码方法有维特比译码算法、MAP、Log-MAP等译码方法，当重叠复用系数K较大时，其计算量和复杂度呈指数率增长，且需要较大的存储量，相应的接收端信号处理复杂度较高。OvXDM (X stands for time T, frequency F, code division C, space S or mixed H, etc.) The common decoding methods are Viterbi decoding algorithm, MAP, Log-MAP and other decoding methods, when the overlap multiplexing coefficient K When it is large, its computational complexity and complexity increase exponentially, and it requires a large amount of storage, and the corresponding signal processing complexity of the receiving end is high.

针对OvXDM***中当重叠复用系数较大时接收端信号处理复杂度较高的问题，目前尚未提出有效的解决方案。Aiming at the problem that the signal processing complexity of the receiving end is high when the overlap multiplexing coefficient is large in the OvXDM system, an effective solution has not been proposed yet.

发明内容Summary of the invention

本发明的一个方面，提供一种接收端信号处理方法，该方法包括：对复用波形矩阵H进行QR分解；根据QR分解的结果和预处理后的信号序列r计算信号序列y，其中，所述预处理后的信号序列r是对接收到的信号序列进行预处理得到的；根据所述信号序列y计算待传输序列X。An aspect of the present invention provides a receiving end signal processing method, the method comprising: performing QR decomposition on a multiplexed waveform matrix H; calculating a signal sequence y according to a result of QR decomposition and a preprocessed signal sequence r, wherein The pre-processed signal sequence r is obtained by pre-processing the received signal sequence; the sequence X to be transmitted is calculated according to the signal sequence y.

本发明的另一个方面，提供一种接收端信号处理装置，包括：分解单元，用于对复用波形矩阵H进行QR分解；第一计算单元，用于根据QR分解的结果和预处理后的信号序列r计算信号序列y，其中，所述预处理后的信号序列r是对接收到的信号序列进行预处理得到的；第二计算单元，用于根据所述信号序列y计算待传输序列X.According to another aspect of the present invention, a receiving end signal processing apparatus includes: a decomposing unit for performing QR decomposition on a multiplexed waveform matrix H; and a first calculating unit for performing a result of QR decomposition and pre-processed The signal sequence r calculates a signal sequence y, wherein the pre-processed signal sequence r is obtained by pre-processing the received signal sequence; the second calculating unit is configured to calculate the sequence X to be transmitted according to the signal sequence y .

实施本发明，通过根据OvXDM***中OvXDM编码特性对数据进行相应译码，避免了现有技术中OvXDM***重叠复用系数较大时接收端信号处理复杂度较高的问题，达到了降低接收端信号处理复杂度的效果。By implementing the invention, the data is correspondingly decoded according to the OvXDM coding characteristic in the OvXDM system, thereby avoiding the problem that the signal processing complexity of the receiving end is high when the overlap multiplexing coefficient of the OvXDM system is large in the prior art, and the receiving end is lowered. The effect of signal processing complexity.

附图说明DRAWINGS

图1a是根据本发明实施例的一种可选的接收端信号处理方法的流程图；FIG. 1a is a flowchart of an optional receiving end signal processing method according to an embodiment of the present invention; FIG.

图1b是根据本发明实施例的一种可选的接收端信号处理方法的流程图；FIG. 1b is a flowchart of an optional receiving end signal processing method according to an embodiment of the present invention; FIG.

图2是根据本发明实施例的OvXDM卷积编码等效模型的示意图；2 is a schematic diagram of an OvXDM convolutional coding equivalent model according to an embodiment of the present invention;

图3是根据本发明实施例的序列正交化示意图；3 is a schematic diagram of sequence orthogonalization according to an embodiment of the present invention;

图4是根据本发明实施例的OvTDM***发送端编码处理框图；4 is a block diagram of a coding process of a transmitting end of an OvTDM system according to an embodiment of the present invention;

图5是根据本发明实施例的OvTDM K路复用波形排列示意图；FIG. 5 is a schematic diagram of an OvTDM K-way multiplexed waveform arrangement according to an embodiment of the present invention; FIG.

图6是根据本发明实施例的OvFDM***发送端编码处理框图；6 is a block diagram of a coding process of a transmitting end of an OvFDM system according to an embodiment of the present invention;

图7是根据本发明实施例的OvFDM K路复用波形排列示意图；7 is a schematic diagram of an OvFDM K-way multiplexed waveform arrangement according to an embodiment of the present invention;

图8是根据本发明实施例的OvTDM接收端处理框图；FIG. 8 is a block diagram of an OvTDM receiving end processing according to an embodiment of the present invention; FIG.

图9是根据本发明实施例的OvFDM接收端处理框图；9 is a block diagram of an OvFDM receiving end processing according to an embodiment of the present invention;

图10a是根据本发明实施例的一种可选的接收端信号处理装置的示意图；FIG. 10a is a schematic diagram of an optional receiving end signal processing apparatus according to an embodiment of the present invention; FIG.

图10b是根据本发明实施例的一种可选的接收端信号处理装置的示意图。Figure 10b is a schematic diagram of an optional receiver signal processing apparatus in accordance with an embodiment of the present invention.

具体实施方式Detailed ways

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。It should be noted that the embodiments in the present application and the features in the embodiments may be combined with each other without conflict. The invention will be described in detail below with reference to the drawings in conjunction with the embodiments.

为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。The technical solutions in the embodiments of the present application are clearly and completely described in the following with reference to the accompanying drawings in the embodiments of the present application. It is an embodiment of the present application, and not all of the embodiments. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present application without departing from the inventive scope shall fall within the scope of the application.

需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本申请的实施例。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。It should be noted that the terms "first", "second" and the like in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or order. It should be understood that the data so used may be interchanged where appropriate to facilitate the embodiments of the present application described herein. In addition, the terms "comprises" and "comprises" and "the" and "the" are intended to cover a non-exclusive inclusion, for example, a process, method, system, product, or device that comprises a series of steps or units is not necessarily limited to Those steps or units may include other steps or units not explicitly listed or inherent to such processes, methods, products or devices.

本发明实施例提供了一种接收端信号处理方法。The embodiment of the invention provides a signal processing method at the receiving end.

图1a是根据本发明实施例的一种可选的接收端信号处理方法的流程图，如图1a所示，该方法包括以下步骤：FIG. 1a is a flowchart of an optional receiving end signal processing method according to an embodiment of the present invention. As shown in FIG. 1a, the method includes the following steps:

步骤S102a：对复用波形矩阵H进行QR分解。Step S102a: QR decomposition is performed on the multiplexed waveform matrix H.

步骤S104a：根据QR分解的结果和预处理后的信号序列r计算信号序列y，其中，预处理后的信号序列r是对接收到的信号序列进行预处理得到的。Step S104a: Calculate the signal sequence y according to the result of the QR decomposition and the pre-processed signal sequence r, wherein the pre-processed signal sequence r is obtained by pre-processing the received signal sequence.

步骤S106a：根据信号序列y计算待传输序列X。Step S106a: Calculate the sequence X to be transmitted according to the signal sequence y.

在本申请实施例中，对复用波形矩阵H进行QR分解，根据QR分解的结果和预处理后的信号序列r计算信号序列y，其中，预处理后的信号序列r是对接收到的信号序列进行预处理得到的，根据信号序列y计算待传输序列X，通过根据OvXDM***中OvXDM编码特性对数据进行相应译码，避免了现有技术中OvXDM***重叠复用系数较大时接收端信号处理复杂度较高的问题，达到了降低接收端信号处理复杂度的效果。In the embodiment of the present application, QR decomposition is performed on the multiplexed waveform matrix H, and the signal sequence y is calculated according to the result of the QR decomposition and the pre-processed signal sequence r, wherein the pre-processed signal sequence r is the received signal The sequence is preprocessed, and the sequence X to be transmitted is calculated according to the signal sequence y, and the data is correspondingly decoded according to the OvXDM coding characteristic in the OvXDM system, thereby avoiding the signal of the receiving end when the overlap multiplexing coefficient of the OvXDM system is large in the prior art. Handling the problem of high complexity has the effect of reducing the complexity of signal processing at the receiving end.

在一个改进的例子中，可以对步骤S102a的复用波形矩阵先进行扩展，然后再进行QR分解；此时对应在步骤S106a中则对扩展后信号序列r _E计算信号序列y。 In a modified example, the multiplexed waveform matrix of step S102a may be first expanded and then QR resolved; at this time, the signal sequence y is calculated for the expanded signal sequence r _E in step S106a.

在一个改进的例子中。可以在步骤S106a之后，计算传输序列X的初步预估值，然后再根据零矩阵计算带传输序列的二次预估值。In a modified example. After step S106a, a preliminary estimate of the transmission sequence X can be calculated, and then a second estimate of the transmitted sequence can be calculated from the zero matrix.

上述的两个改进可以择一进行，也可以通知进行。以下具体通过两个改进同时执行进行说明：The above two improvements can be made alternatively or notified. The following is specifically explained by the simultaneous implementation of two improvements:

图1b是根据本发明另一实施例的一种可选的接收端信号处理方法的流程图，如图1b所示，该方法包括以下步骤：FIG. 1b is a flowchart of an optional receiving end signal processing method according to another embodiment of the present invention. As shown in FIG. 1b, the method includes the following steps:

步骤S102b：对扩展后复用波形矩阵H _E进行QR分解，其中，扩展后复用波形矩阵H _E是对复用波形矩阵H进行扩展得到的； Step S102b: multiplexed after extended QR decomposition matrix H _E waveform, wherein the waveform of the multiplexed extension matrix H _E is multiplexed waveform obtained expanded matrix H;

步骤S104b：根据QR分解的结果和扩展后信号序列r _E计算信号序列y，其中，扩展后信号序列r _E是将预处理后的信号序列r进行扩展得到的，预处理后的信号序列r是对接收到的信号序列进行预处理得到的； Step S104b: Calculate the signal sequence y according to the result of the QR decomposition and the expanded signal sequence r _E , wherein the extended signal sequence r _E is obtained by expanding the pre-processed signal sequence r, and the pre-processed signal sequence r is Pre-processing the received signal sequence;

步骤S106b：根据信号序列y计算待传输序列X的初步预估值；Step S106b: calculating a preliminary estimated value of the sequence X to be transmitted according to the signal sequence y;

步骤S108b：根据置零矩阵计算待传输序列X的二次预估值。Step S108b: Calculate a second estimated value of the sequence X to be transmitted according to the zeroing matrix.

在本申请实施例中，对扩展后复用波形矩阵H _E进行QR分解，根据QR分解的结果和扩展后信号序列r _E计算信号序列y，根据信号序列y计算待传输序列X的初步预估值，根据置零矩阵计算待传输序列X的二次预估值，通过根据OvXDM***中OvXDM编码特性对数据进行相应译码，避免了现有技术中OvXDM***重叠复用系数较大时接收端信号处理复杂度较高的问题，达到了降低接收端信号处理复杂度的效果。 In the embodiment of the present application, QR decomposition is performed on the extended multiplexed waveform matrix H _E , and the signal sequence y is calculated according to the result of the QR decomposition and the extended signal sequence r _E , and a preliminary estimation of the sequence X to be transmitted is calculated according to the signal sequence y. The value is calculated according to the zero-matrix matrix, and the second prediction value of the sequence X to be transmitted is calculated according to the OvXDM coding characteristic in the OvXDM system, thereby avoiding the receiving end of the OvXDM system with large overlap multiplexing coefficient in the prior art. The problem of high signal processing complexity achieves the effect of reducing the complexity of signal processing at the receiving end.

注意：本申请实施例中的待传输序列即为输入信号序列。Note that the sequence to be transmitted in the embodiment of the present application is an input signal sequence.

OvXDM***可表示为重叠时分复用(Overlapped Time Division Multiplexing，OvTDM)***、重叠频分复用(Overlapped Frequency Division Multiplexing，OvFDM)***、重叠码分复用(Overlapped Code Division Multiplexing，OvCDM)***、重叠空分复用(Overlapped Space Division Multiplexing，OvSDM)***、重叠混合复用(Overlapped Hybrid Division Multiplexing，OvHDM)***等，其***等效模型如附图2所示。图3是根据本发明实施例的序列正交化示意图。The OvXDM system can be represented as an Overlapped Time Division Multiplexing (OvTDM) system, an Overlapped Frequency Division Multiplexing (OvFDM) system, an Overlapped Code Division Multiplexing (OvCDM) system, and an overlap. The system equivalent model of the Overlapped Space Division Multiplexing (OvSDM) system and the Overlapped Hybrid Division Multiplexing (OvHDM) system is shown in FIG. 2 . 3 is a schematic diagram of sequence orthogonalization in accordance with an embodiment of the present invention.

以OvTDM***为例，其发送端编码具体处理步骤如下：Take the OvTDM system as an example. The specific processing steps of the sender code are as follows:

(1)首先设计生成发送信号的包络波形h(t)。(1) First, an envelope waveform h(t) that generates a transmission signal is designed.

(2)将(1)中所设计的包络波形h(t)经特定时间移位后，形成其它各个时刻发送信号包络波形h(t-i×ΔT)。(2) After the envelope waveform h(t) designed in (1) is shifted by a specific time, the signal envelope waveform h (t-i × ΔT) is formed at each other time.

(3)将包络波形h(t-i×ΔT)写成矩阵H形式，然后与所要发送的符号向量x相乘，形成发射信号波形。(3) The envelope waveform h(t-i × ΔT) is written in the form of a matrix H, and then multiplied by the symbol vector x to be transmitted to form a transmission signal waveform.

OvTDM***发送端编码处理框图如附图4所示，重叠复用方法遵循平行四边形规则，如附图5所示。The OvTDM system sender encoding process block diagram is shown in Figure 4, and the overlap multiplexing method follows the parallelogram rule, as shown in Figure 5.

以OvFDM***为例，其发送端***编码具体处理步骤如下：Taking the OvFDM system as an example, the specific processing steps of the system coding at the transmitting end are as follows:

(1)首先设计生成发送信号的频谱信号H(f)。(1) First, a spectrum signal H(f) that generates a transmission signal is designed.

(2)将(1)所设计的谱信号H(f)经特定载波频谱间隔ΔB移位后，形成其它各个频谱间隔为ΔB的子载波频谱波形H(f-i×ΔB)。(2) After the spectral signal H(f) designed in (1) is shifted by the specific carrier spectral interval ΔB, the other subcarrier spectral waveforms H (f-i × ΔB) each having a spectral interval of ΔB are formed.

(3)将频谱波形H(f-i×ΔB)写成矩阵H形式，然后与所要发送的符号向量X相乘，形成复调制信号的频谱S(f)。(3) The spectrum waveform H (f-i × ΔB) is written in the form of a matrix H, and then multiplied by the symbol vector X to be transmitted to form a spectrum S(f) of the complex modulated signal.

(4)将(3)生成的复调制信号的频谱进行离散傅氏反变换，最终形成时间域的复调制信号，发送信号可表示为：(4) Performing the discrete Fourier inverse transform on the spectrum of the complex modulated signal generated by (3), and finally forming a complex modulated signal in the time domain, and the transmitted signal can be expressed as:

Signal(t) _TX＝ifft(S(f)) Signal(t) _TX =ifft(S(f))

OvFDM***发送端编码处理框图如附图6所示，重叠复用方法遵循平行四边形规则，如附图7所示。The OvFDM system transmitting end encoding processing block diagram is as shown in FIG. 6, and the overlapping multiplexing method follows the parallelogram rule, as shown in FIG.

OvXDM接收端处理过程：OvXDM receiver processing:

对接收端接收到的信号进行预处理，得到预处理的信号；Pre-processing the signal received by the receiving end to obtain a pre-processed signal;

对所述预处理信号在对应域内按照上述MMSE-QR分解-并行干扰抵消算法进行信号检测，得到输入的信息流。The pre-processed signal is detected in the corresponding domain according to the MMSE-QR decomposition-parallel interference cancellation algorithm to obtain an input information stream.

其中所述预处理过程包括：对接收端接收到的信号进行同步、信道估计、均衡处理等运算。The preprocessing process includes: performing synchronization, channel estimation, equalization processing, and the like on the signal received by the receiving end.

以OvTDM为例，其接收端处理过程如附图8所示，具体步骤如下：Taking OvTDM as an example, the processing of the receiving end is as shown in Figure 8. The specific steps are as follows:

(1)首先对接收信号进行同步，包括载波同步、帧同步、符号时间同步等；(1) First, the received signal is synchronized, including carrier synchronization, frame synchronization, symbol time synchronization, etc.;

(2)按照上述检测算法对预处理后的数据进行相应检测。(2) Perform corresponding detection on the preprocessed data according to the above detection algorithm.

以OvFDM为例，其接收端处理过程如附图9所示，具体步骤如下：Taking OvFDM as an example, the processing of the receiving end is as shown in Figure 9. The specific steps are as follows:

(1)首先对接收信号进行fft(傅里叶变换)运算，使时域信号转换到频域；(1) First performing an Fft (Fourier Transform) operation on the received signal to convert the time domain signal to the frequency domain;

(2)对频域信号进行同步，包括载波同步、帧同步、符号时间同步等；(2) Synchronizing the frequency domain signals, including carrier synchronization, frame synchronization, symbol time synchronization, etc.;

(3)按照上述检测算法对预处理后的数据进行相应检测。(3) Corresponding detection of the preprocessed data according to the above detection algorithm.

可选地，对扩展后复用波形矩阵H _E进行QR分解，根据QR分解的结果和扩展后信号序列r _E计算信号序列y，包括：按照公式H _E＝QR对扩展后复用波形矩阵H _E进行QR分解，其中，Q为(N+L)×L矩阵，Q ^HQ＝I _L，R为L×L的上三角矩阵，矩阵Q分解为N×L的Q ₁和L×L的Q ₂，σI _L＝Q ₂R，

σ ²为噪声功率；根据以下公式计算信号序列y：y＝Q ^Hr _E＝RX-σ ²R ^-HX+Q ₁ ^HN＝RX+η，其中，

为信号序列y中的噪声序列。 Optionally, QR decomposition is performed on the extended multiplexed waveform matrix H _E , and the signal sequence y is calculated according to the result of the QR decomposition and the extended signal sequence r _E , including: expanding the multiplexed waveform matrix H according to the formula H _E =QR pair _E performs QR decomposition, where Q is a (N+L)×L matrix, Q ^H Q=I _L , R is an upper triangular matrix of L×L, and matrix Q is decomposed into N×L Q ₁ and L×L Q ₂ , σI _L = Q ₂ R,

σ ² is the noise power; the signal sequence y is calculated according to the following formula: y = Q ^H r _E = RX - σ ² R - ^H X + Q ₁ ^H N = RX + η, where

Is the noise sequence in the signal sequence y.

首先将复用波形矩阵H扩展为(N+L)×L的H _E矩阵(原复用波形矩阵为N×L矩阵)，预处理后的信号序列r扩展为(N+L)×1的r _E序列，则 First, the multiplexed waveform matrix H is expanded to an (N+L)×L H _E matrix (the original multiplexed waveform matrix is an N×L matrix), and the preprocessed signal sequence r is expanded to (N+L)×1. r _E sequence, then

其中

(σ ²为噪声功率)，I _L为L×L的单位矩阵，0 _L，1为L×1的零向量. among them

(σ ² is the noise power), I _L is the unit matrix of L × L, 0 _{L, 1} is the zero vector of L × 1.

而后对扩展后的H _E矩阵进行QR分解，即： The QR decomposition of the extended H _E matrix is then performed, namely:

其中矩阵Q为(N+L)×L矩阵，并且Q ^HQ＝I _L，R为L×L的上三角矩阵，矩阵Q分解为N×L的Q ₁和L×L的Q ₂，同时满足σI _L＝Q ₂R，进而可得

Where matrix Q is (N+L)×L matrix, and Q ^H Q=I _L , R is an upper triangular matrix of L×L, and matrix Q is decomposed into Q _{1 of} L×L and Q _{2 of} L×L, Satisfy σI _L =Q ₂ R, and then obtain

并且and

对r _E序列进行如下处理，得到处理后的序列y The r _E sequence is processed as follows to obtain the processed sequence y

从上式可看出，该算法相比于QR分解，相比于QR分解，该算法在y的表达式中考虑到噪声的影响，因此，MMSE-QR检测算法能够一定程度上减少噪声的影响，预估结果更精确。It can be seen from the above equation that compared with QR decomposition, the algorithm considers the influence of noise in the expression of y compared to QR decomposition. Therefore, the MMSE-QR detection algorithm can reduce the influence of noise to some extent. The estimated results are more accurate.

可选地，根据信号序列y计算待传输序列X的初步预估值，包括：根据以下公式计算待传输序列X的第k个信号的初步预估值：

其中，y _k＝R _k，k·x _k+η _k+d _k，

R _k，k为矩阵R中第k行、第k列对应的数据，R _k，j为矩阵R中第k行、第j列对应的数据，x _j为待传输序列X中的第j个元素，η _k为信号序列y中的噪声序列η的第k个元素。 Optionally, calculating a preliminary estimated value of the sequence X to be transmitted according to the signal sequence y includes: calculating a preliminary estimated value of the kth signal of the sequence X to be transmitted according to the following formula:

Where y _k =R _k,k ·x _k +η _k +d _k ,

R _k,k is the data corresponding to the kth row and the kth column in the matrix R, R _k,j is the data corresponding to the kth row and the jth column in the matrix R, and x _j is the jth in the sequence X to be transmitted The element, η _{k ,} is the kth element of the noise sequence η in the signal sequence y.

根据处理后序列y进行进一步估计：Further estimation based on the processed sequence y:

处理后的序列的第k个元素y _k可表示为： The kth element y _{k of the} processed sequence can be expressed as:

y _k＝R _k，k·x _k+η _k+d _k y _k =R _k,k ·x _k +η _k +d _k

其中R _k，k为矩阵R中第k行，第k列对应的数据，

R _k，j为矩阵R中第k行，第j列对应的数据，x _j为输入序列X(也即OvXDM编码输入序列)中的第j个元素，d _k独立于上层信号x ₁，x ₂，...，x _k-1，由于为R是上三角矩阵，故可先将最底层信号(第L个信号)解出来： Where R _k,k is the data corresponding to the kth row and the kth column of the matrix R,

R _k,j is the kth row in the matrix R, the data corresponding to the jth column, x _j is the jth element in the input sequence X (that is, the OvXDM coded input sequence), and d _{k is} independent of the upper layer signal x ₁ , x ₂ ,...,x _k-1 , since R is an upper triangular matrix, the lowest signal (the Lth signal) can be solved first:

y _L＝R _L，L·x _L+η _L y _L =R _L,L ·x _L +η _L

对应的估计为：The corresponding estimate is:

对应检测k＝1，2，...，L-1个信号，信号干扰项d _k为： Corresponding to detecting k=1, 2, ..., L-1 signals, the signal interference term d _k is:

对应的估计为：The corresponding estimate is:

可选地，根据置零矩阵计算待传输序列X的二次预估值，包括：根据

计算待传输序列X的二次预估值，其中，G _k＝(H _k ^HH _k) ^-1H _k ^H或G _k＝(H _k ^HH _k+σ ²) ^-1H _k ^H，其中，r _k是将预处理后的信号序列r进行干扰抑制得到的。 Optionally, calculating a second estimated value of the sequence X to be transmitted according to the zeroing matrix, including:

Calculate the second estimate of the sequence X to be transmitted, where G _k = (H _k ^H H _k ) ^-1 H _k ^H or G _k = (H _k ^H H _k + σ ² ) ^-1 H _k ^H , where r _k is obtained by performing interference suppression on the preprocessed signal sequence r.

可选地，在根据置零矩阵计算待传输序列X的二次预估值之前，方法还包括：根据

对预处理后的信号序列r进行干扰抑制，其中，其中，(H) _j表示取H的第j列。 Optionally, before calculating the second estimated value of the sequence X to be transmitted according to the zeroing matrix, the method further includes:

Interference suppression is performed on the preprocessed signal sequence r, where (H) _j represents the jth column of H.

并行干扰抵消算法采用并行处理的方式进行符号间的干扰消除，在输入信号X初始估计值的基础上，恢复各输入信号，在判决信号的过程中，不需进行排序，而是直接进行判决信号。具体做法是：利用检测结果(初始估计值)构造所发送符号的干扰信号估计，在恢复某个输入信号时，都要把其余输入信号的影响作为干扰抵消掉，即在恢复第k个信号时，把第1个，第2个，.....第k-1个，第k+1个，....第L个信号作为干扰抵消掉，得到新的接收向量，然后判决输出第k个信号。该检测算法结合SQRD检测算法，即SQRD-并行干扰抵消算法，其具体检测步骤如下所示：The parallel interference cancellation algorithm uses parallel processing to eliminate interference between symbols. On the basis of the initial estimated value of the input signal X, each input signal is recovered. In the process of determining the signal, no sorting is needed, but the decision signal is directly performed. . The specific method is: constructing the interference signal estimation of the transmitted symbol by using the detection result (initial estimation value), and recovering the influence of the remaining input signals as interference when restoring an input signal, that is, when recovering the kth signal , the first, second, ..... k-1th, k+1th, .... the Lth signal is cancelled as interference, a new reception vector is obtained, and then the output is judged k signals. The detection algorithm combines the SQRD detection algorithm, namely the SQRD-parallel interference cancellation algorithm, and the specific detection steps are as follows:

第一步：根据接收到的信号序列R ₁、复用波形矩阵H，对输入信号X进行初始估计，也即先进行MMSE-QR分解检测估计(如上所述)，得到相应的估值

其中

为输入信号x _k的估计值。接收到的信号经过干扰抑制后的表达式为： The first step is to perform an initial estimation on the input signal X according to the received signal sequence R ₁ and the multiplexed waveform matrix H, that is, first perform MMSE-QR decomposition detection estimation (as described above), and obtain a corresponding estimate.

among them

Is the estimated value of the input signal x _k . The expression of the received signal after interference suppression is:

其中(H) _j表示取H的第j列。从上式可看出在接收到的信号中，将其他所有层的干扰信号都去除了，只留下想要检测的接收信号。 Where (H) _j denotes the jth column of H. It can be seen from the above equation that in the received signal, the interference signals of all other layers are removed, leaving only the received signal that is to be detected.

第二步：计算置零矩阵G _k，可以为迫零检测对应的置零矩阵，也可为最小均方误差检测对应的置零矩阵，即： Step 2: Calculate the zeroing matrix G _k , which can be the zero-setting matrix corresponding to the zero-forcing detection, or the corresponding zero-setting matrix for the minimum mean square error detection, namely:

G _k＝(H _k ^HH _k) ^-1H _k ^H或G _k＝(H _k ^HH _k+σ ²) ^-1H _k ^H G _k =(H _k ^H H _k ) ^-1 H _k ^H or G _k =(H _k ^H H _k +σ ² ) ^-1 H _k ^H

其中H _k表示的是取矩阵H的第k列，σ ²为噪声功率，最后得到检测为： Where H _k represents the kth column of the matrix H, and σ ² is the noise power, and finally the detection is:

将该检测算法替代原OvXDM译码，对应的OvXDM***编码处理过程如下所示：The detection algorithm replaces the original OvXDM decoding, and the corresponding OvXDM system encoding process is as follows:

根据设计参数在调制域内生成包络波形；Generating an envelope waveform in the modulation domain according to design parameters;

根据重叠复用次数将所述包络波形在调制域内按预定的移位间隔进行移位，得到调制域内的各移位包络波形；And shifting the envelope waveform by a predetermined shift interval in the modulation domain according to the number of overlapping multiplexing, to obtain each shift envelope waveform in the modulation domain;

将移位包络波形写成矩阵形式，再与待调制序列中的符号相乘，得到调制域内的复调制包络波形。The shifted envelope waveform is written in a matrix form and multiplied by the symbols in the sequence to be modulated to obtain a complex modulation envelope waveform in the modulation domain.

下面是一个可选的本发明实施例。The following is an optional embodiment of the invention.

根据OvXDM的***特性，首先，假设重叠复用系数为K，复用波形的抽头系数分别定义为[h ₀，h ₁，…，h _K-1]。此时，根据重叠复用关系的卷积特性，若设实信息比特序列长度为L，OvXDM编码后比特序列长度为N，(L+K-1)，则此时，复用波形可用矩阵形式表示为 According to the system characteristics of OvXDM, first, assuming that the overlap multiplexing coefficient is K, the tap coefficients of the multiplexed waveform are defined as [h ₀ , h ₁ , ..., h _K-1 ], respectively. At this time, according to the convolution characteristics of the overlapping multiplexing relationship, if the length of the real information bit sequence is L, and the length of the bit sequence after OvXDM encoding is N, (L+K-1), then the multiplexed waveform can be in the form of a matrix. Expressed as

大小为N×L。

The size is N × L.

设OvXDM编码后的输出向量为Y＝[y ₀，…，y _N-1] ^T，输入向量为X＝[x ₀，…，x _L-1] ^T，OvXDM的编码过程可表示为Y＝HX，即 Let OvXDM encoded output vector be Y=[y ₀ ,...,y _N-1 ] ^T , and the input vector is X=[x ₀ ,...,x _L-1 ] ^T , and the encoding process of OvXDM can be expressed as Y= HX, ie

则此时，接收到的信号序列R ₁可表示为 At this time, the received signal sequence R ₁ can be expressed as

其中，[n ₀，n ₁，…，n _N-1] ^T为白噪声序列。 Where [n ₀ , n ₁ , ..., n _N-1 ] ^T is a white noise sequence.

接收端根据已知的复用波形矩阵H和接收到的信号序列R ₁，进行相应译码。上述接收到的信号序列R ₁与多天线接收序列结构模型相似，均为R ₁＝HX+N，其中X为待传送序列，N为高斯噪声序列，R ₁为接收到的信号序列，不同之处在于矩阵H所代表不同：在多天线***中H表示信道参数矩阵，而在OvXDM***中则表示复用波形矩阵。多天线检测算法包含传统的检测算法，如最小二乘检测算法、最小均方误差检测、最大似然检测、QR分解、双向QR分解、SQRD算法以及串行干扰抵消检测、并行干扰抵消检测等，由于两者结构相似，因而可将检测算法用于对OvXDM***数据进行相应译码。 The receiving end performs corresponding decoding according to the known multiplexed waveform matrix H and the received signal sequence R ₁ . The above received signal sequence R _{1 is} similar to the multi-antenna received sequence structure model, both of which are R ₁ =HX+N, where X is the sequence to be transmitted, N is the Gaussian noise sequence, and R ₁ is the received signal sequence, which is different. The difference is represented by the matrix H: H represents the channel parameter matrix in the multi-antenna system, and the multiplexed waveform matrix is represented in the OvXDM system. The multi-antenna detection algorithm includes traditional detection algorithms such as least squares detection algorithm, minimum mean square error detection, maximum likelihood detection, QR decomposition, bidirectional QR decomposition, SQRD algorithm, serial interference cancellation detection, parallel interference cancellation detection, etc. Because the two structures are similar, the detection algorithm can be used to decode the OvXDM system data accordingly.

可选地，对复用波形矩阵H进行QR分解，包括：分别初始化矩阵Q、矩阵R、和矩阵P，其中，矩阵P是交换矩阵；对初始化后的矩阵Q、矩阵R、和矩阵P进行循环计算，得到最终的矩阵Q、矩阵R、和矩阵P。Optionally, performing QR decomposition on the multiplexed waveform matrix H includes: initializing a matrix Q, a matrix R, and a matrix P, wherein the matrix P is a switching matrix; and performing the initialized matrix Q, the matrix R, and the matrix P Cycle calculations yield the final matrix Q, matrix R, and matrix P.

可选地，对矩阵Q、矩阵R、矩阵P进行初始化，包括：令R＝0，Q＝H，P＝I _L，复用波形矩阵H为N×L的矩阵，矩阵Q为N×L的酉矩阵，且满足Q ^HQ＝I _L，矩阵I _L为L×L的单位矩阵，矩阵R为L×L的上三角矩阵；对初始化后的矩阵Q、矩阵R、矩阵P进行循环计算，得到最终的矩阵Q、矩阵R、矩阵P，包括：在矩阵Q满足第一预设条件的列中查找最小2-范数对应的列；根据查找到的列对矩阵Q中的元素进行归一化和正交化；改变第一预设条件，直至矩阵Q的所有列均完成循环。 Optionally, the matrix Q, the matrix R, and the matrix P are initialized, including: let R=0, Q=H, P=I _L , the multiplexed waveform matrix H is a matrix of N×L, and the matrix Q is N×L酉 matrix, and satisfy Q ^H Q=I _L , matrix I _L is L × L unit matrix, matrix R is L × L upper triangular matrix; cyclic calculation of initialized matrix Q, matrix R, matrix P Obtaining a final matrix Q, a matrix R, and a matrix P, comprising: finding a column corresponding to the minimum 2-norm in a column in which the matrix Q satisfies the first preset condition; and returning the elements in the matrix Q according to the found column Normalization and orthogonalization; changing the first preset condition until all the columns of the matrix Q complete the loop.

可选地，在矩阵Q满足第一预设条件的列中查找最小2-范数对应的列；根据查找到的列对矩阵Q中的元素进行归一化和正交化；改变第一预设条件，直至矩阵Q的所有列均完成循环，包括：步骤1)令i＝1；步骤2)计算Q矩阵中列数l≥i的2-范数，找出其中最小2-范数对应的列数，用k _i表示，即

步骤3)分别交换Q、R、P 矩阵的第k _i列与i列；并且令矩阵R中的主对角线元素r _i，i＝||q _i||；步骤4)根据

对矩阵Q对应的第i列q _i进行归一化；步骤5)根据q _l＝q _l-r _i，lq _i对q _l进行正交化，r _i，l为q _l在q _i方向上的投影系数，其中，i＜l≤L；步骤6)将i+1的值赋值给i；循环执行步骤2)至步骤6)，直至i＝L，得到最终的矩阵Q、矩阵R。 Optionally, searching for a column corresponding to the minimum 2-norm in the column in which the matrix Q satisfies the first preset condition; normalizing and orthogonalizing the elements in the matrix Q according to the found column; changing the first pre- Set the condition until all the columns of the matrix Q complete the loop, including: step 1) let i=1; step 2) calculate the 2-norm of the number of columns l≥i in the Q matrix, find out the minimum 2-norm corresponding The number of columns, denoted by k _i , ie

Step 3) are exchanged Q, R, P K _i of the matrix column and row i; r and causes the main diagonal elements of the matrix R _{_{i, i = || q i ||}} ; Step 4)

Q matrix corresponding to column i q _i is normalized; Step _{_{5) = q l -r i,}} l q i q _l of orthogonalizing according q _{_l,} r _i, _l is the direction q _i q _l The projection coefficient on the above, where i < l ≤ L; step 6) assigns the value of i+1 to i; loops through steps 2) to 6) until i = L, resulting in the final matrix Q, matrix R.

可选地，根据信号序列y计算待传输序列X，包括：根据

计算待传输序列X，其中，y _k＝r _k，k·x _k+η _k+d _k，y _k是信号序列y的第k个元素，

Q(·)是量化判断算子。 Optionally, calculating the sequence X to be transmitted according to the signal sequence y includes: according to

Calculating the sequence X to be transmitted, where y _k =r _k,k ·x _k +η _k +d _k , y _k is the kth element of the signal sequence y,

Q(·) is a quantization judgment operator.

并行干扰抵消检可与传统的检测算法，如迫零检测、最小均方误差检测、QR分解、双向QR分解、SQRD算法、MMSE-QR分解等相结合实现检测过程。本专利主要介绍将MMSE-QR分解-并行干扰检测算法用于OvXDM***数据检测中，其余在此不赘述。Parallel interference cancellation detection can be combined with traditional detection algorithms such as zero-forcing detection, minimum mean square error detection, QR decomposition, two-way QR decomposition, SQRD algorithm, MMSE-QR decomposition, etc. This patent mainly introduces the MMSE-QR decomposition-parallel interference detection algorithm for data detection in the OvXDM system, and the rest is not described here.

首先介绍MMSE-QR分解算法：First introduce the MMSE-QR decomposition algorithm:

(1)首先将复用波形矩阵H扩展为(N+L)×L的H _E矩阵(原复用波形矩阵为N×L矩阵)，将预处理后的信号序列r扩展为(N+L)×1的r _E序列，则 (1) First, the multiplexed waveform matrix H is expanded to an (N+L)×L H _E matrix (the original multiplexed waveform matrix is an N×L matrix), and the preprocessed signal sequence r is expanded to (N+L). ) × 1 r _E sequence, then

其中

(2)而后对扩展后的H _E矩阵进行QR分解，即： (2) Then perform QR decomposition on the extended H _E matrix, namely:

并且and

(3)对r _E序列进行如下处理，得到处理后的序列y (3) Perform the following processing on the r _E sequence to obtain the processed sequence y

从上式可看出，该算法相比于QR分解，MMSE-QR检测算法能够一定程度上减少噪声的影响.It can be seen from the above equation that compared with QR decomposition, the MMSE-QR detection algorithm can reduce the influence of noise to some extent.

(5)根据处理后序列y进行进一步估计：(5) Further estimation based on the processed sequence y:

y _k＝R _k，k·x _k+η _k+d _k y _k =R _k,k ·x _k +η _k +d _k

其中R _k，k为矩阵R中第k行，第k列对应的数据，

R _k，j为矩阵R中第k行，第j列对应的数据，x _j为输入序列X(也即OvXDM编码输入序列)中的第j个元素，d _k独立于上层信号x ₁，x ₂，...，x _k-1，由于为R上三角矩阵，故可先将最底层信号(第L个信号)解出来： Where R _k,k is the data corresponding to the kth row and the kth column of the matrix R,

R _k,j is the kth row in the matrix R, the data corresponding to the jth column, x _j is the jth element in the input sequence X (that is, the OvXDM coded input sequence), and d _{k is} independent of the upper layer signal x ₁ , x ₂ ,...,x _k-1 , since it is a triangular matrix on R, the lowest signal (the Lth signal) can be solved first:

y _L＝R _L，L·x _L+η _L y _L =R _L,L ·x _L +η _L

对应的估计为：The corresponding estimate is:

而后根据发送端采用的调制方式，进行相应解调，而后判决输出。Then, according to the modulation mode adopted by the transmitting end, corresponding demodulation is performed, and then the output is judged.

而并行干扰抵消算法采用并行处理的方式进行符号间的干扰消除，在输入信号X初始估计值的基础上，恢复各输入信号，在判决信号的过程中，不需进行排序，而是直接进行判决信号。具体做法是：利用检测结果(初始估计值)构造所发送符号的干扰信号估计，在恢复某个输入信号时，都要把其余输入信号的影响作为干扰抵消掉，即在恢复第k个信号时，把第1个，第2个，.....第k-1个，第k+1个，....第L个信号作为干扰抵消掉，得到新的接收向量，然后判决输出第k个信号。该检测算法结合SQRD检测算法，即SQRD-并行干扰抵消算法，其具体检测步骤如下所示：The parallel interference cancellation algorithm uses parallel processing to eliminate interference between symbols. On the basis of the initial estimated value of the input signal X, each input signal is recovered. In the process of determining the signal, no sorting is needed, but the decision is directly made. signal. The specific method is: constructing the interference signal estimation of the transmitted symbol by using the detection result (initial estimation value), and recovering the influence of the remaining input signals as interference when restoring an input signal, that is, when recovering the kth signal , the first, second, ..... k-1th, k+1th, .... the Lth signal is cancelled as interference, a new reception vector is obtained, and then the output is judged k signals. The detection algorithm combines the SQRD detection algorithm, namely the SQRD-parallel interference cancellation algorithm, and the specific detection steps are as follows:

其中，

among them,

在另外一些实施例中，将SQRD算法用于OvXDM***数据检测中。具体如下所示：In still other embodiments, the SQRD algorithm is used in OvXDM system data detection. The details are as follows:

SQRD算法首先计算出复用波形矩阵H对应的QR分解中的矩阵Q、R，同时引入一个交换矩阵P，使得矩阵R的上层具有较小的SNR，也即R的对角元素按照从小到大的顺序排列(注意：进行QR分解第k步的SNR _k≈|r _k，k| ²)，而在计算复用波形矩阵H 的QR分解中的矩阵R时，是从上到下一行一行进行的，而计算Q是从左到右一列一列进行的；而后根据矩阵Q、R、P进行检测。假定复用波形矩阵H为N×L的矩阵，则Q为N×L的酉矩阵，且满足Q ^HQ＝I _L，I _L为L×L的单位矩阵，R为L×L的上三角矩阵，具体如下所示： The SQRD algorithm first calculates the matrices Q and R in the QR decomposition corresponding to the multiplexed waveform matrix H, and introduces a switching matrix P, so that the upper layer of the matrix R has a smaller SNR, that is, the diagonal elements of R are from small to large. The order of the order (note: SNR _k ≈ | r _{k, k} | ^{2 of the} _k- th decomposition of the QR decomposition), and when calculating the matrix R in the QR decomposition of the multiplexed waveform matrix H, is performed from top to bottom And the calculation Q is performed from left to right in a column and column; then it is detected according to the matrix Q, R, P. Assuming that the multiplexed waveform matrix H is a matrix of N×L, Q is an N×L 酉 matrix, and satisfies Q ^H Q=I _L , I _L is an L × L unit matrix, and R is an L×L upper triangle. The matrix is as follows:

(1)首先进行初始化：R＝0，Q＝H，P＝I _L (1) Initialize first: R=0, Q=H, P=I _L

(2)然后进行循环计算相应的Q、R、P矩阵：令i＝1 1)计算Q矩阵中列数l≥i的2-范数，找出其中最小2-范数(序列模值)对应的列数，用k _i表示，即

2)而后分别交换Q、R、P矩阵的第k _i列与i列；并且令矩阵R中的主对角线元素r _i，i＝||q _i||；3)接着对矩阵Q对应的第i列q _i进行归一化处理，即

4)对q _l(i＜l≤L)进行正交化，r _i，l为q _l在q _i方向上的投影系数，r _i，l＝q _i ^Hq _l，进而与q _i正交的序列q _l＝q _l-r _i，lq _i(对应图中的e)，如附图3所示； (2) Then calculate the corresponding Q, R, P matrix by loop: let i=1 1) Calculate the 2-norm of the number of columns l≥i in the Q matrix, and find the minimum 2-norm (sequence modulus) The corresponding number of columns, denoted by k _i , ie

2) are then exchanged Q, R, P K _i of the matrix column and row i; r and causes the main diagonal elements of the matrix R _{_{i, i = || q i ||}} ; 3) Next, the corresponding matrix Q The i-th column q _i is normalized, ie

4) Orthogonalize q _l (i<l≤L), r _i,l is the projection coefficient of q _l in the q _i direction, r _i,l =q _i ^H q _l , and then orthogonal to q _i The sequence q _l =q _l -r _i,l q _i (corresponding to e in the figure), as shown in Figure 3;

由于根据误差e与q _i垂直，则 Since it is perpendicular to q _i according to the error e, then

q _i ^H·e＝q _i ^H(q _l-p) q _i ^H ·e=q _i ^H (q _l -p)

＝q _i ^H(q _l-r _i，lq _i) =q _i ^H (q _l -r _i,l q _i )

＝0=0

可得

(其中q _i ^Hq _i＝1)，进而可知q _l在q _i方向上的投影p＝r _i，lq _i，与q _i垂直的e为e＝q _l-p＝q _l-r _i，lq _i，同样也为q _l正交化后对应的序列，表示为q _l＝q _l-r _i，lq _i。 Available

(where q _i ^H q _i =1), and further, the projection of q _l in the q _i direction p=r _i,l q _i , and the e perpendicular to q _i is e=q _l -p=q _l -r _{i , l} q _i , which is also the sequence corresponding to the orthogonalization of q _l , is expressed as q _l =q _l -r _i,l q _i .

5)i＝i+1，再进行(2)中1)-4)的循环，直至矩阵Q的所有列均进行完此循环(i＝L)，最终得到相应的Q、R矩阵。5) i=i+1, and then perform the loop of 1)-4) in (2) until all the columns of matrix Q have completed this cycle (i=L), and finally obtain the corresponding Q and R matrix.

(3)而后进行信号检测：令y＝Q ^Hr，其中r为经过预处理后的接收序列(即上述预处理后的信号序列)，此时y可写为如下形式： (3) Then perform signal detection: Let y=Q ^H r, where r is the pre-processed reception sequence (ie, the pre-processed signal sequence), and y can be written as follows:

其中η _i，i＝1，2，...，L代表噪声部分。

Where η _i , i = 1, 2, ..., L represents the noise portion.

接收矢量第k个元素y _k为： The kth element y _k of the receiving vector is:

y _k＝r _k，k·x _k+η _k+d _k y _k =r _k,k ·x _k +η _k +d _k

于0，则对应的输出为0；而输入数据大于0，则对应的输出为1)。

At 0, the corresponding output is 0; and if the input data is greater than 0, the corresponding output is 1).

最终检测输出序列

其中

Final detection output sequence

among them

本申请实施例还提供了一种接收端信号处理装置，该装置用于执行上述接收端信号处理方法。The embodiment of the present application further provides a receiving end signal processing apparatus, which is used to execute the above-mentioned receiving end signal processing method.

图10a是根据本发明实施例的接收端信号处理装置的示意图，如图10a所示，该装置包括：分解单元10a、第一计算单元20a、第二计算单元30a。FIG. 10a is a schematic diagram of a signal processing apparatus at a receiving end according to an embodiment of the present invention. As shown in FIG. 10a, the apparatus includes a decomposition unit 10a, a first calculation unit 20a, and a second calculation unit 30a.

分解单元10a，用于对复用波形矩阵H进行QR分解。The decomposition unit 10a is configured to perform QR decomposition on the multiplexed waveform matrix H.

第一计算单元20a，用于根据QR分解的结果和预处理后的信号序列r计算信号序列y，其中，预处理后的信号序列r是对接收到的信号序列进行预处理得到的。The first calculating unit 20a is configured to calculate the signal sequence y according to the result of the QR decomposition and the pre-processed signal sequence r, wherein the pre-processed signal sequence r is obtained by pre-processing the received signal sequence.

第二计算单元30a，用于根据信号序列y计算待传输序列X。The second calculation unit 30a is configured to calculate the sequence X to be transmitted according to the signal sequence y.

可选地，分解单元10a包括：初始化子单元、计算子单元。初始化子单元，用于分别初始化矩阵Q、矩阵R、和矩阵P，其中，矩阵P是交换矩阵。计算子单元，用于对初始化后的矩阵Q、矩阵R、和矩阵P进行循环计算，得到最终的矩阵Q、矩阵R、和矩阵P，其中，矩阵Q为酉矩阵，矩阵R为上三角矩阵。Optionally, the decomposition unit 10a includes: an initialization subunit, and a calculation subunit. The initialization subunit is configured to initialize the matrix Q, the matrix R, and the matrix P, respectively, wherein the matrix P is a switching matrix. The calculation subunit is configured to perform cyclic calculation on the initialized matrix Q, the matrix R, and the matrix P to obtain a final matrix Q, a matrix R, and a matrix P, wherein the matrix Q is a unitary matrix and the matrix R is an upper triangular matrix. .

可选地，初始化子单元用于：令R＝0，Q＝H，P＝I _L，复用波形矩阵H为N×L的矩阵，矩阵Q为N×L的矩阵，且满足Q ^HQ＝I _L，矩阵I _L为L×L的单位矩阵，矩阵R为L×L的矩阵。计算子单元用于：在矩阵Q满足第一预设条件的列中查找最小2-范数对应的列；根据查找到的列对矩阵Q中的元素进行归一化和正交化；改变第一预设条件，直至矩阵Q的所有列均完成循环。 Optionally, the initialization subunit is configured to: let R=0, Q=H, P=I _L , the multiplexed waveform matrix H be a matrix of N×L, the matrix Q is a matrix of N×L, and satisfy Q ^H Q =I _L , the matrix I _L is an L × L unit matrix, and the matrix R is an L × L matrix. The calculating subunit is configured to: find a column corresponding to the minimum 2-norm in the column in which the matrix Q satisfies the first preset condition; normalize and orthogonalize the elements in the matrix Q according to the found column; A predetermined condition until all the columns of the matrix Q complete the loop.

可选地，计算子单元用于执行以下步骤：Optionally, the calculation subunit is used to perform the following steps:

步骤1)令i＝1；Step 1) Let i=1;

步骤2)计算Q矩阵中列数l≥i的2-范数，找出其中最小2-范数对应的列数，用k _i表示，即

Step 2) Calculate the 2-norm of the number of columns l ≥ _{i in the} Q matrix, and find the number of columns corresponding to the minimum 2-norm, which is represented by k _i , ie

步骤3)分别交换Q、R、P矩阵的第k _i列与i列；并且令矩阵R中的主对角线元素r _i，i＝||q _i||； K i of the _i-th column and the steps 3) are exchanged Q, R, P matrix; r and causes the main diagonal elements of the matrix R _{_{i, i = || q i ||}} ;

步骤4)根据

对矩阵Q对应的第i列q _i进行归一化； Step 4) According to

Normalizing the i-th column q _i corresponding to the matrix Q;

步骤5)根据q _l＝q _l-r _i，lq _i对q _l进行正交化，r _i，l为q _l在q _i方向上的投影系数，其中，i＜l≤L； Step _{_{5) = q l -r i,}} l q i q _l to be _{orthogonalized,} r _{i, l} is the projection coefficient q _l q _i in the direction according q _l, wherein, i <l≤L;

步骤6)将i+1的值赋值给i；Step 6) assigning the value of i+1 to i;

循环执行步骤2)至步骤6)，直至i＝L，得到最终的矩阵Q、矩阵R。Step 2) to step 6) are performed until i=L, and the final matrix Q and matrix R are obtained.

可选地，第二计算单元30用于：根据

Q(·)是量化判断算子。 Optionally, the second calculating unit 30 is configured to:

Q(·) is a quantization judgment operator.

在一个可实施的改进中，分解单元10还可以对扩展后的复用波形矩阵进行QR分解，此时对应的计算单元对分解后的结果进行处理；In an implementable improvement, the decomposing unit 10 may also perform QR decomposition on the expanded multiplexed waveform matrix, and the corresponding computing unit processes the decomposed result;

在一个可实施的改进中，第二计算单元30待传输序列X的初步预估值；并增加一个第三计算单元，用于根据置零矩阵计算待传输序列X的二次预估值。In an implementable improvement, the second calculation unit 30 is to transmit a preliminary estimate of the sequence X; and a third calculation unit is added for calculating a second estimate of the sequence X to be transmitted according to the zeroing matrix.

如图10b所示，该装置包括：分解单元10b、第一计算单元20b、第二计算单元30b、第三计算单元40b。As shown in FIG. 10b, the apparatus includes a decomposition unit 10b, a first calculation unit 20b, a second calculation unit 30b, and a third calculation unit 40b.

分解单元10b，用于对扩展后复用波形矩阵H _E进行QR分解，其中，扩展后复用波形矩阵H _E是对复用波形矩阵H进行扩展得到的。 The decomposing unit 10b is configured to perform QR decomposition on the expanded multiplexed waveform matrix H _E , wherein the expanded multiplexed waveform matrix H _E is obtained by expanding the multiplexed waveform matrix H.

第一计算单元20b，用于根据QR分解的结果和扩展后信号序列r _E计算信号序列y，其中，扩展后信号序列r _E是将预处理后的信号序列r进行扩展得到的，预处理后的信号序列r是对接收到的信号序列进行预处理得到的。 The first calculating unit 20b is configured to calculate the signal sequence y according to the result of the QR decomposition and the extended signal sequence r _E , wherein the extended signal sequence r _E is obtained by expanding the pre-processed signal sequence r, after pre-processing The signal sequence r is obtained by preprocessing the received signal sequence.

第二计算单元30b，用于根据信号序列y计算待传输序列X的初步预估值。The second calculating unit 30b is configured to calculate a preliminary estimated value of the sequence X to be transmitted according to the signal sequence y.

第三计算单元40b，用于根据置零矩阵计算待传输序列X的二次预估值。The third calculating unit 40b is configured to calculate a second estimated value of the sequence X to be transmitted according to the zeroing matrix.

可选地，QR分解单元按照公式H _E＝QR对扩展后复用波形矩阵H _E进行QR分解，其中，Q为(N+L)×L矩阵，Q ^HQ＝I _L，R为L×L的上三角矩阵，矩阵Q分解为N×L的Q ₁和L×L的Q ₂，σI _L＝Q ₂R，

σ ²为噪声功率；第一计算单元20根据以下公式计算信号序列y：y＝Q ^Hr _E＝RX-σ ²R ^-HX+Q ₁ ^HN＝RX+η，其中，

为信号序列y中的噪声序列。 Optionally, the QR decomposition unit performs QR decomposition on the expanded multiplexed waveform matrix H _E according to the formula H _E =QR, where Q is a (N+L)×L matrix, Q ^H Q=I _L , and R is L× L upper triangular matrix, the matrix Q is decomposed into N × L × L L Q ₁ and the _{_{_{Q 2, σI L = Q 2}}} R,

σ ² is the noise power; the first calculating unit 20 calculates the signal sequence y according to the following formula: y = Q ^H r _E = RX - σ ² R - ^H X + Q ₁ ^H N = RX + η, where

Is the noise sequence in the signal sequence y.

可选地，第二计算单元30根据以下公式计算待传输序列X的第k个信号的初步预估值：

其中，y _k＝R _k，k·x _k+η _k+d _k，

R _k，k为矩阵R中第k行、第k列对应的数据，R _k，j为矩阵R中第k行、第j列对应的数据，x _j为待传输序列X中的第j个元素，η _k为信号序列y中的噪声序列η的第k个元素。 Optionally, the second calculating unit 30 calculates a preliminary estimated value of the kth signal of the sequence X to be transmitted according to the following formula:

Where y _k =R _k,k ·x _k +η _k +d _k ,

可选地，第三计算单元40根据

计算待传输序列X的二次预估值，其中，G _k＝(H _k ^HH _k) ^-1H _k ^H或G _k＝(H _k ^HH _k+σ ²) ^-1H _k ^H，其中，r _k是将预处理后的信号序列r进行干扰抑制得到的。 Optionally, the third calculating unit 40 is configured according to

可选地，装置还包括：干扰抑制单元。干扰抑制单元，用于在第三计算单元40根据置零矩阵计算待传输序列X的二次预估值之前，根据

对预处理后的信号序列r进行干扰抑制，其中，其中，(H) _j表示取H的第j列。 Optionally, the device further comprises: an interference suppression unit. An interference suppression unit, configured to: before the third calculating unit 40 calculates a second estimated value of the sequence X to be transmitted according to the zeroing matrix, according to

所述接收端信号处理装置包括处理器和存储器，上述分解单元10、第一计算单元20、第二计算单元30、第三计算单元40等均作为程序单元存储在存储器中，由处理器执行存储在存储器中的上述程序单元来实现相应的功能。The receiving end signal processing device includes a processor and a memory, and the above-described decomposition unit 10, first calculation unit 20, second calculation unit 30, third calculation unit 40, and the like are all stored as a program unit in a memory, and are executed by the processor. The above described program elements in the memory implement the corresponding functions.

处理器中包含内核，由内核去存储器中调取相应的程序单元。内核可以设置一个或以上，通过调整内核参数来执行接收端信号处理方法。The processor contains a kernel, and the kernel removes the corresponding program unit from the memory. The kernel can be set to one or more, and the receiver signal processing method is performed by adjusting the kernel parameters.

存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)，存储器包括至少一个存储芯片。The memory may include non-persistent memory, random access memory (RAM), and/or non-volatile memory in a computer readable medium, such as read only memory (ROM) or flash memory (flash RAM), the memory including at least one Memory chip.

本发明实施例提供了一种存储介质，其上存储有程序，该程序被处理器执行时实现所述接收端信号处理方法。Embodiments of the present invention provide a storage medium on which a program is stored, and when the program is executed by a processor, the signal processing method of the receiving end is implemented.

本发明实施例提供了一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行所述接收端信号处理方法。An embodiment of the present invention provides a processor, where the processor is configured to run a program, where the program is executed to execute the signal processing method at the receiving end.

本发明实施例提供了一种设备，设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序，处理器执行程序时实现以下步骤：Embodiments of the present invention provide a device including a processor, a memory, and a program stored on the memory and operable on the processor, and the processor implements the following steps when executing the program:

对复用波形矩阵H进行QR分解；根据QR分解的结果和预处理后的信号序列r计算信号序列y，其中，预处理后的信号序列r是对接收到的信号序列进行预处理得到的；根据信号序列y计算待传输序列X。Performing QR decomposition on the multiplexed waveform matrix H; calculating the signal sequence y according to the result of the QR decomposition and the pre-processed signal sequence r, wherein the pre-processed signal sequence r is obtained by pre-processing the received signal sequence; The sequence X to be transmitted is calculated from the signal sequence y.

处理器执行程序时还可以实现以下步骤：分别初始化矩阵Q、矩阵R、和矩阵P，其中，矩阵P是交换矩阵；对初始化后的矩阵Q、矩阵R、和矩阵P进行循环计算，得到最终的矩阵Q、矩阵R、和矩阵P，其中，矩阵Q为酉矩阵，矩阵R为上三角矩阵。The processor may further implement the following steps: initializing the matrix Q, the matrix R, and the matrix P, wherein the matrix P is a switching matrix; performing cyclic calculation on the initialized matrix Q, the matrix R, and the matrix P to obtain a final The matrix Q, the matrix R, and the matrix P, wherein the matrix Q is a unitary matrix and the matrix R is an upper triangular matrix.

处理器执行程序时还可以实现以下步骤：令R＝0，Q＝H，P＝I _L，复用波形矩阵H为N×L的矩阵，矩阵Q为N×L的矩阵，且满足Q ^HQ＝I _L，矩阵I _L为L×L的单位矩阵，矩阵R为L×L的矩阵；在矩阵Q满足第一预设条件的列中查找最小2-范数对应的列；根据查找到的列对矩阵Q中的元素进行归一化和正交化；改变第一预设条件，直至矩阵Q的所有列均完成循环。 When the processor executes the program, the following steps can be implemented: let R=0, Q=H, P=I _L , the multiplexed waveform matrix H be a matrix of N×L, the matrix Q is a matrix of N×L, and satisfy Q ^H Q=I _L , the matrix I _L is an L × L unit matrix, the matrix R is a L×L matrix; the minimum 2-norm corresponding column is found in the column in which the matrix Q satisfies the first preset condition; The columns normalize and orthogonalize the elements in matrix Q; change the first preset condition until all columns of matrix Q complete the loop.

处理器执行程序时还可以实现以下步骤：The following steps can also be implemented when the processor executes the program:

步骤1)令i＝1；Step 1) Let i=1;

步骤4)根据

对矩阵Q对应的第i列q _i进行归一化； Step 4) According to

Normalizing the i-th column q _i corresponding to the matrix Q;

步骤6)将i+1的值赋值给i；Step 6) assigning the value of i+1 to i;

处理器执行程序时还可以实现以下步骤：根据

Q(·)是量化判断算子。 The processor can also implement the following steps when executing the program:

Q(·) is a quantization judgment operator.

此外，处理器执行程序时，还可以执行以下针对扩展后的复用波形矩阵的处理方案：In addition, when the processor executes the program, the following processing scheme for the expanded multiplexed waveform matrix can also be performed:

对扩展后复用波形矩阵H _E进行QR分解，其中，扩展后复用波形矩阵H _E是对复用波形矩阵H进行扩展得到的；根据QR分解的结果和扩展后信号序列r _E计算信号序列y，其中，扩展后信号序列r _E是将预处理后的信号序列r进行扩展得到的，预处理后的信号序列r是对接收到的信号序列进行预处理得到的；根据信号序列y计算待传输序列X的初步预估值；根据置零矩阵计算待传输序列X的二次预估值。 QR decomposition is performed on the expanded multiplexed waveform matrix H _E , wherein the expanded multiplexed waveform matrix H _E is obtained by expanding the multiplexed waveform matrix H; the signal sequence is calculated according to the result of the QR decomposition and the expanded signal sequence r _E y, wherein the extended signal sequence r _E is obtained by expanding the pre-processed signal sequence r, and the pre-processed signal sequence r is obtained by pre-processing the received signal sequence; and calculating according to the signal sequence y The preliminary estimated value of the transmission sequence X; the second estimated value of the sequence X to be transmitted is calculated according to the zeroing matrix.

处理器执行程序时还可以实现以下步骤：按照公式H _E＝QR对扩展后复用波形矩阵H _E进行QR分解，其中，Q为(N+L)×L矩阵，Q ^HQ＝I _L，R为L×L的上三角矩阵，矩阵Q分解为N×L的Q ₁和L×L的Q ₂，σI _L＝Q ₂R，

为信号序列y中的噪声序列。 The processor may also implement the following steps: performing QR decomposition on the expanded multiplexed waveform matrix H _E according to the formula H _E =QR, where Q is a (N+L)×L matrix, Q ^H Q=I _L , R is L × L upper triangular matrix, the matrix Q is decomposed into N × L of Q ₁ and L × L is _{_{_{Q 2, σI L = Q 2}}} R,

Is the noise sequence in the signal sequence y.

处理器执行程序时还可以实现以下步骤：根据以下公式计算待传输序列X的第k个信号的初步预估值：

其中，y _k＝R _k，k·x _k+η _k+d _k，

R _k，k为矩阵R中第k行、第k列对应的数据，R _k，j为矩阵R中第k行、第j列对应的数据，x _j为待传输序列X中的第j个元素，η _k为信号序列y中的噪声序列η的第k个元素。 The processor may also implement the following steps: calculating a preliminary estimate of the kth signal of the sequence X to be transmitted according to the following formula:

Where y _k =R _k,k ·x _k +η _k +d _k ,

处理器执行程序时还可以实现以下步骤：根据

计算待传输序列X的二次预估值，其中，G _k＝(H _k ^HH _k) ^-1H _k ^H或G _k＝(H _k ^HH _k+σ ²) ^-1H _k ^H，其中，r _k是将预处理后的信号序列r进行干扰抑制得到的。 The processor can also implement the following steps when executing the program:

处理器执行程序时还可以实现以下步骤：在根据置零矩阵计算待传输序列X的二次预估值之前，根据

对预处理后的信号序列r进行干扰抑制，其中，其中，(H) _j表示取H的第j列。 The processor may also implement the following steps: before calculating the second estimated value of the sequence X to be transmitted according to the zeroing matrix, according to

本文中的设备可以是服务器、PC、PAD、手机等。The devices in this document can be servers, PCs, PADs, mobile phones, and the like.

本申请还提供了一种计算机程序产品，当在数据处理设备上执行时，适于执行初始化有如下方法步骤的程序：The present application also provides a computer program product, when executed on a data processing device, adapted to perform a process of initializing the method steps as follows:

上述计算机程序产品还可以执行初始化有如下方法步骤的程序：按照公式H _E＝QR对扩展后复用波形矩阵H _E进行QR分解，其中，Q为(N+L)×L矩阵，Q ^HQ＝I _L，R为L×L的上三角矩阵，矩阵Q分解为N×L的Q ₁和L×L的Q ₂，σI _L＝Q ₂R，

为信号序列y中的噪声序列。 The above computer program product may also execute a program for initializing the following steps: QR decomposition of the expanded multiplexed waveform matrix H _E according to the formula H _E = QR, where Q is (N + L) × L matrix, Q ^H Q = I _L, R L × L is a upper triangular matrix, the matrix Q is decomposed into N × L L × L and Q of the _{_{_{Q 1 2, σI L = Q}}} 2 R,

Is the noise sequence in the signal sequence y.

上述计算机程序产品还可以执行初始化有如下方法步骤的程序：根据以下公式计算待传输序列X的第k个信号的初步预估值：

其中，y _k＝R _k，k·x _k+η _k+d _k，

R _k，k为矩阵R中第k行、第k列对应的数据，R _k，j为矩阵R中第k行、第j列对应的数据，x _j为待传输序列X中的第j个元素，η _k为信号序列y中的噪声序列η的第k个元素。 The above computer program product may also execute a program for initializing the following method steps: calculating a preliminary estimate of the kth signal of the sequence X to be transmitted according to the following formula:

Where y _k =R _k,k ·x _k +η _k +d _k ,

上述计算机程序产品还可以执行初始化有如下方法步骤的程序：根据

计算待传输序列X的二次预估值，其中，G _k＝(H _k ^HH _k) ^-1H _k ^H或G _k＝(H _k ^HH _k+σ ²) ^-1H _k ^H，其中，r _k是将预处理后的信号序列r进行干扰抑制得到的。 The above computer program product can also execute a program for initializing the following method steps:

上述计算机程序产品还可以执行初始化有如下方法步骤的程序：在根据置零矩阵计算待传输序列X的二次预估值之前，根据

对预处理后的信号序列r进行干扰抑制，其中，其中，(H) _j表示取H的第j列。 The above computer program product may also execute a program for initializing the following method steps: before calculating the second estimated value of the sequence X to be transmitted according to the zeroing matrix, according to

本领域内的技术人员应明白，本申请的实施例可提供为方法、***、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will appreciate that embodiments of the present application can be provided as a method, system, or computer program product. Thus, the present application can take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment in combination of software and hardware. Moreover, the application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) including computer usable program code.

本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (system), and computer program products according to embodiments of the present application. It will be understood that each flow and/or block of the flowchart illustrations and/or FIG. These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing device to produce a machine for the execution of instructions for execution by a processor of a computer or other programmable data processing device. Means for implementing the functions specified in one or more of the flow or in a block or blocks of the flow chart.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。The computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing device to operate in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture comprising the instruction device. The apparatus implements the functions specified in one or more blocks of a flow or a flow and/or block diagram of the flowchart.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded onto a computer or other programmable data processing device such that a series of operational steps are performed on a computer or other programmable device to produce computer-implemented processing for execution on a computer or other programmable device. The instructions provide steps for implementing the functions specified in one or more of the flow or in a block or blocks of a flow diagram.

在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.

存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。The memory may include non-persistent memory, random access memory (RAM), and/or non-volatile memory in a computer readable medium, such as read only memory (ROM) or flash memory. Memory is an example of a computer readable medium.

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。Computer readable media includes both permanent and non-persistent, removable and non-removable media. Information storage can be implemented by any method or technology. The information can be computer readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read only memory. (ROM), electrically erasable programmable read only memory (EEPROM), flash memory or other memory technology, compact disk read only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, Magnetic tape cartridges, magnetic tape storage or other magnetic storage devices or any other non-transportable media can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary storage of computer readable media, such as modulated data signals and carrier waves.

还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。It is also to be understood that the terms "comprises" or "comprising" or "comprising" or any other variations are intended to encompass a non-exclusive inclusion, such that a process, method, article, Other elements not explicitly listed, or elements that are inherent to such a process, method, commodity, or equipment. An element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in a process, method, article, or device that comprises the element, without further limitation.

本领域技术人员应明白，本申请的实施例可提供为方法、***或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art will appreciate that embodiments of the present application can be provided as a method, system, or computer program product. Thus, the present application can take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment in combination of software and hardware. Moreover, the application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) including computer usable program code.

以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。The above is only an embodiment of the present application and is not intended to limit the application. Various changes and modifications can be made to the present application by those skilled in the art. Any modifications, equivalents, improvements, etc. made within the spirit and scope of the present application are intended to be included within the scope of the appended claims.

Claims

一种接收端信号处理方法，其特征在于，包括：A receiving end signal processing method, comprising:

对复用波形矩阵H进行QR分解；QR decomposition of the multiplexed waveform matrix H;

根据QR分解的结果和预处理后的信号序列r计算信号序列y，其中，所述预处理后的信号序列r是对接收到的信号序列进行预处理得到的；Calculating the signal sequence y according to the result of the QR decomposition and the pre-processed signal sequence r, wherein the pre-processed signal sequence r is obtained by pre-processing the received signal sequence;

根据所述信号序列y计算待传输序列X。The sequence X to be transmitted is calculated from the signal sequence y.
根据权利要求1所述的接收端信号处理方法，其特征在于，The receiving end signal processing method according to claim 1, wherein

对所述复用波形矩阵H进行扩展然后对扩展后复用波形矩阵H _E进行QR分解； Extending the multiplexed waveform matrix H and then performing QR decomposition on the expanded multiplexed waveform matrix H _E ;

预处理后的信号序列r进行扩展后得到扩展后信号序列r _E，根据QR分解的结果和扩展后信号序列r _E计算信号序列y。 The pre-processed signal sequence r is expanded to obtain an extended signal sequence r _E , and the signal sequence y is calculated based on the result of the QR decomposition and the extended signal sequence r _E .
根据权利要求1所述的接收端信号处理方法，其特征在于，The receiving end signal processing method according to claim 1, wherein

根据所述信号序列y计算待传输序列X的初步预估值；Calculating a preliminary estimated value of the sequence X to be transmitted according to the signal sequence y;

根据置零矩阵计算所述待传输序列X的二次预估值。Calculating the second estimated value of the sequence X to be transmitted according to the zeroing matrix.
根据权利要求1所述的接收端信号处理方法，其特征在于，对复用波形矩阵H进行QR分解，包括：The receiving end signal processing method according to claim 1, wherein the QR decomposition of the multiplexed waveform matrix H comprises:

分别初始化矩阵Q、矩阵R、和矩阵P，其中，所述矩阵P是交换矩阵；Initializing matrix Q, matrix R, and matrix P, respectively, wherein the matrix P is a switching matrix;

对初始化后的矩阵Q、矩阵R、和矩阵P进行循环计算，得到最终的矩阵Q、矩阵R、和矩阵P，其中，矩阵Q为酉矩阵，矩阵R为上三角矩阵。The matrix Q, the matrix R, and the matrix P after initialization are cyclically calculated to obtain a final matrix Q, a matrix R, and a matrix P, wherein the matrix Q is a unitary matrix and the matrix R is an upper triangular matrix.
根据权利要求2所述的接收端信号处理方法，其特征在于，对扩展后复用波形矩阵H _E进行QR分解，根据QR分解的结果和所述扩展后信号序列r _E计算信号序列y，包括： The receiving end signal processing method according to claim 2, wherein the extended demultiplexed waveform matrix H _E is QR decomposed, and the signal sequence y is calculated according to the result of the QR decomposition and the expanded signal sequence r _E , including :

按照公式H _E＝QR对所述扩展后复用波形矩阵H _E进行QR分解，其中，Q为(N+L)×L矩阵，Q ^HQ＝I _L，R为L×L的上三角矩阵，矩阵Q分解为N×L的Q ₁和L×L的Q ₂，σI _L＝Q ₂R，
σ ²为噪声功率； The QR decomposition is performed on the extended multiplexed waveform matrix H _E according to the formula H _E = QR, where Q is a (N + L) × L matrix, Q ^H Q = I _L , and R is an upper triangular matrix of L × L , decomposition of a matrix Q N × L of Q ₁ and L × L is _{_{_{Q 2, σI L = Q 2}}} R,
σ ² is the noise power;

根据以下公式计算所述信号序列y：The signal sequence y is calculated according to the following formula:

y＝Q ^Hr _E＝RX-σ ²R ^-HX+Q ₁ ^HN＝RX+η， y=Q ^H r _E =RX-σ ² R ^-H X+Q ₁ ^H N=RX+η,

其中，
为所述信号序列y中的噪声序列。 among them,
Is the noise sequence in the signal sequence y.
根据权利要求3所述的接收端处理信号方法，其特征在于，根据所述信号序列y计算待传输序列X的初步预估值，包括：The receiving end processing signal method according to claim 3, wherein the preliminary estimated value of the sequence X to be transmitted is calculated according to the signal sequence y, comprising:

根据以下公式计算所述待传输序列X的第k个信号的初步预估值：Calculate the preliminary estimate of the kth signal of the sequence X to be transmitted according to the following formula:

其中，y _k＝R _k，k·x _k+η _k+d _k，
R _k，k为矩阵R中第k行、第k列对应的数据，R _k，j为矩阵R中第k行、第j列对应的数据，x _j为所述待传输序列X中的第j个元素，η _k为所述信号序列y中的噪声序列η的第k个元素。7.根据权利要求3所述的接收端信号处理方法，其特征在于，根据置零矩阵计算所述待传输序列X的二次预估值，包括： Where y _k =R _k,k ·x _k +η _k +d _k ,
R _k,k is the data corresponding to the kth row and the kth column in the matrix R, R _k,j is the data corresponding to the kth row and the jth column in the matrix R, and x _j is the number in the sequence X to be transmitted The j elements, η _{k ,} are the _kth element of the noise sequence η in the signal sequence y. The receiving end signal processing method according to claim 3, wherein calculating the second estimated value of the sequence X to be transmitted according to the zeroing matrix comprises:

根据
计算所述待传输序列X的二次预估值，其中，G _k＝(H _k ^HH _k) ^-1H _k ^H或G _k＝(H _k ^HH _k+σ ²) ^-1H _k ^H，其中，r _k是将预处理后的信号序列r进行干扰抑制得到的。 according to
Calculating a quadratic estimate of the sequence X to be transmitted, where G _k =(H _k ^H H _k ) ^-1 H _k ^H or G _k =(H _k ^H H _k +σ ² ) ^-1 H _k ^H Where r _k is obtained by performing interference suppression on the preprocessed signal sequence r.
根据权利要求7所述的接收端信号处理方法，其特征在于，在根据置零矩阵计算所述待传输序列X的二次预估值之前，所述方法还包括：The receiving end signal processing method according to claim 7, wherein before the calculating the second estimated value of the sequence X to be transmitted according to the zeroing matrix, the method further comprises:

根据
对所述预处理后的信号序列r进行干扰抑制，其中，其中，(H) _j表示取H的第j列。 according to
Interference suppression is performed on the preprocessed signal sequence r, wherein (H) _j represents the jth column of H.
根据权利要求4所述的接收端信号处理方法，其特征在于，The receiving end signal processing method according to claim 4, characterized in that

对矩阵Q、矩阵R、矩阵P进行初始化，包括：令R＝0，Q＝H，P＝I _L，所述复用波形矩阵H为N×L的矩阵，矩阵Q为N×L的矩阵，且满足Q ^HQ＝I _L，矩阵I _L为L×L的单位矩阵，矩阵R为L×L的矩阵； The matrix Q, the matrix R, and the matrix P are initialized, including: let R=0, Q=H, P=I _L , the multiplexed waveform matrix H is a matrix of N×L, and the matrix Q is a matrix of N×L And satisfying Q ^H Q=I _L , the matrix I _L is an identity matrix of L×L, and the matrix R is a matrix of L×L;

对初始化后的矩阵Q、矩阵R、矩阵P进行循环计算，得到最终的矩阵Q、矩阵R、矩阵P，包括：在矩阵Q满足第一预设条件的列中查找最小2-范数对应的列；根据查找到的列对矩阵Q中的元素进行归一化和正交化；改变所述第一预设条件，直至矩阵Q的所有列均完成循环。Cycling the initialized matrix Q, the matrix R, and the matrix P to obtain a final matrix Q, a matrix R, and a matrix P, including: finding a minimum 2-norm corresponding in a column in which the matrix Q satisfies the first preset condition Columns; normalize and orthogonalize elements in matrix Q according to the found columns; change the first preset condition until all columns of matrix Q complete the loop.
根据权利要求9所述的接收端信号处理方法，其特征在于，The receiving end signal processing method according to claim 9, wherein

在矩阵Q满足第一预设条件的列中查找最小2-范数对应的列；根据查找到的列对矩阵Q中的元素进行归一化和正交化；改变所述第一预设条件，直至矩阵Q的所有列均完成循环，包括：Finding a column corresponding to the minimum 2-norm in a column in which the matrix Q satisfies the first preset condition; normalizing and orthogonalizing the elements in the matrix Q according to the found column; changing the first preset condition Until all columns of matrix Q complete the loop, including:

步骤1)令i＝1；Step 1) Let i=1;

步骤2)计算Q矩阵中列数l≥i的2-范数，找出其中最小2-范数对应的列数，用k _i表示，即
Step 2) Calculate the 2-norm of the number of columns l ≥ _{i in the} Q matrix, and find the number of columns corresponding to the minimum 2-norm, which is represented by k _i , ie

步骤3)分别交换Q、R、P矩阵的第k _i列与i列；并且令矩阵R中的主对角线元素r _i，i＝||q _i||； K i of the _i-th column and the steps 3) are exchanged Q, R, P matrix; r and causes the main diagonal elements of the matrix R _{_{i, i = || q i ||}} ;

步骤4)根据
对矩阵Q对应的第i列q _i进行归一化； Step 4) According to
Normalizing the i-th column q _i corresponding to the matrix Q;

步骤5)根据q _l＝q _l-r _i，lq _i对q _l进行正交化，r _i，l为q _l在q _i方向上的投影系数，其中，i＜l≤L； Step _{_{5) = q l -r i,}} l q i q _l to be _{orthogonalized,} r _{i, l} is the projection coefficient q _l q _i in the direction according q _l, wherein, i <l≤L;

步骤6)将i+1的值赋值给i；Step 6) assigning the value of i+1 to i;

循环执行步骤2)至步骤6)，直至i＝L，得到最终的矩阵Q、矩阵R。Step 2) to step 6) are performed until i=L, and the final matrix Q and matrix R are obtained.
根据权利要求10所述的接收端信号处理方法，其特征在于，根据所述信号序列y计算待传输序列X，包括：The receiving end signal processing method according to claim 10, wherein calculating the sequence X to be transmitted according to the signal sequence y includes:

根据
k＝L，L-1，..计算所述待传输序列X，其中，
y _k是所述信号序列y的第k个元素，
Q(·)是量化判断算子。 according to
k=L, L-1, .. calculate the sequence X to be transmitted, wherein
y _k is the kth element of the signal sequence y,
Q(·) is a quantization judgment operator.
一种接收端信号处理装置，其特征在于，包括：A receiving end signal processing device, comprising:

分解单元，用于对复用波形矩阵H进行QR分解；a decomposition unit for performing QR decomposition on the multiplexed waveform matrix H;

第一计算单元，用于根据QR分解的结果和预处理后的信号序列r计算信号序列y，其中，所述预处理后的信号序列r是对接收到的信号序列进行预处理得到的；a first calculating unit, configured to calculate a signal sequence y according to the result of the QR decomposition and the pre-processed signal sequence r, wherein the pre-processed signal sequence r is obtained by pre-processing the received signal sequence;

第二计算单元，用于根据所述信号序列y计算待传输序列X。The second calculation unit is configured to calculate the sequence X to be transmitted according to the signal sequence y.
根据权利要求12所述的接收端信号处理装置，其特征在于，The receiving end signal processing device according to claim 12, characterized in that

所述分解单元还用于对扩展后复用波形矩阵H _E进行QR分解，其中，所述扩展后复用波形矩阵H _E是对复用波形矩阵H进行扩展得到的； The decomposing unit is further configured to perform QR decomposition on the extended multiplexed waveform matrix H _E , wherein the extended multiplexed waveform matrix H _E is obtained by expanding the multiplexed waveform matrix H;

所述第一计算单元还用于根据QR分解的结果和扩展后信号序列r _E计算信号序列y，其中，所述扩展后信号序列r _E是将预处理后的信号序列r进行扩展得到的，所述预处理后的信号序列r是对接收到的信号序列进行预处理得到的。 The first calculating unit is further configured to calculate a signal sequence y according to the result of the QR decomposition and the extended signal sequence r _E , wherein the extended signal sequence r _E is obtained by expanding the pre-processed signal sequence r, The pre-processed signal sequence r is obtained by pre-processing the received signal sequence.
根据权利要求12所述的接收端信号处理装置，其特征在于，所述第二计算单元还根据所述信号序列y计算待传输序列X的初步预估值；所述接收端信号处理装置还包括：第三计算单元，用于根据置零矩阵计算所述待传输序列X的二次预估值。The receiving end signal processing apparatus according to claim 12, wherein the second calculating unit further calculates a preliminary estimated value of the sequence X to be transmitted according to the signal sequence y; the receiving end signal processing apparatus further includes And a third calculating unit, configured to calculate a second estimated value of the sequence X to be transmitted according to the zeroing matrix.
根据权利要求13所述的接收端信号处理装置，其特征在于，所述分解单元按照公式H _E＝QR对所述扩展后复用波形矩阵H _E进行QR分解，其中，Q为(N+L)×L矩阵，Q ^HQ＝I _L，R为L×L的上三角矩阵，矩阵Q分解为N×L的Q ₁和L×L的Q ₂，σI _L＝Q ₂R，
σ ²为噪声功率； The receiving end signal processing apparatus according to claim 13, wherein said decomposing unit performs QR decomposition on said expanded multiplexed waveform matrix H _E according to a formula H _E = QR, wherein Q is (N + L) ) × L ^{_{matrix, Q H Q = I L,}} R L × L is a upper triangular matrix, the matrix Q is decomposed into N × L L × L and Q of the _{_{_{Q 1 2, σI L = Q}}} 2 R,
σ ² is the noise power;

所述第一计算单元根据以下公式计算所述信号序列y：The first calculating unit calculates the signal sequence y according to the following formula:

y＝Q ^Hr _E＝RX-σ ²R ^-HX+Q ₁ ^HN＝RX+η，其中，
为所述信号序列y中的噪声序列。 y=Q ^H r _E =RX-σ ² R ^-H X+Q ₁ ^H N=RX+η, where
Is the noise sequence in the signal sequence y.
根据权利要求14所述的接收端信号处理装置，其特征在于，A receiving end signal processing apparatus according to claim 14, wherein:

所述第二计算单元根据以下公式计算所述待传输序列X的第k个信号的初步预估值：
The second calculating unit calculates a preliminary estimated value of the kth signal of the sequence X to be transmitted according to the following formula:

其中，y _k＝R _k，k·x _k+η _k+d _k，
R _k，k为矩阵R中第k行、第k列对应的数据，R _k，j为矩阵R中第k行、第j列对应的数据，x _j为所述待传输序列X中的第j个元素，η _k为所述信号序列y中的噪声序列η的第k个元素。 Where y _k =R _k,k ·x _k +η _k +d _k ,
R _k,k is the data corresponding to the kth row and the kth column in the matrix R, R _k,j is the data corresponding to the kth row and the jth column in the matrix R, and x _j is the number in the sequence X to be transmitted The j elements, η _{k ,} are the _kth element of the noise sequence η in the signal sequence y.
根据权利要求12所述的接收端信号处理装置，其特征在于，所述分解单元包括：The receiving end signal processing apparatus according to claim 12, wherein the decomposition unit comprises:

初始化子单元，用于分别初始化矩阵Q、矩阵R、和矩阵P，其中，所述矩阵P是交换矩阵；Initializing a subunit for initializing a matrix Q, a matrix R, and a matrix P, respectively, wherein the matrix P is a switching matrix;

计算子单元，用于对初始化后的矩阵Q、矩阵R、和矩阵P进行循环计算，得到最终的矩阵Q、矩阵R、和矩阵P，其中，矩阵Q为酉矩阵，矩阵R为上三角矩阵。The calculation subunit is configured to perform cyclic calculation on the initialized matrix Q, the matrix R, and the matrix P to obtain a final matrix Q, a matrix R, and a matrix P, wherein the matrix Q is a unitary matrix and the matrix R is an upper triangular matrix. .
根据权利要求16所述的接收端信号处理装置，其特征在于，A receiving-end signal processing apparatus according to claim 16, wherein:

所述初始化子单元用于：令R＝0，Q＝H，P＝I _L，所述复用波形矩阵H为N×L的矩阵，矩阵Q为N×L的矩阵，且满足Q ^HQ＝I _L，矩阵I _L为L×L的单位矩阵，矩阵R为L×L的矩阵； The initialization subunit is configured to: let R=0, Q=H, P=I _L , the multiplexed waveform matrix H is a matrix of N×L, the matrix Q is a matrix of N×L, and satisfy Q ^H Q =I _L , the matrix I _L is an L × L unit matrix, and the matrix R is an L × L matrix;

所述计算子单元用于：在矩阵Q满足第一预设条件的列中查找最小2-范数对应的列；根据查找到的列对矩阵Q中的元素进行归一化和正交化；改变所述第一预设条件，直至矩阵Q的所有列均完成循环。The calculating subunit is configured to: find a column corresponding to the minimum 2-norm in a column in which the matrix Q satisfies the first preset condition; perform normalization and orthogonalization on the elements in the matrix Q according to the found column; The first predetermined condition is changed until all the columns of the matrix Q complete the loop.
一种存储介质，所述存储介质包括存储的程序，其中，在所述程序运行时控制所述存储介质所在设备执行权利要求1至11任一项所述的接收端信号处理方法。A storage medium comprising a stored program, wherein the device in which the storage medium is located is controlled to perform the receiving end signal processing method according to any one of claims 1 to 11 while the program is running.
一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行权利要求1至11任一项所述的接收端信号处理方法。A processor for executing a program, wherein the program is executed to perform the signal processing method of the receiving end according to any one of claims 1 to 11.