New relaxed stability and stabilization conditions for both discrete and differential linear repetitive processes

The paper develops new results on stability analysis and stabilization of linear repetitive processes. Repetitive processes are a distinct subclass of two-dimensional (2D) systems, whose origins are in the modeling for control of mining and metal rolling operations. The reported systems theory for them has been applied in other areas such iterative learning control, where, uniquely among 2D systems based designs, experimental validation results have been reported. This paper uses a version of the Kalman–Yakubovich–Popov Lemma to develop new less conservative conditions for stability in terms of linear matrix inequalities, with an extension to control law design. Differential and discrete dynamics are analysed in an unified manner, and supporting numerical examples are given.

series of sweeps of a processing tool, which is of finite duration. In the literature, each sweep is often termed a pass and the finite duration is known as the pass length.
An industrial example that highlights the characteristics of repetitive processes is metal rolling operations see, e.g., Rogers et al. (2015), where, in essence, deformation of the workpiece takes place by passing the workpiece between two rolls multiple times, i.e, passes in repetitive process terminology. The objective is to reduce the thickness to a predefined value, where each pass produces a pass profile y k (t), 0 ≤ t ≤ α. In this notation y is the vector or scalar-valued variable under consideration, the non-negative integer k is the pass number and α < ∞ denotes the pass length. Moreover, the pass profile on any pass explicitly contributes to the profile produced on the next pass and so on.
In general, the control problem for repetitive processes is that the sequence of pass profiles {y k } k can contain oscillations that increase in amplitude from pass-to-pass, i.e., in k. Moreover, this undesirable feature cannot be removed by standard control action, e.g., joining successive passes end-to-end to convert the dynamics to those of a standard system, see Rogers et al. (2007) where the deficiencies of this approach are explained. Instead, control design has to be based on regulating the 2-D dynamics that arise from information propagation along the passes and from pass-to-pass.
A stability theory for these processes has been developed (Rogers et al. 2007) that is of the bounded-input bounded-output form based on the pass profile. In particular, a bounded initial pass profile is required to produce a bounded sequence, {y k } k , of pass profiles, where boundedness is defined in terms of the norm on the underlying function space. This theory is based on a general model that includes all linear constant pass length examples as special cases. It allows specific treatment of the boundary conditions at the start of each pass, which must have a very particular and restrictive form for analysis based on the alternative above to produce correct results.
Aside from their use in the analysis and control of physical examples, the repetitive process setting has been used as a basis for solving problems in other areas of control theory, such as iterative learning control (ILC) law design, see, e.g., Rogers et al. (2015) and Paszke et al. (2016)) and iterative algorithms for solving nonlinear dynamic optimal control problems based on the maximum principle (Roberts 2002). In this latter case, a more general form of pass initial conditions is required and the alternative setting discussed above is not applicable. Moreover, in the ILC application designs based on repetitive process theory have been followed through to experimental validation, as in Paszke et al. (2016).
For linear dynamics, the stability theory has been developed the stage where the conditions for this property can be tested using well known tests from standard linear systems theory. However, these tests do not form a general basis for control law design and this has led to the development of Linear Matrix Inequality (LMI) based tests that extend naturally to control law design, where this approach is used in Paszke et al. (2016) with supporting experimental validation. These conditions are sufficient but not necessary and hence there is a level of conservativeness associated with their use.
The main objective of this paper is to develop new stability conditions in terms of LMIs by applying a frequency-discretization technique for differential and discrete linear repetitive processes, i.e., processes where the dynamics along a pass are governed, respectively, by a linear differential or difference equation. The outcome is relaxed stability conditions in terms of LMIs, which can be easily solved via standard numerical software. The main idea supporting the new result in this paper is a significant reduction of the conservatism of the resulting stability condition by dividing the entire frequency domain into sub-intervals and applying the Generalized Kalman-Yakubovich-Popov (GKYP) Lemma to each of them. As a result, some matrix variables are not required to be fixed over the entire frequency range. Also, this property enables control performance objectives to be imposed over specified finite frequency ranges. Hence a reduction in conservatism is possible by avoiding overdesign using the entire frequency range. Additionally, using the Elimination Lemma, standard analysis and design criteria are extended by introducing slack matrix variables. Hence, such conditions are more suitable for stabilizing controllers as products between matrix variables and those describing the process state-space model matrices are eliminated.
Preliminary results on this approach for discrete Li et al. 2015) and differential (Boski et al. 2018;Wang et al. 2017) repetitive processes have been reported. In this paper an approach is developed that enables differential and discrete dynamics to be analysed in unified manner by applying the GKYP and Elimination (Projection) Lemmas. Moreover, the new approach allows for a reduction in conservatism by introducing additional decision variables into final stability conditions in LMI form. Next, the stability conditions are extended to allow control law design and more general cases where the process matrices contain uncertain parameters. Finally, the effectiveness and advantages of the new conditions are demonstrated by two examples.
Throughout this paper the null and identity matrices, respectively, with compatible dimensions are denoted by 0 and I . For a square matrix M, sym{M} denotes M + M T and ρ(·) the spectral radius of its matrix arguments. Furthermore, M 0 (M ≺ 0) means that the symmetric matrix M is positive definite (negative definite). Also, two specific regions of the complex plane are defined: C hp = {λ ∈ C| Re(λ) < 0} (open left half plane) and C uc = {λ ∈ C| |λ| < 1} (an open disc centered at the origin and of unit radius). Finally, ⊗ denotes the Kronecker product, ( ) block entries in symmetric matrices and the superscript * denotes the complex conjugate transpose operation (for both scalars and matrices).
The new results in this paper are formulated in terms of LMIs and hence the following lemmas are useful in transforming non-LMI formulations into a problem subject to LMI constraints, where the first is the Elimination Lemma and the second is a version of the GKYP lemma.
Suppose that Λ(Φ, Ψ ) is a set of complex numbers that represents a certain class of curves in the complex plane defined as where Φ and Ψ are given matrices of dimension 2 × 2. As detailed in Iwasaki and Hara (2005), the imaginary axis and the unit circle are the particular examples of (3) with Ψ = 0 and Φ as and LF, MF and HF denote, respectively, the low, middle and high frequency ranges Table 2 The values of Ψ for discrete dynamics LF MF HF for the imaginary axis, respectively, the unit circle. Furthermore, by appropriately selecting Ψ , the frequency variable λ ∈ Λ(Φ, Ψ ) can be restricted to a certain range of frequencies (finite or semi-finite) and the following result can be established. (i) The frequency domain inequality holds ∀ω ∈ Ω (differential dynamics and λ = jω) or ∀θ ∈ Θ (discrete dynamics and λ = e jθ ) where Ω and Θ are, respectively the frequency ranges given in Tables 1 and 2. (ii) There exist matrices Q 0 and a symmetric matrix P 1 such that Appropriate choices for the matrix in Lemma 2 allows the analysis of the various system properties. Included in these is the bounded realness property, which is related to the stability and finite gain of a linear system. Hence, the matrix is fixed as where P 2 0, = diag{1, −γ 2 } and γ is a given scalar satisfying 0 < γ ≤ 1 and allows the specification of a process gain bound. Also, Lemma 2 provides a dual version of the GKYP lemma, which is exploited in the analysis that follows.

Stability of linear repetitive processes
The linear repetitive processes considered in this paper are of the following form with η representing the differential operator (in case of differential repetitive processes) or one step forward shift operator (in case of discrete repetitive processes). When a differential repetitive process is considered then t is the continuous variable such that t ∈ [0, α], where α is a strictly positive real number, while in case of discrete repetitive process t is a discrete variable such that t ∈ [0, α − 1] and the domain of α is limited to positive integers. Also, x k (t) ∈ R n , u k (t) ∈ R m and y k (t) ∈ R p represent the process state, input and output vectors at time instant t on pass k. The boundary conditions for these processes are the state initial vector on each pass and the initial pass profile. In the analysis of this paper, no loss of generality arises from assuming (for differential processes and with an equivalent statement for discrete processes) that x k+1 (0) = 0, ∀k ≥ 0, and y 0 (t) = f (t), 0 ≤ t ≤ α, where the entries in the vector f (t) ∈ R p are known functions of t.
In repetitive processes described by (8) and the assumed boundary conditions, the influence of the previous pass profile on the dynamics of the next pass arise from the term B 0 y k (t) for the state dynamics and D 0 y k (t) for the pass profile. The stability theory can be applied over the finite and fixed pass length or for all possible pass lengths, where this last case can be analyzed mathematically by considering the case when α → ∞. The first property is termed asymptotic stability and the latter stability along the pass. Moreover, stability along the pass requires the boundedness property to hold for all possible values of the pass length. In contrast, asymptotic stability requires this property for the pass length of the process and is, therefore, a necessary condition for stability along the pass.
Asymptotic stability guarantees that the sequence of pass profiles {y k } k converges as k → ∞ to dynamics described by a standard linear systems state-space model but not that this system is stable. The reason for this fact is the finite pass length, over which duration even an unstable system can only produce a bounded output. Stability along the pass removes this difficulty but there are cases, such as the optimal control application (Roberts 2002) where asymptotic stability suffices (or is all that can be achieved).
The structure of the boundary conditions is area where critical differences exist between repetitive processes and other classes of 2D systems dynamics. In some applications for repetitive processes the initial conditions at the start of each pass are explicit functions of points along the previous pass profile and the structure of these initial conditions play a critical role in stability analysis (Rogers et al. 2007). This form of boundary conditions cannot arise in Roesser and Fornasini Marchesini state-space models for quarter plane causal 2D linear systems and hence no exchange of, say, stability tests is possible (a detailed treatment of this issue is given in Rogers et al. (2007), which also cites the references for these models). In cases where the initial state vector on each pass does not depend on the previous pass profile, the dynamics of a linear repetitive process can be written as a Roesser model. Then, the stability tests and procedures available for the design of stabilizing control laws for the Roesser model can be applied. In particular, Bachelier et al. (2019) interprets the stability along the pass property of a linear repetitive process as the structural stability of the equivalent Roesser model and develops weakly conservative stability and stabilization results based on the polynomial solution parameter dependent LMIs. However, these results are for stabilizability and do not enable design to meet performance specifications over finite frequency ranges.
Suppose that the matrix Φ and the region of the complex plane C A are chosen as in Table 3.
Then the lemma given next establishes conditions for stability along the pass of linear repetitive processes. Table 3 Characterization of Φ and C A Differential process case Discrete process case

