Method and device for calculating weighted value of each path and RAKE receiver
Technical Field
The present invention relates to wireless communication signal receiving technology, and in particular, to a method and an apparatus for calculating weighted values of each path, and a RAKE receiver.
Background
In Code Division Multiplexed (CDMA) spread spectrum systems, the channel bandwidth is much larger than the flat fading bandwidth of the channel. Unlike conventional modulation techniques that require an equalization algorithm to cancel intersymbol interference between adjacent symbols, CDMA spreading codes require good autocorrelation properties in selecting the signal. Thus, the delay spread occurring in the radio channel can be regarded as simply a retransmission of the transmitted signal. Multipath signals may be considered to be uncorrelated if they are delayed relative to each other by more than one chip period.
Since the multipath signals contain available information, the receiver can improve the signal-to-noise ratio of the received signals by combining the multipath signals, and the receiver adopting such signal processing is called a RAKE receiver. Fig. 1 is a schematic structural diagram of a RAKE receiver, which specifically includes the following modules:
a multipath searching and distributing module 110, configured to receive a signal, complete multipath delay capturing and multipath delay distributing on the received signal, and output the processed delay information of each multipath to the descrambling and despreading module 120, the channel estimation module 130, and the weighted value calculation module 140;
a descrambling and despreading module 120, configured to complete descrambling and despreading operations on the chip-level received signal according to the input multipath delay information and scrambling code and spreading code information in the received signal, and output a descrambled and despread symbol-level received signal y to each path merging module 150;
the channel estimation module 130 estimates the channel fading condition of each time delay according to the input multipath time delay information, and outputs the channel fading information corresponding to each multipath time delay to the weighted value calculation module 140 of each path;
the weighted value calculating module 140 of each path calculates the system Gaussian noise power N according to the input multipath delay information, the channel fading information of each path and other external parameters0And receiving the chip level energy E of the useriCalculating a weighting factor ω of each path, and outputting the calculated weighting factor ω to each path merging module 150;
the path merging module 150 uses the formula Y ═ ω for the descrambled and despread path data Y according to the input path weighting factor ωHThe xy is combined, and the signal Y obtained after the combination is output to the demodulation decoding module 160; wherein]HRepresenting a conjugate transpose;
the demodulation decoding module 160 performs corresponding demodulation and decoding operations on the received signal Y, and outputs a signal after demodulation decoding.
In order to meet the increasing rate requirements of users, various technologies with higher and higher peak rates, such as High Speed Downlink Packet Access (HSDPA) and High Speed Uplink Packet Access (HSUPA), are introduced into the existing wireless communication system. The introduction of these new techniques presents significant challenges to conventional RAKE reception techniques. On one hand, the traditional RAKE receiver technology cannot meet the throughput requirement of high-speed service; on the other hand, due to the self-interference characteristic of the CDMA system, it is required that the interference generated by the high-rate service user to other users should be as small as possible, that is, the required signal-to-noise ratio (SNR) of the high-rate service user is reduced as much as possible on the premise of meeting the performance index of the high-rate service. Therefore, it is highly desirable to introduce advanced reception techniques with performance superior to conventional RAKE receivers.
Receivers with superior performance over conventional RAKE receivers include the following:
the enhanced RAKE receiver is improved on the basis of the traditional RAKE receiving technology, and the performance improvement is very limited;
the receiver can theoretically completely eliminate intersymbol interference (ISI) caused by multipath and multiple access interference (MUI) between different users so as to obtain good performance, but has great impact on the existing receiver structure due to great difference between the structure of the receiver and the traditional RAKE receiver;
the general RAKE (G-RAKE) receiver fundamentally considers the dominant factors limiting the performance of the traditional RAKE receiver, not only can achieve the same performance as the equalization receiver, but also has a similar implementation structure as the traditional RAKE receiver.
The structure of the G-RAKE receiver is similar to that of the conventional RAKE receiver shown in fig. 1, except that the specific calculation process of each of the path weight value calculation blocks 140 is different from that of the conventional RAKE receiver.
The calculation of each path weight value of the G-RAKE receiver requires calculation
<math> <mrow> <mi>ω</mi> <mo>=</mo> <msubsup> <mi>R</mi> <mi>u</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mi>h</mi> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
Wherein,
<math> <mrow> <mi>h</mi> <mo>=</mo> <msqrt> <msub> <mi>E</mi> <mn>0</mn> </msub> </msqrt> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
h is the quantity related to each path fading information, a unified Chinese name is not available in the field at present, and h is called each path fading function for convenience of description;
<math> <mrow> <msub> <mi>R</mi> <mi>u</mi> </msub> <mo>=</mo> <msub> <mi>E</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <mi>ISI</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>I</mi> </msub> <msub> <mi>R</mi> <mi>MUI</mi> </msub> <mo>+</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <msup> <mi>n</mi> <mo>′</mo> </msup> </msub> </mrow> </math>
<math> <mrow> <mo>=</mo> <mfrac> <msub> <mi>E</mi> <mn>0</mn> </msub> <msup> <mi>N</mi> <mn>2</mn> </msup> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <munder> <mrow> <mi>l</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mi>l</mi> <mo>≠</mo> <mn>0</mn> </mrow> </munder> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mi>N</mi> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mo>|</mo> <mi>m</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>×</mo> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <mi>iT</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mtext>-iT-</mtext> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mo>+</mo> <mfrac> <msub> <mi>E</mi> <mn>1</mn> </msub> <msup> <mi>N</mi> <mn>2</mn> </msup> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mi>N</mi> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mo>|</mo> <mi>m</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>×</mo> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <mi>iT</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <mi>iT</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>[</mo> <mrow> <mn>1</mn> <mo>-</mo> <mi>δ</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mi>δ</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> </mrow> <mo>]</mo> <mo></mo> </mrow> <mo></mo> </mrow> </math>
in the formula (2) or (3), E0Symbol level energy for the target user, EISymbol level energy, N, for interfering users0The single-sided power spectral density of the system white gaussian noise, the symbol-level energy and the single-sided power spectral density of the system white gaussian noise can be estimated according to a related algorithm, and in a receiver system, the quantities are used as known parameters; n is the spreading factor of the target user, RISIAs a function of intersymbol interference, RMUIAs a function of multiple access interference, TcFor the duration of one spreading chip, T ═ NTc,gl、gqChannel fading factors, τ, for the l-th and q-th paths, respectivelyl、τqIs the multipath propagation delay of the L-th path and the q-th path, L is the total number of effective paths, RpObtaining an autocorrelation function of a transmit pulse shaping filter in advance according to a computational expression of a shaping filter used by a system; d. d1And d2To solve the multipath delay corresponding to the diameter expansion.
According to the above calculation method, the calculation of each path weight value of the G-RAKE receiver as shown in fig. 2 includes the following steps:
step 201: according to the input path fading factor glMultipath position τlAnd symbol level energy E of target user0H is calculated by using the formula (2);
step 202: finding RISI: as can be seen from equation (3):
<math> <mrow> <msub> <mi>R</mi> <mi>ISI</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msup> <mi>N</mi> <mn>2</mn> </msup> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mmultiscripts> <mrow> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <munder> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mi>l</mi> <mo>≠</mo> <mn>0</mn> </mrow> </munder> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <mtext></mtext> </mrow> </mmultiscripts> <msub> <mi>g</mi> <mi>l</mi> </msub> <mrow> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mi>N</mi> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mo>|</mo> <mi>m</mi> <mo>|</mo> <mo>)</mo> </mrow> <mtext>×</mtext> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mtext>-iT-</mtext> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mtext>×</mtext> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <mi>iT</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </mrow> </math>
the step is based on the input fading factors g of each pathlMultipath position taulMultipath time delays d, d corresponding to the solution of the hole enlargement1And d2Andequation (4) calculates the intersymbol interference function RISI;
Step 203: finding RMUI: as can be seen from equation (3):
<math> <mrow> <msub> <mi>R</mi> <mi>MUI</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msup> <mi>M</mi> <mn>2</mn> </msup> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mi>N</mi> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mo>|</mo> <mi>m</mi> <mo>|</mo> <mo>)</mo> </mrow> <mtext>×</mtext> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <mi>iT</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <mi>iT</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mrow> <mrow> <mo>[</mo> <mn>1</mn> <mo>-</mo> <mi>δ</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mi>δ</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>]</mo> </mrow> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo></mo> </mrow> </math>
the step is based on the input fading factors g of each pathlMultipath position taulMultipath time delay d, d corresponding to diameter expanding1And d2And equation (4) calculating a multiple access interference function RMUI;
Step 204: finding Rn′: as can be seen from equation (3):
Rn′=Rp(d1-d2) (6)
this step is based on the autocorrelation function R of the shaping filterpInput multipath time delay d1And d2Calculating the interference function R of the background thermal noise by using the formula (6)n′;
Step 205: according to the input spreading factor N of the target user and the symbol level energy E of the target user0Symbol level energy E of interfering usersISingle-sided power spectral density N of system Gaussian white noise0And the calculation results of step 202 to step 204, calculating to obtain a noise correlation matrix Ru;
Step 206: according to the noise correlation matrix R obtained in step 205uAnd h obtained in step 201, to obtain the path combining weighting factors of the G-RAKE receiver.
However, compared with the conventional RAKE receiver, the processing method for calculating the weighted values of each path adopted by the G-RAKE receiver requires a large amount of calculation overhead, and thus the implementation complexity and cost are high.
Disclosure of Invention
In view of the above, the present invention provides a method for calculating a weighted value of each path, which can effectively reduce the computation of the combining weighted value of each path. The method comprises the following steps:
A. calculating each path fading function h according to the input each path fading factor, the multipath position information and the symbol level energy of the target user;
B. according to the input fading factors of all paths, the multi-path position information and the multi-path time delay corresponding to the de-expanding, the multi-address interference function R is calculatedMUI;
C. According to the autocorrelation function R of the shaping filterpInput multipath time delay, and calculating the interference function R of the background thermal noisen′;
D. Based on the sum E of the inputted chip-level energies of all users independent of the spreading factorc_TSingle-sided power spectral density N of system Gaussian white noise0And the calculation results of the step B and the step C are calculated to obtain a noise correlation matrix Ru;
E. According to the fading functions h of all paths obtained in the step A and the noise correlation matrix R obtained in the step DuAnd calculating to obtain the weighted value omega of each diameter.
Step A, calculating each path fading function h as:
according to the formula <math> <mrow> <mi>h</mi> <mo>=</mo> <msqrt> <msub> <mi>E</mi> <mn>0</mn> </msub> </msqrt> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> </mrow> </math> Calculating each path fading function h, wherein E0Symbol level energy, g, for the target userlIs the channel fading factor, R, of the l pathpIs the autocorrelation function of the transmit pulse shaping filter, d is the multipath time delay corresponding to the de-expansion, τlIs the multipath propagation delay of the l-th path.
Step B calculating multiple access interference function RMUIComprises the following steps:
according to the formula
<math> <mrow> <msub> <mi>R</mi> <mi>MUI</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>i</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>i</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> </mrow> </math> Calculating a multiple access interference function RMUI(ii) a Where N is the spreading factor of the target user, l and q are the labels of the paths, glAnd gqFor each path fading factor, τlAnd τqIs the multipath location information, L is the total number of effective paths, d1And d2To solve the multipath delay corresponding to the diameter expansion.
C, calculating the interference function R of the background thermal noisen′Comprises the following steps: according to the formula Rn′=Rp(d1-d2) Calculating the interference function R of the background thermal noisen′Wherein R ispFor the autocorrelation function of the transmit pulse shaping filter, d1And d2To solve the multipath delay corresponding to the diameter expansion.
Step D calculating the noise correlation matrix RuComprises the following steps:
according to the formula
<math> <mrow> <msub> <mi>R</mi> <mi>u</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>I</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>MUI</mi> </msub> <mo>+</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <msup> <mi>n</mi> <mo>′</mo> </msup> </msub> </mrow> </math>
<math> <mrow> <mo>=</mo> <msub> <mi>E</mi> <mrow> <mi>c</mi> <mo>_</mo> <mi>T</mi> </mrow> </msub> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <mfrac> <mrow> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mi>N</mi> </mfrac> <mo>-</mo> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </math>
Computing a noise correlation matrix Ru(ii) a Where N is the spreading factor of the target user, l and q are the labels of the paths, glAnd gqFor each path fading factor, τlAnd τqIs the multipath location information, L is the total number of effective paths, d1And d2To solve the multipath delay corresponding to the expansion, Ec_TIs the sum of all user chip-level energies, R (n), independent of the spreading factor1,n2) Is n1And n2Of (2), wherein <math> <mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mtext>-[</mtext> <mfrac> <mrow> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> </mrow> <msub> <mi>N</mi> <mi>OVS</mi> </msub> </mfrac> <mtext>]×</mtext> <msub> <mi>N</mi> <mi>OVS</mi> </msub> <mo>,</mo> </mrow> </math> <math> <mrow> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>-</mo> <mo>[</mo> <mfrac> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> </mrow> <msub> <mi>N</mi> <mi>OVS</mi> </msub> </mfrac> <mo>]</mo> <mo>×</mo> <msub> <mi>N</mi> <mi>OVS</mi> </msub> <mo>,</mo> </mrow> </math> NOVSIs the sampling rate.
And E, calculating to obtain the weighted value omega of each diameter as follows: according to the formula <math> <mrow> <mi>ω</mi> <mo>=</mo> <msubsup> <mi>R</mi> <mi>u</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mi>h</mi> </mrow> </math> The weighted value ω of each path is obtained.
Another objective of the present invention is to provide an apparatus for calculating weighted values of each path in a RAKE receiver system, which specifically includes:
each path fading function calculating unit is used for receiving the target user symbol level energy from a descrambling and despreading module of the RAKE receiver, calculating each path fading function h by the multipath fading factor and the multipath position information of the channel estimation module, and inputting each path fading function h obtained by calculation into the weighted value synthesizing unit;
a multiple access interference function calculation unit for receiving multipath fading factor and multipath position information from the channel estimation module of the RAKE receiverAnd calculating a multiple access interference function R from the multipath search assignment module multipath time delaysMUIAnd calculating the multiple access interference function RMUIAn input noise correlation matrix calculation unit;
a thermal noise interference function calculation unit for receiving multipath time delay of multipath search distribution module of RAKE receiver, and autocorrelation function of shaping filter, calculating interference function R of background thermal noisen′And calculating the interference function R of the background thermal noisen′An input noise correlation matrix calculation unit;
a noise correlation matrix calculation unit for calculating a multiple access interference function R from the multiple access interference functionMUIInterference function R of background thermal noise of thermal noise interference function calculation unitn′Sum of all user chip-level energies of the receiver system Ec_TAnd single-sided power spectral density N of system Gaussian white noise0And calculating to obtain a noise correlation matrix RuAnd correlating the resulting noise correlation matrix RuAn input weight value synthesizing unit;
a weight value synthesis unit for synthesizing the weighted values according to the received path fading functions h and the noise correlation matrix RuThe weighted values ω of the respective paths are synthesized and the synthesized weighted values ω of the respective paths are output to the respective path merging modules.
The noise correlation matrix calculation unit calculates to obtain a noise correlation matrix RuThe method comprises the following steps: let the intersymbol interference function be a multiple access interference function RISIEqual to the multiple access interference function RMUIAnd according to an intersymbol interference function RISIMultiple access interference function RMUIBackground thermal noise interference function Rn′Sum of all user chip-level energies of the receiver system Ec_TAnd single-sided power spectral density N of system Gaussian white noise0And calculating to obtain a noise correlation matrix Ru。
The RAKE receiver is a generalized RAKE receiver.
Another objective of the present invention is to provide a RAKE receiver for calculating each path weight value by using the above method, which includes a multipath searching and distributing module, a descrambling and despreading module, a channel estimation module, a path weight value calculating module, a path combining module and a demodulation and decoding module, wherein the path weight value calculating module includes:
each path fading function calculating unit is used for receiving the target user symbol level energy from the descrambling and despreading module, the multipath fading factor and the multipath position information of the channel estimation module, calculating each path fading function h, and inputting each path fading function h obtained by calculation into the weighted value synthesizing unit;
a multiple access interference function calculation unit for receiving the multipath fading factor and multipath position information from the channel estimation module, and the multipath time delay from the multipath search distribution module and calculating the multiple access interference function RMUIAnd calculating the multiple access interference function RMUIAn input noise correlation matrix calculation unit;
a thermal noise interference function calculation unit for receiving the multipath time delay of the multipath search distribution module, and the autocorrelation function of the shaping filter, and calculating the interference function R of the bottom thermal noisen′And calculating the interference function R of the background thermal noisen′An input noise correlation matrix calculation unit;
a noise correlation matrix calculation unit for calculating a multiple access interference function R from the multiple access interference functionMUIInterference function R of background thermal noise of thermal noise interference function calculation unitn′Sum of all user chip-level energies of the receiver system Ec_TAnd single-sided power spectral density N of system Gaussian white noise0And calculating to obtain a noise correlation matrix RuAnd correlating the resulting noise correlation matrix RuAn input weight value synthesizing unit;
a weight value synthesis unit for synthesizing the weighted values according to the received path fading functions h and the noise correlation matrix RuSynthesizing the weighted value omega of each path and outputting to the merging modules of each pathThe combined radial weights ω.
The noise correlation matrix calculation unit calculates to obtain a noise correlation matrix RuThe method comprises the following steps: let the intersymbol interference function be a multiple access interference function RISIEqual to the multiple access interference function RMUIAnd according to an intersymbol interference function RISIMultiple access interference function RMUIBackground thermal noise interference function Rn′Sum of all user chip-level energies of the receiver system Ec_TAnd single-sided power spectral density N of system Gaussian white noise0And calculating to obtain a noise correlation matrix Ru。
The RAKE receiver is a generalized RAKE receiver.
It can be seen from the above technical solutions that, in the process of calculating the weighted value of each path, the method of the present invention does not need to calculate the inter-symbol interference function, so that compared with the G-RAKE receiver, the method of the present invention can effectively reduce the computation of the combining weights of each path, and the quality of the output signal is hardly affected; moreover, the calculation of each radial weighting value is independent of the spreading factor, and users with different spreading factors can share one RAKE receiver without increasing the number of RAKE receivers for users with different spreading factors. By adopting the scheme of the invention, the calculation overhead and the realization complexity can be greatly reduced and the cost can be saved on the premise of outputting the signals with the same quality as the G-RAKE receiver.
Drawings
FIG. 1 is a schematic diagram of a RAKE receiver or G-RAKE receiver architecture;
FIG. 2 is a flow diagram of a prior art implementation of the path weight calculation modules of a G-RAKE receiver;
FIG. 3 is a flow chart of the implementation of the path weight calculation modules of the RAKE receiver according to the present invention;
fig. 4 is a block diagram of the calculation module of each radial weight value of the RAKE receiver according to the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings.
According to the simplified formula (7) proposed by the prior art for the calculation method of the respective path weight values of the G-RAKE receiver:
<math> <mrow> <mfrac> <mn>1</mn> <msup> <mi>N</mi> <mn>2</mn> </msup> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mi>N</mi> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mo>|</mo> <mi>m</mi> <mo>|</mo> <mo>)</mo> </mrow> <mo>×</mo> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <mi>iT</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <mi>iT</mi> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>i</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>i</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
equation (3) can be reduced to:
<math> <mrow> <msub> <mi>R</mi> <mi>u</mi> </msub> <mo>=</mo> <msub> <mi>E</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <mi>ISI</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>I</mi> </msub> <msub> <mi>R</mi> <mi>MUI</mi> </msub> <mo>+</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <msup> <mi>n</mi> <mo>′</mo> </msup> </msub> </mrow> </math>
<math> <mrow> <mo>=</mo> <msub> <mi>E</mi> <mn>0</mn> </msub> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>1</mn> <msup> <mi>N</mi> <mn>2</mn> </msup> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> <mo>-</mo> <mi>N</mi> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>m</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> </mrow> </math>
<math> <mrow> <mo>+</mo> <msub> <mi>E</mi> <mn>1</mn> </msub> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein R (n)1,n2) Is n1And n2The correlation function of (a) is determined,
<math> <mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>-</mo> <mo>[</mo> <mfrac> <mrow> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> </mrow> <msub> <mi>N</mi> <mi>OVS</mi> </msub> </mfrac> <mo>]</mo> <mo>×</mo> <msub> <mi>N</mi> <mi>OVS</mi> </msub> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>-</mo> <mo>[</mo> <mfrac> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> </mrow> <msub> <mi>N</mi> <mi>OVS</mi> </msub> </mfrac> <mo>]</mo> <mo>×</mo> <msub> <mi>N</mi> <mi>OVS</mi> </msub> <mo>.</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
NOVSis the sampling rate of the G-RAKE receiver module.
From equation (3) one can see the matrix RuIn the constituent elements of (1), RISIAnd RMUIThe expression of (c) differs only in the summation of i: the former reduces the summation operation when i is 0 than the latter. Is due to RISIAnd RMUIThis difference in (c) results in a simplified pair R in equation (8)ISIThe calculation of (2) involves a summation operation of the products of the 2N-1 terms. If R is considered in the implementation of G-RAKEISI=RMUIThen equation (8) can be further simplified to:
<math> <mrow> <msub> <mi>R</mi> <mi>u</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>I</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>MUI</mi> </msub> <mo>+</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <msup> <mi>n</mi> <mo>′</mo> </msup> </msub> </mrow> </math>
<math> <mrow> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mn>0</mn> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>I</mi> </msub> <mo>)</mo> </mrow> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mo>=</mo> <msub> <mi>E</mi> <mi>T</mi> </msub> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow> </math>
wherein ET=E0+EIAnd is the symbol level energy. Table 1 shows the computation amount of each time slot of the G-RAKE receiver, where M is the number of iterations for obtaining the weighted value ω, and is generally 1-3; j is the total number of fingers of the G-RAKE receiver,the value is generally J ≈ 2L.
TABLE 1
As can be seen from the data in Table 1, R is determinedISIAnd RMUIOccupies a considerable weight in the total computation of G-RAKE, and therefore R is considered in the G-RAKE implementation processISI=RMUIThe operation amount of the G-RAKE receiver is effectively reduced, and the implementation complexity is reduced.
According to the formula (8), since RISI=RMUICalculating RuTime and energy E0、EIThe matrices multiplied are equal, so that only the total symbol level energy E needs to be knownTWithout distinguishing E therefrom0、EIFurther, the number of parameters required by the G-RAKE receiver can be reduced, and the implementation complexity of the G-RAKE receiver is further reduced.
Further, as can be seen from the formulas (4) and (5), RMUISpreading factor ofCan be extracted, so equation (8) can be rewritten as:
<math> <mrow> <msub> <mi>R</mi> <mi>u</mi> </msub> <mo>=</mo> <msub> <mi>E</mi> <mi>T</mi> </msub> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> </math>
<math> <mrow> <mo>=</mo> <msub> <mi>E</mi> <mrow> <mi>c</mi> <mo>_</mo> <mi>T</mi> </mrow> </msub> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <mfrac> <mrow> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mi>N</mi> </mfrac> <mo>-</mo> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <msub> <mi>N</mi> <mn>0</mn> </msub> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow> </math>
namely RMUIFrom E energy weighting factorTBecome into I.e., the sum of all user chip-level energies independent of the spreading factor, so RuIs independent of the spreading factor. Therefore if R is considered in the G-RAKE implementationISI=RMUIAll users can share the calculation result of one G-RAKE receiver as long as the transmission environment of the users carrying out G-RAKE reception is the same, namely multipath delay and channel fading are the same.
Based on the above analysis, the G-RAKE receiver according to the embodiment of the present invention adopts the structure shown in fig. 1, that is, the structure of the existing G-RAKE receiver or the existing RAKE receiver, except that the implementation process of each radial weight value calculation module in the G-RAKE receiver or the existing RAKE receiver adopts the scheme of the present invention, as specifically shown in fig. 3, the G-RAKE receiver includes the following steps:
step 301: according to the input path fading factor glMultipath position τlAnd symbol level energy E of target user0Calculating each path fading function h by using a formula (2); the step is the same as step 201 in the implementation process of the weight calculation module of the existing G-RAKE receiver;
step 302: finding RMUI: from the formula (12), it can be seen that
<math> <mrow> <msub> <mi>R</mi> <mi>MUI</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>l</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <munderover> <mi>Σ</mi> <mrow> <mi>q</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>L</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>g</mi> <mi>l</mi> </msub> <msubsup> <mi>g</mi> <mi>q</mi> <mo>*</mo> </msubsup> <mo>[</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>i</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>i</mi> <msub> <mi>T</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>R</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>l</mi> </msub> <mo>)</mo> </mrow> <mo>×</mo> <msubsup> <mi>R</mi> <mi>p</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>τ</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> <mo>]</mo> <mo>,</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow> </math>
The step is based on the input fading factors g of each pathlMultipath position taulMultipath time delay d, d corresponding to diameter expanding1And d2And equation (13) calculating a multiple access interference function RMUI;
Step 303: finding Rn′: from the formula (12), it can be seen that
Rn′=Rp(d1-d2) (14)
This step is based on the autocorrelation function R of the shaping filterpInput multipath time delay d1And d2Calculating an interference function R of the background thermal noise according to the formula (14)n′(ii) a Equation (14) is the same as equation (6), so this step is the same as step 204 in the implementation process of the weight calculation module of the existing G-RAKE receiver;
step 304: based on the sum E of the inputted chip-level energies of all users independent of the spreading factorc_TSingle-sided power spectral density N of system Gaussian white noise0And the calculation results of the step 302 and the step 303 are calculated according to the formula (12) to obtain the noise correlation matrix Ru;
Step 305: according to the path fading functions h obtained in step 301 and the noise correlation matrix R obtained in step 304uThe weighted value ω of each path of the G-RAKE receiver is calculated by equation (1).
To implement the above method, the structure and connection relationship of each part of the weighted value calculation module 140 of the RAKE receiver according to the embodiment of the present invention is shown in fig. 4, which specifically includes the following steps:
each path fading function calculating unit 141 for receiving the target user symbol level energy E from the descrambling and despreading module 1200And a multipath fading factor g of the channel estimation module 130lAnd multipath position information taulCalculating each path fading function h according to formula (2) using the received parameters, and inputting each calculated path fading function h to the weighted value synthesizing unit 145;
a multiple access interference function calculation unit 142 for receiving the multipath fading factor g from the channel estimation module 130lAnd multipath position information taulAnd multipath delays d, d of the multipath search assignment module 1101And d2Calculating a multiple access interference function R according to equation (13) using the received parametersMUIAnd calculating the multiple access interference function RMUIAn input noise correlation matrix calculation unit 144;
a thermal noise interference function calculating unit 143 for receiving the multipath time delay d of the multipath search allocating module 1101And d2And the autocorrelation function R of the shaping filterpCalculating the interference function R of the noise floor according to equation (14) using the received parametersn′And calculating the interference function R of the background thermal noisen′An input noise correlation matrix calculation unit 144;
a noise correlation matrix calculation unit 144 for calculating a multiple access interference function R according to the multiple access interference function from the multiple access interference function calculation unit 142MUIInterference function R of background thermal noise of thermal noise interference function calculation unit 143n′Sum of all user chip-level energies of the receiver system Ec_TAnd single-sided power spectral density N of system Gaussian white noise0Calculating to obtain a noise correlation matrix R according to the formula (12)uAnd correlating the resulting noise correlation matrix RuAn input weight value synthesizing unit 145;
a weighted value synthesizing unit 145 for synthesizing the weighted values according to the received path fading functions h and the noise correlation matrix RuThe weighted values ω are synthesized by the formula (1), and the synthesized weighted values ω are output to the path merging modules 150.
Figure 4 shows an embodiment of the invention in which the path weight value calculation blocks have a reduced number of intersymbol interference function calculation units and associated input-output connections compared to prior art G-RAKE receiver path weight value calculation blocks; in the processing method, the expense for calculating the intersymbol interference function is correspondingly saved.
The path weight value calculation module shown in fig. 4 is used as the path weight value calculation module 140 of the RAKE receiver shown in fig. 1, and the connection relationship between the structure and each part of the RAKE receiver shown in fig. 1 is the connection relationship between the structure and each part of the RAKE receiver of the present invention.
Simulation analysis of the G-RAKE receiver proposed by the scheme of the invention shows that R is considered in the implementation process of the G-RAKE receiverISI=RMUIHardly affecting the output result of the G-RAKE receiver. Therefore, R is considered in the implementation of the G-RAKE receiverISI=RMUIOn the premise of not influencing the quality of output signals, the operation amount of the G-RAKE receiver can be effectively reduced, the number of parameters required by the G-RAKE receiver is reduced, and the implementation complexity of the G-RAKE receiver is further reduced; and makes the calculation of the G-RAKE weighted value independent of the spreading factor, users with different spreading factors can share one G-RAKE receiver without adding the number of G-RAKE receivers for users with different spreading factors.
Those skilled in the art should recognize that, since the structural difference between the G-RAKE receiver and the RAKE receiver is mainly due to the difference between the respective radial weight value calculation modules, the method for calculating the respective radial weight values and the apparatus for implementing the calculation method proposed by the present invention can be applied not only to the G-RAKE receiver but also to the RAKE receiver having the same or similar structure.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.