位置を求めるための微分方程式

このQ&Aのポイント
  • 一回微分と二回微分の式から位置を求める方法について説明します。
  • 物理の二次元での空気抵抗がある問題では、x軸では m・d^2/dt^2= -mk(dx/dt)、y軸では -mg-mky(dy/dt)という式が立てられます。
  • 一回微分の式を積分し、初期条件を用いて位置を求めることができます。ただし、y軸の積分方法については説明されていないサイトや資料があるようです。
回答を見る
  • ベストアンサー

一回微分と二回微分の式から位置を求める

物理の二次元での空気抵抗がある問題が出されて式は x軸では m・d^2/dt^2= -mk(dx/dt) y軸では -mg-mky(dy/dt) という式は立てられました。 しかし恥ずかしいことにちゃんと積分をした結果を出すことができません。 m,g,kはそれぞれ定数だとすると x軸は d^2x/dt^2 = -k dx/dt これを一回積分すると dx/dt = Ce^-kt で初期条件より t=0 v=V_xより v = V_x e^-kt そしてこれをさらに積分して x = (-V_x/k) e^-kt +C Cは初期条件よりV_x/k よって x = V_x/k (1-e^-kt) と出せたのはいいのですが yについてがまったくもって導出の仕方がわかりません。 答えは α = -g-k(dy/dt) v = -g/k + (V_y + g/k)e^-kt y = (-g/k)t + (1/k)(V_y + g/k)(1-e^-kt) となっていました。 このyの積分の方法を導出までの計算方法をちゃんと丁寧に説明しているサイトなどがあればお教えいただきたく存じます。 また大変ご面倒をおかけしますが可能ならばこちらにy軸における導出の計算過程をお教えいただければ幸いです。

  • ligase
  • お礼率92% (997/1082)

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

  • ベストアンサー
  • alice_44
  • ベストアンサー率44% (2109/4759)
回答No.3

m(d^2y/dt^2) = -mg - mk(dy/dt) なんですね? α = とか v = とか、変な略し方をするから、式変形を見失うのでは? v = dy/dt と置くならば、dv/dt = - g - kv です。 w = g + kv と置き換えると、(dw/dt)/k = -w で、 x の方程式が解けたのなら、これも w = C e^(-kt) と解けるでしょう。 v の式に戻せば、v = (1/k){C e^(-kt) - g} です。 もし、t = 0 のとき dy/dt = V_y であれば、C = g + k V_y。 dy/dt = -g/k + (g/k + V_y)e^(-kt) となるので、これを積分すれば、 y = Y_0 - (g/k)t + (g/k + V_y)(1/k){1 - e^(-kt)} です。 もし、t = 0 のとき y = 0 であれば、Y_0 = 0。 後半の処理も、x のときと一緒ですね。

ligase
質問者

お礼

ありがとうございます。ちゃんと微分の式で表記して表せるようこころがけます。

その他の回答 (2)

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

おっと, y を忘れた. md^2y/dt^2 = -mg-mk(dy/dt) じゃないとおかしい. で, これは dy/dt に関して変数分離形.

ligase
質問者

お礼

ありがとうございます。

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

確認ですが, 式は md^2/dt^2 = -mg-mky(dy/dt) であってますか? ひょっとして md^2/dt^2 = -mg-mk(dy/dt) だったりしませんか?

ligase
質問者

お礼

md^2/dt^2 = -mg-mk(dy/dt) です。ご指摘ありがとうございます。ご指導よろしくお願い申し上げます。

ligase
質問者

補足

md^2/dt^2 = -mg-mk(dy/dt) です。すみません。おっしゃるとおりで入力間違いです。 md^2/dt^2 = -mg-mk(dy/dt) でお願い申し上げます。