Lemma 3 (Rogers et al. 2007) A linear repetitive process described by (8) is stable along the pass if and only if
The first two conditions in Lemma 3 pose no computational difficulties but the third requires computations for the entire frequency range, i.e., for all points on the imaginary axis of the complex plane for differential processes and the unit circle for discrete processes. Also, the first condition is the necessary and sufficient condition for asymptotic stability. Moreover, the second condition, which governs the dynamics along a pass, is also only a necessary condition for stability along the pass and this property also requires the third condition, i.e., frequency attenuation of the previous pass profile dynamics over the complete frequency range and therefore as k → ∞ where (in the single-input single-output case for ease of presentation) In the remainder of this section the results of Lemmas 1 and 2 are used to develop new LMI based conditions for stability along the pass. In particular, based on the preliminary analysis in  and Paszke and Bachelier (2013), the following result can be established.
Theorem 1 Let γ be a positive scalar satisfying 0 < γ ≤ 1 and Φ chosen as in Table 3. Then a linear repetitive process described by (8) is stable along the pass if there exist compatibly dimensioned P 1 0, P 2 0 such that the LMI Proof Assume that (9) is feasible for some P 1 0, P 2 0. Applying directly the result of the KYP lemma (see Theorems 1 and 2 in Rantzer 1996) gives where G(λ) is defined by iii) in Lemma 3. Moreover,the inequality (10) can be rewritten as where the existence of P 2 0 directly implies that ρ(G(λ)) < 1, ∀λ ∈ Λ(Φ, 0), i.e., feasibility of (9) guarantees that condition iii) of Lemma 3 holds. Furthermore, since P 1 0 and P 2 0, then feasibility of (9) yields Equivalently, the conditions i) and ii) of Lemma 3, respectively, hold and the proof is complete.
Remark 1 In contrast to other classes of 2-D linear systems, the conditions on the matrices D 0 and A in Lemma 3 have physical meanings. The former requires that the initial pass profile sequence {y k (0)} k does not become unbounded with k and the latter that the contribution of the current pass state vector to the pass profile vector is stable. Hence for a given numerical example they should be tested before proceeding to the LMI of (9).

