SOME APPLICATIONS OF THE HILBERT TRANSFORM

Similar documents
U.P.B. Sci. Bull., Series A, Vol. 74, Iss. 4, 2012 ISSN ALTERNATING -GROUPS. Ion ARMEANU 1, Didem OZTURK 2

U.P.B. Sci. Bull., Series A, Vol. 74, Iss. 3, 2012 ISSN SCALAR OPERATORS. Mariana ZAMFIR 1, Ioan BACALU 2

ANOVA IN THE EDUCATIONAL PROCESS

SDP APPROXIMATION OF THE HALF DELAY AND THE DESIGN OF HILBERT PAIRS. Bogdan Dumitrescu

Probability and Statistics for Final Year Engineering Students

ON SOME STATISTICAL INDICATORS OF THE TRAFFIC IN WIRELESS LOCAL AREA NETWORKS

AN ALTERNATIVE TO THE THEORY OF EXTREMA

Gabor wavelet analysis and the fractional Hilbert transform

Hilbert spaces. Let H be a complex vector space. Scalar product, on H : For all x, y, z H and α, β C α x + β y, z = α x, z + β y, z x, y = y, x

Review of Concepts from Fourier & Filtering Theory. Fourier theory for finite sequences. convolution/filtering of infinite sequences filter cascades

Signal interactions Cross correlation, cross spectral coupling and significance testing Centre for Doctoral Training in Healthcare Innovation

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

ENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University

Fourier Series. Spectral Analysis of Periodic Signals

9 Brownian Motion: Construction

ENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University

Vectors in Function Spaces

Shannon meets Wiener II: On MMSE estimation in successive decoding schemes

Recall that any inner product space V has an associated norm defined by

The Continuous-time Fourier

WAVELETS, PROPERTIES OF THE SCALAR FUNCTIONS

ω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the

Deterministic. Deterministic data are those can be described by an explicit mathematical relationship

Rings With Topologies Induced by Spaces of Functions

On Riesz-Fischer sequences and lower frame bounds

Signals and Spectra - Review

A GENERALIZATION OF A CLASSICAL MONTE CARLO ALGORITHM TO ESTIMATE π

Boundedly complete weak-cauchy basic sequences in Banach spaces with the PCP

L p -boundedness of the Hilbert transform

An Introduction to Filterbank Frames

n WEAK AMENABILITY FOR LAU PRODUCT OF BANACH ALGEBRAS

TWO BOUNDARY ELEMENT APPROACHES FOR THE COMPRESSIBLE FLUID FLOW AROUND A NON-LIFTING BODY

Decomposing Bent Functions

SOME NUMERICAL ASPECTS IN MODELING THE LONGITUDINAL IMPACT OF TWO SHAPE MEMORY ALLOY BARS

The Heat Equation John K. Hunter February 15, The heat equation on a circle

On the Hilbert Transform of Wavelets

Valuations. 6.1 Definitions. Chapter 6

A Riesz basis of wavelets and its dual with quintic deficient splines

On Cryptographic Properties of the Cosets of R(1;m)

Equivalent Multivariate Stochastic Processes

Lecture 16: Bessel s Inequality, Parseval s Theorem, Energy convergence

The Feynman-Kac formula

ECE-314 Fall 2012 Review Questions for Midterm Examination II

Machine Learning. A Bayesian and Optimization Perspective. Academic Press, Sergios Theodoridis 1. of Athens, Athens, Greece.

APPROXIMATE ISOMETRIES ON FINITE-DIMENSIONAL NORMED SPACES

ESTIMATES ON THE NUMBER OF LIMIT CYCLES OF A GENERALIZED ABEL EQUATION

Solutions to Problems in Chapter 4

ATTENUATION OF THE ACOUSTIC SCREENS IN CLOSED SPACES

System Identification & Parameter Estimation

Lecture Notes 7 Stationary Random Processes. Strict-Sense and Wide-Sense Stationarity. Autocorrelation Function of a Stationary Process

COMPLEX SIGNALS are used in various areas of signal

