ラプラス変換

数学 において、ラプラス変換は、ピエール=シモン・ラプラス/ l ə ˈ p l ɑː s /)にちなんで名付けられ、変数関数(通常は時間領域)を複素変数関数(複素数値周波数領域、 s領域またはs平面とも呼ばれる)に変換する積分変換です。これらの関数は、時間領域関数には小文字の記号、周波数領域関数には対応する大文字の記号を使用して表記されることがよくあります(例: および⁠ t{\displaystyle t}s{\displaystyle s}×t{\displaystyle x(t)}Xs{\displaystyle X(s)}

この変換は、時間領域における微分積分を、ラプラス領域におけるはるかに簡単な乗算除算に変換するのに有用である(対数が乗算と除算を加算と減算に簡略化するのと同様に)。この変換は科学技術において多くの応用があり、主に常微分方程式積分方程式を代数多項式方程式に簡略化し、畳み込みを乗算に簡略化することで、線形微分方程式[ 1 ]動的システムを解くためのツールとして用いられている。[ 2 ] [ 3 ]

例えば、ラプラス変換によって、単振動子の方程式(フックの法則)は、初期条件を含む代数方程式に変換され、未知関数について解くことができます。解いた後、逆ラプラス変換を用いて元の領域に変換することができます。これは、以下に 示すような表を参照することで容易に行うことができます。 ×t+×t0{\displaystyle x''(t)+kx(t)=0}s2Xss×0×0+Xs0{\displaystyle s^{2}X(s)-sx(0)-x'(0)+kX(s)=0,}×0{\displaystyle x(0)}×0{\displaystyle x'(0)}Xs{\displaystyle X(s).}

ラプラス変換は(適切な関数⁠ ⁠f{\displaystyle f}に対して)積分 によって定義されます。 ここで、sは複素数です。 L{f}s0ftestdt{\displaystyle {\mathcal {L}}\{f\}(s)=\int _{0}^{\infty }f(t)e^{-st}\,dt,}

ラプラス変換は他の多くの変換と関連しています。メリン変換と本質的に同じであり、フーリエ変換と密接に関連しています。フーリエ変換とは異なり、関数のラプラス変換は解析関数であることが多いです。つまり、収束するべき級数を持ち、その係数は元の関数のモーメントを表します。さらに、複素解析の手法、特に線積分は計算を簡素化するために使用できます。

歴史

ラプラス侯爵ピエール=シモン

ラプラス変換は、確率論の研究で同様の変換を使用した数学者天文学者のピエール・シモン・ラプラス侯爵にちなんで名付けられました。[ 4 ] [ 5 ]ラプラスは生成関数の使用について広範囲に渡って著作を残し (1814年)、その結果としてラプラス変換の積分形式が自然に発展しました。[ 6 ]

ラプラスの母関数の使用は、現在Z変換として知られているものと似ており、ニールス・ヘンリク・アーベルが議論した連続変数のケースにはほとんど注意を払っていなかった。[ 7 ]

1744年より、レオンハルト・オイラーは微分方程式の解として形式の積分を研究し 、特にガンマ関数を導入した。[ 8 ]ジョゼフ=ルイ・ラグランジュはオイラーの崇拝者であり、確率密度関数の積分に関する研究において、 ラプラス変換に似た形式の表現を研究した 。 [ 9 ] [ 10 ]zX×e1つの×d× そして zX××d×{\displaystyle z=\int X(x)e^{ax}\,dx\quad {\text{ と }}\quad z=\int X(x)x^{A}\,dx}X×e1つの×1つの×d×{\displaystyle \int X(x)e^{-ax}a^{x}\,dx,}

この種の積分は、1782年に初めてラプラスの注目を集めたようで、彼はオイラーの精神に倣い、積分自体を方程式の解として用いていた。[ 11 ]しかし、1785年にラプラスは、単に積分の形で解を探すのではなく、後に普及する意味で変換を適用し始めたことで、決定的な一歩を踏み出した。彼は、 メリン変換 に似た形の積分を用いて差分方程式全体を変換し、変換された方程式の解を探した。その後、彼はラプラス変換を同じように適用してその特性のいくつかを導き出し、その潜在的な力を理解し始めた。[ 12 ]×sφ×d×{\displaystyle \int x^{s}\varphi (x)\,dx,}

ラプラスはまた、ジョセフ・フーリエフーリエ級数法による拡散方程式の解法は、その解が周期的であるため、限られた空間領域にしか適用できないことを認識していました。1809年、ラプラスはこの変換を用いて、空間に無限に拡散する解を求めました。[ 13 ] 1821年、コーシーはラプラス変換の演算法を開発し、これは現在基礎工学で使用されているのとほぼ同じ方法で線型微分方程式の研究に使用できました。この方法は、世紀の変わり目頃にオリバー・ヘヴィサイドによって普及し、おそらく再発見されました。 [ 14 ]

ベルンハルト・リーマンは1859年の論文『与えられた大きさより小さい素数の個数について』でラプラス変換を用い、その中で反転定理も導出した。リーマンはラプラス変換を用いてリーマンゼータ関数の関数方程式を導出し、彼の手法は今でもヤコビ・シータ関数のモジュラー変換則をこの関数方程式に関連付けるのに用いられている。ヤコビ・シータ関数のモジュラー変換則はポアソン和によって簡単に証明できる。[ 15 ]

20世紀初頭、ヤルマー・メリンはカール・ワイエルシュトラス流の解析学派においてラプラス変換を厳密に研究し、それを微分方程式特殊関数の研究に応用した最初の研究者の一人であった。 [ 16 ] 同じ頃、ヘヴィサイドは演算計算に取り組んでいた。 トーマス・ヨハンネス・スティルチェスは、モーメントに関する研究に関連してラプラス変換の一般化を検討した。この時期の他の貢献者には、マティアス・レルヒ[ 17 ]オリバー・ヘヴィサイドトーマス・ブロムウィッチなどがいた。[ 18 ]

1929年、ヴァネヴァー・ブッシュノーバート・ウィーナーは、フーリエ変換と演算法の両方を適用した電気回路の工学解析のテキストとして演算回路解析を出版し、その中に現代のラプラス変換表の最初の前身の1つを含めました。 1934年、レイモンド・ペイリーノーバート・ウィーナーは、現在ラプラス変換と呼ばれているものについての重要な著作、複素領域におけるフーリエ変換を出版しました(下記参照)。 また1930年代には、ラプラス変換はGHハーディジョン・エデンサー・リトルウッドによるタウバー定理の研究に役立ち、この応用は後にウィダー(1941年)によって詳しく説明され、ウィダーは反転の新しい方法など理論の別の側面を発展させました。 エドワード・チャールズ・ティッチマーシュは、影響力のあるフーリエ積分理論の入門書(1937年)を執筆しました。

この変換が現在(主に工学分野で)広く利用されるようになったのは、第二次世界大戦中およびその直後であり、[ 19 ]以前のヘヴィサイドの演算法に取って代わったためである。ラプラス変換の利点はグスタフ・デーチ[ 20 ]によって強調されており、ラプラス変換という名称は彼に由来すると思われる。

正式な定義

est{\displaystyle \Re (e^{-st})}s領域における様々な複素周波数は次のように表すことができます。軸には純粋な余弦が含まれます。正の軸には減衰余弦が含まれます。負の軸には指数関数的に増加する余弦が含まれます。sσ+ω{\displaystyle (s=\sigma +i\omega ),}eσtコスωt{\displaystyle e^{-\sigma t}\cos(\omega t).}σ0{\displaystyle \sigma =0}σ{\displaystyle \sigma }σ{\displaystyle \sigma }

関数f ( t )のラプラス変換は、すべての実数t ≥ 0に対して定義され、関数F ( s )であり、これは次のように定義される片側変換である。

Fs0ftestdt{\displaystyle F(s)=\int _{0}^{\infty }f(t)e^{-st}\,dt,}   (式1

ここで、sは 実数σωを持つ複素周波数領域パラメータ です。 sσ+ω{\displaystyle s=\sigma +i\omega }

ラプラス変換の別の表記法はL{f}{\displaystyle {\mathcal {L}}\{f\}}Fの代わりに を用いる。[ 3 ]関数記法では となる。これは、特に工学分野では、ダミー変数が関数には現れないことを理解した上で、 と表記されることが多い。 FsL{f}s{\displaystyle F(s)={\mathcal {L}}\{f\}(s)}FsL{ft}{\displaystyle F(s)={\mathcal {L}}\{f(t)\}}t{\displaystyle t}Fs{\displaystyle F(s)}

積分の意味は、対象となる関数の種類によって異なります。積分が存在するための必要条件は、f が[0, ∞)上で局所的に積分可能であることです。無限大で減少する、または指数型( )である局所的に積分可能な関数の場合、積分は(適切な)ルベーグ積分であると理解できます。ただし、多くの用途では、これを条件付き収束する不適切積分と見なす必要があります。さらに一般的には、積分は弱い意味で理解することができ、これについては以下で扱います。 |ft|eB|t|{\displaystyle \vert f(t)\vert \leq Ae^{B\vert t\vert }}

有限ボレル測度μのラプラス変換はルベーグ積分によって定義できる[ 21 ]L{μ}(s)=[0,)estdμ(t).{\displaystyle {\mathcal {L}}\{\mu \}(s)=\int _{[0,\infty )}e^{-st}\,d\mu (t).}

重要な特別なケースとして、μが確率測度、例えばディラックのデルタ関数である場合が挙げられます。演算計算では、測度のラプラス変換は、その測度が確率密度関数fから得られたものとして扱われることがよくあります。その場合、混乱を避けるため、 と書くことがよくあります。 ここで、 0 の下限はの略記です。 L{f}(s)=0f(t)estdt,{\displaystyle {\mathcal {L}}\{f\}(s)=\int _{0^{-}}^{\infty }f(t)e^{-st}\,dt,}limε0+ε.{\displaystyle \lim _{\varepsilon \to 0^{+}}\int _{-\varepsilon }^{\infty }.}

この極限は、0に位置する任意の質点がラプラス変換によって完全に捉えられることを強調しています。ルベーグ積分においてはこのような極限を取る必要はありませんが、ラプラス・スティルチェス変換と関連付けるとより自然に現れます。

両側ラプラス変換

「ラプラス変換」と限定なしに言う場合、通常は片側変換を指します。ラプラス変換は、積分の極限を実軸全体に拡張することで、両側ラプラス変換、または両側ラプラス変換と定義することもできます。その場合、一般的な片側変換は単に両側変換の特殊なケースとなり、変換される関数の定義にはヘビサイドの階段関数が乗算されることが含まれます。

双対ラプラス変換F ( s )は次のように定義されます。

F(s)=estf(t)dt.{\displaystyle F(s)=\int _{-\infty }^{\infty }e^{-st}f(t)\,dt.}   (式2

両側ラプラス変換の代替表記は、 Fではなく⁠ ⁠B{f}{\displaystyle {\mathcal {B}}\{f\}}です。

逆ラプラス変換

二つの可積分関数が同一のラプラス変換を持つのは、ルベーグ測度零点の集合において両者が異なる場合に限られます。これは、変換の値域において逆変換が存在することを意味します。実際、ラプラス変換は可積分関数だけでなく、他の多くの関数空間においても、ある関数空間から別の関数空間への一対一写像です。ただし、その値域を容易に特徴付けることは通常できません。

これが成り立つ典型的な関数空間としては、有界連続関数の空間、L (0, ∞)、あるいはより一般的には(0, ∞)上の緩和超関数の空間が挙げられる。ラプラス変換もまた定義され、緩和超関数の適切な空間に対して単射である。

これらの場合、ラプラス変換の像は収束領域における解析関数の空間に存在します。逆ラプラス変換は、以下の複素積分で与えられ、これは様々な名前(ブロムウィッチ積分フーリエ・メリン積分、メリンの逆公式など) で知られています。

f(t)=L1{F}(t)=12πilimTγiTγ+iTestF(s)ds,{\displaystyle f(t)={\mathcal {L}}^{-1}\{F\}(t)={\frac {1}{2\pi i}}\lim _{T\to \infty }\int _{\gamma -iT}^{\gamma +iT}e^{st}F(s)\,ds,}   (式3

ここでγは実数であり、積分の収束経路はF ( s )の収束領域内にあります。ほとんどの応用では、収束経路は閉じているため、留数定理を使用することができます。逆ラプラス変換の別の公式は、ポストの逆変換公式によって与えられます。ここでの極限は、弱*位相で解釈されます。

実際には、ラプラス変換をテーブルから取得した関数の既知の変換に分解し、検査によって逆変換を構築する方が通常は便利です。

確率論

純粋確率および応用確率において、ラプラス変換は期待値として定義されます。X確率密度関数fを持つ確率変数である場合、 fのラプラス変換は期待値で与えられます 。 ここで、は確率変数期待値です。 L{f}(s)=E[esX],{\displaystyle {\mathcal {L}}\{f\}(s)=\operatorname {E} \left[e^{-sX}\right],}E[r]{\displaystyle \operatorname {E} [r]}r{\displaystyle r}

慣例により、これは確率変数X自体のラプラス変換と呼ばれます。ここで、sを-tに置き換えると Xモーメント生成関数が得られます。ラプラス変換は、マルコフ連鎖などの確率過程初通過時間再生理論など、確率論全般に応用されています。

特に有用なのは、ラプラス変換によって連続確率変数Xの累積分布関数を次のように復元できることである。 [ 22 ]FX(x)=L1{1sE[esX]}(x)=L1{1sL{f}(s)}(x).{\displaystyle F_{X}(x)={\mathcal {L}}^{-1}\left\{{\frac {1}{s}}\operatorname {E} \left[e^{-sX}\right]\right\}(x)={\mathcal {L}}^{-1}\left\{{\frac {1}{s}}{\mathcal {L}}\{f\}(s)\right\}(x).}

代数的構成

ラプラス変換は、正半直線上の関数の畳み込み環に分数体構成を適用することで、純粋に代数的に定義することもできる。結果として得られる抽象作用素の空間はラプラス空間と全く同じであるが、この構成では順方向変換と逆方向変換を明示的に定義する必要がない(収束性の証明に伴う困難を回避するため)。[ 23 ]

収束領域

fが局所的に積分可能な関数(またはより一般的には局所的に有界な変化のボレル測度)である場合、 fのラプラス変換Fsは極限が 存在する限り収束します。 limR0Rf(t)estdt{\displaystyle \lim _{R\to \infty }\int _{0}^{R}f(t)e^{-st}\,dt}

ラプラス変換は、積分が真のルベーグ積分として存在する場合、 絶対収束します 。ラプラス変換は通常、条件付き収束と理解されており、前者の意味では収束しますが、後者の意味で収束しません。 0|f(t)est|dt{\displaystyle \int _{0}^{\infty }\left|f(t)e^{-st}\right|\,dt}

F ( s )が絶対収束する値の集合は、Re( s ) > aまたはRe( s ) ≥ aの形式のいずれかであり、ここでa は、 −∞ ≤ a ≤ ∞を満たす拡張された実定数です(優勢収束定理の結果)。定数aは絶対収束の横座標として知られ、f ( t )の増加動作に依存します。[ 24 ]同様に、両側変換は、 a < Re( s ) < bの形式の帯で絶対収束し、おそらくRe( s ) = aまたはRe( s ) = bの線を含みます。[ 25 ]ラプラス変換が絶対収束するsの値のサブセットは、絶対収束領域、または絶対収束の定義域と呼ばれます。両側の場合、これは絶対収束帯と呼ばれることもあります。ラプラス変換は絶対収束の領域では解析的です。これは、フビニの定理モレラの定理の結果です。

同様に、 F ( s ) が(条件付きまたは絶対的に)収束する値の集合は、条件付き収束領域、または単に収束領域(ROC)と呼ばれます。ラプラス変換がs = s 0で(条件付きで)収束する場合、 Re( s ) > Re( s 0 )を満たすすべてのsに対して自動的に収束します。したがって、収束領域はRe( s ) > aの形をした半平面であり、境界線Re( s ) = a 上のいくつかの点を含む可能性があります。

収束領域Re( s ) > Re( s0 )では、 fのラプラス変換は部分積分によって次のように 表されるF(s)=(ss0)0e(ss0)tβ(t)dt,β(u)=0ues0tf(t)dt.{\displaystyle F(s)=(s-s_{0})\int _{0}^{\infty }e^{-(s-s_{0})t}\beta (t)\,dt,\quad \beta (u)=\int _{0}^{u}e^{-s_{0}t}f(t)\,dt.}

つまり、収束領域において、F ( s ) は、他の関数の絶対収束ラプラス変換として効果的に表現できる。特に、これは解析的である。最も一般的な形では、ラプラス変換は、ある に対してで定義され、絶対値が多項式によってそこで制限される正則関数と、 で支えられ、あるに対してを乗じると緩和分布になる実数直線上の分布との間に一対一の対応を与える。[ 26 ]σR{\displaystyle \sigma \in \mathbb {R} }{sC | Re(s)>σ}{\displaystyle \{s\in \mathbb {C} \ |\ \mathrm {Re} (s)>\sigma \}}[0,){\displaystyle [0,\infty )}eσt{\displaystyle e^{-\sigma t}}σ{\displaystyle \sigma }

fの減衰特性と収束領域内のラプラス変換の特性 との関係に関しては、いくつかのPaley-Wiener の定理があります。

工学応用において、線形時不変(LTI)システムに対応する関数は、すべての有界入力が有界出力を生成する場合、安定である。これは、インパルス応答関数のラプラス変換がRe( s ) ≥ 0の領域で絶対収束することと等価である。結果として、インパルス応答関数のラプラス変換の極が負の実部を持つ限り、LTIシステムは安定である。

この ROC は、システムの因果関係と安定性を知るために使用されます。

性質と定理

ラプラス変換の重要な性質は、時間領域における微分積分を、ラプラス領域におけるsによる乗算と除算に変換することです。したがって、ラプラス変数s は、ラプラス領域における演算子変数とも呼ばれ、微分演算子、または(s −1の場合)積分演算子と呼ばれます。

関数f ( t )g ( t )、およびそれぞれのラプラス変換F ( s )G ( s )が与えられれば、 f(t)=L1{F(s)},g(t)=L1{G(s)},{\displaystyle {\begin{aligned}f(t)&={\mathcal {L}}^{-1}\{F(s)\},\\g(t)&={\mathcal {L}}^{-1}\{G(s)\},\end{aligned}}}

次の表は片側ラプラス変換の性質の一覧である。[ 27 ]

片側ラプラス変換の性質
財産 時間領域 sドメイン コメント
直線性af(t)+bg(t) {\displaystyle af(t)+bg(t)\ }aF(s)+bG(s) {\displaystyle aF(s)+bG(s)\ }積分の基本規則を使用して証明できます。
周波数領域微分 tf(t) {\displaystyle tf(t)\ }F(s) {\displaystyle -F'(s)\ }F ′はsに関するFの 1 次導関数です。
周波数領域一般微分 tnf(t) {\displaystyle t^{n}f(t)\ }(1)nF(n)(s) {\displaystyle (-1)^{n}F^{(n)}(s)\ }より一般的な形式は、F ( s )のn次導関数です。
デリバティブf(t) {\displaystyle f'(t)\ }sF(s)f(0) {\displaystyle sF(s)-f(0^{-})\ }fは微分可能関数であり、その導関数は指数型であると仮定する。これは部分積分によって得られる。
2階微分 f(t) {\displaystyle f''(t)\ }s2F(s)sf(0)f(0) {\textstyle s^{2}F(s)-sf(0^{-})-f'(0^{-})\ }fは2回微分可能であり、2階微分は指数型であると仮定する。次に、 f ′( t )に微分法を適用して計算する。
一般微分 f(n)(t) {\displaystyle f^{(n)}(t)\ }snF(s)k=1nsnkf(k1)(0) {\displaystyle s^{n}F(s)-\sum _{k=1}^{n}s^{n-k}f^{(k-1)}(0^{-})\ }fはn回微分可能であり、 n次導関数は指数型であると仮定する。数学的帰納法により導かれる。
周波数領域積分1tf(t) {\displaystyle {\frac {1}{t}}f(t)\ }sF(σ)dσ {\displaystyle \int _{s}^{\infty }F(\sigma )\,d\sigma \ }これは周波数微分と条件付き収束の性質を利用して推論されます。
時間領域積分 0tf(τ)dτ=(uf)(t){\displaystyle \int _{0}^{t}f(\tau )\,d\tau =(u*f)(t)}1sF(s){\displaystyle {1 \over s}F(s)}u ( t )はヘヴィサイドステップ関数であり、 ( uf )( t )はu ( t )f ( t )畳み込みです。
周波数シフト eatf(t){\displaystyle e^{at}f(t)}F(sa){\displaystyle F(s-a)}
タイムシフト f(ta)u(ta){\displaystyle f(t-a)u(t-a)}

f(t)u(ta) {\displaystyle f(t)u(t-a)\ }

easF(s) {\displaystyle e^{-as}F(s)\ }

easL{f(t+a)}{\displaystyle e^{-as}{\mathcal {L}}\{f(t+a)\}}

a > 0 u ( t )はヘヴィサイドステップ関数である。
時間スケーリング f(at){\displaystyle f(at)}1aF(sa){\displaystyle {\frac {1}{a}}F\left({s \over a}\right)}a > 0
乗算f(t)g(t){\displaystyle f(t)g(t)}12πilimTciTc+iTF(σ)G(sσ)dσ {\displaystyle {\frac {1}{2\pi i}}\lim _{T\to \infty }\int _{c-iT}^{c+iT}F(\sigma )G(s-\sigma )\,d\sigma \ }積分はFの収束領域内に完全に含まれる垂直線Re( σ )= cに沿って行われる。[ 28 ]
畳み込み(fg)(t)=0tf(τ)g(tτ)dτ{\displaystyle (f*g)(t)=\int _{0}^{t}f(\tau )g(t-\tau )\,d\tau }F(s)G(s) {\displaystyle F(s)\cdot G(s)\ }
循環畳み込み(fg)(t)=0Tf(τ)g(tτ)dτ{\displaystyle (f*g)(t)=\int _{0}^{T}f(\tau )g(t-\tau )\,d\tau }F(s)G(s) {\displaystyle F(s)\cdot G(s)\ }周期Tの周期関数の場合。
複雑な活用f(t){\displaystyle f^{*}(t)}F(s){\displaystyle F^{*}(s^{*})}
周期関数f(t){\displaystyle f(t)}11eTs0Testf(t)dt{\displaystyle {1 \over 1-e^{-Ts}}\int _{0}^{T}e^{-st}f(t)\,dt}f ( t )は周期Tの周期関数であり、すべてのt ≥ 0に対してf ( t ) = f ( t + T )が成り立ちます。これは時間シフト特性と等比級数の結果です。
周期的な合計fP(t)=n=0f(tTn){\displaystyle f_{P}(t)=\sum _{n=0}^{\infty }f(t-Tn)}

fP(t)=n=0(1)nf(tTn){\displaystyle f_{P}(t)=\sum _{n=0}^{\infty }(-1)^{n}f(t-Tn)}

FP(s)=11eTsF(s){\displaystyle F_{P}(s)={\frac {1}{1-e^{-Ts}}}F(s)}

FP(s)=11+eTsF(s){\displaystyle F_{P}(s)={\frac {1}{1+e^{-Ts}}}F(s)}

初期値定理
f(0+)=limssF(s).{\displaystyle f(0^{+})=\lim _{s\to \infty }{sF(s)}.}
最終値定理
⁠ ⁠f()=lims0sF(s){\displaystyle f(\infty )=\lim _{s\to 0}{sF(s)}} 、のすべてのが左半平面内にある場合。sF(s){\displaystyle sF(s)}
終値定理は、部分分数分解(またはその他の難しい代数計算)を行うことなく長期的な挙動を導くため有用である。F ( s ) が右手平面に極を持つ場合、または虚軸上に極を持つ場合(例えば、または)、この式の挙動は未定義である。f(t)=et{\displaystyle f(t)=e^{t}}f(t)=sin(t){\displaystyle f(t)=\sin(t)}

べき級数との関係

ラプラス変換は、連続的なべき級数とみなすことができます。[ 29 ] a ( n )が正の整数nの離散関数である場合 、 a ( n )に関連付けられたべき級数は、xが実変数 である級数です ( Z変換を参照)。n上の総和をt上の積分に 置き換えると、べき級数の連続バージョンは次のようになります。 ここで、離散関数a ( n )は連続関数f ( t )に置き換えられます。 n=0a(n)xn{\displaystyle \sum _{n=0}^{\infty }a(n)x^{n}}0f(t)xtdt{\displaystyle \int _{0}^{\infty }f(t)x^{t}\,dt}

べき乗の基数をxからeに変更すると 、0f(t)(elnx)tdt{\displaystyle \int _{0}^{\infty }f(t)\left(e^{\ln {x}}\right)^{t}\,dt}

例えば、すべての有界関数fに対して収束するには、 ln x < 0を満たす必要があります。 s = ln xと代入すると、ラプラス変換が得られます。 0f(t)estdt{\displaystyle \int _{0}^{\infty }f(t)e^{-st}\,dt}

言い換えれば、ラプラス変換は、離散パラメータn が連続パラメータtに置き換えられ、x がe sに置き換えられた、累乗級数の連続アナログです。

べき級数と同様に、⁠ ⁠a(n)=O(ρn){\displaystyle a(n)=O(\rho ^{-n})}ならば、べき級数は⁠ ⁠における解析関数に収束し、 |x|<ρ{\displaystyle \vert x\vert <\rho }⁠ ⁠f(t)=O(eσt){\displaystyle f(t)=O(e^{-\sigma t})}ならば、ラプラス変換は⁠ ⁠(s)>σ{\displaystyle \Re (s)>\sigma }における解析関数に収束します。[ 30 ]

瞬間との関係

これらの量 は関数fのモーメント である。 fの最初のn個のモーメントが絶対収束する場合、積分 の下での繰り返し微分により、 が得られる。 これは確率論において特に重要であり、確率変数X のモーメントは期待値 で与えられるしたがって、関係式は μn=0tnf(t)dt{\displaystyle \mu _{n}=\int _{0}^{\infty }t^{n}f(t)\,dt}(1)n(Lf)(n)(0)=μn.{\displaystyle (-1)^{n}({\mathcal {L}}f)^{(n)}(0)=\mu _{n}.}μn=E[Xn]{\displaystyle \mu _{n}=\operatorname {E} [X^{n}]}μn=(1)ndndsnE[esX](0).{\displaystyle \mu _{n}=(-1)^{n}{\frac {d^{n}}{ds^{n}}}\operatorname {E} \left[e^{-sX}\right](0).}

関数の導関数の変換

ラプラス変換の微分性を利用して関数の導関数を求めると便利な場合が多い。これはラプラス変換の基本式から次のように導出できる。 そして、両側変換の場合は 、 L{f(t)}=0estf(t)dt=[f(t)ests]00estsf(t)dt(by parts)=[f(0)s]+1sL{f(t)},{\displaystyle {\begin{aligned}{\mathcal {L}}\left\{f(t)\right\}&=\int _{0^{-}}^{\infty }e^{-st}f(t)\,dt\\[6pt]&=\left[{\frac {f(t)e^{-st}}{-s}}\right]_{0^{-}}^{\infty }-\int _{0^{-}}^{\infty }{\frac {e^{-st}}{-s}}f'(t)\,dt\quad {\text{(by parts)}}\\[6pt]&=\left[-{\frac {f(0^{-})}{-s}}\right]+{\frac {1}{s}}{\mathcal {L}}\left\{f'(t)\right\},\end{aligned}}}L{f(t)}=sL{f(t)}f(0),{\displaystyle {\mathcal {L}}\{f'(t)\}=s\cdot {\mathcal {L}}\{f(t)\}-f(0^{-}),}L{f(t)}=sestf(t)dt=sL{f(t)}.{\displaystyle {\mathcal {L}}\{f'(t)\}=s\int _{-\infty }^{\infty }e^{-st}f(t)\,dt=s\cdot {\mathcal {L}}\{f(t)\}.}

一般的な結果は 、fのn次導関数を表し、帰納的議論によって確立できます。 L{f(n)(t)}=snL{f(t)}sn1f(0)f(n1)(0),{\displaystyle {\mathcal {L}}\left\{f^{(n)}(t)\right\}=s^{n}\cdot {\mathcal {L}}\{f(t)\}-s^{n-1}f(0^{-})-\cdots -f^{(n-1)}(0^{-}),}f(n){\displaystyle f^{(n)}}

正の実軸上の積分を評価する

ラプラス変換の有用な性質は、 ⁠ の右近傍における の挙動と ⁠ の左近傍におけるの減衰率に関する適切な仮定の下で、次のようになることです。上記の式は部分積分を変形したもので、演算子 と を と ⁠ ⁠ に置き換えたものです 。同等定式化を証明しましょう。 0f(x)g(x)dx=0(Lf)(s)(L1g)(s)ds{\displaystyle \int _{0}^{\infty }f(x)g(x)\,dx=\int _{0}^{\infty }({\mathcal {L}}f)(s)\cdot ({\mathcal {L}}^{-1}g)(s)\,ds}f,g{\displaystyle f,g}0{\displaystyle 0}f,g{\displaystyle f,g}{\displaystyle \infty }ddx{\displaystyle {\frac {d}{dx}}}dx{\displaystyle \int \,dx}L{\displaystyle {\mathcal {L}}}L1{\displaystyle {\mathcal {L}}^{-1}}0(Lf)(x)g(x)dx=0f(s)(Lg)(s)ds.{\displaystyle \int _{0}^{\infty }({\mathcal {L}}f)(x)g(x)\,dx=\int _{0}^{\infty }f(s)({\mathcal {L}}g)(s)\,ds.}

を代入すると、左側は次のように変わります。 ただし、Fubini の定理が成り立つと仮定すると、積分の順序を逆にすると、必要な右側が得られます。 (Lf)(x)=0f(s)esxds{\displaystyle ({\mathcal {L}}f)(x)=\int _{0}^{\infty }f(s)e^{-sx}\,ds}00f(s)g(x)esxdsdx,{\displaystyle \int _{0}^{\infty }\int _{0}^{\infty }f(s)g(x)e^{-sx}\,ds\,dx,}

この方法は、実微積分学の初等的な手法では計算が困難な積分を計算するのに使用できます。例えば、 0sinxxdx=0L(1)(x)sinxdx=01L(sin)(x)dx=0dxx2+1=π2.{\displaystyle \int _{0}^{\infty }{\frac {\sin x}{x}}dx=\int _{0}^{\infty }{\mathcal {L}}(1)(x)\sin xdx=\int _{0}^{\infty }1\cdot {\mathcal {L}}(\sin )(x)dx=\int _{0}^{\infty }{\frac {dx}{x^{2}+1}}={\frac {\pi }{2}}.}

他の変換との関係

ラプラス・スティルチェス変換

関数g  : ℝ → ℝの (一方的) ラプラス・スティールチェス変換は、ルベーグ・スティールチェス積分によって定義されます。{Lg}(s)=0estdg(t) .{\displaystyle \{{\mathcal {L}}^{*}g\}(s)=\int _{0}^{\infty }e^{-st}\,d\,g(t)~.}

関数gは有界変化であると仮定する。gf不定積分であるとすると、 g(x)=0xf(t)dt{\displaystyle g(x)=\int _{0}^{x}f(t)\,d\,t}

とすると、 gのラプラス・スティルチェス変換とfのラプラス変換は一致する。一般に、ラプラス・スティルチェス変換はgに関連付けられたスティルチェス測度のラプラス変換である。したがって、実際には、2つの変換の唯一の違いは、ラプラス変換が測度の密度関数に作用すると考えられるのに対し、ラプラス・スティルチェス変換はその累積分布関数に作用すると考えられるという点である。[ 31 ]

フーリエ変換

を⁠ で支えられた複素数値ルベーグ可積分関数とし、をそのラプラス変換とする。すると、収束領域内では、 関数のフーリエ変換が成り立つ。[ 32 ]f{\displaystyle f}[0,){\displaystyle [0,\infty )}F(s)=Lf(s){\displaystyle F(s)={\mathcal {L}}f(s)}F(σ+iτ)=0f(t)eσteiτtdt,{\displaystyle F(\sigma +i\tau )=\int _{0}^{\infty }f(t)e^{-\sigma t}e^{-i\tau t}\,dt,}f(t)eσt{\displaystyle f(t)e^{-\sigma t}}

実際、フーリエ変換は、(特定の条件下では)双対ラプラス変換の特殊なケースです。主な違いは、関数のフーリエ変換が変数(周波数⁠ ⁠τ{\displaystyle \tau })の複素関数であるのに対し、関数のラプラス変換は複素変数(減衰係数と周波数⁠ )の複素関数であることです。ラプラス変換は通常、 t ≥ 0tの関数の変換に制限されます。この制限の結果として、関数のラプラス変換は変数sの正則関数になります。フーリエ変換とは異なり、分布のラプラス変換は一般に、行儀の良い関数です。複素変数の手法を使用して、ラプラス変換を直接研究することもできます。正則関数として、ラプラス変換はべき級数表現を持ちます。このべき級数は、関数を関数のモーメントの線形重ね合わせとして表現します。この観点は、確率論に応用されています。 σ{\displaystyle \sigma }τ{\displaystyle \tau }

正式には、フーリエ変換は、以下に説明する条件が満たされる場合、 虚数引数s = [ 33 ] [ 34 ]を持つ両側ラプラス変換を評価することと同等である。f^(ω)=F{f(t)}=L{f(t)}|s=iω=F(s)|s=iω=eiωtf(t)dt .{\displaystyle {\begin{aligned}{\hat {f}}(\omega )&={\mathcal {F}}\{f(t)\}\\[4pt]&={\mathcal {L}}\{f(t)\}|_{s=i\omega }=F(s)|_{s=i\omega }\\[4pt]&=\int _{-\infty }^{\infty }e^{-i\omega t}f(t)\,dt~.\end{aligned}}}

このフーリエ変換の規則(フーリエ変換 § その他の規則)では、⁠f^3(ω){\displaystyle {\hat {f}}_{3}(\omega )}の係数が必要である1/逆フーリエ変換について。ラプラス変換とフーリエ変換の関係は、信号または動的システムの周波数スペクトルを決定するためによく使用されます。

上記の関係は、F ( s )の収束領域(ROC)に虚軸σ =0が含まれる場合にのみ有効です。

例えば、関数f ( t ) = cos( ω 0 t )にはラプラス変換F ( s ) = s /( s 2 + ω 0 2 )があり、その ROC はRe( s ) > 0です。s = 0はF ( s )の極である ため、F ( s )s = を代入しても、ディラックのデルタ関数δ ( ω ± ω 0 )に比例する項を含むf ( t ) u ( t )のフーリエ変換は得られません。

しかし、 という関係は、 はるかに弱い条件下でも成り立ちます。例えば、上の例では、 の極限が測度の弱極限として理解されている限り、この関係が成り立ちます(曖昧位相を参照)。境界上の関数のラプラス変換の極限とフーリエ変換を関連付ける一般的な条件は、ペイリー・ウィーナーの定理の形をとります。 limσ0+F(σ+iω)=f^(ω){\displaystyle \lim _{\sigma \to 0^{+}}F(\sigma +i\omega )={\hat {f}}(\omega )}

メリン変換

メリン変換とその逆変換は、単純な変数の変更によって両側ラプラス変換に関連付けられます。

メリン変換で θ = e t と設定すると、両側ラプラス変換が得られます。 G(s)=M{g(θ)}=0θsg(θ)dθθ{\displaystyle G(s)={\mathcal {M}}\{g(\theta )\}=\int _{0}^{\infty }\theta ^{s}g(\theta )\,{\frac {d\theta }{\theta }}}

Z変換

片側 Z 変換は、理想的にサンプリングされた信号のラプラス変換に代入するだけです。 ここで、T = 1/ f sサンプリング間隔(時間単位、たとえば秒)、 f sはサンプリング レート(サンプル/秒またはヘルツ) です。z=defesT,{\displaystyle z{\stackrel {\mathrm {def} }{{}={}}}e^{sT},}

をサンプリングインパルス列(ディラックコムとも呼ばれる)とし、 を 連続時間x ( t )のサンプリング表現とする。ΔT(t) =def n=0δ(tnT){\displaystyle \Delta _{T}(t)\ {\stackrel {\mathrm {def} }{=}}\ \sum _{n=0}^{\infty }\delta (t-nT)}xq(t)=defx(t)ΔT(t)=x(t)n=0δ(tnT)=n=0x(nT)δ(tnT)=n=0x[n]δ(tnT){\displaystyle {\begin{aligned}x_{q}(t)&{\stackrel {\mathrm {def} }{{}={}}}x(t)\Delta _{T}(t)=x(t)\sum _{n=0}^{\infty }\delta (t-nT)\\&=\sum _{n=0}^{\infty }x(nT)\delta (t-nT)=\sum _{n=0}^{\infty }x[n]\delta (t-nT)\end{aligned}}}x[n]=defx(nT) .{\displaystyle x[n]{\stackrel {\mathrm {def} }{{}={}}}x(nT)~.}

サンプリングされた信号x q ( t )のラプラス変換は Xq(s)=0xq(t)estdt=0n=0x[n]δ(tnT)estdt=n=0x[n]0δ(tnT)estdt=n=0x[n]ensT .{\displaystyle {\begin{aligned}X_{q}(s)&=\int _{0^{-}}^{\infty }x_{q}(t)e^{-st}\,dt\\&=\int _{0^{-}}^{\infty }\sum _{n=0}^{\infty }x[n]\delta (t-nT)e^{-st}\,dt\\&=\sum _{n=0}^{\infty }x[n]\int _{0^{-}}^{\infty }\delta (t-nT)e^{-st}\,dt\\&=\sum _{n=0}^{\infty }x[n]e^{-nsT}~.\end{aligned}}}

これは、 ze sT を代入した離散関数x [ n ] の片側Z変換の正確な定義です。 X(z)=n=0x[n]zn{\displaystyle X(z)=\sum _{n=0}^{\infty }x[n]z^{-n}}

最後の2つの式を比較すると、片側Z変換とサンプリングされた信号のラプラス変換の関係が分かります。 Xq(s)=X(z)|z=esT.{\displaystyle X_{q}(s)=X(z){\Big |}_{z=e^{sT}}.}

Z 変換とラプラス変換の類似性は、時間スケール計算の理論でさらに詳しく説明されています。

ボレル変換

ボレル変換 の積分形は、指​​数型整関数f に対するラプラス変換の特殊なケースであり、これは、 ある定数ABに対して成り立つことを意味します。一般化ボレル変換では、指数型ではない関数を変換する際に、指数関数ではなく異なる重み関数を使用することができます。ナクビンの定理は、ボレル変換が適切に定義されるための必要十分条件を示しています。 F(s)=0f(z)eszdz{\displaystyle F(s)=\int _{0}^{\infty }f(z)e^{-sz}\,dz}|f(z)|AeB|z|{\displaystyle |f(z)|\leq Ae^{B|z|}}

基本的な関係

通常のラプラス変換は両側変換の特殊なケースとして表すことができ、両側変換は2つの片側変換の和として表すことができるため、ラプラス変換、フーリエ変換、メリン変換、Z変換の理論は本質的に同じ主題である。しかしながら、これら4つの主要な積分変換にはそれぞれ異なる視点と異なる固有の問題が伴う。

選択されたラプラス変換の表

次の表は、単一変数の多くの一般的な関数のラプラス変換を示しています。[ 35 ] [ 36 ]定義と説明については、表の最後に ある説明を参照してください。

ラプラス変換は線形演算子なので、

  • 和のラプラス変換は、各項のラプラス変換の合計です。L{f(t)+g(t)}=L{f(t)}+L{g(t)}{\displaystyle {\mathcal {L}}\{f(t)+g(t)\}={\mathcal {L}}\{f(t)\}+{\mathcal {L}}\{g(t)\}}
  • 関数の倍数のラプラス変換は、その関数のラプラス変換の倍数です。L{af(t)}=aL{f(t)}{\displaystyle {\mathcal {L}}\{af(t)\}=a{\mathcal {L}}\{f(t)\}}

この線形性と、さまざまな三角関数双曲線関数、複素数 (など) の特性や恒等式を使用すると、定義を直接使用するよりも速く、いくつかのラプラス変換を他のラプラス変換から取得できます。

片側ラプラス変換は、時間領域が非負の実数である関数を入力として受け取ります。そのため、以下の表の時間領域関数はすべてヘビサイドステップ関数u ( t )の倍数になります。

表中の時間遅延τを伴うエントリは、因果的(つまりτ > 0 )であることが要求される。因果システムとは、 t = 0より前のすべての時刻tにおいてインパルス応答h ( t )がゼロとなるシステムである。一般に、因果システムの収束領域は反因果システムの収束領域と同じではない。

選択されたラプラス変換
関数 時間領域 f(t)=L1{F(s)}{\displaystyle f(t)={\mathcal {L}}^{-1}\{F(s)\}}ラプラスs領域 F(s)=L{f(t)}{\displaystyle F(s)={\mathcal {L}}\{f(t)\}}収束領域 参照
単位インパルス δ(t) {\displaystyle \delta (t)\ }1{\displaystyle 1}すべてのs検査
遅延インパルス δ(tτ) {\displaystyle \delta (t-\tau )\ }eτs {\displaystyle e^{-\tau s}\ }すべてのs単位インパルス の時間シフト
単位ステップ u(t) {\displaystyle u(t)\ }1s{\displaystyle {1 \over s}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}単位インパルスを積分する
遅延ユニットステップ u(tτ) {\displaystyle u(t-\tau )\ }1seτs{\displaystyle {\frac {1}{s}}e^{-\tau s}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}単位ステップ の時間シフト
遅延関数と遅延ステップの積 f(tτ)u(tτ){\displaystyle f(t-\tau )u(t-\tau )}esτL{f(t)}{\displaystyle e^{-s\tau }{\mathcal {L}}\{f(t)\}}u-置換、u=tτ{\displaystyle u=t-\tau }
長方形のインパルス u(t)u(tτ){\displaystyle u(t)-u(t-\tau )}1s(1eτs){\displaystyle {\frac {1}{s}}(1-e^{-\tau s})}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}
ランプtu(t) {\displaystyle t\cdot u(t)\ }1s2{\displaystyle {\frac {1}{s^{2}}}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}単位インパルスを2回 ​​積分する
n乗(整数nの場合) tnu(t){\displaystyle t^{n}\cdot u(t)}n!sn+1{\displaystyle {n! \over s^{n+1}}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}n > −1単位ステップをn回 積分する
q乗(複素数qの場合) tqu(t){\displaystyle t^{q}\cdot u(t)}Γ(q+1)sq+1{\displaystyle {\operatorname {\Gamma } (q+1) \over s^{q+1}}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}Re(q)>1{\displaystyle \operatorname {Re} (q)>-1}[ 37 ] [ 38 ]
n乗根 tnu(t){\displaystyle {\sqrt[{n}]{t}}\cdot u(t)}1s1n+1Γ(1n+1){\displaystyle {1 \over s^{{\frac {1}{n}}+1}}\operatorname {\Gamma } \left({\frac {1}{n}}+1\right)}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}上記でq = 1/ nと設定します。
周波数シフトによるntneαtu(t){\displaystyle t^{n}e^{-\alpha t}\cdot u(t)}n!(s+α)n+1{\displaystyle {\frac {n!}{(s+\alpha )^{n+1}}}}Re(s)>α{\displaystyle \operatorname {Re} (s)>-\alpha }単位ステップを積分し、周波数シフトを適用する
周波数シフトによる n乗遅延(tτ)neα(tτ)u(tτ){\displaystyle (t-\tau )^{n}e^{-\alpha (t-\tau )}\cdot u(t-\tau )}n!eτs(s+α)n+1{\displaystyle {\frac {n!\cdot e^{-\tau s}}{(s+\alpha )^{n+1}}}}Re(s)>α{\displaystyle \operatorname {Re} (s)>-\alpha }単位ステップを積分し、周波数シフトを適用し、時間シフトを適用する
指数関数的減衰eαtu(t){\displaystyle e^{-\alpha t}\cdot u(t)}1s+α{\displaystyle {1 \over s+\alpha }}Re(s)>α{\displaystyle \operatorname {Re} (s)>-\alpha }単位ステップ の周波数シフト
両側指数関数的減衰(両側変換のみ) eα|t| {\displaystyle e^{-\alpha |t|}\ }2αα2s2{\displaystyle {2\alpha \over \alpha ^{2}-s^{2}}}α<Re(s)<α{\displaystyle -\alpha <\operatorname {Re} (s)<\alpha }単位ステップ の周波数シフト
指数関数的アプローチ (1eαt)u(t) {\displaystyle (1-e^{-\alpha t})\cdot u(t)\ }αs(s+α){\displaystyle {\frac {\alpha }{s(s+\alpha )}}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}単位ステップマイナス指数関数的減衰
正弦sin(ωt)u(t) {\displaystyle \sin(\omega t)\cdot u(t)\ }ωs2+ω2{\displaystyle {\omega \over s^{2}+\omega ^{2}}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}[ 39 ]
余弦cos(ωt)u(t) {\displaystyle \cos(\omega t)\cdot u(t)\ }ss2+ω2{\displaystyle {s \over s^{2}+\omega ^{2}}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}[ 39 ]
双曲線正弦sinh(αt)u(t) {\displaystyle \sinh(\alpha t)\cdot u(t)\ }αs2α2{\displaystyle {\alpha \over s^{2}-\alpha ^{2}}}Re(s)>|α|{\displaystyle \operatorname {Re} (s)>\left|\alpha \right|}[ 40 ]
双曲線余弦cosh(αt)u(t) {\displaystyle \cosh(\alpha t)\cdot u(t)\ }ss2α2{\displaystyle {s \over s^{2}-\alpha ^{2}}}Re(s)>|α|{\displaystyle \operatorname {Re} (s)>\left|\alpha \right|}[ 40 ]
指数関数的に減少する正弦波 eαtsin(ωt)u(t) {\displaystyle e^{-\alpha t}\sin(\omega t)\cdot u(t)\ }ω(s+α)2+ω2{\displaystyle {\omega \over (s+\alpha )^{2}+\omega ^{2}}}Re(s)>α{\displaystyle \operatorname {Re} (s)>-\alpha }[ 39 ]
指数関数的に減少する余弦波 eαtcos(ωt)u(t) {\displaystyle e^{-\alpha t}\cos(\omega t)\cdot u(t)\ }s+α(s+α)2+ω2{\displaystyle {s+\alpha \over (s+\alpha )^{2}+\omega ^{2}}}Re(s)>α{\displaystyle \operatorname {Re} (s)>-\alpha }[ 39 ]
自然対数ln(t)u(t){\displaystyle \ln(t)\cdot u(t)}1s[ln(s)+γ]{\displaystyle -{1 \over s}\left[\ln(s)+\gamma \right]}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}[ 40 ]
n次の第一種ベッセル関数Jn(ωt)u(t){\displaystyle J_{n}(\omega t)\cdot u(t)}(s2+ω2s)nωns2+ω2{\displaystyle {\frac {\left({\sqrt {s^{2}+\omega ^{2}}}-s\right)^{\!n}}{\omega ^{n}{\sqrt {s^{2}+\omega ^{2}}}}}}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}n > −1[ 41 ]
誤差関数erf(t)u(t){\displaystyle \operatorname {erf} (t)\cdot u(t)}1ses2/4(1erfs2){\displaystyle {\frac {1}{s}}e^{s^{2}/4}\!\left(1-\operatorname {erf} {\frac {s}{2}}\right)}Re(s)>0{\displaystyle \operatorname {Re} (s)>0}[ 41 ]
説明文:

s領域等価回路とインピーダンス

ラプラス変換は回路解析においてよく用いられ、回路要素のs領域への簡単な変換が可能です。回路要素は、位相インピーダンスに非常によく似たインピーダンスに変換できます。

同等のものを以下にまとめます。

s領域等価回路
s領域等価回路

抵抗器は時間領域とs領域で全く同じであることに注意してください。回路要素に初期条件がある場合、ソースが挿入されます。例えば、コンデンサに初期電圧がかかる場合や、インダクタに初期電流が流れる場合、s領域に挿入されたソースはそれを考慮します。

電流源と電圧源の等価物は、上記の表の変換から簡単に導き出されます。

例と応用

ラプラス変換は工学物理学で頻繁に用いられます。線形時不変システムの出力は、その単位インパルス応答を入力信号と畳み込むことで計算できます。この計算をラプラス空間で実行すると、畳み込みは乗算に変換されます。乗算は代数形式であるため、より簡単に解くことができます。詳細については、制御理論を参照してください。ラプラス変換は、多くの関数クラスで可逆です。システムへの入力または出力の簡単な数学的または機能的記述が与えられればラプラス変換は代替の機能的記述を提供し、多くの場合、システムの動作を解析するプロセスや、一連の仕様に基づいて新しいシステムを合成するプロセスを簡素化します。[ 42 ]

ラプラス変換は微分方程式を解くのにも利用でき、機械工学電気工学で広く用いられています。ラプラス変換は線形微分方程式を代数方程式に変換し、代数の正式な規則を用いて解くことができます。そして、逆ラプラス変換を適用することで、元の微分方程式を解くことができます。イギリスの電気技師オリバー・ヘヴィサイドは、ラプラス変換を用いない同様の手法を初めて提案しました。この結果得られた演算法は、ヘヴィサイド計算と呼ばれています。

不定積分の評価

⁠ ⁠L{f(t)}=F(s){\displaystyle {\mathcal {L}}\left\{f(t)\right\}=F(s)}とします。すると(上の表を参照) sL{f(t)t}=s0f(t)testdt=0f(t)estdt=F(s){\displaystyle \partial _{s}{\mathcal {L}}\left\{{\frac {f(t)}{t}}\right\}=\partial _{s}\int _{0}^{\infty }{\frac {f(t)}{t}}e^{-st}\,dt=-\int _{0}^{\infty }f(t)e^{-st}dt=-F(s)}

そこから次のものが得られます: L{f(t)t}=sF(p)dp.{\displaystyle {\mathcal {L}}\left\{{\frac {f(t)}{t}}\right\}=\int _{s}^{\infty }F(p)\,dp.}

極限⁠ ⁠s0{\displaystyle s\rightarrow 0}においては、 極限の交換が正当化され得ることが示される。これは終値定理の結果としてしばしば可能である。交換が正当化されない場合でも、計算は示唆的となることがある。例えば、a ≠ 0 ≠ bの場合、形式的に進めると、 0f(t)tdt=0F(p)dp,{\displaystyle \int _{0}^{\infty }{\frac {f(t)}{t}}\,dt=\int _{0}^{\infty }F(p)\,dp,}0cos(at)cos(bt)tdt=0(pp2+a2pp2+b2)dp=[12lnp2+a2p2+b2]0=12lnb2a2=ln|ba|.{\displaystyle {\begin{aligned}\int _{0}^{\infty }{\frac {\cos(at)-\cos(bt)}{t}}\,dt&=\int _{0}^{\infty }\left({\frac {p}{p^{2}+a^{2}}}-{\frac {p}{p^{2}+b^{2}}}\right)\,dp\\[6pt]&=\left[{\frac {1}{2}}\ln {\frac {p^{2}+a^{2}}{p^{2}+b^{2}}}\right]_{0}^{\infty }={\frac {1}{2}}\ln {\frac {b^{2}}{a^{2}}}=\ln \left|{\frac {b}{a}}\right|.\end{aligned}}}

コンデンサの複素インピーダンス

電気回路理論では、コンデンサに流れる電流は、静電容量と電位の変化率に比例します(SI単位系における式を参照)。これは記号的に、微分方程式で表されます。 ここで、 Cはコンデンサの静電容量、i = i ( t )はコンデンサを流れる電流の時間の関数、v = v ( t )はコンデンサの端子間 電圧の時間の関数です。i=Cdvdt,{\displaystyle i=C{dv \over dt},}

この方程式をラプラス変換すると、 次 の式が得られる 。 I(s)=C(sV(s)V0),{\displaystyle I(s)=C(sV(s)-V_{0}),}I(s)=L{i(t)},V(s)=L{v(t)},{\displaystyle {\begin{aligned}I(s)&={\mathcal {L}}\{i(t)\},\\V(s)&={\mathcal {L}}\{v(t)\},\end{aligned}}}V0=v(0).{\displaystyle V_{0}=v(0).}

V ( s )を解くと、 V(s)=I(s)sC+V0s.{\displaystyle V(s)={I(s) \over sC}+{V_{0} \over s}.}

複素インピーダンスZオーム)の定義は、初期状態V 0をゼロに 保持しながら、複素電圧Vを複素電流Iで割った比です。Z(s)=V(s)I(s)|V0=0.{\displaystyle Z(s)=\left.{V(s) \over I(s)}\right|_{V_{0}=0}.}

この定義と前の式を用いると、次式が得られます。 これはコンデンサの複素インピーダンスの正しい表現です。さらに、ラプラス変換は制御理論において広く応用されています。 Z(s)=1sC,{\displaystyle Z(s)={\frac {1}{sC}},}

インパルス応答

伝達関数を持つ線形時間不変システムを考えるH(s)=1(s+α)(s+β).{\displaystyle H(s)={\frac {1}{(s+\alpha )(s+\beta )}}.}

インパルス応答は、この伝達関数の逆ラプラス変換です。 h(t)=L1{H(s)}.{\displaystyle h(t)={\mathcal {L}}^{-1}\{H(s)\}.}

部分分数展開

この逆変換を評価するために、まず部分分数展開法を用いて H ( s )を展開する。1(s+α)(s+β)=Ps+α+Rs+β.{\displaystyle {\frac {1}{(s+\alpha )(s+\beta )}}={P \over s+\alpha }+{R \over s+\beta }.}

未知定数PRは、伝達関数の対応する極に位置する留数です。各留数は、その特異点が伝達関数全体の形状にどの程度寄与するかを表します。

留数定理によれば、逆ラプラス変換は極とその留数のみに依存します。留数Pを求めるには、方程式の両辺にs + αを掛けて、 1s+β=P+R(s+α)s+β.{\displaystyle {\frac {1}{s+\beta }}=P+{R(s+\alpha ) \over s+\beta }.}

そして、 s = − αとすると、 Rからの寄与は消え、残るのは P=1s+β|s=α=1βα.{\displaystyle P=\left.{1 \over s+\beta }\right|_{s=-\alpha }={1 \over \beta -\alpha }.}

同様に、留数Rは次のように与えられる。 R=1s+α|s=β=1αβ.{\displaystyle R=\left.{1 \over s+\alpha }\right|_{s=-\beta }={1 \over \alpha -\beta }.}

であることに注意し 、 RPを展開したH ( s )の式に 代入すると、 R=1βα=P{\displaystyle R={-1 \over \beta -\alpha }=-P}H(s)=(1βα)(1s+α1s+β).{\displaystyle H(s)=\left({\frac {1}{\beta -\alpha }}\right)\cdot \left({1 \over s+\alpha }-{1 \over s+\beta }\right).}

最後に、線形性特性と指数関数的減少の既知の変換(上記のラプラス変換表項目3を参照)を使用して、 Hsの逆ラプラス変換を実行して、 システムのインパルス応答を 取得 できます。h(t)=L1{H(s)}=1βα(eαteβt),{\displaystyle h(t)={\mathcal {L}}^{-1}\{H(s)\}={\frac {1}{\beta -\alpha }}\left(e^{-\alpha t}-e^{-\beta t}\right),}

畳み込み

畳み込みの性質を利用して、システムが伝達関数1/( s + α )1/( s + β )を持つ一連のフィルタであるかのように同じ結果を得ることができます。つまり、の逆関数 は H(s)=1(s+α)(s+β)=1s+α1s+β{\displaystyle H(s)={\frac {1}{(s+\alpha )(s+\beta )}}={\frac {1}{s+\alpha }}\cdot {\frac {1}{s+\beta }}}L1{1s+α}L1{1s+β}=eαteβt=0teαxeβ(tx)dx=eαteβtβα.{\displaystyle {\mathcal {L}}^{-1}\!\left\{{\frac {1}{s+\alpha }}\right\}*{\mathcal {L}}^{-1}\!\left\{{\frac {1}{s+\beta }}\right\}=e^{-\alpha t}*e^{-\beta t}=\int _{0}^{t}e^{-\alpha x}e^{-\beta (t-x)}\,dx={\frac {e^{-\alpha t}-e^{-\beta t}}{\beta -\alpha }}.}

位相遅延

時間関数 ラプラス変換
sin(ωt+φ){\displaystyle \sin {(\omega t+\varphi )}}ssin(φ)+ωcos(φ)s2+ω2{\displaystyle {\frac {s\sin(\varphi )+\omega \cos(\varphi )}{s^{2}+\omega ^{2}}}}
cos(ωt+φ){\displaystyle \cos {(\omega t+\varphi )}}scos(φ)ωsin(φ)s2+ω2.{\displaystyle {\frac {s\cos(\varphi )-\omega \sin(\varphi )}{s^{2}+\omega ^{2}}}.}

ラプラス変換から始めて、 分数の項を並べ替えることで逆変換を求めます。 X(s)=ssin(φ)+ωcos(φ)s2+ω2{\displaystyle X(s)={\frac {s\sin(\varphi )+\omega \cos(\varphi )}{s^{2}+\omega ^{2}}}}X(s)=ssin(φ)s2+ω2+ωcos(φ)s2+ω2=sin(φ)(ss2+ω2)+cos(φ)(ωs2+ω2).{\displaystyle {\begin{aligned}X(s)&={\frac {s\sin(\varphi )}{s^{2}+\omega ^{2}}}+{\frac {\omega \cos(\varphi )}{s^{2}+\omega ^{2}}}\\&=\sin(\varphi )\left({\frac {s}{s^{2}+\omega ^{2}}}\right)+\cos(\varphi )\left({\frac {\omega }{s^{2}+\omega ^{2}}}\right).\end{aligned}}}

これで、項の逆ラプラス変換を行うことができます。 x(t)=sin(φ)L1{ss2+ω2}+cos(φ)L1{ωs2+ω2}=sin(φ)cos(ωt)+cos(φ)sin(ωt).{\displaystyle {\begin{aligned}x(t)&=\sin(\varphi ){\mathcal {L}}^{-1}\left\{{\frac {s}{s^{2}+\omega ^{2}}}\right\}+\cos(\varphi ){\mathcal {L}}^{-1}\left\{{\frac {\omega }{s^{2}+\omega ^{2}}}\right\}\\&=\sin(\varphi )\cos(\omega t)+\cos(\varphi )\sin(\omega t).\end{aligned}}}

これは単に引数の合計の正弦であり、次のようになります。 x(t)=sin(ωt+φ).{\displaystyle x(t)=\sin(\omega t+\varphi ).}

同様の論理を適用すれば、 L1{scosφωsinφs2+ω2}=cos(ωt+φ).{\displaystyle {\mathcal {L}}^{-1}\left\{{\frac {s\cos \varphi -\omega \sin \varphi }{s^{2}+\omega ^{2}}}\right\}=\cos {(\omega t+\varphi )}.}

統計力学

統計力学では、状態密度のラプラス変換は分配関数を定義する。[ 43 ]すなわち、標準的な分配関数は次のように与えられ 、逆関数は次のように与えられる。 g(E){\displaystyle g(E)}Z(β){\displaystyle Z(\beta )}Z(β)=0eβEg(E)dE{\displaystyle Z(\beta )=\int _{0}^{\infty }e^{-\beta E}g(E)\,dE}g(E)=12πiβ0iβ0+ieβEZ(β)dβ{\displaystyle g(E)={\frac {1}{2\pi i}}\int _{\beta _{0}-i\infty }^{\beta _{0}+i\infty }e^{\beta E}Z(\beta )\,d\beta }

天文スペクトルからの空間構造(時間構造ではない)

ラプラス変換とその逆変換の広く一般的な適用性は、天文学における応用によって実証されています。この応用では、時間領域とスペクトル (周波数領域)を関連付けるのではなく、磁束密度スペクトルが与えられた場合に、点として分解するには遠すぎる天体の無線周波数熱放射源の物質の空間分布に関する情報が提供されます。

物体の特定の特性、例えば球形や一定温度を仮定し、物体のスペクトルに対して逆ラプラス変換を実行することに基づく計算によって、スペクトルと一致する、物体内の物質の分布(中心からの距離の関数としての密度)の唯一の可能なモデルを生成することができる。 [ 44 ]物体の構造に関する独立した情報が利用可能な場合、逆ラプラス変換法は良好な一致を示すことがわかっている。

誕生と死のプロセス

確率でステップが発生するランダムウォークを考えます。[ 45 ]また、時間ステップがパラメータ⁠を持つポアソン過程であると仮定します。すると、ウォークが時間に格子点に存在する確率は、次の式で表されます。 これは積分方程式 系(または微分方程式系と同等)につながります。しかし、これは畳み込み方程式系であるため、ラプラス変換によって次の線型方程式系に変換されます。 これ は 標準的な方法で解くことができます。 {+1,1}{\displaystyle \{+1,-1\}}p,q=1p{\displaystyle p,q=1-p}λ{\displaystyle \lambda }n{\displaystyle n}t{\displaystyle t}Pn(t)=0tλeλ(ts)(pPn1(s)+qPn+1(s))ds(+eλtwhen n=0).{\displaystyle P_{n}(t)=\int _{0}^{t}\lambda e^{-\lambda (t-s)}(pP_{n-1}(s)+qP_{n+1}(s))\,ds\quad (+e^{-\lambda t}\quad {\text{when}}\ n=0).}πn(s)=L(Pn)(s),{\displaystyle \pi _{n}(s)={\mathcal {L}}(P_{n})(s),}πn(s)=λλ+s(pπn1(s)+qπn+1(s))(+1λ+swhen n=0){\displaystyle \pi _{n}(s)={\frac {\lambda }{\lambda +s}}(p\pi _{n-1}(s)+q\pi _{n+1}(s))\quad (+{\frac {1}{\lambda +s}}\quad {\text{when}}\ n=0)}

タウバー理論

上の測度のラプラス変換は次のように与えられる。 ⁠ が小さい場合、指数関数的に減少する積分関数は、領域のより大きな部分集合への測度の集中に対してより敏感になることは 直感的に明らかである。 これをより正確にするために、分布関数を導入する。 正式には、次のような種類の極限が期待される。 タウバー定理は、 としてのラプラス変換の漸近挙動と、としての分布の漸近挙動を関連付ける定理である。 したがって、これらは確率統計の漸近式で重要であり、スペクトル側には推論がより簡単な漸近挙動があることが多い。[ 45 ]μ{\displaystyle \mu }[0,){\displaystyle [0,\infty )}Lμ(s)=0estdμ(t).{\displaystyle {\mathcal {L}}\mu (s)=\int _{0}^{\infty }e^{-st}d\mu (t).}s>0{\displaystyle s>0}μ{\displaystyle \mu }M(t)=μ([0,t)).{\displaystyle M(t)=\mu ([0,t)).}lims0+Lμ(s)=limtM(t).{\displaystyle \lim _{s\to 0^{+}}{\mathcal {L}}\mu (s)=\lim _{t\to \infty }M(t).}s0+{\displaystyle s\to 0^{+}}μ{\displaystyle \mu }t{\displaystyle t\to \infty }

注目すべき2つのタウバー定理は、ハーディ・リトルウッドのタウバー定理ウィーナーのタウバー定理である。ウィーナーの定理は、池原のタウバー定理を一般化したものであり、以下の命題を成す。

A ( x ) をxの非負単調非減少関数とし、 0 ≤ x < ∞で定義されるものとする。ℜ ( s ) > 1 において関数f ( s ) に収束し、ある非負数cに対して、 ℜ( s ) ≥ 1において連続関数 として拡張されるとする。このとき、 x が無限大に向かうときのe x A ( x )の極限は c に等しい。 f(s)=0A(x)exsdx{\displaystyle f(s)=\int _{0}^{\infty }A(x)e^{-xs}\,dx}f(s)cs1{\displaystyle f(s)-{\frac {c}{s-1}}}

この主張は特にリーマンゼータ関数対数微分に適用することができ、素数定理を証明するための非常に簡潔な方法を提供する。[ 46 ]

参照

注記

  1. ^ Lynn, Paul A. (1986)、「ラプラス変換とz変換」、Electronic Signals and Systems、ロンドン:Macmillan Education UK、pp.  225– 272、doi10.1007/978-1-349-18461-3_6ISBN 978-0-333-39164-8ラプラス変換とZ変換はフーリエ変換と密接に関連しています。ラプラス変換はフーリエ変換よりも範囲が広く、自動制御システムを含む連続回路やシステムを記述するためにエンジニアに広く使用されています。
  2. ^ 「微分方程式 - ラプラス変換」Pauls Online Math Notes 、 2020年8月8日閲覧
  3. ^ a b Weisstein, Eric W.、「ラプラス変換」Wolfram MathWorld 、 2020年8月8日閲覧。
  4. ^ディーキン、マイケル・AB(1981)、「ラプラス変換の発展、1737-1937:I.オイラーからスピッツァーまで、1737-1880」正確科学史アーカイブ25(4):343-390ISSN 0003-9519 
  5. ^ "Des Fonctions génératrices" [生成関数について]、Théorie Analytique des Probabilités [分析的確率論] (フランス語) (第 2 版)、パリ、1​​814 年、第 I 章 2-20 章
  6. ^ Jaynes, ET (Edwin T.) (2003),確率論:科学の論理、Bretthorst, G. Larry、ケンブリッジ、イギリス:ケンブリッジ大学出版局、ISBN 0511065892OCLC  57254076
  7. ^アベル、ニールス H. (1820)、「Sur les fonctions génératrices et leurs déterminantes」、āuvres Complètes (フランス語)、vol. II (1839 年出版)、 77 ~ 88ページ 1881年版
  8. ^オイラー 1744オイラー 1753オイラー 1769
  9. ^ラグランジュ 1773
  10. ^グラッタン・ギネス 1997、260ページ
  11. ^グラッタン・ギネス 1997、261ページ
  12. ^グラッタン・ギネス 1997年、261~262ページ
  13. ^グラッタン・ギネス 1997年、262~266ページ
  14. ^ヘヴィサイド、オリバー(2008年1月)微分変換による定積分の解法」電磁気理論第3巻、ロンドン、セクション526、ISBN 9781605206189{{citation}}: CS1 maint: location missing publisher (link)
  15. ^ Edwards, HM (1974), Riemann's Zeta Function , New York: Academic Press, ISBN 0-12-232750-0Zbl  0315.10035
  16. ^ガードナー、マレー F.; バーンズ、ジョン L. (1942)、「ラプラス変換による線形システムの過渡現象の研究」、ニューヨーク:ワイリー付録C
  17. ^ Lerch、Mathias (1903)、「Sur un point de la théorie des fonctions génératrices d'Abel」[逆転公式の証明]、Acta Mathematica (フランス語)、27 : 339–351doi : 10.1007/BF02421315hdl : 10338.dmlcz/501554
  18. ^ Bromwich, Thomas J. (1916)、「動的システムにおける正規座標」ロンドン数学会紀要15 : 401–448doi : 10.1112/plms/s2-15.1.401
  19. ^影響力のある本としては、ガードナー、マレー・F.、バーンズ、ジョン・L.(1942年)、ラプラス変換による線形システムの過渡現象の研究、ニューヨーク:ワイリーなどがある。
  20. ^ Doetsch、Gustav (1937)、Theorie und Anwendung der Laplacesche Transformation [ラプラス変換の理論と応用] (ドイツ語)、ベルリン: Springer1943年の翻訳
  21. ^フェラー 1971、§XIII.1。
  22. ^累積分布関数は確率密度関数の積分です。
  23. ^ミクシンスキ、ヤン(2014年7月14日)、オペレーショナル・カルキュラス、エルゼビア、ISBN 9781483278933
  24. ^ウィダー 1941、第 II 章、§1
  25. ^ウィダー 1941、第 VI 章、§2
  26. ^ Beffa, Federico (2024)、「分布のラプラス変換」弱非線形システム、Cham: Springer Nature Switzerland、pp.  75– 85、doi : 10.1007/978-3-031-40681-2_5ISBN 978-3-031-40680-52026年1月14日取得{{citation}}: CS1 maint: work parameter with ISBN (link)
  27. ^コーン&コーン 1967、226–227 ページ
  28. ^ Bracewell 2000、表14.1、p.385
  29. ^ GhostarchiveWayback Machineにアーカイブ: Mattuck, Arthur(2008年11月7日)「ラプラス変換の起源」YouTube
  30. ^ Widder 1941、38ページ。
  31. ^フェラー 1971、432ページ
  32. ^ローラン・シュワルツ(1966)、『物理科学のための数学』、アディソン・ウェズリー、224ページ。
  33. ^ Titchmarsh, E. (1986) [1948],フーリエ積分理論入門(第2版)、Clarendon Press、p. 6、ISBN 978-0-8284-0324-5
  34. ^タカチ 1953、93ページ
  35. ^ライリー、KF; ホブソン、MP; ベンス、SJ (2010)、物理学と工学のための数学的手法(第3版)、ケンブリッジ大学出版局、p. 455、ISBN 978-0-521-86153-3
  36. ^ Distefano, JJ; Stubberud, AR; Williams, IJ (1995)、「フィードバックシステムと制御」、Schaum's outlines(第2版)、McGraw-Hill、p. 78、ISBN 978-0-07-017052-0
  37. ^ Lipschutz, S.; Spiegel, MR; Liu, J. (2009), Mathematical Handbook of Formulas and Tables , Schaum's Outline Series (3rd ed.), McGraw-Hill, p. 183, ISBN 978-0-07-154855-7– 実数qの場合を示します。
  38. ^ http://mathworld.wolfram.com/LaplaceTransform.html – Wolfram Mathword は複素数qのケースを提供します
  39. ^ a b c dブレイスウェル 1978、227ページ。
  40. ^ a b cウィリアムズ 1973、88ページ。
  41. ^ a bウィリアムズ 1973、89ページ。
  42. ^コーン&コーン 1967、§8.1
  43. ^ RK Pathria; Paul Beal (1996)、統計力学(第2版)、Butterworth-Heinemann、p.  56ISBN 9780750624695
  44. ^ Salem, M.; Seaton, MJ (1974)、「I. 連続スペクトルと輝度等高線」、Monthly Notices of the Royal Astronomical Society167 : 493–510Bibcode : 1974MNRAS.167..493Sdoi : 10.1093/mnras/167.3.493Salem, M. (1974)、「II. 3次元モデル」、Monthly Notices of the Royal Astronomical Society167 : 511–516Bibcode : 1974MNRAS.167..511Sdoi : 10.1093/mnras/167.3.511
  45. ^ a bフェラー『確率論入門』第2巻、pp=479-483
  46. ^池原 誠(1931)、「ランダウの定理の解析的数論への拡張」、数学物理学ジャーナル10 ( 1–4 ): 1– 12、doi : 10.1002/sapm19311011Zbl 0001.12902 

参考文献

モダンな

歴史的

  • オイラー、L. ( 1744)、「De constructione aequationum」[方程式の構築]、オペラ オムニア、第 1 シリーズ (ラテン語)、22 : 150–161
  • オイラー、L. ( 1753)、「Methodus aequationes Differentiales」[微分方程式を解く方法]、オペラ オムニア、第 1 シリーズ (ラテン語)、22 : 181–213
  • Euler, L. (1992) [1769]、「Institutiones calculi integralis, Volume 2」[積分微積分の機関]、Opera Omnia、第 1 シリーズ (ラテン語)、12、バーゼル: Birkhäuser、ISBN 978-3764314743{{citation}}: CS1 maint: work parameter with ISBN (link)第3章~第5章
  • オイラー、レオンハルト(1769)、Institutiones calculiintegris [ Institutions of Integral Calculus ] (ラテン語)、vol. II、パリ: ペトロポリ、ch. 3–5、57–153ページ
  • グラッタン・ギネス、I(1997)「偏微分方程式のラプラス積分解」、ギリスピー、CC(編)、ピエール・シモン・ラプラス1749-1827:精密科学の生涯、プリンストン:プリンストン大学出版局、ISBN 978-0-691-01185-1
  • ラグランジュ、JL (1773)、Mémoire sur l'utilité de la methode、āuvres de Lagrange、vol. 2、171–234 ページ

さらに読む