CN100539456C - Multi-user combined detection method - Google Patents

Abstract

The invention discloses a kind of multi-user combined detection method, comprising: the intermediate sequence according to received signal carries out the estimation of many code channels signal, obtains the aggregate channel response of multiple user signals; According to the aggregate channel response matrix A of described aggregate channel response structure based on the empty code channel number of maximum _Max, described A _MaxBe normal dimension matrix; Utilize described aggregate channel response matrix A based on the empty code channel number of maximum _MaxCarry out joint-detection, obtain the user data symbol of output; According to the empty code channel number of maximum the user data symbol of output is carried out demodulation, obtain joint detection results based on the empty code channel number of maximum; From based on the actual activated code channel data that extract the joint detection results of the empty code channel number of maximum in each maximum empty code channel data, obtain the multi-user association testing result of expectation.Utilize the present invention, can reduce the complexity that processor and application-specific integrated circuit ASIC are realized, improve the joint-detection quality.

Description

Multi-user combined detection method

Technical field

The present invention relates to the mobile communication technology field, be specifically related to a kind of multi-user combined detection method.

Background technology

3-G (Generation Three mobile communication system) adopts CDMA (Wideband Code Division Multiple Access (WCDMA)) technology, performances such as it has that anti-interference, capacity is big, good confidentiality and soft handover, but because the general strict orthogonal not of spreading code, its non-zero cross-correlation coefficient can cause the phase mutual interference between each user, i.e. MAI (multiple access interference).The existence of MAI makes power system capacity be restricted, and brings " near-far interference ", has a strong impact on systematic function.For example, TD-SCDMA (Time Division-Synchronous Code Division Multiple Access) system synthesis has used time division multiple access TDMA, frequency division multiple access FDMA and three kinds of multiplex techniques of code division multiple access CDMA to improve power system capacity.Yet, because the complexity of wireless transmission environment, each user who is in the same time slot of same frequency but adopts different OVSF (orthogonal variable step-length) sign indicating number modulation are when arriving receiving terminal, the orthogonality of its signal is damaged, and because multidiameter delay and other users' interference has caused intersymbol interference ISI and multiple access to disturb MAI.

