Java实现希尔伯特(hilbert)变换

public static void FHT(double[] data, FourierTransform.Direction direction){
        
        int N = data.length;

        
        if (direction == FourierTransform.Direction.Forward){
            
            ComplexNumber[] cdata = new ComplexNumber[N];
            for (int i = 0; i < N; i++)
                    cdata[i] = new ComplexNumber(data[i], 0);

            FourierTransform.FFT(cdata, FourierTransform.Direction.Forward);

            for (int i = 1; i < (N/2); i++)
            {
                    cdata[i].real *= 2.0;
                    cdata[i].imaginary *= 2.0;
            }

           
            for (int i = (N/2)+1; i < N; i++)
            {
                    cdata[i].real = 0.0;
                    cdata[i].imaginary = 0.0;
            }
            FourierTransform.FFT(cdata, FourierTransform.Direction.Backward);

            // Convert back to our initial double array
            for (int i = 0; i < N; i++)
                    data[i] = Math.hypot(cdata[i].real,cdata[i].imaginary);
            
            }
        else {
            FHT(data, FourierTransform.Direction.Forward);

            for (int i = 0; i < data.length; i++)
                data[i] = -data[i];
        }
    }

对应是matlab中的hilbert函数

gitee仓库：HilbertTransform: Java实现希尔伯特(hilbert)变换,请自测，我是通过MATLAB测试的，计算结果能和MATLAB一致