Next Article in Journal
Drone Secure Communication Protocol for Future Sensitive Applications in Military Zone
Previous Article in Journal
Low-Cost Chlorophyll Fluorescence Imaging for Stress Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distributed Mode-Dependent Event-Triggered Passive Filtering for Flexible Manipulator with Semi-Markov Parameters

1
School of Automation and Electrical Engineering, University of Science and Technology Beijing, Beijing 100083, China
2
Institute of Automation, Chinese Academy of Sciences, Beijing 100190, China
3
Huizhou Zhongke Advanced Manufacturing Research Center Co., Ltd., Huizhou 516000, China
*
Author to whom correspondence should be addressed.
Current address: 95 East Zhongguancun Road, Haidian District, Beijing 100000, China.
Sensors 2021, 21(6), 2058; https://doi.org/10.3390/s21062058
Submission received: 15 December 2020 / Revised: 25 February 2021 / Accepted: 26 February 2021 / Published: 15 March 2021
(This article belongs to the Section Sensors and Robotics)

Abstract

:
In this paper, the passive filtering problem of flexible robotic manipulator is investigated over sensor networks in a distributed manner from the control system perspective. The sensor networks are adopted to estimate true states of flexible robotic manipulator. In particular, the semi-Markov model is utilized for flexible manipulators with varying loads in unstructured environment, which is more flexible for practical implementations. Moreover, the new mode-dependent event-triggering mechanism is developed for distributed filter communications. Based on model transformation, sufficient conditions are first established to guarantee prescribed passive performance under disturbances. Then, desired mode-dependent filters are developed with the aid of convex optimization. In the end, several simulations results of a single-link flexible robotic manipulator are provided to verify the usefulness of the developed filtering algorithm.

1. Introduction

In the recent years, robotic manipulators have been more and more widely applied in various fields, such as industrial manufacturing [1], agricultural production [2,3], space exploration [4], etc. Especially, flexible robotic manipulators have attracted a huge attention by scientists and engineers since their distinguishing advantages, which include multi-degree of freedom, changeable structures, higher efficient payload and lower energy consumption [5,6,7]. Meanwhile, the manipulation tasks of robotic manipulators are becoming more and more complex due to practical demands. For instance, one interesting problem is how to achieve stable control of robotic manipulators with payloads of different weights in unstructured working environment. Furthermore, it should be pointed out that in certain tasks the changing of payloads could exhibit random features. Thus, some efforts have been made by applying the hybrid system models. To name a few, in [8], the trajectory control issue of robotic manipulators based on Radial Basis Function (RBF) neural networks is studied by the formulated hybrid switched system model. Furthermore, in [9], multiple models for robotic manipulator are discussed for the adaptive control problem. Besides, the load shifting robot manipulator model is investigated in [10]. It is worth mentioning that Markov jump systems are always used to model complex systems with parameter jumping or different modes, while the jumping modes and the transition probabilities between systems’ modes could be depicted accordingly [11,12,13]. As a result, it is reasonable to describe the flexible robotic manipulators with varying loads by Markov jump models. In addition, it is worth mentioning that the transition probabilities of Markov jump systems are based on fixed exponential distribution, which has ceratin restrictions in practical applications. In fact, the transition probabilities could be time-varying, which gives rises to the researches on so-called semi-Markov jump systems [14,15,16]. As such, semi-Markov jump systems have been intensively studied recently by researchers and various remarkable results on analysis and synthesis problems of semi-Markov jump systems have been reported in the literature. Therefore, it is important to study the flexible robotic manipulators with semi-Markov parameters.
On another active research area, with the substantial development of sensor and network technologies, sensor networks have been extensively studied in recent years, which means that a group of sensors are working collaboratively and changing information with each other. Compared with traditional single sensor, sensor networks can have many working merits, such as better robustness, higher efficiency, lower cost and so on. As one of the fundamental yet significant problems of sensor networks, the state estimation and filtering issues based on sensor networks have been widely studied to estimate the true values of certain states in target systems. In particular, various distributed filtering algorithms over sensor networks have been developed in recent years. Well known examples can be found by distributed H filtering, passivity filtering, dissipative filtering approaches and so on [17,18,19]. For instance, in [20], distributed H filtering is carried out over lossy sensor networks and dissipative filtering with missing measurements is investigated in [21]. Especially, the passivity performance has been proven to be effective for filtering problem from energy input-output perspective [22]. Very recently, by noticing that the common time-triggered strategies of distributed filtering need considerable communication consumption, the novel event-triggered strategies have been proposed to cope with the information exchanges among the sensors [23,24,25]. More precisely, the event-triggered strategies accomplish the signal transmission by certain prescribed event conditions instead of data transmission according to a period of time while a satisfactory task performance can still be ensured. For the distributed filtering of semi-Markov flexible robotic manipulators with diverse modes, the jumping mode information should be used to further reduce the conservatism. However, to the authors’ best knowledge, up to now, there is still little concern on the distributed filtering issues of semi-Markov jump systems within mode-dependent event-triggered framework, which remains an open and challenging problem in this research field.
Motivated by the aforementioned discussions, in this paper, the distributed mode-dependent event-triggered passive filtering algorithm is addressed for semi-Markov flexible manipulators. In comparison of most reported literature, our contributions are made towards the general filtering problem of semi-Markov jump systems with potential applications to flexible robotic manipulators, which can be concluded as follows:
(1) A new distributed mode-dependent event-triggering filtering algorithm is developed over sensor networks, which is based on the formulated semi-Markovian flexible manipulator, which can further extend the event-triggered strategy and can well utilize system mode information. The developed filtering algorithm is with asynchronous sampled-data and is thus more general and applicable in the practical applications.
(2) As an alternative to common H filtering method, the concept of passivity performance is adopted from an energy input-output perspective. By designing the distributed passive filter, the desired passivity performance can be achieved with external disturbance attenuation.
(3) On the basis of convex optimization and Lyapunov–Krasovskii methods, the mode-dependent filter gains and event-triggering parameters are designed simultaneously. At last, the effectiveness of our proposed filtering algorithm is demonstrated via a numerical example.
The remainder structure of our paper would be organized by the following parts. Section 2 formulates the distributed filtering problem of semi-Markovian flexible manipulator while mode-dependent event-triggered filter design is introduced. Section 3 presents the primary design method by proven details. Section 4 demonstrates our proposed algorithm by numerical simulations. In the final section, the paper is concluded and future research prospect is given.
The following notations are standardly used through the paper. R n denotes the n dimensional Euclidean space. Matrix P > 0 implies that P is positive definite. diag { } denotes the block-diagonal matrix. E ( ) represents the expectation operator. L 2 [ 0 , ) represents space of square-integrable vector functions over [ 0 , ) .

