CN101582703B - Method and device for signal channel estimation denoising post-treatment in multiantenna system - Google Patents

Method and device for signal channel estimation denoising post-treatment in multiantenna system Download PDF

Info

Publication number
CN101582703B
CN101582703B CN200810111621A CN200810111621A CN101582703B CN 101582703 B CN101582703 B CN 101582703B CN 200810111621 A CN200810111621 A CN 200810111621A CN 200810111621 A CN200810111621 A CN 200810111621A CN 101582703 B CN101582703 B CN 101582703B
Authority
CN
China
Prior art keywords
channel estimation
user
figuration
estimation results
vector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN200810111621A
Other languages
Chinese (zh)
Other versions
CN101582703A (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.)
China Academy of Telecommunications Technology CATT
Datang Mobile Communications Equipment Co Ltd
Original Assignee
China Academy of Telecommunications Technology CATT
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 China Academy of Telecommunications Technology CATT filed Critical China Academy of Telecommunications Technology CATT
Priority to CN200810111621A priority Critical patent/CN101582703B/en
Publication of CN101582703A publication Critical patent/CN101582703A/en
Application granted granted Critical
Publication of CN101582703B publication Critical patent/CN101582703B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention discloses a method for signal channel estimation denoising post-treatment in a multiantenna system, which comprises the following steps: obtaining a signal channel estimation result * ofthe multiantenna system, wherein the signal channel estimation result * comprises signal channel estimation results * of K users, namely, *=[ *1, *2, ellipsis, and *K], and k=1, 2, ellipsis, and k; obtaining a beamforming combination result *K of the multiantenna signal estimation of a user k and sequentially arranging the beamforming combination results of the users to obtain a row vector *=[*1,*2, ellipsis, and *K]; denoising the row vector * to obtain a denoised signal estimation result *; and according to a tap position set to zero in the denoised signal estimation result *, setting a corresponding position in the signal estimation result * to zero and outputting the signal estimation obtained after zero setting. The invention also discloses a device for signal channel estimation denoising post-treatment in the multiantenna system. The proposal of the invention can fully use the denoising capacity of multiple antennae and effectively improve the accuracy of the signal channel estimation denoising post-treatment.

Description