関連するQ&A

  • 微分方程式で式の変形 空気抵抗を受ける物体の落下

    質量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)}| やっぱり分かりません。教えてください。お願いします。

  • 微分を含む式の積分

    以下の式の積分の方法がわかりません。 ∫{c^2 - (dx/dt)^2}^(1/2) dt 仮に dx/dt=y と置くと、 ∫{c^2 - y^2}^(1/2) dt/dy dy=∫{c^2 - y^2}^(1/2) dt^2/dx dy =∫{c^2 - y^2}^(1/2) dx/y^2 dy となって、わけがわかりません。 どなたかアドバイスお願いします。

  • 物理の微分方程式

    物理の微分方程式 高二です。塾で微分方程式を習ったのですが、さっぱりです。。。。 問 質点が速度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)の答えを教えて下さい とても困っています!部分的でもよいので教えて下さい、お願いします

  • 微分方程式、初期値で微分すると?

    u=u(x,y), v=v(x,y)は(x,y)∈R^2のC^1級関数とし、 (ξ、η)∈R^2を初期値とする初期値問題        dx/dt=u(x,y), x(0)=ξ        dy/dt=v(x,y), y(0)=η の解を x=x(t;ξ,η), y=y(t;ξ,η) とする。 ∂x/∂ξ|[t=0]=1、 ∂/∂t(∂x/∂ξ)|[t=0]=∂u/∂x(ξ,η) を示せ。 って問題なんですけど、初期値(ξ)で偏微分?何のことか全然わかりません。xにt=0代入してからξで微分すれば1にはなりますけどそんなことしたらだめですよね? 助けてください。

  • 微分方程式の問題

    y^2dx+(xy-1)dy=0  解 xy=logy+c の解き方なのですが、同次形で解いていったところ . y=xv (dy/dx=v+xdv/dx)とする。  y^2dx+(xy-1)dy=0 ⇔y^2+(xy-1)dy/dx=0  ←両辺にdxをかける ⇔(y/x)^2+(y/x-1/x^2)dy/dx=0 ←両辺に1/x^2をかける ⇔v^2+(v-1/x^2)(v+xdv/dx)=0  ←y=xv,dy/dxを代入 ここからxとvについてうまくまとめられません。 どなたか教えてください。

  • ラグランジュの方法での位置を微分

    x,y,zがtの関数である位置ベクトル↑rは r↑(x(t),y(t),z(t)), r↑(x,y,z)と書くことができるので 時刻tで微分すると lim(dt→0){r↑(x+dx,y+dy,z+dz)-r↑(x,y,z)}/dt =dr↑/dt=v↑ =(dx/dt,dy/dt,dz/dt) となり速度が導かれますが、 ある流体粒子の位置r↑が ある位置(a,d,c)と任意の時刻tで決まるような関数つまりr↑(a,b,c,t)となる場合 速度は(a,b,c)を固定して偏微分で ∂r↑/dt =lim(dt→0) {r↑(a,b,c,t+dt)-r↑(a,b,c,t)}/dt =∂r↑/dt=v↑ (∂x/dt,∂y/dt,∂z/dt) となるのですか? ラグランジュの方法 https://hitopedia.net/%E3%83%A9%E3%82%B0%E3%83%A9%E3%83%B3%E3%82%B8%E3%83%A5%E3%81%AE%E6%9

  • パラメーター表示による微分の問題

    問題:曲線C:x=cos(t)+tsin(t)、y=sin(t)-tcos(t)に対して、 (1) dy/dxをtの式で表せ。 (2) d^2y/dx^2をtの式で表せ。 (1)はわかるのですが(2)でわからないところがあります。 多分正答はdx/dt、dy/dtを求めたあとdy/dxを求め d^2y/dx^2=d/dt(dy/dx)(dt/dx)という形で解くのだとおもいますが、 d^2y/dt^2を求めたあと d^2y/dx^2=(d^2y/dt^2)(dt/dx)^2と変形して解くと違う解がでてきます。 どうして違う解がでてくろのでしょうか? (dx)^2≠dx^2だからですか? でも d^2y/dx^2=d/dx(dy/dx)という変形ができるから、(dx)^2=dx^2だと思います。 だれかアドバイスください。 後、友達が ∫d(dy/dx)^2=∫(dy/dx)d(dy/dx)=(1/2)(dy/dx)^2+C という積分をしていたのですが、こういう積分もありなんですか? (Cは積分定数)

  • 二階の全微分について

    物理でxyの座標を極座標に変換し加速度を計算するなかで、2階の全微分に困っています。あまり、微分積分は慣れていないので、丁寧に教えていただけると助かります。 http://okwave.jp/qa/q2707943.html でも、同じような質問があります。 一階の全微分はわかりますが、2階の全微分で項が増えるのがわかりません。 具体的には、 Z=f(X,Y), X=g(t) Y=h(t)で、 dZ/dt=(∂Z/∂x)dx/dt+(∂Z/∂y)dy/dt まではよくわかり、これを二階にするときはまず、第1項目(∂Z/∂x)dx/dtが {∂/∂x(∂Z/∂x)dx/dt}dx/dt+{∂/∂y(∂Z/∂x)dx/dt}dy/dt となるだと思うのですが、(∂Z/∂x)d/dt(dx/dt)という項も加わるようです。詳しくその考え方を教えていただけますでしょうか?

  • 二次元の場合の速度の計算の解法

    y・・(yのツードット) ドットは一回微分を示す。つまりこの表記は二回微分とする。 初期条件を(t=0, y=0, y・=V_y0とする) y・・= -g-ky・  (kは正の定数)gは重力加速度 としたとき、これを一回積分すると y・= -g/k + (V_y0 + g/k ) e^-kt となると書いてありましたがなぜそうなるのかが分かりません。 dy/dt = ky が y = Ce^ktというふうになることは変数微分法で理解できているのですが上記の答えを導く解法がわかりません。 ステップバイステップで教えていただければ幸いです。

  • 数||| 微分

    x・dx/dt+y・dy/dt=0の両辺を さらにtで微分したときの式がよく分かりません。 x・dx/dt+y・dy/dt=0 (d^2x/dt^2)+(d^2y/dt^2)=0 d/dt(dx/dt)+d/dt(dy/dt)=0 d/dt(dx/dt+dy/dt)=0 自分でやってみましたが 全く自信がありません...。 分かる方いましたらよろしくお願いします。