計測制御工学 第 6 回講義
複素アナログフィルタの基礎
小林春夫
群馬大学大学院理工学府 電子情報部門
[email protected]
下記から講義使用
出席・講義感想もここから入力してください。
2021
年5
月24
日(
月)
/19 1
欧州のアナログ集積回路 研究教育事情
群馬大学 小林春夫
2020/5/11
Disclaimer:
15年前に訪問した時の情報もあり、
一部は最新情報ではない
/19
欧州地図 2
/19
意外な欧州 3
● 日本は面積も人口も少ない?
欧州で面積&人口が日本より大の国なし
● 欧州は人口少ない、大都市少ない アジアは多い
「アジアは米、欧州は小麦が主食のため」の説 米は何年も備蓄可、小麦は不可
● 「イギリス(ブリテン島)と欧州大陸は別」
の意識が双方にあり (?)
/19
米国と欧州は異なる 4
/19
欧州 2000 年以上の西洋科学の伝統の底力 5
/19
オランダ 6
/19
ベルギー(大学) 7
/19
ベルギー(研究所) 8
/19
イタリア 9
/19
イギリス、アイルランド 10
イギリス: 海外では UK を使う
England, Scotland, Wales, Northern Ireland
の 連合王国United Kingdom
サッカー ワールドカップにも独立
4
チーム出場/19
ドイツ 11
/19
スイス 12
/19
北欧 13
/19
ロシア 14
半導体メーカ技術者より
/19
ポルトガル、イスラエル 15
/19
スペイン バルセロナ市 雑景 16
サグラダファミリア(Sagrada Familia) 聖家族 贖罪教会
/19
フランス 17
/19
フランス パリ 雑景 18
パリ第
6
大学 セーヌ川 エッフェル塔 凱旋門/19
日本企業への売り込み 19
/19
関係ファイル 20
https://kobaweb.ei.st.gunma-u.ac.jp/news/pdf/2019/oldrepo_NapoliUniv.pdf https://kobaweb.ei.st.gunma-u.ac.jp/news/pdf/2019/oldrepo_Europe.pdf
https://kobaweb.ei.st.gunma-u.ac.jp/warehouse/20160722am9IMSTW.pdf
https://kobaweb.ei.st.gunma-u.ac.jp/warehouse/IMSTW20150703.pdf
Complex Signal Processing
in Analog / Mixed-Signal Circuits
Haruo Kobayashi
Minh Tri Tran, Koji Asami Anna Kuwana, Hao San
Division of Electronics and Informatics Gunma University
IPS04 Analog and Power
3rd International Conference on Technology and Social Science
May 8, 2019
Keynote Lecture 03
Full Body
Outline
● Motivation for Complex Signal Processing Research
● RC Polyphase Filter: Transfer Function
● RC Polyphase Filter: Flat Passband Gain Algorithm
● RC Polyphase Filter and Hilbert Filter
● Active Complex Bandpass Filters
● Conclusion
Outline
● Motivation for Complex Signal Processing Research
● RC Polyphase Filter: Transfer Function
● RC Polyphase Filter: Flat Passband Gain Algorithm
● RC Polyphase Filter and Hilbert Filter
● Active Complex Bandpass Filters
● Conclusion
Why My Research for Complex Signal Processing ?
About 15 years ago
at IEEE International Solid-State Circuits Conference San Francisco, CA
The most prestigious conference in IC design
Katholieke Universiteit Leuven (KU Leuven), Belgium
World top research group in analog IC design
Some simple circuit
with curious characteristics presentation
However,
I could not understand its principle
Basics of Complex Signal
2 real signals: I, Q
V signal = I + jQ Complex Signal
V image = I – jQ Image
I = [V signal + V image ]/2
Q = [V signal – V image ]/(2 j)
j ・ j = -1
Basic complex signal processing blocks
Outline
● Motivation for Complex Signal Processing Research
● RC Polyphase Filter: Transfer Function
● RC Polyphase Filter: Flat Passband Gain Algorithm
● RC Polyphase Filter and Hilbert Filter
● Active Complex Bandpass Filters
● Conclusion
Research Goal of First Research
• To establish systematic design and analysis methods of RC polyphase filters.
• As its first step,
to derive explicit transfer functions of the 1st-, 2nd- and 3rd-order
RC polyphase filters.
Features of RC Polyphase Filter
• Its input and output are complex signal.
• Passive RC analog filter
• One of key components in wireless transceiver analog front-end
- I, Q signal generation - Image rejection
• Its explicit transfer function
has not been derived yet.
First-order RC Polyphase Filter
R 1 C 1
I in+
Q in+
I in- Q in -
I out+
Q out+
I out-
Q out-
C 1 C 1
C 1
R 1 R 1 R 1
I: In-Phase, Q: Quadrature-Phase
Differential Complex Input: V in = I in + j Q in
Differential Complex Output: V out = I out + j Q out
I, Q Signal Generation
-2.5 -2 -1.5
-1 -0.5
0 0.5 1 1.5 2 2.5
0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4
voltage [V]
time [us]
I out Q out
1 1
1 C
LO = R
Polyphase Filter
Qin = 0
I in = cos (ω
LOt)
I out = A cos (ω
LOt+θ)
Q out = A sin (ω
LOt+θ)
Single cosine Cosine, Sine signals
Cosine, Sine Signals in Receiver
IF Analog Bandpass
Filter
cos(ω
0t)
I
Q
AD Converter
- sin(ω
0t)
analog digital AD
Converter
Problem when ω LO ≠ 1 /R 1 C 1
-3 -2 -1 0 1 2 3
0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4
voltage [V]
time [us]
Q out I out
1 1
2 C
LO = R
-2.5 -2 -1.5 -1 -0.5 0 0.5
1 1.5
2 2.5
0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4
voltage [V]
time [us]
I out Q out
1 1
1 C
LO = R
2 nd -order RC Polyphase Filter
I out Q out
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4
voltage [V]
time [us]
1 1
2 C
LO = R
The problem of large difference between I out , Q out amplitudes can be alleviated
R1
R1 R1
R1 C1
C1
C1
C1 C2 C2
C2
C2
R2
R2 R2
R2
Iin+
Qin+
Iin- Qin-
Iout+
Qout+
Iout-
Qout-
3 rd -order RC Polyphase Filter
1 1
2 C
LO = R
I out Q out
The amplitude
difference problem is further alleviated.
R1
R1 R1
R1 C1
C1
C1
C1 C2
C2 C2
C2
R2 R2 R2 R2
C3 R3 C3 R3
R3
C3
R3
C3
Iin+
Qin+
Iin- Qin-
Iout+
Qout+
Iout-
Qout-
Pure I, Q Signal Generation
With
3 rd -order harmonics.
Without
3 rd -order harmonics.
Simulation of 3 rd -order Harmonics Rejection
-8 -6 -4 -2 0 2 4 6 8
20 22 24 26 28 30 32 34 36 38 40
voltage [V]
time [ns]
Iout Qout
) (
sin )
sin(
) (
) (
cos )
cos(
) (
3 3
t a
t t
Q
t a
t t
I
LO LO
in
LO LO
in
+
=
+
=
) sin(
) (
) cos(
) (
+
=
+
=
t A
t Q
t A
t I
LO out
LO out
-8 -6 -4 -2 0 2 4 6 8
0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 0.38 0.4
voltage [V]
time [us]
Iin Qin
1 1
3 1
C
LO = R
Image Rejection Filter
Polyphase Filter
I in =
(A + B) cos(ωt)
Q in =
(A - B) sin(ωt)
I out = Acos(ωt) Q out =
Asin(ωt)
t j t
j Be
Ae + − Ae j t
image
Approach (1)
Complex Input: V in = I in + j Q in Complex Output: V out = I out + j Q out Complex Signal Processing
There is NO physical complex signal.
It is only defined mathematically.
Complex Signal Processing is NOT Complex.
Gauss plane
Complex Transfer Function
• Complex Signal Theory
• Complex input
• Complex output
• Complex
Transfer Function
) (
) ) (
(
j V
j j V
G
in
= out
out out
out
in in
in
Q j
I j
V
Q j
I j
V
+
=
+
= ) (
) (
Signals in RC Polyphase Filter
) ( )
( )
(
) ( )
( )
(
) ( )
( )
(
) ( )
( )
(
) ( )
( )
(
) ( )
( )
(
t jQ
t I
t V
t jQ t
I t
V
t Q
t Q
t Q
t I
t I
t I
t Q
t Q
t Q
t I
t I
t I
out out
out
in in
in
out out
out
out out
out
in in
in
in in
in
+
=
+
=
−
=
−
=
−
=
−
=
− +
− +
− +
− +
Differential signal
Complex signal
R 1
C 1
I in+
Q in+
I in-
Q in-
I out+
Q out+
I out-
Q out-
C 1
C 1
R 1
R 1
R 1
C 1
Transfer Function of RC Polyphase Filter
1 2 1
) (
1 ) 1
(
. 1
) 1 (
RC j RC
G
RC j
j RC G
+
= + +
= +
•Transfer Function
•Gain
RC
− 1
Asymmetric
Explanation of I, Q Signal Generation by G 1 ( jω)
𝑄 𝑖𝑛 𝑡 ≡ 0, 𝐼 𝑖𝑛 𝑡 = cos 𝜔𝑡 𝑉 𝑖𝑛 𝑡 = 𝐼 𝑖𝑛 𝑡 + 𝑗 𝑄 𝑖𝑛 𝑡 = cos 𝜔𝑡 = 1
2 [𝑒 𝑗𝜔𝑡 + 𝑒 −𝑗𝜔𝑡 ]
𝐺 1 𝑗𝜔
𝜔= 1 𝑅𝐶
= 0 , 𝐺 1 𝑗𝜔
𝜔= 1 𝑅𝐶
= 2 , ∠𝐺 1 𝑗𝜔 = − 𝜋 4 𝑉 𝑜𝑢𝑡 𝑡 = 1
2 [ 𝐺 1 𝑗𝜔 𝑒 𝑗 𝜔𝑡+∠𝐺 1 𝑗𝜔 + 𝐺 1 −𝑗𝜔 𝑒 𝑗(−𝜔𝑡+∠𝐺 1 (−𝑗𝜔)) ]
= 2
2 cos 𝜔𝑡 − 𝜋
4 + 𝑗 2
2 sin(𝜔𝑡 − 𝜋
4 )
Here 𝑮 𝟏 −𝒋𝝎 𝒆
𝒋(−𝝎𝒕+∠𝑮
𝟏(−𝒋𝝎)) ] = 𝟎
Outline
● Motivation for Complex Signal Processing Research
● RC Polyphase Filter: Transfer Function
● RC Polyphase Filter: Flat Passband Gain Algorithm
● RC Polyphase Filter and Hilbert Filter
● Active Complex Bandpass Filters
● Conclusion
Transfer Function of 2 nd -order RC Polyphase Filter
Transfer Function
Derivation is very complicated, so we used “Mathematica.”
Gain |G2(jω)|
characteristics
𝐺 2 𝑗𝜔 = (1 + 𝜔𝑅 1 𝐶 1 )(1 + 𝜔𝑅 2 𝐶 2 )
1 − 𝜔 2 𝑅 1 𝐶 1 𝑅 2 𝐶 2 + 𝑗𝜔(𝐶 1 𝑅 1 + 𝐶 2 𝑅 2 + 2𝑅 1 𝐶 2 )
Gain of 3 rd -order RC Polyphase Filter
Gain:
Phase:
Gain characteristics
) (
) )) (
( tan(
) (
) (
) ) (
(
3 3 3
2 3
2 3
3 3
j D
j j D
G
j D
j D
j j N
G
R I
I R
−
=
= +
Need for Flat Passband Gain Algorithm
Transfer Function
Gain |G2( jω)|
characteristics
𝐺 2 𝑗𝜔 = (1 + 𝜔𝑅 1 𝐶 1 )(1 + 𝜔𝑅 2 𝐶 2 )
1 − 𝜔 2 𝑅 1 𝐶 1 𝑅 2 𝐶 2 + 𝑗𝜔(𝐶 1 𝑅 1 + 𝐶 2 𝑅 2 + 2𝑅 1 𝐶 2 )
We need flat passband gain
Four Design Parameters
4 parameters : 𝑅 1 , 𝑅 2 , 𝐶 1 , 𝐶 2
𝜔 1 = 1
𝑅 1 𝐶 1 , 𝜔 2 = 1
𝑅 2 𝐶 2 , 𝑋 = 1
𝑅 2 𝐶 1 , 𝑌 = 1 𝑅 1 𝐶 2
4 constraints
Two Constraints from Filter Spec.
● 2 zeros : − 𝜔 1 = − 1
𝑅 1 𝐶 1 , − 𝜔 2 = −1
𝑅 2 𝐶 2
are given from the filter specification.
Proposed Algorithm Uses Third Constraint
● We use the third constraint 𝑋 = 1
𝑅 2 𝐶 1
for passpand gain flattening.
Nyquist Chart of G 2 ( jω)
|G
2(jω
1)|
|G
2(jω
2)|
Gain characteristics |G 2 (jω)| Nyquist chart of G 2 (jω)=X(ω)+j Y(ω)
Y( ω)
X(ω)
|G 2 (jω 1 )|=|G 2 (jω 2 )|
But in general |G 2 (jω 1 )|=|G 2 (jω 2 )|=|G 2 (j√ω 1 ω 2 )|
|G2(j√ω1ω2)|
Our Idea for Flat Passband Gain Algorithm
If we make|G 2 (jω 1 )| =|G 2 (jω 2 )| =|G 2 (𝑗 𝜔 1 𝜔 2 )|, Passband gain becomes flat from ω 1 to ω 2.
Gain characteristics |G 2 (jω)| Nyquist chart of G 2 (jω)=X(ω)+j Y(ω)
|G
2(jω
1)|
|G
2(jω
2)|
|G2(j√ω1ω2)|
Outline
● Motivation for Complex Signal Processing Research
● RC Polyphase Filter: Transfer Function
● RC Polyphase Filter: Flat Passband Gain Algorithm
● RC Polyphase Filter and Hilbert Filter
● Active Complex Bandpass Filters
● Conclusion
Research Objective
Analyze RC polyphase filter
We found that relevance between
RC Polyphase Filter Analog
Complex (I, Q) input
Hilbert Filter Digital
Real part (Vin) input
Vin
Iout+
Qout+
Hilbert Filter
◼ Characteristics
・ Hilbert transform
・ 1 input and 2 outputs
・ It is often implemented in digital filter
Phase
Cosine, Sine Generation with Hilbert Filter
Gain
cos(ωt)+jsin(ωt)
cos(ωt)-jsin(ωt)
cos(ωt)+jsin(ωt)
2 cos(ωt)
cos(ωt) sin(ωt)
x
-ω ω ω component
-ω component
ω component
Hilbert filter
Hilbert Transform
Complex signal from real signal 𝑥 (𝑡) 𝑥 𝑡 → 𝑥 𝑡 + 𝑗𝑦 𝑡
Hilbert transform
𝑦 𝑡 = 1
𝜋 න
−∞
∞ 𝑥(𝜏)
𝑡 − 𝜏 𝑑𝜏 = 𝑥 𝑡 ∗ 1 𝜋𝑡
Impulse response Fourier Transform
ℎ 𝑡 = 1
𝜋𝑡 𝑯 𝝎 = ቊ −𝒋 (𝝎 ≥ 𝟎) 𝒋 (𝝎 < 𝟎)
Frequency characteristic 𝐻 (ω)
𝒀 𝝎 = 𝑯 𝝎 𝑿(𝝎) = ቊ −𝒋𝑿 𝝎 (𝝎 ≥ 𝟎) 𝒋𝑿 𝝎 (𝝎 < 𝟎)
𝐹𝑜𝑢𝑟𝑖𝑒𝑟
David Hilbert 1862-1943
Phase
余弦波の複素平面での解釈
「現世」と「あの世」
ヒルベルトフィルタで負周波数成分のカット
ヒルベルトフィルタで負周波数が見える
1 st order RC Polyphase Filter : Analysis
R1
R1
R1
R1 C1
C1
C1
C1
Iin+
Qin+
Iin-
Qin-
Iout+
Qout+
Iout-
Qout-
Pass band Stop band
: Transfer function
𝐻 1 𝑗𝜔 = 1 + 𝜔𝑅 1 𝐶 1 1 + 𝑗𝜔𝑅 1 𝐶 1
1
1 st order RC Polyphase Filter : Gain and Phase
H 1re − H 1im H 1re + H 1im
π
2 phase lag
π
2 phase lead 𝐻 1𝑟𝑒 𝑗𝜔 = 𝐻 1 𝑗𝜔 + 𝐻 1 ∗ (−𝑗𝜔)
2 = 1
1 + 𝑗𝜔𝑅 1 𝐶 1 𝐻 1𝑖𝑚 𝑗𝜔 = 𝐻 1 𝑗𝜔 − 𝐻 1 ∗ (−𝑗𝜔)
2 = −𝑗 𝜔𝑅 1 𝐶 1 1 + 𝑗𝜔𝑅 1 𝐶 1
𝐻 1 𝑗𝜔 = 𝐻 1𝑟𝑒 𝑗𝜔 + 𝑗𝐻 1𝑖𝑚 (𝑗𝜔)
ー H 1re
ー H 1im
ー H 1
1 st order case Analysis Results
Gain : Hilbert filter only at zero Phase : Completely Hilbert filter
Hilbert filter
RC Polyphase Filter
Results : 2 nd to 4 th RC Polyphase Filter
Gain Phase
2 nd
3 rd
4 th
ー H
1reー H
1imー H
1Analysis Results and Consideration
1 st to 4 th order RC Polyphase Filter Analysis results
Prove for general n-th order case (n = 1, 2, 3, 4, 5, ...)
Gain : Hilbert filter only at zero
Phase : Completely Hilbert filter
Order and Gain
-20 -10 0 10 20
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
相対角周波数 [rad/s] (RC=1/(2π))
|H|
-20 -10 0 10 20
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
相対角周波数 [rad/s] (RC=1/(2π))
|H|
-20 -10 0 10 20
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
相対角周波数 [rad/s] (RC=1/(2π))
|H|
-20 -10 0 10 20
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
相対角周波数 [rad/s] (RC=1/(2π))
|H|
1
st2
nd3
rd4
thThe higher orders,
Increase number of zeros
|𝐻𝑟𝑒| and |𝐻𝑖𝑚| becomes close in wide range
Close to ideal Hilbert transform
Order and Phase
Phase characteristic is not changed
There is always 90 ° phase difference
1
st2
nd3
rd4
thFulfill Hilbert transform in full range
-20 -10 0 10 20
-π π/2 0 π/2 π
相対角周波数 [rad/s] (RC=1/(2π))
∠H
-20 -10 0 10 20
-π π/2 0 π/2 π
相対角周波数 [rad/s] (RC=1/(2π))
∠H
-20 -10 0 10 20
-π π/2 0 π/2 π
相対角周波数 [rad/s] (RC=1/(2π))
∠H
-20 -10 0 10 20
-π π/2 0 π/2 π
相対角周波数 [rad/s] (RC=1/(2π))
∠H
Conclusion
RC polyphase filter is approximation of ideal Hilbert filter for complex input signal
-20 -10 0 10 20
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
相対角周波数 [rad/s] (RC=1/(2π))
|H|
RC Polyphase Filter Hilbert Filter
Phase
Gain Gain Phase
-20 -10 0 10 20
-π π/2 0 π/2 π
相対角周波数 [rad/s] (RC=1/(2π))
∠H
Outline
● Motivation for Complex Signal Processing Research
● RC Polyphase Filter: Transfer Function
● RC Polyphase Filter: Flat Passband Gain Algorithm
● RC Polyphase Filter and Hilbert Filter
● Active Complex Bandpass Filters
● Conclusion
Gm : Transconductance
Input voltage: V in
Output current : I out I out =g m V in
V in g m I out I + I -
I b V in
+ -
dimension of g m
R 1
I out =I + -I -
=g m V in
Complex Bandpass Gm-C Filter
I in
Q in
C
C
g m
- +
V Iout
V Qout g 0
g 0
C s g C
s g
g
jg sc
g jQin Iin
jV V
m
m
Qout Iout
2 0 2 2
0 2
0
+ 2 +
+
−
= +
+
+
+
Gain of Complex Bandpass Gm-C Filter
0 ω
C gm C
g C
gm +
0C g C
gm −
00
1 2 1
g
0
1 g
Center Freq.
C gm
g 0
Q = gm
Gain
Complex Bandpass Active RC Filter
I out+
Q out- I in
I out-
Q out+
Q in
+ -
+ -
g
2g
2g
3g
3g
1g
1C
C
3 2
1