• ベストアンサー

一階微分方程式

この微分方程式の解き方がわかりません。どなたかわかる人がいらしたら、教えてください。 Mdv(t)/dt=-ζv(t)+a*sin(ωt) 初速度をv(0)とおくと、この線形微分方程式の解は、 v(t)=(v(0)+(aω/M)/(ω^2+(ζ/M)^2)exp(-ζt/M)+(a/M)sin(ωt-δ)/√(ω^2+(ζ/M)^2) 公式どおり計算てみましたが、部分積分のところが上手に出来ません。その部分積分は、 v(t)=exp(-ζt/M)[a/M∫exp(ζt/M)*sin(ωt)dt+v(0)] のインテグラルの部分です。

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

  • ベストアンサー
  • qntmphscs
  • ベストアンサー率53% (14/26)
回答No.2

#1です。sin、cosのωが落ちてました。 I=(1/α)*exp(αt)*sin(ωt)-(ω/α)*∫exp(αt)*cos(ωt)dt =(1/α)*exp(αt)*sin(ωt)  -(ω/α)*{(1/α)*exp(αt)*cos(ωt)+(ω/α)*∫exp(αt)*sin(ωt)dt} =(1/α)*exp(αt)*sin(ωt)  -(ω/α^2)*exp(αt)*cos(ωt)-(ω^2/α^2)*∫exp(αt)*sin(ωt)dt =(1/α)*exp(αt)*sin(ωt)  -(ω/α^2)*exp(αt)*cos(ωt)-(ω^2/α^2)*I

その他の回答 (1)

  • qntmphscs
  • ベストアンサー率53% (14/26)
回答No.1

I=∫exp(αt)*sin(ωt)dtと置いて部分積分を2回実行します。 I=(1/α)*exp(αt)*sin(ωt)-(1/α)*∫exp(αt)*cos(ωt)dt =(1/α)*exp(αt)*sin(ωt)  -(1/α)*{(1/α)*exp(αt)*cos(ωt)+(1/α)*∫exp(αt)*sin(ωt)dt} =(1/α)*exp(αt)*sin(ωt)  -(1/α^2)*exp(αt)*cos(ωt)-(1/α^2)*∫exp(αt)*sin(ωt)dt =(1/α)*exp(αt)*sin(ωt)  -(1/α^2)*exp(αt)*cos(ωt)-(1/α^2)*I 両辺に現れたIについて解けば∫exp(αt)*sin(ωt)dtが出ます。 sinとcosのパートを合成すれば片付くように見えますが参考になりますか?

関連するQ&A

  • 微分方程式と積分

    1.次の微分方程式を解け。 (1)y''+2y'+y=3sin2x 同次微分方程式の一般解はu(x)=(C₁+C₂x)exp(-x) と求められるのですが、非同次微分方程式の特殊解u₀(x)が求められません。 どうやって求めればいいのでしょうか。 (2)y''-5y'+6y=x(exp(x)) 非同次微分方程式の特殊解u₀(x)はどうやって求めたらいいのでしょうか。 2.置換積分によって、次の定積分を求めよ。 1.∫[0→π/2] 1/(1+cosx)dx tanx/2=tと置いた後、どうすればいいのでしょうか。 2.∫[0→a] x^2(√a^2-x^2)dx(a>0) x=asintとおくと、dx=acost dt .∫[0→a] x^2(√a^2-x^2)dx=∫[0→π/2] a^2sin^2t*acos^2t dt このあとどうすればいいのでしょうか。 お願いします。

  • 微分方程式

    dx/dt=a^2-x^2 (aは実数の定数) (1)この微分方程式は1階の線形同次・線形非同次・非線形のどれにあてはまるか。 (2)この微分方程式の一般解を変数分離法で求めよ。 考えたことは(1)は非線形だと思いますが、合っていますか? (2)はdx/(x^2-a^2)=-dtと変形し、両辺積分します。  すると、1/(2a)log(|x-a|/|x+a|) = -t + C このあとx=が分からないです。 教えてください。お願いします

  • 微分方程式

    とき方が良くわからないのでわかりやすく回答していただけると嬉しいです。 (1)終端速度を求める問題です mdv/dt=mg-av dv/{(mg/a)-v}=a/m dt 上の式が下の式と同じなのはわかるのですが、左をvの式・右をtの式にするなら mdv/(mg-av^2)=dtで問題ない気がするのですが、なぜここで終わりにしなかったのでしょうか? (2)ばねの問題 m(dx)^2/dt^2 =-kx をとくには dsinwx=wcoswt 、 dcoswx=-wsinwx (wはオメガです) がヒントになると書いてあるのですがこういうのって普通にといていて思いつくものなのでしょうか?それとも一度どこかで聞いて理解し、機械的にこれを利用しているのでしょうか? (3)微分方程式が良くわからないのでとく際のpoint等あれば教えてください よろしくお願いします

  • 未定係数法は一階の線形微分方程式にも使えるのでしょうか? 

    未定係数法は一階の線形微分方程式にも使えるのでしょうか? 一階の線形微分方程式の解き方は dy/dt + p(t)y = g(t) のとき e^∫p(t)dt を両辺にかけて そのあとで両辺を積分してyについて解く と習いました。 そして、未定係数法は2階の線形微分方程式を解く方法の一つとして、 習いました。 ここで疑問に思ったのが、 この未定係数法は一階の線形微分方程式にも使えるのでしょうか? だとしたら下のような手順でよいのでしょうか? 同次式: dy/dt + p(t)y = 0 の一般解を求める (積分定数が残る) 非同次式: dy/dt + p(t)y = g(t) の特殊解を求める (積分定数はない) yの一般解 = 同次式の一般解 + 特殊解 よろしくお願いします。

  • 物理の微分方程式

    物理の微分方程式 高二です。塾で微分方程式を習ったのですが、さっぱりです。。。。 問 質点が速度Vに比例する抵抗力を受けて運動する際、V(t)、X(t)を求めよ。ただし、比例定数をk(>0)とする。 解 ma = mg-kv --1 a = dv/dt --2 v = dx/dt --3 1,2より dv/dt = g-kv/m よって dv/dt = -k/m(v-mg/k) ---4 ←変数分離型     (1/v-mg/k)dv = -k/m dt ----5 ここから積分して、計算して log{v(t)-mg/k} = C-kt/m ----6 (C=log{v(0)-mg/k}) {}は絶対値 そして {v(t)-mg/k} = e^C -e^-kt/m -----7 その後 v(t)=mg/k(1-e^-kt/m)(t≧0) となりました 質問 (1)初期条件ってなんですか?   (2)4→5の過程はなぜやるんですか?変数分離型ってなんですか?   (3)6→7の過程でなぜlogがとれるんですか?   (4)よければx(t)の答えを教えて下さい とても困っています!部分的でもよいので教えて下さい、お願いします

  • 微分回路における微分方程式の解

    微分回路では入力、出力電圧をそれぞれVi、Voとすると dVo/dt + Vo/RC = dVi/dt    (1) という微分方程式が成り立ちます。これを解くと Vo = Vi×exp(-t/RC) になるらしいのですが、(1)の微分方程式を解くことができません。 積分回路の微分方程式 dVo/dt + Vo/RC = Vi/RC は変数分離によって Vo = Vi(1 - exp(-t/RC)) という解が求まったのですが。 (1)の微分方程式の解き方を教えてください。よろしくお願いします。

  • 1階常微分方程式の問題が解けません

    大学の課題で出されたベルヌーイ型の微分方程式がどうしても解けません。 次のような問題です  dx/dt-2(t^2+1)x=-2x^2/t u=x^-1 とおいて同時線形微分方程式にすればいいのかと思ったのですが、 積分が難しくてとてもではないでけど解答の値がでてきませんでした。 解答の値は、  x=±(t^3/4+t/2+c/t)^(1/2) c:任意定数 となっています。 uの置き方をもっと工夫すべきなのか、単に積分計算が出来ていないだけなのか分かりません。 このように数式を書くのが初めてなので分かり辛い書き方ですみません。 お手数ですがよろしくおねがいします。

  • 微分方程式の解き方が分かりません

    この微分方程式が解けません。 dy/dt+y=g(t),y(0)=0のとき g(t)=2(0<=t<=1),0(t<1) 積分範囲で場合分けして解けると思うのですが、y(1)の値が分かれば解けると思います。ちなみに答えは、 y(t)=2(1-exp(-t));(0<=t<=1),2(e-1)exp(-t);(t>1) です。

  • 次の微分方程式を解いてください。

    以下の線形常微分方程式をどなたか解いてください。お願いします。 t=0の時x=0 dx/dt=(cos(t)+sin(t))x/(cos(t)-sin(t))+(cos(t)+sin(t))^2

  • 運動方程式です

    一様な静磁場B=(0.0.Bο)において、質量m、電荷qをもつ荷電粒子の運動方程式は次式で表される。以下の問題に答えよ。 mdv/dt=q(v×B) 運動方程式が次式で表される事を確かめよ。 v={A1sin(qBο/m)t+A2cos(qBο/m)t,A1cos(qBο/m)t-A2sin(qBο/m)t,A3} どのようにやったらこの式が出るのでしょうか。教えてください。