本文實(shí)例講述了Java實(shí)現(xiàn)的n階曲線擬合功能。分享給大家供大家參考，具體如下：
前面一篇文章Java實(shí)現(xiàn)求解一元n次多項式的方法，能解多項式以后，還需要利用那個類，根據(jù)若干采樣點(diǎn)數(shù)據(jù)來對未來數(shù)據(jù)進(jìn)行預(yù)測，擬合的矩陣在上一篇文章中已經(jīng)貼出來了，這里就不說了，本篇主要是如何根據(jù)采樣點(diǎn)來計算系數(shù)矩陣，并計算預(yù)測點(diǎn)的值。
原理很簡單，公式在上一篇文章中也有了，此處直接貼代碼。
其中用到了上一篇文章中寫的類commonAlgorithm.PolynomiaSoluter
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
|
package commonAlgorithm;import commonAlgorithm.PolynomialSoluter;import java.lang.Math;public class LeastSquare { private double[][] matrixA; private double[] arrayB; private double[] factors; private int order; public LeastSquare() { } /* * 實(shí)例化后，計算前，先要輸入?yún)?shù)并生成公式 arrayX為采樣點(diǎn)的x軸坐標(biāo)，按照采樣順序排列 * arrayY為采樣點(diǎn)的y軸坐標(biāo)，按照采樣順序與x一一對應(yīng)排列 order * 為進(jìn)行擬合的階數(shù)。用低階來擬合高階曲線時可能會不準(zhǔn)確，但階數(shù)過高會導(dǎo)致計算緩慢 */ public boolean generateFormula(double[] arrayX, double[] arrayY, int order) { if (arrayX.length != arrayY.length) return false; this.order = order; int len = arrayX.length; // 擬合運(yùn)算中的x矩陣和y矩陣 matrixA = new double[order + 1][order + 1]; arrayB = new double[order + 1]; // 生成y矩陣以及x矩陣中冪<=order的部分 for (int i = 0; i < order + 1; i++) { double sumX = 0; for (int j = 0; j < len; j++) { double tmp = Math.pow(arrayX[j], i); sumX += tmp; arrayB[i] += tmp * arrayY[j]; } for (int j = 0; j <= i; j++) matrixA[j][i - j] = sumX; } // 生成x矩陣中冪>order的部分 for (int i = order + 1; i <= order * 2; i++) { double sumX = 0; for (int j = 0; j < len; j++) sumX += Math.pow(arrayX[j], i); for (int j = i - order; j < order + 1; j++) { matrixA[i - j][j] = sumX; } } // 實(shí)例化PolynomiaSoluter并解方程組，得到各階的系數(shù)序列factors PolynomialSoluter soluter = new PolynomialSoluter(); factors = soluter.getResult(matrixA, arrayB); if (factors == null) return false; else return true; } // 根據(jù)輸入坐標(biāo)，以及系數(shù)序列factors計算指定坐標(biāo)的結(jié)果 public double calculate(double x) { double result = factors[0]; for (int i = 1; i <= order; i++) result += factors[i] * Math.pow(x, i); return result; }} |
希望本文所述對大家java程序設(shè)計有所幫助。
原文鏈接：http://blog.csdn.net/strangerzz/article/details/45250063
