這篇文章給大家介紹使用Java怎么求解一元n次多項(xiàng)式,內(nèi)容非常詳細(xì),感興趣的小伙伴們可以參考借鑒,希望對(duì)大家能有所幫助。
在巴東等地區(qū),都構(gòu)建了全面的區(qū)域性戰(zhàn)略布局,加強(qiáng)發(fā)展的系統(tǒng)性、市場(chǎng)前瞻性、產(chǎn)品創(chuàng)新能力,以專(zhuān)注、極致的服務(wù)理念,為客戶提供成都網(wǎng)站設(shè)計(jì)、網(wǎng)站建設(shè) 網(wǎng)站設(shè)計(jì)制作按需設(shè)計(jì)網(wǎng)站,公司網(wǎng)站建設(shè),企業(yè)網(wǎng)站建設(shè),高端網(wǎng)站設(shè)計(jì),營(yíng)銷(xiāo)型網(wǎng)站建設(shè),外貿(mào)網(wǎng)站制作,巴東網(wǎng)站建設(shè)費(fèi)用合理。
項(xiàng)目需要做趨勢(shì)預(yù)測(cè),采用線性擬合、2階曲線擬合和指數(shù)擬合的算法,各種線性擬合算法寫(xiě)成矩陣大概是這么個(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ù),自己寫(xiě)了個(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; } }
關(guān)于使用Java怎么求解一元n次多項(xiàng)式就分享到這里了,希望以上內(nèi)容可以對(duì)大家有一定的幫助,可以學(xué)到更多知識(shí)。如果覺(jué)得文章不錯(cuò),可以把它分享出去讓更多的人看到。
網(wǎng)站標(biāo)題:使用Java怎么求解一元n次多項(xiàng)式
分享網(wǎng)址:http://chinadenli.net/article10/gogjgo.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供外貿(mào)網(wǎng)站建設(shè)、小程序開(kāi)發(fā)、響應(yīng)式網(wǎng)站、電子商務(wù)、做網(wǎng)站、網(wǎng)站維護(hù)
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如需處理請(qǐng)聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來(lái)源: 創(chuàng)新互聯(lián)