java语言实现一个长度为n_Java语言实现求解一元n次多项式的方法示例[Java代码]...
本文主要向大家介紹了Java語(yǔ)言實(shí)現(xiàn)求解一元n次多項(xiàng)式的方法示例,通過(guò)具體的內(nèi)容向大家展示,希望對(duì)大家學(xué)習(xí)JAVA語(yǔ)言有所幫助。
項(xiàng)目需要做趨勢(shì)預(yù)測(cè),采用線性擬合、2階曲線擬合和指數(shù)擬合的算法,各種線性擬合算法寫成矩陣大概是這么個(gè)形式:
其中x是橫坐標(biāo)采樣值,y是縱坐標(biāo)采樣值,i是采樣點(diǎn)序列號(hào),a是系數(shù),N是采樣點(diǎn)個(gè)數(shù),n是階數(shù),所以線性擬合最后就轉(zhuǎn)成了一個(gè)解高階方程組的問(wèn)題。
不知道有沒(méi)有什么好用的java矩陣運(yùn)算的包,我很不擅長(zhǎng)搜集這種資料,所以只好撿起了已經(jīng)放下多年的線性代數(shù),自己寫了個(gè)java程序用增廣矩陣的算法來(lái)解高階方程組。直接貼代碼好了:
package commonAlgorithm;
public class PolynomialSoluter {
private double[][] matrix;
private double[] result;
private int order;
public PolynomialSoluter() {
}
// 檢查輸入項(xiàng)長(zhǎng)度并生成增廣矩陣
private boolean init(double[][] matrixA, double[] arrayB) {
order = arrayB.length;
if (matrixA.length != order)
return false;
matrix = new double[order][order + 1];
for (int i = 0; i < order; i++) {
if (matrixA[i].length != order)
return false;
for (int j = 0; j < order; j++) {
matrix[i][j] = matrixA[i][j];
}
matrix[i][order] = arrayB[i];
}
result = new double[order];
return true;
}
public double[] getResult(double[][] matrixA, double[] arrayB) {
if (!init(matrixA, arrayB))
return null;
// 高斯消元-正向
for (int i = 0; i < order; i++) {
// 如果當(dāng)前行對(duì)角線項(xiàng)為0則與后面的同列項(xiàng)非0的行交換
if (!swithIfZero(i))
return null;
// 消元
for (int j = i + 1; j < order; j++) {
if (matrix[j][i] == 0)
continue;
double factor = matrix[j][i] / matrix[i][i];
for (int l = i; l < order + 1; l++)
matrix[j][l] = matrix[j][l] - matrix[i][l] * factor;
}
}
// 高斯消元-反向-去掉了冗余計(jì)算
for (int i = order - 1; i >= 0; i--) {
result[i] = matrix[i][order] / matrix[i][i];
for (int j = i - 1; j > -1; j--)
matrix[j][order] = matrix[j][order] - result[i] * matrix[j][i];
}
return result;
}
private boolean swithIfZero(int i) {
if (matrix[i][i] == 0) {
int j = i + 1;
// 找到對(duì)應(yīng)位置非0的列
while (j < order && matrix[j][i] == 0)
j++;
// 若對(duì)應(yīng)位置全為0則無(wú)解
if (j == order)
return false;
else
switchRows(i, j);
}
return true;
}
private void switchRows(int i, int j) {
double[] tmp = matrix[i];
matrix[i] = matrix[j];
matrix[j] = tmp;
}
}
復(fù)制代碼
有更好的算法或者有合適的矩陣運(yùn)算包歡迎交流
PS:這里再為大家推薦幾款計(jì)算工具供大家進(jìn)一步參考借鑒:
在線一元函數(shù)(方程)求解計(jì)算工具:
http://tools.jb51.net/jisuanqi/equ_jisuanqi
科學(xué)計(jì)算器在線使用_高級(jí)計(jì)算器在線計(jì)算:
http://tools.jb51.net/jisuanqi/jsqkexue
在線計(jì)算器_標(biāo)準(zhǔn)計(jì)算器:
http://tools.jb51.net/jisuanqi/jsq
本文由職坐標(biāo)整理并發(fā)布,希望對(duì)同學(xué)們有所幫助。了解更多詳情請(qǐng)關(guān)注編程語(yǔ)言JAVA頻道!
總結(jié)
以上是生活随笔為你收集整理的java语言实现一个长度为n_Java语言实现求解一元n次多项式的方法示例[Java代码]...的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: python入门及日常应用_python
- 下一篇: 民生visa全币种信用卡怎么用最划算