【正文】 ) { double de, ds。 de = 660 * (dzl / 50)*1000。 ds = * dbj。 return ((*((de/ds),10))*)。 } /******* 互阻抗實(shí)部計(jì)算函數(shù) ******/ public static double huzhukangjisuanshibu() { return (*)。 } /******* 互阻抗虛部計(jì)算函數(shù) ******/ public static double huzhukangjisuanxubu(double ddl, double dzl, double d) { double de。 de = 660 * (dzl / 50)*1000。 return (( * ((de / d), 10))*)。 } }2. 復(fù)數(shù)類(lèi) class Complex { public double real。 public double imaginary。 public Complex(double real, double imaginary) { = real。 = imaginary。 } public static Complex operator +(Complex c1, Complex c2) { return new Complex( + , + )。 } public static Complex operator (Complex c1, Complex c2) { return new Complex( , )。 } public static Complex operator *(Complex c1, Complex c2) { return new Complex( * * , * + * )。 } public static Complex operator /(Complex c1, Complex c2) { double c。 c = * + * 。 return new Complex(( * + * )/c, ( * + * )/c)。 } }3. 矩陣計(jì)算類(lèi)部分函數(shù) class Matrix_operation { /// 返回原始矩陣MatrixA的逆矩陣 public static double[,] MInve(double[,] MatrixA) { //獲取矩陣的行列數(shù) int l = (0)。 int k = (1)。 //定義結(jié)果矩陣 double[,] MatrixC = new double[l, l]。 if (l == k) { //定義分塊擴(kuò)展矩陣 double[,] MatrixB = new double[l, 2 * l]。 //初始化分塊擴(kuò)展矩陣l(2*l)為(A|I) for (int i = 0。 i l。 i++) { for (int j = 0。 j 2 * l。 j++) { if (j l) { MatrixB[i, j] = MatrixA[i, j]。 } else { if (j l == i) { MatrixB[i, j] = 1。 } else { MatrixB[i, j] = 0。 } } } } for (int i = 0, j = 0。 i l amp。amp。 j l。 i++, j++) { if (MatrixB[i, j] == 0) { if (i == l 1) { return null。 } int m = i + 1。 for (。 MatrixA[m, j] == 0。 m++) { if (m == l 1) { return null。 } } if (m == l) { return null。 } else { //把i行和m行相加 for (int n = j。 n 2 * l。 n++) { MatrixB[i, n] += MatrixB[m, n]。 } } } double temp = MatrixB[i, j]。 if (temp != 1) { //把i行數(shù)據(jù)，變成以1開(kāi)始的一行數(shù)據(jù) for (int n = j。 n 2 * l。 n++) { if (MatrixB[i, n] != 0) { MatrixB[i, n] /= temp。 } } } //把i行后的所有行的j列變成0 for (int r = l 1。 r i。 r) { temp = MatrixB[r, j]。 for (int s = j。 s 2 * l。 s++) { MatrixB[r, s] = (MatrixB[i, s] * temp)。 } } } //把分塊擴(kuò)展矩陣l(2*l)前l(fā)l轉(zhuǎn)變?yōu)閱挝痪仃? for (int i = l 2。 i = 0。 i) { for (int j = i + 1。 j l。 j++) { if (MatrixB[i, j] != 0) { double tmp = MatrixB[i, j]。 for (int n = j。 n 2 * l。 n++) { MatrixB[i, n] = (tmp * MatrixB[j, n])。 } } } } //返回逆矩陣 for (int i = 0。 i l。 i++) { for (int j = 0。 j l。 j++) { MatrixC[i, j] = MatrixB[i, j + l]。 } } } else { (所計(jì)算的矩陣行列數(shù)不等！)。 } if (MatrixC == null) { (所計(jì)算的矩陣行列式為零！)。 } return MatrixC。 } /// ，其每一元素是參數(shù)矩陣MatrixA和MatrixB對(duì)應(yīng)位置的元素之和 public static double[,] MsAdd(double[,] MatrixA, double[,] MatrixB) { int row1 = (0)。 int col1 = (1)。 int row2 = (0)。 int col2 = (1)。 double[,] MatrixC = new double[row1, col1]。 if (row1 != row2) { (行數(shù)不相等)。 } else if (col1 != col2) {