`Mathematical Problems in EngineeringVolume 2011, Article ID 678575, 25 pageshttp://dx.doi.org/10.1155/2011/678575`
Research Article

## Commutativity of Linear Time-Varying Differential Systems with Nonzero Initial Conditions: A Review and Some New Extensions

1Computer Engineering Department, Halic University, Siracevizler Street, no. 29, 34363 Sisli, Istanbul, Turkey
2Department of Elementary Mathematics Education, Mevlana University, 42003 Konya, Turkey

Received 18 January 2011; Revised 21 June 2011; Accepted 12 July 2011

Copyright © 2011 Muhammet Koksal and Mehmet Emir Koksal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Necessary and sufficient conditions for the commutativity of linear time-varying systems are derived in the case of nonzero initial conditions. It is shown that some commutative class of linear time-varying systems may not commute with arbitrary initial conditions. In this respect, commutativity of Euler differential systems is investigated. Explicit commutativity conditions for the fifth-order systems are solved. New results about the effects of commutativity on system sensitivity and disturbance properties are presented, which is very important for network design and industrial applications where many of the systems are composed of subsystems cooperating one after another in a chain. The results are supported by examples treated either analytically or numerically.

#### 1. Introduction

When two linear time-varying single-input single-output dynamical systems A and B are connected in cascade (or series), the input-output relation of the combined system depends on the parameters of both systems and on which appears first. If both of the connections AB and BA have the same input-output pairs irrespective of the applied input, then we say that these systems are commutative systems; in this case AB and BA are equivalent, that is, AB = BA.

Most of the industrial processes and network designs contain subsystems or subnetworks which are arranged in a cascaded fashion or in a chain, which can be called a serial composition as well. If the change of the order of any two subunits in this chain does not affect the input-output relation of the whole system while it is creating some advantages on the sensitivity, stability, disturbance or noise, or other similar properties, the use of the optimum sequence is inevitable. Therefore, the study of the commutativity conditions is very important from theoretical and as well as practical point of views. And, this paper is aimed at closing the gap in the theory on the related topic.

An exhaustive study of the commutativity of the linear time-varying dynamical systems which can be referred as differential systems, due to their mathematical representations given in Section 2, was introduced in 1988 by the author [1]. This work has been the basic reference for years since it includes the most general necessary and sufficient commutativity conditions for systems of any order but without initial conditions. All the previous results concerning the commutativity conditions were shown to be deduced from the main theorem in that paper. Namely, those for the first-order [2], second-order [35], third-order [6], and fourth-order [7, 8] systems.

In addition to the deductions of the commutativity conditions for low-order systems, the main reference [1] covers some general results in the form of corollaries or theorems of minor importance. For example, Corollary 3.2 states that for the commutativity of two linear systems it is required that either both systems are time-invariant or both systems are time-varying, which is a result originally stated by Marshall [2]. Corollary 3.3 states the commutativity of identical time-varying systems with arbitrary time-invariant forward and feedback path gains, which is originally proved by the author himself for second-order systems [5]. Theorem 5.1 deals with the unconditional commutativity property of Euler Systems in the light of the stated general necessary and sufficient conditions, which is originally treated in the undistributed report [7] and presented in a national conference [8].

Although the exhaustive study [1] included all the previous results reported until its publication, it was not claimed to be the completion of the research on the subject. On the contrary, it forecasted open questions on the commutativity conditions for systems with initial conditions and on the role of the system performance characteristics such as noise, sensitivity, and disturbance, in the cascade connected control or communication networks. Later, some new results on these subjects have appeared in the national [9] and international [10, 11] conference proceedings.

The commutativity conditions with nonzero initial conditions were derived in [9, 10]. It was shown that two initially relaxed commutative systems may not be commutative when they have nonzero initial conditions; further, commutativity with arbitrary initial conditions requires additional constraints on the coefficients of the differential equations describing the systems. The importance of commutativity as far as the system sensitivity is concerned was reported in [9, 11]. It was shown that for the same input-output relation, one system could be much less sensitive to system parameters than the other if the subsystems are commuted.

