WO2017101097A1

WO2017101097A1 - Channel statistical information obtaining method and receiver

Info

Publication number: WO2017101097A1
Application number: PCT/CN2015/097858
Authority: WO
Inventors: 王悦
Original assignee: 华为技术有限公司
Priority date: 2015-12-18
Filing date: 2015-12-18
Publication date: 2017-06-22
Also published as: CN108370283A

Abstract

Disclosed are a channel statistical information obtaining method and a receiver. The method comprises: a receiver receives, during F time slots, received pilot signals comprising pilot signals sent by a transmitter, wherein F is an integer greater than or equal to 1; the receiver calculates, by using the received F received pilot signals, a covariance matrix of received signals of the receiver over a channel for transmitting the pilot signals; the receiver performs signal reconstruction based on compressive covariance sensing by using the covariance matrix, to obtain a reconstruction signal for sparse statistical information of the channel; the receiver obtains final statistical information of the channel by using the reconstruction signal. Embodiments of the present invention can improve the accuracy of channel statistical information.

Description

一种信道统计信息获取方法和接收机Channel statistical information acquisition method and receiver

技术领域Technical field

本发明涉及通信领域，尤其涉及一种信道统计信息获取方法和接收机。The present invention relates to the field of communications, and in particular, to a channel statistics information acquisition method and receiver.

背景技术Background technique

随着通信技术的发展，目前业界中已有提出使用信道统计信息来最优化预编码设计的方案，该方案可以通过信道统计信息来最优化预编码设计以达到最大传输速率或实现最小误符号率。然而，目前都是采用先执行信道估计再求信道统计信息的方法获取最终的信道统计信息。然而，该方案中对信道统计信息估计的准确性会不高。With the development of communication technology, the scheme of using channel statistics to optimize the precoding design has been proposed in the industry. The scheme can optimize the precoding design through channel statistics to achieve the maximum transmission rate or achieve the minimum symbol error rate. . However, at present, the channel statistics information is obtained by first performing channel estimation and then obtaining channel statistics information. However, the accuracy of estimating channel statistics in this scheme may not be high.

发明内容Summary of the invention

本发明实施例提供了一种信道统计信息获取方法和接收机，可以提高信道统计信息的准确性。The embodiment of the invention provides a channel statistics information acquisition method and a receiver, which can improve the accuracy of channel statistics information.

第一方面，本发明实施例提供一种信道统计信息获取方法，包括：In a first aspect, an embodiment of the present invention provides a method for acquiring channel statistics information, including:

接收机在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数；The receiver receives, in F time slots, a received pilot signal including a pilot signal transmitted by a transmitter, wherein the F is an integer greater than or equal to 1;

接收机使用接收到F个接收导频信号计算接收信号的协方差矩阵，其中，该接收信号是指接收机在传输上述导频信号的信道上的接收信号的协方差矩阵；The receiver calculates a covariance matrix of the received signal by using the F received pilot signals, where the received signal refers to a covariance matrix of the received signal of the receiver on the channel transmitting the pilot signal;

接收机使用上述协方差矩阵进行基于压缩协方差感知的信号重建，以获得上述信道的稀疏统计信息的重建信号；The receiver performs signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of the sparse statistical information of the channel;

接收机使用上述重建信号获取上述信道的最终统计信息。The receiver uses the reconstructed signal described above to obtain the final statistical information of the above channel.

该技术方案中，由于使用包括发射机发送的导频信号的接收导频信号求得信号的统计信息，这样相比现有技术求信道统计信息，可以提高信道统计信息的准确性。In the technical solution, since the statistical information of the signal is obtained by using the received pilot signal including the pilot signal transmitted by the transmitter, the accuracy of the channel statistical information can be improved compared to the prior art for obtaining the channel statistical information.

在第一方面的第一种可能的实现方式中，接收机使用接收到F个接收导频信号计算上述接收信号的协方差矩阵的计算方式可以包括：In a first possible implementation of the first aspect, the receiver uses the received F receive pilots The calculation method of calculating the covariance matrix of the received signal by the signal may include:

接收机将接收到F个接收导频信号进行列向量化，以获取这F个接收导频信号的列向量化向量，并使用F个接收导频信号的列向量化向量计算上述接收信号的协方差矩阵。The receiver receives the F received pilot signals for column vectorization to obtain the column vectorization vectors of the F received pilot signals, and calculates the association of the received signals by using the column vectorization vectors of the F received pilot signals. Variance matrix.

该实现方式中，由于使用上述列向量化向量计算上述接收信号的协方差矩阵，而列向量化向量是在计算时相比矩阵计算复杂要低，从而本实施例中，在计算上述接收信号的协方差矩阵时，可以降低计算复杂度。In this implementation, since the covariance matrix of the received signal is calculated by using the above column vectorization vector, the column vectorization vector is computationally less complex than the matrix calculation, so in the embodiment, the received signal is calculated. When the covariance matrix is used, the computational complexity can be reduced.

结合第一方面或者第一方面的第一种可能的实现方式，在第一方面的第二种可能的实现方式中，接收机使用上述协方差矩阵进行基于压缩协方差感知的信号重建，以获得上述信道的稀疏统计信息的重建信号，可以包括：With reference to the first aspect or the first possible implementation manner of the first aspect, in a second possible implementation manner of the first aspect, the receiver uses the covariance matrix to perform signal reconstruction based on compression covariance sensing to obtain The reconstruction signal of the sparse statistical information of the foregoing channel may include:

接收机以上述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪(Orthogonal Matching Pursuit，OMP)算法的输入，并执行该OMP算法，以获得上述信道的稀疏统计信息的重建信号，其中，上述观测矩阵为预先获取的用于进行信号重建的矩阵，上述L为稀疏度。The receiver uses the column vectorization vector, the observation matrix and L of the above-mentioned covariance matrix as input of an Orthogonal Matching Pursuit (OMP) algorithm, and executes the OMP algorithm to obtain a reconstruction signal of the sparse statistical information of the above channel. The above observation matrix is a matrix obtained for signal reconstruction obtained in advance, and the above L is a sparsity degree.

结合第一方面的第二种可能的实现方式，在第一方面的第三种可能的实现方式中，接收机以上述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行该OMP算法，以获得上述信道的稀疏统计信息的重建信号，可以包括：In conjunction with the second possible implementation of the first aspect, in a third possible implementation manner of the first aspect, the receiver uses the column vectorization vector, the observation matrix, and the L of the covariance matrix as the orthogonal matching tracking OMP. The input of the algorithm and performing the OMP algorithm to obtain the reconstructed signal of the sparse statistical information of the channel may include:

接收机计算上述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积，并将内积的绝对值最大的对角线列向量的索引作为本次迭代索引，其中，所述对角线列向量是所述观测矩阵中与所述重建信号对应的矩阵中对角线元素的索引标识对应的列向量，所述i为本次迭代标号，所述第i减1次迭代的残余向量为所述协方差矩阵的列向量化向量经过i减1次迭代后的残余；The receiver calculates the inner product of each diagonal column vector of the above observation matrix and the residual vector of the ith minus one iteration, and indexes the diagonal column vector with the largest absolute value of the inner product as the index of the iteration, wherein The diagonal column vector is a column vector corresponding to an index identifier of a diagonal element in a matrix corresponding to the reconstructed signal in the observation matrix, where i is the current iteration label, and the ith minus 1 The residual vector of the second iteration is a residual of the column vectorization vector of the covariance matrix after i is decremented by one iteration;

接收机判断所述i是否为1，若是，则将所述本次迭代索引更新至待更新集合中，若否，则将所述本次迭代索引和外推元素在所述重建信号的索引更新至所述待更新集合中，其中，当所述本次迭代索引大于历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引减去行索引差的差值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引加上所述行索引差的和值的元素，当所述本次迭代索引小于所述历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引加上行索引差的和值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引减去所述行索引差的差值的元素，所述行索引差为所述本次迭代索引在所述重建信号对应的矩阵中对应的行向量索引与所述历史迭代索引在所述重建信号对应的矩阵中对应的行向量索引之差的绝对值，所述历史迭代索引为本次迭代之前任意一次迭代的迭代索引；The receiver determines whether the i is 1, if yes, updating the current iteration index to the to-be-updated set, and if not, updating the current iteration index and the extrapolated element in the index of the reconstructed signal Up to the to-be-updated set, wherein, when the current iteration index is greater than a historical iteration index, the extrapolation element includes an index of the reconstruction signal equal to a difference between the current iteration index minus a row index difference An element of a value, and an element further including an index of the reconstruction signal equal to the historical iteration index plus a sum of the row index differences, when the current iteration index is smaller than the historical iteration index, The extrapolation element includes an element whose index of the reconstruction signal is equal to the sum of the current iteration index plus the row index difference, and further includes an index at the reconstruction signal equal to the historical iteration index minus the row An element of the difference of the index difference, where the row index difference is that the current row index corresponding to the current iteration index in the matrix corresponding to the reconstruction signal corresponds to the matrix corresponding to the historical iteration index in the reconstructed signal The absolute value of the difference between the row vector indices, which is an iterative index of any iteration before the iteration;

所述接收机将所述待更新集合中的索引所对应的所述观测矩阵中的列向量构成的部分矩阵和所述协方差矩阵的列向量化向量进行基于最小二乘的重建向量估计，以获得第i次的重建信号；The receiver performs a least squares reconstruction vector estimation based on a partial matrix formed by column vectors in the observation matrix corresponding to an index in the to-be-updated set and a column vectorization vector of the covariance matrix. Obtaining the i-th reconstruction signal;

当所述i等于所述L时，所述接收机将所述第i次的重建信号确定为所述信道的稀疏统计信息的重建信号；When the i is equal to the L, the receiver determines the i-th reconstruction signal as a reconstruction signal of the sparse statistical information of the channel;

当所述i小于所述L时，所述接收机将所述协方差矩阵的列向量化向量减去所述部分矩阵和所述第i次的重建信号的乘积之后的结果作为第i次迭代的残余向量，并将i加1，并触发所述计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积的步骤。When the i is smaller than the L, the receiver subtracts the product of the partial matrix and the ith reconstruction signal from the column vectorization vector of the covariance matrix as the ith iteration a residual vector, and i is incremented by one, and triggers the step of calculating the inner product of each diagonal column vector of the observation matrix and the residual vector of the ith minus one iteration.

该实现方式中，通过上述步骤就可以在执行OMP算法时只对重建信号所对应的矩阵的对角线上的元素进行搜索，以及推算出上述外推元素，从而相比传统的OMP算法，该实施例中，迭代的次数远远小于传统的OMP算法中的迭代次数，以减少计算量。In this implementation manner, when the OMP algorithm is executed, only the elements on the diagonal of the matrix corresponding to the reconstructed signal are searched, and the extrapolated elements are derived, thereby comparing the conventional OMP algorithm. In the embodiment, the number of iterations is much smaller than the number of iterations in the conventional OMP algorithm to reduce the amount of calculation.

结合第一方面的上述任意一种可能的实现方式，在第一方面的第四种可能的实现方式中，接收机使用所述重建信号获取所述信道的最终统计信息，可以包括：With reference to any one of the foregoing possible implementations of the first aspect, in a fourth possible implementation manner of the foregoing aspect, the receiving, by the receiver, the final statistics information of the channel by using the reconstructing signal may include:

接收机根据预先获取的稀疏统计信息与最终统计信息的关系信息获取与所述重建信号对应的最终统计信息作为所述信道的最终统计信息。The receiver acquires final statistical information corresponding to the reconstructed signal as final statistical information of the channel according to the relationship information of the pre-acquired sparse statistical information and the final statistical information.

第二方面，本发明实施例提供一种接收机，该接收机被配置实现上述方法的功能，由硬件/软件实现，其硬件/软件包括与上述功能相应的单元。In a second aspect, an embodiment of the present invention provides a receiver configured to implement the functions of the foregoing method, implemented by hardware/software, and the hardware/software includes a unit corresponding to the foregoing functions.

第三方面，本发明实施例提供一种接收机，包括：处理器、网络接口、存储器和通信总线，其中，所述通信总线用于实现所述处理器、网络接口和存储器之间连接通信，所述处理器执行所述存储器中存储的程序用于实现上述方法。In a third aspect, an embodiment of the present invention provides a receiver, including: a processor, a network interface, a memory, and a communication bus, where the communication bus is used to implement connection communication between the processor, the network interface, and the memory, The processor executes a program stored in the memory for implementing the above law.

另外，本发明实施例中，接收机可以直接使用上述稀疏统计信息的重建信号作为上述信道的最终统计信息。这样快速地求得信号的统计信息。In addition, in the embodiment of the present invention, the receiver can directly use the reconstructed signal of the above-mentioned sparse statistical information as the final statistical information of the channel. This quickly finds the statistical information of the signal.

另外，本发明实施例还提供一种存储一个或多个程序的计算机可读存储介质，所述一个或多个程序包括指令，所述指令当被包括屏幕和多个应用程序的接收机执行时使该接收机执行第一方面提供的任意一种实现方式所述的方法。In addition, embodiments of the present invention also provide a computer readable storage medium storing one or more programs, the one or more programs including instructions when executed by a receiver including a screen and a plurality of applications The receiver is caused to perform the method of any one of the implementations provided by the first aspect.

附图说明DRAWINGS

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the description of the prior art will be briefly described below. Obviously, the drawings in the following description are only It is a certain embodiment of the present invention, and other drawings can be obtained from those skilled in the art without any creative work.

图1是本发明实施例提供的信道统计信息获取方法可应用的***架构图；1 is a system architecture diagram of a method for acquiring channel statistics information provided by an embodiment of the present invention;

图2是本发明实施例提供的一种信道统计信息获取方法的流程示意图；2 is a schematic flowchart of a method for acquiring channel statistics information according to an embodiment of the present invention;

图3是本发明实施例提供的另一种信道统计信息获取方法的流程示意图；3 is a schematic flowchart of another method for acquiring channel statistics information according to an embodiment of the present invention;

图4是本发明实施例提供的一种信号重建的流程示意图；4 is a schematic flowchart of signal reconstruction according to an embodiment of the present invention;

图5至图8是本发明实施例提供的效果示意图；5 to FIG. 8 are schematic diagrams of effects provided by an embodiment of the present invention;

图9是本发明实施例提供的一种接收机的结构示意图；FIG. 9 is a schematic structural diagram of a receiver according to an embodiment of the present invention;

图10是本发明实施例提供的另一种接收机的结构示意图；FIG. 10 is a schematic structural diagram of another receiver according to an embodiment of the present disclosure;

图11是本发明实施例提供的另一种接收机的结构示意图。FIG. 11 is a schematic structural diagram of another receiver according to an embodiment of the present invention.

具体实施方式detailed description

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。The technical solutions in the embodiments of the present invention are clearly and completely described in the following with reference to the accompanying drawings in the embodiments of the present invention. It is obvious that the described embodiments are only a part of the embodiments of the present invention, but not all embodiments. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative efforts are within the scope of the present invention.

请参阅图1，图1是本发明实施例提供的信道统计信息获取方法可应用的***架构图，如图1所示，包括发射机11和接收机12，其中：发射机12可以理解为通信***中进行发送信号的设备，例如：基站、接入点(Access Point， AP)等可以与用户设备进行通信的设备，或者上述发射机11还可以是机器对机器(Machine to Machine，M2M)通信中发送信号的设备等其他通信场景中发送信号的设备。接收机12可以理解为通信***中接收信号的设备，例如：用户设备或者M2M通信中接收信号的设备，其中，用户设备可以包括手机、平板电脑、计算机、可穿戴设备或者车载设备等。Referring to FIG. 1, FIG. 1 is a system architecture diagram of a channel statistical information acquisition method according to an embodiment of the present invention. As shown in FIG. 1, a transmitter 11 and a receiver 12 are included, where the transmitter 12 can be understood as a communication. A device that transmits signals in the system, such as a base station, an access point (Access Point, A device such as an AP that can communicate with a user equipment, or the above-described transmitter 11 can also be a device that transmits signals in other communication scenarios such as a device that transmits signals in Machine to Machine (M2M) communication. The receiver 12 can be understood as a device that receives signals in a communication system, such as a user equipment or a device that receives signals in M2M communication, wherein the user equipment can include a mobile phone, a tablet computer, a computer, a wearable device, or an in-vehicle device.

另外，上述***架构可应用的***可以包括毫米波通信***或者多输入多输出(Multiple-Input Multiple-Output，MIMO)通信***，或者可用于毫米波通信技术与MIMO技术结合的通信***。In addition, the system applicable to the above system architecture may include a millimeter wave communication system or a multiple-input multiple-output (MIMO) communication system, or a communication system that can be used in combination with millimeter wave communication technology and MIMO technology.

在上述***架构中发射机11在F个时隙中向接收机12发送导频信号，从而接收机12在该F个时隙中接收到F个包括导频信号的接收导频信号，这样接收机12就可以使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵，以及使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号，从而可以使用所述重建信号获取所述信道的最终统计信息。In the above system architecture, the transmitter 11 transmits pilot signals to the receiver 12 in F time slots, so that the receiver 12 receives F received pilot signals including pilot signals in the F time slots, thus receiving The machine 12 can calculate a covariance matrix of the received signal of the receiver on the channel transmitting the pilot signal by using the F received pilot signals, and perform compression covariance sensing based on the covariance matrix. The signal is reconstructed to obtain a reconstructed signal of the sparse statistical information of the channel, such that the reconstructed signal can be used to obtain final statistical information of the channel.

另外，本发明实施例中，信道的统计信息可以包括信道的协方差矩阵，这里的信道的协方差矩阵可以称作信道协方差矩阵。In addition, in the embodiment of the present invention, the statistical information of the channel may include a covariance matrix of the channel, where the covariance matrix of the channel may be referred to as a channel covariance matrix.

请参阅图2，图2是本发明实施例提供的一种信道统计信息获取方法的流程示意图，如图2所示，包括以下步骤：Referring to FIG. 2, FIG. 2 is a schematic flowchart of a method for acquiring channel statistics information according to an embodiment of the present invention. As shown in FIG. 2, the method includes the following steps:

201、接收机在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数。201. The receiver receives, in F time slots, a received pilot signal including a pilot signal sent by a transmitter, where the F is an integer greater than or equal to 1.

本实施例中，可以是发射机在分配的F个时隙内向接收机发送上述导频信号，这里的发送可以是F个时隙中每个时隙发送一个导频信号，这样步骤201就可以接收到F个上述接收导频信号。这里需要说明的是，由于接收机在接收时，要经过无线信道且信道中存在噪声信号，所以接收机接收到上述接收导频信号不等于上述导频信号。例如：上述导频信号为压缩训练序列，该压缩训练序列可以为根据某一随机分布生成的随机信号，随机信号可以服从贝努力、高斯等随机分布，另外，本实施例中，每个导频信号都可以为相同的压缩训练序列。其中，压缩训练序列的导频信号可以通过如下公式表示： In this embodiment, the transmitter may send the pilot signal to the receiver in the allocated F time slots, where the transmission may be a pilot signal sent in each time slot of the F time slots, so step 201 may F received the above received pilot signals. It should be noted here that since the receiver is going to pass through the wireless channel and there is a noise signal in the channel, the receiver receives the received pilot signal and is not equal to the pilot signal. For example, the pilot signal is a compressed training sequence, and the compressed training sequence may be a random signal generated according to a random distribution, and the random signal may be subjected to a random distribution such as a Bayer effort, Gaussian, etc. In addition, in this embodiment, each pilot is used. The signals can all be the same compression training sequence. Wherein, the pilot signal of the compressed training sequence can be expressed by the following formula:

其中，X为上述导频信号，U_t是N_t×N_t的离散傅立叶变换(Discrete Fourier Transform，DFT)矩阵，()^H()^-1分别表示对矩阵求共轭转置和对矩阵求逆操作，Z是N_t×T的随机信号矩阵可以服从贝努力、高斯等随机分布，N_t是发射天线个数，T表示在一个训练时长内导频信号所占用的时间资源长度，或者可以理解为导频信号的时长。Where X is the above pilot signal, U _t is a Discrete Fourier Transform (DFT) matrix of N _t ×N _t , and () ^H () ^-1 respectively represents conjugate transposition of the matrix and seeking the matrix In the inverse operation, the random signal matrix with Z is N _t ×T can obey the random distribution such as Bell effort and Gaussian, N _t is the number of transmitting antennas, and T is the length of time resources occupied by the pilot signals in a training duration, or It is understood as the duration of the pilot signal.

这样接收机接收到的接收导频信号就可以通过如下公式表示：Thus, the received pilot signal received by the receiver can be expressed by the following formula:

Y＝HX+WY=HX+W

其中，Y是N_t×T的接收导频信号，X是N_t×T的导频信号，H是N_r×N_t的信道矩阵，N_r是接收天线个数，W是N_r×T的加性噪声。Wherein, Y is the received pilot N _t × _T pilot signal, X is N _t × _T a pilot signal, H is the channel matrix of N _r × N _t is, N _r is the number of receiving antennas, W is N _r × T Additive noise.

202、接收机使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。202. The receiver calculates a covariance matrix of the received signal of the receiver on a channel for transmitting the pilot signal by using the F received pilot signals.

在上述F接收导频信号全部接收完毕后，就可以使用这F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。这里的接收信号可以是理解为接收机在上述信道中接收的发射机中发送的任意信号的接收信号。After all of the F received pilot signals are received, the F received pilot signals can be used to calculate a covariance matrix of the received signals of the receiver on the channel on which the pilot signals are transmitted. The received signal herein may be a received signal that is understood to be any signal transmitted by the receiver in the transmitter received in the above channel.

同样的，这里的接收信号可以用如下公式表示：Similarly, the received signal here can be expressed by the following formula:

Y＝HX+WY=HX+W

其中，Y是接收信号，X是发射机发送的信号，H是N_r×N_t的信道矩阵，N_r是接收天线个数，W是N_r×T的加性噪声。Where Y is the received signal, X is the signal transmitted by the transmitter, H is the channel matrix of N _r × N _t , N _r is the number of receiving antennas, and W is additive noise of N _r ×T.

203、接收机使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号。203. The receiver performs signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of sparse statistical information of the channel.

压缩协方差感知，是一种根据观测对象在压缩域的协方差信息来重建其在非压缩域的协方差信息的信号重建技术。Compression covariance sensing is a signal reconstruction technique that reconstructs covariance information in an uncompressed domain based on the covariance information of the observed object in the compressed domain.

204、接收机使用所述重建信号获取所述信道的最终统计信息。204. The receiver uses the reconstructed signal to obtain final statistical information of the channel.

当获取到上述重建信号后，就可以使用该重建信号获取到上述信道的最终统计信息。例如：预先设定好稀疏统计信息与信道的最终统计信息的关系，那么，就可以通过该关系获得最终统计信息。当然，在一些场景中也可以直接使用上述稀疏统计信息的重建信号作为上述信道的最终统计信息。After the reconstructed signal is obtained, the reconstructed signal can be used to obtain the final statistical information of the channel. For example, if the relationship between the sparse statistical information and the final statistical information of the channel is set in advance, the final statistical information can be obtained through the relationship. Of course, in some scenarios you can also make it directly The reconstructed signal of the above sparse statistical information is used as the final statistical information of the above channel.

本实施例中，接收机在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数；所述接收机使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵；所述接收机使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号；所述接收机使用所述重建信号获取所述信道的最终统计信息。由于使用包括发射机发送的导频信号的接收导频信号求得信号的统计信息，这样相比现有技术求信道统计信息，可以提高信道统计信息的准确性。In this embodiment, the receiver receives the received pilot signal including the pilot signal transmitted by the transmitter in F time slots, wherein the F is an integer greater than or equal to 1; the receiver uses the received F Receiving a pilot signal to calculate a covariance matrix of the received signal of the receiver on a channel transmitting the pilot signal; the receiver performing signal reconstruction based on compression covariance sensing using the covariance matrix to obtain Reconstructing the signal of the sparse statistical information of the channel; the receiver uses the reconstructed signal to obtain final statistical information of the channel. Since the statistical information of the signal is obtained by using the received pilot signal including the pilot signal transmitted by the transmitter, the accuracy of the channel statistical information can be improved compared to the prior art for the channel statistical information.

请参阅图3，图3是本发明实施例提供的另一种信道统计信息获取方法的流程示意图，如图3所示，包括以下步骤：Referring to FIG. 3, FIG. 3 is a schematic flowchart of another method for acquiring channel statistics information according to an embodiment of the present invention. As shown in FIG. 3, the method includes the following steps:

301、接收机在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数。301. A receiver receives a received pilot signal including a pilot signal transmitted by a transmitter in F time slots, wherein the F is an integer greater than or equal to 1.

302、将所述接收到F个接收导频信号进行列向量化，以获取所述F个接收导频信号的列向量化向量，并使用所述F个接收导频信号的列向量化向量计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。302. Perform column vectorization on the received F received pilot signals to obtain a column vectorization vector of the F received pilot signals, and calculate a column vectorization vector using the F received pilot signals. The covariance matrix of the received signal on the channel on which the receiver transmits the pilot signal.

本实施例中，步骤302可以是在获取到上述F个列向量化向量后，分别计算这F个列向量化向量中每个列向量化向量与其共轭转置向量的乘积，再将获取的F个乘积的均值作为上述协方差矩阵。In this embodiment, step 302 may be: after obtaining the F column vectorization vectors, calculate the product of each column vectorized vector and its conjugate transpose vector in the F column vectorization vectors, and then obtain the obtained product. The mean of the F products is taken as the above-described covariance matrix.

或者，本实施例中，步骤302可以通过如下公开计算上述协方差矩阵：Alternatively, in this embodiment, step 302 may calculate the covariance matrix by using the following disclosure:

其中，R_y表示上述协方差矩阵，f表示第f个观测时隙，vec()表示对矩阵进行列向量化处理，y是经列向量化后的接收信号向量即y＝vec(Y)，Y表示上述接收导频信号。Where R _y represents the covariance matrix, f represents the fth observation slot, vec() represents the column vectorization process, and y is the received vector vector after the column vectorization, ie y=vec(Y), Y represents the above received pilot signal.

或者，本实施例中，在一些场景中可以直接从上述F个接收导频信号的列向量化向量选择一个列向量化向量作为上述协方差矩阵的列向量化向量，以获取上述协方差矩阵。 Alternatively, in this embodiment, in some scenarios, a column vectorization vector may be directly selected from the column vectorization vectors of the F received pilot signals as a column vectorization vector of the covariance matrix to obtain the covariance matrix.

步骤302由于使用上述列向量化向量计算上述接收信号的协方差矩阵，而列向量化向量是在计算时相比矩阵计算复杂要低，从而本实施例中，在计算上述接收信号的协方差矩阵时，可以降低计算复杂度。Step 302 calculates the covariance matrix of the received signal by using the above column vectorization vector, and the column vectorization vector is lower than the matrix calculation complexity in the calculation, so in the embodiment, the covariance matrix of the received signal is calculated. When you can reduce the computational complexity.

303、接收机使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号。303. The receiver performs signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of the sparse statistical information of the channel.

本实施例中，步骤303可以包括：In this embodiment, step 303 may include:

接收机以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪(Orthogonal Matching Pursuit，OMP)算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号，其中，所述观测矩阵为预先获取的用于进行信号重建的矩阵，所述L为稀疏度。The receiver uses the column vectorization vector of the covariance matrix, the observation matrix, and L as inputs of an Orthogonal Matching Pursuit (OMP) algorithm, and executes the OMP algorithm to obtain sparse statistical information of the channel. Reconstruction signal, wherein the observation matrix is a pre-acquired matrix for signal reconstruction, and the L is a sparsity.

在OMP算法中存在三个输入参数，即观测矩阵、采集向量和稀疏度，当将这三个输入参数输入到OMP算法中，就可以获得一个重建信号，这里OMP算法的特性，该算法为现有的算法，此处不作详细说明。那么，在该实施方式中，以上述协方差矩阵的列向量化向量、观测矩阵和L作为上述OMP算法的输入，从而获取上述信道的稀疏统计信息的重建信号。There are three input parameters in the OMP algorithm, namely the observation matrix, the acquisition vector and the sparsity. When these three input parameters are input into the OMP algorithm, a reconstructed signal can be obtained. Here, the characteristics of the OMP algorithm are Some algorithms are not described in detail here. Then, in this embodiment, the column vectorization vector, the observation matrix, and L of the covariance matrix are used as inputs of the OMP algorithm to acquire the reconstruction signal of the sparse statistical information of the channel.

另外，上述观测矩阵可以是由Z与U_r来获得，其中，Z是N_t×T的随机信号矩阵，N_t是发射天线个数，U_r是N_r×N_r的DFT矩阵，N_r是接收天线个数。例如：上述观测矩阵可以通过如下公开表示：Further, the observation matrix may be obtained by Z and U _r, wherein, Z is a N _t × _T a random signal matrix, N _t is the number of transmit antennas, U _r is a N _{_r} × N _r matrix of the DFT, N _r Is the number of receiving antennas. For example, the above observation matrix can be publicly expressed as follows:

其中，Φ为上述观测矩阵，

表示克罗内克(Kronecker)乘积运算，()^*表示共轭操作。Where Φ is the above observation matrix,

Denotes Kronecker (the Kronecker) product computation, () ^* denotes a conjugate operator.

这里需要说明的是，该实施方式中，可以采用传统的OMP算法获得上述信道的稀疏统计信息的重建信号。但考虑到传统的OMP算法在进行矩阵搜索是矩阵内的每个位置的元素进行搜索，这样计算量会比较大。而在一些场景中信道的统计信息存在稀疏特性，例如：在毫米波技术与MIMO技术结合***的传播环境中，信道的统计信息存在联合稀疏性。这样本实施例中，可以通过如下步骤获取上述信道的稀疏统计信息的重建信号：It should be noted that, in this embodiment, the reconstructed signal of the sparse statistical information of the above channel may be obtained by using a traditional OMP algorithm. However, considering that the traditional OMP algorithm performs a matrix search for searching for elements at each position in the matrix, the amount of calculation will be relatively large. In some scenarios, the statistical information of the channel has sparsity characteristics. For example, in the propagation environment of the system combining millimeter wave technology and MIMO technology, the statistical information of the channel has joint sparsity. In this embodiment, the reconstruction signal of the sparse statistical information of the channel may be obtained by the following steps:

a)接收机计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积，并将内积的绝对值最大的对角线列向量的索引作为本次迭代索引，其中，所述对角线列向量是所述观测矩阵中与所述重建信号对应的矩阵中对角线元素的索引标识对应的列向量，所述i为本次迭代标号，所述第i减1次迭代的残余向量为所述协方差矩阵的列向量化向量经过i减1次迭代后的残余。a) the receiver calculates the respective diagonal column vector of the observation matrix and the ith of the ith minus one iteration An inner product of the remainder vector, and an index of the diagonal column vector having the largest absolute value of the inner product as the current iteration index, wherein the diagonal column vector is the corresponding to the reconstruction signal in the observation matrix The index of the diagonal element in the matrix identifies the corresponding column vector, the i is the iterative label, and the residual vector of the ith minus 1 iteration is the column vectorization vector of the covariance matrix is reduced by 1 times. Remnants after iteration.

