- ベストアンサー
微分方程式のべき級数展開による解法
- べき級数展開を用いて、微分方程式の解を求める方法について説明します。
- 微分方程式の係数を比較することにより、各項の係数の関係式を得ることができます。
- n>=3の場合の係数の関係式は、階乗によって表されます。
- みんなの回答 (8)
- 専門家の回答
質問者が選んだベストアンサー
関連するQ&A
- べき級数で解く微分方程式
次の微分方程式の解を 式(5.1) = y(x) = Σ[i=0,∞] ( a_[i] * x^i ) のべき級数を用いて求めよ。 x (dy/dx) - y = x^k (ただし、kは1以外の自然数) 解答 y を式(5.1)のべき級数で展開し、微分方程式に代入して係数a_iについての関係式を求める。 (1) べき級数展開から次の式を得る。 x Σ[i=0,∞] (i+1)( a_[i+1] * x^i ) - Σ[i=0,∞] ( a_[i] * x^i ) = x^k xの次数ごとに両辺の係数を比較すると、n≠kなるnについて (n-1)a_[n] = 0 となる。 ←疑問点 n≠1 (n≠k) に対して a_[n] = 0 であり、(k-1) * a_[k] = 1より y = 1/(k-1) * x^k を得る。 n=1に対しては、a_[n] = a_[1] ≠ 0でも(n-1) * a_[n] = 0となる。 実際、y = 1/(k-1) * x^k + ax (aは任意の定数) が微分方程式の解となる。 ・・・と本に書いてありますが、「疑問点」のところの比較の方法が分かりません。 まず、i が 0 から n まで変化する過程を自分で計算してみました。 i=0: x * (0+1) a_[0+1] * x^0 - a_[0] * x^0 = a_[1] * x - a_[0] i=1: x * (1+1) a_[1+1] * x^1 - a_[1] * x^1 = 2a_[2] * x^2 - a_[1] * x i=2: x * (2+1) a_[2+1] * x^2 - a_[2] * x^2 = 3a_[3] * x^3 - a_[2] * x^2 : i=n: x * (n+1) a_[n+1] * x^n - a_[n] * x^n = (n+1) a_[n+1] * x^(n+1) - a_[n] * x^n これらを使って「xの次数ごとに両辺の係数を比較する」んですよね。 しかし左辺だけでも、xの次数が1つずつズレていますよね・・・? これらと x^k を具体的にどうやって比較するのでしょうか? x^2ならx^2だけでまとめるんですか? それともx^3とx^2が混ざった形で比較するのですか(どうやってやるのか分かりませんけども)? どうか教えてください。お願いします。
- ベストアンサー
- 数学・算数
- 微分方程式の級数解
次の微分方程式の解を 式(5.1) = y(x) = Σ[i=0,∞] ( a[i] * x^i ) のべき級数を用いて求めよ。 x^2 * (dy/dx) - y = x^2 解答 べき級数展開から次の式を得る。 x^2 * Σ[i=0,∞] (i+1)( a[i+1] * x^i ) - Σ[i=0,∞] ( a[i] * x^i ) = x^2 xの次数ごとに両辺の係数を比較すると、 a[0] = 0 a[1] = 0 a[2] = -1 a[n] = (n-1) a[n-1] (n>=3) なる関係式を得る。これより、n>=3について a[n] = (n-1) ! * a[2] = -(n-1) ! となる。したがって、微分方程式の級数解として y = -x^2 * Σ[i=0,∞] (i+1) ! * x^i ←この式の求め方が分かりません を得る。 ・・・と本に書いてありますが、 y = -x^2 * Σ[i=0,∞] (i+1) ! * x^i の求め方が分かりません。 a[n] = -(n-1) !まで分かっているので、後は代入するだけだと思っていたのですが、やってみると答えが合いません。例えば、 x^2 * Σ[i=0,∞] (i+1)( a[i+1] * x^i ) - Σ[i=0,∞] ( a[i] * x^i ) = x^2 に a[n] = -(n-1) ! a[i] = -(i-1) ! a[i-1] = -(i-2) ! a[i+1] = -i ! など各種取り揃えておいて代入すると x^2 * Σ[i=0,∞] (i+1)( a[i+1] * x^i ) - Σ[i=0,∞] ( a[i] * x^i ) = x^2 x^2 * Σ[i=0,∞] (i+1)( -i ! * x^i ) - Σ[i=0,∞] { -(i-1) ! * x^i } = x^2 (i+1)i ! = (i+1) ! と考えれば x^2 * Σ[i=0,∞] -(i+1) ! * x^i + Σ[i=0,∞] (i-1) ! * x^i = x^2 -x^2 * Σ[i=0,∞] (i+1) ! * x^i + Σ[i=0,∞] (i-1) ! * x^i = x^2 この前半の項が奇しくもこの本の答え y = -x^2 * Σ[i=0,∞] (i+1) ! * x^i と同じになります。 ということは、この後半の項はゼロになるべきということですか?でも、ならないですよね? それとも私の計算が間違っているのでしょうか? どうか正しい解き方を教えてください。お願いします。
- ベストアンサー
- 数学・算数
- 微分方程式の級数解の求め方
微分方程式の級数解の求め方について教えてください。 y' = a^2・y, y(0) = 1 の解が y = f(x) = Σ[n=0→∞]c(n)・x^n であるとします。 この場合に、係数 c10 値と、f(1) の値を求めたいと思います。 以下のように辿ってみましたが、途中でわからなくなりました。 解の式を微分して、 y' = c1 + 2c2・x + 3c3・x^2 + ... 元の方程式を展開すると、 y' = a^2( c0 + c1・x + c2・x^2 + ... ) 両式と y(0) = 1 より、 c1 = c0・a^2 = a^2 2c2・x = c1・x → c2 = c1・x / 2・x = c1 / 2 → a^2 / 2! 3c3・x^2 = c2・x^2 → c3 = c2・x^2 / 3・x^2 = c2 / 3 → a^2 / 3! ゆえに c(n) = a^2 / n! このあと c10 を算出するために上式の a の値は?などとわからなくなりました。 ここまでに誤りがないか、このあとをどうすればよいか、教えていただけないでしょうか。 よろしくお願いします。
- ベストアンサー
- 数学・算数
- ベルヌーイ型の微分方程式を超幾何級数で解きたいです
ベルヌーイ型の微分方程式 y'+P(x)y=Q(x)y^n をガウス形に変換し、超幾何級数を用いて解く方法を教えてください。 よろしくお願いいたします。 P(x)=a/(x-b)、Q(x)=c/(x^2+bx)、n=-1/2 (a, b,c は定数) の場合を解こうとしています。 普通にベルヌーイ型で解こうとすると出来ない積分がでてきます。 そこで、ガウス型に変換し、超幾何級数で解く方法を模索しています。 教科書に載っているやり方で、級数を y=Σ(k=0,∞) ck x^(ρ+k) のように置いても、 また、式を更に微分し常2回微分方程式にしてガウス型への変換を試みていますが、できません。 やり方をご存知の方、いらっしゃいましたらご教示ください。 よろしくお願いいたします。
- 締切済み
- 数学・算数
- 微分方程式の級数解 a[0] * x^n
微分方程式 (d^2 y)/(dx^2) + (1/x) (dy/dx) - (n^2/x^2) y = 0 (x>0) の級数解を、次の問いに従って求めよ。 ただし、n>0とする。 (1) 級数解を y(x) = x^c * Σ[i=0,∞] a[i] * x^i とおいたとき、指数cはどのように求まるか。ただし、a[0] ≠ 0であるとする。 解答 級数解を y(x) = x^c * Σ[i=0,∞] a[i] * x^i とおいて、項別に微分すると dy/dx = x^c * Σ[i=0,∞] (c+i)a[i] * x^(i-1) (d^2 y)/(dx^2) = x^c * Σ[i=0,∞] (c+i)(c+i-1)a[i] * x^(i-2) これを微分方程式に代入して x^c * Σ[i=0,∞] (c+i)(c+i-1)a[i] * x^(i-2) + (1/x) * x^c * Σ[i=0,∞] (c+i)a[i] * x^(i-1) + (n^2/x^2) * x^c * Σ[i=0,∞] a[i] * x^i = 0 x^c * Σ[i=0,∞] (c+i)(c+i-1)a[i] * x^(i-2) + x^c * Σ[i=0,∞] (c+i)a[i] * x^(i-2) + n^2 * x^c * Σ[i=0,∞] a[i] * x^(i-2) = 0 x^c * Σ[i=0,∞] { (c+i)(c+i-1) + (c+i) - n^2 } a[i] * x^(i-2) = 0 x^c * Σ[i=0,∞] [ (c+i) { (c+i-1) + 1} - n^2 ] a[i] * x^(i-2) = 0 x^c * Σ[i=0,∞] { (c+i)(c+i-1+1) - n^2 } a[i] * x^(i-2) = 0 x^c * Σ[i=0,∞] { (c+i)(c+i) - n^2 } a[i] * x^(i-2) = 0 x^c * Σ[i=0,∞] { (c+i)^2 - n^2 } a[i] * x^(i-2) = 0 x^(c-2)の係数について (c^2 - n^2)a[0] = 0 でなければならない。 したがって、a[0] ≠ 0の条件から c^2 = n^2 c = ±n と定まる。 (2) 一般解を級数解で求めよ。 解答 x^(i+c-2) (i=1,2,3,...)の係数について { (c+i)^2 - n^2 } = 0 でなければならない。 c=nのとき、 { (c+i)^2 - n^2 } = (n+i)^2 - n^2 = 2ni + i^2 ≠ 0 であるから、a[i] = 0 (i=1,2,3,...)となる。 すなわち、これに対応する解は a[0] * x^n ←これが分かりません ・・・とまだまだ続くのですが、a[0] * x^nになる理由が分かりません。 自分で考えてみますと、 Σ[i=0,∞] { (c+i)^2 - n^2 } a[i] * x^(i+c-2) = 0 で、(i=1,2,3,...)はすべてa[i] = 0になると言ってるのだから、残るはi=0のみ。 i=0: { (c+0)^2 - n^2 } a[0] * x^(0+c-2) = 0 { c^2 - n^2 } a[0] * x^(c-2) = 0 しかも、c=nなので { n^2 - n^2 } a[0] * x^(n-2) = 0 { 0 } a[0] * x^(n-2) = 0 ・・・x^(n-2)の係数について係数は0という結果になりました。これでいいんですか??? たとえ、{ (c+i)^2 - n^2 } = 2ni + i^2としても、i=0なので0ですよね? このa[0] * x^nはどうやって導いたのでしょうか? 教えてください。お願いします。
- ベストアンサー
- 数学・算数
- 微分方程式
次の微分方程式を解けという問題がわかりません。 y''+4y=sin2x 特性方程式s^2+4=0よりs=±2i(虚数解) 補助方程式の一般解はy=Asin2x+Bcos2x 与方程式の右辺を微分して生じる関数は2sin2x,2cos2xであるが、 これらは上の一般解に含まれている。重複度は2なので、 特殊解を求めるために、 y1=ax^2*sin2xとおく y1'=2a(xsin2x+x^2cos2x) y1''=2a(sin2x+4xcos2x-2x^2sin2x) これらを与方程式に代入すると 2asin2x+8axcos2x-4ax^2sin2x+4ax^2sin2x=sin2x となってしまって解けませんでした。どこを直せばいいでしょうか?
- 締切済み
- 数学・算数
- 微分方程式の解き方が分からず、困っています。
現在、試験に向けて微分方程式の勉強をしているのですが、下記の問題の解き方が分かりません。 教科書を参考に(1)は変数分離系、(2)は同次形、(3)は線形で解こうとしましたが、どの問題も積分するところで複雑な式になってしまい、解けれません。 分かる問題だけでも良いのでアドバイス、解き方を教えてください。よろしくお願いします。 (1)次の微分方程式の一般解を求めよ dy/dx=y^2+1 (2)次の微分方程式の一般解を求めよ y'=(y/x)(log(y/x)+1) (3)次の微分方程式の解でt=0のときx=1の条件を満たすものを求めよ x'cost+xsint=1
- ベストアンサー
- 数学・算数
お礼
ありがとうございます。 求めていたのはこういうのです。正直、両辺を(n-1)!で割ってみると、の部分が少し天下り的な気がしますが、そこはNo.4さんのお礼に書いたような定義で階乗になると当たりをつけていると解釈します。 > b[n]=b[n-1](n≧3) > 初項a[2]の公差0の等差数列 ←nが増えても同じ数字(=a[2]) > 公比1の等比数列の漸化式 ←nが増えても同じ数字(=a[2]) > といってもいいですが, > b[n]=b[n-1]=・・・=b[2]=a[2] 納得です。結局b[n]も初期値のa[2]と同じになる、というのはそうやって表すんですね。覚えておきます。 ありがとうございました!