CN104052691A - MIMO-OFDM system channel estimation method based on compressed sensing - Google Patents
MIMO-OFDM system channel estimation method based on compressed sensing Download PDFInfo
- Publication number
- CN104052691A CN104052691A CN201410312859.6A CN201410312859A CN104052691A CN 104052691 A CN104052691 A CN 104052691A CN 201410312859 A CN201410312859 A CN 201410312859A CN 104052691 A CN104052691 A CN 104052691A
- Authority
- CN
- China
- Prior art keywords
- path
- channel
- vector
- centerdot
- matrix
- 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
Links
Landscapes
- Radio Transmission System (AREA)
Abstract
The invention provides an MIMO-OFDM system channel estimation method based on compressed sensing. The MIMO-OFDM system channel estimation method based on compressed sensing is mainly applied to channel estimation when a receiving terminal is provided with a two-dimensional antenna array. According to the MIMO-OFDM system channel estimation method based on compressed sensing, the time delay, incidence angle and gain of each path of a space channel are estimated in sequence, and channel estimation accuracy can be improved effectively. The MIMO-OFDM system channel estimation method based on compressed sensing comprises the following steps that 1, an initially-estimated value of a channel frequency domain response vector of each pilot frequency sub-carrier is obtained according to the least square criterion; 2, by means of the sparsity of the channel frequency domain response vectors in a time delay domain, the time delay of each path of the channel and an estimated value of a channel time domain response vector of each path of the channel are estimated on the basis of the compressed sensing theory; 3, by means of the sparsity of the channel time domain response vectors in a two-dimensional angle domain, the incidence angle of each path of the channel is estimated on the basis of the compressed sensing theory; 4, the gain coefficient of each path of the channel is estimated according to the least square criterion; 5, estimated values of channel frequency domain responses of all the sub-carriers and antennas are obtained.
Description
Technical field
The present invention relates to the MIMO-OFDM wireless communication system that receiving terminal is equipped with two-dimensional antenna array, particularly the MIMO-OFDM system channel estimation problem based on compressive sensing theory.
Background technology
Multiple-input and multiple-output (MIMO) technology is at transmitting terminal and receiving terminal configuration multiple antenna, by the combination with multiplex technique or diversity technique etc., make full use of the multipath characteristics of scatter channel, in independent, the parallel data flow of space transmission multichannel, in the situation that not increasing system bandwidth, improve exponentially capability of wireless communication system and link reliability, improved the transmission rate of system.Therefore, MIMO technology is subject to extensive concern, is considered to one of important breakthrough of Modern wireless communication and key technology that future wireless system must adopt.In order to meet growing user data demand, antenna configuration number constantly increases, and as extensive mimo system, and two-dimensional antenna array is by the choose reasonable that is this system antenna configuration.
OFDM (OFDM) is a kind of orthogonal multiple carrier modulation technique, has the efficient availability of frequency spectrum.Development along with radio communication, in order to meet people's demand growing to high-rate service, system bandwidth constantly increases, and the frequency selectivity of channel is more outstanding, traditional balancing technique not only complexity is high, and is difficult to eliminate completely the caused intersymbol interference of multidiameter fading channel.Ofdm system converts wideband frequency Selective Fading Channel to a series of arrowbands flat fading channel, effectively solves the problem of intersymbol interference, at aspects such as realizing high speed data transfer, has unique advantage.Therefore, the system architecture based on MIMO and OFDM arises at the historic moment.
Compressed sensing, claims again compression sampling, is a new sampling theory, by the sparse characteristic of exploitation signal, under the condition much smaller than Nyquist sample rate, by stochastical sampling, obtain the discrete sample of signal, then by the perfect reconstruction signal of non-linear algorithm for reconstructing.Due in actual scattering environments, for broadband signal, wireless channel often consists of several main paths, therefore wireless channel can be considered as to a time delay domain condition of sparse channel, utilizes compressive sensing theory to estimate channel response with less number of pilots.Meanwhile, because user is limited to the incident angle number of base station end aerial array, so exploitation channel is in the sparse property of angle domain, can utilize equally compressive sensing theory to obtain the estimation of channel.
Summary of the invention
Technical problem: the invention provides a kind of number of pilots that effectively reduces, there is the MIMO-OFDM system channel estimation method based on compressed sensing of more over-evaluating meter resolution.
Technical scheme: the MIMO-OFDM system channel estimation method based on compressed sensing of the present invention, comprises the following steps:
A) utilize criterion of least squares, obtain the initial estimate of each pilot sub-carrier upper signal channel frequency domain response vector;
B) utilize the initial estimate of the channel frequency domain response vector that described step obtains in a) in the sparse property of time delay domain, based on compressive sensing theory, estimate the estimated value of the time delay in each path of channel and the channel time domain response vector in each path;
C) utilize described step b) in the estimated value of the channel time domain response vector estimated in the sparse property in two dimension angular territory, based on compressive sensing theory, estimate the incidence angle in each path of channel, described incidence angle comprises the vertical angle of pitch and horizontal azimuth;
D) utilize initial estimate, the step b of the channel frequency domain response vector that described step obtains in a)) in time delay, the step c in each path of estimating) incidence angle in each path of estimating, adopt criterion of least squares, estimate the gain coefficient in each path of channel;
E) by described step b) in time delay, the step c in each path of estimating) incidence angle and the steps d in each path of estimating) and in the gain coefficient in estimated each path, the multipath channel models of substitution channel frequency domain response, obtains the estimated value of the channel frequency domain response on all subcarriers, all antennas.
In the preferred version of the inventive method, step b) idiographic flow is:
First, by channel path time delay 0 to uniform discrete between maximum path time delay, obtain minimum time delay interval and be Δ τ and discrete after path number be N
τ, the Fourier transform relation of utilizing channel frequency domain response and channel time domain to respond, is the product of base vector and channel time domain matrix by the channel frequency domain response vector representation on each subcarrier.
Secondly, the initial estimate of the channel frequency domain response vector on all pilot sub-carriers is arranged in to a matrix in order, is designated as channel frequency domain matrix; All base vectors are arranged in to a matrix according to identical order, are designated as projection matrix; Thereby by channel frequency domain Matrix Estimation value representation be:
Wherein,
represent channel frequency domain Matrix Estimation value, Ξ is projection matrix, and A is channel time domain matrix, and W is noise matrix;
Finally, according to formula (12), adopt compressive sensing theory to solve channel time domain matrix, the channel time domain matrix solving is N altogether
τoK, find non-zero row vector wherein, the number L of described non-zero row vector is the channel path sum estimating, l non-zero row vector is the estimated value of the channel time domain response vector in l path, l the corresponding path delay in path is (m-1) Δ τ, wherein, l is the sequence number of non-zero row vector, and m is the channel time domain matrix line number m at l non-zero row vector place.
In the preferred version of the inventive method, step c) idiographic flow is:
First, by vertical angle of pitch uniform discrete, obtain minimum angles and be spaced apart Δ θ, the angle number after discrete is N
θ; By horizontal azimuth uniform discrete, obtain minimum angles and be spaced apart Δ φ, the angle number after discrete is N
φ, utilize the two-dimensional Fourier transform relation between channel time domain response and path incidence angle, the channel time domain response vector estimated value table in each path is shown:
Wherein,
be the channel time domain response vector estimated value in l path, P is projection matrix, γ
lbe the path complex gain vector in l path, w is noise vector;
Then calculate respectively in accordance with the following methods the incidence angle in each path: for l path, utilize formula (18), adopt compressive sensing theory to solve the estimated value that obtains path complex gain vector, according to the peaked position of element in the complex gain vector of path, calculate the incidence angle in l path, described incidence angle comprises the vertical angle of pitch and horizontal azimuth.
In the above-mentioned preferred version of the inventive method, step c) in, every row element arrangement mode and path complex gain vector γ in described projection matrix
lmiddle arrangement of elements mode is identical, is according to the order sequence of the Kronecker product of the vertical angle of pitch and horizontal azimuth institute respective path and obtains;
The method of calculating the incidence angle in l path is: according to the peaked position s of element in the complex gain vector of path, the vertical angle of pitch that obtains l path is
horizontal azimuth is
wherein, n
φ=mod (s, N
φ), n
θ=(s-n
φ)/N
φ+ 1, mod (x, y) represents that x is to y complementation.
In the another kind of preferred version of the inventive method, step c) in, arrangement of elements mode and path complex gain vector γ in described projection matrix
lmiddle arrangement of elements mode is identical, is according to the order of the Kronecker product of horizontal azimuth and vertical angle of pitch institute respective path;
The method of calculating the incidence angle in l path is: according to the peaked position s of element in the complex gain vector of path, the vertical angle of pitch that obtains l path is
horizontal azimuth is
wherein, n
θ=mod (s, N
θ), n
φ=(s-n
θ)/N
θ+ 1.
Beneficial effect: compared with prior art, tool has the following advantages in the present invention:
1, channel estimation methods of the present invention is on the basis of least-squares estimation, further estimated to characterize the parameter (time delay, incidence angle and the gain coefficient that comprise path) of the characteristic of channel, its channel estimation results has significantly and improves compared to least-squares estimation.
2, traditional MIMO-OFDM system channel estimation scheme, as least-squares estimation, number of pilots in an OFDM symbol must be not less than channel time domain length, thereby guarantee that pilot tone is less than coincidence frequency interval at the distribution interval of frequency domain, otherwise can cause that channel estimation results is at the aliasing of time domain.The present invention has utilized the sparse property of channel in time delay domain, and in conjunction with the theory of compressed sensing, number of pilots can be less than channel time domain length, thereby has reduced pilot-frequency expense.
3, traditional incident angle estimation scheme based on subspace algorithm, as MUSIC algorithm, needs a large amount of sample estimate covariance battle arrays.Limited when number of samples, or sample is while having correlation, and estimated performance is very poor, and the number of the incidence angle that can differentiate is less than number of antennas.The present invention utilizes channel in the sparse property in two dimension angular territory, in conjunction with the theory of compressed sensing, estimates incidence angle, not only can effectively reduce observation sample number, and not be subject to the impact of sample correlations and number of antennas, has more wide application.
The present invention adopts the method based on compressed sensing to obtain the estimation of channel parameter, take full advantage of channel in the sparse property in time delay domain and two dimension angular territory, estimated result compares to least-squares estimation tool and increases significantly, reduced pilot-frequency expense simultaneously, for broadband and extensive mimo system, be with a wide range of applications.
Accompanying drawing explanation
Fig. 1 is receiving terminal antenna array configuration schematic diagram.
Embodiment
In order to make those skilled in the art person understand better the present invention, below in conjunction with embodiment and Figure of description, the technical scheme in the present invention is carried out to clear, complete description explanation, should understand these embodiment only for the embodiment of technical solution of the present invention is described, and be not used in, limit the scope of the invention.After having read the present invention, those skilled in the art all fall within to the modification of various equivalents of the present invention and replacement the protection range that the application's claim limits.
Embodiment 1:
1, system model
This programme is considered MIMO-OFDM system.Fig. 1 is receiving terminal antenna configuration schematic diagram, is not general, adopts two-dimensional antenna uniform surface battle array, and wherein horizontal direction has N
hroot antenna, vertical direction has N
vroot antenna, adjacent antenna spacing is d.Time domain multipath channel response matrix is
Wherein, δ (τ) represents Kronecker impulse function, the footpath number that L is multipath channel, τ
l, β
lbe respectively time delay and the gain coefficient in l path.
for the time domain channel matrix of user to aerial array.Suppose that user distance aerial array is enough far away, incident wave is plane wave, and the vertical direction in l footpath and horizontal direction steering vector are respectively
Wherein, θ
l∈ [0, pi/2], φ
l∈ [0, π] is respectively vertical angle of pitch angle and the horizontal azimuth in l footpath, and λ is incident wave wavelength.The multipath channel models of utilizing Fourier transform to obtain channel frequency domain response is:
Wherein, t
sfor the sampling period, N
cfor the number of OFDM symbol subcarrier, H
kbe the channel frequency domain matrix of k subcarrier, the element representative of consumer that its v is capable, h is listed as, to vertical direction v root antenna, the channel frequency domain response of horizontal direction h root antenna on k subcarrier, is expressed as
Wherein, user to vertical direction v root antenna, horizontal direction h root antenna in the channel time domain response in l path is
Suppose that user's number is 1 and is equipped with single antenna, system model is
Y
k=H
kX
k+Z
k (5)
Wherein,
x
krepresent respectively k the reception data on subcarrier and send data, Z
kfor additive white Gaussian noise matrix.To formula (5) equal sign left and right matrix-vector,
y
k=h
kX
k+z
k (6)
Wherein, y
k=vec (Y
k), h
k=vec (H
k), z
k=vec (Z
k), vec () representing matrix vectorization operator.
2, channel initial estimation
Suppose that user's pilot symbol transmitted is X
k, it is y that receiving terminal receives pilot data
k, according to system model formula (6), utilizing criterion of least squares, the initial estimate that obtains each pilot sub-carrier upper signal channel frequency domain response vector is
Suppose that pilot tone is uniformly distributed at frequency domain, wherein pilot interval is k
p, pilot tone number is N
p, pilot set is { 0, k
p..., (N
p-1) k
p.This channel estimation methods is simple, but precision is limited.Next, this programme will estimate time delay, incidence angle and the gain coefficient of the every paths of multipath channel successively based on this initial estimation result, thereby obtain more accurate channel estimation results.
3, time delay is estimated
If channel path maximum delay is τ
max, because every channel path time delay is all at 0~τ
maxbetween random distribution, can be by channel path time delay 0 to τ
maxbetween uniform discrete, discrete interval is Δ τ, the path number after discrete is
wherein,
be expressed as the result that a is rounded up, by formula (3), can obtain being expressed as of channel frequency domain response
Wherein, time delay set Γ's is chosen for
Γ={0,Δτ,2Δτ,…,(N
τ-1)Δτ} (9)
The channel frequency domain response of all antennas on k subcarrier is gathered together, be designated as channel frequency domain response vector, be expressed as
h
k=Au
k (10)
The channel frequency domain vector of k subcarrier wherein
Base vector
Channel time domain matrix
According to formula (7) and formula (10), the initial estimate of channel frequency domain response vector is expressed as
The initial estimate of the channel frequency domain response vector on all pilot sub-carriers is arranged in to channel frequency domain matrix in order
all base vectors are arranged in to projection matrix Ξ according to identical order, thus obtain by channel frequency domain Matrix Estimation value representation, be:
Wherein,
because the most elements of channel time domain matrix A are 0, be a sparse matrix, so formula (12) be the initial estimate of channel frequency domain response vector at the rarefaction representation of time delay domain, based on compressive sensing theory, solve channel time domain matrix A and be
Wherein, μ is regularization parameter, this value choose the precision that affects signal reconstruction, || ||
1, || ||
21-norm and the 2-norm of difference representing matrix.At solved channel time domain matrix
in, if certain a line is not null vector, there is a corresponding time delay path, as channel time domain matrix
m capable be l non-zero row vector, the channel time domain response vector estimated value of l paths is
Corresponding path delay is estimated as
suppose channel time domain matrix
middle non-zero row vector altogether
individual, the channel path estimating adds up to
the corresponding delay set estimating is
4, incidence angle is estimated
Due to vertical angle of pitch random distribution between 0~0.5 π, by vertical angle of pitch uniform discrete, minimum angles is spaced apart Δ θ, and the angle number after discrete is
wherein,
be expressed as the result that a is rounded downwards; Due to horizontal azimuth random distribution between 0~π, by horizontal azimuth uniform discrete, minimum angles is spaced apart Δ φ, and the angle number after discrete is
by (4), channel time domain response can be expressed as
Wherein, angle set Θ, Φ is chosen for
According to the estimated result of formula (14) and formula (16), the channel time domain response vector estimated value table of estimated l paths is shown
Wherein, path complex gain vector
w is noise vector, and P is projection matrix, is expressed as
Wherein,
φ
s=sΔφ,θ
q=qΔθ。Due to path complex gain vector γ
lin to only have an element be not 0, be a sparse vector, so formula (18) is that the estimated value of channel time domain response vector is at the rarefaction representation in two dimension angular territory, based on compressive sensing theory solution path complex gain vector γ
lfor
The estimated path complex gain vector going out
the size of incident angle has been reflected in the residing position of middle element maximum.According to vectorial γ
lthe arrangement mode of middle element, can calculate the vertical angle of pitch and horizontal azimuth.If the path complex gain vector estimating
the peaked position of middle element is s,
(1) if path complex gain vector γ
lthe arrangement mode of middle element is the order of the Kronecker product of the vertical angle of pitch and horizontal azimuth institute respective path, first, for each vertical angle of pitch, by horizontal azimuth order from small to large, arranges path complex gain, and obtaining a length is N
φvector; Then, according to vertical angle of pitch order from small to large, arranging above-mentioned length is N
φvector, obtaining a length is N
θn
φvector, be expressed as
The corresponding vertical angle of pitch of l paths and horizontal azimuth is estimated as
Wherein, n
φ=mod (s, N
φ), n
θ=(s-n
φ)/N
φ+ 1, mod (x, y) represents that x is to y complementation.
(2) if path complex gain vector γ
lthe arrangement mode of middle element is horizontal azimuth and the order of the Kronecker product of vertical angle of pitch institute respective path, first, for each horizontal azimuth, by vertical angle of pitch order from small to large, arranges path complex gain, and obtaining a length is N
θvector; Then, according to horizontal azimuth order from small to large, arranging above-mentioned length is N
θvector, obtaining a length is N
θn
φvector, be expressed as
The corresponding vertical angle of pitch of l paths and horizontal azimuth is estimated as
Wherein, n
θ=mod (s, N
θ), n
φ=(s-n
θ)/N
θ+ 1.
Estimate successively
the set of the incident angle of paths (the vertical angle of pitch and horizontal azimuth) is
5, gain coefficient is estimated
According to the estimated value of the incidence angle in the estimated value of the time delay in the initial estimate of channel frequency domain response, each path, each path, by formula (3) and formula (4), can be obtained
Wherein,
w
v, h, kbe respectively user to initial estimate and the noise of horizontal direction h root antenna, the channel frequency domain response of vertical direction v root antenna on k subcarrier,
be respectively time delay, horizontal azimuth and the vertical angle of pitch in estimated l the path going out, β
lit is the gain coefficient in l path.The initial estimate of the channel frequency domain response on all pilot sub-carriers is arranged in to column vector, according to formula (22),
Wherein,
The path delay and the incidence angle that according to formula (15) and formula (21), obtain, matrix E builds as follows
Wherein,
Adopt criterion of least squares, obtain channel gain coefficient and be
6, channel estimating
According to estimated path delay, incidence angle and gain coefficient, the multipath channel models of substitution channel frequency domain response (2), the estimated value that obtains the channel frequency domain response on all subcarriers and antenna is:
Embodiment 2:
According to compressive sensing theory, can obtain formula (13) and formula (20), adopt optimized algorithm to solve.The present embodiment provides OMP Algorithm for Solving formula (13) and formula (20), is specially:
A. initialization residual
K wherein
pfor pilot interval, N
pfor pilot tone number,
for user is to vertical direction v root antenna, horizontal direction h root antenna is at the initial estimate of the channel frequency domain response of k subcarrier.The set of initialization time delay index
iterations t=1.Make Ψ
(0)for empty matrix.Suppose that multipath number L is known.
B. in the t time iteration, by solving following formula, find the location index in t time delay footpath:
Wherein,
be the j row of projection matrix Ξ, projection matrix Ξ is expressed as
Wherein,
T
sfor the sampling interval.
C. preserve the estimated result of the t time iteration:
D. solve the channel time domain response vector in the path that has completed time delay estimation
E. upgrade residual volume
F. make t=t+1.If t≤L, returns to step b and carries out next iteration; Otherwise, continue.
G. step b~step f has searched out L paths altogether, and the path delay location index set wherein estimating is Λ
(L).For l ∈ Λ arbitrarily
(L), there is a paths, its time delay is
corresponding path gain is
user to the estimated value of the channel time domain response vector of the l paths on all antennas is
Time delay estimates that set sequence is
H. for the estimated L paths going out, search for successively corresponding incidence angle, wherein, the incident angle index value of l paths is
Wherein, p
jj row for projection matrix P.Projection matrix P is provided by following formula
Wherein,
θ
q=(q-1)/c, the minimum antenna distance that d is even antenna plane, λ is carrier wavelength.According to the sortord of row vector element in projection matrix P, can be in the hope of corresponding incidence angle, wherein, every row element arrangement mode and path complex gain vector γ above in projection matrix P
lmiddle arrangement of elements mode is identical, therefore
(1) if path complex gain vector γ
lthe arrangement mode of middle element is the order of the Kronecker product of the vertical angle of pitch and horizontal azimuth institute respective path, first, for each vertical angle of pitch, by horizontal azimuth order from small to large, arranges path complex gain, and obtaining a length is N
φvector; Then, according to vertical angle of pitch order from small to large, arranging above-mentioned length is N
φvector, obtaining a length is N
θn
φvector, be expressed as
The corresponding vertical angle of pitch of l paths and horizontal azimuth is estimated as
Wherein, n
φ=mod (s, N
φ), n
θ=(s-n
φ)/N
φ+ 1, mod (x, y) represents that x is to y complementation.
(2) if the arrangement mode of element is horizontal azimuth and the order of the Kronecker product of vertical angle of pitch institute respective path in path complex gain vector γ l, first, for each horizontal azimuth, by vertical angle of pitch order from small to large, arrange path complex gain, obtaining a length is N
θvector; Then, according to horizontal azimuth order from small to large, arranging above-mentioned length is N
θvector, obtaining a length is N
θn
φvector, be expressed as
The corresponding vertical angle of pitch of l paths and horizontal azimuth is estimated as
Wherein, n
θ=mod (s, N
θ), n
φ=(s-n
θ)/N
θ+ 1.
The set that estimates successively the incident angle (the vertical angle of pitch and horizontal azimuth) of L paths is
Claims (5)
1. the MIMO-OFDM system channel estimation method based on compressed sensing, is characterized in that, the method comprises the following steps:
A) utilize criterion of least squares, obtain the initial estimate of each pilot sub-carrier upper signal channel frequency domain response vector;
B) utilize the initial estimate of the channel frequency domain response vector that described step obtains in a) in the sparse property of time delay domain, based on compressive sensing theory, estimate the estimated value of the time delay in each path of channel and the channel time domain response vector in each path;
C) utilize described step b) in the estimated value of the channel time domain response vector estimated in the sparse property in two dimension angular territory, based on compressive sensing theory, estimate the incidence angle in each path of channel, described incidence angle comprises the vertical angle of pitch and horizontal azimuth;
D) utilize initial estimate, the step b of the channel frequency domain response vector that described step obtains in a)) in time delay, the step c in each path of estimating) incidence angle in each path of estimating, adopt criterion of least squares, estimate the gain coefficient in each path of channel;
E) by described step b) in time delay, the step c in each path of estimating) incidence angle and the steps d in each path of estimating) and in the gain coefficient in estimated each path, the multipath channel models of substitution channel frequency domain response, obtains the estimated value of the channel frequency domain response on all subcarriers and antenna.
2. a kind of MIMO-OFDM system channel estimation method based on compressed sensing according to claim 1, is characterized in that described step b) idiographic flow be:
First, by channel path time delay 0 to uniform discrete between maximum path time delay, obtain minimum time delay interval and be Δ τ and discrete after path number be N
τ, the Fourier transform relation of utilizing channel frequency domain response and channel time domain to respond, is the product of base vector and channel time domain matrix by the channel frequency domain response vector representation on each subcarrier;
Secondly, the initial estimate of the channel frequency domain response vector on all pilot sub-carriers is arranged in to a matrix in order, is designated as channel frequency domain matrix; All base vectors are arranged in to a matrix according to identical order, are designated as projection matrix; Thereby by channel frequency domain Matrix Estimation value representation be:
Wherein,
represent channel frequency domain Matrix Estimation value, Ξ is projection matrix, and A is channel time domain matrix, and W is noise matrix;
Finally, according to formula (12), adopt compressive sensing theory to solve channel time domain matrix, the channel time domain matrix solving is N altogether
τoK, find non-zero row vector wherein, the number L of described non-zero row vector is the channel path sum estimating, l non-zero row vector is the estimated value of the channel time domain response vector in l path, l the corresponding path delay in path is (m-1) Δ τ, wherein, l is the sequence number of non-zero row vector, and m is the channel time domain matrix line number m at l non-zero row vector place.
3. a kind of MIMO-OFDM system channel estimation method based on compressed sensing according to claim 1 and 2, is characterized in that described step c) idiographic flow be:
First, by vertical angle of pitch uniform discrete, obtain minimum angles and be spaced apart Δ θ, the angle number after discrete is N
θ, by horizontal azimuth uniform discrete, obtaining minimum angles and be spaced apart Δ φ, the angle number after discrete is N
φ, utilize the two-dimensional Fourier transform relation between channel time domain response and path incidence angle, the channel time domain response vector estimated value table in each path is shown:
Wherein,
be the channel time domain response vector estimated value in l path, P is projection matrix, γ
lbe the path complex gain vector in l path, w is noise vector;
Then calculate respectively in accordance with the following methods the incidence angle in each path: for l path, utilize formula (18), adopt compressive sensing theory to solve the estimated value that obtains path complex gain vector, according to the peaked position of element in the complex gain vector of path, calculate the incidence angle in l path, described incidence angle comprises the vertical angle of pitch and horizontal azimuth.
4. a kind of MIMO-OFDM system channel estimation method based on compressed sensing according to claim 3, is characterized in that described step c) in, every row element arrangement mode and path complex gain vector γ in described projection matrix P
lmiddle arrangement of elements mode is identical, is according to the order sequence of the Kronecker product of the vertical angle of pitch and horizontal azimuth institute respective path and obtains;
The method of calculating the incidence angle in l path is: according to the peaked position s of element in the complex gain vector of path, the vertical angle of pitch that obtains l path is
horizontal azimuth is
wherein, n
φ=mod (s, N
φ), n
θ=(s-n
φ)/N
φ+ 1, mod (x, y) represents that x is to y complementation.
5. a kind of MIMO-OFDM system channel estimation method based on compressed sensing according to claim 3, is characterized in that described step c) in, every row element arrangement mode and path complex gain vector γ in described projection matrix P
lmiddle arrangement of elements mode is identical, is according to the order sequence of the Kronecker product of horizontal azimuth and vertical angle of pitch institute respective path and obtains;
The method of calculating the incidence angle in l path is: according to the peaked position s of element in the complex gain vector of path, the vertical angle of pitch that obtains l path is
horizontal azimuth is
wherein, n
θ=mod (s, N
θ), n
φ=(s-n
θ)/N
θ+ 1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410312859.6A CN104052691B (en) | 2014-07-02 | 2014-07-02 | MIMO-OFDM system channel estimation method based on compressed sensing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410312859.6A CN104052691B (en) | 2014-07-02 | 2014-07-02 | MIMO-OFDM system channel estimation method based on compressed sensing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104052691A true CN104052691A (en) | 2014-09-17 |
CN104052691B CN104052691B (en) | 2017-02-15 |
Family
ID=51505071
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410312859.6A Active CN104052691B (en) | 2014-07-02 | 2014-07-02 | MIMO-OFDM system channel estimation method based on compressed sensing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104052691B (en) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104539335A (en) * | 2014-12-24 | 2015-04-22 | 无锡北邮感知技术产业研究院有限公司 | Limiting feedback method and device for large-scale antenna system |
CN104869086A (en) * | 2015-05-27 | 2015-08-26 | 东南大学 | MIMO-OFDM communication system downlink channel estimation method based on two-dimensional compression perception and device thereof |
CN104954310A (en) * | 2015-04-24 | 2015-09-30 | 清华大学 | Impulse noise eliminating method and device |
CN106453190A (en) * | 2016-10-10 | 2017-02-22 | 东南大学 | Multi-user sub-carrier index modulation orthogonal frequency-division multiplexing (SIM-OFDM) transmission method |
CN106487432A (en) * | 2016-10-28 | 2017-03-08 | 天津大学 | Based on compressed sensing Massive MIMO safe transmission method |
CN107276924A (en) * | 2016-04-08 | 2017-10-20 | 北京信威通信技术股份有限公司 | The channel estimation methods and device of extensive multiple-input and multiple-output |
CN107872416A (en) * | 2016-09-27 | 2018-04-03 | ***通信有限公司研究院 | The processing method and baseband processing unit of baseband signal |
CN105610477B (en) * | 2016-01-28 | 2018-06-19 | 西南交通大学 | Multiple-input multiple-output system enhancing method of signal multiplexing based on compressed sensing |
CN108242943A (en) * | 2016-12-23 | 2018-07-03 | 上海诺基亚贝尔股份有限公司 | The method and apparatus of precoding is used in communication |
CN108881076A (en) * | 2018-06-30 | 2018-11-23 | 宜春学院 | A kind of compressed sensing based MIMO-FBMC/OQAM system channel estimation method |
CN109005133A (en) * | 2018-07-12 | 2018-12-14 | 南京邮电大学 | Double Sparse multi-path channel models and the channel estimation methods based on this model |
CN109104388A (en) * | 2017-06-20 | 2018-12-28 | 希捷科技有限公司 | The devices, systems, and methods adaptive for regularization parameter |
WO2019010714A1 (en) * | 2017-07-14 | 2019-01-17 | 华为技术有限公司 | Beamforming method and device |
CN107465636B (en) * | 2017-08-21 | 2020-05-01 | 清华大学 | Channel estimation method of millimeter wave large-scale array space-frequency double-broadband system |
CN111385006A (en) * | 2018-12-27 | 2020-07-07 | 财团法人工业技术研究院 | Millimeter wave channel estimation method |
CN111988046A (en) * | 2020-08-19 | 2020-11-24 | 中国电子科技集团公司第三十六研究所 | Nonlinear system correction method and device and electronic equipment |
WO2021223211A1 (en) * | 2020-05-08 | 2021-11-11 | Qualcomm Incorporated | Signaling design for type ii csi-rs spatial domain and frequency domain basis configuration |
CN114039639A (en) * | 2021-11-04 | 2022-02-11 | 电子科技大学 | Combined channel estimation and user detection method in large connection scene |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110261905A1 (en) * | 2010-04-27 | 2011-10-27 | Qualcomm Incorporated | Compressed sensing channel estimation in ofdm communication systems |
CN102891823A (en) * | 2012-10-17 | 2013-01-23 | 东南大学 | Channel estimation method of orthogonal frequency division multi-address mobile communication system |
-
2014
- 2014-07-02 CN CN201410312859.6A patent/CN104052691B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110261905A1 (en) * | 2010-04-27 | 2011-10-27 | Qualcomm Incorporated | Compressed sensing channel estimation in ofdm communication systems |
CN102891823A (en) * | 2012-10-17 | 2013-01-23 | 东南大学 | Channel estimation method of orthogonal frequency division multi-address mobile communication system |
Non-Patent Citations (1)
Title |
---|
GEORG TAUBöCK等: ""Compressive Estimation of Doubly Selective Channels in Multicarrier Systems: Leakage Effects and Sparsity-Enhancing Processing"", 《IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING》 * |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104539335B (en) * | 2014-12-24 | 2017-11-03 | 无锡北邮感知技术产业研究院有限公司 | A kind of limited feedback method and device of extensive antenna system |
CN104539335A (en) * | 2014-12-24 | 2015-04-22 | 无锡北邮感知技术产业研究院有限公司 | Limiting feedback method and device for large-scale antenna system |
CN104954310A (en) * | 2015-04-24 | 2015-09-30 | 清华大学 | Impulse noise eliminating method and device |
CN104954310B (en) * | 2015-04-24 | 2018-12-11 | 清华大学 | Impulse noise removing method and device |
CN104869086A (en) * | 2015-05-27 | 2015-08-26 | 东南大学 | MIMO-OFDM communication system downlink channel estimation method based on two-dimensional compression perception and device thereof |
CN104869086B (en) * | 2015-05-27 | 2017-11-14 | 东南大学 | MIMO ofdm communication system down channels method of estimation, device based on two dimensional compaction perception |
CN105610477B (en) * | 2016-01-28 | 2018-06-19 | 西南交通大学 | Multiple-input multiple-output system enhancing method of signal multiplexing based on compressed sensing |
CN107276924A (en) * | 2016-04-08 | 2017-10-20 | 北京信威通信技术股份有限公司 | The channel estimation methods and device of extensive multiple-input and multiple-output |
CN107872416A (en) * | 2016-09-27 | 2018-04-03 | ***通信有限公司研究院 | The processing method and baseband processing unit of baseband signal |
CN106453190A (en) * | 2016-10-10 | 2017-02-22 | 东南大学 | Multi-user sub-carrier index modulation orthogonal frequency-division multiplexing (SIM-OFDM) transmission method |
CN106487432A (en) * | 2016-10-28 | 2017-03-08 | 天津大学 | Based on compressed sensing Massive MIMO safe transmission method |
CN106487432B (en) * | 2016-10-28 | 2020-01-03 | 天津大学 | Massive MIMO (multiple input multiple output) safe transmission method based on compressed sensing |
CN108242943A (en) * | 2016-12-23 | 2018-07-03 | 上海诺基亚贝尔股份有限公司 | The method and apparatus of precoding is used in communication |
CN109104388A (en) * | 2017-06-20 | 2018-12-28 | 希捷科技有限公司 | The devices, systems, and methods adaptive for regularization parameter |
CN109104388B (en) * | 2017-06-20 | 2021-06-11 | 希捷科技有限公司 | Apparatus, system, and method for regularization parameter adaptation |
CN110870216B (en) * | 2017-07-14 | 2021-06-15 | 华为技术有限公司 | Beam forming method and device |
CN110870216A (en) * | 2017-07-14 | 2020-03-06 | 华为技术有限公司 | Beam forming method and device |
WO2019010714A1 (en) * | 2017-07-14 | 2019-01-17 | 华为技术有限公司 | Beamforming method and device |
US11265054B2 (en) | 2017-07-14 | 2022-03-01 | Huawei Technologies Co., Ltd. | Beamforming method and device |
CN107465636B (en) * | 2017-08-21 | 2020-05-01 | 清华大学 | Channel estimation method of millimeter wave large-scale array space-frequency double-broadband system |
CN108881076A (en) * | 2018-06-30 | 2018-11-23 | 宜春学院 | A kind of compressed sensing based MIMO-FBMC/OQAM system channel estimation method |
CN108881076B (en) * | 2018-06-30 | 2021-05-04 | 宜春学院 | MIMO-FBMC/OQAM system channel estimation method based on compressed sensing |
CN109005133A (en) * | 2018-07-12 | 2018-12-14 | 南京邮电大学 | Double Sparse multi-path channel models and the channel estimation methods based on this model |
CN111385006A (en) * | 2018-12-27 | 2020-07-07 | 财团法人工业技术研究院 | Millimeter wave channel estimation method |
WO2021223211A1 (en) * | 2020-05-08 | 2021-11-11 | Qualcomm Incorporated | Signaling design for type ii csi-rs spatial domain and frequency domain basis configuration |
CN111988046A (en) * | 2020-08-19 | 2020-11-24 | 中国电子科技集团公司第三十六研究所 | Nonlinear system correction method and device and electronic equipment |
CN114039639A (en) * | 2021-11-04 | 2022-02-11 | 电子科技大学 | Combined channel estimation and user detection method in large connection scene |
CN114039639B (en) * | 2021-11-04 | 2023-04-25 | 电子科技大学 | Combined channel estimation and user detection method under large connection scene |
Also Published As
Publication number | Publication date |
---|---|
CN104052691B (en) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104052691A (en) | MIMO-OFDM system channel estimation method based on compressed sensing | |
US8594221B2 (en) | Model-based channel estimator for correlated fading channels and channel estimation method thereof | |
US7385617B2 (en) | Methods for multi-user broadband wireless channel estimation | |
CN107255793B (en) | Array direction finding method and device for broadband OFDM communication signals | |
CN110213185B (en) | Three-dimensional channel parameter estimation method based on atomic norm minimization | |
CN109005133B (en) | Double-sparse multi-path channel model and channel estimation method based on model | |
CN103763223B (en) | Sparse MIMO-OFDM channel estimation method based on space-time correlation of channel | |
CN104869086B (en) | MIMO ofdm communication system down channels method of estimation, device based on two dimensional compaction perception | |
CN102833194B (en) | OFDM (orthogonal frequency division multiplexing) channel estimation method based on symmetrical basis expansion model for quick time-varying channel | |
CN107465636B (en) | Channel estimation method of millimeter wave large-scale array space-frequency double-broadband system | |
CN106506415B (en) | A kind of method of multi-user MIMO-OFDM system channel estimation | |
Zhang et al. | Channel estimation and training design for hybrid multi-carrier mmwave massive MIMO systems: The beamspace ESPRIT approach | |
US11063724B1 (en) | Reduced channel-sounding in MU-MIMO WLANS | |
CN108881076A (en) | A kind of compressed sensing based MIMO-FBMC/OQAM system channel estimation method | |
CN106233685B (en) | The method of hybrid analog-digital simulation digital precode for extensive mimo system | |
Shafin et al. | DoA estimation and RMSE characterization for 3D massive-MIMO/FD-MIMO OFDM system | |
CN103997397A (en) | Space-domain non-uniform pilot design method | |
Uwaechia et al. | Compressed channel estimation for massive MIMO-OFDM systems over doubly selective channels | |
CN101267409A (en) | A MIMO-OFDM dual selective channel tracking method | |
CN114143896A (en) | Large-scale MIMO cross-frequency cooperation robust transmission method | |
Bhoyar et al. | Leaky least mean square (LLMS) algorithm for channel estimation in BPSK-QPSK-PSK MIMO-OFDM system | |
CN114553640B (en) | Cross-frequency-band statistical channel state information estimation method in multi-frequency-band large-scale MIMO system | |
Adeogun et al. | Novel algorithm for prediction of wideband mobile MIMO wireless channels | |
CN113037344A (en) | Multi-user large-scale MIMO multi-source image transmission interference elimination method | |
Wild et al. | Multi-antenna OFDM channel feedback compression exploiting sparsity |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |