Processing math: 74%
A journal of IEEE and CAA , publishes high-quality papers in English on original theoretical/experimental research and development in all areas of automation
X. Xu and  G.-R. Duan,  “High-order fully actuated system models for strict-feedback systems with increasing dimensions,” IEEE/CAA J. Autom. Sinica, vol. 11, no. 12, pp. 2451–2462, Dec. 2024. doi: 10.1109/JAS.2024.124599
Citation: X. Xu and  G.-R. Duan,  “High-order fully actuated system models for strict-feedback systems with increasing dimensions,” IEEE/CAA J. Autom. Sinica, vol. 11, no. 12, pp. 2451–2462, Dec. 2024. doi: 10.1109/JAS.2024.124599

High-Order Fully Actuated System Models for Strict-Feedback Systems W

ith Increasing Dimensions

doi: 10.1109/JAS.2024.124599
Funds:  This paper was partly supported by Shenzhen Key Laboratory of Control Theory and Intelligent Systems (ZDSYS20220330161800001) and the National Natural Science Foundation of China (NSFC) (62350055)
More Information
  • This paper mainly addresses control problems of strict-feedback systems (SFSs) with increasing dimensions. Compared with the commonly-considered SFSs where the subsystems have the same dimension, we aim to handle more complex cases, i.e., the subsystems in the considered SFSs are assumed to have increasing dimensions. By transforming the systems into high-order fully-actuated system (HOFAS) models, the stabilizing controllers can be directly given. Besides first-order SFSs, second-order and high-order SFSs are also considered.

     

  • STRICT-FEEDBACK systems (SFSs) refer to dynamic systems consisting of multiple subsystems in a triangular structure. Control of SFSs has been widely studied for many years with a quantity of interesting results obtained, to name a few, [1]–[5]. A typical SFS is given as follows:

    {˙x1=f1(x1)+G1(x1)x2˙x2=f2(x1,x2)+G2(x1,x2)x3˙xn=fn(x1,x2,,xn)+Gn(x1,x2,,xn)u (1)

    where $ x_i\in\mathbb{R}^{k},\; i = 1,2,\ldots,n, $ are the state vectors, $ u\in\mathbb{R}^k $ is the control input, $ f_i:\mathbb{R}^{ik}\rightarrow\mathbb{R}^k $ and $ G_i:\mathbb{R}^{ik}\rightarrow\mathbb{R}^{k\times k} $ are sufficiently smooth functions. One of the most commonly-used method for handling control problems of such a system is the well-known backstepping method.

    The backstepping technique, pioneered in the 1990s by Petar V. Kokotovick, Miroslav Krstic, and others, originally aimed at crafting stabilizing controls for SFSs [6], [7]. Its fundamental concept involves initiating the design process at the first subsystem already stabilized and then iteratively deriving new controllers to progressively stabilize subsequent subsystems until achieving control over the final external aspect. While this method assumes uniform dimensions for subsystems (i.e., $ x_i,\; i = 1,\ldots,n $), practical systems often feature nonuniform or differing dimensions. Examples include cascades of over-actuated [8] or under-actuated systems [9]. Considering SFSs with nonuniform dimensions enables the development of more general theories, encompassing uniform cases as special instances. However, designing controllers for SFSs with nonuniform dimensions via the traditional backstepping approach proves challenging, if not infeasible. Therefore, addressing the more general scenario-where subsystems possess nonuniform dimensions-requires innovative approaches.

    Fortunately, there are effective approaches available to address control problems in the context of high-order fully actuated systems (HOFAS). HOFAS, also referred to as fully actuated system (FAS) approaches, were originally introduced by G. Duan in a series of papers [10]–[12]. The HOFAS approach offers a distinct perspective for handling control problems by transforming a system into a HOFAS model, which allows for the direct derivation of an exponentially stabilizing controller. Duan has demonstrated that nearly all controllable systems, regardless of linearity or nonlinearity and time-invariance or time-varying nature, can be transformed into HOFASs. The high-order FAS theory has opened up new avenues in control theory and has attracted significant research attention from scholars who have obtained numerous results in this field [13]–[22].

    In [11], it is shown that first-order, second-order and high-order SFSs can all be equivalently transformed into HOFAS models. Although only SFSs with uniform dimension are considered in [11], the HOFAS approach shows its potential in handling control problems of various complicated SFSs. It is thus natural to apply HOFAS approaches to more complicated types of SFSs, one of which is the type of SFSs with increasing dimensions, which motivates this work.

    This paper investigates the control of SFSs with expanding dimensions, introducing HOFAS models tailored for such systems to facilitate the direct derivation of exponentially stabilizing controllers. To our best knowledge, this work is among the first few studies concerning about control of SFSs with increasing dimensions. The contributions of this study can be delineated into two main areas.

    Firstly, we present HOFAS models specifically designed for SFSs with increasing dimensions, building upon prior work wherein various SFSs were transformed into FAS models [11]. While previous research [11] predominantly addressed scenarios where subsystem dimensions remained constant, we extend these findings to encompass more intricate cases, namely SFSs with expanding dimensions. A pivotal technique employed involves the decomposition of variables based on these increasing dimensions. Additionally, whereas existing literature typically focuses on first-order SFSs, we broaden our scope to include second- and higher-order SFSs, enhancing the generality of our approach.

    Secondly, we tackle the control problem associated with SFSs exhibiting increasing dimensions. Leveraging the FAS models tailored for such systems, we derive stabilized controllers directly. A notable advantage of our controllers over those designed using backstepping methods lies in the resulting closed-loop systems, which exhibit linear time-invariance. This feature underscores the efficacy of HOFAS approaches in addressing control challenges posed by more complex SFSs.

    The rest of this work is organized as follows. In Section II, notations are introduced. In Sections III, HOFAS models are proposed for two types of first-order SFSs with increasing dimensions. The matrices $ G_i,\; i = 1,2,\ldots,n $ are first assumed to be constant and then assumed to be matrix functions. Second-and high-order SFSs with increasing dimensions are considered in Section IV. Conclusions are drawn in Section VI.

    Throughout this paper, the following notations are used. The symbol $ I_k $ represents $ k\times k $ identity matrix. $ 0 $ represents zero matrix with appropriate dimension. The symbol $ \triangleq $ means “denoted as” in this work. For two integers $ i,j, $ where $ j\geq i $, the notation $ i\sim j $ represents the series of integers $ i, i+1, \ldots, j $. For a function $ x\in\mathbb{R}^n $, $ x^{(i)} $ denotes its i-th order derivative. Moreover, the following symbol is used:

    x(ij)=(x(i)x(i+1)x(j)). (2)

    For three integers $ i,j,k, $ where $ j\geq i, $ and for $ x_{pk}\in\mathbb{R}^n, \; p = i\sim j $, the following symbols are used:

    xij,k=(xikxi+1,kxjk) (3)
    x(0l(p))pk|p=ij=(x(0l(i))ikx(0l(i+1))i+1,kx(0l(j))jk) (4)

    where $ l:\mathbb{Z}^+\rightarrow\mathbb{Z}^+ $ is an arbitrary function mapping from positive integers into positive integers. In addition, one has the following definition:

    x(0l(p))pk|p=ij,k=12=(x(0l(p))p1|p=ijx(0l(p))p2|p=ij). (5)

    In this subsection, we consider first-order SFSs with increasing dimensions. Consider the following SFS:

    {˙x1=f1(x1)+G1x2˙x2=f2(x1,x2)+G2x3˙xn=fn(x1,x2,xn)+Gnu, (6)

    where $ x_i\in\mathbb{R}^{k_i},\; i = 1,2,\ldots,n $, $ u\in\mathbb{R}^m $, $ G_i\in\mathbb{R}^{k_i\times k_{i+1}},\; i = 1, 2,\ldots, n-1, $ and $ G_n\in\mathbb{R}^{k_n\times m} $. In this work, it is assumed that $ k_1\leq k_2\leq \cdots\leq k_n\leq m. $ Hence, the dimensions of subsystems are increasing. Moreover, the following assumption is needed for constant matrices $ G_i,\; i = 1,2,\ldots,n. $

    Assumption 1: The matrices $ G_i,\; i = 1,2,\ldots,n, $ are all of full-row rank, i.e., $ {\rm rank}(G_i) = k_i,\; i = 1,2,\ldots,n $.

    Assumption 1 ensures that there exist nonsingular matrices $ Q_i, \; i = 1,2,\ldots,n $ and $ G_i, \; i = 1,2,\ldots,n $ such that

    GiQi=[ˉGi0]. (7)

    Based on this fact, let us define $ x_{i+1,1}\in\mathbb{R}^{k_i} $, $ x_{i+1,2}\in\mathbb{R}^{k_{i+1}-k_i} $, $ f_{i1}\in\mathbb{R}^{k_i} $ and $ f_{i2}\in\mathbb{R}^{k_{i+1}-k_i} $ by

    Q1ixi+1=(xi+1,1xi+1,2) (8)
    Q1ifi+1(x1,x2,,xi+1)=(fi+1,1(x1,x2i,1,x2i,2)fi+1,2(x1,x2i,1,x2i,2)),i=1,2,,n1. (9)

    For notation convenience, we also define that $ x_{12} = x_1. $ Meanwhile, define $ u_1\in\mathbb{R}^{k_n} $ and $ u_2\in\mathbb{R}^{m-k_n} $ that

    Q1nu=(u1u2). (10)

    For the first subsystem, it first follows from the definition of $ x_{12} $ and $ x_{21} $ that:

    ˙x12=f1(x1)+[ˉGi0]Q11x2=f1(x12)+ˉG1x21. (11)

    Then, for the second subsystem, it can be obtained that

    [˙x21˙x22]=Q11(f2(x1,x2)+[ˉG20]Q12x3=[f21(x12,x21,x22)f22(x12,x21,x22)]+ˉG2x31. (12)

    After an n-step induction, system (6) can be then transformed into

    {˙x12=f1(x12)+ˉG1x21˙x21=f21(x21,x12,2)+[Ik10]Q11ˉG2x31˙x22=f22(x21,x12,2)+[0Ik2k1]Q11ˉG2x31˙xn1=fn1(x2n,1,x1n,2)+[Ikn0]Q1n1ˉGnu1˙xn2=fn2(x2n,1,x1n,2)+[0Iknkn1]Q1n1ˉGnu1. (13)

    We use the newly obtained system (13) to derive a HOFAS model for the first-order SFS (6). In particular, the following theorem can be obtained.

    Theorem 1: Let Assumption 1 be satisfied, then under the following transformation:

    xi1=ˉG11i1(x(ip)p2|p=1i1Li1(x(0i1p)p2|p=1i1)),i=2,,n (14)

    the SFS (13) can be equivalently transformed into the following HOFAS model:

    x(n+1p)p2|p=1n=Ln(x(0np)p2|p=1n)+ˉG1nu1 (15)

    where $ L_i(\cdot) $ and $ \bar{G}_{1\sim i} $ are given iteratively by

    Li+1(x(0i+1p)p2|p=1i+1)=(˙Li+ˉG1ifi+1,1(x2i+1,1,x1i+1,2)fi+1,2(x2i+1,1,x1i+1,2)) (16)
    ˉG1i+1=[ˉG1i00Iki+1ki]Q1iˉGi+1,i=1,2,,n1. (17)

    with $ L_1(\cdot) = f_1(x_{12}) $ and $ \bar{G}_{1\sim 1} = \bar{G}_1 $.

    Proof: We prove the theorem via the n-step induction.

    Step 1: For the first subsystem, we already have that

    ˙x12=f1(x12)+ˉG1x21=L1(x12)+ˉG11x21. (18)

    Step i ($ 2\leq i\leq n-1) $: Assume that

    x(ip)p2|p=1i1=Li1(x(0i1p)p2|p=1i1)+ˉG1i1xi1 (19)

    then we prove that

    x(i+1p)p2|p=1i=Li(x(0ip)p2|p=1i)+ˉG1ixi+1,1. (20)

    It follows from (19) that:

    xi1=;ˉG11i1(x(ip)p2|p=1i1Li1(x(0i1p)p2|p=1i1)). (21)

    We also recall the following equations from system (13):

    ˙xi1=fi1(x2i,1,x1i,2)+[Iki0]Q1i1ˉGixi+1,1˙xi2=fi2(x2i,1,x1i,2)+[0Ikiki1]Q1i1ˉGixi+1,1. (22)

    Next, taking the derivatives of both sides of the (19) and substituting first equation of (22) into the result, one has that

    x(i+1p)p2|p=1i1=˙Li1+ˉG1i1˙xi1=˙Li1+ˉG1i1(fi1(x2i,1,x1i,2) +[Iki0]Q1i1ˉGixi+1,1)=˙Li1+ˉG1i1fi1(x2i,1,x1i,2) +[ˉG1i10]Q1i1ˉGixi+1,1. (23)

    By combining (23) with the last equation of (22), it can be obtained that

    x(i+1p)p2|p=1i=Li(x(0ip)p2|p=1i)+ˉG1ixi+1,1 (24)

    where

    Li= (˙Li1+ˉG1i1fi1(x2i,1,x1i,2)fi2(x2i,1,x1i,2)) (25)
    ˉG1i=[ˉG1i100Ikiki1]Q1i1ˉGi. (26)

    Step n: Based on the previous $ n-1 $ steps, we can obtain that

    xi1=ˉG11i1(x(ip)p2|p=1i1 Li1(x(0i1p)p2|p=1i1)),i=2,,n1 (27)

    and

    x(np)p2|p=1n1=Ln1(x(0n1p)p2|p=1n1) +ˉG1n1xn1 (28)

    where $ L_i(\cdot) $ and $ \bar{G}_{1\sim i} $ are given iteratively by

    Li+1(x(0i+1p)p2|p=1i+1)=(˙Li+ˉG1ifi+1,1(x2i+1,1,x1i+1,2)fi+1,2(x2i+1,1,x1i+1,2)) (29)
    ˉG1i+1=[ˉG1i00Iki+1ki]Q1iˉGi+1,i=1,2,,n1. (30)

    We recall the following equations from (13):

    ˙xn1=fn1(x2n,1,x1n,2)+[Ikn0]Q1n1ˉGnu1 (31)
    ˙xn2=fn2(x2n,1,x1n,2)+[0Iknkn1]Q1n1ˉGnu1. (32)

    With the similar discussions as that in previous steps, it can be thus obtained that

    x(n+1p)p2|p=1n=Ln(x(0np)p2|p=1n)+ˉG1nu1 (33)

    where

    Ln=(˙Ln1+ˉG1n1fn1(x2n,1,x1n,2)fn2(x2n,1,x1n,2)) (34)
    ˉG1n=[ˉG1n100Iknkn1]Q1n1ˉGn. (35)

    It follows from the definition of $ \bar{G}_{1\sim n} $ that $ \bar{G}_{1\sim n} $ is nonsingular. Therefore, the system (33) is fully-actuated.

    Remark 1: It is noted that in the HOFAS model, only part of the controller u, i.e., $ u_1 $, is used while the rest part $ u_2 $ can be removed. This is mainly because that the system (6) is “over-actuated” when the dimension of controller is greater than the state. Due to this fact, we only need part of the controller to convert the system (6) into a FAS. In fact, most of practical systems are not over-actuated, i.e., $ u_2 $ does not exist in those systems. However, there do exist “over-actuated” practical systems. For example, some real industrial systems are equipped with redundant actuators, such as spacecraft [23], ships [24], and fusion systems [25], for higher reliability and greater flexibility. These practical systems are over-actuated.

    As induced in [11], based on the FAS model (15), the following exponentially stabilizing controller can be directly given:

    u1=ˉG11n(Ln+Ax(0np)p2|p=1n) (36)

    where $ A\in\mathbb{R}^{\sum\limits_{i = 1}^{n}(n+1-i)k_i\times\sum\limits_{i = 1}^{n}(n+1-i)k_i} $, which yields the following linear time-invariant closed-loop system:

    x(n+1p)p2|p=1n=Ax(0np)p2|p=1n. (37)

    In particular, if A is chosen as a block diagonal triangular matrix in the following form:

    [A0n1,1000A0n2,2000A0n] (38)

    where

    A0nk,k=[A0kA1kAnk,k],k=1,2,,n (39)

    with $ A_{pi},\; i = 1,2,\ldots,n, \; 0\leq p\leq n-i, $ being matrices with appropriate dimensions. In this case, the closed-loop system can be written as

    {x(n)12=n1p=0Ap1x(p)12x(n1)22=n2p=0Ap2x(p)22x(1)n2=A0nxn2 (40)

    Remark 2: As is pointed out in [11], the FAS approach has two distinctive advantages compared with backstepping method. First, the FAS approach can always produce a linear time-invariant closed-loop system. Second, the FAS approach can avoid the “differential explosion” problem, which makes the method of backstepping generally not applicable to a SFS with more than 3 or 4 subsystems.

    In the previous subsection, we only assume that $ G_i,\; i = 1,2,\ldots,n $, are constant matrices. In this subsection, we consider the more general case of $ G_i,\; i = 1,2,\ldots,n $, being matrix functions, which is described by the following dynamics:

    {˙x1=f1(x1)+G1(x1)x2˙x2=f2(x1,x2)+G2(x1,x2)x3˙xn=fn(x1,x2,,xn)+Gn(x1,x2,,xn)u (41)

    where $ x_i\in\mathbb{R}^{k_i},\; i = 1,2,\ldots,n $, $ u\in\mathbb{R}^m $. Moreover, $ f_i $ and $ G_i $ are sufficiently smooth functions. In this work, it is assumed that $ k_1\leq k_2\leq \cdots\leq k_n\leq m. $ The following assumption is needed for matrix functions $ G_i(\cdot),\; i = 1,2,\ldots,n. $

    Assumption 2: For any $ x_{1\sim i}\in\mathbb{R}^{k_1+\cdots+k_i} $, the matrix function $ G_i(x_{1\sim i}) $ is of full-row rank, i.e., $ {\rm rank}(G_i) = k_i,\; i = 1, 2,\ldots,n $.

    It follows from Assumption 2 that there exist sufficiently smooth functions $ Q_i: \mathbb{R}^{k_1+\cdots+k_i}\rightarrow\mathbb{R}^{k_{i+1}\times k_{i+1}} $ and $ \bar{G}_i:\mathbb{R}^{k_1+\cdots+k_i} \rightarrow \mathbb{R}^{k_i\times k_{i}} $such that:

    Gi()Qi()=[ˉGi()0] (42)

    where $ Q_i(\cdot) $ and $ \bar{G}_i(\cdot) $ are nonsingular for any $ x_{1\sim i}\in\mathbb{R}^{k_1+\cdots+k_i} $. Define a differential homeomorphism from $ x = \{x_1^T,x_2^T,\ldots, x_n^T\}^T $to $ \bar{x} = \{\bar{x}_1^T,\bar{x}_2^T,\ldots,\bar{x}_n^T\}^T $ as follows:

    ˉx1= x1ˉx2= Q11(x1)x2ˉx3= Q12(x1,x2)x3ˉxn= Q1n1(x1,x2,,xn1)xn. (43)

    Then the inverse mapping from $ \bar{x} $ to x can be written as follows:

    x1= ˉx1x2= Q1(ˉx1)ˉx2ˉQ1(ˉx1)ˉx2x3= Q2(ˉx1,ˉQ1(ˉx1)ˉx2)ˉx3ˉQ2(ˉx1,ˉx2)ˉx3xn= Qn1(ˉx1,ˉQ1(¯x1)ˉx2,,ˉQn2(ˉx1,,ˉxn2)ˉxn1)ˉxn ˉQn1(ˉx1,,ˉxn1)ˉxn. (44)

    Furthermore, define $ \bar{x}_{i1} $, $ \bar{x}_{i2} $ and $ \bar{u}_1 $ as follows:

    {ˉxi1=[Iki10]ˉxi,i=2,,nˉxi2=[0Ikiki1]ˉxi,i=2,,nˉu1=[0Ikn]ˉQ1nu. (45)

    Meanwhile, for notation convenience, we also define that $ \bar{x}_{12} = \bar{x}_1. $ It then follows from (43) and (44) that the following equivalent form of system (41) can be obtained:

    {˙ˉx12=ˉf1(ˉx12)+ˉG1(ˉx12)ˉx21˙ˉx21=ˉf21(ˉx21,ˉx12,2)+[Ik10]ˉQ11ˉG2(ˉx21,ˉx12,2)ˉx3˙ˉx22=ˉf22(ˉx21,ˉx12,2)+[0Ik2k1]ˉQ11ˉG2(ˉx21,ˉx12,2)ˉx3˙ˉxn1=ˉfn1(ˉx2n,1,ˉx1n,2)+[Ikn10]ˉQ1n1ˉGn(ˉx2n,1,ˉx1n,2)ˉu1˙ˉxn2=ˉfn2(ˉx2n,1,ˉx1n,2)+[0Iknkn1]ˉQ1n1ˉGn(ˉx2n,1,ˉx1n,2)ˉu1 (46)

    where

    ˉf1(ˉx12)=f1(ˉx12) (47)

    and for $ i = 2,\ldots, n, $

    ˉfi1(ˉx2i,1,ˉx1i,2)=[Iki10]ˉQ1i1 ×(˙ˉQi1ˉxi+fi(ˉx1,ˉQ1ˉx2,,ˉQi1ˉxi))ˉfi2(ˉx2i,1,ˉx1i,2)=[0Ikiki1]ˉQ1i1 ×(˙ˉQi1ˉxi+fi(ˉx1,ˉQ1ˉx2,,ˉQi1ˉxi)). (48)

    Based on the equivalent system (46), the following HOFAS model can be derived.

    Theorem 2: Let Assumption 2 be satisfied, then under the following transformation:

    ˉxi1=ˉG11i1(ˉx(ip)p2|p=1i1Li1(ˉx(0i1p)p2|p=1i1)),i=2,,n (49)

    the SFS (46) can be equivalently transformed into the following HOFAS model:

    ˉx(0n+1p)p2|p=1n=Ln(ˉx(0np)p2|p=1n)+ˉG1n(ˉx(0np)p2|p=1n)ˉu1 (50)

    where $ L_i(\cdot) $ and $ \bar{G}_{1\sim i}(\cdot) $ are given iteratively by

    Li+1(ˉx(0i+1p)p2|p=1i+1)=(˙Li+˙ˉG1iˉxi+1,1+ˉG1iˉfi+1,1ˉfi+1,2) (51)
    ˉG1i+1(ˉx(0i+1p)p2|p=1i+1)=[ˉG1i00Iki+1ki]ˉQ1iˉGi+1,i=1,2,,n1 (52)

    with $ L_1 = \bar{f}_1(\bar{x}_{12}) $ and $ \bar{G}_{1\sim 1} = \bar{G}_1(\bar{x}_{12}) $.

    Proof: We can also prove the theorem via the n-step induction.

    Step 1: For the first subsystem, we already have

    ˙ˉx12= ˉf1(ˉx12)+ˉG1(ˉx12)ˉx21= L1(ˉx12)+ˉG11(ˉx12)ˉx21. (53)

    Step i ($ 2\leq i\leq n-1) $: Assume that

    ˉx(ip)p2|p=1i1=Li1(ˉx(0i1p)p2|p=1i1)+ˉG1i1(ˉx(0i1p)p2|p=1i1)ˉxi1 (54)

    then we prove that

    ˉx(i+1p)p2|p=1i=Li(ˉx(0ip)p2|p=1i)+ˉG1i(ˉx(0ip)p2|p=1i)ˉxi+1,1. (55)

    It follows from (54) that:

    ˉxi1=ˉG11i1()(ˉx(ip)p2|p=1i1Li1(ˉx(0i1p)p2|p=1i1)). (56)

    We also recall the following equations from system (46):

    ˙ˉxi1= ˉfi1(ˉx2i,1,ˉx1i,2) +[Iki0]ˉQ1i1ˉGi(ˉx2i,1,ˉx1i,2)ˉxi+1,1˙ˉxi2= ˉfi2(ˉx2i,1,ˉx1i,2) +[0Ikiki1]ˉQ1i1ˉGi(ˉx2i,1,ˉx1i,2)ˉxi+1,1. (57)

    Next, taking the derivatives of both sides of the (54) and substituting first equation of (57) into the result, one has that

    ˉx(i+1p)p2|p=1i1= ˙Li1+˙ˉG1i1ˉxi1+ˉG1i1˙ˉxi1= ˙Li1+˙ˉG1i1ˉxi1 +ˉG1i1(ˉfi1+[Iki0]ˉQ1i1ˉGiˉxi+1,1)= ˙Li1+˙ˉG1i1ˉxi1+ˉG1i1ˉfi1 +[ˉG1i10]ˉQ1i1ˉGiˉxi+1,1. (58)

    By combining (58) with the last equation of (57), it can be obtained that

    ˉx(i+1p)p2|p=1i= Li(ˉx(0ip)p2|p=1i) +ˉG1i(ˉx(0ip)p2|p=1i)ˉxi+1,1 (59)

    where

    Li=(˙Li1+˙ˉG1i1ˉxi1+ˉG1i1ˉfi1ˉfi2)ˉG1i=[ˉG1i100Ikiki1]ˉQ1i1ˉGi. (60)

    Step n: Based on the previous $ n-1 $ steps, we can obtain that

    ˉxi1= ˉG11i1(ˉx(ip)p2|p=1i1 Li1(ˉx(0i1p)p2|p=1i1)),i=2,,n1 (61)

    and

    ˉx(np)p2|p=1n1= Ln1(ˉx(0n1p)p2|p=1n1) +ˉG1n1(ˉx(0n1p)p2|p=1n1)ˉxn1 (62)

    where $ L_i(\cdot) $ and $ \bar{G}_{1\sim i}(\cdot) $ are given iteratively by

    Li+1(ˉx(0i+1p)p2|p=1i+1)=(˙Li+˙ˉG1iˉxi+1,1+ˉG1iˉfi+1,1ˉfi+1,2) (63)
    ˉG1i+1(ˉx(0i+1p)p2|p=1i+1)=[ˉG1i00Iki+1ki]ˉQ1iˉGi+1,i=1,2,,n1. (64)

    We recall the following equations from (46):

    ˙ˉxn1= ˉfn1(ˉx2n,1,ˉx1n,2) +[Ikn0]ˉQ1n1ˉGn(ˉx2n,1,ˉx1n,2)ˉu1˙ˉxn2= ˉfn2(ˉx2n,1,ˉx1n,2) +[0Iknkn1]ˉQ1n1ˉGn(ˉx2n,1,ˉx1n,2)ˉu1. (65)

    With the similar discussions as that in previous steps, it can be thus obtained that

    ˉx(0n+1p)p2|p=1n= Ln(ˉx(0np)p2|p=1n) +ˉG1n(ˉx(0np)p2|p=1n)ˉu1 (66)

    where

    Ln=(˙Ln1+˙ˉG1n1ˉxn1+ˉG1n1ˉfn1fn2)ˉG1n=[ˉG1n100Iknkn1]ˉQ1n1ˉGn. (67)

    It follows from the definition of $ \bar{G}_{1\sim n} $ that $ \bar{G}_{1\sim n} $ is nonsingular for any $ \left(\bar{x}_{p2}^{(0\sim n-p)}\Big|_{p = 1\sim n}\right)\in\mathbb{R}^{\sum\limits_{i = 1}^n(n-i)k_i} $. Therefore, the system (66) is a FAS.

    In [11], the author proposes more generalized SFSs of both second-order and high-order. Meanwhile, both second-order and high-order SFSs are converted into FAS models. In this subsection, we further extend the results in [11] to second-order and high-order SFSs with increasing dimensions. Proofs of all theorems in this part are omitted as they can be easily established by applying the induction method introduced in Section III together with the results in [11].

    In this subsection, the following second-order SFS with increasing dimensions is first discussed:

    {¨x1=f1(x(01)1)+G1x2¨x2=f2(x(02)1,x(01)2)+G2x3¨xn1=fn(x(02)i|i=1n2,x(01)n1)+Gn1xn¨xn=fn(x(02)i|i=1n1,x(01)n)+Gnu (68)

    where $ x_i\in\mathbb{R}^{k_i},\; i = 1,2,\ldots,n $, $ u\in\mathbb{R}^m $, $ G_i\in\mathbb{R}^{k_i\times k_{i+1}}, \; i = 1,2,\ldots, n-1, $ and $ G_n\in\mathbb{R}^{k_n\times m} $. We adopt the similar definition for $ x_{i+1,1}\in\mathbb{R}^{k_i} $,$ x_{i+1,2}\in\mathbb{R}^{k_{i+1}-k_i} $, $ f_{i1}\in\mathbb{R}^{k_i} $, $ u_1\in\mathbb{R}^{k_n} $ and $ u_2\in\mathbb{R}^{m-k_n} $ as (8)−(10). Then system (68) can be transformed into the following equivalent form:

    {¨x12=f1(x(01)1)+ˉG1x21¨x21=f21(x(01)21,x(01)12,2)+[Ik10]Q11ˉG2x31¨x22=f22(x(01)21,x(01)12,2)+[0Ik2k1]Q11ˉG2x31¨xn1=fn1(x(02)12,x(02)ij|i=2n1,j=12,x(01)nj|j=12)+[Ikn0]Q1n1ˉGnu1¨xn2=fn2(x(02)12,x(02)ij|i=2n1,j=12,x(01)nj|j=12)+[0Iknkn1]Q1n1ˉGnu1. (69)

    Based on the system (69), we can derive the HOFAS model for the second-order SFSs with increasing dimensions. In particular, we present the following theorem.

    Theorem 3: Let Assumption 1 be satisfied, then under the following transformation:

    xi1=ˉG11i1(x(2i2p)p2|p=1i1Li1(x(02i12p)p2|p=1i1)),i=2,,n (70)

    the SFS (69) can be equivalently transformed into the following HOFAS model:

    x(2n+22p)p2|p=1n=Ln(x(02n+12p)p2|p=1n)+ˉG1nu1 (71)

    where $ L_i(\cdot) $ and $ G_{1\sim i} $ are given iteratively by

    Li+1()=(¨Li+ˉG1ifi+1,1fi+1,2) (72)
    ˉG1i+1=[ˉG1i00Iki+1ki]Q1iˉGi+1,i=1,2,,n1 (73)

    with $ L_1 = f_1\left(x_{12}^{(0\sim 1)}\right) $ and $ \bar{G}_{1\sim 1} = \bar{G}_1 $.

    Next, similar to the first-order SFSs, the more general case where $ G_i,\; i = 1,2,\ldots,n, $ are assumed to be matrix functions can be also considered. In particular, the following second-order SFS is discussed:

    {¨x1=f1(x(01)1)+G1(x(01)1)x2¨x2=f2(x(02)1,x(01)2)+G2(x(02)1,x(01)2)x3¨xn1=fn(x(02)i|i=1n2,x(01)n1)+Gn1(x(02)i|i=1n2,x(01)n1)xn¨xn=fn(x(02)i|i=1n1,x(01)n)+Gn(x(02)i|i=1n1,x(01)n)u (74)

    where $ x_i\in\mathbb{R}^{k_i},\; i = 1,2,\ldots,n $, $ u\in\mathbb{R}^m $. Moreover, $ f_i $ and $ G_i $ are sufficiently smooth functions. In this work, it is assumed that $ k_1\leq k_2\leq \cdots\leq k_n\leq m. $ The following assumption is needed for matrix functions $ G_i,\; i = 1,2,\ldots,n. $

    Assumption 3: For any $\left(x_p^{(0\sim 2)}\Big|_{p = 1\sim i-1},x_{i}^{(0\sim 1)}\right)\in $ $ \mathbb{R}^{3k_1+\cdots+3k_{i-1}+2k_i} $, the matrix function $ G_i(\cdot) $ is of full-row rank, i.e., $ {\rm rank}(G_i) = k_i, \; i = 1,2,\ldots,n $.

    It follows from Assumption 3 that there exist differentiable functions $ Q_i: \mathbb{R}^{3k_1+\cdots+3k_{i-1}+2k_i}\rightarrow\mathbb{R}^{k_{i+1}\times k_{i+1}} $ and $\bar{G}_i: $ $ \mathbb{R}^{3k_1+\cdots+3k_{i-1}+2k_i}\rightarrow\mathbb{R}^{k_i\times k_{i}} $ such that:

    Gi()Qi()=[ˉGi()0] (75)

    where $ Q_i $ and $ \bar{G}_i $ are of full rank. Define the mapping from $ x = \{x_1^T,x_2^T,\ldots,x_n^T\}^T $ to $ \bar{x} = \{\bar{x}_1^T,\bar{x}_2^T,\ldots,\bar{x}_n^T\}^T $ as follows:

    ˉx1= x1ˉx2= Q11(x(01)1)x2ˉx3= Q12(x(02)1,x(01)2)x3ˉxn= Q1n1(x(02)i|i=1n1,x(01)n)xn. (76)

    Then the inverse mapping from $ \bar{x} $ to x can be written as follows:

    x1= ˉx1x2= Q1(ˉx(01)1)ˉx2ˉQ1(ˉx(01)1)ˉx2x3= Q2(ˉx(02)1,(ˉQ1(ˉx(01)1)ˉx2)(01))ˉx3  ˉQ2(ˉx(02)1,ˉx(01)2)ˉx3xn= Qn1((ˉQi1ˉxi)(02)|i=1n2,(ˉQn2ˉxn1)(01))ˉxn ˉQn1(ˉx(0ni)i|i=1n1)ˉxn. (77)

    Define $ \bar{x}_{i1} $, $ \bar{x}_{i2} $ and $ \bar{u}_1 $ as follows:

    {ˉxi1=[Iki10]ˉxi,i=2,,nˉxi2=[0Ikiki1]ˉxi,i=2,,nˉu1=[0Ikn]ˉQ1nu. (78)

    Meanwhile, for notation convenience, we also define that $ \bar{x}_{12} = \bar{x}_1. $ It then follows from (76) and (77) that system (74) can be equivalently transformed into the following form:

    {¨ˉx1=ˉf1(ˉx(01)12)+ˉG1(ˉx(01)12)ˉx21¨ˉx21=ˉf21(ˉx(02)12,ˉx(01)21,ˉx(01)22)+[Ik10]ˉQ11ˉG2(ˉx(02)12,ˉx(01)21,ˉx(01)22)ˉx3¨ˉx22=ˉf22(ˉx(02)12,ˉx(01)21,ˉx(01)22)+[0Ik2k1]ˉQ11ˉG2(ˉx(02)12,ˉx(01)21,ˉx(01)22)ˉx3¨ˉxn1=ˉfn1(ˉx(0n+1)12,ˉx(0ni+2)ij|i=2n1,j=12,ˉx(01)nj|j=12)+[Ikn10]ˉQ1n1×ˉGn(ˉx(0n+1)12,ˉx(0ni+2)ij|i=2n1,j=12,ˉx(01)nj|j=12)ˉu1¨ˉxn2=ˉfn2(ˉx(0n+1)12,ˉx(0ni+2)ij|i=2n1,j=12,ˉx(01)nj|j=12)+[0Iknkn1]ˉQ1n1×ˉGn(ˉx(0n+1)12,ˉx(0ni+2)ij|i=2n1,j=12,ˉx(01)nj|j=12)ˉu1 (79)

    where

    ˉf1(ˉx(01)12)=f1(ˉx(01)12) (80)

    and for $ i = 2,\ldots,n, $

    ˉfi1()=[Iki10]ˉQ1i1(¨ˉQi1ˉxi2˙ˉQi1˙ˉxi+fi)ˉfi2()=[0Ikiki1]ˉQ1i1(¨ˉQi1ˉxi2˙ˉQi1˙ˉxi+fi). (81)

    The FAS model for (79) can be illustrated via the following theorem.

    Theorem 4: Let Assumption 4 be satisfied, then under the following transformation:

    ˉxi1=ˉG11i1(ˉx(2i2p)p2|p=1i1Li1(ˉx(02i12p)p2|p=1i1)),i=2,,n (82)

    the strict-feedback system (79) can be equivalently transformed into the following HOFAS model:

    ˉx(2i+22p)p2|p=1i=Ln(ˉx(02n+12p)p2|p=1n)+ˉG1n(ˉx(02n+12p)p2|p=1n)ˉu1 (83)

    where $ L_i(\cdot) $ and $ \bar{G}_{1\sim i}(\cdot) $ are give iteratively by

    Li+1=(¨Li+¨ˉG1iˉxi+1,1+2˙ˉG1i˙ˉxi+1,1+ˉG1iˉfi+1,1ˉfi+1,2) (84)
    ˉG1i+1=[ˉG1i00Iki+1ki]ˉQ1iˉGi+1,i=1,2,,n1 (85)

    with $ L_1 = \bar{f}_1\left(\bar{x}_{12}^{(0\sim 1)}\right) $ and $ \bar{G}_{1\sim 1} = \bar{G}_1 $.

    In this subsection, a generalized form of SFSs, where each subsystem is high-order, is discussed. The high-order SFSs are firstly considered by Duan in [11]. However, [11] only consider SFSs with uniform dimension. While the dimensions of all subsystems are assumed to be increasing in this work. In particular, we consider the following system:

    {x(m1)1=f1(x(0m11)1)+G1x2x(m2)2=f2(x(0m1)1,x(0m21)2)+G2x2x(mn1)n1=fn(x(0mi)i|i=1n2,x(0mn11)n1)+Gn1xnx(mn)n=fn(x(0mi)i|i=1n1,x(0mn1)n)+Gnu (86)

    where $ x_i\in\mathbb{R}^{k_i},\; i = 1,2,\ldots,n $, $ u\in\mathbb{R}^m $, $ G_i\in\mathbb{R}^{k_i\times k_{i+1}}, \; i = 1, 2,\ldots, n-1, $ and $ G_n\in\mathbb{R}^{k_n\times m} $. We adopt the similar definition for $ x_{i+1,1}\in\mathbb{R}^{k_i} $,$ x_{i+1,2}\in\mathbb{R}^{k_{i+1}-k_i} $, $ f_{i1}\in\mathbb{R}^{k_i} $, $ u_1\in\mathbb{R}^{k_n} $ and $ u_2\in \mathbb{R}^{m-k_n} $ as (8)−(10). Then system (86) can be transformed into the following equivalent form:

    {x(m1)12=f1(x(0m11)12)+ˉG1x21x(m2)21=f21(x(0m1)12,x(0m21)21,x(0m21)22)+[Ik10]Q11ˉG2x31x(m2)22=f22(x(0m1)12,x(0m21)21,x(0m21)22)+[0Ik2k1]Q11ˉG2x31x(mn)n1=fn1(x(0m1)12,x(0mi)ij|i=2n1,j=12,x(0mn1)nj|j=12)+[Ikn0]Q1n1ˉGnu1x(mn)n2=fn2(x(0m1)12,x(0mi)ij|i=2n1,j=12,x(0mn1)nj|j=12)+[0Iknkn1]Q1n1ˉGnu1. (87)

    For notation convenience, we denote that for $ j\geq i\geq 1 $,

    qij=mi+mi+1++mj. (88)

    Based on the equivalent form (87), the HOFAS model for high-order SFS (86) can be obtained, which is given by the following theorem.

    Theorem 5: Let Assumption 1 be satisfied, then under the following transformation:

    xi1=ˉG11i1(x(qpi1)p2|p=1i1Li1(x(0qpi11)p2|p=1i1)),i=2,,n (89)

    the SFS (87) can be equivalently transformed into the following HOFA model:

    x_{p2}^{(q_{p\sim n})}\Big|_{p = 1\sim n} = L_n\left(x_{p2}^{(0\sim q_{p\sim n}-1)}\Big|_{p = 1\sim n}\right)+\bar{G}_{1\sim n}u_1 (90)

    where $ L_i(\cdot) $ and $ \bar{G}_{1\sim i} $ are given iteratively by

    L_{i+1} = \left(\begin{array}{c} {L}_i^{(m_{i+1})}+\bar{G}_{1\sim i}f_{i+1,1} \\ f_{i+1,2}\end{array}\right) (91)
    \bar{G}_{1\sim i+1} = \left[\begin{array}{cc}\bar{G}_{1\sim i}&0\\0&I_{k_{i+1}-k_i}\end{array}\right]Q_i^{-1}\bar{G}_{i+1}, (92)
    i = 1,2,\ldots,n-1 (93)

    with $ L_1 = f_1\left(x_{12}^{(0\sim m_1-1)}\right) $ and $ \bar{G}_{1\sim 1} = \bar{G}_1 $.

    Similar to the first-order and second-order SFSs, we can also consider the following more general high-order SFS with $ G_i,\; i = 1,2,\ldots,n $, being matrix functions:

    \left\{ \begin{aligned} &{x}_1^{(m_1)} = f_1\left(x_1^{(0\sim m_1-1)}\right)+G_1\left(x_1^{(0\sim m_1-1)}\right)x_2 \\& {x}_2^{(m_2)} = f_2\left(x_1^{(0\sim m_1)},x_2^{(0\sim m_2-1)}\right)\\& \;\;\;\; \; \; \; \; \; \; \; \; +G_2\left(x_1^{(0\sim m_1)},x_2^{(0\sim m_2-1)}\right)x_2 \\&\;\;\; \; \; \; \; \; \; \; \; \; \vdots \\& {x}_{n-1}^{(m_{n-1})} = f_n\left(x_i^{(0\sim m_i)}\Big|_{i = 1\sim n-2},x_{n-1}^{(0\sim m_{n-1}-1)}\right)\\& \; \; \; \; \;\;\;\; \; \; \;\;\;\; \; +G_{n-1}\left(x_i^{(0\sim m_i)}\Big|_{i = 1\sim n-2},x_{n-1}^{(0\sim m_{n-1}-1)}\right)x_n\\& {x}_n^{(m_n)} = f_n\left(x_i^{(0\sim m_i)}\Big|_{i = 1\sim n-1},x_{n}^{(0\sim m_{n}-1)}\right)\\& \; \; \; \; \; \;\;\;\; \;\;\; \; \; +G_n\left(x_i^{(0\sim m_i)}\Big|_{i = 1\sim n-1},x_{n}^{(0\sim m_{n}-1)}\right)u\end{aligned}\right. (94)

    where $ x_i\in\mathbb{R}^{k_i},\; i = 1,2,\ldots,n, $ and $ u\in\mathbb{R}^m $. Moreover, $ f_i(\cdot) $ and $ G_i(\cdot) $ are sufficiently smooth functions. In this work, we assume that $ k_1\leq k_2\leq \cdots\leq k_n\leq m. $ The following assumption is needed for system (94).

    Assumption 5: For any $ \left(x_p^{(0\sim m_p)}\Big|_{p = 1\sim i-1},x_{i}^{(0\sim m_{i}-1)}\right)\in $ $\mathbb{R}^{\sum\limits_{p = 1}^{i-1}(m_p+1)k_p+m_ik_i} $, the matrix function $ G_i $ is of full-row rank, i.e., $ {\rm rank}(G_i) = k_i,\; i = 1,2,\ldots,n $.

    It follows from Assumption 5 that there exist sufficiently smooth functions $ Q_i(\cdot) $ and $ \bar{G}_i(\cdot) $ for $ i = 1,2,\ldots,n $ such that:

    G_i(\cdot)Q_i(\cdot) = \left[\begin{array}{cc} \bar{G}_i(\cdot) & 0 \end{array}\right] (95)

    where $ Q_i(\cdot) $ and $ \bar{G}_i(\cdot) $ are nonsingular. Then we define $ \bar{x} = \{\bar{x}_1^T,\bar{x}_2^T,\ldots,\bar{x}_n^T\}^T $ as follows:

    \begin{split} \bar{x}_1 = \ &x_1\\ \bar{x}_2 = \ &Q_1^{-1}\left(x_1^{(0\sim m_1-1)}\right)x_2\\ \bar{x}_3 = \ &Q_2^{-1}\left( x_1^{(0\sim m_1)},x_2^{(0\sim m_2-1)}\right)x_3\\ \vdots\\ \bar{x}_n = \ &Q_{n-1}^{-1} \left(x_i^{(0\sim m_i)}\Big|_{i = 1\sim n-1},x_{n}^{(0\sim m_{n}-1)}\right)x_n. \end{split} (96)

    And the inverse mapping from $ \bar{x} $ to x can be written as follows:

    \begin{split} x_1 = \ &\bar{x}_1\\ x_2 = \ &Q_1\left(\bar{x}_1^{(0\sim m_1-1)}\right)\bar{x}_2\triangleq \bar{Q}_1\left(\bar{x}_1^{(0\sim m_1-1)}\right)\bar{x}_2\\ x_3 = \ &Q_2\left(\bar{x}_1^{(0\sim m_1)},\left(\bar{Q}_1\left(\bar{x}_1^{(0\sim m_1-1)}\right)\bar{x}_2\right)^{(0\sim m_2-1)}\right)\bar{x}_3\\ \ \triangleq\;& \bar{Q}_2\left(\bar{x}_1^{(0\sim m_1+m_2-2)},\bar{x}_2^{(0\sim m_2-1)}\right)\bar{x}_3\\ \vdots\\ x_n = \ &Q_{n-1}\bigg(\left(\bar{Q}_{i-1}\bar{x}_i\right)^{(0\sim m_i)}\Big|_{i = 1\sim n-2} \end{split}
    \begin{split} \ &\left(\bar{Q}_{n-2}\bar{x}_{n-1}\right)^{(0\sim m_{n-1}-1)}\bigg)\bar{x}_n\\ \triangleq \ &\bar{Q}_{n-1}\left(\bar{x}_i^{(0\sim q_{i\sim n-1}-n+i)}\Big|_{i = 1\sim n-1}\right)\bar{x}_n. \end{split} (97)

    Furthermore, define $ \bar{x}_{i1} $, $ \bar{x}_{i2} $ and $ \bar{u}_1 $ as follows:

    \left\{\begin{aligned} &\bar{x}_{i1} = \left[\begin{array}{cc}I_{k_{i-1}}&0\end{array}\right]\bar{x}_i,\; &&i = 2,\ldots,n\\ &\bar{x}_{i2} = \left[\begin{array}{cc}0&I_{k_i-k_{i-1}}\end{array}\right]\bar{x}_i,\; &&i = 2,\ldots,n\\ &\bar{u}_1 = \left[\begin{array}{cc}0&I_{k_n}\end{array}\right]\bar{Q}_n^{-1}u.\end{aligned}\right. (98)

    Meanwhile, for notation convenience, we also define that $ \bar{x}_{12} = \bar{x}_1. $ It then follows from (96) and (97) that the following equivalent form of system (94) can be obtained:

    \left\{ \begin{aligned} \bar{x}_{12}^{(m_1)} = \;&\bar{f}_1\left(\bar{x}_{12}^{(0\sim m_1-1)}\right)+ \bar{G}_1\left(\bar{x}_1^{(0\sim m_1-1)}\right)\bar{x}_{21} \\ \bar{x}_{21}^{(m_2)} =\;& \bar{f}_{21}\left(\bar{x}_{12}^{(0\sim m_1+m_2-1)},\bar{x}_{21}^{(0\sim m_2-1)},\bar{x}_{22}^{(0\sim m_2-1)}\right)\\ & +\left[\begin{array}{cc}I_{k_1}&0\end{array}\right]\bar{Q}_1^{-1} \\ & \times\bar{G}_2\left(\bar{x}_{12}^{(0\sim m_1+m_2-1)},\bar{x}_{21}^{(0\sim m_2-1)},\bar{x}_{22}^{(0\sim m_2-1)}\right)\bar{x}_3 \\ \bar{x}_{22}^{(m_2)} =\;& \bar{f}_{22}\left(\bar{x}_{12}^{(0\sim m_1+m_2-1)},\bar{x}_{21}^{(0\sim m_2-1)},\bar{x}_{22}^{(0\sim m_2-1)}\right) \\ & +\left[\begin{array}{cc}0&I_{k_2-k_1}\end{array}\right]\bar{Q}_1^{-1}\\ & \times\bar{G}_2\left(\bar{x}_{12}^{(0\sim m_1+m_2-1)},\bar{x}_{21}^{(0\sim m_2-1)},\bar{x}_{22}^{(0\sim m_2-1)}\right)\bar{x}_3\\& \vdots \\ \bar{x}_{n1}^{(m_n)} =\; & \bar{f}_{n1}\bigg(\bar{x}_{12}^{(0\sim q_{1\sim n}-n+1)}, \\& \bar{x}_{ij}^{(0\sim q_{i\sim n}-n+i)}\Big|_{i = 2\sim n-1,j = 1\sim2},\bar{x}_{nj}^{(0\sim m_n-1)}\Big|_{j = 1\sim2}\bigg) \\& +\left[\begin{array}{cc}I_{k_{n-1}}&0\end{array}\right]\bar{Q}^{-1}_{n-1}\bar{G}_n\bigg(\bar{x}_{12}^{(0\sim q_{1\sim n}-n+1)}, \\& \bar{x}_{ij}^{(0\sim q_{i\sim n}-n+i)}\Big|_{i = 2\sim n-1,j = 1\sim2},\bar{x}_{nj}^{(0\sim m_n-1)}\Big|_{j = 1\sim2}\bigg)\bar{u}_1, \\& \bar{x}_{n2}^{(m_n)} = \bar{f}_{n2}\bigg(\bar{x}_{12}^{(0\sim q_{1\sim n}-n+1)}, \\ &\bar{x}_{ij}^{(0\sim q_{i\sim n}-n+i)}\Big|_{i = 2\sim n-1,j = 1\sim2},\bar{x}_{nj}^{(0\sim m_n-1)}\Big|_{j = 1\sim2}\bigg) \\ &+\left[\begin{array}{cc}0&I_{k_n-k_{n-1}}\end{array}\right]\bar{Q}^{-1}_{n-1}\bar{G}_n\bigg(\bar{x}_{12}^{(0\sim q_{1\sim n}-n+1)}, \\& \bar{x}_{ij}^{(0\sim q_{i\sim n}-n+i)}\Big|_{i = 2\sim n-1,j = 1\sim2},\bar{x}_{nj}^{(0\sim m_n-1)}\Big|_{j = 1\sim2}\bigg)\bar{u}_1\end{aligned}\right. (99)

    where

    \bar{f}_1 = f_1 (100)

    and for $ i = 2,\ldots,n $

    \bar{f}_{i1} = \left[\begin{array}{cc}I_{k_{i-1}}&0\end{array}\right]\bar{Q}^{-1}_{i-1}\left(-\sum\limits_{k = 1}^{m_i}C_{m_i}^{k}\bar{Q}^{(k)}_{i-1}\bar{x}_i^{(m_i-k)}+f_i\right) (101)
    \bar{f}_{i2} = \left[\begin{array}{cc}0&I_{k_i-k_{i-1}}\end{array}\right]\bar{Q}^{-1}_{i-1}\left(-\sum\limits_{k = 1}^{m_i}C_{m_i}^{k}\bar{Q}^{(k)}_{i-1}\bar{x}_i^{(m_i-k)}+f_i\right) (102)

    with C being the binomial coefficient defined by

    C_{m_i}^k = \frac{m_i!}{(m_i-k)!k!}. (103)

    Based on the equivalent system (99), the following result can be obtained for system (94).

    Theorem 6: Let Assumption 5 be satisfied, then under the following transformation:

    \begin{split} \bar{x}_{i1} = &\bar{G}_{1\sim i-1}^{-1}\bigg( \bar{x}_{p2}^{(q_{p\sim i-1})}\Big|_{p = 1\sim i-1}\\ &-L_{i-1}\left(\bar{x}_{p2}^{(0\sim q_{p\sim i-1}-1)}\Big|_{p = 1\sim i-1}\right)\bigg), \; i = 2,\ldots, n \end{split} (104)

    the SFS (99) can be equivalently transformed into the following HOFAS model:

    \bar{x}_{p2}^{(q_{p\sim n})}\Big|_{p = 1\sim n} = L_n+\bar{G}_{1\sim n}u_1 (105)

    where $ L_i(\cdot) $ and $ \bar{G}_{1\sim i}(\cdot) $ are given iteratively by

    L_{i+1} = \left(\begin{array}{c} {L}_i^{(m_{i+1})}+\displaystyle \sum\limits_{k = 1}^{m_{i+1}}C_{m_{i+1}}^{k}\bar{G}^{(k)}_{1\sim i}\bar{x}_{i1}^{(m_{i+1}-k)}+\bar{G}_{1\sim i}f_{i+1,1} \\ f_{i+1,2}\end{array}\right) (106)
    \bar{G}_{1\sim i+1} = \left[\begin{array}{cc}\bar{G}_{1\sim i}&0\\0&I_{k_{i+1}-k_i}\end{array}\right]Q_i^{-1}\bar{G}_{i+1}, (107)
    i = 1,2,\ldots,n-1 (108)

    with $ L_1 = \bar{f}_1\left(\bar{x}_{12}^{(0\sim m_1-1)}\right) $ and $ \bar{G}_{1\sim 1} = \bar{G}_1\left(\bar{x}_1^{(0\sim m_1-1)}\right) $.

    Remark 3: Different from (43), the mapping (96) is not a differential homeomorphism from x to $ \bar{x} $ as it is not 1-1. However, it can be verified that if the system (105) is stabilized under controller $ u_1 $, then the system (94) is also stabilized. Therefore, the HOFAS model (105) is applicable for control of system (94).

    In this section, we provide two numerical examples to illustrate the effectiveness of our results.

    Consider the following first-order system:

    \qquad\qquad\qquad \dot{x}_1 = 0.1x_1+\left[\begin{array}{cc} 1&1 \end{array}\right]x_2 (109a)
    \qquad\qquad\qquad \dot{x}_2 = x_1x_2+\left[\begin{array}{ccc} 0&1&1\\1&2&2 \end{array}\right]u. (109b)

    It can be seen that system (109) is a first-order SFS with increasing dimensions. By using the method in this paper, we can obtain the following FAS model of system (109):

    \ddot{x}_{12} = 0.1\dot{x}_{12}+x_{12}x_{21}+\left[\begin{array}{cc}0&1\end{array}\right]u_1\setcounter{equation}{110} (110)
    \dot{x}_{22} = x_{12}x_{22}+\left[\begin{array}{cc}1&2\end{array}\right]u_1. (111)

    Then, we can design the following controller:

    u_1 = \left[\begin{array}{cc}-2&1\\1&0\end{array}\right]\left[\begin{array}{c}-2.1\dot{x}_{12}-x_{12}-x_{12}x_{21}\\-x_{12}x_{22}-x_{22}\end{array}\right]. (112)

    The response is then shown by Fig. 1. The effectiveness of our controller can be thus illustrated.

    Figure  1.  Response for system (109) with FAS design.

    Consider the following second-order system:

    \qquad\qquad \ddot{x}_1 =-0.1\dot{x}_1x_1+\left[\begin{array}{cc} 2&3 \end{array}\right]x_2 (113a)
    \qquad\qquad\ddot{x}_2 = \dot{x}_2+\left[\begin{array}{ccc} 0&1&2\\3&2&1 \end{array}\right]u. (113b)

    It can be seen that system (113) is a second-order SFS with increasing dimensions. By using the method in this paper, we can obtain the following FAS model of system (113):

    {x}^{(4)}_{12} = (-0.1\dot{x}_{12}x_{12})^{(2)}+\dot{x_{21}}+\left[\begin{array}{cc}0&1\end{array}\right]u_1\setcounter{equation}{114} (114)
    {x}^{(2)}_{22} = \dot{x}_{22}+\left[\begin{array}{cc}3&2\end{array}\right]u_1. (115)

    Then we can design the following controller:

    \begin{split} u_1 =\;& \left[\begin{array}{cc}-\dfrac{2}{3}& \dfrac{1}{3} \\1&0\end{array}\right]\\ &\times\left[\begin{array}{c}(0.1\dot{x}_{12}x_{12})^{(2)}-4x^{(3)}_{12}-6x^{(2)}_{12}-4x^{(1)}_{12}-4x_{12}-\dot{x_{21}}\\-3\dot{x}_{22}-x_{22}\end{array}\right]. \end{split} (116)

    The response is then shown by Fig. 2. The effectiveness of our controller can be thus illustrated.

    Figure  2.  Response for system (113) with FAS design.

    In this paper, we have provided HOFAS models for SFSs with increasing and increasing dimensions. First-order, second-order and high-order cases have been all considered. We first considered SFSs where $ G_i, \; i = 1,2,\ldots,n, $ are all constant matrices and then extended the results to more general cases where $ G_i,\; i = 1,2,\ldots,n, $ could be matrix functions. A coordinate transformation has been introduced to obtain the HOFAS model. Based on the HOFAS models, the exponentially stabilizing controllers have been directly given. In fact, there is no obstacle for us to extend the results to discrete-time SFSs with increasing dimensions and SFSs with state delays and increasing dimensions. Related proofs can be established by combining the methods of this paper and the existing works. Our future work may concentrate on providing HOFAS models for SFSs with decreasing dimensions.

  • [1]
    M. Krstic and M. Bement, “Nonovershooting control of strict-feedback nonlinear systems,” IEEE Trans. Automatic Control, vol. 51, no. 12, pp. 1938–1943, 2006. doi: 10.1109/TAC.2006.886518
    [2]
    Z. Pan and T. Basar, “Adaptive controller design for tracking and disturbance attenuation in parametric strict-feedback nonlinear systems,” IEEE Trans. Automatic Control, vol. 43, no. 8, pp. 1066–1083, 1998. doi: 10.1109/9.704978
    [3]
    N. Bekiaris-Liberis and M. Krstic, “Stabilization of linear strict-feedback systems with delayed integrators,” Automatica, vol. 46, no. 11, pp. 1902–1910, 2010. doi: 10.1016/j.automatica.2010.07.008
    [4]
    W. Li, X. Yao, and M. Krstic, “Adaptive-gain observer-based stabilization of stochastic strict-feedback systems with sensor uncertainty,” Automatica, vol. 120, p. 109112, 2020. doi: 10.1016/j.automatica.2020.109112
    [5]
    Y. Li, Y. Fan, K. Li, W. Liu, and S. Tong, “Adaptive optimized backstepping control-based rl algorithm for stochastic nonlinear systems with state constraints and its application,” IEEE Trans. Cybernetics, vol. 52, no. 10, pp. 10542–10555, 2021.
    [6]
    P. V. Kokotovic, “The joy of feedback: Nonlinear and adaptive,” IEEE Control Systems Magazine, vol. 12, no. 3, pp. 7–17, 1992. doi: 10.1109/37.165507
    [7]
    M. Krstic, P. V. Kokotovic, and I. Kanellakopoulos, Nonlinear and Adaptive Control Design. John Wiley & Sons, Inc., 1995.
    [8]
    G. Meunier, B. Boulet, and M. Nahon, “Control of an overactuated cable-driven parallel mechanism for a radio telescope application,” IEEE Trans. Control Systems Technology, vol. 17, no. 5, pp. 1043–1054, 2009. doi: 10.1109/TCST.2008.2004812
    [9]
    R. Olfati-Saber, “Cascade normal forms for underactuated mechanical systems,” in Proc. 39th IEEE Conf. Decision and Control (Cat. No. 00CH37187), vol. 3, 2000, pp. 2162–2167.
    [10]
    G. Duan, “High-order fully actuated system approaches: Part I. Models and basic procedure,” Int. Journal of Systems Science, vol. 52, no. 2, pp. 422–435, 2021. doi: 10.1080/00207721.2020.1829167
    [11]
    G. Duan, “High-order fully actuated system approaches: Part II. Generalized strict-feedback systems,” Int. Journal of Systems Science, vol. 52, no. 3, pp. 437–454, 2021. doi: 10.1080/00207721.2020.1829168
    [12]
    G. Duan, “Fully actuated system approaches for continuous-time delay systems: Part 1. Systems with state delays only,” Science China Information Sciences, vol. 66, article 112201, 2023.
    [13]
    D.-W. Zhang, G.-P. Liu, and L. Cao, “Coordinated control of high-order fully actuated multiagent systems and its application: A predictive control strategy,” IEEE/ASME Trans. Mechatronics, vol. 27, no. 6, pp. 4362–4372, 2022. doi: 10.1109/TMECH.2022.3156587
    [14]
    Y. Yu, G.-P. Liu, Y. Huang, and J. M. Guerrero, “Coordinated predictive secondary control for dc microgrids based on high-order fully actuated system approaches,” IEEE Trans. Smart Grid, vol. 15, no. 1, pp. 19–33, 2024.
    [15]
    R. Meng, C. Hua, K. Li, and P. Ning, “Adaptive event-triggered control for uncertain high-order fully actuated system,” IEEE Trans. Circuits and Systems II: Express Briefs, vol. 69, no. 11, pp. 4438–4442, 2022.
    [16]
    R. Dong, C. Hua, K. Li, and R. Meng, “Adaptive fault-tolerant control for high-order fully actuated system with full-state constraints,” Journal of the Franklin Institute, vol. 360, no. 12, pp. 8062–8074, 2023.
    [17]
    M. Cai, X. He, and D. Zhou, “Fault-tolerant tracking control for nonlinear observer-extended high-order fully-actuated systems,” Journal of the Franklin Institute, vol. 360, no. 1, pp. 136–153, 2023. doi: 10.1016/j.jfranklin.2022.11.025
    [18]
    D. Gu and S. Wang, “A high-order fully actuated system approach for a class of nonlinear systems,” Journal of Systems Science and Complexity, vol. 35, no. 2, pp. 714–730, 2022. doi: 10.1007/s11424-022-2041-4
    [19]
    M. Cai, X. He, and D. Zhou, “Low-power fault-tolerant control for nonideal high-order fully actuated systems,” IEEE Trans. Systems, Man, and Cybernetics: Systems, vol. 53, no. 8, pp. 4875–4887, 2023.
    [20]
    D.-W. Zhang, G.-P. Liu, and L. Cao, “Proportional integral predictive control of high-order fully actuated networked multiagent systems with communication delays,” IEEE Trans. Systems, Man, and Cybernetics: Systems, vol. 53, no. 2, pp. 801–812, 2022.
    [21]
    L. Zhang, L. Zhu, and C. Hua, “Practical prescribed time control based on high-order fully actuated system approach for strong interconnected nonlinear systems,” Nonlinear Dynamics, vol. 110, no. 4, pp. 3535–3545, 2022. doi: 10.1007/s11071-022-07820-w
    [22]
    A.-G. Wu, J. Zhang, and Y. Ji, “A fully actuated system approach for stabilization of discrete-time multiple-input nonlinear systems with distinct input delays,” Journal of Systems Science and Complexity, vol. 35, no. 2, pp. 670–687, 2022. doi: 10.1007/s11424-022-2046-z
    [23]
    Q. Hu, B. Li, B. Xiao, and Y. Zhang, Control Allocation for Spacecraft Under Actuator Faults. Springer, 2021.
    [24]
    A. Witkowska and R. Śmierzchalski, “Adaptive dynamic control allocation for dynamic positioning of marine vessel based on backstepping method and sequential quadratic programming,” Ocean Engineering, vol. 163, pp. 570–582, 2018. doi: 10.1016/j.oceaneng.2018.05.061
    [25]
    D. Khimani and M. Patil, “Fault-tolerant sliding mode control for uncertain over-actuated systems,” IFAC Proceedings Volumes, vol. 47, no. 1, pp. 748–753, 2014. doi: 10.3182/20140313-3-IN-3024.00155
  • Related Articles

    [1]Yi Huang, Guo-Ping Liu, Yi Yu, Wenshan Hu. Constrained Networked Predictive Control for Nonlinear Systems Using a High-Order Fully Actuated System Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2025, 12(2): 478-480. doi: 10.1109/JAS.2024.124764
    [2]Xiang Chen, Yujuan Wang, Yongduan Song. Unifying Fixed Time and Prescribed Time Control for Strict-Feedback Nonlinear Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2025, 12(2): 347-355. doi: 10.1109/JAS.2024.124401
    [3]Yuhan Zhang, Zidong Wang, Lei Zou, Yun Chen, Guoping Lu. Ultimately Bounded Output Feedback Control for Networked Nonlinear Systems With Unreliable Communication Channel: A Buffer-Aided Strategy[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(7): 1566-1578. doi: 10.1109/JAS.2024.124314
    [4]Chi Ma, Dianbiao Dong. Finite-time Prescribed Performance Time-Varying Formation Control for Second-Order Multi-Agent Systems With Non-Strict Feedback Based on a Neural Network Observer[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(4): 1039-1050. doi: 10.1109/JAS.2023.123615
    [5]Haijing Wang, Jinzhu Peng, Fangfang Zhang, Yaonan Wang. High-Order Control Barrier Function-Based Safety Control of Constrained Robotic Systems: An Augmented Dynamics Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(12): 2487-2496. doi: 10.1109/JAS.2024.124524
    [6]Jin-Xi Zhang, Kai-Di Xu, Qing-Guo Wang. Prescribed Performance Tracking Control of Time-Delay Nonlinear Systems With Output Constraints[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(7): 1557-1565. doi: 10.1109/JAS.2023.123831
    [7]Zhuwu Shao, Yujuan Wang, Zeqiang Li, Yongduan Song. Dynamic Constraint-Driven Event-Triggered Control of Strict-Feedback Systems Without Max/Min Values on Irregular Constraints[J]. IEEE/CAA Journal of Automatica Sinica, 2024, 11(3): 569-580. doi: 10.1109/JAS.2023.123804
    [8]Maojiao Ye, Qing-Long Han, Lei Ding, Shengyuan Xu. Fully Distributed Nash Equilibrium Seeking for High-Order Players With Actuator Limitations[J]. IEEE/CAA Journal of Automatica Sinica, 2023, 10(6): 1434-1444. doi: 10.1109/JAS.2022.105983
    [9]Hao Liu, Qing-Long Han, Yuzhe Li. Set-Valued State Estimation of Nonlinear Discrete-Time Systems and Its Application to Attack Detection[J]. IEEE/CAA Journal of Automatica Sinica.
    [10]Chaoxu Mu, Yong Zhang, Guangbin Cai, Ruijun Liu, Changyin Sun. A Data-Based Feedback Relearning Algorithm for Uncertain Nonlinear Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2023, 10(5): 1288-1303. doi: 10.1109/JAS.2023.123186
    [11]Kai Zhao, Changyun Wen, Yongduan Song, Frank L. Lewis. Adaptive Uniform Performance Control of Strict-Feedback Nonlinear Systems With Time-Varying Control Gain[J]. IEEE/CAA Journal of Automatica Sinica, 2023, 10(2): 451-461. doi: 10.1109/JAS.2022.106064
    [12]Dong Zhao, Marios M. Polycarpou. Fault Accommodation for a Class of Nonlinear Uncertain Systems With Event-Triggered Input[J]. IEEE/CAA Journal of Automatica Sinica, 2022, 9(2): 235-245. doi: 10.1109/JAS.2021.1004314
    [13]Min Wang, Haotian Shi, Cong Wang. Distributed Cooperative Learning for Discrete-Time Strict-Feedback Multi Agent Systems Over Directed Graphs[J]. IEEE/CAA Journal of Automatica Sinica, 2022, 9(10): 1831-1844. doi: 10.1109/JAS.2022.105542
    [14]Ashish Kumar Jain, Shubhendu Bhasin. Tracking Control of Uncertain Nonlinear Systems With Unknown Constant Input Delay[J]. IEEE/CAA Journal of Automatica Sinica, 2020, 7(2): 420-425. doi: 10.1109/JAS.2019.1911807
    [15]Laura Menini, Corrado Possieri, Antonio Tornambè. Algorithms to Compute the Largest Invariant Set Contained in an Algebraic Set for Continuous-Time and Discrete-Time Nonlinear Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2020, 7(1): 57-69. doi: 10.1109/JAS.2019.1911819
    [16]Daniele Carnevale, Corrado Possieri, Antonio Tornambè. On "Over-Sized" High-Gain Practical Observers for Nonlinear Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2018, 5(3): 691-698. doi: 10.1109/JAS.2018.7511075
    [17]Xinping Guan, Bo Yang, Cailian Chen, Wenbin Dai, Yiyin Wang. A Comprehensive Overview of Cyber-Physical Systems: From Perspective of Feedback System[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(1): 1-14.
    [18]Shanwei Su. Output-feedback Dynamic Surface Control for a Class of Nonlinear Non-minimum Phase Systems[J]. IEEE/CAA Journal of Automatica Sinica, 2016, 3(1): 96-104.
    [19]Qiming Zhao, Hao Xu, Sarangapani Jagannathan. Near Optimal Output Feedback Control of Nonlinear Discrete-time Systems Based on Reinforcement Neural Network Learning[J]. IEEE/CAA Journal of Automatica Sinica, 2014, 1(4): 372-384.
    [20]Tengfei Liu, Zhongping Jiang. Distributed Control of Nonlinear Uncertain Systems: A Cyclic-small-gain Approach[J]. IEEE/CAA Journal of Automatica Sinica, 2014, 1(1): 46-53.

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(2)

    Article Metrics

    Article views (186) PDF downloads(71) Cited by()

    Highlights

    • We present HOFAS models specifically designed for SFSs with increasing dimensions
    • Various SFSs were transformed into FAS models
    • We tackle the control problem associated with SFSs exhibiting increasing dimensions
    • Simulation examples are given with both first and second order cases

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return