# Patent

## Title
Method of double weighting parallel interference cancellation

## Abstract

A method for double weighting parallel interference cancellation in CDMA mobile communication system, integrates ideas of both partial weighting method and weighting method based on Bayes rule. When making a decision on each symbol of the user, calculating reliability coefficient of a decision result on the symbol of the user according to the weighting algorithm based on Bayes rule; weighting regenerated signal of the symbol of the user with the reliability coefficient. During process of MAI estimating and removing, obtaining the MAI on the expected user from weighted regenerated signal of other users in chip level; and then setting a weight value, weighting the MAI with the weight value. Finally, removing weighted MAI from received signal, which means partially removing MAI produced by other users on the expected user. At the same time, the present invention also discloses a double weighting PIC method with simplified algorithm, which transfers weighting in chip level to weighting in symbol level, replacing a Hyperbolic tangent decision with a piecewise linear decision method or a look-up table method.

## Background
FIELD OF THE TECHNOLOGY
The present invention generally relates to multiple user detection (MUD) technique in CDMA mobile communication system, and more particularly, to a method for double weighting parallel interference cancellation in CDMA system.
BACKGROUND OF THE INVENTION
Because of the advantages in great capacity, high service quality and excellent secrecy property, CDMA mobile communication system has become tendency of 3G mobile communication developments. Wherein, multiple user detection (MUD) is one kind of enhanced technique that can overcome limitations on CDMA system capacity produced by multiple access interference (MAI), in order to improve capacity and performance of the CDMA system.
With help of the information of multi-users, MUD technique can jointly detect signals from multi-users, in order to reduce influences of MAI on performance of the receiver, and heighten capacity of the system as much as possible. At present, MUD comprises maximum likelihood sequence detector that is the best detector, linear multiple user detector and interference cancellation multiple user detector that are second best detectors. Wherein, interference cancellation multiple user detection comprises the step of regarding signal from expected user as useful signal, removing and interference of other users from received signals to obtain signal of expected user, and then detecting the signal of the said expected user. In this way, performance of the system can be largely raised.
The interference cancellation multiple user detection can be classified into a serial interference cancellation (SIC) and a parallel interference cancellation (PIC). SIC sorts the user signals in power descending order and cancels interference in serial, the method has a better performance compared with single user detection but with longer time delay, and must sort signals of the users according to their power, which needs heavy calculation and is sensitive to the original signals estimation. PIC removes interference of other users in parallel from received signals, which has the advantages of short time delay and simpler calculation.
Prior PIC methods comprise traditional PIC method, partial PIC method and weighting PIC method that is based on Bayes rule.
Compared with single-user detection, the traditional PIC method can raise system performance in a large extent under high signal-to-noise ratio, but has a lower raise under low signal-to-noise ratio.
Being different from the traditional PIC method which removes MAI influence on expected users completely from received signals, the partial PIC method sets a weight value for each stage of interference cancellation to weight the MAI influence on the expected users, and partially removes the MAI during the interfere cancellation process. The traditional PIC method within Gaussian channel removes all the MAI on the expected users from the received signals, of course, the signal estimation on the expected users under this circumstance is biased; on the other hand, partial PIC method merely removes the MAI partially, which can correct estimation biases on the expected users, making decision results more reliable. At circumstance of low signal-to-noise ratio, partial PIC method has obviously better performance than that of the traditional PIC method.
The patent U.S. Pat. No. 5,418,814 discloses the weighting PIC method that is based on Bayes rule. Although being a weighting method, its weighting principle, which is a symbol-level weighting method that is based on a minimum mean of decision cost, is different from that of the partial PIC method. When creating a decision cost function, the method takes the minimum mean of decision cost as rule, to determine reliability coefficient of the decision result of each symbol, and to make symbol level weighting with the coefficient on the signals regenerated from the symbol. So only part of the interference produced by the symbol of the users is removed during the MAI eliminating process. The method has better performance than that of the traditional PIC method, especially at circumstance of low signal-to-noise ratio, and its improvement in performance is perfectly obvious.
Although both of the two methods above effectively improve performance of traditional PIC method under low signal-to-noise ratio, the extent is very limited.
SUMMARY OF THE INVENTION
It is an object to provide a method for double weighting parallel interference cancellation in CDMA communication system, which improves MUD performance without largely increasing its complexity, especially solves disadvantages of the prior art under circumstances with low signal-to-noise ratio.
A method for double weighting parallel interference cancellation in CDMA system according to the present invention comprises following steps of:
(a) Making a multi-path despreading, channel estimation and a multi-path combining on an input signal of a user by a RAKE receiver;
(b) Making a hard decision on a result of the multi-path combining;
(c) Calculating a reliability coefficient for a result of the hard decision of each symbol by the multi-path combining result and a value of the channel estimation;
(d) Regenerating a weighted signal in chip level for the user by the hard decision result, the reliability coefficient and the channel estimation value;
(e) Estimating a multiple access interference on each user by the weighted regenerated signal in chip level inputted in parallel of each user, and setting a weight value for a interference cancellation in different stages, weighting the multiple access interference on a expected user;
(f) Subtracting the weighted multiple access interference on the expected user from a base band signal of a received signal, and obtaining an output signal of the said expected user, which is also an input signal of the said expected user in the next stage of parallel interference cancellation.
Combining the partial PIC and the weighting PIC based on Bayes rule, bringing forward double weighting PIC method, the method disclosed according to the present invention has advantages of not only weighting algorithm based on Bayes rule, which has minimum decision cost in symbol level, but also partial weighting algorithm, which can make up the statistical signal estimation biases on users. What is more, the method enormously improves gains at circumstance of low signal-to-noise ratio without large increases on calculation workload comparing with the weighting PIC based on Bayes rule, and achieves a large improvement on performance of the system according to the present invention when comparing with partial weighting and weighting methods based on Bayes rule.
It is also a further object to provide a simplified double weighting PIC method, in order to reduce complexity of the algorithm while keeping performance of said double weighting PIC method unchanged.
A further object of the invention is to provide a double weighting PIC method with a simplified algorithm, which reduces complexity of the algorithm while keeping performance of said double weighting PIC method unchanged.
A double weighting PIC method with a simplified algorithm according to the present invention comprises following steps:
(a) Making a multi-path despreading, channel estimation and a multi-path combining on an input signal of a user by a RAKE receiver;
(b) Making a soft decision on a result of the multi-path combining and a value of the channel estimation;
(c) Setting a weight value for the soft decision in different stages, and weighting the soft decision in symbol level;
(d) Regenerating a weighted signal of the user in chip level by the weighted soft decision result and the channel estimation value;
(e) Estimating multiple access interference on each user by the weighted regenerated signal in chip level inputted in parallel of each user;
(f) Subtracting the multiple access interference on an expected user from a base band signal of a received signal, and obtaining an output signal of the expected user, which is also an input signal of the said expected user in the next stage of parallel interference cancellation.
The double weighting PIC method with the simplified algorithm according to the present invention, changing a chip level weighting into symbol level weighting, which can reduce complexity of the algorithm while keeping performance of said double weighting PIC method unchanged.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a schematic block diagram of multiple stages architecture in double weighting PIC receiver.
FIG. 2 is a schematic block diagram of PIC structure in double weighting PIC method.
FIG. 3 is a schematic block diagram of PIC structure of the last stage.
FIG. 4 is a schematic block diagram of PIC structure in simplified double weighting PIC method.