Although the above results are important from both theoretical and application point of views, it is believed that they have not been widely and sufficiently announced at least in an international journal paper. This is due to time limitations of the presentations in the conferences and page limitations in their proceedings which always sound less than journals. Further, some errors in the mentioned literature of major importance but minor reflection need to be corrected. Therefore, the scope of this paper is confined mainly on the explicit commutativity conditions for linear time-varying differential systems with nonzero initial conditions and the effects of commutativity on the system sensitivity. Further, the analytic solution of explicit commutativity conditions for the fifth-order systems is introduced in the literature for the first time. Another original point is the derivation of explicit commutativity conditions for Euler systems which are always commutative without initial conditions. Commutativity of time-invariant linear systems with nonzero initial conditions is also studied as another particular case. Further, new results are obtained for the commutativity of feedback systems with nonzero initial conditions.

This paper, as a survey or a tutorial on the commutativity with nonzero initial conditions, is mainly aimed at presenting the above-mentioned original items of theoretical values and moreover to highlight the importance of commutativity on industrial applications or engineering designs where most of the subsystems are connected in cascade or in a chain.

The rest of the paper is organized as follows: in the following section, the problem of commutativity is defined and for the sake of completeness the main results previously obtained in [1] are summarized as Theorem 2.1 without proof. In Section 3, the commutativity conditions for systems with initial conditions are introduced by Theorem 3.1; some results of minor importance in the above-mentioned literature are correctly stated. Particular results including feedback systems are presented. Section 4 includes the explicit solution of the commutativity conditions for the fifth-order systems. Examples showing the roles of commutativity on system disturbance and system sensitivity to parameter values are presented in Section 6 to show the practical importance of the results. Finally the last section is devoted to conclusions.

#### 2. First Commutativity Conditions

Let the two linear time-varying dynamical systems A and B be represented by th order linear differential equations with variable coefficients , , , respectively, as follows: where and are piecewise continuous mappings from into which is the real line; and are the input and output of the systems and . It is assumed that is also piecewise continuous. Although , that is, A is always assumed an nth order dynamical system, () it is allowed that for , , whilst , which means B is an th order system; is assumed to be equal to unless otherwise stated. It is well known [12] that both systems and hence their interconnections described in the sequel have unique solutions under the cited conditions about the coefficients and unless and , where is the initial time, which is a fact that will not be emphasized in this paper anymore.

For the cascade connection AB shown in Figure 1(a), the output of A is equal to the input of B; so that letting and replacing (2.1b) into (2.1a), the following differential equation is obtained where and are the input and output of the cascaded system AB, respectively; the superscript in the parenthesis denotes ith derivative. Note that B is assumed of degree . Performing the th derivative on the product and rearranging, (2.2) reduces to the following ()th order differential equation describing AB:

Figure 1: Cascade connection of the differential systems described by (2.1a) and (2.1b).

For the connection BA shown in Figure 1(b), the input , , the output , and a similar equation to (2.3a) are obtained with the replacements of , , , by , , , , respectively. The result is

When the two differential equations in (2.3a), (2.3b) are identical except at a finite number of points in every finite subinterval of (these may be the discontinuity points of the coefficients and and their derivatives), the connections AB and BA will have the same output for the same input x if the initial conditions for subsystems A and B are zero. When the coefficients of the derivatives of the same orders are equated, the necessary and sufficient conditions expressed by the following theorem can be obtained [1].

Theorem 2.1 (Koksal 1). The necessary and sufficient conditions for a linear time-varying system of order n (System A) to be commutative with another linear time-varying system of the same or lower order (System B), in the case of zero initial conditions, are that
(i)the coefficients of the differential equation describing the second system must be expressed in terms of the coefficients of the differential equation describing the first system through equations involving properly chosen constants , where (ii)further, with the same constants ’s, the coefficients of the first system must satisfy nonlinear algebraic equations which involve the derivatives of the coefficients as well due to the expression

Although this theorem is apparently stated and originally proved for , that is, for the commutativity of equal order systems, it should be remarked that it can easily be reduced for the case by simply choosing , which results with . Hence, it is valid for systems of any order. Another remark to be noted is that a decreasing order indexing is used for convenience in all of the subscripts and summations in (2.4a) and (2.5a) in accordance with the order of computations from simple to complex; this order corresponds to the decreasing order of derivatives and coefficients in the original equations ((2.1a), (2.1b)).