这里需要说明的是，由于这里是对稀疏统计信息进行重建信号，而且稀疏统计信息对应的矩阵维度与发射天线个数和接收天线个数是相关的，那么，在重建之前就可以确定重建信号对应的矩阵的维度，即可以确定重建信号对应的矩阵的行数和列数。另外，本实施例中，上述重建信号可以为列向量化向量，那么重建信号对应的矩阵就为该重建信号列向量进行矩阵化所得的矩阵。上述对角线元素为重建信号对应的矩阵中对角线上的元素。另外，上述对角线列向量是所述观测矩阵中与所述重建信号对应的矩阵中对角线元素的索引标识对应的列向量可以理解为，上述对角线列向量的列向量索引与上述对角线元素在上述重建信号的行向量索引相同，上述对角线元素的索引标识可以理解为对角线元素在上述重建信号中的行向量索引。这里以重建信号对应的矩阵为9乘9的矩阵进行举例说明，那么，上述重建信号就是81乘1的列向量化向量，而上述观测矩阵可以为包含81列的矩阵。这样在重建信号对应的矩阵的对角线元素在重建信号的行向量索引就为1、11、21、31、41、51、61、71、和81，从而在观测矩阵中的对角线列向量也就包括列向量索引就为1、11、21、31、41、51、61、71、和81的列向量。It should be noted here that since the reconstruction signal is performed on the sparse statistical information, and the matrix dimension corresponding to the sparse statistical information is related to the number of transmitting antennas and the number of receiving antennas, the reconstruction signal can be determined before reconstruction. The dimensions of the matrix can be used to determine the number of rows and columns of the matrix corresponding to the reconstructed signal. In addition, in this embodiment, the reconstructed signal may be a column vectorization vector, and then the matrix corresponding to the reconstructed signal is a matrix obtained by matrixing the reconstructed signal column vector. The above diagonal elements are elements on the diagonal in the matrix corresponding to the reconstructed signal. In addition, the diagonal column vector is a column vector corresponding to the index identifier of the diagonal element in the matrix corresponding to the reconstructed signal in the observation matrix, and the column vector index of the diagonal column vector is as described above. The diagonal elements have the same row vector index of the reconstructed signal, and the index of the diagonal elements can be understood as the row vector index of the diagonal elements in the reconstructed signal. Here, the matrix corresponding to the reconstructed signal is 9 by 9 matrix. Then, the reconstructed signal is an 81 by 1 column vectorization vector, and the above observation matrix may be a matrix containing 81 columns. Thus, the diagonal elements of the matrix corresponding to the reconstructed signal are indexed by the row vectors of the reconstructed signal as 1, 11, 21, 31, 41, 51, 61, 71, and 81, thereby diagonal columns in the observation matrix. The vector also includes column vectors whose column vector indices are 1, 11, 21, 31, 41, 51, 61, 71, and 81.

另外，上述步骤以第一次迭代为例进行说明，在第一次迭代中上述残余向量为上述协方差矩阵的列向量化向量。另外，由于残余向量不断更新，从而每次迭代可以获得不同的列向量索引。In addition, the above steps are described by taking the first iteration as an example. In the first iteration, the residual vector is a column vectorization vector of the covariance matrix. In addition, since the residual vectors are continually updated, different column vector indices can be obtained for each iteration.

另外，这里还需要说明一下，上述本次迭代索引当在上述观测矩阵中时而为列向量索引，当在上述重建信号中时则为行向量索引，因为这里的索引仅是一个序号或者数值。In addition, it should be noted here that the current iteration index is a column vector index when in the above observation matrix, and is a row vector index when in the above reconstructed signal, because the index here is only a sequence number or a value.

b)接收机判断所述i是否为1，若是，则将所述本次迭代索引更新至待更新集合中，若否，则将所述本次迭代索引和外推元素在所述重建信号的索引更新至所述待更新集合中，其中，当所述本次迭代索引大于历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引减去行索引差的差值的元素，以及包括在所述重建信号的索引等于所述历史迭代索引加上所述行索引差的和值的元素，当所述本次迭代索引小于所述历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引加上行索引差的和值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引减去所述行索引差的差值的元素，所述行索引差为所述本次迭代索引在所述重建信号对应的矩阵中对应的行向量索引与所述历史迭代索引在所述重建信号对应的矩阵中对应的行向量索引之差的绝对值，所述历史迭代索引为本次迭代之前任意一次迭代的迭代索引。b) the receiver determines whether the i is 1, if yes, updating the current iteration index to the to-be-updated set, and if not, the current iteration index and the extrapolation element in the reconstructed signal Updating an index to the to-be-updated set, wherein when the current iteration index is greater than a historical iteration index, The extrapolation element includes an element at which an index of the reconstructed signal is equal to a difference of the current iteration index minus a row index difference, and an index included in the reconstructed signal is equal to the historical iteration index plus the row And an element of the sum value of the index difference, when the current iteration index is smaller than the historical iteration index, the extrapolation element includes an index of the reconstruction signal equal to the sum of the current iteration index plus the row index difference And an element further included in the index of the reconstructed signal equal to a difference between the historical iteration index minus the row index difference, where the row index difference is that the current iteration index corresponds to the reconstructed signal The corresponding row vector index in the matrix and the absolute value of the difference between the corresponding row vector index of the historical iteration index in the matrix corresponding to the reconstructed signal, the historical iterative index is an iterative index of any iteration before the iteration .

以当所述本次迭代索引小于所述历史迭代索引时为例，上述外推元素包括在所述重建信号的索引等于所述本次迭代索引减去行索引差的差值的元素可以理解为，包括的外推元素在上述重建信号的索引等于本次迭代索引减去行索引差的差值。同理，上述包括在所述重建信号的索引等于所述历史迭代索引加上所述行索引差的和值的元素可以理解为，包括的外推元素在上述重建信号的索引等于本次迭代索引加上述行索引差的和值。这里还是以重建信号对应的矩阵为9乘9的矩阵进行举例说明，For example, when the current iteration index is smaller than the historical iteration index, the extrapolation element includes an element whose index of the reconstruction signal is equal to the difference between the current iteration index minus the row index difference, which can be understood as The index of the extrapolated element included in the above reconstructed signal is equal to the difference between the current iteration index minus the row index difference. Similarly, the above-mentioned element including the index of the reconstructed signal equal to the historical iteration index plus the sum of the row index differences can be understood as the index of the extrapolated element included in the reconstructed signal is equal to the current iteration index. Add the sum of the above-mentioned row index differences. Here, the matrix in which the matrix corresponding to the reconstructed signal is 9 by 9 is exemplified.

