Next Article in Journal
Stability in the Sense of Hyers–Ulam–Rassias for the Impulsive Volterra Equation
Next Article in Special Issue
Magnetically Suspended Control Sensitive Gyroscope Rotor High-Precision Deflection Decoupling Method using Quantum Neural Network and Fractional-Order Terminal Sliding Mode Control
Previous Article in Journal
Geochemical Characteristics of Deep-Sea Sediments in Different Pacific Ocean Regions: Insights from Fractal Modeling
Previous Article in Special Issue
Improved Results on Delay-Dependent and Order-Dependent Criteria of Fractional-Order Neural Networks with Time Delay Based on Sampled-Data Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Privacy Preservation of Nabla Discrete Fractional-Order Dynamic Systems

1
School of Automation Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China
2
Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China
3
Institute of Electronic and Information Engineering, University of Electronic Science and Technology of China, Dongguan 523808, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(1), 46; https://doi.org/10.3390/fractalfract8010046
Submission received: 4 December 2023 / Revised: 8 January 2024 / Accepted: 10 January 2024 / Published: 11 January 2024

Abstract

:
This article investigates the differential privacy of the initial state for nabla discrete fractional-order dynamic systems. A novel differentially private Gaussian mechanism is developed which enhances the system’s security by injecting random noise into the output state. Since the existence of random noise gives rise to the difficulty of analyzing the nabla discrete fractional-order systems, to cope with this challenge, the observability of nabla discrete fractional-order systems is introduced, establishing a connection between observability and differential privacy of initial values. Based on it, the noise magnitude required for ensuring differential privacy is determined by utilizing the observability Gramian matrix of systems. Furthermore, an optimal Gaussian noise distribution that maximizes algorithmic performance while simultaneously ensuring differential privacy is formulated. Finally, a numerical simulation is provided to validate the effectiveness of the theoretical analysis.

1. Introduction

Fractional calculus is a generalization of classical calculus that extends derivatives and integrals from integer to non-integer orders. Due to their non-locality and hereditary properties, fractional calculus and fractional systems have drawn much attention and have been extensively studied in many fields, including signal processing [1], soft robots [2], medical science [3], economics [4], biology [5], and so on. Notably, its application in modeling various system dynamics and engineering scenarios demonstrates superior efficiency over traditional integer-order methods, particularly within the domain of control systems [6,7,8].
It is worth pointing out that most research on fractional-order systems has mainly focused on continuous-time models. However, with the rapid development of digital computers, discrete fractional calculus may play a more important role than its continuous-time counterparts in today’s digital era. As a result, there has been a growing interest in the study of discrete fractional calculus. For instance, the authors in [9] utilized the Grünwald–Letnikov definition and designed a finite-dimensional controller for discrete fractional-order systems with state disturbance. A modified Mikhailov stability criterion was presented to reduce the computational complexity for a class of non-commensurate discrete-time fractional-order systems in [10]. Recently, significant progress has been made in discrete fractional calculus based on the fractional sum and difference under various definitions [11,12,13] including the Riemann–Liouville and Caputo definitions, respectively. However, the aforementioned literature used delta fractional difference to describe discrete fractional-order systems. This approach limits the general applicability of the results by setting the initial state to zero and rendering the time series to rely on the order of the system. To overcome these limitations, the authors in [14] introduced nabla fractional differences and a discrete nabla analog of the Laplace transform. In [15], the fractional difference of Lyapunov inequalities and a stability analysis of nabla discrete fractional systems were developed. Furthermore, these results were extended to non-linear nabla fractional-order systems [16]. A discrete Mittag–Leffler stability was proposed in [17] to characterize the convergence rule for nabla discrete fractional-order systems and extend the Lyapunov directed method to discrete fractional systems under various definitions. In [18], the authors investigated the asymptotic stability of a non-linear nabla variable-order system. Furthermore, foundational mathematical properties such as the nabla Taylor series and Laplace transform for nabla tempered fractional calculus were explored in [19]. In control systems, the authors designed a discrete-time fractional-order proportional–integral–derivative (FOPID) controller in [20] using finite impulse response FIR filters and model order reduction methods. Additionally, the projective synchronization problem was investigated in [21] for nabla discrete fractional-order neural networks, according to the quaternion theory. By implementing a delay feedback controller, the quasi-uniform synchronization problem for discrete fractional fuzzy neural networks with mixed time delays was first discussed in [22]. The concepts of nabla fractional calculus and nabla fractional order systems have also been extensively studied, particularly in coordination control of multi-agent systems [23,24], as well as in distributed optimization problems [25]. The above research has made great contributions to the development of nabla discrete fractional-order systems, but they do not take into account the issue of privacy leakage within these systems.
Privacy preservation is a crucial aspect of information security, particularly for dynamic systems. In recent years, various methods have been developed for privacy protection, such as k-anonymity [26], l-diversity [27], and secure multiparty computation [28]. However, these approaches typically require a priori knowledge of the attackers’ background information, which may not be available in practice. To address this problem, differential privacy, a method introducing random noise into the data, has become a robust privacy-protection approach. For instance, the authors in [29] provided a comprehensive analysis of differential privacy techniques and demonstrated their application in three practical cases. Recently, differential privacy has been extensively studied in dynamic and cyber–physical systems. The concept of differential privacy was first introduced for dynamic systems in [30], and the authors discussed the Kalman filtering problem with differential privacy requirements. The average privacy-preserving consensus problem of multi-agent systems was investigated in [31], in which the differential privacy of initial values was taken into consideration. The authors in [32] developed a differentially private linear–quadratic control protocol for multi-agent systems and provided quantitative criteria for privacy preservation. Two differentially private gradient methods were proposed for distributed optimization problems in [33] to ensure both privacy and optimization accuracy. For some non-numerical or symbolic data, a novel differential privacy mechanism was proposed in [34] for protecting sensitive symbolic system trajectories and the results were applied to Markov chains. From the point of systems and control theory, the authors in [35] investigated the relation between differential privacy and strong input observability of the systems and designed a privacy-preserving controller for a tracking problem. Furthermore, by assuming that certain initial values are sensitive and the rest are public, a necessary and sufficient condition is established in [36] for preserving the differential privacy and unobservability of linear dynamical systems. However, it should be noted that these studies have primarily focused on integer-order systems where the issue of privacy preservation in nabla fractional-order dynamic systems has not been studied yet.
Inspired by the above discussion, this paper studies privacy protection in nabla discrete fractional-order dynamic systems via differential privacy. The main contributions of this paper are in three aspects:
  • In contrast to the existing privacy-preserving literature that only concentrates on integer-order systems, this paper introduces a differential privacy mechanism for nabla discrete fractional-order dynamic systems. This method introduces Gaussian random noise into the output to prevent sensitive information from being inferred by malicious attackers.
  • A connection between the differential privacy of initial values and the observability of nabla discrete fractional-order systems is provided. Furthermore, the paper proposes a sufficient condition for maintaining differential privacy via the observability Gramian matrix, effectively addressing the analytical challenges introduced by random noise in nabla fractional-order systems.
  • An optimal Gaussian noise distribution is derived from an optimization problem. This distribution ensures the achievement of the desired privacy level while maintaining system performance.
