• ベストアンサー

微分方程式の平衡点の安定性

微分方程式の平衡点の安定性とはどうやったら判別できるのでしょうか？ 例えば、dx/dt＝x(1-x)(1/2-x)という微分方程式については どうやって解けばいいですか？ 下のようなサイトを調べましたが、どうもよく分りません。 http://www4.pf-x.net/~arataka/ode/node7.html

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

  • ベストアンサー
回答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と書き直し、ｘ'の微分方程式に書き直して上記を適用すると、dV/dt<0すなわちx'=0は安定平衡点であることが示されるはずです。

その他の回答 (5)

  • ygv123
  • ベストアンサー率60% (3/5)
回答No.6

ｙｇｖ１２３の者です。（ＡＮｏ１） ＡＮｏ５の回答で解決です。 非常にすっきりしてます。 簡潔に納得できます。 実に簡単、綺麗です。

  • ygv123
  • ベストアンサー率60% (3/5)
回答No.4

ｙｇｖ１２３＝Ｎｏ１ の補足連絡です。 Ｎｏ２・Ｎｏ３の回答でＯＫと思うのですけど。 やはり、 εーδ論法での回答がいるのかな？ 数学ソフトが不具合なので、表示が大変です。 できれば、誰か助けて下さい。 時々覗き、なければ、 そのとき、微力ながら、書くこともあるかも？ ただ、わかります、その残りが？？？？という お困りが、、、 救い天使を私は待つ。

回答No.3

下の回答で、Δt=0.1 としているので、「ExcelのA1セルに0.51、 A2のセルに数式=A1+0.1*A1*(1-A1)*(1/2 -A1) を入力」に訂正します

回答No.2

