微分方程式の平衡点の安定性とはどうやったら判別できるのでしょうか?
例えば、dx/dt=x(1-x)(1/2-x)という微分方程式については
どうやって解けばいいですか?
下のようなサイトを調べましたが、どうもよく分りません。
http://www4.pf-x.net/~arataka/ode/node7.html
No.1
- 回答日時:
まず微分方程式を解く。
dx/dt=x(1-x)(1/2-x)
=x(x-1)(x-!/2)
{1/x(x-1)(1/2-x)}dx=dt
分解する。
{2/x + 2/(x-1) ー4/(x-1/2)}dx=dt
積分する。
lと |が解かりにくいので
logを Logとする。
2Log|x| + 2Log|x-1| -4Log|x-1/2|=t +c
2乗=^2 とする。
Logx^2・(x-1)^2/(x-1/2)^4=t +c
Logをはずす。
x^2・(x-1)^2/(x-1/2)^4=A・e^t
A=e^c
原式を尊重する。
x^2・(1-x)^2/(1/2-x)^4=A・e^t
1/2乗する。
|x|・|1-x|/(1/2-x)^2=B・e^(1/2)・t
B=√A
さてここから本論。
f(p)=0
f(p)=p(1-p)(1/2-p)=0
p=0
p=1
p=1/2
この3ッの平衡点 0 1 1/2
で定義に従って吟味する。
安定 漸近安定 不安定
でそれぞれ調べる。
δ ε での調査です。
あとできますか?
No.2
- 回答日時:
微分方程式の平衡点のことは参考書も枚挙にいとまがないし、サイトもあるので回答する必要はないのではないでしょうか。
Excelで数値実験をしてみると良いと思います。dx/dt=x(1-x)(1/2-x)の平衡点はx=0, 0.5, 1.0でオイラー法によればx(t+Δt) = x(t) + Δt*x(1-x)(1/2-x)
なので例えばΔt=0.1 としてx=0.5 近傍の解を調べるためにExcelのA1セルに0.51 A2のセルに数式=A1+A1*(1-A1)*(1/2 -A1) を入力します。そしてドラッグアンドコピーでこの数式をA3以下の列にコピーします。するとA1のセルを初期値としてオイラー法で解いた結果が得られます。こうして調べてみるとx=0.5は安定な平衡点、x=0, 1.0は不安定平衡点であることがわかります。y=x(1-x)(1/2-x)のグラフを描いてみると、平衡点で、dy/dxの正負が平衡点が安定になるかどうかに対応していることがわかります。数値実験してみてわかるもう一つの事は、この方法では平衡点への収束が遅いことです。平衡点ではdx/dt=0 なので収束が遅くなってしまうのです。そのため数値的に平衡点を求めるためには
http://front.math.ucdavis.edu/0707.1596
にあるような工夫が必要になると思われます
No.4
- 回答日時:
ygv123=No1
の補足連絡です。
No2・No3の回答でOKと思うのですけど。
やはり、
εーδ論法での回答がいるのかな?
数学ソフトが不具合なので、表示が大変です。
できれば、誰か助けて下さい。
時々覗き、なければ、
そのとき、微力ながら、書くこともあるかも?
ただ、わかります、その残りが????という
お困りが、、、
救い天使を私は待つ。
No.5ベストアンサー
- 回答日時:
適当に座標をずらしてx=0が平衡点になるようにしてあるとします。
リアプノフ関数をV(t)=(x(t))^2
で定義します。もしx=0の適当な近傍の中で、すべてのtについてdV/dt<0 であるならば、この近傍内の解はx=0に収束し、x=0が安定平衡点は明らかです(安定でないのはすべてのtについてdV/dt>0 である必要はない)。dx/dt=x(1-x)(1/2-x)のx=0の平衡点を調べてみると、この点の近傍で
dV/dt=2x・(dx/dt)>0
初期値が0より少し大きければ、dx/dt>0なのでオイラー法の次のステップではx=0からさらに大きくなっててしまう。0より少し小さければ、dx/dt<0なのでオイラー法の次のステップではx=0からさらに小さくなっててしまう。つまり安定ではありません。x=1/2を原点に持ってくるように平行移動して、x'=x-1/2と書き直し、x'の微分方程式に書き直して上記を適用すると、dV/dt<0すなわちx'=0は安定平衡点であることが示されるはずです。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 次の微分方程式を解け dx/dt=e^ax+b これがわかりません。詳しく説明して欲しいです 1 2023/05/22 12:35
- 数学 dx/dt = |y| , dy/dt = x (-∞<t<∞) をとけ 1 2022/09/17 09:56
- 数学 dx/dt=x-2y +e^t dy/dt=-3x +2y+1 初期値[1,0] [x,y] この連 3 2023/05/15 18:23
- 数学 dx(t)/dt =rx(t){1-(x(t)/K)} r,Kは正の定数とすると、この微分方程式はラ 1 2022/08/11 16:25
- 物理学 量子力学 球面調和関数 導出 方位角成分 微分方程式の解 2 2022/07/02 13:40
- 数学 微分方程式の問題 1 2023/07/27 12:11
- 数学 x=r・cosθの2回微分 θ=ωtとすると? 5 2022/05/10 23:53
- 数学 微分積分の微分方程式についての問題がわからないです。 2 2022/07/18 17:44
- 数学 微分積分の極限についての問題がわからないです。 1 2023/01/08 13:57
- 数学 積分についてです。 水平投射の微分方程式で、 m×dvx/dt=-Cvx (Cは粘性抵抗) を積分す 1 2022/10/06 19:45
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
積分で1/x^2 はどうなるのでし...
-
e^-2xの積分
-
dy/dx=-x/yの意味が解りません
-
∫1/√x dx 積分せよ 教えて下さい
-
微積分 dの意味
-
自励系のリッカチ方程式 dy/dt ...
-
項の右端につく縦棒の意味を教...
-
積分 Xの-2乗を積分するとどう...
-
∫1/(x^2+1)^2 の不定積分がわ...
-
(1+sinx)/sin2xの積分
-
∫e^cos(x) dx の計算
-
1/X^2の積分ってlogX^2ですか?
-
微分の記号
-
∮a^xdxこれを公式的に導いてほ...
-
dy=dy/dx・dxの求め方
-
積分です
-
写真の赤丸のようになぜ、(d²y/...
-
(x^2-y^2)dy/dx -2xy=0 この微...
-
dy/dx・dxは置換積分を使ってdy?
-
不定積分
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
積分で1/x^2 はどうなるのでし...
-
e^-2xの積分
-
フーリエ級数の問題で、f(x)は...
-
∫1/(x^2+1)^2 の不定積分がわ...
-
積分 Xの-2乗を積分するとどう...
-
∫1/√x dx 積分せよ 教えて下さい
-
項の右端につく縦棒の意味を教...
-
1/X^2の積分ってlogX^2ですか?
-
何が違いますか?
-
微積分 dの意味
-
∫e^cos(x) dx の計算
-
【数学Ⅱ・Ⅲ】微分の問題
-
2次微分の変数変換
-
x/(a^2+x^2)の積分について
-
exp(-ax^2)*cosx の証明
-
∫r/(a^2+r^2)^3/2drの計算の解...
-
x^2 * exp(x^2) dxの不定積分
-
x−1分の2の微分の仕方を教えて...
-
フーリエ変換の問題について
-
台形の任意の高さにおける上辺...
おすすめ情報