OKWAVEのAI「あい」が美容・健康の悩みに最適な回答をご提案!
-PR-
解決
済み

複素積分~フィボナッチ? (その2・その3)

  • すぐに回答を!
  • 質問No.112760
  • 閲覧数106
  • ありがとう数2
  • 気になる数0
  • 回答数4
  • コメント数0

お礼率 53% (212/400)

前回
http://oshiete1.goo.ne.jp/kotaeru.php3?q=110649


『zを複素数とする時、数列x_n(n=0,1,2,..., x_nは実数)に対する変換X(z)を以下のように定義する。
    X(z) = Σ_n=0~∞ x_n z^(-n)
この時以下の問いに答えよ。
(1) |z|>Rの領域において、X(z)は収束するとする。この領域内の原点を含む閉曲線をCとする時、逆変換は
    x_n = (1/2πi) ∫○C X(z) z^(n-1) dz    (∫○CはCを経路とする周回積分記号のつもり。)
となる事を証明せよ。
(2)x_n+2 = x_n+1 + x_n (n=0,1,2,..., x_0=x_1=1)の時、X(z)を求めよ。
(3)前問で求めたX(z)を逆変換する事によって、x_nを求めよ。』

の(1)について解答を頂き解決しました。
ひきづづき(2)(3)に関して教えてください。
通報する
  • 回答数4
  • 気になる
    質問をブックマークします。
    マイページでまとめて確認できます。

質問者が選んだベストアンサー

  • 回答No.2
レベル13

ベストアンサー率 64% (700/1089)

[1] x_n+2 = x_n+1 + x_n
の両辺に z^{-n} をかけて n=0~∞ の和を取る.
x_n のところは X(z) そのもの.
左辺は
[2] Σ_{n=0}^∞ x_{n+2} z^{-n}
   = z^2 Σ_{n=0}^∞ x_{n+2} z^{-n-2}
   = z^2 [ Σ_{m=0}^∞ x_m z^{-m} - x_0 - x_1 z^{-1} ]
   = z^2 [ X(z) - x_0 - x_1 z^{-1} ]
になる.
途中で m=n+2 とおいている.
もともとは m=2~∞ だが,m=0~∞ にして,m=0,1 の分を引いて修正した.
同様にして
[3] Σ_{n=0}^∞ x_{n+2} z^{-n}
   = z [ X(z) - x_0]
[1]に代入整理して,x_0 = x_1 = 1 を使えば
[4]  X(z) = z^2 / (z^2 - z -1)
これが(2)の答.

X(z) がわかったから
[5]  x_n = (1/2πi) ∫○C X(z) z^(n-1) dz
に代入すればよい.
積分路は十分大きな |z| で一周するようにとる(後述).
1位の極 z = α,βが積分路中にあるから,そこでの留数を拾えばよい.
[4]の分母を
[6]  z^2 - z - 1 = (z -α)(z -β)
[7]  α = (√5 + 1)/2,  β=(-√5 + 1)/2
と書いて,
[7]  ∫○C X(z) z^(n-1) dz
    = ∫○C {z^(n+1) / (z -α)(z -β)} dz
    = 2πi [α^{n+1}/(α-β) + β^{n+1}/(β-α)]
あとは,整理して
[8]  x_n = (1/√5) (α^{n+1} - β^{n+1})
が求める答.

十分大きな |z| で一周するように積分路を取るのは,
最初の級数 X(z) の収束条件と関係がある.
"収束半径" (今は逆べき級数だから,|z|>R なら収束,ということになる)
のRがちょうどαになっている.
それは,
http://oshiete1.goo.ne.jp/kotaeru.php3?q=94834
で議論したように,コーシー・アダマールの定理やダランベールの定理によればよい.
したがって,十分大きな |z| で一周とは,
|z|>α で一周ということ.
補足コメント
taropoo

お礼率 53% (212/400)

(1)はなるほどです。
どこからこのような発想が出てくるんでしょうか。
見えている人と見えていない者の差を歴然と感じてしまいます。

ともかく(1)はOKです。

(2)なんですが、
> 1位の極 z = α,βが積分路中にあるから
これはどうしたら分かるのでしょうか?
分母が0になるからαとβが特異点になるのは分かるのですが
それが真性特異点なのか極なのか(極なら何位なのか)除去可能な特異点なのか、その見分け方が分かりません。
ローラン級数に展開すれば良いんでしょうけど、それをせずに手っ取り早く見分ける方法ってありませんか?
ちなみに
z^(n+1) / (z -α)(z -β)
をローラン展開しようと頑張っていますが出来ませんのでそちらの方もよろしくお願いします。
投稿日時 - 2001-08-05 22:44:57
-PR-
-PR-

その他の回答 (全3件)

  • 回答No.1
レベル10

ベストアンサー率 40% (54/135)