EMBODIMENTS OF THE INVENTION
The present invention will be described in more detail hereinafter with reference to the accompanying drawings and embodiments.
The method of the present invention integrates both partial weighting and the weighting based on Bayes rule. When decisions on each symbol of the user are made, a reliability coefficient of decision result of the symbol is calculated according to the weighting algorithm formula based on Bayes rule, and then regenerated signal of the symbol of the user is weighted with the reliability coefficient. During the MAI estimating and removing processes, the MAI on expected users produced by other users is got by the weighted regenerated signal of other users in chip level; after a weight value is set, the obtained MAI is weighted with the weight value. Finally, the weighted MAI is subtracted from received signals, which means partially removing the MAI on expected users produced by other users. Now above method is applied in each stage of PIC architecture.
As shown in FIG. 1, multistage structure in double weighting PIC method is the same with that of traditional PIC method. The first stage PIC structure 1 takes base band signals r(t) of the received signal as input of each user for processing, and the obtained output of each user are input of each user in PIC structure of the next stage; after processing on input of each user in the second stage PIC structure, the obtained output become input of each user in PIC structure of the next stage; stage by stage in this way, after processing on input of each user in PIC structure 2 of the last stage, the obtained output from each user here will be the final result of the multistage PIC structure.
As shown in FIG. 2, PIC structure in double weighting PIC method is the same with that of the weighting PIC method based on Bayes rule. PIC structure of the last stage of the method is same as that of the traditional PIC method, which is shown in FIG. 3.
Now refer to FIG. 1, which shows base band signals r(t) of the received signals enter the first stage PIC structure 1 in parallel mode. As shown in FIG. 2, input signal entering into PIC structure in parallel is transmitted to RAKE receiver 3 of each user respectively. RAKE receiver 3 despreads the input signal first, then makes channel estimation according to the result of despreading, and makes multi-path combining at last, then transmits the multi-path combining result to hard decision maker 4 and decision reliability calculator 7 simultaneously, and sends result of the channel estimation to decision reliability calculator 7 and signal regenerator 5 simultaneously. Hard decision maker 4 makes hard decision for the input, and transmits the result to signal regenerator 5.
In fading channel environment, base band signals r(t) of the received signal can be expressed as:




  
    
      
        
          r
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              
                ∑
                
                  i
                  =
                  1
                
                N
              
              ⁢
              
                
              
              ⁢
              
                
                  ∑
                  
                    l
                    =
                    1
                  
                  L
                
                ⁢
                
                  
                
                ⁢
                
                  
                    a
                    il
                  
                  ⁢
                  
                    
                      S
                      i
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        
                          τ
                          il
                        
                      
                      )
                    
                  
                
              
            
            +
            
              Z
              ⁡
              
                (
                t
                )
              
            
          
          =
          
            
              
                ∑
                
                  i
                  =
                  1
                
                N
              
              ⁢
              
                
              
              ⁢
              
                
                  ∑
                  
                    l
                    =
                    1
                  
                  L
                
                ⁢
                
                  
                
                ⁢
                
                  
                    a
                    il
                  
                  ⁢
                  
                    
                      P
                      i
                    
                  
                  ⁢
                  
                    
                      b
                      i
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        
                          τ
                          il
                        
                      
                      )
                    
                  
                  ⁢
                  
                    
                      c
                      i
                    
                    ⁡
                    
                      (
                      
                        t
                        -
                        
                          τ
                          il
                        
                      
                      )
                    
                  
                
              
            
            +
            
              Z
              ⁡
              
                (
                t
                )
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




Where, r(t) refers to the base band signals of the received signals; ail refers to a channel fading value of the ith user in the lth path, L is a path number; τil refers to a delay of the ith user in the lth path; Si(t) refers to a outbound signal of user i, N refers to total users number; Pi refers to a power of user i; bi(t) refers to a symbol flow of user




  i
  ,
  
    
      
        b
        i
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        ∑
        
          m
          =
          
            -
            ∞
          
        
        ∞
      
      ⁢
      
        
      
      ⁢
      
        
          a
          i
          
            (
            m
            )
          
        
        ⁢
        
          p
          ⁡
          
            (
            
              t
              -
              
                mT
                b
              
            
            )
          
        
      
    
  
  ,
  
    a
    i
    
      (
      m
      )
    
  




refers to the mth symbol of the ith user, p(t) refers to a signal pulse whose cycle is Tb, when there is no hindrance on the deducing process with the algorithm, set p(t) as rectangular pulse (when t∈[0, Tb], p(t)=1; when t∉[0, Tb], p(t)=0); ci(t) refers to a spread-spectrum code of user i; Z(t) refers to channel noise.

In the kth stage PIC structure, input of RAKE receiver 3 of user i is




  
    
      r
      i
      
        (
        k
        )
      
    
    ⁡
    
      (
      t
      )
    
  
  .




When





  
    k
    =
    1
  
  ,
  
    
      
        r
        i
        
          (
          1
          )
        
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        r
        ⁡
        
          (
          t
          )
        
      
      .
    
  




After the RAKE receiver makes multi-path despreading on ri(k)(t), and channel estimation according to the despreading result, multi-path combining is implemented. The lth path despreading result made by RAKE receiver of user i is:




  
    
      
        
          
            
              y
              i
              
                
                  (
                  m
                  )
                
                ⁢
                
                  (
                  k
                  )
                
              
            
            ⁡
            
              (
              l
              )
            
          
          =
          
            
              1
              
                
                  T
                  b
                
              
            
            ⁢
            
              
                ∫
                
                  
                    
                      (
                      
                        m
                        -
                        1
                      
                      )
                    
                    ⁢
                    
                      T
                      b
                    
                  
                  +
                  
                    τ
                    il
                  
                
                
                  
                    mT
                    b
                  
                  +
                  
                    τ
                    il
                  
                
              
              ⁢
              
                
                  ∂
                  
                    
                      r
                      i
                      
                        (
                        k
                        )
                      
                    
                    ⁡
                    
                      (
                      t
                      )
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  
                    c
                    i
                    *
                  
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        τ
                        il
                      
                    
                    )
                  
                
                ⁢
                
                  
                
                ⁢
                
                  d
                  t
                
              
            
          
        
        ⁢
        
          
        
      
    
    
      
        (
        2
        )
      
    
  




Where, l=1, . . . , L.
With the maximal ratio combining, the multi-path combining result made by RAKE receiver will be:




  
    
      
        
          
            Y
            i
            
              
                (
                m
                )
              
              ⁢
              
                (
                k
                )
              
            
          
          =
          
            Re
            ⁢
            
              {
              
                y
                i
                
                  
                    (
                    m
                    )
                  
                  ⁢
                  
                    (
                    k
                    )
                  
                
              
              }
            
          
        
        ⁢
        
          
        
      
    
    
      
        (
        3
        )
      
    
  
  
    
      
        Where
        ,
        
          
            y
            i
            
              
                (
                m
                )
              
              ⁢
              
                (
                k
                )
              
            
          
          =
          
            
              ∑
              
                l
                =
                1
              
              L
            
            ⁢
            
              
            
            ⁢
            
              
                A
                il
                *
              
              ⁢
              
                
                  y
                  i
                  
                    
                      (
                      m
                      )
                    
                    ⁢
                    
                      (
                      k
                      )
                    
                  
                
                ⁡
                
                  (
                  l
                  )
                
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




