Heaviside 演算子法えんざんしほうと変換へんかん

date2026-07-14document_iddoc_35c6837c7649a7c0326f70a219f1fe6ddescriptionHeavisideの演算子法を、微分作用素の多項式・逆作用素・Laplace/Fourier/Z変換との関係から整理する講義である。prerequisites線型作用素方程式 / 多項式作用素と固有値問題 / 作用素カタログ-基礎編 / ラプラス変換の入口type講義content_typelecturestatusactive

mathlinear-operatorheaviside-operator-methodlaplacefourierz-transformlecture

1導入どうにゅう

この講義こうぎでは、Heaviside 演算子法えんざんしほうを「記号きごうを分数ぶんすうのように扱あつかう裏技うらわざ」ではなく、線型作用素方程式せんけいさようそほうていしきlinear operator equation

P (D) y = f

の特殊解とくしゅかいを構成こうせいするための作用素的さようそてきな記法きほうとして読よむことである。

ここで

D = \frac{d}{d t}

は微分作用素びぶんさようそであり、 $P (D)$ は $D$ の多項式作用素たこうしきさようそである。

data/lecture/math/linear-operator/polynomial-operators-and-eigenvalue-problems.lecture.n.md

2なぜ $D$ を記号きごうとして扱あつかうのか

定数係数ていすうけいすう線型せんけい微分方程式びぶんほうていしき