前回からの続きということで、 z^2X(z)-zX(z)-X(z) = (z^2-z-1)X(z) = x_0z^2 + (x_0-x_1)z から X(z) = (x_0z^2 + (x_0+x_1)z)/(z^2-z-1) として、これの逆変換とかいうので求めた場合と X(z) = Σ_n=0~∞ x_n z^(-n) から求めた場合で比較すれば フィボナッチ列の一般項が出てくるの ...続きを読む
前回からの続きということで、
z^2X(z)-zX(z)-X(z)
= (z^2-z-1)X(z)
= x_0z^2 + (x_0-x_1)z
から
X(z) = (x_0z^2 + (x_0+x_1)z)/(z^2-z-1)
として、これの逆変換とかいうので求めた場合と
X(z) = Σ_n=0~∞ x_n z^(-n)
から求めた場合で比較すれば
フィボナッチ列の一般項が出てくるのではないでしょうか?
z^2-z-1が既にそうなりまっせといっている感じですよね。
補足コメント
taropoo

お礼率 53% (212/400)

物分りが悪くて申し訳ありません。
もう少し詳しく教えていただけますでしょうか。
何分社会人ですものですから数学の常識的なものを知りませんものですから。

よろしくお願いいたします。
投稿日時 - 2001-08-02 23:31:40


  • 回答No.3
レベル10

ベストアンサー率 40% (54/135)

siegmundさんが答えられているなぁと思っていたら 自分の回答に間違いがありました。あまり意味はないんですが。 (z^2-z-1)X(z) = x_0z^2 + (x_0-x_1)z から、なぜか X(z) = (x_0z^2 + (x_0+x_1)z)/(z^2-z-1) となってしまっていました。 X(z) = (x_0z^2 + (x_0-x_1)z)/(z^2-z-1) ...続きを読む
siegmundさんが答えられているなぁと思っていたら
自分の回答に間違いがありました。あまり意味はないんですが。

(z^2-z-1)X(z) = x_0z^2 + (x_0-x_1)z
から、なぜか
X(z) = (x_0z^2 + (x_0+x_1)z)/(z^2-z-1)
となってしまっていました。
X(z) = (x_0z^2 + (x_0-x_1)z)/(z^2-z-1)
でした。大変失礼しました。
  • 回答No.4
レベル13

ベストアンサー率 64% (700/1089)

siegmund です. > (1)はなるほどです。 > どこからこのような発想が出てくるんでしょうか。 う~ん,威張るわけじゃないですが,これくらいなら瞬時に見えます. 形としては,フーリエ級数展開と,その係数の積分表示などに似ています. 三角関数などの直交関数で他の関数を展開するなどの経験があれば, 見当がつくんじゃないでしょうか. 私の同僚もほとんど瞬時に見える ...続きを読む
siegmund です.

> (1)はなるほどです。
> どこからこのような発想が出てくるんでしょうか。

う~ん,威張るわけじゃないですが,これくらいなら瞬時に見えます.
形としては,フーリエ級数展開と,その係数の積分表示などに似ています.
三角関数などの直交関数で他の関数を展開するなどの経験があれば,
見当がつくんじゃないでしょうか.
私の同僚もほとんど瞬時に見えると思います.

なぜこういう形の逆べき級数を考えるかは,母関数がもともとのアイデアでしょう.
例えば Hermite (エルミート) 多項式
H_0(x) = 1
H_1(x) = 2x
H_2(x) = 4x^2 - 2
H_3(x) = 8x^3 - 12x
H_4(x) = 16x^4 -48x^2 +12
.........
を扱うとき
  Σ_{n=0}^∞ H_n(x) t^n / n!
という関数(実は,exp(-t^2 + 2tx) になることが知られている)を使うと非常に便利です.
t で展開したときの係数がちょうど H_n(x) になっているというわけです.
これを使うと
  ∫_{-∞}^{∞} H_n(x) H_m(x) exp(-x^2) dx
なんていう面倒そうな積分が簡単に求められます.


> (2)なんですが、
>> 1位の極 z = α,βが積分路中にあるから
> これはどうしたら分かるのでしょうか?

問題の関数は z^(n+1) / (z -α)(z -β) ですね.
z^(n+1) は異常性を与えませんから,問題は 1/(z -α)(z -β) です.
部分分数分解をすればOKです.
  1/(z -α)(z -β) = a/(z -α) + b/(z -β)
とおいて,両辺が等しくなるように a,b を定めれば,
  1/(z -α)(z -β) = [1/(α-β)] [1/(z -α) - 1/(z -β)]
となり,z =α と z =β が1位の極であることがわかります.
留数も見えていますね.

もし,(z+γ) /(z-α)^2 (z-β) なんかでしたら
  (z+γ) /(z-α)^2 (z-β) = (az + b)/(z -α)^2 + c/(z -α) + d/(z -β)
とおいて,a,b,c,d を決めることになります.
お礼コメント
taropoo

お礼率 53% (212/400)

コツとかの問題じゃなくラベルの差ですね。;-P
ありがとうございました。
投稿日時 - 2001-08-08 14:47:41
このQ&Aで解決しましたか?
関連するQ&A
-PR-
-PR-
このQ&Aにこう思った!同じようなことあった!感想や体験を書こう
このQ&Aにはまだコメントがありません。
あなたの思ったこと、知っていることをここにコメントしてみましょう。

その他の関連するQ&A、テーマをキーワードで探す

キーワードでQ&A、テーマを検索する
-PR-
-PR-
-PR-

特集


いま みんなが気になるQ&A

関連するQ&A

-PR-

ピックアップ

-PR-
ページ先頭へ