ラプラス変換

提供: testwiki
ナビゲーションに移動 検索に移動

テンプレート:出典の明記

関数変換

関数解析学において、ラプラス変換(ラプラスへんかん、テンプレート:Lang-en-short)とは、積分で定義される関数空間の間の写像(線型作用素)の一種。関数変換。積分変換の一種。

ラプラス変換の名は18世紀数学者ピエール=シモン・ラプラスにちなむ。

ラプラス変換によりある種の微分積分は積などの代数的な演算に置き換わるため、制御工学などにおいて時間領域の(とくに超越的な)関数を別の領域の(おもに代数的な)関数に変換することにより、計算方法の見通しを良くするための数学的な道具として用いられる。従って、数学の中ではかなり応用寄りの分野である。

フーリエ変換を発展させて、より適用範囲を広げた計算手法である。1899年に電気技師であったオリヴァー・ヘヴィサイドが回路方程式を解くための実用的な演算子を経験則として考案して発表し、後に数学者がその演算子に対し厳密に理論的な裏付けを行った経緯がある。理論的な根拠が曖昧なままで発表されたため、この計算手法に対する懐疑的な声も多かった。この「ヘヴィサイドの演算子」の発表の後に、多くの数学者達により数学的な基盤は1780年の数学者ピエール=シモン・ラプラスの著作にある事が指摘された(この著作においてラプラス変換の公式が頻繁に現れていた)。

フーリエ変換がL^1((-∞,∞))上のゲルファント変換であるのに対しラプラス変換はL^1((0,∞))上のゲルファント変換と説明できる。


これと類似の解法として、より数学的な側面から作られた演算子法がある。こちらは演算子の記号を多項式に見立て、代数的に変形し、公式に基づいて特解を求める方法である。

定義

実数 テンプレート:Math について定義された関数 テンプレート:Mathラプラス変換とは テンプレート:Indent で定義される テンプレート:Mvar の関数 テンプレート:Math のことである。ここで テンプレート:Mvar複素数であり、2 つの実数 テンプレート:Mvar を用いて テンプレート:Math と表すことができる(テンプレート:Mvar虚数単位)。右辺の積分はラプラス積分 テンプレート:En と呼ばれる。これは時間領域から複素平面への写像である。

また、テンプレート:Math として、関数 テンプレート:Math から元の関数 テンプレート:Math を計算することを逆ラプラス変換 テンプレート:En といい、 テンプレート:Indent のように定義されている。ここでcは全ての特異点の実部よりも大きい実数である。右辺の積分はブロムウィッチ積分 テンプレート:En と呼ばれる。これは複素平面から時間領域への写像である。

複素平面を用いたブロムウィッチ積分の解説

これは複素積分となっている。定義通りの積分経路では計算が難しくなるが、閉曲線となるように積分経路を変更して留数を計算することにより簡単に逆ラプラス変換を求める事が可能となる。結果を言えば複素平面上の全ての特異点の留数の総和となる。ここで、テンプレート:Math原関数 テンプレート:Enテンプレート:Math像関数 テンプレート:En という。

ラプラス変換の他の記述の仕方として、次のようなものもある。 テンプレート:Indent

同様に逆ラプラス変換は、次のようにも記述される。 テンプレート:Indent

また、これらの記号を用いた写像 テンプレート:Indent のことも、それぞれラプラス変換逆ラプラス変換と呼ぶ。

普通、ラプラス変換および逆ラプラス変換を行う際には変換表を参照して計算する場合が多いので、前述した定義式にしたがって計算することは少ない。だが場合によっては定義式から計算したほうが簡単なときもある。たとえば逆ラプラス変換をする際に部分分数分解をしなければならない場合、むしろブロムウィッチ積分を計算したほうが早いことも多い。

注:
ラプラス変換は、関数 テンプレート:Math にいったん テンプレート:Math を乗じてからフーリエ変換する操作であると考えることができる(ここで テンプレート:Mathステップ関数である)。
F(s):=F(σ,ω)=θ(t)f(t)eσteiωtdt==s=σ+iω0f(t)estdt

両側ラプラス変換

テンプレート:Main 両側ラプラス変換は積分区間を全実数域へと拡張したもので、以下のように定義される。

F(s)={f(t)}=estf(t)dt

母関数との関係

数列 テンプレート:Mvar の(通常型)母関数

G(an;x)=nanxn

において テンプレート:Math とすると、

G(an;es)=nanesn

となる。 ここで和を積分に変えれば

