• ベストアンサー

微分方程式の一般解

微分方程式の一般解を求める問題なのですが、どうしてもよく分かりません。 y'6+4y'2=40x^3  (ここで'○は微分の回数を示すとします。また、以下ではD=d/dxのことです) 同時方程式(D^6+4D^2)y=D^2(D^4+4)y=D^2{(D^2+2)^2-4D^2}y=D^2(D^2+2+2D)(D^2+2-2D)y=D^2{(D+1)^2+1}{(D-1)^2+1}y=0 の基本解は{1、x、e^(-x)cosx、e^(-x)sinx、e^xcosx、e^xsinx} 次に特殊解Y(x)を求める。 非同次項R(x)=40x^3は同時方程式40D^4y=0の解だから、 特殊解Y(x)の式Yは同時方程式 40D^4・D^2{(D+1)^2+1}{(D-1)^2+1}y=0の解である。 この基本解{1、x、x^2、x^3、x^4、x^5、e^(-x)cosx、e^(-x)sinx、e^xcosx、e^xsinx}から与式の基本解を除いたx^2、x^3、x^4、x^5の一次結合として Y(x)=Ax^2+Bx^+Cx^4+Dx^5とおく。 与式の左辺に代入… と続いていくのですが、どうにもしっくりきません。 答えも y=(20/3)x^3+c1+c2x+e^(x/√2){c3cos(x/√2)+c4sin(x/√2)}+e^(-x/√2){c5cos(x/√2)+c6sin(x/√2)} となり、私の解からでは到底結びつくとは思えないです。 気になるのが 「非同次項R(x)=40x^3は同時方程式40D^4y=0の解だから」 としていますが、本当にこれで良いのか自信もありません。 もし間違えていたら解説をお願いします。 また、他に違うというようなところがあったら指摘してください。 回答、よろしくお願いします。

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

  • ベストアンサー
  • mmky
  • ベストアンサー率28% (681/2420)
回答No.1

参考程度に この場合の解の求め方は、 y'6+4y'2=40x ←(答えが正しいとすると) (D^6+4D^2)y=40x D^2(D^4+4)y=40x (D^4+4)y=∬40xdx (D^4+4)y=(20/3)x^3+c1+c2x とすればよいのでは。 回答は、y=(20/3)x^3+c1+c2x+・・・ の形になりますね。微分方程式は形式によって解の出し方が変わりますので注意しないとね。

その他の回答 (4)

  • mmky
  • ベストアンサー率28% (681/2420)
回答No.5

追伸まで 参考までに特別解の求め方を最後まで書いておきます。 y'6+4y'2=40x ←(この場合で) (D^6+4D^2)y=40x D^2(D^4+4)y=40x (D^4+4)y=∬40xdx (D^4+4)y=(20/3)x^3+c1+c2x η(x)=(1/4){1/(1+D^4/4)}{(20/3)x^3+c1+c2x} =(1/4){1-D^4/4+(D^4/4)^2-・}{(20/3)x3+c1+c2x} =(1/4){(20/3)x^3+c1+c2x } =(5/3)x^3+C1+C2x (註:Dは微分演算子だから級数展開した時、級数の第2項目以降はxが3乗べきまでだから0になる故} それから一般解の求め方の例は、 (D^4+4)y=(20/3)x^3+c1+c2x ここで、 (D^4+4)y=0 と置けば、D=±2i から, D=±(1±i) の4つの虚数根がでますね。 {(1+i)^2=1+2i-1=2i, (1-i)^2=1-2i-1=-2i} だから一般解は公式を使えば、 y=e^x(c3cosx+c4sinx)+e^-x(c5cosx+c6sinx) だから解は、nakaizuさんのご指摘のように、 y=(5/3)x^3+C1+C2x+e^x(c3cosx+c4sinx)+e^-x(c5cosx+c6sinx) の形式になりますね。 そのほか、一般解は例えば、4つの虚数根、±(1±i) を直接使って線形接続すれば、 y=C1e^(1+i)x+C2e^-(1+i)x+C3e^(1-i)x +C3e^-(1-i)x のように書くことも出来ますね。実数部と虚数部に分けて係数をうまくとると公式解のようにsin,cosであらわすことが出来るはずですね。 参考程度に

namisem
質問者

お礼

回答ありがとうございます。 お礼が遅れてしまい申し訳ありません。 細かい解説ありがとうございますm(_ _)m D^2(D^4+4)y=40x (D^4+4)y=∬40xdx この解き方はなるほどですね。他の問題にも使ってみます。 何度も補足ありがとうございました。

  • nakaizu
  • ベストアンサー率48% (203/415)
回答No.4

入力ミスがありました。 y'6+4y'2=40x^3 の解は y=x^5/2+c1+c2x+e^x(c3sin x+c4cos x)+e^(-x)(c5sin x+c6cos x) です。

  • mmky
  • ベストアンサー率28% (681/2420)
回答No.3

#1のmmkyです。#2のnakaizuさんのご指摘どおりですね。 y'6+4y'2=40x の場合、正しくは、 y=(5/3)x^3+c1+c2x+・・・ でしたね。 ごめん。

  • nakaizu
  • ベストアンサー率48% (203/415)
回答No.2

求め方はそれでいいと思いますよ。多少、回りくどいことをしていますが。 ただ、解答が正しいとすると元の方程式は y'6+y'2=40x だったはずです。 y'6+4y'2=40x^3 の解は y=x^5/2+c1+c2x+e^x(c3sin x+c4cos x)+e^(-1)(c5sin x+c6cos x) となります。 問題の取り違えでしょう。

namisem
質問者

お礼

回答ありがとうございます。 お礼が遅くなってしまい申し訳ありません。 う~ん、問題が間違っているみたいですね… どう考えても、上のような答えにはならないですよね 先生に指摘してみたいと思います。 ありがとうございました。

関連するQ&A

  • 微分方程式

    微分方程式の特殊解のおき方がわかりません y"+y=secx 同時微分方程式を解くと y=c1cosx+c2sinx となるところまでできるのですが ここから y(x)=Acosx+Bsinx とおいて計算してもうまくいきません お願いします

  • 1階の線形微分方程式

    1階の線形微分方程式 次の微分方程式の解き方が分かりません。いちおう、自分でもやりましたが、答えを先生が教えてくれないので困っています。さらに(3)はさっぱりです。 (1)y'+2y=6e^x (2)y'+y=sinx (3)xy'-2y=x^3e^x (1),(2)の自分なりで解いてみた答え (1) λ+2=0 λ= -2 よってこの微分方程式の一般解は y1=Ce^-2x ここで、yp=k1*e^x とおいて、ypを微分方程式内に代入をすると、 yp'+2yp=k1*e^x+2k1*e^x=3k1*e^x=6e^x k1=2 y2=2e^x よって y=y1+y2=C*e^-2x+2e^x (2) λ+1=0 λ= -1 よって、求める一般解は y1=Ce^-x ここで、特殊解を考えると yp=L*sinx+M*cosx yp'=L*cosx-M*sinx これを微分方程式に代入して yp'+yp=(L*sinx+M*cosx)+(L*cosx-M*sinx)=(L-M)sinx+(L+M)cosx ここで、 L-M=1 L+M=0 これを解いて L=1/2,M=-1/2 y2=1/2*sinx-1/2*cosx よって、y=y1+y2=Ce^-x+1/2*sinx-1/2*cosx

  • 2階線形微分方程式

    質問なんですが。 微分方程式で y''-2y'+y=(e^x)cosx という問題があるんですが この特殊解を求めるときに y=a(e^x){cosx+(i)sinx} とおいて、これを微分方程式に代入すれば 特殊解がy=-(e^x)cosx となるとなっているのですが。 y=a(e^x){cosx+(i)sinx} とおくというのがよく分かりません。 なんでiがでてくるのかとかも…。 僕は最初 y=a(e^x)cosx+b(e^x)sinxとおいて計算していました。 質問がわかりにくかったらすいません。

  • 2階の非同時線形微分方程式の特殊解

    2階非同時線形微分方程式を解いているのですが、わからない点があるため教えてください。一般解はわかるのですが、特殊解が答と一致しません。どこが間違っているか教えてください。 問1 y''+y'-6y=10e^(2x) 特殊解を求めると y0=ae^(2x)とおくと y0'=2ae^(2x) y0''=4ae^(2x) よって、4ae^(2x)+2ae^(2x)-6(ae^(2x))=10e^(2x) となり、左辺が0になってしまうのですが、どこを直せばいいでしょうか。 答では2xe^2xが特殊解になっています。 問2 y''+y'=x+2 特殊解を求めると y0=ax+b とおくと y0'=a y0''=0 よって、a=x+2 となり、y0=x^2+2xとなったのですが、答えではy0=(1/2)x^2+xとなっています。どこが間違えているか教えてください。 問3 y''+y=5e^xcosx 特殊解を求めると、 y0=e^x(acos+bsinx) y0'=e^x(-asinx+bcosx) y0''=e^x(-acosx-bsinx) よって、 e^x(-acosx-bsinx)+e^x(acos+bsinx)=5e^xcosx となり、問1同様左辺が0になってしまいます。 答では特殊解はe^x(cosx+2sinx)となっています。 問題が多くて申し訳ありませんが、回答お願いします。

  • 連立常微分方程式の一般解について!

    連立常微分方程式の一般解を消去法で求める問題なのですが、解き方を教えていただきたいです。 ・Dy1+y2+y3=0 ・y1+Dy2+y3=sinx-cosx ・y1+y2+Dy3=2sinx ただし、y1=y1(x),y2=y2(x),y3=y3(x),D=d/dx 宜しくお願い致します。

  • フーリエ展開 微分方程式の一般解

    y''+y=f(x)という微分方程式の一般解を求める。 ただし、f(x)=x^2 (-π<x≦π) f(x+2π)=f(x)であるとする。 上記のような問題なのですが、まずf(x)をフーリエ展開すると、f(x)=π^2/3+4Σ(-1)^n/n^2となりました。 この後、係数比較を行うために、yn=Acosnx+Bsinnxとおき、yn'+yn=4(-1)^n/n^2となり、AとBの値を求めることができました。 しかし、この問題の解答はy=c1cosx+c2sinx+(π^2/3)+2xsinx+4Σ[2→∞]{(-1)^n/(1-n^2)n^2}cosnxとなるようで、四つ目の項の2xsinxの出所がよくわかりません。 4Σ[2→∞]{(-1)^n/(1-n^2)n^2}cosnxの部分は、nが2以上のときの場合を表していて、π^2/3はn=0のときの場合を表している。つまり、2xsinxという部分はn=1のときの場合を求めているのではないかというところまで推測できたのですが、何故このような2xsinxという値が出てくるのかわかりません。 n=1のとき、1-n^2が0になってしまうため、別に求めなければいけないというのはなんとなくわかるのですが、上手く2xsinxの値まで辿り着きません。 長くなりましたが、この問題についてわかる方、ご教授お願いします。

  • 線形微分方程式について

    参考書にこのような問題があります、 カッコ内の一組の関数は、与えられた微分方程式の基本解であることをしめせ、(解の確認は直接代入、一次独立性はロンスキアンを用いよ) x^2y''-2xy'+(x^2+2)y=0 [xcosx,xsinx] この問題はまずどのようにして解いていけばよいかわかりません、参考書をみても。 おしえてください。

  • 微分方程式の問題ですが・・・

    y´´-3y´+y=e^x cosx という微分方程式をy=e^x (Acosx+Bsinx)の形で求めよという問題ですが、同次方程式の解と特殊解の解を求めればいいと思うのですが、 特性方程式λ^2 -3λ+1=0で解きます。解の公式で解くとλ=3±√5/2という解がでたのですがあっているのでしょうか?もしあっているとしたら基本解は実数解になるのですが、y=e^x (Acosx+Bsinx)の形で求めよという問ですので基本階は共役複素数解にならないといけないですよね?僕はどこを間違えているのでしょうか?教えてください

  • 微分方程式の解き方

    こんにちは。 現在、問題集で微分方程式の勉強を行っているのですが、下記の問題の解き方が分からず困っています。 y'-y sinx = y^2e^cosx の一般解を求める。 一般解を求めたところ、c=x+esinx となったのですが正しいでしょうか? y'=tany + x secy 変数変換を行うのだと思うのですが、secyの計算方法が分かりません。 少しでもアドバイスを頂ければ幸いです。 よろしくお願い致します。

  • 微分方程式の一般解

    y'+(e^x)y=3e^x という微分方程式があるんですけど 一般解を求めたいんですけど 求めてみたら ∫(3e^x)・e^(e^x)dx+c)e^(-e^x)になりました これってどうやって解くのでしょうか?