markdown
多項式作用素と固有値問題md 4ec75fe
lecture/math/linear-operator/polynomial-operators-and-eigenvalue-problems.lecture.n.md
Download PDF

多項式作用素たこうしきさようそpolynomial operatorと固有値問題こゆうちもんだいeigenvalue problem

mathlinear-operatorpolynomial-operatoreigenvaluelecture

1導入どうにゅう

この講義こうぎでは、特性方程式とくせいほうていしきcharacteristic equationが突然とつぜん現あらわれる公式こうしきではなく、線型作用素せんけいさようそlinear operatorの多項式たこうしきpolynomialを調しらべる問題もんだいから自然しぜんに出でることを確認かくにんする点てんにある。

1 つの線型作用素せんけいさようそlinear operator T:V→V があるとき、多項式たこうしき

p(t)=a0+a1t+⋯+amtm

から

p(T)=a0I+a1T+⋯+amTm

を作つくれる。すると

p(T)x=y

も、前まえの講義こうぎと同おなじ Lx=y 型かたの線型作用素方程式せんけいさようそほうていしきlinear operator equationである。つまり、見みるべき本質ほんしつはやはり

kerp(T),Imp(T)

である。

2用語ようごと定義ていぎ

T:V→V を線型作用素せんけいさようそlinear operatorとする。T0=I、Tk+1=T∘Tk と定義ていぎする。

多項式作用素たこうしきさようそpolynomial operatorとは、多項式たこうしき p(t)=a0+a1t+⋯+amtm に対たいして

p(T)=a0I+a1T+⋯+amTm

で定義ていぎされる線型作用素せんけいさようそlinear operatorである。

固有値こゆうちeigenvalue λ と固有こゆうベクトルeigenvector v は

Tv=λv,v≠0

を満みたす組くみである。

3方針ほうしん

この講義こうぎでは、次つぎの順じゅんで読よむ。

  1. p(T) が再ふたたび線型作用素せんけいさようそlinear operatorであることを証明しょうめいする。
  2. p(T)x=y の可解性かかいせいと一般解いっぱんかいを、核かくkernelと像ぞうimageで読よむ。
  3. Tv=λv なら p(T)v=p(λ)v であることを証明しょうめいする。
  4. 定数ていすう係数けいすうの漸化式ぜんかしきと微分方程式びぶんほうていしきで、なぜ特性とくせい方程式ほうていしきが出でるかを説明せつめいする。

4命題めいだい 1:p(T) は線型作用素せんけいさようそlinear operatorである

4.1主張しゅちょう

T:V→V が線型作用素せんけいさようそlinear operatorなら、任意にんいの多項式たこうしき p(t) に対たいして p(T) も線型作用素せんけいさようそlinear operatorである。

4.2証明しょうめい

まず T が線型せんけいlinearなら、合成ごうせい Tk も線型せんけいlinearである。これは k についての帰納法きのうほうで示しめせる。T0=I は線型せんけいであり、Tk が線型せんけいなら Tk+1=T∘Tk も線型せんけいである。

次つぎに、線型せんけい作用素さようその線型結合せんけいけつごうlinear combinationは線型せんけいである。実際じっさい、S1,…,Sm が線型せんけいで、c1,…,cm がスカラーなら

(∑iciSi)(αx+βy)=∑iciSi(αx+βy)=α∑iciSi(x)+β∑iciSi(y)

である。

したがって

p(T)=a0I+a1T+⋯+amTm

は線型せんけいである。

5命題めいだい 2:p(T)x=y は一般いっぱんの Lx=y と同おなじ

5.1主張しゅちょう

L=p(T) と置おくと、方程式ほうていしき

p(T)x=y

の可解かかい条件じょうけんは

y∈Imp(T)

である。また、解かい xp が 1 つ見みつかれば、解かい全体ぜんたいは

x=xp+h,h∈kerp(T)

である。

5.2証明しょうめい

前まえ講義こうぎの定理ていり 1(可解性かかいせい)と定理ていり 2(一般解いっぱんかい)を L=p(T) に適用てきようするだけである。p(T) は命題めいだい 1 により線型作用素せんけいさようそlinear operatorであるため、Lx=y の一般論いっぱんろんをそのまま使用しようできる。

