数值分析(php兑现)二:线性方程组的两个解法

程序员文章站 2022-04-11 09:17:38

...

数值分析(php实现)二:线性方程组的两个解法
其实这个程序是早就写好的,只不过一直没有时间写上来.
这两个程序都是通过迭代的方法得到线性方程组的解的方法,一个是高斯-塞德尔迭代法,一个是雅可比迭代法

matrix = $array;    }    public function solve(){        $preX = array();        $nowX = array();        $cishu = 17;        $delta = 0.0001;        $matN = count($this->matrix);        for($i = 0;$imatrix[$i][$j] * $nowX[$j]);                    if($j > $i) $sum2 += ($this->matrix[$i][$j] * $preX[$j]);                }                $nowX[$i] = ($this->matrix[$i][$matN] - $sum1 - $sum2)/$this->matrix[$i][$i];                $tempMin = $nowX[$i] > $preX[$i] ? $nowX[$i]-$preX[$i]:$preX[$i]-$nowX[$i];                if($min > $tempMin) $min = $tempMin;            }             $preX = $nowX;            $str = implode(",",$nowX);            echo ($n+1).":($str)"."
";            if($min solve();?>

接下来的是雅可比方法,

matrix = $array;    }    public function solve(){        $preX = array();        $nowX = array();        $cishu = 17;        $matN = count($this->matrix);        for($i = 0;$imatrix[$i][$j] * $preX[$j]);                }                $nowX[$i] = ($this->matrix[$i][$matN] - $sum)/$this->matrix[$i][$i];            }             $preX = $nowX;            $str = implode(",",$nowX);            echo ($n+1).":($str)"."
";        }    }}$a = array(    array(5,2,1,-12),    array(-1,4,2,20),    array(12,-3,10,3));$x = new Yacobi($a);$x->solve();?>

声明：本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系admin@php.cn核实处理。