CN110099016A - A kind of sparse front channel estimation methods of millimeter wave based on deep learning network - Google Patents

A kind of sparse front channel estimation methods of millimeter wave based on deep learning network Download PDF

Info

Publication number
CN110099016A
CN110099016A CN201910397076.5A CN201910397076A CN110099016A CN 110099016 A CN110099016 A CN 110099016A CN 201910397076 A CN201910397076 A CN 201910397076A CN 110099016 A CN110099016 A CN 110099016A
Authority
CN
China
Prior art keywords
matrix
channel
sparse
dimension
network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910397076.5A
Other languages
Chinese (zh)
Other versions
CN110099016B (en
Inventor
许威
张雯惠
徐锦丹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201910397076.5A priority Critical patent/CN110099016B/en
Publication of CN110099016A publication Critical patent/CN110099016A/en
Application granted granted Critical
Publication of CN110099016B publication Critical patent/CN110099016B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • H04L25/0254Channel estimation channel estimation algorithms using neural network algorithms
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2626Arrangements specific to the transmitter only
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Artificial Intelligence (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Power Engineering (AREA)
  • Radio Transmission System (AREA)

Abstract

The invention discloses a kind of sparse front channel estimation methods of millimeter wave based on deep learning network, and using millimeter wave channel sparse characteristic as prior information, the full connection deep neural network of training design is used for millimeter wave surface communications channel estimation.First using full connection phase shifter network, isotropic analog transceiver is designed by configuring each phase shifter uniform phase distribution;Then using the channel sparse information of acquisition and the optimal digital estimator of design as the training data of full connection deep learning network.For the condition of sparse channel under each signal-to-noise ratio, the sparse information of channel is inputted into network, obtains corresponding digital estimator, and then obtain channel estimation results.The condition of sparse channel estimator that the present invention provides can reduce low precision analog-digital converter nonlinear quantization bring error, and use deep learning network implementations, to reduce channel estimation complexity, channel estimation methods that inventive can be optimal on approximation theory.

Description

Millimeter wave sparse array surface channel estimation method based on deep learning network
Technical Field
The invention relates to the field of communication, in particular to a millimeter wave sparse array surface channel estimation method based on a deep learning network.
Background
In recent years, communication technology has been developed in a breakthrough manner and is becoming mature, and the mobile communication industry has been developed rapidly on a global scale. The Multiple Input Multiple Output (MIMO) technology is one of the key technologies for the development of communication technology, and improves the data transmission rate of the system. A plurality of antennas are arranged at a transmitting end and a receiving end of the system, and diversity is formed by the antennas at the transmitting end and the receiving end, so that the stability of the system can be improved. Meanwhile, the number of independent channels between the antennas at the transmitting end and the receiving end is greatly increased, so that the data volume transmitted by the system in unit time is also improved, and the spectrum utilization efficiency of the system is further improved.
Compared with the MIMO technology, the Massive MIMO technology has the advantages that the transmission rate, the energy efficiency, the transmission reliability and the like are greatly improved, the millimeter wave Massive MIMO technology greatly reduces the configuration difficulty of a Massive antenna array surface, and the Massive MIMO technology solves the problem that millimeter wave signals are high in loss and easy to block. In order to reduce power consumption and hardware complexity in the millimeter wave massive MIMO system, a digital-analog hybrid architecture of a small number of radio frequency links may be used.
In order to perform high performance transmission, channel estimation needs to be performed first. Channel estimation of a conventional large-scale multi-antenna technology system is a big challenge in itself, and the use of a low-precision analog-to-digital converter in a hybrid architecture makes channel estimation more difficult while reducing cost and power consumption. Meanwhile, it is also a challenge how to reduce the estimation complexity and pilot overhead while utilizing the sparsity of the millimeter wave channel, and obtain high-precision channel estimation.
Disclosure of Invention
In order to solve the problems, the invention provides a sparse channel estimation method which is used for a low-precision analog-to-digital converter and a broadband mmWave large-scale multi-antenna technology system with a hybrid framework. And training the designed full-connection deep neural network by taking the sparse characteristic of the millimeter wave channel as prior information and taking the selection matrix of the sparse channel and a digital estimator corresponding to the selection matrix as input to obtain the deep neural network suitable for different signal-to-noise ratios, and the deep neural network is used for estimating the communication channel of the millimeter wave array. Firstly, a fully-connected phase shifter network is adopted, and an isotropic analog transceiver is designed by configuring the phase uniform distribution of each phase shifter; and then, the obtained channel sparse information is used as priori knowledge, an optimal digital estimator is designed, and the channel sparse information and the optimal digital estimator are used as training data of the fully-connected deep learning network. For the sparse channel under each signal-to-noise ratio, inputting the sparse information of the channel into the network, and obtaining a corresponding digital estimator so as to obtain a channel estimation result. The sparse channel estimator provided by the invention can reduce the error caused by the nonlinear quantization of the low-precision analog-to-digital converter, and is realized by using a deep learning network, so that the channel estimation complexity is reduced. In a mixed-frame multi-antenna system adopting a low-precision analog-to-digital converter, the performance of the method can approach the theoretically optimal channel estimation method.
In order to achieve the purpose, the invention provides the following technical scheme:
a millimeter wave sparse array surface channel estimation method based on a deep learning network comprises the following steps:
a transmitting end sends a pilot signal, a simulation precoder consisting of fully-connected phase shifters reaches a receiver through a channel, and a receiving end obtains channel estimation by using a simulation estimator, a low-precision analog-to-digital converter quantizer and a deep learning network;
the method is characterized in that: the channel estimation method based on the hardware architecture comprises the following steps:
the method comprises the following steps: base station design simulation pre-coding and sending pilot signal, simulation pre-coding FAmDesigned according to the following formula:
wherein, the [ alpha ], [ beta ]]ijElements representing the ith row and the jth column of the matrix; fAmWith a representation dimension of Nt×NRFtAnalog precoding matrix of, NtRepresenting the number of transmitting antennas, NRFtRepresenting the number of radio frequency links of a transmitting end;representing the phase of the ith row and jth column element of the analog precoding matrix;
step two: receiving terminal optimization deviceAnalog estimator WAm,WAmCalculated according to the following formula:
wherein, WAmWith a representation dimension of Nr×NRFrAnalog estimation matrix of, NrRepresenting the number of receiving antennas, NRFrRepresenting the number of radio frequency links of a receiving end;representing the phase of the ith row and jth column element of the simulation estimation matrix;
step three: the receiving end obtains sparse characteristic information of the channel, namely the position of a non-zero element in a sparse channel matrix, through a compressed sensing technology, such as a classical Orthogonal Matching Pursuit (OMP) technology, and the sparse characteristic information is used as prior information of subsequent channel estimation;
step four: all possible selection matrixes P corresponding to sparse channel characteristic information by the receiving endv[k]Is stored and a selection matrix P is selected according to eachv[k]Designing corresponding digital estimation matrix
Step five: selecting the matrix Pv[k]And a digital estimation matrix corresponding theretoAnd as a training set of the fully-connected neural network, obtaining the fully-connected feedforward network suitable for each signal-to-noise ratio through neural network training. Within the coherent time, the receiving end selects the matrix P according to the corresponding selection matrix of the sparse channelv[k]Obtaining a digital estimation matrix of the network output using the network structureDetecting the received signal by using the data estimation matrix to obtain a channel estimation value
Further, the matrix P is selected in the fourth stepv[k]Expressed by the following formula:
wherein e isπ(i)(π(i)∈{1,2,...,NrNt}) represents a dimension of NrNtThe pi (i) -th element of x 1 is a vector of 1 and the remaining elements are 0. N is a radical ofvRepresenting a dimension of N on the k-th subcarrierrNtVector component h of x 1 channel vector projected onto angular domainv[k]Number of medium non-zero elements.
For non-zero channel elements in the channel matrix, the selection matrix P of the channelv[k]Has N possible forms, for all possible sets Pv1[k],Pv2[k],...,PvN[k]Each Pvi[k]The possibility of occurrence of (i ═ 1.., N) is allNamely, it is
Wherein N is the total number of the possible numbers,c represents a combination number formula. N is a radical ofvRepresenting a dimension of N on the k-th subcarrierrNtVector component h of x 1 channel vector projected onto angular domainv[k]Number of medium non-zero elements.
Further, projecting the channel to the virtual angle domain to obtain a channel vector component hv[k],hv[k]Calculated according to the following formula:
wherein A istWith a representation dimension of Nt×NtA transmit dictionary matrix composed of transmit front response vectors,representation matrix AtConjugation of (1);represents the kronecker product; a. therWith a representation dimension of Nr×NrA receiving dictionary matrix composed of the receiving array surface response vectors; h [ k ]]Denotes the dimension N on the k sub-carrierr×NtVec (H [ k ]) of]) Representation matrix H [ k ]]Vectorization of (2).
Wherein A istExpressed by the following formula:
wherein,(p∈{1,2,...,Nt}) represents a dimension of NtX 1 of the transmitted wavefront response vector, whereWherein N istP × Q, P denotes the number of antennas on the horizontal axis of the transmitting antenna array, and Q denotes the number of antennas on the vertical axis of the transmitting antenna array;
Arexpressed by the following formula:
wherein,(q∈{1,2,...,Nrdenotes the dimension NrX 1, where,wherein N isrI × J, I denotes the number of antennas on the horizontal axis of the reception antenna array, and J denotes the number of antennas on the vertical axis of the reception antenna array;
further, the optimal number estimation matrix in the fourth stepUsing the channel estimation minimum mean square error criterion,it can be calculated as follows:
wherein K belongs to {1, 2.,. K } represents the kth subcarrier, and K represents the total number of subcarriers;representing N on the k sub-carrierRFr×NvM denotes the number of channel uses within the coherence time, i.e., the number of channel estimates ηbRepresents a distortion factor related to the quantization bit number b of an analog-to-digital converter (ADC); omega k]Denotes the dimension on the k sub-carrier as MNRFr×NvThe measurement matrix of (2); omegaH[k]Represents the matrix omega k]The conjugate transpose of (1);represents the large scale fading coefficient of the channel;with a representation dimension of Nv×NvThe identity matrix of (1);representing the variance of each element of the equivalent noise vector,represents an Additive White Gaussian Noise (AWGN) variance; p denotes a transmit pilot power.
Further, the channel estimation times M are expressed by the following formula:
indicating rounding up.
The measurement matrix Ω [ k ] is represented by the following formula:
Ω[k]=Φ[k]ΨPv[k],
wherein the dimension is NRFr×NrNtPilot correlation matrix ofsm[k](M ∈ {1, 2.., M }) denotes that the dimension at the time of the mth training is NRFrX 1 transmit pilot vector.
Where Ψ represents the dimension NrNt×NrNtThe spatial transformation matrix, Ψ, is represented by the following equation:
wherein A istWith a representation dimension of Nt×NtA transmit dictionary matrix composed of transmit array surface response vectors, ArWith a representation dimension of Nr×NrA receive dictionary matrix composed of receive front response vectors.
Wherein, Pv[k]With a representation dimension of NrNt×NvThe selection matrix of (2).
Further, in step fiveThe method is obtained through a full-connection deep learning network with forward feedback according to the following steps:
(1) for a known number of non-zero elements NvOf the sparse channel, the channel matrix H [ k ]]Denotes the dimension N on the k sub-carrierr×NtUnder each signal-to-noise ratio, the receiving end calculates all possible P according to the sparse characteristic informationv[k]Formed covariance matrix C k]=ΩH[k]Ω[k]Wherein Ω [ k ]]=Φ[k]ΨPv[k]. At the same time, according to the method as claimed in claim 4, each P is calculatedv[k]Formed covariance matrix C k]Corresponding toDesigning a new numerical estimation matrixAnd will all possible C k]And correspondingStored in a database.
(2) Data are extracted from the database and divided into two groups, namely training data and testing data. Carrying out complex value splitting operation on the training data, and carrying out training set C [ k ]]And W'D[k]Split into a real matrix CR[k]、And an imaginary matrix CI[k]、Two parts.
(3) Then C is addedR[k]、And CI[k]、Performing matrix vectorization operation to obtain dimension Nv 2X 1 column vector cR[k]、And cI[k]、C is toR[k]、As input and training targets of the real part deep learning network, cI[k]、As input to the imaginary deep learning network and as a training target.
(4) Two deep learning fully-connected networks are constructed, and the network structures are the same. All are two layers of forward feedback fully-connected neural networks, the number of neurons in the first layer is N, whereinMeanwhile, a bias (bias) connection is set in the first layer, and the transfer function of the first layer is set as a softmax function; the first layer output is connected to the second layer neurons, the number of neurons is N, and the bias (bias) connection is also set at the second layer.
The softmax function is defined as follows:
wherein, the output of each layer of fully-connected network is as follows:
wherein W and b represent parameters of a fully-connected neural network, yiAnd biThe i-th element, x, representing y and bjDenotes the j-th element of x, Wi,jRepresents an element with position (i, j) in W.
(5) The fully-connected neural network is trained, and the deep learning network is tested by taking complex value splitting and matrix vectorization of test data as input, so that a stable network structure under each signal-to-noise ratio is obtained.
(6) In the coherent time, for the sparse channel under each signal-to-noise ratio, a channel covariance matrix C [ k ] containing sparse information of the channel]Input network to obtain corresponding digital estimator
Further, the channel estimation value in step fiveObtained according to the following formula:
wherein (C)HA conjugate transpose operation representing a matrix; ()-1Representing the inversion operation of the matrix.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1) for the broadband channel estimation, OFDM modulation is adopted in the invention, so that frequency domain channel estimation is carried out on each narrowband subcarrier, the invention takes the millimeter wave channel sparse characteristic as prior information, and takes a selection matrix of a sparse channel and a digital estimator corresponding to the selection matrix as input to train a designed full-connection deep neural network, so as to obtain the deep neural network suitable for different signal-to-noise ratios and used for millimeter wave array surface communication channel estimation; 2) the scheme adopts a fully-connected phase shifter network, and designs an isotropic analog transceiver by configuring the phase uniform distribution of each phase shifter; and designing an optimal digital estimator by taking the obtained channel sparse characteristic information as prior knowledge, and taking the optimal digital estimator and the channel sparse characteristic information as training data of the fully-connected deep learning network. For the sparse channel under each signal-to-noise ratio, inputting the sparse information of the channel into a network, and obtaining a corresponding digital estimator so as to obtain a channel estimation result; 3) the sparse channel estimator provided by the invention can reduce the error caused by the nonlinear quantization of the low-precision analog-to-digital converter, and is realized by using a deep learning network, so that the channel estimation complexity is reduced. In a mixed-frame multi-antenna system adopting a low-precision analog-to-digital converter, the performance of the method can approach the theoretically optimal channel estimation method.
Drawings
FIG. 1 is a block diagram of the system of the present invention;
fig. 2 is a graph of Normalized Mean Square Error (NMSE) of channel estimation as a function of signal-to-noise ratio (SNR) for an 8 x 4 millimeter wave antenna array, quantized by a 4-bit analog-to-digital converter at the receiving end;
Detailed Description
The technical solutions provided by the present invention will be described in detail below with reference to specific examples, and it should be understood that the following specific embodiments are only illustrative of the present invention and are not intended to limit the scope of the present invention.
The invention provides a millimeter wave sparse array surface channel estimation method based on a deep learning network, which is characterized in that millimeter wave channel sparse characteristics are used as prior information, a selection matrix of a sparse channel and a corresponding digital estimator are used as input to train a designed full-connection deep neural network, and the deep neural network suitable for different signal-to-noise ratios is obtained and used for millimeter wave array surface communication channel estimation.
Firstly, a fully-connected phase shifter network is adopted, and an isotropic analog transceiver is designed by configuring the phase uniform distribution of each phase shifter; and then, the obtained channel sparse information is used as priori knowledge, an optimal digital estimator is designed, and the channel sparse information and the optimal digital estimator are used as training data of the fully-connected deep learning network. For the sparse channel under each signal-to-noise ratio, inputting the sparse information of the channel into the network, and obtaining a corresponding digital estimator so as to obtain a channel estimation result. The sparse channel estimator provided by the invention can reduce the error caused by the nonlinear quantization of a low-precision analog-to-digital converter (ADC), and is realized by using a deep learning network, so that the channel estimation complexity is reduced. In a mixed-frame multi-antenna system adopting a low-precision analog-to-digital converter, the performance of the method can approach the theoretically optimal channel estimation method.
As shown in fig. 1, under each signal-to-noise ratio, channel prior sparse information is obtained through a compressed sensing technology, and an optimal digital estimator under a corresponding minimum mean square error is designed according to all possible sparse channels and is used as a training input of a deep learning full-connection network. Through training and testing the deep learning network, the network which can be used under each signal-to-noise ratio is obtained and used for channel estimation. In a hybrid architecture, a transmitting end sends a pilot signal, an analog precoder consisting of fully-connected phase shifters reaches a receiver through a channel, and sparse information of the channel is input into a network for sparse channels under various signal-to-noise ratios through an analog estimator and a low-precision analog-to-digital converter quantizer, so that a corresponding digital estimator can be obtained, and a channel estimation result is obtained. The sparse channel estimator provided by the invention can reduce the error caused by the nonlinear quantization of the low-precision analog-to-digital converter, and is realized by using a deep learning network, so that the channel estimation complexity is reduced. In a mixed-frame multi-antenna system adopting a low-precision analog-to-digital converter, the performance of the method can approach the theoretically optimal channel estimation method.
The channel estimation method provided by the invention comprises the following steps:
the method comprises the following steps: base station design simulation pre-coding and sending pilot signal, simulation pre-coding FAmDesigned according to the following formula:
wherein, the [ alpha ], [ beta ]]ijElements representing the ith row and the jth column of the matrix; fAmWith a representation dimension of Nt×NRFtAnalog precoding matrix of, NtRepresenting the number of transmitting antennas, NRFtRepresenting the number of radio frequency links of a transmitting end;representing the phase of the ith row and jth column element of the analog precoding matrix;
step two: receiving end optimization design simulation estimator WAm,WAmCalculated according to the following formula:
wherein, WAmWith a representation dimension of Nr×NRFrAnalog estimation matrix of, NrRepresenting the number of receiving antennas, NRFrRepresenting the number of radio frequency links of a receiving end;representing a matrix of simulated estimatesThe phase of the ith row and jth column element;
step three: the receiving end obtains sparse characteristic information of the channel, namely the position of a non-zero element in a sparse channel matrix, through a compressed sensing technology, such as a classical Orthogonal Matching Pursuit (OMP) technology, and the sparse characteristic information is used as prior information of subsequent channel estimation.
Step four: all possible selection matrixes P corresponding to sparse channel characteristic information by the receiving endv[k]Is stored and a selection matrix P is selected according to eachv[k]Designing corresponding digital estimation matrix
Step five: selecting the matrix Pv[k]And a digital estimation matrix corresponding theretoAnd as a training set of the fully-connected neural network, obtaining the fully-connected feedforward network suitable for each signal-to-noise ratio through neural network training. Within the coherent time, the receiving end selects the matrix P according to the corresponding selection matrix of the sparse channelv[k]Obtaining a digital estimation matrix of the network output using the network structureDetecting the received signal by using the data estimation matrix to obtain a channel estimation value
Further, the matrix P is selected in the fourth stepv[k]Expressed by the following formula:
wherein e isπ(i)(π(i)∈{1,2,...,NrNt}) represents a dimension of NrNtThe pi (i) -th element of x 1 is a vector of 1 and the remaining elements are 0. N is a radical ofvRepresenting a dimension of N on the k-th subcarrierrNtVector component h of x 1 channel vector projected onto angular domainv[k]Number of medium non-zero elements.
For non-zero channel elements in the channel matrix, the selection matrix P of the channelv[k]Has N possible forms, for all possible sets Pv1[k],Pv2[k],...,PvN[k]Each Pvi[k]The possibility of occurrence of (i ═ 1.., N) is allNamely, it is
Wherein N is the total number of the possible numbers,c represents a combination number formula. N is a radical ofvRepresenting a dimension of N on the k-th subcarrierrNtVector component h of x 1 channel vector projected onto angular domainv[k]Number of medium non-zero elements.
Further, projecting the channel to the virtual angle domain to obtain a channel vector component hv[k],hv[k]Calculated according to the following formula:
wherein A istWith a representation dimension of Nt×NtA transmit dictionary matrix composed of transmit front response vectors,representation matrix AtConjugation of (1);represents the kronecker product; a. therWith a representation dimension of Nr×NrA receiving dictionary matrix composed of the receiving array surface response vectors; h [ k ]]Denotes the dimension N on the k sub-carrierr×NtVec (H [ k ]) of]) Representation matrix H [ k ]]Vectorization of (2).
Wherein A istExpressed by the following formula:
wherein,(p∈{1,2,...,Nt}) represents a dimension of NtX 1 of the transmitted wavefront response vector, whereWherein N istP denotes the number of antennas on the horizontal axis of the transmitting antenna array, and Q denotes the number of antennas on the vertical axis of the transmitting antenna array.
ArExpressed by the following formula:
wherein,(q∈{1,2,...,Nrdenotes the dimension NrX 1, where,wherein N isrI × J, I tableThe number of antennas on the horizontal axis of the receiving antenna array is shown, and J represents the number of antennas on the vertical axis of the receiving antenna array.
Further, the optimal number estimation matrix in the fourth stepUsing the channel estimation minimum mean square error criterion,it can be calculated as follows:
wherein K belongs to {1, 2.,. K } represents the kth subcarrier, and K represents the total number of subcarriers;representing N on the k sub-carrierRFr×NvM denotes the number of channel uses within the coherence time, i.e., the number of channel estimates ηbRepresents a distortion factor related to the quantization bit number b of an analog-to-digital converter (ADC); omega k]Denotes the dimension on the k sub-carrier as MNRFr×NvThe measurement matrix of (2); omegaH[k]Represents the matrix omega k]The conjugate transpose of (1);represents the large scale fading coefficient of the channel;with a representation dimension of Nv×NvThe identity matrix of (1);representing the variance of each element of the equivalent noise vector,to representAdditive White Gaussian Noise (AWGN) variance; p denotes a transmit pilot power.
Further, the channel estimation times M are expressed by the following formula:
indicating rounding up.
The measurement matrix Ω [ k ] is represented by the following formula:
Ω[k]=Φ[k]ΨPv[k],
wherein the dimension is NRFr×NrNtPilot correlation matrix ofsm[k](M ∈ {1, 2.., M }) denotes that the dimension at the time of the mth training is NRFrX 1 transmit pilot vector.
Where Ψ represents the dimension NrNt×NrNtThe spatial transformation matrix, Ψ, is represented by the following equation:
wherein A istWith a representation dimension of Nt×NtA transmit dictionary matrix composed of transmit array surface response vectors, ArWith a representation dimension of Nr×NrA receive dictionary matrix composed of receive front response vectors.
Wherein, Pv[k]With a representation dimension of NrNt×NvThe selection matrix of (2).
Further, in step fiveBy feed-forward as followsThe fully connected deep learning network yields:
(1) for a known number of non-zero elements NvOf the sparse channel, the channel matrix H [ k ]]Denotes the dimension N on the k sub-carrierr×NtUnder each signal-to-noise ratio, the receiving end calculates all possible P according to the sparse characteristic informationv[k]Formed covariance matrix C k]=ΩH[k]Ω[k]Wherein Ω [ k ]]=Φ[k]ΨPv[k]. At the same time, according to the method as claimed in claim 4, each P is calculatedv[k]Formed covariance matrix C k]Corresponding toDesigning a new numerical estimation matrixAnd will all possible C k]And correspondingStored in a database.
(2) Data are extracted from the database and divided into two groups, namely training data and testing data. Carrying out complex value splitting operation on the training data, and carrying out training set C [ k ]]And W'D[k]Split into a real matrix CR[k]、And an imaginary matrix CI[k]、Two parts.
(3) Then C is addedR[k]、And CI[k]、Performing matrix vectorization operation to obtain dimension Nv 2X 1 column vector cR[k]、And cI[k]、C is toR[k]、As input and training targets of the real part deep learning network, cI[k]、As input to the imaginary deep learning network and as a training target.
(4) Two deep learning fully-connected networks are constructed, and the network structures are the same. All are two layers of forward feedback fully-connected neural networks, the number of neurons in the first layer is N, whereinMeanwhile, a bias (bias) connection is set in the first layer, and the transfer function of the first layer is set as a softmax function; the first layer output is connected to the second layer neurons, the number of neurons is N, and the bias (bias) connection is also set at the second layer.
The softmax function is defined as follows:
wherein, the output of each layer of fully-connected network is as follows:
wherein W and b represent parameters of a fully-connected neural network, yiAnd biThe i-th element, x, representing y and bjDenotes the j-th element of x, Wi,jRepresents an element with position (i, j) in W.
(5) The fully-connected neural network is trained, and the deep learning network is tested by taking complex value splitting and matrix vectorization of test data as input, so that a stable network structure under each signal-to-noise ratio is obtained.
(6) In the coherent time, for the sparse channel under each signal-to-noise ratio, a channel covariance matrix C [ k ] containing sparse information of the channel]Input network to obtain corresponding digital estimator
Further, the channel estimation value in step fiveObtained according to the following formula:
wherein (C)HA conjugate transpose operation representing a matrix; ()-1Representing the inversion operation of the matrix.
As shown in fig. 2, under sparse channel, the channel estimation method proposed by the present invention is extremely close to the theoretically optimal Minimum Mean Square Error (MMSE) estimator, especially under low signal-to-noise ratio (SNR). Under a sparse channel, the invention provides the deep learning network, so that the precision of channel estimation is further improved and the complexity is reduced.
The technical means disclosed in the invention scheme are not limited to the technical means disclosed in the above embodiments, but also include the technical scheme formed by any combination of the above technical features. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principle of the present invention, and such improvements and modifications are also considered to be within the scope of the present invention.

