The present disclosure relates generally to the field of signal processing. More particularly, the present disclosure relates to high-performance filter bank channelizers.
The M-path polyphase analysis filter bank channelizer is quite a remarkable digital signal processing technology. In its simplest realization, the maximally decimated filter bank, the bank outputs M baseband time series from translated spectral spans with bandwidth and sample rate fs/M from M spectral bands centered at integer multiples of fs/M. Modifications to the channelizer are many and include offsets of channelizer center frequencies, non-maximal decimation from M-to-1 to M/2-to-1 or to 3M/4-to-1 along with various post channelization signal conditioning options.
While filter bank channelizers are of significant use and importance in the digital signal processing field, it would be beneficial to improve the performance of such channelizers by reducing the amount of computational processing that must be carried out by digital signal processors and other devices on which such channelizers are implemented. Doing so would significantly improve the performance and speed of such channelizers. Accordingly, what would be desirable are high-performance filter bank channelizers which address the foregoing, and other, needs.
The present disclosure relates to high-performance filter bank channelizers. In one embodiment, a high-performance channelizer is provided which includes a digital direct synthesis (DDS) module generating a heterodyne signal; a mixer in communication with DDS module and mixing the heterodyne signal with an input signal; and an M-path channelizer in communication with the mixer, the M-path channelizer processing an output signal of the mixer to generate a plurality of output channels, wherein the heterodyne signal shifts an input spectrum of the input signal by an offset between input and output centers. The heterodyne signal operates at a high input sample rate.
In another embodiment, a high-performance channelizer is provided which includes an input commutator receiving and commutating an input signal; an M-path polyphaser filter in communication with the commutator; and an M-path inverse discrete Fourier transform module processing outputs of the polyphaser filter, wherein the M-path polyphaser filter introduces a plurality of phase rotations in a time domain resulting in a workload reduction for a processor on which the channelizer is implemented. The plurality of phase rotations can be inserted at a rate of 1/30th of an input rate of the channelizer.
In another embodiment, a resampling channelizer is provided and includes a frequency division multiplex (FDM) commutator receiving and commutating an FDM input signal; an M/2-path input data buffer in communication with the FDM commutator; an M-path polyphaser filter in communication with the input data buffer; a circular output buffer in communication with the M-path polyphaser filter; an M-point inverse fast Fourier (IFFT) module in communication with the circular output buffer; and a time division multiplex (TDM) commutator in communication with the M-point IFFT module and generating a TDM output signal, wherein the M-path polyphaser filter is operated at a sample rate above fs/M.
In still another embodiment, a half-band filter is provided, and includes an upper filter path including even indices of a low-pass filter; a lower filter path including even symmetric filter coeeficients; a switch in communication with the upper and lower filter paths and switching an input signal between the upper and lower filter paths; and a mixer in communication with the upper and lower filter paths and mixing outputs of the upper and lower filter paths.
In yet another embodiment, a cascaded half-band filter is provided, and includes an input commutator receiving and commutating an input signal; a first M-path filter in communication with the input commutator; a first M-point circular buffer in communication with the first M-path filter; a first M-point inverse fast Fourier transform (IFFT) module in communication with the first M-point circular buffer; a second M-point IFFT module in communication with the first IFFT module; a second M-point circular buffer in communication with the second M-point IFFT module; a second M-path filter in communication with the second M-point circular buffer; and an output commutator in communication with the second M-path filter and generating an output signal, wherein the input commutator, the first M-path filter, the first M-point circular buffer, and the first M-point IFFT module form an analysis channelizer, and the second M-point IFFT module, the second M-point circular buffer, the second M-path filter, and the output commutator form a synthesis channelizer, the analysis channelizer cascaded with the synthesis channelizer.
The foregoing features of the invention will be apparent from the following Detailed Description of the Invention, taken in connection with the accompanying drawings, in which:
The present disclosure relates to high-performance filter bank channelizers, as described in detail below in connection with
In its most common incarnation, a polyphase down sampling channelizer simultaneously down converts and down samples M equally spaced, fixed bandwidth signals.
The time delay response of each path filter aligns the time origin of their sampled data sequences formed at their outputs to a single common output time origin. This task is accomplished by the all-pass characteristics of the M-path partitioned filter that apply the required differential time delay to the individual input time series. Finally, the IFFT block performs the equivalent of a beam-forming operation; the coherent summation of the time aligned signals at each output port with selected phase profile. Note that the channel spacing, the channel bandwidth, and the sample rate are all fs/M. This form of the channelizer is called a maximally decimated filter bank.
As a multi-channel channelizer in which we extract and separate adjacent channels, the signal bandwidth must be less than the channel spacing. Under this condition, there is spectral gap between the input channel bands. The gap is required for the channel filters to have a non-zero transition bandwidth between the channel bands. Discussed below are the signal bandwidths and channel spacing, as well as the required filter characteristics. Also discussed below is the option of increasing the transition bandwidth of the channelizer and follow the channelizer with filters that form the desired narrow transition bandwidth. Because these filters operate at the reduced output sample rate, their reduced length and clock speed offer significant implementation advantages.
The expression to determine the number of taps in a finite impulse response (FIR) filter is shown in Equation (2) where fs is sample rate, Δf is transition bandwidth and K(A) is a parameter proportional to A, the out of band attenuation level. The estimate from Equation (2) sets the filter length to 3273 taps. When designed by the FIRPM algorithm the estimate proved to be very good. We adjusted the filter length to be 1 less than the closest multiple of 30, a filter of length 3269 taps which met the design specifications. When partitioned into the 30 arms of the 30 path polyphase filter we find each arm contains 109 taps:
We can implement the 30 path filters of the channelizer with each path operating at 1/30th of the input rate 720 MHz/30 or 24 MHz, a comfortable speed for an FPGA implementation. We simulated the design in MATLAB, and the spectral responses of the prototype filter are shown in
The rotator sequence is periodic in twice the length of the output vectors formed by the polyphase filter. Noting the sign change at the midpoint of the rotator vector we apply the rotators to the filter output in the same way the lower half of the butterfly of a radix-2 FFT forms its sum. We form the weighted sum of the even indexed data vectors and the weighted sum of the odd indexed data vectors and apply the complex rotator weights to their difference. Alternate vector outputs have to be sign changed to keep the channel spectral bin center at direct current (DC) rather than at the half sample rate.
The modified form of the channelizer incorporating the half band width frequency offset discussed above is shown in
We now consider a modification to the channelizer that holds promise of reduced workload but still meets the design requirements. We recognize that the workload in the channelizer is dominated by the large number of coefficients in the polyphase filter partition. As commented on earlier, in (2), this number is large because the ratio of sample rate to transition bandwidth is large. We can reduce the channelizer computational workload if we increase the transition bandwidth. This would reduce the filter length but would result in a filter that doesn't meet the design requirements. Our response to this problem is we use a second filter applied to the channelizer output to form the reduced transition bandwidth at reduced cost because of its reduced sample rate.
If we increase the transition bandwidth of the channelizer filter we will also have to increase the output sample rate of the channelizer. We present a modified version of the Nyquist theorem. The Nyquist criterion tells us that the sample rate should exceed the two-sided bandwidth, leaving the question: “By how much?” The channelizers disclosed herein answer that question. As shown in Equation (4) below, we should exceed the signal's two-sided BW by the transition BW of the anti-alias filters.
Nyquist: fs>2 Sided BW
harris: fs=2 sided BW+Filter Transition BW
The excess bandwidth typically increases the filter sample rate (fs) by 10% to 20%. In the modern era, we raise the sample rate to accommodate a significant increase in transition BW and then use a follow-up DSP filters to reduce the bandwidth and sample rate to the desired lower values.
When we run the M path polyphase filter bank at rates above fs/M, the architecture changes and the channelizer is known as non-maximally decimated filter bank. We have a few options for the amount to raise the output sample rate. One common and easy to implement options is to double the sample rate from fs/M to 2 fs/M. For our particular example we raise the output sample rate from fs/30 or 720/30 or 24 MHz to 720/15 or 48 MHz by delivering 15 samples to the channelizer and form 30 output samples for every 15 input samples. We could have selected some other ratio which increased the sample rate by a smaller amount such as 720/20 or 36 MHz by delivering 20 samples to the channelizer and form 30 output samples for every 20 input samples. In the first case we would increase the sample rate by 100% and in the second we would increase the sample rate by 50%. The option space for the first case is quite wide. Whichever choice we would finally select, we would be sure to include the internal frequency shift option in the resampled channelizer. The common configuration of the M/2-to-1 resampling channelizer is shown in
We designed the 30 path channelizer to meet the specifications indicated in
We now have the task of designing the cascade filter that will reduce the transition bandwidth of channelized time series to the desired 0.5 MHz. It would be nice if that filter be configured to reduce the sample rate from 48 MHz to 24 MHz. The first filter option that comes to mind is a true half-band finite impulse response (FIR) filter. Access to this option is why we selected 48 MHz for the channelizer output sample rate. The qualifier, true for the half band filter, is that we want the design to have zero value on alternate output samples. We can achieve that goal with a windowed sinc series or with the half-band technique that uses the FIRPM algorithm to design the filter's odd index non-zero weights and inserting the even index zeros and center tap. The former design is characterized by non-uniform pass band and stop band ripple levels while the later has the equal ripple response of the standard FIRPM design.
We selected the half-band technique and designed the half band filter to operate at the 48 MHz sample rate with the desired 0.5 MHz transition BW and the 50 dB stopband levels. The filter length require to meet these requirements is 233 taps and the spectral characteristics of the filter is shown in
A second option for the half-band filter is a linear phase all-pass infinite impulse response (IIR) filter. We designed and simulated the IIR version of the half band filter which has a 2-to-1 down sampling option implementation very similar to that of
A third option for the half band filter is the cascade of a pair of analysis and synthesis channelizers. The analysis synthesizer partition the input spectrum into a set of reduced sample rate baseband channels. The prototype filter in the analysis channelizer is Nyquist filter designed so adjacent channels cross at their-6 dB levels. The synthesis channelizer performs the perfect reconstruction of these baseband channels as they are aliased up to their original center frequencies by the M/2 up-sampling process. The bandwidth filtering option is performed by the binary mask between the analysis and synthesis banks. The stop band corresponds to the channels not participating in the assembly of the super channel formed by the channels passed to the synthesis process from the output of the analysis process. This architecture is shown in
We designed the 40 path filter with 6-taps per path and synthesized the half band filter with 20 selected channels in an offset bin channelizer. The spectral characteristics are shown in FIG. 14. The computational workload for the synthesized filter is 12 multiplies per input sample, 12 multiplies per output sample for the path filters and 10 multiplies per input sample for the two 40 point IFFTs. When the output 40 path filter is replaced with a 20-path filter to reduce the output sample rate the workload for the output channelizes is cut in half for a combined workload of 25 multiplies per input sample.
It is noted that the various channelizers and filters disclosed herein could be implemented using any suitable processor such as an application-specific integrated circuit (ASIC), a digital signal processor (DSP), a field-programmable gate array (ASIC), a microprocessor, or as software executed by a general-purpose processor. It is additionally noted that the channelizers and filters could be implemented in a radiofrequency transceiver, which could include, but is not limited to, a cellular transceiver (e.g., base station or mobile device supporting one or more communications protocols such as 3GPP, 4G, 5G, etc.), a satellite transceiver (e.g., an earth station or a satellite in space), a wireless networking transceiver (e.g., a WiFi base station or WiFi-enabled device), a short-range (e.g., Bluetooth) transceiver, or any other radiofrequency transceiver.
Advantageously, the channelizers and filters disclosed herein meet a severe set of specifications, such as operation at high sample rate fs with specifications that lead to very long filter lengths. Fortunately, the M-path polyphase channelizer performs M-to-1 down sampling to each of the M-paths. This means that each path operates at the reduced sample rate fs/M. Two filters can be implemented, one operating at the high input sample rate and one operating at the lower output sample rate. In this process, the first filter, with a wider transition bandwidth, reduces the bandwidth and sample rate.
Having thus described the system and method in detail, it is to be understood that the foregoing description is not intended to limit the spirit or scope thereof. It will be understood that the embodiments of the present disclosure described herein are merely exemplary and that a person skilled in the art can make any variations and modification without departing from the spirit and scope of the disclosure. All such variations and modifications, including those discussed above, are intended to be included within the scope of the disclosure. What is desired to be protected by Letters Patent is set forth in the following claims.
This application is a continuation of U.S. patent application Ser. No. 17/559,834 filed Dec. 22, 2021, which is related and claims priority to provisional patent application Ser. No. 63/129,980, filed on Dec. 23, 2020, and entitled High-Performance Filter Bank Channelizers, the entire contents of which is expressly incorporated by reference herein.
Number | Date | Country | |
---|---|---|---|
63129980 | Dec 2020 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 17559834 | Dec 2021 | US |
Child | 18611548 | US |