تحول هيلبرت

تحول هلبرت، بالأحمر، لـموجة مربعة، بالأزرق.

تحوّل هلبرت Hilbert transform هي عملية خطية تستخدم في الرياضيات وعمليات معالجة الإشارة. تحويل هيلبرت سمي بعدما قدم ديفيد هيلبرت حلا لحالة خاصة من معضلة ريمان وهيلبرت

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

جدول تحولات هلبرت مختارة

الإشارة
تحول هلبرت[fn 1]

[fn 2]
[fn 2]

(see Dawson function)
Sinc function
Rectangular function
Dirac delta function
Characteristic Function
ملاحظات
  1. ^ Some authors (e.g., Bracewell) use our −H as their definition of the forward transform. A consequence is that the right column of this table would be negated.
  2. ^ أ ب The Hilbert transform of the sin and cos functions can be defined in a distributional sense, if there is a concern that the integral defining them is otherwise conditionally convergent. In the periodic setting this result holds without any difficulty.

An extensive table of Hilbert transforms is available (King 2009). Note that the Hilbert transform of a constant is zero.


تحول هلبرت المتقطع

Figure 1: Filter whose frequency response is bandlimited to about 95% of the Nyquist frequency
Figure 2: Hilbert transform filter with a highpass frequency response
Figure 3.
Figure 4. The Hilbert transform of cos(ωt) is sin(ωt). This figure shows sin(ωt) and two approximate Hilbert transforms computed by the MATLAB library function, hilbert(·)
Figure 5. Discrete Hilbert transforms of a cosine function, using piecewise convolution

For a discrete function,    with discrete-time Fourier transform (DTFT),    and discrete Hilbert transform    the DTFT of    in the region  −π < ω <  π is given by:

The inverse DTFT, using the convolution theorem, is:

حيث:

which is an infinite impulse response (IIR). When the convolution is performed numerically, an FIR approximation is substituted for h[n], as shown in Figure 1. An FIR filter with an odd number of anti-symmetric coefficients is called Type III, which inherently exhibits responses of zero magnitude at frequencies 0 and Nyquist, resulting in this case in a bandpass filter shape. A Type IV design (even number of anti-symmetric coefficients) is shown in Figure 2. Since the magnitude response at Nyquist does not drop out, it approximates an ideal Hilbert transformer a little better than the odd-tap filter. However:

  • A typical (i.e. properly filtered and sampled) u[n] sequence has no useful components at the Nyquist frequency.
  • The Type IV impulse response requires a ½ sample shift in the h[n] sequence. That causes the zero-valued coefficients to become non-zero, as seen in Figure 2. So a Type III design is potentially twice as efficient as Type IV.
  • The group delay of a Type III design is an integer number of samples, which facilitates aligning with to create an analytic signal. The group delay of Type IV is halfway between two samples.

The MATLAB function, hilbert(u,N), convolves a u[n] sequence with the periodic summation:[1][2]

and returns one cycle (N samples) of the periodic result in the imaginary part of a complex-valued output sequence. The convolution is implemented in the frequency domain as the product of the array    with samples of the −i•sgn(ω) distribution (whose real and imaginary components are all just 0 or ±1). Figure 3 compares a half-cycle of hN[n] with an equivalent length portion of h[n]. Given an FIR approximation for   denoted by    substituting    for the −i•sgn(ω) samples results in an FIR version of the convolution.

The real part of the output sequence is the original input sequence, so that the complex output is an analytic representation of u[n]. When the input is a segment of a pure cosine, the resulting convolution for two different values of N is depicted in Figure 4 (red and blue plots). Edge effects prevent the result from being a pure sine function (green plot). Since hN[n] is not an FIR sequence, the theoretical extent of the effects is the entire output sequence. But the differences from a sine function diminish with distance from the edges. Parameter N is the output sequence length. If it exceeds the length of the input sequence, the input is modified by appending zero-valued elements. In most cases, that reduces the magnitude of the differences. But their duration is dominated by the inherent rise and fall times of the h[n] impulse response.

An appreciation for the edge effects is important when a method called overlap-save is used to perform the convolution on a long u[n] sequence. Segments of length N are convolved with the periodic function:

When the duration of non-zero values of    is M < N, the output sequence includes N − M + 1 samples of    M-1 outputs are discarded from each block of N, and the input blocks are overlapped by that amount to prevent gaps.

Figure 5 is an example of using both the IIR hilbert() function and the FIR approximation. In the example, a sine function is created by computing the Discrete Hilbert transform of a cosine function, which was processed in four overlapping segments, and pieced back together. As the FIR result (blue) shows, the distortions apparent in the IIR result (red) are not caused by the difference between h[n] and hN[n] (green and red in Fig 3). The fact that hN[n] is tapered (windowed) is actually helpful in this context. The real problem is that it's not windowed enough. Effectively, M = N, whereas the overlap-save method needs M < N.

انظر أيضا

مراجع

وصلات خارجية

  • Derivation of the boundedness of the Hilbert transform
  • Mathworld Hilbert transform — Contains a table of transforms
  • Analytic Signals and Hilbert Transform Filters
  • Eric W. Weisstein, Titchmarsh theorem at MathWorld.
  • Johansson, Mathias. "The Hilbert transform" (PDF). Archived from the original (PDF) on 2012-02-05. a student level summary of the Hilbert transformation.
  • "GS256 Lecture 3: Hilbert Transformation" (PDF). Archived from the original (PDF) on 2012-02-27. an entry level introduction to Hilbert transformation.
  1. ^ see Convolution Theorem
  2. ^ For the summed result obtained slightly differently see http://www.rle.mit.edu/dspg/documents/HilbertComplete.pdf eq. (17), (18), and unlabeled eq. below (18).
الكلمات الدالة: