- 締切済み
3元連立微分方程式の安定条件を教えて下さい
定数係数の2元連立微分方程式の安定条件は、ヤコビアン行列JのトレースtrJ<0, det(J)>0, (trJ)^2>4det(J)というのは知っています。では、3元連立微分方程式の場合、安定・不安定はどのように判定すればよいのか教えて頂けませんでしょうか。あるいは、参考文献または、URLでも結構です。
- nihonsumire
- お礼率92% (301/324)
- 経済学・経営学
- 回答数2
- ありがとう数13
- みんなの回答 （2）
- 専門家の回答
みんなの回答
- kiyos06
- ベストアンサー率82% (64/78)
>鞍点などの判定が出来るのでしょうか。 1)固有値が、実数部0で重根になると鞍点かも。 1.1)それに近い状態でも、鞍点もどきでしょう。 1.2)多分、近接根でも。 1.3)ただし、重根 -> 鞍点とは限らない。
- kiyos06
- ベストアンサー率82% (64/78)
1)2元の時は、固有値計算で、(λ -a11)(λ -a22) +a12 a21 =0の2根について実数部が負になればよい。 1.1)2実根と共役根の場合分けと、実数部確認で可能 1.2)定数が4つなので、detとtrで何とかなる。 2)3元の時は、3次方程式になって、3根の実部が負になることが必要 2.1)det,tr等だけでは多分無理 2.2)3次方程式の各係数を求めて、ラウス・フルビッツに持ち込む? 2.3)YahooやGoogleで「系 安定性 ラウス」「系 安定性 フルビッツ」を検索する。
お礼
ありがとうございます。やはり、ラウス・フルビッツ判定条件を使うのですね。現在、学習中です。ただ、鞍点などの判定が出来るのでしょうか。
関連するQ&A
- 三元連立微分方程式の解き方教えてください。
dx/dt=(A-C)/A*y*z dy/dt=((C-A)x-L)/A*z dz/dt=L/C*y この連立方程式なのですが，x,y,z以外は全て定数です。 まずは，ひとつの変数のみの方程式を作ろうと思いzだけの式を作りました。 z*z'''-z''*z'+((A-C)/A)^2*z^3*z'=0 という式になりました（'は微分回数です。） そして更に左辺は(z*z'')のtで微分した形に近いんじゃないかと思い解いていったのですが・・・・ここらへんからもうさっぱりでどうすればよいのかわからなくなってしまいました。もしよろしければこういう三元の連立で変数がyzみたいに積の項がある方程式のよい解き方を教えてください。パソコンで計算してみたらx,yは円の運動をしました。
- ベストアンサー
- 数学・算数
- 連立微分方程式について
連立微分方程式を勉強しているのですが、少しわからないところがあるので今回質問させて頂きました。 添付データの一行目の数式の連立微分方程式を解きたいのですが、 テキストの回答欄には 「二行目と三行目は任意のt,sに対して可換であるため、 四行目のようになる」 という説明書きがあるのですが、意味がよくわかりません…。 それがどう四行目につながっているのかも、意味不明です。。。 五行目では、exp(2t)をくくりだしているのはわかるのですが、 二個目の行列の中の-costが-cost/exp(2t)でないのはどうしてでしょうか？ かなり初歩的なところで躓いていると思いますので、お手数ですが詳しく教えて頂ければ嬉しいです。 よろしくおねがいしますm(_ _)m
- 締切済み
- 数学・算数
- 連立微分方程式の解き方について
Aの連立微分方程式 {y'₁(x)=4y₁(x)+2y₂(x) {y'₂(x)=2y₁(x)+3y₂(x) 初期条件 y₁(0) = 2, y₂(0) = 1 ↑ Aの連立微分方程式からまず行列を書くので |4 2| |2 3| ↑ と書きます。しかしここから固有値を求め、さらに固有値に対する固有ベクトル計算するとどうしても計算できなく、計算サイト(WolframAlpha)を使っても√を含む値が出てしまい、計算できなくなってしまいます。 問題を記載するにあたってタイピングミス等はありません。 誰か、わかる人教えてもらえないでしょうか？ マジで解けなくてほんとに困っています。
- 締切済み
- 数学・算数
- 連立微分方程式
点P(x,y)は連立微分方程式 dx/dt=y dy/dt=-x を満たすものとする。t=0で原点以外の点から出発した点P(x,y)は、tが増加するにつれてどのようにふるまうか述べよ。図を用いてもよい。 この問題の解き方がよく分かりません。 連立微分方程式について、色々な文献を見てみたのですが、どうもいまいちです。 上の連立方程式を2つともdt=のかたちにして、dx/y=dy/-xという式にし、変数を分離して両辺を積分して・・・すると、x^2+y^2=Cという式に なりました。 円の方程式っぽいです。 でも、tは消えてしまい・・・ よく分からなくなってきました。 そもそもここまでの解き方も自分は間違っているのでしょうか？？ ご意見やヒント、解答ヨロシクお願いしますm(_ _)m
- 締切済み
- 数学・算数
- 連立微分方程式の解き方について質問です
『１階線形微分方程式系 X'=2X-Y , Y'=5X-4Y を考える。 (１)行列Ａ＝（２ －１）（５ －４）とするとき、 Ａ＝aＩ+Ｂ 、 trB=０ をみたす定数a 及び 行列Ｂを求めよ。ただし、Ｉは二次単位行列である。 (２)exp(tA)を計算せよ。 (３）微分方程式の解を求めよ。』 (１)はできたのですが、(２)がわかりません……。 どうすれば、eを行列で累乗できるのしょうか。 (１)(２)と(３)のつながりもわかりません……。 よろしくお願いします。
- ベストアンサー
- 数学・算数
- 連立微分方程式の問題
x, u, v,を実数，a, τを実定数とする。次の連立微分方程式を解いてu, vを求めよ。 式は添付画像をご参照ください。 という問題です。 vを消去してuの微分方程式に書き換えたところ d4y/dx4(u)+(4a^2)(d4y/dx4)(u)=0 という式が得られてこの式を解くことができなくて... 私は間違っているのかそれとも別のやり方でやるべきですか。 この連立微分方程式の解き方をご存知の方がいらっしゃいましたら、ご指導お願いします。
- ベストアンサー
- 数学・算数
- 連立微分方程式の解き方がわかりません
連立微分方程式の解き方がわかりません ｆ(x,Pn)=φ(x,Pn)+φ(x, -Pn) と、ｘで偏微分した ｆ'(x,Pn)=φ'(x,Pn)+φ'(x, -Pn)=λｆ(x,Pn) λは、実定数 の連立微分方程式 です。 （要は、固有値問題） 整理すると、 φ'(x,Pn)+φ'(x, -Pn)=λφ(x,Pn)+λφ(x, -Pn) これを、以下と置ければ、解けるのですが、 φ'(x,Pn)=λφ(x,Pn) φ'(x, -Pn)=λφ(x, -Pn) そう置いて、いいものかどうか、わかりません。
- ベストアンサー
- 数学・算数
- 連立1次方程式の構成の問題について
3つの変数(a,b,c)を未知とする、連立1次方程式があり、さらに、(b,c,d)を未知数とするもう1つの連立1次方程式があります。未知数としてb,cの部分は重なっていますが、今のままでは２つの連立1次方程式は全く別ものとなっています。しかし、やはりb,cは共通だと考えた場合、今度は(a,b,c,d)を未知数とする4×4の連立方程式を構成して解くということになります。 つまり、2組の連立1次方程式(３元）から１つの4元の連立1次方程式を作るということになります。この場合、4元連立1次方程式を作る方法は唯一であるはずなのですが、どのように考えたらいいでしょうか。解き方としては逆行列などを作用させて...と考えます。変な行列を作ってしまたら（例えば1つの行が(0,0,0,0)とか）になると逆行列が作れず唯一の解が出ないと思います。でもちょっと考えたら（作り方をまちがえたら？）そうなってしまう可能性があります。このように２つの連立方程式からちょっと大きな別の連立方程式を作って逆行列で解く方法について教えて頂きたいのですが。 なお、3×3は行列式は非ゼロであり、至極無理のないものを考えています。変なマトリックス（係数が10^(-8)とか）は全く想定していません。 最終的にはプログラム化していくことを考えていますが、今はその前段階の考え方についてお尋ねします。 よろしくお願いします。
- 締切済み
- 数学・算数
お礼
ありがとうございます。浅学非才のため、kiyos06 さんのコメントは、理解できてません。もう少し学習をしてみます。