The rest of this paper can be organized as follows: Section 2 introduces some notations and basic concepts of nabla fractional calculus. Section 3 formulates the problem statement for nabla discrete fractional systems. Section 4 analyzes the differential privacy for nabla fractional order systems. Section 5 presents a numerical example, and Section 6 concludes the paper.

2. Preliminaries

2.1. Notations

N , R and Z ( Z + ) denotes the natural number set, real number set, and integers (positive integers) set, respectively. For any x 1 , , x m R n , we denote [ x 1 ; ; x m ] as a vector [ x 1 T x m T ] T R n m . For a vector x R n , denote | | x | | p as its p-norm. For a matrix A R n × n , denote | | A | | 2 as its 2-norm, and its maximum and minimum eigenvalues are denoted by λ m a x ( A ) and λ m i n ( A ) , respectively. Further, A 0 means that A is symmetric and positive-definite. For any a , b R , N a + 1 = { a + 1 , , } . For a random variable ω , ω N n ( μ , Σ ) means that ω has a non-degenerate multivariate Gaussian distribution with mean μ R n and covariance Σ 0 , and its distribution has the following probability density:
p ( ω ; μ , Σ ) = 1 ( 2 π ) n d e t ( Σ ) 1 / 2 e | ω μ | 2 Σ 1 / 2 .
The so-called Q-function is defined by Q ( ω ) = 1 2 π ω e v 2 2 d v , where Q ( ω ) < 1 / 2 for ω > 0 , and R ( ϵ , δ ) = ( Q 1 ( δ ) + ( Q 1 ( δ ) ) 2 + 2 ϵ ) / 2 ϵ .

2.2. Discrete Fractional Calculus

Some basic concepts and lemmas about discrete fractional calculus are presented.
Definition 1.
(Gamma function) For any k R { , 2 , 1 , 0 } , the gamma function is defined by
Γ ( k ) = 0 e s s k 1 d s .
Definition 2.
(Rising function) The rising function is defined by
t α ¯ = Γ ( k + α ) k , k R { , 2 , 1 , 0 } .
Definition 3.
For f : N a + 1 n R , its n-th integer-order backward difference is defined by
n f ( k ) = i = 0 n ( 1 ) i n i f ( k i ) , k N a + 1 ,
where n Z + , N a + 1 = { a + 1 , a + 2 , } , a R and n i Γ ( p + 1 ) Γ ( q + 1 ) Γ ( p q + 1 ) .
Definition 4.
For f : N a + 1 n R , its α-th sum with the Riemann–Liouville definition is defined as
a R k α f ( k ) j = 0 k a 1 ( 1 ) j α j f ( k j ) , k N a + 1 ,
where α ( n 1 , n ) , n Z + , a R .
Then, Equation (4) equals to the following expression
a R k α f ( k ) = j = 0 k a 1 ( j + 1 ) α 1 ¯ Γ ( α ) f ( k j ) = j = a + 1 k ( k j + 1 ) α 1 ¯ Γ ( α ) f ( j ) .
Definition 5.
For f : N a + 1 n R , the nabla Riemann–Liouville and Caputo fractional differences are defined as
a R k α f ( k ) = a R k α n n f ( k ) , k N a + 1 ,
a C k α f ( k ) = n a R k α n f ( k ) , k N a + 1 ,
where α ( n 1 , n ) , n Z + , a R .

3. Problem Statement

Consider the following linear nabla discrete fractional-order systems with the Riemann–Liouville definition
0 R k α x ( k ) = A x ( k ) + B u ( k ) , k N 1 , y ( k ) = C x ( k ) + D u ( k ) ,
where 0 < α < 1 , y ( k ) R p is the output, and u ( k ) R m is the input. A R n × n , B R n × m , C R p × n , and D R p × m are the known constant matrices, respectively. In addition, x ( k ) R n is the pseudo-state, which can be calculated using the following formula [14]
x ( k ) = F A , α ( k ) x 0 + s = 1 k F A , α ( k s + 1 ) B u ( s ) ,
where x ( 0 ) = x 0 is the initial state and F A , α ( k ) = n = 0 k A n ( k n + 1 ) ( n + 1 ) α 1 ¯ Γ ( ( n + 1 ) α ) .
In this paper, we will focus on a finite dataset, which is consistent with studying the properties of a nabla dynamical system over a finite time. Given y ( 0 ) , u ( 0 ) , and a finite time L Z + , we denote the output sequence Y L : = [ y ( 0 ) ; ; y ( L ) ] R ( L + 1 ) p , and the input sequence U L : = [ u ( 0 ) ; ; u ( L ) ] R ( L + 1 ) m . Then, the output sequence Y L R ( L + 1 ) p of (8) can be written as
Y L = O L x 0 + H L U L ,
where O L R ( L + 1 ) p × n and H L R ( L + 1 ) p × ( L + 1 ) m are
O L = C C F ( 1 ) C F ( 2 ) C F ( L ) T ,
H L = D 0 0 0 C F ( 1 ) B + D 0 C F ( 2 ) B C F ( 1 ) B + D 0 0 C F ( L ) B C F ( L 1 ) B C F ( 1 ) B + D ,
F ( L ) = F A , α ( L ) = n = 0 L ( L n + 1 ) α n + α 1 ¯ Γ ( α n + 1 ) A n .
Throughout this paper, we consider the initial values x ( 0 ) = x 0 of system (8) to be privacy-sensitive.

4. Differential Privacy for Nabla Discrete Fractional-Order Systems

In this section, we study the differential privacy of nabla discrete fractional-order systems from three aspects. First, we analyze differential privacy based on a Gaussian mechanism that perturbs system outputs with random noise. Second, we give rank conditions for the observability of nabla discrete fractional-order systems and investigate the connection between observability and differential privacy. In particular, we characterize how the observability Gramian matrix impacts privacy guarantees. Finally, we formulate and solve an optimization problem to determine the minimum output noise quantity required to satisfy a specified level of differential privacy.

4.1. Differential Privacy with Output Noise

To prevent x 0 from being inferred from the published output data, random noise w ( k ) R p is added to the output such that y w , k = y ( k ) + w ( k ) , and we denote the noise sequence W L : = [ w ( 0 ) ; ; w ( L ) ] R ( L + 1 ) p . Thus, the mapping H : R n R ( L + 1 ) p from the initial value x 0 to the output sequence Y w , K can be described as
Y w , L = H ( x 0 ) : = O L x 0 + H L U L + W L .
Next, some definitions of differential privacy are given.
Definition 6.
Given c > 0 and r Z + , a pair of input data x 0 , x 0 R n is said to be c-adjacent and is denoted by A d j r c if
| | x 0 x 0 | | r c .
Then, we define the differential privacy of the mechanism (14).
Definition 7.
Given ϵ > 0 and δ 0 , the mechanism (14) is said to be ( ϵ , δ ) -differentially private at a finite time instant L Z + , if, for all S r a n g e ( H ) ,
P [ H ( x 0 ) S ] e ϵ P [ H ( x 0 ) S ] + δ .
Theorem 1.
The mechanism (14) induced by W L N ( L + 1 ) p ( μ , Σ w ) is ( ϵ , δ ) -differentially private at a finite time L Z + for A d j 2 c with ϵ > 0 and 0 δ 1 , if the covariance matrix Σ w 0 is chosen such that
λ m a x 1 / 2 ( O Σ , L ) c R ( ϵ , δ ) ,
where O Σ , L is defined by
O Σ , L = O L T Σ w 1 O L .
Proof. 
For any given ( x 0 , x 0 ) A d j 2 c , define
p L = ( L + 1 ) p , ζ = Y L Y L , Y L = O L x 0 + H L U L , Y L = O L x 0 + H L U L .
For any ϵ > 0 and S ( R p L ) , it follows that
P ( Y L + W L S ) = R p L 1 S ( Y L + w ) p W L ( w ) d w = S p W L ( z Y L ζ ) d z = 1 ( 2 π ) p K d e t ( Σ w ) S e x p | | z Y L ζ | | Σ w 1 2 2 d z .
Next, standard computation yields
| | z Y L ζ | | Σ w 1 2 2 = | | z Y L | | Σ w 1 2 2 + 2 ( z Y L ) T Σ w 1 ζ | | ζ | | Σ w 1 2 2 = ϵ | | z Y L | | Σ w 1 2 2 + 2 ( z Y L ) T Σ w 1 ζ | | ζ | | Σ w 1 2 2 ϵ 2 .
Inspired by these equalities, we divide S into S = ( S S ˜ ( a ¯ ) ) S ˜ ( a ¯ ) with
S ˜ ( a ¯ ) = { z S : 2 ( z Y L ) T Σ w 1 ζ | | ζ | | Σ w 1 2 2 ϵ 0 } .
We apply the following computation to the most right-hand side of (20)
S e x p | | z Y L ζ | | Σ w 1 2 2 d z e ϵ S S ˜ ( a ¯ ) e x p | | z Y L | | Σ w 1 2 2 d z + S ˜ ( a ¯ ) e x p | | z Y L | | Σ w 1 2 2 × 2 ( z Y L ) T Σ w 1 ζ | | ζ | | Σ w 1 2 2 d z e ϵ S e x p | | z Y L | | Σ w 1 2 2 d z + R p L e x p | | z Y L | | Σ w 1 2 2 1 S ˜ ( a ¯ ) d z .
From (16) we obtain, for any ϵ > 0 , S B ( R q L ) ; thus, we have
P ( Y L + W L S ) e ϵ P ( Y L + W L S ) + P U ϵ | | ζ | | Σ w 1 | | ζ | | Σ w 1 2 .
and U N 1 ( 0 , 1 ) . Then, the mechanism is ( ϵ , δ ) differential private if Q ( ϵ / | | ζ | | Σ w 1 | | ζ | | Σ w 1 / 2 ) δ , i.e.,
| | ζ | | Σ w 1 1 R ( ϵ , δ ) ,
for any ( x 0 , x 0 ) A d j 2 c . The inequality (25) holds if (17) is satisfied because
| | ζ | | Σ w 1 = | | O L ( x 0 x 0 ) | | Σ w 1 c λ m a x 1 / 2 ( O Σ , L ) .
From (18), λ m a x 1 / 2 ( O Σ w , L ) is the 2-induced matrix norm of Σ w 1 / 2 O L , denoted by | | Σ w 1 / 2 O L | | 2 . This can be bounded as
λ m a x 1 / 2 ( O Σ w , L ) = | | Σ w 1 / 2 O L | | | | Σ w | | 2 | | O L | | 2 = λ m i n 1 / 2 ( Σ w ) λ m a x 1 / 2 ( O ( L + 1 ) q , L ) .
and consequently,
λ m a x 1 / 2 ( O Σ w , L ) λ m i n 1 / 2 ( Σ w ) λ m a x 1 / 2 ( O ( L + 1 ) q , L ) .
Therefore, for any given c , ϵ > 0 and 0 < δ < 1 / 2 , one can make the Gaussian mechanism ( ϵ , δ ) -differentially private if one makes the minimum eigenvalue of the covariance matrix Σ w sufficiently large such that
λ m i n 1 / 2 ( Σ w ) c λ m a x 1 / 2 ( O ( L + 1 ) q , L ) R ( ϵ , δ ) .
Then, together with (28) and (29), we can obtain (17). □
Remark 1.
In the special case with an independent and identically distributed(i.i.d.) Gaussian noise, where Σ = σ 2 I ( L + 1 ) q , σ > 0 , (29) can be written as
σ c λ m a x 1 / 2 ( O ( L + 1 ) q , L ) R ( ϵ , δ ) ,
to make the mechanism ( ϵ , δ ) -differentially private.

4.2. Connection with Observability

