java dsp 高通滤波

以下是使用Java实现高通滤波的示例代码：

import org.apache.commons.math3.filter.DefaultMeasurementModel;
import org.apache.commons.math3.filter.DefaultProcessModel;
import org.apache.commons.math3.filter.KalmanFilter;
import org.apache.commons.math3.filter.MeasurementModel;
import org.apache.commons.math3.filter.ProcessModel;

public class HighPassFilter {
    public static void main(String[] args) {
        // 定义高通滤波器的参数
        double dt = 0.1; // 采样时间间隔
        double cutoffFrequency = 1.0; // 截止频率

        // 创建过程模型
        ProcessModel processModel = new DefaultProcessModel(
                new double[]{1.0}, // 状态转移矩阵A
                null, // 控制输入矩阵B
                null, // 过程噪声矩阵Q
                null, // 初始状态矩阵x0
                null // 初始状态协方差矩阵P0
        );

        // 创建测量模型
        MeasurementModel measurementModel = new DefaultMeasurementModel(
                new double[]{1.0}, // 观测矩阵H
                null // 观测噪声矩阵R
        );

        // 创建卡尔曼滤波器
        KalmanFilter kalmanFilter = new KalmanFilter(processModel, measurementModel);

        // 设置高通滤波器的截止频率
        double alpha = Math.exp(-2.0 * Math.PI * cutoffFrequency * dt);
        kalmanFilter.setAlpha(alpha);

        // 输入观测数据
        double[] measurements = {1.0, 2.0, 3.0, 4.0, 5.0};

        // 进行滤波
        for (double measurement : measurements) {
            kalmanFilter.predict();
            kalmanFilter.correct(new double[]{measurement});
            double[] filteredState = kalmanFilter.getStateEstimation();
            System.out.println("Filtered value: " + filteredState[0]);
        }
    }
}

请注意，上述示例使用了Apache Commons Math库中的KalmanFilter类来实现高通滤波。你需要在项目中添加对该库的依赖。