CN109150783A - A kind of channel estimation methods and device - Google Patents
A kind of channel estimation methods and device Download PDFInfo
- Publication number
- CN109150783A CN109150783A CN201710508459.6A CN201710508459A CN109150783A CN 109150783 A CN109150783 A CN 109150783A CN 201710508459 A CN201710508459 A CN 201710508459A CN 109150783 A CN109150783 A CN 109150783A
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- virtual carrier
- channel estimation
- value
- channel
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
- H04L27/2695—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2668—Details of algorithms
- H04L27/2673—Details of algorithms characterised by synchronisation parameters
- H04L27/2675—Pilot or known symbols
Abstract
This application discloses a kind of channel estimation methods and devices to avoid the generation of time domain mirror image diameter when using DFT channel estimation under multi-user installation scene, improves channel estimation accuracy.The frequency domain received data of symbol where a kind of channel estimation methods provided by the present application comprise determining that the demodulated reference signal in each time slot, and the basic sequence in sequence is sent using the demodulated reference signal of multi-user installation, initial channel estimation is carried out, the channel estimation value of frequency domain is obtained；Virtual carrier, and the increased virtual carrier in two sides of the channel estimation value in the frequency domain are determined using virtual carrier generator matrix；The wherein virtual carrier generator matrix is that the positional relationship of the channel estimation value according to virtual carrier Yu the frequency domain generates；Using the channel estimation value for having the frequency domain of virtual carrier is increased in two sides, the joint time domain channel estimated value of multi-user installation is determined.
Description
Technical field
This application involves field of communication technology more particularly to a kind of channel estimation methods and device.
Background technique
The base station long term evolution (Long Term Evolution, LTE) physically ink Shared Channel (Physical
Uplink Shared Channel, PUSCH) receive symbol level equalizer when need using uplink (Up Link, UL) demodulation ginseng
Examine the channel estimation of signal (Demodulation Reference Signal, DMRS) progress.In multi-user's multiple-input, multiple-output
Multiple user equipment (User when (Most In Most Out, MIMO) (Multiple Users MIMO, MU-MIMO) scene
Equipment, UE) it is multiplexed same running time-frequency resource, the DMRS sequence of different UE is distinguish by frequency domain phase rotation, is being passed
In the channel estimation based on discrete Fourier transform (Discrete Fourier Transform, DFT) method of system, receiver
There are the characteristics that cyclic shift progress UE differentiation in time domain using the CIR of multiple UE and is respectively completed channel estimation.
In uplink assignment resource block (Resource Block, RB) number less (narrowband), the channel estimation based on DFT exists
Time resolution reduce, time-domain windowed noise suppression decreased effectiveness, at the same the non-integer point sampling problem of time domain multipath at narrowband more
To be obvious, caused by the leakage of time domain multipath energy it is also more serious.
To solve this problem, the prior art mainly adopts the modes for increasing virtual carrier to improve time resolution, a kind of increase
Virtual carrier way is to increase frequency domain mirror image, and concrete operations are in inverse discrete Fourier transform (Inverse Discrete
Fourier Transform, IDFT) the symmetrical mirror image of preoperative channel estimation in frequency domain increase (IDFT points increase at this time),
Time domain after the completion of IDFT completes multi-path resolved and adding window noise suppression, carries out DFT operation after Time Domain Processing, frequency domain data later is gone
Final channel estimation in frequency domain is obtained after mirror image.
In conclusion the CIR of more UE is arranged on the multiple positions of time domain in MU-MIMO scene, observed in frequency domain visible
Actual sub-carriers and virtual carrier intersection phase difference are big compared with SU-MIMO scene, are using the channel estimation side based on frequency domain mirror image
When method, there are SPA sudden phase anomalies between mirrored portion and original signal part, it will lead in frequency domain symbol that phase is discontinuous (or to be construed to
Destroy the orthogonality of IDFT), it is viewed as false diameter (mirror image diameter) occur in time domain, destroys the suppression of original multi-path resolved and adding window
It makes an uproar rule, increases the mean square error (Mean Squared Error, MSE) of the channel estimation based on DFT.To avoid this from asking
Topic, needs to carry out the Multipath searching of this UE after the CIR cyclic shift of each UE to time domain axis edge, this mode leads to each UE points
Not Wan Cheng IDFT/DFT operation, computation complexity is high.
Summary of the invention
The embodiment of the present application provides a kind of channel estimation methods and device, to use under multi-user installation scene
When DFT channel estimation, the generation of time domain mirror image diameter is avoided, improves channel estimation accuracy.
A kind of channel estimation methods provided by the embodiments of the present application include:
The frequency domain received data of symbol where determining the demodulated reference signal in each time slot, and use multi-user installation
Demodulated reference signal sends the basic sequence in sequence, carries out initial channel estimation, obtains the channel estimation value of frequency domain；
Virtual carrier, and the increased imaginary load in two sides of the channel estimation value in the frequency domain are determined using virtual carrier generator matrix
Wave；The wherein virtual carrier generator matrix is that the positional relationship of the channel estimation value according to virtual carrier Yu the frequency domain generates；
Using the channel estimation value for having the frequency domain of virtual carrier is increased in two sides, the joint time domain channel of multi-user installation is determined
Estimated value.
In this method provided by the embodiments of the present application, the generating mode of virtual carrier, i.e. virtual carrier generator matrix are improved, is
It is arranged in the arrangement mode generation of the two sides of the channel estimation value of frequency domain according to virtual carrier, is determined using virtual carrier generator matrix
Virtual carrier, and the increased virtual carrier in two sides of the channel estimation value in the frequency domain utilize the frequency for having virtual carrier in two sides increase
The channel estimation value in domain determines the joint time domain channel estimated value of multi-user installation, therefore by improving virtual carrier and original son
The phase continuity of carrier wave reduces channel estimation to avoid the generation of time domain mirror image diameter under multi-user installation scene
Error.
Optionally, the virtual carrier generator matrix W is generated using following formula:
Wherein:
V=0...K-1,12N+K...12N+2K-1
M=0...12N-1
N indicates the number for the resource block that the ascending resource of multiple user equipment multiplexings is included, and K is indicated in the frequency domain
The increased imaginary load wave number in every side of channel estimation value, p indicate that the frequency domain signal-to-noise ratio of virtual carrier application, I indicate unit diagonal matrix.
Optionally, using following formula, virtual carrier is determined using the virtual carrier generator matrix
Wherein
Wherein, describedIndicate the channel estimation value of the frequency domain,WithRespectively indicate the letter needed in the frequency domain
The increased virtual carrier in the two sides of road estimated value.
Optionally, using the channel estimation value for having the frequency domain of virtual carrier is increased in two sides, the joint of multi-user installation is determined
Time domain channel estimated value, specifically includes:
Calculate the cyclic shift value of the demodulated reference signal of each user equipmentWherein j=0...U-1, U indicate user
Number of devices；
Calculate so thatThe value range of maximum α, α are as follows: 0 to 12N-
1；
Calculated α completes the phase rotation of H, the channel estimation value H of the frequency domain after obtaining phase rotationr,
Middle H indicates that two sides increase the channel estimation value for having the frequency domain of virtual carrier,
Utilize HrDetermine the joint time domain channel estimated value of multi-user installation.
By analyzing the cyclic shift value of multiple UE in the embodiment of the present application, additive phase is introduced in frequency domain before IDFT
Rotation, so that the midpoint on time domain axis is avoided in the CIR arrangement of multi-user installation, when this operation can be arranged in avoid CIR
The channel estimation results bigger error of UE near the midpoint of domain.
Optionally, H is utilizedrThe joint time domain channel estimated value for determining multi-user installation, specifically includes:
To the HrAdding window is carried out, H is obtainedw；
Calculate HwCorresponding time domain channel shock response；
Multi-path location, adding window noise suppression are calculated separately to each user equipment in time domain channel shock response, obtained multiple
The joint time domain channel estimated value of user equipment；
Equipment for each user calculates inverse cyclic shift value, and after removing noise figure outside window, utilizes the inverse circulation
Shift value, removes the additional cyclic shift value of the time domain channel estimated value of the user equipment, and is converted to frequency domain, finally obtains this
The channel estimation in frequency domain value of user equipment.
In the embodiment of the present application, window adding in frequency domain operation is introduced before IDFT is operated, can reduce the introducing due to virtual carrier
Bring time domain energy leakage problem, it is ensured that the raising of the accuracy of channel estimation.
Optionally, equipment for each user determines inverse cyclic shift value using following formula
Wherein j=0...U-1.
A kind of channel estimating apparatus provided by the embodiments of the present application, comprising:
First unit for the frequency domain received data of symbol where determining the demodulated reference signal in each time slot, and makes
The basic sequence in sequence is sent with the demodulated reference signal of multi-user installation, initial channel estimation is carried out, obtains the channel of frequency domain
Estimated value；
Second unit, for determining virtual carrier using virtual carrier generator matrix, and in the channel estimation value of the frequency domain
The increased virtual carrier in two sides；The wherein virtual carrier generator matrix is the channel estimation value according to virtual carrier and the frequency domain
What positional relationship generated；
Third unit, for determining multi-user installation using the channel estimation value for having the frequency domain of virtual carrier is increased in two sides
Joint time domain channel estimated value.
Optionally, the second unit generates the virtual carrier generator matrix W using following formula:
Wherein:
V=0...K-1,12N+K...12N+2K-1
M=0...12N-1
N indicates the number for the resource block that the ascending resource of multiple user equipment multiplexings is included, and K is indicated in the frequency domain
The increased imaginary load wave number in every side of channel estimation value, p indicate that the frequency domain signal-to-noise ratio of virtual carrier application, I indicate unit diagonal matrix.
Optionally, the second unit uses following formula, determines virtual carrier using the virtual carrier generator matrix
Wherein
Wherein, describedIndicate the channel estimation value of the frequency domain,WithRespectively indicate the letter needed in the frequency domain
The increased virtual carrier in the two sides of road estimated value.
Optionally, the third unit, is specifically used for:
Calculate the cyclic shift value of the demodulated reference signal of each user equipmentWherein j=0...U-1, U indicate user
Number of devices；
Calculate so thatThe value range of maximum α, α are as follows: 0 to 12N-
1；
Calculated α completes the phase rotation of H, the channel estimation value H of the frequency domain after obtaining phase rotationr,
Middle H indicates that two sides increase the channel estimation value for having the frequency domain of virtual carrier,
Utilize HrDetermine the joint time domain channel estimated value of multi-user installation.
Optionally, the third unit utilizes HrThe joint time domain channel estimated value for determining multi-user installation, specifically includes:
To the HrAdding window is carried out, H is obtainedw；
Calculate HwCorresponding time domain channel shock response；
Multi-path location, adding window noise suppression are calculated separately to each user equipment in time domain channel shock response, obtained multiple
The joint time domain channel estimated value of user equipment；
Equipment for each user calculates inverse cyclic shift value, and after removing noise figure outside window, utilizes the inverse circulation
Shift value, removes the additional cyclic shift value of the time domain channel estimated value of the user equipment, and is converted to frequency domain, finally obtains this
The channel estimation in frequency domain value of user equipment.
Optionally, third unit equipment for each user determines inverse cyclic shift value using following formula
Wherein j=0...U-1.
Another embodiment of the application provides a kind of calculating equipment comprising memory and processor, wherein the storage
Device is for storing program instruction, and the processor is for calling the program instruction stored in the memory, according to the journey of acquisition
Sequence executes any of the above-described kind of method.
Another embodiment of the application provides a kind of computer storage medium, and the computer-readable recording medium storage has
Computer executable instructions, the computer executable instructions are for making the computer execute any of the above-described kind of method.
Detailed description of the invention
In order to more clearly explain the technical solutions in the embodiments of the present application, make required in being described below to embodiment
Attached drawing is briefly introduced, it should be apparent that, the drawings in the following description are only some examples of the present application, for this
For the those of ordinary skill in field, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is a kind of overall procedure schematic diagram of channel estimation methods provided by the embodiments of the present application；
Fig. 2 is a kind of idiographic flow schematic diagram of channel estimation methods provided by the embodiments of the present application；
Fig. 3 is a kind of structural schematic diagram of channel estimating apparatus provided by the embodiments of the present application；
Fig. 4 is a kind of structural schematic diagram for calculating equipment provided by the embodiments of the present application.
Specific embodiment
The embodiment of the present application provides a kind of channel estimation methods and device, to use under multi-user installation scene
When DFT channel estimation, the generation of time domain mirror image diameter is avoided, improves channel estimation accuracy.
Long term evolution (Long Term Evolution, LTE) ascending physical signal uplink shared channel (Physical
Uplink Shared Channel, PUSCH) channel estimation, at multi-user's multiple-input, multiple-output (Most In Most Out, MIMO)
Under (Multiple Users MIMO, MU-MIMO) scene, when user equipment (User Equipment, UE) scheduled distribution
When RB number is less, the conventional channel estimation methods error based on DFT is larger.Had using the method that mirror-image fashion increases virtual carrier
Conducive to reduction evaluated error, but this mode for increasing virtual carrier can bring the energy leakage of multipath, companion under multi-user scene
With there is mirror image diameter, channel estimating performance is caused to deteriorate (not having this problem under single user scene).In order to adapt to multi-user
Scene, mirror image discrete Fourier transform (Discrete Fourier Transform, DFT) channel estimation needs when MU-MIMO
Time-frequency conversion is carried out by user, time-frequency conversion number increases several times than conventional DFT channel estimation, and complexity is very high.This Shen
Please embodiment propose a kind of new improvement project, the appearance of mirror image diameter can be avoided while increasing time resolution, reduction
Original diameter energy leakage so that multiple UE of MU-MIMO channel impulse response (Channel Impulse Response,
CIR) processing can pass through inverse discrete Fourier transform (Inverse Discrete Fourier Transform, an IDFT) behaviour
It is unified after work to complete, it avoids carrying out IDFT operation to each UE, reduces computation complexity.
Specifically, for the channel estimation problems under MU-MIMO scene, the embodiment of the present application improves the generation of virtual carrier
Mode improves the phase continuity of virtual carrier and original subcarrier, the generation of time domain mirror image diameter is avoided under multi-user scene, drops
Low channel estimation errors.The cyclic shift parameter that multiple UE are analyzed in the embodiment of the present application, it is attached in frequency domain introducing before IDFT
Add phase rotation, so that the midpoint (phase discontinuously obvious region) on time domain axis is avoided in the CIR arrangement of multi-user.This behaviour
Make that the UE channel estimation results bigger error near time domain midpoint can be arranged in avoid CIR.In addition, being introduced before IDFT is operated
Window adding in frequency domain operation is reduced since the introducing bring time domain energy of virtual carrier is revealed.
True carrier wave described in the embodiment of the present application is the frequency domain position that data, pilot tone are really carried in OFDM, imaginary load
Wave is the frequency domain position for not carrying data, pilot tone.
Technical solution provided by the embodiments of the present application is illustrated with reference to the accompanying drawing.
Referring to Fig. 1, a kind of channel estimation methods provided by the embodiments of the present application include:
S101, determine demodulated reference signal in each time slot where symbol frequency domain received data, and use multi-user
The demodulated reference signal of equipment sends the basic sequence in sequence, carries out initial channel estimation, obtains the channel estimation value of frequency domain；
S102, virtual carrier is determined using virtual carrier generator matrix, and increase in the two sides of the channel estimation value of the frequency domain
Virtual carrier；The wherein virtual carrier generator matrix is the positional relationship of the channel estimation value according to virtual carrier Yu the frequency domain
It generates；
S103, using two sides increase have virtual carrier frequency domain channel estimation value, when determining the joint of multi-user installation
Domain channel estimation value.
In this method provided by the embodiments of the present application, the generating mode of virtual carrier, i.e. virtual carrier generator matrix are improved, is
It is arranged in the arrangement mode generation of the two sides of the channel estimation value of frequency domain according to virtual carrier, is determined using virtual carrier generator matrix
Virtual carrier, and the increased virtual carrier in two sides of the channel estimation value in the frequency domain utilize the frequency for having virtual carrier in two sides increase
The channel estimation value in domain determines the joint time domain channel estimated value of multi-user installation, therefore by improving virtual carrier and original son
The phase continuity of carrier wave reduces channel estimation to avoid the generation of time domain mirror image diameter under multi-user installation scene
Error.
The illustration of a specific embodiment is given below.
Referring to fig. 2, a kind of optional detailed process provided by the embodiments of the present application is as follows:
S201: under MU-MIMO scene, multiple uplink UE multiplexed resources are N number of resource block (RB), and base station side completes radio frequency
(Radio Frequency, RF) multiple antennas receives, by digital front-end (Digital Front End, DFE) or quick Fu
The frequency-region signal of DMRS symbol in each time slot is obtained after leaf transformation (Fast Fourier Transformation, FFT) processing
Yi, wherein i=0...12N-1, wherein i indicates sub-carrier indices.
S202: the basic sequence X in sequence is sent using the DMRS of more UEi, wherein i=0...12N-1 is completed based on most
The initial channel estimation of small square law (Least Square, LS) method, obtains the frequency domain representation of channel estimationSpecifically:Wherein, i=0...12N-1.
S203: defining virtual carrier arrangement mode is to increaseTwo sides, whereinFor vector expression, increaseTwo sides are
Respectively increase some elements end to end in vector, it is K that every side, which increases imaginary load wave number, and defining Q point DFT, (Q point DFT means that DFT calculates rule
Mould is Q point, i.e. FQSize be Q*Q) be transformed to FQ(m, n), wherein F indicates the canonical matrix of DFT transform, output frequency domain index
Collection is combined into m, and input time domain indexed set is combined into n, for example, output frequency domain index complete or collected works be 0,1,,, N-1, indexed set here
Closing m may be { 0,1,2,3 }.The frequency domain signal-to-noise ratio for defining virtual carrier application is p, to keep virtual carrier signal-to-noise ratio, compared with primary carrier
Height uses lower default signal-to-noise ratio (for example taking snr=0db) herein.Virtual carrier generator matrix W is calculated using following formula:
Wherein:
V=0...K-1,12N+K...12N+2K-1
M=0...12N-1
Wherein, W dimension is 2K*12N, and virtual carrier generator matrix W can be precalculated, and saves as constant look-up table.
F is standard FFT matrix, is square matrix, F12NThe as square matrix of 12N*12N, other similarly, details are not described herein.
Wherein, size (L1) indicate L1Length.
I indicates unit diagonal matrix.
As it can be seen that improving the generating mode of virtual carrier by this step, the optimization of virtual carrier generator matrix minimizes target
Are as follows: before and after increasing virtual carrier, the mean square error of original sub-carrier channels estimation is reduced, when can reduce under MU-MIMO scene
The generation of domain mirror image diameter reduces channel estimation errors.
S204: virtual carrier is generated using virtual carrier generator matrix W:WhereinRespectively willWithIt is appended toTwo sides (i.e. vectorHead and the tail), obtain:
S205: definition UE number is U, and the DMRS cyclic shift value of each UE obtainsWherein j=0...U-1；It defines attached
Add shifting value α, attempts to calculate α, so thatMaximum, the value range of α are as follows: 0 to
12N-1.And calculated α completes the phase rotation of H, obtains:
Wherein, the DMRS cyclic shift value for calculating each UE can be realized using the prior art as defined in LTE standard, such as be adopted
The DMRS cyclic shift value of each UE is calculated with following formula:It is specifically no longer superfluous
It states.
This step S205 makes the time domain channel response of multiple UE be distributed the midpoint as far as possible far from time domain axis, so that multiple UE
Domain channel response is minimized in true carrier wave and virtual carrier intersection phase change.
This step analyzes the cyclic shift parameter of multiple UE, introduces additive phase rotation in frequency domain before IDFT, so that
The midpoint on entire time domain axis is avoided in the CIR arrangement of multi-user, can be missed to avoid the channel estimation with special displacement UE
Difference is bigger than normal.
S206: the time domain leakage in order to reduce channel response after subsequent DFT, the frequency domain channel obtained in the previous step
Respond HrUpper adding window Wh(Hamming window), concrete operations are as follows: Hw=Hr.*Wh。
This step introduces window adding in frequency domain operation before IDFT is operated, and is not necessarily to anti-windowing operation, window adding in frequency domain after DFT operation
The Gibbs' effect at edge in result in time domain can be reduced, therefore is reduced since the introducing bring time domain energy of virtual carrier is let out
Dew.Because window adding in frequency domain crossover position concentrates on virtual carrier, without increasing reversed window operation after subsequent DFT operation.
S207: the time domain channel shock response of more UE is calculated using IDFT:
S208: following operation is completed to each UE in time domain channel shock response: calculating multi-path location, adding window noise suppression.
Later, the joint time domain channel estimation h of more UE is obtained, this operation once completes the channel estimation of all UE in time domain, specifically,
IDFT operates to obtain the time domain impulse response of all UE, then calculates multi-path location to each UE and adding window noise suppression operates,
Do not have to repeatedly progress IDFT thus to operate.
S209: being directed to each UE, completes inverse cyclic shift (inverse process of step S205) in time domain:Wherein j=0...U-1, this operation in simultaneously remove additional cycles displacement.
S210: being directed to each UE, completes DFT operation after removing noise figure outside window, the frequency domain channel for obtaining current UE is estimated
MeterWherein j=0...U-1, concrete operations areWhereinIt indicates h cyclic shift
It completes, is avoided to each UE as it can be seen that the embodiment of the present application is unified after an IDFT is operated to the CIR processing of more UE
IDFT operation is carried out, computation complexity is reduced.
Device provided by the embodiments of the present application is described below.
Correspond to the above method ground, referring to Fig. 3, a kind of channel estimating apparatus provided by the embodiments of the present application, comprising:
First unit 11, for the frequency domain received data of symbol where determining the demodulated reference signal in each time slot, and
The basic sequence in sequence is sent using the demodulated reference signal of multi-user installation, initial channel estimation is carried out, obtains the letter of frequency domain
Road estimated value；
Second unit 12, for determining virtual carrier using virtual carrier generator matrix, and in the channel estimation value of the frequency domain
The increased virtual carrier in two sides；The wherein virtual carrier generator matrix is the channel estimation value according to virtual carrier and the frequency domain
Positional relationship generate；
Third unit 13, for determining that multi-user sets using the channel estimation value for having the frequency domain of virtual carrier is increased in two sides
Standby joint time domain channel estimated value.
Optionally, the second unit generates the virtual carrier generator matrix W using following formula:
Wherein:
V=0...K-1,12N+K...12N+2K-1
M=0...12N-1
N indicates the number for the resource block that the ascending resource of multiple user equipment multiplexings is included, and K is indicated in the frequency domain
The increased imaginary load wave number in every side of channel estimation value, p indicate that the frequency domain signal-to-noise ratio of virtual carrier application, I indicate unit diagonal matrix.
Optionally, the second unit uses following formula, determines virtual carrier using the virtual carrier generator matrix
Wherein
Wherein, describedIndicate the channel estimation value of the frequency domain,WithRespectively indicate the letter needed in the frequency domain
The increased virtual carrier in the two sides of road estimated value.
Optionally, the third unit, is specifically used for:
Calculate the cyclic shift value of the demodulated reference signal of each user equipmentWherein j=0...U-1, U indicate user
Number of devices；
Calculate so thatThe value range of maximum α, α are as follows: 0 to 12N-
1；
Calculated α completes the phase rotation of H, the channel estimation value H of the frequency domain after obtaining phase rotationr,Wherein H indicates that two sides increase the channel estimation value for having the frequency domain of virtual carrier,
Utilize HrDetermine the joint time domain channel estimated value of multi-user installation.
Optionally, the third unit utilizes HrThe joint time domain channel estimated value for determining multi-user installation, specifically includes:
To the HrCarry out adding window Wh, obtain Hw, Hw=Hr.*Wh；
Calculate HwCorresponding time domain channel shock response
Multi-path location, adding window noise suppression are calculated separately to each user equipment in time domain channel shock response, obtained multiple
The joint time domain channel estimated value h of user equipment；
Equipment for each user calculates inverse cyclic shift valueWherein j=0...U-1,
And after removing noise figure outside window, using the inverse cyclic shift value, the additional of the time domain channel estimated value of the user equipment is removed
Cyclic shift value, and frequency domain is converted to, finally obtain the channel estimation in frequency domain value of the user equipment
Optionally, third unit equipment for each user determines inverse cyclic shift value using following formula
Wherein j=0...U-1.
Other than apparatus above, the embodiment of the present application provides a kind of calculating equipment, which is specifically as follows table
Face computer, portable computer, smart phone, tablet computer, personal digital assistant (Personal Digital
Assistant, PDA) etc..As shown in figure 4, the calculating equipment may include central processing unit (Center Processing
Unit, CPU), memory, input-output apparatus etc., input equipment may include keyboard, mouse, touch screen etc., output equipment
It may include display equipment, such as liquid crystal display (Liquid Crystal Display, LCD), cathode-ray tube (Cathode
Ray Tube, CRT) etc..
Memory may include read-only memory (ROM) and random access memory (RAM), and provide storage to processor
The program instruction and data stored in device.In the embodiment of the present application, memory can be used for storing the embodiment of the present application offer
Channel estimation methods program.
Processor is by the program instruction for calling memory to store, and processor according to the program instruction of acquisition for executing:
The frequency domain received data of symbol where determining the demodulated reference signal in each time slot, and use the demodulation reference of multi-user installation
Signal sends the basic sequence in sequence, carries out initial channel estimation, obtains the channel estimation value of frequency domain；Square is generated using virtual carrier
Battle array determines virtual carrier, and the increased virtual carrier in two sides of the channel estimation value in the frequency domain；Wherein the virtual carrier generates square
Battle array, is that the positional relationship of the channel estimation value according to virtual carrier Yu the frequency domain generates；There is virtual carrier using increasing in two sides
Frequency domain channel estimation value, determine the joint time domain channel estimated value of multi-user installation.
The embodiment of the present application provides a kind of computer storage medium, for being stored as calculating used in above-mentioned calculating equipment
Machine program instruction, it includes the programs for executing above-mentioned channel estimation methods.
The computer storage medium can be any usable medium or data storage device that computer can access, packet
Include but be not limited to magnetic storage (such as floppy disk, hard disk, tape, magneto-optic disk (MO) etc.), optical memory (such as CD, DVD,
BD, HVD etc.) and semiconductor memory (such as it is ROM, EPROM, EEPROM, nonvolatile memory (NAND FLASH), solid
State hard disk (SSD)) etc..
In conclusion technical solution provided by the embodiments of the present application, improves the generating mode of virtual carrier, i.e. virtual carrier is raw
At matrix, it is that the positional relationship of the channel estimation value according to virtual carrier Yu the frequency domain generates, utilizes virtual carrier generator matrix
It determines virtual carrier, and the increased virtual carrier in two sides of the channel estimation value in the frequency domain, has virtual carrier using increasing in two sides
Frequency domain channel estimation value, determine the joint time domain channel estimated value of multi-user installation, therefore by improving virtual carrier and former
There is the phase continuity of subcarrier, to avoid the generation of time domain mirror image diameter under multi-user installation scene, reduces channel
Evaluated error.Further, it by analyzing the cyclic shift value of multiple UE in the embodiment of the present application, is introduced before IDFT in frequency domain
Additive phase rotation, so that the midpoint on time domain axis is avoided in the CIR arrangement of multi-user installation, this operation can be to avoid CIR
It is arranged in the channel estimation results bigger error of the UE near time domain midpoint.Further, it in the embodiment of the present application, is operated in IDFT
Preceding introducing window adding in frequency domain operation, can reduce the introducing bring time domain energy leakage problem due to virtual carrier, it is ensured that letter
The raising of the accuracy of road estimation.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more,
The shape for the computer program product implemented in usable storage medium (including but not limited to magnetic disk storage and optical memory etc.)
Formula.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Obviously, those skilled in the art can carry out various modification and variations without departing from the essence of the application to the application
Mind and range.In this way, if these modifications and variations of the application belong to the range of the claim of this application and its equivalent technologies
Within, then the application is also intended to include these modifications and variations.
Claims (14)
1. a kind of channel estimation methods, which is characterized in that this method comprises:
The frequency domain received data of symbol where determining the demodulated reference signal in each time slot, and use the demodulation of multi-user installation
Reference signal sends the basic sequence in sequence, carries out initial channel estimation, obtains the channel estimation value of frequency domain；
Virtual carrier, and the increased virtual carrier in two sides of the channel estimation value in the frequency domain are determined using virtual carrier generator matrix；
The wherein virtual carrier generator matrix is that the positional relationship of the channel estimation value according to virtual carrier Yu the frequency domain generates；
Using the channel estimation value for having the frequency domain of virtual carrier is increased in two sides, the joint time domain channel estimation of multi-user installation is determined
Value.
2. the method according to claim 1, wherein generating the virtual carrier generator matrix W using following formula:
Wherein:
V=0...K-1,12N+K...12N+2K-1
M=0...12N-1
N indicates the number for the resource block that the ascending resource of multiple user equipment multiplexings is included, and K indicates the channel in the frequency domain
The increased imaginary load wave number in every side of estimated value, p indicate that the frequency domain signal-to-noise ratio of virtual carrier application, I indicate unit diagonal matrix.
3. according to the method described in claim 2, it is characterized in that, utilizing the virtual carrier generator matrix using following formula
Determine virtual carrier
Wherein
Wherein, describedIndicate the channel estimation value of the frequency domain,WithIt respectively indicates and needs to estimate in the channel of the frequency domain
The increased virtual carrier in the two sides of evaluation.
4. according to the method described in claim 3, it is characterized in that, thering is the channel of the frequency domain of virtual carrier to estimate using increasing in two sides
Evaluation determines the joint time domain channel estimated value of multi-user installation, specifically includes:
Calculate the cyclic shift value of the demodulated reference signal of each user equipmentWherein j=0...U-1, U indicate user equipment
Number；
Calculate so thatThe value range of maximum α, α are as follows: 0 to 12N-1；
Calculated α completes the phase rotation of H, the channel estimation value H of the frequency domain after obtaining phase rotationr, wherein H table
Show that two sides increase the channel estimation value for having the frequency domain of virtual carrier,
Utilize HrDetermine the joint time domain channel estimated value of multi-user installation.
5. according to the method described in claim 4, it is characterized in that, utilizing HrDetermine that the joint time domain channel of multi-user installation is estimated
Evaluation specifically includes:
To the HrAdding window is carried out, H is obtainedw；
Calculate HwCorresponding time domain channel shock response；
Multi-path location, adding window noise suppression are calculated separately to each user equipment in time domain channel shock response, obtain multiple users
The joint time domain channel estimated value of equipment；
Equipment for each user calculates inverse cyclic shift value, and after removing noise figure outside window, utilizes the inverse cyclic shift
Value, removes the additional cyclic shift value of the time domain channel estimated value of the user equipment, and be converted to frequency domain, finally obtains the user
The channel estimation in frequency domain value of equipment.
6. according to the method described in claim 5, it is characterized in that, equipment for each user, is determined inverse using following formula
Cyclic shift value
Wherein j=0...U-1.
7. a kind of channel estimating apparatus characterized by comprising
First unit, for the frequency domain received data of symbol where determining the demodulated reference signal in each time slot, and using more
The demodulated reference signal of user equipment sends the basic sequence in sequence, carries out initial channel estimation, obtains the channel estimation of frequency domain
Value；
Second unit, for determining virtual carrier, and the two sides of the channel estimation value in the frequency domain using virtual carrier generator matrix
Increased virtual carrier；The wherein virtual carrier generator matrix is the channel estimation value that the frequency domain is arranged according to virtual carrier
What the arrangement mode of two sides generated；
Third unit, for determining the connection of multi-user installation using the channel estimation value for having the frequency domain of virtual carrier is increased in two sides
Close time domain channel estimated value.
8. device according to claim 7, which is characterized in that the second unit generates the imaginary load using following formula
Wave generator matrix W:
Wherein:
V=0...K-1,12N+K...12N+2K-1
M=0...12N-1
N indicates the number for the resource block that the ascending resource of multiple user equipment multiplexings is included, and K indicates the channel in the frequency domain
The increased imaginary load wave number in every side of estimated value, p indicate that the frequency domain signal-to-noise ratio of virtual carrier application, I indicate unit diagonal matrix.
9. device according to claim 8, which is characterized in that the second unit uses following formula, utilizes the void
Carrier wave generator matrix determines virtual carrier
Wherein
Wherein, describedIndicate the channel estimation value of the frequency domain,WithIt respectively indicates and needs to estimate in the channel of the frequency domain
The increased virtual carrier in the two sides of evaluation.
10. device according to claim 9, which is characterized in that the third unit is specifically used for:
Calculate the cyclic shift value of the demodulated reference signal of each user equipmentWherein j=0...U-1, U indicate user equipment
Number；
Calculate so thatThe value range of maximum α, α are as follows: 0 to 12N-1；
Calculated α completes the phase rotation of H, the channel estimation value H of the frequency domain after obtaining phase rotationr, wherein H table
Show that two sides increase the channel estimation value for having the frequency domain of virtual carrier,
Utilize HrDetermine the joint time domain channel estimated value of multi-user installation.
11. device according to claim 10, which is characterized in that the third unit utilizes HrDetermine multi-user installation
Joint time domain channel estimated value, specifically includes:
To the HrAdding window is carried out, H is obtainedw；
Calculate HwCorresponding time domain channel shock response；
Multi-path location, adding window noise suppression are calculated separately to each user equipment in time domain channel shock response, obtain multiple users
The joint time domain channel estimated value of equipment；
Equipment for each user calculates inverse cyclic shift value, and after removing noise figure outside window, utilizes the inverse cyclic shift
Value, removes the additional cyclic shift value of the time domain channel estimated value of the user equipment, and be converted to frequency domain, finally obtains the user
The channel estimation in frequency domain value of equipment.
12. device according to claim 11, which is characterized in that third unit equipment for each user uses
Following formula determines inverse cyclic shift value
Wherein j=0...U-1.
13. a kind of calculating equipment characterized by comprising
Memory, for storing program instruction；
Processor executes according to the program of acquisition for calling the program instruction stored in the memory: determining each time slot
In demodulated reference signal where symbol frequency domain received data, and using multi-user installation demodulated reference signal send sequence
In basic sequence, carry out initial channel estimation, obtain the channel estimation value of frequency domain；Imaginary load is determined using virtual carrier generator matrix
Wave, and the increased virtual carrier in two sides of the channel estimation value in the frequency domain；The wherein virtual carrier generator matrix, is according to void
What the positional relationship of the channel estimation value of carrier wave and the frequency domain generated；Utilize the channel for the frequency domain for having virtual carrier in two sides increase
Estimated value determines the joint time domain channel estimated value of multi-user installation.
14. a kind of computer storage medium, which is characterized in that the computer-readable recording medium storage has computer executable
Instruction, the computer executable instructions are for making the computer perform claim require 1 to 6 described in any item methods.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710508459.6A CN109150783B (en) | 2017-06-28 | 2017-06-28 | Channel estimation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710508459.6A CN109150783B (en) | 2017-06-28 | 2017-06-28 | Channel estimation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109150783A true CN109150783A (en) | 2019-01-04 |
| CN109150783B CN109150783B (en) | 2020-12-04 |
Family
ID=64803468
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710508459.6A Active CN109150783B (en) | 2017-06-28 | 2017-06-28 | Channel estimation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109150783B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112995084A (en) * | 2021-02-07 | 2021-06-18 | 比科奇微电子(杭州)有限公司 | Signal processing method and processing device |
| WO2022141588A1 (en) * | 2020-12-31 | 2022-07-07 | 北京小米移动软件有限公司 | Communication method and apparatus |
| CN115842701A (en) * | 2023-02-22 | 2023-03-24 | 南京创芯慧联技术有限公司 | Channel estimation method, device, communication equipment and storage medium |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007124761A1 (en) * | 2006-04-27 | 2007-11-08 | Telecom Italia S.P.A. | Frequency domain channel estimation in a single carrier frequency division multiple access system |
| WO2008035836A1 (en) * | 2006-09-21 | 2008-03-27 | Electronics And Telecommunications Research Institute | Time synchronization method and frequency offset estimation method using the same in ofdm network |
| CN101577692A (en) * | 2009-06-15 | 2009-11-11 | 北京邮电大学 | Channel estimating method of orthogonal frequency division multiplexing system and device thereof |
| CN101778063A (en) * | 2010-03-18 | 2010-07-14 | 展讯通信（上海）有限公司 | Channel estimation method and device thereof |
| CN101883074A (en) * | 2010-06-29 | 2010-11-10 | 北京邮电大学 | Cyclic prefix (CP) and virtual carrier based blind frequency offset estimation method in OFDM (Orthogonal Frequency Division Multiplexing) system |
| CN101997807A (en) * | 2009-08-31 | 2011-03-30 | 大唐移动通信设备有限公司 | Channel estimation method and device |
| CN102006255A (en) * | 2009-09-03 | 2011-04-06 | 卓胜微电子（上海）有限公司 | Estimation method of frequency deviation of integral multiple of subcarriers in orthogonal frequency division multiplexing (OFDM) system |
| CN102611656A (en) * | 2012-03-21 | 2012-07-25 | 武汉邮电科学研究院 | Enhanced channel estimation method and enhanced channel estimation device suitable for uplink of LTE (long term evolution) system |
| CN102694763A (en) * | 2012-05-31 | 2012-09-26 | 重庆邮电大学 | Method for assessing integer frequency offset of TD-LTE system |
| CN103283199A (en) * | 2010-10-25 | 2013-09-04 | 爱立信(中国)通信有限公司 | Method and arrangement in wireless communications system |
| CN104639472A (en) * | 2013-11-12 | 2015-05-20 | 普天信息技术研究院有限公司 | Uplink MU-MIMO (multi-user-multiple input multiple output) channel estimation method |
| US20150341101A1 (en) * | 2014-05-22 | 2015-11-26 | Lg Electronics Inc. | Method and apparatus for channel estimation in wireless communication system |
| CN106161316A (en) * | 2015-04-08 | 2016-11-23 | 北京三星通信技术研究有限公司 | Pilot frequency sequence reference signal sending, receiving method and device thereof |
-
2017
- 2017-06-28 CN CN201710508459.6A patent/CN109150783B/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007124761A1 (en) * | 2006-04-27 | 2007-11-08 | Telecom Italia S.P.A. | Frequency domain channel estimation in a single carrier frequency division multiple access system |
| WO2008035836A1 (en) * | 2006-09-21 | 2008-03-27 | Electronics And Telecommunications Research Institute | Time synchronization method and frequency offset estimation method using the same in ofdm network |
| US20100027723A1 (en) * | 2006-09-21 | 2010-02-04 | Hee-Wook Kim | Time synchronization method and frequency offset estimation method using the same in ofdm network |
| CN101577692A (en) * | 2009-06-15 | 2009-11-11 | 北京邮电大学 | Channel estimating method of orthogonal frequency division multiplexing system and device thereof |
| CN101997807A (en) * | 2009-08-31 | 2011-03-30 | 大唐移动通信设备有限公司 | Channel estimation method and device |
| CN102006255A (en) * | 2009-09-03 | 2011-04-06 | 卓胜微电子（上海）有限公司 | Estimation method of frequency deviation of integral multiple of subcarriers in orthogonal frequency division multiplexing (OFDM) system |
| CN101778063A (en) * | 2010-03-18 | 2010-07-14 | 展讯通信（上海）有限公司 | Channel estimation method and device thereof |
| CN101883074A (en) * | 2010-06-29 | 2010-11-10 | 北京邮电大学 | Cyclic prefix (CP) and virtual carrier based blind frequency offset estimation method in OFDM (Orthogonal Frequency Division Multiplexing) system |
| CN103283199A (en) * | 2010-10-25 | 2013-09-04 | 爱立信(中国)通信有限公司 | Method and arrangement in wireless communications system |
| CN102611656A (en) * | 2012-03-21 | 2012-07-25 | 武汉邮电科学研究院 | Enhanced channel estimation method and enhanced channel estimation device suitable for uplink of LTE (long term evolution) system |
| CN102694763A (en) * | 2012-05-31 | 2012-09-26 | 重庆邮电大学 | Method for assessing integer frequency offset of TD-LTE system |
| CN104639472A (en) * | 2013-11-12 | 2015-05-20 | 普天信息技术研究院有限公司 | Uplink MU-MIMO (multi-user-multiple input multiple output) channel estimation method |
| US20150341101A1 (en) * | 2014-05-22 | 2015-11-26 | Lg Electronics Inc. | Method and apparatus for channel estimation in wireless communication system |
| CN106161316A (en) * | 2015-04-08 | 2016-11-23 | 北京三星通信技术研究有限公司 | Pilot frequency sequence reference signal sending, receiving method and device thereof |
Non-Patent Citations (3)
| Title |
|---|
| 卓维娜: ""MIMo-OFDM系统信道估计和预编码技术研究"", 《中国优秀硕士学位论文全文数据库（信息科技辑）》 * |
| 吕喆等: ""一种LTE上行多用户MIMO信道估计算法"", 《电视技术》 * |
| 王亚林: ""LTE-A系统中信道估计技术研究"", 《中国优秀硕士学位论文全文数据库（信息科技辑）》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022141588A1 (en) * | 2020-12-31 | 2022-07-07 | 北京小米移动软件有限公司 | Communication method and apparatus |
| CN112995084A (en) * | 2021-02-07 | 2021-06-18 | 比科奇微电子(杭州)有限公司 | Signal processing method and processing device |
| CN115842701A (en) * | 2023-02-22 | 2023-03-24 | 南京创芯慧联技术有限公司 | Channel estimation method, device, communication equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109150783B (en) | 2020-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11096069B2 (en) | Apparatus and method for hybrid multiple access wireless communication system | |
| US10148463B2 (en) | Method and device for estimating frequency offset | |
| CN103379078B (en) | A kind of method and apparatus of frequency domain equalization detection | |
| CN109150783A (en) | A kind of channel estimation methods and device | |
| CN103384226A (en) | Method and device for frequency domain equalization detection | |
| CN104348763A (en) | Channel measuring method for large-scale antennas, and user terminal | |
| CN103916229A (en) | Interference suppression method and device | |
| US10581651B2 (en) | Apparatus and method for transmitting/receiving signal in multi-carrier system | |
| CN103167549A (en) | Measuring reference signal processing method and device | |
| CN106170936A (en) | The method and apparatus determining signal to noise ratio in radio communication | |
| JP2014183582A (en) | Noise plus interference space covariance matrix determination device and interference rejection combining device | |
| CN102546483A (en) | Method for broadband common-frequency interference noise estimation and interference suppression and corresponding system | |
| CN115834303B (en) | Adaptive frequency domain channel estimation method, device, communication equipment and storage medium | |
| CN105471778B (en) | A kind of signal detecting method and device | |
| JP5788088B2 (en) | Receiving apparatus and receiving method | |
| KR102375186B1 (en) | Apparatus and method for performing channel decoding operation in communication system | |
| CN114465864A (en) | Phase noise compensation method, phase noise compensation device, terminal, storage medium, and program product | |
| CN102480440B (en) | Channel estimation methods under the co-channel interference environment of a kind of broadband and system | |
| CN114826832A (en) | Channel estimation method, neural network training method, device and equipment | |
| CN112398769A (en) | Modulation mode detection method and device, electronic equipment and readable storage medium | |
| WO2019205111A1 (en) | Method for merging data, apparatus, and device | |
| US20150207588A1 (en) | Receiving device and receiving method, and computer program | |
| CN110881200B (en) | Signal transmission method and system based on non-orthogonal multiple access and electronic equipment | |
| CN117118791B (en) | Frequency offset estimation method, device, computer equipment and storage medium | |
| CN103378883A (en) | Method and device for equalizing frequency domain under condition of uplink MU-MIMO |
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 |