Next Article in Journal
Synthesis and Computer Study of Population Dynamics Controlled Models Using Methods of Numerical Optimization, Stochastization and Machine Learning
Next Article in Special Issue
Dynamical Analysis of Discrete-Time Two-Predators One-Prey Lotka–Volterra Model
Previous Article in Journal
Enhancing Gain for UWB Antennas Using FSS: A Systematic Review
Previous Article in Special Issue
Analysis of Fractional-Order Nonlinear Dynamic Systems with General Analytic Kernels: Lyapunov Stability and Inequalities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hierarchical Quasi-Fractional Gradient Descent Method for Parameter Estimation of Nonlinear ARX Systems Using Key Term Separation Principle

by
Naveed Ishtiaq Chaudhary
1,
Muhammad Asif Zahoor Raja
1,
Zeshan Aslam Khan
2,
Khalid Mehmood Cheema
3,* and
Ahmad H. Milyani
4
1
Future Technology Research Center, National Yunlin University of Science and Technology, 123 University Road, Section 3, Douliou 64002, Taiwan
2
Department of Electrical Engineering, International Islamic University, Islamabad 44000, Pakistan
3
School of Electrical Engineering, Southeast University, Nanjing 210096, China
4
Department of Electrical and Computer Engineering, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(24), 3302; https://doi.org/10.3390/math9243302
Submission received: 29 November 2021 / Revised: 14 December 2021 / Accepted: 16 December 2021 / Published: 18 December 2021
(This article belongs to the Special Issue Fractional Calculus and Nonlinear Systems)

Abstract

:
Recently, a quasi-fractional order gradient descent (QFGD) algorithm was proposed and successfully applied to solve system identification problem. The QFGD suffers from the overparameterization problem and results in estimating the redundant parameters instead of identifying only the actual parameters of the system. This study develops a novel hierarchical QFDS (HQFGD) algorithm by introducing the concepts of hierarchical identification principle and key term separation idea. The proposed HQFGD is effectively applied to solve the parameter estimation problem of input nonlinear autoregressive with exogeneous noise (INARX) system. A detailed investigation about the performance of HQFGD is conducted under different disturbance conditions considering different fractional orders and learning rate variations. The simulation results validate the better performance of the HQFGD over the standard counterpart in terms of estimation accuracy, convergence speed and robustness.

1. Introduction

Fractional calculus has emerged as an important tool to effectively model a variety of systems or processes that arise in the various applications relating to physics, engineering and applied sciences [1,2,3,4,5,6,7,8]. Over the last decade, a trend of developing novel adaptive algorithms based on strong foundations of fractional calculus can be seen [9,10,11,12,13,14,15]. The idea of incorporating the fractional derivative into the conventional iterative procedure of the least mean square (LMS) algorithm was first introduced in 2009 [16]. The fractional derivative-based LMS, i.e., F-LMS, was then exploited to solve communication [17,18,19], vibration rejection [20], control [21], recommender systems [22,23], power signal modelling [24], signal processing [25] and time series prediction problems [26]. The variants of the F-LMS were proposed to solve different identification and parameter estimation problems [27,28,29,30,31]. Cheng et al. presented the concept of a variable initial value and proposed an innovative F-LMS (IF-LMS) [32], and Chaudhary et al. applied the IF-LMS for the parameter estimation of power signals [33]. Chen et al. compared the performance of different fractional order gradient methods [34]. Wei. et al. suggested a modification to handle a long memory property of fractional calculus and proposed a generalized fractional order gradient method (GFGM) [35]. Liu et al. generalized the GFGM to solve convex optimization problem with high dimensions, and proposed a quasi-fractional gradient (QFGD) method with theoretical analysis and convergence proof [36]. Furthermore, they also applied QFGD to solve system identification problem through overparameterization and estimated the redundant parameters instead of estimating only the actual parameters of the system. In order to avoid the problem of overparameterization, we propose hierarchical QFGD (HQFGD) by integrating the hierarchical identification principle and key term separation idea with QFGD for efficient parameter estimation. The hierarchical identification principle decomposes the system into multiple subsystems, thus reducing the overall dimensions of the system, and the key term separation technique allows one to avoid the overparameterization issue [37,38,39]. In this study, we propose HQFGD method for efficient parameter estimation of input nonlinear autoregressive exogenous noise (INARX) systems. The INARX system belongs to a class of block-oriented nonlinear systems that are widely used to model a variety of nonlinearities [40,41,42,43]. The main contributions of this study in terms of salient features are summarized as:
  • A novel hierarchical quasi-fractional gradient descent, HQFGD, algorithm is presented by integrating the hierarchical identification theory and key term separation technique with the quasi-fractional gradient method.
  • The hierarchical identification procedure decomposes the system into subsystems, thus reducing the overall dimensions of the system, and the key term separation technique allows one to avoid the overparameterization issue.
  • The accuracy and robustness of the proposed HQFGD is established through effective parameter estimation of input nonlinear autoregressive exogenous noise, INARX, systems under different disturbance conditions, fractional orders and learning rate variations.
  • The comparison with the standard counterpart validates the efficacy of the proposed HQFGD scheme in terms of convergence sped and estimation accuracy.
The remainder of this paper is structured as follows: Section 2 gives the description of the INARX system. Section 3 describes the HGD approach and Section 4 provides the design of the hierarchical QFGD for INARX systems. Section 5 discusses the results of numerical experimentation and Section 6 presents the conclusions with some future works.

2. Nonlinear ARX System Model

The governing mathematical relation of the INARX system with input x, output y, and exogeneous noise n is written as [44,45]:
o ( t ) = V ( z ) W ( z ) x ¯ ( t ) + 1 W ( z ) n ( t ) ,
where W ( z ) , V ( z ) and x ¯ ( t ) are defined as
W ( z ) = 1 + w 1 z 1 + w 2 z 2 + , , + w n w z n w ,
V ( z ) = r 0 + v 1 z 1 + v 2 z 2 + , , + v n v z n v ,
x ¯ ( t ) = u 1 γ 1 ( x ( t ) ) + u 2 γ 2 ( x ( t ) ) + + u p γ p ( x ( t ) ) .
Rewriting (1) after rearranging,
o ( t ) = ( 1 W ( z ) ) o ( t ) + V ( z ) x ¯ ( t ) + n ( t ) .
Putting Expressions (2)–(4) into (5), let v 0 = 1 , and considering x ¯ ( t ) as a key term, apply key term separation:
o ( t ) = k = 1 n w w k o ( t k ) + k = 0 n v v k x ¯ ( t k ) + n ( t ) = k = 1 n w w k o ( t k ) + v 0 x ¯ ( t ) + k = 1 n v v k x ¯ ( t k ) + n ( t ) = k = 1 n w w k o ( t k ) + k = 1 n v v k x ¯ ( t k ) + i = 1 p u i γ i ( x ( t ) ) + n ( t )
The information vector is defined as:
Ω w ( t ) = [ o ( t 1 ) , o ( t 2 ) , , o ( t n w ) ] T n w Ω v ( t ) = [ x ¯ ( t 1 ) , x ¯ ( t 2 ) , , x ¯ ( t n v ) ] T n v γ ( t ) = [ γ 1 ( u ( t ) ) , γ 2 ( u ( t ) ) , , γ p ( u ( t ) ) ] T p
and the corresponding parameter vectors are:
w = [ w 1 , w 2 , , w n w ] T n w v = [ v 1 , v 2 , , v n v ] T n v u = [ u 1 , u 2 , , u p ] T p
Putting (7) and (8) into (6), the output of the INARX system using key term separation [37] is written as:
o ( t ) = Ω w T ( t ) w + Ω v T ( t ) v + γ T ( t ) u + n ( t ) .

3. Hierarchical Gradient Descent Method

The brief overview of hierarchical identification based gradient descent (HGD) is presented in this section. The intermediate variables are defined as:
o w ( t ) = o ( t ) Ω v T ( t ) v γ T ( t ) u ,
o v ( t ) = o ( t ) Ω w T ( t ) w γ T ( t ) u ,
o u ( t ) = o ( t ) Ω w T ( t ) w Ω v T ( t ) v ,
Decomposing (9) into three sub-systems:
o w ( t ) = Ω w T ( t ) w + n ( t ) ,
o v ( t ) = Ω v T ( t ) v + n ( t ) ,
o u ( t ) = γ T ( t ) u + n ( t ) .
From (13)–(15), the cost functions are defined as:
J w ( w ) = ( o w ( t ) Ω w T ( t ) w ) 2 ,
J v ( v ) = ( o v ( t ) Ω v T ( t ) v ) 2 ,
J u ( u ) = ( o u ( t ) γ T ( t ) u ) 2 .
The following update rules are obtained by minimizing cost functions (16)–(18) through the gradient descent approach:
w ^ ( t ) = w ^ ( t 1 ) + α 1 Ω w ( t ) [ o w ( t ) Ω w T ( t ) w ^ ( t 1 ) ] = w ^ ( t 1 ) + α 1 Ω w ( t ) [ o ( t ) Ω v T ( t ) v ( t 1 ) γ T ( t ) u ^ ( t 1 ) Ω w T ( t ) w ^ ( t 1 ) ]
v ^ ( t ) = v ^ ( t 1 ) + α 2 Ω v ( t ) [ o r ( t ) Ω v T ( t ) v ^ ( t 1 ) ] = v ^ ( t 1 ) + α 2 Ω v ( t ) [ o ( t ) Ω w T ( t ) w ^ ( t 1 ) γ T ( t ) u ^ ( t 1 ) Ω v T ( t ) v ^ ( t 1 ) ]
u ^ ( t ) = u ^ ( t 1 ) + α 3 γ ( t ) [ o u ( t ) γ T ( t ) u ^ ( t 1 ) ] = u ^ ( t 1 ) + α 3 γ ( t ) [ o ( t ) Ω w T ( t ) w ^ ( t 1 ) Ω v T ( t ) v ^ ( t 1 ) γ T ( t ) u ^ ( t 1 ) ]
There are unmeasurable inner terms in Ω v ( t ) of (19)–(21), whose estimation is defined as:
Ω ^ v ( t ) = [ x ¯ ^ ( t 1 ) , x ¯ ^ ( t 2 ) , , x ¯ ^ ( t n v ) ] T ,
x ¯ ^ ( t ) = i = 1 p u ^ i ( t ) γ i ( x ( t ) ) = γ T ( t ) u ^ ( t ) .
Replacing Ω v ( t ) with Ω ^ v ( t ) in (19)–(21) gives the HGD for INARX identification:
w ^ ( t ) = w ^ ( t 1 ) + α 1 Ω w ( t ) [ o ( t ) Ω ^ v T ( t ) v ( t 1 ) γ T ( t ) u ^ ( t 1 ) Ω w T ( t ) w ^ ( t 1 ) ] ,
v ^ ( t ) = v ^ ( t 1 ) + α 2 Ω v ( t ) [ o ( t ) Ω w T ( t ) w ^ ( t 1 ) γ T ( t ) u ^ ( t 1 ) Ω ^ v T ( t ) v ^ ( t 1 ) ] ,
u ^ ( t ) = u ^ ( t 1 ) + α 3 γ ( t ) [ o ( t ) Ω w T ( t ) w ^ ( t 1 ) Ω ^ v T ( t ) v ^ ( t 1 ) γ T ( t ) u ^ ( t 1 ) ] .

4. Hierarchical Quasi Fractional Gradient Descent Method

In this section, the design of the proposed hierarchical quasi fractional gradient descent HQFGD is presented. The general iterative update rule of QFGD for a multi-variable function f ( s ) is written as reported in the recent literature [36]:
s ( t ) = s ( t 1 ) α Γ ( 2 μ ) f ( s ) | t = t 1 [ ( s ( t 1 ) s ( t 2 ) + ε ) 1 μ ] ,
where 0 < μ < 2 is the fractional order, α is the learning rate, ε is a vector of small numbers, is the conventional first order gradient, · denotes the 2 norms of a vector, and represents the Hadamard product. The search direction in QFGD is produced through the Hadamard product between the direction of the fractional gradient and another vector on the basis of the fractional gradient direction, as described in [36]. This helps in avoiding zigzag during gradient descent search by rotating the fractional gradient direction. Thus, f ( s ) | t = t 1 [ ( s ( t 1 ) s ( t 2 ) + ε ) 1 μ ] is called the QFGD direction. Further details about QFGD can be seen in [36].
Minimizing the cost functions (16)–(18) by using (27), the iterative update relations of the QFGD for w , v , u are written as:
w ^ ( t ) = w ^ ( t 1 ) + α 1 Ω w ( t ) Γ ( 2 μ ) [ o w ( t ) Ω w T ( t ) w ^ ( t 1 ) ] [ ( w ^ ( t 1 ) w ^ ( t 2 ) + ε ) 1 μ ] ,
v ^ ( t ) = v ^ ( t 1 ) + α 2 Ω v ( t ) Γ ( 2 μ ) [ o v ( t ) Ω v T ( t ) v ^ ( t 1 ) ] [ ( v ^ ( t 1 ) v ^ ( t 2 ) + ε ) 1 μ ] ,
u ^ ( t ) = u ^ ( t 1 ) + α 3 γ ( t ) Γ ( 2 μ ) [ o u ( t ) γ ( t ) u ^ ( t 1 ) ] [ ( u ^ ( t 1 ) u ^ ( t 2 ) + ε ) 1 μ ] .
Using (10)–(12) in (28)–(30), respectively:
w ^ ( t ) = w ^ ( t 1 ) + α 1 Ω w ( t ) e w ( t ) Γ ( 2 μ ) [ ( w ^ ( t 1 ) w ^ ( t 2 ) + ε ) 1 μ ] ,
v ^ ( t ) = v ^ ( t 1 ) + α 2 Ω v ( t ) e v ( t ) Γ ( 2 μ ) [ ( v ^ ( t 1 ) v ^ ( t 2 ) + ε ) 1 μ ] ,
u ^ ( t ) = u ^ ( t 1 ) + α 3 γ ( t ) e u ( t ) Γ ( 2 μ ) [ ( u ^ ( t 1 ) u ^ ( t 2 ) + ε ) 1 μ ] ,
where e w , e v and e u are defined as:
e w ( t ) = [ o ( t ) γ T ( t ) u ^ ( t 1 ) Ω v T ( t ) v ^ ( t 1 ) Ω w T ( t ) w ^ ( t 1 ) ] ,
e v ( t ) = [ o ( t ) Ω w T ( t ) w ^ ( t 1 ) γ T ( t ) u ^ ( t 1 ) Ω v T ( t ) v ^ ( t 1 ) ] ,
e u ( t ) = [ o ( t ) Ω w T ( t ) w ^ ( t 1 ) Ω v T ( t ) v ^ ( t 1 ) γ T ( t ) u ^ ( t 1 ) ] .
Using (22)–(23) to replace Ω v ( t ) with estimate Ω ^ v ( t ) in (31)–(33) gives the hierarchical QFGD method for parameter estimation of INARX systems:
w ^ ( t ) = w ^ ( t 1 ) + α 1 Ω w ( t ) e w ( t ) Γ ( 2 μ ) [ ( w ^ ( t 1 ) w ^ ( t 2 ) + ε ) 1 μ ] e w ( t ) = [ o ( t ) γ T ( t ) u ^ ( t 1 ) Ω ^ v T ( t ) v ^ ( t 1 ) Ω w T ( t ) w ^ ( t 1 ) ]
v ^ ( t ) = v ^ ( t 1 ) + α 2 Ω v ( t ) e v ( t ) Γ ( 2 μ ) [ ( v ^ ( t 1 ) v ^ ( t 2 ) + ε ) 1 μ ] e v ( t ) = [ o ( t ) Ω w T ( t ) w ^ ( t 1 ) γ T ( t ) u ^ ( t 1 ) Ω ^ v T ( t ) v ^ ( t 1 ) ]
u ^ ( t ) = u ^ ( t 1 ) + α 3 γ ( t ) e u ( t ) Γ ( 2 μ ) [ ( u ^ ( t 1 ) u ^ ( t 2 ) + ε ) 1 μ ] e u ( t ) = [ o ( t ) Ω w T ( t ) w ^ ( t 1 ) Ω ^ v T ( t ) v ^ ( t 1 ) γ T ( t ) u ^ ( t 1 ) ]
Putting μ = 1 into (34)–(36) reduces the HQFGD to conventional HGD (24)–(26). The overall graphical flow of the study is presented in Figure 1.

5. Results and Discussion

In this section, results of the proposed HQFGD are presented for parameter estimation of INARX system, and detail analyses is conducted regarding estimation accuracy, convergence speed and robustness of the scheme. Considering the following INARX system for numerical experimentation [37]:
W ( z ) = 1 + w 1 z 1 + w 2 z 2 = 1 + 1.35 z 1 + 0.75 z 2 ,
V ( z ) = v 1 z 1 + v 2 z 2 = 0.25 z 1 + z 2 ,
x ¯ ( t ) = u 1 γ 1 ( x ( t ) ) + u 2 γ 2 ( x ( t ) ) = 0.5 x ( t ) + 0.9 x 2 ( t ) .
The parameters to be estimated are:
Φ = [ w , v , u ] T = [ w 1 , w 2 v 1 , v 2 u 1 , u 2 ] T = [ 1.35 , 0.75 , 0.25 , 1 , 0.5 , 0.9 ] T
The following evaluation metrics are developed to measure the estimation performance of the HQFGD for INARX
Fitness = Φ Φ ^ Φ ,
MSE = mean ( Φ Φ ^ ) 2 ,
where · denotes the 2 norms of a vector. The simulations are conducted in Matlab, where the input signal with length is randomly generated with zero mean and unit variance, while the disturbance/noise signal is generated with normal distribution and constant variance. The considered input-output data length is 30,000. The performance of the HQFGD is assessed for two learning rates ( α = 0.0005, 0.001), five fractional orders ( μ = 0.7, 0.8, 0.9, 1.0, 1.1, 1.2) and three disturbance levels ( σ = 0.02, 0.09, 0.2).
The fitness iterative plots of the HQFGD based on the learning rate parameter are presented in Figure 2 for μ = 0.8, 1.0, 1.2 and σ = 0.02, 0.2. It is seen that the convergence speed of the HQFGD is faster for higher value of α but at the cost of more steady state mis adjustments, while the smaller value of α is good in steady state dynamics with smooth convergence curve but suffers from slow convergence speed. Thus, the small value of the learning rate ( α = 0.0005) is selected, and the tabular presentation of the iterative fitness values is given in Supplementary Tables S1–S3 for σ = 0.02, 0.09 and 0.2, respectively.
The fitness iterative plots of the HQFGD based on the fractional order variations are presented in Figure 3 for σ = 0.02, 0.09 and 0.2. The results indicate that the HQFGD is convergent for all considered fractional orders with relatively better convergence speed in the cases of μ = 0.8 and 0.9. The presented results further indicate that the HQFGD is robust against disturbances, with little decrease in the accuracy level for high noise variance.
The parameter estimation accuracy of the HQFGD is presented through parameter iterative plots and is given in Figure 4 for μ = 0.8, 1, 1.2, σ = 0.02 and α = 0.0005. The results show that the HQFGD accurately estimated the parameters of INARX system for all fractional order variations with relatively better convergence speed in the case of μ = 0.8.
The performance of the HQFGD for parameter estimation of INARX system is further assessed through MSE-based evaluation metrics, and the results are presented in Table 1 and Table 2 for α = 0.0005 and 0.001, respectively, in the case of all fractional orders and disturbance levels. These MSE results further verify the fitness-based results, showing that the HQFGD is accurate and convergent for all fractional orders. The HQFGD is relatively better regarding convergence speed than the conventional HGD for μ = 0.8 and 0.9, while the HQFGD is relatively better than the HGD with regard to final accuracy for μ = 1.1 and 1.2.

6. Conclusions

The following are the main conclusions:
  • A novel design of hierarchical quasi-fractional gradient descent, HQFGD, is presented for effective parameter estimation of input nonlinear autoregressive systems with exogeneous disturbance, i.e., INARX systems.
  • The HQFGD is developed by incorporating the hierarchical identification principle and key term separation idea into the structure of QFGD. The hierarchical identification procedure decomposes the INARX system into different subsystems and reduces the computational complexity of the conventional counterpart.
  • The HQFGD effectively estimates the parameters of the INARX system by considering only the actual system parameters and avoiding estimation of the redundant parameters due to overparameterization problems caused by the product of the cross terms.
  • The HQFGD is accurate, robust, and convergent in comparison with the standard counterpart for parameter estimation of INARX systems.
  • The HQFGD is relatively better with regard to convergence speed than the conventional HGD for μ = 0.8 and 0.9, while the HQFGD is relatively better than the HGD as regards final accuracy for μ = 1.1 and 1.2.
In future, one may exploit the proposed methodology for solving complex estimation problems [46,47] with hysteresis/dead zone nonlinearities [48,49] and investigate the recently introduced fractional derivative definitions for designing novel fractional gradient methods [50,51,52] to solve complex optimization problems [53,54,55].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/math9243302/s1. Results of fitness adaptation of the HQFGD in tabular form are given as a supplementary material in Supplementary Tables S1–S3 for 0.02, 0.09 and 0.2 noise levels, respectively.

Author Contributions

Conceptualization, N.I.C. and Z.A.K.; methodology, N.I.C. and M.A.Z.R.; software, N.I.C.; validation, M.A.Z.R. and Z.A.K.; resources, K.M.C. and A.H.M.; writing—original draft preparation, N.I.C.; writing—review and editing, Z.A.K. and M.A.Z.R.; project administration, K.M.C. and A.H.M.; funding acquisition, K.M.C. and A.H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge the kind help and great support of the late J.A. Tenreiro Machado in developing the new paradigm of fractional order adaptive algorithms by exploring the applications of fractional calculus in adaptive methods. He motivated us to explore the applications of novel fractional order algorithms to solve different parameter estimation and nonlinear identification problems.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sabatier, J.; Agrawal, O.P.; Machado, J.T. Advances in Fractional Calculus, 1st ed.; Springer: Dordrecht, The Netherlands, 2007; Volume 4. [Google Scholar]
  2. Machado, J.T.; Kiryakova, V.; Mainardi, F. Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1140–1153. [Google Scholar] [CrossRef] [Green Version]
  3. Tejado, I.; Pérez, E.; Valério, D. Fractional calculus in economic growth modelling of the group of seven. Fract. Calc. Appl. Anal. 2019, 22, 139–157. [Google Scholar] [CrossRef]
  4. Rashid, S.; Hammouch, Z.; Aydi, H.; Ahmad, A.G.; Alsharif, A.M. Novel computations of the time-fractional Fisher’s model via generalized fractional integral operators by means of the Elzaki transform. Fractal Fract. 2021, 5, 94. [Google Scholar] [CrossRef]
  5. Sun, H.; Zhang, Y.; Baleanu, D.; Chen, W.; Chen, Y. A new collection of real world applications of fractional calculus in science and engineering. Commun. Nonlinear Sci. Numer. Simul. 2018, 64, 213–231. [Google Scholar] [CrossRef]
  6. Masood, Z. Fractional Dynamics of Stuxnet Virus Propagation in Industrial Control Systems. Mathematics 2021, 9, 2160. [Google Scholar] [CrossRef]
  7. Valério, D.; Ortigueira, M.D.; Tenreiro, J.M.; Lopes, A.M. Continuous-time fractional linear systems: Steady-state responses. In Volume 6 Applications in Control; De Gruyter: Berlin, Germany, 2019; pp. 149–174. [Google Scholar]
  8. Christ, L.F.; Valério, D.; Coelho, R.M.; Vinga, S. Models of bone metastases and therapy using fractional derivatives. J. Appl. Nonlinear Dyn. 2018, 7, 81–94. [Google Scholar] [CrossRef]
  9. Pires, E.S.; Machado, J.T.; de Moura Oliveira, P.B.; Cunha, J.B.; Mendes, L. Particle swarm optimization with fractional-order velocity. Nonlinear Dyn. 2010, 61, 295–301. [Google Scholar] [CrossRef] [Green Version]
  10. Mousavi, Y.; Alfi, A. Fractional calculus-based firefly algorithm applied to parameter estimation of chaotic systems. Chaos Solitons Fractals 2018, 114, 202–215. [Google Scholar] [CrossRef]
  11. Ray, P.K.; Paital, S.R.; Mohanty, A.; Foo, Y.E.; Krishnan, A.; Gooi, H.B.; Amaratunga, G.A. A hybrid firefly-swarm optimized fractional order interval type-2 fuzzy PID-PSS for transient stability improvement. IEEE Trans. Ind. Appl. 2019, 55, 6486–6498. [Google Scholar] [CrossRef]
  12. Muhammad, Y.; Khan, R.; Raja, M.A.Z.; Ullah, F.; Chaudhary, N.I.; He, Y. Design of fractional swarm intelligent computing with entropy evolution for optimal power flow problems. IEEE Access 2020, 8, 111401–111419. [Google Scholar] [CrossRef]
  13. Yousri, D.; Mirjalili, S. Fractional-order cuckoo search algorithm for parameter identification of the fractional-order chaotic, chaotic with noise and hyper-chaotic financial systems. Eng. Appl. Artif. Intell. 2020, 92, 103662. [Google Scholar] [CrossRef]
  14. Kiani-B, A.; Fallahi, K.; Pariz, N.; Leung, H. A chaotic secure communication scheme using fractional chaotic systems based on an extended fractional Kalman filter. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 863–879. [Google Scholar] [CrossRef]
  15. Wang, X.; Zhang, F.; Kung, H.T.; Johnson, V.C.; Latif, A. Extracting soil salinization information with a fractional-order filtering algorithm and grid-search support vector machine (GS-SVM) model. Int. J. Remote Sens. 2020, 41, 953–973. [Google Scholar] [CrossRef]
  16. Zahoor, R.M.A.; Qureshi, I.M. A modified least mean square algorithm using fractional derivative and its application to system identification. Eur. J. Sci. Res. 2009, 35, 14–21. [Google Scholar]
  17. Khan, A.A.; Shah, S.M.; Raja, M.A.Z.; Chaudhary, N.I.; He, Y.; Machado, J.A.T. Fractional LMS and NLMS Algorithms for Line Echo Cancellation. Arab. J. Sci. Eng. 2021, 46, 9385–9398. [Google Scholar] [CrossRef]
  18. Raja, M.A.Z.; Akhtar, R.; Chaudhary, N.I.; Zhiyu, Z.; Khan, Q.; Rehman, A.U.; Zaman, F. A new computing paradigm for the optimization of parameters in adaptive beamforming using fractional processing. Eur. Phys. J. Plus 2019, 134, 275. [Google Scholar] [CrossRef]
  19. Shah, S.M.; Samar, R.; Raja, M.A.Z.; Dynamics, N. Fractional-order algorithms for tracking Rayleigh fading channels. Nonlinear Dyn. 2018, 92, 1243–1259. [Google Scholar] [CrossRef]
  20. Yin, W.; Wei, Y.; Liu, T.; Wang, Y. A novel orthogonalized fractional order filtered-x normalized least mean squares algorithm for feedforward vibration rejection. Mech. Syst. Signal Process. 2019, 119, 138–154. [Google Scholar] [CrossRef]
  21. Chaudhary, N.I.; Ahmed, M.; Khan, Z.A.; Zubair, S.; Raja, M.A.Z.; Dedovic, N. Design of normalized fractional adaptive algorithms for parameter estimation of control autoregressive autoregressive systems. Appl. Math. Model. 2018, 55, 698–715. [Google Scholar] [CrossRef]
  22. Khan, Z.A.; Chaudhary, N.I.; Zubair, S. Fractional stochastic gradient descent for recommender systems. Electron. Mark. 2019, 29, 275–285. [Google Scholar] [CrossRef]
  23. Khan, Z.A.; Zubair, S.; Chaudhary, N.I.; Raja, M.A.Z.; Khan, F.A.; Dedovic, N. Design of normalized fractional SGD computing paradigm for recommender systems. Neural Comput. Appl. 2020, 32, 10245–10262. [Google Scholar] [CrossRef]
  24. Chaudhary, N.I.; Aslam, M.S.; Baleanu, D.; Raja, M.A.Z. Design of sign fractional optimization paradigms for parameter estimation of nonlinear Hammerstein systems. Neural Comput. Appl. 2020, 32, 8381–8399. [Google Scholar] [CrossRef]
  25. Chaudhary, N.I.; Zubair, S.; Raja, M.A.Z. A new computing approach for power signal modeling using fractional adaptive algorithms. ISA Trans. 2017, 68, 189–202. [Google Scholar] [CrossRef]
  26. Shoaib, B.; Qureshi, I.M. Adaptive step-size modified fractional least mean square algorithm for chaotic time series prediction. Chin. Phys. B 2014, 23, 050503. [Google Scholar] [CrossRef]
  27. Cheng, S.; Wei, Y.; Sheng, D.; Chen, Y.; Wang, Y. Identification for Hammerstein nonlinear ARMAX systems based on multi-innovation fractional order stochastic gradient. Signal Process. 2018, 142, 1–10. [Google Scholar] [CrossRef]
  28. Chaudhary, N.I.; Raja MA, Z.; He, Y.; Khan, Z.A.; Machado, J.T. Design of multi innovation fractional LMS algorithm for parameter estimation of input nonlinear control autoregressive systems. Appl. Math. Model. 2021, 93, 412–425. [Google Scholar] [CrossRef]
  29. Aslam, M.S.; Chaudhary, N.I.; Raja, M.A.Z. A sliding-window approximation-based fractional adaptive strategy for Hammerstein nonlinear ARMAX systems. Nonlinear Dyn. 2017, 87, 519–533. [Google Scholar] [CrossRef]
  30. Zubair, S.; Chaudhary, N.I.; Khan, Z.A.; Wang, W. Momentum fractional LMS for power signal parameter estimation. Signal Process. 2018, 142, 441–449. [Google Scholar] [CrossRef]
  31. Chaudhary, N.I.; Zubair, S.; Aslam, M.S.; Raja, M.A.Z.; Machado, J.T. Design of momentum fractional LMS for Hammerstein nonlinear system identification with application to electrically stimulated muscle model. Eur. Phys. J. Plus 2019, 134, 407. [Google Scholar] [CrossRef]
  32. Cheng, S.; Wei, Y.; Chen, Y.; Li, Y.; Wang, Y. An innovative fractional order LMS based on variable initial value and gradient order. Signal Process. 2017, 133, 260–269. [Google Scholar] [CrossRef]
  33. Chaudhary, N.I.; Latif, R.; Raja, M.A.Z.; Machado, J.T. An innovative fractional order LMS algorithm for power signal parameter estimation. Appl. Math. Model. 2020, 83, 703–718. [Google Scholar] [CrossRef]
  34. Chen, Y.; Gao, Q.; Wei, Y.; Wang, Y. Study on fractional order gradient methods. Appl. Math. Comput. 2017, 314, 310–321. [Google Scholar] [CrossRef]
  35. Wei, Y.; Kang, Y.; Yin, W.; Wang, Y. Generalization of the gradient method with fractional order gradient direction. J. Frankl. Inst. 2020, 357, 2514–2532. [Google Scholar] [CrossRef]
  36. Liu, J.; Zhai, R.; Liu, Y.; Li, W.; Wang, B.; Huang, L. A quasi fractional order gradient descent method with adaptive stepsize and its application in system identification. Appl. Math. Comput. 2021, 393, 125797. [Google Scholar] [CrossRef]
  37. Chen, H.; Xiao, Y.; Ding, F. Hierarchical gradient parameter estimation algorithm for Hammerstein nonlinear systems using the key term separation principle. Appl. Math. Comput. 2014, 247, 1202–1210. [Google Scholar] [CrossRef]
  38. Ding, F.; Chen, H.; Xu, L.; Dai, J.; Li, Q.; Hayat, T. A hierarchical least squares identification algorithm for Hammerstein nonlinear systems using the key term separation. J. Frankl. Inst. 2018, 355, 3737–3752. [Google Scholar] [CrossRef]
  39. Ding, F.; Ma, H.; Pan, J.; Yang, E. Hierarchical gradient-and least squares-based iterative algorithms for input nonlinear output-error systems using the key term separation. J. Frankl. Inst. 2021, 358, 5113–5135. [Google Scholar] [CrossRef]
  40. Giri, F.; Bai, E.W. Block-Oriented Nonlinear System Identification, 1st ed.; Springer: London, UK, 2010; Volume 1. [Google Scholar]
  41. Billings, S.A. Nonlinear System Identification: NARMAX Methods in the Time, Frequency, and Spatio-Temporal Domains; John Wiley & Sons: Hoboken, NJ, USA, 2013; Volume 1. [Google Scholar]
  42. Schoukens, J.; Ljung, L. Nonlinear system identification: A user-oriented road map. IEEE Control Syst. Mag. 2019, 39, 28–99. [Google Scholar]
  43. Le, F.; Markovsky, I.; Freeman, C.T.; Rogers, E. Recursive identification of Hammerstein systems with application to electrically stimulated muscle. Control Eng. Pract. 2012, 20, 386–396. [Google Scholar] [CrossRef] [Green Version]
  44. Mehmood, A.; Zameer, A.; Chaudhary, N.I.; Raja, M.A.Z. Backtracking search heuristics for identification of electrical muscle stimulation models using Hammerstein structure. Appl. Soft Comput. 2019, 84, 105705. [Google Scholar] [CrossRef]
  45. Mehmood, A.; Zameer, A.; Chaudhary, N.I.; Ling, S.H.; Raja, M.A.Z. Design of meta-heuristic computing paradigms for Hammerstein identification systems in electrically stimulated muscle models. Neural Comput. Appl. 2020, 32, 12469–12497. [Google Scholar] [CrossRef]
  46. Mehmood, A.; Chaudhary, N.I.; Zameer, A.; Raja, M.A.Z. Novel computing paradigms for parameter estimation in power signal models. Neural Comput. Appl. 2020, 32, 6253–6282. [Google Scholar] [CrossRef]
  47. Mehmood, A.; Shi, P.; Raja, M.A.Z.; Zameer, A.; Chaudhary, N.I. Design of backtracking search heuristics for parameter estimation of power signals. Neural Comput. Appl. 2020, 33, 1479–1496. [Google Scholar] [CrossRef]
  48. Prasad, V.; Mehta, U. Modeling and parametric identification of Hammerstein systems with time delay and asymmetric dead-zones using fractional differential equations. Mech. Syst. Signal Process. 2022, 167, 108568. [Google Scholar] [CrossRef]
  49. Prasad, V.; Kothari, K.; Mehta, U. Parametric identification of nonlinear fractional Hammerstein models. Fractal Fract. 2020, 4, 2. [Google Scholar] [CrossRef] [Green Version]
  50. Atangana, A.; Shafiq, A. Differential and integral operators with constant fractional order and variable fractional dimension. Chaos Solitons Fractals 2019, 127, 226–243. [Google Scholar] [CrossRef]
  51. Ghanbari, B.; Atangana, A. A new application of fractional Atangana–Baleanu derivatives: Designing ABC-fractional masks in image processing. Phys. A Stat. Mech. Its Appl. 2020, 542, 123516. [Google Scholar] [CrossRef]
  52. Atangana, A. Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination? Chaos Solitons Fractals 2020, 136, 109860. [Google Scholar] [CrossRef]
  53. Shi, L.; Wang, X.; Hou, H. Research on Optimization of Array Honeypot Defense Strategies Based on Evolutionary Game Theory. Mathematics 2021, 9, 805. [Google Scholar] [CrossRef]
  54. Posypkin, M.; Khamisov, O. Automatic Convexity Deduction for Efficient Function’s Range Bounding. Mathematics 2021, 9, 134. [Google Scholar] [CrossRef]
  55. Lera, D.; Posypkin, M.; Sergeyev, Y.D. Space-filling curves for numerical approximation and visualization of solutions to systems of nonlinear inequalities with applications in robotics. Appl. Math. Comput. 2021, 390, 125660. [Google Scholar] [CrossRef]