若第一次迭代中沿对角线搜索选出的观测矩阵的列向量索引为21，那么，将索引21添加到上述待更新集合中，第一次迭代中不执行对角线外推算；第二次迭代中沿对角线搜索出的观测矩阵的列向量索引为51，将索引51作为第二迭代索引，而索引21和索引51分别在重建信号对应的矩阵对应元素(3,3)和元素(6,6)，而元素(3,3)和元素(6,6)的行索引差为3，那么，就可以外推出51减3等于48的索引，以及21加3等于24的索引，其中，索引48和索引24在重建信号对应的矩阵中对应元素(3,6)和元素(6,3)；第三次迭代中沿对角线搜索出的观测矩阵的列向量索引为81，那么，将索引81作为第三迭代索引，而索引21、索引51和索引81分别在重建信号对应的矩阵对应元素(3,3)、元素(6,6)和元素(9,9)，其中，索引21、索引51为历史迭代索引，且元素(3,3)与元素(9,9)的行索引差为6，而元素(6,6)与元素(9,9)的行索引差为3，那么第三次迭代就可以外推出81减3等于78的索引，以及51加3等于54的索引，以及外推出81减6等于75的索引，以及21加6等于27的索引，其中，索引 78、索引54、索引75和索引27在重建信号对应的矩阵中对应元素(6,9)、元素(9,6)、元素(3,9)和元素(9,3)。If the column vector index of the selected observation matrix along the diagonal search in the first iteration is 21, then the index 21 is added to the above-mentioned to-be-updated set, and the diagonal extrapolation is not performed in the first iteration; In the second iteration, the column vector index of the observation matrix searched along the diagonal line is 51, the index 51 is taken as the second iteration index, and the index 21 and the index 51 are respectively corresponding to the matrix corresponding elements (3, 3) and the corresponding signals of the reconstructed signal. The element (6,6), and the element (3,3) and the element (6,6) have a row index difference of 3, then, it is possible to extrapolate the index of 51 minus 3 equal to 48, and the index of 21 plus 3 equals 24 Wherein the index 48 and the index 24 correspond to the element (3, 6) and the element (6, 3) in the matrix corresponding to the reconstructed signal; the column vector index of the observation matrix searched along the diagonal in the third iteration is 81 Then, the index 81 is used as the third iteration index, and the index 21, the index 51, and the index 81 respectively correspond to the matrix corresponding elements (3, 3), elements (6, 6), and elements (9, 9) corresponding to the reconstruction signal. Wherein, the index 21 and the index 51 are historical iteration indexes, and the row index difference between the element (3, 3) and the element (9, 9) is 6, and the element (6, 6) and the element ( 9,9) the row index difference is 3, then the third iteration can be extrapolated 81 minus 3 equal to 78 index, and 51 plus 3 is equal to 54 index, and the extrapolation 81 minus 6 is equal to 75 index, and 21 Add 6 to the index of 27, where the index 78. Index 54, index 75, and index 27 correspond to elements (6, 9), elements (9, 6), elements (3, 9), and elements (9, 3) in the matrix corresponding to the reconstructed signal.

c)接收机将所述待更新集合中的索引所对应的所述观测矩阵中的列向量构成的部分矩阵和所述协方差矩阵的列向量化向量进行基于最小二乘的重建向量估计，以获得第i次的重建信号。c) the receiver performs a least squares reconstruction vector estimation based on a partial matrix composed of column vectors in the observation matrix corresponding to the index in the set to be updated and a column vectorization vector of the covariance matrix Obtain the i-th reconstruction signal.

通过步骤a)和b)就可以获取包括索引的待更新集合，从而可以使用该待更新集合内的索引所对应的所述观测矩阵中的列向量构成部分矩阵，再将该部分矩阵与协方差矩阵的列向量化向量进行基于最小二乘的重建向量估计，以获得第i次的重建信号。这里的部分矩阵可以是将上述待更新集合中的索引作为部分矩阵的列向量索引，即构造包括索引为上述待更新集合中的索引的列向量的部分矩阵。The set to be updated including the index can be obtained through steps a) and b), so that the partial matrix can be formed by using the column vector in the observation matrix corresponding to the index in the to-be-updated set, and then the partial matrix and the covariance can be formed. The column vectorization vector of the matrix performs a least squares reconstruction vector estimation to obtain the i-th reconstruction signal. The partial matrix here may be a column vector index that takes an index in the above-mentioned set to be updated as a partial matrix, that is, constructs a partial matrix including a column vector whose index is an index in the above-mentioned set to be updated.

d)当所述i等于所述L时，所述接收机将所述第i次的重建信号确定为所述信道的稀疏统计信息的重建信号。d) when the i is equal to the L, the receiver determines the ith reconstruction signal as a reconstruction signal of the sparse statistical information of the channel.

e)当所述i小于所述L时，所述接收机将所述协方差矩阵的列向量化向量减去所述部分矩阵和所述第i次的重建信号的乘积之后的结果作为第i次迭代的残余向量，并将i加1，并触发所述计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积的步骤。e) when the i is smaller than the L, the receiver subtracts the product of the partial matrix and the ith reconstruction signal from the column vectorization vector of the covariance matrix as the ith The residual vector of the next iteration, and increments i, and triggers the step of calculating the inner product of each diagonal column vector of the observation matrix and the residual vector of the ith minus one iteration.

通过上述a)至e)的步骤就可以在执行OMP算法时只对重建信号所对应的矩阵的对角线上的元素进行搜索，以及推算出上述外推元素，从而相比传统的OMP算法，该实施例中，迭代的次数远远小于传统的OMP算法中的迭代次数，以减少计算量。另外，上述介绍的OMP算法可以定义为基于对角线搜索的OMP(Diagonal Search Orthogonal Matching Pursuit，DS-OMP)算法。Through the steps a) to e) above, it is possible to search only the elements on the diagonal of the matrix corresponding to the reconstructed signal when performing the OMP algorithm, and to derive the above extrapolated elements, thereby comparing the conventional OMP algorithm. In this embodiment, the number of iterations is much smaller than the number of iterations in the conventional OMP algorithm to reduce the amount of computation. In addition, the OMP algorithm described above may be defined as an OMP (Diagonal Search Orthogonal Matching Pursuit, DS-OMP) algorithm based on diagonal search.

304、接收机使用所述重建信号获取所述信道的最终统计信息。304. The receiver uses the reconstructed signal to obtain final statistical information of the channel.

本实施例中，可以预先获取稀疏统计信息与最终统计信息的关系信息，这样步骤304可以包括：In this embodiment, the relationship information between the sparse statistical information and the final statistical information may be obtained in advance, and the step 304 may include:

例如：稀疏统计信息与最终统计信息存在如下关系： For example, sparse statistics have the following relationship with the final statistics:

其中，上述

为上述最终统计信息，

为上述稀疏统计信息的重建信号且满足

即可以理解为将向量

矩阵化可得到

U_t是N_t×N_t的DFT矩阵，U_r是N_r×N_r的DFT矩阵，

表示克罗内克(Kronecker)乘积运算，()^*表示共轭操作。Among them, the above

For the above final statistics,

Reconstructing the signal for the above sparse statistical information and satisfying

Can be understood as a vector

Matrixization is available

U _t is a DFT matrix of N _t ×N _t , and U _r is a DFT matrix of N _r ×N _r ,

Represents the Kronecker product operation, and () ^* indicates the conjugate operation.

通过上述公式可以是接收机可以是在执行完基于压缩协方差感知的信号重建算法后，接收机根据信号重建结果来获得信道的统计信息

再将向量

矩阵化为

再获取上述信道的最终统计信息

The above formula may be that the receiver may obtain the channel statistical information according to the signal reconstruction result after performing the signal reconstruction algorithm based on compression covariance sensing.

Vector

Matrixization

Obtain the final statistics of the above channels.

其中，上述

的说明如下：Among them, the above

The description is as follows:

在理想场景或者在无噪声场景下接收机接收的接收信号的协方差矩阵的列向量r_y可以变换成如下公式表示：The column vector r _{y of the} covariance matrix of the received signal received by the receiver in an ideal scene or in a noiseless scene can be transformed into the following formula:

其中，U_r是N_r×N_r的DFT矩阵，

表示Kronecker乘积运算，()^*表示共轭操作，R_h表示信道的协方差矩阵：而R_h又可以通过如下公式表示：Where U _r is a DFT matrix of N _r ×N _r ,

Represents the Kronecker product operation, () ^* represents the conjugate operation, and R _h represents the covariance matrix of the channel: and R _h can be expressed by the following formula:

R_h＝E{hh^H}R _h =E{hh ^H }

其中，E{}表示求数学期望，h＝vec(H)，H为信道矩阵，且H又可以经信道分解表示为：Where E{} denotes the mathematical expectation, h=vec(H), H is the channel matrix, and H can be expressed by channel decomposition as:

其中，H_v是一个稀疏的矩阵，这是因为在一些传播环境中，例如：毫米波技术与MIMO技术结合的传播环境中，由于有效散射体的个数是有限的，使得有效多径信道的个数L远小于信道的维度，即信道稀疏度是L。Where H _v is a sparse matrix, because in some propagation environments, such as the millimeter wave technology combined with MIMO technology, the number of effective scatterers is limited, making the effective multipath channel The number L is much smaller than the dimension of the channel, that is, the channel sparsity is L.

另外，根据信道矩阵H的分解表示和克罗内克(Kronecker)乘积的展开式，信道协方差矩阵R_h可以进一步表示为：In addition, according to the decomposition representation of the channel matrix H and the expansion of the Kronecker product, the channel covariance matrix _Rh can be further expressed as:

其中，

h_v＝vec(H_v)。 among them,

h _v =vec(H _v ).

由此可见，由于信道协方差矩阵

包含了多个时隙内稀疏信道h_v＝vec(H_v)的统计信息，所以各个时隙内的信道的独立稀疏性通过

表现为了多时隙信道间的联合稀疏性。Thus, due to the channel covariance matrix

Contains statistical information of the sparse channel h _v =vec(H _v ) in multiple time slots, so the independent sparsity of the channels in each time slot passes

Expressed as joint sparsity between multi-slot channels.

将上式R_h带入Bring the above formula R _h into

进一步有：Further there are:

另外，进一步简化数学符号，令

而且

则上式可以表示为

In addition, further simplify the mathematical symbols, so

and

Then the above formula can be expressed as

需要说明的是，上述推导仅是公式之间的推导，接收机在实际应用中不需要进行上述推导，可以获取到如下公开并进行使用：It should be noted that the above derivation is only a derivation between formulas, and the receiver does not need to perform the above derivation in practical applications, and can be obtained and disclosed as follows:

另外，需要说明的是，本实施例中计算信道的最终统计信息时并不限定使用上述公式进行计算，也可以通过大量实验数据得到上述重建信号与最终统计信息之间的关系，再使用该关系获取到信道的最终统计信息等。In addition, it should be noted that, in the present embodiment, when calculating the final statistical information of the channel, the calculation is not limited by using the above formula, and the relationship between the reconstructed signal and the final statistical information may also be obtained through a large amount of experimental data, and then the relationship is used. Get the final statistics of the channel, etc.

下面以图4为例对接收机使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号的步骤进行举例说明：The following takes FIG. 4 as an example to illustrate the steps of the receiver using the covariance matrix to perform signal reconstruction based on compression covariance sensing to obtain a reconstruction signal of the sparse statistical information of the channel:

且在该实例接收机进行基于压缩协方差感知的信号重建就转化为了如下数学问题：And the signal reconstruction based on compression covariance sensing in the example receiver is transformed into the following mathematical problem:

其中，arg min表示求解令目标函数最小的

作为重建结果，|| ||₁|| ||₂分别表示向量的1范数和2范数，其中，1范数是对联合稀疏性的部分体现，而2范数是对加性噪声的约束，λ表示拉格朗日参数。Where arg min indicates that the solution makes the objective function the smallest

As a result of the reconstruction, || || ₁ || || ₂ represent the 1 norm and 2 norm of the vector, respectively, where the 1 norm is part of the joint sparsity, and the 2 norm is the additive noise. Constraint, λ represents the Lagrangian parameter.

该实例具体通过图4所示的步骤求出上述

即获得到上述稀疏统计信息的重建信号，如图4所示，包括：This example is specifically obtained by the steps shown in FIG.

That is, the reconstruction signal to the above sparse statistical information is obtained, as shown in FIG. 4, including:

1、输入。其中，该步骤可以是输入r_y,Φ,L，其中，L为稀疏度，或者理解为迭代次数，或者理解为上述

所对应的矩阵

中对角线上的非零元素的个数。1. Enter. Wherein, the step may be inputting r _y , Φ, L, where L is a sparsity degree, or is understood as the number of iterations, or is understood as the above

Corresponding matrix

The number of non-zero elements on the diagonal.

2、初始化。该步骤中可以是将迭代标号初始化i＝0，以及重建向量初始化r⁽⁰⁾＝0(其中加粗的零代表元素值全为零的零向量)，以及将残余向量初始化u⁽⁰⁾＝r_y，以及集合初始化

其中Γ,Λ分别表示待更新集合以及对角线搜索空间。2, initialization. In this step, the iteration label may be initialized with i=0, and the reconstruction vector initialization r ⁽⁰⁾ =0 (where the thickened zero represents a zero vector whose element values are all zero), and the residual vector is initialized u ⁽⁰⁾ = r _y , and set initialization

Where Γ, Λ respectively represent the set to be updated and the diagonal search space.

3、定义对角线搜索空间。该步骤可以是将后续搜索空间限定在信道协方差矩阵的对角线上的各个元素位置，而不是整个信道协方差矩阵的全部元素位置。即只对信道协方差矩阵的对角线上的非零元素进行搜索。上述信道协方差矩阵为上述

对应的矩阵

或者上述

进行矩阵化后的矩阵

3. Define the diagonal search space. This step may be to limit subsequent search spaces to individual element locations on the diagonal of the channel covariance matrix, rather than all element locations of the entire channel covariance matrix. That is, only non-zero elements on the diagonal of the channel covariance matrix are searched. The above channel covariance matrix is as described above

Corresponding matrix

Or the above

Matrix after matrixing

4、沿对角线搜索。该步骤中在每一次迭代中，首先更新迭代标号，即i＝i+1，然后在已定义的对角线搜索空间内，寻找矩阵Φ的列向量索引，该索引编号t⁽ⁱ⁾(其中，上脚标⁽ⁱ⁾表示当前的迭代标号)为4. Search along the diagonal. In each iteration, in this iteration, the iteration label is first updated, ie i=i+1, and then in the defined diagonal search space, the column vector index of the matrix Φ is searched, the index number t ⁽ⁱ⁾ (where , the footer ⁽ⁱ⁾ indicates the current iteration label)

其中，arg max|<u^(i-1),Φ_.,j>|表示搜索令上一次迭代更新的残余向量u和观测矩阵Φ的列向量的内积的绝对值最大的列向量索引，而j∈Λ表示搜索范围仅仅限定在对角线的元素位置上。Where arg max|<u ^(i-1) , Φ _{., j} >| represents the column vector index with the largest absolute value of the inner product of the column vector of the last iteration update and the column vector of the observation matrix Φ. j∈Λ indicates that the search range is only limited to the position of the elements of the diagonal.

通过步骤可以实现计算矩阵Φ的各个对角线列向量与第i减1次迭代的残余向量的内积，并将内积的绝对值最大的对角线列向量的索引作为本次迭代索引。Through the steps, the inner product of each diagonal line vector of the calculation matrix Φ and the residual vector of the ith minus one iteration can be realized, and the index of the diagonal column vector with the largest absolute value of the inner product is used as the index of the current iteration.

5、判断当前迭代是否是首次迭代，若是则直接进行步骤7，否则执行步骤6。5. Determine whether the current iteration is the first iteration. If yes, proceed directly to step 7, otherwise go to step 6.

6、进行对角线外推算。对角线外推算是指，根据之前搜索到的对角线元素位置和当前搜索到的对角线元素位置来获得对应对角线外的元素位置，这一操作是利用了信道协方差矩阵的埃尔米特结构特性，对角线上的非零元素反映了稀疏信道的自相关特性，而对应对角线外的非零元素反映了稀疏信道的互相关特性。举个简单的对角线外推算的例子以作示例，在信道协方差矩阵为9乘9的矩阵中，若第一次迭代出在上述矩阵Φ中的列向量索引21，而索引21对应于协方差矩阵(上述

进行矩阵化后的矩阵)中的元素(3,3)，即第一次迭代沿对角线搜索选出的协方差矩阵的对角线元素位置为(3,3)，即搜索出在协方差矩阵中(3,3)的元素为非零元素，第一次迭代中不执行对角线外推算；第二次迭代出在上述矩阵Φ中的列向量索引51，而索引51对应于协方差矩阵中的元素(6,6)，即第二次迭代中沿对角线搜索选出的协方差矩阵的对角线元素位置为(6,6)，第二次对角线外推算出的协方差矩阵的对角线外元素位置为(3,6)(6,3)，而元素(3,6)和(6,3)在上述

中的列向量索引为48和24；第三次迭代出在上述矩阵Φ中的列向量索引81，而索引81对应于协方差矩阵中的元素(9,9)，即第三次迭代中沿对角线搜索选出的协方差矩阵的对角线元素位置为(9,9)，第三次对角线外推算出的协方差矩阵的对角线外元素位置为(3,9)(9,3)(6,9)(9,6)，而元素(3,9)(9,3)(6,9)(9,6)在上述

中的列向量索引为75、27、78和54，若有后续迭代以此类推。6. Perform diagonal extrapolation. Diagonal extrapolation refers to obtaining the position of the element outside the diagonal line based on the previously searched diagonal element position and the currently searched diagonal element position. This operation utilizes the channel covariance matrix. Hermite's structural characteristics, the non-zero elements on the diagonal reflect the autocorrelation properties of the sparse channel, while the non-zero elements outside the diagonal reflect the cross-correlation properties of the sparse channel. As an example, a simple diagonal extrapolation algorithm is used. In a matrix with a channel covariance matrix of 9 by 9, if the first iteration is performed, the column vector index 21 in the matrix Φ is obtained, and the index 21 corresponds to Covariance matrix

The element (3, 3) in the matrix after matrixing, that is, the position of the diagonal element of the covariance matrix selected by the diagonal search in the first iteration is (3, 3), that is, the search is in the association. The elements of (3,3) in the variance matrix are non-zero elements, and the diagonal extrapolation is not performed in the first iteration; the second iteration produces the column vector index 51 in the above matrix Φ, and the index 51 corresponds to the association The elements (6, 6) in the variance matrix, that is, the diagonal element positions of the covariance matrix selected along the diagonal search in the second iteration are (6, 6), and the second diagonal extrapolation is calculated. The position of the diagonal element of the covariance matrix is (3,6) (6,3), while the elements (3,6) and (6,3) are above

The column vector indices in the index are 48 and 24; the third iteration produces the column vector index 81 in the above matrix Φ, and the index 81 corresponds to the element in the covariance matrix (9, 9), that is, in the third iteration The position of the diagonal element of the covariance matrix selected by the diagonal search is (9, 9), and the position of the diagonal element of the covariance matrix calculated by the third diagonal extrapolation is (3, 9) ( 9,3)(6,9)(9,6), and the elements (3,9)(9,3)(6,9)(9,6) are above

The column vector indices in the are 75, 27, 78, and 54, if there are subsequent iterations and so on.

7、集合更新。该步骤可以是将本次迭代中对角线搜索和对角线外推算获得的索引更新到待更新集合Γ中。7, collection update. This step may be to update the index obtained by the diagonal search and the diagonal extrapolation in this iteration into the to-be-updated set.

8、最小二乘估计。该步骤可以是根据更新后的集合Γ中所对应的列向量索引所构成的部分矩阵

和r_y来进行基于最小二乘的重建向量估计，具体通过如下公式进行最小二乘的重建向量估计8. Least squares estimation. The step may be a partial matrix formed by the column vector index corresponding to the updated set Γ

And r _y to perform the least squares reconstruction vector estimation, and the least squares reconstruction vector estimation is performed by the following formula

其中，下脚标Γ⁽ⁱ⁾和.,Γ⁽ⁱ⁾分别表示向量的对应元素位置和矩阵的对应列向量位置，

表示对矩阵

求伪逆运算。另外，这里的需要说明的是，由于上述

为列向量化向量，那么，在上述

中的索引就是在上述

中的位置。Wherein, the lower footmarks ⁽ⁱ⁾ and ., Γ ⁽ⁱ⁾ respectively represent the corresponding element position of the vector and the corresponding column vector position of the matrix,

Representation matrix

Pseudo-inverse operation. In addition, the need to explain here is due to the above

Vectorizing the vector for the column, then, above

The index in is in the above

The location in .

9、残余更新。该步骤可以是将当前迭代的最后一步，即对残余向量u⁽ⁱ⁾进行残余更新，例如：通过如下公式表示：9. Residual updates. This step may be the last step of the current iteration, ie the residual update of the residual vector u ⁽ⁱ⁾ , for example: by the following formula:

10、输出。该步骤是将最后一次(第L次)迭代中的最小二乘重建向量作为最终重建结果

并输出。10. Output. This step is to use the least squares reconstruction vector in the last (Lth) iteration as the final reconstruction result.

And output.

通过上述步骤就可以获取到上述

接收机根据信号重建结果来获得信道的统计信息，即将向量

矩阵化为

再根据下式获得对信道协方差矩阵的最终估计Through the above steps, you can get the above

The receiver obtains the statistical information of the channel according to the signal reconstruction result, that is, the vector

Matrixization

Then obtain the final estimate of the channel covariance matrix according to the following formula.

本实施例中，在图2所示的实施例的基础上增加了多种可选的实施方式，且都可以实现提高信道统计信息的准确性。In this embodiment, a plurality of optional implementation manners are added on the basis of the embodiment shown in FIG. 2, and the accuracy of the channel statistics information can be improved.

下面请参阅图5至图8，图5至图8是本发明实施例与估计信道统计信息的现有技术进行对比的实验数据图表。Referring to FIG. 5 to FIG. 8, FIG. 5 to FIG. 8 are diagrams of experimental data in comparison with the prior art for estimating channel statistical information according to an embodiment of the present invention.

其中，图5至图8中Pr表示正确信道统计信息估计的概率，SNR表示信噪比，对比现有技术(虚线方块)，本发明实施例提供的技术方案(实线三角)的性能更优。例如：在SNR为10时，现有技术的Pr为0.4左右，而本发明实施例提供的技术方案在SNR为10时，Pr可以达到0.7。In FIG. 5 to FIG. 8 , Pr represents the probability of correct channel statistical information estimation, and SNR represents a signal to noise ratio. Compared with the prior art (dashed square), the technical solution (solid triangle) provided by the embodiment of the present invention has better performance. . For example, when the SNR is 10, the Pr of the prior art is about 0.4, and the technical solution provided by the embodiment of the present invention can achieve a Pr of 0.7 when the SNR is 10.

另外，在图6中L表示稀疏度，通过图6可知，随着信道稀疏性的变差(即“L”增加，信道从特别稀疏变为一般稀疏)，本发明实施例提供的技术方案的性能变化差异比较平缓，且性能变化远远小于现有技术的性能衰减程度，换言之，本发明实施例提供的技术方案更具更好的稀疏信道环境鲁棒性。In addition, in FIG. 6, L represents the sparsity, and as shown in FIG. 6, as the channel sparsity deteriorates (ie, the "L" increases, the channel changes from a particularly sparse to a generally sparse), the technical solution provided by the embodiment of the present invention The performance variation is relatively flat, and the performance change is far less than the performance degradation degree of the prior art. In other words, the technical solution provided by the embodiment of the present invention has better robustness of the sparse channel environment.

在图7和图8中，T表示在每个时隙内导频信号所占用的时间资源长度，F表示导频信号的观测时隙的个数，或者两者也可以理解为估计信道的统计信息所使用的观测开销。通过图7和图8可知，在获得相同估计准确性的前提下，本发明实施例提供的技术方案所用的观测开销比现有技术方案所用的观测开销更少，不论是每个时隙内的时间资源长度“T”，还是观测时隙的个数“F”都比现有技术方案所用的开销更少。In FIG. 7 and FIG. 8, T denotes the time resource length occupied by the pilot signal in each slot, F denotes the number of observation slots of the pilot signal, or both can be understood as the estimated channel statistics. The observation overhead used by the information. It can be seen from FIG. 7 and FIG. 8 that, under the premise that the same estimation accuracy is obtained, the observation overhead used by the embodiments of the present invention is less than the observation overhead used in the prior art scheme, regardless of the time slot. The time resource length "T", or the number of observation slots "F", is less than the overhead used in prior art solutions.

