締切済み

JavaScriptについての質問です。

2016/06/13 23:35

ガウスの消去法による解の求め方をプログラムしたいのですが、上手くいかなかったので質問しました。以下のプログラムから実行はできるのですが、解が求められていませんでした。どこを変えればきちんとできるのでしょうか、回答お願いします。 var A=[[1,-1,1],[1,2,0],[2,0,3]],b=[5,1,9]; WScript.Echo("ガウスの消去法の解:\n"+Gauss(A,b)); function Gauss(A,b) { var s,x=new Array(A.length),M=new Array(A.length); for(var i=k+1;i<A.length-1;k++) { for(var i=k+1;i<A.length;i++) { M[i]=new Array(A.length); M[i][k]=(A[i][k])/(A[k][k]); for(var j=k+1;j<A.length ;j++) A[i][j]=A[i][j]-(M[i][k]*A[k][j]); b[i]=b[i]-(M[i][k]*b[k]); } } for(var i=A.lengh-1;i>=0;i--) { s=0; for(var k=i+1;k<A.length;k++) s+=A[i][k]*x[k]; x[i]=(1/A[i][i])*(b[i]-s); } return x; }

tsukibashi
お礼率0% (0/19)

JavaScript
回答数7
ありがとう数0

みんなの回答 （7）
専門家の回答

みんなの回答

amanojaku1
ベストアンサー率54% (265/488)

2016/06/19 11:58 回答No.7

> 「htm」ファイルにコピペして下さい。拡張子が「htm」、「html」のファイルの事です。例えば「Gauss.htm」など。

ログインすると、全ての回答が全文表示されます。

amanojaku1
ベストアンサー率54% (265/488)

2016/06/16 23:14 回答No.6

> for( t = MA.length-1; 0<=t; t--){ 無駄がありました、下記のように変更して下さい。 for( t = MA.length-1; 1<=t; t--){

ログインすると、全ての回答が全文表示されます。

amanojaku1
ベストアンサー率54% (265/488)

2016/06/16 22:48 回答No.5

> (ガウスの消去法にしろ、ガウス・ジョルダン法にしろ)向かって左から袈裟懸けになる行列の対角成分：1にしたい対角要素がゼロの場合は計算不能向かって左から袈裟懸けになる行列の対角成分：1にしたい対角要素がゼロにならないように行を交換する手法があるようです。

ログインすると、全ての回答が全文表示されます。

amanojaku1
ベストアンサー率54% (265/488)

2016/06/16 22:28 回答No.4

> GJ = GaussJordan(MA,MB); > document.write('GaussJordan(MA,MB); '); > function GaussJordan(MA,MB){ 「GaussJordan」を「Gauss」に変更して下さい。

ログインすると、全ての回答が全文表示されます。

amanojaku1
ベストアンサー率54% (265/488)

2016/06/16 22:12 回答No.3

これはガウス(Gauss)の消去法です(^_^;こちらの方が簡単でしたね) (ガウスの消去法にしろ、ガウス・ジョルダン法にしろ)向かって左から袈裟懸けになる行列の対角成分：1にしたい対角要素がゼロの場合は計算不能となり、ゼロでないにしても非常に小さい値の場合、誤差が大きくなってしまうと言う問題があるようです。「htm」ファイルにコピペして下さい。 <html> <head> <meta http-equiv="Content-Type" content="Text/HTML; charset=Shift_JIS"> <title>test</title> </head> <body bgcolor="#8f9fff" text="#000000" link="#2f2fff" vlink="#d700d7" alink="#00ffff" > <script type="text/javascript">  </script> </body> </html>

ログインすると、全ての回答が全文表示されます。

amanojaku1
ベストアンサー率54% (265/488)

2016/06/16 20:53 回答No.2

> function GaussJordan(MA,MB){ の次に下記を追加して下さい。 var i,j,t,tg;

ログインすると、全ての回答が全文表示されます。

amanojaku1
ベストアンサー率54% (265/488)

2016/06/16 20:46 回答No.1

これはガウス・ジョルダン(Gauss-Jordan)法です。 (ガウスの消去法にしろ、ガウス・ジョルダン法にしろ)向かって左から袈裟懸けになる行列の対角成分：1にしたい対角要素がゼロの場合は計算不能となり、ゼロでないにしても非常に小さい値の場合、誤差が大きくなってしまうと言う問題があるようです。「htm」ファイルにコピペして下さい。 <html> <head> <meta http-equiv="Content-Type" content="Text/HTML; charset=Shift_JIS"> <title>test</title> </head> <body bgcolor="#8f9fff" text="#000000" link="#2f2fff" vlink="#d700d7" alink="#00ffff" > <script type="text/javascript">  </script> </body> </html>

ログインすると、全ての回答が全文表示されます。