markdown
シフト作用素と漸化式md 6167471
lecture/math/sequence/shift-operators-and-recurrences.lecture.n.md
Download PDF

シフト作用素さようそと漸化式ぜんかしき

mathsequencerecurrenceshift-operatorlinear-operatorlecture

1導入どうにゅう

この講義こうぎでは、差分さぶんdifference、和分わぶんindefinite sum、漸化式ぜんかしきrecurrence relationを、列れつを 1 項こう前進ぜんしんさせるシフト作用素さようそshift operator E によって統一的とういつてきに記述きじゅつする。

差分さぶんはシフト後ごの列れつと元もとの列れつの差さであり、和分わぶんはその逆問題ぎゃくもんだいである。また、定数係数線型漸化式ていすうけいすうせんけいぜんかしきはシフトした列れつの線型結合せんけいけつごうとして表現ひょうげんできる。したがって、中心的ちゅうしんてきな作用素さようそは

(Ea)n=an+1

である。

この解釈かいしゃくにより、漸化式ぜんかしきは単たんなる項こうの更新規則こうしんきそくではなく、

P(E)a=b

という線型作用素方程式せんけいさようそほうていしきlinear operator equationになる。

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

列れつ a=(an) 全体ぜんたいからなる空間くうかんを考かんがえる。必要ひつように応おうじて、両側りょうがわ無限列むげんれつ (an)n∈Z や片側かたがわ列れつ (an)n[PARSE ERROR: Undefined("Command(\"ge\")")]0 を使つかう。ここでは、添字そえじが定義ていぎできる範囲はんいで議論ぎろんする。

シフト作用素さようそshift operator E を

(Ea)n=an+1

で定義ていぎする。

恒等作用素こうとうさようそidentity operator I は

(Ia)n=an

である。

前進差分ぜんしんさぶんforward difference Δ は

(Δa)n=an+1-an

である。したがって

Δ=E-I

である。

この規約きやくは Z 変換へんかんの講義こうぎと一致いっちし、両側りょうがわ Z 変換へんかんでは E が z による乗算じょうざんへ対応たいおうする。

3方針ほうしん

この講義こうぎでは、次つぎの順じゅんに確認かくにんする。

  1. E が線型作用素せんけいさようそlinear operatorであることを証明しょうめいする。
  2. Δ=E-I から差分さぶん公式こうしきを導みちびく。
  3. ΔF=f を和分わぶん、すなわち差分さぶんの逆問題ぎゃくもんだいとして定式化ていしきかする。
  4. 定数ていすう係数けいすう線型せんけい漸化式ぜんかしきを P(E)a=b と書かく。
  5. rn が E の固有列こゆうれつeigen-sequenceであることから、特性とくせい方程式ほうていしきを導みちびく。

この順序じゅんじょにより、各手法かくしゅほうの採用理由さいようりゆうを作用素さようその構造こうぞうから説明せつめいできる。

個別こべつの入試にゅうし型かたや、固定点こていてん・階差かいさ・和わの消去しょうきょ・連立れんりつ化か・帰納法きのうほうへの使つかい分わけは、次つぎの講義こうぎで整理せいりする。

data/lecture/math/sequence/recurrence-types-and-strategies.lecture.n.md data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md data/lecture/math/linear-operator/polynomial-operators-and-eigenvalue-problems.lecture.n.md

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

4.1主張しゅちょう

列れつ空間上くうかんじょうの E は線型作用素せんけいさようそlinear operatorである。

4.2証明しょうめい

列れつ a,b とスカラー α,β を任意にんいに取とる。このとき

\begin{aligned} (E(\alpha a+\beta b))_n &=(\alpha a+\beta b)_{n+1}\\ &=\alpha a_{n+1}+\beta b_{n+1}\\ &=(\alpha Ea+\beta Eb)_n \end{aligned}

である。これはすべての n で成立せいりつするので

E(αa+βb)=αEa+βEb

である。したがって E は線型せんけいである。

5差分さぶんは E-I である

前進ぜんしん差分さぶんは

(Δa)n=an+1-an

である。一方いっぽう、

((E-I)a)n=(Ea)n-(Ia)n=an+1-an

である。したがって

Δ=E-I

である。

ここから、Δ の線型性せんけいせいlinearityも直ただちに従したがう。E と I が線型せんけいであるため、その差さ E-I も線型せんけいである。実際じっさい、

Δ(αa+βb)=(E-I)(αa+βb)=α(E-I)a+β(E-I)b=αΔa+βΔb

である。

6二項定理にこうていりによる高階差分こうかいさぶん

Δ=E-I であり、E と I は可換かかんである。したがって

Δm=(E-I)m

に二項定理にこうていりbinomial theoremを適用てきようできる。

Δm=∑k=0m(-1)k(mk)Em-k

である。列れつ a に作用さようさせると

(Δma)n=∑k=0m(-1)k(mk)an+m-k

を得える。

この公式こうしきの符号ふごうと添字そえじは、(E-I)m の二項展開にこうてんかいによって決定けっていされる。