微分方程式の平衡点のことは参考書も枚挙にいとまがないし、サイトもあるので回答する必要はないのではないでしょうか。Excelで数値実験をしてみると良いと思います。dx/dt＝x(1-x)(1/2-x)の平衡点はx=0, 0.5, 1.0でオイラー法によれば 　x(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 にあるような工夫が必要になると思われます

  • ygv123
  • ベストアンサー率60% (3/5)
回答No.1

まず微分方程式を解く。 ｄｘ／ｄｔ＝ｘ（１－ｘ）（１／２－ｘ） ＝ｘ（ｘ－１）（ｘ－！／２） ｛１／ｘ（ｘ－１）（１／２－ｘ）｝ｄｘ＝ｄｔ 分解する。 ｛２／ｘ　＋　２／（ｘ－１）　ー４／（ｘ－１／２）｝ｄｘ＝ｄｔ 積分する。 ｌと　｜が解かりにくいので ｌｏｇを　Ｌｏｇとする。 ２Ｌｏｇ｜ｘ｜　＋　２Ｌｏｇ｜ｘ－１｜　－４Ｌｏｇ｜ｘ－１／２｜＝ｔ　＋ｃ ２乗＝＾２　　とする。 Ｌｏｇｘ＾２・（ｘ－１）＾２／（ｘ－１／２）＾４＝ｔ　＋ｃ Ｌｏｇをはずす。 ｘ＾２・（ｘ－１）＾２／（ｘ－１／２）＾４＝Ａ・ｅ＾ｔ Ａ＝ｅ＾ｃ 原式を尊重する。 ｘ＾２・（１－ｘ）＾２／（１／２－ｘ）＾４＝Ａ・ｅ＾ｔ １／２乗する。 ｜ｘ｜・｜１－ｘ｜／（１／２－ｘ）＾２＝Ｂ・ｅ＾（１／２）・ｔ Ｂ＝√Ａ さてここから本論。 ｆ（ｐ）＝０ ｆ（ｐ）＝ｐ（１－ｐ）（１／２－ｐ）＝０ ｐ＝０ ｐ＝１ ｐ＝１／２ この３ッの平衡点　０　　１　　　１／２ で定義に従って吟味する。 安定　　漸近安定　　不安定 でそれぞれ調べる。 δ　ε　での調査です。 あとできますか？ 　

buchar
質問者

補足

回答ありがとうございます。 でもすいません、残りのδ、εでの調査がわからないです・・。

  1. カテゴリ
  2. 学問・教育
  3. 数学・算数

関連するQ&A

  • 微分方程式での平衡点

    微分方程式でx>0　y>0における平衡点を求める問題なのですが dx/dt＝x^2+2y^2-9　という問題で 単純な問題とですぐ解けるといわれましたが、 2変数になると良く分かりません、単純にdx/dtが0になる点のことでいいのでしょうか？　その場合x=1、y=2が答えになると思います。 深く考えすぎかもしれませんがアドバイスお願いします

  • 振り子の問題について　（平衡点など

    d^2x/dt^2 + a*dx/dt +sin(x)=0 この問題は非線形微分方程式です dx/dt=y dy/dt=-a*y -sin(x) とおいて連立させて平衡点を求めれば （0,0）（0,π）？になるとお思います この微分方程式を平衡点のまわりで線形化させたいのですが解法が全然分かりません アドバイスをいただけるか参考サイトを教えていただけけないでしょうか

  • 微分方程式

    こんにちは。微分方程式の問題が解けなくて困っています。 次のx(t)に関する微分方程式 d^2x/dt^2=-1/x^2 ただし初期条件はt=0でx=X0(x0>0）,dx/dt=√２であるとする。 問題　与式の両辺にdx/dtを乗じて積分することにより、初期条件を満たすxについての1階微分方程式をもとめよ。 必要ならば、公式d/dt(dx/dt)^2=2*(dx/dt)*(d^2x/dt^2) 少し問題の書き方がおかしいかもしれませんが（微分の書き方）どなたかお願いします。 自分なりにといたのですが 与式の両辺にdx/dtをかけて dx/dt(d^2x/dt^2)=-1/x^2*(dx/dt） 与えられた公式をつかい (1/2)*d/dt*(dx/dt)^2=-dx/dt*(1/x^2) ∫(1/2)*d/dt*(dx/dt)^2=-∫dx/dt*(1/x^2)　????? と与えられたヒント通りにしてそこからどうしたらいいのかわからなくなってしまいました・・・

  • 微分方程式

    こんにちは。微分方程式の問題が解けなくて困っています。 次のx(t)に関する微分方程式 d^2x/dt^2=-1/x^2 ただし初期条件はt=0でx=X0(x0>0）,dx/dt=√２であるとする。 (1)　与式の両辺にdx/dtを乗じて積分することにより、初期条件を満たすxについての1階微分方程式をもとめよ。 必要ならば、公式d/dt(dx/dt)^2=2*(dx/dt)*(d^2x/dt^2) (2)0<x0<1のときt(t≧0)餓変化した場合のx(t)の最大値を求めよ。 (1)は与式の両辺にdx/dtをかけて dx/dt(d^2x/dt^2)=-1/x^2*(dx/dt） 与えられた公式をつかい (1/2)*d/dt*(dx/dt)^2=-dx/dt*(1/x^2) (1/2)*d/dx*(dx/dt)^2=-(1/x^2) 両辺xで積分すると (dx/dt)^2=2/x+2(1-1/X0)(初期条件より） (2) は dt/dxが0すなわち1/xが-(1-1/X0）のときかとおもったのですが よくわからないです。 どなたかおねがいします。。

  • 微分方程式について

    微分方程式について。 yやdy/dxの形ならば解けるのですが ちょっと変わった形になると解けずに困っております。 回答お願いします。 1 未知関数x(t),y(t)に関する微分方程式 x´(t)=y(t), y´(t)=-x(t)を 初期条件x(0)=a, y(0)=bの下で解け。 2 x=x(t)を変数tのC^∞級関数とする。 このとき、 d^2x/dt^2 +(dx/dt)^2 -4=0 を解け。 3 tの関数x(t)が次の微分方程式を満たすとする x´+x^2+a(t)x+b(t)=0 ただしx´=dx/dtである。 ・x(t)=u´(t)/u(t)のとき、関数u(t)の満たす微分方程式を求めよ。 ・微分方程式 x´=x(1-x)の一般解を求めよ。 長いですが回答お願いします

  • 不動点の安定性

    微分方程式 dx/dt=x(1-x)(1/2-x)を前進オイラー法で解く場合のすべての 不動点の安定性について調べよ。という問題なのですが、誰か 解法を教えて下さい。よろしくお願いします。

  • 非線形微分方程式の特異点

    次の連立微分方程式について回答お願いします dx/dt = -x+y dx/dt = 2x+1-e^y 1,特異点を求める 2,特異点の近傍でe^yを一次の項まで近似し、方程式を線形化 3,線形化した方程式をベクトル表示 4,特性方程式（固有方程式）を解き、固有値を求める 5,この連立方程式の特異点のタイプは何か

  • 微分方程式の解き方

    dx/dt = x - (x + y)(x^2+y^2)^(1/2) dy/dt = y - (x - y)(x^2+y^2)^(1/2) という微分方程式があります。 この方程式の解を厳密に求めることはできないようですが、 (x^2+y^2)^(1/2) = r x = r cosθ y = r sinθ と置くことにより、上記の微分方程式の答えが、 dθ/dt = r dr/dt = r(1-r) を満たすことが分かるそうです。 ところで、上の微分方程式からどうやってこれを導くのでしょうか？勘でしょうか？

  • 微分方程式を解く問題が分かりません。

    微分方程式を解く問題が分かりません。 次の微分方程式が解けません。 {(d^2)x}/{d(t^2)}+2ε(dx/dt)+(ω^2)x=0 ただしε<ωとする。また初期条件をt=0でx=0、dx/dｔでv0とする。 が解けません。x=e^(αt)とおいて解いていくようなのですが・・・。 よろしくお願いします。

  • 微分方程式についてです

    dy/dx = -y の微分方程式で ,x=10tと置いた場合． dx/dt = 10 ∴dt/dx = 1/10 を使って， dy/dx = (dy/dt)(dt/dx) から， -y =(dy/dt)(1/10) ∴dy/dt = -10y とするのはいいのでしょうか． 1次近似なので，x=10t は分割が小さくなっただけのような気がするんですけど．

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する