この点てんが重要じゅうようである。p(T) という記号きごうが複雑ふくざつでも、方程式ほうていしきの型かたは変かわっていない。可解性かかいせいは像ぞうimage、自由度じゆうどは核かくkernelで決きまる。

6命題めいだい 3:固有こゆうベクトルeigenvector上うえで p(T) は p(λ) になる

6.1主張しゅちょう

Tv=λv、v≠0 とする。このとき任意にんいの多項式たこうしき p について

p(T)v=p(λ)v

である。

6.2証明しょうめい

まず Tkv=λkv を k について帰納法きのうほうで示しめす。k=0 では T0v=Iv=v=λ0v である。Tkv=λkv と仮定かていすると

Tk+1v=T(Tkv)=T(λkv)=λkTv=λkλv=λk+1v

である。

したがって

\begin{aligned} p(T)v &=(a_0I+a_1T+\cdots+a_mT^m)v\\ &=a_0v+a_1Tv+\cdots+a_mT^mv\\ &=a_0v+a_1\lambda v+\cdots+a_m\lambda^m v\\ &=p(\lambda)v \end{aligned}

である。

7特性方程式とくせいほうていしきが出でる理由りゆう

命題めいだい 3 から、もし p(λ)=0 なら

p(T)v=p(λ)v=0

である。したがって

v∈kerp(T)

である。つまり、p(T)x=0 の解かいを探さがすとき、T の固有こゆうベクトルeigenvectorを使つかうと、作用素さようそ方程式ほうていしきが数かずの方程式ほうていしき

p(λ)=0

へ落おちる。

これが特性方程式とくせいほうていしきcharacteristic equationの抽象的ちゅうしょうてきな理由りゆうである。特性とくせい方程式ほうていしきは「たまたま試ためす式しき」ではなく、固有値問題こゆうちもんだいeigenvalue problemで作用素さようそをスカラーへ還元かんげんするための式しきである。

8例れい 1:漸化式ぜんかしきrecurrence relationとシフト作用素さようそshift operator

列れつ a=(an) に対たいして、シフト作用素さようそshift operator E を

(Ea)n=an+1

で定義ていぎする。定数ていすう係数けいすう線型せんけい漸化式ぜんかしき

an+2-c1an+1-c2an=0

は

(E2-c1E-c2I)a=0

と書かける。ここで

p(t)=t2-c1t-c2

と置おけば、これは

p(E)a=0

である。

指数しすう型かたの列れつ rn は

E(rn)=rn+1=rrn

を満みたすので、E の固有列こゆうれつeigen-sequenceである。よって

p(E)(rn)=p(r)rn

である。したがって rn が解かいになる条件じょうけんは

p(r)=0

である。これが定数ていすう係数けいすう線型せんけい漸化式ぜんかしきの特性とくせい方程式ほうていしきである。

data/lecture/math/sequence/shift-operators-and-recurrences.lecture.n.md

9例れい 2:微分作用素びぶんさようそdifferential operatorと定数ていすう係数けいすう ODE

滑なめらかな関数かんすうの空間くうかんで、微分作用素びぶんさようそdifferential operator D を

Df=f′

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

ay′'+by′+cy=0

は

(aD2+bD+cI)y=0

と書かける。

指数関数しすうかんすう erx は

D(erx)=rerx

を満みたすので、D の固有関数こゆうかんすうeigenfunctionである。よって

(aD2+bD+cI)erx=(ar2+br+c)erx

である。したがって erx が解かいになる条件じょうけんは

ar2+br+c=0

である。これは二階にかい定数ていすう係数けいすう微分方程式びぶんほうていしきの特性とくせい方程式ほうていしきである。

data/lecture/math/differential-equations/second-order-linear-constant-coefficient-odes.lecture.n.md

10同伴行列どうはんぎょうれつとの関係かんけい

モニック多項式たこうしき

p(t)=tn+cn-1tn-1+⋯+c1t+c0

