位置を求めるための微分方程式
- 一回微分と二回微分の式から位置を求める方法について説明します。
- 物理の二次元での空気抵抗がある問題では、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)
- 数学・算数
- 回答数3
- ありがとう数3
- みんなの回答 (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 のときと一緒ですね。
その他の回答 (2)
- Tacosan
- ベストアンサー率23% (3656/15482)
おっと, y を忘れた. md^2y/dt^2 = -mg-mk(dy/dt) じゃないとおかしい. で, これは dy/dt に関して変数分離形.
お礼
ありがとうございます。
- Tacosan
- ベストアンサー率23% (3656/15482)
確認ですが, 式は md^2/dt^2 = -mg-mky(dy/dt) であってますか? ひょっとして md^2/dt^2 = -mg-mk(dy/dt) だったりしませんか?
お礼
md^2/dt^2 = -mg-mk(dy/dt) です。ご指摘ありがとうございます。ご指導よろしくお願い申し上げます。
補足
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)}| やっぱり分かりません。教えてください。お願いします。
- ベストアンサー
- 数学・算数
- 物理の微分方程式
物理の微分方程式 高二です。塾で微分方程式を習ったのですが、さっぱりです。。。。 問 質点が速度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にはなりますけどそんなことしたらだめですよね? 助けてください。
- ベストアンサー
- 数学・算数
- ラグランジュの方法での位置を微分
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というふうになることは変数微分法で理解できているのですが上記の答えを導く解法がわかりません。 ステップバイステップで教えていただければ幸いです。
- ベストアンサー
- 数学・算数
お礼
ありがとうございます。ちゃんと微分の式で表記して表せるようこころがけます。