G(at;es)=atestdt

となり、関数 テンプレート:Mvar のラプラス変換と一致する。この意味においてラプラス変換は母関数の「連続版」とみなすことができる。 こうした理由により、母関数とラプラス変換は同種の性質を満たすことがある。たとえば母関数の性質

G(an;x)G(bn;x)=G(an*bn;x)

はラプラス変換の性質

[f](s)[g](s)=[f*g](s)

に対応する。ここで テンプレート:Math畳み込み積。

性質

ラプラス変換と逆ラプラス変換は互いに他の逆変換である。 テンプレート:Indentここで、テンプレート:Mvar恒等変換を表わす。

線型性

ラプラス変換は線型性を持ち、したがって特に重ね合わせの原理 を用いて計算することが可能である。ラプラス変換が線型性を持つとは、任意の関数 テンプレート:Math に対して テンプレート:Indent が成り立つということである。ただし、テンプレート:Mathテンプレート:Mvar に関係しない定数。逆ラプラス変換も同様に線形性を持ち、 テンプレート:Indent

が成り立つ。したがって、与えられた関数を部分分数分解できるとき、各因子がラプラス変換の表にあるものに合致すれば、その変換が求められる。

相似性

テンプレート:Math のとき、 テンプレート:Indent が成立する。

微分式

時間 テンプレート:Mvar に関する導関数のラプラス変換は多項式の差となって現れる。実際に、一階の導関数をラプラス変換すると以下のように テンプレート:Math(元の式に テンプレート:Math を代入した値)が現れる。

テンプレート:Indent

また、二階導関数の場合は テンプレート:Math に加え、テンプレート:Math における微分係数 テンプレート:Math が現れる。

テンプレート:Indent

これを繰り返すと、一般の テンプレート:Mvar 階の導関数のラプラス変換は以下のようになる。

テンプレート:Indent

積分式

テンプレート:Indent

畳み込み

関数の畳み込みはラプラス変換で積(値ごとの積)に写される。 テンプレート:Indent

これは、テンプレート:Math かつ テンプレート:Indent ならば テンプレート:Indent と書くこともできる。

初期値の定理・最終値の定理

ラプラス変換の原関数の初期値(テンプレート:Math での値)や最終値(テンプレート:Math における極限値)を表す初期値の定理 テンプレート:En および最終値の定理 テンプレート:En と呼ばれる公式が以下のような式によって与えられる。

初期値の定理
テンプレート:Mvar の関数 テンプレート:Mathテンプレート:Math で連続ならば
f(0)=limt0f(t)=limssF(s)
が成り立つ。特に、テンプレート:Mvar が微分可能なときは部分積分により容易に証明できる。
最終値の定理
テンプレート:Mvar の関数 テンプレート:Mathテンプレート:Math で収束するなら
f()=limtf(t)=lims0sF(s)sΔ0
が成り立つ。ただし、テンプレート:Mathテンプレート:Math を含む角領域である。

性質一覧表

  • 表中の凡例
u(t) : ヘビサイド関数
(f*g)(t) : テンプレート:Mvarテンプレート:Mvar畳み込み
f(t) : テンプレート:Mathテンプレート:Math 階微分
f(n)(t) : テンプレート:Mathテンプレート:Mvar 階微分
片側ラプラス変換の性質(その 1)
性質 原関数
f(t)=1{F(s)}
('テンプレート:Mvar' 領域 / 時間領域)
像関数
F(s)={f(t)}
('テンプレート:Mvar' 領域 / 周波数領域)
備考
線形性 af(t)+bg(t) aF(s)+bG(s)
相似性 f(at) 1aF(sa) ただし、テンプレート:Math
移動 f(ta)u(ta) easF(s)
f(t+λ) eλs{F(s)0λestf(t)dt} 移動第 2 則
ただし、テンプレート:Math
テンプレート:Math微分 f(t) sF(s)f(0) ただし、テンプレート:Mvarテンプレート:Math微分可能とする。
テンプレート:Math 階微分 f(t) s2F(s)sf(0)f(0) ただし、テンプレート:Mvarテンプレート:Math 階微分可能とする。
テンプレート:Mvar 階微分 f(n)(t) snF(s)sn1f(0)f(n1)(0) ただし、テンプレート:Mvarテンプレート:Mvar 階微分可能とする。
積分 0tf(τ)dτ=(u*f)(t) 1sF(s)
1(n1)!0t(tq)n1f(q)dq 1snF(s) ただし、テンプレート:Math
0t0τn10τ1f(τ)dτdτ1dτn1 1snF(s)
畳み込み (f*g)(t)=0tf(u)g(tu)du F(s)G(s)
周期関数 f(t)=f(t+T) 11eTs0Testf(t)dt テンプレート:Math は周期 テンプレート:Mvar の周期関数。
片側ラプラス変換の性質(その 2)
性質 像関数
F(s)={f(t)}
('テンプレート:Mvar' 領域 / 周波数領域)
原関数
f(t)=1{F(s)}
('テンプレート:Mvar' 領域 / 時間領域)
備考
移動 F(sa) eatf(t)
テンプレート:Math 階微分 F(s) tf(t) ただし、テンプレート:Mvarテンプレート:Math 階微分可能とする。
テンプレート:Math 階微分 F(s) t2f(t) ただし、テンプレート:Mvarテンプレート:Math 階微分可能とする。
テンプレート:Mvar 階微分 F(n)(s) (t)nf(t) ただし、テンプレート:Mvarテンプレート:Mvar 階微分可能とする。
積分 sF(σ)dσ f(t)t
sσn1σ1F(σ)dσdσ1dσn1 f(t)tn
畳み込み (F*G)(s)=12πicic+iF(σ)G(sσ)dσ f(t)g(t)