而且，由于本发明实施例提供的技术方案只需进行一次信号重建操作，这样本发明实施例中计算的复杂度就比较低，且本发明实施例中还可以采用DS-OMP算法，而该算法中只在对角线上进行搜索，从而可以更进一步降低计算杂度，由此带来的好处是可大幅降低***的硬件复杂度和功耗。Moreover, since the technical solution provided by the embodiment of the present invention only needs to perform a signal reconstruction operation, the complexity of the calculation in the embodiment of the present invention is relatively low, and the DS-OMP algorithm can also be used in the embodiment of the present invention, and the algorithm Only search on the diagonal, so you can further reduce the meter The benefit of this is that the hardware complexity and power consumption of the system can be greatly reduced.

下面为本发明装置实施例，本发明装置实施例用于执行本发明方法实施例一至二实现的方法，为了便于说明，仅示出了与本发明实施例相关的部分，具体技术细节未揭示的，请参照本发明实施例一和实施例二。The following is a device embodiment of the present invention. The device embodiment of the present invention is used to perform the method for implementing the first to second embodiments of the present invention. For the convenience of description, only parts related to the embodiment of the present invention are shown, and the specific technical details are not disclosed. Please refer to Embodiment 1 and Embodiment 2 of the present invention.

请参阅图9，图9是本发明实施例提供的一种接收机的结构示意图，如图9所示，包括：接收单元91、计算单元92、重建单元93和获取单元94，其中：Referring to FIG. 9, FIG. 9 is a schematic structural diagram of a receiver according to an embodiment of the present invention. As shown in FIG. 9, the method includes: a receiving unit 91, a calculating unit 92, a reconstructing unit 93, and an obtaining unit 94, where:

接收单元91，用于在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数。The receiving unit 91 is configured to receive, in the F time slots, a received pilot signal including a pilot signal sent by a transmitter, where the F is an integer greater than or equal to 1.

所述计算单元92，用于使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。The calculating unit 92 is configured to calculate a covariance matrix of the received signal of the receiver on a channel for transmitting the pilot signal by using the F received pilot signals.

所述重建单元93，用于使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号。The reconstruction unit 93 is configured to perform signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of the sparse statistical information of the channel.

所述获取单元94，用于使用所述重建信号获取所述信道的最终统计信息。The obtaining unit 94 is configured to obtain final statistical information of the channel by using the reconstruction signal.

请参阅图10，图10是本发明实施例提供的另一种接收机的结构示意图，如图10所示，包括：接收单元101、计算单元102、重建单元103和获取单元104，其中：Referring to FIG. 10, FIG. 10 is a schematic structural diagram of another receiver according to an embodiment of the present invention. As shown in FIG. 10, the method includes: a receiving unit 101, a calculating unit 102, a reconstructing unit 103, and an obtaining unit 104, where:

接收单元101，用于在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数。The receiving unit 101 is configured to receive, in the F time slots, the receiving of the pilot signal sent by the transmitter. a pilot signal, wherein the F is an integer greater than or equal to one.

计算单元102，用于计算单元用于将所述接收到F个接收导频信号进行列向量化，以获取所述F个接收导频信号的列向量化向量，并使用所述F个接收导频信号的列向量化向量计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。The calculating unit 102 is configured to perform, by the computing unit, performing column vectorization on the received F received pilot signals, to obtain a column vectorization vector of the F received pilot signals, and using the F receiving guides A column vectorization vector of the frequency signal calculates a covariance matrix of the received signal of the receiver on the channel transmitting the pilot signal.

重建单元103，用于使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号。The reconstruction unit 103 is configured to perform signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of the sparse statistical information of the channel.

本实施例中，重建单元103可以用于以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号，其中，所述观测矩阵为预先获取的用于进行信号重建的矩阵，所述L为稀疏度。In this embodiment, the reconstruction unit 103 may be configured to use the column vectorization vector, the observation matrix, and L of the covariance matrix as an input of the orthogonal matching tracking OMP algorithm, and execute the OMP algorithm to obtain the channel. A reconstruction signal of sparse statistical information, wherein the observation matrix is a pre-acquired matrix for performing signal reconstruction, and the L is a sparsity degree.

该实施方式中，重建单元103可以包括：In this embodiment, the reconstruction unit 103 may include:

第一计算子单元1031，用于计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积，并将内积的绝对值最大的对角线列向量的索引作为本次迭代索引，其中，所述对角线列向量是所述观测矩阵中与所述重建信号对应的矩阵中对角线元素的索引标识对应的列向量，所述i为本次迭代标号，所述第i减1次迭代的残余向量为所述协方差矩阵的列向量化向量经过i减1次迭代后的残余；a first calculating sub-unit 1031, configured to calculate an inner product of each diagonal line vector of the observation matrix and a residual vector of the ith minus one iteration, and to maximize a diagonal column vector of an absolute value of the inner product An index is used as the current iteration index, wherein the diagonal column vector is a column vector corresponding to an index identifier of a diagonal element in a matrix corresponding to the reconstructed signal in the observation matrix, and the i is an iteration a label, the residual vector of the ith minus one iteration is a residual of a column vectorization vector of the covariance matrix after i is decremented by one iteration;

判断单元1032，用于判断所述i是否为1，若是，则将所述本次迭代索引更新至待更新集合中，若否，则将所述本次迭代索引和外推元素在所述重建信号的索引更新至所述待更新集合中，其中，当所述本次迭代索引大于历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引减去行索引差的差值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引加上所述行索引差的和值的元素，当所述本次迭代索引小于所述历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引加上行索引差的和值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引减去所述行索引差的差值的元素，所述行索引差为所述本次迭代索引在所述重建信号对应的矩阵中对应的行向量索引与所述历史迭代索引在所述重建信号对应的矩阵中对应的行向量索引之差的绝对值，所述历史迭代索引为本次迭代之前任意一次迭代的迭代索引；The determining unit 1032 is configured to determine whether the i is 1, and if yes, update the current iteration index to the to-be-updated set, if not, the current iteration index and the extrapolation element in the reconstruction Updating an index of the signal to the to-be-updated set, wherein when the current iteration index is greater than a historical iteration index, the extrapolation element includes an index of the reconstruction signal equal to the current iteration index minus a row An element of a difference value of the index difference, and an element further including an index of the reconstruction signal equal to the historical iteration index plus a sum of the row index differences, when the current iteration index is smaller than the historical iteration index And the extrapolating element includes an element whose index of the reconstructed signal is equal to the sum of the current iteration index plus the row index difference, and further includes an index of the reconstructed signal equal to the historical iteration index minus An element of a difference value of the row index difference, where the row index difference is a row row index corresponding to the current iteration index in a matrix corresponding to the reconstruction signal, and the historical iteration index is in the Built absolute difference of the index of the row vector in the matrix corresponding to the signal corresponding to the present sub-iteration index history Iterative index of any iteration before iteration;

最小二乘单元1033，用于将所述待更新集合中的索引所对应的所述观测矩阵中的列向量构成的部分矩阵和所述协方差矩阵的列向量化向量进行基于最小二乘的重建向量估计，以获得第i次的重建信号；a least squares unit 1033, configured to perform a least squares reconstruction based on a partial matrix composed of column vectors in the observation matrix corresponding to an index in the to-be-updated set and a column vectorization vector of the covariance matrix Vector estimation to obtain the i-th reconstruction signal;

确定单元1034，用于当所述i等于所述L时，将所述第i次的重建信号确定为所述信道的稀疏统计信息的重建信号；a determining unit 1034, configured to determine, when the i is equal to the L, the reconstructed signal of the ith time as a reconstruction signal of the sparse statistical information of the channel;

第二计算子单元1035，用于当所述i小于所述L时，将所述协方差矩阵的列向量化向量减去所述部分矩阵和所述第i次的重建信号的乘积之后的结果作为第i次迭代的残余向量，并将i加1，并触发所述计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积的操作。a second calculating sub-unit 1035, configured to subtract, after the i is smaller than the L, a result of subtracting a product of the partial matrix and the ith reconstruction signal from a column vectorization vector of the covariance matrix As the residual vector of the ith iteration, i is incremented by 1, and the operation of calculating the inner product of each diagonal column vector of the observation matrix and the residual vector of the ith minus one iteration is triggered.

获取单元104，用于使用所述重建信号获取所述信道的最终统计信息。The obtaining unit 104 is configured to obtain final statistical information of the channel by using the reconstruction signal.

本实施例中，获取单元104可以用于根据预先获取的稀疏统计信息与最终统计信息的关系信息获取与所述重建信号对应的最终统计信息作为所述信道的最终统计信息。In this embodiment, the obtaining unit 104 may be configured to obtain, according to the relationship information of the pre-acquired sparse statistical information and the final statistical information, final statistical information corresponding to the reconstructed signal as final statistical information of the channel.

本实施例中，在图9所示的实施例的基础上增加了多种可选的实施方式，且都可以实现提高信道统计信息的准确性。In this embodiment, a plurality of optional implementation manners are added on the basis of the embodiment shown in FIG. 9, and the accuracy of the channel statistics information can be improved.

请参阅图11，图11是本发明实施例提供的另一种接收机的结构示意图，如图11所示，包括：处理器111、网络接口112、存储器113和通信总线114，其中，所述通信总线114用于实现所述处理器111、网络接口112和存储器113之间连接通信，所述处理器111执行所述存储器113中存储的程序用于实现以下方法：Referring to FIG. 11, FIG. 11 is a schematic structural diagram of another receiver according to an embodiment of the present invention. As shown in FIG. 11, the method includes: a processor 111, a network interface 112, a memory 113, and a communication bus 114. The communication bus 114 is configured to implement connection communication between the processor 111, the network interface 112, and the memory 113, and the processor 111 executes a program stored in the memory 113 for implementing the following method:

在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数；Receiving, in F time slots, a received pilot signal including a pilot signal transmitted by a transmitter, wherein the F is an integer greater than or equal to 1;

使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵；Calculating a covariance matrix of the received signal of the receiver on a channel transmitting the pilot signal using the received F received pilot signals;

使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号；Performing signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of sparse statistical information of the channel;

使用所述重建信号获取所述信道的最终统计信息。 The final statistical information of the channel is obtained using the reconstruction signal.

本实施例中，处理器111执行的使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵的程序，可以包括：In this embodiment, the program executed by the processor 111 to calculate the covariance matrix of the received signal of the receiver on the channel for transmitting the pilot signal by using the F received pilot signals may include:

将所述接收到F个接收导频信号进行列向量化，以获取所述F个接收导频信号的列向量化向量，并使用所述F个接收导频信号的列向量化向量计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。Performing column vectorization on the received F received pilot signals to obtain a column vectorization vector of the F received pilot signals, and calculating the column vectorization vector using the F received pilot signals A covariance matrix of received signals on a channel on which the receiver transmits the pilot signal.

本实施例中，处理器111执行的使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号的程序，可以包括：In this embodiment, the program executed by the processor 111 to perform signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of the sparse statistical information of the channel may include:

以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号，其中，所述观测矩阵为预先获取的用于进行信号重建的矩阵，所述L为稀疏度。Performing a column vectorization vector, an observation matrix, and L of the covariance matrix as an input of an orthogonal matching tracking OMP algorithm, and executing the OMP algorithm to obtain a reconstruction signal of sparse statistical information of the channel, where The observation matrix is a pre-acquired matrix for signal reconstruction, and the L is sparsity.

本实施例中，处理器111执行的以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号的程序，可以包括：In this embodiment, the column vectorization vector, the observation matrix, and L of the covariance matrix performed by the processor 111 are used as input of the orthogonal matching tracking OMP algorithm, and the OMP algorithm is executed to obtain the sparseness of the channel. The program for reconstructing signals of statistical information may include:

计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积，并将内积的绝对值最大的对角线列向量的索引作为本次迭代索引，其中，所述对角线列向量是所述观测矩阵中与所述重建信号对应的矩阵中对角线元素的索引标识对应的列向量，所述i为本次迭代标号，所述第i减1次迭代的残余向量为所述协方差矩阵的列向量化向量经过i减1次迭代后的残余；Calculating an inner product of each diagonal line vector of the observation matrix and a residual vector of the ith minus one iteration, and using an index of a diagonal column vector having the largest absolute value of the inner product as the index of the iteration, wherein The diagonal column vector is a column vector corresponding to an index identifier of a diagonal element in a matrix corresponding to the reconstruction signal in the observation matrix, where i is the current iteration label, and the ith is reduced by 1 time The residual vector of the iteration is a residual of the column vectorization vector of the covariance matrix after i is decremented by one iteration;