2. Preliminaries and Problem Formulation

2.1. Flexible Manipulator with Semi-Markov Parameters

Since flexible robotic manipulators always manipulate varying loads in unstructured environment, it is flexible and reasonable to model the parameter changing features caused by varying loads for practical implementations. For some applications, the load shifting can be modeled by semi-Markov chains, such that the parameter changing of flexible robotic manipulators can be conducted by semi-Markov chains accordingly.
In order to describe the semi-Markov chain, give a probability space ( O , F , P ) , where O denotes sample space, F represents δ -algebra of subsets of O and P stands for probability measure on F . Let { σ ( t ) , t 0 } denote a continuous-time discrete-state semi-Markov process on taking values in a finite set S = { 1 , , N } . Furthermore, the transition probability matrix Π : = ( π i j ( h ) ) , h > 0 , i , j I is defined in Equation (1)
Pr ( σ ( t + h ) = j | σ ( t ) = i ) = π i j ( h ) h + o ( h ) , i j , 1 + π i i ( h ) h + o ( h ) , i = j ,
where lim ( o ( h ) / h ) = 0 with o ( h ) denoting weak infinitesimal generator, π i j ( h ) 0 , i j , is the transition rate from mode i at time t to mode j at time t + h , satisfying π i i ( h ) = j = 1 , j i N π i j ( h ) , i S .
Consider the following flexible-link manipulator with the semi-Markov jumping dynamics and σ ( t ) denotes the semi-Markov chain in ( O , F , P ) , which is defined in Equation (2):
θ ˙ m = ω m , ω ˙ m = k J m ( σ ( t ) ) ( θ l θ m ) B J m ( σ ( t ) ) ω m + 1 J m ( σ ( t ) ) u , θ ˙ l = ω l , ω ˙ l = k J l ( σ ( t ) ) ( θ l θ m ) m ( σ ( t ) ) g h J l ( σ ( t ) ) sin ( θ l ) ,
where J m ( σ ( t ) ) denotes the inertia of motor, J l ( σ ( t ) ) denotes the inertia of link, θ m denotes the angular rotation of motor, θ l denotes the angular position of link, ω m denotes the angular velocity of motor, k denotes the joint elastic constant, m ( σ ( t ) ) denotes the link mass, l denotes the link length, g is the gravity constant, B is the viscosity, ω l denotes the angular velocity of link. More details of the parameters can be found in the literature [26].
Furthermore, by assuming certain u that can stabilize the manipulator and taking into account the external disturbance w 1 ( t ) L 2 [ 0 , ) , the following state-space nonlinear model can be obtained in Equation (3):
x ˙ ( t ) = A ( σ ( t ) ) x ( t ) + g ( σ ( t ) , x ( t ) ) + B ( σ ( t ) ) w 1 ( t ) z ( t ) = E ( σ ( t ) ) x ( t )
where x ( t ) = [ θ m , ω m , θ 1 , ω 1 ] T represents formulated manipulator system state, z ( t ) represents the output to be estimated over sensor network, w ( t ) denotes the disturbance and A ( σ ( t ) ) , B ( σ ( t ) ) , E ( σ ( t ) ) are defined in Equations (4)–(7)
A ( σ ( t ) ) = 0 1 0 0 k J m ( σ ( t ) ) B J m ( σ ( t ) ) k J m ( σ ( t ) ) 0 0 0 0 1 k J l ( σ ( t ) ) 0 k J l ( σ ( t ) ) 0 ,
g ( σ ( t ) , x ( t ) ) = 0 0 0 m ( σ ( t ) ) g h J 1 ( σ ( t ) ) sin ( θ 1 ) ,
B ( σ ( t ) ) = 0 1 J m ( σ ( t ) ) 0 0 ,
E ( σ ( t ) ) = 1 0 0 0 0 1 0 0 .
Moreover, it follows that g ( σ ( t ) , x ( t ) ) is a nonlinear function and can satisfy the Lipschitz nonlinear conditions with ι ( σ ( t ) ) = m ( σ ( t ) ) g h J 1 .
Remark 1.
The formulated state-space model can be generally utilized for nonlinear semi-Markov jump systems with Lipschitz conditions. Under this context, the distributed filtering problem can be further solved by convex optimization design.

2.2. Distributed Filter Design over Sensor Networks

