GB2425236A - An analog equaliser characterised by iterative means arranged in operation to generate an estimate of marginal posterior expectations for received bit values - Google Patents

Abstract

Prior art analog solutions to MIMO detection have the drawback that the number of transistors used increases exponentially in proportion to the number of receiver channels. The current invention aims to find a lower complexity solution to analog equalisation for applications such as MIMO detection and mass storage readers. The analog equaliser includes at least one analogue processing block (APB) arranged in operation to iteratively generate an estimate of marginal posterior expectations for received bit values. A further independent claim states that the equalisation method comprises the step of passing a plurality of log-likelihood marginal posterior expectations to an APB and generating in the APB a revised estimate of the log-likelihood marginal posterior expectations using coordinate descent optimisation.

Description

Apparatus and Method of Egualisation The invention relates to an apparatus

and method of equalisation, and in particular to an apparatus and method of equalisation for MIMO decoding and reading of recordable media with reduced component complexity.

In modern high-speed wireless communications networks, multipath signal propagation is an increasingly significant problem. In traditional wireless communication, a transmit antenna emits an electromagnetic (EM) signal to a receive antenna over an interve ing space. However, any obstructions to the signal within that space scatter the EM signal, resulting in copies of the signal reaching the receive antenna at different times and at different intensities via different paths; an effect known as channel spread. In a digital signal, channel spread results in an overlap between successive received bits, and this reduces the confidence in any given bit value received.

To increase bit transmission rates requires shorter bit representations. Consequently, the overlap caused by the same channel spread correspondingly increases, making disambiguation of the received bit stream more difficult. Therefore in high-speed wireless networks, there is a need to mitigate the effect of channel spread.

One approach is multiple input, multiple output (MIMO) communication, wherein multiple transmitter and receiver aerials are used. MIMO systems improve communications robustness by providing multiple, path-independent copies of the transmitted data. This is typically achieved by use of space-time coding techniques, for example Alamouti orthogonal space-time block coding (see S. M. Alamouti, A Simple Transmit Diversity Technique for Wireless Communications, IEEE Journal on Select Areas in Communications, vol. 16, no. 8, Oct. 1998). The result is a set of received signals in which path induced interference differs for each copy of the data, simplifying disambiguation of the common and disparate signal components.

However, MIMO decoding is not a trivial task. Typical detectors use digital signal processing (DSP) methods to decode the MIMO signal; this may involve multiple sampling of each candidate bit signal for each MIMO receiver, and calculating and aggregating bit value probabilities for each sample. These steps incur large computational costs relative to the actual bit rate. The computational costs in turn carry a corresponding power cost that is significant in portable MIMO devices, and can cause a processor bottleneck that limits throughput in high data rate applications. This problem also occurs in other applications where a receive signal is equalised to estimate the source signal, such as surface reading in magnetic storage media.

Recently, an alternative method of MIMO detection has been proposed using analog circuitry rather than digital signal processing (see Piechocki, R. J., Garrido, J., McNamara, D., and McGreen, J., Analogue MIMO detector: The Concept and Initial Results', IEEE First International Symposium on Wireless Communications Systems, Mauritius, 20 22nd Sept. 2004).

Advantageously, analog circuitry does not require sampling of the incoming signal, and can operate directly on the soft' (probabilistic) values observed by the receivers.

Moreover, the circuitry can be constructed to operate in parallel on the multiple receiver channels.

In consequence, equivalent detector processing can be performed several orders of magnitude more quickly than by the DSP equivalent, whilst simultaneously requiring less power.

However, the analog solution to MIMO detection proposed in Piechocki et. al. above has the drawback that the number of transistors used increases exponentially in proportion to the number of receiver channels. Consequently it is desirable to find a lower complexity solution to analog equalisation for applications such as MIMO detection and mass storage readers.

Accordingly, aspects of the present invention seek to mitigate, alleviate or eliminate the above-mentioned problem.

In one aspect of the present invention, an analog equaliser comprises at least a first analog processing block arranged in operation to generate successively improved estimates of marginal posterior expectations (MPEs) for received bit values.

In one configuration of the above aspect, successive estimates of the MPEs are obtained using a coordinate descent minimisation means.

In another configuration of the above aspect, the analog equaliser updates the MPEs by a temperature factor.

In another configuration of the above aspect, in use a single analog processing block feeds its own outputs back to form its own inputs in successive iterative cycles.

In another configuration of the above aspect, the analog equaliser comprises a plurality of K sets of K-I Rk'.tanh calculation circuits, k=l, . . .,K, where K is the number of MPE estimates, means to sum each of the K sets of K-I Rk'.tanh calculation outputs, means to subtract each said sum from a respective filtered signal Zk, k=l,. . .,K, and means to scale an output signal for each MPE estimate in inverse proportion to a mean

field annealing factor T.

In another aspect of the present invention, an ASIC comprises an analog equaliser as described herein.

In an aspect of the present invention, a multiple input, multiple output detector comprises an analog equaliser as described herein.

In an aspect of the present invention, a wireless communications device comprises an analog equaliser as described herein.

In an aspect of the present invention, a bulk storage device comprises an analog equaliser as described herein.

In another aspect of the present invention, a method of equalisation comprisses the step of passing a plurality of log-likelihood marginal posterior expectations to an analog processing block (APB), the APB in turn generating a revised estimate of the log- likelihood marginal posterior expectations using coordinate descent optimisation.

Although the present invention has been described hereinabove with reference to a number of separate aspects, in accordance with the present invention any aspect of the present invention described previously can be used in conjunction with any other aspect of the present invention.

Embodiments of the present invention will now be described by way of example with reference to the accompanying drawings, in which: Figure IA is a schematic diagram of a single analog processing block in accordan e with an embodiment of the present invention.

Figure 1 B is a schematic diagram of an analog MIMO detector comprising a plurality of analog processing blocks, in accordance with an embodiment of the present invention.

Figure 2 is a MOS transconductance amplifier, for use in an embodiment of the present invention.

Figure 3 is an absolute value circuit, for use in an embodiment of the present invention.

Figure 4 is a Gilbert multiplier, for use in an embodiment of the present invention.

Figure 5 is a tanh calculator circuit in accordance with an embodiment of the pres invention.

Figure 6 is a tanh calculator circuit in accordance with an embodiment of the present invention.

Figure 7 is a graph showing the response of a tanh calculator block in accordance with an embodiment of the present invention, with input voltage on the x-axis and output current on the y-axis.

Figure 8 is a subtraction circuit, for use in an embodiment of the present invention.

Figure 9 is a transresistance circuit, for use in an embodiment of the present invention.

Figure 10 is a voltage shifter, for use in an embodiment of the present invention.

Figure 11 is a graph showing the DC response of a transresistance circuit, with input current on the x-axis and output voltage on the y-axis.

Figure 12 is a current to voltage converter, for use in an embodiment of the present invention.

Figure 13 illustrates sharing an absolute value circuit between corresponding tanh calculator circuits in accordance with an embodiment of the invention.

Figure 14 illustrates the arrangement of circuits to estimate one of three received bits in an analog MIMO detector in accordance with an embodiment of the present invention.

Figure 15 is a schematic diagram of an analog equaliser arranged in conjunction with a channel detector, in accordance with an embodiment of the present invention Figure 16 is a schematic diagram of a single recurrent analog processg block in accordance with an embodiment of the present invention.

Figures 1 7A and 1 7B are graphs illustrating marginal posterior expectation estimates for a recurrent analog processing block and a plurality of analog processing blocks respectively, where time is on the x-axis, and estimate values are expressed as voltages on the y-axis.

An analog MIMO detector is disclosed. In the following description, a number of specific details are presented in order to provide a thorough understanding of embodiments of the present invention. It will be apparent to a person skilled in the art, however, that these specific details need not be employed to practice the present invention.

Theoretical background

A model of MIMO signal reception well known in the art is y=Hx+n (1) where x is a transmitted data vector from NT transmit antennas, y is the received vccor, H is the matrix of channels between each transmit and receive antennas and n is a vector of independent Gaussian noise distributions.

Applying a channel-matched filter HH, the signal model becomes z=H'y=H'Hx+ H11nRx+v (2) where R = H"H and v = H"n N(O,oR).

Thus the likelihood of a value in filtered signal z is f(zx, R, cr2) exP[__(z - Rx' R' (z - Rx)]. (3) j2ro-JRI Assuming uniform priors, then f(xz,R,a,)x f(zx,R,c,). (4) The marginal posterior expectations mk =E{xkzk} (5) are then computed, and the sign taken, to derive a binary value.

It can be further shown that mk E{tanh[zk_RkkxkJJ}. (6) (e.g. see Fabricius, T and Winther, 0, Correcting the bias of subtractive interference cancellation in CDMA - Advanced mean field theory', submitted to IEEE trans.

Information Theory, and soft published at http://isp.imm.dtu. dk/staff/winther/ fabricius.ieeeinformationtheory.pdf, as at 15/03/2005).

However, it is difficult to evaluate the expectations in equation (6), as the tanh function is asymptotic in character.

Description of new approach

Notably, a mean field approach to the problem may be taken by approximating the expectations (tanh(...)) as tanh((...)). Applying this approximation to equation (6) gives mk = E{tanhHzk - RkkxkJJ} a,, k!=k tanh[E{_(zk - Rkkxk.J}J k!=k mk = tanh[_[zk _RkkmkJJ. (7) k'!=k Thus, applying a mean field approach to the model of marginal posterior expectations (hereinafter MPE5) generates a recursive solution for mk. Iterating equation (7) corresponds to a coordinate descent minimisation of the cost function of the free energy, thus guaranteeing that at least a local minimum value is found.

To improve the chance of finding a global minimum in the descent minimisation problem, mean field annealing may be applied. In mean field annealing, a variabh temperature Tis used in place of the noise Q2,, in equation (7). Thus an annealing schedule may be used wherein following an initial and comparatively high value of T, subsequent iterations of equation (7) use decreasing values back toward or even below the original noise floor (See Fabricius, T and Winther, 0, ibid). Note that such mean field annealing does not result in probabilistic gradient ascent dependent on T, as in simulated annealing.

The inventors of the present invention have appreciated that a suitable structure m A be constructed in analog circuitry to iterate toward a solution for mk in the fashion described above, so providing the building blocks for a lower complexity analog equaliser for MIMO detection based upon approximating mk rather than providing an exact calculation.

Applying mean field annealing to equation (7) gives a general form for mk as; mk =atanh[--zk 1?kmkJJ () T k'!=k where temperature T lowers with each re-estimation of mk, and a is a constant.

In principle, the tanh function is readily implementable using analog circuits; for example, the well-known Gilbert multiplier circuit inherently computes tanh for its inputs, normally relying on the lineansing approximation tanh(x) x for small values to perform its multiplier operation.

In analog calculations, means such as mk are generally computed as differential currents.

However, in a recursive problem such as that of equation (8), distributing the new values of one mk iteration at the input to a further iteration would require k-2 copies of each current using current mirrors resulting in a large number of components.

To reduce this problem, the posterior expectations mk are re-expressed as log-likelihood ratios using an inverse hyperbolic tangent: Lk = 2 tanh' (9) Substituting equation (9) into (8) gives Lk -aRk tanh(LJJ (10) The output of (10) may be expressed using a differential voltage, rather than a differential current. This allows the MPE estimates of one iterative stage to be directly connected to the inputs of another with no need for copying.

Overview of the analog process Referring now to Figures 1A and IB, in an embodiment of the present invention, the analog MIMO detector employs a plurality of analog processing blocks (APBs) (102, 104, 108, collectively referenced as 110), each implementing an iteration stage of equation 10.

APB 104 is taken as an example in Figure 1 B and discussed in functional terms below.

For clarity, implementation issues are deferred to specific sections further on.

A plurality of K log likelihood ratio MPEs Lk, k=1, ...K, are input 210 to the APB 104 as voltages. These inputs are distributed to a plurality of tanh calculators (collectively referenced 240). For each klh bit whose MPE is being estimated, there are K-i tanh circuits. All but the ICJh input is passed to the respective K-i tanh calculators in each of these K sets.

Thus, for example, the plurality of K-i tanh calculators 243 associated with L3 will receive all Lk inputs other than L3 itself. In addition, each set of K-i tanh calculators receives the relevant elements of the channel estimation matrix R. The output of each individual tanh calculator is, tanh&J for one value of k for its corresponding Lk.

The outputs of the K-i tanh functions for each Lk are then summed, and subtracted from z 260. The output is then a partial calculation Zk - a Rk tanhLJ of equation (10) k!=k 2 for each Lk.

The value of Trelevant to the current iteration stage is then used in scaling 280 an output voltage representing the updated estimate for each Lk 220, which may then be used as input for the next iteration, or taken as the final estimate.

The implementation of the functional elements disclosed above will now be discussed in detail.

Tanh function (1) - nonlinearity and multiplication Referring now to Figure 2, in an embodiment of the present invention the analog circuit to perform the hyperbolic tangent and multiplication comprises a simple differential pair with a mirror to subtract the output currents, i.e. a basic transconductance amplifier.

The output current of this circuit is given by out = b tanh[3LJ.

Therefore, an output current Rk tanhL) can be obtained by representing log likelihood ratios Lk as differential voltages (normalised to the thermal voltage), and defining the bias current b as proportional to the input value Rk' (i.e. lb = IrefRk) from the channel estimator.

Advantageously, the plurality of output currents may be summed together without the need for any additional transistors according to Kirchoff's current law, to provide partial calculation a Rk tanhLJ of equation (10).

k'!=k 2 However, whilst the bias current to the above circuit must be positive, the input vahes of Rk' may be positive or negative. This discrepancy must be resolved.

Tanh function (2) - absolute values and sign management Referring now to Figure 3, in a first embodiment of the present invention, to overcome the problem of bias current sign a translinear circuit is provided to obtain the absolute value of Rk'.

The circuit of Figure 3 comprises two translinear loops. The first ioop uses transistors M2-M1-M3-M4, so that 2'1 = 1314. Substituting the values of the currents gives: sma1Ijsmo11 = outl (IOutl - 1,). The second translinear loop is formed by M4 and M5, forming a simple current mirror so that oui2 = 14 = oU! - ln* These equations can be solved to obtain the output currents in terms of the input currents: outi = ( + + out? = ( gn + 4Iaii + where small is a small bias current compared to rn* This circuit then behaves like an absolute value circuit that for positive input currents obtains I 0, and for negative input currents obtains outI 0, out 2 in Obtaining the absolute value of Rk' allows the transconductance amplifier to operate as desired. However, the correct sign of the input value Rk' must still be preserved for subsequent reincorporation.

In an embodiment of the present invention, the sign is obtained from the voltage difference between the gates of M2 and M4 of the tanh block. If I, >0, then the current flows through M3; otherwise, the current flows through M4. Thus depending on the sign of the voltage at the gate of M4 will either be close to Vdd or to ground. A differential voltage representing the sign of the input signal is obtained when the voltage at the gate of M4 is compared with a reference, for example the gate of M2.

To apply this sign to obtain the desired expression one may use a Gilbert multiplier, as depicted in Figure 4.

The output current of the multiplier is I, = small tanh, tanh J. Therefore, the combined circuit uses the absolute value of the current representing the value Rk' and a log-likelihood ratio voltage as inputs, and further applies the sign of the bias current, to obtain the final following expression: oul = rof sgn(R) tanh[fJ = IrejR. tanh(-J (11) Note that using a four quadrant Gilbert multiplier may suggest the use of Rk' as an input voltage, without the need for using an absolute value process. However, simulations have shown that the values of matrix R are typically too big for the tanh(x) x approximation to hold for multiplications.

The combined tanh function circuit comprising the transconductance and absolutc values circuits disclosed above can be seen in Figure 5.

However, in a preferred embodiment of the present invention, the sign of Rk' may be more efficiently managed using the tanh function circuit illustrated in Figure 6.

Here, transistors M5 to M9 form an absolute value circuit, and by virtue of the use of NMOS rather than PMOS transistors these form a current sink rather than a current source. Consequently, currents Ij and Iay2 pass into this sink and so the absolute value circuit can connect directly with the tanh block formed by transistors Ml to M4.

This preferred embodiment operates in a similar manner to the embodiment described previously. A first translinear loop formed by M5 and M6 acts in a clockwise direction, and a second translinear loop formed by M7 and M8 acts in an anticlockwise direction.

Consequently the current through M5 and M6 is forced to be a small current small and the currents through M7 and M8 are Iay2 and (Iav2+Iin) respectively. This arrangement gives current equation small = av2 (av2 + I,), for which the solution is r r AT 2 av2 - in V in 11 small Transistors M8 and M9 form a current mirror so that the current through M9 is also (Iav2+hn). Therefore, Jay] is expressed by avi = in + av2 = in + + 4lsmaii2) For Iin>>Ismall, this circuit therefore behaves like an absolute value circuit that for positive input currents produces api, av2 0 and for negative input currents produces api 0 av2 1inL Similarly, the sign also feeds directly into the tanh calculation as a function of whether Jay] or av2 is non-zero, as opposing differential pairs in the tanh block receive the input current depending upon the original sign.

The differential voltage applied to both differential pairs is the same as in the previous embodiment, but with the terminals swapped, i.e, +AV is applied to one and -AV to the other, inverting the polarity between them. Finally, output currents of the differential pairs are summed appropriately (positive with negative and negative with positive), where positive and negative refer to the terminal of zW applied to each branch; for example in the case of Figure 6, the current that comes from the first pair from the Vii branch is summed with the current coming from the second pair from the V12 branch.

This arrangement preserves the sign information, and so the final current = J,, tanh I\V is obtained by subtracting 1o2 from ol. The n' in the preceding expression is a relative scaling factor reflecting the use of MOS transistors instead of BJTs.

Note that as there is a plurality (k-i) of tanh calculations to perform prior to sumniation for the kt' MPE, it may be advantageous to keep 1o2 and Jj separate, sum all 102s and J is, and then perform the subtraction to obtain an aggregate tanh calculation. This would then only require one current mirror, but the currents involved would be much higher than for individual subtractions and so performance could be affected. The e two alternative approaches would therefore be a matter for selection by a person skilled in the art, dependant upon the current levels used in their specific detector hardware.

For the embodiment of either Figure 5 or Figure 6, the response of the circuit for different voltage and current inputs is shown in Figure 7, demonstrating Rk'-scaled tanh function outputs.

In an alternative embodiment of the present invention, the whole circuit of Figure 6 is implemented using PMOS transistors and the absolute value circuit of Figure 3, but now changing the differential pairs to use PMOS transistors. In this configuration, both parts of the circuit work as current sources and so direct connection between them would similarly be possible.

Summation As noted previously, for either of the above embodiments, the Out for each of the plurality of Rk' k' !=k tanh calculations are then simple to sum using Kirchoff's law, as prescribed in equation 10.

Subtraction from zk Referring now to Figure 8, a subtractor circuit as illustrated that may be used for the first embodiment of the tanh circuit as seen in Figure 5. This subtractor circuit acts as a current mirror to subtract currents. It takes six currents with duplicate values for zk and for Rk' (as split by the absolute value circuit), and then sums the positive and negative parts before finally obtaining the difference between the two resulting currents.

In the preferred embodiment of the tanh circuit as seen in figure 6, subtraction from Zk may be more simply achieved by swapping the input voltage terminals to produce = -I, tanhi AV * This is then added' to h using Kirchofs law, for example 2nVT) during summation of the k-i tanh calculations.

In either case, the output current represents the partial calculation Zk aR1 tanhLi-J of equation 10.

k'!=k 2 Formatting the result (1) - current to voltage conversion and annealing The result of the above processes gives an interim a-posteriori log-likelihood ratio for the received symbol mk as seen in equation 10, and expressed as a current: = 1rej L = ref [zk - tanh, where ref is a chosen reference current used in k!=k 2 mapping the values of Rk and Zk within the circuits, i.e. Iz=Iref*zk and Irlref*Rk.

However, recall that the tanh calculation circuit described above takes a differential voltage representation of the log-likelihood ratios as its inputs. Therefore, in order to pass the result of the above process to a subsequent iterative stage of the equaliser, it is necessary to reexpress the output current as a differential voltage normalised by the thermal voltage. During this step, the mean field annealing factor 2a/T may also be included.

To this end, a circuit is needed that linearly converts a current to a voltage in the form A V0 = VT ref Typically, current-to-voltage conversion is by means of diode-connected transistors that have a logarithmic response, i.e. = V log-k-. ref

In an embodiment of the present invention, a circuit similar to the absolute value cir'uit of Figure 3 may be used to provide current-tovoltage conversion; as noted previously in relation to this circuit, it only acts as an absolute value circuit when bias current sma// is small compared to I,,,. However, when this condition is not met, the circuit acts as a current to voltage converter, albeit with a response that follows an inverse hyperbolic tangent relationship rather than a truly linear one. Renaming small as Ib: !( +J4I +) AV1=V1og 1 (- +J4I +,,) + /J

-

- T og In =2Vtanh b (12) Fortunately, however, like tanh the inverse is approximately linear for small values and tanh'(x)x when x is less than about 2. Therefore providing that I,fl<2Ib, such a circuit can be used.

Figure 9 illustrates a transresistance circuit, incorporating the above absolute value circuit acting as a current-to-voltage converter, to implement equation 12. Recalling that a current = I, . L currently gives the log-likelihood ratio within the APB, one obtains a voltage representation given by

LL

= 2VT tanh1 VT (13) 4+1f4/J As noted above, such a circuit only acts approximately linearly up to some maximum input value ImaxM*Ib, where b is the bias current. In the case of the circuit of Figure 9, this value Mis approximately 2.

Therefore in order to linearly obtain log likelihood voltages up to some value Lm, t is necessary to ensure that Lmax =M or thereabouts. This is achieved by setting = Linax --. The transresistance block consequently outputs AV V =V rn out' T 7 bias max ref Thus for the circuit of Figure 9, the bias current must be set to umax = Lmax Then, for example, if one sets the maximum Lm to 10, b = rej and the output of the circuit will give the value AV01 V In an embodiment of the present invention, the annealing factor 2a/T is also incorporated into the current to voltage converter scaling calculation. This may bc achieved simply by substituting I,, = L with I,, = 2-L,,iax Formatting the result (2) - voltage scaling Having obtained a differential voltage output using the circuit described above, it will therefore be necessary to amplify it by a restorative factor. In the above-enumerated example this would be a factor of 5 to obtain the desired value for the log-likelihood ratio.

Referring now also to Figure 10, in an embodiment of the present invention, the amplification can be performed with a MOS level shifter, as seen in Figure 10 and as incorporated into the transresistance circuit of Figure 9.

Providing that the transistors are working in strong inversion, the circuit has a transfer function of the type AV0 / AV where W is the width of the transistors, L is the length, and 1 and 2 refer to the bottom and top pair of transistors in Figure 10 respectively. Because transistors are typically fabricated with similar lengths, in order to amplify the voltage by a factor of 5 therefore requires the bottom pair of transistors to be 25 times wider than the top pair.

In an embodiment of the present invention, the annealing factor 2a/T can be included during amplification. Note that because the value of T varies as a position of each analog processing block in within the equaliser's iterative chain, further including the annealing factor 2a/T means that the amplification, and therefore the relative width of the bottom and top transistors in the circuit of Figure 9, will need to be selected to reflect the value of Tin each analog processing block.

Thus for example, given an equaliser in which a=0.5, Lmax=l 0 and K=1, and wherein T varied from 10 to 0.1 over the annealing schedule, then the amplification factor in the first analog processing block would be 2*0. 5*10/(1*10)=1, whilst the amplification factor in the last would be 2*0. 5*10/(1*0.1)=100. Referring now to Figure 11, the substantially linear current-to-voltage

relationship of the transresistance circuit incorporating the MOS level shifter is shown in the graph, with input current on the x-axis and output voltage on the y-axis.

The person skilled in the art will appreciate that alternative arrangements are possible to implement the transresistance circuit, and these are envisaged within the scope of the present invention. Similarly, the person skilled in the art will appreciate that other current to voltage converters are possible, and these are likewise envisaged within the scope of the present invention. For example, referring now to Figure 12, in an alternative embodiment of the present invention, the circuit of Figure 12 may be employed as a current to voltage converter. In this circuit, Ml and M2 form a current mirror so that 2 = A Moreover, it can be seen that I = A + rn Currents Ij and 2 are then forced to sum 4; + b* Solving this system of equations for Ii and 2 gives: = ! (1, + i,,,) and 2 = (ii, - i,,).

By virtue of the aforementioned logarithmic relationship, the output voltage is

I

AV V. log[ = V. log'b = v.io[ : IhI = 2V tanhIJ.

Thus this circuit also provides acceptable linearity, but for a smaller range (M O.5) than the previously disclosed adapted absolute value circuit of Figure 3 (M 2).

Considerations for multiple analog processing blocks In one analog processing block, the output of matrix R=R, ij=1,..,K from the channel estimator provides the required inputs to the mk bit estimations in the form of the kth row of the matrix (excepting Rkk as it is not needed).

However, for a MIMO detector comprising an iterative chain of discrete analog processing blocks, each block requires a copy of R. Recall, however, that generating copies using current mirrors assumes that the current is positive, whilst values of R may be positive or negative. In an embodiment of the present invention, values of R are passed through an absolute value circuit to obtain two value streams (one positive current and one null current, which is which being dependent on the sign of the input value as detailed previously). The positive currents may then be copied.

Advantageously therefore, one would only need one absolute value circuit for each corresponding tanh calculator in the chain of analog processing blocks. Such an arrangement is depicted in Figure 13.

In Figure 13, one sees that the absolute value circuit of Figure 3, using PMOS transistors, can be directly input to a tanh circuit as the current minors change the direction of the currents, allowing direct connection. Consequently the tanh calculator depicted is that of Figure 6 minus the absolute value circuit shown therein, i.e. only transistors Ml to M4, directly connected to the copies of the absolute value currents.

In summary, Figure 14 provides a schematic diagram illustrating how an analogue processor block (APB) according to an embodiment of the present invention would be arranged to decode bit I of 3 bits (i.e. K=3) using process described herein. The individual sections of the APB are labelled as follows: AV: The absolute value circuit of Figure 3; MIRR: A current minor; TANT-I: The tanh circuit of figure 6, without the absolute value elements (i.e. just transistors Ml to M4); TRANSR: The transresistance circuit of Figure 9, and; SUBTR: The subtraction circuit of Figure 8.

Relation between numbers of receiver channels and transistors The total number of transistors needed, based on an embodiment of the invention, is in the order of Trans=N(K(K-1).l8+K.l3) (14) where N is the required number of iterative stages and K is the number of bits to decode, based on the number of transmit antennas and the modulation scheme.

Thus, advantageously, the number of transistors used is proportional only to the square of the transmitted bits, rather than being exponential as is the case in the prior art.

Thus for this example, a quadrature phase shift key (QPSK) system with 10 transmit antennas (K=20), and decoding using 10 stages (N8=10), would need 71,000 transistors.

This number of transistors may still be considered too many for some applications, however, so additional means may be considered to reduce the number of stages used.

Variant embodiments of the analog MIMO detector Firstly, referring now to Figure 15, a channel detector is shown arranged in operation to pass estimates of Lk, k=1, . . . ,K to the analog equaliser. In an embodiment of the pr sent invention, the inclusion of a channel detector is expected to significantly reduce the number of iterative stages required to 2 or 3 in the example above. Consequently, the number of transistors used would fall to 14,000 or 21,000 respectively.

The analog equaliser by itself starts with an initial log likelihood value of zero for each Lk and then iterates towards a solution. Consequently, the provision of an initial estimate for each Lk from a channel detector avoids the need to iterate from zero to an approximate value, so reducing the number of iteration stages needed. Optionally, as shown in Figure 15 the output of the equaliser may then be passed back to the channel detector.

Referring now to Figure 16, in an embodiment of the present invention a single analog processing block (APB) comprises a variable gain amplifier instead of the MOS Ic vl shifter of Figure 10. In this arrangement, the outputs of the APB are fed back to its inputs in continuous time, while the amplification gain and temperature T are also altered continuously as appropriate.

The use of the variable gain amplifier allows the APB to act as successive iterative stages, as T changes value (recall that the MOS level shifter in each chained APB required different transistors to accommodate different values of 7).

Referring to Figures 1 7A and 1 7B, simulations suggest that such a feedback mechanism with a variable gain amplifier is stable. Figure 1 7A shows a single APB stage iterating using feedback, while Figure 1 7B shows a multiple APB stage feed forward version.

Both are shown to converge similarly on the same randomly chosen input values.

Advantageously, however, in the single-stage case the number of transistors used would only be approximately 7,000.

In an alternative embodiment, instead of a variable gain amplifier, the APB switches between different MOS level shifters according to a timing scheme, or according to the presence of a subsequent input voltage to the APB, so ratcheting through a set of fixed amplifications corresponding to each iteration of mk, and then outputting the result of the final level shifter.

It will be clear to a person skilled in the art that a combination of feed forward and feedback APBs may be used in assembling iterative stages. For example, a set of three feed forward stages may be used, wherein the output of the third stage feeds back to the input of the first, with variable gain amplifiers or switched level shifters set appropriately.

It will similarly be clear to a person skilled in the art that although the iterative stages described herein have been implemented using both bipolar (BJT) and CMOS devices, a fully CMOS implementation is possible by polarising the transistors in weak inversion where the voltage-current relationship is exponential as in the case of BJTs.

It will be clear to a person skilled in the art that currents for vector z may be copied in a similar fashion to R. It will be clear to a person skilled in the art that the analog equaliser described herein may comprise a discrete entity, for example an ASIC, or plurality of entities, for example separate analog processing blocks. Similarly it will be clear to a person skilled in the art that the equaliser may form part of an analog MIMO detector, or equaliser for a reader of a magnetic storage medium. A more general device may be adapted to incorporate the analog equaliser, such as an entertainment device for games or streaming media, a laptop or PDA, or a hard drive. Alternatively, the equaliser may be a functionally separable component such as in a plug-in circuit board, or a peripheral such as a PCMIA card.

It will be understood that the analog equaliser disclosed herein provides one or more of the following advantages: i. high speed detection for applications such as MIMO and disk reading; ii. innate computation of soft values; iii. iterative convergence toward an estimate of the MPEs; iv. reduced complexity (a sub-exponential relationship between transistors and receiver channels), and; v. a plurality of architectures to meet different component count restrictions.

Claims (23)

1. CLAIMS: 1. An analog equaliser including an estimation means, said

   estimation means comprising at least a first analog processing block (hereinafter APB) is arranged in operation to iteratively generate an estimate of marginal posterior expectations (hereinafter MPEs) for received bit values.
2. 2. An analog equaliser according to claim 1 wherein the estimation means further comprises scaling means to update MPEs according to a mean field annealing factor.
3. 3. An analog equaliser according to any one of the preceding claims wherein the estimation means further comprises coordinate descent minimisation means for obtaining the MPE estimates.
4. 4. An analog equaliser according to any one of the preceding claims compris ig at least two APBs operably connected in succession, and arranged in operation such that iteration is achieved by successive reestimation of the MPEs by succesive respective APBs).
5. 5. An analog equaliser according to claim 4 wherein each successive APB is arranged in operation to apply a lower mean field annealing factor than the preceding APB.
6. 6. An analog equaliser according to any one of the preceding claims and comprising a processing chain of APBs, wherein the outputs of the last in said processing chain of APBs are arranged in use so as to feed back to the inputs of the first APB in the chain, whilst a mean field annealing factor of each APB is lowered accordingly.
7. 7. An analog equaliser according to any one of claims 1 to 3 wherein the outputs of a single APB are arranged in use so as to feed back to its own inputs, whilst a mean

   field annealing factor is lowered accordingly.
8. 8. An analog equaliser according to any one of the preceding claims compris i circuitry arranged in operation to perform a calculation of the general form Ak - - b Ck,, tanh(4JJ given input values for A, B and c
9. 9. An analog equaliser according to claim 8 comprising circuitry operable to calculate 2a1 Lk =_-t\zk -aRktanh(__J k'!=k where L are log-likelihood ratios of the MPEs, T is an annealing factor, z is a signal model and R is a channel cross-correlation matrix, for k=1,. .,K bits.
10. 10. An analog equaliser according to any one of the preceding claims, comprising an analog processing block (APB) characterised by; a plurality of K sets of K-i Rk'.tanh calculation circuits, k=1,. . .,K, where K is the number of MPE estimates; means for summing each of the K sets of K-i Rk'.tanh calculation outputs; means for subtracting each said sum from a respective filtered signal Zk, k=l,...,K,and; means for scaling an output signal for each MPE estimate in inverse proportion

    to a mean field annealing factor T.
11. 11. An analog equaliser according to claim 10 wherein each tanh calculation circuit comprises a transconductance amplifier, an absolute value generator, and a Gilbert multiplier.
12. 12. An analog equaliser according to claim 10 wherein each tanh calculation cruit comprises an absolute value circuit constructed with PMOS transistors forming two translinear ioops, operably coupled to a transconductance amplifier.
13. 13. An analog equaliser according to claim 10 and comprising a plurality of APBs, wherein corresponding tanh calculation circuits within each APB share one absolute value circuit.
14. 14. An analog equaliser according to claim 10 further comprising a transresistance circuit operable to substantially linearly convert in input current to an output voltage.
15. 15. An ASIC comprising an analog equaliser in accordance with any one of the preceding claims.
16. 16. A multiple input, multiple output (MIMO) detector comprising an analog equaliser in accordance with any one of the preceding claims.
17. 17. A mobile communications device comprising an analog equaliser in accordance with any one of any one of the preceding claims.
18. 18. A mobile communications device according to claim 17 wherein the mobi1e communications device is any one of; i. a plug-in circuit board; ii. a PCMIA card; iii. aPDA; iv. a laptop computer, and; v. an entertainment device.
19. 19. A magnetic data storage device comprising an analog equaliser in accordance with any one of claims ito 15.
20. 20. A method of equalisation comprising the step of passing a plurality of log- likelihood marginal posterior expectations to an analog processing block (APB), and generating in said APB a revised estimate of the loglikelihood marginal posterioi expectations using coordinate descent optimisation.
21. 21. The method of claim 20 wherein subsequent re-estimations are performed by corresponding subsequent APBs.
22. 22. An analog equaliser substantially according to any of the illustrated embodiments of the invention, with reference to the accompanying drawings.
23. 23. A method of equalisation substantially according to any of the illustrated embodiments of the invention, with reference to the accompanying drawings.