New LMI-based characterization of stability along the pass
Alternative LMI characterizations of the last two results can be obtained by introducing auxiliary slack matrix variables. These in turn, can be used to overcome non-convexity arising in control law design or reduce the level of conservativeness in many control problems for repetitive processes, such as robust control. To proceed, partition the matrix Φ (see (4)) into scalar elements as and introduce the following matrices of compatible dimensions where F 1 , F 2 and F 3 are additional compatibly dimensioned slack matrix variables. Then the theorem established next gives new LMI-based conditions for stability along the pass of linear repetitive processes.
Theorem 2 Let γ be a positive scalar satisfying 0 < γ ≤ 1 and Φ chosen as in Table 3. Then a linear repetitive process described by (8) is stable along the pass if there exist compatibly dimensioned matrices P 1 0, Proof Assume that the LMI defined in (12) has a feasible solution. Then this LMI can be rewritten as (1) with Then application of Lemma 1 enables (12) to be equivalently rewritten as (2). Also since Λ ⊥ = 0, the LMI in (12) reduces to the second inequality in (2), i.e., where by construction the matrix ⊥ is Furthermore, after some routine matrix manipulations the inequality (14) can be rewritten as where and by Lemma 1, feasibility of (15) implies that the inequality Finally, this last inequality is equivalent to (9) and by Theorem 1 stability along the pass is ensured.

Remark 2
The design conditions of Theorem 2 are LMIs that can be easily and effectively solved via numerical software. In addition, the scalar parameter γ can be minimized to compute the minimal process gain.

Remark 3
The LMI variables in (12) could be selected real symmetric instead of complex Hermitian without introducing conservatism, see Pipeleers and Vandenberghe (2011) for further details.

Remark 4
Revisiting the conditions in Theorem 2, it is found that Lemma 1 can again be invoked to reduce the number of slack matrix variables without adding conservatism to the solution. Specifically, the slack matrix variables can be constrained as W 2 = β 1 W 1 where β 1 is an arbitrary positive scalar. In particular, in (13) set Then choose Λ T ⊥ = [β 1 I − I 0] and check if the first inequality in (2) holds. In particular, under the above choice of Λ and W , the first inequality in (2) becomes Next, assign and by Lemma 1 the inequality (16) is also equivalent to Finally, the above inequality can be written as −2β 1 P 1 ≺ 0 for differential process case and holds for any β 1 > 0 and P 1 0 or as (β 2 1 − 1)P 1 ≺ 0 for discrete process case and holds for any |β 1 | < 1 and P 1 0.

State feedback based control of linear repetitive processes
In this section, the stability results of the previous section are extended to the case of controlled processes. The 2D system's structure of these processes requires that any control law has structure of feedback action on the current pass (stabilization and/or assignment of the state matrix eigenvalues to satisfy the along the pass specifications) and feedforward action from the previous pass (convergence of the pass-to-pass dynamics). One such control law has the form where K 1 and K 2 are matrices to be found. Application of this control law results in the controlled process state-space model where the meaning of η is the same as in (8). Furthermore, by defining it follows immediately that the results of the previous section can be directly used to develop a condition suitable for computing control law gains K 1 and K 2 . In particular, the following result is obtained by applying Theorem 2 to (18).
Theorem 3 Let γ be a positive scalar satisfying 0 < γ ≤ 1 and Φ chosen as in Table 3. Suppose also that the control law (17) is applied to a repetitive process described by (8).
Then the resulting controlled process is stable along the pass if there exist matrices P 1 0, is feasible. In addition, if the above LMIs are feasible, the corresponding matrices in the control law (17) are given by Proof Assume that the LMI (20) holds. Also, it follows immediately that the feasibility of (20) implies that W 1 is non-singular. Also setting Y = K , W T 1 leads to the interpretation of the condition of Theorem 2 for (18) and hence stability along the pass for the controlled process. This completes the proof.
Remark 5 This last result requires the following constraints on the multipliers W 1 = W 2 and W 3 = [0 I ]W 1 and introduces conservativeness but still some additional freedom is introduced by the slack matrix variables F 1 , F 2 , F 3 . A possible way to reduce conservatism is to apply the constraints where β 2 and β 3 are some scalars. Hence, line searches on these scalars are required to make the LMI less conservative (in terms of performance provided by γ ).

Less conservative conditions for stability along the pass
The results of Theorems 1, 2 and 3 can introduce a significant degree of conservativeness since it has been obtained by keeping P 2 (and some other matrix variables) constant and independent of frequency. This means that P 2 must work for entire imaginary axis or the unit circle, i.e., ∀λ ∈ Λ(Φ, 0). One possible way to overcome this problem is to extend the results developed in previous sections of this paper to allow piecewise constant matrix variables defined over union of segments on the entire imaginary axis or the unit circle in the complex plane.
These segments must be contiguous to ensure continuity of λ over Λ(Φ, 0) and this can be achieved by dividing the entire frequency range into a number of finite or semi-finite ranges. Also, specific choices of a pair of Φ define (semi-) finite frequency ranges as shown in Tables 1 and 2 and correspond to segments of the imaginary axis (differential processes) or unit circle (discrete processes). Moreover, using such frequency partitioning enables specifications to be imposed over frequency ranges relevant to a specific application of repetitive processes. Furthermore, to allow for more flexibility in practical design, frequency specifications in different ranges can be introduced without augmenting the plant with frequency dependent scalings or weights.
To establish the main result of this section, consider the partitioning where N ≥ 2. Then by different choices of Ψ (i.e. Ψ 1 ,Ψ 2 ,. . ., Ψ N ) the specific ranges of the complex curves are characterized. Moreover, both continuous-time (frequency variable λ on the imaginary axis) and discrete-time cases (λ on the unit circle) can be treated by a different choice of Φ-see (4).
In what follows, the partitioning (22) implies that the frequency variable λ is defined over unions of finite or semi-finite frequency ranges. In particular, for the continuous-time setting Table 4 The values of Ω i for differential dynamics Where 2 ≤ r ≤ N − 1 and ω cr = (ω r +ω r −1 ) 2 Table 5 The values of Θ i for discrete dynamics Where 2 ≤ r ≤ N − 1 and θ dr = θ r −θ r −1 2 , θ cr = θ r +θ r −1 2 and λ = jω, (22) imposes the following frequency partitioning where Ω 1 = {ω : {|ω| < ω 1 }}, Ω r = {ω : {ω r −1 ≤ |ω| ≤ ω r }}, Ω N = {ω : {|ω| > ω N }} and 2 ≤ r ≤ N − 1. In the discrete-time case, λ = e jθ and the following partitioning results where and 2 ≤ r ≤ N − 1. Ω i and Θ i are the frequency ranges given below in Tables 4 and 5. Each matrix Ω i or Θ i is of dimension 2×2 and consists of three different elements (scalars) which depend on the chosen frequency interval. In particular, the following partitioning of Ψ is enforced The values of τ , υ and ς for specified choices of frequency intervals are also shown in Tables 4  and 5. Also, the application of piecewise constant matrices P 1i , P 2i and Q i ∀i = 1, . . . , N , is possible and hence the LMI feasibility problem of the next theorem is a novel condition for stability along the pass.
Theorem 4 Let γ be a positive scalar satisfying 0 < γ ≤ 1 and Φ chosen as in Table 3. Also, suppose that the entire frequency range is arbitrarily divided into N different frequency intervals as given in (22)-(24) and the matrix Ψ chosen as in Tables 4 or 5. Then a linear repetitive process described by (8) is stable along the pass if there exist P 1i 0, P 2i 0, Proof Suppose that there exist matrices P 1i 0, P 2i 0 and Q i 0 for i = 1, . . . , N , such that the LMIs (25) hold for i = 1, . . . , N . Then by simple matrix manipulations on each LMI in (25) it follows that Finally, the proof is completed applying Lemma 2 to each frequency interval with the choice of piecewise constant matrices P 1i 0, P 2i 0, Q i 0, ∀i = 1, . . . , N .
By analogy with the derivation of Theorem 3, the above result can be directly extended to the case of controlled processes (18) where the control law has the structure given in (17). Therefore, with the notation used above, the following result extends Theorem 4 to the controlled process.
Theorem 5 Let γ be a positive scalar satisfying 0 < γ ≤ 1, Φ is as in Table 3 and Ψ chosen as in Tables 4 or 5. Suppose also that the control law (17) is applied to a repetitive process described by (8) and the entire frequency range is arbitrarily divided into N different frequency intervals as given in (22). Then the controlled process is stable along the pass if there exist matrices P 1i 0, P 2i 0, are feasible for i = 1, . . . , N . In addition, if these LMIs are feasible, the required control law matrices K 1 and K 2 of (17) can be calculated using (21).

