- ベストアンサー
連立方程式(行列)
- みんなの回答 （2）
- 専門家の回答
質問者が選んだベストアンサー
＞言語化してもらえると…ってのはさすがに頼りすぎですかね?ｗ ですねぇ・・・。(笑) よっぽど切羽詰った状況なら手を貸してもいいですが、 別にそこまで大したものでもないでしょう？ まずは、自分で考えてみて下さい。どうしてもダメだったり、 機能拡張したいけど方法が分からないとかあれば、また書き込 んでおいて下さい。 ＞逆行列の方はもうちょっと検索知してみますね。 n行1列のベクトルが、n行n列の単位行列になるだけで本質的な 違いは何一つありません。結局の所、単に掃き出し法をしてる だけですよ。 頑張ってみてください。
その他の回答 (1)
- POTATO_XP
- ベストアンサー率10% (24/230)
連立方程式の解を求めるアルゴリズムとしては、直接法、反復法１、反復法２と数種類あります。 ■直接法 ガウスの消去法、LU分解法 ■反復法１ ヤコビ法、ガウス・ザイデル法 ■反復法２ 最急降下法、共役勾配法、ICCG法 etc... 検索して、調べてみてください。参考までに、ガウスの消去法の手順を示しておきます。 n=0から初める、n=0～N-1の範囲で以下を繰り返す 1).n行n列成分の値でn行成分の値を全て割る。 m=0から始める、m=0～N-1の範囲で以下を繰り返す、但しm=nとなる場合は処理をとばす 2).m行n列成分の値で倍したn行成分でm行成分を全て引く。 注意する事は、1).の処理でn行n列成分０の場合の特殊処理を入れる事。あと、n行n列成分の値を割る場合、この値が極力大きな値となる様に工夫した方が精度が良くなります。実用的な所だと100×100以上の規模の計算も必要になる場合もあり、誤差が累積するとの理由で直接法はあまり用いられず、反復法の出番になります。 以上、ご参考まで。
お礼
消去法の手順を言葉で説明するのに困ってたので助かりました。 逆行列の方はもうちょっと検索知してみますね。 言語化してもらえると…ってのはさすがに頼りすぎですかね?ｗ
関連するQ&A
- n≠mの場合の連立一次方程式の求解について
連立1次方程式AX=Bと解についてお尋ねします。 この問題はいろんなケースがあり、Aの行列がnxmでn=m, n>m n<mで分類されます。あと行列式がゼロかどうかとかです。n≠mを考えるのですがtA(=Aの転置）を右からかけるとtAAX=tABとなり、CX=D(C=tAA(←対称な正方行列), D=tAB）となるので、det(C)がゼロでなければ解が求まることになります。CX=Dの解は求解可能なら求め方はいくつもありそうです。しかし、考えてみるとこの解Xはどういう解になるのでしょうか。n>mの場合、式が未知数より多いのでデータ解析での最小二乗法の解になるのでしょうか。一方でn<mだと解は求まらないはずなのでこれはdet(C)はゼロになるのでしょうか。 今一度n>mの場合ですが、式が多いのですが、2つの式を足して1つにするなどすると、式を減らすことができると思います。そのようにしてmxmの正方行列（対象は期待できないかも）として求解することは可能ですね。そのような解はどういう解になるのでしょうか。有限要素法による求解がこれに相当するように見えるのですが。 初等的な線形代数には載っていないと思いますが、現実的な問題に取り組むと初めから正方行列ではない場合が多いと思いますが。n≠mの解の一般的な考え方なのだろうと思いますが。 よろしくお願いします。
- ベストアンサー
- 数学・算数
- 逆行列が存在しない連立一次方程式
Ax = b Aはm×n行列 m≠n xはn×1 bはm×1 という式があるとき、 x = (At・A)^-1・b At： Aの転置 と変形すると一般化逆行列を求めなくてもxが求められるように思うのですが、この式変形は間違ってますか？
- ベストアンサー
- 数学・算数
- 連立方程式についてのご質問
連立方程式についてのご質問です。 X_sは、N*1ベクトル A_ijは、N*N行列 Bは=(1,....,1)'N*1ベクトル |A_11 A_12 ・・・ A_1K | |X_1| |B | |A_12 A_22 ・ |＊| ・ |＝| ・| | ・ ・ ・ | | ・ | | ・| |A_K1 A_k2 ・・・ A_KK | |X_k| |B | (| |は行列を表しています。) 上記式をXについて解きたいのですが、各変数がベクトル・行列表記なのでよく分かりません。 X=～。の形で解を求めることができるのでしょうか？ (やはりN>K,N=K,N<Kの場合によって違う・・？) それと、具体的な数値解はどのように求めるのでしょう？ A_ij,X_s,Bがスカラーである場合は、ガウスの消去法、吐き出し法などで計算することができます。 上記のように各要素がベクトル・行列表記の場合はどのように計算するのでしょうか。 未熟な質問で申し訳ございませんが、よろしくお願い致します。
- 締切済み
- 数学・算数
- 連立線形方程式の数値解
Ａ＝｜－２ ０ －１ ８｜ ｂ＝｜８｜ ｜ ０ １ ０ ０｜ ｜９｜ ｜ ７ －３ －１ ０｜ ｜２｜ ｜ ０ －２ ４ １｜ ｜９｜ 上のＡ、ｂについて連立線形方程式Ａｘ＝ｂの数値解を解く問題で、この行列式のままだとヤコビやガウス＝ザイデル法では解けない（収束しない）ので、これを収束するように修正したいのですが色々試してみてもできず、ここで聞くに至りました。 収束させるためには具体的にどうしたらいいのか教えて下さい。よろしくお願いします。
- ベストアンサー
- 数学・算数
- 行列・非同次連立一次方程式
「非同次連立方程式の解が一意に決まるとき、係数行列と拡大行列の階数の関係と、その理由を説明せよ」という問題について、やや曖昧なところがあるので質問します。 まず階数の関係については、係数行列をm行n列としたとき、（係数行列の階数）＝（拡大行列の階数）＝nとなった場合に、解が一意に定まるという解答だと思うのですが、これで正しいでしょうか？ 続いて、その理由について、係数行列がn列ということは、未知数がn個あり、階数がnということは式がn個並ぶはずなので、n個の未知数にn個の式で一意解になるというのは、なんとなくわかるのですが、やや説明が曖昧な気がして、悩んでいます。 もし、これよりも明解な理由があれば、ご教授願いたいです。よろしくお願いします。
- ベストアンサー
- 数学・算数
- 行列について教えてください。
行列についていくつか質問させて下さい。 １）(A+B)^-1の結果 （^-1は逆行列です。自分ではそのままA^-1+B^-1だ と思うのですが・・・） ２）置換行列について 例えば、線形方程式Ax=b(A:正方行列,x:求める解)を解く問題でA^=PAQ,x^=Q^-1x,b^=Pbと置換行列P,Qを用いてその後の処理（行列分割）を行うと、本で読んだのですがなぜ、置換行列を使っておきかえるのでしょうか？ 置換しなくてもその後の処理（行列分割）はできると思います。 一つでもいいので解かる方よろしくお願いします。
- ベストアンサー
- 数学・算数
- 連立一次方程式を解くプログラムについて
数値計算の本を見たら必ず載っている連立1次方程式の解法ですが、どのようなタイプの行列でも解くことができるものにはどのような解法があるでしょうか。もちろん、解くことができる範囲でということではあります。その意味でガウスの消去法（ピボット付）になるでしょうか。ガウスの消去法は解き方に基本的な制約はないですね。一方、共役勾配法の説明を見ると"対称正定値行列の場合、..."となっており、その範囲でしか考えていないということでしょうか。そうなるとかなり絞られることになってしまいます。任意の行列は変換して対称正定値に変換できる、ということでもないと思いますが。 有限要素法に関連した連立方程式解法についても書籍1冊分の解説とかありそうですが。高速化のために長い解説があったとしても前提によって使える範囲が狭いものが多いように思えるのですが。よろしくお願いします。
- 締切済み
- 数学・算数
- 不等号の制約つきの連立一次方程式の解き方
連立一次方程式に対して、不等号のついた制約条件があった場合、どのように解けばよいのでしょうか。制約条件がないときは、一般的には、ガウスの消去法などで、行列の演算（Ａｘ＝Ｂの形）をすればよいと思うのですが、それに対して、その中の変数に条件、例えば、x1>0,x2>0,,,,などがついたときの解き方がわかりません。週末の間に解決したいので、すぐに回答が欲しいです。よろしくお願い致します。
- ベストアンサー
- 数学・算数
お礼
わざわざ私の戯言(ｗ)にまでつきあっていただいてありがとうございます。 あとは自力でがんばります! もし多少違ってもそれはそれで勉強になりますしねｗ