Moving+orizon)ault(stimation of a+ybrid6ystem using the6witched
ARX-Laguerre0odel
IBTISSEM BENGHARAT, ABDELKADER MBAREK
Research Laboratory of Automatic Signal and Image Processing,
National School of Engineers of Monastir,
University of Monastir, 5019,
TUNISIA
n
Abstract: Using moving horizon Fault estimation, we propose in this article a fault estimation of a linear hybrid
system using a reduced complexity SARX-Laguerre model. The linear hybrid systems are approximated by a
SARX-Laguerre model (Switched ARX-Laguerre).This latter is obtained by expanding a discrete time SARX
model parameters on Laguerre orthonormal bases. The resulting model ensures an efficient complexity reduction
with respect to the classical SARX model. This parametric complexity reduction is still subject to an optimal
choice of the Laguerre poles defining Laguerre bases. Therefore, we propose in this paper to identify the pa-
rameters of the SARX-Laguerre model by using a recursive algorithm to identify the Fourier coefficients and a
metaheuristic algorithm to identify the poles. The proposed model is built from the system input / output obser-
vations and is used to develop a fault estimation scheme. The scheme of the fault estimation based on the Moving
Horizon fault Estimation (MHE). The proposed fault estimation using moving horizon procedure is tested on
numerical simulation.
Key-Words: - Hybrid system, Switching models, Orthogonal bases, SARX-Laguerre model, Moving Horizon,
Fault Estimation.
Received: April 15, 2024. Revised: August 19, 2024. Accepted: October 14, 2024. Published: November 28, 2024.
1 Introduction
Modern technological systems always present an in-
crease on complexity and become more and more sus-
ceptible to faults. Therefore, in the last two decades
considerable research has been focused on model
based fault estimation, [1], [2], [3], [4], [5], [6], [7],
[8], [9]. The most important work in fault estimation
uses the moving horizon faults estimation (MHE) to
optimize faults. The MHE method is based on an es-
timate over a horizon using the past input and output
information and with this data, it is able to estimate
the current and past states based on a model of the
system. This estimator is formulated to minimize the
quadratic estimation error between the system output
and the model output over a horizon, [10].
The goal of this work is to synthesize a moving
horizon fault estimation (MHE) for linear hybrid sys-
tems using a reduced complexity switching model. A
hybrid system is a dynamic system that explicitly and
simultaneously involves continuous and discrete be-
haviors. These systems are classically consisting of
continuous processes interacting with or supervised
by discrete processes. Switched systems form an im-
portant class of hybrid systems, [11], [12], [13], [14],i
which consist of a finite number of subsystems and a
switching rule indicating the active subsystem at each
instant, [15]. The MHE is used to estimate the fault
actuator for a class of linear hybrid system known
as switching system by solving an optimization prob-
lem over a past time horizon under certain constraints.
The reduced complexity model used in the fault esti-
mation is the SARX-Laguerre obtained after the ex-
pansion of a Switched ARX model (SARX) on La-
guerre orthonormal bases by filtering the input and
the output of every submodel on two independent La-
guerre bases, [16], [17].
In last recent years, the use of orthogonal basis
in modelling both of linear systems, [18], [19], and
nonlinear systems, [16], [20], [21], has experienced
a very great development due to its capacity in the
reduction of parametric complexity. This modelling
is applied by filtering the input of the linear model
as FIR model, [22], or ARX model, and nonlinear
model as Volterra model by decomposing each ker-
nel on Laguerre basis, [23]. This input filtering guar-
antees a parametric reduction which can be signifi-
cant when the considered system is linear with a dom-
inant first order dynamic. However, in the case of
scattered poles and an oscillating system the Laguerre
model requires a huge number of Laguerre functions
and then many numbers of parameters to represent
systems with various representative modes. To over-
come this drawback, the study, [24], propose an al-
ternative solution to represent any complex dynamic
linear systems with a reduced parameter complexity
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
109
Volume 15, 2024
model using the Laguerre bases. This idea is based
on the expansion of the linear model on two indepen-
dent Laguerre bases by filtering the input and the out-
put of the linear model as ARX model. Due to its
high efficiency in reducing the parameter complexity,
this approach was extended to the case of linear mul-
tiple model approach were a ARX-Laguerre multiple
model approach is proposed by [16], and to the case
of Multiple Input Multiple Output (MIMO) systems,
where a ARX-Laguerre MIMO model is proposed by
[25]. Thus, to reduce the number of parameters of the
SARX-Laguerre model an optimal choice of the poles
is necessary. In this paper we propose to use a meta-
heuristic algorithm to optimize the SARX-Laguerre
poles and using recursive procedure to identify the
Fourier coefficients.
This paper is organized as follows: in section 2,
we present the recursive representation of the SARX-
Laguerre model obtained after the expansion of ev-
ery submodel on two Laguerre bases. In section 3,
firstly we focus on the parametric identification of
the SARX-Laguerre model where we propose a recur-
sive algorithm to estimate the Fourier coefficients so
that the regularized square error will be minimized.
Secondly, we use the genetic algorithm to optimize
the poles. In section 4, we present the horizon fault
estimation scheme of a linear hybrid system using
SARX-Laguerre model. Finally, the proposed MHE
using SARX-Laguerre model scheme is validated by
a numerical example.
2 SARX-Laguerre0odel
After the expansion of the strictly causal discrete time
linear hybrid systems known as Switched Auto Re-
gressive eXogenous (SARX) models on Laguerre or-
thonormal independent bases in the goal to reduce the
parameters complexity, we obtain the flowing recur-
sive representation of the new model entitled SARX-
Laguerre model.
X(k+ 1) = AλkX(k) + bλk
yy(k) + bλk
uu(k)
y(k) = (Cλk)TX(k)(1)
with, the state vector X(k)and the parameters
vector Cλkare a Nλk
a+Nλk
bdimensional vectors,
where Nλk
aand Nλk
bare the truncated orders, defined
as:
X(k) =
x0,y(k)
.
.
.
xNλk
a1,y(k)
x0,u(k)
.
.
.
xNλk
b1,u(k)
, Cλk=
gλk
0,a
.
.
.
gλk
Nλk
a1,a
gλk
0,b
.
.
.
gλk
Nλk
b1,b
(2)
and Aλkis a square matrix of dimension Nλk
a+
Nλk
bdefined as:
Aλk=Aλk
y0Nλk
a,Nλk
b
0Nλk
b,Nλk
aAλk
u(3)
where 0Nλk
a,Nλk
b
and 0Nλk
b,Nλk
aare two null matrices
of dimensions (Nλk
a×Nλk
b)and (Nλk
b×Nλk
a)re-
spectively, and Aλk
yand Aλk
uare two square matrices
of dimension Nλk
aand Nλk
brespectively. The matrix
Aλk
yis given, for r=t= 1, . . . , Na, by the following
relation:
Aλk
y(r, t) =
ξλk
aif r =t
(ξλk
a)(rt1)(1 (ξλk
a)2)if r > t
0if r < t
,
(4)
and the matrix Aλk
uis given, for r=t= 1, . . . , Nb,
as:
Aλk
u(r, t) =
ξλk
bif r =t
(ξλk
b)(rt1)(1 (ξλk
b)2)if r > t
0if r < t
(5)
The vectors bλk
yand bλk
uare of dimension Nλk
a+Nλk
b
and given by:
bλk
y=1(ξλk
a)2
1
ξλk
a
(ξλk
a)2
.
.
.
(ξλk
a)Na1
0Nb,1
(6)
bλk
u=1(ξλk
b)2
0Na,1
1
ξλk
b
(ξλk
b)2
.
.
.
i(ξλk
b)Nb1
(7)
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
110
Volume 15, 2024
with 0Na,1and 0Nb,1are two null vectors of dimension
Naand Nbrespectively.
We note that the SARX-Laguerre model is char-
acterized by (2 s)poles to be optimized and (s N)
Fourier Coefficients to be identified, where (N=
Nλk
a+Nλk
b). In this paper, we use a recursive proce-
dure to identify the Fourier coefficients and a meta-
heuristic algorithm to optimize the poles.
3 Identification of SARX-Laguerre
0odel
3.1 Recursive,dentification of Fourier
&oefficients
In this section, we propose a recursive method to iden-
tify the Fourier coefficients of the SARX-Laguerre
model. The proposed method consists in the mini-
mization of the regularized square error Jλh
reg (h)given
as:
Jλh
reg (h) =
h
k=1 ym(k)(Cλk)TX(k)2+
α
Nλh
a1
n=0 gλh1
n,a (h1) gλh
n,a(h)2+
Nλh
b1
n=0 gλh1
n,b (h1) gλh
n,b(h)2
,
(8)
where ym(k)is the measured output at time instant
kcorrespond to the output of the λkth subsystem,
gλh1
n,a (h1) and gλh1
n,b (h1) are the Fourier coeffi-
cients at time instant (h1) and gλh
n,a(h)and gλh
n,b(h)
are the Fourier coefficients at time instant (h), for
λh, λh1[1, . . . , s]. The parameter (α > 0)
acts on the importance given either to the first part or
to the second part of the criterion Jλh
reg (h). The cost
function Jλh
reg (h)given by (8) can be written in matrix
form as:
Jλh
reg (h) = Ym
hYλh
h2+αCλh1
h1Cλh
h2(9)
where the vectors Ym
hand Yλh
hcontains, for (k=
1, . . . , h), the λkth subsystem output and the λkth
submodel output respectively. Cλh1
h1and Cλh
hare
the parameter vectors of the λhth ARX-Laguerre sub-
model at time instants h1and h, where Cλh
his
defined as:
Cλh
h=gλh
0,a(h), . . . , gλh
Nλh
a1,a(h),
gλh
0,b(h), . . . , gλh
Nλh
b1,b(h)T(10)
The vector Yλh
hcan be written as:
Yλh
h=XhCλh
h(11)
where X(h)is a matrix regrouping the state vector
X(k)given by (2) of the λkth ARX-Laguerre sub-
model for (k= 1, . . . , h)determined recursively as:
Xh=Xh1
X(h)T(12)
Then, the cost function Jλh
reg (h)given by (9) can be
rewritten as:
Jλh
reg (h) = Ym
hXhCλh
hTYm
hXhCλh
h
+Cλh1
h1Cλh
hTαCλh1
h1Cλh
h(13)
Since the criterion (13) is quadratic in Cλh
h, its min-
imization with respect to this vector gives a global
minimum such as:
Jλh
reg (h)
Cλh
h
=
2XT
hYm
hXT
hXhCλh
hT
2αCλh1
h1Cλh
h
(14)
The vector
Cλh
hgrouping the estimated parameters
of the SARX-Laguerre model, for λk {1, . . . , s},
is calculated by the following relation:
Cλh
h=XT
hXh+αI1XT
hYm
h+α
Cλh1
h1
(15)
where Iis a (Nλk
a+Nλk
b)identity matrix.
3.2 SARX-Laguerre3oles2ptimization
using*enetic$lgorithms
Genetic Algorithms are considered as function opti-
mizers inspired by the evolution theory, and which
are applied to a broad range of optimization problems.
The implementation of a genetic algorithm requires
an initial population randomly generated. The opti-
mality of each solution from the population is quanti-
fied by computing a criterion named "fitness". Rank-
ing the fitness values assigned to the solutions al-
lows Genetic Algorithms to eliminate the bad ones
and make copies of the best ones, so that the popula-
tion size remains unchanged. This operator is named
"selection". From two good solutions "Parents", the
next operator named "Cross Over" gives rise to two
"Children" which are closer than "Parents" to the op-
timal solution. Changes then these new individuals
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
111
Volume 15, 2024
during the "Mutation" operator preventing the algo-
rithm to converge to a local optimum.
In this paper, we use the Genetic Algorithm to
optimize the poles of the SARX-Laguerre model.
Where, in the case of modeling linear hybrid system
decomposed to (s)subsystem, we are led to optimize
(2 s)poles and to identify (N S)Fourier coefficients,
where (N=Nλk
a+Nλk
b). To optimize the poles we
propose to use the genetic algorithm and to identify
the poles we use the recursive algorithm given in the
previous subsection. The proposed identification pro-
cedure of the SARX-Laguerre mode is illustrated by
Figure 1. The optimization begins by generating m
initial populations, each population is a set of (2 s)
poles. The best poles are used to compute the outputs
of the SARX-Laguerre model used in term to evaluate
the fitness in the next step. The aim of the optimiza-
tion is to select the best poles that minimize the Nor-
malized Mean Square Error "NMSE" between the
system outputs and the SARX-Laguerre model out-
put. The fitness function associated with the output is
given by:I
NMSE =M
k=1[ym(k)y(k)]2
M
k=1[ym(k)]2(16)
where ym(k)is the measured output at time instant k
and y(k)is the output model at time instant kgiven
by relation (1).
( )
m
y k
( )u k
( )y k
+
-
Switching hybrid
system ࡺࡹࡿࡱFitness
Poles optimization using
Genetic Algorithm
SARX-Laguerre
model
Fourier coefficients
identification
k
l
ǣǣሺ࢑ሻǣǣ
Fig. 1: Identification procedure of SARX-Laguerre
model
4 Moving Horizon Fault Estimation
based on SARX-Laguerre Model
In the following, we propose to use the moving hori-
zon fault estimation (MHE) based on the SARX-
Laguerre model to solve the fault estimation of the lin-
ear hybrid systems. The MHE is used to estimate the
actuator faults from the error between the estimated
output y(k)and the system output y(k). The SARX-
Laguerre model taking account of actuator faults can
be written from (1) as:
X(k+ 1) = AλkX(k) + bλk
yy(k)+
bλk
u(u(k) + f(k)),
y(k) = (Cλk)T
X(k),
(17)
where
X(k)and y(k)are the estimated state vector
and the estimated output model respectively, when
f(k)is the actuator faults. The estimate actuators
faults
f(k)is obtained by an online minimization
at every sample time of the quadratic cost function
Jf(k)defined as follows:
Jf(k) =
k
j=kNf+1
(y(j)ym(j))2(18)
with, Nfis the fault estimation horizon, y(j)is the es-
timated model output when the actuator fault is con-
sidered as given by (17) and ym(j)is the measured
system output. The quadratic criterion Jf(k)given
by (18) can be rewritten as:
Jf(k) =
k
j=kNf+1 (Cλj)TAλj
X(j1)+
bλj
yy(j1) + bλj
u(u(j1) +
f(j1))ym(j)2,
(19)
where
fis the fault actuator to be optimized using the
MHE method. The MHE method is formulated by
minimizing the criterion Jf(k)given by (19) at every
sample time. The criterion Jf(k)can be written in
matrix form as:
Jf(k) =
YNf(k)Ym,Nf(k)
2(20)
where the vector Ym,Nf(k)is determined from the
measured output system as follows:
Ym,Nf(k) = [ym(kNf+ 1), . . . , ym(k)]T(21)
and the vector
YNf(k)is determined from the model
output y(k)given by (17) as follows:
YNf(k) = y(kNf+ 1), . . . , y(k)T(22)
From relation (17), the vector
YNf(k)can be written
in matrix form as:
YNf(k) = CAX(k) + CByYNf(k)+
CBuUNf(k) + CBuFNf(k)(23)
with :
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
112
Volume 15, 2024
The matrices C,A,Byand Buare defined as:
C=
(Cλ(kNf+1) )T01,N . . . 01,N
01,N (Cλ(kNf+2) )T. . . 01,N
.
.
..
.
.....
.
.
01,N 01,N . . . (Cλk)T
(24)
A=
Aλ(kNf+1) 01,N . . . 01,N
0N,N Aλ(kNf+2) . . . 0N,N
.
.
..
.
.....
.
.
0N,N 0N,N . . . Aλk
(25)
By=
bλ(kNf+1)
y0(N,1) ... 0(N,1)
0(N,1) bλ(kNf+2)
y... 0(N,1)
.
.
..
.
.....
.
.
0(N,1) 0(N,1) . . . bλk
y
(26)
Bu=
bλ(kNf+1)
u0(N,1) . . . 0(N,1)
0(N,1) bλ(kNf+2)
u. . . 0(N,1)
.
.
..
.
.....
.
.
0(N,1) 0(N,1) . . . bλk
u
(27)
The vectors YNf(k)and UNf(k)are respectively
the vector of the output model without actuator
fault and the vector of the input system for (j=
kNf, k Nf+ 1, . . . , k)defined as follows:
YNf(k) =
y(kNf)
y(kNf+ 1)
.
.
.
y(k1)
,
UNf(k) =
u(kNf)
u(kNf+ 1)
.
.
.
u(k1)
(28)
The vector FNf(k)is regroup the actuator fault
for (j=kNf, k Nf+ 1, . . . , k)to be opti-
mized by minimizing the criterion Jf(k), defined
as:
FNf(k) =
f(kNf)
fNf(kNf+ 1)
.
.
.
fNf(k1)
(29)
The vector X(k)regroup the state vector X(j)
for (j=kNf+ 1, k Nf+ 2, . . . , k)
X(k) =
X(kNf)
X(kNf+ 1)
.
.
.
X(k1)
(30)
Then, the criterion Jf(k)given by (20) can be
written as:
Jf(k) =
CAX(k) + CByYNf(k)+
CBuUNf(k) + CBu
FNf(k)Ym,Nf(k)
2
(31)
The criterion (31) can be rewritten as:
Jf(k) = CAX(k) + CByYNf(k)+
CBuUNf(k)Ym,Nf(k)2+
CBu
FNf(k)2+ 2 CAX(k) + CByYNf(k)+
CBuUNf(k)Ym,Nf(k)CBu
FNf(k).
(32)
The criterion Jf(k)given by (32) is a quadratic in
terms of the actuator fault to be optimized
FNf(k).
In the case where constraints on the actuator fault
are taken into account, the optimized actuator fault is
given as follows:
F(k1) = min
F(k1)Sf
Jf(k)(33)
where the criterion Jf(k)is given by relation (32) and
F(k1) is defined as follows:
F(k1) =
f(kNf+ 1) . . .
f(k1) T
(34)
with Sfis the admissible set of actuator fault con-
straints defined as:
Sf=
F(k)/Γ
F(k1) V(35)
where Γand Vare defined as:
Γ = [ INfINf]T, V = [ Fmax Fmin ]T
(36)
with INfis a Nfdimensional identity matrix and
Fmin and Fmax are Nfdimensional vectors defined
as follows:
Fmin = [ fm, . . . , fm], Fmax = [ fM, . . . , fM]
(37)
fmand fMare the bounds of actuator fault that can
be chosen from the physical constraints.
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
113
Volume 15, 2024
5 Numerical Example
To validate the proposed moving horizon fault estima-
tion of hybrid systems using SARX-Laguerre model,
we consider a linear hybrid system composed of three
linear submodels. The system is given by a switched
ARX model (SARX) as follows:
y(k) = a1(λk)y(k1) + a2(λk)y(k2)+
a3(λk)y(k3) + b1(λk)u(k1)+
b2(λk)u(k2) + b3(λk)u(k3) + ε(t)
(38)
with λk {1,2,3}and the parameters a1,a2,a3,b1,
b2and b3are given as follows:
a1= [ 0.4 0.2 0.55 ]
a2= [ 0.25 0.35 0.2]
a3= [ 0.25 0.20.15 ]
b1= [ 0.85 1.15 0.5]
b2= [ 0.8 0.5 0.6]
b3= [ 0.25 0.5 0.4]
(39)
The SARX model corresponding to (38) contains 18
parameters. In the goal to reduce the parameters
number, we model this system bay a SARX-laguerre
where we decompose every ARX submodel on two
Laguerre bases with a truncated order (Nλk
a=Nλk
b=
2) for (λk= 1,2,3). Then, the obtained SARX-
Laguerre model is characterised by NL= 12 Fourier
coefficients to be identified and 6poles to be opti-
mized. To obtain the SARX-Laguerre model, we ap-
plied the proposed identification procedure using a
genetic algorithm to optimise the poles and a recur-
sive method to identify the Fourier coefficients. To
generate the identification data we use an excitation
input u(t)for M= 2000 observations as given in
Figure 2. In Figure 3 we draw the evolution of the
arbitrary switching signal and in Figure 4 we plot the
evolution of the SARX model output. For α= 0.5,
the optimal values of SARX-Laguerre poles and the
Fourier coefficients are given in Table 1. To vali-
Table 1.SARX-Laguerre parameters
λkPoles Fourier coefficients
ξλk
aξλk
b
1 -0.7478 0.5203 C1= [0.0541,0.0218,1.1794,0.0007]
2 0.0202 0.4113 C2= [0.3895,0.0016,1.2869,0.0030]
3 -0.0580 0.4414 C3= [0.5733,0.0029,1.3113,0.0023]
date the capability of the moving horizon fault estima-
tion of hybrid systems using SARX-Laguerre model
in fault estimation, we applied to the SARX system
given by (38) a constant input u= 1 and we added
to the input a constant fault at every case of switch-
ing signal as given in Figure 5. By considering con-
straints to the fault where fm= 0 and fM= 0.85 and
by applying the proposed fault estimation procedure,
we obtain the optimal fault as given in Figure 6.
0 500 1000 1500 2000
Number of iterations
0.5
0.6
0.7
0.8
0.9
u(k)
Input u(k)
Fig. 2: Evolution of input signal
0 500 1000 1500
Number of iterations
1
1.5
2
2.5
3
k
k
Fig. 3: Evolution of the arbitrary switching signal
6 Conclusion
In this paper, we propose a fault estimation scheme for
a class of linear hybrid systems using moving hori-
zon fault estimation. Linear hybrid systems are ap-
proximated by a reduced complexity SARX-Laguerre
model obtained after expanding the SARX model on
independent Laguerre orthonormal bases. To identify
the SARX-Laguerre model, we propose in this paper
to identify the Fourier coefficients using a recursive
algorithm and a metaheuristic algorithm to optimize
the poles. The identified model is used to develop
an online moving horizon fault estimation algorithm.
This algorithm can be used to develop a fault tolerant
adaptive control algorithm for nonlinear hybrid sys-
tems. The proposed fault estimation scheme is eval-
uated on a numerical example and the performances
are assessed in fault estimation.
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
114
Volume 15, 2024
Fig. 4: Evolution of the SARX model output
0
0.5
1
1.5
2
2.5
3
Switching signal
Added fault
0 50 100 150 200 250 300
Number of iterations
0
0.1
0.2
0.3
0.4
0.5
0.6
Added fault
Optimal fault
Fig. 6: Added fault and optimal fault
[3] A. Hasan, M. Tahavori, and H. S. Midtiby,
``Model-based fault diagnosis algorithms for
robotic systems,'' IEEE Access, vol. 11, pp.
2250--2258, 2023. [Online]. Available: https:
//doi.org/10.1109/ACCESS.2022.3233672
[4] C.-C. Peng and Y.-H. Cheng, ``Data driven
based modeling and fault detection for the mat-
lab/simulink turbofan engine: An ARX model
approach,'' in 2022 IEEE Conference on Control
Technology and Applications (CCTA). IEEE,
2022, pp. 498--503. [Online]. Available: https:
//doi.org/10.1109/CCTA49430.2022.9966098
[5] D. Lee and J. Jeong, ``Bearing fault detection
based on few-shot learning in siamese network,''
WSEAS Transactions on Systems, vol. 21, 2022.
[Online]. Available: https://doi.org/10.37394/
23202.2022.21.30
[6] S. A. Herve, Y. T. A. Aurelien, H. S. I.
Judith, and N. Leandre, ``Networked iterative
learning fault diagnosis algorithm for systems
with sensor random packet losses, time-varying
delays, limited communication and actuator
failure : Application to the hydroturbine
governing system,'' WSEAS TRANSACTIONS
on SYSTEMS and CONTROL, vol. 16, 2021.
[Online]. Available: https://doi.org/10.37394/
23203.2021.16.20
[7] S. X. Ding, ``Optimal fault detection and
estimation: a unified scheme and least squares
solutions,'' IFAC-PapersOnLine, vol. 51, no. 24,
pp. 465--472, 2018. [Online]. Available: https:
//doi.org/10.1016/j.ifacol.2018.09.618
Fig. 5: Switching signal and added fault
References:
[1] M. L. C. Peixoto, P. H. S. Coutinho, A. T.
Nguyen, T. M. Guerra, and R. M. Palhares,
``Fault estimation for nonlinear parameter-
varying time-delayed systems,'' Applied Mathe-
matics and Computation, vol. 465, 2024. [On-
line]. Available: https://doi.org/10.1016/j.amc.
2023.128405
[2] A. Prosvirin, F. Piltan, and J. Kim, ``Hy-
brid rubbing fault identification using a deep
learning-based observation technique,'' IEEE
Transactions on Neural Networks and Learn-
ing Systems, 11 2021. [Online]. Available:
https://doi.org/10.1109/TNNLS.2020.3027160
0 50 100 150 200 250 300
Number of iterations
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
115
Volume 15, 2024
[8] A. Mellit, G. Tina, and S. Kalogirou,
``Fault detection and diagnosis methods for
photovoltaic systems: A review,'' Renew-
able and Sustainable Energy Reviews, vol. 91,
pp. 1--17, 2018. [Online]. Available: https:
//doi.org/10.1016/j.rser.2018.03.062
[9] C. P. Raluca Rob, Gelu Ovidiu Tirian, ``In-
telligent signal processing application for real-
time faults detection,'' WSEAS Transactions on
Electronics, vol. 7, pp. 92--101, 2016. [On-
line]. Available: https://wseas.com/journals/
electronics/2016/a245817-089.pdf
[10] B. M. Lima, M. M. Morato, P. R. C. Mendes,
and J. E. Normey-Rico, ``Moving horizon
estimation of faults in renewable microgrids,''
IFAC-PapersOnLine, vol. 52, pp. 311--316,
2019. [Online]. Available: https://doi.org/10.
1016/j.ifacol.2019.06.080
[11] D. Krokavec and A. Filasová, ``Interval
observers for switched metzler-takagi-sugeno
fuzzy systems,'' WSEAS Transactions on Sys-
tems and Control, vol. 18, 2023. [Online]. Avail-
able: https://doi.org/10.37394/23203.2023.18.
10
[12] H. Torabi and H. Keshvri-Khor, ``Identifi-
cation of switched arx systems using an
iterative weighted least squares algorithm,''
International Journal of Industrial Electronics,
Control and Optimization, vol. 5, no. 1, pp.
23--30, 2022. [Online]. Available: https:
//doi.org/10.22111/ieco.2021.38734.1358
[13] K. K. Monfared, Y. Neyshabouri, A. Mire-
mad, S. Ahmadi, and H. Iman-Eini, ``Opti-
mal switching-sequence-based model predic-
tive control for a hybrid multilevel stat-
com,'' IEEE Transactions on Industrial Elec-
tronics, vol. 69, no. 10, pp. 9952--9960, 2022.
[Online]. Available: https://doi.org/10.1109/
TIE.2022.3146592
[14] G. Maślak and P. Orłowski, ``Microgrid opera-
tion optimization using hybrid system modeling
and switched model predictive control,'' Ener-
gies, vol. 15, no. 3, 2022. [Online]. Available:
https://www.mdpi.com/1996-1073/15/3/833
[15] E. P. van Horssen, D. Antunes, and M. Heemels,
``Switched LQG control for linear systems
with multiple sensing methods,'' Automat-
ica, vol. 103, pp. 217--229, 2019. [Online].
Available: https://doi.org/10.1016/j.automatica.
2019.01.036
[16] A. Mbarek and K. Bouzrara, ``Fault tolerant
control for MIMO nonlinear systems via mpc
based on mimo arx-laguerre multiple mod-
els,'' Mathematical Problems in Engineering,
vol. 2019, pp. 1--26, 2019. [Online]. Avail-
able: https://doi.org/10.1155/2019/9012182
[17] A. Sameh, M. Abdelkader, G. Tarek, and
R. José, ``Optimal multimodel representation
by laguerre filters applied to a communicating
two tank system,'' Journal of Systems Science
and Complexity, vol. 31, no. 3, pp. 621--646,
2018. [Online]. Available: https://doi.org/10.
1007/s11424-017-6047-2
[18] C. B. Njima and T. Garna, ``PIO output
fault diagnosis by ARX-laguerre model ap-
plied to 2nd order electrical system,'' IEEE
Access, vol. 8, pp. 83 052--83 061, 2020.
[Online]. Available: https://doi.org/10.1109/
ACCESS.2020.2990696
[19] S. Maraoui and K. Bouzrara, ``Paramet-
ric identification of linear systems us-
ing autoregressive with eXogenous input
model expansion on meixner-like orthonor-
mal bases,'' Transactions of the Institute of Mea-
surement and Control, vol. 42, no. 2, pp.
306--321, January 2020. [Online]. Available:
https://doi.org/10.1177/0142331219868621
[20] R. Singh and M. Sznaier, ``On identifica-
tion of nonlinear arx models with sparsity
in regressors and basis functions,'' IFAC-
PapersOnLine, vol. 54, no. 7, pp. 720--725,
2021, 19th IFAC Symposium on System Iden-
tification SYSID 2021. [Online]. Available:
https://doi.org/10.1016/j.ifacol.2021.08.446
[21] X. Song, Y. Li, T. Huang, and R. Chan,
``A sparse multiwavelet-based generalized la-
guerre–volterra model for identifying time-
varying neural dynamics from spiking activi-
ties,'' Entropy, vol. 19, no. 8, 2017. [Online].
Available: https://doi.org/10.3390/e19080425
[22] I. bengharat, A. Mbarek, and H. Dhouibi, ``Ex-
pansion of the switched FIR model parameters
on laguerre bases,'' in 2019 International Con-
ference on Signal, Control and Communication
(SCC), 2019, pp. 93--98. [Online]. Avail-
able: https://doi.org/10.1109/SCC47175.2019.
9116185
[23] X. Wang, T. Nguyen, and J. E. Schutt-Ainé,
``Laguerre–volterra feed-forward neural net-
work for modeling pam-4 high-speed links,''
IEEE Transactions on Components, Packaging
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
116
Volume 15, 2024
and Manufacturing Technology, vol. 10, no. 12,
pp. 2061--2071, 2020. [Online]. Available:
https://doi.org/10.1109/TCPMT.2020.3039486
[24] K. Bouzrara, T. Garna, J. Ragot, and
H. Messaoud, ``Decomposition of an ARX
model on Laguerre orthonormal bases,'' ISA
Transactions, vol. 51, no. 6, pp. 848-
-860, 2012. [Online]. Available: https:
//doi.org/10.1016/j.isatra.2012.06.005
[25] A. El Anes, K. Bouzrara, T. Garna,
and H. Messaoud, ``Expansion of MIMO
ARX model on laguerre orthonormal bases,''
in Electrical Engineering and Software Appli-
cations ICEESA, 2013. [Online]. Available:
https://doi.org/10.1109/ICEESA.2013.6578403
Contribution of Individual Authors to the
Creation of a Scientific Article (Ghostwriting
Policy)
The authors equally contributed to the present re-
search, at all stages from the formulation of the prob-
lem to the final findings and solution.
Sources of Funding for Research Presented in a
Scientific Article or Scientific Article Itself
No funding was received for conducting this study.
Conflicts of Interest
The authors have no conflicts of interest to declare.
Creative Commons Attribution License 4.0
(Attribution 4.0 International , CC BY 4.0)
This article is published under the terms of the
Creative Commons Attribution License 4.0
https://creativecommons.org/licenses/by/4.0/deed.en
_US
WSEAS TRANSACTIONS on ELECTRONICS
DOI: 10.37394/232017.2024.15.13
Ibtissem Bengharat, Abdelkader Mbarek
E-ISSN: 2415-1513
117
Volume 15, 2024