新しく発言する  EXIT  インデックスへ

ガウスの消去法連立1次方程式の解Ax=b


  ガウスの消去法 連立1次方程式の解 Ax=b 山中和義 2007/11/25 10:21:25 
  !ガウス・ジョルダン(Gauss-Jordan)法N元... 山中和義 2007/11/25 10:22:50 
  │└FullBASICなら逆行列による計算がお手軽です... 山中和義 2007/11/25 10:28:23 
Re: !ガウス・ジョルダン(Gauss-Jordan)法N元...  返事を書く  ノートメニュー
山中和義 <drdlxujciw> 2007/11/25 10:28:23
FullBASICなら逆行列による計算がお手軽です。
行列計算が可能ですから、記述が簡潔になります。


!逆行列による N元連立1次方程式の解 Ax=b

OPTION ARITHMETIC rational !有理数(分数)

! 3*x1 -4*x2 +2*x3 = 0
! 2*x1 +5*x2 +3*x3 = 6
!-2*x1 +3*x2 -2*x3 = 1

!連立方程式を記述する
LET N=3 !変数(元)の数

DATA 3,-4, 2 !左辺
DATA 2, 5, 3
DATA -2, 3,-2

DIM A(N,N) !係数のみ配列で記録する
MAT READ A

DATA 0 !右辺
DATA 6
DATA 1
DIM b(N)
MAT READ b

DIM x(N) !解

!行列Aの逆行列をA'とする。Ax=bに左からかけると、A'Ax=A'b ∴x=A'b

DIM Ai(N,N) !逆行列を求める
MAT Ai=INV(A)
MAT PRINT Ai !確認

MAT x=Ai*b !左からかける

MAT PRINT x !結果を表示する

END
  ヤコビ反復法(Jacobi) 山中和義 2008/09/29 10:28:16 
  │└ガウス・ザイデル反復法(Gauss-Seidel) 山中和義 2008/09/29 10:32:12 
  !共役勾配法(conjugategradientmethod)N元... 山中和義 2008/10/06 20:58:41 

 インデックスへ  EXIT
新規発言を反映させるにはブラウザの更新ボタンを押してください。