Background
In a mobile communication system, the motion of a mobile station causes doppler shift and channel quality change, and high requirements are put on technologies such as channel estimation and link adaptation in the mobile communication system, so the motion velocity of the mobile station needs to be measured in the mobile communication system to switch operating modes according to different motion velocities.
The signal received by the mobile station from the base station will change in phase and amplitude due to the motion of the mobile station, and therefore, the change in the velocity of the mobile station can be measured using the change. For the measurement of the moving speed of the mobile station, there are many companies and research institutions that have conducted extensive research, and there are two main types of measurement methods: one is to determine the moving speed of the mobile station by using the correlation between the frequency of the change of the best cell and the speed of the mobile station; the other is to determine the doppler shift and thus the moving velocity of the mobile station by using the correlation between the time domain or frequency domain variation characteristics of the received signal and the doppler shift. Two types of measurement methods are described in detail below:
in the first category, cell signal quality measurement methods are preferred.
If the moving speed of the mobile station is fast, the switching of the adjacent cells is inevitably carried out, so the moving speed corresponding to the mobile station can be obtained by recording the switching condition of the mobile station. However, in the mobile communication system adopting the hierarchical cell structure, the moving speed of the mobile station is determined based on the number of changes of the microcell base station having the best signal quality measured by the mobile station per unit time, rather than the actual number of handovers. Therefore, in the specific measurement step, a threshold of the number of changes corresponding to each speed needs to be preset, and the movement speed of the mobile station can be determined by comparing the number of changes of the femtocell base station with the best signal quality with the threshold.
The method needs to utilize the indication information of the high layer, simultaneously monitors the signal quality of a plurality of cells, records the change condition of the adjacent micro cell with the best signal quality in real time, has a complex realization process, and has a long period for measuring the movement speed because the cell is switched for a long time.
And the second type, a measurement method according to the received signal variation.
The movement of the mobile station causes the change of the channel information, and further causes the change of the amplitude and phase of the received signal and the channel information, and the change of the movement speed of the mobile station can be obtained by using the change characteristic. The following two methods are commonly used:
a. the first zero-crossing measuring method comprises the following steps:
the first zero-crossing measurement is to estimate the moving speed of the mobile station by using the relation between the channel-related statistical characteristics and the Doppler frequency shift caused by the movement of the mobile station.
According to the classical statistical theory, the time domain correlation property of the channel information satisfies the following formula:
wherein τ represents a time interval, f
dDenotes the maximum Doppler shift, J
0(. cndot.) represents a first class of zeroth order Bessel functions,
representing the noise power spectral density.
Bessel function J0(x) The value of x at the first zero crossing is a fixed value: x is 2.405. If appropriate τ is selected to make J0(x) Equal to zero, the maximum doppler shift can be found:
<math><mrow><msub><mi>f</mi><mi>d</mi></msub><mo>=</mo><mfrac><mn>2.405</mn><mrow><mn>2</mn><mi>πτ</mi></mrow></mfrac></mrow></math>
thus, by calculating
(τ) the maximum doppler shift can be calculated and the velocity of the mobile station can be measured.
The first zero-crossing measurement method needs to firstly obtain correct time-domain impulse response information of the mobile channel, so that the requirement on mobile channel estimation is high, the estimation precision of the method is greatly influenced by the estimation result of the mobile channel, in addition, the selection of the time interval tau is sensitive for finding the first zero point, when the motion speed is different, the tau value during zero crossing is changed, and the complexity of the system is increased.
b. Doppler power spectrum measurement method:
for the doppler power spectral characteristics of a mobile channel, it is generally assumed theoretically that it satisfies the classical doppler power spectral profile as shown in fig. 1, with the largest doppler shifts-Fd and Fd on both sides of the vertical axis. Therefore, the maximum Doppler frequency shift of the mobile channel can be estimated by using the Doppler power spectrum of the mobile channel, and the motion speed of the mobile station can be further obtained.
In the method, firstly, a mobile channel parameter of a time domain is obtained, then Fourier transformation is carried out on the mobile channel parameter to obtain a Doppler power spectrum of the mobile channel, then a catastrophe point on the power spectrum is searched to obtain the maximum Doppler frequency shift, and further the motion speed of the mobile station is obtained.
The premise of the Doppler power spectrum measurement method is as follows: supposing that the doppler power spectrum of the mobile signal meets the distribution diagram of the classical doppler power spectrum, and in practical situations, the doppler power spectrum of the mobile channel has a large difference from the distribution diagram of the classical doppler power spectrum, and the discontinuity point in the power spectrum is not equal to the actual maximum doppler frequency shift, so that the measurement error is large; meanwhile, the pilot spacing must meet the requirement of the minimum sampling rate of the doppler spectrum.
As can be seen from the above description of the method for measuring the moving speed of the mobile station in the prior art, the conventional measuring method is complex in implementation and is not easy to implement.
Disclosure of Invention
In view of the above, the main objective of the present invention is to provide a method and an apparatus for measuring the moving speed of a mobile station, which can simply and effectively measure the moving speed of the mobile station.
In order to achieve the purpose, the technical scheme of the invention is realized as follows:
the invention provides a method for measuring the moving speed of a mobile station, which determines the corresponding relation between the relevant included angle value of a channel vector and the moving speed of the mobile station in a set time interval, and the method also comprises the following steps:
A. obtaining channel information estimation values of pilot subcarriers in a downlink subframe, and respectively calculating relevant included angle values of corresponding channel vectors of the selected pilot subcarriers in corresponding time intervals;
B. and determining the movement speed of the mobile station according to the calculation result and the corresponding relation.
Wherein, when only two rows of pilot frequency sub-carriers are selected in each downlink sub-frame, the determining the moving speed of the mobile station in step B is: directly determining the movement speed of the mobile station according to the calculation result and the corresponding relation;
when more than two pilot frequency sub-carriers are selected in each downlink sub-frame, the motion speed of the mobile station determined in the step B is: the moving speed of the mobile station corresponding to each calculation result is determined, and then the final moving speed of the mobile station is calculated.
The determining the corresponding relationship is specifically as follows: and determining the corresponding relation according to the actual transmission environment or simulation estimation.
And the calculation result is the relevant included angle value in the step A.
Further comprising the following steps between the step A and the step B:
a1, continuously counting the relevant included angle values of the channel vectors corresponding to the pilot subcarriers with the same time of each pilot subcarrier in the step A and the corresponding time intervals in a certain number of downlink subframes, and respectively and correspondingly calculating the mean value of the relevant included angles in each time interval according to the obtained relevant included angle values;
correspondingly, the calculation result is the mean value of the correlation included angles.
The method further comprises the following steps between the step A1 and the step B:
a2, respectively eliminating the relevant included angle values which deviate from the relevant included angle mean value corresponding to the relevant included angle value by a certain ratio, and respectively recalculating each new relevant included angle mean value;
correspondingly, the calculation result is the new mean value of the correlation included angle.
The corresponding channel vectors are: and two rows of channel vectors corresponding to the pilot frequency sub-carriers corresponding to each time interval.
The invention also provides a device for measuring the moving speed of the mobile station, which comprises: a corresponding relation obtaining module, a related included angle mean value calculating module and a motion speed determining module, wherein,
a corresponding relation obtaining module for determining the corresponding relation between the relevant included angle value of the channel vector and the moving speed of the mobile station in each set time interval and sending the determined corresponding relation to the moving speed determining module;
the device comprises a correlation included angle mean value calculation module, a motion speed determination module and a channel vector calculation module, wherein the correlation included angle mean value calculation module is used for obtaining a channel information estimation value of a pilot frequency subcarrier in each downlink subframe, calculating a correlation included angle value of a channel vector according to the estimation value and sending a calculation result to the motion speed determination module;
and the movement speed determining module is used for determining the movement speed of the mobile station according to the corresponding relation and the calculation result respectively sent by the corresponding relation obtaining module and the related included angle mean value calculating module.
Wherein the calculation result is: the relevant clip angle value.
The correlation angle mean calculation module is further configured to: continuously counting relevant included angle values corresponding to a certain number of downlink subframes, and calculating a relevant included angle mean value according to the obtained relevant included angle values;
correspondingly, the calculation result is: the mean value of the correlation angle.
The correlation angle mean calculation module is further configured to: respectively eliminating related included angle values deviating from the corresponding related included angle mean value by a certain ratio in the related included angle values, and respectively recalculating each new related included angle mean value;
correspondingly, the calculation result is: new mean correlation angle.
The method and the device for measuring the moving speed of the mobile station provided by the invention utilize the pilot frequency subcarrier in the downlink subframe to calculate the channel information estimation value and then calculate the relevant included angle value of the channel vector, thereby obtaining the moving speed of the mobile station, and the method and the device are simple to realize and do not need complex software and hardware equipment. The invention continuously counts the relevant included angle values corresponding to a plurality of downlink sub-frames, or calculates the relevant included angle values corresponding to a plurality of rows of pilot frequency sub-carriers in each downlink sub-frame, and calculates the mean value of the relevant included angles, thereby increasing the objectivity and reliability of measuring the movement speed of the mobile station. After the mean value of the relevant included angles is calculated, the relevant included angle values deviating from the mean value of the relevant included angles greatly are removed, and then the mean value of the relevant included angles is recalculated, so that the reliability of the measurement of the moving speed of the mobile station can be further enhanced.
Detailed Description
The basic idea of the invention is: and calculating the related included angle values of two channel vectors at a certain time interval on the pilot frequency sub-carrier in the downlink sub-frame, and determining the movement speed of the mobile station according to the related included angle values.
The basic principle of the invention is as follows:
in an Orthogonal Frequency Division Multiplexing (OFDM) system, the correlation function of channel information values on a certain pilot subcarrier can be expressed as:
wherein f isdFor maximum Doppler frequency offset, J0(. cndot.) is a first class of zeroth order Bessel function. For actual channel information estimates, there is often a noise signal, so the above equation can be written as:
wherein,
for noise power spectral density, H is the channel information value, k is the kth time, τ is the time interval, f
dMaximum Doppler shift, J
0(. cndot.) is a first class of zeroth order Bessel function,
is the noise power spectral density.
To overcome the effect of noise, the phase of the channel vector at both times k and k + τ may be usedThe average value of the pinch angle values is used to smooth the compensation. If the channel information values at adjacent pilot symbols are shorter in time interval, the included angle is fdAnd the pilot frequency time interval tau, and in a linear relation, the Doppler frequency shift of the mobile station can be obtained by solving the relevant included angle value, and then the motion speed of the mobile station is solved. The included angle is expressed as follows:
<math><mrow><mi>θ</mi><mo>=</mo><mi>arg</mi><mrow><mo>(</mo><mi>corr</mi><mrow><mo>(</mo><msub><mover><mi>H</mi><mo>~</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>H</mi><mo>~</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>τ</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>
if limiting fdThe value of x tau is such that the value of the relevant included angle is 0, pi]In the range of (1), wherein the above condition is easily satisfied in an actual mobile communication system, the magnitude of the doppler shift can be determined by using the cumulative mean of the correlation included angle values, and the moving speed of the mobile station can be obtained.
The following detailed description is made by way of specific embodiments in conjunction with the accompanying drawings.
Fig. 2 is a flowchart illustrating a method for measuring a moving speed of a mobile station using two rows of pilot signals according to the present invention, as shown in fig. 2, the method includes:
step 201: according to the actual transmission environment or simulation estimation, the corresponding relation between the relevant included angle value of the channel vector and the moving speed of the mobile station under a certain fixed time interval tau is obtained.
The time interval τ may be selected according to the requirement of measurement accuracy, but is not limited, when a certain time interval τ is selected, the time interval τ selected in the following step cannot be changed.
For the selection of the time interval τ, two rows of pilot subcarriers may be selected first, and then the time interval τ is obtained, where k and k + τ in step 203 are two time instants corresponding to the pilot subcarriers selected in this step.
For example, in a 3GPP Long Term Evolution (LTE) Time Division Duplex (TDD) Type2 system, it is assumed that a channel model is a global system for mobile communications (GSM) TU channel, each subframe is 5ms, each subframe includes 7 slots, only one slot in one subframe is used for a downlink channel, two columns of OFDM symbols in a downlink common pilot frequency domain have pilot symbols inserted, a pilot symbol time domain is separated by six OFDM symbols, and a signal-to-noise ratio is 10 dB. Then, a curve diagram of the mobile station moving speed and the relevant included angle value shown in fig. 3 is obtained through simulation statistics, and then the corresponding relation of the mobile station moving speed and the relevant included angle value shown in table 1 is obtained according to fig. 3:
angle (degree) |
[0-26] |
[26-42] |
[42-80] |
Speed (kilometer/hour) |
Low speed (0-125) |
Medium speed (75-230) |
High speed (170- |
TABLE 1
The above is only an example of the LTE TDD Type2 system, and the method and apparatus of the present invention are not limited to the LTE TDD Type2 system, and can be applied to TDD and Frequency Division Duplex (FDD) systems.
Specifically, how to obtain the corresponding relationship between the moving speed of the mobile station and the related included angle value according to the actual transmission environment or simulation estimation belongs to the known technology, and is not described herein again.
Step 202: and obtaining the channel information estimation value of the pilot frequency subcarrier in the downlink subframe.
Here, the obtained signal values of the pilot subcarriers may be expressed as:
Y=Hx+n
where H is the channel information value, x is the transmitted pilot symbol value, and n is noise.
In the downlink sub-frame, the signal value Y of the pilot frequency sub-carrier is divided by the transmitted pilot frequency symbol x, so as to obtain the channel information estimation value of the pilot frequency sub-carrier Wherein the influence of the noise n is neglected.
Step 203: and according to the obtained channel information estimation value, taking the time as k and the time interval tau, and calculating the related included angle value of the channel vector on the pilot frequency subcarrier at the corresponding k time and the k + tau time.
Wherein the channel vector corresponds to the channel information estimate. The correlation angle values refer to: the value of the correlation clip between the channel vector at time k and the channel vector at time k + τ.
The associated clip angle value of the two channel vectors is calculated according to the following formula:
<math><mrow><mi>θ</mi><mo>=</mo><mi>arg</mi><mrow><mo>(</mo><mi>corr</mi><mrow><mo>(</mo><msub><mover><mi>H</mi><mo>~</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><msub><mover><mi>H</mi><mo>~</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>τ</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>
where τ is the time interval and k denotes the kth time.
There is no limitation on the selection of the time k, and the selection can be performed at will, but once a certain time k is selected, the pilot subcarriers at the same time k must be selected when the selection of the pilot subcarriers in the next frame is performed in step 204.
Step 204: and continuously counting the corresponding correlation included angle values of the plurality of downlink subframes, and calculating a mean value of the correlation included angles according to the obtained plurality of correlation included angle values.
There are several methods for processing the statistical correlation pinch value, for example: the related included angle value of the current statistics can be partially covered with the related included angle value of the previous statistics, specifically, the related included angle mean value is calculated by setting the related included angle value corresponding to 100 frames of statistics each time, when the related included angle value corresponding to 100 frames of statistics is counted, the related included angle mean value is calculated, when 150 frames of statistics are counted, if the related included angle mean value needs to be calculated, the related included angle value corresponding to 100 frames to 150 frames is covered with the related included angle value corresponding to 0 to 100 frames of previous statistics, and therefore the related included angle mean value is calculated by still using the related included angle value corresponding to 100 frames; or, only the corresponding related included angle values of 100 frames may be stored, and after counting a new related included angle value, the stored corresponding related included angle values of 100 frames are dynamically updated; or, when the storage space for storing the correlation included angle value is large, the correlation included angle values of more frames can be stored, so as to meet the requirement of performing correlation included angle mean calculation by counting the correlation included angle values of more frames, such as: the calculation of the mean value of the correlation included angle is carried out by counting the correlation included angle values of 100 frames each time, and the stored correlation included angle value can be far larger than 100 frames, so that when the corresponding correlation included angle values of 200 frames or more need to be counted for measuring the movement speed of the mobile station at a certain time, and the calculation of the mean value of the correlation included angle is carried out, enough correlation included angle values can still be provided.
Continuously counting relevant included angle values corresponding to a plurality of downlink subframes, namely, circularly executing steps 202-203, and obtaining the relevant included angle value corresponding to one downlink subframe every time when executing; and then, calculating the mean value of the related included angles according to the counted related included angle values corresponding to the plurality of downlink subframes.
The number of the correlation included angle values that need to be calculated as the mean value of the correlation included angles depends on the actually required measurement accuracy and the signal-to-noise ratio of the current received signal, and there is no specific upper limit to the number of the correlation included angle values that need to be counted.
In addition, the time intervals τ used in the calculation of the correlation pinch value are all the time intervals τ selected in step 201.
Step 205: and eliminating the correlation included angle value which deviates from the correlation included angle mean value and is larger, and then recalculating a new correlation included angle mean value.
The determination of the data to be rejected depends on the actually required measurement accuracy, and for example, it may be set as follows: and eliminating the related included angle value deviating from the related included angle mean value by more than 20 percent of the related included angle mean value.
Step 206: and determining the motion speed of the mobile station according to the obtained new correlation included angle mean value by referring to the obtained corresponding relation shown in the table 1.
The determining the moving speed of the mobile station specifically comprises: and determining which angle interval of the corresponding relation the newly calculated average value of the related angles is in, and then the moving speed of the mobile station is the speed interval corresponding to the angle interval.
In the calculation process of the moving speed of the mobile station shown in fig. 2, steps 204 and 205 are optional steps.
When step 204 is removed, step 205 must be removed accordingly, and at this point, in step 206, the velocity of the mobile station is determined based on the correspondence and the associated included angle value calculated in step 202.
Step 204 may be retained when step 205 is eliminated, and in this case, the motion speed of the mobile station is determined in step 206 according to the correspondence and the mean correlation angle calculated in step 204.
In the process of calculating the moving speed of the mobile station in fig. 2, the correlation included angle value, the correlation included angle mean value, and the recalculated correlation included angle mean value may be collectively referred to as a calculation result.
In addition, the above-mentioned method flow for measuring the moving speed of the mobile station once is only given, and when the moving speed of the mobile station is continuously measured, it is only necessary to return to the step 202 after the step 206 is executed each time, and details are not described here.
Fig. 2 is a schematic flow chart of a method for determining a moving speed of a mobile station by using two rows of pilot signals, and further, in order to measure the moving speed of the mobile station more accurately, the moving speed of the mobile station may be obtained by solving a correlation included angle value of a corresponding channel vector by using the multiple rows of pilot signals. Wherein each time corresponds to a row of pilot subcarriers.
When the mobile station movement speed is determined by using the multi-column pilot signals, namely, the relevant included angle values of the corresponding channel vectors of the multi-column pilot subcarriers in the set time interval are respectively calculated in each downlink subframe. Correspondingly, on the basis of steps 201-206, as shown in fig. 4, the method for determining the moving speed of the mobile station by using the multi-column pilot channels comprises the following steps:
step 401: according to the actual transmission environment or simulation estimation, the corresponding relation between the relevant included angle value of the channel vector and the movement speed of the mobile station in the set time interval of every two rows of pilot frequency sub-carriers is obtained.
Each two rows of pilot subcarriers are set with a time interval, and the time intervals can be different, so that the corresponding relation between the relevant included angle value of the channel vector and the moving speed of the mobile station in each set time interval needs to be obtained respectively.
The method for specifically obtaining the corresponding relationship in each time interval is the same as that in step 201.
In this step, as in step 201, a certain row of pilot subcarriers may be selected first, a time interval between every two pilot subcarriers is calculated, then a part or all of the calculated time interval is selected as the selected time interval, and two rows of pilot subcarriers corresponding to each selected time interval are determined. When selecting the time intervals, the time intervals with larger numerical difference should be selected. When the time interval τ is selected, the corresponding pilot subcarrier is determined, and in this case, the calculation of the correlation included angle value is directly performed without selecting the pilot subcarrier in step 403.
Step 402: as in step 202.
Step 403: and respectively calculating the pilot frequency sub-carrier of each selected certain time in the current downlink sub-frame and the relevant included angle value of the corresponding channel vector of the calculated pilot frequency sub-carrier in a set time interval.
Wherein, the calculation of each related included angle value is the same as the calculation method of the related included angle value in the step 203.
Step 404: and counting multiple rows of pilot frequency sub-carriers in multiple downlink subframes and corresponding related included angle values in a set time interval, and then respectively calculating the mean value of the corresponding related included angles in each time interval.
The specific statistical method and the calculation method of the mean value of the correlation included angles are the same as those in step 204, and the differences are only that: the mean values of different correlation angles over a plurality of time intervals need to be repeatedly calculated.
Step 405: and correspondingly eliminating the relevant included angle value corresponding to each time interval respectively, and recalculating a new relevant included angle mean value corresponding to each time interval.
The method for removing and recalculating the new mean value of the correlation angle is the same as the method in step 205, and the difference is only that: the corresponding new mean correlation angles for a number of different time intervals need to be repeatedly calculated.
Step 406: and correspondingly determining the movement speed of the mobile station according to the new correlation included angle mean value corresponding to each time interval by respectively taking the corresponding relation corresponding to each time interval as the basis.
The specific processing method is the same as step 206, and the differences are only that: the velocity interval corresponding to the new mean correlation angle in each time interval needs to be determined repeatedly.
Step 407: the final moving velocity of the mobile station is calculated from the moving velocities of the plurality of mobile stations obtained in step 406.
Since each new correlation angle mean value can obtain a velocity interval according to the corresponding relationship, and the velocity intervals are not necessarily completely the same, it is necessary to perform certain processing on the obtained multiple velocity intervals, so as to obtain the final velocity interval corresponding to the moving velocity of the mobile station.
The specific method can be as follows: arranging corresponding speed intervals in the order from small to large according to the time intervals in the obtained speed intervals; then, the intersection calculation of the speed intervals is carried out: comparing whether the first speed interval and the second speed interval have intersection, if not, taking the second speed interval as an effective speed interval, and if so, taking the speed interval of the intersection as the effective speed interval; and then comparing whether the effective speed interval and the third speed interval have an intersection, if not, taking the third speed interval as the effective speed interval, if so, taking the speed interval of the intersection as the effective speed interval, and using the judgment processing method of the effective speed interval and the third speed interval to process the subsequent speed interval, and so on, until the final speed interval is compared with the effective speed interval, obtaining the final effective speed interval, wherein the final effective speed interval is the speed interval corresponding to the movement speed of the mobile station.
The above method for processing a plurality of speed sections is given, but the method is not limited to the above method for processing the speed sections, and in practical applications, the above processing method may be adaptively changed, or other processing methods for the sections may be used.
Here, steps 404 and 405 are optional steps.
When step 404 is removed, step 405 must be removed accordingly, and in step 406, the moving speed of the mobile station is determined according to the corresponding relationship of each time interval and the relevant included angle value corresponding to each time interval calculated in step 402.
When step 405 is eliminated, step 404 may be retained, and in this case, the motion speed of the mobile station is correspondingly determined in step 406 according to the correspondence and the mean correlation angle calculated in step 404.
In the process of calculating the moving speed of the mobile station in fig. 4, the correlation included angle value, the correlation included angle mean value, and the new correlation included angle mean value may be collectively referred to as a calculation result.
Fig. 5 is a schematic structural diagram of an apparatus for measuring a moving speed of a mobile station according to the present invention, as shown in fig. 5, the apparatus includes: a correspondence obtaining module 510, a correlation angle mean calculating module 520, and a motion velocity determining module 530, wherein,
a corresponding relation obtaining module 510, configured to obtain a corresponding relation between a relevant included angle value of the channel vector and a moving speed of the mobile station according to an actual transmission environment or simulation estimation, and send the obtained corresponding relation to the moving speed determining module 530.
When multiple rows of pilot signals are used, the correspondence obtaining module 510 needs to obtain correspondences in multiple fixed time intervals; when two rows of pilot signals are used, the corresponding relation obtaining module 510 only needs to obtain the corresponding relation in the time interval in which the pilot signals are set.
A mean value calculation module 520 of the correlation included angle, configured to obtain an estimated value of channel information on a pilot subcarrier in a downlink subframe, and calculate a correlation included angle value of a channel vector according to the estimated value, when a certain number of correlation included angles of the downlink subframe are obtained, calculate a mean value of the correlation included angle, compare each correlation included angle value with a mean value of the corresponding correlation included angle, remove a correlation included angle value that deviates from the mean value of the correlation included angle by a relatively large amount, and then send a new mean value of the correlation included angle, which is recalculated, to the motion speed determination module 530.
When multiple rows of pilot subcarriers are taken from each downlink subframe, the calculated correlation included angle value, the calculated correlation included angle mean value and the calculated new correlation included angle mean value are all multiple and respectively correspond to each two rows of pilot subcarriers and the time interval thereof.
The mean value of the correlation included angle calculating module 520 may also calculate only the correlation included angle value, and directly send the correlation included angle value to the motion speed determining module 530 without performing subsequent calculation of the mean value of the correlation included angle; the mean value of the correlation angle calculation module 520 may also calculate only the mean value of the correlation angle and send the mean value of the correlation angle to the motion speed determination module 530 without performing the subsequent steps of removing the correlation included angle value and recalculating the mean value of the correlation included angle.
Here, the correlation included angle value, the correlation included angle mean value, and the new correlation included angle mean value are collectively referred to as a calculation result.
A motion speed determining module 530, configured to determine a motion speed of the mobile station according to the corresponding relationship and the calculation result respectively sent by the corresponding relationship obtaining module 510 and the related included angle mean calculating module 520.
Wherein, when only two rows of pilot subcarriers are selected in each downlink subframe, the determining the motion speed of the mobile station is as follows: directly obtaining the movement speed of the mobile station according to the calculation result and the corresponding relation; when a plurality of columns of pilot subcarriers are selected in each downlink subframe, the movement speed of the mobile station is determined as follows: firstly, respectively obtaining the movement speed of the mobile station according to each calculation result and the corresponding relation; since each calculation result corresponds to a corresponding relationship and further corresponds to the movement velocity of one mobile station, after the movement velocities of the mobile stations are obtained, the obtained movement velocities of the plurality of mobile stations need to be calculated to obtain the final movement velocity of the mobile station.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention.