This application claims priority to Chinese Patent Application No. 201110422765.0, filed on Dec. 16, 2011, which is hereby incorporated by reference in its entirety.
The present invention relates to the field of communications, and in particular, to a multipath searching method and a multipath searcher.
A multipath searcher is an important module of an RAKE receiver in a wideband code division multiple access (Wideband Code Division Multiple Access, WCDMA) system, and its performance directly influences performance and implementation complexity of the RAKE receiver. A function of multipath searching (or time delay estimation) is to obtain signal energy distribution characteristics at different time delay points through a matching filter, identify multipath positions having greater energy, and distribute their time amounts to different receiving paths of the RAKE receiver. The multipath searching provides information such as a multipath position and multipath strength of a wireless channel for subsequent synchronization, channel estimation, Rake reception, channel equalization, measurement, and so on.
An existing multipath searching method includes calculating steps in the following order: performing correlation and coherent accumulation on received signals and reference sequences; then performing time delay estimation; performing incoherent accumulation after a value obtained through the time delay estimation is converted into an absolute value and squared; and finally, selecting a maximum value and outputting a final time delay estimation result, where the time delay estimation specifically is:
1. Calculate correlation values of the received signals and reference signals to find a position of a global maximum value.
2. Take M discrete values in a range near the global maximum value, and perform serial estimation on each of the taken values, that is, first l−1 maximum values are multiplied by a raised cosine (Raised Cosine, RC) waveform, and subtracted from an original correlated waveform, then an lth maximum value is found, and finally L maximum value positions in each group of M groups are obtained.
3. Calculate a sum of L reconstructed waveforms, compare mean-square errors between the original waveform and the reconstructed waveforms, and select maximum value positions of the group with a minimum mean-square error as an estimated output of a multipath delay.
When multipath waveforms are reconstructed, that is, when a value of the original waveform is multiplied by an RC waveform, an ideal multipath waveform is a multiplied product of a gain of the path and the RC waveform. However, in the foregoing multipath searching method provided in the prior art, when multipath waveforms are reconstructed, a used value of the original waveform is a result of superimposing all multipath waveforms, and because influences of other multipath waveforms are included when the multipath waveforms are reconstructed, a result finally obtained has low accuracy, resulting in that multipaths with an interval less than one chip cannot be distinguished, and weak multipaths are difficult to be detected.
Embodiments of the present invention provide a multipath searching method and a multipath searcher, to improve accuracy of delay estimation, thereby improving performance of a receiver.
An embodiment of the present invention provides a multipath searching method, where the method includes: performing processing on received signals ri(t) and pilot sequences pi(t) to obtain upsampling signals y(τ);
subtracting contribution values of a current multipath waveform to other multipath waveforms by using a raised cosine function with a peak position set to zero and according to a global maximum value of |y(τ)|2 and a total of M sampling points {circumflex over (τ)}l(m) near the global maximum value, to obtain an upsampling signal yL(τ) having a minimum residual in M groups of upsampling signals yL(m)(τ), where M is an integer greater than or equal to 1, and |y(τ)|2 is a squared absolute value of y(τ); and
obtaining a time delay power spectrum of the upsampling signal yL(τ) having the minimum residual.
An embodiment of the present invention provides a multipath searcher, where the apparatus includes: a pre-processing module, configured to perform processing on received signals ri(t) and pilot sequences pi(t) to obtain upsampling signals y(τ);
a minimum residual obtaining module, configured to subtract contribution values of a current multipath waveform to other multipath waveforms by using a raised cosine function with a peak position set to zero and according to a global maximum value of |y(τ)|2 and M sampling points {circumflex over (τ)}l(m) near the global maximum value, to obtain an upsampling signal yL(τ) having a minimum residual in M groups of upsampling signals yL(m)(τ), where M is an integer greater than or equal to 1, and |y(τ)|2 is a squared absolute value of y(τ); and
a time delay power spectrum obtaining module, configured to obtain a time delay power spectrum of the upsampling signal yL(τ) having the minimum residual.
It maybe known from the foregoing embodiments of the present invention that, the contribution values of the current multipath waveform to the other multipath waveforms may be subtracted by using the raised cosine function with the peak position set to zero and according to the global maximum value of the squared absolute value |y(τ)|2 of the upsampling signal y(τ) and the M sampling points {circumflex over (τ)}l(m) near the global maximum value, to obtain the upsampling signal yL(τ) having the minimum residual in the M groups of upsampling signals yL(m)(τ), and further obtain the time delay power spectrum of the upsampling signal yL(τ) having the minimum residual. Therefore, compared with the prior art, in the method provided by the embodiments of the present invention, when multipath waveforms are reconstructed, an influence of the current multipath waveform on other multipaths is eliminated, so that multipaths with an interval less than one chip may be better distinguished, and weak multipaths may also be better detected, thereby improving the performance of the receiver.
To illustrate the technical solutions in the embodiments of the present invention more clearly, accompanying drawings required for describing the prior art or the embodiments are the briefly introduced in the following. Apparently, the accompanying drawings in the following description are merely some embodiments of the present invention, and persons of skill in the art may further obtain other drawings from these accompanying drawings.
The technical solutions in the embodiments of the present invention are described clearly and completely in the following with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the embodiments to be described are merely part rather than all of the embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by persons skilled in the art shall fall within the protection scope of the present invention.
S101: Perform processing on received signals ri(t) and pilot sequences pi(t) to obtain upsampling signals y(τ).
In the embodiment of the present invention, first, a correlator may be used to perform a correlation operation on N received signals ri(t) and pilot sequences pi(t), where the correlation operation is Σi=0N−1ri(t+τ) pi*(t) and pi*(t) is a complex conjugate of the pilot sequences pi(t). Then, coherent accumulation is performed on a result of the correlation operation Σi=0N−1ri(t+τ) pi*(t) on I pilot symbols to obtain yc(τ)=Σi=0I−1yi(τ), where I is a positive integer; and a coherent accumulator may be used to perform an operation of the coherent accumulation.
Next, upsampling (upsampling), for example, K-multiple) upsampling is performed on Σi=I−1yi(τ) obtained through the coherent accumulation, that is, K−1 zero values are inserted between every two sampling points. Afterward, filtering is performed on upsampled values obtained through the K-multiple upsampling to obtain y(τ). In the embodiment of the present invention, an interpolation filter may be used to perform filtering on the upsampled values.
S102: Subtract contribution values of a current multipath waveform to other multipath waveforms by using a raised cosine function with a peak position set to zero and according to a global maximum value of |y(τ)|2 and a total of M sampling points {circumflex over (τ)}l(m) near the global maximum value, to obtain an upsampling signal yL(τ) having a minimum residual in M groups of upsampling signals yL(m)(τ), where M is an integer greater than or equal to 1, and |y(τ)|2 is a squared absolute value of y(τ).
|y(τ)|2 is the squared absolute value of y(τ) obtained after the upsampling and the filtering are performed on Σi=0I−1yi(τ) in step S101. If a value of |y(τ)|2 is obtained, it is not difficult to obtain the global maximum value of |y(τ)|2, thereby acquiring a position of the global maximum value of |y(τ)|2. In the embodiment of the present invention, the position of the global maximum value of |y(τ)|2 is represented by {circumflex over (τ)}g.
Although the global maximum value of |y(τ)|2 and its position can be found through the foregoing method, taking influences of other multipath waveforms into consideration, the foregoing maximum value and its position are not very accurate. Therefore, in the embodiment of the present invention, M sampling points {circumflex over (τ)}l(m) are taken at and near the position {circumflex over (τ)}g which is of the global maximum value of |y(τ)|2 and found in advance, M groups of yL(m)(τ) are obtained through an iteration operation and residuals are calculated, and then the M residuals are compared to obtain the upsampling signal yL(τ) having the minimum residual in the M groups of yL(m)(τ). Specifically, the iteration operation is as follows:
S1021: Take M sampling points in [{circumflex over (τ)}g−Δτ, {circumflex over (τ)}g+Δτ.
S1022: Perform the iteration operation on each sampling point {circumflex over (τ)}l(m) of the M sampling points according to iteration formulas yl(m)(τ)=yl−1(m)(τ)−yl−1(m)({circumflex over (τ)}l−1(m))R(τ−{circumflex over (τ)}l−1(m)) and
where l takes a value of 2, 3, . . . , L, m takes a value of 1, 2, 3, . . . , M, and L and M may take empirical values of 3 and 8 respectively.
For each sampling point {circumflex over (τ)}l(m), one group yl(m)(τ) may be obtained through one iteration. According to a value of l, L yl(m)(τ) may be obtained through one iteration for each sampling point {circumflex over (τ)}l(m). In this way, when iterations for all the sampling points are finished, it is equivalent to performing M×L iterations, and taking a last iteration result yL(m)(τ) of each group of the M groups.
S1023: Calculate a residual for the last iteration result yL(m)(τ) of each group according to a formula
Finally, residuals in all the groups are compared to obtain the upsampling signal yL(τ) having the minimum residual in the M groups of yL(m)(τ).
It should be noted that, in the embodiment of the present invention, R(τ) in yl(m)(τ)=yl−1(m)(τ)−yl−1(m)({circumflex over (τ)}l−1(m))R(τ−{circumflex over (τ)}l−1(m)) is a standard raised cosine function (Raised Cosine, RC) with a peak position set to zero, that is,
where sin(.) and cos(.) represent sine and cosine functions respectively, β represents a rolloff factor, and T represents a symbol interval. Taking a WCDMA system as an example, β=0.22, and T−1/3.84e6 s (second).
S103: Obtain a time delay power spectrum of the upsampling signal yL(τ) having the minimum residual.
In the embodiment of the present invention, the time delay power spectrum of the upsampling signal yL(τ) having the minimum residual may be obtained. Specifically, incoherent accumulation is performed on zj(τ)=|yL(τ)|2 on J time slots according to z(τ)=Σj=0J−1zj(τ), where J is an integer greater than or equal to 1. Filtering is performed on z(τ), for example, α filtering is performed according to vn(τ)=ανn−1(τ)+(1−α)z(τ) to obtain the time delay power spectrum vn(τ), where α has a value range of 0<α<1. According to common sense, a maximum value of the time delay power spectrum vn(τ) may be found on the time delay power spectrum vn(τ), thereby obtaining a position of the maximum value of the time delay power spectrum vn(τ).
It maybe known from the foregoing multipath searching method provided by the embodiment of the present invention that, the contribution values of the current multipath waveform to the other multipath waveforms may be subtracted by using the raised cosine function with the peak position set to zero and according to the global maximum value of the squared absolute value |y(τ)|2 of the upsampling signal y(τ) and the total of M sampling points {circumflex over (τ)}l(m) near the global maximum value, to obtain the M groups of upsampling signals yL(m)(τ) and compare their residuals, and further obtain the time delay power spectrum of the upsampling signal y(τ) having the minimum residual. Therefore, compared with the prior art, in the method provided by the embodiment of the present invention, when multipath waveforms are reconstructed, an influence of the current multipath waveform on other multipaths is eliminated, so that multipaths with an interval less than one chip may be better distinguished, and weak multipaths may also be better detected, thereby improving performance of a receiver.
the pre-processing module 201 is configured to perform processing on received signals ri(t) and pilot sequences pi(t) to obtain upsampling signals y(τ);
the minimum residual obtaining module 202 is configured to subtract contribution values of a current multipath waveform to other multipath waveforms by using a raised cosine function with a peak position set to zero and according to a global maximum value of |y(τ)|2 and M sampling points {circumflex over (τ)}l(m) near the global maximum value, to obtain an upsampling signal yL(τ) having a minimum residual in M groups of upsampling signals yL(m)(τ), where M is an integer greater than or equal to 1, and |y(τ)|2 is a squared absolute value of y(τ); and
the time delay power spectrum obtaining module 203 is configured to obtain a time delay power spectrum of the upsampling signal yL(τ) having the minimum residual.
It should be noted that, in the foregoing implementation manner of the multipath searcher, dividing of each functional module is merely for exemplary illustration, and the foregoing functions may be assigned to different functional modules for implementation according to needs in practical applications, for example, configuration requirements of corresponding hardware or consideration of convenience of software implementation, that is, an internal structure of the multipath searcher may be divided into different functional modules to implement all or part of the functions described above. Moreover, in practical applications, corresponding functional modules in this embodiment may be implemented through corresponding hardware, and may also be accomplished through corresponding software executed by corresponding hardware. For example, the pre-processing module may be hardware having a function of executing the performing processing on received signals ri(t) and pilot sequences pi(t) to obtain upsampling signals y(τ), such as a preprocessor, and may also be a general processor or another hardware device capable of executing corresponding computer programs to accomplish the foregoing function. For another example, the minimum residual obtaining module may be hardware having a function of executing the subtracting contribution values of a current multipath waveform from the other multipath waveforms by using a raised cosine function with a peak position set to zero and according to a global maximum value of a|y(τ)|2 and M sampling points {circumflex over (τ)}l(m) near the global maximum value, to obtain an upsampling signal yL(τ) having a minimum residual in M groups of upsampling signals yL(m)(τ), such as a minimum residual obtaining device, and may also be a general processor or another hardware device capable of executing corresponding computer programs to accomplish the foregoing function (the principle described above is applicable to all the embodiments provided by this specification).
The pre-processing module 201 shown in
the correlator 301 is configured to perform correlation on the received signals ri(t) and the pilot sequences pi(t), to obtain yi(τ)=Σi=0N−1ri(t+τ) pi*(t), where pi*(t) represents a complex conjugate of the pilot sequences pi(t);
the coherent accumulator 302 is configured to perform coherent accumulation on yi(τ) on I pilot symbols to obtain yc(τ)=Σi=0I−1yi(τ), where I is a positive integer; and
the sampling and filtering unit 303 is configured to perform upsampling and filtering on yc(τ) to obtain y(τ).
The minimum residual obtaining module 202 shown in
the maximum value searching submodule 401 is configured to search |y(τ)|2 for a position {circumflex over (τ)}g of the global maximum value of |y(τ)|2;
the iteration submodule 402 is configured to take a total of M sampling points {circumflex over (τ)}l(m) at the position {circumflex over (τ)}g and near the position {circumflex over (τ)}g, obtain the M groups of upsampling signals yL(m)(τ) through an iteration operation, and calculate residuals of the M groups of yL(m)(τ); and
the comparing submodule 403 is configured to compare the M residuals, to obtain the upsampling signal yL(τ) having the minimum residual in the M groups of yL(m)(τ).
The iteration submodule 402 shown in
the sampling unit 501 is configured to take M sampling points in [{circumflex over (τ)}g−Δτ, {circumflex over (τ)}g+Δτ];
the iteration unit 502 is configured to perform an iteration operation on each sampling point {circumflex over (τ)}l(m) of the M sampling points according to formulas and yl(m)(τ)=yl−1(m)(τ)−yl−1(m)({circumflex over (τ)}l−1(m))R(τ−{circumflex over (τ)}l−1(m)) and
to obtain the M groups of upsampling signals yL(m)(τ), where R(τ) represents a raised cosine wave function with a peak position set to zero, l takes a value of 2, 3, . . . , L, and m takes a value of 1, 2, 3, . . . , M; and the residual calculating unit 503 is configured to calculate a residual for each group of upsampling signals yL(m)(τ) according to a formula
The time delay power spectrum obtaining module 203 shown in
the incoherent accumulator 601 is configured to perform incoherent accumulation on zj(τ)=|yL(τ)|2 on J time slots according to z(τ)=Σj=0J−1zj(τ), where J is an integer greater than or equal to 1; and
the filter 602 is configured to perform filtering on z(τ) according to νn(τ)=ανn−1(τ)+(1−α)z(τ) to obtain a time delay power spectrum νn(τ), where α has a value range of 0<α<1.
It should be noted that, because content such as information exchange between the modules/units in the foregoing apparatus and execution processes is based on a same conception as the method embodiments of the present invention, and brings same technical effects as the method embodiments of the present invention, reference may be made to the description in the method embodiments of the present invention for details, which are not repeatedly described here.
Persons of ordinary skill in the art may understand that all or part of the steps in various methods in the foregoing embodiments may be accomplished by a program instructing relevant hardware, for example, one or multiple or all of the following methods:
performing processing on received signals ri(t) and pilot sequences pi(t) to obtain upsampling signals y(τ);
subtracting contribution values of a current multipath waveform to other multipath waveforms by using a raised cosine function with a peak position set to zero and according to a global maximum value of |y(τ)|2 and a total of M sampling points {circumflex over (τ)}l(m) near the global maximum value, to obtain an upsampling signal yL(τ) having a minimum residual in M groups of upsampling signals yL(m)(τ), where M is an integer greater than or equal to 1, and |y(τ)|2 is a squared absolute value of y(τ); and
obtaining a time delay power spectrum of the upsampling signal yL(τ) having the minimum residual.
Persons of ordinary skill in the art may understand that all or part of the steps of the methods in the foregoing embodiments may be accomplished by a program instructing relevant hardware. The program may be stored in a computer readable storage medium, and the storage medium may include: a read only memory (ROM, Read Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk, an optical disk, and so on.
The multipath searching method and the multipath searcher provided by the embodiments of the present invention are introduced in detail in the foregoing. In this specification, specific examples are used for illustrating principles and implementation manners of the present invention. The foregoing description of the embodiments are merely used to help understanding of the method of the present invention and its core idea. Meanwhile, persons of ordinary skill in the art may make modifications to the specific implementation manners and application scopes according to the idea of the present invention. In conclusion, the content of the specification should not be construed as a limitation to the present invention.
Number | Date | Country | Kind |
---|---|---|---|
201110422765.0 | Dec 2011 | CN | national |