電気回路学 Ⅱ
コミュニケーションネットワークコース
5
セメ山田 博仁
ラプラス変換における初期条件の扱
1.
キャパシタの初期電荷い
2.
コイルの初期電流C
i(t)
q(0)
i(t) L
i(0)
t t
t
q i t dt
dt C t i dt
t C i
dt t C i
t
v
0 00
1 ( 0 ) ( )
) ( )
1 ( )
1 ( ) (
+
の意味 q i t dt sC q I sC s
t C
v
t( 0 ) ( )
) ( )
0 1 (
)
(
£
0
£
これをラプラス変換すると、
v(t)
v(t)
流れる電流と両端の電圧との関係は、
初期電荷
q(0)
により発生する 電圧がv(t)
と同一方向なら+
、逆なら ‒
dt t L di t
v ( )
) (
) 0 ( )
( )
( t sLI s Li
v
£
これをラプラス変換すると、
流れる電流と両端の電圧との関係は、
初期電流
i(0)
がi(t)
と同一方向なら ‒、逆なら +過渡関数波
過渡関数波とは
?
単位ステップや単位インパルスを、時間微分或いは積分した関数で表され る一連の波形を過渡関数波と呼ぶ。
単位ステップと単位インパルス
0 t 1
a
0 a t a
1
1( )
1( )
1( 0 ) 1 0 1 u
0( t )
s s u
t u dt s
t
du
£ ££
図
(a)
の波形を時間t
で微分すると図(b)
の波形を得る。(a)
(b)
a → 0
の極限を考えると、図(a)
の波形は単位ステップ
u
–1(t)
となり、図(b)
の波形は単位インパル スu
0(t)
となる。) ) (
(
0
1
u t
dt
t
du
即ち、過渡関数波
単位ダブレット
0 2a t a
1
a
(a) (b)
a → 0
(d) 0 t a 1
2a a a
1
a
2
)
~ (
1
t u t
図
(a)
の三角波を時間微分すると、図(b)
のような正および負の方形波が続 いて現れる波形となる。これを で表せば、 の時間積分は0
となるが、 の1
次モーメント を考えると、図(d)
のように その時間積分は − 1 となることが分かる。そこで、a → 0
の極限を考えて、 を考えると、図
(c)
のように高さは無限に高く、幅が 無限に小さい正と負のインパルスが、t = 0
の時刻に同時に存在する波形と なる。これを単位ダブレットと呼び、その1
次モーメントは − 1となる。
(c) 0 t
2
1 a
2a a
)
~ (
1
t u
2
1
a
)
~ (
1
t
u ~ ( )
1
t
u ~ ( )
1
t u )
~ (
1
t u t
) ( )
~ (
1
1
t u t
u
0 t
)
1
( t +∞ u
–∞
単位ダブレット
1
)
1
(
tu t dt
s
dt t t du
u
( )
)
(
01 £
£
また、単位ダブレットは単位インパルスを時間微分した ものであるから、そのラプラス変換は、 となる。
過渡関数波
高次の特異波形
単位インパルス
u
0(t)
をk
回微分した特異な関数をu
k(t)
で表す。それ は、正負のインパルスが時刻t = 0
に同時にk + 1
個 発生する波形であ る。
k n k k dt n
t u
t
n k k! ) 1 ( ) 0
(
u
k( t ) s
kまた、ラプラス変換は、 となる。£ その
k
次モーメントは、であり、有限確定値をとる。
0 t
)
1
( t +∞ u
–∞
単位ダブレット
t
で微分t
で微分0 t
)
2
( t +∞ u
–∞
0 t
)
3
( t +∞ u
–∞
t = 0
で同時過渡関数波
単位ランプ
0 t 1
u
–1(t)
単位インパルス
u
0(t)
をk
回積分して得られる関数をu
–k(t)
で表す。1
回積分したものは、図(a)
の単位ステップu
–1(t)
で、2
回積分したものは 図(b)
に示すように、時刻t = 0
から直線的に増加する波形であり、3
回 積分したものは図(c)
に示すように、時刻t = 0
から放物線的に増加する 波形となる 。これら一群の関数を単位ランプと呼ぶ。(a)
単位ステップu
-1(t) (b)
単位半無限ランプu
–2(t) 0 t
1
u
–2(t)
1
(c)
単位放物線ランプu
–3(t) 0 t
1
u
–3(t)
1 ,
2 , 1 ),
)! ( 1 ) (
(
11
u t k
k t t
u
k k
過渡関数波
単位ランプのラプラス変換は、
u
k( t )
£
0tu
k1( t ) dt
£
0t
0tu
0( t )( dt )
k s
k£ となる。
例
5.3.1
例5.3.2
f(t)
がt = a
で連続なら、) ( )
( )
0
( t a f t dt f a
u
の関係が成り立つ。
u
-1(t)sinωt u
-1(t) sinωt
時刻
t = 0
に突然現れる正弦波単位インパルス
u
0(t)
を用いて過渡関数波
繰り返す波形のラプラス変換
−∞
から時刻t = 0
までf(t) = 0
で、t > 0
では周期T
をもって同じ波形 が繰り返されるようなとき、その波形f(t)
を、0 < t < T
の1
周期の間でのみ
f(t)
に等しく、それ以外の全ての時刻t
では0
になる波形f
0(t)
をもって表せば、
f ( t ) f
0( t ) f
0( t T ) f
0( t 2 T )
となる。従って、ラプラス変換F(s)
は、sT
sT sT
sT sT
st st
st st
e s F
s F e
e e
s F e
s F s
F
dt e T t
f dt
e T t f dt
e t f dt
e t f s
F
1
) (
) ( ) 1
( )
( )
( )
(
) 2 ( )
( )
( )
( )
(
0
0 2
2 0
0 0
0 0
0 0
0 0
0
ただし、
F
0( s )
0f
0( t ) e
stdt
0Tf ( t ) e
stdt
F
0(s)
を、ウェイデリッチによる定常ラプラス変換と呼ぶ。あるいは、
である。e
sTt t f
f
1
) ) (
(
£ 0£
f
0(t)
0 T 2T 3T t
f(t)
過渡関数波
例
5.4.1
図に示すように、
t < 0
で0 、 t >
0
では方形波が繰り返すような波形 のラプラス変換F(s)
は、
T t t
t T t
f
, 2 0 0
0 2 1
)
0
(
として、0 t 1
2
T
TT
2
3
2TT
2 5
2 20 2
0
0 0
0
1 1
) ( )
(
sT T
T st st
st
e
s s
dt e e dt
e t f s
F
従って、
2 2
1 1 1
) 1
(
TsT s sT
e e s
s s e F
f
0(t)
過渡関数波
展開定理
F(s)
のラプラス逆変換を求めるにあたり、 F(s) を部分分数に展開し、展開式の各項についてラプラス逆変換するのが便利。例えば、
at
at
e
e a a a a
s a
s a
a s s a a
s s
1 1 1 1 1 1 1 1 1 1 1
) (
1
1 1 11 £ £ £
£
(1) F(s)
が1
位の極のみからなるときn n
s s
C s
s C s
s C s
s s C
F
3 3 2
2 1
)
1(
と書ける。ここで、
s
j(j = 1, 2, ‥‥ , n)
はF(s)
の1
位の極であり、C
j(j = 1, 2, ‥‥ , n)
は 極s
j の留数である。
sj
j s
j
s s F s
C ( )
従って、
nj s s
st j
n
j
t s
j j
j
s s F s e
e C s
F t
f
1 1
1
( ) ( )
)
(
£t s j
e
js
s
1
1
£ より、 F(s) のラプラス逆変換は、
となる。
過渡関数波
展開定理
(2) F(s)
が2
位以上の極をもつとき ( )
)
(
11 1 2
1 13 1
1 12 1
11 1
1 1
1
F s
s s
C s
s C s
s C s
s s C
F
k k k k
と書ける。ここで、
F
1(s)
はもはや、s
1 に極を持たない有理関数であり、C
1j(j = 1, 2,
‥‥ , k
1)
は定数である。従って、
( )
!
)
! ( 1
! 2
! 2
! ) 1
(
1 1
1 1
1
1 1 1
1 1 2
2 1 1
2 12
1 1 11
1 1 1
1 1 1
1 1
1
s j F
k e C t
s F e
t C t C
k C C t k
C t t
f
k
j
t s j k j
t s k k
k k
k
£
£となる。
F(s)
を部分分数に展開(s = s
1 でローラン展開)
して、s
1 以外の極s
2, s
3, ‥‥ , s
n についても、F
1(s)
について行う。回路網関数
微積分方程式とラプラス変換
静止の状態にある
(
全ての初期条件を0
とした)
回路に励振を加えたとき、応答のラプラス変換と励振のラプラス変換との比を回路網関数という。
(
応答のラプラス変換) = (
回路網関数)×(
励振のラプラス変換)
全ての初期条件を
0 (i(0) = 0, q(0) = 0)
と置いてしまえば、励振e(t)
および 応答i(t)
は、それぞれのラプラス変換E(s)
およびI(s) 、 (
ほぼ)1 、 1
対 応する。従って、e(t)
およびi(t)
で考える代わりにE(s)
およびI(s)
で考 えて、これらラプラス変換したものも励振および応答と呼んでいる。また、全ての初期条件を
0
としたとき、Z(s) = E(s)/I(s)
をインピーダンス関数、
Y(s) = I(s)/E(s)
をアドミタンス関数と呼んでいたが、より一般的には次のように定義する。
回路網関数
(network function)
の代わりにシステム関数(system function)
、伝達関数
(transfer function)
などと呼ぶこともある。回路網関数は、対象としている回路網の構造を与えれば一意に定まる。
回路網関数
回路網関数には、
(
電圧)/(
電流)
を表すインピーダンス関数(impedance function) 、 (
電流)/(
電圧)
を表すアドミタンス関数(admittance function)
、 さらにまた励振と応答が同じ節点対(
端子対)
で測られるとき駆動点関数(driving-point function)
、異なる節点対(
端子対)
で測られるとき伝達関数と呼ばれる。
V
1(s)
回路網V
2(s)
I
1(s) I
2(s)
例えば以下の回路網において、
V
1(s)/I
1(s), V
2(s)/I
2(s)
は、駆動点インピーダンス関数I
1(s)/V
1(s), I
2(s)/V
2(s)
は、駆動点アドミタンス関数V
1(s)/I
2(s), V
2(s)/I
1(s)
は、伝達インピーダンス関数I
1(s)/V
2(s), I
2(s)/V
1(s)
は、伝達アドミタンス関数 である。複素記号演算との関 係
回路網関数を
H(s)
とすると、s
をjω
で置き換えたH(jω)
は、複素記号 演算で得られるインピーダンスやアドミタンスと一致する。即ち、ラプラス 変換による演算で定義される回路網関数H(s)
は、複素記号演算で定義され る回路網関数H(jω)
を拡張したもので、s ↔ jω
で相互に置き換わる。
dt j s
s j
dt d
e A e
A t
A
j t st1 1
ˆ ) ˆ
cos(
(
ラプラス変換)
(
初期条件の導入) E(s ) ( (Y(s)E(s))
代数演算) I(s) (
ラプラス逆変換)
(
微分方程式の標準的解法)
e(t ) i(t )
周波数域解析 時間域解析
微分方程式またはラプラス変換による時間域解析のプロセスを下図に示す。
線形電気回路の時間域解析
時間域解析と周波数域解 析
線形電気回路の解析にラプラス変換を用いれば、初期条件も導入しながら、
多くの関数のラプラス変換がラプラス変換表を用いて機械的に行える。その 後、
s
関数の代数演算によって応答のラプラス変換を求め、ラプラス変換表 を用いてラプラス逆変換を行えば、時間応答が求められる。時間域解析と周波数域解析
時間
t
の関数としての励振e(t)
に対する回路網の応答i(t)
を求めること を、時間域解析(time domain analysis)
と呼んでいる。これに対し、e(t),
i(t)
のラプラス変換E(s)
とI(s)
の関係を求めることを周波数域解析(frequency domain analysis)
と言う。周波数域解析では、一般的には初期条件を考慮しない。全ての初期条件を
0
として扱う。ラプラス変換による演算法は、ヘビサイド
(Oliver Heaviside)
によって導入 されたヘビサイドの演算子法を数学的に明確にする過程で変形されたもの。周波数域解析における重ねの
周波数域解析における重ね合わせの理と初期条件
理
周波数域解析で初期条件を扱う必要のある場合には、初期条件に関連した 項を強制振動項と同格に扱い、強制振動の一成分であると考える。
RLC
直列回路を例に見てみると、回路方程式のラプラス変換は、 ( ) ( 0 ) 1 ( ) ( 0 )
) ( )
( I s q
i sC s
sI L s RI s
E
で表されるから、
) 0 ( )
( 1 )
(
) 0 ) (
( 1 )
(
) ( )
( 1 ) (
3 2
1
Li s
sC I R
sL
sC s q
sC I R
sL
s E s
sC I R
sL
が成り立つ。
) ( ) ( )
( )
(
2 31
s I s I s I s
I
即ち、
( 0 ) , ( 0 ) ),
( Li
sC s q
E
を各々独立した励振と見なした場合の応答)
( ), ( ),
(
2 31