a y^{'}' + b y^{'} + c y = f (t)

は

(a D^{2} + b D + c I) y = f (t)

と書かける。 $P (D) = a D^{2} + b D + c I$ と置おけば

P (D) y = f

である。

この形かたちにする理由りゆうは、微分びぶんを何回なんかいも展開てんかいするのではなく、「未知関数みちかんすう $y$ に作用さようする 1 つの線型作用素せんけいさようそ」として見みたいからである。すると一般論いっぱんろんにより、同次解どうじかいは $\ker P (D)$ 、特殊解とくしゅかいは $P (D) y = f$ を満みたす 1 本ほんの解かいとして整理せいりできる。

3形式的けいしきてきな逆作用素ぎゃくさようそ

Heaviside 演算子法えんざんしほうでは、形式的けいしきてきに

y = \frac{1}{P (D)} f

のように書かくことがある。これは「 $P (D)$ の逆作用素ぎゃくさようそを $f$ に作用さようさせる」という意味いみであり、普通ふつうの数かずの割わり算ざんではない。

注意ちゅういすべき点てんは、 $P (D)$ が核かくを持もつ場合ばあい、逆作用素ぎゃくさようそは全体ぜんたいでは一意いちいに定義ていぎできないことである。したがって、この記法きほうは特殊解とくしゅかいを 1 つ作つくる道具どうぐとして読よむのが安全あんぜんである。一般解いっぱんかいは、その後あとに $\ker P (D)$ を足たして得える。

4命題めいだい：指数入力しすうにゅうりょくでは $P (D)$ が $P (λ)$ になる

4.1主張しゅちょう

$P$ を多項式たこうしきとする。このとき

P (D) e^{λ t} = P (λ) e^{λ t}

である。

4.2証明しょうめい

$D (e^{λ t}) = λ e^{λ t}$ である。したがって帰納法きのうほうにより

D^{k} e^{λ t} = λ^{k} e^{λ t}

である。 $P (s) = a_{0} + a_{1} s + \dots + a_{m} s^{m}$ と書かけば

\begin{aligned} P(D)e^{\lambda t} &=(a_0I+a_1D+\cdots+a_mD^m)e^{\lambda t}\\ &=(a_0+a_1\lambda+\cdots+a_m\lambda^m)e^{\lambda t}\\ &=P(\lambda)e^{\lambda t} \end{aligned}\begin{aligned} P(D)e^{\lambda t} &=(a_0I+a_1D+\cdots+a_mD^m)e^{\lambda t}\\ &=(a_0+a_1\lambda+\cdots+a_m\lambda^m)e^{\lambda t}\\ &=P(\lambda)e^{\lambda t} \end{aligned}

である。□

5特殊解とくしゅかいの作つくり方かた

もし $P (λ) \neq 0$ なら、

P (D) (\frac{e^{λ t}}{P (λ)}) = e^{λ t}

である。したがって

P (D) y = e^{λ t}

の特殊解とくしゅかいとして

y_{p} = \frac{e^{λ t}}{P (λ)}

を取とれる。

右辺うへんが $\cos ω t$ や $\sin ω t$ の場合ばあいは、 $e^{i ω t}$ を使つかって複素数ふくそすうで計算けいさんし、最後さいごに実部じつぶや虚部きょぶを取とる。

6共鳴きょうめいする場合ばあい

$P (λ) = 0$ の場合ばあい、 $e^{λ t}$ は同次方程式どうじほうていしき $P (D) y = 0$ の解かいである。したがって

\frac{e^{λ t}}{P (λ)}

という式しきは使つかえない。未定係数法みていけいすうほうで $t$ を掛かける補正ほせいが必要ひつようになるのは、この重複ちょうふくを避さけるためである。

言いい換かえると、 $e^{λ t}$ をそのまま特殊解候補とくしゅかいこうほにしても、 $P (λ) = 0$ なら

P (D) e^{λ t} = 0

となり、右辺うへんの $e^{λ t}$ を作つくれない。これは候補こうほが同次解空間どうじかいくうかんに入はいっている、という意味いみである。

$λ$ が $P$ の $m$ 重根じゅうこんなら、候補こうほを

t^{m} e^{λ t}

の形かたちまで広ひろげる。この $t^{m}$ は暗記あんきする飾かざりではない。同次解どうじかいと重複ちょうふくする候補こうほから離はなれ、 $P (D)$ を作用さようさせたときに右辺うへんを作つくれる方向ほうこうへ移うつるための補正ほせいである。

7Laplace 変換へんかんとの関係かんけい

Heaviside 演算子法えんざんしほうの形式計算けいしきけいさんは、Laplace 変換へんかんで正当化せいとうかしやすい。初期値しょきちが 0 の場合ばあい、

[PARSE ERROR: Undefined("Command(\"mathcal\")")] L {D y} = s Y (s)

に近ちかい形かたちになるので、

[PARSE ERROR: Undefined("Command(\"mathcal\")")] L {P (D) y} = P (s) Y (s)

となる。したがって

P (D) y = f

は

P (s) Y (s) = F (s)

へ変かわり、

Y (s) = \frac{F (s)}{P (s)}

と解とける。

初期値しょきちが 0 でない場合ばあいは、 $[PARSE ERROR: Undefined("Command(\"mathcal\")")] L {y^{'}} = s Y (s) - y (0)$ のように初期値項しょきちこうが現あらわれる。したがって、Heaviside 演算子法えんざんしほうだけで初期値しょきちを無視むししてよいわけではない。Laplace 変換へんかんは、この初期値しょきちの扱あつかいを明示めいじする。

data/lecture/math/analysis/introduction-to-laplace-transform.lecture.n.md

8変換へんかんの対応表たいおうひょう

変換へんかん	対象たいしょう	中心作用素ちゅうしんさようそ	変換後へんかんご	向むく問題もんだい
Laplace transform	連続時間れんぞくじかん $t [PARSE ERROR: Undefined("Command(\"ge\")")] 0$	微分作用素びぶんさようそ $D$	$D \mapsto s$ と初期値項しょきちこう	初期値問題しょきちもんだい、ステップ入力にゅうりょく、デルタ入力にゅうりょく
Fourier transform	全時間ぜんじかん・周波数解析しゅうはすうかいせき	微分作用素びぶんさようそ $D$	$D \mapsto i ω$	波なみ、振動しんどう、PDE、畳たたみ込こみ
Z-transform	離散時間列りさんじかんれつ $x [n]$	シフト作用素さようそ $E$ 、遅延作用素ちえんさようそ $R$	両側りょうがわで $E \mapsto z$ 、 $R \mapsto z^{- 1}$	漸化式ぜんかしき、差分方程式さぶんほうていしき、デジタルフィルタ

Fourier 変換へんかんが効きく理由りゆうは、 $e^{i ω t}$ が微分作用素びぶんさようその固有関数こゆうかんすうだからである。Z 変換へんかんが効きく理由りゆうは、 $z^{- n}$ の重おもみが添字そえじのシフトを $z$ や $z^{- 1}$ の掛かけ算ざんへ変かえるからである。

data/lecture/math/analysis/introduction-to-fourier-transform.lecture.n.md data/lecture/math/analysis/introduction-to-z-transform.lecture.n.md

9畳たたみ込こみと逆作用素ぎゃくさようそ

線型時不変系せんけいじふへんけいでは、入力にゅうりょく $f$ に対たいする出力しゅつりょくが

y = h * f

と表あらわされることがある。ここで $h$ はインパルス応答おうとうである。変換へんかんを使つかうと

Y = H F

となるので、畳たたみ込こみ作用素さようそは掛かけ算ざんに変かわる。

この見方みかたでは、 $[PARSE ERROR: Undefined("Command(\"dfrac\")")] 1 P (D)$ は「 $P (D)$ の逆ぎゃくに対応たいおうするインパルス応答おうとうと畳たたみ込こむ作用素さようそ」として理解りかいできる。つまり、形式的けいしきてきな逆作用素ぎゃくさようそは、変換へんかんや畳たたみ込こみを通つうじて具体化ぐたいかされる。

data/lecture/math/differential-equations/step-functions-delta-functions-and-convolution.lecture.n.md

10どこまで成立せいりつするか

Heaviside 演算子法えんざんしほうは、定数係数ていすうけいすう線型せんけい微分方程式びぶんほうていしきでは強力きょうりょくである。しかし、係数けいすうが変数へんすうに依存いぞんすると、 $D$ と乗法作用素じょうほうさようそが一般いっぱんに可換かかんでないため、多項式たこうしきのように扱あつかえない。

また、 $P (D)^{- 1}$ は核かくがあると一意いちいではない。境界条件きょうかいじょうけんや初期条件しょきじょうけんを含ふくめると、どの逆ぎゃくを選えらぶかが問題もんだいになる。したがって、形式計算けいしきけいさんの後あとには、元もとの微分方程式びぶんほうていしきと条件じょうけんへ代入だいにゅうして検算けんざんする。

11一言ひとことでいうと

Heaviside 演算子法えんざんしほうは、 $P (D) y = f$ を作用素方程式さようそほうていしきとして見みて、形式的けいしきてきな逆作用素ぎゃくさようそで特殊解とくしゅかいを作つくる方法ほうほうである。Laplace・Fourier・Z 変換へんかんは、この作用素さようそを代数的だいすうてきな掛かけ算ざんへ移うつす座標変換ざひょうへんかんとして働はたらく。