It is evident that if the system (8) is observable, then the initial value x 0 can be determined from the output. This suggests an inherent relationship between the initial value exposure risk and the observability properties of linear nabla discrete fractional-order systems. Consequently, this subsection aims to further investigate the connection between observability and the protection of initial-value data via differential privacy.
First, we introduce the following definitions.
Definition 8.
The initial state x 0 of system (8) is said to be observable in a finite time L Z + if it can be uniquely determined from the observation of output sequence Y L . Otherwise, the system of (8) is said to be unobservable.
Definition 9.
The initial values x 0 are unobservable, if, for any pair of c-adjacent initial state x 0 , x 0 R n , the following equation
P [ H ( x 0 ) S ] = P [ H ( x 0 ) S ] ,
holds for all S r a n g e ( H ) .
Next, we give some rank conditions of observability for linear nabla fractional-order systems.
Theorem 2.
The following statements are equivalent.
1 
The initial states x 0 of system (8) is observable in a finite time L Z + .
2 
The observability matrix O L has rank n.
3 
The observability Gramian matrix R O ( k ) = s = 1 k F A , α T ( s ) C T C F A , α ( s ) has rank n.
Proof. 
(1)⇔(2)
(Sufficiency). According to Equation (9), we have
y ( k ) = C x ( k ) + D u ( k ) = C [ F A , α ( k ) x 0 + s = 1 k F A , α ( k s + 1 ) B u ( s ) ] + D u ( k ) .
Then, we have
C F A , α ( k ) x 0 = y ( k ) C s = 1 k F A , α ( k s + 1 ) B u ( s ) D u ( k ) = y ˜ ( k ) .
According to (33), we have
C C F ( 1 ) C F ( L ) x 0 = y ˜ ( 0 ) y ˜ ( 1 ) y ˜ ( L ) .
Now, suppose r a n k ( O L ) = n . Then, given u ( k ) and a known y ( k ) , x 0 can be determined uniquely and x 0 R n . This means the system (8) is observable.
(Necessity). Suppose the system is observable, then x 0 can be uniquely determined from the system of Equation (34). The system yields a unique solution for x 0 , which means there exists at least n rows in O L that are linearly independent. Consequently, this implies that the rank of O L is n. The proof is complete.
(1)⇔(3)
(Sufficiency). Suppose the matrix R O ( k ) has rank n. Multiplying both sides of (33) by F A , α T ( s ) C T and taking summation over the discrete time interval, we have
s = 1 L F A , α T ( s ) C T C F A , α ( s ) x 0 = s = 1 L F A , α T ( s ) C T y ˜ ( s ) R O ( k ) x 0 = s = 1 L F A , α T ( s ) C T y ˜ ( s ) .
Since rank ( R O ( k ) ) = n , the matrix is invertible and
x 0 = R O ( k ) 1 s = 1 L F A , α T ( s ) C T y ˜ ( s ) .
Thus, the system (8) is observable.
(Necessity). The contraction method is used in the proof. Suppose rank( R O ( k ) ) < n , then there exists a non-zero vector ξ R n such that R O ( k ) ξ = 0 . Then, we obtain
0 = ξ T R O ( k ) ξ = s = 1 L ξ F A , α T ( s ) C T C F A , α ( s ) ξ = s = 1 L | | C F A , α ( s ) ξ | | 2 2 ,
which implies C F A , α ( s ) ξ = 0 for all L > 0 . Thus, x ( 0 ) = x 0 + ξ yields the same response for the system (8) as x ( 0 ) = x 0 , which contradicts the assumption that the system (8) is completely observable. Therefore, r a n k ( R O ( k ) ) = n . □
Then, we have the following theorems for differential privacy from the observability perspective.
Theorem 3.
Then, the mechanism (14) induced by W L N ( L + 1 ) p ( μ , Σ w ) is ( ϵ , δ ) -differentially private at a finite time L Z + for A d j 2 c with ϵ > 0 and 0 δ 1 , if the covariance matrix Σ w 0 is chosen such that
λ m i n 1 / 2 ( Σ w ) c ( λ m a x 1 / 2 ( R O ( k ) ) ) R ( ϵ , δ ) .
Proof. 
It holds that
| O L ( x 0 x 0 ) | Σ w 1 λ m a x 1 / 2 ( Σ w 1 ) ( | O L ( x 0 x 0 ) | 2 ) c λ m a x 1 / 2 ( Σ w 1 ) ( λ m a x 1 / 2 ( R O ( k ) ) ) .
Since 1 / λ m a x ( Σ w 1 ) = λ m i n ( Σ ) , from (25), we can obtain (38). □
Remark 2.
From (38), it can be inferred that the Gaussian mechanism’s differential privacy is defined by the maximum eigenvalue of the observability Gramian matrix. This implies that the less input observable system can be made highly differentially private with only a minimal amount of noise added.
Theorem 4.
The initial values x 0 of system (8) are unobservable if and only if
r a n k ( O L ) < n .
Proof. 
(Sufficiency). We suppose (40) holds. Thus, there exists non-zero vector ρ k e r ( O L ) such that O L ρ = 0 , then we obtain O L x 0 = O L ( x 0 + ρ ) , which, together with (31), implies
P [ H ( x 0 ) S ] = P [ H ( x 0 + ρ ) S ] ,
for all H r a n g e ( H ) . Thus, this completes the proof by Definition 9.
(Necessity). To prove the necessity, we use the contraction method. With the unobservability of x 0 , we suppose that (40) does not hold, i.e., r a n k ( O L ) = n . Then, for all x 0 , x 0 , we have O L ( x 0 x 0 ) 0 , yielding O L x 0 O L x 0 . Thus, from the definition of H in (10), it is clear that there exists H r a n g e ( H ) such that
P [ H ( x 0 ) S ] P [ H ( x 0 ) S ] .
This contradicts the fact that the x 0 of the system (8) is unobservable. Therefore, it can be concluded that (40) holds, which completes the proof. □
Remark 3.
In terms of Definitions 7 and 9, a pair of ( x 0 , x 0 ) that satisfies (31) would yield (16) with ϵ = δ = 0 , implying perfect privacy. However, it cannot be conclusively said that the unobservability of x 0 directly infers its ( 0 , 0 ) -differential privacy. Conversely, it is evident that for ϵ = δ = 0 , the ( ϵ , δ ) -differential privacy of x 0 indeed implies its unobservability.

4.3. Minimization of the Noise

Building upon Equation (18), in this subsection, we determine the optimal Gaussian noise with the minimum energy that fulfills the sufficient condition (17) for differential privacy as established in Theorem 1. Specifically, we formulate an optimization problem wherein the covariance matrix of the Gaussian noise, denoted as Σ w 0 , is chosen as the decision variable.
Problem 1:
minimize Σ w 0 T r ( Σ w ) ,
subject to λ m a x ( O L T Σ w 1 O L ) 1 c 2 R 2 ( ϵ , δ ) .
Under certain specific assumptions and conditions, the solution to Problem 1 can be derived as follows:
Theorem 5.
Assume that O L is full rank, so the optimal solution to problem 1 is Σ w * = c 2 R 2 ( ϵ , δ ) O L O L T .
Proof. 
Denote Φ = c R ( ϵ , δ ) O L , so that Σ w * = Φ Φ T . From (44), we have I Φ T Σ w Φ 0 . According to the Schur complement, we have
Σ w Φ Φ T I 0 .
From (45), it implies that Σ w Σ w * . Consequently, for all Σ w that satisfy the constraint, T r ( Σ w ) T r ( Σ w * ) . Therefore, we can find that Σ w * is the optimal solution with the minimal trace for Problem 1. □

