Channel impulse response estimator for a viterbi equalizer

Abstract: A Viterbi equalizer can provide an accurate equalizing characteristic. This Viterbi equalizer comprises a synchronizing signal data detecting circuit for detecting a synchronizing signal data portion from a reception signal data series, a transmission line characteristic estimating circuit for modelling an impulse response between a transmitter and a receiver by comparing a synchronizing signal data detected by the synchronizing signal data detecting circuit with a reference signal, and a decoding circuit for decoding a transmission data series by using a Viterbi algorithm on the basis of a transmission model obtained from the transmission line characteristic estimating circuit.

Description: BACKGROUND OF THE INVENTION Field of the Invention: The present invention relates to a Viterbi equalizer in a receiver of a transmission system in accordance with the precharacterizing part of claim 1. Such a Viterbi equalizer is known from second IEE National Conference on Telecommunications, York, GB, April 2 to 5, 1989, pages 61-66, IEE London, 1989; XP-41168; L.B. Lopes "Performance of Viterbi Equalizers for the GSM System". Description of the Prior Art: Heretofore, the development of digital mobile telephones has been advanced so far in U.S., European countries and Japan. A mobile communication system such as the mobile telephone or the like is affected by a so-called multipath because a highrise building or the like is located between a high-speed mobile station such as automobile and a base station so that a transmission characteristic between the base station and the mobile station is deteriorated considerably. Therefore, a data transmission with less error is impossible. In addition, such an equivalent transmission characteristic is fluctuated with time. In the above-mentioned mobile communication system, in order to realize a reception with less error, an equalizing technique for correcting the above transmission characteristic becomes indispensable. Japanese laid-open patent publication No. 4-183041 describes as a conventional equalizing technique a Viterbi equalizer in which transmission data is decoded by using a transmission characteristic provided between the base station and the mobile station on the basis of a maximum likelihood sequence estimation. More specifically, according to an example of the prior art, in a Viterbi decoder in which a maximum likelihood reception is carried out following a time fluctuation of a characteristic of a transmission line, there is provided a Viterbi decoder which comprises a register for storing a plurality of sample values of a reception signal, a transmission line characteristic calculating circuit for receiving the plurality of sample values from the above register and estimating respective transmission line characteristics relative to a plurality of sequences, a branch metric calculating circuit for calculating distances of the plurality of samples by calculating hypothetical signal reception points relative to the respective sequences from the transmission line characteristics calculated by the transmission line characteristic calculating circuit and a Viterbi processor supplied with an output of the branch metric calculating circuit and judging a received signal on the basis of a Viterbi algorithm. According to this conventional Viterbi decoder, however, since a plurality of sample values of the received signal are compared, it cannot be expected that the transmission characteristic between a transmitter and a receiver is estimated accurately. In the Viterbi equalizer which is known from the above-mentioned article of L.B. Lopes "Performance of Viterbi Equalizers for the GSM System" as corresponding to the precharacterizing part of claim 1, the reception signal data series produced at the output of a low pass filter is fed as a whole, i.e. directly to a channel estimation and synchronizing unit. Said known Viterbi equalizer described by L.B. Lopes does not include a "separate" synchronizing signal data detecting means for separately detecting a synchronizing signal data portion from the reception signal data series, because the estimation procedure carried out in the channel estimation and synchronizing unit of the known Viterbi equalizer involves essentially the convolution of the received signal with the stored central 16 bits of the sequence. There is no comparison of the synchronizing signal data portion with a reference signal. EP-A-0 496 152 A3 describes a Viterbi equalizer having a variable number of states. A received radio signal is converted to the base band and coupled to a measurement unit which serves to examine those received signal. For measuring multipath interference by comparing a received signal with what that received signal is expected to be, the measurement unit captures a sample of a received data spanning an expected time of arrival of a training sequence, which is embedded in the received data signal, correlates this sample against the expected form of the training sequence to produce samples which represent a complex channel response and applies an input signal, which is provided from said produced samples to a windowing unit. The windowing unit serves for adapting the number of states of the Viterbi algorithm in dependence upon the multipath interference duration as measured on the basis of the training sequence. The function of the windowing unit of the known Viterbi equalizer is however very different from the function of a transmission line characteristic estimation means, as it is used in the present Viterbi equalizer, which models an impulse response between a transmitter and a receiver. OBJECTS AND SUMMARY OF THE INVENTION Therefore, it is an object of the present invention to provide an improved Viterbi equalizer in which the aforesaid shortcomings and disadvantages encountered with the prior art can be eliminated. More specifically, it is another object of the present invention to provide a Viterbi equalizer in which a transmission line characteristic can be estimated by using a synchronizing signal pattern. It is a further object of the present invention to provide a Viterbi equalizer in which a transmission line characteristic between a transmitter and a receiver can be accurately estimated by using a method of least squares. According to an essential aspect, the present invention provides a Viterbi equalizer in a receiver of a transmission system, said Viterbi equalizer comprising: transmission line characteristic estimating means for estimating a transmission line response characteristic on the basis of a reception signal data series, and decoding means for decoding the reception signal data by using a Viterbi algorithm on the basis of information provided by said transmission line characteristic estimating means, characterized in that synchronizing signal data detecting means are provided for detecting a synchronizing signal data portion from said reception signal data series and supplying said detected synchronizing signal data portion to said transmission line characteristic estimating means, said transmission line characteristic estimating means being arranged for modelling an impulse response between a transmitter and said receiver by comparing said synchronizing signal data portion detected and supplied by said synchronizing signal data detecting means with a reference signal; and said decoding means is decoding the reception signal data series on the basis of the modelled impulse response obtained from said transmission line characteristic estimating means. In an advantageous development of the Viterbi equalizer according to the present invention, said synchronizing signal data detecting means detects the synchronizing signal data portion by correlating the reception data with a synchronizing signal pattern previously stored. According to an advantageous development, said transmission line characteristic estimating means of said Viterbi equalizer employs a method of least squares when the modelled impulse response between the transmitter and the receiver is obtained by comparing said synchronizing signal data with said reference signal, wherein said reference signal is a reference signal of said detected synchronizing signal pattern. Another advantageous development of said Viterbi equalizer provides that said transmission line characteristic estimating means estimates a transmission line characteristic by correlating a synchronizing signal actually received with a known synchronizing signal pattern within the receiver in advance. According to a further development, said decoding means of said Viterbi equalizer of the invention calculates a branch metric of said Viterbi algorithm on the basis of said modelled impulse response obtained from said transmission line characteristic estimating means. According to a further advantageous development, the Viterbi equalizer of the invention provides the decoding means as calculating a branch metric of said Viterbi algorithm on the basis of said modelled impulse response obtained from said transmission line characteristic estimating means and further calculates a path metric using said calculated branch metric. The above and other objects, features, and advantages of the present invention will become apparent from the following detailed description of illustrative embodiments thereof to be read in conjunction with the accompanying drawings, in which like reference numerals are used to identify the same or similar parts in the several views. BRIEF DESCRIPTION OF THE DRAWINGS Fig. 1 is a block diagram showing a fundamental circuit arrangement of a Viterbi equalizer according to an embodiment of the present invention; Figs. 2A and 2B are diagrams used to explain arrangements of frame and slot of a communication channel of a so-called GSM (group special mobile) system, respectively; Fig. 3 is a diagram showing an example of a data series used in the present invention; Fig. 4 is diagram used to explain the present invention; Fig. 5 is a block diagram showing a generalized transmission line equivalent model; Fig. 6 is a block diagram showing a specific transmission line equivalent model; Fig. 7 is a diagram showing a trellis expression; Fig. 8 is a block diagram showing a logical unit of the Viterbi equalizer; FIG. 9 is a diagram used to explain a metric calculation and a remaining path; FIG. 10 is a flowchart to which references will be made in explaining operation of a Viterbi equalizer according to a first embodiment of the present invention; and FIG. 11 is a flowchart to which references will be made in explaining operation of the Viterbi equalizer according to a second embodiment of the present invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS A Viterbi equalizer according to embodiments of the present invention will hereinafter be described with reference to the drawings. FIG. 1 shows a fundamental circuit arrangement of a Viterbi equalizer according to an embodiment of the present invention. In this embodiment, the Viterbi equalizer shown in FIG. 1 is applied to a GSM (group special mobile) system that is adopted by mobile telephone used in European countries. As shown in FIG. 1, a reception signal supplied to an input terminal 1 is supplied to a branch metric calculating circuit 21 which constructs a Viterbi estimation unit 2 and also the reception signal is supplied to a synchronizing (sync.) signal data detection unit 3. Sync. signal data from the sync. signal data detection unit 3 is supplied to a transmission line characteristic estimation unit 4. A communication channel from the base station to the mobile station (automobile) according to the GSM system adopted in European countries has a frame arrangement shown in FIGS. 2A and 2B. Since each time slot is added at its central portion with a synchronizing signal pattern (SYNC pattern) having a well-known pattern as shown in FIG. 2B and then transmitted, the transmission line characteristic estimation unit 4 estimates an impulse response (hereinafter referred to as a channel response) interposed between the transmitter and the receiver by utilizing this sync. pattern. The GSM system uses a GMSK (Gaussian Minimum Shift Keying) modulation system. In this case, a signal of the high frequency transmission system is converted into a baseband signal by a demodulator and therefore the present invention will hereinafter be described as the signal processing in the base band. In the GSM system, there are designated in advance 8 kinds of data series and one of those data series is illustrated in FIG. 3. A conventional procedure for modellng the channel response by utilizing such sync. pattern will be described. Now, the case that a channel response is presented as shown in FIG. 4 will be described by way of example (this channel response is unknown in actual practice). In FIG. 4, the unit of the time base direction is equal to a symbol transmission interval. The sync. pattern in FIG. 4 is the sync. pattern of FIG. 3. Sync. data that will be received when the reception signal passes the transmission system having such channel response is expressed by the following equation (1): where y i is the reception signal, x i is the sync. pattern, and h i is the channel response. These values are those that were sampled by the symbol time interval T, respectively. When a reception signal corresponding to a sync. pattern portion is calculated in accordance with the equation (1), then there is obtained an output signal shown in FIG. 4. On the receiver side, known informations are the sync. pattern x i and the reception signal y i . In the conventional processing procedure of modelling in the transmission line characteristic estimation unit 4, the sync. data portion is detected by correlating the reception signal with the sync. pattern. Then, a cross-correlation function r j between the sync. signal data portion and the sync. signal pattern is calculated by the following equation (2): Then, a normalization is carried out by using a maximum value of this cross-correlation function r j . Cross-correlation functions thus calculated are illustrated in FIG. 4. The channel response is estimated by this cross-correlation function and an estimated value is supplied to the branch metric calculating circuit 21. After the channel response was estimated, transmission data series is decoded by using a Viterbi algorithm. FIG. 5 shows a generalized transmission line equivalent model. An example of a transmission line equivalent model in which the generalized transmission line equivalent model of FIG. 5 is practically arranged under the condition that the channel response length is limited will be described with reference to FIG. 6. If the transmission line equivalent model is arranged as shown in FIG. 6, such a transmission line equivalent model can be regarded as a convolutional encoder in which a constraint length is 4 and a coding ratio r is 1/1. This convolutional encoder is different from an ordinary convolutional encoder in that an adder 71 carries out a linear operation, symbols that are input to shift registers T 0 , T 1 , T 2 and T 3 are two values of <+1> and <-1> and respective outputs of the shift registers T 0 , T 1 , T 2 and T 3 are weighted by weighting amounts corresponding to channel responses h -1 , h 0 , h +1 and h +2 and which are then added together by the adder 71. A symbol G that is transmitted when the transmission line equivalent model is arranged as shown in FIG. 6, is expressed by the following equation (3): where <T j > represents the content stored in a register T j . FIG. 7 is a trellis diagram showing the transition of the internal state of the transmission line in the transmission equivalent model shown in FIG. 6. Three alphabets corresponding to respective state nodes S i in FIG. 7 represent internal states of the shift register in each time slot. The shift register takes values of <+1> and <-1> and these values are expressed as H and L for simplicity. Modifying the ordinary trellis diagram, FIG. 7 shows the state that a transition shown by a solid line occurs when an information input symbol <-1> is input and a transition shown by a broken line occurs when an information input symbol <+1> is input. Reception signal data Y k is input to the branch metric calculating circuit 21, which then calculates a likelihood concerning the transition. While some calculations for calculating a likelihood have been proposed so far, a Humming distance that is the most general evaluation measure in the Viterbi decoder is applied in a broad sense. A branch metric in time slot t(k) is calculated by the following equation (4): b (k, S i → S n ) = | Y k - G k | where Y k is the reception signal data and G k is the symbol transmitted from the equivalent transmission line model. G k takes the value calculated by the foregoing equation (3). The branch metric obtained in the branch metric calculating circuit 21 is supplied to an ACS (Add Compare Select) circuit 22. The ACS circuit 22 is comprised of an adder, a comparator and a selector. The ACS circuit 22 adds the branch metric and a path metric of one time slot before stored in a path metric memory circuit 23 and selects a smaller value as a remaining path. The path metric is a value which results from adding branch metrics in the remaining path. An output signal of the ACS circuit 22 is supplied through a normalization circuit 24 to the path metric memory circuit 23, and the output signal of the ACS circuit 22 is supplied to a maximum likelihood path detection circuit 25. The maximum likelihood path detection circuit 25 detects a path having a minimum path metric value and outputs the content corresponding to the detected path and stored in a path memory 26 as decoded data. The path memory 26 estimates and stores therein information bit series. FIG. 8 shows a logical unit which constructs the Viterbi equalizer. In FIG. 8, calculated amounts express the following contents, respectively. P(k-1, S i ) : path metric having a remaining path which reaches the state node S i in time slot t(k-1) P(k-1, S j ) : path metric having a remaining path which reaches the state node S j in time slot t(k-1) b(k, S i → S n ) : branch metric corresponding to a transition from the state node S i to the state node S n in time slot t(k) b(k, S j → S n ) : branch metric corresponding to a transition from the state node S j to the state node S n in time slot t(k) M(k-1, S i ) : path memory owned by the remaining path reaching the state node S i in time slot t(k-1) M(k-1, S j ) : path memory owned by the remaining path reaching the state node S k in time slot (k-1) <-1>, <+1> : information symbols that are estimated to be transmitted in time slot t(k) P(k, S n ) : path metric owned by the remaining path reaching the state node S n in time slot t(k) M(k, S n ) : path memory owned by the remaining path reaching the state node S n in time slot t(k) If the constraint length is taken as k, then there exist 2 k-1 states so that the logical units shown in FIG. 8 corresponding to the states 2 k-1 are required fundamentally. The scale of the path metric memory circuit 23 is reduced by providing the normalization circuit 24 similarly to the block arrangement of the Viterbi equalizer shown in FIG. 1. Also, an overflow that takes place when the path metric is calculated can be avoided. As a practical normalization processing, the minimum value of the path metric is detected and the detected value is subtracted from each path metric amount. Accordingly, there exist 2 k-1 remaining paths thus selected similarly to the 2 k-1 states. In each time slot, an operation for selecting the remaining path and an operation for updating a path metric corresponding to the path and the path memory 26 are repeated. It is known that, if the above operations are carried out for a sufficiently long time, then the remaining paths are merged to the same path. FIG. 9 shows such a state. A length in which the path is merged from the latest processing time point is called an end path length. How to update the path memory of FIG. 8 is determined depending on the respective states. For example, in the logical unit of "LLL", <-1> is determined. In the logical unit "HLL", <+1> is determined. Respective values will hereinafter be determined similarly. In the maximum likelihood judgement, the path having the minimum path metric is detected and the content of the path memory corresponding to the detected path is output as an information symbol of a timing point advanced by an amount of the end path length (about three to four times the constraint length in general). A signal processing flow of this Viterbi equalizer will be described with reference to a flowchart forming FIG. 10. Referring to FIG. 10, following the start of operation, when the reception signal data Y k is supplied to the input terminal 1, the sync. pattern is detected at step S1. Then, the processing proceeds to step S2, whereat a cross-correlation function of the sync. pattern of the reception signal data Y k and the sync. pattern previously stored is calculated by the transmission line characteristic estimation unit 4 to thereby estimate the channel response (step S3). Then, the branch metric calculating circuit 21 calculates the branch metric in step S4, and the processing proceeds to step S5, where the calculation of N'th state is started. In step S6, an address of a state-1 that is one time slot before is set. Then, the processing proceeds to step S7, where the path metric stored in the path metric memory circuit 23 at its address thus set is read out. Then, the read out path metric and the branch metric calculated at step S4 are added by the ACS circuit 22 and the added output is stored in a register P1 in step S8. In step S9, an address of a state-2 of one previous time slot is set. Then, the processing proceeds to step S10, where a path metric stored in the path metric memory circuit 23 at its set address is read out. In the next step S11, this path metric is added to the branch metric calculated at step S4 by the ACS circuit 22, and the added output is stored in a register P2. In the next steps S12 and S13, the values stored in the registers P1, P2 are compared and selected by the ACS circuit 22. Then the value thus selected is output at step S14. In step S15, the path metric memory circuit 23 is updated by this selected value, and the path memory 26 is also updated by this selected value at step S16. In the next decision step S17, it is determined whether or not the processing from step S5 to S16 is carried out 2 k-1 times. If a NO is output at decision step S17, then the processing returns to step S5. If on the other hand a YES is output at decision step S17, then the processing proceeds to step S18, where a path having the minimum path metric value is detected by the maximum likelihood path detection circuit 25. In the next step S19, a normalization is effected by subtracting the minimum value of the path metric from each of the path metric amounts. In step S20, an address of a maximum likelihood path is set by the maximum likelihood path detection circuit 25. In the next step S21, the content of the path memory 26 is output as decoded data. Comparing the channel response and the cross-correlation function r j shown in FIG. 4, it can be confirmed that the channel response can be estimated with a certain accuracy, whereas it is unavoidable that a "false impulse response" that should not appear also is detected. A cause for this lies in an autocorrelation function a j of the sync. pattern. This autocorrelation function a j is expressed by the following equation (5): Autocorrelation functions thus calculated are illustrated in FIG. 4. Study of FIG. 4 reveals that there exist several peaks having a considerably large level in addition to the main peak. There is then the risk that an accuracy required when the channel response is estimated is deteriorated by the above-mentioned peaks. According to another embodiment of the present invention, since the sync. data is used as a reference signal and the impulse response between the transmitter and the receiver is modelled by using the least squares method, the impulse response between the transmitter and the receiver can be uniquely modelled. The model thus determined is a model whose error becomes minimum in a sense of the least squares estimation, which as a result can provide a satisfactory equalizing characteristic. In accordance with the second embodiment of the present invention, the transmission line characteristic estimation unit 4 shown in FIG. 1 is arranged as follows. The transmission line equivalent model is formed as the channel response as shown in FIG. 5. If so, the signal y i that is predicted to be received is expressed by the aforementioned equation (1). On the other hand, expressing the signal Y i that is actually received, an error ε concerning the i'th symbol is expressed by the following equation (6): ε i = y i - Y i Calculating a square sum E of this error ε i , we have: Then, an impulse series h n is determined so as to minimize this error E. In this embodiment, the least squares method is applied so that the above-mentioned equation (7) is partially differentiated with respect to the impulse series h n as: Substituting n = -k, -(k-1), ..., 0, ..., +(k-1), +k into the above equation (8) yields the following simultaneous equations: The coefficient matrixes in the above simultaneous equations (9) become symmetrical matrixes so that calculations for each element need not be carried out for all elements. Further, in order to solve the above simultaneous equations, the coefficient matrixes are resolved in an LU fashion. The transmission line characteristic estimation unit 4 according to this embodiment can determine the channel response with high accuracy by the above-mentioned means. A rest of the arrangements is formed similarly to that of the Viterbi equalizer shown in FIG. 1 and therefore need not be described in detail. The signal processing flow of the Viterbi equalizer according to this embodiment will be described with reference to FIG. 11. Referring to FIG. 11, following the start of operation, when the reception signal data Y k is supplied to the input terminal 1, the sync. pattern portion is detected at step S31. The detection of the sync. signal pattern portion is carried out by correlating the reception data Y k with the sync. signal pattern previously stored. Then, in step S32, an impulse response between a transmitter and a receiver is modelled in the transmission line characteristic estimating unit 4 by using a method of least squares on the basis of a reference signal of the sync. signal pattern portion thus detected. Then the channel response is estimated in step S33. Then, the branch metric calculating circuit 21 calculates the branch metric in step S34, and the processing proceeds to step S35, where the calculation of N'th state is started. In step S36, an address of a state-1 that is one time slot before is set. Then the processing proceeds to step S37, where the path metric stored in the path metric memory circuit 23 at its address thus set is read out. Then the read out path metric and the branch metric calculated at step S34 are added by the ACS circuit 22, and the added output is stored in a register P1 in step S38. In step S39, an address of a state-2 of one time slot before is set. Then, the processing proceeds to step S40, where a path metric stored in the path metric memory circuit 23 at its set address is read out. In the next step S41, this path metric is added to the branch metric calculated at step S34 by the ACS circuit 22, and the added output is stored in a register P2. In the next steps S42 and S43, the values stored in the registers P1, P2 are compared and selected by the ACS circuit 22. Then the value thus selected is output at step S44. In step S45, the path metric memory circuit 23 is updated by this selected value, and the path memory 26 is also updated by this selected value at step S46. In the next decision step S47, it is determined whether or not the processing from step S35 to S46 is carried out 2 k-1 times. If a NO is output at decision step S47, then the processing returns to step S35. If on the other hand a YES is output at decision step S47, then the processing proceeds to step S48, where a path having the minimum path metric value is detected by the maximum likelihood path detection circuit 25. In the next step S49, a normalization is effected by subtracting the minimum value of the path metric from each of the path metric amounts. In step S50, an address of a maximum likelihood path is set by the maximum likelihood path detection circuit 25. In the next step S51, the content of the path memory 26 is output as decoded data. According to the first embodiment, since the impulse response between the transmitter and the receiver is modelled by using the sync. signal pattern portion as the reference signal, there is then the advantage that the impulse response between the transmitter and the receiver can be modelled uniquely. Furthermore, according to the second embodiment of the present invention, since the transmission line characteristic is estimated by using the method of least squares in the transmission model and the error of the model becomes minimum, there is then the advantage that a satisfactory equalizing characteristic can be obtained. Having described preferred embodiments of the invention with reference to the accompanying drawings, it is to be understood that the invention is not limited to those precise embodiments and that various changes and modifications could be effected therein by one skilled in the art without departing from the scope of the invention as defined in the appended claims.

