応用数学 - ラプラス変換

提供: MochiuWiki : SUSE, EC, PCB

概要

ラプラス変換(Laplace transform)は、t0 で定義された関数 f(t) から無限積分を用いて、新しい関数 F(s) を作り出す積分変換である。
ラプラス変換は、微分方程式を代数演算に帰着させて解くための方法を与える。

ラプラス変換を使用すると、ある種の微分方程式に関しては、ややこしい積分を使用せずに代数的に解くことが可能になる。

ラプラス変換の名前は、ラプラス(Pierre-SimonLaplace, 1749 - 1827)に因む。
ラプラスは微分方程式の解法としてラプラス変換を導入したわけではない。

「ラプラス変換を使用した微分方程式の解法」は、ヘビサイド(Oliver Heaviside, 1850 -1925)によるものである。
ヘビサイドは演算子を使用した微分方程式の解法を発見した。
この演算子の数学的基礎付けが他の数学者達によって行われ、ラプラス変換との関係が明らかにされた。


ラプラス変換の定義

定義: f(t)t0 で定義された関数とする。
この時、以下のsの関数 F(s)f(t) のラプラス変換と言い、[f(t)] または L[f]と書く。

F(s)=[f(t)]

F(s)=0+f(t)estdt

f(t) を原関数(primitive function)、F(s) を像関数(image function)または f(t) のラプラス変換と呼ぶ。
F(s) のsは実数または複素数である。


f(t)t>0 で定義されていて、t=0 では定義されていない場合でも、f(0)の値を適当に決めて t0 で定義されているとしてよい。

f(t) から F(s) に変換する式の説明
f(t)est を乗算して f(t)est を作る時、s は定数とみなしてよい。
次に、この関数を 0+ の範囲で積分する。(積分は変数tで行うので、sはまだ定数とみなしてよい)
無限積分は定積分なので、積分の結果、変数tは消える。(ここで定数と見なしていたsは残る)
これを、sの関数 F(s) とする。

ラプラス変換の定義式の F(s) は、右辺の無限積分が収束する場合のみ定義される。
右辺の無限積分はsの値によって収束する場合と収束しない場合が考えられる。

ラプラス変換は、変数tの関数の集合(原関数の集合)から、変数sの関数の集合(像関数の集合)への変換規則を表している。


ラプラス変換の定義に関する例題

ラプラス変換の定義に従って、関数 f(t)=1 のラプラス変換 [f] を求めよ。

[f]=0+f(t)estdt=0+estdt


sの値により場合分けする。

s=0 の時

[f]=0+e0t˙dt=0+dt=limβ+0βdt=limβ+[t]0β=limβ+β=+


s0 の時