判断所述i是否为1，若是，则将所述本次迭代索引更新至待更新集合中，若否，则将所述本次迭代索引和外推元素在所述重建信号的索引更新至所述待更新集合中，其中，当所述本次迭代索引大于历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引减去行索引差的差值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引加上所述行索引差的和值的元素，当所述本次迭代索引小于所述历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引加上行索引差的和值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引减去所述行索引差的差值的元素，所述行索引差为所述本次迭代索引在所述重建信号对应的矩阵中对应的行向量索引与所述历史迭代索引在所述重建信号对应的矩阵中对应的行向量索引之差的绝对值，所述历史迭代索引为本次迭代之前任意一次迭代的迭代索引；Determining whether the i is 1, if yes, updating the current iteration index to the to-be-updated set, and if not, updating the index of the current iteration index and the extrapolation element to the reconstructed signal to the In the update set, wherein, when the current iteration index is greater than the historical iteration index, the extrapolation element includes an index of the reconstruction signal equal to the difference between the current iteration index minus the row index difference An element, and an element further including an index of the reconstruction signal equal to the historical iteration index plus a sum of the row index differences, when the current iteration index is smaller than the historical iteration index, the extrapolation An element includes an element at which an index of the reconstructed signal is equal to a sum value of the current iteration index plus a row index difference, and further comprising an index at the reconstructed signal equal to the historical iteration index minus the row An element of the difference value of the difference, wherein the row index difference is that the current row index corresponding to the current iteration index in the matrix corresponding to the reconstructed signal corresponds to the matrix corresponding to the historical iteration index in the reconstructed signal The absolute value of the difference between the row vector indices, which is an iterative index of any iteration before the iteration;

将所述待更新集合中的索引所对应的所述观测矩阵中的列向量构成的部分矩阵和所述协方差矩阵的列向量化向量进行基于最小二乘的重建向量估计，以获得第i次的重建信号；And performing a least squares reconstruction vector estimation on the partial matrix formed by the column vectors in the observation matrix corresponding to the index in the to-be-updated set and the column vectorization vector of the covariance matrix to obtain the i-th time Reconstruction signal;

当所述i等于所述L时，将所述第i次的重建信号确定为所述信道的稀疏统计信息的重建信号；When the i is equal to the L, determining the ith reconstruction signal as a reconstruction signal of the sparse statistical information of the channel;

当所述i小于所述L时，将所述协方差矩阵的列向量化向量减去所述部分矩阵和所述第i次的重建信号的乘积之后的结果作为第i次迭代的残余向量，并将i加1，并触发所述计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积的步骤。When the i is smaller than the L, the result of subtracting the product of the partial matrix and the ith reconstruction signal from the column vectorization vector of the covariance matrix is used as the residual vector of the ith iteration. And i is incremented by one, and the step of calculating the inner product of each diagonal column vector of the observation matrix and the residual vector of the ith minus one iteration is triggered.

本实施例中，处理器111执行的使用所述重建信号获取所述信道的最终统计信息的程序，可以包括：In this embodiment, the program executed by the processor 111 to obtain the final statistical information of the channel by using the reconstruction signal may include:

根据预先获取的稀疏统计信息与最终统计信息的关系信息获取与所述重建信号对应的最终统计信息作为所述信道的最终统计信息。And obtaining, according to the relationship information of the pre-acquired sparse statistical information and the final statistical information, final statistical information corresponding to the reconstructed signal as final statistical information of the channel.

另外，本发明实施例还提供一种存储一个或多个程序的计算机可读存储介质，所述一个或多个程序包括指令，所述指令当被包括屏幕和多个应用程序的接收机执行时使该接收机执行本发明实施例提供的任意一种实现方式所述的方法。 In addition, embodiments of the present invention also provide a computer readable storage medium storing one or more programs, the one or more programs including instructions when executed by a receiver including a screen and a plurality of applications The receiver is configured to perform the method described in any one of the implementation manners provided by the embodiments of the present invention.

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存取存储器(Random Access Memory，简称RAM)等。One of ordinary skill in the art can understand that all or part of the process of implementing the foregoing embodiments can be completed by a computer program to instruct related hardware, and the program can be stored in a computer readable storage medium. When executed, the flow of an embodiment of the methods as described above may be included. The storage medium may be a magnetic disk, an optical disk, a read-only memory (ROM), or a random access memory (RAM).

以上所揭露的仅为本发明较佳实施例而已，当然不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。 The above is only the preferred embodiment of the present invention, and the scope of the present invention is not limited thereto, and thus equivalent changes made in the claims of the present invention are still within the scope of the present invention.

Claims

一种信道统计信息获取方法，其特征在于，包括：A method for acquiring channel statistics information, comprising:

接收机在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数；The receiver receives, in F time slots, a received pilot signal including a pilot signal transmitted by a transmitter, wherein the F is an integer greater than or equal to 1;

所述接收机使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵；The receiver calculates a covariance matrix of the received signal of the receiver on a channel transmitting the pilot signal by using the F received pilot signals;

所述接收机使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号；The receiver performs signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of sparse statistical information of the channel;

所述接收机使用所述重建信号获取所述信道的最终统计信息。The receiver uses the reconstructed signal to obtain final statistical information for the channel.
如权利要求1所述的方法，其特征在于，所述接收机使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵，包括：The method of claim 1 wherein said receiver calculates a covariance matrix of received signals of said receiver on a channel transmitting said pilot signal using said received received pilot signals, including :

所述接收机将所述接收到F个接收导频信号进行列向量化，以获取所述F个接收导频信号的列向量化向量，并使用所述F个接收导频信号的列向量化向量计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。The receiver performs column vectorization on the received F received pilot signals to obtain a column vectorization vector of the F received pilot signals, and uses the column vectorization of the F received pilot signals. A vector calculates a covariance matrix of the received signal of the receiver on a channel on which the pilot signal is transmitted.
如权利要求1或2所述的方法，其特征在于，所述接收机使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号，包括：The method according to claim 1 or 2, wherein the receiver performs signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of the sparse statistical information of the channel, including:

所述接收机以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号，其中，所述观测矩阵为预先获取的用于进行信号重建的矩阵，所述L为稀疏度。The receiver uses the column vectorization vector, the observation matrix, and L of the covariance matrix as an input of the orthogonal matching tracking OMP algorithm, and executes the OMP algorithm to obtain a reconstruction signal of the sparse statistical information of the channel, The observation matrix is a pre-acquired matrix for performing signal reconstruction, and the L is a sparsity degree.
如权利要求3所述的方法，其特征在于，所述接收机以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号，包括：The method according to claim 3, wherein said receiver performs a column vectorization vector, an observation matrix, and L of said covariance matrix as an input of an orthogonal matching tracking OMP algorithm, and executes said OMP algorithm, Reconstructing signals for obtaining sparse statistical information of the channel, including:

所述接收机计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积，并将内积的绝对值最大的对角线列向量的索引作为本次迭代索引，其中，所述对角线列向量是所述观测矩阵中与所述重建信号对应的矩阵中对角线元素的索引标识对应的列向量，所述i为本次迭代标号，所述第i减1 次迭代的残余向量为所述协方差矩阵的列向量化向量经过i减1次迭代后的残余；The receiver calculates an inner product of each diagonal line vector of the observation matrix and a residual vector of the ith minus one iteration, and uses an index of a diagonal column vector having the largest absolute value of the inner product as the current iteration An index, wherein the diagonal column vector is a column vector corresponding to an index identifier of a diagonal element in a matrix corresponding to the reconstruction signal in the observation matrix, where the i is an iteration label, i minus 1 The residual vector of the second iteration is a residual of the column vectorization vector of the covariance matrix after i is decremented by one iteration;

所述接收机判断所述i是否为1，若是，则将所述本次迭代索引更新至待更新集合中，若否，则将所述本次迭代索引和外推元素在所述重建信号的索引更新至所述待更新集合中，其中，当所述本次迭代索引大于历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引减去行索引差的差值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引加上所述行索引差的和值的元素，当所述本次迭代索引小于所述历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引加上行索引差的和值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引减去所述行索引差的差值的元素，所述行索引差为所述本次迭代索引在所述重建信号对应的矩阵中对应的行向量索引与所述历史迭代索引在所述重建信号对应的矩阵中对应的行向量索引之差的绝对值，所述历史迭代索引为本次迭代之前任意一次迭代的迭代索引；Determining, by the receiver, whether the i is 1, if yes, updating the current iteration index to the to-be-updated set, and if not, the current iteration index and the extrapolating element in the reconstructed signal Indexing is updated to the to-be-updated set, wherein when the current iteration index is greater than a historical iteration index, the extrapolation element includes an index of the reconstruction signal equal to the current iteration index minus a row index difference An element of the difference value, and an element further including an index of the reconstruction signal equal to the historical iteration index plus a sum of the row index differences, when the current iteration index is smaller than the historical iteration index, The extrapolation element includes an element whose index of the reconstruction signal is equal to a sum value of the current iteration index plus a row index difference, and further includes an index at the reconstruction signal equal to the historical iteration index minus the An element of a difference value of a row index difference, wherein the row index difference is a row row index corresponding to the current iteration index in a matrix corresponding to the reconstruction signal, and the historical iteration index is in the reconstruction letter The absolute difference matrix corresponding to the index of the corresponding row vector, prior to the present iteration history any iteration index of the next iteration iteration index;

所述接收机将所述待更新集合中的索引所对应的所述观测矩阵中的列向量构成的部分矩阵和所述协方差矩阵的列向量化向量进行基于最小二乘的重建向量估计，以获得第i次的重建信号；The receiver performs a least squares reconstruction vector estimation based on a partial matrix formed by column vectors in the observation matrix corresponding to an index in the to-be-updated set and a column vectorization vector of the covariance matrix. Obtaining the i-th reconstruction signal;

当所述i等于所述L时，所述接收机将所述第i次的重建信号确定为所述信道的稀疏统计信息的重建信号；When the i is equal to the L, the receiver determines the i-th reconstruction signal as a reconstruction signal of the sparse statistical information of the channel;

当所述i小于所述L时，所述接收机将所述协方差矩阵的列向量化向量减去所述部分矩阵和所述第i次的重建信号的乘积之后的结果作为第i次迭代的残余向量，并将i加1，并触发所述计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积的步骤。When the i is smaller than the L, the receiver subtracts the product of the partial matrix and the ith reconstruction signal from the column vectorization vector of the covariance matrix as the ith iteration a residual vector, and i is incremented by one, and triggers the step of calculating the inner product of each diagonal column vector of the observation matrix and the residual vector of the ith minus one iteration.
如权利要求1-4中任一项所述的方法，其特征在于，所述接收机使用所述重建信号获取所述信道的最终统计信息，包括：The method according to any one of claims 1 to 4, wherein the receiver uses the reconstructed signal to obtain final statistical information of the channel, including:

所述接收机根据预先获取的稀疏统计信息与最终统计信息的关系信息获取与所述重建信号对应的最终统计信息作为所述信道的最终统计信息。The receiver acquires final statistical information corresponding to the reconstructed signal as final statistical information of the channel according to the relationship information between the pre-acquired sparse statistical information and the final statistical information.
一种接收机，其特征在于，包括：接收单元、计算单元、重建单元和获取单元，其中： A receiver, comprising: a receiving unit, a calculating unit, a reconstructing unit, and an acquiring unit, wherein:

所述接收单元，用于在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数；The receiving unit is configured to receive, in F time slots, a received pilot signal including a pilot signal sent by a transmitter, where the F is an integer greater than or equal to 1;

所述计算单元，用于使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵；The calculating unit is configured to calculate, by using the F received pilot signals, a covariance matrix of the received signal of the receiver on a channel for transmitting the pilot signal;

所述重建单元，用于使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号；The reconstructing unit is configured to perform signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of sparse statistical information of the channel;

所述获取单元，用于使用所述重建信号获取所述信道的最终统计信息。The acquiring unit is configured to obtain final statistical information of the channel by using the reconstruction signal.
如权利要求6所述的接收机，其特征在于，所述计算单元用于将所述接收到F个接收导频信号进行列向量化，以获取所述F个接收导频信号的列向量化向量，并使用所述F个接收导频信号的列向量化向量计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。The receiver according to claim 6, wherein said calculating unit is configured to perform column vectorization on said received F received pilot signals to obtain column vectorization of said F received pilot signals a vector, and using a column vectorization vector of the F received pilot signals to calculate a covariance matrix of the received signal of the receiver on a channel on which the pilot signal is transmitted.
如权利要求6或7所述的接收机，其特征在于，所述重建单元用于以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号，其中，所述观测矩阵为预先获取的用于进行信号重建的矩阵，所述L为稀疏度。The receiver according to claim 6 or 7, wherein said reconstruction unit is configured to perform a column vectorization vector, an observation matrix, and an L of the covariance matrix as an input of an orthogonal matching tracking OMP algorithm, and execute And the OMP algorithm is configured to obtain a reconstruction signal of the sparse statistical information of the channel, where the observation matrix is a pre-acquired matrix for performing signal reconstruction, and the L is a sparsity degree.
如权利要求8所述的接收机，其特征在于，所述重建单元包括：The receiver of claim 8 wherein said reconstruction unit comprises:

