与えられた連立微分方程式について、以下の問いに答えます。 (1) 連立微分方程式を $\frac{d\mathbf{x}}{dt} = A\mathbf{x}$ の形で表現する。 (2) 行列 $A$ の固有値と固有ベクトルを求める。 (3) 連立微分方程式の一般解を求める。

応用数学微分方程式線形代数固有値固有ベクトル連立微分方程式
2025/6/1

1. 問題の内容

与えられた連立微分方程式について、以下の問いに答えます。
(1) 連立微分方程式を dxdt=Ax\frac{d\mathbf{x}}{dt} = A\mathbf{x} の形で表現する。
(2) 行列 AA の固有値と固有ベクトルを求める。
(3) 連立微分方程式の一般解を求める。

2. 解き方の手順

(1) dxdt=Ax\frac{d\mathbf{x}}{dt} = A\mathbf{x} の形に表現する。
与えられた連立微分方程式は以下の通りです。
dx1dt=x1x2x3\frac{dx_1}{dt} = x_1 - x_2 - x_3
dx2dt=x1+x2x3\frac{dx_2}{dt} = x_1 + x_2 - x_3
dx3dt=2x1x2\frac{dx_3}{dt} = 2x_1 - x_2
したがって、AAx\mathbf{x} は以下のようになります。
x=(x1x2x3)\mathbf{x} = \begin{pmatrix} x_1 \\ x_2 \\ x_3 \end{pmatrix}
A=(111111210)A = \begin{pmatrix} 1 & -1 & -1 \\ 1 & 1 & -1 \\ 2 & -1 & 0 \end{pmatrix}
(2) 行列 AA の固有値と固有ベクトルを求める。
固有方程式は AλI=0|A - \lambda I| = 0 で与えられます。
AλI=1λ1111λ121λ=0|A - \lambda I| = \begin{vmatrix} 1-\lambda & -1 & -1 \\ 1 & 1-\lambda & -1 \\ 2 & -1 & -\lambda \end{vmatrix} = 0
(1λ)((1λ)(λ)1)+(1)(λ+2)1(12(1λ))=0(1-\lambda)((1-\lambda)(-\lambda) - 1) + (-1)(-\lambda + 2) - 1(-1 - 2(1-\lambda)) = 0
(1λ)(λ+λ21)+λ2+1+22λ=0(1-\lambda)(-\lambda+\lambda^2 - 1) + \lambda - 2 + 1 + 2 - 2\lambda = 0
λ+λ21+λ2λ3+λ+λ2+1+22λ=0-\lambda+\lambda^2-1+\lambda^2-\lambda^3+\lambda + \lambda - 2 + 1 + 2 - 2\lambda = 0
λ3+2λ22λ=0-\lambda^3 + 2\lambda^2 - 2\lambda = 0
λ(λ22λ+2)=0-\lambda(\lambda^2 - 2\lambda + 2) = 0
したがって、固有値は λ1=0,λ2=1+i,λ3=1i\lambda_1 = 0, \lambda_2 = 1 + i, \lambda_3 = 1 - i です。
固有ベクトルを求める。
λ1=0\lambda_1 = 0 のとき、(A0I)v=0(A - 0I)\mathbf{v} = 0
(111111210)(v1v2v3)=(000)\begin{pmatrix} 1 & -1 & -1 \\ 1 & 1 & -1 \\ 2 & -1 & 0 \end{pmatrix} \begin{pmatrix} v_1 \\ v_2 \\ v_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}
v1v2v3=0v_1 - v_2 - v_3 = 0
v1+v2v3=0v_1 + v_2 - v_3 = 0
2v1v2=02v_1 - v_2 = 0
v1=12v2v_1 = \frac{1}{2} v_2
12v2v2v3=0    12v2=v3\frac{1}{2} v_2 - v_2 - v_3 = 0 \implies -\frac{1}{2} v_2 = v_3
v1=(121)\mathbf{v}_1 = \begin{pmatrix} 1 \\ 2 \\ -1 \end{pmatrix}
λ2=1+i\lambda_2 = 1 + i のとき、(A(1+i)I)v=0(A - (1+i)I)\mathbf{v} = 0
(i111i1211i)(v1v2v3)=(000)\begin{pmatrix} -i & -1 & -1 \\ 1 & -i & -1 \\ 2 & -1 & -1-i \end{pmatrix} \begin{pmatrix} v_1 \\ v_2 \\ v_3 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}
iv1v2v3=0-iv_1 - v_2 - v_3 = 0
v1iv2v3=0v_1 - iv_2 - v_3 = 0
2v1v2(1+i)v3=02v_1 - v_2 - (1+i)v_3 = 0
第一式と第二式から iv1v2=v1iv2-iv_1 - v_2 = v_1 - iv_2
v2=1+i1iv1=iv1v_2 = \frac{1+i}{1-i} v_1 = i v_1
iv1iv1v3=0    v3=2iv1-iv_1 - iv_1 - v_3 = 0 \implies v_3 = -2iv_1
v2=(1i2i)\mathbf{v}_2 = \begin{pmatrix} 1 \\ i \\ -2i \end{pmatrix}
λ3=1i\lambda_3 = 1 - i のとき、v3=v2=(1i2i)\mathbf{v}_3 = \overline{\mathbf{v}_2} = \begin{pmatrix} 1 \\ -i \\ 2i \end{pmatrix}
(3) 連立微分方程式の一般解を求める。
一般解は x(t)=c1eλ1tv1+c2eλ2tv2+c3eλ3tv3\mathbf{x}(t) = c_1 e^{\lambda_1 t} \mathbf{v}_1 + c_2 e^{\lambda_2 t} \mathbf{v}_2 + c_3 e^{\lambda_3 t} \mathbf{v}_3 で与えられます。
x(t)=c1(121)+c2e(1+i)t(1i2i)+c3e(1i)t(1i2i)\mathbf{x}(t) = c_1 \begin{pmatrix} 1 \\ 2 \\ -1 \end{pmatrix} + c_2 e^{(1+i)t} \begin{pmatrix} 1 \\ i \\ -2i \end{pmatrix} + c_3 e^{(1-i)t} \begin{pmatrix} 1 \\ -i \\ 2i \end{pmatrix}
x(t)=c1(121)+c2et(cost+isint)(1i2i)+c3et(costisint)(1i2i)\mathbf{x}(t) = c_1 \begin{pmatrix} 1 \\ 2 \\ -1 \end{pmatrix} + c_2 e^{t}(\cos t + i\sin t) \begin{pmatrix} 1 \\ i \\ -2i \end{pmatrix} + c_3 e^{t}(\cos t - i\sin t) \begin{pmatrix} 1 \\ -i \\ 2i \end{pmatrix}
x1(t)=c1+et(c2(cost+isint)+c3(costisint))x_1(t) = c_1 + e^t (c_2 (\cos t + i\sin t) + c_3 (\cos t - i\sin t))
x2(t)=2c1+et(c2(icostsint)+c3(icostsint))x_2(t) = 2c_1 + e^t (c_2 (i\cos t - \sin t) + c_3 (-i\cos t - \sin t))
x3(t)=c1+et(c2(2icost+2sint)+c3(2icost+2sint))x_3(t) = -c_1 + e^t (c_2 (-2i\cos t + 2\sin t) + c_3 (2i\cos t + 2\sin t))
別の定数を用いて表現すると
x1(t)=c1+et(C2cost+C3sint)x_1(t) = c_1 + e^t (C_2 \cos t + C_3 \sin t)
x2(t)=2c1+et(C3cost+C2sint)x_2(t) = 2c_1 + e^t (-C_3 \cos t + C_2 \sin t)
x3(t)=c1+et(2C3cost2C2sint)x_3(t) = -c_1 + e^t (2C_3 \cos t - 2C_2 \sin t)

3. 最終的な答え

(1) dxdt=(111111210)x\frac{d\mathbf{x}}{dt} = \begin{pmatrix} 1 & -1 & -1 \\ 1 & 1 & -1 \\ 2 & -1 & 0 \end{pmatrix} \mathbf{x}
(2) 固有値: λ1=0,λ2=1+i,λ3=1i\lambda_1 = 0, \lambda_2 = 1+i, \lambda_3 = 1-i
固有ベクトル: v1=(121),v2=(1i2i),v3=(1i2i)\mathbf{v}_1 = \begin{pmatrix} 1 \\ 2 \\ -1 \end{pmatrix}, \mathbf{v}_2 = \begin{pmatrix} 1 \\ i \\ -2i \end{pmatrix}, \mathbf{v}_3 = \begin{pmatrix} 1 \\ -i \\ 2i \end{pmatrix}
(3) 一般解:
x1(t)=c1+et(C2cost+C3sint)x_1(t) = c_1 + e^t (C_2 \cos t + C_3 \sin t)
x2(t)=2c1+et(C3cost+C2sint)x_2(t) = 2c_1 + e^t (-C_3 \cos t + C_2 \sin t)
x3(t)=c1+et(2C3cost2C2sint)x_3(t) = -c_1 + e^t (2C_3 \cos t - 2C_2 \sin t)
ここで、c1,C2,C3c_1, C_2, C_3 は任意の定数です。

