markdown
差分方程式の基本md bcc718a
lecture/math/sequence/difference-equation-basics.lecture.n.md
Download PDF

差分方程式さぶんほうていしきの基本きほん

date2026-07-14document_iddoc_5c082dc1f3665a541ded29f583c8d4bbdescription差分方程式を前進差分とシフト作用素 E-I によって定式化し、一次方程式の解法と線型作用素方程式との対応を説明する。prerequisites数列と漸化式 / 一次漸化式 / シフト作用素と漸化式type講義content_typelecturestatusactiverelateddata/lecture/math/sequence/sequences-portal.lecture.n.md / data/lecture/math/sequence/first-order-recurrences.lecture.n.md / data/lecture/math/sequence/shift-operators-and-recurrences.lecture.n.md / data/lecture/math/sequence/recurrence-types-and-strategies.lecture.n.md / data/lecture/math/calculus/introduction-to-differential-equations.lecture.n.md / data/lecture/math/analysis/introduction-to-z-transform.lecture.n.md / data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md
mathsequenceundergraduatelecture

1導入どうにゅう

この講義こうぎでは、差分方程式さぶんほうていしきが離散的りさんてきな変化へんかを数列すうれつの項こうの関係かんけいとして記述きじゅつする方程式ほうていしきであることを説明せつめいする。連続的れんぞくてきな変化へんかを記述きじゅつする微分方程式びぶんほうていしきとの構造的対応こうぞうてきたいおうも明確めいかくにする。

漸化式ぜんかしきを離散時間りさんじかんの方程式ほうていしきとして定式化ていしきかし、微分方程式びぶんほうていしきとの構造的対応こうぞうてきたいおうを説明せつめいする。ここでは an+1-an を離散的りさんてきな変化量へんかりょうとして解釈かいしゃくする。

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

差分方程式さぶんほうていしきDifference equation とは、離散的りさんてきな複数ふくすうの添字そえじにおける未知列みちれつの値あたいの関係かんけいを表あらわす方程式ほうていしきである。高階こうかいの方程式ほうていしきは an+2 と an のような非隣接項ひりんせつこうを含ふくみうるが、この講義こうぎでは隣接りんせつする 2 項こうの一次方程式いちじほうていしきを主おもに扱あつかう。

前進差分ぜんしんさぶんForward difference とは、

Δan=an+1-an

で定義ていぎされる差さである。

3方針ほうしん

まず an+1-an を離散変化量りさんへんかりょうとして採用さいようする理由りゆうを説明せつめいする。次つぎに、一次差分方程式いちじさぶんほうていしきを変形へんけいし、等比数列とうひすうれつと望遠和ぼうえんわによる解法かいほうを導出どうしゅつする。

data/lecture/math/sequence/first-order-recurrences.lecture.n.md data/lecture/math/sequence/shift-operators-and-recurrences.lecture.n.md data/lecture/math/calculus/introduction-to-differential-equations.lecture.n.md

4直感的ちょっかんてきな説明せつめい

微分方程式びぶんほうていしきにおける y′(x) が連続的れんぞくてきな局所変化率きょくしょへんかりつを表あらわすのに対たいし、差分方程式さぶんほうていしきにおける an+1-an は添字そえじの 1 単位たんいの増加ぞうかに対応たいおうする増分ぞうぶんを表あらわす。

したがって、差分方程式さぶんほうていしきは、時間じかんや反復回数はんぷくかいすうが整数値せいすうちをとる現象げんしょうの記述きじゅつに適てきする。

5厳密げんみつな説明せつめい

5.11. 差分さぶんを使つかう理由りゆう

数列すうれつの変数へんすう n は 1, 2, 3, ... という離散値りさんちをとる。したがって微分びぶん

dadn

をそのまま考かんがえるより、

Δan=an+1-an

によって変化量へんかりょうを記述きじゅつすることが自然しぜんである。

5.22. 最もっとも基本きほんの形かたち

an+1=ran

の場合ばあい、

an=a1rn-1