5. Numerical Example

In this section, a numerical example was provided to validate the effectiveness of the proposed algorithm.
Consider a feedback system in Figure 1, where the dynamic of plant P is given as
P : k α x p ( k ) = A p x p ( k ) + B p u p ( k ) , y p ( k ) = C p x p ( k ) ,
where x p R n p , u p R m p , and y p R q p are the state, input, and output, respectively, and A p , B p , C p are known matrices.
The dynamic of controller C is given as
C : k α x c ( k ) = A c x c ( k ) + B c e ( k ) , u p ( k ) = C c x c ( k ) ,
where x c R n c is the state of the controller, y r ( k ) denotes the reference signal, and e ( k ) = y r ( k ) y p ( k ) is the tracking error signal. Then, the composite system consisting of the plant (46) and controller (47) is
k α x ^ ( k ) = A ^ x ^ ( k ) + B ^ y r ( k ) , y p ( k ) = C ^ x ^ ( k ) , e ( k ) = y r ( k ) C ^ x ^ ( k ) ,
where x ^ ( k ) = x p ( k ) x c ( k ) , A ^ = A p B p C c B c C p A c , B ^ = 0 B c , C ^ = C p 0 .
We assume that the attackers with access to the output y p ( k ) try to estimate the initial state of plant x p ( 0 ) . The control objective is to minimize the tracking error e ( k ) and preserve the privacy of the initial state of plant x p ( 0 ) . In order to protect the sensitive information of x p ( 0 ) , we consider adding output noise w ( k ) to y p ( k ) such that y p ( k ) = C ^ x ^ ( k ) + w ( k ) , and the feedback system with output noise is shown in Figure 1.
The dynamics of system (46) and (47) are given by A p = 1.2 0.5 1 0 , B p = [ 0.3 , 0 ] T , C p = [ 0.2 , 0 ] , A c = 1 1 0 0.1 , B c = [ 0 , 1 ] T , C c = [ 1.5 , 0 ] , and α = 0.8 . We design the output noise to make the initial state of plant x p ( 0 ) differentially private for L = 30 , c = 2 , ϵ = 0.2 , and δ = 0.1 , and the distribution of W L is given by N ( 0 , Σ w ) . We then obtain R ( ϵ , δ ) = 0.2711 .
In what follows, we compare the following three cases (Table 1).
Figure 2 depicts the system output y p ( k ) and reference signal y r ( k ) for the three cases. The output error is the smallest in the noise-free case since the trajectory of y r ( k ) is fully utilized, increasing the potential for information leakage. In the case of i.i.d. noise, which has a uniform frequency distribution, there is an implication of an excessive out-of-band noise addition beyond the minimum required level. Consequently, this leads to a significant deterioration in tracking performance. In contrast, the optimized minimum noise strikes a balance, suppressing error fluctuations while preserving some tracking capability.

6. Conclusions

In this paper, we have studied a privacy-preserving problem of sensitive initial states for nabla discrete fractional-order dynamic systems. We developed a Gaussian mechanism by introducing output noise with ( ϵ , δ ) -differential privacy preserved. We analyzed the observability of nabla discrete fractional-order dynamic systems and derived its connection between differential initial-value privacy. Specifically, the differential privacy level of initial states can be characterized by the maximum eigenvalue of the observability Gramian matrix. Moreover, we formulated an optimization problem to derive the minimum energy random noise that satisfies the differential privacy requirement. In future works, the differential privacy approaches will be explored to develop a privacy-preserving controller for nabla fractional-order systems and the considered method will be utilized for nabla fractional-order multi-agent systems.

Author Contributions

Conceptualization, J.M.; methodology, J.M.; software, J.M.; validation, J.M.; formal analysis, J.M.; investigation, J.M.; resources, J.M.; data curation, J.M.; writing—original draft preparation, J.M.; writing—review and editing, J.M. and Z.P.; visualization, J.M.; supervision, J.H. and Z.P.; project administration, J.H.; funding acquisition, J.H. and Z.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 62203089, Grant 62103341, Grant 61473061, and Grant 12271083, in part by the Project funded by China Postdoctoral Science Foundation under Grant 2021M700695, in part by the Sichuan Science and Technology Program under Grant 2022NSFSC0890, and in part by the Guangdong Basic and Applied Basic Research Foundation under Grant 2022A1515110135.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Yamni, M.; Daoui, A.; El Ogri, O.; Karmouni, H.; Sayyouri, M.; Qjidaa, H.; Flusser, J. Fractional Charlier moments for image reconstruction and image watermarking. Signal Process. 2020, 171, 107509. [Google Scholar] [CrossRef]
  2. Wu, J.; Xu, Z.; Zhang, Y.; Su, C.Y.; Wang, Y. Modeling and tracking control of dielectric elastomer actuators based on fractional calculus. ISA Trans. 2023, 138, 687–695. [Google Scholar] [CrossRef] [PubMed]
  3. Traver, J.E.; Nuevo-Gallardo, C.; Tejado, I.; Fernández-Portales, J.; Ortega-Morán, J.F.; Pagador, J.B.; Vinagre, B.M. Cardiovascular circulatory system and left carotid model: A fractional approach to disease modeling. Fractal Fract. 2022, 6, 64. [Google Scholar] [CrossRef]
  4. Yousefpour, A.; Jahanshahi, H.; Munoz-Pacheco, J.M.; Bekiros, S.; Wei, Z. A fractional-order hyper-chaotic economic system with transient chaos. Chaos Solitons Fractals 2020, 130, 109400. [Google Scholar] [CrossRef]
  5. Yin, Y.; Guo, J.; Peng, G.; Yu, X.; Kong, Y. Fractal operators and fractional dynamics with 1/2 order in Biological Systems. Fractal Fract. 2022, 6, 378. [Google Scholar] [CrossRef]
  6. Chen, P.; Luo, Y.; Zheng, W.; Gao, Z.; Chen, Y. Fractional order active disturbance rejection control with the idea of cascaded fractional order integrator equivalence. ISA Trans. 2021, 114, 359–369. [Google Scholar] [CrossRef]
  7. Belkhatir, Z.; Laleg-Kirati, T.M. High-order sliding mode observer for fractional commensurate linear systems with unknown input. Automatica 2017, 82, 209–217. [Google Scholar] [CrossRef]
  8. Qi, F.; Qu, J.; Chai, Y.; Chen, L.; Lopes, A.M. Synchronization of incommensurate fractional-order chaotic systems based on linear feedback control. Fractal Fract. 2022, 6, 221. [Google Scholar] [CrossRef]
  9. Alessandretti, A.; Pequito, S.; Pappas, G.J.; Aguiar, A.P. Finite-dimensional control of linear discrete-time fractional-order systems. Automatica 2020, 115, 108512. [Google Scholar] [CrossRef]
  10. Stanisławski, R.; Latawiec, K.J. A modified Mikhailov stability criterion for a class of discrete-time noncommensurate fractional-order systems. Commun. Nonlinear Sci. Numer. Simul. 2021, 96, 105697. [Google Scholar] [CrossRef]
  11. Baleanu, D.; Wu, G.C.; Bai, Y.R.; Chen, F.L. Stability analysis of Caputo–like discrete fractional systems. Commun. Nonlinear Sci. Numer. Simul. 2017, 48, 520–530. [Google Scholar] [CrossRef]
  12. Wei, Y.; Su, N.; Zhao, L.; Cao, J. LMI based stability condition for delta fractional order system with sector approximation. Chaos Solitons Fractals 2023, 174, 113816. [Google Scholar] [CrossRef]
  13. Zhu, Z.; Lu, J.G. LMI-based robust stability analysis of discrete-time fractional-order systems with interval uncertainties. IEEE Trans. Circuits Syst. I-Regul. Pap. 2021, 68, 1671–1680. [Google Scholar] [CrossRef]
  14. Atici, F.M.; Eloe, P.W. Linear systems of fractional nabla difference equations. Rocky Mt. J. Math. 2011, 41, 353–370. [Google Scholar] [CrossRef]
  15. Wei, Y.; Chen, Y.; Liu, T.; Wang, Y. Lyapunov functions for nabla discrete fractional order systems. ISA Trans. 2019, 88, 82–90. [Google Scholar] [CrossRef]
  16. Wei, Y.; Zhao, L.; Wei, Y.; Cao, J. Lyapunov theorem for stability analysis of nonlinear nabla fractional order systems. Commun. Nonlinear Sci. Numer. Simul. 2023, 126, 107443. [Google Scholar] [CrossRef]
  17. Wei, Y.; Wei, Y.; Chen, Y.; Wang, Y. Mittag–Leffler stability of nabla discrete fractional-order dynamic systems. Nonlinear Dyn. 2020, 101, 407–417. [Google Scholar] [CrossRef]
  18. Hioual, A.; Ouannas, A.; Grassi, G.; Oussaeif, T.E. Nonlinear nabla variable-order fractional discrete systems: Asymptotic stability and application to neural networks. J. Comput. Appl. Math. 2023, 423, 114939. [Google Scholar] [CrossRef]
  19. Wei, Y.; Zhao, L.; Zhao, X.; Cao, J. Enhancing the Mathematical Theory of Nabla Tempered Fractional Calculus: Several Useful Equations. Fractal Fract. 2023, 7, 330. [Google Scholar] [CrossRef]
  20. Stanisławski, R.; Rydel, M.; Li, Z. A new reduced-order implementation of discrete-time fractional-order pid controller. IEEE Access 2022, 10, 17417–17429. [Google Scholar] [CrossRef]
  21. Zhao, M.; Li, H.L.; Zhang, L.; Hu, C.; Jiang, H. Quasi-projective synchronization of discrete-time fractional-order quaternion-valued neural networks. J. Frankl. Inst. 2023, 360, 3263–3279. [Google Scholar] [CrossRef]
  22. Gao, P.; Zhang, H.; Ye, R.; Stamova, I.; Cao, J. Quasi-uniform synchronization of fractional fuzzy discrete-time delayed neural networks via delayed feedback control design. Commun. Nonlinear Sci. Numer. Simul. 2023, 126, 107507. [Google Scholar] [CrossRef]
  23. Ma, J.; Hu, J.; Zhao, Y.; Ghosh, B.K. Leader-Following Consensus Control of Nabla Discrete Fractional Order Multi-Agent Systems. IFAC-PapersOnLine 2020, 53, 2897–2902. [Google Scholar] [CrossRef]
  24. Yuan, X.; Mo, L.; Yu, Y.; Ren, G. Containment control of fractional discrete-time multi-agent systems with nonconvex constraints. Appl. Math. Comput. 2021, 409, 126378. [Google Scholar] [CrossRef]
  25. Hong, X.; Zeng, Y.; Zhou, S.; Wei, Y. Nabla fractional distributed optimization algorithm with directed communication topology. In Proceedings of the 2023 6th International Symposium on Autonomous Systems (ISAS), Nanjing, China, 23–25 June 2023; pp. 1–6. [Google Scholar]
  26. Wang, J.; Cai, Z.; Yu, J. Achieving personalized k-anonymity-based content privacy for autonomous vehicles in CPS. IEEE Trans. Ind. Inform. 2019, 16, 4242–4251. [Google Scholar] [CrossRef]
  27. Ashkouti, F.; Khamforoosh, K.; Sheikhahmadi, A. DI-Mondrian: Distributed improved Mondrian for satisfaction of the L-diversity privacy model using Apache Spark. Inf. Sci. 2021, 546, 1–24. [Google Scholar] [CrossRef]
  28. Cramer, R.; Damgård, I.B.; Nielsen, J.B. Secure Multiparty Computation; Cambridge University Press: Cambridge, UK, 2015. [Google Scholar]
  29. Dwork, C. Differential privacy: A survey of results. In Proceedings of the International Conference on Theory and Applications of Models of Computation, Xi’an, China, 25–29 April 2008; Springer: Berlin/Heidelberg, Germany, 2008; pp. 1–19. [Google Scholar]
  30. Le Ny, J.; Pappas, G.J. Differentially private filtering. IEEE Trans. Automat. Contr. 2013, 59, 341–354. [Google Scholar] [CrossRef]
  31. Liu, X.K.; Zhang, J.F.; Wang, J. Differentially private consensus algorithm for continuous-time heterogeneous multi-agent systems. Automatica 2020, 122, 109283. [Google Scholar] [CrossRef]
  32. Yazdani, K.; Jones, A.; Leahy, K.; Hale, M. Differentially private LQ control. IEEE Trans. Automat. Contr. 2022, 68, 1061–1068. [Google Scholar] [CrossRef]
  33. Wang, Y.; Nedić, A. Tailoring gradient methods for differentially-private distributed optimization. IEEE Trans. Automat. Contr. 2023, 1–16. [Google Scholar] [CrossRef]
  34. Chen, B.; Leahy, K.; Jones, A.; Hale, M. Differential privacy for symbolic systems with application to Markov Chains. Automatica 2023, 152, 110908. [Google Scholar] [CrossRef]
  35. Kawano, Y.; Cao, M. Design of privacy-preserving dynamic controllers. IEEE Trans. Automat. Contr. 2020, 65, 3863–3878. [Google Scholar] [CrossRef]
  36. Wang, L.; Manchester, I.R.; Trumpf, J.; Shi, G. Differential initial-value privacy and observability of linear dynamical systems. Automatica 2023, 148, 110722. [Google Scholar] [CrossRef]
