微分方程式の解法と疑問点
- 微分方程式dv/dt=-(k/m)vを解く過程で、疑問点があります。
- 疑問点１: e^(-C1)=0なのか、疑問点２: なぜvが正の値をとるのか
- 疑問点３: v0は何を表すか、わかりません。
- ベストアンサー
微分方程式 1/vと初期条件
質量mの質点が、速度に比例する抵抗のみの作用のもとで運動しているとき、微分方程式dv/dt=-(k/m)vが成り立つという。v(0)=0としてこの微分方程式を解け。という問題で3点わからない箇所があります。 自分の解き方では、-(k/m)=aとして、dt/dv=1/av 両辺をvで積分して、 ∫(dt/dv)dv=1/a∫(1/v)dv 、t=(1/a)(log|v|+C1) C1は積分定数、at-C1=log|v|、 |v|=e^(at-C1)、v=±e^(at)*e^(-C1) よってv(t)=±e^(at)*e^(-C1) 、t=0とv(0)=0を代入して 0=±e^(0)*e^(-C1)、0=±1*e^(-C1)よりe^(-C1)=0という矛盾がおこりました。これが１つ目のわからない点です。 二つの目のわからない点は、vが正の値をとるようになっていることです。 (1/a)∫(1/v)dv=(1/a)(logv+C1)のようになっていることです。 本の解答では、dt/dv=-(m/k)1/v、 -(k/m)dt/dv=1/vだから、 -(k/m)t=0+∫(v0→v)dv/v=log(v/v0)、よってv=v0e^{-(k/m)t}となっています。 v0=v(0)としたら、本の答えがv=0になってしまうので、それはないと思いますが、３つ目としてv0は何を表すかわかりません。 どなたか、なぜe^(-C1)=0になるのか、なぜvが正の値をとるか、v0は何かを教えてくださいおねがいします。
- situmonn9876
- お礼率91% (646/703)
- 数学・算数
- 回答数4
- ありがとう数4
- みんなの回答 （4）
- 専門家の回答
質問者が選んだベストアンサー
> 分母のvが０でないという条件だけでは、微分方程式を解けない そんなことは言っていない。 dt/dv=1/av と変形した段階で、「vで除算」してますが、「v(0) = 0」だったらその時「0で除算している」から、その段階で実は変でしょ？ということを言っています。だから「e^(-C1)=0という矛盾がおこ」る。 私も多分「v(0) = v0」の間違いだとは思っていますが、その場合でも v0 = 0だったら結局同じ問題点を議論しないといけません。
その他の回答 (3)
- tmppassenger
- ベストアンサー率76% (285/372)
> 不定積分でt=(1/a)(log|v|+C1) C1は積分定数 としたあと、 > v(0)もvも正、v(0)もvも負と場合分けして、絶対値記号をはずし、 > vを求めても間違ってはいないのでしょうか？ 勿論全く間違ってないですが、結局「積分定数はどうやって決めるか？」という事を考えた時、実は「最初から（積分範囲をきちんと考えた上で）定積分で考えれば自ら導出される」、という事です。
お礼
定積分で微分方程式を解く方法も、身に付けておきたくなりました。 お返事ありがとうございます。
- jcpmutura
- ベストアンサー率84% (311/366)
本当に問題に 「 v(0)=0としてこの微分方程式を解け 」 と書いてあるならば、問題が間違っています 解答が 「 v=v0e^{-(k/m)t} 」 ならば 問題は 「 v(0)=v0としてこの微分方程式を解け 」 でなければなりません vは正の値ではありません (v/v0)が正の値をとるようになっているのです
お礼
問題文の訂正ありがとうございます。すっきりしました。
- tmppassenger
- ベストアンサー率76% (285/372)
この手の微分方程式を解く時は本当はもっと慎重にならないといけないが、大抵そのように書いてない事が多い。 (1) 一般に0に対する除算は定義できないことを念頭に、 > -(k/m)=aとして、dt/dv=1/av と変形した時の、vが「分母」に来ていることの正当性をきちんと議論する事（つまり、分母が非零の時と、零になってしまうので除算出来ない場合をきちんと分け、その両者の解が「接続」することがあるかを議論しないといけない）。 (2) vが正の領域でも負の領域でも、∫(v0→v)dv/v=log(v/v0) でしょう？よく考えてみましょう。積分範囲が0を跨いではいけないことに注意。 (3) > -(k/m)dt/dv=1/vだから、-(k/m)t=0+∫(v0→v)dv/v=log(v/v0) と書いた時、dtの積分が t に化けていますが、この時dtをどの範囲で積分したのかを考えれば答えはでるはず。分からなければ「置換積分」とは何だったかを確認すること。 最後に、ここでも「vが分母にある」ことの議論をきちんとする必要があることに留意。
お礼
分母のvが０でないという条件だけでは、微分方程式を解けないとの指摘ありがとうございます。
関連するQ&A
- 物理の微分方程式
物理の微分方程式 高二です。塾で微分方程式を習ったのですが、さっぱりです。。。。 問 質点が速度Vに比例する抵抗力を受けて運動する際、V（ｔ）、X（ｔ）を求めよ。ただし、比例定数をｋ（>０）とする。 解 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≧０) となりました 質問 (1)初期条件ってなんですか？ (2)4→5の過程はなぜやるんですか？変数分離型ってなんですか？ (3)6→7の過程でなぜlogがとれるんですか？ (4)よければx(t)の答えを教えて下さい とても困っています！部分的でもよいので教えて下さい、お願いします
- ベストアンサー
- 物理学
- 微分方程式で式の変形 空気抵抗を受ける物体の落下
質量mの物体が速度の2乗に比例する空気の抵抗を受けながら落下する問題を考えよう。 鉛直上向きにy軸をとり、時間をtとすると、速度はdy/dtで表される。重力加速度の大きさをg, 抵抗力を係数をkとすると、運動方程式は次のようになる。 m (d^2y)/(dt^2) = -mg + k(dy/dt)^2 この方程式にはyが含まれていない。 速度を v = dy/dt とおけば、(d^2y)/(dt^2) = dv/dt であるから、運動方程式(2.19)は次のようにvについての1階の微分方程式に帰着される。 dv/dt = -g + k/m v^2 (2.20) この微分方程式は、次のように変数分離形の1階常微分方程式であり 1/ (v^2 - mg/k) dv/dt = k/m 両辺をtで積分すると 1/{2√(mg/k)} ∫[1/{(v-√(mg/k)} - 1/{(v+√(mg/k)}] dv = k/m ∫dt log |{v-√(mg/k)}/{v+√(mg/k)}| = 2√(kg/m)t + C (2.21) ここで、t=0 で v=0 として物体の落下だけを考えることにすると、y軸は鉛直上向きを正の方向としているので t>0 では v=dy/dt<0 dv/dt<0 となる。 したがって、式(2.20)から0<-v<√(mg/k)であることがわかり、式(2.21)からvは次のようになる。 v=dy/dt = -√(mg/k) * [1-e^{-2√(kg/m)t - C}] / [1+e^{-2√(kg/m)t - C}] (2.22) ・・・と本に書いてあるんですが、どうやってこの(2.22)を導き出したのかが分かりません。 勘でやってみますと、 log |{v-√(mg/k)}/{v+√(mg/k)}| = 2√(kg/m)t + C (2.21) の両辺でeをとって e^[log |{v-√(mg/k)}/{v+√(mg/k)}|] = e^{2√(kg/m)t + C} |{v-√(mg/k)}/{v+√(mg/k)}| = e^{2√(kg/m)t + C} |{v-√(mg/k)}| = e^{2√(kg/m)t + C} * |{v+√(mg/k)}| やっぱり分かりません。教えてください。お願いします。
- ベストアンサー
- 数学・算数
- 速度と微分方程式
速度と微分方程式 １．最初、２０ｍ/sの早さで走っていた自動車が一定の加速度で速さを増し、４．０秒後には４０ｍ/sになった。 （１）加速度一定の条件を用いて、早さに関する微分方程式を書け。 （２）初期条件を書け。 （３） （１）の微分方程式を（２）の初期条件を元に解け。 （４）速さの時間に関する式を求めよ。 （５） （４）を基に微分方程式を書け。これを解いて、この間に進んだ距離を求めよ。 (1) a=dv/dt ∫dv=a∫dt v+c=at+c v=at+c (2) 初期条件 t=0のときv=20m/s (3) 初期条件より v=5t+20 （３）までの答えはこれでいいでしょうか？ （４）と（５）がわからないのでどなたかよろしくお願いします。
- 締切済み
- 物理学
- 積分後m/k分合わない 空気抵抗を受ける物体の落下
本には v=dy/dt = -√(mg/k) * [1-e^{-2√(kg/m)t - C}] / [1+e^{-2√(kg/m)t - C}] (2.22) ここで、t=0 のときに v=0 であるから C=0 と定まる。 さらに、式(2.22)をtで積分するとyが次のように求まる。 y = -√(mg/k)t - m/k log( 1+e^2√(kg/m)t ) + C' (2.23) ・・・と書いてあります。 これを自力でやってみました。 「t=0 のときに v=0 であるから C=0 と定まる」とあるので、式(2.22)はCを消して実質 v=dy/dt = -√(mg/k) * [1-e^{-2√(kg/m)t}] / [1+e^{-2√(kg/m)t}] (2.22)' になります。これをtで積分すると、 ∫(dy/dt) dt = -√(mg/k) ∫[ [1-e^{-2√(kg/m)t}] / [1+e^{-2√(kg/m)t}] ] dt ∫dy = -√(mg/k) ∫[ [1+e^{-2√(kg/m)t} -2e^{-2√(kg/m)t}] / [1+e^{-2√(kg/m)t}] ] dt y = -√(mg/k) ∫[ 1 - [2e^{-2√(kg/m)t}] / [1+e^{-2√(kg/m)t}] ]dt y = -√(mg/k) [∫dt - 2∫[ [e^{-2√(kg/m)t}] / [1+e^{-2√(kg/m)t}] ]dt ] y = -√(mg/k)t + 2√(mg/k)∫[ [e^{-2√(kg/m)t}] / [1+e^{-2√(kg/m)t}] ]dt ここで分母を微分すると、[1+e^{-2√(kg/m)t}]' = -2√(kg/m)・e^{-2√(kg/m)t} ということで、後ろの項は、分子が-2√(kg/m)・e^{-2√(kg/m)t}であれば、積分すると晴れて log | [1+e^{-2√(kg/m)t}] |にできます。 ちょうど積分記号∫の前に+ 2√(mg/k)がありますので取り入れて y = -√(mg/k)t -∫[ [2√(mg/k)・e^{-2√(kg/m)t}] / [1+e^{-2√(kg/m)t}] ]dt y = -√(mg/k)t - log | 1+e^{-2√(kg/m)t} | + C' ・・・あれ？ 式(2.23)の m/k はどこから降ってきたのでしょうか？ どこか計算を抜かしていますでしょうか？ どうか教えてください。お願いします。
- ベストアンサー
- 数学・算数
- I = (V_i - V)/Rをtについて微分
共振回路の各素子の電流・電圧の式に対する式は以下のようになる。 V_i = RI + V ・・・ (1) I = C dV/dt + 1/L ∫ V dt ・・・ (2) 式(2)の両辺をtに関して微分すると、 dI/dt = C d^2V/dt^2 + V/L ・・・ (3) また、式(1)から、I = (V_i - V)/R だから、式(3)に代入すると、 1/R (dV_i/dt - dV/dt) = C d^2V/dt^2 + V/L ・・・この先も続くのですが、ここで I = (V_i - V)/R を代入すると、なぜこうなるのかが分かりません。 右辺は何も変わりはありませんね。 I = (V_i - V)/R を微分すると dI/dt = 1/R (dV_i/dt - dV/dt) になり、式(3)の左辺と同じなので代入してああなった、ということですか？ だとしたら、その微分の仕方が分かりません。 dV_i/dt と dV/dt はなんとなく分かるのですが 1/R って微分しても変わらないんでしたっけ？ (x^n)' = nx^(n-1) (1/R) = {R^(-1)}' = (-1)R^(-1-1) = -R^(-2) ??? あれ～、何故か定積分の式なら覚えてたりするんですけどね・・・。 tが文字の中に含まれてたら解けそうな気がするんですけど・・・あれあれ？ 記憶喪失してます。どうかお願いします。
- ベストアンサー
- 数学・算数
- 運動方程式の微分積分の計算
運動方程式の微分積分の計算方法がわかりません。詳しく教えてもらえると嬉しいです。よろしく、お願いします。以下はテキストの抜粋です。 m・dv/dt = F(r) 両辺に速度 v=dr/dt をかけると mv・dv/dt = F(r)・dr/dt となる。ここで、 v・dv/dt = d/dt(1/2v^2) ← この式変形が、分かりません。1/2も不明です。 と変形できるので、上の式は d/dt [1/2 mv^2(t)] = F・dr(t)/dt
- ベストアンサー
- 物理学
お礼
お礼コメントの間違いを指摘してくださりありがとうございます。 v(0)もvがとる値の一つと考えて、v(0)=0のとき0で割ってしまう異常になる というふうに考えました。
補足
話はそれますが、よかったらお返事ください。途中の計算のしかたについてですが、自分のように不定積分でt=(1/a)(log|v|+C1) C1は積分定数 としたあと、v(0)もvも正、v(0)もvも負と場合分けして、絶対値記号をはずし、vを求めても間違ってはいないのでしょうか？