For the sensor network, a direct graph G = V , E , A is presented to describe the communication topology of N sensor nodes, where V = v 1 , v 2 , , v N and E V × V stand for the set of nodes and edges, respectively, A = a i j R N × N represent the adjacency matrix with a i i = 0 for any i. A is associated with the edges of G are positive, i.e., a i j > 0 ε i j E . In addition, the corresponding Laplacian matrix of G is defined as L .
As depicted in Figure 1, it is supposed that all the sampler of sensors are time-driven with sampling instants t k and zero-order-holders are event-driven. The sampling period is set by h k = t k + 1 t k d ¯ . The measured output of sensor m can be obtained in Equation (8):
y m ( t ) = C m ( σ ( t ) ) x ( t ) + D m ( σ ( t ) ) w 2 ( t ) , m = 1 , 2 , , N ,
where y m ( t ) represents the measured output of sensor m, C m ( σ ( t ) ) R q and D m ( σ ( t ) ) R m are known constant matrices, w 2 ( t ) stands for the disturbance over sensors belonging to L 2 [ 0 , ) . Moreover, the event-triggered generator with event condition is deployed to transmit the latest data.
Consequently, denote x ^ m ( t ) as the estimated value of x ( t ) by sensor m and the sensor estimation error y ˜ m ( t ) can be obtained in Equation (9):
y ˜ m ( t ) = y m ( t ) C m ( σ ( t ) ) x ^ m ( t ) .
Then, y ˜ m ( t k ) is sampled at each t k and y ˜ m ( t δ m ) is transmitted according to the network topology when the event-triggering condition is satisfied. Besides, the broadcasting instant t δ + 1 m of sensor m satisfies the following event condition in Equation (10):
t δ + 1 i = min t k > t δ m t k y ˜ m ( t k ) y ˜ m ( t δ m ) κ m ( σ ( t ) ) y ˜ m ( t δ m ) ,
where 0 < κ m ( σ ( t ) ) < 1 denotes the threshold parameter.
Remark 2.
It can be found that our proposed event-triggered strategy can lead to asynchronous local information exchanges between the sensors, which is more practical and applicable for distributed sensor networks.
As such, the distributed mode-dependent filter can be designed in Equation (11):
x ^ m . ( t ) = A ( σ ( t ) ) x ^ m ( t ) + g ( σ ( t ) , x ^ m ( t ) ) + F m ( σ ( t ) ) y ˜ m ( t δ m ) + K m ( σ ( t ) ) m = 1 N a m n y ˜ n t δ m ( t ) j y ˜ m ( t δ m ) , z ^ m ( t ) = E ( σ ( t ) ) x ^ m ( t ) , t t δ m , t δ + 1 m , m j ( t ) = arg min p t t p m | t t p m .
where z ^ m ( t ) denotes the estimate of z ( t ) , F m ( σ ( t ) ) and K m ( σ ( t ) ) denotes mode-dependent filter gains to be designed.
By defining the sampled-data error in Equation (12):
ϵ m ( s k ) = y ˜ m t k y ˜ m t δ m , t δ m t k < t δ + 1 m ,
and dividing the interval t δ m , t δ + 1 m with t k = t δ m t δ + 1 m h k + 1 t k , t k + 1 , it can be derived in Equation (13) that
x ^ m . ( t ) = A ( σ ( t ) ) x ^ m ( t ) + g ( σ ( t ) , x ^ m ( t ) ) + F m ( σ ( t ) ) y ˜ m ( s k ) + K m ( σ ( t ) ) n = 1 N a m n y ˜ n t k y ˜ m ( t k ) F m ( σ ( t ) ) ϵ m ( t k ) K m ( σ ( t ) ) j = 1 N a m n ϵ n t k ϵ m ( t k ) , z ^ m ( t ) = E ( σ ( t ) ) x ^ m ( t ) , t t k , t k + 1 .
As a result, by letting e m ( t ) = x ( t ) x ^ m ( t ) and z ˜ m ( t ) = z ( t ) z ^ m ( t ) , the filtering error dynamics can be obtained in Equation (14):
e ˙ m ( t ) = A ( σ ( t ) ) e m ( t ) + g ( σ ( t ) , x ( t ) ) g ( σ ( t ) , x ^ m ( t ) ) F m ( σ ( t ) ) C m ( σ ( t ) ) e m ( t k ) K m ( σ ( t ) ) n = 1 N a m n C n ( σ ( t ) ) e n ( t k ) C m ( σ ( t ) ) e m ( t k ) + F m ( σ ( t ) ) ϵ m ( t k ) + K m ( σ ( t ) ) m = 1 N a m n ϵ j t k ϵ i ( t k ) + B ( σ ( t ) ) w 1 ( t ) F m ( σ ( t ) ) D m ( σ ( t ) ) w 2 ( t k ) K m ( σ ( t ) ) n = 1 N a m n D n ( σ ( t ) ) v 2 t k D m ( σ ( t ) ) v 2 t k , z ˜ m ( t ) = E ( σ ( t ) ) e m ( t ) , t t k , t k + 1 ,
which can be further rewritten in Equation (15):
e ˙ ( t ) = A ¯ ( σ ( t ) ) e ( t ) + G ( σ ( t ) , e ( t ) ) ( F ( σ ( t ) ) K ( σ ( t ) ) L ( σ ( t ) ) ) C ( σ ( t ) ) e ( t k ) + ( F ( σ ( t ) ) K ( σ ( t ) ) L ( σ ( t ) ) ) ϵ ( t k ) + B ¯ ( σ ( t ) ) w 1 ( t ) ( F ( σ ( t ) ) K ( σ ( t ) ) L ( σ ( t ) ) ) D ¯ ( σ ( t ) ) w 2 t k , z ˜ ( t ) = E ¯ ( σ ( t ) ) e ( t ) , t t k , t k + 1
where
e ( t ) = e 1 T ( t ) , e 2 T ( t ) , , e N T ( t ) T , z ˜ ( t ) = z ˜ 1 T ( t ) , z ˜ 2 T ( t ) , , z ˜ N T ( t ) T , ϵ ( t ) = ϵ 1 T ( t ) , ϵ 2 T ( t ) , , ϵ N T ( t ) T , A ¯ ( σ ( t ) ) = I N A ( σ ( t ) ) , G ( σ ( t ) , e ( t ) ) = [ ( g ( σ ( t ) , x ( t ) ) g ( σ ( t ) , x ^ 1 ( t ) ) ) T , ( g ( σ ( t ) , x ( t ) ) g ( σ ( t ) , x ^ 2 ( t ) ) ) T , ( g ( σ ( t ) , x ( t ) ) g ( σ ( t ) , x ^ N ( t ) ) ) T ] T , C ( σ ( t ) ) = d i a g C i ( σ ( t ) ) , F ( σ ( t ) ) = d i a g F i ( σ ( t ) ) , K ( σ ( t ) ) = d i a g K i ( σ ( t ) ) , L ( σ ( t ) ) = L ( σ ( t ) ) I , B ¯ ( σ ( t ) ) = 1 B ( σ ( t ) ) , D ¯ ( σ ( t ) ) = D 1 T ( σ ( t ) ) , D 2 T ( σ ( t ) ) , , D N T ( σ ( t ) ) T , E ¯ ( σ ( t ) ) = I E ( σ ( t ) ) .
Moreover, it can be verified that G ( σ ( t ) ) T ( e ( t ) ) G ( σ ( t ) ) ( e ( t ) ) Γ ( σ ( t ) ) e T ( t ) e ( t ) , Γ ( σ ( t ) ) = d i a g ι ( σ ( t ) ) , ι ( σ ( t ) ) , , ι ( σ ( t ) ) .

2.3. Filtering Objective

For denoting simplicity, denote mode σ ( t ) as i index and employ the novel input-delay method. Then, the overall filtering error dynamics can be deduced in Equation (16):
e ˙ ( t ) = A ¯ i e ( t ) + G i ( e ( t ) ) ( F i K i L i ) C i e ( t d ( t ) ) + ( F i K i L i ) ϵ ( t d ( t ) ) + B ¯ i w 1 ( t ) ( F i K i L i ) D ¯ i w 2 t d ( t ) , z ˜ ( t ) = E ¯ i e ( t ) , t t k , t k + 1
where d ( t ) = t t k , t t k , t k + 1 with 0 d ( t ) < d ¯ .
Before proceeding further, the following passivity performance index in Equation (17) is introduced for the distributed filters.
Definition 1.
Under the zero-initial condition, the modified passivity performance γ is said to be achieved in the mean-square sense, if it holds that
2 E { 0 T z ˜ T ( t ) S ( w 1 ( t ) + w 2 t d ( t ) ) d t } γ 0 T w 1 T ( t ) w 1 ( t ) + w 2 T t d ( t ) w 2 t d ( t ) d t ,
where matrix S is with appropriate dimension.
The control objective is to design the desired filter gains K i and F i for all filters, such that the passivity performance can be achieved according to Definition 1.
The following useful lemma is given for later derivations.
Lemma 1.
Ref. [27] For any matrix M > 0 , scalars τ > 0 , τ ( t ) satisfying 0 τ ( t ) τ , vector function x ˙ ( t ) : [ τ , 0 ] R n such that the concerned integrations are well defined in Equation (18), then
τ t τ t x ˙ T ( s ) M x ˙ ( s ) d s ζ T ( t ) U ζ ( t ) ,
where
ζ ( t ) = [ x T ( t ) , x T ( t τ ( t ) ) , x T ( t τ ) ] T , U = M M 0 * 2 M M * * M .

3. Distributed Filtering Analysis and Filter Gain Design

In this section, the mode-dependent distributed filters are designed with proven details.
Theorem 1.
For given h k and mode-dependent filter gains K i and F i , the passivity filtering can be achieved according to Definition 1, if there exist mode-dependent matrices P i > 0 , mode-dependent parameters κ i , and matrices Q > 0 , R > 0 , such that Ξ i , k < 0 , for all i N , k = 1 , 2 , , K , where
Ξ i , k : = Ξ 1 i , k Ξ 2 i , k * Ξ 3 i , k , Ξ 1 i , k : = 2 P i A ¯ i + Q R + Γ i + j = 1 N π i j , k P j P i F i C i + P i K i L i C i + R 0 P i F i P i K i L i * 2 R + C T ( Λ I q ) C R C T ( Λ I q ) * * Q R 0 * * * ( Λ I N ) I q , Ξ 2 i , k : = P i P i B ¯ i E ¯ i T P i F i D ¯ i + P i K i L i D ¯ i E ¯ i T S d ¯ A ¯ i T R 0 0 0 d ¯ C i T F i T R + d ¯ C i T L i T K i T R 0 0 0 0 0 0 0 d ¯ F i T R d ¯ L i T K i T R , Ξ 3 i , k : = I 0 0 d ¯ R * γ I 0 d ¯ B ¯ i T R * * γ I d ¯ D ¯ i T F i T R + d ¯ D ¯ i T L i T K i T R * * * R .
Proof. 
For each mode i, the mode-dependent Lyapunov-Krasovskii functionals are constructed in Equation (19):
V ( i , t ) = V 1 ( i , t ) + V 2 ( i , t ) + V 3 ( i , t ) ,
where
V 1 ( i , t ) = e T ( t ) P i e ( t ) , V 2 ( i , t ) = t d ¯ t e T ( φ ) Q e ( φ ) d φ , V 3 ( i , t ) = d ¯ d ¯ 0 t + φ t e ˙ T ( η ) R e ˙ ( η ) d η d φ .
The weak infinitesimal operator L of V ( i , t ) is defined in Equation (20):
L V ( i , t ) : = lim Δ 0 1 Δ { E { V ( σ ( t + Δ ) , t + Δ ) | σ ( t ) = i } V ( i , t ) } ,
with
lim Δ 0 1 Δ G i ( h + Δ ) G i ( h ) 1 G i ( h ) = 0 , lim Δ 0 1 Δ 1 G i ( h + Δ ) 1 G i ( h ) = 1 , lim Δ 0 1 Δ q i j ( G i ( h ) G i ( h + Δ ) ) Δ ( 1 G i ( h ) ) = π i j ( h ) ,
where h is the elapsed time at mode i, G i ( h ) represents the cumulative distribution function of the sojourn time, q i j denotes the probability intensity jumping from mode i to mode j and π i j ( h ) : = q i j π i ( h ) .
As a result, it can be deduced that
L V 1 ( i , t ) = lim 0 1 [ j = 1 , j i N Pr { σ t + = j | σ ( t ) = i } e T t + P j e t + + Pr σ t + = i σ ( t ) = i e T t + P i e t + e T t P i e t ] lim 0 1 [ j = 1 , j i N q i j G i ( h + ) G i ( h ) 1 G i ( h ) e T t + P j e t + + G i ( h + ) G i ( h ) 1 G i ( h ) e T t + P i e t + ] e T t P i e t = e ˙ T ( t ) P i e ( t ) + e T ( t ) P i e ˙ ( t ) + j = 1 N π i j ( h ) e T ( t ) P j e ( t ) = 2 e T ( t ) P i e ˙ ( t ) + j = 1 N π i j ( h ) e T ( t ) P j e ( t ) = 2 e T ( t ) P i ( A ¯ i e ( t ) + G i ( e ( t ) ) ( F i K i L i ) C i e ( t d ( t ) ) + ( F i K i L i ) ϵ ( t d ( t ) ) + B ¯ i w 1 ( t ) ( F i K i L i ) D ¯ i w 2 t d ( t ) ) + j = 1 N π i j ( h ) e T ( t ) P j e ( t )
L V 2 ( i , t ) = e T ( t ) Q e ( t ) e T ( t d ¯ ) Q e ( t d ¯ ) ,
L V 3 ( i , t ) = d ¯ d ¯ 0 e ˙ T ( t ) R e ˙ ( t ) d ¯ d ¯ 0 e ˙ T ( t + φ ) R e ˙ ( t + φ ) d φ = d ¯ 2 e ˙ T ( t ) R e ˙ ( t ) d ¯ t d ¯ t e ˙ T ( φ ) R e ˙ ( φ ) d φ .
In light of Lemma 1 as stated in Equation (18), one can obtain that
d ¯ t d ¯ t e ˙ T ( φ ) R e ˙ ( φ ) d φ e T ( t ) e T ( t d ( t ) ) e T ( t d ¯ ) T R R 0 * 2 R R * * R e T ( t ) e T ( t d ( t ) ) e T ( t d ¯ )
Moreover, from the Equation (16), one has
d ¯ 2 e ˙ T ( t ) R e ˙ ( t ) = η T ( t ) d ¯ A ¯ i T d ¯ C i T ( F i K i L i ) T 0 d ¯ ( F i K i L i ) T d ¯ I d ¯ B ¯ i T d ¯ D ¯ i T ( F i K i L i ) T R d ¯ A ¯ i T d ¯ C i T ( F i K i L i ) T 0 d ¯ ( F i K i L i ) T d ¯ I d ¯ B ¯ i T d ¯ D ¯ i T ( F i K i L i ) T T η ( t )
where
η ( t ) = e T ( t ) e T ( t d k ( t ) ) e T ( t d ¯ ) ϵ T ( t d ( t ) ) G i T ( e ( t ) ) w 1 T ( t ) w 2 T t d ( t ) T
From the definition of nonlinear function G i as stated in Equation (1), it yields that
G i T ( e ( t ) ) G i ( e ( t ) ) Γ i e T ( t ) e ( t ) ,
which means that
Γ i e T ( t ) e ( t ) G i T ( e ( t ) ) G i ( e ( t ) ) 0 .
The event-triggering function in Equation (10) implies that
m = 1 N ϵ m T ( t k ) ϵ m ( t k ) m = 1 N κ i 2 ( y ˜ m ( t k ) ϵ i ( t k ) ) T ( y ˜ m ( t k ) ϵ m ( t k ) )
which can lead to
e ( t d ( t ) ) ϵ ( t d ( t ) ) T C i T ( Λ i I ) C i C i T ( Λ i I ) * ( Λ i I ) I e ( t d ( t ) ) ϵ ( t d ( t ) ) 0 ,
where Λ i = d i a g { κ 1 2 , κ 2 2 , , κ N 2 } .
Then, it can be obtained by Schur complement that
L V ( i , t ) + Γ i e T ( t ) e ( t ) G i T ( e ( t ) ) G i ( e ( t ) ) 2 e T ( t ) E ¯ i T S ( w 1 ( t ) + w 2 t d ( t ) ) γ w 1 T ( t ) w 1 ( t ) γ w 2 T t d ( t ) w 2 t d ( t ) + e ( t d ( t ) ) ϵ ( t d ( t ) ) T C i T ( Λ i I ) C i C i T ( Λ i I ) * ( Λ i I ) I e ( t d ( t ) ) ϵ ( t d ( t ) ) η T ( t ) Ξ ˜ i η ( t ) ,
where
Ξ ˜ i : = Ξ ˜ 1 i Ξ ˜ 2 i * Ξ ˜ 3 i , Ξ ˜ 1 i : = 2 P i A ¯ i + Q R + Γ i + j = 1 N π i j ( h ) P j P i ( F i K i L i ) C i + R 0 * 2 R + C T ( Λ I q ) C R * * Q R , Ξ ˜ 2 i : = P i ( F i K i L i ) P i P i B ¯ i 2 E ¯ i T P i ( F i K i L i ) D ¯ i 2 E ¯ i T S d ¯ A ¯ i T C T ( Λ I q ) 0 0 0 d ¯ C i T ( F i K i L i ) T 0 0 0 0 0 , Ξ ˜ 3 i : = ( Λ I N ) I q 0 0 0 d ¯ ( F i K i L i ) T * I 0 0 d ¯ I * * γ I 0 d ¯ B ¯ i T * * * γ I d ¯ D ¯ i T ( F i K i L i ) T * * * * R 1 .
By performing congruent transformation to Ξ ˜ i by d i a g { I , I , I , I , I , I , I , R } , Ξ i can be obtained. Therefore, it can be verified that when Ξ i < 0 holds, the following inequality holds by integrating between 0 and T that
2 E { 0 T z ˜ T ( t ) S ( w 1 ( t ) + w 2 t d ( t ) ) d t } + γ 0 T w 1 T ( t ) w 1 ( t ) + w 2 T t d ( t ) w 2 t d ( t ) d t 0 .
Finally, by taking into account the time-varying dwell time h ( t ) [28], one has π i j ( h ) = 1 K λ k π i j , k , 1 K λ k = 1 , λ k 0 . This implies that if Ξ i , k < 0 holds, the passivity performance can be achieved according to Definition 1 in the mean-square sense, which completes the proof. □
Theorem 2.
For given h k , the passivity filtering can be achieved according to Definition 1, if there exist mode-dependent matrices P i > 0 , F ˜ i and K ˜ i , mode-dependent parameters κ i , and matrices Q > 0 , R > 0 , such that Θ i , k < 0 , for all i N , k = 1 , 2 , , K , where
Θ i , k : = Θ 1 i , k Θ 2 i , k * Θ 3 i , k , Θ 1 i , k : = 2 P i A ¯ i + Q R + Γ i + j = 1 N π i j , k P j F ˜ i C i + K ˜ i L i C i + R 0 F ˜ i K ˜ i L i * 2 R + C T ( Λ I q ) C R C T ( Λ I q ) * * Q R 0 * * * ( Λ I N ) I q , Θ 2 i , k : = P i P i B ¯ i E ¯ i T S F ˜ i D ¯ i + K ˜ i L i D ¯ i E ¯ i T d ¯ A ¯ i T P i 0 0 0 d ¯ C i T F ˜ i T + d ¯ C i T L i T K ˜ i T 0 0 0 0 0 0 0 d ¯ F ˜ i T d ¯ L i T K ˜ i T , Θ 3 i , k : = I 0 0 d ¯ P i * γ I 0 d ¯ B ¯ i T P i * * γ I d ¯ D ¯ i T F ˜ i T + d ¯ D ¯ i T L i T K ˜ i T * * * R 2 P i .
With the above feasible solutions, the mode-dependent filter gains K i and F i can be calculated by
F i = P i 1 F ˜ i , K i = P i 1 K ˜ i .
Proof. 
Based on matrix transformation, the proof can follow directly from Theorem 1. □
Remark 3.
It is worth mentioning that the derived conditions are in the form of strict LMIs, which can be easily solved by Matlab or other mathematical convex optimization softwares. Once the feasible solution is solved, then the corresponding filter gains can be obtained.

4. Illustrative Example

In this section, the simulation results are provided to validate the effectiveness of our derived results.
For the illustrative example, a numerical simulation is carried out in Matlab 2017a with numerical parameters.
Consider the flexible manipulator with following semi-Markov parameters according to system (3):
A 1 = 0 1 0 0 50 1.25 50 0 0 0 0 1 10 0 10 0 , g 1 ( x ( t ) ) = 0 0 0 5 sin ( θ 1 ) , B 1 = 0 5 0 0 , E 1 = 1 0 0 0 0 1 0 0 .
and
A 2 = 0 1 0 0 25 0.625 25 0 0 0 0 1 5 0 5 0 , g 2 ( x ( t ) ) = 0 0 0 4 sin ( θ 1 ) , B 2 = 0 2.5 0 0 , E 2 = 1 0 0 0 0 1 0 0 .
The sensor network consisted of 4 sensors are use with following parameters according to system (8) and its communication topology is depicted in Figure 2,
C 11 = 1.2 0 0 0 0 1.2 0 0 , D 11 = 0.5 0.5 , C 21 = 1.1 0 0 0 0 1.1 0 0 , D 21 = 0.6 0.6 , C 31 = 1.3 0 0 0 0 1.3 0 0 , D 31 = 0.4 0.4 , C 41 = 1.4 0 0 0 0 1.4 0 0 , D 41 = 0.3 0.3 ,
and
C 12 = 0.9 0 0 0 0 0.9 0 0 , D 12 = 0.2 0.2 , C 22 = 0.8 0 0 0 0 0.8 0 0 , D 22 = 0.4 0.3 , C 32 = 0.7 0 0 0 0 0.7 0 0 , D 32 = 0.6 0.4 , C 42 = 0.6 0 0 0 0 0.6 0 0 , D 42 = 0.7 0.6 ,
where the Laplacian can obtained as follows:
L = 2 1 1 0 0 2 1 1 0 0 1 1 1 0 0 1
Based on the above directed communication topology, sensor nodes 1–4 can exchange information accordingly. It is noticed that all sensors are with different parameters to verify the general applicability of our filtering design. In practical applications, the adopted sensor networks can also be chosen with same parameters.
In the simulation, the sampling period is randomly set by 0.1 s and 0.2 s, such that d ¯ = 0.2 s. The transition rates are supposed to be π 11 ( h ) [ 1.6 , 1.4 ] and π 22 ( h ) [ 1.9 , 1.1 ] , which implies that π 11 , 1 = 1.4 , π 11 , 2 = 1.6 , π 22 , 1 = 1.1 and π 22 , 2 = 1.9 with K = 2 . Moreover, the mode-dependent event-triggered scalars are supposed to be κ 1 = 0.1 and κ 2 = 0.2 . The passivity performance index is chosen by γ = 10 and the disturbances are set by w 1 ( t ) = 0.1 s i n ( t ) and w 2 ( t ) = 0.1 c o s ( t ) . With these parameters, the desired mode-dependent filter gains can be obtained by
K 11 = 0.0631 0.0022 0.2574 0.0358 0.0315 0.0021 0.0361 0.0063 , K 21 = 0.0611 0.0020 0.2371 0.0285 0.0299 0.0014 0.0321 0.0050 , K 31 = 0.0884 0.0024 0.2714 0.0296 0.0409 0.0016 0.0338 0.0056 , K 41 = 0.0988 0.0024 0.2703 0.0397 0.0437 0.0027 0.0336 0.0076 ,
and
K 12 = 0.0236 0.0022 0.1970 0.0037 0.0069 0.0007 0.0318 0.0019 , K 22 = 0.0114 0.0060 0.2303 0.0010 0.0149 0.0006 0.0354 0.0017 , K 32 = 0.0248 0.0108 0.2609 0.0118 0.0156 0.0029 0.0395 0.0005 , K 42 = 0.0471 0.0007 0.1247 0.0021 0.0021 0.0010 0.0203 0.0020 ,
and
F 11 = 0.0940 0.0068 0.9188 0.3561 0.0043 0.0158 0.1647 0.0651 , F 21 = 0.0761 0.0081 0.9721 0.3893 0.0110 0.0175 0.1716 0.0710 , F 31 = 0.1059 0.0044 0.7079 0.3739 0.0093 0.0172 0.1277 0.0674 , F 41 = 0.1148 0.0029 0.5616 0.3428 0.0202 0.0150 0.1046 0.0615 ,
and
F 12 = 1.0483 0.0050 0.7521 0.5298 0.3404 0.0511 0.2321 0.0889 , F 22 = 1.0430 0.0045 1.0643 0.5668 0.3113 0.0520 0.2824 0.0958 , F 32 = 1.0082 0.0081 1.1269 0.6160 0.2850 0.0549 0.2896 0.1038 , F 42 = 1.0100 0.0047 1.3132 0.6969 0.2634 0.0643 0.3202 0.1174 .
By random initial conditions for sensor nodes, Figure 3, Figure 4, Figure 5 and Figure 6 show the state trajectories of the filtering errors with disturbances. It can be seen that all the filtering errors can converge to zeros despite of system mode jumping and external disturbances, which means that the designed mode-dependent filters can effectively estimate the true state of flexible manipulator according to Definition 1. This also implies that the filtering errors e m ( t ) = x ( t ) x ^ m ( t ) for all sensors can be mean-square stable when there is no disturbance, which can verify our derived results. Figure 7, Figure 8, Figure 9 and Figure 10 depict event triggering instants and release interval of the sensors, where the release intervals are larger and the signal transmission among the sensors is event-triggered instead of traditional time-triggered schemes. It can be found that the developed event-triggered strategy can considerably decrease the numbers of communications compared with tradition time-triggered schemes (time-varying sampling periods of 0.1 s and 0.2 s). This shows the considerable advantages on decreasing signal transmissions among the sensor networks, where the event-triggered instants are larger than the tradition time-triggered instants. Thus, the simulation results can support our developed filter designs.