The above theorem expresses the necessary and sufficient conditions for systems A and B in the case of zero initial conditions for the differential systems in ((2.1a), (2.1b)). Although the initial clue for the proof is given above, the full proof is exhaustive and given in [1] which also discusses the consequences of the theorem. Therefore, it is not included in this paper. Instead, the case of nonzero initial conditions is considered next.

#### 3. Second Commutativity Conditions

Let the systems in ((2.1a), (2.1b)) be commutative when they are initially relaxed, that is, they satisfy the conditions of Theorem 2.1. Consider now the nonzero initial conditions expressed bywhere and express the vector of the first and the remaining n-m initial values of and it is derivatives according to the partitioning in (3.1a). With these initial conditions, the commutativity requirements are expressed by the following theorem.

Theorem 3.1 (Koksal 2). The necessary and sufficient conditions for the commutativity of two linear time-varying systems A and B which are described by the differential equations in (2.1a) and (2.1b), and with the initial conditions in (3.1a) and (3.1b), respectively, are (i)the conditions of Theorem 2.1 on the coefficients (2.4a), (2.5a)) be satisfied,(ii)the initial conditions satisfy the following at the initial time : or equivalently, where the matrices A1, A2, B1, B2 are defined in the proof (see (3.8a)–(3.8d)).

Proof. Having satisfied the conditions of Theorem 2.1, the differential equations in (2.3a) and (2.3b) will have the same coefficients. For the same input, they will have identical solutions if their initial conditions are also equal.
For the connection AB, , and for BA, ; hence to have the same initial conditions for in (2.3a) and (2.3b), must be satisfied at the initial time . On the other hand, using for AB and for BA in (2.1b) and (2.1a),are obtained for the connections AB and BA, respectively. Using the formula the kth derivative of a product can be obtained from (3.4a) and (3.4b), respectively.
The last two equations can be written in matrix form aswhere the entries of the coefficient matrix are defined byThe above expressions for the entries of the coefficient matrices in ((3.7a), (3.7b) imply that , , , are matrices of orders , , , , respectively. Further, the square matrices and are lower triangular matrices with nonzero diagonal elements and , respectively; therefore, they have well-defined inverses. Hence, for can be computed from (3.7b), and for can be computed from (3.7a). If these computed values are inserted in (3.3), this equation can be rearranged in matrix from as in (3.2a) or equivalently (3.2b).
The linear time varying systems A and B are commutative for all inputs if and only if (1) AB and BA are described by identical differential equations and (2) this differential equation has the same initial conditions for AB and BA at the initial time . Hence, in addition to the conditions of Theorem 2.1, the conditions obtained in ((3.2a), (3.2b)) for the equivalent initial conditions of (2.3a) and (2.3b) constitute the necessary and sufficient conditions of the commutatitivity of the systems A and B. Equation (2.3b) must be satisfied by the initial conditions of A and B. Hence this completes the proof.

This theorem is stated improperly in the present literature [9, 10], which formulated that the second condition above as the coefficient matrix in (3.2b) be identically zero, which is an impossible situation unless and system B is an identity.

We now consider the consequences of Theorem 3.1 in the form of corollaries. First consider the commutativity of A and B for all arbitrary initial conditions. This is a case which we call absolute commutativity. Even the first condition about the coefficients is satisfied, the commutativity for all arbitrary initial conditions is not possible if . This is obvious from the first block of equations in (3.2a). If , ((3.2a), (3.2b)) imply that ; that is, for the commutativity of all initial conditions. Equation (3.8d) gives that ; hence , which means that system B is an identity. This result can be expressed by the following corollary.

Corollary 3.2 (Absolute Commutativity). There are no time-varying dynamical systems which are absolutely commutative at any time . The only absolutely commutative system of a time-varying system of order is the identity. Note that the corollary holds for time-invariant systems as well.

When the case of the special values of the initial conditions is considered for A and B at a particular initial time (conditional commutativity), the first rows and the remaining rows of (3.2b) can be expressed as follows:

Corollary 3.3 (Conditional Commutativity). If the systems A and B of Theorem 3.1 are commutative with zero initial conditions, the necessary and sufficient conditions for their commutativity under nonzero initial conditions at time t0 as well are (i)the first m initial values of A expressed in (3.1a) must be identical with those of B expressed in (3.1b);(ii)the initial condition vector of A, that is, in (3.1a) must be in the null space of the value of the coefficient matrix in (3.9b) at the initial time to .

Corollary 3.4 (Conditional Commutativity of Systems of Equal Orders). If A and B are commutative systems with the same order , for their commutativity with nonzero initial conditions at the initial time , it is necessary and sufficient that both of these systems must have the same initial condition vectors in the null space of matrix at .

This corollary follows from (3.9a), (3.9b)) with .