A kind of signal channel estimation denoising post-treatment method and apparatus of multiaerial system
Technical field
The present invention relates to the channel estimation technique field, particularly a kind of signal channel estimation denoising post-treatment method and apparatus of multiaerial system.
Background technology
In a lot of communication systems, in order to improve transmission reliability, the channel impulse response of the channel that need obtain passing through in the signals transmission, i.e. channel estimating.For example, in the TD-SCDMA system, in order to realize joint-detection, wave beam forming or other measuring process, channel estimation process is essential.In a lot of communication systems, channel estimating often realizes through the mode of sending training sequence.In data detection process, can suppress noise and interference targetedly by channel information owing to knowing, improve data detection performance.
But the training itself that is used as channel estimating also can receive the influence of noise and interference, makes channel estimating not accurate enough, in addition, much adopts in the system of blind Channel Estimation, because channel estimation methods is limit, channel estimating is also very inaccurate.This will be to process deleterious impact such as the Sequence Detection of later use channel estimating, wave beam forming, measurements.
For this reason, after obtaining channel estimating, often need carry out some special noise reduction process, be called signal channel estimation denoising post-treatment.A kind of method relatively more commonly used is exactly the noise gate reprocessing.The main contents of noise gate reprocessing comprise: estimate the noise power in the channel estimating; Set a thresholding according to estimated noise power; The channel estimation value that is lower than this thresholding is thought noise, and only keeps the channel estimating tap (directly) greater than this thresholding.For multiaerial system, often each antenna independently calculates channel estimation results, and each antenna carries out this noise gate noise reduction process respectively; Perhaps at first utilize a plurality of antennas to estimate respectively to obtain averaging after the noise power; Obtain an average noise power; Then power and average by antenna is asked in each tap of the channel estimating that obtains on each antenna respectively, the result after average is carried out threshold processing.Describe in detail through mathematic(al) representation below:
If antenna element k aThe channel impulse response that last estimation obtains can be written as:
h ‾ ( k a ) = [ h k a ( 1 ) , h k a ( 2 ) , . . . , h k a ( L ) ] T , k a = 1 , . . . , K a - - - ( 1 )
Here L is a training sequence length.K aBe the antenna element number.
At first, through the channel estimation results on each antenna, estimate to obtain the noise σ on each antenna Ka 2, k a=1...K aAnd it asked on average obtain average noise:
σ ‾ 2 = 1 K a Σ k a = 1 K a σ k a 2 - - - ( 2 )
Then to obtain on each antenna the channel estimating structure ask power with, promptly obtain:
P ( k a ) = [ | h k a ( 1 ) | 2 , | h k a ( 2 ) | 2 , . . . , | h k a ( L ) | 2 ] T , k a = 1 , . . . , K a - - - ( 3 )
P ‾ = 1 K a Σ k a = 1 K a p ( k a ) - - - ( 4 )
According to sequence each antenna in the formula (1) is carried out threshold processing then; Promptly set through
Figure S2008101116211D00026
thresholding
Figure S2008101116211D00027
in
Figure S2008101116211D00028
less than the tap position of this thresholding, to tap zero setting corresponding in the formula (1).
Existing signal channel estimation denoising post-treatment method is under the high s/n ratio condition, and the estimation of noise is a kind of efficient ways more accurately under the condition.But under the low signal-to-noise ratio condition, some footpath will be submerged in the noise owing to power is less, is abandoned thereby be taken as noise, will cause channel estimating inaccurate like this.In addition, for the situation that many antennas receive, simple here power addition, the not noise reduction capability of good use smart antenna.
Summary of the invention
In view of this, the objective of the invention is to, propose a kind of signal channel estimation denoising post-treatment method and apparatus of multiaerial system, can utilize the noise reduction capability of smart antenna that channel estimation results is carried out noise reduction, improve the accuracy of channel estimating.
The signal channel estimation denoising post-treatment method of the multiaerial system that the embodiment of the invention proposes comprises the steps:
Try to achieve the channel estimation results of multiaerial system
Figure S2008101116211D00029
Said channel estimation results
Figure S2008101116211D000210
The channel estimation results that comprises K user
Figure S2008101116211D00031
Promptly h ^ = [ h ^ 1 , h ^ 2 , . . . , h ^ K ] , k = 1,2 , · · · , K ;
Try to achieve the figuration amalgamation result of the multi-antenna channel estimation of user k
Figure S2008101116211D00033
Each user's figuration amalgamation result is arranged in order, obtains a row vector h ~ = [ h ~ 1 , h ~ 2 , . . . , h ~ K ] ;
carries out noise reduction process to said row vector, obtains the channel estimation results
Figure S2008101116211D00036
behind the noise reduction
Tap position according to zero setting in the channel estimation results behind the said noise reduction
Figure S2008101116211D00037
; With corresponding also zero setting of position in the channel estimation results
Figure S2008101116211D00038
, and the channel estimation results after the output zero setting processing.
The signal channel estimation denoising post-treatment device of the multiaerial system that the embodiment of the invention proposes comprises:
Channel estimation module is used for each antenna to multiaerial system and carries out channel estimating and obtain channel estimation results that the said channel estimation results of channel estimation results
Figure S2008101116211D00039
comprises K user promptly
h ^ = [ h 1 ^ , h 2 ^ , . . . , h K ^ ] , k = 1,2 , · · · , K ;
Figuration merges module; Be used to construct figuration combined vector or the figuration that each user's multi-antenna channel estimates and merge matrix; And carry out figuration with the channel estimation results that said figuration combined vector or figuration merge each user that matrix obtains said channel estimation module and merge, and each user's amalgamation result is lined up a row vector;
The noise reduction post-processing module; The row vector that said figuration merging module is obtained carries out the noise reduction reprocessing; Obtain channel estimation results
Figure S2008101116211D000314
behind the noise reduction to
Figure S2008101116211D000315
in the tap position of zero setting; With corresponding also zero setting of position in the channel estimation results
Figure S2008101116211D000316
, and the channel estimation results after the output zero setting processing.
Can find out from above technical scheme; Merge through the multi-antenna channel estimated result being carried out figuration; According to the figuration amalgamation result channel estimation results is carried out noise reduction, this figuration merging process is equivalent to a kind of effect of space filtering of multi-beam, thereby has played noise reduction and the effect that suppresses to disturb.The present invention program can make full use of the noise reduction capability of many antennas, has promoted the accuracy of signal channel estimation denoising post-treatment preferably, thereby has effectively promoted systematic function, in the system that utilizes many antennas and channel estimating, has very important realization and is worth.
Description of drawings
Fig. 1 realizes the flow chart of the signal channel estimation denoising post-treatment of multiaerial system for the embodiment of the invention;
Fig. 2 is the waveform sketch map of the equivalent wave beam of the figuration merging in the embodiment of the invention;
Fig. 3 improves the waveform sketch map of the equivalent broad beam of back figuration merging for embodiment of the invention scheme;
Fig. 4 improves the waveform sketch map of equivalence two wave beams of back figuration merging for embodiment of the invention scheme;
Fig. 5 is the modular structure figure of the signal channel estimation denoising post-treatment device of the multiaerial system of embodiment of the invention proposition;
Fig. 6 is the map as a result that embodiment of the invention scheme and prior art are carried out emulation.
Embodiment
The present invention program's the noise reduction capability of to the effect that utilizing many antennas is carried out noise reduction process to channel estimation results, thereby finds useful footpath more accurately, the elevator system performance.The present invention program mainly is applicable to little spacing antenna array, like uniform circular array, and even linear array etc.
For making the object of the invention, technical scheme and advantage clearer, the present invention is done further to set forth in detail below in conjunction with specific embodiment.
The flow process of the signal channel estimation denoising post-treatment of embodiment of the invention realization multiaerial system is as shown in Figure 1, comprises the steps:
Step 101: try to achieve each antenna channel estimated result.The channel estimation methods that adopts can be a prior art, as carrying out channel estimating through training sequence or carrying out blind Channel Estimation etc.Channel estimation results is designated as:
Figure S2008101116211D00041
Wherein Represent the channel estimation results that each antenna element is total,
Figure S2008101116211D00043
Expression is to k user's channel estimation results, h ^ = [ h ^ 1 , h ^ 2 , . . . , h ^ K ] ;
Figure S2008101116211D00045
Represent k aThe channel estimation results of individual antenna element, h ^ = h ^ 1 h ^ 2 · · · h ^ K a ;
Figure S2008101116211D00047
Represent k aIn the channel estimation results of individual antenna element, belong to k user's component, h ^ k a = [ h ^ 1 k a , h ^ 2 k a , . . . , h ^ K k a ] .
Step 102: channel estimation results is carried out the noise reduction first time, and the noise-reduction method that is adopted here can be a prior art, like the method for mentioning in the background technology.If current each the antenna channel estimated result after the noise reduction process is:
h = h 1 h 2 · · · h K a - - - ( 5 )
And:
h k a = [ h 1 k a , h 2 k a , . . . , h K k a ] , h k k a = [ h k ( k a , 1 ) , h k ( k a , 2 ) , . . . , h k ( k a , W ) ] . - - - ( 6 )
Wherein, h k aRepresent k aChannel estimation results behind the individual antenna element noise reduction, h k k aFor belonging to k user's component in the channel estimation results behind this noise reduction.
Can each antenna channel estimated result of k user be written as:
h k = h k 1 h k 2 · · · h k k a - - - ( 7 )
Then have: h=[h 1, h 2..., h K] (8)
Here k aExpression antenna subscript, K representes total number of users.W representes that each user's window is long, and it has represented certain multidiameter delay.
Step 103: the figuration combined vector of trying to achieve each user's multi-antenna channel estimation w k = [ w 1 k , w 2 k , . . . , w K a k ] . Here the figuration combined vector w of user k kCan be according to the channel estimation results h in the formula (7) kTry to achieve.The method for solving here is same as existing beamforming algorithm, promptly at first through (7), tries to achieve the space covariance matrix of user k R k = h k h k H , Then through finding the solution R kEigenvalue of maximum characteristic of correspondence vector, said characteristic vector is said user's figuration combined vector w k, this is commonly called eigen beam figuration method (EBB algorithm), perhaps passes through R k(Direction Arrive Wave DOA) estimates, pass through arrival bearing's estimated result of user k then, and antenna obtains said user's figuration combined vector w corresponding to the array response vector of this direction to carry out arrival bearing user k, this is commonly called fixed beam figuration method.
Owing to use the system of little spacing antenna array to tend to realize that DOA estimates and descending wave beam forming, thereby can directly take in the weight coefficient of wave beam forming and be used as the figuration combined vector w here kRealize that like this amount of calculation does not increase.
Step 104: the channel estimation results
Figure S2008101116211D00061
before the noise reduction carried out figuration and merged the first time that step 101 is obtained, and was about to the channel estimation results that the figuration combined vector multiply by respective user; And it is vectorial that each user's figuration amalgamation result is lined up a row.Promptly have:
h ~ k = w k h ^ k - - - ( 9 )
h ~ = [ h ~ 1 , h ~ 2 , . . . , h ~ K ] - - - ( 10 )
Theoretical according to smart antenna, the signal to noise ratio of the merging channel estimation results in the formula (10) is estimated with respect to each antenna channel in the formula (6)
Figure S2008101116211D00064
Middle signal to noise ratio is wanted high about 10log (K a) (dB), therefore, it is can Billy accurately a lot of with the channel estimation results in the formula (6) to adopt in the formula (10) channel estimation results to carry out noise reduction process.
Step 105: the row vector to obtaining after the figuration merging carries out noise reduction process according to the single antenna noise reduction algorithm; Channel estimation results
Figure S2008101116211D00065
the noise reduction process here adopts existing noise-reduction method after obtaining noise reduction, repeats no more.
Step 106: each antenna channel estimated result before the noise reduction
Figure S2008101116211D00067
is carried out noise reduction process according to
Figure S2008101116211D00066
; The tap position of zero setting in promptly to
Figure S2008101116211D00068
, the middle also zero setting of relevant position with
Figure S2008101116211D00069
.In fact exactly will merge the tap zero setting of thinking noise in the channel estimating of back, can write down the zero setting position, the tap zero setting of correspondence position during original channel is estimated.Channel estimation results after the noise reduction process is as final channel estimation results.
Here mainly be the noise reduction and inhibition interference performance that has utilized smart antenna, promoted the signal to noise ratio of channel estimating, thereby promoted the right probability of useful grade, and then promoted systematic function.More serious for outer area interference, under the condition that perhaps signal to noise ratio is lower, the method is here compared existing method and is had improvement in performance clearly.
For step 102; Noise reduction mainly is to have promoted the accuracy that figuration merging weights calculate in the step 103 for the first time, if amount of calculation is too big, can skip; Directly utilize the original channel that obtains in the step 101 to estimate to carry out figuration and merge weights calculating, also can obtain preferable performance.
For in the step 103, the figuration combined vector can directly be taken the weight vector in down beam shaping.For wherein estimating earlier DOA, obtain the method for figuration weight vector again through DOA, then directly the system of getting is used to locate or the DOA estimated result of other functions calculates, thereby saves amount of calculation greatly.So, the present invention is actually with few amount of calculation lifting and has brought bigger systematic function to promote.
Formed a wave beam because the figuration here merges to be equivalent to like Fig. 2, wherein, disk representation space angle, from 0 ° of angle to ± 180 °, disk diameter to the size that gains of numeric representation.Curve representation is along footpath, the space of all directions, is main lobe 201 along the comparatively thick of 90 ° of directions wherein, and what other was less is secondary lobe 202.This wave beam has played certain space filtering effect, thereby has promoted signal to noise ratio.But main lobe 201 width are limited, thereby, when the environmental diffusion is very big, has a small amount of secondary lobe 202 and drop on beyond the main lobe scope, this moment, filtering also possibly given by space filtering in these useful footpaths.Thereby, can adopt following method to obtain the figuration combined vector in the step 103.Promptly calculate the broad beam figuration weight coefficient of all directions in advance; And be stored in the system; In step 103, obtain then after the DOA estimation of user k; According to arrival bearing's estimated result through user k, search the broad beam figuration weight coefficient of all directions of being preserved, get corresponding figuration weight coefficient as figuration combined vector w k
Through above-mentioned processing, in fact be exactly the form that the wave beam among Fig. 2 is become Fig. 3, wherein main lobe 301 becomes more roomy, and secondary lobe 302 is then less comparatively speaking.Through adopting the method for broad beam, increased the scope of application of the present invention, make that it also can the obtained performance gain under the big environment in angular spread angle.Here beamwidth can be according to the residing environmental selection in base station.
In addition, under the bigger condition of angular spread, possibly have a plurality of arrival bearings corresponding to each user's signal, like this, for step 103, step 104 can also be carried out according to following realization:
Step 103: utilize each cell channel estimated result, obtain figuration combined vector corresponding to each each arrival bearing of user:
w ( k ) ( k d ) = [ w k 1 , w k 2 , . . . , w k K a ] , k = 1 . . . K , k d = 1 . . . K d k - - - ( 11 )
And the structure figuration merges matrix:
A ( k ) = [ w ( k ) ( 1 ) T , w ( k ) ( 2 ) T , . . . , w ( k ) ( K d k ) T ] T , k = 1 . . . K - - - ( 12 )
Here k d = 1 . . . K d k All K user's k in the expression sub-district u dIndividual arrival bearing, K d kAll arrival bearings of k=1...K user of expression expression.
The channel estimation results that can obtain according to step 101 is found the solution and is obtained the figuration combined vector w of user k=1...K corresponding to each arrival bearing (k) (k d )
The method for solving here mainly contains two kinds:
(1) method of characteristic is promptly at first tried to achieve the space covariance matrix of user k=1...K R ( k ) = h ( k ) h ( k ) H , k = 1 . . . K , Try to achieve R then (k)All characteristic values
Figure S2008101116211D00082
And characteristic vector
Figure S2008101116211D00083
Characteristic value here and characteristic vector are corresponding one by one, and be theoretical according to smart antenna, the characteristic value here
Figure S2008101116211D00084
In maximum some values be associated with several effective arrival bearings, several minimum values are associated with noise, thereby can set a thresholding Γ, and think greater than the characteristic value of this thresholding corresponding an arrival bearing, otherwise its correspondence noise.This several features value characteristic of correspondence vector is our figuration combined vector to be asked.
Get:
Figure S2008101116211D00085
is as the figuration combined vector.In fact, the effective vectorial number when here is for the moment, promptly deteriorates to the figuration combined vector that eigen beam figuration method recited above (EBB algorithm) is tried to achieve.
(2) DOA method, this method are at first wanted effective arrival bearing of estimating user k=1...K.Be made as
Figure S2008101116211D00086
Then the pairing array response vector of each effective arrival bearing is the figuration combined vector w corresponding to each arrival bearing (k) (k d ), can be with figuration combined vector w according to formula (12) (k) (k d )Merge, construct figuration and merge matrix A (k)Here a plurality of effective arrival bearing who needs estimating user k=1...K needs the DOA method of estimation of high-resolution for the less environment of angular spread, also need try to achieve user's space covariance matrix as the one of which earlier R k = h k h k H , k = 1 . . . K . The DOA method of estimation here can be existing any DOA method of estimation, like the Bartlett spectral method, and the Capon least variance method; The MUSIC algorithm; ESPRIT algorithm etc., effectively the DOA number also can calculate through existent method, as the information standard (Akaike ' s information criterion; AIC) criterion or minimum description length (Rissanen ' s minimumdescription length, MDL) criterion etc.Obviously, when effective arrival bearing's number is one, then deteriorate to the figuration vector that fixed beam figuration method above-mentioned (GOB algorithm) obtains.
Step 104: the result
Figure S2008101116211D00088
before the noise reduction first time of the channel estimating in the formula (7) is carried out figuration merge, and each user's amalgamation result is lined up a row vector.Promptly have:
h ~ k = A ( k ) h ^ k - - - ( 9 )
h ~ = [ h ~ 1 , h ~ 2 , . . . , h ~ K ] - - - ( 10 )
The process that this figuration merges has just been accomplished a kind of effect of space filtering of multi-beam; Thereby the effect of having played noise reduction and having suppressed to disturb; The figuration merging process here is equivalent to and has utilized multi-beam figuration directional diagram as shown in Figure 4 to carry out space filtering (is example with two wave beams), and wherein main lobe 401 has two different directions.The method of formation broad beam shown in Figure 3 can be regarded as a special case of multi-beam figuration merging here.
The embodiment of the invention proposes a kind of signal channel estimation denoising post-treatment device that is applicable to multiaerial system, and its modular structure is as shown in Figure 5, comprising:
Channel estimation module 501 is used for each antenna to multiaerial system and carries out channel estimating and obtain channel estimation results
Figure S2008101116211D00091
Wherein
Figure S2008101116211D00092
Represent the channel estimation results that each antenna element is total,
Figure S2008101116211D00093
Expression is to k user's channel estimation results,
Figure S2008101116211D00094
Represent k aThe channel estimation results of individual antenna element,
Figure S2008101116211D00095
Represent k aIn the channel estimation results of individual antenna element, belong to k user's component.
Preliminary noise reduction module 502 is used for the channel estimation results that said channel estimation module 501 obtains is carried out the noise reduction first time, and the channel estimation results h that incites somebody to action behind the noise reduction first time is input to figuration merging module 503.Also can omit this module in the practical application, make amount of calculation further reduce.
Figuration merges module 503; Be used to construct figuration combined vector or the figuration that each user's multi-antenna channel estimates and merge matrix; And carry out figuration with the channel estimation results that said figuration combined vector or figuration merge each user that matrix obtains said channel estimation module 501 and merge, and each user's amalgamation result is lined up a row vector.
Noise reduction post-processing module 504; Be used for the row vector that said figuration merging module 503 obtains is carried out the noise reduction reprocessing; Channel estimation results
Figure S2008101116211D00097
carries out noise reduction process according to
Figure S2008101116211D00098
to each antenna channel estimated result before the noise reduction
Figure S2008101116211D00099
after obtaining noise reduction; The tap position of zero setting in promptly to
Figure S2008101116211D000910
, the middle also zero setting of relevant position with .Channel estimation results after the noise reduction process is as final channel estimation results
In the TD-SCDMA system, utilize as above scheme to carry out emulation, simulated conditions is following:
1) CASE3 channel condition, single sub-district 12.2K is professional, single user and 6 users.
2) 120 degree sectors, even linear array, 6 antennas, 1000 data blocks.
3) utilize original estimation result to carry out user DOA and estimate that method of estimation adopts the BartLett spectral method, promptly tries to achieve the average Rxx corresponding to W space covariance matrix of the individual tap of W (W=128/Kcell) in certain user's the window.Try to achieve the Bartlett spectrum of Rxx subsequently, think that its spectrum maximum power place is the subscriber signal arrival bearing.Here trying to achieve Bartlett spectral corner degree step-length (resolution) is 2 degree.
4) user's incident footpath angular spread obedience
Figure S2008101116211D00101
goes up the Laplace distributed random variable.
Here β=11.44.
The demodulation performance contrast that emulation obtains is as shown in Figure 6, and wherein ordinate representes to transmit Block Error Rate (BLER), and abscissa Ior/Ioc representes signal to noise ratio, and unit is a decibel (dB).It is thus clear that, utilize 6 antenna linear arrays, utilize algorithm of the present invention under amount of calculation increase situation seldom, under single user's condition, signal to noise ratio can obtain the 0.5dB gain, and under 6 user's conditions, signal to noise ratio can obtain the gain greater than 1dB.
The above is merely preferred embodiment of the present invention, not in order to restriction the present invention, all any modifications of within spirit of the present invention and principle, being done, is equal to and replaces and improvement etc., all should be included within protection scope of the present invention.