Proof
The result follows from straightforward application of the same steps as in proof of Theorem 3 where the matrix variables P 1 , P 2 , F 1 , F 2 , F 3 are directly replaced by their frequency range dependent form P 1i , P 2i , F 1i , F 2i , F 3i and hence details are omitted.
Remark 6 These last two results enable different performance specifications to be imposed in different frequency ranges. In particular, it is reasonable from the practical point of view to impose different control performance by choosing different process gains in specified frequency ranges, e.g., to ensure that G(λ) meets the specification where γ i ∈ (0, 1] ∀i = 1, . . . , N .

Robustness
This section further extends the results obtained so far to processes with uncertain model data, i.e., uncertain parameters are present in the process matrices. In the following, we will consider two types of parameter uncertainties: norm-bounded and polytopic uncertainties.

Norm-bounded uncertainty case
It is assumed that the deviation of the process parameters from their nominal values is normbounded. Specifically, the uncertainty associated with the state-space model (8) is of the form where the matrices A, B, B 0 , C, D and D 0 represent the nominal process dynamics of (8) and are assumed to be time and pass number invariant. The matrices ΔA(t), ΔB(t), ΔB 0 (t), ΔC(t), ΔD(t) and ΔD 0 (t) denote time-varying uncertainties where E 1 , E 2 , H 1 , H 2 and H 3 are given real constant matrices with appropriate dimensions and δ(t) is an unknown and time-varying perturbation satisfying δ The proofs of the results in this section make use of the following result.
Lemma 4 (Petersen 1987) Given matrices X , Y , Z = Z T , δ(t) of compatible dimensions, then for all δ(t) satisfying δ T (t)δ(t) I if, and only if, there exists ε > 0 such that Introduce the notation Then, based on Theorem 4, the following results are obtained for stability along the pass in the presence of the parametric uncertainties in the state-space model matrices described by (30).
Theorem 6 Consider an uncertain linear repetitive process described by the version of (8) with uncertainty described by (30). Also, let γ be a positive scalar satisfying 0 < γ ≤ 1, Φ is as in Table 3 and Ψ chosen as in Tables 4 or 5. Also, suppose that the entire frequency range is arbitrarily divided into N different frequency intervals as given in (22)-(24). Then a linear repetitive process described by (8) and (30) is robustly stable along the pass if there exist compatibly dimensioned matrices P 1i 0, P 2i 0, Q i 0, W 1i , W 2i , W 3i , F 1i , F 2i , F 3i and a positive scalar 1 such that the following LMIs ⎡ are feasible for i = 1, . . . , N , where Proof Suppose that the LMIs (31) are feasible for i = 1, . . . , N . Then application of Schur's complement formula to (31) gives

and by Lemma 4 the last inequality is feasible if and only if
The last inequality is (25) applied to the uncertainty case. Moreover, the LMI of (31) can be obtained by employing the same steps used in the absence of uncertainty and hence a version of (4) applied to the plant with uncertainty is obtained. Finally, by the result of Theorem 3, feasibility of (31) ensures that a linear repetitive process of the form (8) and (30) is robustly stable along the pass.
The following result extends Theorem 5 to the case when norm bounded uncertainty is present. This is achieved by applying the result of Lemma 4 to Theorem 5 for uncertain processes.
Theorem 7 Consider an uncertain linear repetitive process described by the version of (8) with uncertainty of the form (30) present and with the control law (17) applied. Also, let γ be a positive scalar satisfying 0 < γ ≤ 1, Φ is as in Table 3 and Ψ chosen as in Tables 4 or 5. Suppose that the entire frequency range is arbitrarily divided into N different frequency intervals as given in (22)-(24). Then, the controlled process is robustly stable along the pass if there exist P 1i 0, P 2i 0, Q i 0, F 1i , F 2i , F 3i , Y and a positive scalar 2 such that the following LMIs are feasible for i = 1, . . . , N , where E 1 and H 1 are as in (31) and If these LMIs are feasible, the corresponding matrices in the control law (17) can be computed using (21).

