Combined interference cancellation with FEC decoding for high spectral efficiency satellite communications

Abstract

A combined interference cancellation communication system employing forward error code (FEC) decoding for high spectral efficiency satellite communications. The disclosed system enables efficient utilization of available bandwidth through overlapping adjacent channels. A receiver is used to receive a waveform having data information and noise information. A filter bank is coupled with the receiver to receive and filter waveform and output channel information received by the receiver. The channel information received includes a combination of data signals and adjacent channel interference signals. A demodulator is provided to provide an estimation signal representative of an estimation of at least one parameter of the channel information. Soft-input and soft-output decoders are provided to receive the channel information in order to calculate and estimate interference values based on the estimation signal. Additionally, the soft-input and soft-output interference canceler is provided for receiving the output channel information and the estimated interference value calculated from the decoders in order to provide a data signal based on substantially without interference the channel information and estimated interference value. Thus, a substantially more accurate data signal is provided. Typically, such systems employ remote ground terminals, e.g., VSAT, which are used for communicating via a geosynchronous satellite from a remote location to a central hub station or other remote locations. A particular advantage of the disclosed systems is their relatively low site cost and small earth station size.

Description

BACKGROUND OF THE INVENTION

[0002] 1. Field of the Invention

[0003] The present invention relates generally to a system and method for performing noise cancellation in narrowband satellite communication systems. More particularly, the present invention relates to performing noise cancellation for estimated signal parameters on a demodulated signal.

[0004] 2. Description of the Related Art

[0005] Efficient use of available bandwidth in wireless, such as satellite, communications applications is a problem of paramount importance. An example of such a narrow band satellite includes very small aperture terminal (VSAT) systems. VSAT systems use compact earth stations that are installed at one or more customer's premises to provide links among the premises over a wide coverage area. Typically, in such systems, remote ground terminals are used for communicating via a geosynchronous satellite from a remote location to a central hub station or other remote locations. The central hub station communicates with multiple remote ground terminals. VSAT systems are used to handle customer network requirements, from small retail sites up to major regional offices, and can support two-way data, voice, multi-media, and other types of data. A particular advantage of these systems is their relatively low site cost and small earth-station size.

[0006] In wireless systems, multiple users share the same bandwidth. Channel sharing through fixed-allocation, demand assigned or random-allocation modes is known as multiple access. Two of the more commonly known basic multiple-access techniques include time division multiple access (TDMA) and code division multiple access (CDMA).

[0007] VSAT type systems have traditionally implemented TDMA using time division multiplexed (TDM) mode. Such systems generally are used for low speed (300 bps to 19,200 bps) data communications such as credit card processing and verification, point-of-sale inventory control and general business data connectivity. A typical TDM/TDMA network, when implemented in a star topology (FIG. 1), uses a large satellite hub system that manages all network terminal access and routing. Data is transmitted to and from the hub in short bursts on satellite channels that are shared with a number of other VSAT terminals. The hub communicates with these VSAT terminals over a higher speed outbound TDM satellite carrier. The terminals transmit back to the hub on assigned inbound carriers using TDM protocols. Such a combination enables a predetermined number of slots in time each second that each terminal can transmit data. In addition, more or less time can dynamically be assigned to the terminals based upon each terminal's individual requirements.

[0008] In contrast, in a CDMA type system a user's station signal is multiplied by a unique spreading code at a high speed to be spread in a wide frequency band. Thereafter, the signal is transmitted to a transmission path. In a receiving side, the signal that was multiplexed by the spreading code is subjected to a despreading process to detect a desired signal. Signal detection is based on a unique spreading code assigned to a user's station. If despreading is carried out with reference to a particular code used to spread a transmission signal, a user's station signal is correctly reproduced.

[0009] Regardless of the access technique used, increased efficiency and lower cost is a primary goal. Accordingly, efficiencies in bandwidth may be realized using techniques such as crowding of adjacent channels, frequency re-use, and increasing of data rates, generally resulting in an increased amount of data traveling through the limited amount of available bandwidth. Unfortunately, however, such techniques introduce a significant amount of interference which must be canceled. Combined multi-user detection and decoding is believed to have the potential to improve performance to match that of an interference-free system. However, most development heretofore has been in CDMA based systems. In particular, it is known that the optimum receiver for CDMA system employing Forward Error Control (FEC) coding combines the trellises of both the multi-user detector and the FEC code. However, the complexity of such a receiver is exponential in the product of the number of users and the constraint length of the code. This complexity makes the use of the optimal detector prohibitive for even small systems.

[0010] Copending U.S. Patent Application entitled “Combined Interference Cancellation with FEC Decoding for High Spectral Efficiency Satellite Communications” Ser. No. 09/436,670, filed Nov. 10, 1999 discloses a method and apparatus for a low complexity cancellation scheme in narrow band type satellite applications that allows for efficient utilization of available bandwidth by eliminating the interference resulting from the aggressive channel crowding.

[0011] However, there is a need for performing iterative noise cancellation on signal parameters such as frequency, timing phase and signal strength. The invention should preferably use a demodulator to achieve this result.

SUMMARY OF THE INVENTION

[0012] Briefly, the present invention relates to a satellite communications system and method for achieving efficient utilization of available bandwidth for satellite applications such as fixed wireless, mobile satellite systems and other narrow-band type applications. A soft decision-feedback scheme is used iteratively in combination with Forward Error Correction (FEC) decoding for interference cancellation to enable efficient use of the available bandwidth through aggressive crowding of adjacent channels.

[0013] In a first embodiment of the present invention, a multiple channel decoding receiver is provided. The multiple channel receiver includes a matched-filter bank that is used to receive signals and provide initial estimates of data. The estimates are fed to an interference canceler, which provides channel information to a demodulator. The demodulator outputs a signal representative of an estimation of at least one parameter of the channel information. The estimated signal is then provided to a decoder which calculates an estimated interference value based on the estimation signal. The estimated interference value comprises soft interference estimates which are fed back into the interference canceler. The soft estimates of the interfering signals are subtracted from the matched-filter outputs in the interference canceler to generate new, refined, soft approximations of the data. These refined soft estimates are then fed to the channel demodulators for the next iteration. The process is repeated iteratively until the channel interference reaches acceptable levels.

BRIEF DESCRIPTION OF THE DRAWINGS

[0014] The invention is pointed out with particularity in the appended claims. However, other features of the invention will become more apparent, and the invention will be best understood by referring to the following detailed description in conjunction with the accompanying drawings, in which:

[0015] FIG. 1 illustrates a VSAT system in a star topology;

[0016] FIG. 2 is a schematic diagram of an exemplary ACI transmitter model in accordance with an embodiment of the present invention;

[0017] FIG. 3 is a block diagram of a first exemplary embodiment of an ACI receiver model in accordance with the present invention;

[0018] FIG. 4 illustrates a spectral view of a signal model for known systems;

[0019] FIG. 5 illustrates a spectral view of a signal model in accordance with an embodiment of the present invention;

[0020] FIG. 6 is a block diagram of a multi-channel receiver which combines interference cancellation with forward error correction (FEC) decoding in accordance with the present invention;

[0021] FIG. 7 is a block diagram of a single-channel receiver which combines interference cancellation with forward error correction (FEC) decoding in accordance with the present invention;

[0022] FIGS. 8 a and 8b are graphs illustrating the performance of the proposed iterative multi-channel receiver at channel spacing of 0.75TS−1 with 4 states and 16 states convolutional code, respectively; and

[0023] FIG. 9 is a graph illustrating the performance of the proposed iterative single-channel receiver at channel spacing of 0.75TS1 with 4 states convolutional code.

DETAILED DESCRIPTION OF THE INVENTION

[0024] The present invention relates to a satellite communications system and method for achieving efficient utilization of available bandwidth for satellite applications. In particular, a soft decision-feedback scheme is used iteratively for interference cancellation in combination with FEC decoding to enable efficient use of the available bandwidth using techniques such as crowding of adjacent channels, frequency re-use, and increasing the data rates. A particular advantage of such a system is the ability to eliminate interference, such as adjacent channel interference (ACI), that may be introduced during, for example, channel crowding, thereby resulting in a higher spectral efficiency. For example, the present embodiment enables a satellite system to operate at a bandwidth efficiency level of 2.66 bits-per-second/Hz with minimum additional energy requirement in the signal-to-noise ratio range of interest using only a four state FEC code. This corresponds to an approximately 55% improvement in spectral utilization over current systems that employ similar modulation techniques. This improvement is expected to increase when more efficient FEC codes are used.

[0025] Turning now to FIG. 2, there is illustrated an exemplary VSAT system communicating in a TDMA mode. It is to be noted that the principals embodied in the present embodiment may also be applied to other satellite communications systems as well. As shown in the figure, the VSAT system, such as available from Hughes Network Systems, includes a central hub station 102 that controls one or more earth stations 104A-104B located on customers' premises. The earth stations 104A-104B and the central hub station 102 communicate with each other using a geosynchronous satellite 106. Each of the earth stations 104A-104B has a receiver 108A-108B for receiving and decoding signals received from the satellite 106 and transmitters 110A-110B for transmitting data to the satellite 106. The hub, or base station, station 102 similarly includes a receiver 112 for receiving and decoding signals received from the satellite 106 and a transmitter 114 for transmitting data to the satellite 106.

[0026] Turning to FIG. 3, an exemplary ACI transmitter model 120 is shown which may be used in the earth station transmitters 11OA-11OB and the hub transmitter 114. The transmitter 120 receives data from a first source 122A to an Mth source 122C. Converters 124A-124C convert the data from a binary phase shift keying (BPSK) signal to a quadrature phase shift keying (QPSK) signal. The resultant frequency domain pulse 126 is interleaved and transmitted as signal s(t) 128. The signal s(t) 128 which models the situation of adjacent channel interference caused by signal crowding, consists of the signal in noise as

r (t)=s(t)+n(t)   (1)

[0027] The n(t) is the standard additive white Gaussian noise (AWGN) with single-sided power spectral density level of N0 (Watts/Hz). The signals s(t) models the situation of ACI in which there are M adjacent data sources that are identical and independent. Each source transmits a QPSK signal at the rate of TS−1 with an arbitrary unit-energy pulse shape p(t). The signal is described in complex form as

s(t)=Re {{tilde over (s)}(t)ej2nfct}  (2)

[0028] where ƒc is the carrier frequency and {tilde over (s)}(t) is the baseband complex envelope of the signal and is mathematically expressed as 1$\begin{array}{cc}\tilde{s}\left(t\right)=\sum _{m=1}^M\sum _{k=-\infty }^{\infty }\left[a_{m,k}p\left(t-{\mathrm{kT}}_s-{\epsilon }_m{T}_s\right)e^{j\left(2n\Delta f_{m}t+\theta m\right)}\right]& \left(3\right)\end{array}$

[0029] The data streams {am,i; m=1, . . . , M} consist of QPSK symbols taking on the equi-likely values (±1±j) and are statistically related as 2$\begin{array}{cc}E\left\{a_{m,i}^{*}{a}_{n,j}\right\}=\{\begin{array}{cc}2,& m=n,i=j\\ 0,& \mathrm{otherwise}\end{array}& \left(4\right)\end{array}$

[0030] Assuming without loss of generality that M is an odd integer and that the desired center, or (M+1)/2-th, conveys the desired data and that the other signals are viewed as being adjacent-channel interferers ((M−1)/2 ones on either side), namely Δf1=<. . . <Δf(m+1)/2<. . . <ΔfM. The present embodiment is used when the channel spacing values, compared with that of the center channel, are small enough to cause large amounts of overlap in the spectra. Furthermore, in practical systems, these channels are equally spaced in frequency, say by Δf. In terms of the above, then 3$\begin{array}{cc}\Delta f_m=\left(m-\frac{M+1}{2}\right)\Delta f;m=1,2\dots ,M& \left(5\right)\end{array}$

[0031] It is to be noted that M should not be interpreted as being the number of channels in the entire available bandwidth. Instead, it is the number of channels that the receiver wishes to process jointly to announce a decision regarding the desired data stream. Although the use of guard bands is known when separating channels, their use consumes a non-trivial amount of bandwidth, thereby decreasing spectral efficiency. Therefore, in the absence of an installed guard band, the outermost or “edge” channels, i.e., 1st and Mth, will always have interference. However, as described further below, the present embodiment does not require that the receiver compensate for these edge channels. For example, M is chosen to be seven. As such, the receiver processes seven channels in the presence of two additional signals. Note that in the presence of guard band at the edges of the M channels, the receiver will jointly receive all the M channels. This feature is particularly useful for the base-station which is interested in decoding more than one channel.

[0032] Two criteria of performance are considered. The first is the bandwidth efficiency, η, in bits-per-second/Hz defined as the ratio of the data rate to the bandwidth used. This quantity is obtained in this case as 4$\begin{array}{cc}\eta \stackrel{\Delta }{=}\frac{2}{\Delta f{T}_s}& \left(6\right)\end{array}$

[0033] The second is the energy efficiency defined as the signal-to-noise ration per bit required to achieve a specific bit error probability Pb(E) of the desired channel. Other measures of performance may also be used, such as symbol error probability and word error probability.

[0034] It can be seen that the interference-free performance is 5$\begin{array}{cc}{P}_b\left(E\right)=Q\left(\sqrt{\frac{2E_b}{N_0}}\right)& \left(7\right)\end{array}$

[0035] where Eb is the average bit energy and Q(x) is defined as the Gaussian probability integral 6$\begin{array}{cc}Q\left(x\right)\stackrel{\Delta }{=}\underset{\chi }{\overset{\infty }{\int }}\frac{1}{\sqrt{2n}}{e}^{-\frac{y^2}{2}}dy& \left(8\right)\end{array}$

[0036] The performance described by the above equations is used as a benchmark to quantify the ability of the proposed receivers to suppress interference.

[0037] From statistical theory of hypothesis testing, it is known that the optimum receiver is the one that minimizes sequence error probability and is derived from implementing the average likelihood-ration function (ALF). The likelihood function conditional on knowing a given signal is AWGN is then obtained as 7$\begin{array}{cc}\bigwedge \left(r\left(t\right)\right)=\exp \left(-\frac{1}{N_0}\underset{-\infty }{\stackrel{\infty }{\int }}|\tilde{r}\left(t\right)-\tilde{s}\left(t\right){|}^{2}dt\right)& \left(9\right)\end{array}$

[0038] where r(t) is the baseband complex envelope of the received waveform. Thus, 8$\begin{array}{cc}\begin{array}{c}\ln \left(\bigwedge \left(\tilde{r}\left(t\right)\right)\right)=-\frac{1}{N_0}\underset{-\infty }{\stackrel{\infty }{\int }}|\tilde{r}\left(t\right)\left[\sum _{m=1}^M\sum _{k=-\infty }^{\infty }{a}_{m,k}p\left(t-{\mathrm{kT}}_s-{\epsilon }_m{T}_s\right)e^{j\left(2n\Delta f_{m}t+\theta m\right)}\right]{|}^{2}dt-\\ =\frac{1}{N_0}\underset{-\infty }{\stackrel{\infty }{\int }}|\tilde{r}\left(t\right){|}^{2}dt+\frac{1}{N_0}2\mathrm{Re}\left\{\underset{-\infty }{\stackrel{\infty }{\int }}{\tilde{r}\left(t\right)\left[\sum _{m=1}^M\sum _{k=-\infty }^{\infty }{a}_{m,k}p\left(t-{\mathrm{kT}}_s-{\epsilon }_m{T}_s\right)e^{j\left(2n\Delta f_{m}t+\theta m\right)}\right]}^{*}\right\}dt-\\ \frac{1}{N_0}\underset{-\infty }{\stackrel{\infty }{\int }}|\sum _{m=1}^M\sum _{k=-\infty }^{\infty }{a}_{m,k}p\left(t-{\mathrm{kT}}_s-{\epsilon }_m{T}_s\right)e^{j\left(2n\Delta f_{m}t+\theta m\right)}{|}^{2}dt\end{array}& \left(10\right)\end{array}$

[0039] By absorbing terms that are independent of the sought symbols and after some algebraic manipulations, maximizing the ALF over the information symbols is equivalent to maximizing the quantity J(a1,a2, . . . aM) of the metric, defined as 9$\begin{array}{cc}\begin{array}{c}j\left(a_1,a_2,\dots ,a_m\right)\stackrel{\Delta }{=}2Re\{\sum _{k=-\infty }^{\infty }\left[\sum _{m=1}^M{a}_{m,k}^{*}{x}_m\left(\left(k+{\epsilon }_m\right)T_s\right]\right\}-\\ \sum _{k=-\infty }^{\infty }\sum _{l=-\infty }^{\infty }\left[\sum _{m=1}^M\sum _{n=1}^M{a}_{m,k}{C}_{m,n}\left(\left(k+{\epsilon }_m{T}_{s^{\prime }}*l+{\epsilon }_n\right)T_s\right)a_{n,l}^{*}\right]\end{array}& \left(11\right)\end{array}$

[0040] The above expression shows that {xm((k+ε)Ts); m=1, . . . , M} is a set of sufficient statistics which consists of an exemplary bank of matched filters 132, as shown in FIGS. 4 and 5, matched to the modulating signal in each channel, then sampled at the symbol rate of Ts−1. More generally defined as 10$\begin{array}{cc}{x}_m\left(t\right)\stackrel{\Delta }{=}\underset{-\infty }{\stackrel{\infty }{\int }}r\left(a\right)e^{-j\left(2n\Delta f_{m}a+\theta m\right)}p*\left(a-t\right)da\mathrm{and}\mathrm{also}& \left(12\right)\\ C_{l,j}\left(t_1,t_2\right)\stackrel{\Delta }{=}\left[\int p*\left(a\right)p\left(a+t_2-t_1\right)e^{-\mathrm{j2n}(\Delta f_j-\Delta f_l)a}da\right]x\exp (-j\left(2n\left(\Delta f_j-\Delta f_l\right)t_2+\left({\theta }_j-{\theta }_l\right)\right)& \left(13\right)\end{array}$

[0041] The function of the optimal rule, or the maximum likelihood sequence estimation receiver 136, is to determine the sequence of information symbols (a1,a2, . . . aM) that maximizes the metric shown above. If there are N symbols in a frame, then the most straightforward way of implementing the optimum receiver requires 4MN computations of the metric. However, this procedure can be implemented in the most efficient way by generalizing the modified Viterbi Algorithm (VA) of G. Ungerboeck, “Adaptive Maximum-Likelihood Receiver for Carrier-Modulated Data Transmission Systems,” IEEE Transactions on Communications, pp. 624-636, May 1974.

[0042] Accordingly, the metric shown above can be made recursive by the relation 11$\begin{array}{cc}\begin{array}{c}{J}_p\left(a_1,p,a_2,p,\dots ,a_{M,p}\right)=J_{p-1}(a_1,p-1,a_2,p-1,\dots ,a_M,p-\\ 1)+\mathrm{Re}\{\sum _{m=1}^M{a}_{m,p}^{*}[2x_m\left(\left(p+{\epsilon }_m\right)T_s\right)-\\ \sum _{n=1}^M{C}_{n,m}\left(\left(p+{\epsilon }_m\right)T_s\right),\left(p+{\epsilon }_m\right)T_s)a_{n,p}-\\ 2\sum _{n=1}^M\sum _{k\le p-1}{C}_{n,m}(\left(k+{\epsilon }_n\right)T_{s^{\prime }}(p+\\ {\epsilon }_m)T_s)a_{n,k}]\}\end{array}& \left(14\right)\end{array}$

[0043] Equivalently, the interference channel whose impulse response spans L symbols can be viewed as a finite-state discrete-time machine where the state at discrete time i is defined as

Sl{double overscore (Δ)}(a1,l−1, . . . ,al−L;al,l−1, . . . ,a2,l−L; . . .; aM,l−1, . . . ,aM,l−L)   (15)

[0044] The VA then tracks the paths through the trellis and provides the solution to the problem of maximum-likelihood estimate of the state sequence. Thus, it is clear that the trellis has a maximum of 4ML states. Note that the efficiency of this modified VA stems from the fact that maximizing the likelihood function requires computing N4ML instead of 4MN metrics, wherein L is typically much smaller than N. Thus, even-though the complexity remains exponential in the number of channels, making the optimum rule computationally intensive, the complexity becomes independent of N. It is to be noted that reduced-complexity versions of the vector VA, which use decision-feedback on a per-survivor basis, may also be used. It is worth noting that this receiver does not include the FEC decoding which is done separately. This will result in a loss in performance. The reason for this is the huge complexity of the optimum joint receiver. In this invention we present a low complexity receiver for joint demodulation and channel decoding that achieves very close performance to the optimum receiver.

[0045] Before providing a suboptimal cancellation structures, the effective channel that appears at the output of the matched-filter bank, {xm(t);m=1,2,. . . ,M}, is characterized to yield 12$\begin{array}{cc}{x}_m\left(t\right)=\left[\sum _{n=1}^M\sum _{i=-\infty }^{\infty }{a}_{n,i}{C}_{n,m}\left(\left(i+{\epsilon }_n\right)T_{s^{\prime }}t\right)\right]+n_m\left(t\right)& \left(16\right)\end{array}$

[0046] From the above equation it is clear that the equivalent lowpass interference channel is described by the previously defined impulse response Cl,j(t1,t2), which can be equivalently represented in the frequency domain as 13$\begin{array}{cc}{C}_{l,j}\left(t_l,t_2\right)\stackrel{\Delta }{=}\left[\int P\left(f\right)P*\left(f-\left(\Delta f_j-\Delta f_l\right)\right)e^{\mathrm{j2}nf\left(t_2-t_l\right)}df\right]x\exp (-j\left(2n\left(\Delta f_j-\Delta f_l\right)t_2+\left({\theta }_j-{\theta }_l\right)\right)& \left(17\right)\end{array}$

[0047] This function represents the effective channel impulse response at the output of the jth matched-filter when excited by the Ith data source. It consists of the cascade of the pulse-shaping filter and the complex multiplier at the transmitter side, the channel, and the matched-filter at the receiver. It is to be noted that the impulse response in this case is time-varying, a condition that results from the presence of complex-exponential multipliers (or frequency shifters) in the system. As the channel spacing is increased, the magnitude of the impulse response decreases but its duration is increased, resulting in an equivalent channel with larger memory span. From above,

nm(t){double overscore (Δ)}∫ñ(a)p*(a−t)e−j(2nΔƒma+θm)da   (18)

[0048] wherein, the {nm(t);m=1,2,. . . ,M} is a set of zero-mean complex Gaussian random processes with covariance

E{n*l(t)nj(t′)}=NOClj(t,t′)   (19)

[0049] Several assumptions are made to simplify notation. First, in regards to the pulse shaping, it is assumed that the pulse p(t) selected satisfies the Nyquist criterion of zero inter-symbol interference. This criterion is expressed in time as 14$\begin{array}{cc}\stackrel{\infty }{\underset{-\infty }{\int }}p*\left(t\right)p\left(t+{\mathrm{nT}}_s\right)dt=\{\begin{array}{cc}1,& n=0\\ 0,& \mathrm{otherwise}\end{array}& \left(20\right)\end{array}$

[0050] or in frequency as 15$\begin{array}{cc}\frac{1}{T_s}\sum _{n=-\infty }^{\infty }|P\left(f+\frac{n}{T_s}\right){|}^2=1& \left(21\right)\end{array}$

[0051] This states that the aliased or folded version of the auto-power spectrum associated with the selected pulse must be flat. (The aliased version is what results when replicating the function at multiples of the symbol rate Ts−1.) Note that the root-raised cosine pulse, which is a practical and bandwidth-efficient choice, is defined in the frequency domain as 16$\begin{array}{cc}P\left(f\right)=\{\begin{array}{cc}1,& 0\le \left|f\right|\le \left(1-\beta \right)/2T_s\\ \sqrt{\frac{1}{2}\left[1-\sin n{T}_s\left(f-\frac{1}{2T_s}\right)/\beta \right]},& \left(1-\beta \right)/2T_2\le \left(1+\beta \right)/2T_s\\ 0,& \mathrm{otherwise}\end{array}& \left(22\right)\end{array}$

[0052] where β is the roll-off parameter. Second, the spectral overlap of these channels does not exceed 50%. This, along with practical values of the roll-off parameter, indicates that the ACI on a given channel results from one adjacent interferer on either side. Third, it is assumed that the ACI extends over a finite time interval spanning L symbols. The actual value of L is directly related to the amount of spectral overlap that exists between the channels. From basic principles of Fourier transforms, the value of L, which can be thought of as the memory of the interference channel, is larger for smaller overlap.

[0053] Fourth, it is assumed that the receiver is able to maintain phase coherence and time synchronism. For the time synchronism situation, a situation representative of a satellite down-link application, it is assumed that the relative time delays are zero. Based on the above set of assumptions, the matched-filter statistic is described as 17$\begin{array}{cc}\begin{array}{c}{x}_m\left({\mathrm{kT}}_s\right)=a_{m,k}+\sum _{i=-L}^L{a}_{m-1,k-i}{C}_{m-1,m}\left(\left(k-i\right)T_{s^{\prime }}{\mathrm{kT}}_s\right)+\\ \sum _{i=-L}^L{a}_{m+1,k-i}{C}_{m+1,m}\left(\left(k-i\right)T_{s^{\prime }}{\mathrm{kT}}_s\right)+n_m\left({\mathrm{kT}}_s\right)\end{array}& \left(23\right)\end{array}$

[0054] The first term on the right-hand side of the above equation is the desired information symbol; the second term is the ACI contribution from the left channel; while the third term is the ACI contribution from the right channel. The ACI is determined by the symbol-spaced samples of the cross-correlation between transmit and receive filters. The set {nm(kTs);m=1,2,. . . ,M} has elements that are statistically correlated across different m's but independent for a specific m.

[0055] Combined multi-user detection and decoding has received considerable attention recently with its potential to improve performance to match that of the single-user system. However, most of the work done has focused on spread-spectrum CDMA systems. It has been shown that the optimum receiver, for a CDMA system employing FEC coding, combines the trellises of both the multi-user detector and the FEC code. The complexity of this receiver is exponential in the product of the number of users and the constraint length of the code, making the implementation of the optimal detector prohibitive for even small systems.

[0056] FIG. 4 which illustrates a spectral view of a signal model for known systems highlights the inefficiencies of the current art. As previously discussed above, the bandwidth is limited and an ideal scenario is to pack as many channels as possible into the limited available bandwidth. However, limitations such as interference between channels limit the number of channels that can be packed into the available bandwidth.

[0057] FIG. 5 illustrates a spectral view of a signal model in accordance with an embodiment of the present invention. Specifically, FIG. 5 depicts a 40% improvement in bandwidth efficiency over prior art systems. Channels are packed closer together in the limited available bandwidth resulting in greater bandwidth capacity. By decoding additional channels, an even greater improvement can be achieved.

[0058] Recently, a new powerful class of concatenated convolutional codes was proposed which use parallel concatenation of two (or more) recursive systematic convolutional codes (constituent codes) fed by two information sequences, of which the second is obtained from the first through the interposition of a long interleaver. One of the key factors contributing to the remarkable performance of this coding scheme is the elegant iterative Soft-Input Soft-Output (SISO) decoding structure whose performance was shown, via simulation, to approach that of the maximum-likelihood decoding, at signal-to-noise ratios very close to the Shannon limit, with much less complexity. This decoder is based on iteratively decoding the component codes and passing the so-called extrinsic information, which is a part of the component decoder soft output, to the next decoding stage. The impressive performance, achieved by this iterative decoding architecture, has encouraged several researchers to consider applying the same principle in the other sub-modules of the receivers.

[0059] In this section, several strategies for joint iterative interference cancellation and FEC decoding, which provide a tradeoff between performance and complexity, are investigated. We will restrict ourselves to binary convolutional codes; the extension to higher order codes is straightforward. The proposed algorithms function independently of the type of channel decoders in the sense that any SISO decoders can be used (for example, MAP, Log-MAP, or SOVA). The choice of a particular decoder is determined by the allowable complexity at the receiver. Based on the Gray mapping assumption, each QPSK symbol is

a m,k =a m,k I +ja m,k Q   (24)

[0060] where am,kI, and am,kQ ε{−1,1}. Without loss of generality, we will derive the different detection rules for the in-phase k-th binary symbol associated with the m-th channel.

[0061] Denote,

A+{double overscore (Δ)}[am−1,k−L, . . . ,am−1,k+, +1,a m,kQ,am+1,k−L, . . . , am+1,k+L:aj,l ε{ejπi/4},am,kQ ε{−1,1}]  (25)

[0062] The set A− can be similarly defined. Using the independence assumption, which is justified by the interleaving, the maximum-a-posteriori (MAP) detector is given by 18$\begin{array}{cc}{L}_{k,m}^I=\log {\sum }_{\underset{\_}{a}\in A^{+}}\frac{P\left({\chi }_{m,k}|\underset{\_}{a}\right){\prod }_{j\in \left\{m-1,m+1\right\},k-L\le l\le k+L}P\left(a_{j,l}\right)P\left(a_{k,m}^Q\right)}{{\sum }_{\underset{\_}{a}\in A^{-}}P\left({\chi }_{m,k}|\underset{\_}{a}\right){\prod }_{j\in \left\{m-1,m+1\right\},k-L\le l\le k+L}P\left(a_{j,l}\right)P\left(a_{k,m}^Q\right)}& \left(26\right)\end{array}$

[0063] where Lk,mI is the updated log-likelihood ratio, P(χm,k|a) is the conditional Gaussian distribution of the matched filter output as per (26). In (26), P(aj,l) and P(ak,mQ) are obtained from the soft outputs of the previous iteration as follows

P (aj,l=aj,lI+jaj,lQ)=P(aj,lI)P(aj,lQ)   (27)

[0064] 19 $$ P ⁡ ( a j , l I = 1 ) = e λ j , l I 1 + e λ j , l I ( 28 ) P ⁡ ( a j , l I = - 1 ) = 1 1 + e λ j , l I ( 29 )

[0065] where λj,lI is the output log-likelihood ratio of the previous iteration. It is clear that results similar to (28), (29) hold for P(aj,lQ).

[0066] The MAP detector requires a complexity of the order O(2×42(2L+1)) which can be prohibitive for practical applications. Therefore, in the following, a lower complexity detection rule based on the MMSE principle is developed. The iterative MAP detection rule has been proposed for CDMA signals.

[0067] This scheme, depicted in FIG. 6, uses the soft information supplied by the M single-user decoders to calculate the optimum, feed-forward and feed-back, filter weights after each iteration. FIG. 6 shows the multi-channel receiver 130 which combines interference cancellation with forward error cancellation (FEC) decoding using a matched filter 132 and channel estimation 136 as inputs to an interference canceler (IC) 134. Multiple channels M are provided from the interference canceler 134 to demodulators 138a . . . 138b which provide outputs to optional equalizers 140a . . . 140b. The equalizers 140a . . . 140b provide outputs to deinterleaver 142a . . . 142b, which provide outputs to SISO decoders 144a . . . 144b and interleavers 146a . . . 146b for outputs 1 through M. The interleaver output 146a and 146b are fed back to the interference canceler 134 in a closed loop.

[0068] The received signal is processed at the match filter 132 where interference associated with the received signal is suppressed or the transmitted pulse of the signal is matched. Specifically, the out of band interference is suppressed. The suppressed signal is then provided to the IC 134 where a portion of the interference is canceled.

[0069] The suppressed signal is provided to the demodulators 138a through 138b, where timing, phase, frequency and signal strength estimation is performed on the signal. The estimated parameters are used by the IC 134 to regenerate the interference signal. The interference suppressed signal is used at the successive iteration by the demodulators 138a, 138b to re-estimate the parameters again. With each iteration, the demodulators 138a, 138b encounter less interference.

[0070] The interference can not practically be eliminated completely. Each iteration reduces the amount of interference. However, there are trade offs. For each iteration, the reduction of the amount of interference becomes less. In addition, different parameters can have different convergence rates. For example, the amplitude convergence can gain satisfactory accuracy in the first few iterations while the phase and timing estimation can require further iterations. To accommodate this situation the demodulators 138a, 138b functions can be terminated earlier than that of the decoders 144a, 144b and IC's 134 functions. In other words, the estimation is discontinued at the demodulators 138a, 138b for the parameters that converge.

[0071] In order to limit the number of iterations performed on adjacent channels, in one embodiment of the present invention, the multiple channels M are grouped into carrier groups. Within each carrier group, interference cancellation techniques are performed to reduce the spacing between channels within the carrier groups. However, sufficient spacing is utilized between carrier groups so that interference cancellation is not required. In other words, sufficient spacing is provided between the carrier groups so that minimal interference occurs between carrier groups.

[0072] If the multiple channels M were not grouped into carrier groups, the possibility of continuous interference cancellation could occur. That is, in large systems having a number of channels, multiple occurrences of interference could occur. The interference cancellation process would continue the process for each channel. By having the multiple channels M grouped into carrier groups, the number of iterations is limited by the carrier group size.

[0073] In another embodiment of the present invention, each carrier group is further divided into two subgroups—an even and odd group. Demodulation/decoding can be performed first in one subgroup. When the receiver processes the second sub group, the interference information gained from processing the first subgroup can be utilized for cancelling out part of the interference in the second subgroup. In this manner, the second subgroup can benefit from an extra iteration without having to perform that extra iteration.

[0074] Equalizers 140a, 140b are optional and equalize the received signal from demodulators 138a, 138b. In one embodiment of the invention, equalizers 140a, 140b can be used where there is severe inter-symbol channel interference. The equalizers 140a, 140b are depicted as being connected between demodulators 138a, 138b and decoders 144a, 144b. However, the equalizers 144a, 144b can be inserted after and/or in parallel with the demodulators 138a, 138b and can be included with the IC 134. Thus, the equalizers 140a, 140b can benefit from the reduced interference level and the reliability information from the soft output decoder 144a, 144b.

[0075] However, it is not just the equalizers 140a, 140b which benefit from the reduced interference level and the reliability information from the soft output of the decoders 144a, 144b, the demodulators 138a, 138b also benefit. The decoders 144a, 144b provide reliable hard decisions for the tracking of loops for the demodulators 138, 138b. For example, the decoders 144a, 144b provide information on what's being transmitted e.g., a zero or one. Further, the decoders 144a, 144b also provides soft information on the probability of the information being correct e.g., the information being transmitted has a probability of 0.7 of being a one. This provides reliability to demodulators 138a, 138b.

[0076] Although the present invention is discussed in terms of linear systems, it will be appreciated by those skilled in the art that the present invention is applicable to nonlinear systems. The present invention can be applied to any FEC based transmitter system where the transmitter characteristics are known or can be estimated by the receiver.

[0077] FIG. 7 on the other hand shows the single-channel embodiment of a receiver 148 which combines interference cancellation and FEC decoding employing the match filter bank 132 with a minimum means square error (MMSE) transversal filter 150 which provides a single channel to a deinterleaver 154, SISO decoder 156, and interleaver 158 which provides input to calculate feed forward and feed back coefficients at 160. The feed forward and feed back coefficients calculation 160 also receives the channel estimation 152 to provide feed forward and feed back coefficients to the MMSE transversal filter 150. For this we extend the algorithms for CDMA signals to the current narrow band TDMA application. Let χ be a [M(2L+1)×1] complex vector of the matched filter bank outputs from the k−L)-th to the (k+L)-th samples. Then χ can be written as 20$\begin{array}{cc}\underset{\_}{\chi }=\sum _{j=l}^M\sum _{l=k-2L}^{k+2L}{\underset{\_}{C}}_{j,l}{a}_{j,l}+\underset{\_}{n}& \left(30\right)\end{array}$

[0078] where, for example,

C j,l ={square root}{square root over (Eb,j)}[Cj,l((l+εj)Ts,(k−L+εl)Ts), . . . , Cj,M((l+εj)Ts,(k+L+εM)Ts)]T

[0079] is the ACI associated with aj,l, and n is the Gaussian noise vector. Let, RI, the interference correlation matrix, be defined as

RI{double overscore (Δ)}[Cl,k−2L. . . Cm,k−1Cm,k+1. . . CM,k+2L]T   (31)

[0080] and

a {double overscore (Δ)}[a1,k−2L, . . . ,am,k−1,a m,k+1, . . . ,aM,k+2L]T   (32)

[0081] then

χ =Cm,kam,k+RIa+n  (33)

[0082] Now, the updated decoder input is calculated from

ym,k=cƒTχ+cb   (34)

Lm,kI=Re{ym,k}  (35)

Lm,kQ=Im{ym,k}  (36)

[0083] where cƒ is the [M(2L+1)×1] feed-forward coefficients vector, cb is the feed-back coefficient. Restricting the filter to have a single feed-back coefficient, rather than a vector, should not result in a loss of degrees of freedom. The coefficients cƒ, cb are obtained through minimizing the MSE between the data symbol and its estimate, given by 21$\begin{array}{cc}\begin{array}{c}\mathrm{MSE}=E\left[|y_{m,k}-a_{m,k}{|}^2\right]\\ E\left[|{\underset{\_}{c}}_f^T\underset{\_}{X}+c_b-a_{m,k}{|}^2\right]\\ E\left[|{\underset{\_}{c}}_f^T\left({\underset{\_}{C}}_{m,k}{a}_{m,k}+R^I\underset{\_}{a}+\underset{\_}{n}\right)=c_b-a_{m,k}{|}^2\right]\end{array}& \left(37\right)\end{array}$

[0084] It may be shown that the optimum solution based on MMSE must satisfy the following relations

E[c ƒ T (Cm,kam,k+RIa+n)+cb]=0   (38)

E[ (cƒT(Cm,kam,k+RIa+n)+cb−am,k)χH]=0   (39)

[0085] Note that the relation (38) ensures that the output of the MMSE filter is unbiased, while (39) is a direct application of the orthogonality principle. Solving (38), (39), we obtain the following results for the feed-forward and feed-back filter coefficients

c b =−c ƒ T R I E[a]   (40)

c ƒ T =C m,k H (A+B+Rn−CCH)−1   (41)

[0086] where, by definition, we have

A{double overscore (Δ)}Cm,kCm,kH   (42)

B{double overscore (Δ)}RIE[aaH]RIH   (43)

C{double overscore (Δ)}RI E[a]  (44)

[0087] and Rn is the [M(2L+1)×M(2L+1)] noise covariance matrix which may be constructed using a component-wise relation. In (43)-(44), the E[a] and E[aaH] values are obtained from the following component-wise relations

E (aj,l)=E(aj,lI)+jE(aj,lQ)   (45)

[0088] 22 $$ E ⁡ ( a j , l I ) = ⅇ λ j , l I - 1 ⅇ λ j , l I + 1 ( 46 )

E(aj,laj,l*)=1   (47)

E(aj,lai,n)=E(aj,l)E(ai,j)   (48)

[0089] where (48) follows from the independence assumption. In the first decoding iteration, we select E[a]=0. The feed-forward filter coefficients vector, cƒ, in this iteration is given by the MMSE equations derived in (41), and the feedback coefficient cb=0. After each iteration, the E[a] values are re-calculated using the decoder's soft outputs. The E[a] values are then used to generate the new set of filter coefficients as described. In the asymptotic case where |E[a]|=1, the receiver is equivalent to the subtractive interference canceler. This is expected, since |E[a]|=1 means that previous iteration decisions, for the interference, are error-free. In this case, the subtractive interference canceler is the optimum solution. The two main sources of complexity in this algorithm are the matrix inversion operation required in (41), and the need for M SISO channel decoders. In the following, we will investigate how to lower the computational complexity by proposing a soft interference cancellation algorithm that does not require a matrix inversion operation. In addition, the performance of the MMSE SISO detector will be studied assuming single-channel decoding.

[0090] Based on (40), (41), ym,k can be written as

y m,k =c ƒ T (χ−RIE[a])   (49)

[0091] By observing that the matrix inversion operation is only required to compute cƒT, the following approximation is proposed

c ƒ T =e (m−1)(2L+1)+L+1 T   (50)

[0092] where em−1)(2L+1)+L+1T=[0, . . . 0,1(m−1)(2L+1)+L+1,0, . . . , 0]. Hence,

y m,k =χm,k−(RIE[a])m,k   (51)

[0093] The complexity of this algorithm is a linear function of the product of the number of interfering users and the interference memory (i.e., O(2(L+1))). This algorithm can be regarded as a soft subtractive interference canceler. This is so as the decoder's soft outputs are used to calculate estimates of the transmitted symbols, E[a]; the estimates of the transmitted symbols and the interference cross-correlation matrix, RI, are used to generate updated estimates of the interference signals, at the output of the matched filter. The interference estimates are then subtracted from the matched filter output xm,k resulting in the next decoding iteration input. It is interesting to note that this MMSE-based development results in a scheme similar to the one proposed under a different derivation for CDMA signals.

[0094] Although the iterative SISO MMSE algorithm was developed assuming the use of M channel decoders, it is straightforward to modify the algorithm to be used with M matched filters but with an arbitrary number of decoders that is smaller than M. For any undecoded user d, we have

E[ad,l]=0 for any l   (52)

[0095] after any iteration. It is particularly interesting to consider the case of single-channel decoding which is displayed in FIG. 7. In this case, the soft information is fed-back to cancel the ISI, appearing at the matched filter bank outputs, and the feed-forward MMSE filter coefficients are used to suppress the ACI. It is easy to see that, assuming error-free feedback, this algorithm is capable of suppressing the ACI, except for the edge channels effect, asymptotically.

[0096] Monte-Carlo simulations are implemented to evaluate the bit error rate performance and demonstrate the effectiveness of the proposed solutions that combine interference cancellation and FEC decoding. The FEC encoding considered is the optimum rate ½ convolutional code with 4 and 16 states and the decoder uses the soft output Viterbi algorithm (SOVA). The channel spacing used is Δƒ=0.75Ts−1 (Hz) which corresponds to a spectral efficiency level of 2.6671 bps/Hz or an improvement of 55% compared to the current state-of-the-art.

[0097] FIGS. 8 a and 8b compare the performance of the soft interference cancellation scheme, the conventional receiver, and the interference-free system. The receiver processes seven channels jointly in a presence of a total of nine QPSK sources. It is clear that the performance of the proposed iterative decoding and interference cancellation algorithm is better than the conventional receiver and very close to the interference-free system, with a difference of less than 0.5 dB using four iterations when the input SNR is about 4 dB. It is also noted that as the SNR increases, the performance gap between the iterative algorithm and the conventional receiver increases while it diminishes more compared with the interference-free system.

[0098] FIG. 9 includes the performance achieved by the iterative MMSE algorithm, assuming single-channel coding, the conventional receiver, the feed-forward MMSE receiver, and the interference-free system. The receiver has three matched filters and it was assumed that only three channels are transmitting simultaneously (i.e., neglecting the edge effect). It is clear that the proposed algorithm provides considerable gain in performance compared to the conventional receiver and the feed-forward MMSE receiver. However, the difference in performance between the interference-free bound and the single-channel decoding algorithm is between 1.5-2 dBs which may be unacceptable in some cases. The performance can be improved upon by decoding more channels, offering a tradeoff between performance quality and receiver computational load.

[0099] The present invention and its performance have been described primarily in association with the customary AWGN channel. It is however clear from concepts in this patent to extend it to the case of fading channel. This can be done, for example, by estimating the channel directly and incorporating this estimate into the reconstruction of the interference for subsequent cancellation.

[0100] In order to prevent the synchronizing sequence from being obscured by the adjacent channels, a system level mechanism is needed to ensure that only one user is transmitting during its known symbol period.

[0101] Obviously, many modifications and variations of the present invention are possible in light of the above teachings. Thus, it is to be understood that, within the scope of the appended claims, the invention may be practiced otherwise than as specifically described above.

Claims

1. A system for enabling efficient utilization of available bandwidth through overlapping adjacent channels comprising: a receiver, for receiving a waveform having data information and noise information, a filter bank, adapted to receive and filter said waveform and output channel information, said channel information including a combination of data signals and adjacent channel interference signals; at least one demodulator, adapted to output an estimation signal representative of an estimation of at least one parameter of said channel information; at least one decoder, adapted to calculate an estimated interference value based on said estimation signal; and an interference canceler, adapted to estimate a data signal substantially without interference based on said output channel information and said estimated interference value.

2. The system of claim 1 further comprising: one or more equalizers adapted to equalize the estimation signals from said at least one demodulator.

3. The system of claim 2 wherein said equalizers are in parallel with said at least one demodulator.

4. The system of claim 2 wherein said equalizers are in series with said at least one demodulator.

5. The system of claim 1 wherein said signal parameter comprises at least one of a frequency parameter for determining a frequency value; a timing parameter for determining a timing value; a phase parameter for determining a phase value; and a signal strength parameter for determining a signal strength value.

6. The system of claim 1 wherein said channels comprise carrier groups.

7. The system of claim 6 wherein said carrier groups comprise odd channels.

8. The system of claim 6 wherein said carrier groups comprise even channels.

9. The system of claim I wherein said relatively more accurate estimated data signal is fed back into said interference canceler for a predetermined number of iterations.

10. The system of claim 1 wherein said interference canceler is designed based on the minimum means square error criterion (MMSE).

11. The system of claim 1 wherein said interference canceler is equipped with feed-back coefficients to subtract the estimated interference and feed-forward coefficients to suppress the residual interference.

12. The system of claim 11 wherein the feed-forward and feed-back coefficients of the interference canceler are optimized in every iteration using the feed-back information, from said at least one decoder.

13. The system of claim 12 wherein, in another embodiment, the same feed-forward coefficients (matched filter coefficients) are used in all iterations to reduce the complexity involved in the optimization process.

14. The system of claim 1 wherein, in another embodiment, the interference canceler is designed using the maximum-a-posteriori (MAP) rule.

15. The system of claim 1 wherein said at least one decoder provides soft information to said interference canceler.

16. The system of claim 1 wherein said at least one decoder provides hard information to said interference canceler.

17. A method for enabling efficient utilization of available bandwidth through overlapping adjacent channels comprising: receiving a waveform having data information and noise information, receiving and filtering said waveform and output channel information, said channel information including a combination of data signals and adjacent channel interference signals; estimating an output signal representative of an estimation of at least one parameter of said channel information via at least one demodulator adapted to receive said channel information; calculating an estimated interference value based on said output estimation signal via at least one decoder adapted to receive said output estimated signal; and estimating a data signal substantially without interference based on said channel information and said estimated interference value via an interference canceler adapted to receive said signals to produce a more accurate data signal.

18. The method of claim 17, wherein soft-input/soft-output decoders are used to obtain the estimates of the data.

19. The method of claim 17, wherein soft-input/hard-output decoders are used to obtain the estimates of the data.

20. The method of claim 18, wherein the soft-input/soft-output decoders comprise at least one of a Maximum a-posteriori (MAP) algorithm; a Log-MAP algorithm; and a Soft-output Viterbi (SOVA) algorithm.

21. The method of claim 17, wherein said channels comprise carrier groups.

22. The method of claim 21 wherein said carrier groups comprise odd channels.

23. The method of claim 21 wherein said carrier groups comprise even channels.

24. The method of claim 18 further including the step of calculating subsequent estimated interference signals from said relatively more accurate data signals.

25. The method of claim 24 further including the step of feeding back said relatively more accurate data signals into said calculating step a predetermined number of times.

26. The method of claim 24 further including the step of repeating said feed back step to output increasingly accurate estimated data signals.

27. The method of claim 21 further including the step of: equalizing said channel information with one or more equalizers adapted to receive said channel information from said demodulators.

28. The method of claim 17 wherein said signal parameter comprises at least one of a frequency parameter for determining a frequency value; a timing parameter for determining a timing value; a phase parameter for determining a phase value; and a signal strength parameter for determining a signal strength value.