Ail is an estimation value of ail√{square root over (Pi)}, ail is a channel fading value in the lth path of the ith user, Pi refers to a power of user i.
The above multi-path combining results of the RAKE receiver can be expressed as




  
    
      Y
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        
          μ
          i
        
        ⁢
        
          a
          i
          
            (
            m
            )
          
        
      
      +
      
        n
        i
      
    
  
  ,




wherein, ni is Gaussian White Noise, which obeys normal distribution N(0, σi2), σ2 is a noise power of





  
    n
    i
  
  ;
  
    a
    i
    
      (
      m
      )
    
  




refers to the mth symbol of user i, whose value is +1 or −1, μi is a real number relating to channel fading.

The decision result for the mth symbol of the ith user is:




  
    
      
        
          
            a
            ^
          
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        =
        
          sgn
          ⁢
          
            {
            
              Y
              i
              
                
                  (
                  m
                  )
                
                ⁢
                
                  (
                  k
                  )
                
              
            
            }
          
        
      
    
    
      
        (
        5
        )
      
    
  




The decision reliability calculator calculates the reliability coefficient of the decision result made by hard decision maker through two input signals, and transmits the reliability coefficient to signal regenerator 5. The reliability coefficient calculating formula of âi(m)(k) is expressed as:




  
    
      
        
          f
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        =
        
          tanh
          ⁢
          
            {
            
              w
              ⁢
              
                
                  
                    μ
                    i
                  
                  ⁢
                  
                    
                    
                      Y
                      i
                      
                        
                          (
                          m
                          )
                        
                        ⁢
                        
                          (
                          k
                          )
                        
                      
                    
                    
                  
                
                
                  σ
                  i
                  2
                
              
            
            }
          
        
      
    
    
      
        (
        6
        )
      
    
  




Where, w is a positive real number, which is used to make up errors in noise power estimation. The w values in different signal-to-noise ratios can be determined via experiments, wherein w=1 at higher signal-to-noise ratios.
After obtaining regenerated signals of the user from three input signals, the signal regenerator transmits the regenerated signals to device 6 for MAI estimation and partial interference cancellation. Weighted regenerated signal in chip level of user i is:




  
    
      
        
          
            g
            i
            
              (
              k
              )
            
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∑
            
              l
              =
              1
            
            L
          
          ⁢
          
            
          
          ⁢
          
            
              A
              il
            
            ⁢
            
              
                ∑
                
                  n
                  =
                  
                    -
                    ∞
                  
                
                ∞
              
              ⁢
              
                
              
              ⁢
              
                
                  f
                  i
                  
                    
                      (
                      n
                      )
                    
                    ⁢
                    
                      (
                      k
                      )
                    
                  
                
                ⁢
                
                  
                    a
                    ^
                  
                  i
                  
                    
                      (
                      n
                      )
                    
                    ⁢
                    
                      (
                      k
                      )
                    
                  
                
                ⁢
                
                  p
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        nT
                        b
                      
                      -
                      
                        τ
                        il
                      
                    
                    )
                  
                
                ⁢
                
                  
                    c
                    i
                  
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        τ
                        il
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




Where, p(t) is a signal pulse whose cycle is Tb.
As shown in FIG. 2, the base band signals r(t) of the received signals are also transmitted into MAI estimation and partial interference cancellation device 6, which estimates MAI on each user according to the regenerated signals of each user inputted in parallel; the resulted signal after removing MAI on certain user from the base band signals r(t) of the received signals will be the output of the said user in the present stage PIC structure and the input of the said user in next stage of PIC structure.
In the kth stage PIC structure, MAI estimation of user i will be:




  
    
      
        
          
            I
            ^
          
          i
          
            (
            k
            )
          
        
        =
        
          
            ∑
            
              
                j
                =
                1
              
              ,
              
                j
                ≠
                i
              
            
            N
          
          ⁢
          
            
          
          ⁢
          
            
              g
              j
              
                (
                k
                )
              
            
            ⁡
            
              (
              t
              )
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




Suppose weight value of the kth stage PIC structure is p(k); the MAI from formula (8) can be weighted and then interference cancellation is made to obtain output ri(k+1)(t) of user i in the kth stage PIC structure according to formula (9). ri(k+1)(t) is the input of RAKE receiver of user i in PIC structure of the next stage.




  
    
      
        
          
            r
            i
            
              (
              
                k
                +
                1
              
              )
            
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            r
            ⁡
            
              (
              t
              )
            
          
          -
          
            
              p
              
                (
                k
                )
              
            
            ⁢
            
              
                I
                ^
              
              i
              
                (
                k
                )
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




Different weight value p(k) can be set up for interference cancellations in different stages, and preferably set p(1)<p(2) . . . <p(k) . . . <p(S), wherein k is the kth stage of the interference cancellation, S refers to a stage number of the interference cancellation.
Same process will be done for parallel input signals in PIC structure of the next stage. The process goes in this way stage by stage, and when it goes PIC structure of the last stage, the parallel input signals will be transmitted respectively into RAKE receivers 3 of each user. After multi-path dispreading, channel estimation and multi-path combining on the input signal made by RAKE receiver of the user, soft output of the user can be got. Said soft output of each user will be final result of each user in multistage PIC structure. In other words, multi-path dispreading process according to formula (2), and multi-path combining according to formula (3), (4) are made on the input signal of user i. And the soft output of user i obtained through multi-path combining will be final result of user i in the multistage PIC structure. In the receiver, soft output of the user will be transmitted to its decoder for decoding.
Integrating the partial PIC and the weighting PIC based on Bayes rule, the present invention provides the double weighting PIC method, which has advantages of not only the weighting algorithm based on Bayes rule, which has advantages of Minimum decision Cost in symbol level, but also advantages of the partial weighting algorithm, which can make up the statistical signal estimation biases on users. At the same time, the method has better performance than that of weighting PIC based on Bayes rule, which can improve gains at circumstance of signal-to-noise ratio without obvious increases in calculation workload and provide a largely improved performance compared with both partial weighting and weighting based on Bayes rule.
Through analysis on calculation formula of said double weighting PIC method, it is noticed that, formula (6) comprises the Hyperbolic tangent operation, which is hard to implemented with prior hardware, and needs to develop new arithmetic chips, so the cost is much higher; moreover,




  
    p
    
      (
      k
      )
    
  
  ⁢
  
    
      I
      ^
    
    i
    
      (
      k
      )
    
  




in formula (9) is chip level multiplication with heavy calculation workload.

In order to solve above disadvantages, the present invention further provides a double weighting PIC method with simplified algorithm, whose multistage architecture is shown in FIG. 1, PIC structure of the simplified algorithm is illustrated in FIG. 4, PIC structure of the simplified algorithm in the last stage is shown in FIG. 3.
Here is one of embodiments with simplified algorithm according to the present invention:
As shown in FIG. 1, base band signals r(t) of the received signals enter in parallel to the first stage PIC structure 1 in FIG. 1. Referring to FIG. 4, the input signals r(t) entered the PIC structure are transmitted into RAKE receiver 3 of each user respectively. RAKE receiver 3 implements channel estimation after dispreading the input signal, and completes multi-path combining at least. RAKE receiver 3 transmits the multi-path combining result to soft decision maker 8, and then at the same time sends the channel estimation result to soft decision maker 8 and signal regenerator 5.
In the kth stage PIC structure, multi-path combining result of user i can be presented by:




  
    
      
        
          Y
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        =
        
          
            
              μ
              i
            
            ⁢
            
              a
              i
              
                (
                m
                )
              
            
          
          +
          
            n
            i
          
        
      
    
    
      
        (
        10
        )
      
    
  