「応用数学」の関連問題

直径 $d = 20 \text{ mm}$、長さ $l = 400 \text{ mm}$ の円形断面軸の一端が壁に固定されている。軸端に $T = 300 \text{ Nm}$ のトルクを作用さ...

力学材料力学ねじりトルク極断面二次モーメント横弾性係数単位変換
2025/6/6

ある船が川を $60 km$ 上るのに $5$ 時間、下るのに $3$ 時間かかった。このとき、以下の2つの問いに答える。 (1) この船の静水時の速さを求めなさい。 (2) この川の流れの速さを求め...

速度距離連立方程式文章問題
2025/6/6

2種類の財 $x$ と $y$ があり、効用関数が $u(x, y) = x^{\frac{1}{7}}y^{\frac{6}{7}}$ で与えられています。財 $x$ の価格を $p_x > 0$、...

経済学効用関数最適化ラグランジュ乗数法ミクロ経済学
2025/6/6

2つの財 $x$ と $y$ があり、効用関数が $u(x, y) = x^{\frac{1}{7}}y^{\frac{6}{7}}$ で与えられています。各財の価格は $p_x > 0$、$p_y ...

最適化効用関数ラグランジュ乗数法経済学
2025/6/6

この問題は、効用最大化問題を解くものです。所得$m$、x財の価格$p_x$、y財の価格$p_y$が与えられたとき、それぞれの効用関数$u(x,y)$のもとで、最適な消費計画$(x, y)$を求める問題...

効用最大化ラグランジュ乗数法経済学偏微分
2025/6/6

効用関数 $u(x, y) = xy$ のもとで、x財の価格が $p_x > 0$、y財の価格が $p_y > 0$、所得が $m > 0$ であるときの最適消費プラン (x, y) を求める問題です...

最適化効用関数ラグランジュ乗数法経済学
2025/6/6

$L(x, y, \lambda) = x^\alpha y^{1-\alpha} + \lambda(M - p_x x - p_y y)$ ここで $\lambda$ はラグランジュ乗数で...

経済学ミクロ経済学効用関数需要関数ラグランジュ乗数
2025/6/6

与えられた制約条件の下で、関数を最大化する最適化問題を解きます。 (1) $\max_{x,y} xy$ subject to $x+y-2=0$ (2) $\max_{x,y} x^3y^2$ su...

最適化制約付き最適化ラグランジュの未定乗数法微分最大値
2025/6/6

(1) 制約条件 $x + y - 2 = 0$ の下で、$xy$ を最大化する問題。 (2) 制約条件 $x + 2y - 10 = 0$ の下で、$x^3y^2$ を最大化する問題。

最大化制約条件微分最適化ラグランジュの未定乗数法
2025/6/6

地球と太陽の距離が $1.5 \times 10^{11} \mathrm{m}$ であり、光の速度が $3.0 \times 10^{8} \mathrm{m/s}$ であるとき、光が太陽から地球ま...

物理距離速度時間指数
2025/6/6