5. Conclusions

This paper discusses distributed filters design of flexible manipulator with semi-Markov parameters based on sensor networks. Moreover, the passivity performance is adopted to cope with the external disturbances of manipulator and sensor networks. By developing mode-dependent event-triggered schemes to achieve information exchanges, the distributed filtering can be accomplished in an asynchronous framework. By employing mode-dependent Lyapunov–Krasovskii approach, sufficient filtering conditions can be deduced and desired filter gains are designed, such that the passivity performance can be achieved. The correctness of our design method is finally demonstrated by a numerical example. In the future study, an interesting issue would be extending our current results to the cases with sampled semi-Markov processes, which means that the observed jumping modes may be different from the true system modes, which is more complex but more practical in the real world applications.

Author Contributions

Conceptualization; methodology; Y.J.; writing–original draft preparation, Y.J.; review and editing, S.Z.; funding acquisition, S.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Key Research and Development Program of China 2020AAA0105900, Huizhou Science and Technology Project 2020SB0106006, and Zhongkai Science and Technology Project 20201218A002.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. George Thuruthel, T.; Ansari, Y.; Falotico, E.; Laschi, C. Control strategies for soft robotic manipulators: A survey. Soft Robot. 2018, 5, 149–163. [Google Scholar] [CrossRef] [PubMed]
  2. Yu, S.; Yu, X.; Shirinzadeh, B.; Man, Z. Continuous finite-time control for robotic manipulators with terminal sliding mode. Automatica 2005, 41, 1957–1964. [Google Scholar] [CrossRef]
  3. Edan, Y.; Miles, G.E. Systems engineering of agricultural robot design. IEEE Trans. Syst. Man Cybern. 1994, 24, 1259–1265. [Google Scholar] [CrossRef]
  4. Moosavian, S.A.A.; Papadopoulos, E. Free-flying robots in space: An overview of dynamics modeling, planning and control. Robotica 2007, 25, 537. [Google Scholar] [CrossRef] [Green Version]
  5. Chen, W. Dynamic modeling of multi-link flexible robotic manipulators. Comput. Struct. 2001, 79, 183–195. [Google Scholar] [CrossRef]
  6. Chaoui, H.; Sicard, P.; Gueaieb, W. ANN-based adaptive control of robotic manipulators with friction and joint elasticity. IEEE Trans. Ind. Electron. 2009, 56, 3174–3187. [Google Scholar] [CrossRef]
  7. Kiang, C.T.; Spowage, A.; Yoong, C.K. Review of control and sensor system of flexible manipulator. J. Intell. Robot. Syst. 2015, 77, 187–213. [Google Scholar] [CrossRef]
  8. Yu, L.; Fei, S.; Huang, J.; Gao, Y. Trajectory switching control of robotic manipulators based on RBF neural networks. Circuits Syst. Signal Process. 2014, 33, 1119–1133. [Google Scholar] [CrossRef]
  9. Ciliz, M.K.; Narendra, K.S. Adaptive control of robotic manipulators using multiple models and switching. Int. J. Robot. Res. 1996, 15, 592–610. [Google Scholar] [CrossRef]
  10. Wang, X.; Zhao, J. Autonomous switched control of load shifting robot manipulators. IEEE Trans. Ind. Electron. 2017, 64, 7161–7170. [Google Scholar] [CrossRef]
  11. De Farias, D.P.; Geromel, J.C.; Do Val, J.B.; Costa, O.L.V. Output feedback control of Markov jump linear systems in continuous-time. IEEE Trans. Autom. Control. 2000, 45, 944–949. [Google Scholar] [CrossRef]
  12. Shen, H.; Men, Y.; Wu, Z.G.; Cao, J.; Lu, G. Network-based quantized control for fuzzy singularly perturbed semi-Markov jump systems and its application. IEEE Trans. Circuits Syst. I Regul. Pap. 2018, 66, 1130–1140. [Google Scholar] [CrossRef]
  13. Li, H.; Shi, P.; Yao, D. Adaptive sliding-mode control of Markov jump nonlinear systems with actuator faults. IEEE Trans. Autom. Control. 2016, 62, 1933–1939. [Google Scholar] [CrossRef]
  14. Qi, W.; Zong, G.; Karimi, H.R. Sliding mode control for nonlinear stochastic singular semi-Markov jump systems. IEEE Trans. Autom. Control. 2019, 65, 361–368. [Google Scholar] [CrossRef]
  15. Shen, H.; Dai, M.; Yan, H.; Park, J.H. Quantized output feedback control for stochastic semi-Markov jump systems with unreliable links. IEEE Trans. Circuits Syst. II Express Briefs 2018, 65, 1998–2002. [Google Scholar] [CrossRef]
  16. Zong, G.; Ren, H. Guaranteed cost finite-time control for semi-Markov jump systems with event-triggered scheme and quantization input. Int. J. Robust Nonlinear Control. 2019, 29, 5251–5273. [Google Scholar] [CrossRef]
  17. Wu, L.; Shi, P.; Gao, H.; Wang, C. H filtering for 2D Markovian jump systems. Automatica 2008, 44, 1849–1858. [Google Scholar] [CrossRef]
  18. Sathishkumar, M.; Sakthivel, R.; Kwon, O.; Kaviarasan, B. Finite-time mixed H and passive filtering for Takagi–Sugeno fuzzy nonhomogeneous Markovian jump systems. Int. J. Syst. Sci. 2017, 48, 1416–1427. [Google Scholar] [CrossRef]
  19. Mathiyalagan, K.; Park, J.H.; Sakthivel, R. Robust reliable dissipative filtering for networked control systems with sensor failure. IET Signal Process. 2014, 8, 809–822. [Google Scholar] [CrossRef]
  20. Dong, H.; Wang, Z.; Gao, H. Distributed H filtering for a class of Markovian jump nonlinear time-delay systems over lossy sensor networks. IEEE Trans. Ind. Electron. 2012, 60, 4665–4672. [Google Scholar] [CrossRef]
  21. Zhang, M.; Shen, C.; Wu, Z.G.; Zhang, D. Dissipative filtering for switched fuzzy systems with missing measurements. IEEE Trans. Cybern. 2019, 50, 1931–1940. [Google Scholar] [CrossRef] [PubMed]
  22. Ding, D.; Wang, Z.; Shen, B. Recent advances on distributed filtering for stochastic systems over sensor networks. Int. J. Gen. Syst. 2014, 43, 372–386. [Google Scholar] [CrossRef]
  23. Zhang, X.M.; Han, Q.L.; Zhang, B.L. An overview and deep investigation on sampled-data-based event-triggered control and filtering for networked systems. IEEE Trans. Ind. Inform. 2016, 13, 4–16. [Google Scholar] [CrossRef]
  24. Tian, E.; Wang, Z.; Zou, L.; Yue, D. Probabilistic-constrained filtering for a class of nonlinear systems with improved static event-triggered communication. Int. J. Robust Nonlinear Control. 2019, 29, 1484–1498. [Google Scholar] [CrossRef]
  25. Meng, X.; Chen, T. Event triggered robust filter design for discrete-time systems. IET Control. Theory Appl. 2014, 8, 104–113. [Google Scholar] [CrossRef]
  26. Sahab, A.R.; Modabbernia, M.R. Backstepping method for a single-link flexible-joint manipulator using genetic algorithm. Int. J. Innov. Comput. Inf. Control 2011, 7, 4161–4170. [Google Scholar]
  27. Park, P.; Ko, J.W.; Jeong, C. Reciprocally convex approach to stability of systems with time-varying delays. Automatica 2011, 47, 235–238. [Google Scholar] [CrossRef]
  28. Ma, C.; Wu, W.; Li, Y. Distributed mode-dependent state estimation for semi-Markovian jumping neural networks via sampled data. Int. J. Syst. Sci. 2019, 50, 216–230. [Google Scholar] [CrossRef]