Where, ni is Gaussian White Noise, which is subject to normal distribution




  
    N
    ⁡
    
      (
      
        0
        ,
        
          σ
          i
          2
        
      
      )
    
  
  ;
  
    a
    i
    
      (
      m
      )
    
  




is the mth symbol of user i, whose value is +1 or −1, μi is a real number relating to channel fading.

Soft decision maker 8 makes soft decision on the input signal, the soft decision on multi-path combining result in RAKE receiver of user i is




  
    
      ξ
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        f
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
      ⁢
      
        
          a
          ^
        
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
    
  
  ,




and





  
    f
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  
  ⁢
  
    
      a
      ^
    
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  




satisfies:





  
    
      
        
          
            f
            i
            
              
                (
                m
                )
              
              ⁢
              
                (
                k
                )
              
            
          
          ⁢
          
            
              a
              ^
            
            i
            
              
                (
                m
                )
              
              ⁢
              
                (
                k
                )
              
            
          
        
        =
        
          tanh
          ⁢
          
            {
            
              w
              ⁢
              
                
                  
                    μ
                    i
                  
                  ⁢
                  
                    Y
                    i
                    
                      
                        (
                        m
                        )
                      
                      ⁢
                      
                        (
                        k
                        )
                      
                    
                  
                
                
                  σ
                  
                    i
                    2
                  
                
              
            
            }
          
        
      
    
    
      
        (
        11
        )
      
    
  




Where,




  
    a
    ^
  
  i
  
    
      (
      m
      )
    
    ⁢
    
      (
      k
      )
    
  




is the decision result on the mth symbol of user i,





  
    
      
        a
        ^
      
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      sgn
      ⁢
      
        {
        
          Y
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        }
      
    
  
  ,
  
    f
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  




is reliability coefficient of





  
    
      a
      ^
    
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  
  .




Replacing the Hyberbolic tangent decision of said double weighting PIC method with a piecewise linear decision, setting a piecewise linear decision function as L(x), the method of replacing the Hyperbolic tangent decision with the piecewise linear decision is to apply the piecewise linear decision function L(x) to approach the Hyperbolic tangent function tanh(x). Deducting process of the piecewise linear decision function is as follows.
{circle around (1)}. Define the piecewise linear decision function L(x)
Because the Hyperbolic tangent function is an odd function: tanh(−x)=−tanh(x);, define

L(−x)=−L(x).

{circle around (2)}. Determine a threshold value θ
The Hyperbolic tangent function has following characteristics: when x→∞, tanh(x)→1; Therefore, set the threshold value θ>0, when x>θ, set L(x)=1;
{circle around (3)}. Determine a linear parameter Q
When 0≦x≦θ, equally divide field [0, θ] into Q small intervals, the qth small interval is




  
    [
    
      
        x
        
          q
          -
          1
        
      
      ,
      
        x
        q
      
    
    ]
  
  ,
  
    
      x
      q
    
    =
    
      
        q
        ⁢
        
          
        
        ⁢
        θ
      
      Q
    
  
  ,
  
    
      x
      0
    
    =
    0
  
  ,
  
    
      
        x
        Q
      
      =
      θ
    
    ;
  




{circle around (4)}. The expression of L(x) within the qth small interval is as following.
In the interval [xq−1, xq], define L(x) as a line segment connecting between point Cq and Dq. Wherein, coordinate of point Cq is Cq=(xq−1, tanh(xq−1)), coordinate of point Dq is Dq=(xq, tanh(xq)). And approximate the tanh(x) curve within the interval [xq−1, xq] with said line segment CqDq. The equation of the line segment CqDq is:




  
    
      
        
          
            L
            q
          
          ⁡
          
            (
            x
            )
          
        
        =
        
          
            tanh
            ⁡
            
              (
              
                x
                
                  q
                  -
                  1
                
              
              )
            
          
          +
          
            
              
                
                  tanh
                  ⁡
                  
                    (
                    
                      x
                      q
                    
                    )
                  
                
                -
                
                  tanh
                  ⁡
                  
                    (
                    
                      x
                      
                        q
                        -
                        1
                      
                    
                    )
                  
                
              
              
                
                  x
                  q
                
                -
                
                  x
                  
                    q
                    -
                    1
                  
                
              
            
            ⁢
            
              (
              
                x
                -
                
                  x
                  
                    q
                    -
                    1
                  
                
              
              )
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




{circle around (5)}. With L(−x)=−L(x), the expression of L(x) in the interval [−xq, −xq−] is:

L(x)=−Lq(−x)

{circle around (6)}. The expression of the piecewise linear decision function L(x) is:




  
    
      
        
          L
          ⁡
          
            (
            x
            )
          
        
        =
        
          {
          
            
              
                
                  1
                  ,
                  
                    x
                    >
                    θ
                  
                
              
            
            
              
                
                  
                    
                      L
                      q
                    
                    ⁡
                    
                      (
                      x
                      )
                    
                  
                  ,
                  
                    x
                    ∈
                    
                      [
                      
                        
                          x
                          
                            q
                            -
                            1
                          
                        
                        ,
                        
                          x
                          q
                        
                      
                      ]
                    
                  
                
              
            
            
              
                
                  
                    -
                    
                      
                        L
                        q
                      
                      ⁡
                      
                        (
                        
                          -
                          x
                        
                        )
                      
                    
                  
                  ,
                  
                    x
                    ∈
                    
                      [
                      
                        
                          -
                          
                            x
                            q
                          
                        
                        ,
                        
                          -
                          
                            x
                            
                              q
                              -
                              1
                            
                          
                        
                      
                      ]
                    
                  
                
              
            
            
              
                
                  
                    -
                    1
                  
                  ,
                  
                    x
                    <
                    
                      -
                      θ
                    
                  
                
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




Soft decision maker 8 transmits the soft decision result to soft decision weighting device 9, which weights the soft decision result with formula (14), and sends the weighted result to signal regenerator 5. The weighting operation is symbol level.




  
    
      
        
          ρ
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        =
        
          
            ξ
            i
            
              
                (
                m
                )
              
              ⁢
              
                (
                k
                )
              
            
          
          ·
          
            p
            
              (
              k
              )
            
          
        
      
    
    
      
        (
        14
        )
      
    
  




Signal regenerator 5 obtains the regenerated signal from two input signals according to following formula, and transmits the regenerated signal to MAI estimation and partial interference cancellation device 6, the weighted regenerated signal in chip level of user i can be expressed with following formula:




  
    
      
        
          
            g
            i
            
              (
              k
              )
            
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∑
            
              l
              =
              1
            
            L
          
          ⁢
          
            
          
          ⁢
          
            
              A
              il
            
            ⁢
            
              
                ∑
                
                  n
                  =
                  
                    -
                    ∞
                  
                
                ∞
              
              ⁢
              
                
              
              ⁢
              
                
                  ρ
                  i
                  
                    
                      (
                      n
                      )
                    
                    ⁢
                    
                      (
                      k
                      )
                    
                  
                
                ⁢
                
                  
                    a
                    ^
                  
                  i
                  
                    
                      (
                      n
                      )
                    
                    ⁢
                    
                      (
                      k
                      )
                    
                  
                
                ⁢
                
                  p
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        nT
                        b
                      
                      -
                      
                        τ
                        il
                      
                    
                    )
                  
                
                ⁢
                
                  
                    c
                    i
                  
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        τ
                        il
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




As shown in FIG. 4, base band signals r(t) of the received signals are also transmitted to MAI estimation and partial interference cancellation device 6, which estimates MAI on each user through the regenerated signal of each user inputted in parallel, in the kth stage PIC algorithm, estimation of MAI on user i is:




  
    
      
        
          I
          i
          
            (
            k
            )
          
        
        =
        
          
            ∑
            
              
                j
                =
                1
              
              ,
              
                j
                ≠
                i
              
            
            N
          
          ⁢
          
            
          
          ⁢
          
            
              g
              j
              
                (
                k
                )
              
            
            ⁡
            
              (
              t
              )
            
          
        
      
    
    
      
        (
        16
        )
      
    
  




After MAI on user i is calculated, the method removes MAI on user i from base band signals r(t) of the received signals, and implements interference cancellation on MAI according to formula (17):




  
    
      
        
          
            r
            i
            
              (
              
                k
                +
                1
              
              )
            
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            r
            ⁡
            
              (
              t
              )
            
          
          -
          
            
              I
              ^
            
            i
            
              (
              k
              )
            
          
        
      
    
    
      
        (
        17
        )
      
    
  








  
    r
    i
    
      (
      
        k
        +
        1
      
      )
    
  
  ⁡
  
    (
    t
    )
  




is the output signal of user i in the kth stage PIC structure, and the input signal into RAKE receiver of user i in PIC structure of the next stage.

The signal that is obtained after removing MAI on the user from base band signals r(t) of the received signals will become output signal of the said user in the present stage PIC structure, as well as input signal of the said user in next stage PIC structure. In the next stage of PIC structure, the parallel input signals will be processed in the same way. The operation is carried out stage by stage in this way; when being processed in PIC structure of the last stage, as shown in FIG. 3, the parallel input signals are transmitted into RAKE receiver 3 of each user separately, which implements despreading, channel estimation and multi-path combining on the input signal. The soft output from each user here will be the final result of the multistage PIC structure. In the receiver, the soft output of the user is transmitted to decoder of the said user for decoding.
Another embodiment of simplified algorithm according to the present invention is as following.
Now refer to FIG. 1. The base band signals r(t) of the received signals are transmitted into the first stage PIC structure 1 in parallel. As shown in FIG. 4, said input signals r(t) are respectively transmitted into RAKE receiver 3 of each user, which despreads the input signal first and makes channel estimation, and last implements multi-path combining. RAKE receiver 3 transmits the multi-path combining result to soft decision maker 8, and sends channel estimation result to soft decision maker 8 and signal regenerator 5 at the same time. In the kth stage PIC structure, the multi-path combining result of user i can be expressed as:




  
    
      
        
          Y
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        =
        
          
            
              μ
              i
            
            ⁢
            
              a
              i
              
                (
                m
                )
              
            
          
          +
          
            n
            i
          
        
      
    
    
      
        (
        18
        )
      
    
  




ni is Gaussian White Noise, which is subject to normal distribution




  
    N
    ⁡
    
      (
      
        0
        ,
        
          σ
          i
          2
        
      
      )
    
  
  ;
  
    a
    i
    
      (
      m
      )
    
  




is the mth symbol of user i, whose value is +1 or −1. μi is a real number that is related with channel fading.

Soft decision maker 8 makes soft decision on the input signal, the soft decision on the multi-path combining result of RAKE receiver of user i is




  
    
      ξ
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        f
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
      ⁢
      
        
          a
          ^
        
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
    
  
  ,




where





  
    f
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  
  ⁢
  
    
      a
      ^
    
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  




satisfies following formula:





  
    
      
        
          
            f
            i
            
              
                (
                m
                )
              
              ⁢
              
                (
                k
                )
              
            
          
          ⁢
          
            
              a
              ^
            
            i
            
              
                (
                m
                )
              
              ⁢
              
                (
                k
                )
              
            
          
        
        =
        
          tanh
          ⁢
          
            {
            
              w
              ⁢
              
                
                  
                    μ
                    i
                  
                  ⁢
                  
                    Y
                    i
                    
                      
                        (
                        m
                        )
                      
                      ⁢
                      
                        (
                        k
                        )
                      
                    
                  
                
                
                  σ
                  i
                  2
                
              
            
            }
          
        
      
    
    
      
        (
        19
        )
      
    
  




Where,




  
    a
    ^
  
  i
  
    
      (
      m
      )
    
    ⁢
    
      (
      k
      )
    
  




is a decision result of the mth symbol of user i,





  
    
      
        a
        ^
      
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      sgn
      ⁢
      
        {
        
          Y
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        }
      
    
  
  ,
  
    f
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  




is a reliability coefficient of





  
    
      a
      ^
    
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  
  .




When replacing the Hyperbolic tangent decision in the double weighting PIC method with a look-up table method, setting decision function of the look-up table method as T(x), the method of replacing the Hyperbolic tangent decision with the look-up table method is to approach the Hyperbolic tangent function tanh(x) with the decision function T(x) of the look-up table method, whose deducing process is as follows.
{circle around (1)}. Define decision function T(x) of the look-up table method
Because the Hyperbolic tangent function is an odd function: tanh(−x)=−tanh(x), define

T(−x)=−T(x);

{circle around (2)}. Determine a threshold value θ
Because the Hyperbolic tangent function has the characteristics of: when x→∞, tanh(x)→1; set the threshold value θ>0 according to the present invention, when x>θ, set T(x)=1;
{circle around (3)}. Determine a linear parameter Q
When 0≦x≦θ, divide the field [0, θ] into Q small intervals equally, where the qth small interval is




  
    [
    
      
        x
        
          q
          -
          1
        
      
      ,
      
        x
        q
      
    
    ]
  
  ,
  
    
      x
      q
    
    =
    
      
        q
        ⁢
        
          
        
        ⁢
        θ
      
      Q
    
  
  ,
  
    
      x
      0
    
    =
    0
  
  ,
  
    
      
        x
        Q
      
      =
      θ
    
    ;
  




{circle around (4)}. Expression of the T(x) in small interval q is as following.
In the small interval [xq−1, xq], set midpoint of the small interval as




  
    
      
        x
        
          q
          -
          1
        
      
      +
      
        x
        q
      
    
    2
  
  ,




and define T(x) as follows:





  
    
      
        
          T
          ⁡
          
            (
            x
            )
          
        
        =
        
          tanh
          ⁡
          
            (
            
              
                
                  x
                  
                    q
                    -
                    1
                  
                
                +
                
                  x
                  q
                
              
              2
            
            )
          
        
      
    
    
      
        (
        20
        )
      
    
  




{circle around (5)}. By means of T(−x)=−T(x), expression of the T(x) in the interval [−θ, 0] can be obtained.
{circle around (6)}. Expression of the decision function T(x) in look-up table method is:




  
    
      
        
          1
          ,
          
            x
            >
            θ
          
        
        ⁢
        
          

        
        ⁢
        
          
            T
            ⁡
            
              (
              x
              )
            
          
          =
          
            {
            
              
                
                  
                    
                      
                        
                          tanh
                          ⁡
                          
                            (
                            
                              
                                
                                  x
                                  
                                    q
                                    -
                                    1
                                  
                                
                                +
                                
                                  x
                                  q
                                
                              
                              2
                            
                            )
                          
                        
                        ,
                        
                          x
                          ∈
                          
                            [
                            
                              
                                x
                                
                                  q
                                  -
                                  1
                                
                              
                              ,
                              
                                x
                                q
                              
                            
                            ]
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          -
                          
                            tanh
                            ⁡
                            
                              (
                              
                                
                                  
                                    x
                                    
                                      q
                                      -
                                      1
                                    
                                  
                                  +
                                  
                                    x
                                    q
                                  
                                
                                2
                              
                              )
                            
                          
                        
                        ,
                        
                          x
                          ∈
                          
                            [
                            
                              
                                -
                                
                                  x
                                  q
                                
                              
                              ,
                              
                                -
                                
                                  x
                                  
                                    q
                                    -
                                    1
                                  
                                
                              
                            
                            ]
                          
                        
                      
                    
                  
                
                ⁢
                
                  

                
                -
                1
              
              ,
              
                x
                <
                
                  -
                  θ
                
              
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




Soft decision maker 8 transmits the soft decision result to soft decision weighting device 9, which weights the soft decision result in symbol level with formula (22), and sends the weighted result to signal regenerator 5.




  
    
      
        
          ρ
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        =
        
          
            ξ
            i
            
              
                (
                m
                )
              
              ⁢
              
                (
                k
                )
              
            
          
          ·
          
            p
            
              (
              k
              )
            
          
        
      
    
    
      
        (
        22
        )
      
    
  




According to the formula below, signal regenerator 5 obtains the regenerated signal from two input signals, and transmits the regenerated signals to MAI estimation and partial interference cancellation device 6; the weighted regenerated signal of user i in chip level can be expressed as:




  
    
      
        
          
            g
            i
            
              (
              k
              )
            
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∑
            
              l
              =
              1
            
            L
          
          ⁢
          
            
          
          ⁢
          
            
              A
              il
            
            ⁢
            
              
                ∑
                
                  n
                  =
                  
                    -
                    ∞
                  
                
                ∞
              
              ⁢
              
                
              
              ⁢
              
                
                  ρ
                  i
                  
                    
                      (
                      n
                      )
                    
                    ⁢
                    
                      (
                      k
                      )
                    
                  
                
                ⁢
                
                  p
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        nT
                        b
                      
                      -
                      
                        τ
                        il
                      
                    
                    )
                  
                
                ⁢
                
                  
                    c
                    i
                  
                  ⁡
                  
                    (
                    
                      t
                      -
                      
                        τ
                        il
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        23
        )
      
    
  




As shown in FIG. 4, base band signals r(t) of the received signals are also transmitted into MAI estimation and partial interference cancellation device 6, which estimates MAI on each user based on the regenerated signal inputted in parallel of each user, in the kth stage PIC algorithm, the estimations of MAI on user i is:




  
    
      
        
          
            I
            ^
          
          i
          
            (
            k
            )
          
        
        =
        
          
            ∑
            
              
                j
                =
                1
              
              ,
              
                j
                ≠
                i
              
            
            N
          
          ⁢
          
            
          
          ⁢
          
            
              g
              j
              
                (
                k
                )
              
            
            ⁡
            
              (
              t
              )
            
          
        
      
    
    
      
        (
        24
        )
      
    
  




After MAI on user i is calculated, the method removes MAI on user i from base band signals r(t) of the received signals, and implements interference cancellation on MAI according to formula (25):




  
    
      
        
          
            r
            i
            
              (
              
                k
                +
                1
              
              )
            
          
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            r
            ⁡
            
              (
              t
              )
            
          
          -
          
            
              I
              ^
            
            i
            
              (
              k
              )
            
          
        
      
    
    
      
        (
        25
        )
      
    
  








  
    r
    i
    
      (
      
        k
        +
        1
      
      )
    
  
  ⁡
  
    (
    t
    )
  




is the output signal of user i in the kth stage PIC structure, and as well as input signal of RAKE receiver of user i in the next stage PIC structure.

The signal that is obtained after removing MAI on the user from base band signals r(t) of the received signals will become the output signal of the said user in the present stage PIC structure, as well as the input signal of the said user in next stage of PIC structure. In PIC structure of the next stage, the parallel input signals will be processed in the same way. The operation is carried out stage by stage in this way; when processed in PIC structure of the last stage, as shown in FIG. 3, the parallel input signals are transmitted into RAKE receiver 3 of each user separately, which implements despreading, channel estimation and multi-path combining on the input signal. The soft output from each user here will be the final result of the multistage PIC structure. In the receiver, the soft output of the user is transmitted to its decoder for decoding.

## Claims

1. A double weighting parallel interference cancellation method that can be used in a CDMA mobile communication system, comprises:
(a) Making a multi-path despreading, channel estimation and a multi-path combining on an input signal of a user by a RAKE receiver;
(b) Making a hard decision on a result of the multi-path combining;
(c) Calculating a reliability coefficient for a result of the hard decision of each symbol by the multi-path combining result and a value of the channel estimation;
(d) Regenerating a weighted signal in chip level for the user by the hard decision result, the reliability coefficient and the channel estimation value;
(e) Estimating a multiple access interference on each user by the weighted regenerated signal in chip level inputted in parallel of each user, and setting a weight value for a interference cancellation in different stages, weighting the multiple access interference on a expected user;
(f) Subtracting the weighted multiple access interference on the expected user from a base band signal of a received signal, and obtaining an output signal of the said expected user, which is also an input signal of the said expected user in the next stage parallel interference cancellation.

2. The double weighting parallel interference cancellation method according to claim 1, comprises, if the said input signal of user i in the kth stage parallel interference cancellation is expressed as ri(k)(t),
In step (a), calculating the said multi-path combining result of user i by a formula





  
    
      Y
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      Re
      ⁢
      
        {
        
          y
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        }
      
    
  
  ,



 which can be expressed as





  
    
      Y
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        
          μ
          i
        
        ⁢
        
          a
          i
          
            (
            m
            )
          
        
      
      +
      
        n
        i
      
    
  
  ,



 where, ni is a Gaussian White Noise, which is subject to a normal distribution of





  
    N
    ⁡
    
      (
      
        0
        ,
        
          σ
          i
          2
        
      
      )
    
  
  ,
  
    σ
    2
  



 is a noise power of





  
    n
    i
  
  ,
  
    a
    i
    
      (
      m
      )
    
  



 is the mth symbol of user i, whose value is +1 or −1, μi is a real number that is related with channel fading;
In step (b), calculating the said hard decision result of the mth symbol of user i by a formula





  
    
      
        a
        ^
      
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      sgn
      ⁡
      
        (
        
          Y
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        )
      
    
  
  ;



In step (c), calculating the said reliability coefficient of the said decision result on the mth symbol of user i by





  
    
      f
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      tanh
      ⁢
      
        {
        
          w
          ⁢
          
            
              
                μ
                i
              
              ⁢
              
                
                
                  Y
                  i
                  
                    
                      (
                      m
                      )
                    
                    ⁢
                    
                      (
                      k
                      )
                    
                  
                
                
              
            
            
              σ
              i
              2
            
          
        
        }
      
    
  
  ,



 where, w is a positive real number;
In step (d), calculating the said weighted regenerated signal of user i in chip level by a formula





  
    
      
        g
        i
        
          (
          k
          )
        
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        ∑
        
          l
          =
          1
        
        L
      
      ⁢
      
        
      
      ⁢
      
        
          A
          il
        
        ⁢
        
          
            ∑
            
              n
              =
              
                -
                ∞
              
            
            ∞
          
          ⁢
          
            
          
          ⁢
          
            
              f
              i
              
                
                  (
                  n
                  )
                
                ⁢
                
                  (
                  k
                  )
                
              
            
            ⁢
            
              
                a
                ^
              
              i
              
                
                  (
                  n
                  )
                
                ⁢
                
                  (
                  k
                  )
                
              
            
            ⁢
            
              p
              ⁡
              
                (
                
                  t
                  -
                  
                    nT
                    b
                  
                  -
                  
                    τ
                    il
                  
                
                )
              
            
            ⁢
            
              
                c
                i
              
              ⁡
              
                (
                
                  t
                  -
                  
                    τ
                    il
                  
                
                )
              
            
          
        
      
    
  
  ,



 where, Ail is a estimation value of ail√{square root over (Pi)}, ail refers to a channel fading value of user i in the path l, pi is a power of user i;
In step (e), calculating the said multiple access interference estimation of user i by a formula





  
    
      I
      i
      
        (
        k
        )
      
    
    =
    
      
        ∑
        
          
            j
            =
            1
          
          ,
          
            j
            ≠
            i
          
        
        N
      
      ⁢
      
        
      
      ⁢
      
        
          g
          j
          
            (
            k
            )
          
        
        ⁡
        
          (
          t
          )
        
      
    
  
  ;



Assuming the said weight value of the kth stage parallel interference cancellation is p(k), then calculating the said output signal of user i in step (f) by a formula





  
    
      
        r
        i
        
          (
          
            k
            +
            1
          
          )
        
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        r
        ⁡
        
          (
          t
          )
        
      
      -
      
        
          p
          
            (
            k
            )
          
        
        ⁢
        
          
            I
            ^
          
          i
          
            (
            k
            )
          
        
      
    
  
  ,



 where, r(t) is the said base band signal of the said received signal.

3. The double weighting parallel interference cancellation method according to claim 2, wherein step (c) of calculating the said reliability coefficient of the said decision result on the mth symbol of user i comprises, replacing the said Hyperbolic tangent decision with a piecewise linear decision, which means approaching a Hyperbolic tangent function tanh(x) with a piecewise linear decision function L(x), a expression of the piecewise linear decision function L(x) is:



  
    1
    ,
    
      x
      >
      θ
    
  




  
    
      L
      ⁡
      
        (
        x
        )
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    
                      
                        
                          
                            L
                            q
                          
                          ⁡
                          
                            (
                            x
                            )
                          
                        
                        ,
                        
                          x
                          ∈
                          
                            [
                            
                              
                                x
                                
                                  q
                                  -
                                  1
                                
                              
                              ,
                              
                                x
                                q
                              
                            
                            ]
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          -
                          
                            
                              L
                              q
                            
                            ⁡
                            
                              (
                              
                                -
                                x
                              
                              )
                            
                          
                        
                        ,
                        
                          x
                          ∈
                          
                            [
                            
                              
                                -
                                
                                  x
                                  q
                                
                              
                              ,
                              
                                -
                                
                                  x
                                  
                                    q
                                    -
                                    1
                                  
                                
                              
                            
                            ]
                          
                        
                      
                    
                  
                
              
            
          
          .
          
            

          
          ⁢
          
            -
            1
          
        
        ,
        
          x
          <
          
            -
            θ

