数学 において、ラプラス変換は、 ピエール=シモン・ラプラス ()にちなんで名付けられ、実 変数関数 (通常は 、時間領域 )を複素変数 関数 (複素数値周波数領域、 s 領域 またはs 平面 とも呼ばれる)に変換する積分変換 です。これらの関数は、時間領域関数には小文字の記号、周波数領域関数には対応する大文字の記号を使用して表記されることがよくあります(例: および ) 。 t {\displaystyle t} s {\displaystyle s} × ( t ) {\displaystyle x(t)} X ( s ) {\displaystyle X(s)}
この変換は、時間領域における微分 と積分を 、ラプラス領域におけるはるかに簡単な乗算 と除算に変換するのに有用である( 対数が 乗算と除算を加算と減算に簡略化するのと同様に)。この変換は科学 技術において多くの応用があり、主に常微分方程式 と積分方程式を 代数 多項式方程式 に簡略化し、畳み込みを 乗算 に簡略化することで、線形微分方程式 [ 1 ] や動的システムを 解くためのツールとして用いられている。[ 2 ] [ 3 ]
例えば、ラプラス変換によって、単振動子 の方程式(フックの法則 )は、初期条件 と を含む代数方程式に変換され、未知関数について解くことができます。解いた後、逆ラプラス変換を用いて元の領域に変換することができます。これは、以下に 示すような表を参照することで容易に行うことができます。 × ″ ( t ) + け × ( t ) = 0 {\displaystyle x''(t)+kx(t)=0} s 2 X ( s ) − s × ( 0 ) − × ′ ( 0 ) + け X ( s ) = 0 、 {\displaystyle s^{2}X(s)-sx(0)-x'(0)+kX(s)=0,} × ( 0 ) {\displaystyle x(0)} × ′ ( 0 ) {\displaystyle x'(0)} X ( s ) 。 {\displaystyle X(s).}
ラプラス変換は(適切な関数 f {\displaystyle f} に対して)積分 によって定義されます。 ここで、sは 複素数 です。 L { f } ( s ) = ∫ 0 ∞ f ( t ) e − s t d t 、 {\displaystyle {\mathcal {L}}\{f\}(s)=\int _{0}^{\infty }f(t)e^{-st}\,dt,}
ラプラス変換は他の多くの変換と関連しています。メリン変換 と本質的に同じであり、フーリエ変換 と密接に関連しています。フーリエ変換とは異なり、関数のラプラス変換は解析関数であることが多いです。つまり、収束するべき 級数 を持ち、その係数は元の関数のモーメントを表します。さらに、 複素解析 の手法、特に線積分は 計算を簡素化するために使用できます。
歴史 ラプラス侯爵ピエール=シモン ラプラス変換は、確率論の 研究で同様の変換を使用した数学者 で天文学者の ピエール・シモン・ラプラス侯爵 にちなんで名付けられました。[ 4 ] [ 5 ] ラプラスは生成関数 の使用について広範囲に渡って著作を残し (1814年)、その結果としてラプラス変換の積分形式が自然に発展しました。[ 6 ]
ラプラスの母関数の使用は、現在Z変換 として知られているものと似ており、ニールス・ヘンリク・アーベル が議論した連続変数の ケースにはほとんど注意を払っていなかった。[ 7 ]
1744年より、レオンハルト・オイラーは 微分方程式の解として形式の積分を研究し 、特にガンマ関数 を導入した。[ 8 ] ジョゼフ=ルイ・ラグランジュはオイラーの崇拝者であり、 確率密度関数の 積分に関する研究において、 ラプラス変換に似た形式の表現を研究した 。 [ 9 ] [ 10 ] z = ∫ X ( × ) e 1つの × d × そして z = ∫ X ( × ) × あ d × {\displaystyle z=\int X(x)e^{ax}\,dx\quad {\text{ と }}\quad z=\int X(x)x^{A}\,dx} ∫ X ( × ) e − 1つの × 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 ] によって強調されており、ラプラス変換という名称は彼に由来すると思われる。
ℜ ( e − s t ) {\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 ) であり、これは次のように定義される片側変換である。
F ( s ) = ∫ 0 ∞ f ( t ) e − s t d t 、 {\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 ] 関数記法 では となる。これは、特に工学分野では、ダミー変数が関数 には現れないことを理解した上で、 と表記されることが多い。 F ( s ) = L { f } ( s ) {\displaystyle F(s)={\mathcal {L}}\{f\}(s)} F ( s ) = L { f ( t ) } {\displaystyle F(s)={\mathcal {L}}\{f(t)\}} t {\displaystyle t} F ( s ) {\displaystyle F(s)}
積分の意味は、対象となる関数の種類によって異なります。積分が存在するための必要条件は、f が [0, ∞) 上で局所的に積分 可能であることです。無限大で減少する、または指数型 ( )である局所的に積分可能な関数の場合、積分は(適切な)ルベーグ積分 であると理解できます。ただし、多くの用途では、これを∞ で条件付き収束する 不適切積分 と見なす必要があります。さらに一般的には、積分は弱い意味 で理解することができ、これについては以下で扱います。 | f ( t ) | ≤ あ e B | t | {\displaystyle \vert f(t)\vert \leq Ae^{B\vert t\vert }}
有限ボレル測度 μ のラプラス変換はルベーグ積分によって定義できる[ 21 ] L { μ } ( s ) = ∫ [ 0 , ∞ ) e − s t d μ ( t ) . {\displaystyle {\mathcal {L}}\{\mu \}(s)=\int _{[0,\infty )}e^{-st}\,d\mu (t).}
重要な特別なケースとして、μが 確率測度 、例えばディラックのデルタ関数 である場合が挙げられます。演算計算 では、測度のラプラス変換は、その測度が確率密度関数f から得られたものとして扱われることがよくあります。その場合、混乱を避けるため、 と書くことがよくあります。 ここで、 0 − の下限はの略記です。 L { f } ( s ) = ∫ 0 − ∞ f ( t ) e − s t d t , {\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 ) = ∫ − ∞ ∞ e − s t f ( t ) d t . {\displaystyle F(s)=\int _{-\infty }^{\infty }e^{-st}f(t)\,dt.} (式2 )
両側ラプラス変換の代替表記は、 F ではなく B { f } {\displaystyle {\mathcal {B}}\{f\}} です。
二つの可積分関数が同一のラプラス変換を持つのは、ルベーグ測度 零点の集合において両者が異なる場合に限られます。これは、変換の値域において逆変換が存在することを意味します。実際、ラプラス変換は可積分関数だけでなく、他の多くの関数空間においても、ある関数空間から別の関数空間への一対一写像 です。ただし、その値域を容易に特徴付けることは通常できません。
これが成り立つ典型的な関数空間としては、有界連続関数の空間、L ∞ (0, ∞) 、あるいはより一般的には(0, ∞) 上の緩和超関数の空間 が挙げられる。ラプラス変換もまた定義され、緩和超関数の適切な空間に対して単射である。
これらの場合、ラプラス変換の像は収束領域 における解析関数 の空間に存在します。逆ラプラス変換は、以下の複素積分で与えられ、これは様々な名前( ブロムウィッチ積分 、フーリエ・メリン積分、 メリンの逆公式 など) で知られています。
f ( t ) = L − 1 { F } ( t ) = 1 2 π i lim T → ∞ ∫ γ − i T γ + i T e s t F ( s ) d s , {\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 [ e − s X ] , {\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 ] F X ( x ) = L − 1 { 1 s E [ e − s X ] } ( x ) = L − 1 { 1 s L { 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 のラプラス変換F (s ) は極限が 存在する限り収束します。 lim R → ∞ ∫ 0 R f ( t ) e − s t d t {\displaystyle \lim _{R\to \infty }\int _{0}^{R}f(t)e^{-st}\,dt}
ラプラス変換は、積分が真のルベーグ積分として存在する場合、 絶対収束します 。ラプラス変換は通常、条件付き収束 と理解されており、前者の意味では収束しますが、後者の意味で収束しません。 ∫ 0 ∞ | f ( t ) e − s t | d t {\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 ) = ( s − s 0 ) ∫ 0 ∞ e − ( s − s 0 ) t β ( t ) d t , β ( u ) = ∫ 0 u e − s 0 t f ( t ) d t . {\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} } { s ∈ C | R e ( 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 ) = L − 1 { F ( s ) } , g ( t ) = L − 1 { G ( s ) } , {\displaystyle {\begin{aligned}f(t)&={\mathcal {L}}^{-1}\{F(s)\},\\g(t)&={\mathcal {L}}^{-1}\{G(s)\},\end{aligned}}}
次の表は片側ラプラス変換の性質の一覧である。[ 27 ]
初期値定理 f ( 0 + ) = lim s → ∞ s F ( s ) . {\displaystyle f(0^{+})=\lim _{s\to \infty }{sF(s)}.} 最終値定理 f ( ∞ ) = lim s → 0 s F ( s ) {\displaystyle f(\infty )=\lim _{s\to 0}{sF(s)}} 、のすべての極 が左半平面内にある場合。s F ( s ) {\displaystyle sF(s)} 終値定理は、 部分分数 分解(またはその他の難しい代数計算)を行うことなく長期的な挙動を導くため有用である。F ( s ) が 右手平面に極を持つ場合、または虚軸上に極を持つ場合(例えば、または )、この式の挙動は未定義である。f ( t ) = e t {\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 = 0 ∞ a ( n ) x n {\displaystyle \sum _{n=0}^{\infty }a(n)x^{n}} ∫ 0 ∞ f ( t ) x t d t {\displaystyle \int _{0}^{\infty }f(t)x^{t}\,dt}
べき乗の基数をxから e に変更すると 、∫ 0 ∞ f ( t ) ( e ln x ) t d t {\displaystyle \int _{0}^{\infty }f(t)\left(e^{\ln {x}}\right)^{t}\,dt}
例えば、すべての有界関数f に対して収束するには、 ln x < 0 を満たす必要があります。 − s = ln x と代入すると、ラプラス変換が得られます。 ∫ 0 ∞ f ( t ) e − s t d t {\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 } における解析関数に収束します。
瞬間との関係 これらの量 は関数f のモーメント である。 f の最初のn 個のモーメントが絶対収束する場合、積分 の下での繰り返し微分 により 、 が得られる。 これは確率論において特に重要であり、確率変数X のモーメントは期待値 で与えられる。 したがって、関係式は μ n = ∫ 0 ∞ t n f ( t ) d t {\displaystyle \mu _{n}=\int _{0}^{\infty }t^{n}f(t)\,dt} ( − 1 ) n ( L f ) ( n ) ( 0 ) = μ n . {\displaystyle (-1)^{n}({\mathcal {L}}f)^{(n)}(0)=\mu _{n}.} μ n = E [ X n ] {\displaystyle \mu _{n}=\operatorname {E} [X^{n}]} μ n = ( − 1 ) n d n d s n E [ e − s X ] ( 0 ) . {\displaystyle \mu _{n}=(-1)^{n}{\frac {d^{n}}{ds^{n}}}\operatorname {E} \left[e^{-sX}\right](0).}
ラプラス変換の微分性を利用して関数の導関数を求めると便利な場合が多い。これはラプラス変換の基本式から次のように導出できる。 そして、両側変換の場合は 、 L { f ( t ) } = ∫ 0 − ∞ e − s t f ( t ) d t = [ f ( t ) e − s t − s ] 0 − ∞ − ∫ 0 − ∞ e − s t − s f ′ ( t ) d t (by parts) = [ − f ( 0 − ) − s ] + 1 s L { 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 ) } = s ⋅ L { f ( t ) } − f ( 0 − ) , {\displaystyle {\mathcal {L}}\{f'(t)\}=s\cdot {\mathcal {L}}\{f(t)\}-f(0^{-}),} L { f ′ ( t ) } = s ∫ − ∞ ∞ e − s t f ( t ) d t = s ⋅ L { 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 ) } = s n ⋅ L { f ( t ) } − s n − 1 f ( 0 − ) − ⋯ − f ( n − 1 ) ( 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)}}
正の実軸上の積分を評価する ラプラス変換の有用な性質は、 の右近傍における の挙動と の左近傍におけるの減衰率に関する適切な仮定の下で、次のようになることです。上記の式は部分積分を変形したもので、演算子 と を と に置き換えたものです 。同等の 定式 化を証明しましょう。 ∫ 0 ∞ f ( x ) g ( x ) d x = ∫ 0 ∞ ( L f ) ( s ) ⋅ ( L − 1 g ) ( s ) d s {\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 } d d x {\displaystyle {\frac {d}{dx}}} ∫ d x {\displaystyle \int \,dx} L {\displaystyle {\mathcal {L}}} L − 1 {\displaystyle {\mathcal {L}}^{-1}} ∫ 0 ∞ ( L f ) ( x ) g ( x ) d x = ∫ 0 ∞ f ( s ) ( L g ) ( s ) d s . {\displaystyle \int _{0}^{\infty }({\mathcal {L}}f)(x)g(x)\,dx=\int _{0}^{\infty }f(s)({\mathcal {L}}g)(s)\,ds.}
を代入すると、左側は次のように変わります。 ただし、Fubini の定理が成り立つと仮定すると、積分の順序を逆にすると、必要な右側が得られます。 ( L f ) ( x ) = ∫ 0 ∞ f ( s ) e − s x d s {\displaystyle ({\mathcal {L}}f)(x)=\int _{0}^{\infty }f(s)e^{-sx}\,ds} ∫ 0 ∞ ∫ 0 ∞ f ( s ) g ( x ) e − s x d s d x , {\displaystyle \int _{0}^{\infty }\int _{0}^{\infty }f(s)g(x)e^{-sx}\,ds\,dx,}
この方法は、実微積分学の初等的な手法では計算が困難な積分を計算するのに使用できます。例えば、 ∫ 0 ∞ sin x x d x = ∫ 0 ∞ L ( 1 ) ( x ) sin x d x = ∫ 0 ∞ 1 ⋅ L ( sin ) ( x ) d x = ∫ 0 ∞ d x x 2 + 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 : ℝ → ℝの (一方的) ラプラス・スティールチェス変換は、ルベーグ・スティールチェス積分 によって定義されます。{ L ∗ g } ( s ) = ∫ 0 ∞ e − s t d g ( t ) . {\displaystyle \{{\mathcal {L}}^{*}g\}(s)=\int _{0}^{\infty }e^{-st}\,d\,g(t)~.}
関数g は有界変化 であると仮定する。gが f の不定積分 であるとすると、 g ( x ) = ∫ 0 x f ( t ) d t {\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 ) = L f ( s ) {\displaystyle F(s)={\mathcal {L}}f(s)} F ( σ + i τ ) = ∫ 0 ∞ f ( t ) e − σ t e − i τ t d t , {\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 ≥ 0 のt の関数の変換に制限されます。この制限の結果として、関数のラプラス変換は変数s の正則関数 になります。フーリエ変換とは異なり、分布のラプラス変換は一般に 、行儀の良い 関数です。複素変数の手法を使用して、ラプラス変換を直接研究することもできます。正則関数として、ラプラス変換はべき級数表現を持ちます。このべき級数は、関数を関数の モーメント の線形重ね合わせとして表現します。この観点は、確率論に応用されています。 σ {\displaystyle \sigma } τ {\displaystyle \tau }
正式には、フーリエ変換は、以下に説明する条件が満たされる場合、 虚数引数s = iω [ 33 ] [ 34 ]を持つ両側ラプラス変換を評価することと同等である。 f ^ ( ω ) = F { f ( t ) } = L { f ( t ) } | s = i ω = F ( s ) | s = i ω = ∫ − ∞ ∞ e − i ω t f ( t ) d t . {\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 / 2π 逆フーリエ変換について。ラプラス変換とフーリエ変換の関係は、信号 または動的システムの周波数スペクトル を決定するためによく使用されます。
上記の関係は、F ( s ) の収束領域(ROC)に虚軸σ =0が含まれる 場合にのみ 有効です。
例えば、関数f ( t ) = cos( ω 0 t ) にはラプラス変換F ( s ) = s /( s 2 + ω 0 2 ) があり、その ROC はRe( s ) > 0 です。s = iω 0は F ( s ) の極である ため、F ( s ) にs = iω を代入しても、ディラックのデルタ関数 δ ( ω ± ω 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 ∞ θ s g ( θ ) d θ θ {\displaystyle G(s)={\mathcal {M}}\{g(\theta )\}=\int _{0}^{\infty }\theta ^{s}g(\theta )\,{\frac {d\theta }{\theta }}}
片側 Z 変換は、理想的にサンプリングされた信号のラプラス変換に代入するだけです。 ここで、T = 1/ f s はサンプリング間隔 (時間単位、たとえば秒)、 f s はサンプリング レート (サンプル/秒 またはヘルツ ) です。z = d e f e s T , {\displaystyle z{\stackrel {\mathrm {def} }{{}={}}}e^{sT},}
をサンプリングインパルス列(ディラックコム とも呼ばれる)とし、 を 連続時間x ( t )のサンプリング表現とする。 Δ T ( t ) = d e f ∑ n = 0 ∞ δ ( t − n T ) {\displaystyle \Delta _{T}(t)\ {\stackrel {\mathrm {def} }{=}}\ \sum _{n=0}^{\infty }\delta (t-nT)} x q ( t ) = d e f x ( t ) Δ T ( t ) = x ( t ) ∑ n = 0 ∞ δ ( t − n T ) = ∑ n = 0 ∞ x ( n T ) δ ( t − n T ) = ∑ n = 0 ∞ x [ n ] δ ( t − n T ) {\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 ] = d e f x ( n T ) . {\displaystyle x[n]{\stackrel {\mathrm {def} }{{}={}}}x(nT)~.}
サンプリングされた信号x q ( t ) のラプラス変換は X q ( s ) = ∫ 0 − ∞ x q ( t ) e − s t d t = ∫ 0 − ∞ ∑ n = 0 ∞ x [ n ] δ ( t − n T ) e − s t d t = ∑ n = 0 ∞ x [ n ] ∫ 0 − ∞ δ ( t − n T ) e − s t d t = ∑ n = 0 ∞ x [ n ] e − n s T . {\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}}}
これは、 z → e sT を代入した離散関数x [ n ] の片側Z変換の正確な定義です。 X ( z ) = ∑ n = 0 ∞ x [ n ] z − n {\displaystyle X(z)=\sum _{n=0}^{\infty }x[n]z^{-n}}
最後の2つの式を比較すると、片側Z変換とサンプリングされた信号のラプラス変換の関係が分かります。 X q ( s ) = X ( z ) | z = e s T . {\displaystyle X_{q}(s)=X(z){\Big |}_{z=e^{sT}}.}
Z 変換とラプラス変換の類似性は、時間スケール計算 の理論でさらに詳しく説明されています。
ボレル変換 の積分形は、指数型整関数 f に対するラプラス変換の特殊なケースであり、これは、 ある定数A とB に対して成り立つことを意味します。一般化ボレル変換では、指数型ではない関数を変換する際に、指数関数ではなく異なる重み関数を使用することができます。ナクビンの定理は 、ボレル変換が適切に定義されるための必要十分条件を示しています。 F ( s ) = ∫ 0 ∞ f ( z ) e − s z d z {\displaystyle F(s)=\int _{0}^{\infty }f(z)e^{-sz}\,dz} | f ( z ) | ≤ A e B | 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 { a f ( t ) } = a L { f ( t ) } {\displaystyle {\mathcal {L}}\{af(t)\}=a{\mathcal {L}}\{f(t)\}} この線形性と、さまざまな三角関数 、双曲線関数 、複素数 (など) の特性や恒等式を使用すると、定義を直接使用するよりも速く、いくつかのラプラス変換を他のラプラス変換から取得できます。
片側ラプラス変換は、時間領域が非負の 実数である関数を入力として受け取ります。そのため、以下の表の時間領域関数はすべてヘビサイドステップ関数 u ( t ) の倍数になります。
表中の時間遅延τ を伴うエントリは、因果的 (つまりτ > 0 )であることが要求される。因果システムとは、 t = 0 より前のすべての時刻tにおいて インパルス応答 h ( t ) がゼロとなるシステムである。一般に、因果システムの収束領域は反因果システム の収束領域と同じではない。
s 領域等価回路とインピーダンスラプラス変換は回路解析 においてよく用いられ、回路要素のs 領域への簡単な変換が可能です。回路要素は、位相 インピーダンスに非常によく似たインピーダンス に変換できます。
同等のものを以下にまとめます。
s 領域等価回路抵抗器は時間領域とs 領域で全く同じであることに注意してください。回路要素に初期条件がある場合、ソースが挿入されます。例えば、コンデンサに初期電圧がかかる場合や、インダクタに初期電流が流れる場合、s 領域に挿入されたソースはそれを考慮します。
電流源と電圧源の等価物は、上記の表の変換から簡単に導き出されます。
例と応用 ラプラス変換は工学 や物理学で頻繁に用いられます。 線形時不変システム の出力は、その単位インパルス応答を入力信号と畳み込むことで計算できます。この計算をラプラス空間で実行すると、畳み込みは乗算に変換されます。乗算は代数形式であるため、より簡単に解くことができます。詳細については、制御理論 を参照してください。ラプラス変換は、多くの関数クラスで可逆です。システムへの入力または出力の簡単な数学的または機能的記述が与えられれば、 ラプラス変換は代替の機能的記述を提供し、多くの場合、システムの動作を解析するプロセスや、一連の仕様に基づいて新しいシステムを合成するプロセスを簡素化します。[ 42 ]
ラプラス変換は微分方程式を解くのにも利用でき、機械工学 や電気工学 で広く用いられています。ラプラス変換は線形微分方程式を代数方程式に変換し、代数の正式な規則を用いて解くことができます。そして、逆ラプラス変換を適用することで、元の微分方程式を解くことができます。イギリスの電気技師オリバー・ヘヴィサイドは 、ラプラス変換を用いない同様の手法を初めて提案しました。この結果得られた演算法は、ヘヴィサイド計算と呼ばれています。
不定積分の評価 L { f ( t ) } = F ( s ) {\displaystyle {\mathcal {L}}\left\{f(t)\right\}=F(s)} とします。すると(上の表を参照) ∂ s L { f ( t ) t } = ∂ s ∫ 0 ∞ f ( t ) t e − s t d t = − ∫ 0 ∞ f ( t ) e − s t d t = − 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 } = ∫ s ∞ F ( p ) d p . {\displaystyle {\mathcal {L}}\left\{{\frac {f(t)}{t}}\right\}=\int _{s}^{\infty }F(p)\,dp.}
極限 s → 0 {\displaystyle s\rightarrow 0} においては、 極限の交換が正当化され得ることが示される。これは終値定理 の結果としてしばしば可能である。交換が正当化されない場合でも、計算は示唆的となることがある。例えば、a ≠ 0 ≠ b の場合、形式的に進めると、 ∫ 0 ∞ f ( t ) t d t = ∫ 0 ∞ F ( p ) d p , {\displaystyle \int _{0}^{\infty }{\frac {f(t)}{t}}\,dt=\int _{0}^{\infty }F(p)\,dp,} ∫ 0 ∞ cos ( a t ) − cos ( b t ) t d t = ∫ 0 ∞ ( p p 2 + a 2 − p p 2 + b 2 ) d p = [ 1 2 ln p 2 + a 2 p 2 + b 2 ] 0 ∞ = 1 2 ln b 2 a 2 = ln | b a | . {\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 = C d v d t , {\displaystyle i=C{dv \over dt},}
この方程式をラプラス変換すると、 次 の式が得られる 。 I ( s ) = C ( s V ( s ) − V 0 ) , {\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}}} V 0 = v ( 0 ) . {\displaystyle V_{0}=v(0).}
V ( s ) を解くと、 V ( s ) = I ( s ) s C + V 0 s . {\displaystyle V(s)={I(s) \over sC}+{V_{0} \over s}.}
複素インピーダンスZ (オーム )の定義は、初期状態V 0を ゼロに 保持しながら、複素電圧Vを複素電流 I で割った比です。Z ( s ) = V ( s ) I ( s ) | V 0 = 0 . {\displaystyle Z(s)=\left.{V(s) \over I(s)}\right|_{V_{0}=0}.}
この定義と前の式を用いると、次式が得られます。 これはコンデンサの複素インピーダンスの正しい表現です。さらに、ラプラス変換は制御理論において広く応用されています。 Z ( s ) = 1 s C , {\displaystyle Z(s)={\frac {1}{sC}},}
インパルス応答 伝達関数 を持つ線形時間不変システムを考えるH ( s ) = 1 ( s + α ) ( s + β ) . {\displaystyle H(s)={\frac {1}{(s+\alpha )(s+\beta )}}.}
インパルス応答 は、この伝達関数の逆ラプラス変換です。 h ( t ) = L − 1 { H ( s ) } . {\displaystyle h(t)={\mathcal {L}}^{-1}\{H(s)\}.}
部分分数展開 この逆変換を評価するために、まず部分分数展開法を用いて H ( s ) を展開する。1 ( s + α ) ( s + β ) = P s + α + R s + β . {\displaystyle {\frac {1}{(s+\alpha )(s+\beta )}}={P \over s+\alpha }+{R \over s+\beta }.}
未知定数P とR は、伝達関数の対応する極に位置する留数 です。各留数は、その特異点 が伝達関数全体の形状にどの程度寄与するかを表します。
留数定理 によれば、逆ラプラス変換は極とその留数のみに依存します。留数P を求めるには、方程式の両辺にs + α を掛けて、 1 s + β = P + R ( s + α ) s + β . {\displaystyle {\frac {1}{s+\beta }}=P+{R(s+\alpha ) \over s+\beta }.}
そして、 s = − α とすると、 R からの寄与は消え、残るのは P = 1 s + β | s = − α = 1 β − α . {\displaystyle P=\left.{1 \over s+\beta }\right|_{s=-\alpha }={1 \over \beta -\alpha }.}
同様に、留数R は次のように与えられる。 R = 1 s + α | s = − β = 1 α − β . {\displaystyle R=\left.{1 \over s+\alpha }\right|_{s=-\beta }={1 \over \alpha -\beta }.}
であることに注意し 、 R とPを展開した H ( s ) の式に 代入すると、 R = − 1 β − α = − P {\displaystyle R={-1 \over \beta -\alpha }=-P} H ( s ) = ( 1 β − α ) ⋅ ( 1 s + α − 1 s + β ) . {\displaystyle H(s)=\left({\frac {1}{\beta -\alpha }}\right)\cdot \left({1 \over s+\alpha }-{1 \over s+\beta }\right).}
最後に、線形性特性と指数関数的減少の既知の変換(上記のラプラス変換表 の項目 #3を参照)を使用して、 H ( s ) の逆ラプラス変換を実行して、 システムのインパルス応答を 取得 できます。h ( t ) = L − 1 { H ( s ) } = 1 β − α ( e − α t − e − β 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 + β ) = 1 s + α ⋅ 1 s + β {\displaystyle H(s)={\frac {1}{(s+\alpha )(s+\beta )}}={\frac {1}{s+\alpha }}\cdot {\frac {1}{s+\beta }}} L − 1 { 1 s + α } ∗ L − 1 { 1 s + β } = e − α t ∗ e − β t = ∫ 0 t e − α x e − β ( t − x ) d x = e − α t − e − β 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 )}} s sin ( φ ) + ω cos ( φ ) s 2 + ω 2 {\displaystyle {\frac {s\sin(\varphi )+\omega \cos(\varphi )}{s^{2}+\omega ^{2}}}} cos ( ω t + φ ) {\displaystyle \cos {(\omega t+\varphi )}} s cos ( φ ) − ω sin ( φ ) s 2 + ω 2 . {\displaystyle {\frac {s\cos(\varphi )-\omega \sin(\varphi )}{s^{2}+\omega ^{2}}}.}
ラプラス変換から始めて、 分数の項を並べ替えることで逆変換を求めます。 X ( s ) = s sin ( φ ) + ω cos ( φ ) s 2 + ω 2 {\displaystyle X(s)={\frac {s\sin(\varphi )+\omega \cos(\varphi )}{s^{2}+\omega ^{2}}}} X ( s ) = s sin ( φ ) s 2 + ω 2 + ω cos ( φ ) s 2 + ω 2 = sin ( φ ) ( s s 2 + ω 2 ) + cos ( φ ) ( ω s 2 + ω 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 ( φ ) L − 1 { s s 2 + ω 2 } + cos ( φ ) L − 1 { ω s 2 + ω 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 ).}
同様の論理を適用すれば、 L − 1 { s cos φ − ω sin φ s 2 + ω 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 ( β ) = ∫ 0 ∞ e − β E g ( E ) d E {\displaystyle Z(\beta )=\int _{0}^{\infty }e^{-\beta E}g(E)\,dE} g ( E ) = 1 2 π i ∫ β 0 − i ∞ β 0 + i ∞ e β E Z ( β ) 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 = 1 − p {\displaystyle p,q=1-p} λ {\displaystyle \lambda } n {\displaystyle n} t {\displaystyle t} P n ( t ) = ∫ 0 t λ e − λ ( t − s ) ( p P n − 1 ( s ) + q P n + 1 ( s ) ) d s ( + e − λ t when 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 ( P n ) ( s ) , {\displaystyle \pi _{n}(s)={\mathcal {L}}(P_{n})(s),} π n ( s ) = λ λ + s ( p π n − 1 ( s ) + q π n + 1 ( s ) ) ( + 1 λ + s when 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 ) = ∫ 0 ∞ e − s t d μ ( 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)).} lim s → 0 + L μ ( s ) = lim t → ∞ M ( t ) . {\displaystyle \lim _{s\to 0^{+}}{\mathcal {L}}\mu (s)=\lim _{t\to \infty }M(t).} s → 0 + {\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 ) = ∫ 0 ∞ A ( x ) e − x s d x {\displaystyle f(s)=\int _{0}^{\infty }A(x)e^{-xs}\,dx} f ( s ) − c s − 1 {\displaystyle f(s)-{\frac {c}{s-1}}}
この主張は特にリーマンゼータ関数 の対数微分 に適用することができ、素数定理 を証明するための非常に簡潔な方法を提供する。[ 46 ]
参照
注記 ^ Lynn, Paul A. (1986)、「ラプラス変換とz 変換」、Electronic Signals and Systems 、ロンドン:Macmillan Education UK、pp. 225– 272、doi :10.1007/978-1-349-18461-3_6 、ISBN 978-0-333-39164-8 ラプラス変換とZ変換はフーリエ変換と密接に関連しています。ラプラス変換はフーリエ変換よりも範囲が広く、自動制御システムを含む連続回路やシステムを記述するためにエンジニアに広く使用されています。 ^ 「微分方程式 - ラプラス変換」 、 Pauls Online Math Notes 、 2020年8月8日 閲覧 ^ a b Weisstein, Eric W.、 「ラプラス変換」 、 Wolfram MathWorld 、 2020年8月8日 閲覧。 ^ ディーキン、マイケル・AB(1981)、 「ラプラス変換の発展、1737-1937:I.オイラーからスピッツァーまで、1737-1880」 、 正確科学史アーカイブ 、 25 (4): 343-390 、 ISSN 0003-9519 ^ "Des Fonctions génératrices" [生成関数について]、 Théorie Analytique des Probabilités [ 分析的確率論 ] (フランス語) (第 2 版)、パリ、1814 年、第 I 章 2-20 章 ^ Jaynes, ET (Edwin T.) (2003), 確率論:科学の論理 、Bretthorst, G. Larry、ケンブリッジ、イギリス:ケンブリッジ大学出版局、 ISBN 0511065892 、OCLC 57254076 ^ アベル、ニールス H. (1820)、「Sur les fonctions génératrices et leurs déterminantes」、 āuvres Complètes (フランス語)、vol. II (1839 年出版)、 77 ~ 88 ページ 1881年版 ^ オイラー 1744 、オイラー 1753 、オイラー 1769 ^ ラグランジュ 1773 ^ グラッタン・ギネス 1997、260 ページ^ グラッタン・ギネス 1997、261 ページ^ グラッタン・ギネス 1997年 、261~262ページ^ グラッタン・ギネス 1997年 、262~266ページ^ ヘヴィサイド、オリバー (2008年1月) 「 微分変換による定積分の解法」 電磁気理論 第3巻、ロンドン、セクション526、 ISBN 9781605206189 {{citation }}: CS1 maint: location missing publisher (link )^ Edwards, HM (1974), Riemann's Zeta Function , New York: Academic Press, ISBN 0-12-232750-0 、Zbl 0315.10035 ^ ガードナー、マレー F.; バーンズ、ジョン L. (1942)、 「ラプラス変換による線形システムの過渡現象の研究」 、ニューヨーク:ワイリー 付録C^ Lerch、Mathias (1903)、「Sur un point de la théorie des fonctions génératrices d'Abel」[逆転公式の証明]、 Acta Mathematica (フランス語)、 27 : 339–351 、 doi : 10.1007/BF02421315 、 hdl : 10338.dmlcz/501554 ^ Bromwich, Thomas J. (1916)、 「動的システムにおける正規座標」 、 ロンドン数学会紀要 、 15 : 401–448 、 doi : 10.1112/plms/s2-15.1.401 ^ 影響力のある本としては、ガードナー、マレー・F.、バーンズ、ジョン・L.(1942年)、 ラプラス変換による線形システムの過渡現象の研究 、ニューヨーク:ワイリーなどがある。 ^ Doetsch、Gustav (1937)、 Theorie und Anwendung der Laplacesche Transformation [ ラプラス変換の理論と応用 ] (ドイツ語)、ベルリン: Springer 1943年の翻訳^ フェラー 1971 、§XIII.1。^ 累積分布関数は確率密度関数の積分です。 ^ ミクシンスキ、ヤン(2014年7月14日)、 オペレーショナル・カルキュラス 、エルゼビア、 ISBN 9781483278933 ^ ウィダー 1941 、第 II 章、§1^ ウィダー 1941 、第 VI 章、§2^ Beffa, Federico (2024)、 「分布のラプラス変換」 、 弱非線形システム 、Cham: Springer Nature Switzerland、pp. 75– 85、 doi : 10.1007/978-3-031-40681-2_5 、 ISBN 978-3-031-40680-5 、2026年1月14日取得 {{citation }}: CS1 maint: work parameter with ISBN (link )^ コーン&コーン 1967 、226–227 ページ^ Bracewell 2000 、表14.1、p.385^ Ghostarchive と Wayback Machine にアーカイブ: Mattuck, Arthur(2008年11月7日) 「ラプラス変換の起源」 、 YouTube ^ フェラー 1971、432 ページ^ ローラン・シュワルツ (1966)、 『物理科学のための数学』 、アディソン・ウェズリー 、224ページ。^ Titchmarsh, E. (1986) [1948], フーリエ積分理論入門 (第2版)、 Clarendon Press 、p. 6、 ISBN 978-0-8284-0324-5 ^ タカチ 1953、93 ページ^ ライリー、KF; ホブソン、MP; ベンス、SJ (2010)、 物理学と工学のための数学的手法 (第3版)、ケンブリッジ大学出版局、p. 455、 ISBN 978-0-521-86153-3 ^ Distefano, JJ; Stubberud, AR; Williams, IJ (1995)、 「フィードバックシステムと制御」 、Schaum's outlines(第2版)、McGraw-Hill、p. 78、 ISBN 978-0-07-017052-0 ^ 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 の場合を示します。^ http://mathworld.wolfram.com/LaplaceTransform.html – Wolfram Mathword は複素数qのケースを提供します ^ コーン&コーン 1967 、§8.1^ RK Pathria; Paul Beal (1996)、 統計力学 (第2版)、Butterworth-Heinemann、p. 56 、 ISBN 9780750624695 ^ Salem, M.; Seaton, MJ (1974)、「I. 連続スペクトルと輝度等高線」、 Monthly Notices of the Royal Astronomical Society 、 167 : 493–510 、 Bibcode : 1974MNRAS.167..493S 、 doi : 10.1093/mnras/167.3.493 、Salem, M. (1974)、「II. 3次元モデル」、 Monthly Notices of the Royal Astronomical Society 、 167 : 511–516 、 Bibcode : 1974MNRAS.167..511S 、 doi : 10.1093/mnras/167.3.511 ^ a b フェラー 『確率論入門』第2巻、pp=479-483 ^ 池原 誠 (1931)、「ランダウの定理の解析的数論への拡張」、 数学物理学ジャーナル 、 10 ( 1–4 ): 1– 12、 doi : 10.1002/sapm19311011 、 Zbl 0001.12902
参考文献
モダンな ブレイスウェル、ロナルド・N.(1978)『フーリエ変換とその応用 (第2版)』マグロウヒル工学社、ISBN 978-0-07-007013-4 Bracewell, RN (2000) 『フーリエ変換とその応用』 (第3版)、ボストン:McGraw-Hill、ISBN 978-0-07-116043-8 フェラー、ウィリアム (1971)『確率論とその応用入門』第2巻 、第2版、ニューヨーク:ジョン・ワイリー・アンド・サンズ 、MR 0270403 Korn, GA; Korn, TM (1967)、『科学者と技術者のための数学ハンドブック (第2版)』、McGraw-Hill Companies、ISBN 978-0-07-035370-1 Widder, David Vernon (1941), The Laplace Transform , Princeton Mathematical Series, v. 6, Princeton University Press , MR 0005923 ウィリアムズ、J.(1973)、ラプラス変換 、問題解決者、ジョージ・アレン&アンウィン、ISBN 978-0-04-512021-5 Takacs, J. (1953)、「Fourier amplitudok meghatarozasa Operatorszamitassal」、Magyar Hiradastechnika (ハンガリー語)、IV ( 7– 8): 93– 96
歴史的 オイラー、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 ページ
さらに読む
外部リンク