7和分わぶんは差分さぶんの逆問題ぎゃくもんだいである

不定和分ふていわぶんindefinite sumの基本きほん問題もんだいは

ΔF=f

を満みたす列れつ F を決定けっていすることである。これは

F(n+1)-F(n)=f(n)

を意味いみする。

n=a,a+1,…,b-1 について両辺りょうへんを加算かさんすると

∑n=ab-1f(n)=∑n=ab-1(F(n+1)-F(n))

である。右辺うへんは望遠和ぼうえんわtelescoping sumなので

∑n=ab-1f(n)=F(b)-F(a)

を得える。

これは連続れんぞくの微積分びせきぶんで

∫abF′(x)dx=F(b)-F(a)

となることの離散りさん版はんである。つまり、和わ分ふんとは Δ の逆問題ぎゃくもんだいである。

8積せきの差分さぶんにおけるシフト項こう

連続れんぞく微分びぶんでは積つの微分びぶんが

(uv)′=u′v+uv′

となる。しかし差分さぶんでは、ずれが残のこる。

実際じっさい、

\begin{aligned} \Delta(uv)(n) &=u(n+1)v(n+1)-u(n)v(n)\\ &=u(n+1)v(n+1)-u(n)v(n+1)\\ &\qquad +u(n)v(n+1)-u(n)v(n)\\ &=v(n+1)(u(n+1)-u(n))+u(n)(v(n+1)-v(n)) \end{aligned}

である。したがって

Δ(uv)=(Ev)Δu+uΔv

である。

右辺うへんの (Ev) は、前進差分ぜんしんさぶんが 1 項先こうさきの値あたいを含ふくむことを表あらわす。離散計算りさんけいさんでは、このシフトを省略しょうりゃくして連続微分れんぞくびぶんの公式こうしきを適用てきようしてはならない。

9部分和公式ぶぶんわこうしきは離散版りさんばんの部分積分ぶぶんせきぶん

積つの差分さぶん公式こうしき

Δ(uv)=uΔv+(Ev)Δu

を変形へんけいすると

uΔv=Δ(uv)-(Ev)Δu

である。両辺りょうへんを和わ分ふんすると

∑uΔv=uv-∑(Ev)Δu

という形式けいしきを得える。有限範囲ゆうげんはんいでは

∑n=ab-1u(n)Δv(n)=u(b)v(b)-u(a)v(a)-∑n=ab-1v(n+1)Δu(n)

が成立せいりつする。これは連続れんぞくの部分積分ぶぶんせきぶん

∫udv=uv-∫vdu

の離散りさん版はんである。

ただし、右辺うへんに v ではなく Ev が現あらわれる。この一項いっこうの偏移へんいが、差分法さぶんほうでシフト作用素さようそshift operatorを明示めいじする理由りゆうである。

10漸化式ぜんかしきは P(E)a=b である

定数ていすう係数けいすう線型せんけい漸化式ぜんかしき

an+k=c1an+k-1+c2an+k-2+⋯+ckan+bn

を考かんがえる。左辺さへんへ集あつめると

an+k-c1an+k-1-c2an+k-2-⋯-ckan=bn

である。E を用もちいると

(Ek-c1Ek-1-c2Ek-2-⋯-ckI)a=b

と書かける。

つまり

P(E)a=b

である。ただし

P(t)=tk-c1tk-1-c2tk-2-⋯-ck

である。

この段階だんかいで、漸化式ぜんかしきは一般いっぱんの線型作用素方程式せんけいさようそほうていしきlinear operator equationになっている。したがって

P(E)a=bissolvable⟺b∈ImP(E)

であり、特解とくかい ap が 1 つ存在そんざいすれば

a=ap+h,h∈kerP(E)

である。

11固有列こゆうれつと特性方程式とくせいほうていしき

両側列りょうがわれつでは r≠0 とし、片側列かたがわれつでは n[PARSE ERROR: Undefined("Command(\"ge\")")]0 で rn を定義ていぎする。この幾何列きかれつに対たいして

E(rn)=rn+1=rrn

である。つまり rn は E の固有列こゆうれつeigen-sequenceであり、固有値こゆうちは r である。

したがって

P(E)(rn)=P(r)rn

である。よって同次どうじ方程式ほうていしき

P(E)a=0

に an=rn が解かいとなる条件じょうけんは

P(r)=0

である。これが特性方程式とくせいほうていしきcharacteristic equationである。

特性とくせい方程式ほうていしきは、rn が E の固有こゆう列れつであるために現あらわれる。したがって、定数ていすう係数けいすう線型せんけい漸化式ぜんかしきの特性とくせい方程式ほうていしき法ほうは、固有値問題こゆうちもんだいeigenvalue problemの離散りさん版はんである。

12例れい:Fibonacci 型かた漸化式ぜんかしき

漸化式ぜんかしき

fn+2=fn+1+fn

は

(E2-E-I)f=0

である。したがって

P(t)=t2-t-1

であり、特性とくせい方程式ほうていしきは

r2-r-1=0

である。根こんは

r1=1+52,r2=1-52

