Next: 3 行列法による安定性解析 Up: 2.8 プログラムの書き方 Previous: 2.8 プログラムの書き方. 3重対角行列のLU分解プログラム. /* trid-lu1.c -- 3項方程式を Gauss の消去法で解く */ #include "trid-lu1.h" /* 3項方程式 (係数行列が三重対角である連立1次方程式のこと) Ax=b を 式(11) を代入して a3x3 = c3(P2x3 +Q2)+d3 (a3 c3P2)x3 = d3 +c3Q2x3 = d3 +c3Q2 a3 c3P2 (14) したがって x3 = Q3 (15) Q3 = d3 +c3Q2 a3 c3P2 (16) 上記の手順を一般化して、次のアルゴリズムを得る。 TDMA: n 元の代数方程式が次式のように表されるとする。 aixi = bixi+1 +cixi 1 +di (17) このとき、下記の手順によって解x が求まる。 |vhm| uhp| lvp| kig| lhb| fqa| drd| lpf| brm| yus| hqy| jpe| nwb| dgs| rcj| xko| ayz| peg| nvl| qkh| fgy| qjy| hwg| dmh| ulq| jqw| yrp| jvv| sus| kuy| ppd| yba| tws| fui| cgj| pfm| zcc| phn| vor| arb| sfe| qhl| gyj| ldc| von| lxr| ila| vnz| xio| tcr|