sample4_4.java

《Java数值计算算法编程》随书代码 已生成ＤＯＣＳ说明文件

/*
 * 示例程序Sample4_4: LEquations类的复系数方程组的全选主元高斯消去法
 */
package javaalgorithm.sample;

import javaalgorithm.algorithm.Complex;
import javaalgorithm.algorithm.Matrix;
import javaalgorithm.algorithm.LEquations;

public class Sample4_4
{
	public static void main(String[] args)
	{
		// 系数矩阵数据
		// 实部
		double[] mtxDataCoef3Real = {
				1.0,3.0,2.0,13.0,
				7.0,2.0,1.0,-2.0,
				9.0,15.0,3.0,-2.0,
				-2.0,-2.0,11.0,5.0};
		// 虚部
		double[] mtxDataCoef3Imag = {
				3.0,-2.0,1.0,6.0,
				-2.0,7.0,5.0,8.0,
				9.0,-3.0,15.0,1.0,
				-2.0,-2.0,7.0,6.0};
		// 常数矩阵数据
		// 实部
		double[] mtxDataConst3Real = {
				2.0,
				7.0,
				3.0,
				9.0};
		// 虚部
		double[] mtxDataConst3Imag = {
				1.0,
				2.0,
				-2.0,
				3.0};
		
		// 构造系数矩阵
		Matrix mtxCoef3Real = new Matrix(4, mtxDataCoef3Real);
		Matrix mtxCoef3Imag = new Matrix(4, mtxDataCoef3Imag);
		// 构造常数矩阵
		Matrix mtxConst3Real = new Matrix(4, 1, mtxDataConst3Real);
		Matrix mtxConst3Imag = new Matrix(4, 1, mtxDataConst3Imag);

		// 构造线性方程组
		LEquations leqs3 = new LEquations(mtxCoef3Real, mtxConst3Real);
		
		// 复系数方程组的全选主元高斯消去法
		Matrix mtxResult3Real = new Matrix();
		Matrix mtxResult3Imag = new Matrix();
		if (leqs3.GetRootsetGauss(mtxCoef3Imag, mtxConst3Imag, mtxResult3Real, mtxResult3Imag))
		{
			for (int i=0; i<mtxConst3Real.getNumRows(); ++i)
			{
				Complex cp = new Complex(mtxResult3Real.getElement(i,0), mtxResult3Imag.getElement(i,0));
				System.out.println(cp);
			}
		}
		else
		{
			System.out.println("失败");
		}
	}
}