Claims (7)

1. A millimeter wave sparse array surface channel estimation method based on a deep learning network is characterized by comprising the following steps:
the method comprises the following steps: base station design simulation pre-coding and sending pilot signal, simulation pre-coding FAmDesigned according to the following formula:
wherein, the [ alpha ], [ beta ]]ijI-th row, i-th row of the representation matrixElements of column j; fAmWith a representation dimension of Nt×NRFtAnalog precoding matrix of, NtRepresenting the number of transmitting antennas, NRFtRepresenting the number of radio frequency links of a transmitting end;representing the phase of the ith row and jth column element of the analog precoding matrix;
step two: receiving end optimization design simulation estimator WAm,WAmCalculated according to the following formula:
wherein, WAmWith a representation dimension of Nr×NRFrAnalog estimation matrix of, NrRepresenting the number of receiving antennas, NRFrRepresenting the number of radio frequency links of a receiving end;representing the phase of the ith row and jth column element of the simulation estimation matrix;
step three: the receiving end obtains sparse characteristic information of the channel through a compressed sensing technology, namely the position of a non-zero element in a sparse channel matrix, and the sparse characteristic information is used as prior information of subsequent channel estimation;
step four: all possible selection matrixes P corresponding to sparse channel characteristic information by the receiving endv[k]Is stored and a selection matrix P is selected according to eachv[k]Designing corresponding digital estimation matrix
Step five: selecting the matrix Pv[k]And a digital estimation matrix corresponding theretoAs a training set of fully connected neural networks, via the neural networkThe network training obtains a fully-connected forward feedback network suitable for each signal-to-noise ratio, and the receiving end selects a matrix P according to the sparse channelv[k]Obtaining a digital estimation matrix of the network output using the network structureDetecting the received signal by using the data estimation matrix to obtain a channel estimation value
2. The deep learning network-based millimeter wave front channel estimation method according to claim 1, wherein in the fourth step, a compressed sensing technique is used to obtain sparse feature information of the channel, and a corresponding selection matrix P is usedv[k]Expressed by the following formula:
wherein e isπ(i)(π(i)∈{1,2,…,NrNt}) represents a dimension of NrNtX 1. pi (i) th vector with 1 element and 0 elements, NvRepresenting a dimension of N on the k-th subcarrierrNtVector component h of x 1 channel vector projected onto angular domainv[k]The number of the non-zero elements in the group,
for non-zero channel elements in the channel matrix, the selection matrix P of the channelv[k]Has N possible forms, for all possible sets Pv1[k],Pv2[k],…,PvN[k]Each Pvi[k]The possibility of occurrence of (i ═ 1, …, N) is allNamely, it is
Wherein N is the total number of the possible numbers,c represents a combination number formula. N is a radical ofvRepresenting a dimension of N on the k-th subcarrierrNtVector component h of x 1 channel vector projected onto angular domainv[k]Number of medium non-zero elements.
3. The millimeter wave sparse front channel estimation method based on the deep learning network as claimed in claim 2, wherein the channel vector component h is obtained by projecting the channel onto a virtual angle domainv[k],hv[k]Calculated according to the following formula:
wherein A istWith a representation dimension of Nt×NtA transmit dictionary matrix composed of transmit front response vectors,representation matrix AtConjugation of (1);represents the kronecker product; a. therWith a representation dimension of Nr×NrA receiving dictionary matrix composed of the receiving array surface response vectors; h [ k ]]Denotes the dimension N on the k sub-carrierr×NtVec (H [ k ]) of]) Representation matrix H [ k ]]Vectorization of (2);
wherein A istExpressed by the following formula:
wherein,with a representation dimension of NtX 1 of the transmitted wavefront response vector, whereWherein N istP × Q, P denotes the number of antennas on the horizontal axis of the transmitting antenna array, and Q denotes the number of antennas on the vertical axis of the transmitting antenna array;
Arexpressed by the following formula:
wherein,with a representation dimension of NrX 1, where,wherein N isrI denotes the number of antennas on the horizontal axis of the reception antenna array, and J denotes the number of antennas on the vertical axis of the reception antenna array.
4. The deep learning network-based millimeter wave front channel estimation method according to claim 3, wherein the optimal number estimation matrix in the fourth stepUsing the channel estimation minimum mean square error criterion,calculated according to the following formula:
wherein, k ∈ {1,2, …, K denotes the kth subcarrier, K denotes the total number of subcarriers;representing N on the k sub-carrierRFr×NvM denotes the number of channel uses, i.e. the number of channel estimates, within the coherence time ηbRepresents a distortion factor related to the quantization bit number b of an analog-to-digital converter (ADC); omega k]Denotes the dimension on the k sub-carrier as MNRFr×NvThe measurement matrix of (2); omegaH[k]Represents the matrix omega k]The conjugate transpose of (1);represents the large scale fading coefficient of the channel;with a representation dimension of Nv×NvThe identity matrix of (1);representing the variance of each element of the equivalent noise vector,represents an Additive White Gaussian Noise (AWGN) variance; p denotes a transmit pilot power.
5. The deep learning network-based millimeter wave front channel estimation method according to claim 4, wherein the channel estimation times M are expressed by the following formula:
represents rounding up;
the measurement matrix Ω [ k ] is represented by the following formula:
Ω[k]=Φ[k]ΨPv[k],
wherein the dimension is NRFr×NrNtPilot correlation matrix ofsm[k](M ∈ {1,2, …, M }) denotes the dimension N at the mth trainingRFrX 1 transmit pilot vector.
Where Ψ represents the dimension NrNt×NrNtThe spatial transformation matrix, Ψ, is represented by the following equation:
wherein A istWith a representation dimension of Nt×NtA transmit dictionary matrix composed of transmit array surface response vectors, ArWith a representation dimension of Nr×NrA receiving dictionary matrix composed of the receiving array surface response vectors;
wherein, Pv[k]With a representation dimension of NrNt×NvThe selection matrix of (2).
6. The deep learning network-based millimeter wave front channel estimation method according to claim 1, wherein in step fiveThe method is obtained through a full-connection network with forward feedback, and comprises the following steps:
(1) for a known number of non-zero elements NvOf the sparse channel, the channel matrix H [ k ]]Denotes the dimension N on the k sub-carrierr×NtUnder each signal-to-noise ratio, the receiving end calculates all possible P according to the sparse characteristic informationv[k]Formed covariance matrix C k]=ΩH[k]Ω[k]Wherein Ω [ k ]]=Φ[k]ΨPv[k]. A method, meter according to claim 4Calculate each Pv[k]Formed covariance matrix C k]Corresponding toDesigning a new numerical estimation matrixAnd will all possible C k]And correspondingStoring the data in a database;
(2) extracting data from database, dividing the data into two groups of training data and test data, carrying out complex value splitting operation on the training data, and collecting training set C [ k ]]And W'D[k]Split into a real matrix CR[k]、And an imaginary matrix CI[k]、Two parts;
(3) then C is addedR[k]、And CI[k]、Performing matrix vectorization operation to obtain dimension Nv 2X 1 column vector cR[k]、And cI[k]、C is toR[k]、As input and training targets of the real part deep learning network, cI[k]、As input and training targets for the imaginary deep learning network;
(4) two deep learning fully-connected networks are constructed, the network structures are the same, the two deep learning fully-connected networks are both two layers of forward feedback fully-connected neural networks, the number of neurons in the first layer is N, whereinMeanwhile, a bias (bias) connection is set in the first layer, and the transfer function of the first layer is set as a softmax function; connecting the first layer output with the second layer neuron, wherein the number of the neurons is N, and setting bias (bias) connection on the second layer;
the softmax function is defined as follows:
wherein, the output of each layer of fully-connected network is as follows:
wherein W and b represent parameters of a fully-connected neural network, yiAnd biThe i-th element, x, representing y and bjDenotes the j-th element of x, Wi,jRepresents an element with position (i, j) in W;
(5) training the fully-connected neural network, and testing the deep learning network by taking complex value splitting and matrix vectorization of test data as input, thereby obtaining a stable network structure under each signal-to-noise ratio;
(6) in the coherent time, for the sparse channel under each signal-to-noise ratio, a channel covariance matrix C [ k ] containing sparse information of the channel]The input of the network is carried out,corresponding digital estimators can be obtained
7. The deep learning network-based millimeter wave front channel estimation method according to claim 1, wherein the channel estimation value in step (5)Calculated according to the following formula:
wherein (C)HA conjugate transpose operation representing a matrix; ()-1Representing the inversion operation of the matrix.
CN201910397076.5A 2019-05-14 2019-05-14 Millimeter wave sparse array surface channel estimation method based on deep learning network Active CN110099016B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910397076.5A CN110099016B (en) 2019-05-14 2019-05-14 Millimeter wave sparse array surface channel estimation method based on deep learning network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910397076.5A CN110099016B (en) 2019-05-14 2019-05-14 Millimeter wave sparse array surface channel estimation method based on deep learning network

Publications (2)

Publication Number Publication Date
CN110099016A true CN110099016A (en) 2019-08-06
CN110099016B CN110099016B (en) 2022-05-31

Family

ID=67447898

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910397076.5A Active CN110099016B (en) 2019-05-14 2019-05-14 Millimeter wave sparse array surface channel estimation method based on deep learning network

Country Status (1)

Country Link
CN (1) CN110099016B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110519188A (en) * 2019-08-20 2019-11-29 电子科技大学 A kind of compressed sensing based multi-user's time-varying millimeter wave channel estimation methods
CN110535500A (en) * 2019-09-03 2019-12-03 电子科技大学 A kind of millimeter wave MIMO mixed-beam forming optimization method based on deep learning
CN111713035A (en) * 2020-04-07 2020-09-25 东莞理工学院 MIMO multi-antenna signal transmission and detection technology based on artificial intelligence
CN111835351A (en) * 2020-06-01 2020-10-27 东南大学 Deep learning-based quantitative signal reconstruction method and arrival angle estimation method
CN112187323A (en) * 2020-09-29 2021-01-05 国网江苏省电力有限公司丹阳市供电分公司 IRS-based large-scale MIMO (multiple input multiple output) cascade channel estimation method under mixed low-precision architecture
CN112615801A (en) * 2020-12-16 2021-04-06 西安交通大学 Channel estimation method, medium, and apparatus based on compressed sensing and deep learning
CN113193893A (en) * 2021-04-30 2021-07-30 东南大学 Millimeter wave large-scale MIMO intelligent hybrid beam forming design method
CN113972939A (en) * 2021-09-09 2022-01-25 浙江大学 Antenna system precoding method and device based on double time scales and deep learning
CN114584448A (en) * 2022-02-16 2022-06-03 山东大学 SM-OFDM signal grouping detection method based on deep neural network
CN114759956A (en) * 2022-04-13 2022-07-15 东南大学 Single-bit ADC uplink multi-user MIMO deep unfolding pre-coding method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107566303A (en) * 2017-07-27 2018-01-09 东华大学 A kind of millimeter wave channel estimation methods based on Bayes's compressed sensing
US9924880B2 (en) * 2015-02-11 2018-03-27 Samsung Electronics Co., Ltd. RF doppler bio-signal sensor for continuous heart rate variability and blood pressure monitoring
WO2018093411A2 (en) * 2016-11-16 2018-05-24 Intel IP Corporation Station (sta), access point (ap) and methods for beam refinement and related signaling
CN108881074A (en) * 2018-05-08 2018-11-23 东南大学 Broadband millimeter-wave channel estimation methods under a kind of low precision mixed architecture
CN109005133A (en) * 2018-07-12 2018-12-14 南京邮电大学 Double Sparse multi-path channel models and the channel estimation methods based on this model
WO2019014264A1 (en) * 2017-07-12 2019-01-17 Intel IP Corporation Enhanced frame format for single user 2.16 ghz channel wireless communications
CN109743268A (en) * 2018-12-06 2019-05-10 东南大学 Millimeter wave channel estimation and compression method based on deep neural network

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9924880B2 (en) * 2015-02-11 2018-03-27 Samsung Electronics Co., Ltd. RF doppler bio-signal sensor for continuous heart rate variability and blood pressure monitoring
WO2018093411A2 (en) * 2016-11-16 2018-05-24 Intel IP Corporation Station (sta), access point (ap) and methods for beam refinement and related signaling
WO2019014264A1 (en) * 2017-07-12 2019-01-17 Intel IP Corporation Enhanced frame format for single user 2.16 ghz channel wireless communications
CN107566303A (en) * 2017-07-27 2018-01-09 东华大学 A kind of millimeter wave channel estimation methods based on Bayes's compressed sensing
CN108881074A (en) * 2018-05-08 2018-11-23 东南大学 Broadband millimeter-wave channel estimation methods under a kind of low precision mixed architecture
CN109005133A (en) * 2018-07-12 2018-12-14 南京邮电大学 Double Sparse multi-path channel models and the channel estimation methods based on this model
CN109743268A (en) * 2018-12-06 2019-05-10 东南大学 Millimeter wave channel estimation and compression method based on deep neural network

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YUCHENG WANG: "Wideband mmWave Channel Estimation for Hybrid Massive MIMO With Low-Precision ADCs", 《IEEE》 *
杜若宸: "基于毫米波的5G通信定位一体化系统", 《通信设计与应用》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110519188A (en) * 2019-08-20 2019-11-29 电子科技大学 A kind of compressed sensing based multi-user's time-varying millimeter wave channel estimation methods
CN110535500A (en) * 2019-09-03 2019-12-03 电子科技大学 A kind of millimeter wave MIMO mixed-beam forming optimization method based on deep learning
CN110535500B (en) * 2019-09-03 2021-09-21 电子科技大学 Millimeter wave MIMO mixed beam forming optimization method based on deep learning
CN111713035A (en) * 2020-04-07 2020-09-25 东莞理工学院 MIMO multi-antenna signal transmission and detection technology based on artificial intelligence
CN111835351B (en) * 2020-06-01 2023-07-28 东南大学 Quantized signal reconstruction method and arrival angle estimation method based on deep learning
CN111835351A (en) * 2020-06-01 2020-10-27 东南大学 Deep learning-based quantitative signal reconstruction method and arrival angle estimation method
CN112187323A (en) * 2020-09-29 2021-01-05 国网江苏省电力有限公司丹阳市供电分公司 IRS-based large-scale MIMO (multiple input multiple output) cascade channel estimation method under mixed low-precision architecture
CN112615801A (en) * 2020-12-16 2021-04-06 西安交通大学 Channel estimation method, medium, and apparatus based on compressed sensing and deep learning
CN113193893A (en) * 2021-04-30 2021-07-30 东南大学 Millimeter wave large-scale MIMO intelligent hybrid beam forming design method
CN113972939A (en) * 2021-09-09 2022-01-25 浙江大学 Antenna system precoding method and device based on double time scales and deep learning
CN113972939B (en) * 2021-09-09 2022-07-12 浙江大学 Antenna system precoding method and device based on double time scales and deep learning
WO2023035736A1 (en) * 2021-09-09 2023-03-16 浙江大学 Antenna system precoding method and apparatus based on two time scales and deep learning
CN114584448A (en) * 2022-02-16 2022-06-03 山东大学 SM-OFDM signal grouping detection method based on deep neural network
CN114584448B (en) * 2022-02-16 2023-09-22 山东大学 SM-OFDM signal grouping detection method based on deep neural network
CN114759956A (en) * 2022-04-13 2022-07-15 东南大学 Single-bit ADC uplink multi-user MIMO deep unfolding pre-coding method
CN114759956B (en) * 2022-04-13 2024-01-30 东南大学 Single-bit ADC uplink multi-user MIMO deep spreading precoding method

Also Published As

Publication number Publication date
CN110099016B (en) 2022-05-31

Similar Documents

Publication Publication Date Title
CN110099016B (en) Millimeter wave sparse array surface channel estimation method based on deep learning network
CN110099017B (en) Channel estimation method of hybrid quantization system based on deep neural network
Hu et al. Super-resolution channel estimation for mmWave massive MIMO with hybrid precoding
CN108933745B (en) Broadband channel estimation method based on super-resolution angle and time delay estimation
CN112737985B (en) Large-scale MIMO channel joint estimation and feedback method based on deep learning
CN107018099B (en) Time-varying channel estimation method for millimeter wave multi-user MIMO system
CN110855585B (en) Channel estimation method of millimeter wave large-scale MIMO system
CN106559367A (en) MIMO ofdm system millimeter wave channel estimation methods based on low-rank tensor resolution
CN107135024A (en) A kind of mixed-beam figuration Iterative Design method of low complex degree
CN108881074B (en) Broadband millimeter wave channel estimation method under low-precision hybrid architecture
CN110401476A (en) A kind of codebook-based millimetre-wave attenuator multi-user parallel beam training method
CN109257309A (en) A kind of high performance extensive MIMO downlink transmission channel estimation method
CN109347529B (en) Channel estimation and hybrid beam forming method for resisting non-ideality of phase shifter
CN114567525B (en) Channel estimation method and device
Rodríguez-Fernández et al. A frequency-domain approach to wideband channel estimation in millimeter wave systems
CN112448742A (en) Hybrid precoding method based on convolutional neural network under non-uniform quantization
CN110719127B (en) Millimeter wave MIMO system beam forming method with constant modulus constraint
CN112769462A (en) Millimeter wave MIMO broadband channel estimation method based on joint parameter learning
Zamanipour A survey on deep-learning based techniques for modeling and estimation of massivemimo channels
Shoukath Ali et al. Time domain channel estimation for time and frequency selective millimeter wave MIMO hybrid architectures: sparse Bayesian learning-based Kalman filter
CN112398513A (en) Beam forming method of massive MIMO system
CN110086734B (en) Channel tracking and recovery method in time-varying large-scale MIMO system
CN116405077A (en) Large-scale MIMO hybrid beamforming method based on deep learning
CN110213186A (en) The angle of arrival and signal-noise ratio estimation method of ESPRIT based on mixed structure
Zheng et al. Joint position, orientation and channel estimation in hybrid mmWave MIMO systems

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant