Medical, Pharma, Engineering, Science, Technology and Business

**Edris Mohsen ^{*}, Lyndon J Brown and Jie Chen**

Department of Electrical and Computer Engineering Western University, London, Ontario, Canada

- *Corresponding Author:
- Edris Mohsen

Department of Electrical and Computer Engineering Western University

London, Ontario, Canada

**Tel:**+1 519-661-2111

**E-mail:**[email protected]

**Received date:** July 14, 2017; **Accepted date:** July 25, 2017; **Published date:** July 27, 2017

**Citation: **Mohsen E, Brown LJ, Chen J (2017) A Real time Alternative to the Hilbert Huang Transform Based on Internal Model Principle. J Electr Electron Syst 6: 234. doi: 10.4172/2332-0796.1000234

**Copyright:** © 2017 Mohsen E, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Electrical & Electronic Systems

This article presents a new tuning approach for an adaptive internal-model-principle based signal identification algorithm whose computational costs are low enough to allow a realtime implementation. The algorithm allows an instantaneous Fourier decomposition of non-stationary signals that have a strongly predictable component. The algorithm is implemented as a feedback loop resulting in a closed loop system with a frequency response of a bandpass filter with notches at the frequencies of the Fourier decomposition. This is achieved through real time selection of the coefficients of the transfer functions in the feedback loop. Previously these coefficients were selected by solving a large set of coupled linear equations. Rules for explicitly solving for these parameters are given that only involve evaluating frequency responses at the frequencies of the instantaneous Fourier decomposition. This allows realtime implementation on a low cost lap top with sampling rates up to 10 kHz.

Internal model principle; Frequency identification; Adaptive multiple notch filters; Periodic disturbance; State variables; Bandpass filter; Instantaneous Fourier decomposition

In this article, we are interested in the problem of identifying signals of the following form

(1)

(2)

and n(t) is measurement noise. These are signals that are the sums of n periodic components with each component composed of mi harmonics. The periods, the harmonic amplitudes and relative phases can vary slowly in time. By identification, we mean determining the values *ω _{i}, Ä�_{ij}* and

Several techniques have been developed in the literature to solve this problem. The most traditional technique is the fast Fourier transform. Newer techniques include wavelet analysis. These approaches suffer from not allowing continuous estimations of the frequencies and have difficult trade-offs between time and frequency resolutions. Other approaches are based on the use of adaptive notch filters [1] and output regulation [2]. A new approach that has been widely applied is the Hilbert Huang Transform (HHT) [3]. Control engineers treat similar problems where exact tracking of reference signals or rejection of disturbances is required. Approaches that accomplish this include repetitive controllers [4] and adaptive feed-forward cancellation (AFC) [5]. The repetitive controller is based on a fundamental control theory principle called the internal model principle (IMP). This principle was presented by Francis and Wonham and states that the output error can be driven asymptotically to zero by placing a model of exogenous signals in a stable feedback loop [6]. Unfortunately small errors in this model can lead to significant degradation in the performance of internal model principle controllers. This problem of uncertainty in the signal model can be overcome with adaptive controllers [7]. In achieving asymptotically perfect rejection of disturbances it is inherent that the disturbance is completely identified. Thus, these types of controllers can be turned into signal processing algorithms by replacing the process to be controlled with tuning functions [8].

Unfortunately, to successfully implement this algorithm requires being able to tune a stable feedback control loop for the entire range of possible frequencies in the model given by equation (1). Fortunately, it has been shown that in the signal processing framework, the simplest tuning solution, i.e. selecting all of the gains to be one, is guaranteed to be stable. This algorithm has been successfully applied to the problem of the repeatable disturbances seen in disk drive head control [9]. Unfortunately, by resorting to this simple tuning approach, there is no control over the dynamics and noise rejection characteristics of the algorithm.

When the frequencies are known a priori, the report [10] shows how the dynamics of the algorithm can be completely specified. Unfortunately this article requires solving a set more than coupled linear equations which are a function of the signal’s frequencies. Unless the sample rate is less than 1Hz this will not be feasible to do each sample. This article shows how these parameters can be explicitly solved by simply evaluating some frequency response functions at certain frequencies.

In Section II, an instantaneous Fourier decomposition (IFD) algorithm [11] that is similar in approach to the HHT is presented. In Section III an updated formula for calculating the instantaneous frequencies are given. In Section IV, the new realtime tuned algorithm is presented. In Section V, the ability of the proposed algorithm to identify the periodic signal with uncertain frequencies is demonstrated. Conclusions are drawn in Section VI.

A preliminary version of this article was presented at the 30th annual IEEE Canadian Conference on Electrical and Computer Engineering (IEEE 2017 CCECE) in Windsor [12].

The HHT proceeds from the realization that the Hilbert transform gives a mathematically precise definition of instantaneous frequency that agrees with our intuitive understanding when applied to narrowband signals. In this narrowband case, the instantaneous frequency can be approximated as the derivative of the angle of the narrowband signal and times the quadrature of that signal where the quadrature can be approximated by either a scaled version of the derivative or integral of the signal. The HHT uses an empirical method to break down signals into narrowband signals. This empirical method is numerically intensive and not compatible with a realtime implementation.

Our algorithm uses the same approximations to estimate the instantaneous frequencies as the HHT but uses an alternative, notch filter based approach that simultaneous calculates the quadrature signals and decomposes the signal into narrow band signals. The structure of the adaptive instantaneous frequency decomposition is shown in **Figure 1**, where G(s) is a tuning function.

Each of the transfer functions IMi,j are an internal model for a sinusoid of frequency . When the model frequencies and the signal frequencies match, i.e., and the closed loop system is stable, each *u _{ij}* will be a single sinusoidal and meet the HHT definition of an intrinsic function. The basic algorithm is the state space based implementation of the internal models given by

(3)

where i=1,2,……, n and j=1,2,……, m_{i}. A_{ij}, is expressed as follows

(4)

This is taken from [11] with minor modifications to fit the signal model that was given in equations (1) and (2). The gains *K _{1ij}, K_{2ij}* have been moved to the input vector from the output vector so that adjustments in their value do not directly change u, i.e., a bumpless transfer. Consequently, the responses at

(5)

(6)

i.e., the first state is the sinusoidal component of the original signal and the second state is its quadrature. While the states are time varying, when the signal parameters are time varying, when the signal parameters are time invariant is time invariant as is

Since the state variables *x _{1i1}* and

(7)

can be used to update the frequency estimates.

Thus a quasi-periodic signal can be decomposed into a sum of narrow band signals, {*u _{ij}*}={

(8)

(9)

In ref. [11], it is establish for sufficiently small K_{ai} the algorithm is locally exponentially stable when G(s) and the *K _{1ij}, K_{2ij}* are chosen so that the feedback loop in

**Off-line tuning**

As with ref. [10], we satisfy the above stability assumption by designing the closed loop system to incorporate a bandpass filter with notch filter. Let a 2^{nd} order desirable bandpass filter be given by

(10)

We choose the controller parameters to be such that the transfer function from d to e is

(11)

where are small real numbers, and are the notches frequency. The presence of the numerator of the second term is a fundamental consequence of the internal model principle. Therefore, the ability of the algorithm to improve noise rejection is achieved.

An analysis of **Figure 1** gives

(12)

Where,

And

(13)

Note in equations (11,12) Ð� represents and Æ© represents . The terms *Y _{kl}* are the product of all the terms except the

**On-line frequency identification**

Now the crucial question is how to choose *G(s)* and *K _{1jk}, K_{2jk}* and implement the algorithm without needing to solve a set or 2nt+4 linear equations. It can be seen that all of the terms in the denominator except the term containing

(14)

This generates 2 complex and complementary conjugate equations with 2 unknowns, i.e., the real part of either equation gives *K _{1jk}* and the imaginary gives

Extending this to the following product

and equating the following sets , we get

**Linear dependency of equations:** When when i ≠ k then the equations to be solved become linearly dependent. With our solution technique this is reflected in the fact that the denominator of equation (12) will be zero when we substitute in and it will not be possible to calculate two pairs of internal model gains. Further, while it is theoretical possible to solve when the frequencies are extremely close, we get solutions that lead to unstable results because of numerical stability issues. To solve this problem, while calculating the controller gains, we drop the approximately redundant internal model when the frequencies become close, i.e., within 0.1%. After calculating the controller gains, the two redundant models are each assigned half of the gain. That is when jωi=lωk, we drop Internal model *IM _{l,k}* from the design stage. Let be the calculated controllers gains. Then and. It should be noted that the threshold for HHT to distinguish between close frequencies is 10%.

In this particular section, the effectiveness of our real time implementation of our proposed adaptive algorithm is verified via simulation. The model configuration parameters that are used with the matlab/simulink (R2016) environment are as follows: Solver ode5 (Dormand-prince) selection with fundamental sample time is 0.0025 s. Therefore, the sampling rate in our case is selected to be 400Hz, then the Nyquist frequency is 200 Hz. The code generation with C language and tool chain (Microsoft visual C++ 2012 V11.1 n-make 164-bit windows). All random numbers ware zero mean.

Our signal to be identified was produced by summing the outputs of two copies of the model shown in **Figure 2**. The feedback loop containing the pure delay is called a repetitive controller and is capable of producing any periodic disturbance with period T. The value T was an integrated band limited white disturbance. The frequency cutoff of this noise was 20 rad/s and the variance was 0.5. The initial conditions for both fundamental frequencies are 4.2 and 5 Hz. The disturbance input to the repetitive controller causes the amplitudes and relative phases to vary slowly with time as well. This random signal was band limited to 50 Hz and had variance 0.1. Additional measurement noise was added to the sums of these two signals. This noise was band limited to 50 Hz and had a variance of 0.1. The low pass filter had a cutoff frequency of 100 rad/s concentrating the energy in the harmonics to below the 4th and third harmonic, respectively though signal was present in all harmonics up to the Nyquist frequency.

The frequency adaption gains were chosen as K_{a}=1.95 or with frequency 7.5% to 10% of the fundamental frequencies (**Table 1**). The closed loop transfer function was chosen to be a second order Chebyshev band-pass filter with 1 dB band-pass ripple, and low and high band-pass frequencies are 1 and 50 Hz, respectively. So the bandpass filter transfer function is given by

(15)

For and and are small real numbers ε_{ij} = 0.1 the coefficients of the simple tuning function are given in the **Table 2**.

d1=b1 × 10^{04} |
c1 |
c2 × 10^{04} |
c3 × 10^{05} |
c4 × 10^{06} |

3.2624 | 200 | 3.897 | 2.369 | 1.403 |

**Table 1:** Bandpass filter parameters.

b1 × 10^{04} |
a1 |
a2 × 10^{04} |
a3 × 10^{05} |
a4 × 10^{06} |

3.2624 | 373.4367 | 8.6845 | 2.8156 | 1.4027 |

**Table 2:** Values of simple tuning function *G(S)* (b_{1}; a_{1}; a_{2}; a_{3} and a_{4}).

Under these conditions, a 50 s Matlab simulation could be performed in under 5 s. The identified frequencies are shown in **Figure 3**. We can see good identification and tracking of the fundamental frequencies. **Figures 4** and **5** show a close up of the actual outputs of the signal generators and the identified signals (very good tracking of amplitude and relative phases). The first component has significant DC which we have not attempted to identify. In particular, there is no way to distinguish and hence identify the DC content of the two true signals. Again we show good matches and thus we are able to identify these periodic signals in real time.

To get an overview of the signal frequency content and the accuracy of the identified models, the FFT transforms of the signal to be identified and the error signal are shown in **Figure 6**. It can be seen that most of harmonics have been identified although there is a huge DC component in both the signal and (*e*) in the proposed algorithm, which is as anticipated as we did not attempt to identify it.

**Figure 7** displays the quasi-periodic signals *d _{i}(t)*, the identified signal

It can be shown that the choice for controller parameters of *G(s)*=1, *K _{1ij}*=0 and

The authors would like to thank Libyan government and Libyan Ministry of higher education and scientific research for their fund and support for his research despite the hard environment that they work at due to unstable situation in Libya. In addition we would like to acknowledge the research funding provided by Western University.

This article has shown the instantaneous Fourier decomposition algorithm that is based on the orthogonal state variables of an internal model principle controller. First we examined how we implemented this algorithm off-line by matching the coefficients of nominators and denominators in both equations (11,12). Second, and the main contribution in this article is to develop a means of calculating the controller parameters that has a lower computational burden such that it can be successfully implemented in realtime. As a result of our work in section IV-B, the schema has been successfully implemented online after solving for the issue of overlapping harmonics from different signal components. One of the models is removed from the design process eliminating the dependent equations. The associated controller gains are distributed equally in the implemented controller. Thus, the instantaneous Fourier decomposition is implemented in real time, the frequency is identified with high speed of convergence and the predictable disturbance is identified as well as the system stability is guaranteed.

Our future work will be conducted to identify the uncertain frequencies of periodic signals and eliminate periodic disturbances in discrete state space form.

- Regalia PA (1991) An improved lattice-based adaptive IIR notch filter. IEEE Transactions on Signal Processing 39: 2124-2128.
- Kim H, Shim H, Jo NH (2014) Adaptive add-on output regulator for rejection of sinusoidal disturbances and application to optical disc drives. IEEE Transactions on Industrial Electronics 61: 5490-5499.
- Huang NE, Shen Z, Long SR, Wu MC, Shih HH (1998) The empirical mode decomposition and the hilbert spectrum for nonlinear and non-stationary time series analysis. Proceedings of the Royal Society A 454: 903-995.
- Hara S, Yamamoto Y, Omata T, Nakano M (1988) Repetitive control system: a new type servo system for periodic exogenous signals. IEEE Transactions on Automatic Control 33: 659-668.
- Bodson M, Sacks A, Khosla P (1992) Harmonic generation in adaptive feed forward cancellation schemes. Proceedings of the 31st IEEE Conference on Decision and Control 2: 1261-1266.
- Francist BA, Wonham WM (1976) The internal model principle of control theory. Elsevier 12: 457-465.
- Brown LJ, Zhang Q (2001) Control for canceling periodic disturbances with uncertain frequency. Proceedings of the 40th IEEE Conference on Decision and Control 5: 4909-4914.
- Kim W, Kim H, Chung CC, Tomizuka M (2011) Adaptive output regulation for the rejection of a periodic disturbance with an unknown frequency. IEEE Transactions on Control Systems Technology 19: 1296-1304.
- Nagashima M, Usui K, Kobayashi M (2007) Rejection of unknown periodic disturbances in magnetic hard disk drives. IEEE Transactions on Magnetics 43: 3774-3778.
- Zhang Q, Brown LJ (2003) Designing of adaptive bandpass filter with adjustable notch for frequency demodulation. Proceedings of the American Control Conference 4: 2931-2936.
- Brown LJ, Zhang Q (2003) Identification of periodic signals with uncertain frequency. IEEE Transactions on Signal Processing 51: 1538-1545.
- Mohsen E, Brown LJ (2017) Realtime implementation of an internal-model-principle signal identifier. IEEE 30th Canadian Conference on Electrical and Computer Engineering (CCECE).

Select your language of interest to view the total content in your interested language

- Advantages of E-Business
- Antenna Design
- Applied Statistics
- Applied Thermodynamics
- Biochemical Thermodynamics
- Business & Management
- Capacitor
- Catalyst
- Circuits
- Coal Consumption
- Cognitive Systems Engineering
- Commercialization of New Techniques
- Computation theory
- Computer Networking
- Conditioning Monitoring
- Contingency Analysis
- Controllers
- Cordless
- Current
- Cyber Physical System
- Cyber Security
- Datamining
- Design and Microfabrication
- Digital Signal Processing
- Directional Microphone
- Dynamic Stabilizer
- Dynamical System
- Electrical Engineering
- Electricity
- Electromagnetism
- Entrepreneurship
- Femtosecond Laser
- Fibre Optical Sensors
- Finance management
- Global optimization
- Health care management
- Image Recognition
- Industrial Crystallization
- Information Science
- Information Visualization
- Infrared
- Integrated Circuit
- Internet computing
- Laser Diode
- Laser Spectrum
- Logistics
- Machines
- Management Cybernetics
- Manufacturing system
- Materials Management
- Mathematical optimization
- Mechatronics and Robotics
- Microphone Array
- Modeling and Simulation
- Modern Thermodynamics
- Molecular Catalysis
- Molecular thermodynamics
- Multimedia Network
- Nano Photonics
- Nanomaterials in Catalysis
- Network
- Nonlinear Dynamics
- Nonlinear Optics
- Nucleophile Catalysis
- Omnidirectional Microphone
- Operations Research
- Optical Fiber
- Optical Switch
- Optics
- Organometallic Catalysis
- Oscillations
- Oscillator
- Oxidation Catalysis
- Parallel Processing
- Photon
- Photonic Band Gap
- Photonic Crystal
- Photonic Crystal Fiber
- Photonics
- Power Analysis
- Power System
- Power System Stabilizers
- Process Engineering
- Production and Operations Management
- Proline Catalysis
- Protein Folding Thermodynamics
- Protein Thermodynamics
- Quantum Dot
- Quantum Well
- Quaternary Logics
- Radar Imaging
- Radio Waves
- Refractive Index
- Reliability engineering
- Robot Audition
- Satellite Communication
- Semiconductor
- Semiconductor Laser
- Sensor
- Signal Crystal
- Signal Processing
- Simulator
- Single Mode
- Smart Grid
- Soft Computing
- Soft Computing and Computational Intelligent
- Software design
- Spatial Database
- Spectroscopy and Catalysis
- Spectrum
- Stochastic control
- Supply Chain Management
- Technologies Management
- Technology
- Telecommunications Journal
- Thermodynamics Databases and Analysis
- Thermodynamics Material Science
- Thermodynamics Physics
- Thermodynamics experiments updates
- Two Photonabsorption
- Ubiquitous Computing
- Wave Length
- Wave Theory
- Web Developer
- Wireless Journals
- Wireless Sensor

- Total views:
**207** - [From(publication date):

July-2017 - Oct 19, 2017] - Breakdown by view type
- HTML page views :
**172** - PDF downloads :
**35**

Peer Reviewed Journals

International Conferences 2017-18