制御工学や電気回路の解析を学び始めると、必ずと言っていいほど登場するのがラプラス変換です。
その中でも「微分のラプラス変換」は、微分方程式を代数方程式に変換するための要となる考え方です。
L{f'(t)}という記号を見て、何を意味しているのか分からず戸惑った方も多いのではないでしょうか。
本記事では、微分のラプラス変換の公式とその性質について、できるだけ丁寧に解説していきます。
s領域という概念や、微分定理の証明、さらには微分方程式への応用まで幅広く取り上げます。
数式だけでなく具体例も交えながら進めていきますので、初めて学ぶ方でも理解しやすい内容になっているはずです。
それではまず、微分のラプラス変換の結論となる公式から見ていきましょう。
微分のラプラス変換の結論、L{f'(t)}=sF(s)−f(0)という公式
それではまず、微分のラプラス変換に対する結論について解説していきます。
結論から言うと、関数f(t)の一階微分f'(t)のラプラス変換は、次の式で表されます。
L{f'(t)} = sF(s) − f(0)
ここでF(s)はf(t)のラプラス変換であり、sはラプラス変換における複素変数です。
f(0)は時間tがゼロのときの初期値を意味しています。
微分という時間領域の演算が、s領域では単純にsを掛け算するだけの操作に置き換わるという点が、この公式の最大のポイントでしょう。
この性質のおかげで、微分方程式を解く際に厄介な微分演算を回避できるのです。
公式が成り立つ理由の概要
なぜsを掛けるだけで微分が表現できるのか、疑問に思う方も多いはずです。
これはラプラス変換の定義式に対して部分積分を適用することで導かれます。
部分積分の過程でf(0)という初期値の項が自然に現れてくる点が特徴的です。
詳しい証明については後の章で改めて取り上げます。
高階微分への拡張公式
二階微分や三階微分についても、同様の考え方で公式が成立します。
L{f”(t)} = s²F(s) − sf(0) − f'(0)
L{f”'(t)} = s³F(s) − s²f(0) − sf'(0) − f”(0)
次数が上がるごとに、初期値の項が増えていく規則性に気づくでしょう。
n階微分の一般式は以下のようにまとめられます。
L{f⁽ⁿ⁾(t)} = sⁿF(s) − sⁿ⁻¹f(0) − sⁿ⁻²f'(0) − ⋯ − f⁽ⁿ⁻¹⁾(0)
この一般式さえ覚えておけば、どんな階数の微分方程式にも対応できます。
公式が実務で持つ意味
単なる数式の変形に見えるかもしれませんが、この公式には工学的に大きな意味があります。
電気回路のコイルやコンデンサの微分方程式を、初期条件込みでそのまま代数方程式に変換できるのです。
ステップ応答やインパルス応答の解析でも、この性質が土台になっています。
初期値を式の中に自然に組み込めるという点が、他の変換手法にはない強みでしょう。
ラプラス変換の基礎をおさらい
続いては、微分のラプラス変換を理解するための前提知識を確認していきます。
基礎が曖昧なままだと、公式を暗記するだけで終わってしまいがちです。
ここでしっかりと土台を固めておきましょう。
ラプラス変換の定義式
ラプラス変換とは、時間関数f(t)を複素変数sの関数F(s)に変換する操作のことです。
F(s) = L{f(t)} = ∫₀^∞ f(t)e⁻ˢᵗ dt
この積分はtがゼロから無限大までの範囲で計算されます。
指数関数e⁻ˢᵗを掛けて積分することで、時間領域の情報がs領域に集約されるイメージです。
制御工学ではこのs領域での計算のしやすさが重宝されています。
s領域とは何か
s領域とは、時間tを変数とする時間領域に対して、複素変数sを変数とする世界を指します。
時間領域では微分方程式として表現される問題が、s領域では代数方程式として扱えるようになります。
これはちょうど、時間軸での複雑な波形を、周波数のような別の視点から見直す感覚に近いでしょう。
sは一般にσ+jωという複素数で表され、実部と虚部それぞれが減衰と振動の性質を担っています。
この抽象的な変換こそが、微分方程式を簡単に解くための鍵になっているのです。
代表的な関数のラプラス変換表
実際の計算では、よく使う関数の変換結果を一覧として持っておくと非常に便利です。
| 時間関数 f(t) | ラプラス変換 F(s) |
|---|---|
| 1(単位ステップ) | 1/s |
| t | 1/s² |
| tⁿ | n!/sⁿ⁺¹ |
| e^(at) | 1/(s−a) |
| sin(ωt) | ω/(s²+ω²) |
| cos(ωt) | s/(s²+ω²) |
この表を頭に入れておくと、微分定理と組み合わせた計算がぐっとスムーズになります。
特にe^(at)やsin、cosは微分方程式の解として頻出するため、覚えておいて損はないでしょう。
微分のラプラス変換の証明
続いては、なぜあの公式が成り立つのか、その証明過程を確認していきます。
証明を理解しておくと、公式を丸暗記するよりも応用が効きやすくなります。
部分積分による導出過程
まずは定義式にf'(t)を代入するところから始めます。
L{f'(t)} = ∫₀^∞ f'(t)e⁻ˢᵗ dt
この積分に対して部分積分を適用します。
u=e⁻ˢᵗ、dv=f'(t)dtと置くと、部分積分の公式∫u dv = [uv] − ∫v duが使えます。
∫₀^∞ f'(t)e⁻ˢᵗ dt = [f(t)e⁻ˢᵗ]₀^∞ + s∫₀^∞ f(t)e⁻ˢᵗ dt
右辺の第一項は、tが無限大でe⁻ˢᵗがゼロに収束すると仮定すると、−f(0)だけが残ります。
右辺の第二項は、まさにsF(s)そのものです。
これらをまとめると、L{f'(t)} = sF(s) − f(0)という公式が導かれます。
部分積分という基本的な積分技法だけで、この重要な公式が得られるというのは興味深い点でしょう。
収束のための条件
先ほどの証明では、tが無限大に近づくときにf(t)e⁻ˢᵗがゼロに収束することを前提にしていました。
この条件が成り立たない関数に対しては、公式をそのまま適用できません。
一般に、f(t)が指数位数の関数であれば、sの実部を十分大きく取ることで収束条件を満たせます。
指数位数とは、|f(t)|がMe^(at)以下に抑えられるという性質のことです。
この条件は工学で扱う多くの信号に対して自然に満たされるため、過度に心配する必要はないでしょう。
高階微分公式の証明の流れ
二階微分の公式は、一階微分の公式を二回適用することで導けます。
f”(t)を(f'(t))’と考え、まずg(t)=f'(t)としてL{g'(t)} = sL{g(t)} − g(0)を適用します。
そこにL{f'(t)} = sF(s) − f(0)を代入すると、L{f”(t)} = s²F(s) − sf(0) − f'(0)という式が得られます。
この操作を繰り返すことで、n階微分に対する一般公式も帰納的に証明できるのです。
数学的帰納法の考え方を使えば、どんな階数でも同じロジックで説明がつくでしょう。
微分方程式への応用
続いては、微分のラプラス変換を実際の微分方程式にどう使うのかを確認していきます。
公式だけを知っていても、使い方が分からなければ意味がありません。
常微分方程式を解く基本手順
ラプラス変換を用いて微分方程式を解く手順は、大きく三つのステップに分かれます。
| 手順 | 内容 |
|---|---|
| ステップ1 | 微分方程式の両辺をラプラス変換し、代数方程式に変換する |
| ステップ2 | F(s)について代数的に解く |
| ステップ3 | 逆ラプラス変換を行い、時間領域の解f(t)に戻す |
この三段階のプロセスさえ押さえておけば、多くの線形微分方程式に対応できます。
特に定数係数の微分方程式であれば、比較的機械的に解を求められるでしょう。
初期値問題との相性の良さ
微分のラプラス変換公式には初期値f(0)が自動的に組み込まれるという点が、初期値問題を解くうえで大きな利点になります。
通常の解法では、一般解を求めたあとに初期条件を代入して積分定数を決定するという手順が必要です。
一方、ラプラス変換を使えば、変換の段階で初期条件がすでに式に反映されています。
そのため、特殊解や一般解を分けて考える必要がなく、一気に初期値問題の解にたどり着けるのです。
この効率の良さこそが、制御工学でラプラス変換が愛用される理由の一つでしょう。
具体例で計算の流れを確認
実際に簡単な例で確認してみましょう。
次の微分方程式を初期条件y(0)=1のもとで解きます。
y'(t) + 2y(t) = 0
両辺をラプラス変換すると、L{y'(t)} + 2L{y(t)} = 0となります。
微分定理を適用すると、sY(s) − y(0) + 2Y(s) = 0という式になります。
y(0)=1を代入すると、sY(s) − 1 + 2Y(s) = 0です。
これを整理すると、Y(s)(s+2) = 1、つまりY(s) = 1/(s+2)が得られます。
先ほどの変換表からL⁻¹{1/(s−a)} = e^(at)であることを思い出すと、y(t) = e⁻²ᵗという解が導かれます。
実際にこの解を元の微分方程式に代入すると、条件を満たしていることが確認できるでしょう。
積分のラプラス変換との違い
続いては、微分のラプラス変換と対をなす積分のラプラス変換との違いを確認していきます。
両者の関係を理解しておくと、s領域での操作がより直感的につかめるようになります。
積分定理の公式
関数f(t)の積分∫₀ᵗ f(τ)dτに対するラプラス変換は、次のように表されます。
L{∫₀ᵗ f(τ)dτ} = F(s)/s
微分がsを掛ける操作であったのに対し、積分はsで割る操作に対応しています。
この対称性は偶然ではなく、微分と積分が互いに逆演算であることの自然な帰結でしょう。
微分と積分の対応関係
時間領域での微分と積分、そしてs領域でのsの掛け算と割り算がきれいに対応しています。
| 時間領域の操作 | s領域での対応 |
|---|---|
| 微分 d/dt | sを掛ける(初期値を引く) |
| 積分 ∫dt | sで割る |
この対応関係のおかげで、微分方程式だけでなく積分方程式もラプラス変換で扱えるのです。
微分と積分がs領域では単なる四則演算に還元されるという事実は、ラプラス変換の最大の魅力と言っても過言ではないでしょう。
使い分けのポイント
実際の問題では、微分方程式なのか積分方程式なのかを見極めることが第一歩になります。
回路解析では、コイルの電圧降下は電流の微分で表され、コンデンサの電圧は電流の積分で表されます。
そのため、L形やRC回路のようにインダクタンスとキャパシタンスが混在する回路では、微分定理と積分定理を両方使う場面が出てくるでしょう。
どちらの定理を使うべきか迷ったときは、まず対象の物理量が時間の何を表しているのかを整理することをおすすめします。
計算時の注意点とよくある間違い
続いては、微分のラプラス変換を計算する際に陥りやすいミスについて確認していきます。
公式自体はシンプルですが、実際に手を動かすと意外なところでつまずくものです。
初期値の代入忘れ
最も多いミスは、f(0)やf'(0)といった初期値の項を書き忘れてしまうことです。
特に高階微分になると項数が増えるため、うっかり一つ落としてしまうケースが目立ちます。
n階微分の公式では初期値の項がn個現れるという点を意識しておくと、抜け漏れを防ぎやすくなるでしょう。
計算のあとに一度、次数と項の数が一致しているか見直す習慣をつけておくと安心です。
収束条件の見落とし
公式を機械的に適用するあまり、収束条件を確認せずに計算を進めてしまう方も少なくありません。
指数関数的に急激に増大する関数の場合、sの範囲によっては変換自体が定義できないこともあります。
実務上は工学的に妥当な信号を扱うことがほとんどなので過度に神経質になる必要はありませんが、理論的な背景として頭の片隅に置いておくとよいでしょう。
逆変換の際に起こりやすいミス
Y(s)を求めたあとの逆ラプラス変換で、部分分数分解を誤ってしまうケースもよく見られます。
分母が二次式や三次式になると、部分分数分解の係数決定でケアレスミスが起きやすくなります。
逆変換前には必ず分母を因数分解し、変換表と照らし合わせながら一つずつ丁寧に対応させることが大切です。
焦って計算を進めるよりも、途中式を省略せずに確認しながら進めるほうが結果的に早く正確な答えにたどり着けるでしょう。
まとめ
今回は微分のラプラス変換について、公式とその性質を中心に解説してきました。
結論として、一階微分に対してはL{f'(t)} = sF(s) − f(0)という公式が成り立ちます。
高階微分についても、sのべき乗と初期値の項が規則的に増えていく一般公式が存在します。
この性質のおかげで、時間領域では厄介だった微分方程式を、s領域では代数方程式として簡単に解けるようになるのです。
証明の過程では部分積分が使われており、収束条件など理論的な背景も押さえておくとより理解が深まります。
積分のラプラス変換と対比させることで、微分と積分がs領域でどう対応しているのかもクリアになったのではないでしょうか。
初期値の代入忘れや収束条件の見落とし、逆変換時のミスといった注意点も、事前に知っておけば十分に回避できます。
ぜひ本記事の内容を参考に、微分のラプラス変換を使った微分方程式の解法にチャレンジしてみてください。