第一计算子单元，用于计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积，并将内积的绝对值最大的对角线列向量的索引作为本次迭代索引，其中，所述对角线列向量是所述观测矩阵中与所述重建信号对应的矩阵中对角线元素的索引标识对应的列向量，所述i为本次迭代标号，所述第i减1次迭代的残余向量为所述协方差矩阵的列向量化向量经过i减1次迭代后的残余；a first calculation subunit, configured to calculate an inner product of each diagonal line vector of the observation matrix and an inner product of the residual vector of the ith minus one iteration, and to index the diagonal column vector with the largest absolute value of the inner product As the current iteration index, the diagonal column vector is a column vector corresponding to an index identifier of a diagonal element in a matrix corresponding to the reconstruction signal in the observation matrix, and the i is an iteration label And the residual vector of the ith minus one iteration is a residual of the column vectorization vector of the covariance matrix after i is decremented by one iteration;

判断单元，用于判断所述i是否为1，若是，则将所述本次迭代索引更新至待更新集合中，若否，则将所述本次迭代索引和外推元素在所述重建信号的索引更新至所述待更新集合中，其中，当所述本次迭代索引大于历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引减去行索引差的差值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引加上所述行索引差的和值的元素，当所述本次迭代索引小于所述历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引加上行索引差的和值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引减去所述行索引差的差值的元素，所述行索引差为所述本次迭代索引在所述重建信号对应的矩阵中对应的行向量索引与所述历史迭代索引在所述重建信号对应的矩阵中对应的行向量索引之差的绝对值，所述历史迭代索引为本次迭代之前任意一次迭代的迭代索引；a determining unit, configured to determine whether the i is 1, if yes, updating the current iteration index to the to-be-updated set, and if not, the current iteration index and the extrapolating element in the reconstructing signal The index is updated to the to-be-updated set, wherein when the current iteration index is greater than the historical iteration index, the extrapolation element includes an index of the reconstruction signal equal to the current iteration index minus a row index An element of the difference difference, and an element further included in the index of the reconstruction signal equal to the historical iteration index plus the sum of the row index differences, when the current iteration index is smaller than the historical iteration index The extrapolation element includes an index at the reconstruction signal equal to the current iteration index plus a line An element of the sum of the sum values, and an element further including an index of the reconstructed signal equal to a difference between the historical iteration index minus the row index difference, the row index difference being the current iteration index An absolute value of a difference between a corresponding row vector index in the matrix corresponding to the reconstructed signal and a corresponding row vector index of the historical iteration index in the matrix corresponding to the reconstructed signal, the historical iterative index being arbitrary before the iteration Iterative index of one iteration;

最小二乘单元，用于将所述待更新集合中的索引所对应的所述观测矩阵中的列向量构成的部分矩阵和所述协方差矩阵的列向量化向量进行基于最小二乘的重建向量估计，以获得第i次的重建信号；a least squares unit, configured to perform a least squares reconstruction vector based on a partial matrix composed of column vectors in the observation matrix corresponding to an index in the to-be-updated set and a column vectorization vector of the covariance matrix Estimating to obtain the ith reconstruction signal;

确定单元，用于当所述i等于所述L时，将所述第i次的重建信号确定为所述信道的稀疏统计信息的重建信号；a determining unit, configured to determine the i-th reconstruction signal as a reconstruction signal of the sparse statistical information of the channel when the i is equal to the L;

第二计算子单元，用于当所述i小于所述L时，将所述协方差矩阵的列向量化向量减去所述部分矩阵和所述第i次的重建信号的乘积之后的结果作为第i次迭代的残余向量，并将i加1，并触发所述计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积的操作。a second calculating subunit, configured to: when the i is smaller than the L, subtract a result of subtracting the product of the partial matrix and the ith reconstruction signal from a column vectorization vector of the covariance matrix The residual vector of the i-th iteration, and i is incremented by one, and triggers the operation of calculating the inner product of each diagonal column vector of the observation matrix and the residual vector of the ith minus one iteration.
如权利要求6-9中任一项所述的接收机，其特征在于，所述获取单元用于根据预先获取的稀疏统计信息与最终统计信息的关系信息获取与所述重建信号对应的最终统计信息作为所述信道的最终统计信息。The receiver according to any one of claims 6 to 9, wherein the obtaining unit is configured to obtain a final statistic corresponding to the reconstructed signal according to the relationship information between the pre-acquired sparse statistical information and the final statistical information. Information is used as the final statistical information of the channel.
一种接收机，其特征在于，包括：处理器、网络接口、存储器和通信总线，其中，所述通信总线用于实现所述处理器、网络接口和存储器之间连接通信，所述处理器执行所述存储器中存储的程序用于实现以下方法：A receiver, comprising: a processor, a network interface, a memory, and a communication bus, wherein the communication bus is configured to implement connection communication between the processor, a network interface, and a memory, the processor executing The program stored in the memory is used to implement the following methods:

在F个时隙中接收包括发射机发送的导频信号的接收导频信号，其中，所述F为大于或者等于1的整数；Receiving, in F time slots, a received pilot signal including a pilot signal transmitted by a transmitter, wherein the F is an integer greater than or equal to 1;

使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵；Calculating a covariance matrix of the received signal of the receiver on a channel transmitting the pilot signal using the received F received pilot signals;

使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号；Performing signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of sparse statistical information of the channel;

使用所述重建信号获取所述信道的最终统计信息。The final statistical information of the channel is obtained using the reconstruction signal.
如权利要求11所述的接收机，其特征在于，所述处理器执行的使用接收到F个接收导频信号计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵的程序，包括：The receiver according to claim 11, wherein said processor performs the use of receiving F received pilot signals to calculate a connection of said receiver on a channel for transmitting said pilot signal The procedure for receiving the covariance matrix of the signal, including:

将所述接收到F个接收导频信号进行列向量化，以获取所述F个接收导频信号的列向量化向量，并使用所述F个接收导频信号的列向量化向量计算所述接收机在传输所述导频信号的信道上的接收信号的协方差矩阵。Performing column vectorization on the received F received pilot signals to obtain a column vectorization vector of the F received pilot signals, and calculating the column vectorization vector using the F received pilot signals A covariance matrix of received signals on a channel on which the receiver transmits the pilot signal.
如权利要求11或12所述的接收机，其特征在于，所述处理器执行的使用所述协方差矩阵进行基于压缩协方差感知的信号重建，以获得所述信道的稀疏统计信息的重建信号的程序，包括：The receiver according to claim 11 or 12, wherein the processor performs a signal reconstruction based on compression covariance sensing using the covariance matrix to obtain a reconstruction signal of sparse statistical information of the channel. Programs, including:

以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号，其中，所述观测矩阵为预先获取的用于进行信号重建的矩阵，所述L为稀疏度。Performing a column vectorization vector, an observation matrix, and L of the covariance matrix as an input of an orthogonal matching tracking OMP algorithm, and executing the OMP algorithm to obtain a reconstruction signal of sparse statistical information of the channel, where The observation matrix is a pre-acquired matrix for signal reconstruction, and the L is sparsity.
如权利要求13所述的接收机，其特征在于，所述处理器执行的以所述协方差矩阵的列向量化向量、观测矩阵和L作为正交匹配追踪OMP算法的输入，并执行所述OMP算法，以获得所述信道的稀疏统计信息的重建信号的程序，包括：The receiver according to claim 13, wherein said processor performs a column vectorization vector, an observation matrix, and L of said covariance matrix as an input of an orthogonal matching tracking OMP algorithm, and performs said The OMP algorithm, the program for obtaining a reconstructed signal of the sparse statistical information of the channel, comprising:

计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积，并将内积的绝对值最大的对角线列向量的索引作为本次迭代索引，其中，所述对角线列向量是所述观测矩阵中与所述重建信号对应的矩阵中对角线元素的索引标识对应的列向量，所述i为本次迭代标号，所述第i减1次迭代的残余向量为所述协方差矩阵的列向量化向量经过i减1次迭代后的残余；Calculating an inner product of each diagonal line vector of the observation matrix and a residual vector of the ith minus one iteration, and using an index of a diagonal column vector having the largest absolute value of the inner product as the index of the iteration, wherein The diagonal column vector is a column vector corresponding to an index identifier of a diagonal element in a matrix corresponding to the reconstruction signal in the observation matrix, where i is the current iteration label, and the ith is reduced by 1 time The residual vector of the iteration is a residual of the column vectorization vector of the covariance matrix after i is decremented by one iteration;

判断所述i是否为1，若是，则将所述本次迭代索引更新至待更新集合中，若否，则将所述本次迭代索引和外推元素在所述重建信号的索引更新至所述待更新集合中，其中，当所述本次迭代索引大于历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引减去行索引差的差值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引加上所述行索引差的和值的元素，当所述本次迭代索引小于所述历史迭代索引时，所述外推元素包括在所述重建信号的索引等于所述本次迭代索引加上行索引差的和值的元素，以及还包括在所述重建信号的索引等于所述历史迭代索引减去所述行索引差的差值的元素，所述行索引差为所述本次迭代索引在所述重建信号对应的矩阵中对应的行向量索引与所述历史迭代索引在所述重建信号对应的矩阵中对应的行向量索引之差的绝对值，所述历史迭代索引为本次迭代之前任意一次迭代的迭代索引；Determining whether the i is 1, if yes, updating the current iteration index to the to-be-updated set, and if not, updating the index of the current iteration index and the extrapolation element to the reconstructed signal to the In the update set, wherein, when the current iteration index is greater than the historical iteration index, the extrapolation element includes an index of the reconstruction signal equal to the difference between the current iteration index minus the row index difference An element, and an element further including an index of the reconstruction signal equal to the historical iteration index plus a sum of the row index differences, when the current iteration index is smaller than the historical iteration index, the extrapolation An element includes an element having an index of the reconstruction signal equal to a sum value of the current iteration index plus a row index difference, and further comprising an index at the reconstruction signal equal to the historical iteration index minus the row index difference An element of the difference, the row index difference being the corresponding one of the current iteration index corresponding to the reconstruction signal An absolute value of a difference between a corresponding row vector index in the matrix and a corresponding row vector index of the historical iteration index in the matrix corresponding to the reconstructed signal, the historical iteration index being an iterative index of any one iteration before the iteration;

将所述待更新集合中的索引所对应的所述观测矩阵中的列向量构成的部分矩阵和所述协方差矩阵的列向量化向量进行基于最小二乘的重建向量估计，以获得第i次的重建信号；And performing a least squares reconstruction vector estimation on the partial matrix formed by the column vectors in the observation matrix corresponding to the index in the to-be-updated set and the column vectorization vector of the covariance matrix to obtain the i-th time Reconstruction signal;

当所述i等于所述L时，将所述第i次的重建信号确定为所述信道的稀疏统计信息的重建信号；When the i is equal to the L, determining the ith reconstruction signal as a reconstruction signal of the sparse statistical information of the channel;

当所述i小于所述L时，将所述协方差矩阵的列向量化向量减去所述部分矩阵和所述第i次的重建信号的乘积之后的结果作为第i次迭代的残余向量，并将i加1，并触发所述计算所述观测矩阵的各个对角线列向量与第i减1次迭代的残余向量的内积的步骤。When the i is smaller than the L, the result of subtracting the product of the partial matrix and the ith reconstruction signal from the column vectorization vector of the covariance matrix is used as the residual vector of the ith iteration. And i is incremented by one, and the step of calculating the inner product of each diagonal column vector of the observation matrix and the residual vector of the ith minus one iteration is triggered.
如权利要求11-14中任一项所述的接收机，其特征在于，所述处理器执行的使用所述重建信号获取所述信道的最终统计信息的程序，包括：The receiver according to any one of claims 11 to 14, wherein the program executed by the processor to obtain final statistical information of the channel using the reconstruction signal comprises:

根据预先获取的稀疏统计信息与最终统计信息的关系信息获取与所述重建信号对应的最终统计信息作为所述信道的最终统计信息。 And obtaining, according to the relationship information of the pre-acquired sparse statistical information and the final statistical information, final statistical information corresponding to the reconstructed signal as final statistical information of the channel.