#### 4. Reduction for Special Cases

##### 4.1. Time-Invariant Systems with Nonzero Initial Conditions

If the coefficients of linear systems in (2.1a) and (2.1b) are constant, A and B become time-invariant. Since time-invariant linear systems with zero initial conditions are commutative, they commute [1]. Although primitive, this can be proved by considering the Laplace transform of (2.1a) and (2.1b) with zero initial conditions and showing that AB and BA have the same transfer function where and are the Laplace Transforms of the input and the output of the connections AB and BA. An alternate way would be to show that the conditions of Theorem 2.1 are satisfied. This necessitates the reduction of the formulas in (2.4a), (2.4b) and (2.5a), (2.5b)) for constant ’s. Since the fact to be shown is obvious from the above discussion, this approach is omitted.

The commutativity of time-invariant systems with nonzero initial conditions, however is not straightforward. In this case, Corollaries 3.3 and 3.4 are valid for all initial times , which is due to the time invariance of the coefficients. For the last reason, the expressions in (3.8a)–(3.8d) for the entries of the coefficient matrices in ((3.2a), (3.2b)), ((3.7a), (3.7b), (3.9a), (3.9b)) simplify to All the other entries of , , , which are not defined above are zero. Therefore, in addition to the lower triangular matrices and , and become upper triangular (trapezoidal if ).

##### 4.2. Euler Systems with Nonzero Initial Conditions

It is well known that linear time-varying systems described by Euler-type differential equations [13] are always commutative when the initial conditions are zero [1, 7, 8]. A detailed proof in [1] is given based on the impulse response representations of the individual systems as well as that of the interconnection. That paper also shows that the cascade connection of any two Euler systems is another Euler system, and further, all the commutative pairs of an Euler system are of the same type, that is, Euler-type. Although the proof can be done by showing that any Euler system A of order and another one B of the same or lower order satisfy the conditions of Theorem 2.1, the attempt here will be directly to the case of initial conditions.

It is shown by a simple example in [10] that the commutativity of Euler systems with arbitrary nonzero initial conditions is not straightforward as in the case of time-invariant systems. Here in the following, the problem is treated on a general base as an application of the results obtained in the previous section.

The differential systems A and B described by (2.1a) and (2.1b), respectively, become Euler-type systems if the time-varying coefficients are set to Using the formula for the th derivative of the th power, that is, all the previous results about the commutativity of Euler systems with nonzero initial conditions, in particular Theorem 3.1, Corollaries 3.3 and 3.4, can be derived with slight modifications. Since , only for , the initial time can be any real number except 0. The mentioned modifications are on the expressions in (3.8a)–(3.8d) for the entries of the coefficient matrices , , , which are used in the statement of the theorem and its mentioned corollaries. The Euler systems (3.8a)–(3.8d) reduce to the following:Note that the entries , , , of the matrices , , , of orders , , , , which are not defined in (4.5a),(4.5b), (4.5c), (4.5d), respectively, are all zero. Hence as in the time-invariant case, and are upper triangular.

Examples to explicit forms of the coefficient matrices used to express the constraints on the initial conditions for commutativity will be given in Section 6.

##### 4.3. Feedback Systems

Due to their importance in automatic control theory, commutativity of feedback systems is expressed by a corollary.

Corollary 4.1. Any two versions A and B of a linear time-varying system S simply obtained by constant forward and feedback path gains and , and , respectively, are commutative with any set of initial conditions of S if and only if where