[f]=0+estdt=limβ+0βestdt=limβ+[ests]0β=limβ+esβ1s=limβ+1esβs={+(s<0 の と き )1s(s>0 の と き )


以上より、s>0 の時のみ L[f] は収束し、次式となる。

[1]=1s(s>0)



ラプラス変換の存在条件

定理:
関数 f(t)t0 において区分的に連続であるとする。
この時、f(t) に対し、以下を満たす定数Mとαが存在すれば、s>α であるsについてラプラス変換 L[f] が存在する。
|f(t)|Meαt

説明:
この定理の条件を満たす関数を指数α位の関数と言う。
f(t) が指数α位のとき、ラプラス変換の定義式は s>α の範囲で収束して、ラプラス変換が存在する。



基本的な関数のラプラス変換公式

下表のような原関数 f(t) と像関数 F(s) の対応関係を変換表として使用すれば、ラプラス変換を機械的に実行できる。

基本的な変換表
f(t) F(s) 条件
1 1s s>0
t 1s2 s>0
tn n!sn+1 s>0,n=0,1,2,
tα Γ(α+1)sα+1 α>1,s>0
eat 1sa s>α
sinat as2+a2 s>0,a:実定数
cosat ss2a2 s>0,a:実定数
単位階段関数 U(ta) eass s>0,a>0
デルタ関数 δ(t) 1 <s<


法則
f(t) F(s) 条件
線形法則
sinhat
as2a2 s>a0
線形法則
coshat
ss2a2 s>a0
相似法則
f(at)
1aF(sa) s>αa
第1移動法則
eλtsinat
a(sλ)2+a2 s>λ
第1移動法則
eλtcosat
sλ(sλ)2+a2 s>λ
第2移動法則
U(tλ)f(tλ)
eλsF(s) s>a


※備考
sinhat=eateat2
coshat=eat+eat2



ラプラス変換公式を使用した例題

ラプラス変換公式を使用して、以下の関数のラプラス変換を求めよ。

(1) [t2]
(2) [t]
(3) [sinωt]
(4) [cost3]
(5) [et]


  • (1)の解答
    [t2]=2!s(2+1)=2s3(s>0)
  • (2)の解答
    [t]=[t12]=Γ(12+1)s12+1=12Γ(12)s32=π2ss(s>0)
    参考: Γ(12)=π,Γ(α+1)=αΓ(α)(α>0)
  • (3)の解答
    [sin2t]=2s2+22=2s2+4(s>0)
  • (4)の解答
    [cost3]=ss2+(13)2=ss2+19=9s9s2+1(s>0)
  • (5)の解答
    [et]=1s(1)=1s+1(s>1)



線形法則

定理:
[f(t)][g(t)] がそれぞれ s>α1,s>α2 において存在する時、
任意の定数 a,bについて以下が成立する。

[af(t)+bf(t)]=a[f(t)]+b[g(t)](s>max(α1,α2))


例題:
線形法則とラプラス変換公式を使用して、以下のラプラス変換を求めよ。
[3et+2e2t]

解答:
線形法則より以下が成立する。
[3et+2e2t]=3[et]+2[e2t]

したがって、変換公式を使用すると以下が求められる。
 与 式 =3[et]+2[e2t]=31s+1+21s2=3s+1+2s2

sの範囲は s>1 かつ s>2 より、 s>2 となる。



相似法則

定理:
[f(t)]=F(s)(s>α) の時、正の定数aについて以下が成立する。
[f(at)]=1aF(sa)(s>αa)

説明:
f(t) の変数tをa倍した f(at) に対応するラプラス変換は、 f(t) のラプラス変換 F(s) の変数を 1a 倍して、さらにそれを 1a 倍したものになる。

例題:
相似法則とラプラス変換の公式 [et]=1s1(s>1) を使用して、以下の値を求めよ。
[et2]

解答:
この問題は、F(s)=1s1,a=12 の場合である。
[et2]=112F(s12)=21s121=22s1

sの範囲は、s12>1 より、 s>12 となる。



第1移動法則

定理:
[f(t)]=F(s)(s>a) の時、以下が成立する。
[eλtf(t)]=F(sλ)(s>a+λ)

説明:
f(t) に指数関数 eλt を乗算するということは、 f(t) のラプラス変換 F(s)λ だけ平行移動することに対応する。

例題:
第1移動法則およびラプラス変換の公式 [t]=1s2 を使用して、以下の値を求めよ。
[te3t]

解答:
第1移動法則において、 f(t)=t,λ=3 とする。
[t]=1s2 より、
[te3t]=[e3tt]=1(s3)2

ここで、 s3>0 より、 s>3



第2移動法則

定理:
[f(t)]=F(s)(s>a) の時、正の定数 λ について以下が成立する。
[U(tλ)f(tλ)]=eλsF(s)(s>a)

説明:
f(t)λ だけ平行移動させることは、ラプラス変換した側では F(s) に指数関数 eλs を乗算することに対応する。
U(tλ)f(tλ) のグラフは、tλ の時は f(tλ)t<λ の時は 0 となる。


例題:
第2移動法則を用いて以下の値を求めよ。
[U(t2)(t2)2]

解答:
U(t2)(t2)2 は第2移動法則において、 λ=2,f(t)=t2 の場合に対応する。
第2移動法則およびラプラス変換の公式 [tn]=n!sn+1(s>0) を使用して計算する。

[U(t2)(t2)2]=e2s[t2]=e2s2!s2+1(第2移動法則)=2s3e2s(s>0)(変換公式)



微分法則

定理1:
f(x)t0 で定義された連続な関数で、以下を満たす定数Mとαが存在し、微分可能であるとする。
|f(t)|Meαt

この時、df(t)dt が区分的に連続であれば、s>α の時 [df(t)dt] が存在し、以下が成立する。
[df(t)dt]=s[f(t)]f(0)(s>α)

定理1の説明:
f(t) を微分することは、 f(t) のラプラス変換 F(s) をs倍して f(0) を減算することに対応する。
原関数を微分することは、像関数をs倍した後、定数 f(0) を減算することに対応する。


定理2:
t0 において、 f(t),df(t)dt,,d(n1)f(t)dt(n1) は連続な関数で、以下を満たす定数Mとαが存在すると仮定する。
|f(i)(t)|Meαt(i=0,1,2,,n1)

この時、 t0 において dnf(t)dtn が区分的に連続であれば、 s>α の時、 [dnf(t)dtn] が存在し、以下が成立する。
[f(n)(t)]=sn[f(t)]sn1f(0)sd(n2)f(0)dt(n2)d(n1)f(0)dt(n1)(s>α)

定理2の説明:
原関数の世界で微分することは、像関数の世界では「sの多項式を作る」ことに対応する。
この定理の式はラプラス変換を使用して微分方程式を解く際に重要な役割を果たす。


微分方程式を解く際に必要となる公式:
f(t),df(t)dt,s が上記の定理の条件を満たしているとする。
この時、以下が成立する。
(1) [df(t)dt]=s[f(t)]f(0)
(2) [d2f(t)dt2]=s2[f(t)]sf(0)df(0)dt


例題:
f(t)=tet について以下を求めよ。
(1) df(t)dt
(2) f(0)
(3) [f(t)] (微分法則を使用する)

解答:
(1)
積の微分公式より、 df(t)dt=et+tet

(2)
t=0 の時、 f(0)=0×e0=0

(3)
n=1 の微分法則により、 [df(t)dt]=s[f(t)]f(0)
これに、 f(t),df(t)dt,f(0) を代入する。

すると、 [et+tet]=s[tet]0 となる。

[et+tet]=s[tet]0[et]+L[tet]=s[tet] 線 形 性 よ り (s1)[tet]=[et](s1)[tet]=1s1(s>1)[tet]=1(s1)2(s>1)



積分法則

定理:
f(t)t0 において区分的に連続であり、ある定数M, αに対して、 |f(t)|Meαt が成立すると仮定する。
この時、以下が成立する。
[0tf(u)du]=1s[f(t)](s>max(0,α))

説明:
原関数を積分することは、像関数をsで除算することに対応する。

参考:
原関数の微分は、像関数をs倍すること(その後に定数を減算すること)に対応する。


例題:
積分法則を使用して以下の値を求めよ。
[0tsin(3u)du]

解答:
積分法則より、
[0tsin(3u)du]=1s[sin(3t)]=1s3s2+32=3s(s2+9)(s>0)

参考:
[sin(at)]=as2+a2(s>0)



tn積法則

定理:
f(t)は、 t0 において区分的に連続で指数α位の関数とする。
この時、 [tnf(t)](n=1,2,) が存在し、 [f(t)]=F(s) とおくとき、以下が成立する。
[tnf(t)]=(1)ndnF(s)dsn(s>α,n=1,2,3,)

説明:
原関数の世界で「関数を tn 倍する」ことは、像関数の世界では「関数を n 回微分して (1)n 倍する」ことに対応する。

tn 積法則の具体例:
f(t) は、先の定理の仮定を満たしているとする。
この時、以下が成立する。
(1) n= 正 整 数  の場合:
[tnf(t)]=(1)ndnF(s)dsn

(2) n=1 の場合:
[tf(t)]=dF(s)ds

(3) n=2 の場合:
[t2f(t)]=d2F(s)ds2


例題:
tn 積法則を使用して以下の値を求めよ。
[te2t]

解答:
f(t)=e2t とおくと、
[f(t)]=L[e2t]=1s2=(s2)1(s>2)

tn 積法則の n=1 の場合を適用すると以下が得られる。
[te2t]=dF(s)ds=dds(s2)1

ここで、合成関数の微分公式を使用して、
 与 式 =(1)(s2)2=1(s2)2(s>2)



合成積

定義:
f(t),g(t)t0 において区分的に連続であるとする。
この時、以下の関数 (f*g)(t)f(t)g(t) の合成積(convolution)あるいは畳み込みという。

(f*g)(t)=0tf(τ)g(tτ)dτ


合成積の性質
合成積は、以下の法則を満たす。
(1) (f*g)=(g*f) (交換法則)
(2) (f*g)*h=f*(g*h) (結合法則)


  • 合成積 (f*g)(t)=0tf(τ)g(tτ)dτ の直観的意味
    時刻 t のとき、耳に到達する音の波形を h(t) とする。
    h(t) は、その時より τ だけ前の時刻に出た音の波形 g(tτ)f(τ) だけ増幅(あるいは減少)させた波形 f(τ)g(tτ) の影響を受ける。
    時刻 t においては、 τ0 から t までの全ての影響を受けるので、それらを全部重ねあわせた波形 (f*g)(t) が耳に到達する。
    したがって、耳に到達する音の波形は、 h(t)=(f*g)(t) となる。


定理:
f(t),g(t)t0 で区分的に連続な関数であり、かつ指数α位の関数とする。
[f(t)]=F(s),[g(t)]=G(s) とする時、 [(f*g)(t)] が存在し、以下が成立する。

[(f*g)(t)]=F(s)G(s)(s>α)

注意:
f(t)g(t) の(合成積ではない普通の)積のラプラス変換 [f(t)g(t)] については、一般に [f(t)g(t)][f(t)][g(t)] は一致しない。
すなわち、 [f(t)g(t)]F(s)G(s)


例題:
f(t)=t,g(t)=sin(t) について、合成積のラプラス変換公式を使用して [(f*g)(t)] を求めよ。

解答:
ラプラス変換の基本公式より、
[t]=1s2(s>0)
[sin(t)]=1s2+1(s>0)

合成積のラプラス変換公式より、
[(f*g)(t)]=[f(t)][g(t)]=[t][sin(t)]=1s21s2+1(s>0)=1s2(s2+1)(s>0)



ラプラス変換による微分方程式の解法



ラプラス変換の応用

ラプラス変換は、微分方程式を解くために用いられる。
また、システムの安定性や過渡解の解析に用いられる。

主な応用分野

  • 古典的自動制御: 伝達関数を求めるための計算手法。
  • 電気回路: 過渡現象を表す微分方程式の解法。



他の積分変換との比較

フーリエ変換

フーリエ変換(Fourier transform)とラプラス変換では積分範囲と係数部分が異なる。
また、ラプラス変換の変数sに対応する変数が iω (i: 虚数単位)になる。

フーリエ変換は、ラプラス変換とは異なり、周波数分析に使用される。

フーリエ変換は、スイッチをON/OFFした後、時間が十分経過した後の電圧・電流の波形の分析に使用される。
ラプラス変換は、過渡状態の電圧・電流の分析に使用される。(過渡状態とは、スイッチをON/OFFした直後の状態である)

Z変換

Z変換(Z-transform)は、両側ラプラス変換を離散化したものである。
Z変換は、デジタル信号処理で用いられる。

変換を離散化することによりコンピュータで効率よく扱うことができるようになる。

離散フーリエ変換

離散フーリエ変換(DFT: Discrete Fouriertransform)は、フーリエ変換を離散化したものである。
離散フーリエ変換を高速に計算するためのアルゴリズムとして高速フーリエ変換(FFT: Fast Fouriertransform)がある。
これらの技術はデジタル信号処理で使用される。