US20060258294A1 - Mobile station speed estimation - Google Patents
Mobile station speed estimation Download PDFInfo
- Publication number
- US20060258294A1 US20060258294A1 US11/460,427 US46042706A US2006258294A1 US 20060258294 A1 US20060258294 A1 US 20060258294A1 US 46042706 A US46042706 A US 46042706A US 2006258294 A1 US2006258294 A1 US 2006258294A1
- Authority
- US
- United States
- Prior art keywords
- envelope
- mobile station
- speed
- signal
- estimate
- 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.)
- Abandoned
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/02—Systems for determining distance or velocity not using reflection or reradiation using radio waves
- G01S11/10—Systems for determining distance or velocity not using reflection or reradiation using radio waves using Doppler effect
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/01—Reducing phase shift
Definitions
- the present invention relates to the estimation of mobile stations speed using Doppler spread estimation.
- a radio signal from a transmitter to a receiver will travel along different propagation paths as the signal is scattered by obstacles, such as houses and other objects. This leads to signals received with different time delays, so called multipath propagation.
- Doppler shifts When the receiver starts moving Doppler shifts will be introduced.
- Doppler spread or Doppler spectrum, is generated.
- the Doppler spread can also be used to infer the scalar speed of a receiver relative to local radio reflectors. This scalar speed is not the radial velocity that one can determine by measuring a Doppler shift. In most cases, the speed will correspond to the speed of the receiver, such as in the case of a receiver in a moving car surrounded by stationary street furniture.
- a method of estimating the speed of a mobile station in a MIMO, SIMO or MISO wireless communication system comprising:
- the radio signals may all be on the same carrier frequency, for instance where MIMO, SIMO or MISO techniques are employed to improve spectrum utilisation efficiency.
- MIMO, SIMO or MISO techniques are employed to improve spectrum utilisation efficiency.
- the separation of signals on different carriers can be achieved by simply filtering the received signals.
- the transfer function of a path from a transmitter to a moving receiver is subject to three amplitude effects, namely path loss, slow fading and fast fading. It is the fast fading which is associated with Doppler spread.
- the fast fading operates at time scales many orders of magnitude smaller than path loss and slow fading. Accordingly, for a signal transmitted at a constant or slowly varying power, the derivative of the envelope of the received signal will be determined predominantly by the fast fading, i.e. Doppler spread, and the path loss and slow fading effects can be disregarded.
- Modulation effects can be compensated for by demodulating the received signal with a signal corresponding to the modulating signal which can be predetermined, for example by means of reference codes transmitted at predetermined times.
- the path transfer function can be estimated by techniques, other than by demodulating the received signal with the modulating signal, such as that described in Wu, J. et al., “Blind Channel Estimation Based on Subspace for Multicarrier CDMA”, Vehicular Technology Conference, 2001. VTC 2001 Spring. IEEE VTS 53rd, Volume: 4, 2001.
- the received signal's envelope will be proportional to the transfer function of the path.
- said values for the derivatives of said envelopes are derived by low-pass filtering an envelope signal representing the corresponding path transfer function envelope to band limit it and filtering the band-limited envelope signal using an IIR or FIR filter. More preferably, said envelope signals comprise respective sequences of samples representing the corresponding path transfer function envelope.
- the computing of said estimates of the Doppler spreads comprises determining the variances of said derivative values. More preferably, the computing of said estimates of the Doppler spreads comprises determining a value indicative of the received powers of said radio signals. Still more preferably, the Doppler spread estimates are each calculated by determining the square root of the result of dividing twice the corresponding variance by the value indicative of the received power of the respective radio signal.
- the speed of the mobile station may calculated in accordance with the formula: speed ⁇ f d spread f c
- f sdpread is the Doppler spread of the signal having the greatest energy during calculation of said derivatives and f c is the carrier frequency of said radio signal.
- the speed of the mobile station may be calculated as the mean, a linearly or non-linearly weighted mean or the median of speeds calculated for a plurality of said signals in accordance with the formula: speed ⁇ f d spread f c
- d spread is the Doppler spread of a signal and f c is the carrier frequency of said signal.
- the speed estimation method may perform processing applicable to different speed ranges, e.g a low speed range and a high speed range, in parallel and then select the more appropriate speed value as the final speed estimate.
- mobile stations including processing means configured, e.g. by software, for performing methods according to the present invention.
- FIG. 1 is a block diagram of a mobile station
- FIG. 2 is a flowchart illustrating a first speed estimating program according to the present invention
- FIG. 3 is a flowchart illustrating a signal conditioning process of the program of FIG. 2 ;
- FIG. 4 is a flowchart illustrating a speed estimating process of the program of FIG. 2 ;
- FIG. 5 is a flowchart illustrating a speed estimate selection process of the program of FIG. 2 ;
- FIG. 6 illustrates a muli-input/multi-output (MIMO) radio communication system
- FIG. 7 is a block diagram of another mobile station
- FIG. 8 is a flowchart illustrating a second speed estimating program according to the present invention.
- FIG. 9 is a flowchart illustrating a signal conditioning process of the program of FIG. 8 ;
- FIG. 10 is a flowchart illustrating a speed estimating process of the program of FIG. 8 ;
- FIG. 11 is a flowchart illustrating a speed estimate selection process of the program of FIG. 8 ;
- FIG. 12 is a flowchart illustrating a third speed estimating program according to the present invention.
- FIG. 13 is a flowchart illustrating a signal conditioning process of the program of FIG. 12 ;
- FIG. 14 is a flowchart illustrating a speed estimating process of the program of FIG. 12 ;
- FIG. 15 is a flowchart illustrating a speed estimate selection process of the program of FIG. 12 ;
- a WCDMA mobile station comprises an antenna 1 , an rf subsystem 2 , a baseband DSP (digital signal processor) subsystem 3 , an analogue audio subsystem 4 , a loudspeaker 5 , a microphone 6 , a controller 7 , a liquid crystal display 8 , a keypad 9 , memory 10 , a battery 11 and a power supply circuit 12 .
- DSP digital signal processor
- the rf subsystem 2 contains if and rf circuits of the mobile telephone's transmitter and receiver and a frequency synthesizer for tuning the mobile station's transmitter and receiver.
- the antenna 1 is coupled to the rf subsystem 2 for the reception and transmission of radio waves.
- the baseband DSP subsystem 3 is coupled to the rf subsystem 2 to receive baseband signals therefrom and for sending baseband modulation signals thereto.
- the baseband DSP subsystems 3 includes codec and RAKE functionality, which are well-known in the art, and is programmed for estimating a received signal Doppler spread and the speed of the mobile station.
- the analogue audio subsystem 4 is coupled to the baseband DSP subsystem 3 and receives demodulated audio therefrom.
- the analogue audio subsystem 4 amplifies the demodulated audio and applies it to the loudspeaker 5 .
- Acoustic signals, detected by the microphone 6 are pre- amplified by the analogue audio subsystem 4 and sent to the baseband DSP subsystem 3 for coding.
- the controller 7 controls the operation of the mobile telephone. It is coupled to the rf subsystem 2 for supplying tuning instructions to the frequency synthesizer and to the baseband DSP subsystem 3 for supplying control data and management data for transmission.
- the controller 7 operates according to a program stored in the memory 10 .
- the memory 10 is shown separately from the controller 7 . However, it may be integrated with the controller 7 .
- the display device 8 is connected to the controller 7 for receiving control data and the keypad 9 is connected to the controller 7 for supplying user input data signals thereto.
- the battery 1 is connected to the power supply circuit 12 which provides regulated power at the various voltages used by the components of the mobile telephone.
- the controller 7 is programmed to control the mobile station for speech and data communication and with application programs, e.g. a WAP browser, which make use of the mobile station's data communication capabilities.
- application programs e.g. a WAP browser
- the DSP subsystem 3 is programmed to perform calculations appropriate for low and high speeds in parallel and then select the appropriate result.
- the speed estimation program of the DSP subsystem 3 receives a stream of in-phase and quadrature baseband signal components and provides time-domain samples of the transfer function (H) of the signal paths having the various delays handled by the RAKE processing.
- these transfer function samples are obtained by demodulating the received signal for reference symbols with values representing the transmitted signal.
- Signal conditioning is performed on the transfer function samples (step p 1 a and p 1 b ) using algorithms adapted for low and high speed ranges respectively.
- the low and high speed ranges overlap in the present example. However, the upper edge of the lower range may simply meet the lower edge of the upper range.
- the conditioned signals are then used to produce speed estimates (steps p 2 a and p 2 b ) and the appropriate speed estimate is selected (step p 3 ).
- the signal conditioning processes p 1 a , p 1 b loop through the RAKE fingers which process transfer function signals for respective different path delays.
- the mobile station does not employ a RAKE system, such as where the mobile station is not a CODMA device, there will be no need for this loop.
- step s 1 it is first determined whether the RAKE finger to be processed is locked (step s 1 ).
- a “locked” finger is one whose allocated lifetime has expired. If the finger is locked at step s 1 , processing for the current finger is skipped.
- step s 2 the state of the finger, i.e. deallocated or allocated to a radio path, is determined (step s 2 ).
- step s 3 it is determined whether the finger's lifetime has expired. If the finger's lifetime has expired at step s 3 , the finger is locked (step s 4 ) and further processing for it is skipped. However, if the finger's lifetime has not expired at step s 3 , all variables and parameters for the current finger are reset (step s 5 ).
- step s 6 it is determined whether a valid reference symbol has been detected.
- the reference symbols occur regularly in signals transmitted to the mobile station and have known and values. Therefore, values associated with like signals are used for calculation of the Doppler spread and speed.
- the power at which the reference signals are transmitted is constant. However, if the power at which they are transmitted varies sufficiently quickly to affect the determined derivative value for the envelope, the magnitude of the envelope can be normalised using information about their transmission power. This information may be provided in control channels, for example.
- in-phase (I) and quadrature (Q) component values for the transfer function are obtained (step s 7 ).
- I and Q values are used to calculated the magnitude of the transfer function's envelope (r) (step s 8 ).
- the previous value of transfer function's envelope magnitude is used to avoid gaps in the sequence of envelope magnitude values (step s 9 ).
- Continuous invocation mode is the normal mode of operation with an envelope value be obtained for each slot, i.e. regularly and frequently.
- the speed estimation program may not be run in some slots. For instance the program may not be called during the transmission gap in compressed mode.
- the new envelope magnitude value is fed into a low-pass filtering process (step s 1 ).
- the low-pass filtering is provided by implementing a 3 rd order Butterworth IIR filter with a cutoff at 500 Hz.
- the filtering process for the low speed thread also performs downsampling by a factor of 2.
- the bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered.
- 500 Hz corresponds approximately to a speed of 125 km/h at 2.17 GHz.
- the filtering process s 11 is skipped and the output of the filtering process s 11 is replaced with the current envelope magnitude value (step s 12 ).
- N is 750.
- step s 14 If all of the fingers have been processed at step s 14 , the signal conditioning is complete. Otherwise, the process returns to step s 1 .
- the bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered.
- 1 kHz corresponds approximately to a speed of 250 km/h at 2.17 GHz.
- the results of the low speed and high speed signal conditioning process p 1 a , p 1 b are supplied to the low speed and high speed speed estimating processes p 2 a , p 2 b respectively.
- an initialisation routine sets a “finger valid flag” to true and a “maximum energy” variable to ⁇ 1. After initialisation, the process enters a loop so that all of the RAKE fingers are processed in turn.
- step s 2 it is determined whether the current finger's lifetime is greater than or equal to a threshold (step s 2 ).
- the threshold is set to identify fingers what are sufficiently old to provide sufficiently good data for reliable speed estimation. If the current finger's lifetime is less than the threshold, the process moves on to the next finger.
- a value b1 and an energy related value for the current finger are calculated (step s 23 ).
- the energy related value for the current finger is calculated by multiplying ⁇ circumflex over (b) ⁇ 0 by the current finger's age.
- step s 24 it is determined whether the current finger has the largest energy related value of the fingers processed so far (step s 24 ) by comparing its energy with the “maximum energy” variable. If this is not the case, the process moves on to the next finger. However, if this is the case, the “finger valid flag” is set to true to indicate that data is now available for speed estimation, the “maximum energy” variable is set to the energy of the current finger and an index, identifying the current finger, is stored (step s 25 ).
- the high speed speed estimating process p 2 b is identical except for the valid value ranges used at steps s 30 and s 32 .
- the speed selection process p 3 determines whether both speed estimate process have produced reliable estimates (step s 41 ). If one or both have not produced reliable estimates, as indicated by their respective reliability flags, no speed estimate is selected, a speed estimate reliability flag is set to false (step s 42 ) and the process terminates.
- step s 44 It is then determined whether the compensated low speed speed estimate is less than a threshold located within the overlap between the low and high speed ranges (step s 44 ). If the low speed speed estimate is less than the threshold, it is determined whether the estimate is less than 0 (step s 45 ) and, if so, the final speed estimate is set to 0 (step s 46 ). If the compensated low speed speed estimate is not less than 0 at step s 45 , the final speed estimate is set to the compensated low speed speed estimate (step s 47 ).
- the final speed estimate is set to the compensated high speed speed estimate (step s 48 ).
- the speed estimate reliability flag is set to true.
- speed estimate reliability flag If the speed estimate reliability flag is true, a speed estimate will now be available for use by the mobile station and can be reported to a fixed network node where it may be used to control handovers, for example.
- a MIMO WCDMA mobile communication system comprises at least one base station 21 and at least one mobile station 22 .
- the base station has two antennas 23 , 24 and the mobile station 22 also has two antennas 25 , 26 .
- the base station 21 concurrently transmits different data from each antenna 23 , 24 using the same carrier frequency.
- the mobile station 22 uses both of its antennas 25 , 26 to receive transmissions from the base station 21 . Consequently, there are four spatially distinct channels between the base station's antennas 23 , 24 and the mobile station's antennas.
- a first spatial channel comprises the paths from the first base station antenna 23 to the first mobile station antenna 25 , which consist of a direct path and paths involving reflections from first and second buildings 27 , 28 .
- a second spatial channel comprises the paths from the first base station antenna 23 to the second mobile station antenna 26 , which consist of a direct path and paths involving reflections from first and second buildings 27 , 28 .
- a third spatial channel comprises the paths from the second base station antenna 24 to the first mobile station antenna 25 , which consist of a direct path and paths involving reflections from first and second buildings 27 , 28 .
- a fourth spatial channel comprises the paths from the second base station antenna 24 to the second mobile station antenna 26 , which consist of a direct path and paths involving reflections from first and second buildings 27 , 28 .
- the mobile station 22 additionally comprises an rf subsystem 32 , a baseband DSP (digital signal processor) subsystem 33 , an analogue audio subsystem 34 , a loudspeaker 35 , a microphone 36 , a controller 37 , a liquid crystal display 38 , a keypad 39 , memory 40 , a battery 41 and a power supply circuit 42 .
- a baseband DSP digital signal processor
- the rf subsystem 32 contains if and rf circuits of the mobile station's transmitter and receiver and a frequency synthesizer for tuning the mobile station's transmitter and receiver.
- the antennas 25 , 26 are coupled to the rf subsystem 2 for the reception and transmission of radio waves.
- the baseband DSP subsystem 33 is coupled to the rf subsystem 32 to receive baseband signals therefrom and for sending baseband modulation signals thereto.
- the baseband DSP subsystems 33 includes codec and RAKE functionality, which are well-known in the art, and is programmed for estimating a received signal Doppler spread and the speed of the mobile station 22 .
- the analogue audio subsystem 34 is coupled to the baseband DSP subsystem 33 and receives demodulated audio therefrom.
- the analogue audio subsystem 34 amplifies the demodulated audio and applies it to the loudspeaker 35 .
- Acoustic signals, detected by the microphone 36 are pre-amplified by the analogue audio subsystem 34 and sent to the baseband DSP subsystem 33 for coding.
- the controller 37 controls the operation of the mobile station 22 . It is coupled to the rf subsystem 32 for supplying tuning instructions to the frequency synthesizer and to the baseband DSP subsystem 33 for supplying control data and management data for transmission.
- the controller 37 operates according to a program stored in the memory 40 .
- the memory 40 is shown separately from the controller 37 . However, it maybe integrated with the controller 37 .
- the display device 38 is connected to the controller 37 for receiving control data and the keypad 39 is connected to the controller 7 for supplying user input data signals thereto.
- the battery 41 is connected to the power supply circuit 42 which provides regulated power at the various voltages used by the components of the mobile station 22 .
- the controller 37 is programmed to control the mobile station 22 for speech and data communication and with application programs, e.g. a WAP browser, which make use of the mobile station's data communication capabilities.
- application programs e.g. a WAP browser
- each of the mobile station's antennas 25 , 26 receives the signals transmitted from each of the base station's antennas 23 , 24 , it is necessary for the transfer functions of the various paths, followed by the signals in the four spatial channels, to be known in order to separate the transmitted signals in the mobile station 22 . Consequently, the DSP subsystem 33 is programmed to obtain channel estimates for the each of the paths. The details of this process are not material to the present invention. However, the channel estimates may be obtained by means of reference symbols in the transmitted signals. In such an example, different reference symbols are used for each transmitted signal. This enables the contributions from the different base station antennas 23 , 24 to be separated.
- the different path delays means that the reference codes are detected, e.g. by autocorrelation, at different times depending on the path followed. Consequently, the contributions of signals from the same base station antenna 23 , 24 but which followed different paths can also be separated.
- the transfer functions of the various paths can be determined once the received signals have been separated according to source antenna and path.
- the DSP subsystem 33 is programmed to perform speed estimation calculations appropriate for low and high speeds in parallel and then select the appropriate result.
- the speed estimation program of the DSP subsystem 33 receives a stream of in-phase and quadrature baseband signal components and provides time-domain samples of the transfer functions (H) of the signal paths from both base station antennas 23 , 24 having the various delays handled by the RAKE processing.
- these transfer function samples are obtained by demodulating the received signals for reference symbols with values representing the transmitted signals.
- Signal conditioning is performed on the transfer function samples (step p 11 a and p 11 b ) using algorithms adapted for low and high speed ranges respectively.
- the low and high speed ranges overlap in the present example. However, the upper edge of the lower range may simply meet the lower edge of the upper range.
- the conditioned signals are then used to produce speed estimates (steps p 12 a and p 12 b ) and the appropriate speed estimate is selected (step p 13 ).
- the signal conditioning processes p 11 a , p 11 b loop through the RAKE fingers of both RAKE processes, i.e. the RAKE process for the signals from the first base station antenna 23 and the RAKE process for the signals from the second base station antenna 24 , which process transfer function signals for respective different paths.
- RAKE process for the signals from the first base station antenna 23
- RAKE process for the signals from the second base station antenna 24 which process transfer function signals for respective different paths.
- the mobile station does not employ a RAKE system, such as where the mobile station is not a CDMA device, there will be no need for this loop.
- step s 101 it is first determined whether the RAKE finger to be processed is locked.
- a “locked” finger is one whose allocated lifetime has expired. If the finger is locked at step s 101 , processing for the current finger is skipped.
- step s 101 the state of the finger, i.e. deallocated or allocated to a radio path, is determined (step s 102 ).
- step s 103 it is determined whether the finger's lifetime has expired. If the finger's lifetime has expired at step s 103 , the finger is locked (step s 104 ) and further processing for it is skipped. However, if the finger's lifetime has not expired at step s 103 , all variables and parameters for the current finger are reset (step s 105 ).
- step s 106 it is determined whether a valid reference symbol has been detected.
- the reference symbols occur regularly in signals transmitted to the mobile station and have known values. Therefore, values associated with like signals are used for calculation of the Doppler spread and speed.
- the power at which the reference signals are transmitted is constant. However, if the power at which they are transmitted varies sufficiently quickly to affect the determined derivative value for the envelope, the magnitude of the envelope can be normalised using information about their transmission power. This information may be provided in control channels, for example.
- in-phase (I) and quadrature (Q) component values for the transfer function are obtained (step s 107 ).
- I and Q values are used to calculated the magnitude of the transfer function's envelope (r) (step s 108 ).
- the previous value of transfer function's envelope magnitude is used to avoid gaps in the sequence of envelope magnitude values (step s 109 ).
- Continuous invocation mode is the normal mode of operation with an envelope value be obtained for each slot, i.e. regularly and frequently.
- the speed estimation program may not be run in some slots. For instance the program may not be called during the transmission gap in compressed mode.
- the new envelope magnitude value is fed into a low-pass filtering process (step s 111 ).
- the low-pass filtering is provided by implementing a 3 rd order Butterworth IIR filter with a cutoff at 500 Hz.
- the filtering process for the low speed thread also performs downsampling by a factor of 2.
- the bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered.
- 500 Hz corresponds approximately to a speed of 125 km/h at 2.17 GHz.
- the filtering process s 111 is skipped and the output of the filtering process s 111 is replaced with the current envelope magnitude value (step s 112 ).
- N is 750.
- step s 114 If all of the fingers have been processed at step s 114 , the signal conditioning is complete. Otherwise, the process returns to step s 101 .
- the bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered.
- 1 kHz corresponds approximately to a speed of 250 km/h at 2.17 GHz.
- the results of the low speed and high speed signal conditioning process p 11 a , p 11 b are supplied to the low speed and high speed speed estimating processes p 12 a , p 12 b respectively.
- an initialisation routine sets a “finger valid flag” to true and a “maximum energy” variable to ⁇ 1. After initialisation, the process enters a loop so that all of the RAKE fingers are processed in turn.
- step s 122 it is determined whether the current finger's lifetime is greater than or equal to a threshold.
- the threshold is set to identify fingers what are sufficiently old to provide sufficiently good data for reliable speed estimation. If the current finger's lifetime is less than the threshold, the process moves on to the next finger.
- a value ⁇ circumflex over (b) ⁇ 0 and an energy related value for the current finger are calculated (step s 123 ).
- the energy related value for the current finger is calculated by multiplying 10 by the current finger's age.
- step s 124 it is determined whether the current finger has the largest energy related value of the fingers processed so far (step s 124 ) by comparing its energy with the “maximum energy” variable. If this is not the case, the process moves on to the next finger. However, if this is the case, the “finger valid flag” is set to true to indicate that data is now available for speed estimation, the “maximum energy” variable is set to the energy of the current finger and an index, identifying the current finger, is stored (step s 125 ).
- the high speed speed estimating process p 12 b is identical except for the valid value ranges used at steps s 130 and s 132 .
- the speed selection process p 13 determines whether both speed estimate process have produced reliable estimates (step s 141 ). If one or both have not produced reliable estimates, as indicated by their respective reliability flags, no speed estimate is selected, a speed estimate reliability flag is set to false (step s 142 ) and the process terminates.
- step s 144 It is then determined whether the compensated low speed speed estimate is less than a threshold located within the overlap between the low and high speed ranges (step s 144 ). If the low speed speed estimate is less than the threshold, it is determined whether the estimate is less than 0 (step s 145 ) and, if so, the final speed estimate is set to 0 (step s 146 ). If the compensated low speed speed estimate is not less than 0 at step s 145 , the final speed estimate is set to the compensated low speed speed estimate (step s 147 ).
- the final speed estimate is set to the compensated high speed speed estimate (step s 148 ).
- the speed estimate reliability flag is set to true.
- speed estimate reliability flag If the speed estimate reliability flag is true, a speed estimate will now be available for use by the mobile station and can be reported to a fixed network node where it may be used to control handovers, for example.
- the speed estimate is determined from the signal received with the highest energy.
- the mobile station 22 is physically as described with reference to FIG. 7 but differs in the RAKE processing and the program for calculating the speed estimate.
- the RAKE processing differs in that a separate RAKE process is provided for each antenna combination. That is, a first RAKE process processes signals from the first base station antenna 23 which are received at the first mobile station antenna 25 , a second RAKE process processes signals from the second base station antenna 24 which are received at the first mobile station antenna 25 , a third RAKE process processes signals from the first base station antenna 23 which are received at the second mobile station antenna 26 and a fourth RAKE process processes signals from the second base station antenna 24 which are received at the second mobile station antenna 26 .
- the speed estimation program of the DSP subsystem 33 receives a stream of in-phase and quadrature baseband signal components and provides time-domain samples of the transfer functions (H) of the signal paths from both base station antennas 23 , 24 having the various delays handled by the four RAKE processes.
- these transfer function samples are obtained by demodulating the received signals for reference symbols with values representing the transmitted signals.
- Signal conditioning is performed on the transfer function samples (step p 21 a and p 21 b ) using algorithms adapted for low and high speed ranges respectively.
- the low and high speed ranges overlap in the present example. However, the upper edge of the lower range may simply meet the lower edge of the upper range.
- the conditioned signals are then used to produce speed estimates (steps p 22 a and p 22 b ) and the appropriate speed estimate is selected (step p 23 ).
- the signal conditioning processes p 21 a , p 21 b loop through the RAKE fingers of each RAKE process in four parallel threads.
- step s 201 it is first determined whether the RAKE finger to be processed is locked.
- a “locked” finger is one whose allocated lifetime has expired. If the finger is locked at step s 201 , processing for the current finger is skipped.
- the state of the finger i.e. deallocated or allocated to a radio path, is determined (step s 202 ).
- step s 203 it is determined whether the finger's lifetime has expired. If the finger's lifetime has expired at step s 203 , the finger is locked (step s 204 ) and further processing for it is skipped. However, if the finger's lifetime has not expired at step s 203 , all variables and parameters for the current finger are reset (step s 205 ).
- step s 206 it is determined whether a valid reference symbol has been detected.
- the reference symbols occur regularly in signals transmitted to the mobile station and have known values. Therefore, values associated with like signals are used for calculation of the Doppler spread and speed.
- the power at which the reference signals are transmitted is constant. However, if the power at which they are transmitted varies sufficiently quickly to affect the determined derivative value for the envelope, the magnitude of the envelope can be normalised using information about their transmission power. This information may be provided in control channels, for example.
- in-phase (I) and quadrature (Q) component values for the transfer function are obtained (step s 207 ).
- I and Q values are used to calculated the is magnitude of the transfer function's envelope (r) (step s 208 ).
- the previous value of transfer function's envelope magnitude is used to avoid gaps in the sequence of envelope magnitude values (step s 209 ).
- Continuous invocation mode is the normal mode of operation with an envelope value be obtained for each slot, i.e. regularly and frequently.
- the speed estimation program may not be run in some slots. For instance the program may not be called during the transmission gap in compressed mode.
- the new envelope magnitude value is fed into a low-pass filtering process (step s 211 ).
- the low-pass filtering is provided by implementing a 3 rd order Butterworth IIR filter with a cutoff at 500 Hz.
- the filtering process for the low speed thread also performs downsampling by a factor of 2.
- the bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered.
- 500 Hz corresponds approximately to a speed of 125 km/h at 2.17 GHz.
- the filtering process s 211 is skipped and the output of the filtering process s 211 is replaced with the current envelope magnitude value (step s 212 ).
- N is 750.
- step s 214 If all of the fingers have been processed at step s 214 , the signal conditioning is complete. Otherwise, the process returns to step s 201 .
- the bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered.
- 1 kHz corresponds approximately to a speed of 250 km/h at 2.17 GHz.
- the four results of the low speed and high speed signal conditioning processes p 21 a , p 21 b are supplied to respective threads of the low speed and high speed speed estimating processes p 22 a , p 22 b.
- an initialisation routine sets a “finger valid flag” to true and a “maximum energy” variable to ⁇ 1. After initialisation, the thread enters a loop so that all of the RAKE fingers of the associated RAKE process are processed in turn.
- step s 222 it is determined whether the current finger's lifetime is greater than or equal to a threshold.
- the threshold is set to identify fingers what are sufficiently old to provide sufficiently good data for reliable speed estimation. If the current finger's lifetime is less than the threshold, the process moves on to the next finger.
- the energy related value for the current finger is calculated by multiplying ⁇ circumflex over (b) ⁇ 0 by the current finger's age.
- step s 223 it is determined whether the current finger has the largest energy related value of the fingers processed so far (step s 224 ) by comparing its energy with the “maximum energy” variable. If this is not the case, the process moves on to the next finger. However, if this is the case, the “finger valid flag” is set to true to indicate that data is now available for speed estimation, the “maximum energy” variable is set to the energy of the current finger and an index, identifying the current finger, is stored (step s 225 ).
- a mean or a linearly or non-linearly weighted sum, weighted according to E ⁇ r 2 ⁇ , of the reliable speed values generated by the threads, if any, is calculated to produce the final low speed speed estimate.
- the speed value is generated by the thread corresponding to the maximum value of E ⁇ r 2 ⁇ .
- the median of the, if any, reliable speed value can be taken a the final low speed speed estimate. If at least one thread produced a reliable estimate, a low speed reliability flag is set to true. Otherwise, the low speed reliability flag is set to false.
- the high speed speed estimating process p 22 b is identical except for the valid value ranges used at steps s 230 and s 232 and the setting of a high speed reliability flag.
- the speed selection process p 23 determines whether both speed estimate process have produced reliable estimates by inspecting the low and high speed reliability flags (step s 241 ). If one or both have not produced reliable estimates, as indicated by their respective reliability flags, no speed estimate is selected, a speed estimate reliability flag is set to false (step s 242 ) and the process terminates.
- step s 244 It is then determined whether the compensated low speed speed estimate is less than a threshold located within the overlap between the low and high speed ranges (step s 244 ). If the low speed speed estimate is less than the threshold, it is determined whether the estimate is less than 0 (step s 245 ) and, if so, the final speed estimate is set to 0 (step s 246 ). If the compensated low speed speed estimate is not less than 0 at step s 245 , the final speed estimate is set to the compensated low speed speed estimate (step s 247 ).
- the final speed estimate is set to the compensated high speed speed estimate (step s 248 ).
- the speed estimate reliability flag is set to true.
- speed estimate reliability flag If the speed estimate reliability flag is true, a speed estimate will now be available for use by the mobile station and can be reported to a fixed network node where it may be used to control handovers, for example.
- neural networks are used to generate the final low speed and high speed speed estimates from the value produced by the four threads in each of the low speed and high speed estimation processes p 22 a , p 22 b .
- the neural networks are trained from simulations sweeping over relevant signal to noise ratios and various radio channel configurations.
- both the base station and the mobile station have two antennas.
- other numbers of antennas at the base and mobile station may be used.
- the present invention can also be applied in SIMO (Single Input Multi Output) and MISO (Multi Input Single Output) systems where the same problem of multiple concurrent signals on the same carrier arises.
- SIMO Single Input Multi Output
- MISO Multi Input Single Output
- speed estimates are determined for two speed ranges. However, calculations may be performed for more or fewer speed ranges as required. More ranges give better performance but at the cost of increased complexity. Also, higher order low-pass filters give better performance but again at the cost of increased complexity.
- the Doppler spread and speed calculations are performed in a mobile station.
- these calculations could be performed at a fixed network node using signal strength reports from a mobile station.
- the calculations could be performed in part by a mobile station, e.g. the signal conditioning processes, and in part by a fixed network node, e.g. the speed estimation and selection processes.
Abstract
The Doppler spreads of the received signals in a MIMO, SIMO or MISO mobile communication system are estimated from the derivatives of the path transfer function envelopes for the received signals. The estimated Doppler spreads is used to estimate the speed of the mobile station.
Description
- This is a continuation-in-part application of application Ser. No. 10/073,909, filed Feb. 14th, 2002.
- The present invention relates to the estimation of mobile stations speed using Doppler spread estimation.
- A radio signal from a transmitter to a receiver will travel along different propagation paths as the signal is scattered by obstacles, such as houses and other objects. This leads to signals received with different time delays, so called multipath propagation. When the receiver starts moving Doppler shifts will be introduced. When several paths arrive at the same time delay a Doppler spread, or Doppler spectrum, is generated.
- Knowledge of the Doppler spread has been found to be useful for enhancing the operation of receivers as described, for example, in Morelli, M. et al., “Further Results in Carrier Frequency Estimation for Transmissions Over Flat Fading Channels”, IEEE Communications Letters vol. 2, no. 12, December 1998. The Doppler spread can also be used to infer the scalar speed of a receiver relative to local radio reflectors. This scalar speed is not the radial velocity that one can determine by measuring a Doppler shift. In most cases, the speed will correspond to the speed of the receiver, such as in the case of a receiver in a moving car surrounded by stationary street furniture.
- The present inventors have established that the Doppler spread associated with a transmission path can be estimated reliably from the derivative of the transmission path's transfer function. According to the present invention, there is provided a method of estimating the speed of a mobile station in a MIMO, SIMO or MISO wireless communication system, the method comprising:
-
- receiving a plurality of radio signals;
- deriving a value for the derivative of the envelope of the path transfer function for each radio signal;
- computing an estimate of the Doppler spread of each of said radio signals from said derivative values; and
- deriving a value for the speed of said mobile station from said Doppler spread estimates.
- The radio signals may all be on the same carrier frequency, for instance where MIMO, SIMO or MISO techniques are employed to improve spectrum utilisation efficiency. The separation of signals on different carriers can be achieved by simply filtering the received signals.
- The transfer function of a path from a transmitter to a moving receiver is subject to three amplitude effects, namely path loss, slow fading and fast fading. It is the fast fading which is associated with Doppler spread. The fast fading operates at time scales many orders of magnitude smaller than path loss and slow fading. Accordingly, for a signal transmitted at a constant or slowly varying power, the derivative of the envelope of the received signal will be determined predominantly by the fast fading, i.e. Doppler spread, and the path loss and slow fading effects can be disregarded. Modulation effects can be compensated for by demodulating the received signal with a signal corresponding to the modulating signal which can be predetermined, for example by means of reference codes transmitted at predetermined times.
- However, the path transfer function can be estimated by techniques, other than by demodulating the received signal with the modulating signal, such as that described in Wu, J. et al., “Blind Channel Estimation Based on Subspace for Multicarrier CDMA”, Vehicular Technology Conference, 2001. VTC 2001 Spring. IEEE VTS 53rd, Volume: 4, 2001.
- In the minimal case of an unmodulated carrier, the received signal's envelope will be proportional to the transfer function of the path.
- Preferably, therefore, said values for the derivatives of said envelopes are derived by low-pass filtering an envelope signal representing the corresponding path transfer function envelope to band limit it and filtering the band-limited envelope signal using an IIR or FIR filter. More preferably, said envelope signals comprise respective sequences of samples representing the corresponding path transfer function envelope.
- Preferably, the computing of said estimates of the Doppler spreads comprises determining the variances of said derivative values. More preferably, the computing of said estimates of the Doppler spreads comprises determining a value indicative of the received powers of said radio signals. Still more preferably, the Doppler spread estimates are each calculated by determining the square root of the result of dividing twice the corresponding variance by the value indicative of the received power of the respective radio signal.
- Preferably, the Doppler spread estimates are each calculated in accordance with the formula:
- where r is the magnitude of the corresponding radio signal.
- The speed of the mobile station may calculated in accordance with the formula:
- where fsdpread is the Doppler spread of the signal having the greatest energy during calculation of said derivatives and fc is the carrier frequency of said radio signal.
- Alternatively, the speed of the mobile station may be calculated as the mean, a linearly or non-linearly weighted mean or the median of speeds calculated for a plurality of said signals in accordance with the formula:
- where dspread is the Doppler spread of a signal and fc is the carrier frequency of said signal.
- The speed estimation method may perform processing applicable to different speed ranges, e.g a low speed range and a high speed range, in parallel and then select the more appropriate speed value as the final speed estimate.
- According to the present invention, there are also provided mobile stations including processing means configured, e.g. by software, for performing methods according to the present invention.
-
FIG. 1 is a block diagram of a mobile station; -
FIG. 2 is a flowchart illustrating a first speed estimating program according to the present invention; -
FIG. 3 is a flowchart illustrating a signal conditioning process of the program ofFIG. 2 ; -
FIG. 4 is a flowchart illustrating a speed estimating process of the program ofFIG. 2 ; -
FIG. 5 is a flowchart illustrating a speed estimate selection process of the program ofFIG. 2 ; -
FIG. 6 illustrates a muli-input/multi-output (MIMO) radio communication system; -
FIG. 7 is a block diagram of another mobile station; -
FIG. 8 is a flowchart illustrating a second speed estimating program according to the present invention; -
FIG. 9 is a flowchart illustrating a signal conditioning process of the program ofFIG. 8 ; -
FIG. 10 is a flowchart illustrating a speed estimating process of the program ofFIG. 8 ; -
FIG. 11 is a flowchart illustrating a speed estimate selection process of the program ofFIG. 8 ; -
FIG. 12 is a flowchart illustrating a third speed estimating program according to the present invention; -
FIG. 13 is a flowchart illustrating a signal conditioning process of the program ofFIG. 12 ; -
FIG. 14 is a flowchart illustrating a speed estimating process of the program ofFIG. 12 ; -
FIG. 15 is a flowchart illustrating a speed estimate selection process of the program ofFIG. 12 ; - Preferred embodiments of the present invention will now be described, by way of example, with reference to the accompanying drawings.
- Referring to
FIG. 1 , a WCDMA mobile station comprises anantenna 1, an rf subsystem 2, a baseband DSP (digital signal processor)subsystem 3, ananalogue audio subsystem 4, aloudspeaker 5, amicrophone 6, a controller 7, aliquid crystal display 8, akeypad 9,memory 10, abattery 11 and apower supply circuit 12. - The rf subsystem 2 contains if and rf circuits of the mobile telephone's transmitter and receiver and a frequency synthesizer for tuning the mobile station's transmitter and receiver. The
antenna 1 is coupled to the rf subsystem 2 for the reception and transmission of radio waves. - The
baseband DSP subsystem 3 is coupled to the rf subsystem 2 to receive baseband signals therefrom and for sending baseband modulation signals thereto. Thebaseband DSP subsystems 3 includes codec and RAKE functionality, which are well-known in the art, and is programmed for estimating a received signal Doppler spread and the speed of the mobile station. - The
analogue audio subsystem 4 is coupled to thebaseband DSP subsystem 3 and receives demodulated audio therefrom. Theanalogue audio subsystem 4 amplifies the demodulated audio and applies it to theloudspeaker 5. Acoustic signals, detected by themicrophone 6, are pre- amplified by theanalogue audio subsystem 4 and sent to thebaseband DSP subsystem 3 for coding. - The controller 7 controls the operation of the mobile telephone. It is coupled to the rf subsystem 2 for supplying tuning instructions to the frequency synthesizer and to the
baseband DSP subsystem 3 for supplying control data and management data for transmission. The controller 7 operates according to a program stored in thememory 10. Thememory 10 is shown separately from the controller 7. However, it may be integrated with the controller 7. - The
display device 8 is connected to the controller 7 for receiving control data and thekeypad 9 is connected to the controller 7 for supplying user input data signals thereto. - The
battery 1 is connected to thepower supply circuit 12 which provides regulated power at the various voltages used by the components of the mobile telephone. - The controller 7 is programmed to control the mobile station for speech and data communication and with application programs, e.g. a WAP browser, which make use of the mobile station's data communication capabilities.
- In order to estimate the speed of the mobile station, the
DSP subsystem 3 is programmed to perform calculations appropriate for low and high speeds in parallel and then select the appropriate result. - Referring to
FIG. 2 , the speed estimation program of theDSP subsystem 3 receives a stream of in-phase and quadrature baseband signal components and provides time-domain samples of the transfer function (H) of the signal paths having the various delays handled by the RAKE processing. In the present example, these transfer function samples are obtained by demodulating the received signal for reference symbols with values representing the transmitted signal. - Signal conditioning is performed on the transfer function samples (step p1 a and p1 b) using algorithms adapted for low and high speed ranges respectively. The low and high speed ranges overlap in the present example. However, the upper edge of the lower range may simply meet the lower edge of the upper range. The conditioned signals are then used to produce speed estimates (steps p2 a and p2 b) and the appropriate speed estimate is selected (step p3).
- The signal conditioning processes p1 a, p1 b loop through the RAKE fingers which process transfer function signals for respective different path delays. Of course, if the mobile station does not employ a RAKE system, such as where the mobile station is not a CODMA device, there will be no need for this loop.
- Referring to
FIG. 3 , in the low speed signal conditioning estimation process p1 a, it is first determined whether the RAKE finger to be processed is locked (step s1). A “locked” finger is one whose allocated lifetime has expired. If the finger is locked at step s1, processing for the current finger is skipped. - However, if the current finger is not locked at step s1, the state of the finger, i.e. deallocated or allocated to a radio path, is determined (step s2).
- In the case of the current finger being in the deallocated state at step s2, it is determined whether the finger's lifetime has expired (step s3). If the finger's lifetime has expired at step s3, the finger is locked (step s4) and further processing for it is skipped. However, if the finger's lifetime has not expired at step s3, all variables and parameters for the current finger are reset (step s5).
- In the case of the current finger being allocated at step s2 and following step s5, it is determined whether a valid reference symbol has been detected (step s6). The reference symbols occur regularly in signals transmitted to the mobile station and have known and values. Therefore, values associated with like signals are used for calculation of the Doppler spread and speed. In the present case, the power at which the reference signals are transmitted is constant. However, if the power at which they are transmitted varies sufficiently quickly to affect the determined derivative value for the envelope, the magnitude of the envelope can be normalised using information about their transmission power. This information may be provided in control channels, for example.
- If a valid reference symbol has been detected, in-phase (I) and quadrature (Q) component values for the transfer function are obtained (step s7). When the aforementioned I and Q values have been obtained, they are used to calculated the magnitude of the transfer function's envelope (r) (step s8).
- However, if a valid reference symbol has not been detected, the previous value of transfer function's envelope magnitude is used to avoid gaps in the sequence of envelope magnitude values (step s9).
- When the envelope magnitude value has been obtained, it is determined whether continuous invocation mode is being used (step s10). Continuous invocation mode is the normal mode of operation with an envelope value be obtained for each slot, i.e. regularly and frequently. However, under some circumstances, the speed estimation program may not be run in some slots. For instance the program may not be called during the transmission gap in compressed mode.
- If operation is continuous, the new envelope magnitude value is fed into a low-pass filtering process (step s1). The low-pass filtering is provided by implementing a 3rd order Butterworth IIR filter with a cutoff at 500 Hz. The filtering process for the low speed thread also performs downsampling by a factor of 2. The bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered. Thus, 500 Hz corresponds approximately to a speed of 125 km/h at 2.17 GHz.
- If the operation is discontinuous, the filtering process s11 is skipped and the output of the filtering process s11 is replaced with the current envelope magnitude value (step s12).
- The result of the filtering process s11 or the current envelope magnitude, as the case may be, is used to calculate the square of the envelope magnitude (rt 2), the derivative of the envelope magnitude ({dot over (r)}t) and the square of the derivative of the envelope magnitude ({dot over (r)}t 2). Since, the envelope magnitudes are time-spaced samples, the derivative is approximated in the low speed thread using:
where ˆindicates an approximate or estimated value and ΔT is the sampling interval. It can be seen that this is a two-tap FIR filter. - These values are used to update accumulated values
for the square of the envelope magnitude, the approximate derivative of the envelope magnitude and the square of the derivative of the envelope magnitude, where N is an observation window length in terms of sample (step s13). In the present example, N is 750. - If all of the fingers have been processed at step s14, the signal conditioning is complete. Otherwise, the process returns to step s1.
- The high speed signal conditioning p1 b is the same as the low speed signal conditioning except that the low-pass filter has a cut-off at 1 kHz rather than 500 kz and the envelope magnitude derivative {dot over (rt)} is approximated using:
- which is a 3-tap FIR filter having a zero coefficient for the middle tap.
- The bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered. Thus, 1 kHz corresponds approximately to a speed of 250 km/h at 2.17 GHz.
- The results of the low speed and high speed signal conditioning process p1 a, p1 b are supplied to the low speed and high speed speed estimating processes p2 a, p2 b respectively.
- Referring to
FIG. 4 , at the start of the low speed speed estimating process p2 a, an initialisation routine (step s21) sets a “finger valid flag” to true and a “maximum energy” variable to −1. After initialisation, the process enters a loop so that all of the RAKE fingers are processed in turn. - At the start of the loop, it is determined whether the current finger's lifetime is greater than or equal to a threshold (step s2). The threshold is set to identify fingers what are sufficiently old to provide sufficiently good data for reliable speed estimation. If the current finger's lifetime is less than the threshold, the process moves on to the next finger.
- If at step s22, the current finger is sufficiently mature, a value b1 and an energy related value for the current finger are calculated (step s23). b is indicative of the power of the current finger and is calculated using:
where E{r2} is the mean ofFIG. 3 . - The energy related value for the current finger is calculated by multiplying {circumflex over (b)}0 by the current finger's age.
- Following step s23, it is determined whether the current finger has the largest energy related value of the fingers processed so far (step s24) by comparing its energy with the “maximum energy” variable. If this is not the case, the process moves on to the next finger. However, if this is the case, the “finger valid flag” is set to true to indicate that data is now available for speed estimation, the “maximum energy” variable is set to the energy of the current finger and an index, identifying the current finger, is stored (step s25).
- When all of the fingers have been processed (step s26), it is determined whether data suitable for speed estimation has been obtained by checking the state of the “finger valid flag” (step s27). If the “finger valid flag” is false, no further calculations are performed and the result of the low speed speed estimating process p2 a is noted as being unreliable by setting a reliability flag to false. However, if the “finger valid flag” is true, the Doppler spread {circumflex over (f)}d spread is estimated for the last finger whose index was stored at step s25 using:
where E{{dot over ({circumflex over (r)})}2} is the mean of
calculated at step s1I ofFIG. 3 and E{{dot over ({circumflex over (r)})}} is the mean of
calculated at step s11 ofFIG. 3 (step s29). - The newly calculated Doppler spread estimate is checked to determine whether it falls within a range of valid values (step s30). If the Doppler spread estimate is outside the valid values range, the process proceeds to step s28. However, if it is within the valid values range, it is used to computer a speed estimate for the mobile station using:
where c is the velocity of light in free space and fc is the carrier frequency of the signal being received. It is determined whether the newly estimated speed value is within a valid range, i.e. below 500 km/h in the present example. (step s32). If not, the process moves to step s28. If the speed value is valid, the reliability flag is set to true. - The high speed speed estimating process p2 b is identical except for the valid value ranges used at steps s30 and s32.
- Referring to
FIG. 5 , when the speed estimation processes p2 a, p2 b have been completed, the speed selection process p3 determines whether both speed estimate process have produced reliable estimates (step s41). If one or both have not produced reliable estimates, as indicated by their respective reliability flags, no speed estimate is selected, a speed estimate reliability flag is set to false (step s42) and the process terminates. - If both speed estimation processes p2 a, p2 b produced valid speed estimates, the estimates are compensated for systematic errors introduced by the low-pass filtering and the speed estimation algorithm for low speeds (step s43). The compensation algorithm is a standard linear mapping of the form:
{circumflex over (ν)} compensated=α({circumflex over (ν)}+β)
and in the present example α=1.7 and β=−12 for the low speed speed estimate and α=1.1 and β=−20 for the high speed speed estimate. - It is then determined whether the compensated low speed speed estimate is less than a threshold located within the overlap between the low and high speed ranges (step s44). If the low speed speed estimate is less than the threshold, it is determined whether the estimate is less than 0 (step s45) and, if so, the final speed estimate is set to 0 (step s46). If the compensated low speed speed estimate is not less than 0 at step s45, the final speed estimate is set to the compensated low speed speed estimate (step s47).
- If the compensated low speed speed estimate is not below the threshold at step s44, the final speed estimate is set to the compensated high speed speed estimate (step s48).
- Following steps s46, s47 and s48, the speed estimate reliability flag is set to true.
- If the speed estimate reliability flag is true, a speed estimate will now be available for use by the mobile station and can be reported to a fixed network node where it may be used to control handovers, for example.
- It will be appreciated that many modifications may be made to the preferred embodiment described above. For example, additional loops may be introduced into the signal conditioning and speed estimating processes for iterating through different signal sources in the case of transmit diversity being used.
- Referring to
FIG. 6 , a MIMO WCDMA mobile communication system comprises at least onebase station 21 and at least onemobile station 22. The base station has twoantennas mobile station 22 also has twoantennas - The
base station 21 concurrently transmits different data from eachantenna mobile station 22 uses both of itsantennas base station 21. Consequently, there are four spatially distinct channels between the base station'santennas FIG. 6 , a first spatial channel comprises the paths from the firstbase station antenna 23 to the firstmobile station antenna 25, which consist of a direct path and paths involving reflections from first andsecond buildings base station antenna 23 to the secondmobile station antenna 26, which consist of a direct path and paths involving reflections from first andsecond buildings base station antenna 24 to the firstmobile station antenna 25, which consist of a direct path and paths involving reflections from first andsecond buildings base station antenna 24 to the secondmobile station antenna 26, which consist of a direct path and paths involving reflections from first andsecond buildings - Referring to
FIG. 7 , themobile station 22 additionally comprises anrf subsystem 32, a baseband DSP (digital signal processor)subsystem 33, ananalogue audio subsystem 34, aloudspeaker 35, amicrophone 36, acontroller 37, aliquid crystal display 38, akeypad 39,memory 40, abattery 41 and apower supply circuit 42. - The
rf subsystem 32 contains if and rf circuits of the mobile station's transmitter and receiver and a frequency synthesizer for tuning the mobile station's transmitter and receiver. Theantennas - The
baseband DSP subsystem 33 is coupled to therf subsystem 32 to receive baseband signals therefrom and for sending baseband modulation signals thereto. Thebaseband DSP subsystems 33 includes codec and RAKE functionality, which are well-known in the art, and is programmed for estimating a received signal Doppler spread and the speed of themobile station 22. - The
analogue audio subsystem 34 is coupled to thebaseband DSP subsystem 33 and receives demodulated audio therefrom. Theanalogue audio subsystem 34 amplifies the demodulated audio and applies it to theloudspeaker 35. Acoustic signals, detected by themicrophone 36, are pre-amplified by theanalogue audio subsystem 34 and sent to thebaseband DSP subsystem 33 for coding. - The
controller 37 controls the operation of themobile station 22. It is coupled to therf subsystem 32 for supplying tuning instructions to the frequency synthesizer and to thebaseband DSP subsystem 33 for supplying control data and management data for transmission. Thecontroller 37 operates according to a program stored in thememory 40. Thememory 40 is shown separately from thecontroller 37. However, it maybe integrated with thecontroller 37. - The
display device 38 is connected to thecontroller 37 for receiving control data and thekeypad 39 is connected to the controller 7 for supplying user input data signals thereto. - The
battery 41 is connected to thepower supply circuit 42 which provides regulated power at the various voltages used by the components of themobile station 22. - The
controller 37 is programmed to control themobile station 22 for speech and data communication and with application programs, e.g. a WAP browser, which make use of the mobile station's data communication capabilities. - Since each of the mobile station's
antennas antennas mobile station 22. Consequently, theDSP subsystem 33 is programmed to obtain channel estimates for the each of the paths. The details of this process are not material to the present invention. However, the channel estimates may be obtained by means of reference symbols in the transmitted signals. In such an example, different reference symbols are used for each transmitted signal. This enables the contributions from the differentbase station antennas base station antenna - Since the reference symbols are transmitted with known characteristics, the transfer functions of the various paths can be determined once the received signals have been separated according to source antenna and path.
- In order to estimate the speed of the mobile station, the
DSP subsystem 33 is programmed to perform speed estimation calculations appropriate for low and high speeds in parallel and then select the appropriate result. - Referring to
FIG. 8 , the speed estimation program of theDSP subsystem 33 receives a stream of in-phase and quadrature baseband signal components and provides time-domain samples of the transfer functions (H) of the signal paths from bothbase station antennas - Signal conditioning is performed on the transfer function samples (step p11 a and p11 b) using algorithms adapted for low and high speed ranges respectively. The low and high speed ranges overlap in the present example. However, the upper edge of the lower range may simply meet the lower edge of the upper range. The conditioned signals are then used to produce speed estimates (steps p12 a and p12 b) and the appropriate speed estimate is selected (step p13).
- The signal conditioning processes p11 a, p11 b loop through the RAKE fingers of both RAKE processes, i.e. the RAKE process for the signals from the first
base station antenna 23 and the RAKE process for the signals from the secondbase station antenna 24, which process transfer function signals for respective different paths. Of course, if the mobile station does not employ a RAKE system, such as where the mobile station is not a CDMA device, there will be no need for this loop. - Referring to
FIG. 9 , in the low speed signal conditioning estimation process p11 a, it is first determined whether the RAKE finger to be processed is locked (step s101). A “locked” finger is one whose allocated lifetime has expired. If the finger is locked at step s101, processing for the current finger is skipped. - However, if the current finger is not locked at step s101, the state of the finger, i.e. deallocated or allocated to a radio path, is determined (step s102).
- In the case of the current finger being in the deallocated state at step s102, it is determined whether the finger's lifetime has expired (step s103). If the finger's lifetime has expired at step s103, the finger is locked (step s104) and further processing for it is skipped. However, if the finger's lifetime has not expired at step s103, all variables and parameters for the current finger are reset (step s105).
- In the case of the current finger being allocated at step s102 and following step s105, it is determined whether a valid reference symbol has been detected (step s106). The reference symbols occur regularly in signals transmitted to the mobile station and have known values. Therefore, values associated with like signals are used for calculation of the Doppler spread and speed. In the present case, the power at which the reference signals are transmitted is constant. However, if the power at which they are transmitted varies sufficiently quickly to affect the determined derivative value for the envelope, the magnitude of the envelope can be normalised using information about their transmission power. This information may be provided in control channels, for example.
- If a valid reference symbol has been detected, in-phase (I) and quadrature (Q) component values for the transfer function are obtained (step s107). When the aforementioned I and Q values have been obtained, they are used to calculated the magnitude of the transfer function's envelope (r) (step s108).
- However, if a valid reference symbol has not been detected, the previous value of transfer function's envelope magnitude is used to avoid gaps in the sequence of envelope magnitude values (step s109).
- When the envelope magnitude value has been obtained, it is determined whether continuous invocation mode is being used (step s110). Continuous invocation mode is the normal mode of operation with an envelope value be obtained for each slot, i.e. regularly and frequently. However, under some circumstances, the speed estimation program may not be run in some slots. For instance the program may not be called during the transmission gap in compressed mode.
- If operation is continuous, the new envelope magnitude value is fed into a low-pass filtering process (step s111). The low-pass filtering is provided by implementing a 3rd order Butterworth IIR filter with a cutoff at 500 Hz. The filtering process for the low speed thread also performs downsampling by a factor of 2. The bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered. Thus, 500 Hz corresponds approximately to a speed of 125 km/h at 2.17 GHz.
- If the operation is discontinuous, the filtering process s111 is skipped and the output of the filtering process s111 is replaced with the current envelope magnitude value (step s112).
- The result of the filtering process s111 or the current envelope magnitude, as the case maybe, is used to calculate the square of the envelope magnitude (r1 2), the derivative of the envelope magnitude ({dot over (r)}t) and the square of the derivative of the envelope magnitude ({dot over (r)}t 2). Since, the envelope magnitudes are time-spaced samples, the derivative is approximated in the low speed thread using:
where ˆindicates an approximate or estimated value and ΔT is the sampling interval. It can be seen that this is a two-tap FIR filter. - These values are used to update accumulated values
for the square of the envelope magnitude, the approximate derivative of the envelope magnitude and the square of the derivative of the envelope magnitude, where N is an observation window length in terms of sample (step s113). In the present example, N is 750. - If all of the fingers have been processed at step s114, the signal conditioning is complete. Otherwise, the process returns to step s101.
- The high speed signal conditioning p11 b is the same as the low speed signal conditioning except that the low-pass filter has a cut-off at 1 kHz rather than 500 kHz and the envelope magnitude derivative {dot over (r)}t is approximated using:
which is a 3-tap FIR filter having a zero coefficient for the middle tap. - The bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered. Thus, 1 kHz corresponds approximately to a speed of 250 km/h at 2.17 GHz.
- The results of the low speed and high speed signal conditioning process p11 a, p11 b are supplied to the low speed and high speed speed estimating processes p12 a, p12 b respectively.
- Referring to
FIG. 10 , at the start of the low speed estimating process p12 a, an initialisation routine (step s121) sets a “finger valid flag” to true and a “maximum energy” variable to −1. After initialisation, the process enters a loop so that all of the RAKE fingers are processed in turn. - At the start of the loop, it is determined whether the current finger's lifetime is greater than or equal to a threshold (step s122). The threshold is set to identify fingers what are sufficiently old to provide sufficiently good data for reliable speed estimation. If the current finger's lifetime is less than the threshold, the process moves on to the next finger.
- If at step s122, the current finger is sufficiently mature, a value {circumflex over (b)}0 and an energy related value for the current finger are calculated (step s123). {circumflex over (b)}0. is indicative of the power of the current finger and is calculated using:
where E{r2} is the mean of
calculated at step s111 inFIG. 9 . - The energy related value for the current finger is calculated by multiplying 10 by the current finger's age.
- Following step s123, it is determined whether the current finger has the largest energy related value of the fingers processed so far (step s124) by comparing its energy with the “maximum energy” variable. If this is not the case, the process moves on to the next finger. However, if this is the case, the “finger valid flag” is set to true to indicate that data is now available for speed estimation, the “maximum energy” variable is set to the energy of the current finger and an index, identifying the current finger, is stored (step s125).
- When all of the fingers have been processed (step s126), it is determined whether data suitable for speed estimation has been obtained by checking the state of the “finger valid flag” (step s127). If the “finger valid flag” is false, no further calculations are performed and the result of the low speed estimating process p12 a is noted as being unreliable by setting a reliability flag to false. However, if the “finger valid flag” is true, the Doppler spread {circumflex over (r)}d spread is estimated for the last finger whose index was stored at step s125 using:
where E{{dot over ({circumflex over (r)})}2} is the mean of
calculated at step s111 ofFIG. 9 and E{{dot over ({circumflex over (r)})}} is the mean of
calculated at step s111 ofFIG. 9 (step s129). - The newly calculated Doppler spread estimate is checked to determine whether it falls within a range of valid values (step s130). If the Doppler spread estimate is outside the valid values range, the process proceeds to step s128. However, if it is within the valid values range, it is used to computer a speed estimate for the mobile station using:
where c is the velocity of light in free space and f is the carrier frequency of the signal being received. It is determined whether the newly estimated speed value is within a valid range, i.e. below 500 km/h in the present example. (step s132). If not, the process moves to step s128. If the speed value is valid, the reliability flag is set to true. - The high speed speed estimating process p12 b is identical except for the valid value ranges used at steps s130 and s132.
- Referring to
FIG. 11 , when the speed estimation processes p12 a, p12 b have been completed, the speed selection process p13 determines whether both speed estimate process have produced reliable estimates (step s141). If one or both have not produced reliable estimates, as indicated by their respective reliability flags, no speed estimate is selected, a speed estimate reliability flag is set to false (step s142) and the process terminates. - If both speed estimation processes p12 a, p12 b produced valid speed estimates, the estimates are compensated for systematic errors introduced by the low-pass filtering and the speed estimation algorithm for low speeds (step s143). The compensation algorithm is a standard linear mapping of the form:
{circumflex over (ν)}compensated=α({circumflex over (ν)}+β)
and in the present example α=1.7 and β=−12 for the low speed speed estimate and α=1.1 and β=−20 for the high speed speed estimate. - It is then determined whether the compensated low speed speed estimate is less than a threshold located within the overlap between the low and high speed ranges (step s144). If the low speed speed estimate is less than the threshold, it is determined whether the estimate is less than 0 (step s145) and, if so, the final speed estimate is set to 0 (step s146). If the compensated low speed speed estimate is not less than 0 at step s145, the final speed estimate is set to the compensated low speed speed estimate (step s147).
- If the compensated low speed speed estimate is not below the threshold at step s144, the final speed estimate is set to the compensated high speed speed estimate (step s148).
- Following steps s146, s147 and s148, the speed estimate reliability flag is set to true.
- If the speed estimate reliability flag is true, a speed estimate will now be available for use by the mobile station and can be reported to a fixed network node where it may be used to control handovers, for example.
- Thus, the speed estimate is determined from the signal received with the highest energy.
- In another MIMO system embodiment, the
mobile station 22 is physically as described with reference toFIG. 7 but differs in the RAKE processing and the program for calculating the speed estimate. - The RAKE processing differs in that a separate RAKE process is provided for each antenna combination. That is, a first RAKE process processes signals from the first
base station antenna 23 which are received at the firstmobile station antenna 25, a second RAKE process processes signals from the secondbase station antenna 24 which are received at the firstmobile station antenna 25, a third RAKE process processes signals from the firstbase station antenna 23 which are received at the secondmobile station antenna 26 and a fourth RAKE process processes signals from the secondbase station antenna 24 which are received at the secondmobile station antenna 26. - Referring to
FIG. 12 , the speed estimation program of theDSP subsystem 33 receives a stream of in-phase and quadrature baseband signal components and provides time-domain samples of the transfer functions (H) of the signal paths from bothbase station antennas - Signal conditioning is performed on the transfer function samples (step p21 a and p21 b) using algorithms adapted for low and high speed ranges respectively. The low and high speed ranges overlap in the present example. However, the upper edge of the lower range may simply meet the lower edge of the upper range. The conditioned signals are then used to produce speed estimates (steps p22 a and p22 b) and the appropriate speed estimate is selected (step p23).
- The signal conditioning processes p21 a, p21 b loop through the RAKE fingers of each RAKE process in four parallel threads.
- Referring to
FIG. 13 , in each thread of the low speed signal conditioning estimation process p21 a, it is first determined whether the RAKE finger to be processed is locked (step s201). A “locked” finger is one whose allocated lifetime has expired. If the finger is locked at step s201, processing for the current finger is skipped. - However, if the current finger is not locked at step s201, the state of the finger, i.e. deallocated or allocated to a radio path, is determined (step s202).
- In the case of the current finger being in the deallocated state at step s202, it is determined whether the finger's lifetime has expired (step s203). If the finger's lifetime has expired at step s203, the finger is locked (step s204) and further processing for it is skipped. However, if the finger's lifetime has not expired at step s203, all variables and parameters for the current finger are reset (step s205).
- In the case of the current finger being allocated at step s202 and following step s205, it is determined whether a valid reference symbol has been detected (step s206). The reference symbols occur regularly in signals transmitted to the mobile station and have known values. Therefore, values associated with like signals are used for calculation of the Doppler spread and speed. In the present case, the power at which the reference signals are transmitted is constant. However, if the power at which they are transmitted varies sufficiently quickly to affect the determined derivative value for the envelope, the magnitude of the envelope can be normalised using information about their transmission power. This information may be provided in control channels, for example.
- If a valid reference symbol has been detected, in-phase (I) and quadrature (Q) component values for the transfer function are obtained (step s207). When the aforementioned I and Q values have been obtained, they are used to calculated the is magnitude of the transfer function's envelope (r) (step s208).
- However, if a valid reference symbol has not been detected, the previous value of transfer function's envelope magnitude is used to avoid gaps in the sequence of envelope magnitude values (step s209).
- When the envelope magnitude value has been obtained, it is determined whether continuous invocation mode is being used (step s210). Continuous invocation mode is the normal mode of operation with an envelope value be obtained for each slot, i.e. regularly and frequently. However, under some circumstances, the speed estimation program may not be run in some slots. For instance the program may not be called during the transmission gap in compressed mode.
- If operation is continuous, the new envelope magnitude value is fed into a low-pass filtering process (step s211). The low-pass filtering is provided by implementing a 3rd order Butterworth IIR filter with a cutoff at 500 Hz. The filtering process for the low speed thread also performs downsampling by a factor of 2. The bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered. Thus, 500 Hz corresponds approximately to a speed of 125 km/h at 2.17 GHz.
- If the operation is discontinuous, the filtering process s211 is skipped and the output of the filtering process s211 is replaced with the current envelope magnitude value (step s212).
- The result of the filtering process s211 or the current envelope magnitude, as the case may be, is used to calculate the square of the envelope magnitude (rt 2), the derivative of the envelope magnitude ({dot over (r)}t) and the square of the derivative of the envelope magnitude ({dot over (r)}t 2). Since, the envelope magnitudes are time-spaced samples, the derivative is approximated in the low speed thread using:
where ˆindicates an approximate or estimated value and ΔT is the sampling interval. It can be seen that this is a two-tap FIR filter. - These values are used to update accumulated values
for the square of the envelope magnitude, the approximate derivative of the envelope magnitude and the square of the derivative of the envelope magnitude, where N is an observation window length in terms of sample (step s213). In the present example, N is 750. - If all of the fingers have been processed at step s214, the signal conditioning is complete. Otherwise, the process returns to step s201.
- The high speed signal conditioning p21 b is the same as the low speed signal conditioning except that the low-pass filter has a cut-off at 1 kHz rather than 500 Hz and the envelope magnitude derivative {dot over (r)}t is approximated using:
which is a 3-tap FIR filter having a zero coefficient for the middle tap. - The bandwidth of the filter is set to twice the maximum Doppler shift that can occur in the speed range covered. Thus, 1 kHz corresponds approximately to a speed of 250 km/h at 2.17 GHz.
- The four results of the low speed and high speed signal conditioning processes p21 a, p21 b are supplied to respective threads of the low speed and high speed speed estimating processes p22 a, p22 b.
- Referring to
FIG. 14 , at the start of each thread of the low speed estimating process p22 a, an initialisation routine (step s221) sets a “finger valid flag” to true and a “maximum energy” variable to −1. After initialisation, the thread enters a loop so that all of the RAKE fingers of the associated RAKE process are processed in turn. - At the start of the thread's loop, it is determined whether the current finger's lifetime is greater than or equal to a threshold (step s222). The threshold is set to identify fingers what are sufficiently old to provide sufficiently good data for reliable speed estimation. If the current finger's lifetime is less than the threshold, the process moves on to the next finger.
- If at step s222, the current finger is sufficiently mature, a value {circumflex over (b)}0 and an energy related value for the current finger are calculated (step s223). {circumflex over (b)}0 is indicative of the power of the current finger and is calculated using:
where E{r2} is the mean of
calculated at step s211 inFIG. 13 . - The energy related value for the current finger is calculated by multiplying {circumflex over (b)}0 by the current finger's age.
- Following step s223, it is determined whether the current finger has the largest energy related value of the fingers processed so far (step s224) by comparing its energy with the “maximum energy” variable. If this is not the case, the process moves on to the next finger. However, if this is the case, the “finger valid flag” is set to true to indicate that data is now available for speed estimation, the “maximum energy” variable is set to the energy of the current finger and an index, identifying the current finger, is stored (step s225).
- When all of the fingers have been processed (step s226), it is determined whether data suitable for speed estimation has been obtained by checking the state of the “finger valid flag” (step s227). If the “finger valid flag” is false, no further calculations are performed and the result of the low speed estimating process p22 a is noted as being unreliable by setting a reliability flag to false. However, if the “finger valid flag” is true, the Doppler spread {circumflex over (f)}d spread is estimated for the last finger whose index was stored at step s225 using:
where E{{dot over ({circumflex over (r)})}2} is the mean of
calculated at step s211 ofFIG. 9 and E{{dot over ({circumflex over (r)})}} is the mean of
calculated at step s211 ofFIG. 13 (step s129). - The newly calculated Doppler spread estimate is checked to determine whether it falls within a range of valid values (step s230). If the Doppler spread estimate is outside the valid values range, the process proceeds to step s228. However, if it is within the valid values range, it is used to computer a speed estimate for the mobile station using:
where c is the velocity of light in free space and f is the carrier frequency of the signal being received. It is determined whether the newly estimated speed value is within a valid range, i.e. below 500 km/h in the present example, (step s232). If not, the process moves to step s228. If the speed value is valid, the reliability flag is set to true. - When the four threads have terminated, a mean or a linearly or non-linearly weighted sum, weighted according to E{r2}, of the reliable speed values generated by the threads, if any, is calculated to produce the final low speed speed estimate. Alternatively the speed value is generated by the thread corresponding to the maximum value of E{r2}. Alternatively, the median of the, if any, reliable speed value can be taken a the final low speed speed estimate. If at least one thread produced a reliable estimate, a low speed reliability flag is set to true. Otherwise, the low speed reliability flag is set to false.
- The high speed speed estimating process p22 b is identical except for the valid value ranges used at steps s230 and s232 and the setting of a high speed reliability flag.
- Referring to
FIG. 15 , when the speed estimation processes p22 a, p22 b have been completed, the speed selection process p23 determines whether both speed estimate process have produced reliable estimates by inspecting the low and high speed reliability flags (step s241). If one or both have not produced reliable estimates, as indicated by their respective reliability flags, no speed estimate is selected, a speed estimate reliability flag is set to false (step s242) and the process terminates. - If both speed estimation processes p22 a, p22 b produced valid speed estimates, the estimates are compensated for systematic errors introduced by the low-pass filtering and the speed estimation algorithm for low speeds (step s243). The compensation algorithm is a standard linear mapping of the form:
νcompensated=α({circumflex over (ν)}+β)
and in the present example α=1.7 and β=−12 for the low speed speed estimate and α=1.1 and β=−20 for the high speed speed estimate. - It is then determined whether the compensated low speed speed estimate is less than a threshold located within the overlap between the low and high speed ranges (step s244). If the low speed speed estimate is less than the threshold, it is determined whether the estimate is less than 0 (step s245) and, if so, the final speed estimate is set to 0 (step s246). If the compensated low speed speed estimate is not less than 0 at step s245, the final speed estimate is set to the compensated low speed speed estimate (step s247).
- If the compensated low speed speed estimate is not below the threshold at step s244, the final speed estimate is set to the compensated high speed speed estimate (step s248).
- Following steps s246, s247 and s248, the speed estimate reliability flag is set to true.
- If the speed estimate reliability flag is true, a speed estimate will now be available for use by the mobile station and can be reported to a fixed network node where it may be used to control handovers, for example.
- In another MIMO system embodiment, neural networks are used to generate the final low speed and high speed speed estimates from the value produced by the four threads in each of the low speed and high speed estimation processes p22 a, p22 b. The neural networks are trained from simulations sweeping over relevant signal to noise ratios and various radio channel configurations.
- In the foregoing MIMO system embodiments, both the base station and the mobile station have two antennas. However, other numbers of antennas at the base and mobile station may be used.
- The present invention can also be applied in SIMO (Single Input Multi Output) and MISO (Multi Input Single Output) systems where the same problem of multiple concurrent signals on the same carrier arises. In each case, there are multiple paths and one path is selected by a process substantially as described above to provide the Doppler estimate.
- In the foregoing, speed estimates are determined for two speed ranges. However, calculations may be performed for more or fewer speed ranges as required. More ranges give better performance but at the cost of increased complexity. Also, higher order low-pass filters give better performance but again at the cost of increased complexity.
- In the foregoing, the Doppler spread and speed calculations are performed in a mobile station. However, these calculations could be performed at a fixed network node using signal strength reports from a mobile station. Alternatively, the calculations could be performed in part by a mobile station, e.g. the signal conditioning processes, and in part by a fixed network node, e.g. the speed estimation and selection processes.
Claims (50)
1-48. (canceled)
49. A method of estimating the Doppler spread of a radio signal, the method comprising:
receiving a radio signal;
deriving values for a derivative of an envelope of a path transfer function for said radio signal; and
computing an estimate of the Doppler spread of said radio signal from the derivative values.
50. A method according to claim 49 , wherein said values for the derivative of said envelope are derived by low-pass filtering an envelope signal representing the path transfer function envelope to band limit the envelope signal and filtering the band-limited envelope signal using a Finite Impulse Response (FIR) filter.
51. A method according to claim 50 , wherein said envelope signal comprises a sequence of samples representing the path transfer function envelope.
52. A method according to claim 49 , wherein the computing of said estimate of the Doppler spread comprises determining a variance of the derivative values.
53. A method according to claim 52 , wherein the computing of said estimate of the Doppler spread comprises determining a value indicative of the received power of said radio signal.
54. A method according to claim 53 , wherein the Doppler spread estimate is calculated by determining the square root of the result of dividing twice said variance by said value indicative of the received power of the radio signal.
55. A method according to claim 49 , wherein the Doppler spread estimate is calculated in accordance with the formula:
where r is the magnitude of the radio signal.
56. A method of estimating the speed of a mobile station in a wireless communication system, the method comprising:
receiving a radio signal;
deriving values for a derivative of an envelope of a path transfer function for said radio signal;
computing an estimate of the Doppler spread of said radio signal from the derivative values; and
deriving a value for the speed of said mobile station from said Doppler spread estimate.
57. A method according to claim 56 , wherein said values for the derivative of said envelope are derived by low-pass filtering an envelope signal representing the path transfer function envelope to band limit the envelope signal and filtering the band-limited envelope signal using a Finite Impulse Response (FIR) filter.
58. A method according to claim 57 , wherein said envelope signal comprises a sequence of samples representing the path transfer function envelope.
59. A method according to claim 57 , wherein the computing of said estimate of the Doppler spread comprises determining a variance of the derivative values.
60. A method according to claim 59 , wherein the computing of said estimate of the Doppler spread comprises determining a value indicative of the received power of said radio signal.
61. A method according to claim 60 , wherein the Doppler spread estimate is calculated by determining the square root of the result of dividing twice said variance by said value indicative of the received power of the radio signal.
62. A method according to claim 56 , wherein the Doppler spread estimate is cacluated in accordance with the formula:
where r is the magnitude of the radio signal.
63. A method according to claim 56 , wherein the speed of the mobile station is calculated in accordance with the formula:
where fdspread is the Doppler spread and fc is the carrier frequency of said radio signal.
64. A method of estimating the speed of a mobile station in a wireless communication system, the method comprising:
receiving a radio signal;
deriving first and second values for a derivative of an envelope of said radio signal;
computing first and second estimates of the Doppler spread of said radio signal from the derivative values; and
deriving a value for the speed of said mobile station from said Doppler spread estimates.
65. A method according to claim 64 , wherein said first value for the derivative of the envelope of said radio signal is derived by low-pass filtering an envelope signal representing the path transfer function envelope to band limit the envelope signal and filtering the band-limited envelope signal using a Finite Impulse Response (FIR) filter and said second value for the derivative of the envelope of said radio signal is derived by low-pass filtering an envelope signal representing the path transfer function envelope to band limit the envelope signal and filtering the band-limited envelope signal using an FIR filter, the first value being derived using a low-pass filter characteristic having a lower cut-off frequency than that of the low-pass filter characteristic used for deriving said second value.
66. A method according to claim 65 , wherein the first value for the derivative of the envelope of said radio signal is further derived by downsampling the output of the low-pass filtering.
67. A method according to claim 65 , wherein said envelope signal comprises a sequence of samples representing the path transfer function envelope.
68. A method according to claim 65 , wherein the computing of said first and second estimates of the Doppler spread comprises determining variances of said first and second derivative values, respectively.
69. A method according to claim 68 , wherein the computing of each of said estimates of the Doppler spread comprises determining a value indicative of the received power of said radio signal.
70. A method according to claim 69 , wherein the first and second Doppler spread estimates are calculated by determining the square root of the result of dividing twice the respective variance by said value indicative of the received power of the radio signal.
71. A method according to claim 64 , wherein the Doppler spread estimates are calculated in accordance with the formula:
where r is the magnitude of the radio signal.
72. A method according to claim 64 , comprising deriving first and second speed estimates values from the first and second Doppler spreads, respectively, and selecting the first or second speed estimate value in dependence on the magnitude of the first or second speed estimate value to provide said speed estimate.
73. A method according to claim 72 , wherein the first and second speed estimate values are calculated in accordance with the formula:
where fdspread is a Doppler spread and fc is the carrier frequency of said radio signal.
74. A mobile station of a wireless communication network, the mobile station including radio receiver means for receiving a radio signal and outputting a baseband signal and processing means for processing said baseband signal, wherein the processing means is configured for estimating the Doppler spread of a radio signal by deriving values for a derivative of an envelope of a path transfer function for a radio signal, received by said radio receiver means, and computing an estimate of the Doppler spread of said radio signal from the derivative values.
75. A mobile station according to claim 74 , wherein said processing means is configured such that said values for the derivative of said envelope are derived by low-pass filtering an envelope signal representing the path transfer function envelope to band limit the envelope signal and filtering the band-limited envelope signal using a Finite Impulse Response (FIR) filter.
76. A mobile station according to claim 75 , wherein the radio receiver means outputs said baseband signal as a plurality of samples.
77. A mobile station according to claim 74 , wherein the processing means is configured such that the computing of said estimate of the Doppler spread comprises determining a variance of the derivative values.
78. A mobile station according to claim 77 , wherein the processing means is configured such that the computing of said estimate of the Doppler spread comprises determining a value indicative of the received power of said radio signal.
79. A mobile station according to claim 78 , wherein the processing means is configured such that the Doppler spread estimate is calculated by determining the square root of the result of dividing twice said variance by said value indicative of the received power of the radio signal.
80. A mobile station according to claim 74 , wherein the processing means is configured such that the Doppler spread estimate is calculated in accordance with the formula:
where r is the magnitude of the radio signal.
81. A mobile station of a wireless communication network, the mobile station including radio receiver means for receiving a radio signal and outputting a baseband signal and processing means for processing said baseband signal, wherein the processing means is configured for estimating the speed of the mobile station by deriving values for a derivative of a path transfer function envelope for a radio signal, received by said radio receiver means, computing an estimate of the Doppler spread of said radio signal from the derivative values and deriving a value for the speed of the mobile station from said Doppler spread estimate.
82. A mobile station according to claim 81 , wherein said processing means is configured such that said values for the derivative of said envelope are derived by low-pass filtering an envelope signal representing the path transfer function envelope to band limit the envelope signal and filtering the band-limited envelope signal using a Finite Impulse Response (FIR) filter.
83. A mobile station according to claim 82 , wherein the radio receiver means outputs said baseband signal as a plurality of samples.
84. A mobile station according to claim 81 , wherein the processing means is configured such that the computing of said estimate of the Doppler spread comprises determining a variance of the derivative values.
85. A mobile station according to claim 84 , wherein the processing means is configured such that the computing of said estimate of the Doppler spread comprises determining a value indicative of the received power of said radio signal.
86. A mobile station according to claim 85 , wherein the processing means is configured such that the Doppler spread estimate is calculated by determining the square root of the result of dividing twice said variance by said value indicative of the received power of the radio signal.
87. A mobile station according to claim 81 , wherein the processing means is configured such that the Doppler spread estimate is calculated in accordance with the formula:
where r is the magnitude of the radio signal.
88. A mobile station of a wireless communication network, the mobile station including radio receiver means for receiving a radio signal and outputting a baseband signal and processing means for processing said baseband signal, wherein the processing means is configured for estimating the speed of the mobile station by deriving first and second values for a derivative of an envelope of a path transfer function envelope of a radio signal, received by said radio receiver means, computing first and second estimates of the Doppler spread of said radio signal from the derivative values and deriving a value for the speed of said mobile station from said Doppler spread estimates.
89. A mobile station according to claim 88 , wherein the processing means is configured such that said first value for the derivative of said envelope is derived by low-pass filtering an envelope signal representing the path transfer function envelope to band limit the envelope signal and filtering the band-limited envelope signal using a Finite Impulse Response (FIR) filter and said second value for the derivative of said envelope is derived by low-pass filtering an envelope signal representing the path transfer function envelope to band limit the envelope signal and filtering the band-limited envelope signal using an FIR filter, the first value being derived using a low-pass filter characteristic having a lower cut-off frequency than that of the low-pass filter characteristic used for deriving said second value.
90. A mobile station according to claim 89 , wherein the processing means is configured such that the first value for the derivative of the envelope of said radio signal is further derived by downsampling the output of the low-pass filtering.
91. A mobile station according to claim 89 , wherein the radio receiver means outputs the baseband signal as a sequence of samples.
92. A mobile station according to claim 89 , wherein the processing means is configured such that the computing of said first and second estimates of the Doppler spread comprises determining variances of said first and second derivative values, respectively.
93. A mobile station according to claim 92 , wherein the processing means is configured such that the computing of each of said estimates of the Doppler spread comprises determining a value indicative of the received power of said radio signal.
94. A mobile station according to claim 93 , wherein the processing means is configured such that the first and second Doppler spread estimates are calculated by determining the square root of the result of dividing twice the respective variance by said value indicative of the received power of the radio signal.
95. A mobile station according to claim 88 , wherein the processing means is configured such that the Doppler spread estimates are calculated in accordance with the formula:
where r is the magnitude of the radio signal.
96. A mobile station according to claim 88 , wherein the processing means is configured for deriving first and second speed estimate values from the first and second Doppler spreads, respectively, and selecting the first or second speed estimate value in dependence on the magnitude of the first or second speed estimate value to provide said speed estimate.
97. A mobile station according to claim 96 , wherein the processing means is configured such that the first and second speed estimate values are calculated in accordance with the formula:
where fdspread is the Doppler spread and fc is the carrier frequency of said radio signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/460,427 US20060258294A1 (en) | 2002-02-14 | 2006-07-27 | Mobile station speed estimation |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/073,909 US20030153274A1 (en) | 2002-02-14 | 2002-02-14 | Doppler spread estimation |
US10/639,531 US7218934B2 (en) | 2002-02-14 | 2003-08-13 | Mobile station speed estimation |
US11/460,427 US20060258294A1 (en) | 2002-02-14 | 2006-07-27 | Mobile station speed estimation |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/639,531 Division US7218934B2 (en) | 2002-02-14 | 2003-08-13 | Mobile station speed estimation |
Publications (1)
Publication Number | Publication Date |
---|---|
US20060258294A1 true US20060258294A1 (en) | 2006-11-16 |
Family
ID=37419771
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/639,531 Expired - Fee Related US7218934B2 (en) | 2002-02-14 | 2003-08-13 | Mobile station speed estimation |
US11/460,427 Abandoned US20060258294A1 (en) | 2002-02-14 | 2006-07-27 | Mobile station speed estimation |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/639,531 Expired - Fee Related US7218934B2 (en) | 2002-02-14 | 2003-08-13 | Mobile station speed estimation |
Country Status (1)
Country | Link |
---|---|
US (2) | US7218934B2 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090224931A1 (en) * | 2008-03-06 | 2009-09-10 | Research In Motion Limited | Safety for Mobile Device Users While Driving |
EP2103956A1 (en) | 2008-03-18 | 2009-09-23 | Research In Motion Limited | Estimation of the speed of a mobile device |
US20110075773A1 (en) * | 2009-09-25 | 2011-03-31 | Icera Inc. | Method and system for processing a signal |
CN108490427A (en) * | 2018-02-07 | 2018-09-04 | 浙江大学 | A kind of moving target indoor positioning and real-time tracing method |
Families Citing this family (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7218934B2 (en) * | 2002-02-14 | 2007-05-15 | Nokia Corporation | Mobile station speed estimation |
TW579636B (en) * | 2002-10-25 | 2004-03-11 | Benq Corp | Method and system for estimating movement speed of mobile phone |
US11309943B2 (en) | 2004-04-02 | 2022-04-19 | Rearden, Llc | System and methods for planned evolution and obsolescence of multiuser spectrum |
US9826537B2 (en) | 2004-04-02 | 2017-11-21 | Rearden, Llc | System and method for managing inter-cluster handoff of clients which traverse multiple DIDO clusters |
US10985811B2 (en) | 2004-04-02 | 2021-04-20 | Rearden, Llc | System and method for distributed antenna wireless communications |
US8654815B1 (en) | 2004-04-02 | 2014-02-18 | Rearden, Llc | System and method for distributed antenna wireless communications |
US10749582B2 (en) | 2004-04-02 | 2020-08-18 | Rearden, Llc | Systems and methods to coordinate transmissions in distributed wireless systems via user clustering |
US10425134B2 (en) | 2004-04-02 | 2019-09-24 | Rearden, Llc | System and methods for planned evolution and obsolescence of multiuser spectrum |
US10200094B2 (en) * | 2004-04-02 | 2019-02-05 | Rearden, Llc | Interference management, handoff, power control and link adaptation in distributed-input distributed-output (DIDO) communication systems |
US10277290B2 (en) | 2004-04-02 | 2019-04-30 | Rearden, Llc | Systems and methods to exploit areas of coherence in wireless systems |
US10187133B2 (en) | 2004-04-02 | 2019-01-22 | Rearden, Llc | System and method for power control and antenna grouping in a distributed-input-distributed-output (DIDO) network |
US11394436B2 (en) | 2004-04-02 | 2022-07-19 | Rearden, Llc | System and method for distributed antenna wireless communications |
US10886979B2 (en) | 2004-04-02 | 2021-01-05 | Rearden, Llc | System and method for link adaptation in DIDO multicarrier systems |
US9819403B2 (en) * | 2004-04-02 | 2017-11-14 | Rearden, Llc | System and method for managing handoff of a client between different distributed-input-distributed-output (DIDO) networks based on detected velocity of the client |
US11451275B2 (en) | 2004-04-02 | 2022-09-20 | Rearden, Llc | System and method for distributed antenna wireless communications |
US9312929B2 (en) | 2004-04-02 | 2016-04-12 | Rearden, Llc | System and methods to compensate for Doppler effects in multi-user (MU) multiple antenna systems (MAS) |
US8542763B2 (en) | 2004-04-02 | 2013-09-24 | Rearden, Llc | Systems and methods to coordinate transmissions in distributed wireless systems via user clustering |
US9685997B2 (en) | 2007-08-20 | 2017-06-20 | Rearden, Llc | Systems and methods to enhance spatial diversity in distributed-input distributed-output wireless systems |
EP1908195A1 (en) * | 2005-07-04 | 2008-04-09 | Nokia Corporation | Speed detection method in communication system, receiver, network element and processor |
US9097783B2 (en) | 2006-04-28 | 2015-08-04 | Telecommunication Systems, Inc. | System and method for positioning using hybrid spectral compression and cross correlation signal processing |
WO2008038520A1 (en) * | 2006-09-25 | 2008-04-03 | Nec Corporation | Mobile communication terminal and method for detecting mobile speed of the same |
US8989155B2 (en) | 2007-08-20 | 2015-03-24 | Rearden, Llc | Systems and methods for wireless backhaul in distributed-input distributed-output wireless systems |
US8160503B2 (en) * | 2009-02-27 | 2012-04-17 | Research In Motion Limited | Method and system for characterizing a radio channel of a wireless network using variability of synchronization |
US20120120942A1 (en) * | 2009-06-19 | 2012-05-17 | Sharp Kabushiki Kaisha | Wireless communication system, transmitter and wireless communication method |
CA2814009A1 (en) * | 2010-10-08 | 2012-04-12 | Loctronix Corporation | Doppler aided inertial navigation |
US11050468B2 (en) | 2014-04-16 | 2021-06-29 | Rearden, Llc | Systems and methods for mitigating interference within actively used spectrum |
US11189917B2 (en) | 2014-04-16 | 2021-11-30 | Rearden, Llc | Systems and methods for distributing radioheads |
US10194346B2 (en) | 2012-11-26 | 2019-01-29 | Rearden, Llc | Systems and methods for exploiting inter-cell multiplexing gain in wireless cellular systems via distributed input distributed output technology |
US11190947B2 (en) | 2014-04-16 | 2021-11-30 | Rearden, Llc | Systems and methods for concurrent spectrum usage within actively used spectrum |
US10164698B2 (en) | 2013-03-12 | 2018-12-25 | Rearden, Llc | Systems and methods for exploiting inter-cell multiplexing gain in wireless cellular systems via distributed input distributed output technology |
US9973246B2 (en) | 2013-03-12 | 2018-05-15 | Rearden, Llc | Systems and methods for exploiting inter-cell multiplexing gain in wireless cellular systems via distributed input distributed output technology |
US10488535B2 (en) | 2013-03-12 | 2019-11-26 | Rearden, Llc | Apparatus and method for capturing still images and video using diffraction coded imaging techniques |
US9923657B2 (en) | 2013-03-12 | 2018-03-20 | Rearden, Llc | Systems and methods for exploiting inter-cell multiplexing gain in wireless cellular systems via distributed input distributed output technology |
US10547358B2 (en) | 2013-03-15 | 2020-01-28 | Rearden, Llc | Systems and methods for radio frequency calibration exploiting channel reciprocity in distributed input distributed output wireless communications |
JP6221405B2 (en) * | 2013-06-26 | 2017-11-01 | 富士通株式会社 | Mobile terminal, speed estimation method, speed estimation program |
US11290162B2 (en) | 2014-04-16 | 2022-03-29 | Rearden, Llc | Systems and methods for mitigating interference within actively used spectrum |
US11516051B2 (en) | 2018-03-06 | 2022-11-29 | Samsung Electronics Co., Ltd. | Method and apparatus for AI-based UE speed estimation using uplink SRS measurements |
US11496339B2 (en) * | 2020-04-03 | 2022-11-08 | Samsung Electronics Co., Ltd. | Doppler spread estimation based on supervised learning |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3713154A (en) * | 1960-08-20 | 1973-01-23 | Telefunken Patent | Radar |
US5396253A (en) * | 1990-07-25 | 1995-03-07 | British Telecommunications Plc | Speed estimation |
US20020042279A1 (en) * | 2000-08-31 | 2002-04-11 | Alcatel | Receiver device for a mobile radiocommunication unit employing a speed estimator |
US6449489B1 (en) * | 1999-09-14 | 2002-09-10 | Lucent Technologies Inc. | Method and apparatus for adaptive adjustments of user probe signal |
US20020172307A1 (en) * | 2001-03-27 | 2002-11-21 | David Sandberg | Method and apparatus for estimating doppler spread |
US6505053B1 (en) * | 1999-11-05 | 2003-01-07 | At&T Corp. | Method for sinusoidal modeling and prediction of fast fading processes |
US6542745B1 (en) * | 1999-02-08 | 2003-04-01 | Mitsubishi Denki Kabushiki Kaisha | Method of estimating the speed of relative movement of a transmitter and a receiver, in communication with one another, of a telecommunication system |
US6563861B1 (en) * | 1999-03-22 | 2003-05-13 | Ericsson, Inc. | Doppler spread estimation system |
US6636574B2 (en) * | 2001-05-31 | 2003-10-21 | Motorola, Inc. | Doppler spread/velocity estimation in mobile wireless communication devices and methods therefor |
US6680969B1 (en) * | 1999-03-22 | 2004-01-20 | Ericsson, Inc. | Methods for estimating doppler spreads including autocorrelation function hypotheses and related systems and receivers |
US6697640B1 (en) * | 1998-11-11 | 2004-02-24 | Nokia Corporation | Method and apparatus for directional radio communication |
US6880969B2 (en) * | 2001-03-23 | 2005-04-19 | Denso Corporation | Temperature sensor and production method thereof |
US7054394B2 (en) * | 2002-02-18 | 2006-05-30 | Telefonaktiebolaget Lm Ericsson (Publ) | Doppler shift and spread estimation method and apparatus |
US7107056B2 (en) * | 2002-10-25 | 2006-09-12 | Benq Corporation | Method and system for estimating movement speed of mobile unit |
US7218934B2 (en) * | 2002-02-14 | 2007-05-15 | Nokia Corporation | Mobile station speed estimation |
-
2003
- 2003-08-13 US US10/639,531 patent/US7218934B2/en not_active Expired - Fee Related
-
2006
- 2006-07-27 US US11/460,427 patent/US20060258294A1/en not_active Abandoned
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3713154A (en) * | 1960-08-20 | 1973-01-23 | Telefunken Patent | Radar |
US5396253A (en) * | 1990-07-25 | 1995-03-07 | British Telecommunications Plc | Speed estimation |
US6697640B1 (en) * | 1998-11-11 | 2004-02-24 | Nokia Corporation | Method and apparatus for directional radio communication |
US6542745B1 (en) * | 1999-02-08 | 2003-04-01 | Mitsubishi Denki Kabushiki Kaisha | Method of estimating the speed of relative movement of a transmitter and a receiver, in communication with one another, of a telecommunication system |
US6680969B1 (en) * | 1999-03-22 | 2004-01-20 | Ericsson, Inc. | Methods for estimating doppler spreads including autocorrelation function hypotheses and related systems and receivers |
US6563861B1 (en) * | 1999-03-22 | 2003-05-13 | Ericsson, Inc. | Doppler spread estimation system |
US6449489B1 (en) * | 1999-09-14 | 2002-09-10 | Lucent Technologies Inc. | Method and apparatus for adaptive adjustments of user probe signal |
US6505053B1 (en) * | 1999-11-05 | 2003-01-07 | At&T Corp. | Method for sinusoidal modeling and prediction of fast fading processes |
US20020042279A1 (en) * | 2000-08-31 | 2002-04-11 | Alcatel | Receiver device for a mobile radiocommunication unit employing a speed estimator |
US6880969B2 (en) * | 2001-03-23 | 2005-04-19 | Denso Corporation | Temperature sensor and production method thereof |
US20020172307A1 (en) * | 2001-03-27 | 2002-11-21 | David Sandberg | Method and apparatus for estimating doppler spread |
US6636574B2 (en) * | 2001-05-31 | 2003-10-21 | Motorola, Inc. | Doppler spread/velocity estimation in mobile wireless communication devices and methods therefor |
US7218934B2 (en) * | 2002-02-14 | 2007-05-15 | Nokia Corporation | Mobile station speed estimation |
US7054394B2 (en) * | 2002-02-18 | 2006-05-30 | Telefonaktiebolaget Lm Ericsson (Publ) | Doppler shift and spread estimation method and apparatus |
US7107056B2 (en) * | 2002-10-25 | 2006-09-12 | Benq Corporation | Method and system for estimating movement speed of mobile unit |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090224931A1 (en) * | 2008-03-06 | 2009-09-10 | Research In Motion Limited | Safety for Mobile Device Users While Driving |
US7898428B2 (en) | 2008-03-06 | 2011-03-01 | Research In Motion Limited | Safety for mobile device users while driving |
EP2103956A1 (en) | 2008-03-18 | 2009-09-23 | Research In Motion Limited | Estimation of the speed of a mobile device |
US20090240464A1 (en) * | 2008-03-18 | 2009-09-24 | Research In Motion Limited | Estimation of the Speed of a Mobile Device |
US7895013B2 (en) | 2008-03-18 | 2011-02-22 | Research In Motion Limited | Estimation of the speed of a mobile device |
US20110075773A1 (en) * | 2009-09-25 | 2011-03-31 | Icera Inc. | Method and system for processing a signal |
US9083446B2 (en) * | 2009-09-25 | 2015-07-14 | Icera Inc. | Method and system for processing a signal |
CN108490427A (en) * | 2018-02-07 | 2018-09-04 | 浙江大学 | A kind of moving target indoor positioning and real-time tracing method |
Also Published As
Publication number | Publication date |
---|---|
US20040097197A1 (en) | 2004-05-20 |
US7218934B2 (en) | 2007-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7218934B2 (en) | Mobile station speed estimation | |
CA2224934C (en) | Rake receiver architecture for a mobile demodulator used in a cdma communication system | |
US7376176B2 (en) | Method and equipment for regulating dynamically an average area of a channel estimation | |
US7515902B2 (en) | Velocity estimation apparatus in mobile communication environments | |
US6928274B2 (en) | Receiver device for a mobile radiocommunication unit employing a speed estimator | |
EP1210803B1 (en) | Determination of data rate, based on power spectral density estimates | |
JPH08307382A (en) | Cdma rake receiver having subchip resolution | |
EP0808031B1 (en) | Spread spectrum multi-path demodulator | |
JP2002520982A (en) | Adaptive setting method of path selection threshold for DS-CDMA receiver | |
EP1212871B1 (en) | Estimation, in a rake receiver, of the rate of change of the channel | |
EP1600792A2 (en) | Velocity estimation apparatus and method | |
JPWO2004068749A1 (en) | Fading frequency estimation device | |
US20120033768A1 (en) | Non-coherent detection apparatus and method for ieee 802.15.4 lr-wpan bpsk receiver | |
US20030003923A1 (en) | Equipment for the calculation of mobile handset position | |
US8290495B2 (en) | Apparatus and method for velocity estimation in mobile communication system | |
US6118806A (en) | Signal synthesis method and apparatus under diversity reception | |
CN1111976C (en) | TSTD transmitter for limiting transmission power of antenna and controlling method thereof for base station in mobile communication system | |
JP2001036952A (en) | Cdma receiving device in cdma mobile communication system and reception signal power measuring method | |
CN101283558A (en) | Communication channel estimation | |
CN1146135C (en) | Method and apparatus for configuring Rake receiver | |
US7352799B2 (en) | Channel estimation in spread spectrum system | |
EP1443672B1 (en) | Path searching circuit, method and program in a CDMA communication system | |
US20030153274A1 (en) | Doppler spread estimation | |
US20030050017A1 (en) | Process and device for estimating the speed of movement of a mobile terminal, in particular a cellular mobile telephone | |
EP2485417A1 (en) | Method and device for selecting noise path |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |