• 締切済み

あってますか??

d^2y/dx^2-dy/dx-6y=0,y(0)=1,y`(0)=0のとき、dy/dx=zとおくと、dz/dx=z+6yになる。yを求めよ。っていう問題で、 d/dx=Dっておいて(D^2-D-6)y=0 だから(D-3)(D+2)y=0 よって(D-3)y=0または(D+2)y=0 (D-3)y=0はdy/dx-3y=0でdy/dx=3y だからdy/y=3dx 両辺積分すると、logy=3x+C (Cは積分定数) y=e^(3x+C) M=e^Cとおくと y=Me^3x 同様にy=Ne^(-2x) あわせて、y=Me^3x+Ne^(-2x) y(0)=1,y`(0)=0より M=2/5,N=3/5 よってy=2/5e^3x+3/5e^(-2x) あってますか?? 変なところあったら教えてください>< 問題文の方法使ってないような気がするんですけど だめなんでしょうか?? 連立微分方程式を使うやり方がわかる人は 教えてください!!!

みんなの回答

  • takamm
  • ベストアンサー率40% (4/10)
回答No.2

  d^2y/dx^2-dy/dx-6y=0を変形して、   d/dx(dy/dx+αy)=β(dy/dx+αy) となるとします。これを展開すると、   d^2y/dx^2+d(α-β)y/dx-αβy=0 これより、問題の式と各項の係数を比較して、   α-β=-1   -αβ=-6 これを解いて、(α,β)=(2,3),(-3,-2) (i) (α,β)=(2,3)のとき、   d/dx(dy/dx+2y)=3(dy/dx+2y) ここで、dy/dx+2y=zとおきます   dz/dx=3z   z=Ae^(3x) これより、   dy/dx+2y=Ae^(3x) x=0のとき、y(0)=1,y'(0)=0だから、それぞれを代入して、A=2が得られる。よって、   dy/dx+2y=2e^(3x) ―(*) (ii)(α,β)=(-3,-2)のとき、(i)と同様にして、   dy/dx-3y=-3e^(-2x)  ―(**) (i),(ii)より、得られた2式より、   y=(2/5)e^(3x)+(3/5)e^(-2x) (∵(*)-(**)) 一応、微分方程式を連立で解いてますが、どうでしょうか?なぜ、dy/dx=zと置きたいのかは分かりませんが、上のようにおけば、解きやすいかと思います。dy/dx=zというのは、問題文に含まれていたのでしょうか?まず、そこが正しいか確認したほうがよいかと思います。dy/dx=zとおいて、解けないといってるわけではないのですが、自分としてはdy/dx=zとおいても解き方は思いつかないです。

  • Tacosan
  • ベストアンサー率23% (3656/15482)
回答No.1

「連立微分方程式を使う方法」が何を指しているかはしりませんが, 例えば dy/dx = z dz/dx = z+6y を行列で書くと (d/dx) (y, z) = ((0, 1), (6, 1)) (y, z) の形になる (ベクトル, 行列の形は雰囲気で察してください) から右辺の行列を A とおくと (y, z) = e^(Ax) (y(0), z(0)) とできるんですが.... いいのかなぁ, これで.

関連するQ&A

  • 偏微分方程式の解き方

    x(y-z) (∂z/∂x) + y(z-x) (∂z/∂y) = z(x-y) この微分方程式を解く問題で、解答を見ても、理解できない部分があるため、質問させていただきます。 ~解答~ 補助方程式 dx/x(y-z) = dy/y(z-x) = dz/z(x-y) これより、 (1/x)dx / (y-z) = (1/y)dy / (z-x) = (1/z)dz / (x-y) と変形できます。 ここまでは分かるのですが、 これに加比の理を適用すると、 d(logx+logy+logz) / 0 = ((1/x)dx+(1/y)dy+(1/z)dz) / ((y-z)+(z-x)+(x-y)) = (1/x)dx / (y-z) = (1/y)dy / (z-x) = (1/z)dz / (x-y) ↑ここの1つ目のイコールが何故、成り立つのかが理解できません。 d(logx+logy+logz)を計算したら、1/x + 1/y +1/z になってしまわないでしょうか? 逆に、積分してlogになったのだとしても、dが残る理由が理解できません。 よろしくお願いします。 一応、続きも書いておきます。 ここで、d(logx+logy+logz) / 0 より、 d(logx+logy+logz) / 0 = d(logxyz) = 0 よって、logxyz = C' ゆえに、xyz = C (積分定数) このあと、もう1つの解を出して、一般解とします。

  • 微分の変形

    条件として x=e^z dx/dz=e^z=x y'=dy/dx=dy/dz*dz/dx=1/x*dy/dz なのですが、 y"=d^2y/dx^2=d/dx*(1/x)*dy/dz+1/x*d/dx*(dy/dz) の変形の理由がわかりません。 どなたかご教授お願い申し上げます。

  • 完全形でない3変数関数の微分方程式の解法

    全微分方程式A(x,y,z)dx+B(x,y,z)dy+C(x,y,z)dz=0がある。この式をPとおく。ここで、ベクトル値関数f=[A,B,C]とおき、f・(rotf)=0となるならばPは積分可能でその一般解は下記の手順により求まる。 手順1:Pについてdz=0とすると、Adx+Bdy=0となる。この式をQとおく。これが(∂A/∂y)=(∂B/∂x)を満たすとき、また満たさないときは積分因子μをかけることによりこのQの一般解ξ(x,y,z)=E (Eは定数)が得られる。 手順2:Pの両辺にλをかけたものの一般解を求める。するとλAdx=(∂ξ/∂x)となる。これから、λの値を求める。 手順3:ξの全微分はdξ=(∂ξ/∂x)dx+(∂ξ/∂y)dy+(∂ξ/∂z)dzとなり、このうち(∂ξ/∂x)dx+(∂ξ/∂y)dyはλAdx+λBdyとなるが、最後の(∂ξ/∂z)dzだけはλRdzとなるかは不明である。 dξ=(∂ξ/∂x)dx+(∂ξ/∂y)dy+(∂ξ/∂z)dzと(∂ξ/∂x)dx+(∂ξ/∂y)dy=λAdx+λBdyより、λAdx+λBdy=dξ-(∂ξ/∂z)dzとなる。 するとPの両辺にλをかけた式は、λAdx+λBdy+λCdz=dξ+{λC-(∂ξ/∂z)}dz=0となる。 ここで、λC-(∂ξ/∂z)=ηとおくと、λAdx+λBdy+λCdz=dξ+ηdz=0となり、2変数の全微分方程式dξ+ηdz=0が得られる。この解が結局全微分方程式Pの一般解となる。 ここで質問です。 手順1でdz=0とした式Adx+Bdy=0 (∂A/∂y)=(∂B/∂x)、またはμAdx+μBdy=0 (∂μA/∂y)=(∂μB/∂x)を解くとこの一般解、ξ(x,y,z)=Eが得られ、この関数ξの全微分はdξ=(∂ξ/∂x)dx+(∂ξ/∂y)dy=Adx+Bdy=0、またはdξ=(∂ξ/∂x)dx+(∂ξ/∂y)dy=μAdx+μBdy=0になるのが分かります。 手順2,3でλAdx+λBdy+λCdz=0という式が出てきますが、これはλをかける事により完全形になっていると思われます。しかしなぜλAdx=(∂ξ/∂x)となるのかが分かりません。ξはAdx+Bdy=0の解として現れる関数なので、λAdx+λBdy+λCdz=0を満たす関数は別にあり、例えばこれをσとすると、この関数の全微分はdσ=(∂σ/∂x)dx+(∂σ/∂y)dy+(∂σ/∂z)dz=λAdx+λBdy+λCdz=0となり、λAdx=(∂σ/∂x)dxとなるのではないのでしょうか? それともこの関数σがξと一致すると仮定しているのでしょうか? それからもう1つ気になるのですが、手順3で「最後の(∂ξ/∂z)dzだけはλRdzとなるかは不明である。」とありますが、これもよく意味が分かりません。なぜ(∂ξ/∂z)dzだけλRdzとはなるか分からないのでしょうか? おそらく私が根本的に間違っていると思いますので、詳しい方教えてください。お願いします。

  • 対数微分法の過程

    y=a^xについて、両辺の自然対数をとると logy=xloga この両辺をxについて微分すると d/dy(logy)dy/dx=d/dx(xloga) 1/ydy/dx=loga となるようなのですが、両辺をxで微分するとどうしてそうなるのかわかりません。教えて下さい 特にlogyをyで微分して、それをyをxについて微分したものにかけているのがよくわかりません

  • 変数を置換える微分方程式について

    お世話になります、 以下の考え方で問題がないかご教授願います。 (どうして大学のテキストは解答がついてないのでしょうか、成否の確認ができません…) 問題:y´=y^2/(xy-x^2) 最初は変数分離で一瞬で解けると思ったのですが、分母がネックになりました。 右辺を1変数に置き換えられないか考えて、分子分母に1/(xy)を乗じます。 y^2/(xy-x^2)=(y/x)/(1-x/y) において、 z=(y/x)としてz^(-1)=(x/y),dz/dy=1/xよりdy=x dz ∴ z/(1-z^(-1))=(x・dz)/dx 両辺を逆数化して (1-z^(-1))/z=dx/(x・dz) ∫(1-z^(-1))/z dz=∫1/x dx    log|z|+z^(-1)+C=log|x| <Cは積分定数です> つまり、log x=log z+z^(-1)+C log x=log z+log e^z^(-1)+log e^C log x=log{ze^z^(-1)e^C} x=ze^z^(-1)C´ <e^C=C´とする> x=(y/x)e^(x/y)C´ 一応、微分は無くなったので正解となるのでしょうか? お世話になります。

  • 対数微分法について

    例えばy=sinx^xなどという関数は両辺自然対数をとりますよね そのとき、左辺はlogyとなり 「両辺xについて微分したとき」左辺はy'/yとなりますが 「xについて微分なのになぜyがxの関数かのように微分されているのですか?」 考えられたことは、logyを微分したら、d(logy)/dy×(dy/dx)でlogy/dxと同じことになるので、d(logy)/dyは1/yですよね。ということは・・・?dy/dxはy'ということでしょうか?けどyっていうのはxという文字を含んでいませんよね・・・。 合成関数みたいな感じでしょうか・・・?合成関数って微分したら中身をさらに微分するけど・・・ y'ってやるとyの中身は・・・? などと混乱してしまいました。 アドバイスお願いします。

  • 一階線形微分について

    y' - 2y/x = (x^2) cos3xについて、解きなさいという問題なのですが、 同次形の場合について解いてみて、 dy/dx - 2y/x = 0 1dx/y = 2dx/x 両辺を積分 logy = 2∫1dx/x + C = logx^2 + C y = x^2 + e^C e^C = C(x) y = x^2 + C(x) となったのですが、答えがe^□にならないといけないのでしょうか。非同次形を続けて解いてみたら可笑しな値になりました。いまいち一階線形微分について理解出来ていないので、ご指導頂けたら有り難いです。

  • ベクトル場と線積分

    電磁気の問題で大変困っています。 一様なベクトル場、E→=(Eo,0,0)があり、次の二つの異なる経路C1とC2について、点A(a,0,0)から点B(0,a,0)まで線積分せよ。 C1:X+Y=a、Z=0 C2:X^2+Y^2=a^2,Z=0(X,Y≧0) という問題です。 C1の場合、自分は、 ∫_c1E→・(dx,dy,dz)=∫_c1(Eo,0,0)・(dx,dy,dz)=∫Eodx(xはa→0) =[Eox]=-Eoa で一応問題の答えはこれであっているのですが、このような計算過程でいいのでしょうか? 次のC2では、自分は、 ∫_c2E→・(dx,dy,dz)=∫_c2 (Eo,0,0)・(dx,dy,dz) ここで、x=acosθ、y=asinθとおき、 dx=-asinθdθ、dy=acosθdθ θの変化は、0→π/2で、 上式=∫Eo・(-asinθdθ)=[Eoa・cosθ]=Eoa としました。 この問題の答えは、Eoaπ/2となると思うのですが。。。   確実にC2の場合は自分の計算は間違っていると思います。 どなたか説明して頂けませんか? お願いします。 同じくC1についても間違いのある点をご指摘して頂ければ幸いです。 お願いします。

  • dy/dx (y+1)を積分して(y+1)^2?

    次の微分方程式の一般解を求めよ。 (1+y) (d^2y)/(dx^2) + (dy/dx)^2 = 0 dy/dx = p とおくと、      (1+y)p (dp)/(dy) + p^2 = 0 となり、      (i) (1+y) (dp)/(dy) + p = 0      (ii) p = 0 の2通りが考えられる。 (i)の場合      1/p (dp)/(dy) + 1/(1+y) = 0 の両辺をyで積分して      log |p(y+1)| = C_1 つまり、      dy/dx (y+1) = C_1 両辺をxで積分して、      (y+1)^2 = C_1x + C_2     ←? という解を得る。 ・・・と本に書いてあります。しかし、 「両辺をxで積分して」の計算は間違ってないですか? 自分が計算すると、      dy/dx (y+1) = C_1      ∫ (y+1) dy/dx dx = C_1∫dx      ∫ (y+1) dy = C_1∫dx      ∫y dy + ∫1 dy = C_1∫dx      y^2/2 + y = C_1x + C_2 になります。 積分して(y+1)^2になるなら、元々は2(y+1)じゃないといけないですよね、きっと。 ということで、どなたか検算をお願いします。

  • ベッセルの方程式の問題の解き方が分かりません

     次のベッセルの方程式の問題の解き方が分かりません。  数学に詳しい方、よろしければご教示願えないでしょうか。 問題は、  ベッセルの方程式に帰着できるさまざまな方程式がある。示されている置換を 使って、次の微分方程式の一般解を求めよ。 4*x^2*y" + 4*x*y' + (x - ν^2)*y = 0 (√x = z)  このように解いてみました。  ベッセルの微分方程式は、 x^2*y" + x*y' + (x^2 - ν^2)*y = 0 で、  一般解は、 y(x) = A*Jν(x) + B*Yν(x) ここで、A と Bは任意定数、Jν(x)は第1種ベッセル関数、Yν(x)は第2種ベッセル 関数。 √x = z より、 dz/dx = 1 / (2*√x) y'とy"は、 y' = dy/dx = (dy/dz)*(dz/dx) = (dy/dz)/(2*√x) y" = d^2y/dx^2 = (d/dx)*(dy/dx) = (d/dz)/(2*√x)*(dy/dz)/(2*√x) = (d^2y/dz^2)/(4*x) ゆえに、 4*x^2*y" + 4*x*y' + (x - ν^2)*y = 4*x^2*(d^2y/dz^2)/(4*x) + 4*x*(dy/dz)/(2*√x) + (x - ν^2)*y = x*(d^2y/dz^2) + 2*√x*(dy/dz) + (x - ν^2)*y = z^2*(d^2y/dz^2) + 2*z*(dy/dz) + (z^2 - ν^2)*y = 0 となって、第 2項目が z*(dy/dz) にならず、2*z*(dy/dz) になってしまいます。  本の回答をみると、 A*Jν(√x) + B*Yν(√x) となっているので、問題の微分方程式を、 z^2*(d^2y/dz^2) + z*(dy/dz) + (z^2 - ν^2)*y = 0 に変形したのだと思いますが、どのようにすれば良いのでしょうか ?  同様に下記の問題も、 x^2*y" + x*y' + 4*(x^4 - ν^2)*y = 0 (x^2 = z) 同じ解き方をしたため、 z^2*(d^2y/dz^2) + z*(dy/dz) + (z^2 - ν^2)*y = 0 に変形できませんでした。  なにとぞよろしくお願いします。