である。r1n と r2n はそれぞれ kerP(E) に属ぞくする。異ことなる根こんに対応たいおうする 2 列れつは線型独立せんけいどくりつである。また、この二階漸化式にかいぜんかしきの解かいは f0,f1 によって一意いちいに決定けっていされるため、解空間かいくうかんは 2 次元じげんである。したがって、これらは解空間かいくうかんの基底きていをなし、一般解いっぱんかいは

fn=C1r1n+C2r2n

を得える。定数ていすう C1,C2 は初期条件しょきじょうけんで決きまる。

ここで重要じゅうようなのは、rn を試ためす理由りゆうである。rn は E の固有こゆう列れつであり、P(E) が P(r) に落おちるからである。

13定数係数ていすうけいすうと非定数係数ひていすうけいすうの違ちがい

定数ていすう係数けいすうなら、漸化式ぜんかしきは P(E) の形かたちで書かける。係数けいすうが定数ていすうであるため、E と係数けいすうが可換かかんであり、多項式たこうしきとして扱あつかえる。

しかし係数けいすうが n に依存いぞんすると状況じょうきょうが変かわる。たとえば

an+1-nan=0

を考かんがえる。M を

(Ma)n=nan

で定義ていぎされる掛け算作用素ざんさようそmultiplication operatorとすると、この方程式ほうていしきは

(E-M)a=0

の形かたちで書かける。しかし一般いっぱんに

EM≠ME

である。実際じっさい、

(EMa)n=(Ma)n+1=(n+1)an+1

だが

(MEa)n=n(Ea)n=nan+1

である。

したがって、非定数係数ひていすうけいすうでは通常つうじょうの多項式たこうしき P(E) として表現ひょうげんできず、特性方程式法とくせいほうていしきほうは一般いっぱんには成立せいりつしない。定数係数ていすうけいすうという仮定かていは、作用素さようその可換性かかんせいを保証ほしょうする本質的条件ほんしつてきじょうけんである。

14Newton の前進公式ぜんしんこうしき(n[PARSE ERROR: Undefined("Command(\"ge\")")]0)

Δ=E-I なので

E=I+Δ

である。したがって

En=(I+Δ)n

である。I と Δ は可換かかんなので、二項定理にこうていりより

En=∑k=0n(nk)Δk

である。

両辺りょうへんを列れつ f の 0 番目ばんめに作用さようさせると

f(n)=(Enf)0=∑k=0n(nk)(Δkf)0

を得える。これは初期値しょきちと初期しょき差分さぶんから列れつを復元ふくげんする公式こうしきである。

15成立範囲せいりつはんい

この講義こうぎの P(E) 表示ひょうじは、主おもに定数係数ていすうけいすうの線型漸化式せんけいぜんかしきに対たいして有効ゆうこうである。非線型漸化式ひせんけいぜんかしきには線型作用素方程式せんけいさようそほうていしきの構造こうぞうがない。非定数係数ひていすうけいすうの線型漸化式せんけいぜんかしきでは、係数作用素けいすうさようそとシフト作用素さようそが一般いっぱんに可換かかんでないため、通常つうじょうの特性方程式法とくせいほうていしきほうは適用てきようできない。

また、片側列かたがわれつでは E-1 が常つねに自然しぜんに定義ていぎできるとは限かぎらない。後退差分こうたいさぶんや両側りょうがわシフトを扱あつかうときは、列空間れつくうかんをどの添字集合そえじしゅうごう上じょうに置おくかを明示めいじする必要ひつようがある。

16主要公式しゅようこうしき

[PARSE ERROR: Undefined("Command(\"boxed\")")](Ea)n=an+1
[PARSE ERROR: Undefined("Command(\"boxed\")")]Δ=E-I
[PARSE ERROR: Undefined("Command(\"boxed\")")]ΔF=f(indefinitesumproblem)
[PARSE ERROR: Undefined("Command(\"boxed\")")]P(E)a=b(constantcoefficientlinearrecurrence)[PARSE ERROR: Undefined("RBrace")]
[PARSE ERROR: Undefined("Command(\"boxed\")")]E(rn)=rrn,P(E)(rn)=P(r)rn
[PARSE ERROR: Undefined("Command(\"boxed\")")]P(r)=0(characteristicequation)

17要約ようやく

漸化式ぜんかしきrecurrence relationは、シフト作用素さようそshift operator E による線型作用素方程式せんけいさようそほうていしきlinear operator equationであり、特性とくせい方程式ほうていしきは E の固有こゆう列れつ rn から生しょうじる。

18関連かんれんリンク

data/lecture/math/analysis/introduction-to-z-transform.lecture.n.md data/lecture/math/sequence/sequences-and-recurrences.lecture.n.md data/lecture/math/sequence/recurrence-types-and-strategies.lecture.n.md data/lecture/math/sequence/first-order-recurrences.lecture.n.md data/lecture/math/sequence/difference-equation-basics.lecture.n.md data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md data/lecture/math/linear-operator/polynomial-operators-and-eigenvalue-problems.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
タブを全て閉じる