Field of the Invention
This invention relates in general to methods for digital filter bank analysis and synthesis and in particular to the methods for multiplexing and demultiplexing of frequency division multiple access (FDMA) signals.
Background of the Invention
Digital filter bank and spectrum analysis and synthesis are extremely important in a variety of fields, including voice and data communication. The IRIDIUM.TM. Satellite Cellular Communications System, for example, calls for use of an L-band modem link using a 24-channel frequency division multiple access (FDMA) signal. The FDMA signal must be channelized to provide for proper interfacing to single-channel modulators and demodulators. At the transmitter, the baseband signals must be filtered and translated to the appropriate frequencies for summation into an overall FDMA signal (i.e., multiplexed). In the receiver, the various channels of the FDMA input signal must be separated by filtering and translated to baseband (i.e., demultiplexed) for input to a bank of demodulators.
Taking a brute force approach to generating and separating FDMA signals requires a channel filter and local oscillator for each FDMA channel. This approach, however, results in a large size and weight for a large number of channels, which is intolerable for a space-borne system.
Various techniques have been described for digital filter bank and spectrum analysis and synthesis. The text Multirate Digital Signal Processing, by Ronald E. Crochiere and Lawrence R. Rabiner, published in 1983 by Prentice-Hall, Inc., Englewood Cliffs, N.J., develops such techniques for several cases, including a complex filter bank and single-sideband filter bank. The Crochiere and Rabiner technique is not applied to teach a hardware-efficient method for generating and separating FDMA signals, however, such as are required for the IRIDIUM.TM. System.
Thus, what is needed is a method for FDMA signal multiplexing and demultiplexing which reduces filter computation and storage requirements, frequency translation computations, and minimizes size and weight. Specifically, the method should allow the transmission and reception of real signals while manipulating complex baseband signals (I and Q data of a quadrature phase shift key (QPSK) signal).
Summary of the Invention
Accordingly, it is an advantage of the present invention to provide a new and improved method for FDMA signal multiplexing and demultiplexing. It is a further advantage of the present invention to provide such a method to reduce filter computation, storage requirements, and frequency translation computations. It is still a further advantage of the present invention to be conductive to a hardware implementation to minimize size and weight.
To achieve these advantages, a method for multiplexing a series of frequency division multiple access (FDMA) signals to produce a multiplexed output signal is contemplated which includes the steps of providing a plurality of complex channels corresponding to the series of FDMA signals. Each of the series of FDMA signals is offset relative to each of the plurality of complex channels. Each of the series of FDMA signals is filtered to be centered within and occupy one-half of an available bandwidth of each of the plurality of complex channels, leaving a series of "off" channels between each pair of successive FDMA signals. The real part of the multiplexed output signal is processed such that a plurality of higher frequency signals of the series of FDMA signals image into the series of "off" channels.
The method for demultiplexing a multiplexed input signal into a sequence of frequency division multiple access (FDMA) signals includes the step of providing a demultiplexing channelization pattern comprising a plurality of complex channels in a complex channel spectrum. The multiplexed input signal contained in a plurality of input channels is aligned with the demultiplexing channelization pattern, wherein each of the plurality of input channels are one-half a bandwidth of each of the plurality of complex channels. The plurality of input channels contain a sequence of FDMA signals separated by a reverse sequence of images of the FDMA signals. The plurality of input channels is overlaid onto the plurality of complex channels such that consecutive even-numbered input channels correspond to consecutive complex channels of the plurality of complex channels and the sequence of FDMA signals can be extracted from the plurality of complex channels.
The above and other features and advantages of the present invention will be better understood from the following detailed description taken in conjunction with the accompanying drawings.
Brief Description of the Drawings
In FIG. 1, there is shown a diagram of a complex output spectrum from an FDMA synthesizer (multiplexer).
In FIG. 2, there is shown a diagram of a real output spectrum from an FDMA multiplexer.
In FIG. 3, there is shown a diagram of analyzer (demultiplexer) channel extractions from a real input spectrum.
In FIG. 4, there is shown a schematic diagram of a Generalized Discrete Fourier Transform (GDFT) polyphase multiplexer which operates in accordance with a preferred method of the invention.
In FIG. 5, there is shown a schematic diagram of a GDFT polyphase demultiplexer which operates in accordance with a preferred method of the invention.
Detailed Description of the Prior Art
The text Multirate Digital Signal Processing, by Ronald E. Crochiere and Lawrence R. Rabiner, referred to above and herein incorporated by reference, provides a useful framework and notation for mathematically describing the present method for multiplexing and demultiplexing.
Standard decimator and interpolator structures, using a direct structure and employing a low pass filter, are not efficient because a portion of the decimator output values are ignored and a portion of the interpolator calculations involve multiplication by zero. A more efficient method uses a polyphase network structure in which an original filter response has been apportioned equally across branches to create various branch filters. That is, where N is the total number of coefficients in the prototype filter and M is the interpolation or decimation factor, each branch has N/M coefficients.
The equations describing the input and output relationships of a conventional multiplexer/demultiplexer operation are: ##EQU1## where x is the input signal, n is dependent on the sampling process, K is the number of channels, k=0,1, . . . , K-1 and denotes the signal channels, X.sub.k denotes the kth channel signal, f is a synthesis (multiplexing) filter function, M is a sampling rate reduction factor, and W.sub.k is e.sup.2.pi.j/k.
The polyphase structure described above is described by the following equations: ##EQU2## where .rho.=0, 1, 2, . . . , M-1
The overbar notation indicates a clockwise commutator implementation versus a counterclockwise implementation. x.sub..rho. with the "hat" represents the polyphase output signals, q.sub..rho. and p.sub..rho. are the polyphase branches for the bandpass filters, X.sub.k are the input signals in the synthesizer (i.e., multiplexer) and W.sub.M =e.sup.2.pi.j/M.
Both equations 3 and 6 include an operation in the form of a Discrete Fourier Transform (DFT). Examination of equations 3 and 6 also reveals that the total number of filter calculations required for the polyphase multiplexer/demultiplexer has been reduced by M.
Detailed Description of the Preferred Embodiment
The actual implementation of the polyphase multiplexing and demultiplexing method for the IRIDIUM.TM. System involves an FDMA signal format of 32 channels. Only the middle 24 channels are used, however, due to Intermediate Frequency (IF) filter attenuation present in the analog portion of the transceiver. In the multiplexer, it was desired to generate 32 real channels within the allowable bandwidth (0 to .pi., normalized, for a real signal).
The multiplexing/demultiplexing method described assumes complex channels, in which the normalized bandwidth encompasses 0 to 2.pi.. But when the real part of the signal is extracted for transmission, the frequencies from .pi. to 2.pi. will alias into the band from 0 to .pi.. Therefore, it was required to generate 32 complex channels in such a way so that when the real part of the signal was taken, 32 unique real channels would result.
In FIG. 1, there is shown a complex output spectrum from a multiplexer operating in accordance with the method described herein. Within the spectrum of 0 to 2.pi., 32 adjacent complex channels, labelled Ch. 0, Ch. 1, . . . , Ch. 31, are defined. It can be seen that the channels are defined with a one-fourth channel offset, s, of the complex channel edge with respect to 0 (i.e., an even stacking scheme). This offset is critical to generating unique channels. Because of this frequency offset, the equations describing the multiplexer and demultiplexer operation exhibit the form of a Generalized Discrete Fourier Transform (GDFT), as will be seen below.
Also shown in FIG. 1 are 32 individual single FDMA channels within each complex channel, and the signal waveforms associated with each signal FDMA channel. Only half of the available bandwidth of each complex channel is used by the individual FDMA channels, which is controlled by the response of the filters used. Each of the individual FDMA channels is centered within its corresponding complex channel. This creates what are essentially "off" channels in between successive complex channels. The "off" channels are of the same bandwidth as the single FDMA channels and separate successive single FDMA channels. When the real part of the output signal is taken, the higher frequency channels "fold back" in and fill in these "off" channels, creating 32 real contiguous channels in the region from 0 to .pi., as shown in FIG. 2.
In FIG. 2, there is shown the real output spectrum using a multiplexing method in accordance with the present invention. The individual FDMA signals originally present in the 32 complex channels of FIG. 1 remain in FIG. 2. However, because of the aliasing and imaging of the FDMA channel signals symmetrically about the point .pi. on the spectrum when the real part of the output spectrum is taken, individual FDMA channel signals occur a second time in the spectrum. The choice of the quarter complex channel offset described above results in the images of the individual FDMA channel signals occurring exactly between each of the original FDMA channel signals, without overlap, in the previous "off" channels.
In FIG. 2, the individual FDMA channel signal labelled with channel number 0 in FIG. 2 occurs both immediately to the right of 0 and immediately to the left of 2.pi. in the spectrum. Moving to the right along the channel spectrum in FIG. 2, in the previously "off" channel between original FDMA channel signals 0 and 1, the image of channel signal 31 occurs. Similarly, in the previously "off" channel between original channel signals 1 and 2, the image of channel signal 30 occurs, and so on. The aliasing and imaging creates a symmetric individual FDMA signal spectrum with respect to the point .pi. on the spectrum when the real part of the output signal is taken. The real output spectrum from 0 to .pi. is repeated in the image of the real output from .pi. to 2.pi..
In FIG. 3, there is shown a diagram of analyzer (demultiplexer) channel extractions from a real input spectrum. The upper portion of FIG. 3 shows the real output spectrum and image of the real output spectrum from FIG. 2. The lower portion of FIG. 3 illustrates the complex channels which correspond to the real output spectrum and image of the real output spectrum channels of FIG. 2. The length along the normalized spectrum between successive dashed vertical lines in the lower portion indicates the defined width of each complex channel making up the demultiplexer. The complex channel width is twice the width of each individual FDMA channel from the real output and image of the real output in the upper portion of FIG. 3. The complex channel width is offset one-quarter channel, as was the case in the multiplexer of FIG. 1, so that each single FDMA channel is centered within its corresponding complex channel. The channel 0 FDMA signal thus starts at 0 on the demultiplexer spectrum.
The demultiplexer channel extractions shown in the lower portion of FIG. 3 correspond to every other FDMA channel signal from the upper portion of FIG. 3. When the actual input channels are overlaid on the demultiplexer channelization pattern, the demultiplexer "skips" the odd channels and captures only the even channels in the region from 0 to .pi.. The odd channels are not lost, however, since the demultiplexer will actually capture the images of the odd channels in the region from .pi. to 2.pi..
As an example, channels 20 and 30 in extracted form at the lower portion of FIG. 3 are shown as corresponding to the FDMA channel signals 20 and 30, respectively, of the image of the real signal at the top of FIG. 3. The demultiplexer extractions of the FDMA channel signals 0 through 16 occur from the FDMA channel real signal (the 0 to .pi. portion of the upper spectrum in FIG. 3). The demultiplier extractions of the FDMA channel signals corresponding to FDMA channel signals 17 through 31 come from the image of the real FDMA signals (the .pi. to 2.pi. portion of the upper spectrum of FIG. 3).
Modifying the multiplexer and demultiplexer equations to reflect the above described method including changes in channel stacking and translation from complex to real signals leads to the following equations: ##EQU3## where X.sub.k.sup.GDFT is the kth channel signal in the synthesizer (multiplexer) and .rho.=0, 1, . . . , M-1. Also, ##EQU4## for k=0, 1, . . . , M-1. The DFT operation in both equations contains a one-quarter channel frequency offset yielding a GDFT form.
In FIG. 4, there is shown a schematic diagram of a GDFT polyphase multiplexer 10 which operates in accordance with a preferred method of the invention. The GDFT polyphase multiplexer 10 has the following inputs:
which are input, respectively, into multipliers 12, 14, 16, and 18. The ellipses indicate that only four representatives of the m inputs are shown. In addition, the signal j.sup.m is also input into each of m multipliers (only multipliers 12, 14, 16, and 18 are shown in FIG. 4). The output of each of multipliers 12, 14, 16, and 18 form the 0, 1, k, and m-1 inputs, respectively, of the m inputs for inverse FFT 20. Inverse FFT 20 outputs 0, 1, p, and m-1 correspond to the representative inputs above and are input into multipliers 22, 24, 26, and 28, respectively. along with the respective inputs:
The output of multipliers 22, 24, 26, and 28 are input into polyphase branch filters (of which polyphase branch filters 30, 32, 34, and 36, respectively, are representatives) providing filtering functions:
respectively. The outputs of polyphase branch filters 30, 32, 34, and 36 produce the respective signals:
The polyphase branch filter outputs are then switched through operator 38, and the real portion of the polyphase branch filter outputs is taken by real function 40. The net result is multiplexed output signal x(n).
In FIG. 5, there is shown a schematic diagram of a GDFT polyphase demultiplexer 42 which operates in accordance with a preferred method of the invention. The GDFT polyphase demultiplexer 42 receives the multiplexed input signal x(n) into operator 44, which separates the multiplexed input signal x(n) into m inputs. Representative of the m inputs, the inputs x.sub.o (m), x.sub.1 (m), x.sub.p (m), and x.sub.m-1 (m) are shown in FIG. 5 as input, respectively, into polyphase branch filters 46, 48, 50, and 52. Polyphase branch filters 46, 48, 50, and 52 perform the following filtering operations, respectively:
The outputs of polyphase branch filters 46, 48, 50, and 52 are input into multipliers 54, 56, 58, and 60, respectively, which also receive the respective inputs:
Multipliers 54, 56, 58, and 60 produce outputs labelled 0, 1, p, and M-1, respectively, into FFT 62.
The FFT 62 output signals above are input into multipliers 64, 66, 68, and 70, respectively. Each of multipliers 64, 66, 68, and 70 also receives a j.sup.-m input. The output of multipliers 64, 66, 68, and 70 (and the remainder of the multipliers not shown in FIG. 5 which make up the m channels) are signals representing the demultiplexed version of the original input signal labelled 0, 1, k, and M-1, and represent the following signals:
Thus, a method for FDMA signal multiplexing and demultiplexing has been described which overcomes specific problems and accomplishes certain advantages relative to prior art methods and mechanisms. The improvements over known technology are significant. The channelization scheme modifies and combines the complex and single-sideband approaches to polyphase filter banks to allow the transmission and reception of real signals while manipulating complex baseband signals (I and Q data of a QPSK signal). The polyphase method described reduces the filter computation and storage requirements by a factor M equal to the number of FDMA channels and reduces the frequency translation computation by a factor of log.sub.2 M/M due to the use of the FFT. The multiplexing and demultiplexing method is conducive to a hardware efficient implementation for generating and separating FDMA signals. Such hardware efficiency is of paramount importance on the IRIDIUM.TM. satellite cellular communications system.
Thus, there has also been provided, in accordance with an embodiment of the invention, a method for FDMA multiplexing and demultiplexing that fully satisfies the aims and advantages set forth above. While the invention has been described in conjunction with a specific method, many alternatives, modifications, and variations will be apparent to those of ordinary skill in the art in light of the foregoing description. Accordingly, the invention is intended to embrace all such alternatives, modifications, and variations as fall within the spirit and broad scope of the appended claims.