4. The double weighting parallel interference cancellation method according to claim 2, wherein step (c) of calculating the said reliability coefficient of the said decision result on the mth symbol of user i comprises, replacing the said Hyperbolic tangent decision with a look-up table method, which means approaching a Hyperbolic tangent function tanh(x) with a decision function T(x) of the look-up table method, a expression of the said decision function T(x) of the look-up table method is:




  1
  ,
  
    
      x
      >
      
        θ
        ⁢
        
          

        
        ⁢
        
          T
          ⁡
          
            (
            x
            )
          
        
      
    
    =
    
      {
      
        
          
            
              
                
                  
                    
                      
                        
                          tanh
                          ⁡
                          
                            (
                            
                              
                                
                                  x
                                  
                                    q
                                    -
                                    1
                                  
                                
                                +
                                
                                  x
                                  q
                                
                              
                              2
                            
                            )
                          
                        
                        ,
                        
                          x
                          ∈
                          
                            [
                            
                              
                                x
                                
                                  q
                                  -
                                  1
                                
                              
                              ,
                              
                                x
                                q
                              
                            
                            ]
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          -
                          
                            tanh
                            ⁡
                            
                              (
                              
                                
                                  
                                    x
                                    
                                      q
                                      -
                                      1
                                    
                                  
                                  +
                                  
                                    x
                                    q
                                  
                                
                                2
                              
                              )
                            
                          
                        
                        ,
                        
                          x
                          ∈
                          
                            [
                            
                              
                                -
                                
                                  x
                                  q
                                
                              
                              ,
                              
                                -
                                
                                  x
                                  
                                    q
                                    -
                                    1
                                  
                                
                              
                            
                            ]
                          
                        
                      
                    
                  
                
              
            
          
          .
          
            

          
          ⁢
          
            -
            1
          
        
        ,
        
          x
          <
          
            -
            θ

5. The double weighting parallel interference cancellation method according to claim 1, wherein step (e) comprises, setting the different weight value p(k) for the interference cancellation in different stages, wherein p(1)<p(2) . . . <p(k) . . . <p(S), where, k is the kth stage of the interference cancellation, and S is a stage number of the interference cancellation.

6. A double weighting parallel interference cancellation method that can be used in CDMA mobile communication system, comprises:
(a) Making a multi-path despreading, channel estimation and a multi-path combining on an input signal of a user by a RAKE receiver;
(b) Making a soft decision on a result of the multi-path combining and a value of the channel estimation;
(c) Setting a weight value for the soft decision in different stages, and weighting the soft decision in symbol level;
(d) Regenerating a weighted signal of the user in chip level by the weighted soft decision result and the channel estimation value;
(e) Estimating multiple access interference on each user by the weighted regenerated signal in chip level inputted in parallel of each user;
(f) Subtracting the multiple access interference on an expected user from a base band signal of a received signal, and obtaining an output signal of the expected user, which is also an input signal of the said expected user in the next stage of parallel interference cancellation.

7. The double weighting parallel interference cancellation method according to claim 6, comprises, if the said input signal of user i in the kth stage parallel interference cancellation is expressed as




  
    
      r
      i
      
        (
        k
        )
      
    
    ⁡
    
      (
      t
      )
    
  
  ,



In step (a), calculating the said multi-path combining result of user i by a formula





  
    
      Y
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      Re
      ⁢
      
        {
        
          y
          i
          
            
              (
              m
              )
            
            ⁢
            
              (
              k
              )
            
          
        
        }
      
    
  
  ,



 which can be expressed as





  
    
      Y
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        
          μ
          i
        
        ⁢
        
          a
          i
          
            (
            m
            )
          
        
      
      +
      
        n
        i
      
    
  
  ,



 where, ni is a Gaussian White Noise, which is subject to a normal distribution of





  
    N
    ⁡
    
      (
      
        0
        ,
        
          σ
          i
          2
        
      
      )
    
  
  ,
  
    σ
    2
  



 is a noise power of





  
    n
    i
  
  ⁢
  
    a
    i
    
      (
      m
      )
    
  



 is the mth symbol of user i, whose value is +1 or −1, μi is a real number that is related with channel fading;
In step (b), expressing the said soft decision of the said multi-path combining result for user i as





  
    
      ξ
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        f
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
      ⁢
      
        
          a
          ^
        
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
    
  
  ,



 calculating the said soft decision by a formula





  
    
      
        f
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
      ⁢
      
        
          a
          ^
        
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
    
    =
    
      tanh
      ⁢
      
        {
        
          w
          ⁢
          
            
              
                μ
                i
              
              ⁢
              
                Y
                i
                
                  
                    (
                    m
                    )
                  
                  ⁢
                  
                    (
                    k
                    )
                  
                
              
            
            
              σ
              i
              2
            
          
        
        }
      
    
  
  ,



 where, W is a positive real number,





  
    a
    ^
  
  i
  
    
      (
      m
      )
    
    ⁢
    
      (
      k
      )
    
  



 is a soft decision result for the mth symbol of user





  i
  ,
  
    f
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  



 is a reliability coefficient of





  
    
      a
      ^
    
    i
    
      
        (
        m
        )
      
      ⁢
      
        (
        k
        )
      
    
  
  ;



Assuming the said weight value of the said soft decision in the kth stage is p(k), and then in step (c) calculating the said weighted soft decision of user i in symbol level by a formula





  
    
      ρ
      i
      
        
          (
          m
          )
        
        ⁢
        
          (
          k
          )
        
      
    
    =
    
      
        ζ
        i
        
          
            (
            m
            )
          
          ⁢
          
            (
            k
            )
          
        
      
      ·
      
        p
        
          (
          k
          )
        
      
    
  
  ;



In step (d), calculating the said weighed regenerated signal in chip level of user i by a formula





  
    
      
        g
        i
        
          (
          k
          )
        
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        ∑
        
          l
          =
          1
        
        L
      
      ⁢
      
        
      
      ⁢
      
        
          A
          il
        
        ⁢
        
          
            ∑
            
              n
              =
              
                -
                ∞
              
            
            ∞
          
          ⁢
          
            
          
          ⁢
          
            
              ρ
              i
              
                
                  (
                  n
                  )
                
                ⁢
                
                  (
                  k
                  )
                
              
            
            ⁢
            
              
                a
                ^
              
              i
              
                
                  (
                  n
                  )
                
                ⁢
                
                  (
                  k
                  )
                
              
            
            ⁢
            
              p
              ⁡
              
                (
                
                  t
                  -
                  
                    nT
                    b
                  
                  -
                  
                    τ
                    il
                  
                
                )
              
            
            ⁢
            
              
                c
                i
              
              ⁡
              
                (
                
                  t
                  -
                  
                    τ
                    il
                  
                
                )
              
            
          
        
      
    
  
  ,



 where, Ail is a estimation value of ail√{square root over (Pi)}, ail refers to a channel fading value of user i in the path l, and pi is a power of user i;
In step (e), calculating the said multiple access interference on user i by a formula





  
    
      
        I
        ^
      
      i
      
        (
        k
        )
      
    
    =
    
      
        ∑
        
          
            j
            =
            1
          
          ,
          
            j
            ≠
            i
          
        
        N
      
      ⁢
      
        
      
      ⁢
      
        
          g
          j
          
            (
            k
            )
          
        
        ⁡
        
          (
          t
          )
        
      
    
  
  ;



In step (f), calculating the said output signal of user i by a formula





  
    
      
        r
        i
        
          (
          
            k
            +
            1
          
          )
        
      
      ⁡
      
        (
        t
        )
      
    
    =
    
      
        r
        ⁡
        
          (
          t
          )
        
      
      -
      
        
          I
          ^
        
        i
        
          (
          k
          )
        
      
    
  
  ,



 where, r(t) is the said base band signal of the said received signal.

8. The double weighting parallel interference cancellation method according to claim 7, wherein step (b) of making the said soft decision of the said multi-path combining result for user i comprises, replacing the said Hyperbolic tangent decision with a piecewise linear decision, which means approaching a Hyperbolic tangent function tanh(x) with a piecewise linear decision function L(x), a expression of the piecewise linear decision function L(x) is:




  
    L
    ⁡
    
      (
      x
      )
    
  
  =
  
    {
    
      
        
          
            
              1
              ,
              
                x
                >
                θ
              
            
          
        
        
          
            
              
                
                  L
                  q
                
                ⁡
                
                  (
                  x
                  )
                
              
              ,
              
                x
                ∈
                
                  [
                  
                    
                      x
                      
                        q
                        -
                        1
                      
                    
                    ,
                    
                      x
                      q
                    
                  
                  ]
                
              
            
          
        
        
          
            
              
                -
                
                  
                    L
                    q
                  
                  ⁡
                  
                    (
                    
                      -
                      x
                    
                    )
                  
                
              
              ,
              
                x
                ∈
                
                  [
                  
                    
                      -
                      
                        x
                        q
                      
                    
                    ,
                    
                      -
                      
                        x
                        
                          q
                          -
                          1
                        
                      
                    
                  
                  ]
                
              
            
          
        
        
          
            
              
                -
                1
              
              ,
              
                x
                <
                
                  -
                  θ
                
              
            
          
        
      
      .

9. The double weighting parallel interference cancellation method according to claim 7, wherein step (b) of making the said soft decision of the said multi-path combining result for user i comprises, replacing the said Hyperbolic tangent decision with a look-up table method, which means approaching a Hyperbolic tangent function tanh(x) with a decision function T(x) of the look-up table method, a expression of the said decision function T(x) of the said look-up table method is:




  
    T
    ⁡
    
      (
      x
      )
    
  
  =
  
    {
    
      
        
          
            
              1
              ,
              
                x
                >
                θ
              
            
          
        
        
          
            
              
                tanh
                ⁡
                
                  (
                  
                    
                      
                        x
                        
                          q
                          -
                          1
                        
                      
                      +
                      
                        x
                        q
                      
                    
                    2
                  
                  )
                
              
              ,
              
                x
                ∈
                
                  [
                  
                    
                      x
                      
                        q
                        -
                        1
                      
                    
                    ,
                    
                      x
                      q
                    
                  
                  ]
                
              
            
          
        
        
          
            
              
                -
                
                  tanh
                  ⁡
                  
                    (
                    
                      
                        
                          x
                          
                            q
                            -
                            1
                          
                        
                        +
                        
                          x
                          q
                        
                      
                      2
                    
                    )
                  
                
              
              ,
              
                x
                ∈
                
                  [
                  
                    
                      -
                      
                        x
                        q
                      
                    
                    ,
                    
                      x
                      
                        q
                        -
                        1
                      
                    
                  
                  ]
                
              
            
          
        
        
          
            
              
                -
                1
              
              ,
              
                x
                <
                θ
              
            
          
        
      
      .

10. The double weighting parallel interference cancellation method according to claim 6, wherein step (c) comprises, setting the different weight value p(k) for the soft decision in different stages, wherein p(1)<p(2) . . . <p(k) . . . <p(S), where, k is the kth stage of the interference cancellation, and S is a stage number of the interference cancellation.