である。これは等比数列とうひすうれつである。この関係式かんけいしきは、添字そえじが 1 増加ぞうかするごとに直前ちょくぜんの項こうを r 倍ばいする規則きそくを表あらわす。したがって、

a2=ra1,a3=ra2=r2a1,[PARSE ERROR: Undefined("Command(\"dots\")")]

となり、帰納的きのうてきに

an=a1rn-1

を得える。

5.33. 一次いちじ非同次ひどうじ差分方程式さぶんほうていしき

an+1=ran+b

を n≥1 で対象たいしょうとし、初期値しょきちを a1 とする。まず定数解ていすうかい α を求もとめると

α=rα+b

したがって、

α=b1-r(r≠1)

である。この定数解ていすうかいを、更新写像こうしんしゃぞう a↦ra+b によって不変ふへんな平衡点へいこうてんequilibrium pointという。平衡点へいこうてんからの偏差へんさを新あらたな列れつとすると、定数項ていすうこう b を消去しょうきょして同次方程式どうじほうていしきへ変換へんかんできる。

そこで

cn=an-α

とおくと

cn+1=an+1-α=ran+b-α=r(an-α)=rcn

となる。したがって

cn=Crn-1,C=c1=a1-α

で、

an=α+(a1-α)rn-1

である。すなわち、平衡点へいこうてんを減算げんざんすることにより、非同次方程式ひどうじほうていしきを同次方程式どうじほうていしきへ変換へんかんできる。

ただし r=1 のときは

an+1=an+b

であるから、

an+1-an=b

であり、各かく反復はんぷくにおける増分ぞうぶんは b である。実際じっさい、

a2-a1=b,a3-a2=b,[PARSE ERROR: Undefined("Command(\"dots\")")],an-an-1=b

を加算かさんすると中間項ちゅうかんこうが相殺そうさいされ、

an-a1=(n-1)b

となる。したがって、

an=a1+(n-1)b

となる。したがって、r≠1 の公式こうしきは r=1 に適用てきようできず、r=1 の場合ばあいには独立どくりつの導出どうしゅつを要ようする。

5.44. E-I による作用素表示さようそひょうじ

添字そえじを 1 増加ぞうかさせるシフト作用素さようそshift operator E を

(Ea)n=an+1

と定義ていぎする。このとき、前進差分ぜんしんさぶんは

Δa=(E-I)a

と表現ひょうげんされる。

差分さぶん方程式ほうていしき

Δa=f

は、線型作用素方程式せんけいさようそほうていしきlinear operator equation

(E-I)a=f

である。一般いっぱんの線型作用素方程式せんけいさようそほうていしきと同様どうように、1 つの特解とくかい a(p) が得えられれば、すべての解かいは

a=a(p)+h,(E-I)h=0

で与あたえられる。ここで (E-I)h=0 は hn+1=hn を意味いみするため、連続れんぞくする整数添字せいすうそえじの領域上りょういきじょうで h は定数列ていすうれつである。したがって、不定和分ふていわぶんにおいて積分定数せきぶんていすうに対応たいおうする対象たいしょうは、差分作用素さぶんさようその核かくである。

data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md

6数学的すうがくてきな解釈かいしゃく

6.1行列ぎょうれつによる解釈かいしゃく

一次いちじの差分方程式さぶんほうていしき

an+1=ran+b

は、定数座標ていすうざひょう 1 を付加ふかして

(an+11)=(rb01)(an1)

と表現ひょうげんでき、漸化式ぜんかしきは行列反復ぎょうれつはんぷくとして解釈かいしゃくされる。

この表現ひょうげんでは、漸化式ぜんかしきの反復適用はんぷくてきようが行列ぎょうれつの累乗るいじょうに対応たいおうする。定数係数線型同次高階漸化式ていすうけいすうせんけいどうじこうかいぜんかしきも、連続れんぞくする項こうを状態じょうたいベクトルにまとめることにより、伴随行列はんずいぎょうれつによる一次系いちじけいへ変換へんかんできる。この場合ばあい、伴随行列はんずいぎょうれつの対角化可能性たいかくかかのうせいと Jordan 構造こうぞうが解かいの形式けいしきを規定きていする。