Figure 1. Illustration of distributed filtering on sensor network.
Figure 1. Illustration of distributed filtering on sensor network.
Sensors 21 02058 g001
Figure 2. The topology of sensor network.
Figure 2. The topology of sensor network.
Sensors 21 02058 g002
Figure 3. Filtering errors of sensor 1.
Figure 3. Filtering errors of sensor 1.
Sensors 21 02058 g003
Figure 4. Filtering errors of sensor 2.
Figure 4. Filtering errors of sensor 2.
Sensors 21 02058 g004
Figure 5. Filtering errors of sensor 3.
Figure 5. Filtering errors of sensor 3.
Sensors 21 02058 g005
Figure 6. Filtering errors of sensor 4.
Figure 6. Filtering errors of sensor 4.
Sensors 21 02058 g006
Figure 7. Event-triggered instants of sensor 1.
Figure 7. Event-triggered instants of sensor 1.
Sensors 21 02058 g007
Figure 8. Event-triggered instants of sensor 2.
Figure 8. Event-triggered instants of sensor 2.
Sensors 21 02058 g008
Figure 9. Event-triggered instants of sensor 3.
Figure 9. Event-triggered instants of sensor 3.
Sensors 21 02058 g009
Figure 10. Event-triggered instants of sensor 4.
Figure 10. Event-triggered instants of sensor 4.
Sensors 21 02058 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ji, Y.; Zheng, S. Distributed Mode-Dependent Event-Triggered Passive Filtering for Flexible Manipulator with Semi-Markov Parameters. Sensors 2021, 21, 2058. https://doi.org/10.3390/s21062058

AMA Style

Ji Y, Zheng S. Distributed Mode-Dependent Event-Triggered Passive Filtering for Flexible Manipulator with Semi-Markov Parameters. Sensors. 2021; 21(6):2058. https://doi.org/10.3390/s21062058

Chicago/Turabian Style

Ji, Yidao, and Suiwu Zheng. 2021. "Distributed Mode-Dependent Event-Triggered Passive Filtering for Flexible Manipulator with Semi-Markov Parameters" Sensors 21, no. 6: 2058. https://doi.org/10.3390/s21062058

APA Style

Ji, Y., & Zheng, S. (2021). Distributed Mode-Dependent Event-Triggered Passive Filtering for Flexible Manipulator with Semi-Markov Parameters. Sensors, 21(6), 2058. https://doi.org/10.3390/s21062058

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