学步园

矩阵变换解一元N次方程的源代码，这是这次我们做仿真电学实验软件Team里一个人写的不过他不太喜欢注释，可能看得要辛苦点
using System;<br /> using System.Collections.Generic;<br /> using System.Text;</p> <p>namespace Object<br /> {<br /> public class Matrix<br /> {<br /> public Matrix()<br /> { }</p> <p> public Matrix(double[,] matrix)<br /> {<br /> this.matrix = matrix;<br /> }<br /> double[,] matrix = new double[1, 1];//定义数组<br /> public double[,] GetMatrix//获取设置数组<br /> {<br /> get { return matrix; }<br /> set { matrix = value; }<br /> }<br /> public void Account(double[,] matrix, double[] mb, int number)//消除前排主项<br /> {<br /> for (int i = 0; i < number; i++)<br /> {<br /> for (int j = i + 1; j < number; j++)<br /> {<br /> double temp;<br /> if (matrix[j, i] != 0)<br /> {<br /> temp = matrix[i, i] / matrix[j, i];<br /> for (int k = 0; k < number; k++)<br /> {<br /> matrix[j, k] = matrix[i, k] - matrix[j, k] * temp;//消元<br /> }<br /> mb[j] = mb[i] - mb[j] * temp;<br /> }</p> <p> }<br /> } Console.WriteLine(mb[2].ToString());<br /> for (int i = 0; i < number; i++)<br /> {<br /> for (int j = 0; j < number; j++)<br /> {<br /> Console.Write(matrix[i, j].ToString() + " ");<br /> }<br /> Console.Write(mb[i].ToString());<br /> Console.WriteLine();<br /> }<br /> Console.ReadLine();<br /> }<br /> public void ReAccont(double[,] matrix, double[] mb, int number)////消除后排主项<br /> {<br /> for (int i = number - 1; i >= 0; i--)<br /> {</p> <p> for (int j = i - 1; j >= 0; j--)<br /> {<br /> double temp;<br /> if (matrix[j, i] != 0)<br /> {<br /> temp = matrix[i, i] / matrix[j, i];<br /> for (int k = 0; k < number; k++)<br /> {<br /> matrix[j, k] = matrix[i, k] - matrix[j, k] * temp;//消元 </p> <p> }<br /> mb[j] = mb[i] - mb[j] * temp;<br /> }<br /> }<br /> }<br /> for (int i = 0; i < number; i++) //输出消元后的矩阵<br /> {<br /> for (int j = 0; j < number; j++)<br /> {<br /> Console.Write(matrix[i, j].ToString() + " ");<br /> }<br /> Console.Write(mb[i].ToString());<br /> Console.WriteLine();<br /> }<br /> Console.ReadLine();<br /> for (int i = 0; i < number; i++)<br /> {<br /> if (matrix[i, i] != 0)<br /> {<br /> mb[i] = mb[i] / matrix[i, i];<br /> matrix[i, i] = 1;<br /> }<br /> else<br /> {<br /> mb[i] = 0;<br /> }<br /> }<br /> for (int i = 0; i < number; i++) //高斯-约当消去法<br /> {<br /> for (int j = 0; j < number; j++)<br /> {<br /> Console.Write(matrix[i, j].ToString() + " /t ");<br /> }<br /> Console.Write(mb[i].ToString());<br /> Console.WriteLine();<br /> }<br /> Console.ReadLine();<br /> for (int i = 0; i < number; i++) //输出最后结果<br /> {<br /> Console.WriteLine("x" + i.ToString() + "=" + mb[i].ToString());<br /> }<br /> Console.ReadLine();<br /> }<br /> }<br /> class Program<br /> {<br /> static void Main(string[] args)<br /> {<br /> double[,] mat = new double[,] {<br /> { 2, 2, 3, 2, 3, 8, 4, 6, 9, 10 },<br /> { 4, 7, 7, 12, 55, 6, 9, 8, 4, 65 },<br /> { -2, 4, 5, 5, 4, 6, 8, 4, 6, 9 },<br /> { 0,0,5,3,4,-4,-8,9,6,7},<br /> {1,5,6,8,9,7,4,6,5,2},<br /> {5,8,6,7,11,45,6,25,5,4},<br /> {5,9,8,9,77,8,16,19,20,16},<br /> {8,9.2,65,6.3,5,6,8,4,6,7},<br /> {4,5,8,12,20,65,46,35,10,1},<br /> {0,0,0,0,0,4,5,6,8,9}<br /> };<br /> double[] ma = new double[] { 3, 1, -7, 4, 6, 8, 9, 5, 12, 4 };<br /> Matrix m = new Matrix(mat);</p> <p> m.Account(mat, ma, 10);<br /> m.ReAccont(mat, ma, 10);</p> <p> }<br /> }<br /> }