に対たいして、同伴行列どうはんぎょうれつcompanion matrixは、p(t) を特性多項式とくせいたこうしきcharacteristic polynomialとして持もつように作つくられる行列ぎょうれつである。

同伴どうはん行列ぎょうれつが重要じゅうようなのは、抽象的ちゅうしょうてきな多項式たこうしき p(t) を「ある線型せんけい作用素さようその固有値問題こゆうちもんだい」として実現じつげんする標準ひょうじゅんモデルになるからである。漸化式ぜんかしきでは、状態じょうたいベクトル

(an+k-1an+k-2⋮an)

を 1 歩ほ進すすめる行列ぎょうれつが同伴どうはん行列型ぎょうれつがたになる。すると、漸化式ぜんかしきの特性とくせい方程式ほうていしきと同伴どうはん行列ぎょうれつの固有値問題こゆうちもんだいが一致いっちする。

ここで混同こんどうしてはならないのは、同伴どうはん行列ぎょうれつは「多項式たこうしき作用素さようそ p(T) の一般いっぱん定義ていぎ」ではなく、特定とくていの多項式たこうしきを行列ぎょうれつの特性とくせい多項式たこうしきとして実現じつげんする具体的ぐたいてきな模型もけいだという点てんである。一般論いっぱんろんは p(T) と kerp(T) であり、同伴どうはん行列ぎょうれつはそれを座標ざひょうで実現じつげんする代表例だいひょうれいである。

data/lecture/math/linear-algebra/companion-matrix-basics.lecture.n.md

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

p(T) の定義ていぎは、T が同おなじ空間くうかん V から V へ戻もどる線型作用素せんけいさようそlinear operatorであるときに成立せいりつする。T:V→W で V と W が異ことなる場合ばあい、一般いっぱんには T2 が定義ていぎできないため、多項式たこうしき p(T) はそのまま作つくれない。

また、p(λ)=0 から v∈kerp(T) は従したがうが、kerp(T) のすべてが単純たんじゅんな固有こゆうベクトルだけで張はられるとは限かぎらない。重根じゅうこんや対角化たいかくか不能ふのうな場合ばあいには、一般化いっぱんか固有こゆうベクトルgeneralized eigenvectorや Jordan 構造こうぞうが現あらわれる。この点てんは、一般化いっぱんか前まえの固有こゆうベクトルだけの議論ぎろんと、一般化いっぱんか後のちの一般化固有空間いっぱんかこゆうくうかんgeneralized eigenspaceの議論ぎろんを分わけて扱あつかう必要ひつようがある。

12最終形さいしゅうけい

[PARSE ERROR: Undefined("Command(\"boxed\")")]p(T)=a0I+a1T+⋯+amTm
[PARSE ERROR: Undefined("Command(\"boxed\")")]p(T)x=yissolvable⟺y∈Imp(T)
[PARSE ERROR: Undefined("Command(\"boxed\")")]p(T)xp=y⇒{x∣p(T)x=y}=xp+kerp(T)
[PARSE ERROR: Undefined("Command(\"boxed\")")]Tv=λv⇒p(T)v=p(λ)v
[PARSE ERROR: Undefined("Command(\"boxed\")")]p(λ)=0⇒v∈kerp(T)

13一言ひとことでいうと

特性方程式とくせいほうていしきcharacteristic equationは、線型作用素せんけいさようそlinear operatorの多項式たこうしきpolynomial p(T) を固有こゆうベクトルeigenvector上うえで数かずの方程式ほうていしき p(λ)=0 に変換へんかんしたものである。

14演習えんしゅうリンク

data/exercise/math/linear-operator/linear-operator-equations.exercise.n.md

15関連かんれんリンク

data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md data/lecture/math/linear-algebra/eigenvalues-and-eigenvectors.lecture.n.md data/lecture/math/linear-algebra/minimal-polynomial-basics.lecture.n.md data/lecture/math/linear-algebra/companion-matrix-basics.lecture.n.md data/lecture/math/sequence/shift-operators-and-recurrences.lecture.n.md
raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる