java离散点数据积分的思路

积分一般是针对函数进行，现在只有离散点，需要将离散点通过插值或者拟合的方式变成函数。一重积分本质上就是求面积（函数的曲线和坐标轴围成的面积），最简单的方法是沿着坐标横轴将面积分成很多段梯形，用这些梯形的面积来代表需要求取的面积。

本例涉及到了对离散数据进行插值或者拟合，在下面的案例中并不会介绍这方面的内容，请自行百度或google.

先给出求取任何2个点之间积分的函数

/**
     * 计算任意2个点之间的积分（或者叫面积）
     * a和b之间的距离不要太大，如果太大的话，函数中变量n的值也需要相应变大，本例中n的值为10
     * 本例只准备用该函数计算2个相邻采样点之间的积分，所以n设置的比较小
     * @param samplingFun 需要积分的函数，通过拟合或者插值将原始的离散数据拟合成的函数
     * @param a 积分起始点
     * @param b 积分结束点
     * @return 积分的结果
     */
    private static double integralInner(SamplingFun samplingFun, double a, double b) {
       int n = 10;//每一段之间分为10个梯形
       double s = 0;
       double h = (b - a) / n;
       for(int i = 0; i < n; i++){
           s += (samplingFun.value(a + i * h) + samplingFun.value(a + (i + 1) * h)) / 2.0;
       }
       return s * h;
    }

 

对采样点的集合进行积分的函数

/**
     * 对采样数据的集合进行积分的函数，梯形法,实测是可用的
     * @param points 采样点的集合
     * @param samplingInterval 采样间隔
     */
    public static void integral2(List<SamplingPoint> points, int samplingInterval){

        //可用，且时间复杂度不高，因为会保存已经算过的面积
        double integral = 0;//保存已经计算过的面积
        SamplingFun samplingFun = new SamplingFun(samplingInterval, points);//通过插值的方式得到的函数
        for(int i = 0; i < points.size() - 1; i++){
            //0和1之间的面积属于1
            //1和2之间的面积 加上 已经算过的0和1之间的面积 等于0到2之间的面积，属于2
            //2和3之间的面积 加上 已经算过的0和2之间的面积 等于0到3之间的面积，属于3
            //最后一个是1022和1023之间的面积 加上已经算过的0和1022之间的面积 等于0和1023之间的面积，属于1023
            double integralRes = integralInner(samplingFun, i * samplingInterval * 1.0 / 1000000.0,
                    (i + 1) * samplingInterval * 1.0 / 1000000.0);
            points.get(i + 1).setRealY((float) (integral + integralRes));
            integral += integralRes;
        }
        points.get(0).setRealY(0);
    }