A i A j = 0. A i B j = δ i,j. A i B j = 1. A i. B i B j = 0 B j

arxiv: v1 [math.oa] 2 Mar 2014

Complex Analysis Qualifying Exam Solutions

The Hilbert Transform

The Hilbert transform

TIME-FREQUENCY ANALYSIS AND HARMONIC GAUSSIAN FUNCTIONS

= F (b) F (a) F (x i ) F (x i+1 ). a x 0 x 1 x n b i

ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process

Homework If the inverse T 1 of a closed linear operator exists, show that T 1 is a closed linear operator.

Outline of Fourier Series: Math 201B

COMPUTATION OF FOURIER TRANSFORMS FOR NOISY B FOR NOISY BANDLIMITED SIGNALS

A NOTE ON ENRICHED CATEGORIES

Math 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.

AN APPROACH TO THE NONLINEAR LOCAL PROBLEMS IN MECHANICAL STRUCTURES

Lecture 27 Frequency Response 2

Continuous Time Signal Analysis: the Fourier Transform. Lathi Chapter 4

Blind separation of instantaneous mixtures of dependent sources

Review of Probability

Question Paper Code : AEC11T02

Fourier Analysis and Power Spectral Density

NYQUIST STABILITY FOR HYSTERESIS SWITCHING MODE CONTROLLERS

A brief introduction to trace class operators

ON THE BOUNDEDNESS BEHAVIOR OF THE SPECTRAL FACTORIZATION IN THE WIENER ALGEBRA FOR FIR DATA

Optimum Sampling Vectors for Wiener Filter Noise Reduction

FUNCTIONAL ANALYSIS LECTURE NOTES: ADJOINTS IN HILBERT SPACES

WE begin by briefly reviewing the fundamentals of the dual-tree transform. The transform involves a pair of

Stochastic Process II Dr.-Ing. Sudchai Boonto

Hilbert Spaces. Contents

Nonlinear Stationary Subdivision

Applied and Computational Harmonic Analysis

ECE-340, Spring 2015 Review Questions

Let R be the line parameterized by x. Let f be a complex function on R that is integrable. The Fourier transform ˆf = F f is. e ikx f(x) dx. (1.

Representation of Signals & Systems

LECTURE-20 : ISOLATED SINGULARITIES. A punctured domain is an open set with a point removed. For p Ω, we use the notation

So reconstruction requires inverting the frame operator which is often difficult or impossible in practice. It follows that for all ϕ H we have

MULTIPLEXING AND DEMULTIPLEXING FRAME PAIRS

Czechoslovak Mathematical Journal

General Inner Product and The Fourier Series

NOTES ON FRAMES. Damir Bakić University of Zagreb. June 6, 2017

Study of some equivalence classes of primes

SMALL COMBINATORIAL CARDINAL CHARACTERISTICS AND THEOREMS OF EGOROV AND BLUMBERG

1 Mathematical Preliminaries

FORMAL GROUPS OF CERTAIN Q-CURVES OVER QUADRATIC FIELDS

Decompositions of frames and a new frame identity

Spectral Analysis. Jesús Fernández-Villaverde University of Pennsylvania

Modeling with Itô Stochastic Differential Equations

EAS 305 Random Processes Viewgraph 1 of 10. Random Processes

arxiv: v1 [math.ca] 6 Feb 2015

An Investigation of 3D Dual-Tree Wavelet Transform for Video Coding

Transcription:

U.P.B. Sci. Bull. Series A, Vol. 71, Iss. 3, 2009 ISSN 1223-7027 SOME APPLICATIONS OF THE HILBERT TRANSFORM Ştefania Constantinescu 1 În acest articol, se dau unele proprietăţi ale transformării Hilbert, cu interpretări în teoria matematică a semnalelor aleatoare. Este studiată reprezentarea semnalelor analitice ca o generalizare a formei complexe a semnalelor sinusoidale, un rol special fiind acordat proprietăţilor spectrale, legate de conceptul de densitate de putere. Este de asemenea analizat un exemplu practic din domeniul circuitelor electronice. In this paper we recall some properties of the Hilbert transform interpreted in terms of the mathematical theory of the random signals. We study the representation of the analitical signals as a generalization of the complex form of the sinusoidal signals, a special part being focused on the spectral properties, in connection with the concept of power density. It is also presented an example taken from practice of the electronic circuits. Keywords: Fourier transform, Hilbert transform, analytical representation. 1. Introduction By making use of the Hilbert transform, one can describe some constructions used in the signal theory in mathematical rigorous terms. In this paper, one recalls the main facts regarding the Hilbert transform and one analyzes some concepts like that of analytic associated signal, convolution filters, spectral densities of stationary random signal. 2. On the Hilbert transform Let be a function from. For any, one can consider the convolution, which belongs to. If the Fourier image is bounded, then the linear operator : is well defined and bounded; indeed, if, then. and if (0 fixed), then by the Parseval relation,. and since is dense in, this inequality extends to. 1 Drd., Dept. of Mathematics, University POLITEHNICA of Bucharest, Romania

36 Ştefania Constantinescu Take now (0. This function is not integrable on ; instead, its primitive is and one can take its associated distribution and then, the distribution, defined by 1.. principal value for any test-function. One knows that. 1, so can be identified with the distribution [2]. For any function : from, one defines its Hilbert image. (1) Explicitly, for any,.. 1. Since., from (1) one gets:... (2) Because, it follows that is a function which belongs to. PROPOSITION 1. Suppose that : belongs to and 0 for any 0. If,, then is even and is odd; moreover,. (3) Proof. Indeed, since 0,, then, for any ; hence. But, therefore. Thus, and it remains to separate the real and imaginary parts. One knows that the Hilbert operator :, 1 1 is linear and continuous. Moreover, for any,, that is. Particularly, is bijective and. From the relation (2), one obtains another definition of the Hilbert transform:.. One also knows the following table:

Some applications of the Hilbert transform 37 sin cos cos sin 1 According to prop. 1, one gets: and. Such a pair of Hilbert images can also be obtained in another way. Namely, let : be a function from and put 1 cos, 1 sin. By the Fourier representation formula, 1 2 cos cos sin, for any real. Consider now the conjugate-function, defined by sin cos. Then sin lim sin lim lim lim... Thus, and similarly,. NOTE. One can also define the discrete Hilbert transform. Namely, consider the sequence,, defined by if is odd and 0 for even. For any sequence,, such that 0 for 0 and, its discrete Hilbert image is, where (, in analogy with (1)). If and put,, and and it founds again (3); [5]

38 Ştefania Constantinescu 3.. Some representations of the signals It is useful to introduce some terminology from Signal Processing. Recall that any function : from is also called a signal with finite energy; its norm is the energy of. Any linear continuous operator : is also called a continual filter. If there is a function (called weight) such that for any, then one says that is a convolution filter, with the transfer function. By definition, if : belongs to, one can consider a new signal, defined by,, (4) called the analytic signal associated to [1], [4]. Mathematically, this terminology is justified by the fact that under some conditions, is the restriction of a complex holomorphic ( analytic) function. Example. For cos, one gets cos sin and for sin, sincos. By extension, one defines the distribution. PROPOSITION 2. The filter, is a convolution filter, with the transfer function 2, 0 (5) 0, 0 Proof. For any, we have 1 1 1, hence. Thus, is a convolution filter. Its transfer function is 1 1. 1, for any real 0. COROLLARY. For any and any real,., i.e... NOTE. The passage from to is a generalization of the complex form of the sinusoidal signals. By (4), the imaginary part of is just the Hilbert image of and conversely, if we know, then the signal can be directly recovered, since.

Some applications of the Hilbert transform 39 Fix 0. For any signal : from, one can consider,,, defined by. ;, for ant real t. NOTE. If :, is a convolution filter with the transfer function, then by proposition 2, the filter. ( 0 fixed) will have the transfer function 2, 0. 0, 0 Thus, if is symmetrical around (i.e. null outside an interval centered at ), then is even; the filter is said band-pass and low-pass. Therefore, by considering the Hilbert transform, the band-pass filters can be replaced by low-pass ones. From (4), it follows that. cos sin Hence, from and, one can recover ; conversely,. cos.sin and sin.cos. These formulas are different ways of signal representations, which use the Hilbert transforms. In what follows, we will estimate their spectral densities and mean powers. 4. Spectral densities estimations Let Ω,, be a probability field and, be a stationary random signal ( random process) relative to this field; for any, is supposed to be a random variable from Ω, hence having means and variances. Then one can define its autocorrelative function., for any real ; (6) this is independent in, one to the stationarity hypothesis. Obviously, 0 and is even. The Fourier image of, is called the power spectral density of the random signal. By the Fourier inversion formula,. (7) Since is even, 2 is real. cos and for any real,

40 Ştefania Constantinescu PROPOSITION 3. For any stationary random signal, the power spectral density of the analytic signal is 4, 0. (8) 0, 0 Proof. Generally, if :,, is a convolution filter with the weight hence, then it holds the following relation for the autocorrelations of inputs/outputs: and for spectral densities,, for any real ; here is the transfer function of [1]. Applying this for the filter from prop.2, it follows that. and apply the proposition 2. In analogy to proposition 3, we will give the power spectral densities for, and. PROPOSITION 4. Fix 0 and a stationary random signal in Ω. Then a) ; b) and if is symmetrical around, i.e. null outside an interval centered in, this equals. Proof. a)........ Then.. b) On the other hand,....... By direct computation,, hence and.., therefore, 2 2.. If is symmetrical around, then by a), is symmetrical around 0, hence an even function. Then has only real values and consequently,. 0 for any and 2. Example. Fix 0, 0 and suppose that, is an ideal random signal with the power spectral density,, 0, In this case, one says that is a white noise with the frequency band width equal to 2B. By prop. 4, we can explicit the mean powers for,, and. Namely,

