No.7ベストアンサー
- 回答日時:
x”(t)=A・x(t)と直接解いたほうが簡単ですが
システマチックでなく発展性がなく運に左右される解きかたなので
コンピュータ計算に向いているy’(t)=C・y(t)によって解きましょう
というより
y(t)=exp(C・t)・y(0)という風に既に解けているのですね
exp(C・t)を簡単に表現するかどうかは好みの問題と言い張ることができる
何しろ2次のベクトル微分方程式を1次ベクトル微分方程式にするやり方は普遍てきなのだから
x”(t)=A・x(t)だとAをいちいち対角かしないと行けないが
y’(t)=C・y(t)はもうすでに解けているのだからほうっておいても言いのです
子の回答を持っていって文句を言う先生を論破してください
nubouさんの回答を参考に友達と考えた結果、
ようやく理解することができました。
多くの御回答ほんとうにありがとうございました。
No.6
- 回答日時:
正方行列が湯に足り行列によって対角かされるための必要十分条件はその正方行列が正規行列であることである
という低利があります
C=
[0 E]
[A 0]
でありAは対象行列だから
C・C^T=C^T・C=
[E 0]
[0 A^2]
でありAは正規行列になります
従ってCは対角行列に相似なのです
これはさすがに猿にも人間にも難しいジョルダンの標準形を持ち出さなくても子の問題を解けるということを意味します
No.5
- 回答日時:
私の質問の回答は
A=
[-2 0 0]
[1 -2 1]
[0 1 -2]
でなく
A=
[-2 1 0]
[1 -2 1]
[0 1 -2]
です
書いた後間違いに気がついたので理解を確認するために宿題にさせてもらいました
さてジョルダンの標準形ですがCの子誘致がすべて異なればCは対角行列に相似です
またたとえCの固有多項式が多重根を持っていてもCが独立な固有ベクトルを6つ持てばCは対角行列に相似です
任意の行列は対角行列に相似だとは限らないがジョルダンの標準形には相似です
No.4
- 回答日時:
忙しかったプログラミングも一服したので暇つぶしを兼ねて出血サービスしちゃいましょう
x1(t)を左の失点の右への変異とし
x2(t)を中の失点の右への変異とし
x3(t)を右の失点の右への変異とし
x(t)=[x1(t),x2(t),x3(t)]^Tとし
m1/k=αとし
m2/k=βとし
m3/k=γとし
P=
[α 0 0]
[0 β 0]
[0 0 γ]
A=
[-2 0 0]
[1 -2 1]
[0 1 -2]
とし
B= P^(-1)・Aとすれば
x(t)''=B・x(t)
さらに
y(t)=[x(t)^T,x1(t)’,x2(t)’,x3(t)’]^T
とし
Eを3次正方行列とし
C=
[0 E]
[B 0]
とすれば
y(t)’=C・y(t)
解いて
y(t)=exp(C・t)・y(0)
Jを適当な6次ジョルダンの標準形行列とし
Qを適当な6次正則行列として
J=Q^(-1)・C・Qであるならば
y(t)=Q・exp(J・t)・Q^(-1)・y(0)
となる
exp(J・t)はきわめて簡単になるので後は猿でもできます
さて記述に間違いがあります
補足に書いてください
この回答への補足
実は解答の指針というのが決められていて、
1.3つの質点の運動方程式を立てる
2.運動方程式に現れる係数の作る対称行列を対角化する直交行列を求める。そのために行列の固有値を求め、対応する規格化した固有ベクトルから直交行列を作る。
3.対角化されて求まった基準座標の従う運動方程式を書き下し、各振動の方程式から一般解を求める。
というものです。
2.3.がこの問題の難しい部分なのでしょうが、肝心の運動方程式を立てるところが・・・さっぱりなんです。
せっかく回答いただいたのに、『ジョルダンの標準形』というものを習っていないものでよく分からないのです。本当にスイマセン。
No.3
- 回答日時:
失礼しました
ちょん簿です
x1(t)を左の失点の右への変異とし
x2(t)を中の失点の右への変異とし
x3(t)を右の失点の右への変異とし
x(t)=[x1(t),x2(t),x3(t)]^Tとし
m1/k=αとし
m2/k=βとし
m3/k=γとし
P=
[α 0 0]
[0 β 0]
[0 0 γ]
A=
[-2 0 0]
[1 -2 1]
[0 1 -2]
とすれば
x(t)''=P^(-1)・A・x(t)
なので
P^(-1)・Aを対角化またはジョルダンの標準形に変換すればすぐにもとまりますね
No.2
- 回答日時:
x1(t)を左の失点の右への変異とし
x2(t)を中の失点の右への変異とし
x3(t)を右の失点の右への変異とし
x(t)=[x1(t),x2(t),x3(t)]^Tとし
m1/k=αとし
m2/k=βとし
m3/k=γとし
P=
[α 0 0]
[0 β 0]
[0 0 γ]
A=
[-2 0 0]
[1 -2 1]
[0 1 -2]
とすれば
x(t)’=P^(-1)・A・x(t)
なので
P^(-1)・Aを対角化またはジョルダンの標準形に変換すればすぐにもとまりますね
No.1
- 回答日時:
解析力学を知っているなら3つの質点の平衡状態からのずれをu_1、u_2、u_3とでもしてラグランジアンLを書いてオイラー・ラグランジュ方程式から運動方程式を求めるのが間違いないかと。
(ポテンシャルを Uとすれば
U=1/2 k u_1^2+1/2 k (u_2-u_1)^2+1/2 k (u^3-u_2)^2+1/2 k u_3^2)
運動方程式を書いたらそれを行列の形で書いて線形代数の処方箋に従って行列を
対角化すれば解けます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 物理の単振動の問題で分からない所を教えてください 1 2023/05/10 20:59
- 物理学 高1力学の運動量の問題です。問題を一通り解いたのですが、行き詰まってしまったのでご回答頂ければ嬉しい 3 2022/06/29 11:20
- 物理学 写真の図は単振動の動きを段階的に表したものです。 (加速度=a、力=F、ばね定数=k、物体の質量=m 11 2022/08/24 21:57
- 物理学 写真の図は単振動の動きを段階的に表したものです。 (加速度=a、力=kx、ばね定数=k、物体の質量= 8 2022/08/24 23:39
- 物理学 半径r、質量Mの半円板を円板平面平行に円の中心を回転軸として微小振幅で振り子運動をさせる。このときに 4 2023/08/10 14:08
- 数学 数学微分方程式の問題です。次に書く問題を教えて欲しいです。上端を固定された長さlの棒の先に質量mの質 2 2022/04/29 21:27
- 物理学 物理力学の問題を教えてください 2 2022/07/21 15:18
- 物理学 高1物理です。一通り解いてみたのですが、力学的エネルギーの計算で行き詰まってしまいました。解法を示し 2 2022/07/11 00:47
- 物理学 量子力学 生成消滅演算子 2 2022/08/04 23:17
- 物理学 振り子運動について 2 2023/01/08 16:34
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
2物体の運動を重心系で考えると...
-
高校物理基礎で、変位と位置の...
-
高1力学の運動量の問題です。問...
-
物理なんですけど、変位=x座標...
-
スパン方向とはどの方向ですか?
-
運動方程式の時間反転
-
ラディアル方向・タンジェンシ...
-
座標変換に関して質問です。参...
-
「陽に含まない」について
-
変位座標から角速度・角度の求め方
-
質量m 半径aの一様な円環の慣性...
-
LEDの色温度と色座標について
-
遠心力とコリオリ力について
-
平行軸の定理について
-
力学の問題です。(空気抵抗あり...
-
新幹線での座席回転について
-
波数ベクトル
-
PDF-XChange Viewerで、回転し...
-
角パイ・単管パイプの耐荷重を...
-
図・図形の回転ハンドルが出ない
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
2物体の運動を重心系で考えると...
-
質量m 半径aの一様な円環の慣性...
-
電磁気学でわからないところが...
-
スパン方向とはどの方向ですか?
-
径方向?放射方向?
-
高校物理基礎で、変位と位置の...
-
ばね定数の異なる3つのバネと...
-
段差を乗り越えるのに必要なト...
-
「陽に含まない」について
-
2つのバネに挟まれた物体の振動...
-
2物体の慣性モーメント
-
ラディアル方向・タンジェンシ...
-
表式ってなんですか?数学用語?
-
力学の問題です。(1)しかわかり...
-
すべり面とすべり方向
-
SPEC 時間を止めて、銃で放たれ...
-
特殊相対性理論におけるトンネ...
-
変位座標から角速度・角度の求め方
-
エネルギー積分とは? 運動方程...
-
循環座標(解析力学)
おすすめ情報