Claims (11)

1. the signal channel estimation denoising post-treatment method of a multiaerial system is characterized in that, comprises the steps:
The channel estimation results
Figure FSB00000824776800013
i.e. that the said channel estimation results of channel estimation results
Figure FSB00000824776800011
that obtains multiaerial system comprises K user is k=1; 2;, K;
The figuration amalgamation result
Figure FSB00000824776800015
that obtains the multi-antenna channel estimation of user k is arranged in order each user's figuration amalgamation result, obtains a row vector
Figure FSB00000824776800016
Figure FSB00000824776800017
carries out noise reduction process to said row vector, obtains the channel estimation results
Figure FSB00000824776800018
behind the noise reduction
Tap position according to zero setting in the channel estimation results behind the said noise reduction
Figure FSB00000824776800019
; With corresponding also zero setting of position in the channel estimation results
Figure FSB000008247768000110
, and the channel estimation results after the output zero setting processing.
2. the signal channel estimation denoising post-treatment method of multiaerial system according to claim 1; It is characterized in that the figuration amalgamation result that the multi-antenna channel of the said user of obtaining k is estimated comprising:
Channel estimation results according to user k Try to achieve the figuration combined vector w of the multi-antenna channel estimation of said user k k, and said figuration combined vector multiply by said user's channel estimation results, obtain said user's figuration amalgamation result
Figure 2008101116211100001DEST_PATH_IMAGE002
3. the signal channel estimation denoising post-treatment method of multiaerial system according to claim 2 is characterized in that, the figuration combined vector w that said multi-antenna channel of trying to achieve said user is estimated kFor:
Channel estimation results according to said user k
Figure FSB000008247768000114
Try to achieve said user's space covariance matrix
Figure FSB000008247768000115
Then through finding the solution R kEigenvalue of maximum characteristic of correspondence vector, said characteristic vector is said user's figuration combined vector w k
4. the signal channel estimation denoising post-treatment method of multiaerial system according to claim 2 is characterized in that, said channel estimation results according to user k
Figure FSB000008247768000116
Try to achieve the figuration combined vector w of the multi-antenna channel estimation of said user k kFor:
Channel estimation results according to said user k
Figure FSB000008247768000117
Try to achieve said user's space covariance matrix
Figure FSB000008247768000118
According to said space covariance matrix R kCarry out arrival bearing user DOA and estimate to obtain arrival bearing's estimated result of user k; The figuration combined vector w that obtains said user corresponding to the array response vector of this direction according to arrival bearing's estimated result and the antenna of said user k k
5. the signal channel estimation denoising post-treatment method of multiaerial system according to claim 2 is characterized in that, said channel estimation results according to user k
Figure FSB00000824776800021
Try to achieve the figuration combined vector w of said user's multi-antenna channel estimation kComprise:
Calculate and preserve the broad beam figuration weight coefficient of all directions in advance;
Channel estimation results according to said user k
Figure FSB00000824776800022
Try to achieve said user's space covariance matrix
Figure FSB00000824776800023
According to said space covariance matrix R kCarry out arrival bearing user DOA and estimate to obtain arrival bearing's estimated result of user k, the arrival bearing's estimated result according to user k checks in the figuration combined vector w of corresponding figuration weight coefficient as said user then k
6. the signal channel estimation denoising post-treatment method of multiaerial system according to claim 1; It is characterized in that the said channel estimation results
Figure FSB00000824776800024
that obtains multiaerial system further comprises afterwards:
Figure FSB00000824776800025
carries out the noise reduction first time to said channel estimation results, obtains the channel estimation results h after the noise reduction process for the first time;
The figuration amalgamation result
Figure FSB00000824776800026
that the multi-antenna channel of the then said user of obtaining k is estimated is:
The channel estimation results h of the first time after the noise reduction process according to user k k, try to achieve the figuration amalgamation result that the multi-antenna channel of said user k is estimated
Figure FSB00000824776800027
7. the signal channel estimation denoising post-treatment method of multiaerial system according to claim 1; It is characterized in that the figuration amalgamation result that the multi-antenna channel of the said user of obtaining k is estimated comprising:
Channel estimation results according to said user k
Figure FSB00000824776800029
Try to achieve said user's space covariance matrix Try to achieve R kAll characteristic values and each characteristic value characteristic of correspondence vector;
With said characteristic value and predefined thresholding relatively, will be greater than the characteristic value characteristic of correspondence vector of thresholding as the figuration combined vector;
Merge matrix according to said figuration combined vector structure figuration; The channel estimation results
Figure FSB000008247768000211
that figuration is merged matrix and user k multiplies each other, and obtains the figuration amalgamation result
Figure FSB000008247768000212
of the multi-antenna channel estimation of said user k
8. the signal channel estimation denoising post-treatment method of multiaerial system according to claim 1; It is characterized in that the figuration amalgamation result
Figure FSB00000824776800031
that the multi-antenna channel of the said user of obtaining k is estimated is:
To be used for user k is carried out the figuration combined vector w of the weight vector of down beam shaping as said user's multi-antenna channel estimation k, and said figuration combined vector multiply by said user's channel estimation results, obtain said user's figuration amalgamation result
Figure FSB00000824776800032
9. the signal channel estimation denoising post-treatment method of multiaerial system according to claim 1; It is characterized in that the figuration amalgamation result
Figure FSB00000824776800033
that the multi-antenna channel of the said user of obtaining k is estimated comprising:
A plurality of effective arrival bearing of estimating user k=1...K; The pairing array response vector of each effective arrival bearing; Be figuration combined vector corresponding to each arrival bearing; Merge matrix according to said figuration combined vector structure figuration; The channel estimation results
Figure FSB00000824776800034
that figuration is merged matrix and user k multiplies each other, and obtains the figuration amalgamation result
Figure FSB00000824776800035
of the multi-antenna channel estimation of said user k
10. the signal channel estimation denoising post-treatment device of a multiaerial system; It is characterized in that; Comprise: channel estimation module; Being used for each antenna to multiaerial system carries out channel estimating and obtains channel estimation results
Figure FSB00000824776800038
i.e. that the said channel estimation results of channel estimation results
Figure FSB00000824776800036
Figure FSB00000824776800037
comprises K user
Figure FSB00000824776800039
k=1; 2;, K;
Figuration merges module; Be used to construct figuration combined vector or the figuration that each user's multi-antenna channel estimates and merge matrix; And carry out figuration with the channel estimation results
Figure FSB000008247768000310
that said figuration combined vector or figuration merge each user that matrix obtains said channel estimation module and merge, and each user's amalgamation result is lined up a row vector;
The noise reduction post-processing module; The row vector that said figuration merging module is obtained carries out the noise reduction reprocessing; Obtain channel estimation results
Figure FSB000008247768000311
behind the noise reduction to
Figure FSB000008247768000312
in the tap position of zero setting; With corresponding also zero setting of position in the channel estimation results
Figure FSB000008247768000313
, and the channel estimation results after the output zero setting processing.
11. the signal channel estimation denoising post-treatment device of multiaerial system according to claim 10 is characterized in that said device further comprises
Preliminary noise reduction module, the channel estimation results that is used for said channel estimation module is obtained carries out first noise reduction, and the channel estimation results behind the first noise reduction is input to figuration merging module.
CN200810111621A 2008-05-15 2008-05-15 Method and device for signal channel estimation denoising post-treatment in multiantenna system Active CN101582703B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200810111621A CN101582703B (en) 2008-05-15 2008-05-15 Method and device for signal channel estimation denoising post-treatment in multiantenna system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200810111621A CN101582703B (en) 2008-05-15 2008-05-15 Method and device for signal channel estimation denoising post-treatment in multiantenna system

