Issue with DSP.Util.hilbert

Specific Domains Modelling & Simulations
1

I was looking to get the envelope of a signal and looking at the example at Discrete-time analytic signal using Hilbert transform - MATLAB hilbert and checking this in octave using the signal package from forge I am getting different results from Matlab/Octave to Julia.

Ex (Matlab):

xr = [1 2 3 4];
x = hilbert(xr)
x = 1×4 complex

   1.0000 + 1.0000i   2.0000 - 1.0000i   3.0000 - 1.0000i   4.0000 + 1.0000i

Ex (Julia):

julia> xr = [1 2 3 4];

julia> x = Util.hilbert(xr)
1×4 Array{Complex{Float64},2}:
 1.0+0.0im  2.0+0.0im  3.0+0.0im  4.0+0.0im

Any guessing to whats going on? I was using DSP.jl before the 1.0v release of Julia and it was fine.

2

It looks like passing an integer array is the problem.
I really think that’s buggy behavior, might be worth filing an issue.

julia> x = [ 1.0, 2.0, 3.0, 4.0 ]
4-element Array{Float64,1}:
 1.0
 2.0
 3.0
 4.0

julia> hilbert(x)
4-element Array{Complex{Float64},1}:
 1.0 + 1.0im
 2.0 - 1.0im
 3.0 - 1.0im
 4.0 + 1.0im