行列の連立方程式 (mathematica)

締切済

質問者：krclight417
質問日時：2012/11/19 21:08
回答数：5件

mathematica4.1を使用して，行列で表した連立方程式を計算したいと思っています．例えば

A= a11 a12 a13
a21 a22 a23
a31 a32 a33
B=b11 b12 b13
b21 b22 b23
b31 b32 b33

C,Dも同様に定義して
A .x+B .y=1
C .x+D .y=0
x=x1 y=y1
x2 y2
x3 y3
のような連立方程式のx,yについて解きたいのですが，どのようにmathematicaで表現すればよいかが分かりません．
どなたか教えてください．お願いします．
当然ではありますが，A,B,C,Dはすべて既知としています．

通報する

この質問への回答は締め切られました。

質問の本文を隠す

回答 (5件)

最新から表示
回答順に表示

No.5

回答者： kazuma1956
回答日時：2013/10/12 01:24

Mathematicaでの回答例

(* 以下の様に定義します *)
(* 例えば、Aの(i,j)成分は、a[i,j]となります *)
A = Array[a, {3, 3}];
B = Array[b, {3, 3}];
CC = Array[c, {3, 3}];
DD = Array[d, {3, 3}];
X = Array[x, {3}];
Y = Array[y, {3}];

(* このまま、
Solve[{A.X+B.Y={1,1,1},C.X+D.Y=={0,0,0}},Join[X,Y]]
を実行したいところですが、結構時間がかかるようです。*)

(* {{A,B},{C,D}}を結合した大きな行列を作ります。*)
M = ArrayFlatten[{{A, B}, {CC, DD}}];

(* M.{x[1],x[2],x[3],y[1],y[2],y[3]}={1,1,1,0,0,0}
を解けばよいのですが、これも時間がかかります *)

(*そこで、LinearSolveを使います。これは短時間で処理が済みますが、出力
が膨大になるようです *)
(* このような計算のとき、Det[M] !=0 という前提で、計算されるので、ご注意ください *)

LinearSolve[M][{1,1,1,0,0,0}]

- 0
- 件

通報する

No.4

回答者： 178-tall
回答日時：2012/11/21 19:26

Mathematica の使用経験はありませんが、たとえば、

(Linear Algebra パッケージ？) MatrixManipulation パッケージに
BlockMatrix 関数があるようです。

　　　

参考URL：http://mathematica.stackexchange.com/questions/7 …

- 0
- 件

通報する

この回答へのお礼

遅くなりまして申し訳ありません．回答有り難うございます．
しかし，mathematicaの経験が浅いため，このBlockMatrix関数をどう適用すればよいのか…
勉強してみます．

通報する

お礼日時：2012/11/26 16:12

No.3

回答者： alice_44
回答日時：2012/11/20 13:54

行列を並べて作った行列を「ブロック行列」とか「区分行列」とか言うが、

ブロック行列の逆行列について、↓のような資料があった。
http://ja.wikipedia.org/wiki/%E5%8C%BA%E5%88%86% …

あるいは、スカラーの 2 元 2 連立方程式と同じように、
C x + D y = 0 を y = … の形に解いて A x + B y = 1 へ代入する
などでもよいのかもしれない。
A, B, C, D の中にひとつでも正則な行列があれば、このやり方がスムースにいく。

たとえ四つとも非正則であっても、連立一次方程式の不能の場合、不定の場合を処理する
方法を知っていれば、det D = 0 の時に C x + D y = 0 を y = … と扱うことが
できると思う。解の一意性を判定して不定解を正しく扱う部分を、プログラムで自動的に
やらせるには、ゴタゴタした場合分けをコーディングしなければならないだろうけど。

# やっぱり、6 連立方程式に展開してしまったほうが
# 楽なように見えるなあ。