Proof. Let the system be described by Then the feedback systems A and B shown in Figure 2 are described by (2.1a) and (2.1b), respectively, with Since these coefficients can be related as in (2.4a) withand (2.5a) is satisfied due to (4.9b), conditions of Theorem 2.1 are satisfied. Hence, A and B are always commutative with zero initial conditions.
For the case of nonzero initial conditions, first note that so that A and B have the same initial conditions as S; so Corollary 3.4 is satisfied with . Then express bi’s in terms of ai’s by using (4.8) and (4.9a), (4.9b) as follows: Using these expressions in (3.8c), (3.8d) to evaluate the coefficient matrices , appearing in (3.10), the following results are obtained: Finally inserting these expressions in (3.10), is obtained. Since is nonsingular and , , the common initial condition vector can be arbitrary chosen if the scalar coefficient of in (4.12) is zero; that is, Otherwise, A and B cannot commute with any nonzero initial conditions; hence this completes the proof.

Figure 2: Feedback systems obtained from S by constant forward and feedback path gains.

#### 5. Explicit Formulas for the Fifth-Order Systems

The explicit commutativity conditions are important from two points of views. The first is the ease for the applicability of the use of the commutativity conditions stated in Theorems 2.1 and 3.1. The second is their contribution on the development of the commutativity theory. In fact, the starting point of the commutativity theory considered in the sense of this paper is the explicit commutativity conditions suggested by Marshall for the first-order systems [2]. Later the explicit commutativity conditions for second-order systems [35] and those for third-order systems [6] constituted the base of the general commutativity conditions originally presented in [6]. Later the explicit conditions of commutativity for forth-order systems appeared in the undistributed work [7] and presented in [8]. All of these results were summarized in the main reference [1].

Here as in the following theorem, the explicit commutativity conditions for fifth-order linear time-varying systems are presented for the first time.

Theorem 5.1 (Koksal 3). The necessary and sufficient conditions for a fifth-order linear time-varying system of order n (System A) to be commutative with another linear time-varying system of the fifth- or lower-order (System B), in the case of zero initial conditions, are that
(i)the coefficients of the differential equation describing the second system must be expressed in terms of the coefficients of the differential equation describing the first system through 6 equations involving 6 properly chosen constants , , , , , , where (ii)further, with the same constants ’s, the coefficients of the first system must satisfy 4 nonlinear algebraic equations where

Proof. When the general conditions of Theorem 2.1 are applied for the commutativity of fifth-order systems, (5.1a) results simply as the matrix form of (2.4a) for . The explicit solutions given in (5.1b) are obtained by evaluating the integrals in (2.4b) with , for each successively in the order given in (5.1b). However, it is very exhaustive to derive the explicit expressions for given in this equation. The derivation details are wearying and too lengthy to include in the proof, though they involve an attentive treatment of the ordinary integration methods. But the verification is straight foreword by inserting them in (2.4b) for . Finally (5.2a), (5.2b) follows as the direct result of the second condition (2.5a), (2.5b)) of Theorem 2.1 for and it is included for the sake of completeness of the theorem.
The above equations in (5.1b) explicitly represent the coefficients ‘s of the same or lower order commutative systems in terms of the coefficients ‘s of a fifth-order system. Using these equations in (3.8a)–(3.8d) and inserting the results in ((3.2a), (3.2b)), ((3.7a), (3.7b), (3.9a), (3.9b)), the commutativity conditions for systems with initial conditions are also expressed explicitly. Although it appears that these explicit results are inevitably lengthy, ((5.1a), (5.1b) and all the subsequent equations obtained from them greatly simplify if some of the coefficients ’s are constant or completely zero.
Another remark about shortening all the apparent results (not only for fifth-order systems) obtained in this paper is the assumption ; that is, the normalization of the coefficients ’s dividing by . In this case, all the derivatives of will be zero, and all its integer or fractional powers will be 1. But this is experienced complicating the expressions for the derivatives of the remaining coefficients due to the existence of the denominator term generated by normalization.

#### 6. Effects of Commutativity on System Sensitivity and Disturbance

The importance of commutativity for system sensitivity was first reported by the author [9, 11]. It was shown that when the order of connections of a second-order system and its first-order commutative pair was changed, the impulse response and hence the input-output transfer characteristics could be much more stable to parameter changes in one case than the other. This fact is very important in systems which use cascade connection of subsystems, which is usually the case in most of the telecommunication, control, and other practical applications. The effect of noise or disturbance interfering at an intermediate stage of a system may also depend on the order of connections of the subsystems. Therefore, the convenient order must be preferred when cascading the systems which need to be commutative to function the same input-output relation.

To illustrate the use of the derived formulas for the commutativity of fifth-order systems, the first example is chosen accordingly for the study of sensitivity and disturbance in conjunction with the commutativity.

Example 6.1. Let A be a fifth-order system described by Consider the search for the fifth- or lower-order commutative pairs of A. Computing ’s expressed in (5.1b) and inserting them in (5.1a), the following matrix equation is obtained for the fifth- or lower-order commutative systems B: In view of (5.2a), (5.2b), the second commutativity condition of Theorem 2.1 implies where ’s are computed using (2.5b). Note that the entries above the second cross diagonal are not computed since they are not needed to arrive the conclusion that . With this information, all the commutative pairs of A are described by Since , this case is obviously a special application of Corollary  3 in [1]. In fact, system B described by (6.4) has a block-diagram representation shown in Figure 3; that is, B is obtained from A by an arbitrary forward and feedback path gains and , respectively.
Consider now the commutativity with nonzero initial conditions. For the system A in (6.1) and the system B in (6.4), and Corollary 3.3 necessitates A and B have the same initial conditions in the null space of which is the coefficient matrix in (3.9b).
Obviously, the null space is empty for all initial times unless , and for , A and B are commutative for any set of initial conditions.
The above conclusions are verified by Simulink. The results are plotted in Figure 4. In this figure, the curve I shows the equal responses of AB and BA to a unit step for , with zero initial conditions. For curves II-AB and BA, equal nonzero initial conditions (, , , all others are zero) are used with the same values of input and c0, c5; obviously A and B are not commutative since . In Figure 4(b), I shows the step response of both AB and BA with the equal initial conditions when , , that is, ; hence the commutativity conditions hold. When the same coefficients but unequal initial conditions are used ( is changed to 1), AB and BA yield different responses as shown by the curves II-AB, BA.
To show the sensitivities of AB and BA to parameter values and initial conditions, consider the commutative case , , , , for which the unit step response was shown in Figure 4(b) (I). If the commutativity condition is spoiled slightly by changing to −1.2, then the responses of AB and BA change as shown in Figure 4(c) by the curves I-AB, BA. Obviously, the connection AB is much more sensitive to than BA. On the other hand, when is reduced to 0.60, then AB and BA give the responses shown in Figure 4(c) by II-AB, BA; AB is still more sensitive than BA for this case.
To test the disturbance by an external noise, a unit step is injected in to the connection between A and B when both systems are initially relaxed and unforced with , . Under the same conditions, a saw-tooth wave with period 0.2 and ranging from –18.5 to 20 is applied as the second disturbance. The results are shown in Figure 4(d) by the curves I-AB, BA and II-AB, BA, respectively. Comparing the results with the zero response, it is seen that the output of AB is much more effected than that of BA by both of the disturbances.

Figure 3: The only commutative class of systems of A considered in Example 6.1.
Figure 4: Simulation results for Example 6.1.

Example 6.2. Consider the second-order time-invariant linear systems where A is chosen unstable to emphasize on the effects of commutativity on sensitivity and disturbance. For the commutativity of these systems with nonzero initial conditions as well, Corollary 3.4 requires The coefficient matrix in (6.8) is the one given in (3.10) and obtained from it by using formulas in (4.2) for the given systems in this example. To have commutativity with nonzero initial conditions, it is necessary that this coefficient matrix must be singular which is achieved if Further, (6.8) requires Obviously, if , then and the systems are commutative for any set of equal initial conditions; that is, and can be chosen arbitrarily. If , using (6.9), (6.10) is reduced to
It is important to note that all of these commutativity results can be equally obtained by taking the Laplace transform of (6.6) with the general nonzero initial conditions for both systems A and B; and then computing and equating the -domain input-output relations for the connections AB and BA. The result is where represent the initial values at time . Equating the appropriate coefficients of and , results equivalent to the commutativity conditions in (6.7), (6.8) are obtained. However, this elementary process is not valid for time varying systems for which Laplace transform can apply only under specific conditions, and it is not as systematic and easily applicable as given in this paper for higher-order time-invariant systems.
The simulation results for the case , with a unit step input are shown in Figure 5(a). For and , , both connections yield the response I which is coincident with the plot of the analytic result computed from (6.12) or (6.13) by setting and inserting the given numerical values in; for and , the responses for AB and BA are still equal as noted by II; when b1 is changed to –2.5 in case II, the commutativity condition in (6.9) is spoiled and different responses are obtained as indicated by III-AB and III-BA for AB and BA. Further results are shown in Figure 5 for the case ; with , (6.9) implies ; and with , (6.11) implies . For and , , it can be shown by using (6.12) or (6.13) that an input cancels the pole in (6.12), (6.13)) at , and hence it prevents unboundedness and produces an output This analytical solution is shown in Figure 5(b) by I; in the same figure, II indicates the equal outputs of both AB and BA obtained by Simulink. A slight difference originates from the deficiency of Simulink due to the particular case in concern, that is, exact pole cancellation cannot be achieved numerically, so that simulation results get unbounded as . To test the sensitivity to the initial conditions, is changed to 1.01; the result is that the output of BA blows up almost two times before that of AB, as seen from the curves III-AB and III-BA in Figure 5(b), respectively; hence BA is more sensitive then AB to . In Figure 5(c), the responses to are compared for AB and BA; I is obtained for AB and BA with the same commutativity conditions valid for Figure 5(b) (I and II); II-AB and II-BA are obtained for AB and BA when is increased from to 10. Obviously, AB is almost insensitive to this change and BA is still more sensitive then AB. To show the disturbance effects, under the zero input conditions of the case in Figure 5(b) (I, II), a step of magnitude of 0.1 is injected at time between A and B. The output without disturbance (I-AB, BA) and the output with disturbance (II-AB, II-BA) are shown in Figure 5(d); obviously, AB is almost insensitive to this disturbance and BA is much more affected then AB.

Figure 5: Simulation and analytical results for Example 6.2.

Example 6.3. This example is preciously chosen so that the systems are both time varying and they have analytic solutions. A is defined as The coefficients of this system satisfy the first set of commutativity conditions expressed by (2.5a), and with (2.4a) all of its commutative pairs are expressed by Further, if A and B are required to be commutative with nonzero initial conditions as well, Corollary 3.3 requiresat the initial time . For the commutativity with nonzero initial conditions, the coefficient matrix in (6.18b) must be singular, that is, its determinant must be zero and this is achieved if and only if With this condition, (6.18b) implies the relation between equal initial conditions of systems A and B, so that they are commutative in case of nonzero initial conditions as well.
For the analytic solution of differential equations describing systems A and B, these equations are first written in state-space form, then the successive decomposition method is used to decompose the state-coefficient -matrix [13, 14]. A second decomposition is used to evaluate the state transition matrix [15]. Knowing the state transition matrices, the general solutions with arbitrary initial conditions at are obtained for systems A and B by the usual methods. The computed results are expressed for as where . Obviously both systems are strongly stable due to the terms in the exponents.
Using the input-output relations for the connections AB and BA depicted in Figure 1, the analytic solutions are obtained for any input and any set of initial conditions. The derivation details and the results are not given here for space limitations. Comparison of the like terms in the solutions reveals that the responses with zero initial conditions are equal. Further, to have equivalence with nonzero initial conditions, it is verified that (6.19a), (6.19b) must hold, and A and B must have identical initial conditions satisfying (6.19b).
The simulation results and some analytical solutions for the example are shown in Figure 6. In this figure, equal initial conditions are used for A and B in all cases; it is set that , , as to satisfy the commutativity conditions. In Figure 6(a), the input is 100 sin t; I is for , so that (6.19b) is satisfied; hence AB and BA give the same response. However, when is set to 0, then (6.19b) is spoiled and AB and BA give different responses as shown by the curves II-AB and II-BA. Observe that the connection BA is more sensitive to this change than AB. In Figure 6(b), the input is a square wave of magnitude 100 and period 1, , ; since all the commutativity conditions hold AB, and BA yield the same responses I-AB, BA; however, when the coefficient in (6.17b) is omitted, the first commutativity condition gets violated and the connections give different responses indicated by II-AB and II-BA. To see the disturbance (or noise) effects, the input is set to 10, , . Commutativity conditions hold and both connections yield the response shown in Figure 6(c) by the curve I-AB, BA. But when a saw tooth disturbance signal (with period 0.5, maximum and minimum values ) is injected between the connections, BA is much more affected than AB as the responses II-AB and II-BA show. Finally in Figure 6(d), the curve I shows the output (scale 100 : 1) of AB and BA with zero input and equal initial conditions , ; the results for AB, BA, and the analytic solution are all coincident. When is increased by an amount of , the deviations and analytically expressed by and the simulation results are plotted for (curves II-AB, II-BA). Analytical and simulation results are checked to be coincident. Obviously, the connection BA is more sensitive than AB to the initial condition relevant to the derivative.