この行列表現ぎょうれつひょうげんにより、一次いちじだけでなく二階にかいおよび高階こうかいの漸化式ぜんかしきも共通きょうつうの枠組わくぐみによって整理せいりできる。

6.2解析学的かいせきがくてきな解釈かいしゃく

an+1-an=kan

は

Δan=kan

であり、これは微分方程式びぶんほうていしき

y′=ky

の離散的対応物りさんてきたいおうぶつである。連続系れんぞくけいの解かいが指数関数しすうかんすうとなるのに対たいし、離散系りさんけいの解かいは等比数列とうひすうれつとなる。

6.3複素関数ふくそかんすうによる解釈かいしゃく

差分方程式さぶんほうていしきに現あらわれる rn は、複素指数関数ふくそしすうかんすうによって増減ぞうげんと振動しんどうを統一的とういつてきに表現ひょうげんできる。r=ρeiθ と極形式きょくけいしきで表あらわすと

rn=ρneinθ

である。ここで

einθ=cos(nθ)+isin(nθ)

である。したがって、振動しんどうを伴ともなう増減ぞうげんは、単一たんいつの複素指数型ふくそしすうがたによって記述きじゅつできる。

この表現ひょうげんは倍率ばいりつと回転かいてんを同時どうじに記述きじゅつする。特とくに二階差分方程式にかいさぶんほうていしきが共役複素根きょうやくふくそこんをもつ場合ばあい、実数解じっすうかいが三角関数さんかくかんすうによって表現ひょうげんされる。

母関数ぼかんすうgenerating functionとは、列れつを係数けいすうにもつ形式的けいしきてきまたは収束しゅうそくする冪級数べききゅうすうである。ここでは添字そえじを 0 から開始かいしするように再添字付さいそえじづけし、A(w)=∑n[PARSE ERROR: Undefined("Command(\"ge\")")]0anwn とする。漸化式ぜんかしきに wn を乗じょうじて総和そうわを取とると、シフトは w の代数因子だいすういんしと初期項しょきこうへ変換へんかんされる。w=z-1 と置おけば、これは片側かたがわ Z 変換へんかんの冪級数表示べききゅうすうひょうじである。

7適用方針てきようほうしん

  • 隣接項りんせつこう an+1,an の関係かんけいを差分方程式さぶんほうていしきとして定式化ていしきかする。
  • an+1=ran+b では、平衡点へいこうてんの減算げんざんによる同次化どうじかの可否かひを最初さいしょに検討けんとうする。
  • 微分方程式びぶんほうていしきとの構造的対応こうぞうてきたいおうを利用りようして、増加ぞうか、減衰げんすい、振動しんどうを解釈かいしゃくする。

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

ここでは一次いちじの基本形きほんけいだけを扱あつかった。また an+1=ran+b の解法かいほうは r≠1 と r=1 で場合分ばあいわけが必要ひつようである。高階こうかいの差分方程式さぶんほうていしきや非線形ひせんけいの差分方程式さぶんほうていしきでは、別べつの手法しゅほうが必要ひつようである。

9基本公式きほんこうしき

[PARSE ERROR: Undefined("Command(\"boxed\")")]Δan=an+1-an
[PARSE ERROR: Undefined("Command(\"boxed\")")]an+1=ran+b⇒an=b1-r+(a1-b1-r)rn-1(n[PARSE ERROR: Undefined("Command(\"ge\")")]1,r≠1)
[PARSE ERROR: Undefined("Command(\"boxed\")")]an+1=an+b⇒an=a1+(n-1)b(n[PARSE ERROR: Undefined("Command(\"ge\")")]1)

10要約ようやく

  • 差分方程式さぶんほうていしきは、数列すうれつの変化へんかを離散的りさんてきな変化量へんかりょうによって記述きじゅつする方程式ほうていしきである。
  • シフト作用素さようそによる表示ひょうじでは Δ=E-I であり、一般解いっぱんかいは特解とくかいと定数列ていすうれつの和わである。

11関連かんれんリンク

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