A kind of speed measurement method and equipment based on channel estimating
Technical field
The present invention relates to communication technical field, relate in particular to a kind of speed measurement method and equipment based on channel estimating.
Background technology
Current Computational Method of Velocity Measurement comprises: correlation class algorithm and level crossing rate (lcr) method.
Speed-measuring method 1: correlation class methods.
ACF(auto-correlation function) method: for the channel of Rayleigh fading, derive normalized time domain auto-correlation function expression formula as follows:
(1)
represent 0 rank Bessel function of the first kind,
, be the time interval for calculating correlated channels response; Based on estimated channel response, the auto-correlation of the channel response of statistics on certain hour interval, releases corresponding speed by formula (1) is counter.
For Rice channel, because arrival bearing is not equally distributed, and be subject to the impact of Rice factor K, need to revise formula (1), and can not directly use, do not launch in detail at this.
Cov(Covariance, auto-covariance) method: it is defined as:
(2)
represent channel coefficients envelope square, T
1for the time-domain sampling time interval of channel, the theoretical expression of this definition amount is:
(3)
, be similar to ACF, can utilize the corresponding relation between statistic and speed, the result of estimating based on channel response obtains statistic, and then obtains the estimation of speed.
Speed-measuring method 2: based on the method for level crossing rate (lcr), Doppler spread can cause signal to rise and fall in time domain, every mobile half wavelength distance generally, signal amplitude is deep fade once; The envelope of measuring channel response exceedes within a certain period of time the number of times of a certain thresholding or obtains the number of times of maximum or become the minimizing time from maximum.For example, by level decline number of times Le in the statistical unit time, can estimate speed; Suppose carrier frequency fc, the light velocity is c, can estimate speed v=c/fc*Le.
It should be noted that, in above-mentioned correlation class algorithm and level crossing rate (lcr) method, do not consider the characteristic (as PDP(time domain power delay spectrum) of actual channel etc.), but theoretical expression based on certain, therefore its rate accuracy that can reach is limited.
In order to improve rate accuracy, in prior art, propose a kind ofly to do poor method based on channel response and (can consider the difference of interchannel, to improve rate accuracy), in its scope capable of measuring speed, linear between velocity function and speed, both corresponding relations are pointed out by the concrete statistics under various channels, and its defined speed statistical function is:
(4)
Wherein, k represents subcarrier number, and n represents the n time sampling of channel,
represent the n time sampling of k subcarrier of channel response.Channel without decline and the Rayleigh fading of single footpath under, can prove in theory
linear with Doppler frequency, so linear with speed.In actual applications, channel response is produced by channel estimating conventionally, and because the method has been utilized the difference (i.e. sampling before and after variation sizes) of channel, therefore can be described as based on channel response and change big or small method.
Realizing in process of the present invention, inventor finds at least to exist in prior art following problem:
In level crossing rate (lcr) method, cannot go out Le by accurate statistics.Signal, owing to being subject to noise and channel effect, is observed and can be found a lot of burrs from time domain; Therefore first to carry out denoising, deburring processing to signal, otherwise cannot accurate statistics level decline number of times.In addition, the statistical of level crossing rate (lcr) also has impact to velocity estimation accuracy rate, and the precision of this algorithm is not high, and needs under the low speed the timing statistics of growing.
In correlation class algorithm, inapplicable under Rice channel, need in other cases to utilize Rice's factor correction algorithm, but Rice's factor K is difficult for determining.Due to Bezier Non-monotonic function, for accurate estimating speed, need to ensure
; And under high speed Doppler spread
larger,
value need very littlely just can carry out the estimation of Doppler spread, thereby make correlation class algorithm under high-speed case, be subject to larger restriction.In addition, correlation class algorithm also can affect the precision of algorithm.
For change big or small method based on channel response, under different channel scenes, between velocity function and the simulation curve of speed, there are differences, and velocity function and speed are linear relationship in theory, but under different channels, to there is larger performance difference, and in prior art, ignore this difference, thereby cause the precision of channel estimating to be subject to larger performance loss.
Summary of the invention
The embodiment of the present invention provides a kind of speed measurement method and equipment based on channel estimating, to improve rate accuracy.
In order to achieve the above object, the embodiment of the present invention provides a kind of speed measurement method based on channel estimating, comprising:
Acquisition channel response is estimated, and is chosen the channel response of time domain most powerful path in described channel response estimation;
Obtain velocity information according to the channel response of described time domain most powerful path.
The embodiment of the present invention provides a kind of speed measuring equipment based on channel estimating, comprising:
Choose module, estimate for obtaining channel response, and choose the channel response of time domain most powerful path in described channel response estimation;
Obtain module, for obtaining velocity information according to the channel response of described time domain most powerful path.
Compared with prior art, the embodiment of the present invention at least has the following advantages:
Can effectively solve Computational Method of Velocity Measurement based on the channel variation size velocity function under different channels and add up the problem of dispersing, significantly increase the convergence of speed statistical function under different channels.
Brief description of the drawings
In order to be illustrated more clearly in technical scheme of the present invention, below the accompanying drawing of required use during embodiment is described is briefly described, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, do not paying under the prerequisite of creative work, can also obtain according to these accompanying drawings other accompanying drawing.
Fig. 1 is a kind of speed measurement method schematic flow sheet based on channel estimating that the embodiment of the present invention one provides;
Fig. 2 is a kind of speed measuring equipment structural representation based on channel estimating that the embodiment of the present invention two provides.
Embodiment
Below in conjunction with the accompanying drawing in the present invention, the technical scheme in the present invention is clearly and completely described, obviously, described embodiment is only a part of embodiment of the present invention, instead of whole embodiment.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtaining under creative work prerequisite, belong to the scope of protection of the invention.
Embodiment mono-
The embodiment of the present invention one provides a kind of speed measurement method based on channel estimating, and as shown in Figure 1, the method comprises the following steps:
Step 101, obtains channel response and estimates.Wherein, can obtain channel response based on existing system flow process and estimate, in this process embodiment of the present invention, repeat no more.
Step 102, choose channel response estimate in the channel response of time domain most powerful path.
In the embodiment of the present invention, the process of choosing the channel response of time domain most powerful path in channel response estimation comprises: in the time that channel response is estimated as time domain channel response, choose the channel response of time domain most powerful path in time domain channel response; Or,
In the time that channel response is estimated as domain channel response, domain channel response is converted to time domain channel response, and chooses the channel response of time domain most powerful path in the time domain channel response after conversion.
Further, in choosing the channel response of time domain most powerful path in time domain channel response or choosing the time domain channel response after conversion in the process of the channel response of time domain most powerful path, can choose the strongest time delay footpath of power in the time domain channel response after time domain channel response or conversion, and re-construct the channel response of time domain most powerful path by the strongest time delay footpath of power.
Step 103, obtains velocity information according to the channel response of time domain most powerful path.
In the embodiment of the present invention, the process that obtains velocity information according to the channel response of time domain most powerful path comprises: the speed statistical function that calculates channel according to the channel response of time domain most powerful path; And according to the corresponding relation of speed statistical function and predefined speed statistical function and velocity information, obtain velocity information (being the estimation of speed or the interval of velocity estimation).
Further, calculating according to the channel response of time domain most powerful path in the process of speed statistical function of channel, the channel response of time domain most powerful path is time domain channel response, and the speed statistical function that calculates channel is time domain speed statistical function; Or,
The channel response of time domain most powerful path is converted to domain channel response, and the speed statistical function that calculates channel is frequency domain speed statistical function.
It should be noted that, the technical scheme that the embodiment of the present invention provides can be applicable in LTE system, for LTE system being carried out to up testing the speed or descending testing the speed; Certainly, also can be applied in other system and test the speed, it is applicable to the system of all uses based on channel response variation Computational Method of Velocity Measurement, in the embodiment of the present invention, no longer application scenarios is described.
The technical scheme embodiment of the present invention being provided below in conjunction with LTE system is further detailed, based on PUSCH(Physical Uplink Shared Channel, Physical Uplink Shared Channel) the up mode of testing the speed comprise the following steps:
Step 1, obtain channel response based on PUSCH channel (different scenes this channel can different) and estimate.
In the present embodiment, if this channel response is estimated as time domain channel response, can suppose that the time domain channel response (being initial channel estimation) of two time slot Slot of a certain user on scheduling subframe i is h
i1and h (n)
i2(n), n=1 ..., N.
If this channel response is estimated as domain channel response, the domain channel response (being initial channel estimation) that can suppose two time slot Slot of a certain user on scheduling subframe i is H
i1and H (k)
i2(k), k represents subcarrier number, and value is k=1 ..., K, and K equals the sub-carrier number that PUSCH occupies.
Step 2, choose channel response estimate in the channel response of time domain most powerful path.
Mode one: in the time that channel response is estimated as time domain channel response, choose the time delay footpath that in time domain channel response, power is the strongest, and re-construct the channel response of time domain most powerful path by the strongest time delay footpath of power.
Concrete, choose the time delay footpath that power is the strongest and be:
with
, and determine the strongest time delay footpath n of power obtaining on two Slot
0=n
1or n
2; Or the strongest time delay footpath of power of choosing two Slot upper signal channels estimations is:
.
The channel response that re-constructs time domain most powerful path by the strongest time delay footpath of power is:
。
In the embodiment of the present invention, choosing the time delay footpath (n that power is the strongest
1and n
2) time, common two Slot upper signal channels are estimated the strongest identical (n in time delay footpath of power obtaining
1=n
2=n
0), if different, can some slot being as the criterion (is n
0=n
2or n
0=n
1); Or directly obtain n
0(
), to ensure an only having most powerful path.In the present embodiment, be by time delay the footpath the strongest power of choosing
footpath.
Mode two: in the time that channel response is estimated as domain channel response, domain channel response is converted to time domain channel response, and choose the strongest time delay footpath of power in the time domain channel response after conversion, re-construct the channel response of time domain most powerful path by the strongest time delay footpath of power.
Concrete, domain channel response is converted to time domain channel response and comprises: by H
i1(k) carry out IFFT(Fourier inversion) obtain h to time domain
i1(n), and by H
i2(k) carry out IFFT and obtain h to time domain
i2(n), n=1 ..., N, and N<=K.
Choosing the time delay footpath that power is the strongest is:
with
, and determine the strongest time delay footpath n of power obtaining on two Slot
0=n
1or n
2; Or the strongest time delay footpath of power of choosing two Slot upper signal channels estimations is:
.
The channel response that re-constructs time domain most powerful path by the strongest time delay footpath of power is:
。
In the embodiment of the present invention, choosing the time delay footpath (n that power is the strongest
1and n
2) time, common two Slot upper signal channels are estimated the strongest identical (n in time delay footpath of power obtaining
1=n
2=n
0), if different, can be as the criterion by some slot; Or directly obtain n
0(
), to ensure an only having most powerful path.In the present embodiment, be by time delay the footpath the strongest power of choosing
footpath.
Step 3, calculate the speed statistical function of channel according to the channel response of time domain most powerful path (re-constructing the channel response of time domain most powerful path).
Mode one: because the channel response of time domain most powerful path is time domain channel response, the speed statistical function that can directly calculate channel is time domain speed statistical function.
Time domain speed statistical function is:
。
Mode two: because the channel response of time domain most powerful path is time domain channel response, the channel response of time domain most powerful path can be converted to domain channel response, and the speed statistical function that calculates channel is frequency domain speed statistical function.Domain channel response after conversion is: H
i1' (k), H
i2' (k), k=1 ..., K.
Frequency domain speed statistical function is:
。
Step 4, according to the corresponding relation of speed statistical function and predefined speed statistical function and velocity information, obtain velocity information (being the estimation of speed or the interval of velocity estimation).
Concrete, in the embodiment of the present invention, can be to I(I>=1) velocity function in individual subframe averages, and obtains the estimated value of velocity function.
Afterwards, according to the corresponding relation of velocity function and speed, obtain the estimation of this velocity estimation or velocity interval.
Afterwards, can average velocity estimation result repeatedly, obtain the estimation of present speed.
Said process is to describe in the up mode of testing the speed of PUSCH, and below the descending mode of testing the speed of LTE is described.
Compared with up, LTE is descending can use CRS(Common Reference Signal, public reference signal) channel estimation results carry out velocity estimation, it can use the channel estimation results in CRS channel estimation results or the part bandwidth on full bandwidth, and the time interval of the channel of work difference can be in a subframe or adjacent sub-frame or several subframe.
In sum, existingly change in big or small method different concrete interchannels based on channel response, there is larger difference in its defined velocity function and rate curve, if do not processed, has directly reduced the precision of its velocity estimation, and process, can there be two kinds of thinkings:
The one, extract the characteristic of channel that determines curve, for different curves is formulated different velocity functions and the corresponding relation of speed, and need in actual applications to extract the corresponding characteristic of actual channel, select its corresponding corresponding relation; The 2nd, adopt certain processing convergence of function under different channels that gather way.
In the embodiment of the present invention, provide a kind of constringent method under increase channel, by channel response is processed, extracted most powerful path, and carry out the statistics of velocity function based on channel response corresponding to this most powerful path, can significantly strengthen the convergence of channel.And do not increasing the precision that under corresponding relation curve, raising is tested the speed, it has larger terseness.
And can effectively solve Computational Method of Velocity Measurement based on the channel variation size velocity function under different channels and add up the problem of dispersing, significantly increase the convergence of speed statistical function under different channels.
Embodiment bis-
Based on the inventive concept same with said method, a kind of speed measuring equipment based on channel estimating is also provided in the embodiment of the present invention, as shown in Figure 2, this equipment comprises:
Choose module 11, estimate for obtaining channel response, and choose the channel response of time domain most powerful path in described channel response estimation;
Obtain module 12, for obtaining velocity information according to the channel response of described time domain most powerful path.
The described module 11 of choosing, specifically in the time that described channel response is estimated as time domain channel response, chooses the channel response of time domain most powerful path in described time domain channel response; Or,
In the time that described channel response is estimated as domain channel response, described domain channel response is converted to time domain channel response, and chooses the channel response of time domain most powerful path in the time domain channel response after conversion.
The described module 11 of choosing, is further used for choosing the strongest time delay footpath of power in the time domain channel response after time domain channel response or conversion, and re-constructs the channel response of time domain most powerful path by the strongest time delay footpath of described power.
The time domain channel response of supposing two time slot Slot of a certain user on scheduling subframe i is h
i1and h (n)
i2(n), described in choose module 11, be further used for choosing the time delay footpath that power is the strongest and be:
with
, and determine the strongest time delay footpath n of power obtaining on two Slot
0=n
1or n
2; Or the strongest time delay footpath of power of choosing two Slot upper signal channels estimations is:
;
The channel response that re-constructs time domain most powerful path by the strongest time delay footpath of described power is:
。
The domain channel response of supposing two time slot Slot of a certain user on scheduling subframe i is H
i1and H (k)
i2(k), k represents subcarrier number, described in choose module 11, be further used for H
i1(k) carry out Fourier inversion IFFT and obtain h to time domain
i1(n), and by H
i2(k) carry out IFFT and obtain h to time domain
i2(n), n=1 ..., N;
Choosing the time delay footpath that power is the strongest is:
with
, and determine the strongest time delay footpath n of power obtaining on two Slot
0=n
1or n
2; Or the strongest time delay footpath of power of choosing two Slot upper signal channels estimations is:
;
The channel response that re-constructs time domain most powerful path by the strongest time delay footpath of described power is:
。
Described acquisition module 12, specifically for calculating the speed statistical function of channel according to the channel response of described time domain most powerful path; According to the corresponding relation of described speed statistical function and predefined speed statistical function and velocity information, obtain velocity information.
Described acquisition module 12, is further used in the time that the channel response of described time domain most powerful path is time domain channel response, and the speed statistical function that calculates channel is time domain speed statistical function; Or,
The channel response of described time domain most powerful path is converted to domain channel response, and the speed statistical function that calculates channel is frequency domain speed statistical function.
The channel response of described time domain most powerful path comprises:
,
The domain channel response after conversion is: H
i1' (k), H
i2' (k);
Described time domain speed statistical function is:
,
Described frequency domain speed statistical function is:
。
Wherein, the modules of apparatus of the present invention can be integrated in one, and also can separate deployment.Above-mentioned module can be merged into a module, also can further split into multiple submodules.
Through the above description of the embodiments, those skilled in the art can be well understood to the mode that the present invention can add essential general hardware platform by software and realize, and can certainly pass through hardware, but in a lot of situation, the former is better execution mode.Based on such understanding, the part that technical scheme of the present invention contributes to prior art in essence in other words can embody with the form of software product, this computer software product is stored in a storage medium, comprise that some instructions (can be personal computers in order to make a computer equipment, server, or the network equipment etc.) carry out the method described in each embodiment of the present invention.
It will be appreciated by those skilled in the art that accompanying drawing is the schematic diagram of a preferred embodiment, the module in accompanying drawing or flow process might not be that enforcement the present invention is necessary.
It will be appreciated by those skilled in the art that the module in the device in embodiment can be distributed in the device of embodiment according to embodiment description, also can carry out respective change and be arranged in the one or more devices that are different from the present embodiment.The module of above-described embodiment can be merged into a module, also can further split into multiple submodules.
The invention described above embodiment sequence number, just to describing, does not represent the quality of embodiment.
Disclosed is above only several specific embodiment of the present invention, and still, the present invention is not limited thereto, and the changes that any person skilled in the art can think of all should fall into protection scope of the present invention.