Figure 6: Simulation and analytical results for Example 6.3.

#### 7. Conclusions

Commutativity conditions of two linear time-varying systems with nonzero initial conditions are thoroughly investigated in this paper. The derived conditions can be used to test the commutativity without solving the differential equations describing the system. Although equivalent conditions can be obtained in another way comparing the analytic solutions of the interconnections, very few classes of time-varying systems have explicit analytic solutions only. Therefore, the method of finding explicit solutions is not general and easily applicable as the method presented in this paper since it does not require finding the analytic solutions. This method can be used for all linear time varying systems as well.

The explicit commutativity conditions for fifth-order systems are obtained by an enormous effort and presented in this paper. Moreover, the commutativity is studied in view of system sensitivity and disturbance; it is shown that for commutative pairs, one order of connection might be more advantages than the other, which is an important fact from the applications point of view.

Most of the results are verified by examples treated analytically and/or by Simulink. It is seen that the analytical solutions and the results of Simulink are almost same with an exception where there is a perfect pole cancellation which cannot be achieved numerically by Simulink.

#### References

1. M. Koksal, “An exhaustive study on the commutativity of time varying systems,” International Journal of Control, vol. 47, no. 5, pp. 1521–1537, 1988.
2. E. Marshall, “Commutativity of time varying systems,” Electronics Letters, vol. 13, no. 18, pp. 539–540, 1977.
3. M. Koksal, “Commutativity of second order time varying systems,” International Journal of Control, vol. 36, no. 3, pp. 541–544, 1982.
4. S. V. Saleh, “Comments on ‘Commutativity of second-order time-varying systems’,” International Journal of Control, vol. 37, pp. 1195–1 196, 1983.
5. M. Koksal, “Corrections on ‘commutativity of second order time varying systems’,” International Journal of Control, vol. 38, no. 1, pp. 273–274, 1983.
6. M. Koksal, “General conditions for commutativity of time-varying systems,” in Proceedings of the (IASTED) International Association of Science and Technology for Development Conference, Telecommunication and Control (TELECON '84), pp. 225–229, Halkidiki, Greece, 1984.
7. M Koksal, A Survey on the Commutativity of Time-Varying Systems, Middle East Technical University, Gaziantep, Turkey, 1985.
8. M. Koksal, “Commutativity of 4th order systems and Euler systems,” in Proceedings of the National Congress of Electrical Engineering, Adana, Turkey, 1985, Paper no:BI-6.
9. M. Koksal, “Effects of initial conditions on the commutativity and effects of commutativity on the system sensitivity,” in Proceedings of the 2nd National Electrical Engineering Conference, pp. 570–573, Ankara, Turkey, 1987.
10. M. Koksal, “Effects of nonzero initial conditions on the commutativity of linear time-varying systems,” in Proceedings of the International Conference on Modelling and Simulation (AMSE '88), vol. 1, pp. 49–55, AMSE Press, Istanbul, Turkey, 1988.
11. M. Koksal, “Effects of commutativity on system sensitivity,” in Proceedings of the 6th International Symposium on Networks, Systems and Signal Processing (ISYNT2 '89), pp. 61–62, Zagreb, Yogoslavia, 1989.
12. C. A. Desoer, Notes for a Second Course on Linear Systems, Van Nostrand Rheinhold, New York, NY, USA, 1970.
13. M. Y. Wu, “A successive decomposition method for the solution of linear time-varying systems,” International Journal of Control, vol. 33, no. 1, pp. 181–186, 1981.
14. M. Y. Wu, “On the commutative class of linear time-varying systems,” International Journal of Control, vol. 23, no. 3, pp. 433–444, 1976.
15. M. Y. Wu and M. Horowitz, “On solution, stability and transformation of linear time-varying systems,” International Journal of Control, vol. 22, no. 2, pp. 169–180, 1975.