Proof
The result follows from a straightforward application the same steps as in proof of Theorem 6 for a controlled process and hence the details are omitted.

Polytopic uncertainty case
An alternative way of dealing with uncertain processes is to assume that the repetitive process matrices are not exactly known but reside within a given polytope. Specifically, this means that all process matrices are actually dependent on a real parameter vector ξ with a polytopic dependency: where It means that the various matrices A j , B 0 j , B j , C j , D 0 j , D j define the vertices of the polytope P in which the actual process matrices A, B 0 , B, C, D 0 and D lie. As known, this clearly removes the requirement for keeping the same matrix variables over the entire parametric uncertainty domain and given frequency range. In particular, some matrix variables can be made parameter-dependent (linearly dependent on ξ ) due to their appearance free from products with the repetitive process matrices. Simply, we assume that matrix functions P 1i (ξ ), , which are to be determined, are expected to have the Theorem 8 Let γ be a positive scalar satisfying 0 < γ ≤ 1 and Φ chosen as in Table 3. Also, suppose that the entire frequency range is arbitrarily divided into N different frequency intervals as given in (22)-(24) and the matrix Ψ chosen as in Tables 4 or 5. Then an uncertain linear repetitive process with process model matrices described by (33) and (34) is robustly stable along the pass if there exist P 1i j 0, P 2i j 0, Q i j 0, F 1i j , F 2i j , F 3i j , W 1i , W 2i and W 3i , such that, the LMIs ⎡ The following result is a linearly parameter-dependent approach to controller design for processes with polytopic uncertainty.
Theorem 9 Let γ be a positive scalar satisfying 0 < γ ≤ 1, Φ is as in Table 3 and Ψ chosen as in Tables 4 or 5. Suppose also that the control law (17) is applied to an uncertain linear repetitive process (33) and the entire frequency range is arbitrarily divided into N different frequency intervals as given in (22). Then the controlled process is robustly stable along the pass if there exist matrices P 1i j 0, P 2i j 0, are feasible for i = 1, . . . , N and j = 1, . . . , M. In addition, if these LMIs are feasible, the required control law matrices K 1 and K 2 of (17) can be calculated using (21).
Proof The proof can be carried out by employing similar lines to that of the proofs of Theorem 5 and 8. Hence the details are omitted.

Simulation based case study
To justify the effectiveness and feasibility of the formulated results, two numerical examples are presented. The first example is formulated to show the reduced conservatism compared with Li et al. (2015) for stability along the pass of discrete repetitive processes. Next, the second example illustrates the results of controller design procedure for differential repetitive processes compared with Rogers et al. (2007) and its extended form in Maniarski et al. (2020).

Example 1
Let the process matrices in (8) be give as follows (as noted, the below example is adopted from Li et al. (2015)) Solving the LMI condition proposed in Theorem 2 by applying Matlab LMI Toolbox routines, gives that the LMI (12) is feasible and the minimum γ is 0.9758. By testing the stability condition in Li et al. (2015) (formulated in Theorem 1 of Li et al. 2015) over entire frequency range, it follows that the result in Li et al. (2015) does not support feasibility for γ ≤ 1 and hence fails to decide the stability along the pass of the considered process. Thus, the effectiveness of the proposed condition is clearly demonstrated since our result is less restrictive compared with Li et al. (2015) for entire frequency range. The stability along the pass for the process in Example 1 can be confirmed with Li et al. (2015) only in the cases when the entire frequency range is partitioned uniformly into 4 subintervals or partitioned non-uniformly into 2 subintervals ([0, 0.5] ∪ [0.5, π]). See also Fig. 1.