Some applications of the Hilbert transform 41 0 1 2 1 2 1.. Similarly,. and also,.. An application Consider an electronic circuit like that from figure 1, where the input signal is the tension, which can be regarded as a stationary random signal with the power spectral density. Denote by the tension at the output. By Kirchhoff law,, where the current intensity is.. Thus for any real. By applying the Fourier transform, it follows that., hence the transfer function of the circuit, regarded as a filter, will be 1 1. If, 0, where, it follows that 0, 0. Moreover,, hence. For, is neglectable and so will be. The previous circuit is an example of low-pass filter. 5. Conclusions Fig. 1. In 1 we have recalled some properties of the Hilbert transform. In 2 we analyzed the notion of analytic signal associated to a given signal and we introduced a modulation around a fixed frequency. This assures a correspondence between band-pass and low-pass filters. The filter plays a special part (prop. 2, 3). Another concept which is analyzed in 3 is that of power spectral density; this permits to estimate the mean powers of

42 Ştefania Constantinescu some signals and to study the behavior of some filters (prop. 4). In this work, we presented some concepts taken from Electrical Engineering and we gave some significant examples. R E F E R E N C E S [1] V.Croitoru, O. Stănăşilă, Probabilităţi aplicate şi semnale aleatoare (Applied probabilities and random signals), Ed. Matrix-Rom, 2007 [2] C. Gasquet, P. Witomski, Analyse de Fourier et applications, Masson, 1990 [3] N.G. Kingsbury, A dual-tree complex wavelet transform, Proc. IEEE Conf. Image Processing, Vancouver, Sept. 2000 [4] S.L. Marple, Computing the discrete-time analytic signal, IEEE Trans. of Signal Processing, vol. 47, nr. 9, 2600-2603, 2002 [5] I.W.Selesnick, Hilbert pairs of wavelet bases, Internet, NSF Career, 987452, 2004 [6] S. Tolwinski, The Hilbert transform and data analysis, RIG Projct 2007