Claims: A Viterbi equalizer in a receiver of a transmission system, said Viterbi equalizer comprising: transmission line characteristic estimating means (4) for estimating a transmission line response characteristic on the basis of a reception signal data series, and decoding means (2) for decoding the reception signal data by using a Viterbi algorithm on the basis of information provided by said transmission line characteristic estimating means (4), characterized in that synchronizing signal data detecting means (3) are provided for detecting a synchronizing signal data portion (SYNC) from said reception signal data series and supplying said detected synchronizing signal data portion (SYNC) to said transmission line characteristic estimating means (4), said transmission line characteristic estimating means being arranged for modelling an impulse response between a transmitter and said receiver by comparing said synchronizing signal data portion (SYNC) detected and supplied by said synchronizing signal data detecting means (3) with a reference signal; and said decoding means is decoding the reception signal data series on the basis of the modelled impulse response obtained from said transmission line characteristic estimating means (4). The Viterbi equalizer according to claim 1, wherein said synchronizing signal data detecting means detects the synchronizing signal data portion by correlating the reception data with a synchronizing signal pattern previously stored. The Viterbi equalizer according to claim 2, wherein said transmission line characteristic estimating means (4) employs a method of least squares when the modelled impulse response between the transmitter and the receiver is obtained by comparing said synchronizing signal data with said reference signal, said reference signal being a reference signal of said detected synchronizing signal pattern. The Viterbi equalizer according to claim 1, wherein said transmission line characteristic estimating means (4) estimates a transmission line characteristic by correlating a synchronizing signal actually received with a known synchronizing signal pattern within a receiver in advance. The Viterbi equalizer according to claim 1, wherein said decoding means (2) calculates a branch metric of said Viterbi algorithm on the basis of said modelled impulse response obtained from said transmission line characteristic estimating means (4). The Viterbi equalizer according to claim 1, wherein said decoding means (2) calculates a branch metric of said Viterbi algorithm on the basis of said modelled impulse response obtained from said transmission line characteristic estimating means (4) and further calculates a path metric by using said calculated branch metric.

Inventors: SATO TERUO, SATO, TERUO, SATO, TERUO, C/O SONY CORPORATION
Applicants: SONY CORP, SONY CORPORATION
Application year: 1993
Main classification: H04L25
IPC codes: H04L25, H04B3, H03H21, H03H15, H04B1, H03H17