怎样用c语言编程解多元方程组,求用C#编程实现：求一个多元一次方程组解

满意答案

yujiabian

2013.11.14

采纳率：54%    等级：12

已帮助：17454人

高斯消元 再代回原方程求解

算法如下

_rows是指方程组行数;_cols是列数

void gaoss(double[,]a)//高斯消元求未知数X，

{

string print = "";

int L = _rows - 1;

int i, j, l, n, m, k = 0;

double[] temp1 = new double[_rows];

/*第一个do-while是将增广矩阵消成上三角形式*/

do

{

n = 0;

for (l = k; l < L; l++)

temp1[n++] = a[l + 1, k] / a[k, k];

for (m = 0, i = k + 1; i < _rows; i++, m++)

{

for (j = k; j < _cols; j++)

a[i, j] -= temp1[m] * a[k, j];

}

k++;

} while (k < _rows);

///*第二个do-while是将矩阵消成对角形式，并且重新给k赋值,最后只剩下对角线和最后一列的数，其它都为0*/

k = L - 1;

do

{

n = 0;

for (l = k; l >= 0; l--)

temp1[n++] = a[k - l, k + 1] / a[k + 1, k + 1];

for (m = 0, i = k; i >= 0; i--, m++)

{

for (j = k; j < _cols; j++)

a[k - i, j] -= temp1[m] * a[k + 1, j];

}

k--;

} while (k >= 0);

/*下一个for是解方程组*/

for (i = 0; i < _rows; i++)

{

double value = a[i, _rows] / a[i, i];

print += "X"+(i+1)+"="+value+" ";

}

MessageBox.Show(print, "方程的解为：");

}

10分享举报