CN1357181A - Method and system for fast maximum posteriori decoding - Google Patents

Method and system for fast maximum posteriori decoding Download PDF

Info

Publication number
CN1357181A
CN1357181A CN99814657.9A CN99814657A CN1357181A CN 1357181 A CN1357181 A CN 1357181A CN 99814657 A CN99814657 A CN 99814657A CN 1357181 A CN1357181 A CN 1357181A
Authority
CN
China
Prior art keywords
measured
branch transition
likelihood ratio
facing state
status quantitavie
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN99814657.9A
Other languages
Chinese (zh)
Other versions
CN1143471C (en
Inventor
J·莱策尔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Clastres LLC
Telefonaktiebolaget LM Ericsson AB
Original Assignee
Telefonaktiebolaget LM Ericsson AB
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US09/252,028 external-priority patent/US6343368B1/en
Application filed by Telefonaktiebolaget LM Ericsson AB filed Critical Telefonaktiebolaget LM Ericsson AB
Publication of CN1357181A publication Critical patent/CN1357181A/en
Application granted granted Critical
Publication of CN1143471C publication Critical patent/CN1143471C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/004Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045Arrangements at the receiver end
    • H04L1/0054Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/37Decoding methods or techniques, not specific to the particular type of coding provided for in groups H03M13/03 - H03M13/35
    • H03M13/39Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes
    • H03M13/3905Maximum a posteriori probability [MAP] decoding or approximations thereof based on trellis or lattice decoding, e.g. forward-backward algorithm, log-MAP decoding, max-log-MAP decoding
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/004Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045Arrangements at the receiver end
    • H04L1/0055MAP-decoding

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Theoretical Computer Science (AREA)
  • Error Detection And Correction (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Detection And Correction Of Errors (AREA)
  • Detection And Prevention Of Errors In Transmission (AREA)
  • Reduction Or Emphasis Of Bandwidth Of Signals (AREA)

Abstract

Processing delay associated with MAP decoding (and variants thereof) is reduced by processing forward and reverse state metrics in parallel. Further reductions in processing complexity can be achieved by omitting redundant operations and multiplexing processing units functionality.

Description

The method and system of fast maximum posteriori decoding
Background
The present invention relates generally to the fault processing in the field of wireless communications, more specifically to the signal that uses the decoding of maximum a posteriori probability (MAP) decoding technique with the error correcting code transmission.
The explosive increase of the development of business communication system, particularly cellular radiotelephone system forces the system designer seeks increases power system capacity under situation about communication quality not being reduced to outside the tolerable thresholding of user method.A kind of technology that reaches these targets relates to wherein using analog-modulated that the system transition that data are added on the carrier wave is added to system on the carrier wave for wherein using digital modulation with data.
In radio digital communication system, the standardization air interfaces specify most systems parameter, comprise modulation type, burst format, communication protocol etc.For example ETSI European Telecommunications Standards Institute (BTSI) has stipulated a kind of global system for mobile communications (GSM) standard, it use that Guassian Minimum Shift Keying GMSK (GMSK) modulation scheme character rate uses that time division multiple access inserts that (TDMA) controls on as radio frequency (RF) physical channel of 271ksps or link, the communicating by letter of speech and data message.In the U.S., telecommunications industry association (TIA) has published many interim standards, as IS-54 and IS-136, defined the digital AMPSAdvanced Mobile Phone Service (D-AMPS) of various version, a kind of use difference quadrature phase shift keying (DQPSK) modulation scheme is at the tdma system of RF link.
Tdma system is subdivided into one or more RF channels with available band, and each channel further is divided into the physical channel of a plurality of corresponding tdma frame time slots again.Logic channel has specified the physical channel of modulation and coding to form by one or several.In these systems, travelling carriage is by sending and receive the base station communication of digital information pulse train and a plurality of dispersions on up link and down link RF channel.
At present the growth number of the travelling carriage that uses has produced the demand of multi-voice and data channel more in the cellular telecommunication system.The base station becomes at interval nearer as a result, and making increases in interference adjacent or that operate in the nearer sub-district at interval between the travelling carriage on the same frequency.Although digital technology provides the how useful channel on the given frequency spectrum, also need to reduce and disturb, or more particularly will increase the ratio (being carrier wave-interference (C/I) ratio) of carrier signal strength and interference.
For various communication services are provided, need user's bit rate of corresponding minimum.For example to speech and/or data service, corresponding speech quality of user's bit rate and/or data throughout, user's bit rate is high more, and speech quality is good more and/or data throughout is big more.Total user's bit rate is determined by the combination technique of the speech coding of selecting, error correction coding (chnnel coding), modulation scheme is connected assignable resource with each number (as coding, time slot, frequency etc.).
Concrete care is error correction or chnnel coding in this standard.Digital communication system adopts various technology to handle the wrong information that receives.In general, these technology comprise the technology of the information that those help receivers to correct a mistake to receive, and can make as forward error correction (FEC) technology and those to receive the technology that wrong information resends receiver, as automatic repeat requests technology.The FEC technology comprises for example data convolution or the block encoding before modulation, and these codings are used to allow to correct by noise and disturb the error code that causes.The FEC coding relates to the data bit that the coded-bit that uses greater number is represented some.The common like this encoding rate by them is indicated convolution code as 1/2 and 1/3, and wherein lower code check provides stronger error protection, but then user's bit rate is lower for given channel bit rate.
At receiver side, the sequence that receives is decoded so that allow further deal with data.Certainly the channel of transmission signals makes this decoding processing have more challenge thereon.This challenging diagram is seen Fig. 1.Here, symbols streams u is input to convolution coder 10.Encoded symbols streams x goes up transmission and is received by the form of receiving entity with symbols streams y at a transmission channel 12 (for example an air interface or Wireline), and this symbols streams has reflected the influence of channel to transmission signals.MAP decoder 14 is interpreted as output symbol stream  with the symbols streams y that receives, and this is the estimation of original symbol stream u.Attention is discussed in order to simplify this, has ignored the description of various other processing well known in the art, and for example the modulation and demodulation of carrying out respectively at transmitter and receiver side are handled.
Decoder 14 can use the MAP decoding algorithm, this algorithm is also referred to as by the maximal posterior probability algorithm of symbol or bcjr algorithm, it merges the disclosure content here as a reference by Bahl, Cocke, Jelinek and Raviv proposition in March, 1974, IEEE Trans.Inform.Theory volume IT-20 284-287 page or leaf " minimized the optimal decoding of the uniform enconding of symbol error rate ".Also proposed the various mutation of MAP algorithm, these mutation have more commercial flexibility from the viewpoint of using.For example, merge the disclosure content here as a reference at the Log-Map of in March, 1996 people such as Hagenauer in the works " iterative decoding of binary block and convolution code " of IEEE Trans.Inform.Theory volume IT-42 429-445 page or leaf (carry out handle at log-domain) and Max-Log-Map (simple version of Log-Map).
Because with respect to not too complicated decoder (as the Viterbi decoder), the calculation of complex of MAP decoder and only have quite little performance to increase is so adopted the less relatively of MAP decoder from it since the seventies occurs under most conditions.But along with the continuous development of error correction coding, to the interest of MAP decoder owing to rekindled near the agricultural limit earlier.The professional in present technique field will understand energy per bit and the monolateral noise density (E that the first agricultural limit has determined to be used for the minimum of reliable transmission b/ N o) ratio.Although farming has proved this limit at the forties latter stage earlier, also do not develop the error correction coding that this theoretical possibility performance can be provided.But occurred a class what is called " turbo (turbine) " sign indicating number recently, they more approach the operation on the formerly agricultural limit.Because the arrival of these and other novel error correction codings, the MAP decoder has been got back to leading position.
The detailed computing of MAP algorithm is described below, and the formula of reviewing here can be understood measuring of relating in the MAP algorithm better.It is that a decoder of realizing the MAP algorithm in direct mode simply is provided that of MAP decoder may use.That is to say that to the sequence with N symbol that receives, the first step is calculated and stored all branch transition and measures with forward-facing state and measure.Calculated N branch transition measure measure with N forward-facing state after, calculate all N back to status quantitavie.At last, calculate the back behind status quantitavie, calculated likelihood ratio.Even but this direct MAP decoding processing realizes that related processing delay is also quite high when using powerful processor to carry out these calculating.Present communication needs more and more higher bit rate, promptly be need be more and more littler the processing delay tolerance limit.
Reducing of finding in the literature finds in the WO98/20617 disclosure of exercise question for " the soft-decision output decoder of the code word of decoding convolutional encoding " by being intended that of the processing delay related with using the MAP algorithm, and merging is as a reference especially here for its disclosure.Here, described the Log-Map decoder, wherein the viterbi decoder of first " general " is at initial condition t 0Beginning also provides a plurality of forward-facing states to measure α for each state in each time interval on length of window 2L.The viterbi decoder of second " general " also is provided, and it is at the second time t 2LBeginning also provides a plurality of backs to status quantitavie β for each state in each time interval.Then processor is carried out each state and is asked max calculation twice, uses forward-facing state to measure to first state, to second state use the back to status quantitavie and between to first and second states path use the branch transition state.
Although this MAP realization of decoding realizes reducing processing delay with respect to direct MAP algorithm, its complexity is still higher relatively and say so inferior good after second " general " viterbi decoder uses initial estimation to calculate when beginning at every turn on the meaning of status quantitavie.Therefore be desirable to provide the added complexity with the processing delay that reduces, limited (if there is) and do not need to take the MAP decoder of the inferior good mutation of algorithm to realize making it possible to strong error correction coding is carried out efficient coding, described strong error correction coding is considered for next generation communication system.
Summary
These and other shortcomings and the restriction of the method and system of traditional transmission information can be overcome according to the present invention, processing delay (as with respect to direct applied one approximate 2 factor) can be reduced to decoder and the decoding technique that status quantitavie calculates (RSMC) parallel processing by using a kind of wherein forward-facing state measure calculation (FSMC) and back.Implementation complexity also reduces with respect to for example WO98/20617 and can further reduce time delay by approval processing sequential and the elimination unnecessary extra process unit (as branch transition metric computation unit and likelihood ratio computing unit) related with various tasks.
The accompanying drawing summary
Read following detailed in conjunction with the accompanying drawings, these and other purposes of the present invention, feature and advantage will become more obvious, wherein:
Fig. 1 is that explanation is relating to the general block diagram that sends and receive a signal on the channel that uses error correction coding;
Fig. 2 is the block diagram according to the demonstration decoder of exemplary embodiment of the present invention;
Fig. 3 is the flow chart of describing according to the demonstration decode procedure of exemplary embodiment of the present invention;
Fig. 4 is the sequential chart of timing aspect of operation of the demonstration program of the demonstration decoder of key diagram 2 and Fig. 3;
Fig. 5 is another sequential chart of timing aspect of operation of the demonstration program of the demonstration decoder of key diagram 2 and Fig. 3;
Fig. 6 is the block diagram according to the decoder of another exemplary embodiment of the present invention; And
Fig. 7 is the block diagram that can use exemplary wireless communications systems of the present invention.
Describe in detail
The unrestricted purpose for explanation has in the following description been delivered concrete details such as concrete circuit, circuit element, technology etc. so that thoroughly understand the present invention.But to the professional in present technique field obviously the present invention can realize with other embodiments that are different from detail.In other example, will ignore the detailed description of known method, equipment and circuit for not fuzzy description of the invention.
Should provide some wherein to consider environment more of the present invention so that be given in order to begin this discussion about the additional detail of MAP decoding algorithm and demonstration mutation such as LogMAP algorithm and Max-Log-MAP algorithm.Use some variablees in the following discussion, be defined as follows:
S k: the state of k node in the grating texture
Measure calculation relates to probability, and its symbol is as follows: A, B and C presentation of events, then to give a definition effectively:
P (A): the probability of incident A
(A, B): incident A and B unite the probability of generation to P
(C): incident A, B and C unite the probability of generation to P for A, B
P (A|B): the conditional probability of incident A when given incident B takes place
(A, B|C): incident A and B united the conditional probability of generation to P when given incident C took place
(A|B, C): given incident B and C unite the conditional probability MAP algorithm that incident A takes place when taking place to P
The MAP algorithm comprises following 4 steps:
1. branch transition measure calculation (BTMC): from the symbol y that receives at channel 12 kBranch Computed shifts measures γ k(S K-1, S k):
γ k(S k-1,S k)=P(y k,S k|S k-1)
(1)
=P(y k|S k-1,S k)·P(S k|S k-1)
2. forward-facing state measure calculation (FSMC): initially to measure α 0(S 0) the recursive calculation forward-facing state measures α k(S k), initially measure initial condition definition according to convolution coder: α k ( S k ) = Σ S k - 1 α k - 1 ( S k - 1 ) · Y k ( S k - 1 , S k ) - - - ( 2 )
3. (RSMC) calculated to status quantitavie in the back: initially to measure β N(S N) after the recursive calculation to status quantitavie β k(S k), initially measure end-state definition according to convolution coder: β k - 1 ( S k - 1 ) = Σ S k β k ( S k ) · Y k ( S k - 1 , S k ) - - - ( 3 )
4. likelihood ratio is calculated (LRC): calculate soft output ∧ k, the likelihood ratio u of itself and given y kOf equal value.Note molecule and denominator with respectively at u k=+1 and u kThose states of=-1 are to last: Λ k = P ( u k = + 1 | y ) P ( u k = - 1 | y ) = Σ S k - 1 , S k , u k = + 1 P ( S k - 1 , S k , y ) Σ S k - 1 , S k , u k = - 1 P ( S k - 1 , S k , y ) (4) = Σ S k - 1 , S k , u k = + 1 α k - 1 ( S k - 1 ) · Y k ( S k - 1 , S k ) · β k ( S k ) Σ S k - 1 , S k , u k = - 1 α k - 1 ( S k - 1 ) · Y k ( S k - 1 , S k ) · β k ( S k ) The Log-MAP algorithm
If take the logarithm and the symbol a=ln α of all parameters in the formula (1-4), b=ln β, c=ln γ and L=ln Δ, then the formula of BTMC, FSMC, RSMC and LRC becomes 4 steps in following formula (1 '-4 ') and the processing Log-MAP algorithm:
1.BTMC:
c k(S k-1,S k)=ln(P(y k|S k-1,S k))+ln(P(S k|S k-1))(1′)
2. initially measure and be a 0(S 0)=ln (α 0(S 0)) FSMC: a k ( S k ) = ln ( Σ S k - n e a k - 1 ( S k - 1 ) + c k ( S k - 1 , S k ) ) - - - ( 2 ′ )
3. initially measure and be b 0(S 0)=ln (β 0(S 0)) RSMC: b k - 1 ( S k - 1 ) = ln ( Σ S k e b k ( S k ) + c k ( S k - 1 , S k ) ) - - - ( 3 ′ )
4.LRC: L k = ln ( Σ S k - 1 , S k , u k = + 1 e a k - 1 ( S k - 1 ) + c k ( S k - 1 , S k ) + b k ( S k ) ) -
(4′) - ln ( Σ S k - 1 , S k , u k = - 1 e a k - 1 ( S k - 1 ) + c k ( S k - 1 , S k ) + b k ( S k ) ) The Max-Log-MAP algorithm
If use approximation: ln ( Σ i e x i ) = Max i ( x i ) - - - ( 5 )
That formula (1 '-4 ') becomes is following (1 "-4 "), the Log-MAP algorithm becomes the Max-Log-MAP algorithm. and four steps of then carrying out the Max-Log-MAP algorithm are as follows:
1.BTMC:
c k(S k-1,S k)=ln(P(y k|S k-1,S k))+ln(P(S k|S k-1))(1″)
2. initially measure and be a 0(S 0)=ln (α 0(S 0)) FSMC: a k ( S k ) = Max S k - 1 ( a k - 1 ( S k - 1 ) + c k ( S k - 1 , S k ) ) - - - ( 2 ′ ′ )
3. initially measure and be b 0(S 0)=1n (β 0(S 0)) RSMC: b k - 1 ( S k - 1 ) = Max S k ( b k ( S k ) + c k ( S k - 1 , S k ) ) - - - ( 3 ′ ′ )
4.LRC: L k = Max S k - 1 , S k , u k = + 1 ( a k - 1 ( S k - 1 ) + c k ( S k - 1 , S k ) + b k ( S k ) ) -
(4″) - Max S k - 1 , S k , u k = - 1 ( a k - 1 ( S k - 1 ) + c k ( S k - 1 , S k ) + b k ( S k ) )
The same with the many algorithms that find in the document, the MAP theory is interesting, but theoretical realization has proposed realistic problem to the design engineer.In MAP algorithm and above-mentioned mutation thereof, a such problem relates to processing delay.The demonstration of having described relative WO 98/20617 above realizes that this realization wants to reduce the processing delay related with the MAP algorithm.Exemplary embodiment of the present invention provides the another kind of realization that also reduces processing delay, and computation complexity is less and estimated value is not used to status quantitavie calculating in each back.
Figure 2 shows that the block diagram of an exemplary embodiment of the present invention.Because the present invention can be applicable to original MAP algorithm, and Log-MPA and Max-Log-MAP (with other any mutation) algorithm, signal name and reference symbol among the figure that the exemplary embodiment is here provided are described below.To the effective signal name of above-mentioned agreement of MAP algorithm on every signal line with top or with first signal with reference to giving out.After the MAP signal name or under square brackets in provide effective another signal name of the naming convention related with Log-MAP algorithm or Max-Log-MAP algorithm.For example referring to the upper left corner of Fig. 2, the holding wire related with input above BTMC unit 20 has a MAP signal with reference to P (y k| S K-1, S k) and a Log-MAP (or Max-Log-MAP) signal with reference to ln (P (y k| S K-1, S k)).
The demonstration decoder of Fig. 2 has the different functional block of a plurality of following work according to the present invention.BTM is that the mode that the concrete version of the MAP algorithm that uses with decoder in BTMC unit 20 and 22 is determined is calculated.If for example adopt the MAP algorithm, then branch metric calculates according to formula (1), if adopt the Log-MAP algorithm, then according to formula (1 ') if or adopt the Max-Log-MAP algorithm, then calculate according to formula (1 ").
Similarly, calculating forward-facing state in the mode that depends on the MAP mutation that decoder uses in FSMC unit 24 measures.Just, if adopt the MAP algorithm to calculate FSM, if adopt the Log-MAP algorithm to calculate, if use the Max-Log-MAP algorithm to calculate according to formula (2 ") according to formula (2 ') according to formula (2).In RSMC unit 26, calculate the back similarly to status quantitavie, if adopt the MAP algorithm, if adopt the Log-MAP algorithm to calculate, if use the Max-Log-MAP algorithm to calculate according to formula (3 ") according to formula (3 ') according to formula (3).
In LRC unit 28 and 30, use branch metric, the forward-facing state in unit 20,22,24 and 26, determined to measure and reverse state measure calculation likelihood ratio.In particular,, then calculate likelihood ratio,, then calculate,, then calculate according to formula (4 ") if adopt the Max-Log-MAP algorithm according to formula (4 ') if adopt the Log-MAP algorithm according to formula (4) if adopt the MAP algorithm.
Fig. 2 also comprises time delay unit 32 and 34 and memory cell 36 and 38.These unit are used to cushion and control the measurement value of delivering to LRC unit 28 and 30 by timing and control unit 40.For avoiding fuzzy this figure, ignored from the holding wire of other functional blocks of timing and control unit 40 to Fig. 2.But the professional in present technique field will understand the timing (I/O) of each square frame in timing and control unit 40 control charts 2, can calculate and transmit measurement value and regularly and between each functional block of control unit 40 and Fig. 2 provide holding wire in the mode that describes below of for example Fig. 3-5 like this.
In order to understand the working method according to demonstration MAP decoder of the present invention better, Fig. 3 provides the flow chart of an explanation demonstration work decoding.Here receive a sequence deciphering in step 50.This sequence that receives is passable, for example, handles (as down-conversion, demodulation etc.) earlier with the variety of way that the professional in present technique field knows better at receiver.Next in step 52, the initial condition that forward-facing state is measured and the back will be used to the status quantitavie recurrence of unit 24 and 26 is set respectively.Determine the initial condition of FSMC unit from the initial condition (fix typically, so receiver being known) of convolution coder.If end-state is known,, then determine the initial condition of RSMC unit from the end-state of convolution coder as when convolution coder uses the convolution process that stops.Otherwise for example when convolution coder uses hangover convolution process, can use the estimation of end-state.
For example, consider a system, wherein the convolution coder of transmitter has i.e. 4 kinds of different possible the states of 3 constraint length, and wherein receiver adopts the Log-MAP algorithm of the realization according to the present invention.If the initial condition of encoder is for example S 0=0, then initial FSM can be set to a 0(S 0=0)=0 and a 0(S 0=1)=a 0(S 0=2)=a 0(S 0=3)=-∞.If end-state is known as for example S N=2, then Initial R SM can be set to b N(S N=2)=0 and b N(S N=0)=b N(S N=1)=b 0(S 0=3)=-∞.Receiver is not known the end-state of convolution coder else if, and then Initial R SM can be set to b N(S N=0)=b N(S N=1)=b N(S N=2)=b N(S N=3)=0.
Get back to Fig. 3 now, as a reference affix k carries out initialization in step 54 pair decode procedure.This affix can for example corresponding time step based on character rate and/or sampling rate selection.Carry out initialization with initial FSM and RSM respectively in this moment FSMC unit 24 and RSMC 26 in step 52 storage.In step 56, increase affix k and decode procedure and begin.The demonstration decode procedure of describing in this embodiment of the present invention is by a periodic duty, and wherein forward-facing state is measured and be back to the status quantitavie parallel computation, has reduced decoding delay like this.Processing method is odd number or even number and slightly different according to symbolic number N in the sequence that receives.
First in circulation promptly comprises step 56-62, and the branch transition of calculating and storing receiving sequence is measured, forward-facing state is measured and a back part to status quantitavie.In particular, (if N is an even number) or k=N/2-1/2 (if N is an odd number) from step-length K=1 to step-length k=N/2 calculate affix and are the RSM (step 60) that the BTM of k and N-k+1 (step 58), FSM (step 60) that affix is k and affix are N-k.FSM and RSM are stored in respectively in unit 36 and 38.In this, as determining that in declaring piece 62 decode procedure has had enough metric data and begun to calculate likelihood ratio, and remaining is measured to be calculated so that process expanded to and comprises that likelihood ratio calculates.
If odd number symbol N is arranged, then only calculate one time likelihood ratio during the iteration in the first time of Extended Cyclic in receiving sequence.In step 64,, use α here like this if k=N/2+1/2 (promptly having only N is that odd number is just for true) handling process forwards step 66 to K-1(S K-1), γ k(S K-1, S k) and B k(S k) to calculate affix be the LR of N/2+1/2.The flow process circulation turns back to step 56 then increases affix k, calculates and stores one group and newly measure and calculate then extra likelihood ratio.
If on the other hand N be even number or for the first time after the iteration when N is odd number, flow process forwards piece 68 to from piece 64, wherein through two likelihood ratios of the each calculating of circulation.Exactly,, calculate and the storage affix is the BTM of k and N-k+1 up to step-length k=N from step-length k=N/2+1 (if N is an even number) or k=N/2+3/2 (if N is an odd number), affix is that the FSM of k, RSM and the affix that affix is N-k are the LR of K and N-k+1.When decoder arrived step-length k=N, process 70 was come out along "No" branch from Decision Block.Then can be from the sequence after the decoder output decoding in step 72.
As previously mentioned, timing and control unit 40 (Fig. 2) are controlled calculating of measuring and the demonstration decode procedure that the transmission value is described with reference to figure 3 with acquisition between various computing units.Because forward-facing state is measured and the back be that first sign-on from receiving sequence calculates then to status quantitavie is that last sign-on from receiving sequence calculates to the status quantitavie parallel processing and because forward-facing state is measured, therefore for some time at interval, therebetween receiving sequence is calculated and measured, but do not calculate likelihood ratio.As mentioned above, in case affix arrives N/2, just can begin to calculate likelihood ratio.The sequential of exemplary embodiment of the present invention can use the sequential chart of Figure 4 and 5 to see clearlyer, wherein with top Fig. 2 and the same naming convention of 3 uses.
Fig. 4 explanation is when N is even number and measure timing with the likelihood ratio compute associations.Wherein can see at moment k=0, initial forward status quantitavie of FSMC 24 outputs, and RSMC26 therewith exports an initial back to status quantitavie simultaneously.But LRC 28 and 30 does not provide likelihood ratio constantly at this.When affix k increased, FSMC 24 and RSMC 26 continuation output forward directions and back were to measuring.When k=N/2+1, LRC unit 28 and 30 all uses previous forward direction that calculates and store and back to begin to export likelihood ratio to status quantitavie.This will proceed to k=N always.
Fig. 5 explanation is when N is odd number and measure timing with the likelihood ratio compute associations.Wherein can see, again, at moment k=0, initial forward status quantitavie of FSMC 24 outputs, and RSMC 26 therewith exports an initial back to status quantitavie simultaneously.But LRC 28 and 30 does not provide likelihood ratio constantly at this.When affix k increased, FSMC 24 and RSMC 26 continuation output forward directions and back were to measuring.When k=N/2+1, LRC unit 30 first likelihood ratios of output.When k=N/2+3/2, LRC unit 28 and 30 all uses previous forward direction that calculates and store and back to begin to export likelihood ratio to status quantitavie then.This will proceed to k=N always.
Calculating in the likelihood ratio according to aforementioned formula, can see being used for determining that forward-facing state is measured and repeat the back during the calculating section of status quantitavie is calculating likelihood ratio.Exactly, forward-facing state measure with the branch status quantitavie multiply each other (or addition) and the back is common computing (as relatively formula (2) and formula (4)) to multiply each other (or addition) of status quantitavie and branch's status quantitavie.According to another exemplary embodiment of the present invention shown in Figure 6, these common computings only need be carried out once like this, promptly during measure calculation rather than in measure calculation and likelihood ratio computing interval.In this drawing, common reference number is reused in the unit that the decoder with Fig. 2 is had same input and output, and the signal naming convention also is the same.So only FSMC 80 has different reference numbers and different output with RSMC 82.Be further noted that and in the decoder of Fig. 6, ignored time delay unit 32 and 34 that this is because LRC unit 28 and 30 has directly used FSM and RSM respectively.
The professional that other of these exemplary embodiment change the present technique field also is obvious.If for example BTMC unit 20 and 22 processing delay are less than or equal to half of processing delay of FSMC unit and RSMC unit respectively, then can from the decoder Fig. 2 and 6, remove a BTMC unit.For example can provide single BTMC unit to come alternately to provide a branch transition to measure to FSMC unit and RSMC unit.Can provide a multiplexer (not shown) between FSMC unit and RSMC unit, to switch its output in the output of BTMC unit.Certainly, also need need need between the unit that the end Branch Computed of sequence grating texture is measured, switch from unit that the beginning Branch Computed of sequence grating texture is measured and those to the input of single BTMC at those.
Similarly, be less than or equal to half processing delay of the processing delay related, then also can remove one of them LRC unit with FSMC unit and RSMC unit if LRC unit 28 and 30 has respectively.In this embodiment, can provide a demodulation multiplexer (not shown) to come alternately input with the single LRC of being outputted to of FSMC and RSMC unit.According to another exemplary embodiment of the present invention, can use the channelization treatment technology to make up FSMC and RSMC unit further to reduce complexity according to decoder of the present invention.
Decoding technique according to the present invention may be used on any communication system and/or environment.But as mentioned above, these MAP decoding techniques (with its mutation) can find concrete application in wireless communications environment, transmission channel 12 is that an air interface and encoder 10 and decoder 14 are parts (or vice versa) of base station and travelling carriage in this environment.For some versatilities are provided, Fig. 7 and following description provide some can realize the general discussion of a kind of like this demonstration system of the present invention.
Fig. 7 represents a block diagram that comprises the demonstration cellular mobile radiotelephone system of an exemplary base station 110 and travelling carriage 120.The base station comprises control and processing unit 130, and it is connected to MSC 140, the latter and then be connected to the PSTN (not shown) again.The general features of this cellular radiotelephone system is being known in the art, as authorize the United States Patent (USP) 5th of people's such as Wejke exercise question for " the adjacent auxiliary switching in the cellular communication system ", 175, No. 867 and authorize the United States Patent (USP) 5th of people's such as Raith exercise question for " a kind of method that in wireless communication system, communicates ", 603, describe in 081, the disclosure of these two patents merges in this application as a reference.
A plurality of voice/datas (promptly professional) channel is handled by a Traffic Channel transceiver 150 in base station 110, and transceiver is by control and processing unit 130 controls.Also have each base station to comprise a control channel transceiver 160, it can handle more than one control channel.Control channel transceiver 160 is by control and processing unit 130 controls.Control channel transceiver 160 on the control channel of base station or sub-district to the travelling carriage broadcast control information that locks onto this control channel.Should be appreciated that transceiver 150 and 160 can realize with individual equipment, as professional and control transceiver 170, and transceiver needs not be special use to the distribution of antenna.As the part of the signal processing of wherein carrying out, base station 110 can comprise a decoder, as mentioned above, is used to remove the error correction coding of the signal correction connection that transmits with travelling carriage 120 on control channel or Traffic Channel.
Travelling carriage 120 is received in the information of broadcasting on the control channel in its business and control channel transceiver 170.Again, receive the part of information processing as this, travelling carriage 120 can comprise the decoder that above-mentioned relatively any previous exemplary embodiment is described.Then, the control channel information (this information comprises the feature of the candidate cell that travelling carriage will lock) that processing unit 180 assessments receive, and which sub-district definite travelling carriage should be locked in.In case distributed a Traffic Channel, the information that processing unit 180 is also deciphered and demodulation receives on Traffic Channel for travelling carriage 120.
Although only describe the present invention in detail with reference to several exemplary embodiment, the professional in the present technique field will understand not departing under the situation of the present invention can carry out various modifications.For example above-mentioned processing can be used to decipher the received signal sequence through turbine coding, wherein will calculate likelihood ratio for a sequence before the value of symbol that receives is carried out any hard decision.So as can be seen the present invention can be applied to equally wherein with hard symbol judgement in series, before carrying out hard symbol judgement or even do not carry out hard decision and produce the decoding technique of soft information.Therefore the present invention only wants that by following the claim that comprises all its equivalents is defined.

Claims (25)

1. the method for the symbol sebolic addressing that receives of a decoding comprises step:
Calculating first and second branch transition measures;
Measuring definite forward-facing state based on described first branch transition measures;
With determine step that described forward-facing state is measured concurrently, measure based on described second branch transition and to determine that the back is to status quantitavie;
Based on described branch transition measure, described forward-facing state is measured and described back to the combination calculation likelihood ratio of status quantitavie;
Use described likelihood ratio to decipher the symbol sebolic addressing of described reception.
2. the process of claim 1 wherein that described branch transition is measured, described forward-facing state is measured, described back to status quantitavie and described likelihood ratio according to the MAP algorithm computation.
3. the process of claim 1 wherein that described branch transition is measured, described forward-facing state is measured, described back to status quantitavie and described likelihood ratio according to the Log-MAP algorithm computation.
4. the process of claim 1 wherein that described branch transition is measured, described forward-facing state is measured, described back to status quantitavie and described likelihood ratio according to the Max-Log-MAP algorithm computation.
5. the process of claim 1 wherein that the step that described first and second branch transition of described calculating are measured also comprises step:
Described first and second branch transition of parallel computation are measured.
6. the process of claim 1 wherein that the step that described first and second branch transition of described calculating are measured also comprises step:
Alternately calculating described first and second branch transition measures.
7. the process of claim 1 wherein that the step of the described likelihood ratio of described calculating also comprises step:
Two likelihood ratios of parallel computation.
8. the process of claim 1 wherein that the step of the described likelihood ratio of described calculating also comprises step:
Alternately calculate:
Based on first likelihood ratio that described forward-facing state one of is measured, described back one of is measured to one of status quantitavie and described first branch transition; And
Based on second likelihood ratio that described forward-facing state one of is measured, described back one of is measured to one of status quantitavie and described second branch transition.
9. the process of claim 1 wherein and describedly determine that described forward direction and back use the channelization treatment technology to carry out to the parallel step of status quantitavie.
10. the method for claim 1 also comprises step:
Storing described forward-facing state measures with described back to status quantitavie; And
After measuring and be described, behind status quantitavie, calculates the described forward-facing state of having stored predetermined number described likelihood ratio.
11. the method for claim 10 also comprises step:
Storing described first and second branch transition measures.
12. the process of claim 1 wherein that described decoding step also comprises step:
The value of symbol of the symbol sebolic addressing that hard decision receives after having calculated all likelihood ratios of described sequence.
13. decipher N method that receives the sequence of symbol and comprise step for one kind:
(a) initialization affix k;
(b) increase described affix k;
(c) calculating first branch transition based on the sequence of described reception measures;
(d) calculating second branch transition based on the sequence of described reception measures;
(e) measuring definite forward-facing state based on described first branch transition measures;
(f) the described forward-facing state of storage is measured;
(g) with determine step that described forward-facing state is measured concurrently, measure based on described second branch transition and to determine that the back is to status quantitavie;
(h) the described back of storage is to status quantitavie;
(i) repeating step (b)-(h) is up to k≤N/2, measures with described back based on described forward-facing state then to begin to calculate likelihood ratio to status quantitavie; And
(j) use described likelihood ratio to decipher the symbol sebolic addressing of described reception.
14. the method for claim 13, wherein said branch transition is measured, described forward-facing state measure described back to status quantitavie and described likelihood ratio according to the MAP algorithm computation.
15. the method for claim 13, wherein said branch transition is measured, described forward-facing state is measured with described likelihood ratio according to the Log-MAP algorithm computation.
16. the method for claim 13, wherein said branch transition is measured, described forward-facing state measure described back to status quantitavie and described likelihood ratio according to the Max-Log-MAP algorithm computation.
17. the method for claim 13, the step that described first and second branch transition of wherein said calculating are measured also comprises step:
Described first and second branch transition of parallel computation are measured.
18. the method for claim 13, the step that described first and second branch transition of wherein said calculating are measured also comprises step:
Using same processing unit alternately to calculate described first and second branch transition measures.
19. the method for claim 13, the step of the described likelihood ratio of wherein said calculating also comprises step:
Two likelihood ratios of parallel computation.
20. the method for claim 13, the step of the described likelihood ratio of wherein said calculating also comprises step:
Alternately calculate:
Based on first likelihood ratio that described forward-facing state one of is measured, described back one of is measured to one of status quantitavie and described first branch transition; And
Based on second likelihood ratio that described forward-facing state one of is measured, described back one of is measured to one of status quantitavie and described second branch transition.
21. the method for claim 13, wherein said definite described forward direction and back use the channelization treatment technology to carry out to the parallel step of status quantitavie.
22. the method for claim 13, wherein said decoding step also comprises step:
The value of symbol of the symbol sebolic addressing that hard decision receives after having calculated all likelihood ratios of described sequence.
23. the transceiver of a received signal comprises:
The reception processing unit of the sequence that reception will be deciphered; With
A MAP error correcting deocder from back to status quantitavie that decipher described sequence by parallel computation forward direction and.
24. the transceiver of claim 23, wherein said MAP error correcting deocder use described forward direction of Log-MAP algorithm computation and back to status quantitavie.
25. the transceiver of claim 23, wherein said MAP error correcting deocder use described forward direction of Max-Log-MAP algorithm computation and back to status quantitavie.
CNB998146579A 1998-12-18 1999-12-13 Method and system for fast maximum posteriori decoding Expired - Fee Related CN1143471C (en)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
IB9802068 1998-12-18
WOPCT/IB98/02068 1998-12-18
US09/252,028 US6343368B1 (en) 1998-12-18 1999-02-18 Method and system for fast maximum a posteriori decoding
US09/252,028 1999-02-18
US09/252028 1999-02-18

Publications (2)

Publication Number Publication Date
CN1357181A true CN1357181A (en) 2002-07-03
CN1143471C CN1143471C (en) 2004-03-24

Family

ID=26318724

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB998146579A Expired - Fee Related CN1143471C (en) 1998-12-18 1999-12-13 Method and system for fast maximum posteriori decoding

Country Status (7)

Country Link
EP (1) EP1142183B1 (en)
JP (1) JP2002533991A (en)
CN (1) CN1143471C (en)
AT (1) ATE242945T1 (en)
AU (1) AU758505B2 (en)
DE (1) DE69908820T2 (en)
WO (1) WO2000038366A1 (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2360858B (en) * 2000-03-20 2004-08-18 Motorola Inc High-speed maximum a posteriori (MAP) architecture with optimized memory size and power consumption
JP3613134B2 (en) 2000-05-12 2005-01-26 日本電気株式会社 High speed turbo decoder
US7127664B2 (en) * 2000-09-18 2006-10-24 Lucent Technologies Inc. Reconfigurable architecture for decoding telecommunications signals
US6865710B2 (en) * 2000-09-18 2005-03-08 Lucent Technologies Inc. Butterfly processor for telecommunications
ATE323341T1 (en) * 2001-06-21 2006-04-15 Cit Alcatel LOG MAP DECODING
US7116732B2 (en) 2001-06-21 2006-10-03 Alcatel Method and apparatus for decoding a bit sequence
US6718504B1 (en) 2002-06-05 2004-04-06 Arc International Method and apparatus for implementing a data processor adapted for turbo decoding
US7055102B2 (en) * 2002-12-06 2006-05-30 Sandbridge Technologies, Inc. Turbo decoder using parallel processing
KR101282568B1 (en) 2006-08-08 2013-07-04 삼성전자주식회사 Method and apparatus for calculating log likeligddo ratio in communication system
JP6188106B2 (en) 2014-09-24 2017-08-30 株式会社日立国際電気 MIMO receiving apparatus and MIMO receiving method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0766736A (en) * 1993-08-25 1995-03-10 Sony Corp Viterbi decoding device
DE4437984A1 (en) * 1994-10-25 1996-08-14 Philips Patentverwaltung Transmission system with soft output decoding
JP2798123B2 (en) * 1995-11-17 1998-09-17 日本電気株式会社 Viterbi decoding device
KR100195745B1 (en) * 1996-08-23 1999-06-15 전주범 Add compare selecter of vitervi decoder

Also Published As

Publication number Publication date
EP1142183B1 (en) 2003-06-11
AU758505B2 (en) 2003-03-20
JP2002533991A (en) 2002-10-08
EP1142183A1 (en) 2001-10-10
DE69908820T2 (en) 2004-04-15
DE69908820D1 (en) 2003-07-17
CN1143471C (en) 2004-03-24
AU1978000A (en) 2000-07-12
WO2000038366A1 (en) 2000-06-29
WO2000038366A8 (en) 2001-12-20
ATE242945T1 (en) 2003-06-15

Similar Documents

Publication Publication Date Title
CN1178397C (en) Soft decision output decoder for decoding convolutionally encoded codewords
CN100355201C (en) Reduced soft output information packet selection
CN1155160C (en) Method and apparatus for transmitting and receiving
US6343368B1 (en) Method and system for fast maximum a posteriori decoding
CN1178399C (en) Highly parallel MAP decoder
US7227908B2 (en) Information transmission apparatus, information transmission method, information reception apparatus, and information reception method
US6798852B2 (en) Methods and systems for extracting a joint probability from a map decision device and processing a signal using the joint probability information
US7159168B2 (en) Iterative decoding with likelihood weighting
CN1205791C (en) Systems and methods for receiving modulated signal containing encoded and unencoded bits using multi-pass demodulation
CN1104795C (en) Interference mitigation by joint detection of cochannel signals
US5633881A (en) Trellis encoder and decoder based upon punctured rate 1/2 convolutional codes
US6728323B1 (en) Baseband processors, mobile terminals, base stations and methods and systems for decoding a punctured coded received signal using estimates of punctured bits
CN1233139C (en) Effective band width cascade connected lattice code modem and decording method
EP2750299A1 (en) Interference cancellation in a multi-user receiver
CN1143471C (en) Method and system for fast maximum posteriori decoding
EP1820276A2 (en) Turbo decoder with stake heritage for data block redundant version decoding
CN1211931C (en) Memory architecture for MAP decoder
US20030018941A1 (en) Method and apparatus for demodulation
JP4675312B2 (en) Encoding device, decoding device, transmitter, and receiver
US7634703B2 (en) Linear approximation of the max* operation for log-map decoding
EP1564894B1 (en) Method and device for sequence estimation
CN1148006C (en) Method and decoder for decoding turbo code
JP2004248217A (en) Receiving apparatus and method

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C06 Publication
PB01 Publication
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: CLUSTER CO., LTD.

Free format text: FORMER OWNER: TELEFONAKTIEBOLAGET LM ERICSSON (SE) S-126 25 STOCKHOLM, SWEDEN

Effective date: 20150209

Owner name: AOPU DISI CELLULAR TECHNOLOGY CO., LTD.

Free format text: FORMER OWNER: CLUSTER CO., LTD.

Effective date: 20150209

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20150209

Address after: Texas, USA

Patentee after: Telefonaktiebolaget LM Ericsson (publ)

Address before: Delaware

Patentee before: Clastres LLC

Effective date of registration: 20150209

Address after: Delaware

Patentee after: Clastres LLC

Address before: Stockholm

Patentee before: Telefonaktiebolaget LM Ericsson

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20040324

Termination date: 20181213