Disturb and intersymbol interference in order to eliminate multiple access, in the TD-SCDMA system, base station and terminal have all adopted associated detection technique.Associated detection technique is on the basis of traditional detection technology, make full use of the prior information that causes all subscriber signals that MAI disturbs and multipath thereof (as the subscriber channel sign indicating number that knows and training sequence, each user's channel estimating etc.), the separation of subscriber signal being used as the joint detection process that is mutually related of a unification finishes, thereby has good interference free performance, reduced the requirement of system to the power control precision, therefore can effectively utilize the link frequency spectrum resource more, improve power system capacity significantly.

The specific implementation method of joint detection algorithm has multiple, roughly is divided into three major types such as nonlinear algorithm, linear algorithm and decision feedback algorithms.Wherein, linear combined detection algorithm mainly contains two kinds of ZF-BLE (ZF data block linear equalization) and MMSE-BLE (Minimum Mean Square Error data block linear equalization); Decision feedback algorithms mainly contains two kinds of ZF-BDFE (ZF data block decision feedback equalization) and MMSE-BDFE (Minimum Mean Square Error data block decision feedback equalization).

Wherein, the Mathematical Modeling formula of ZF-BLE is:

${\hat{d}}_{c,\mathrm{ZF}-\mathrm{BLE}}={\left(A^H{R}_n^{-1}A\right)}^{-1}{A}^H{R}_n^{-1}e,---\left(1\right)$

The Mathematical Modeling formula of MMSE-BLE is:

${\hat{d}}_{c,\mathrm{MMSE}-\mathrm{BLE}}={(A^H{R}_n^{-1}A+R_d^{-1})}^{-1}{A}^H{R}_n^{-1}e,---\left(2\right)$

In above formula, A is a channel combination response matrix,

Be the interchannel noise correlation inverse matrix, e is the reception data,

Be the inverse matrix of intersymbol noise correlation matrix,

Be the output symbol of ZF-BLE,

Output symbol for MMSE-BLE.

The Mathematical Modeling of ZF-BDFE is that ZF-BLE adds decision-feedback; The Mathematical Modeling of MMSE-BDFE is that MMSE-BLE adds decision-feedback.

At present, the implementation of TD-SCDMA base station physical layer segment multi-user association detection algorithm mainly is based on variable actual activated code channel number, promptly carries out different processing according to different actual activated code channel numbers.

With the example that is embodied as of ZF-BLE, present implementation promptly is to find the solution following linear equation:

$A^H{R}_n^{-1}A{\hat{d}}_{c,\mathrm{ZF}-\mathrm{BLE}}=A^H{R}_n^{-1}e---\left(3\right)$

Wherein, A is the aggregate channel response matrix, is (L _Data+ W _FL-1) * K _VcThe change of N dimension matrix, its structure be as shown in Figure 1:

K among Fig. 1 _VcFor reality activates empty code channel number (parameter), W _FLBe channel estimation window long (constant) that N is a symbolic number, Q _MaxBe maximum spreading gain, L _Data=N * Q _Max

So-called empty code channel is the combination of the channel code of different spreading factors being regarded as the virtual channel sign indicating number of several maximum spreading factors in implementation algorithm, thereby when algorithm is realized, a code channel is decomposed into a plurality of empty code channels, the virtual channel sign indicating number by-1,0,1 three kind of combinations of values form.As, the maximum spreading factor of the channel code that uses in system is 16, then spreading factor is that 8 channel code can regard that 2 spreading factors are the combination of 16 virtual channel sign indicating number as, thus, 1 spreading factor be 8 code channel to be decomposed into 2 spreading factors be 16 empty code channel, if reality has been used 3 channel codes, 1 spreading factor is 16,1 spreading factor is that 8,1 spreading factors are 4, then K _Vc=7.Actual activated code channel exists corresponding relation with the empty code channel number of actual activation.

Being the inverse matrix of interchannel noise correlation matrix, is (L _Data+ W _FL-1) * (L _Data+ W _FL-1) normal dimension symmetric positive definite matrix; E is (L for receiving data _Data+ W _FL-1) normal dimensional vector * 1.

Find the solution the joint detection results that above-mentioned linear equation (3) obtains the multi-user Its structure as shown in Figure 2.This result is carried out demodulation, can obtain each user's received signal.

At present, this implementation based on variable actual activated code channel number be applied to current two kinds of different frameworks processor (based on pipelined architecture with based on the parallel computation framework) or ASIC (application-specific integrated circuit (ASIC)) different shortcomings is arranged.

For DSP based on pipelined architecture, usually compiler is that the optimization of the code of constant is better than the code that cycle-index is a variable for cycle-index, and because actual activated code channel number is variable, cycle-index is a variable, therefore it is relatively poor to optimize effect, brings the increase of size of code and the problems such as increase of time expense thus.Simultaneously, because cycle-index does not fixedly cause the unfixing of united detection processing time delay, therefore can bring delay variation.

For processor and ASIC circuit based on the parallel computation framework, because actual activated code channel number is variable, but causing the parallel processing number is parameter rather than fixed constant, and expects that based on the processor and the ASIC circuit of parallel computation framework the parallel processing number is a fixed constant.Processing demands variable and line number not only can be brought problems such as design difficulty is big, algorithm is realized complexity, size of code is big, and complicated implementation also can be brought problems such as time overhead is big, delay variation is big.

Summary of the invention

The purpose of this invention is to provide a kind of multi-user combined detection method, be unfavorable for that based on the detection mode of variable actual activated code channel number processor and ASIC realize and optimize and make restituted signal produce the shortcoming of delay variation to overcome prior art, reduce the complexity that processor and ASIC realize that multi-user association detects, improve signal detection quality.

For this reason, the invention provides following technical scheme:

A kind of multi-user combined detection method comprises step:

A, carry out many code channels signal according to the intermediate sequence of received signal and estimate, obtain the aggregate channel response of multiple user signals;

B, according to the aggregate channel response matrix A of described aggregate channel response structure based on the empty code channel number of maximum _Max, described A _MaxBe normal dimension matrix;

C, utilize described aggregate channel response matrix A based on the empty code channel number of maximum _MaxCarry out joint-detection, obtain the user data symbol of output;

D, according to the empty code channel number of maximum to output user data symbol carry out demodulation, obtain joint detection results based on the empty code channel number of maximum;

E, from based on the actual activated code channel data that extract the joint detection results of the empty code channel number of maximum in each maximum empty code channel data, obtain the multi-user association testing result of expectation.

Described aggregate channel response matrix A based on the empty code channel number of maximum _MaxBe (L _Data+ W _FL-1) * K _MaxThe normal dimension matrix of N; Wherein, K _MaxBe the empty code channel number of maximum, W _FLLong for channel estimation window, N is a symbolic number, L _Data=N * Q _Max, Q _MaxBe maximum spreading gain.

Described step B comprises:

Obtain each submatrix in the aggregate channel response matrix;

Valid data in each submatrix are rearranged, be specially: the preceding K that the data of activated code channel aggregate channel impulse response correspondence in each submatrix is placed on this submatrix _VcOn the individual position, all the other positions fill in 0, described K _VcFor reality activates empty code channel number;

Form aggregate channel response matrix A by each submatrix after rearranging based on the empty code channel number of maximum _Max

Described step C is specially:

C1, employing ZF data block linear equalizing method are carried out joint-detection to multiple user signals; Perhaps

C2, employing Minimum Mean Square Error data block linear equalizing method are carried out joint-detection to multiple user signals; Perhaps

C3, employing ZF data block decision feedback equalization method are carried out joint-detection to multiple user signals; Perhaps

C4, employing Minimum Mean Square Error data block decision feedback equalization method are carried out joint-detection to multiple user signals.

Described step C1 is specially:

Find the solution linear equation: ${A_{\max }}^H{R}_n^{-1}{A}_{\max }{\hat{d}\′}_{c,\mathrm{ZF}-\mathrm{BLE}}={A_{\max }}^H{R}_n^{-1}e,$ The user data symbol that obtains exporting

Wherein,

Expression interchannel noise correlation inverse matrix is (L _Data+ W _FL-1) * (L _Data+ W _FL-1) normal dimension symmetric positive definite matrix, e is (L for receiving data _Data+ W _FL-1) normal dimensional vector * 1,

K _MaxThe normal dimensional vector of N * 1.

Described step C2 is specially:

Find the solution linear equation: $({A_{\max }}^H{R}_n^{-1}{A}_{\max }+R_d^{-I}){\hat{d}\′}_{c,\mathrm{MMSE}-\mathrm{BLE}}={A_{\max }}^H{R}_n^{-1}e,$ The user data symbol that obtains exporting

Wherein,

Expression interchannel noise correlation inverse matrix,

The inverse matrix of expression intersymbol noise correlation matrix is K _MaxN * K _MaxThe normal dimension symmetric positive definite matrix of N, e is (L for receiving data _Data+ W _FL-1) normal dimensional vector * 1.

By above technical scheme provided by the invention as can be seen, the present invention utilize multi-user association detect in maximum code channel number be the characteristics of fixed numbers and multi-user association detection algorithm itself, in the implementation procedure that multi-user association detects based on the maximum code channel number of fixing rather than based on variable activated code channel number.Aggregate channel by multiple user signals responds the aggregate channel response matrix of structure based on the empty code channel number of maximum, and utilize it to carry out joint-detection, from joint detection results, extract the incoming symbol of the valid data of actual activated code channel correspondence then as decoding based on the empty code channel number of maximum, realize simple, reliable.

Apply the present invention to processor and ASIC circuit based on the parallel computation framework, because the aggregate channel response matrix based on the empty code channel number of maximum is normal dimension matrix, therefore only need to design at finding the solution normal dimensional linear equation, can reduce design difficulty greatly, reduce the complexity that realizes, reduce the code expense, and can maximally utilise computation capability, thereby reduce time overhead based on the processor or the ASIC circuit of parallel computation framework at the design of finding the solution normal dimensional linear equation.In addition, the processing time of basic fixed also helps reducing delay variation, and simple relatively design helps increasing the stability of program.

Equally, apply the present invention to processor, mean cycle-index for fixing it owing to find the solution normal dimensional linear equation, and cycle-index is a fixed value based on pipelined architecture, better program optimization effect must be brought, therefore size of code and time expense can be reduced.Therefore the fixing stability that has guaranteed the united detection processing time delay substantially of cycle-index reduced delay variation, improved the detection quality.

Description of drawings

Fig. 1 is the structural representation of aggregate channel response matrix A and submatrix thereof in the ZF-BLE Mathematical Modeling;

Fig. 2 is a joint detection results of utilizing the ZF-BLE method to obtain in the prior art;

Fig. 3 is the baseband model of TD-SCDMA system;

Fig. 4 is the burst structure schematic diagram of TD-SCDMA system;

Fig. 5 is the realization flow figure of the inventive method;

Fig. 6 is the aggregate channel response matrix A in the inventive method _MaxAnd the structural representation of submatrix;

Fig. 7 is the correlation matrix in the inventive method

And the structural representation of submatrix;

Fig. 8 is the correlation matrix in the inventive method

And the structural representation of submatrix;

The 9th, the correlation matrix in the inventive method

And the structural representation of submatrix;

Figure 10 is based on the user data symbol of maximum empty code channel number output

Structural representation.

Embodiment

Core of the present invention be to utilize multi-user association detect in maximum code channel number be the characteristics of fixed numbers and multi-user association detection algorithm itself, in the implementation procedure that multi-user association detects based on the maximum code channel number of fixing rather than based on variable activated code channel number.Intermediate sequence according to received signal carries out the estimation of many code channels signal, utilize the aggregate channel response matrix of aggregate channel response structure based on the empty code channel number of maximum, this matrix is normal dimension matrix, utilize and to carry out joint-detection by normal dimension matrix, from joint detection results, extract the incoming symbol of the valid data of actual activated code channel correspondence then as decoding based on the empty code channel number of maximum.

In order to make those skilled in the art person understand the present invention program better, the present invention is described in further detail below in conjunction with drawings and embodiments.

At first, with reference to Fig. 3, Fig. 3 shows the baseband model of TD-SCDMA system:

Wherein, d ^(k)Be the emission data of code channel k, c ^(k)Be the spreading code of code channel k, h ^(k)Be the channel impulse response of code channel k, c ^(k)And h ^(k)Combination can regard aggregate channel sign indicating number b as ^(k), each user's data finally is superimposed, and then is superimposed with The noise, forms received signal e.

The structure of the basic physical channel of TD-SCDMA system is as shown in Figure 4: be made up of the identical data block (352 chips) of former and later two sizes, training sequence Miamble (144 chips) placed in the middle and the GP (protection at interval) of 16 chips.Burst is carried by time slot, burst block of a time slot carrying.Joint-detection in the synchronous time division system is by known Midamble is carried out channel estimating, many code channels combined signal handled carried out.

According to above-mentioned baseband model, total received signal is: e=Ad+n.

Wherein, d is the data symbol sequence of emission, and A is the matrix relevant with channel impulse response h with spreading code c, and n is a noise, and channel impulse response h can estimate to draw according to the training sequence code in the received signal (Midamble).Can estimate the d sequence thus.

According to code channel activating situation and channel impulse response, can obtain the aggregate channel response of multiple user signals.In the present invention, utilize the aggregate channel response matrix A of this aggregate channel response structure based on maximum code channel number _Max, in construction process, fill in 0 value on the position with non-activated code channel correspondence, then according to A _MaxMultiple user signals is carried out joint-detection, though comprised the data of non-activated code channel in the joint detection results that obtains like this, be 0 value owing to what fill in, therefore, can from this result, extract the valid data of actual activated code channel correspondence, the testing result that obtains expecting.

With reference to Fig. 5, Fig. 5 shows the realization flow of the inventive method, may further comprise the steps:

Step 501: the maximum empty code channel number that obtains system.

The present technique field personnel know, in the TD-SCDMA system, also will carry out data-modulated and band spectrum modulation through the digital data stream after the physical channel mapping, data-modulated is that continuous several bits are mapped as a symbol, and the complex symbol after the data-modulated carries out band spectrum modulation again.Band spectrum modulation mainly is divided into spread spectrum and two steps of scrambling.At first with spreading code to the data-signal spread spectrum, its spread spectrum coefficient is between 1-16, scrambled code then is added to scrambler in the signal behind the spread spectrum.The effect of spreading code is the different user that is used for distinguishing in the same time slot, is called channel code again.In CDMA (code division multiple access) system, can use a plurality of channel codes.

When carrying out joint-detection,, regard the channel code of different spreading factors the combination of the virtual channel sign indicating number of several maximum spreading factors as, thereby a code channel can be decomposed into a plurality of empty code channels in order to simplify the complexity of calculating.

Step 502: the maximum empty code channel number according to system carries out joint-detection to it, obtains the joint detection results based on the empty code channel number of maximum.

At first, according to the intermediate sequence of received signal, calculate the aggregate channel response of each subscriber signal.

The present technique field personnel know, the training sequence in the burst structure is used to carry out channel estimating, measurement etc.In same sub-district, the intermediate code that different user adopted in the same time slot is produced after cyclic shift by a basic intermediate code.

The front is mentioned, and aggregate channel response matrix A is the matrix relevant with channel impulse response h with spreading code c, and channel impulse response h can estimate to draw according to the training sequence code in the received signal (Midamble).

Therefore, can obtain each user's channel impulse response h,, obtain the aggregate channel response of multiple user signals then according to code channel activating situation and each user's channel impulse response according to the training sequence in each user's burst structure in the received signal.In the present invention, utilize the aggregate channel response matrix A of this aggregate channel response structure based on the empty code channel number of maximum _Max, this matrix is (L _Data+ W _FL-1) * K _MaxThe normal dimension matrix of N; Wherein, K _MaxBe the empty code channel number of maximum, W _FLLong for channel estimation window, N is a symbolic number, L _Data=N * Q _Max, Q _MaxBe maximum spreading gain.In construction process, fill in 0 value on the position with non-activated code channel correspondence.

Then, utilization is based on the aggregate channel response matrix A of the empty code channel number of maximum _MaxCarry out joint-detection, obtain the user data symbol of output.When carrying out joint-detection, can adopt ZF-BLE (ZF data block linear equalization), MMSE-BLE (Minimum Mean Square Error data block linear equalization), ZF-BDFE (ZF data block decision feedback equalization), MMSE-BDFE methods such as (Minimum Mean Square Error data block decision feedback equalizations).After the user data symbol that obtains exporting, according to the empty code channel number of maximum the user data symbol of output is carried out demodulation again, can obtain joint detection results based on the empty code channel number of maximum.

Step 503:, obtain the multi-user association testing result of expectation from based on extracting actual activated code channel data the joint detection results of the empty code channel number of maximum.

The present invention is applicable to the realization based on multi-user association detection algorithms such as ZF-BLE, ZF-BDFE, MMSE-BLE, MMSE-BDFE, is that the invention will be further described for example below with ZF-BLE.

The ZF-BLE algorithm mainly based on formula noted earlier (1), is promptly separated linear equation:

$A^H{R}_n^{-1}A{\hat{d}}_{c,\mathrm{ZF}-\mathrm{BLE}}=A^H{R}_n^{-1}e.$

If K _MaxBe the empty code channel number of maximum (constant).

The implementation of prior art is for to activate empty code channel number K based on reality _Vc(parameter) finds the solution above-mentioned linear equation, promptly finds the solution K _VcThe change dimensional linear equation of * N (N is a symbolic number) is to obtain expected result ${\hat{d}}_{c,\mathrm{ZF}-\mathrm{BLE}};$

The present invention then is based on K _MaxAbove-mentioned linear equation is found the solution, promptly find the solution K _MaxThe normal dimensional linear equation of * N (N is a symbolic number) is again with the gained result

Carry out to obtain expected result after the simple process

If A _MaxBe (L _Data+ W _FL-1) * K _MaxThe normal dimension matrix of N is made of N sub-piece B and matrix element 0, as shown in Figure 6:

Wherein, sub-piece B is (Q _Max+ W _FL-1) * K _MaxNormal dimension matrix, valid data are at its preceding K _VcRow, back (K _Max-K _Vc) columns is according to filling out 0.

Be (L _Data+ W _FL-1) * (L _Data+ W _FL-1) normal dimension symmetric positive definite matrix, so Be K _MaxN * (L _Data+ W _FL-1) normal dimension matrix, as shown in Figure 7:

Wherein, sub-piece B _ArgjBe K _Max* (L _Data+ W _FL-1) normal dimension matrix.

Then

Be K _MaxN * K _MaxThe normal dimension positive definite symmetric matrices of N, as shown in Figure 8:

Wherein, sub-piece H _iBe K _Max* K _MaxThe normal dimension matrix of N.

With H _i Back (K on the diagonal _Max-K _Vc) individual data fill out any non-zero values, λ ₀, λ ₁... λ _n, n=K _Max-K _Vc-1, λ _i≠ 0, make

Become non-singular matrix.

Because receiving data e is (L _Data+ W _FL-1) normal dimensional vector * 1, therefore,

Be K _MaxThe normal dimensional vector of N * 1, as shown in Figure 9.

Thus, separate linear equation:

${A_{\max }}^H{R}_n^{-1}{A}_{\max }{\hat{d}\′}_{c,\mathrm{ZF}-\mathrm{BLE}}={A_{\max }}^H{R}_n^{-1}e---\left(4\right)$

Can try to achieve user data symbol based on the empty code channel number output of maximum Be K _MaxThe normal dimensional vector of N * 1, its structure as shown in figure 10.

Contrast K shown in Figure 10 _MaxThe normal dimensional vector of N * 1

With K shown in Figure 2 _VcThe change dimensional vector of N * 1 Can see

In comprised In all data

I=0,1...N-1, therefore can from

In every K _MaxK before the extracting data _VcIndividual data are desired multi-user association testing result

MMSE-BLE is mainly based on formula (2).When handling, need find the solution following linear equation based on the empty code channel number of maximum:

$({A_{\max }}^H{R}_n^{-1}{A}_{\max }+R_d^{-1}){\hat{d}}_{c,\mathrm{MMSE}-\mathrm{BLE}}={A_{\max }}^H{R}_n^{-1}e,---\left(5\right)$

Wherein,

Being the inverse matrix of intersymbol noise correlation matrix, is K _MaxN * K _MaxThe positive definite symmetrical matrix of N.Therefore, according to above description, be equally applicable to the MMSE-BLE algorithm based on the implementation of the empty code channel number of maximum for ZF-BLE.

ZF-BDFE is that ZF-BLE adds decision-feedback, and MMSE-BLE is that MMSE-BLE adds decision-feedback, therefore according to foregoing description, applies the present invention to ZF-BDFE and MMSE-BLE, and it realizes that principle is similar, does not repeat them here.

As seen, the present invention carries out multi-user association based on the empty code channel number of maximum and detects in mobile radio system, and in implementation procedure, each matrix of generation is normal dimension matrix, is convenient to the realization of processor and ASIC circuit.Apply the present invention to processor and ASIC circuit based on the parallel computation framework, owing to only need to design at finding the solution normal dimensional linear equation, therefore can reduce the complexity of design difficulty and realization greatly, reduce the code expense, and can maximally utilise computation capability, thereby reduce time overhead based on the processor or the ASIC circuit of parallel computation framework at the design of finding the solution normal dimensional linear equation.In addition, the processing time of basic fixed also helps reducing delay variation, and simple relatively design helps increasing the stability of program.Apply the present invention to processor, mean cycle-index for fixing it owing to find the solution normal dimensional linear equation, and cycle-index is a fixed value, must bring better program optimization effect, therefore reduces size of code, reduces time overhead based on pipelined architecture.Therefore the fixing stability that has guaranteed the united detection processing time delay substantially of cycle-index reduced delay variation.

Though described the present invention by embodiment, those of ordinary skills know, the present invention has many distortion and variation and do not break away from spirit of the present invention, wish that appended claim comprises these distortion and variation and do not break away from spirit of the present invention.

Claims (6)

1, a kind of multi-user combined detection method is characterized in that, comprises step:

A, carry out many code channels signal according to the intermediate sequence of received signal and estimate, obtain the aggregate channel response of multiple user signals;

B, according to the aggregate channel response matrix A of described aggregate channel response structure based on the empty code channel number of maximum _Max, described A _MaxBe normal dimension matrix;

C, utilize described aggregate channel response matrix A based on the empty code channel number of maximum _MaxCarry out joint-detection, obtain the user data symbol of output;

D, according to the empty code channel number of maximum to output user data symbol carry out demodulation, obtain joint detection results based on the empty code channel number of maximum;

E, from based on the actual activated code channel data that extract the joint detection results of the empty code channel number of maximum in each maximum empty code channel data, obtain the multi-user association testing result of expectation.

2, method according to claim 1 is characterized in that, described aggregate channel response matrix A based on the empty code channel number of maximum _MaxBe (L _Data+ W _FL-1) * K _MaxThe normal dimension matrix of N; Wherein, K _MaxBe the empty code channel number of maximum, W _FLLong for channel estimation window, N is a symbolic number, L _Data=N * Q _Max, Q _MaxBe maximum spreading gain.

3, method according to claim 1 is characterized in that, described step B comprises:

Obtain each submatrix in the aggregate channel response matrix;

Valid data in each submatrix are rearranged, be specially: the preceding K that the data of activated code channel aggregate channel impulse response correspondence in each submatrix is placed on this submatrix _VcOn the individual position, all the other positions fill in 0, described K _VcFor reality activates empty code channel number;

Form aggregate channel response matrix A by each submatrix after rearranging based on the empty code channel number of maximum _Max

4, method according to claim 3 is characterized in that, described step C is specially:

C1, employing ZF data block linear equalizing method are carried out joint-detection to multiple user signals; Perhaps

C2, employing Minimum Mean Square Error data block linear equalizing method are carried out joint-detection to multiple user signals; Perhaps

C3, employing ZF data block decision feedback equalization method are carried out joint-detection to multiple user signals; Perhaps

C4, employing Minimum Mean Square Error data block decision feedback equalization method are carried out joint-detection to multiple user signals.

5, method according to claim 4 is characterized in that, described step C1 is specially: find the solution linear equation: ${A_{\max }}^H{R}_n^{-1}{A}_{\max }{\hat{d}\′}_{c,\mathrm{ZF}-\mathrm{BLE}}={A_{\max }}^H{R}_n^{-1}e,$ The user data symbol that obtains exporting Wherein, Expression interchannel noise correlation inverse matrix is (L _Data+ W _FL-1) * (L _Data+ W _FL-1) normal dimension symmetric positive definite matrix, e is (L for receiving data _Data+ W _FL-1) normal dimensional vector * 1,

Be K _MaxThe normal dimensional vector of N * 1.

6, method according to claim 4 is characterized in that, described step C2 is specially:

Find the solution linear equation: ${(A_{\max }}^H{R}_n^{-1}{A}_{\max }+R_d^{-1}{\widehat{)d}\′}_{c,\mathrm{MMSE}-\mathrm{BLE}}={A_{\max }}^H{R}_n^{-1}e,$ The user data symbol that obtains exporting Wherein,

Expression interchannel noise correlation inverse matrix,

The inverse matrix of expression intersymbol noise correlation matrix is K _MaxN * K _MaxThe normal dimension symmetric positive definite matrix of N, e is (L for receiving data _Data+ W _FL-1) normal dimensional vector * 1.

Images