importClass(java.text.DecimalFormat);
let df = new DecimalFormat("0.##");
/***
@return 计算的结果
*/
function mathDeterminantCalculation(value) {
// 当矩阵的行数大于2时
for (let i = 0; i < value.length; i++) {
// 检查数组对角线位置的数值是否是0,如果是零则对该数组进行调换,查找到一行不为0的进行调换
if (value[i][i] == 0) {
value = changeDeterminantNoZero(value, i, i);
}
for (let j = 0; j < i; j++) {
// 让开始处理的行的首位为0处理为三角形式
// 如果要处理的列为0则和自己调换一下位置,这样就省去了计算
if (value[i][j] == 0) {
continue;
}
// 如果要是要处理的行是0则和上面的一行进行调换
if (value[j][j] == 0) {
let temp = value[i];
value[i] = value[i - 1];
value[i - 1] = temp;
continue;
}
let ratio = -(value[i][j] / value[j][j]);
value[i] = addValue(value[i], value[j], ratio);
}
}
return value;
}
/***
@return 将i行i列之前数字置换为0后的新的行
/
function addValue(currentRow,frontRow,
ratio){
for (let i = 0; i < currentRow.length; i++) {
currentRow[i] += frontRow[i] ratio;
currentRow[i] =parseFloat(currentRow[i]);
//( df.format(currentRow[i]));
//Double.parseDouble(df.format(currentRow[i]));
}
return currentRow;
}
/**
要判断的列
*/
function changeDeterminantNoZero(determinant,
line, row){
for (let j = line; j < determinant.length; j++) {
// 进行行调换
if (determinant[j][row] != 0) {
let temp = determinant[line];
determinant[line] = determinant[j];
determinant[j] = temp;
return determinant;
}
}
return determinant;
}
/**
@return 增广矩阵
*/
function transferMatrix(coefficient,
value) {
let temp = new Array(coefficient.length);
if (coefficient.length != value.length) {
return temp;
}
// 将方程值添加到系数矩阵中
for (let i = 0; i < coefficient.length; i++) {
temp[i] =new Array(coefficient[0].length + 1);
for (let j = 0; j < coefficient[0].length; j++) {
temp[i][j] = coefficient[i][j];
}
}
for (let i = 0; i < value.length; i++) {
temp[i][temp[i].length - 1] = value[i];
}
return temp;
}
/**
@return 非零行的个数
*/
function effectiveMatrix(value) {
for (let i = value.length - 1; i > -1; i--) {
for (let j = 0; j < value[i].length; j++) {
if (value[i][j] != 0) {
return i + 1;
}
}
}
return 0;
}
/**
@return 方程组的解
/
function calculationResult(mathMatrix) {
// 有解时方程组的个数等于方程组的未知数
let result = new Array(mathMatrix.length);
log(mathMatrix);
for (let i = mathMatrix.length - 1; i > -1; i--) {
let temp = 0;
for (let j = mathMatrix[i].length; j > 0; j--) {
// 第一个为方程的解,需要将解赋值给临时变量
if (mathMatrix[i][j - 1] != 0) {
if (j == mathMatrix[i].length) {
temp = mathMatrix[i][j - 1];
} else if ((j - 1 > -1 )&& ((result[j - 1])!=undefined)) {
temp -= mathMatrix[i][j - 1] result[j - 1];
//log((j - 1 > -1 )&&(typeof(result[j - 1])!=undefined)+";j="+j+";undefined="+(typeof(result[j - 1])!=undefined));
} else {
result[i] = temp / mathMatrix[i][j - 1];
continue;
}
}
}
}
return result;
}
function main() {
/
// 方程的未知数的个数
int n = 3;
// 系数矩阵
double[][] test = { { 2, 3, 1 }, { 1,1, 1 }, { 1, 2, -1 } };
// 方程的解
double[] value = { 11, 6, 2 };
方程组的解为x1=1.0 方程组的解为x2=2.0 方程组的解为x3=3.0
*/
/*
// 方程的未知数的个数
int n = 4;
// 系数矩阵
double[][] test = { { 3, 1, -1, 1 },{ 1, -1, 1,2 }, {2,1,2,-1},{ 1,0, 2, 1 } };
// 方程的解
double[] value ={ -3, 4, 7,6 };
方程组的解为x1=1.0 方程组的解为x2=-2.0 方程组的解为x3=3.0 方程组的解为x4=-1.0
/
/
// 方程的未知数的个数
int n = 4;
// 系数矩阵
double[][] test = { { 1, -3, 4, -5 },{ 1, -1, -1,2 }, {1,2,0,5},{ 2,-1, 3, -2 } };
// 方程的解
double[] value ={ 0, 0, 0,0 };
方程组有零解
*/
// 方程的未知数的个数
let n = 5;
// 系数矩阵
let test = [[2,1, 1,1,1 ],[ 1, 2, 1,1,1 ], [1,1,3,1,1],[ 1,1,1,4,1 ],[1,1,1,1,5]];
// 方程的解
let value =[ 4,5,9,0,-5 ];
n=3;
value=[11,0,-2];
test=[[3,1,-2],[1,1,-1],[-1,-1,-3]]
try {
// 转换成增广矩阵并进行初等行变化
let mathMatrix = mathDeterminantCalculation(transferMatrix(
test, value));
// 找出非零行的个数
let checkMatrixRow = effectiveMatrix(mathMatrix);
// 根据未知数的个数和方程组非零行的个数来判断当前方程组的解的情况
if (n > checkMatrixRow) {
toastLog("未知数有" + n + "个,消元法后获取的阶梯方程组有"
checkMatrixRow + "个方程,等于未知数个数,所以该方程有解");
let result = calculationResult(mathMatrix);
for (let i = 0; i < result.length; i++) {
toastLog("方程组的解为x" + (i + 1) + "="+
df.format(result[i]));
}
}
} catch (e) {
// TODO Auto-generated catch block
log(e.name + ": " + e.message);
}
}
main();
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。