Figure 1. Feedback system with output noise.
Figure 1. Feedback system with output noise.
Fractalfract 08 00046 g001
Figure 2. Reference signal y r ( k ) and plant output y p ( k ) under three random noises.
Figure 2. Reference signal y r ( k ) and plant output y p ( k ) under three random noises.
Fractalfract 08 00046 g002
Table 1. Comparison of different random noises.
Table 1. Comparison of different random noises.
Noise-Freei.i.d. NoiseThe Minimum Noise
W L = 0 W l N ( 0 , σ i . i . d 2 I ) W L N ( 0 , Σ w * )
T r ( Σ w ) = 0 T r ( σ i . i . d 2 I ) = 34.545 T r ( Σ w * ) = 2.2784
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ma, J.; Hu, J.; Peng, Z. Privacy Preservation of Nabla Discrete Fractional-Order Dynamic Systems. Fractal Fract. 2024, 8, 46. https://doi.org/10.3390/fractalfract8010046

AMA Style

Ma J, Hu J, Peng Z. Privacy Preservation of Nabla Discrete Fractional-Order Dynamic Systems. Fractal and Fractional. 2024; 8(1):46. https://doi.org/10.3390/fractalfract8010046

Chicago/Turabian Style

Ma, Jiayue, Jiangping Hu, and Zhinan Peng. 2024. "Privacy Preservation of Nabla Discrete Fractional-Order Dynamic Systems" Fractal and Fractional 8, no. 1: 46. https://doi.org/10.3390/fractalfract8010046

APA Style

Ma, J., Hu, J., & Peng, Z. (2024). Privacy Preservation of Nabla Discrete Fractional-Order Dynamic Systems. Fractal and Fractional, 8(1), 46. https://doi.org/10.3390/fractalfract8010046

Article Metrics

Back to TopTop