変換表

変換表 原関数
f(t)=1{F(s)}
'テンプレート:Mvar' 領域 / 時間領域
像関数
F(s)={f(t)}
'テンプレート:Mvar' 領域 / 周波数領域
収束域
単位インパルス δ(t) 1 alls
単位ステップ関数 u(t) 1s Re{s}>0
ランプ関数 tu(t) 1s2 Re{s}>0
テンプレート:Mvar
テンプレート:Mvar は整数)
tnn!u(t) 1sn+1 Re{s}>0
(n>1)
テンプレート:Mvar
テンプレート:Mvar は複素数)
tqΓ(q+1)u(t) 1sq+1 Re{s}>0
(Re{q}>1)
テンプレート:Mvar 乗根 tnu(t)=t1/nu(t) 1s1+1/nΓ(1+1n) Re{s}>0
指数減衰 eαtu(t) 1s+α Re{s}>α
テンプレート:Mvar 乗の指数減衰 tnn!eαtu(t) 1(s+α)n+1 Re{s}>α
理想遅延 δ(tτ) eτs
遅延付き単位ステップ関数 u(tτ) 1seτs Re{s}>0
遅延付き テンプレート:Mvar 乗の指数減衰 (tτ)nn!eα(tτ)u(tτ) 1(s+α)n+1eτs Re{s}>α
指数関数的接近 (1eαt)u(t) αs(s+α) Re{s}>0
正弦関数 sin(ωt)u(t) ωs2+ω2 Re{s}>0
余弦関数 cos(ωt)u(t) ss2+ω2 Re{s}>0
双曲線正弦関数
(ハイパボリックサイン)
sinh(αt)u(t) αs2α2 Re{s}>|α|
双曲線余弦関数
(ハイパボリックコサイン)
cosh(αt)u(t) ss2α2 Re{s}>|α|
正弦波の指数減衰 eαtsin(ωt)u(t) ω(s+α)2+ω2 Re{s}>α
余弦波の指数減衰 eαtcos(ωt)u(t) s+α(s+α)2+ω2 Re{s}>α
自然対数 ln(tt0)u(t) 1s[ln(t0s)+γ] Re{s}>0
第 1 種ベッセル関数 Jn(ωt)u(t) ωn(s+s2+ω2)ns2+ω2 Re{s}>0
(n>1)
第 1 種変形ベッセル関数 In(ωt)u(t) ωn(s+s2ω2)ns2ω2 Re{s}>|ω|
第 2 種ベッセル関数
(次数が テンプレート:Math の場合)
Y0(αt)u(t) 2sinh1(s/α)πs2+α2 Re{s}>0
第 2 種変形ベッセル関数
(次数が テンプレート:Math の場合)
K0(αt)u(t) cos1(s/α)α2s2  Re{s}<|α| 
誤差関数 erf(t)u(t) es2/4(1erf(s/2))s Re{s}>0
凡例

テンプレート:Col-begin テンプレート:Col-break

テンプレート:Col-break

テンプレート:Col-end

関連文献

関連項目

テンプレート:Div col

テンプレート:Div col end

外部リンク

テンプレート:制御理論 テンプレート:Normdaten