Figure 1. Graphical abstract of the study.
Figure 1. Graphical abstract of the study.
Mathematics 09 03302 g001
Figure 2. Learning curves of the HQFGD for different step size variations. (a) μ = 0.8, σ = 0.02; (b) μ = 1, σ = 0.02; (c) μ = 1.2, σ = 0.02; (d) μ = 0.8, σ = 0.2; (e) μ = 1, σ = 0.2; (f) μ = 1.2, σ = 0.2.
Figure 2. Learning curves of the HQFGD for different step size variations. (a) μ = 0.8, σ = 0.02; (b) μ = 1, σ = 0.02; (c) μ = 1.2, σ = 0.02; (d) μ = 0.8, σ = 0.2; (e) μ = 1, σ = 0.2; (f) μ = 1.2, σ = 0.2.
Mathematics 09 03302 g002aMathematics 09 03302 g002b
Figure 3. Learning curves of the HQFGD for different fractional orders. (a) σ = 0.02; (b) σ = 0.09; (c) σ = 0.2.
Figure 3. Learning curves of the HQFGD for different fractional orders. (a) σ = 0.02; (b) σ = 0.09; (c) σ = 0.2.
Mathematics 09 03302 g003aMathematics 09 03302 g003b
Figure 4. Parameters learning curves of the HQFGD for noise level 0.02. (a) Estimate of w 1 ; (b) Estimate of w 2 ; (c) Estimate of v 1 ; (d) Estimate of v 2 ; (e) Estimate of u 1 ; (f) Estimate of u 2 .
Figure 4. Parameters learning curves of the HQFGD for noise level 0.02. (a) Estimate of w 1 ; (b) Estimate of w 2 ; (c) Estimate of v 1 ; (d) Estimate of v 2 ; (e) Estimate of u 1 ; (f) Estimate of u 2 .
Mathematics 09 03302 g004aMathematics 09 03302 g004b
Table 1. MSE results of HQFGD for different noise levels for α = 0.0005.
Table 1. MSE results of HQFGD for different noise levels for α = 0.0005.
Noise μ w 1 w 2 v 1 v 2 u 1 u 2 MSE
0.020.71.34990.75030.25041.00040.49980.90049.71 × 10−8
0.81.34980.75030.25041.00040.49980.90049.01 × 10−8
0.91.34980.75030.25031.00040.49990.90048.01 × 10−8
1.01.34980.75020.25021.00040.49990.90046.92 × 10−8
1.11.34970.75020.25001.00040.49990.90036.75 × 10−8
1.21.34960.75000.24961.00060.49990.90031.24 × 10−7
0.090.71.34960.75110.25161.00120.49940.90141.28 × 10−6
0.81.34960.75110.25161.00120.49940.90131.21 × 10−6
0.91.34960.75100.25151.00110.49950.90131.10 × 10−6
1.01.34960.75100.25131.00110.49950.90129.53 × 10−7
1.11.34950.75080.25101.00110.49950.90127.77 × 10−7
1.21.34950.75070.25051.00110.49960.90116.10 × 10−7
0.20.71.34860.75310.25461.00360.49830.90381.05 × 10−5
0.81.34870.75300.25441.00350.49840.90379.87 × 10−6
0.91.34870.75280.25421.00330.49840.90369.07 × 10−6
1.01.34880.75270.25391.00310.49850.90358.05 × 10−6
1.11.34880.75250.25351.00290.49860.90336.81 × 10−6
1.21.34890.75230.25281.00280.49880.90305.00 × 10−6
Φ 1.35000.75000.25001.00000.50000.90000
Table 2. MSE results of HQFGD for different noise levels for α = 0.001.
Table 2. MSE results of HQFGD for different noise levels for α = 0.001.
Noise μ w 1 w 2 v 1 v 2 u 1 u 2 MSE
0.020.71.34980.75070.25091.00100.49970.90064.46 × 10−7
0.81.34980.75070.25091.00100.49970.90064.27 × 10−7
0.91.34980.75060.25091.00090.49970.90054.00 × 10−7
1.01.34980.75060.25081.00090.49980.90053.64 × 10−7
1.11.34980.75060.25081.00080.49980.90053.20 × 10−7
1.21.34980.75050.25071.00070.49980.90052.70 × 10−7
0.090.71.34940.75240.25301.00360.49910.90205.41 × 10−6
0.81.34940.75230.25301.00350.49910.90195.19 × 10−6
0.91.34940.75220.25291.00330.49910.90194.85 × 10−6
1.01.34940.75210.25281.00310.49910.90194.41 × 10−6
1.11.34940.75190.25271.00280.49910.90183.88 × 10−6
1.21.34940.75180.25261.00240.49920.90173.27 × 10−6
0.20.71.34820.75670.25801.01060.49730.90554.36 × 10−5
0.81.34810.75620.25801.01030.49740.90544.14 × 10−5
0.91.34800.75590.25781.00980.49740.90543.86 × 10−5
1.01.34800.75560.25761.00910.49740.90523.51 × 10−5
1.11.34800.75520.25731.00830.49750.90513.08 × 10−5
1.21.34810.75470.25691.00730.49760.90492.59 × 10−5
Φ 1.35000.75000.25001.00000.50000.90000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chaudhary, N.I.; Raja, M.A.Z.; Khan, Z.A.; Cheema, K.M.; Milyani, A.H. Hierarchical Quasi-Fractional Gradient Descent Method for Parameter Estimation of Nonlinear ARX Systems Using Key Term Separation Principle. Mathematics 2021, 9, 3302. https://doi.org/10.3390/math9243302

AMA Style

Chaudhary NI, Raja MAZ, Khan ZA, Cheema KM, Milyani AH. Hierarchical Quasi-Fractional Gradient Descent Method for Parameter Estimation of Nonlinear ARX Systems Using Key Term Separation Principle. Mathematics. 2021; 9(24):3302. https://doi.org/10.3390/math9243302

Chicago/Turabian Style

Chaudhary, Naveed Ishtiaq, Muhammad Asif Zahoor Raja, Zeshan Aslam Khan, Khalid Mehmood Cheema, and Ahmad H. Milyani. 2021. "Hierarchical Quasi-Fractional Gradient Descent Method for Parameter Estimation of Nonlinear ARX Systems Using Key Term Separation Principle" Mathematics 9, no. 24: 3302. https://doi.org/10.3390/math9243302

APA Style

Chaudhary, N. I., Raja, M. A. Z., Khan, Z. A., Cheema, K. M., & Milyani, A. H. (2021). Hierarchical Quasi-Fractional Gradient Descent Method for Parameter Estimation of Nonlinear ARX Systems Using Key Term Separation Principle. Mathematics, 9(24), 3302. https://doi.org/10.3390/math9243302

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop