Delay-Probability-Distribution-Dependent H ∞ FIR Filtering Design with Envelope Constraints

This paper studies the problem of H ∞ finite-impulse response (FIR) filtering design of time-delay system. The time-delay considered here is time-varying meanwhile with a certain stochastic characteristic, and the probability of delay distribution is assumed to be known. Furthermore, the requirement of pulse-shape is also considered in filter design. Employing the information about the size and probability distribution of delay, a delay-probability-distribution-dependent criterion is proposed for the filtering error system. Based on a Lyapunov-Krasovskii functional, a set of linear matrix inequalities (LMIs) are formulated to solve the problem. At last, a numerical example is used to demonstrate the effectiveness of the filter design approach proposed in the paper.


Introduction
In the studies about filtering problem, one most significant approach frequently applied in the past decades is Kalman filtering, the main idea of which is to minimize the variance of the estimation error assuming considered system dynamics to be exactly known and the external disturbances to be stationary Gaussian noises with known statistical properties [1,2].However, in many practical engineering applications, the statistical details about external noise are not available [3][4][5][6][7][8].In these cases, many approaches are introduced to improve systems' robustness, such as H ∞ , H 2 , and mixed H ∞ /H 2 filtering [2,[9][10][11][12][13].In this paper, the H ∞ filtering approach is utilized.
On the other hand, time-delays are frequently encountered in practical engineering systems, such as manufacturing systems, power systems, and networked control systems [14][15][16][17].Existence of delay makes the analysis and synthesis of systems a much more difficult task; meanwhile it is also the source of instability and poor performance in many cases [13,[18][19][20].The main approaches to solve delay problems can be classified into delay-dependent approach and delay-independent approach.It has been shown in [21,22] that the results obtained using delay-dependent approaches are generally less conservative than the delay-independent approaches ones [23].Acknowledging this fact, the delaydependent approach is applied in this paper.
In fact, the variation of delay may often stick to some probability distribution in spite of its varying and underivable property [24,25].Furthermore, in many real systems such as networked control systems, the time-varying delay may have some abrupt burst, leading to very large delay with a very small probability [26].In this sense, the discussion about time-delay should not only depend on its size but also on its probability distribution.In this paper, a new filter design approach and new stability criteria for the filtering error system taking the stochastic characteristic of time-varying delay into account is proposed.
While an H ∞ optimal filter can catch the frequencydomain property, the time-domain constraints such as envelope constraints or bounds on signals cannot be handled by this frequency-domain approach [27].Among various timedomain specifications, envelope constraints, which make requirement on the pulse-shape, have significant applications in many practical engineering systems, such as communication systems, radar, sonar systems, and signal processing systems [28][29][30][31].For instance, in deconvolution filtering and data channel equalization problems, it is extremely important to achieve a desired pulse-shape through designing an appropriate filter [27].
Therefore, aiming at incorporating both frequencydomain and time-domain constraints into the problem, we intend to design a filter satisfying the H ∞ performance and subject to envelope constraints in outputs.Meanwhile, timevarying delays with certain stochastic characteristics in the transmission channel are also taken into account.With the proposed filter design approach, a more general condition of time-varying delay problem can be solved.As in most situations, although detailed and exact information about delay cannot be achieved, the delay's probability distribution characteristics can be predicted or observed relatively easily.Once the probability information is gotten, the filter design approach can be developed.
In this paper, based on a Lyapunov-Krasovskii functional, we first present an H ∞ optimal solution to the design of a finite-impulse response (FIR) filter using information about the range of time-varying delay and its probability distribution.Then, the envelope constraints are taken into consideration.The resultant filter is called an H ∞ optimal Envelope-Constrained FIR (ECFIR) filter.We obtain the solution via solving an LMI optimization problem.At last, a numerical example is presented to illustrate the effectiveness of the proposed filtering design approach.

Problem Formulation and Preliminaries
Consider a filtering system shown in Figure 1, where Σ  represents a linear dynamic system with state-space realization given by where   () ∈ R   is the model state vector, () ∈ R   is the input signal, () ∈ R   is the source signal generated by the model, and   ,   ,   are known constant matrices with appropriate dimensions.Then the output () is transmitted through a channel with time-varying delay modeled by where   () ∈ R Obviously, it can be seen from the definition that  ∈ Ω 1 is equal to the occurrence of event () ∈ [0,  1 ) and  ∈ Ω 2 means that the event () ∈ [ 1 ,  2 ] occurs.Therefore, a stochastic variable () can be defined as where According to Assumptions 1 and 2, the system model described by (2) can be rewritten as At the receiving end, we are interested in designing a linear filter with state-realization as follows: where   () ∈ R   is the filter state vector, ŝ(), is the estimated signal of source signal () and   ,   ,   ,   have the following form: The transfer function of the filter is given by where (0), (1), . .., and (  ) are parameters to be determined.Define the filtering error as () = ()−ŝ().Then, via augmenting the models Σ  and Σ  , the filtering error system is given as follows: where Before giving the main results, we need following definitions at first.Definition 4. For a given function (()), its stochastic difference operator is defined as Definition 5 (see [32]).The filtering error system in ( 10) is said to be stochastically stable if for any initial condition   (0) and zero exogenous noise   (), there exists a positive definite  independent of   (0), such that the following condition is satisfied: Definition 6. System ( 10) is said to be stochastically stable with an H ∞ norm bound , if the following conditions hold.
(2) For all nonzero   () ∈  2 [0, ∞) and under zero initial conditions, the following inequality holds: Now, with the definitions above, we present the objective of this paper.
Given the filtering system shown in Figure 1, we are interested in designing a filter in the form of ( 7)- (8) such that (a) the filtering error system (10) is asymptotically stable in the stochastic sense; (b) the filtering error system (10) possesses a minimized H ∞ performance level ; (c) a time-domain envelope constraint is imposed on the output signal ŝ() as follows: where () and () are the lower and upper bounds of the time-domain mask, respectively.

Main Results
In this section, based on the Lyapunov-Krasovskii stability theorem, a delay-probability-distribution-dependent approach is proposed to solve the H ∞ FIR filter design problem subject to envelope constraints described in (15).First, a stability criterion for the filtering error system described in (10) is proposed.Then the envelope constraints are taken into consideration.An H ∞ optimal ECFIR filter design approach is given at last.
Proof.First, define a Lyapunov-Krasovskii functional as follows: where and  =   > 0,  1 =   1 > 0,  2 =   2 > 0,  1 =   1 > 0, and  2 =   2 > 0 are Lyapunov matrices to be determined.Then using the stochastic difference operator defined in (12), we obtain Using the Jensen inequality [33], the following expressions are obtained: Thus, we have Thus, we obtain where By Schur complement, it can be concluded from (17) that Υ < 0. By similar lines as in [32], the stochastic stability can be guaranteed if condition (17) holds.Then, define the performance index as follows:

Mathematical Problems in Engineering
Considering the fact that () ≥ 0, under the zero initial condition, we have Thus,  < 0 is equal to where Through applying Schur complement, it is shown that (Θ + Υ) < 0 can be guaranteed by condition (17).That is to say, once ( 17) is satisfied, the H ∞ performance can be guaranteed to be less than .Thus, the proof is completed.At this point, the second desired property of the system will be considered, which is the envelope constraints demand.First, some notations are introduced [34]: . . .

𝑓 (𝑛
where  is an  × (  + 1) matrix,  =  +   + 1, is a given signal, and Are, respectively, the upper and lower bounds.Therefore, the constraint of ( 15) is equal to where diag(•) denotes a conversion from a vertical vector to a diagonal matrix.
Based on Theorem 7 and (33), we can establish another theorem to determine the filter that satisfies the envelope constraint meanwhile possessing optimal H ∞ performance.Theorem 8.An H ∞ optimal filter of the form (7)-( 8) satisfying envelope constraint in (15) can be obtained by solving the following LMI optimization problem: where Ξ is defined in (17).

An Illustrative Example
In this section, an example is given to support the filter design method proposed in the paper.Consider a filtering system as shown in Figure 1.The parameters for Σ  are given by The parameters for the delay channel Σ  are given by Using Theorem 8, the H ∞ optimal filter is obtained via using the LMI toolbox of MATLAB with   chosen to be 5.The The expected envelope constraints and () (the output of Σ  ) corresponding to a particular case where input signal  () is chosen to be unit impulse signal are shown in Figure 2. The transmitted signal () through Σ  which is generated with no noise added is also given in the figure.The filter output ŝ() and filtering error () are given in Figures 3 and  4, respectively.
Furthermore, to illustrate the performance of the designed filter, we add the disturbance signal V() chosen as white noise with mean of zero and variance of 1 × 10 −3 into the system.The resultant filter output and filtering error are shown in Figures 5 and 6, respectively.It is shown that the designed filter is effective.

Conclusions
In this paper, we have solved the filtering design problem of time-delay system.The time-delay considered here is timevarying meanwhile with a certain stochastic characteristic, and the probability of delay distribution is assumed to be known.Furthermore, the envelope constraints are also considered in the process of filtering design.The delaydistribution-dependent criterion is formed for the filtering error system, employing the information about not only the size of delay but also its probability distribution.A set of linear matrix inequalities (LMIs) are formulated to solve the problem.Through solving the LMI optimization problem, the H ∞ performance is minimized and pulse-shape demand imposed by envelope constraints is satisfied.Finally, an illustrative example is presented to demonstrate the effectiveness of the filtering design approach.For future research directions, extending the filter design approach proposed in this paper to networked control systems and distributed systems is an interesting issue.Besides, more general filter    design approaches considering delays in different forms with different characteristics also deserve further investigation.

Figure 3 :
Figure 3: Output of the filter without disturbance.

Figure 5 :
Figure 5: Output of the filter with disturbance.
,   ,   ,   ,   ,   are all known constant system matrices with appropriate dimensions.As is shown in (2), the source signal () suffers from influence of time-varying delay () and disturbance from the environment represented by V().The output of transmission channel is (), which is also the input signal of the filter.We are going to use the corrupted signal () to reconstruct original source signal.