Publications (2)

Publication Number Publication Date
CN101582703A CN101582703A (en) 2009-11-18
CN101582703B true CN101582703B (en) 2012-09-12

Family

ID=41364696

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810111621A Active CN101582703B (en) 2008-05-15 2008-05-15 Method and device for signal channel estimation denoising post-treatment in multiantenna system

Country Status (1)

Country Link
CN (1) CN101582703B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102638425B (en) * 2012-03-13 2015-08-12 京信通信系统(中国)有限公司 The channel estimation methods of long term evolution downlink channel and device
CN108282423B (en) * 2018-01-03 2021-01-19 深圳市极致汇仪科技有限公司 Noise reduction method and system for MIMO channel estimation

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1855770A (en) * 2005-04-28 2006-11-01 上海原动力通信科技有限公司 Beam shaping method for inhibiting interferes
CN1870610A (en) * 2005-05-25 2006-11-29 上海原动力通信科技有限公司 Channel estimating method of array antenna

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1855770A (en) * 2005-04-28 2006-11-01 上海原动力通信科技有限公司 Beam shaping method for inhibiting interferes
CN1870610A (en) * 2005-05-25 2006-11-29 上海原动力通信科技有限公司 Channel estimating method of array antenna

Also Published As

Publication number Publication date
CN101582703A (en) 2009-11-18

