• ベストアンサー

三元連立一階方程式系の一般解の分類について教えてください。

DX=AX という行列で表されたしきです。 演算子Dはtで一階微分することです。 Xは(x(t),y(t),z(t))の転置行列です。 A=(aij) これは行列Aの固有値を出して x'=λ1*x y'=λ2*y z'=λ3*z としてとくとよいのですか?

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

  • ベストアンサー
noname#108554
noname#108554
回答No.1

いいんじゃないでしょうか、だいたい。 解いた後、逆変換を忘れないでください。

from006
質問者

お礼

どうもありがとうございました。逆変換もしっかりとできました。

関連するQ&A

  • 連立2階微分方程式

    2階微分方程式の連立が解けません…(つд⊂;) できれば解き方のヒントをいただけると助かります。 yとzがxの関数で,k>0のとき (d^2y/dx^2)+k^2*z(x)=0 (d^2z/dx^2)+k^2*y(x)=0

  • 2階線形微分方程式の解法

    受験生です。 ずっと昔の大学入試問題です。答えがなく、悩んでいます。 d^2y/dx^2-(a+b)(dy/dx)+aby=0 (ただしdyやdxは微分演算子です) なのですが、おそらくa=bとa≠bとで分けるのだと思いますが、 両者ともどのようにして解けばよいのか分かりません。 解だけは載っていまして、 y=A*exp(ax)+B*exp(bx)とy=(Ax+B)*exp(ax) でした。 とりあえず私はa=bのときをやってみまして、 (d/dx-a)^2y=0と形式的に書き直して、 (d/dx-a)(dy/dx-ay)=0 とし、 (dy/dx-ay)=zとおいて ・(d/dx-a)z=0 ・(dy/dx-ay)=z を満たす解を探そうとしました。 上の方の式は直ぐにz=exp(ax)と出ましたが、 これを下の式に代入した後が分からなくなってしまいました。 勉強した範囲では、一階や二階の微分方程式の解は 一般解と特殊階の和で表せるということでしたが、それを元に 考えてみてもここから進みません。 質問は、以上の行き詰ってしまった所から先の解法と、 もうひとつの解であるy=A*exp(ax)+B*exp(bx)の導出方法です。 詳しい方、ご教授お願いできませんか。

  • 連立微分方程式 行列

    この問題がどうしても分からないので数学のできる方お願いします。 Q.次の連立微分方程式を解け   x'(t)= x(t) -y(t)+4z(t)   y'(t)=3x(t)+2y(t) -z(t)   z'(t)=2x(t) +y(t) -(t) この問題の前の問題で A=(1 -1 4) (3 2 -1) (2 1 -1)←行列 の固有値、固有ベクトル、Aの対角化、自然数nに対するA^n を求めたのですがここで使用するのでしょうか?

  • 連立微分方程式

    物理の勉強をしていたところ、その中の問題の一つを解く過程で連立微分方程式が出てきてしまいました。 解き方がいまいちわからないので、ご鞭撻の程よろしくお願いいたします。 以下に問題と自分なりの解答を記載します。 微分は全て’を用います。また、全て時間で微分しています。大文字は全て定数です。 また、時間微分演算子をtで表します。 例:x’=dx/dt   t=d/dt    t^2=d^2/dt^2 ------------------------------------------------------------- 問題  Ax''=-Bx+Cy'    -(1)    Ay''=-By-Cx'   -(2) --------------------------------------------------------------- 解答 (1)+(2):  A(x''+y'')+B(x+y)=C(y'-x') →   (At^2+B)(x+y)=Ct(y-x)    -(3) 同様に(1)-(2):  (At^2+B)(x-y)=Ct(y+x)    -(4) (3)×(y-x)-(4)×(x+y) :   0=Ct((y+x)^2+(y-x)^2) ∴ (y+x)^2+(y-x)^2=D 従って、y^2+x^2=E これより、この物理運動はx-y平面で円運動になる。 となったのですが、解き方が完全に思いつきの我流ですので、 間違ってるとしか思えません。(そもそも演算子をバラしていいのかどうか・・・。) 正しい解き方をご教授ください。 よろしくお願いいたします。

  • 微分方程式の一般解を求めたいです。

    dy/dx = (a+by)(c(x)+d(x)y) ここで、a,bは定数、c(x),d(x)はxの区間Iで連続とする。 (1)この微分方程式は、変数変換y = 1/b(1/z - a)により次の線形微分方程式に変換されるという。 dz/dx = f(x)z + g(x) をf(x),g(x)をa,b,c(x),d(x)を用いて表せ。 ********************************************* これはf(x) = ad(x) - bc(x) g(x) = -d(x) として答えがでました。 ********************************************* (2)a = b = 1,c(x) = x + 2/x , d(x) = xとするとき、微分方程式の一般解を求めよ。 dz/dx = -2z/x -x という式になると思うんですけど一般解をどう導き出していいのか分かりません。よろしくお願いします。

  • 1階線形偏微分方程式の一般解

    数学のことでちょっと皆様のお知恵を拝借いたしたく質問します。 次の偏微分方程式の一般解の求め方を教えてください。 ∂T(x,t)/∂t + (q(t)/S)(∂T(x,t)/∂x) = c(T_w(x,t) - T(x,t)) c,S:定数 僕の所有する参考書によるとこの種の方程式は ラグランジュの偏微分方程式と呼ばれていて、 ちょっとだけ一般解の求め方が書いてありました。 しかし、どうしても一般解にたどりつけません。 その方法とは、偏微分方程式 P(x,y,z)(∂z/∂x) + Q(x,y,z)(∂z/∂y) = R(x,y,z) に対して連立補助方程式 dx/P = dy/Q = dz/R を解いた解を f(x,y,z) = a, g(x,y,z) = b (a,bは積分定数) とする。φを任意の関数として、一般解は  φ(f,g) = 0 である。 という解法です。しかし、T_wが邪魔でうまくいかないです。 詳しい参考書を手に入れようにも近くに本屋がないのでお手上げです。 どなたかご教授お願いしますm( _ _ )m

  • 行列の連立微分方程式

    Aという行列に対して、dy/dt=Ay y(0)=(c0、c1、c2)^T【←転置】という連立微分方程式を解け。 ただしAは α 0 0 0 β 1 0 0 β という3×3の行列であるとする。 という問題なのですが、おそらくy=exp(tA)y(0)という形にして計算するのではないかと思うのですが、ここのexp(tA)を求める過程について質問です。 この問題を解く前に、 e^(αt) 0 0 0 e^(βt) te^(βt) 0 0 e^(βt) という3×3行列Xについて、 X^(-1)dX/dt=A という計算をさせられているので、おそらくこの式を利用して解くのだと思いますが、どのように利用したらいいのか見当がつきません。 exp(tA)=X になるのではないかという予想はあるのですが、どうしたらそうなるのか、過程等わかる方がいましたら、ご教授お願いします。

  • 三元連立微分方程式の解き方教えてください。

    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は円の運動をしました。

  • 連立微分方程式

    x=x(t), y=y(t)をtの関数として、次の連立微分方程式を考える。    dx/dt=2x+y dy/dt=x+2y (1)z=x+y, w=x-yとおいて、z,wについての微分方程式に書き換えなさい。 z+w=2x, z-w=2yつまり、x=(z+w)/2, y=(z-w)/2 これを、x,yにするだけでよいのでしょうか?

  • 微分方程式の一階化について質問です。

    微分方程式の一階化について質問です。 x''+ax'+a^2x=0という微分方程式('は時間微分)を一階化して行列表現せよ。 という問題に出会いました。 あくまでこれは問題の出だし部分でその後、固有値を求めたりなど色々続くのですが 正直、初めて見る問題で第一問目が一番分かりません。 数学に詳しい方、どうかよろしくお願い致します!!