New Directions in Channelized Receivers and Transmitters
fred harris
2-December 2011
1
New Directions in Channelized Receivers and Transmitters fred - - PowerPoint PPT Presentation
New Directions in Channelized Receivers and Transmitters fred harris 2-December 2011 1 Motivation For Using Multirate Filters 2 Processing Task: Obtain Digital Samples of Complex Envelope Residing at Frequency f C Analog Digital R
fred harris
2-December 2011
1
2
Analog Digital Multi-Channel FDM Input S ignal S ingle-Channel Base banded Output S ignal
R eceiver
Input Spectrum
f fC
S
elected Narrow B and S ignal
3
U tah
4
f f f f
L
ass Filter fs
e
1 2 3 3 4 4 2 1
Analog Signal Processing
5
f f f fs/2
1 2 3
.... ....
L
ass Filter fs
1 2 3 4
6
L
ass Filter fs/M M:1
e
5 6
f f f f fs/2 fs/2 fs/2 fs/M
fs
3 4 5 6
.... .... .... .... .... .... .... ....
L
ass Filter fs
1 2 3 4
~ ~ ~ ~
Ignoring Good Advice!
7
ADC M:1 s(t) s(n) h(n) r(n) r(nM) s(n) e -j n e -j n L O CLK L OWP AS S FIL T E R
Digital Down Converter (DDC)
8
fs/2 fs/2 fs/2 fs/M
f f f f CHANNE L OF INT E R E S T T R ANS LAT E D S P E CT R U M F IL T E R E D S P E CT R U M S P E CT R AL R E P L ICAT E S AT DOWN-S AMP L E D R AT E OU T P U T DIGIT AL F IL T E R R E S P ONS E INP U T ANAL OG F IL T E R R E S P ONS E
.... ....
9
10
Signal and Filter are at Different Frequencies Which One to Move?
f f f f CHANNE L OF INT E R E S T T R ANS LAT E D S P E CT R U M T R ANS LAT E D F IL T E R F IL T E R R E S P ONS E
Second Option First Option
11
B and P ass F ilter fs/M M:1
jn e e
5 6
f f f f fs/2 fs/2 fs/M fs/M
fs
3 4 5 6
.... .... .... .... .... .... .... ....
L
ass Filter fs
1 2 3 4
~ ~ ~ ~
No Spectral Image
12
ADC M:1 s(t) h(n) s(n) r(n) r(n) r(nM) e -j n e j n e j n L O CLK BANDP AS S FIL T E R
Up-Convert Filter, Filter Signal at IF, Down Convert Output of Filter
13
( )
j n j n k k j n j k k j n j n
Down-Convert Signal at Input to Low-Pass Filter Down-Convert Signal at Output of Band-Pass Filter Up-Convert Low Pass Filter To Become Complex Band-Pass Filter
14
Not Finished: Moving Down Converter from Input to Output Replaces 2-Multipliers (Complex Scalar) with 4-Multipliers (Complex Product)
15
ADC M:1 s(t) h(n) s(n) r(nM) r(n) r(nM) e -j Mn e j n e j n e j Mn L O CL K BANDP AS S FIL T E R
Only Down Convert the Samples we Intend to Keep! Let the Resampler Alias the Center Frequency to Baseband 0 Aliases to M0 Mod(2) 0 Aliases to 0 if M0 = k2
16
REORDERED FUNDAMENTAL OPERATION
fs/2 fs/2 fs/M fs/M
f f f f CHANNE L OF INT E R E S T T R ANS LAT E D F IL T E R F IL T E R E D S P E CT R U M AL IAS E D R E P L ICA T E S A T DOWN-S AMP L E D R AT E T R ANS LAT E D R E P L ICAT E S AT DOWN-S AMP L E D R AT E INP U T ANAL OG F IL T E R R E S P ONS E
.... .... .... ....
17
Successive Transformations Turn Sampled Data Version of Edwin Armstrong’s Heterodyne Receiver to Tuned Radio Frequency (TRF) Receiver and then to Aliased TRF Receiver.
Digital
B and-P ass M-to-1
e -j kn
H(Z e )
k Digital L
ass M-to-1
H(Z
)
e-j kn
Digital
B and-P ass
M-to-1
H(Z e )
k
k k
M θ = k 2π 2π
Any Multiple of Output Sample Rate Aliases to Baseband Armstrong Nyquist Equivalency Theorem Digital
B and-P ass
M-to-1
e
H(Z e )
k
18
19
x(n) x(n) (n) y(n) h(n) w(n) w(n) r(n) r(n) R (f)= H(f)G(f) h(n) H(f) h(n) g(n) G(f) g(n) f f f x(n) x(n) (n) v(n) g(n) w(n) w(n) r(n) r(n) R (f)= G(f)H(f) h(n) H(f) h(n) g(n) G(f) g(n) f f f
20
X X Y W W L
1
L
2
L = L L
3 1 2
X X Y W W L
1
L
2
L = L L
3 1 2
Filter R esample R esampled F ilter
21
22
2
1 1 2 2 1 2 1 2 3 3 2 2 3 3 2 2 3 3 2 2 1 2 2 1 2 3 2 3
3
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) where ( ) ( ) ( )
k k k k k k k k k k
y n x n k h k dy n y n k g k x n k k h k g k x n k h k k g k x n k h k k g k x n k f k f n h n k g k
X(Z) X(Z) X(Z) Y (Z) Y (Z) y(n) y(n) Y (Z )
M
Y (Z )
M
Y (Z )
M
Z
M
y(nM) y(nM) y(nM) x(nM+ r) x(n) x(n) x(n) H(Z) M:1 M:1 M:1 Z H(Z )
M r
Z H(Z)
r
r r
23
Coefficient Assignment of Low-Pass Polyphase Partition
Extract Delays T
Non-Zero Coefficient
For M-to-1 resample start at Index r and Increment by M For 3-to-1 resample start at index r and increment by 3
C0 C0 C0 C3 C3 C3 C4 C4 C4 C5 C5 C5 C6 C6 C6 C7 C7 C7 C8 C8 C8 C9 C9 C9 C10 C10 C10 C11 C11 C11 C1 C1 C1 C2 C2 C2
This mapping from 1-D to 2-D is used by Cooley-Tukey FFT. Polyphase Filters and CT
24
H( Z) x(n) y(n) y(nM) M:1
1
( ) ( )
N n n
H Z h n Z
1 1 ( )
( ) ( )
M N r nM r n
H Z h r nM Z
1 1
( ) ( )
M N r nM r n
H Z Z h r nM Z
Z Z Z Z
H ( Z )
M
H ( Z )
1
M
H ( Z )
2
M
H ( Z )
M-2
M
H ( Z )
M-1
M
.... .... ....
x(n) y(n) y(nM) M:1
M-Path Partition Supports M-to-1 Down Sample Also Supports Rational Ratio M-to-Q and M-to-Q/P Down Sample!
25
M Delays 1 Delay M Delays 1 Delay
Input Clock, T Output Clock, MT
Z Z
M:1 M:1 M-Units of Delay at Input Rate Same as 1-Unit of Delay at Output Rate H(Z ) H(Z )
M:1 M:1
Z Z Z Z
H ( Z ) H ( Z )
1
H ( Z )
2
H ( Z )
M-2
H ( Z )
M-1
.... .... ....
x(n) y(nM,k) M:1 M:1 M:1 M:1 M:1
H ( Z ) H ( Z )
1
H ( Z )
2
H ( Z )
M-2
H ( Z )
M-1
.... ....
x(n) y(nM)
8-tap Coefficient B ank S elect
IF S T AGE RF S T AGE IF S T AGE DDS DDS S
P
Clock Clock CARR IER P L L T IMING P LL DAC ADC DIGIT AL L OW
AS S DIGIT AL L OW
AS S DIGIT AL L OW
AS S DIGIT AL L OW
AS S DIGIT AL L OW
AS S DIGIT AL L OW
AS S DIGIT AL L OW
AS S DIGIT AL L OW
AS S S hape & U psample Matched Filter Interpolate Decimate I-Q T able Detect
cos( n) cos( n) sin( n) sin( n)
Modulator Demodulator
f Analog Analog
Analog
Oscillator P A R F Carrier Gain Control Oscillator Carrier W aveform VGA L NA
R F Carrier
IF
Channel Channel Bits Bits
Transmitter Process, Up-Sample and Up-Convert Receiver Process, Down Convert and Down-Sample
Modulator Raises Sample Rate & Applies Heterodyne at High Output Sample Rate! De-Modulator Applies Heterodynes at High Input Rate & then Reduces
DDS DIGIT AL L OW
AS S DIGIT AL L OW
AS S
M:1 M:1
DDS CIC CIC M:1 2:1 2:1 2:1 2:1 M:1 S cale Factor CIC Correction CIC Correction Half Band Half Band Half Band Half Band
z z z z-1 z-1 z-1
z-1
z z-1
M:1 Integrators Derivative Filters
DDS DIGIT AL L OW
AS S DIGIT AL L OW
AS S
M:1 M:1
… …
… … … …
0 0 1 1 M-1 M-1
DDS DDS
8-tap 8-tap Coefficient
Filter Number 1 2 3 4 5 6 7 8 9 10 Total Number Taps 3 3 3 3 7 7 7 7 11 19 70 Operations Per Filter 2‐A 2‐Shifts 2‐A 2‐Shifts 2‐A 2‐Shifts 2‐A 2‐Shifts 4‐A 2‐Mult 4‐A 2‐Mult 4‐A 2‐Mult 4‐A 2‐Mult 6‐A 3‐Mult 10‐A 5‐Mult ___ Adds Ref to Input 2 2/2 2/4 2/8 4/16 4/32 4/64 4/128 6/256 10/512 4.26 Mult Ref to Input 2/16 2/32 2/64 2/128 3/256 5/512 0.27
2 4 6 8 10 12 14 16 18
0.2 0.4 0.6 0.8 1
5 10 15 20 25 30 35 40 45
0.5 1 1.5 2
10 20 30 40 50 60 70 80 90
1 2 3 4
20 40 60 80 100 120 140 160 180 200
2 4 6 8
50 100 150 200 250 300 350 400 15 10 5 5 10 15
Impulse and Frequency Response of Last Stage Referred to Earlier Stages
Filter Number 1 2 3 4 5 6 7 8 9 10 Total Number Taps 1 1 1 1 3 3 3 3 3 4 23 Operations Per Filter 3‐A 1‐Mult 3‐A 1‐Mult 3‐A 1‐Mult 3‐A 1‐Mult 7‐A 3‐Mult 7‐A 3‐Mult 7‐A 3‐Mult 7‐A 3‐Mult 7‐A 3‐Mult 9‐A 4‐Mult ___ Adds Ref to Input 3/2 3/4 3/8 3/16 7/32 7/64 7/128 7/256 7/512 9/1024 3.25 Mult Ref to Input 1/2 1/4 1/8 1/16 3/32 3/64 3/128 3/256 3/512 4/1024 1.12
Impulse and Frequency Response of Last Stage Referred to Earlier Stages
5 10 15 20 25 30
0.2 0.4 0.6 0.8 1
10 20 30 40 50 60
0.5 1 1.5 2
20 40 60 80 100 120
1 2 3 4
50 100 150 200
2 4 6 8
100 150 200 250 300 350 400 450 15 10 5 5 10 15
2 4 6 8 10 12 14 16 18 20
Impulse Response
0.2 0.4 0.6 0.8 1
Frequency Response
0 8 0 6 0 4 0 2 0 2 0 4 0 6 0 8 1
Group Delay
Impulse, Frequency, & Group Delay Response of 2-Path Linear Phase, Recursive Half-Band Filter
h1 h1 h5 h5 h7 h7 h9 h9 1 h3 h3
1 L
as s L
as s High-P ass High-P ass Half Band F ilter: h(n) h(2n) h(2n+ 1) 2-to-1
f
0.5 0.25
L P HP
2-Point DFT
4-Point IFFT Extracts Signal Component From One-of-Four Selected Nyquist Zones
Half Band Filters Centered on Cardinal Directions Each Reduces BW 2-to-1 and Reduces Sample Rate 2-to-1
2 1 3 L
ass Filter P
Hilbert Filter High-P ass Filter Neg Freq Hilbert Filter
)
1 1 1 1 H (Z )
2
H (Z )
1
2
Z H (Z )
2
2
Z H (Z )
3
2
4-P
IFFT
4-Path Polyphase Filter Path-0 Not Used 2.5-Multiplies per Input 4-Phase Rotators fsk/4: {c0 c1 c2 c3} fs0/4: {1 1 1 1} fs1/4: {1 j -1 -j} fs2/4: {1 -1 1 -1} fs3/4: {1 -j -1 j} 2-to-1 Down-Sample
h1 h1 h5 h5 h13 h13 h9 h9 1 h17 h17
c 1 c 2 c 3
Bands Centered
(DC and fS/2) Alias To DC When Down-Sampled 2-to-1 Bands Centered
(+fS/4 and –fS/4) Alias To fs/2 When Down-Sampled 2-to-1
Showing Alias Free Pass, Transition, and Aliased Bands
Alias Free P ass Band Alias Free P ass Band T ransition Bandwidth T ransition Bandwidth Folded B andwidth Due to 2-to-1 Down S ample L
ass Band-0 P
ass B and-1 Negative F req-P ass Band-3 High-P ass Band-2
Any Narrowband Signal Must Reside in One of the 4 Alias Free Band Intervals. The Alias Free Band Intervals Overlap!
1 2
1 2
1 2
1 2
1 2
1 2
1 2
1 2
0.5
Frequency Response Log Mag (dB)
0.5
Frequency Response
Log Mag (dB)
ctra of Signal Aliased to Different Sampled Data Frequencies in Successive 2-to-1 Sample Rate Reductions.
Channelizer 4-P ath 2-to-1 Down S ample Channelizer 4-P ath 2-to-1 Down S ample Channelizer 4-P ath 2-to-1 Down S ample Channelizer 4-P ath 2-to-1 Down S ample
...... ......
Filter Channel S elect DDS
5 10 15 20 0.5 1 Impulse Response: Half Band Filter
2 2 2 2.5 [1 ] 2 4 8 1 1 1 2.5 2.5 [1 ] 2 4 8 2.5 3 7.5 N
RL I-Q I-Q I-Q
5 10 15 20 25 30 35 40 45 100
Frequency / MHz Received Signal Spectrum 17.9 17.92 17.94 17.96 17.98 18 18.02 18.04 18.06 18.08 18.1 100
F / MH Zoom-in
0.2
0.1 0.2 0.3 Frequency/MHz Stage 6
0.05 0.1 0.15
Frequency/MHz Stage 7
0.05
Frequency/MHz Stage 8
0.02 0.04
Frequency/MHz Stage 9
+0.56 = -0.44 0.28 3
= 0.14
ampled Data Frequency Locations on Successive Aliases
0.5 1
0.5 1
0.5 1
0.5 1
0.5 1
0.5 1
0.01 0.02 0.03 0.04
Frequency / MHz M a g n itu d e / d B Post Processing--Input Signal Spectrum (ch # 600)
0.01 0.02 0.03 0.04
Frequency / MHz M a g n itu d e / d B Post Processing--Down Convert (ch # 600)
M a g n itu d e / d B Post Processing--Filtering (ch # 600)
A 375-to-1 down-sample: 90 MHz to 240 kHz with a 30 kHz output BW 80 dB dynamic range. Require 6 CIC stages. The gain of each stage is 375: Gain of 6 stages becomes (375)6 or 2.8 .1015 or 52 bits growth in the CIC integrators. With 16-bit input data integrator bit width is 16+52 or 68. Six integrators in both I & Q paths would be circulating 816 bits per input sample which if converted to the 16-bit width required of the arithmetic in the half-band filters proves to be same number
Number of operations for the I-Q half band filter chain is on the
represents a workload 1/6 of the CIC chain. The efficient cascade CIC filter chain can be replaced with an even more efficient cascade four-path half band filter chain.
0.1 0.2 0.3 0.4 0.5 8 2
Group Delay: Two Path, 4-Coefficient, Linear Phase 2-Path Filter
0.5
0.1 0.2 0.3 0.4 0.5 9 8 1
Group Delay: Detail
0.1 0.2 0.3 0.4 0.5
Frequency response
Channelizer 4-P ath 2-to-1 Down S ample Channelizer 4-P ath 2-to-1 Down S ample Channelizer 4-P ath 2-to-1 Down S ample Channelizer 4-P ath 2-to-1 Down S ample
...... ......
Filter Channel S elect DDS
2 2 2 2.0 [1 ] 2 4 8 1 1 1 2.0 2.0 [1 ] 2 4 8 2.0 3 6.0 N
RL I-Q I-Q I-Q
5 10 15 20 25 30
Impulse Response, Two-Path, 4-Coefficient, Linear Phase IIR
2 1 1
( ) ( )
M N j k r nM M r n
G Z e Z h r nM Z
1 1
( ) ( ) ( )( ) ( )
k k k
N N j n j j n n n n
G Z h n e Z h n e Z H e Z
1 1 ( ) ( )
k
M N j r nM r nM r n
1 1
( ) ( )
k k
M N j j nM r nM r n
G Z e Z h r nM e Z
k k
θ = k 2π 2π θ = k M
Modulation Theorem of Z-Transform
Z Z Z Z
H ( Z )
M
H ( Z )
1
M
H ( Z )
2
M
H ( Z )
M-2
M
H ( Z )
M-1
M
.... .... ....
x(n) y(n) y(nM) M:1 e e e e e j k0 j k1 j k2 j k(M-2) j k(M-1) 2 2 2 2 2
M M M M M
Z Z Z Z
H ( Z ) H ( Z )
1
H ( Z )
2
H ( Z )
M-2
H ( Z )
M-1
.... .... ....
x(n) y(nM,k) M:1 M:1 M:1 M:1 M:1 e e e e e j k0 j k1 j k2 j k(M-2) j k(M-1) 2 2 2 2 2
M M M M M
We Reduce Sample Rate M-to-1 Prior to Reducing Bandwidth (Nyquist is Raising His Eyebrows!) We Intentionally Alias the Spectrum. (Were you Paying Attention in school when they discussed the importance of anti-aliasing filters?) M-fold Aliasing! M-Unknowns! M-Paths supply M-Equations We can the separate Aliases!
Z Z Z Z
H ( Z ) H ( Z )
1
H ( Z )
2
H ( Z )
M-2
H ( Z )
M-1
.... .... ....
x(n) y(nM,k) M:1 M:1 M:1 M:1 M:1 e e e e e j k0 j k1 j k2 j k(M-2) j k(M-1) 2 2 2 2 2
M M M M M
W t Ph S i f f l ibl
Polyphase Partition with Commutator Replacing the “r” Delays in the “r-th” Path
H ( Z ) H ( Z )
1
H ( Z )
2
H ( Z )
M-2
H ( Z )
M-1
.... ....
x(n) y(nM,k) e e e e e j k0 j k1 j k2 j k(M-2) j k(M-1) 2 2 2 2 2
M M M M M
Note: We don’t assign Phase Spinners to Select Desired Center Frequency Till after Down Sampling And Path Processing This Means that The Processing for every Channel is the same till the Phase Spinner No longer LTI, Filter now has M-Different Impulse Responses! Now LTV or PTV Filter.
Digital
B and-P ass
M-to-1
H(Z e )
k
Digital L
ass M-to-1
H(Z
)
e-j kn
Rather than selecting center frequency at input and reduce sample rate at output, we reverse the order, reduce sample rate at input and select center frequency at output. We perform arithmetic
M-P ath Digital P
M-to-1
H(Z )
r
e
M rk
n n-1 n-2 n-3 n-4 n-5 n-6 n-7 n-8 n-9 n n-1 n-2 n-3 n-4 n-5 n-6 n-7 n-8 n n-1 n-2 n-3 n-4 n-5 n-6 n-7 n-8 n+ 6 n+ 5 n+ 4 n+ 3 n+ 2 n+ 1 n-10 n-11 n-12 n-13 n-14
— — — — — —
n n-1 n-2 n-3 n-4 n-5 n n-1 n-2 n-3 n-4 n-5 n n-1 n-2 n-3 n-4 n-5 n-12 n-13 n-14 n-15 n-16 n-17 n-12 n-13 n-14 n-15 n-16 n-17 n-12 n-13 n-14 n-15 n-16 n-17 n+ 6 n+ 5 n+ 4 n+ 3 n+ 2 n+ 1 n-6 n-7 n-8 n-9 n-10 n-11 n-6 n-7 n-8 n-9 n-10 n-11 n-6 n-7 n-8 n-9 n-10 n-11
— — — — — — — — — — — —
ADC s(t) h(r+ nM) h(M-1+ nM) h(0+ nM) h(1+ nM) s(n) r(nM) r(nM,k) CL K B ANDP AS S FIL T E R P OL Y P HAS E P AR T IT ION L OW P AS S FIL T E R P OL Y P HAS E P AR T IT ION P HAS E R OT AT OR S AL IAS E D HE T E R ODY NE e j r k
M
... ... ..... .....
.. this is very stuff....
Impulse Response of 6-Path Polyphase Partition Prior to 6-to-1 Resampling
Frequency Response of 6-Path Polyphase Partition Prior to 6-to-1 Resampling
Phase Response of 6-Path Polyphase Partition Prior to 6-to-1 Resampling
Overlay Phase Response of 6-Path Polyphase Partition Prior to 6-to-1 Resampling
Nyquist Zone +1 Nyquist Zone +2 Nyquist Zone
Nyquist Zone
Nyquist Zone Phase Aligned 2/6 Phase Shifts
Phase Shifts
Phase Shifts 2/3 Phase Shifts
3-D Paddle-Wheel Phase Profiles, 6-Path Partition Prior to 6-to-1 Resampling
Overlay 3-D Paddle-Wheel Phase Profiles, 6-Path Partition Prior to 6-to-1 Resampling
Overlay 3-D Paddle-Wheel Phase Profiles, Showing Phase Shifts in +1 Nyquist Zone
Overlay 3-D Paddle-Wheel Phase Profiles, Phase Shifted to Align Phases in +1 Nyquist Zone
PolyChanDemo
H(Z ) H (Z )
y(nM,k) y(nM,k) y (nM) y(n,k)
M-to-1 M-to-1
1 1 2 2 2 2 M M
x(n) x(n) x(n)
M M j nk j nk 2 2
r r S tandard DDC P
2-P
Filters 1-P
Filter
ide Input Down Conversion to Output of Filter Where it anishes Due to Down Sampling. Rotators in Filter Factor Out and are Applied to Path Outputs Rather than to Coefficients. Advantage: Real sequence is made complex at output of Filter Rather than at Input to Filter
… …
… … … …
0 0 1 1 M-1 M-1
DDS
H (Z ) H (Z ) H (Z ) H (Z )
1 2 M-1
. . . . . .
x(n) y(nM,k)
e
M j 1k2
e
M j 0k2
e
M j 2k2
e
M j (M-1)k2
Nyquist Rate for Filter is 200 kHz + 200 kHz = 400 kHz or fs/50
200 kHz 80 dB 0.1 dB 200 kHz 20 MHz
f
/Octave 400 T ap FIR Filter 20 MHz Input S ample R ate 20 MHz Output S ample R ate
Polyphase Partition of Low‐Pass Filter
…
0 1 2 49 48
50-to-1 400 T aps 8 T aps 20 kHz 20 MHz 400 kHz 400 kHz P
L
ass Filter
Cascade Polyphase Filter Down‐Sampling and Up‐Sampling
… …
… … … …
0 0 1 1 2 2 49 49 48 48
50-to-1 1-to-50 400 T aps 400 T aps 20 MHz 20 MHz 20 MHz 20 MHz 400 kHz 400 kHz P
L
ass Filter P
L
ass Filter 8-taps 8-taps
8-tap 8-tap 20 MHz 20 MHz 400 kHz Coeffic ient B ank Coeffic ient B ank S elec t S elec t
400 T ap FIR Filter 20 MHz Input S ample R ate 20 MHz Output S ample R ate
16‐Ops/Input 60‐Ops/Input
400-T ap L
Filter 20MHz 20MHz 20MHz 400 kHz 20MHz 8-T ap F ilter 8-T ap F ilter
Coefficient Bank Coefficient Bank S tate Machine S elect S elect White Box White Box
H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z )
1 1 2 2 M-1 M-1
...... ......
x(n) y(nM) y(n)
f f f
H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z )
1 1 2 2 M-1 M-1
...... ......
x(n) y(nM,k) y(n,k)
e
M j 1k2
e
M j 1k2
e
M j 0k2
e
M j 0k2
e
M j 2k2
e
M j 2k2
e
M j (M-1)k2
e
M j (M-1)k2
f f f
H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z ) H (Z )
1 1 2 2 M-1 M-1
...... ......
y(nM,k )
1
y(nM,k )
2
y(n,k )+ y(n+ k )
1 2
e
M j rk12
e
M j rk22
e
M j rk12 e M j rk22
f
Workload for Multiple M-Path Filters
hen k > Log2(M)/2 Build all channels and discard the channels you don’t need! M=16, Log2(16)/2 = 2: thus if you want 2 or more, Build them all! M=128, Log2(128)/2 = 3.5: thus if you want 4 or more, Build them all! M=1024, Log2(1024)/2 = 5: thus if you want 5 or more, Build them all!
M-Channel Channelizer: Resampled M-Path Narrowband Filter Channels Alias to Baseband: Phase Aligned Sums Separate Aliases: rk Performed at Low Output Rate Rather Than at High Input Rate. One Input Filter Services M-Output Channels
fs h (n) h (n)
2
h (n)= h(r+ nM)
r
P
P artition h (n)
M-2
h (n)
1
h (n)
3
h (n)
M-1
FDM T DM
M-PNT IFF T
.. this is very stuff....
H(Z ) H(Z ) H (Z )
y(nM,k ) y(nM,k ) y(nM,k ) y(nM k ) y (nM) y(n,k ) y(n,k )
e e
M-to-1 M-to-1 M-to-1
1 1 1 2 2 2 2 2 2 2 M M
x(n) x(n) x(n)
e e
M M M M j nk j nk j nk j nk 2 2 2 2
r r S tandard DDC P
1 1 1 1 1 2 2 2 2 2
2-P
Filters 2-P
Filters 1-P
Filter
y(m) x(m) x(n)
6 6 3 3 3 9 9 1 1 2 2 1 7 7 4 4 4 10 10 20 20 2 8 8 5 5 5 25 25 15 15 m n m
fs fs fs f f f 5 5
.... .... .... .... .... ....
F ilter R ejected S pectrum R ejected S pectrum
m) m) n)
6 3 3 9 1 2 1 7 4 4 10 20 2 8 5 5 25 15 m n m
fs fs fs f f f 5 5
.... .... ....
F ilter R ejected S pectrum R ejected S pectrum
1
( ) ( )
N n n
H Z h n Z
C0 C3 C4 C5 C6 C7 C8 C9 C12 C10 C13 C11 C14 C1 C2
1-to-5
1 1 ( )
( ) ( )
N M M r nM r n
H Z h r nM Z
C0 C3 C4 C5 C6 C7 C8 C12 C9 C13 C10 C14 C11 C1 C2
1-to-5
Polyphase Partition of Resampling Filter
1 1
( ) ( )
N M M r nM r n
H Z Z h r nM Z
C0 C3 C4 C5 C6 C7 C8 C12 C9 C13 C10 C14 C11 C1 C2
1-to-5
C0 C3 C4 C5 C6 C7 C8 C12 C9 C13 C10 C14 C11 C1 C2
1-to-5
Factor Delays and Rearrange
M Delays 1 Delay M Delays 1 Delay
Input Clock, T Output Clock, MT Z H(Z ) Z H(Z )
M:1 M:1 M:1 M:1 x(n) y(m) y(m) x(n)
C0 C3 C4 C5 C6 C7 C8 C12 C9 C13 C10 C14 C11 C1 C2
1-to-5 1-to-5 1-to-5 1-to-5 1-to-5
C0 C3 C4 C5 C6 C7 C8 C12 C9 C13 C10 C14 C11 C1 C2
Replace Up Samplers, Delays, and Summer with M-Port Output Commutator
2 1
( ) ( )
N j k n n M n
G Z h n e Z
1 2 1 ( ) ( ) 1 2 2 1 1 2 1
( ) ( ) ( ) ( ) ( ) ( )
N M M j r nM k r nM M r n N M M j rk j nM r nM M M r n N M M j rk r nM M r n
G Z h r nM e Z G Z Z e h r nM e Z G Z Z e h r nM Z
Spin The Delays, Don’t Touch the M-Path Partitioned Weights
fs fB
W
fs M f
F0 F1 F2 FM-2 FM-1
M-Point IFFT Supplies Phase Spinners to Form Up Converters to all Multiples of Input Sample Rate
All Output Channels Centered on Multiples of Input Sample Rate Example: Multiples of 6-MHz
h (5n) h (5n) h (5n)
2
h (5n)
2
h (5n)
4
h (5n)
4
h ( 5n)
1
h ( 5n)
1
h ( 5n)
3
h ( 5n)
3
x(n) y(m,k) y(m,k) e e e
j j j n 2 2 5 5 kr kr
Input Offset Frequency
Input Offset Observed at Output
y(m) x(m) x(n) x(n)
3 1 4 2 5 m n n m
fs fs fs fs f f f f 5 5 F ilter R ejected S pectrum R ejected S pectrum
Input Spectrum Shifted Spectrum Up-Sampled Spectrum Distorted Spectrum
y(m) x(m) x(n) x(n)
3 1 4 2 5 m n n m
fs fs fs fs f f f f 5 5 2 F ilters R ejected S pectrum R ejected S pectrum
Input Spectrum Shifted Spectrum Up-Sampled Spectrum T wo Filter Spectra
fs h (n) h (n)
2
h (n)= h(r+ nM)
r
P
P artition h (n)
M-2
h (n)
1
h (n)
3
h (n)
M-1
FDM T DM
M-P NT FF T
Critically Sampled fS=fC Nyquist Sampling fS=2fC QRT Nyquist Filter fS=2fC
QRT Nyquist Filter fS=4fC Nyquist Sampling fS=4fC
Various Filter-Channelizer Configurations
f f f f f f f f f f fs fs fs fs fs f = fs
f = fs
f = fs
f = fs
f = fs
f = fs
B W
f = fs
B W
f = fs
B W
f = 2fs
B W
f = fs
B W
f = fs
B W
f = fs
B W
f = fs
B W
ch(k-1) ch(k-1) ch(k-1) ch(k-1) ch(k-1) ch(k+ 1) ch(k+ 1) ch(k+ 1) ch(k+ 1) ch(k+ 1) ch(k) ch(k) ch(k) ch(k) ch(k) 0.8 2.0 2.0 4.0 1.0 1.0 1.0 1.0 2.0 0.1 dB 0.1 dB 3.0 or 6.0 dB 6.0 dB 0.1 dB
f f
B W B W S ample R ate= 2 Channel S pacing S ample R ate = Channel S pacing Channel S pac ing Channel S pac ing Aliased T ransition B andwidth
(M 1)
H ( Z )
2
H ( Z )
1 2
H ( Z )
2 2
H ( Z )
2
H ( Z )
2
H ( Z )
2 2
.... ....
y(n ) y(n)
M 2
Z Z Z Z Z
.... .... .... .... ....
....
M M M M M M 2 2 2 2 2 2M 2 :1
+
( 1)+
( - 1)(M 1)
H ( Z )
2
H ( Z )
1
2
H ( Z )
2
2
H ( Z )
2
H ( Z )
2
H ( Z )
2
.... ....
y(n )
M 2
Z Z Z Z Z
.... .... .... .... ....
....
M
M M M M M M
2
2 2 2 2 2 2
:1
M 2 :1 M 2 :1 M 2 :1 M 2 :1 M 2 :1 M 2 :1
1)
+
( 1)
+
( -1)
( )
Path Filters: Polynomials in ZM Converted to Polynomials in Z2
H ( Z )
2
H ( Z )
1
2
H ( Z )
2
2
H ( Z )
2
H ( Z )
2
H ( Z )
2
H ( Z )
2
.... ....
y(n )
M 2
Z Z Z Z Z Z
.... .... .... ....
....
M
M M M M M
2
2 2 2 2
:1
M 2 :1 M 2 :1 M 2 :1 M 2 :1 M 2 :1 M 2 :1
( 1)
+
( -1) ( )
Use Noble Identity to Pull M/2‐to‐1 Resampler Through Delays in Lower Half of Paths
Z Z Z
H ( Z ) H ( Z )
1
H ( Z )
2
H ( Z )
( -1)
H ( Z )
( )
H ( Z )
M-2
H ( Z )
M 1
.... .... .... .... ....
x(n) y(n )
M 2
M 2 M 2 2 2 2 2 2 2 2
Replace Delays and M/2‐to‐1 Resamplers with Dual Input M/2 Path Commutator
H (Z )
k 2
Z H (Z )
2 (k+ M/2)
h(k) h(k+M/2) h(k+M) h(k+3 ) h(k+2M) h(k+5 ) h(k+3M) h(k+7 ) h(k+4M) h(k+9 ) M 2 M 2 M 2 M 2
fs f fs fs fs fs fs fs fs fs M M M M M M M M M 2 2
esample M/2-to-1 R esample
..... ..... ..... .....
f f
0 5T T 1.5T T T 0.5T 0.5T mT (m+1)T
fs 2 fs/M M Path Poly Phase Filter in Z M-PNT IFFT
.....
flg=0
flg=0
flg=1 f l g = 1
Circular Buffer
2
1 1 2 2 3 M/2-1 M-1 M/2 M/2+ 1 M-2 M-1
M M M
ath P
M-P
F DM
T DM
M-P ath Input Data B uffer
uffer
S tate E ngine
1 1 2 2 3 M/2-1 M-1 M/2 M/2+ 1 M-2 M-1
M M M
ath P
M-P
FDM T DM
M-P ath Input Data B uffer
uffer
S tate E ngine
M P
IF FT M P
IF FT M P ath F ilter M P ath F ilter
nM/2 P
S hift M-P
Circular Buffer nM/2 P
S hift M-P
Circular Buffer
x x x x
Frequency Domain Filtering With Cascade M/2‐to‐1 Analysis and 1‐to‐M/2 Synthesis Channelizers
Impulse response and Frequency Response 25‐Enabled Ports: 2.4 MHz Bandwidth
.... .... .... .... .... .... .... .... .... .... .... ....
S elector
Input Channel Configuration M M M M N N
Anal. Anal. N-1 B w0 P P
N-1
B WN-1 2-to-M, M-Channel S ynthesis Channelizer 2fs Mfs
1 2 M/2-1 M/2 M/2+ 1 M-1
.... …. M-P
ath P
M-P
M-P ath Input Data B uffer
Circular Output B
uffer
S tate E ngine
Analysis Channleizers
1 2 M/2-1 M/2 M/2+ 1 M-1
....
.... .... .... .... .... .... .... .... .... .... .... ....
…. M-P
ath P
M-P
S elector M-P ath Input Data B uffer
uffer
S tate E ngine M M M M N N
S ynth. S ynth. N-1 B w0 P P
N-1
B WN-1 M/2-to-1, M-Channel Channelizer fs
10
40 FLTR 0
10
40 FLTR 1
10
40 FLTR 2
10
40 FLTR 3
10
40 FLTR 4
10
40 FLTR 5
10
40 FLTR 6
10
40 FLTR 7
10
40 FLTR 8
10
40 FLTR 9
10
40 FLTR 10
10
40 FLTR 11
10
40 FLTR 12
10
40 FLTR 13
10
40 FLTR 14
10
40 FLTR 15
10
40 FLTR 16
10
40 FLTR 17
10
40 FLTR 18
10
40 FLTR 19
10
40 FLTR 20
10
40 FLTR 21
10
40 FLTR 22
10
40 FLTR 23
10
40 FLTR 24
10
40 FLTR 25
10
40 FLTR 26
10
40 FLTR 27
10
40 FLTR 28
10
40 FLTR 29
10
40 FLTR 30
10
40 FLTR 31
10
40 FLTR 32
10
40 FLTR 33
10
40 FLTR 34
10
40 FLTR 35
10
40 FLTR 36
10
40 FLTR 37
10
40 FLTR 38
10
40 FLTR 39
40 FLTR 40
40 FLTR 41
40 FLTR 42
40 FLTR 43
40 FLTR 44
40 FLTR 45
40 FLTR 46
40 FLTR 47
20
40 Sig 3 Magnitude (dB)
20
40 Sig 6
20
40 Sig 8
20
40 Sig 4 Frequency (MHz) Magnitude (dB)
20
40 Sig 5 Frequency (MHz)
20
40 Sig 7 Frequency (MHz)
10 20 30
40 Sig 1 F (MH ) Magnitude (dB)
10 20 30
40 Sig 2 F (MH )
.. ..
.. .. .. .. .. ..
M P
IF FT M P
IF FT M P ath F ilter M P ath F ilter
nM/2 P
S hift M-P
Circular Buffer nM/2 P
S hift M-P
Circular Buffer
HERE POINTED HAIR BOSS. THIS REPORT EXPLAINS HOW A SMALL FREQUENCY OFFSET AT THE INPUT SAMPLE RATE IS CONVERTED TO THE SAME FREQUENCY OFFSET FROM THE CHANNEL CENTER FREQUENCY AT THE HIGH OUTPUT SAMPLE RATE.