Similar Documents

Publication Publication Date Title
CN101482610B (en) Self-adapting special interference restraint technology for phased array radar
KR101768587B1 (en) Covariance matrix estimation method for reducing nonstationary clutter and heterogeneity clutter
CN106443721A (en) Rapid convergence type anti-interference method of satellite navigation antenna array
CN102175989A (en) Method for measuring incoherently distributed signal two-dimensional DOA (direction of arrival)
CN107092007A (en) A kind of Wave arrival direction estimating method of virtual second order array extension
CN102830387A (en) Data preprocessing based covariance matrix orthogonalization wave-beam forming method
CN105137437A (en) Target detection method based on spatial domain phase variance weighting
CN111948634A (en) Target detection method and device based on covariance matrix reconstruction under interference condition
CN101436893A (en) Beam size enlargement apparatus and method for restraining interference of intelligent antenna
CN106443594A (en) Radar antenna array steady beam forming method based on sparse constraint
CN107124216A (en) A kind of Capon robust adaptive beamforming method and system for array error
CN104931942B (en) A kind of method for eliminating radar multi-path jamming
CN105204008A (en) Adaptive antenna wave beam forming nulling widening method based on covariance matrix extension
CN103941234B (en) A kind of Beam-former method for designing for multipath signal
CN110824414A (en) Device and method for estimating angle of arrival
CN104360355A (en) Anti-interference method and device
CN107290732A (en) A kind of single base MIMO radar direction-finding method of quantum huge explosion
CN109861933A (en) A kind of millimeter wave mimo channel estimation method based on MUSIC algorithm and precoding
CN101582703B (en) Method and device for signal channel estimation denoising post-treatment in multiantenna system
CN101615943B (en) Method for estimating elevation angle of intelligent antenna multi-subarray system
CN110212951B (en) Large-scale MIMO channel estimation method based on Butler matrix
CN107167776A (en) The adaptive beam-forming algorithm compensated based on subspace
CN104716982A (en) Robust anti-interference processing method and device of spread frequency system
CN101588190B (en) Method and device for channel estimation, and corresponding base station system
CN101594322B (en) Multi-antenna channel estimation method and device as well as corresponding base station system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: INST OF TELECOMMUNICATION SCIENCE AND TECHNOLGOY

Free format text: FORMER OWNER: DATANG MOBILE COMMUNICATION EQUIPMENT CO., LTD.

Effective date: 20110621

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100083 NO. 29, XUEYUAN ROAD, HAIDIAN DISTRICT, BEIJING TO: 100191 NO. 40, XUEYUAN ROAD, HAIDIAN DISTRICT, BEIJING

TA01 Transfer of patent application right

Effective date of registration: 20110621

Address after: 100191 Haidian District, Xueyuan Road, No. 40,

Applicant after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Address before: 100083 Haidian District, Xueyuan Road, No. 29,

Applicant before: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.

C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee after: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

TR01 Transfer of patent right

Effective date of registration: 20210603

Address after: 100085 1st floor, building 1, yard 5, Shangdi East Road, Haidian District, Beijing

Patentee after: DATANG MOBILE COMMUNICATIONS EQUIPMENT Co.,Ltd.

Address before: 100191 No. 40, Haidian District, Beijing, Xueyuan Road

Patentee before: CHINA ACADEMY OF TELECOMMUNICATIONS TECHNOLOGY

TR01 Transfer of patent right