Example 2
As an illustrative example of the new results in this paper consider the linearized differential linear repetitive process model that represents the metal rolling process given in Bochniak Here λ 1 = 600[N /m] is the stiffness of the adjustment mechanism spring, λ 2 = 2000 [N/m] is the hardness of the metal strip and M = 100 [kg] denotes the lumped mass of the roll-gap adjusting mechanism. In this case, the nominal state-space model matrices in (8)  Suppose that the parameter M is uncertain and possibly time-varying, and its values are assumed to be within the interval [90, 112.5]. Then this case corresponds to the scenario when the matrices E 1 , E 2 , H 1 , H 2 and H 3 are The plots of ρ(G( jω)) given in Fig. 2 show that the nominal and uncertain repetitive processes (for the extremes of the uncertain element ranges, i.e. δ(t) = 1) are unstable along the pass. Clearly, the condition iii) of Lemma 3 is not satisfied for these processes and hence the LMIs in (12)  Hence ω 1 = 1.7[rad/s], ω 2 = 2.29[rad/s] and ω 3 = 3[rad/s]. However, there is no systematic way to determine what is the best choice of the frequency sub-intervals but our intention is to  Fig. 3 Plot of ρ(G C ( jω)) for the models: nominal and two with the extreme uncertainties keep resonance and antiresonance frequencies within separate sub-intervals since different performance (provided by γ i -see Remark 6) may be achieved in specified frequency ranges (sub-intervals). Solving the LMI condition obtained in Theorem 7 for γ = 1 and using (21) we obtain the following control law matrices K 1 = 1 , 0 3 · [−1.2527 1.4887], K 2 = 5 , 51.9732 The controlled uncertain process is stable along the pass as seen in Fig. 3 where the spectral radius of the controlled process is below 1 for the entire frequency range, i.e., ρ(G C ( jω)), Note that the existing results on control law design for uncertain differential linear repetitive process in Rogers et al. (2007) and its extended form Maniarski et al. (2020) are feasible for this example as demonstrated in Fig. 4. However, new design procedure provides better performance for low frequencies. In addition, the controlled dynamics were simulated over 20 passes where the initial pass profile is chosen as y 0 (t) = sin(2πt/2 , 0) + 0.5 sin(4πt/2 , 0) [mm] for 0 ≤ t ≤ 20. Also, it was hypothetically assumed that the sheet metal thickness should be finally close to zero, hence the quality of the control over each pass profile can be measured by the RMSE (Root Mean Squared Error -computed along each pass) and shown in Fig. 5. This also confirms robust stability along the pass of the controlled process and better performance obtained by the new stabilization design since the convergence rate is higher than the alternatives.

Conclusions and future works
This paper has developed new results on the stability and stabilization for linear repetitive processes. Both discrete and differential dynamics of a given repetitive process have been considered and novel conditions for stability along the pass and control law design have been obtained. The major benefit of these new results is that stability along the pass tests can be directly extended to the control law design algorithms for nominal and uncertain models. Also, the advantage over current results in this last aspect is the avoidance of product terms between the state-space model matrices and some Lyapunov/LMI decision matrices. This decoupling has been achieved through the use of slack matrix variables and therefore a reduction in the conservatism of the analysis and design tests is possible. Finally, the numerical examples to demonstrate the effectiveness of the new results have been given.
Future research should include a detailed investigation into the H ∞ robust performance in the presence of external bounded non-repetitive disturbances in both the state and pass profile vectors. Another area is design of the more complex controller structures, e.g. a dynamic controller, and their experimental validation. Moreover, it is worth noting that reducing the conservatism of conditions for stability along the pass recently gained increasing attention, and significant progress has been made based on interpreting the stability along the pass of a linear repetitive process as structural stability of 2D Roesser model. Therefore, it is also promising to develop robust control results by extending the methods of Bachelier et al. (2019). Wojciech Paszke received his M.Sc. and Ph.D. degrees, both in electrical engineering, from the Technical University of Zielona Gora, in 2000 and 2005, respectively. Between 2008 and 2010 he was affiliated to Eindhoven University of Technology, The Netherlands, where he has been a control systems expert on high precision positioning of electron microscope. Currently, he is affiliated to the Institute of Automation, Electronic and Electrical Engineering at the University of Zielona Gora, Poland. His research is focused on multidimensional (nD) systems, repetitive processes, iterative learning control schemes and convex optimization in robust control problems.
Eric Rogers is currently a Professor of control systems theory and design with the School of Electronics and Computer Science, University of Southampton, Southampton, U.K. His major research interests are in iterative learning control theory and applications in both engineering and healthcare, multidimensional systems theory and applications, long-range autonomy for autonomous underwater vehicles, and flow control. Dr. Rogers is also the Editor-in-Chief of the International Journal of Control and Multidimensional Systems and Signal Processing. He has an extensive record of advisory work for funding agencies and industries in U.K. and elsewhere.