The relation between phase-type representation and positive system realization in both the discrete and continuous time is discussed. Using the Perron-Frobenius theorem of nonnegative matrix theory, a transformation from positive realization to phase-type realization is derived under the excitability condition. In order to explain the connection, some useful properties and characteristics such as irreducibility, excitability, transparency, and order reduction for positive realization and phase-type representation are discussed. In addition, the connection between the phase-type renewal process and the feedback positive system is discussed in the stabilization concept.
1. Introduction
Positive system problems have been developed in applications areas such as biological models, production systems, and economic applications. The realization problem for positive system has extensively been considered in many research papers as [1–5]. Many research activities and applications have been devoted to the field of phase-type distributions [6–8]. The positive representation problem finding a Markov chain associated with phase-type distribution has considerable connections with the positive realization problem in the control theory [9].
We will discuss the relationship between phase-type representation and positive realization by using the Perron-Frobenius theorem introduced in [10–13]. The Perron-Frobenius theorem is an important concept for the study of positive systems. For example, the Perron-Frobenius theorem can be used to derive a transformation from positive realization to phase-type representation. We use tools and results within the broad research area of nonnegative matrix theory, which enable us to explore the characteristics and properties of a Metzler matrix and nonnegative matrix. Metzler matrices are replaced by Z-matrices, in particular, by the M-matrices introduced in [10–13].
The connection between phase-type and positive realization has restrictively been proved in irreducible representation cases by remarking that it can be easily simplified to an irreducible case by discarding some states [9]. Under the irreducible assumption, it is proven that the positive realization can be transformed into a phase-type representation [9]. We will show that a positive realization normalized by a positive number can be transformed into a phase-type representation. We modify the correspondence between positive realizations and phase-type representations under more general assumptions. We use excitable systems as a subclass of the positive systems introduced in [2, 14]. However, the phase-type representation has a benefit that the number of free parameters in the representation can be reduced, compared with the general positive realization.
We will discuss the properties and characteristics, such as irreducibility, excitability, transparency, and stabilization introduced in [3, 8, 15–18]. Excitability and transparency are similar to the reachability and observability of positive linear systems [17, 18]. There exist unreachable and unobservable positive states that are excitable and transparent [2, 14]. The properties of excitability and transparency are discussed furthermore. We will demonstrate how to discard some unnecessary states when a representation is not irreducible. The common connection between the phase-type renewal process introduced in [7, 8] and the state feedback control introduced in [15, 16] will be handled. In addition, the relation and the common characteristics between discrete phase-type distributions and discrete-time positive systems are discussed in a similar manner to that applied to the continuous case.
An outline of the paper is as follows. In Section 2, some relevant background materials, including definitions and preliminary results are provided. Section 3 addresses the relationship between phase-type distributions and positive realizations in the continuous-time domain and discusses the common properties and characteristics, such as irreducibility, excitability, transparency, and stabilization. Section 4 discusses the relationship and the common characteristics between discrete phase-type and discrete-time positive systems in a similar manner to the continuous case. Finally, the conclusions follow.
2. Preliminaries
Before proceeding, we introduce some basic notations. An n×n nonnegative matrix A is denoted by A≥0 if its entries are nonnegative and at least one entry is positive. A=[aij] is defined by a strict positive matrix (i.e., A>0) if all entries aij>0. The associated directed graph, G(A), consists of n vertices {v1,…,vn}, where aij≠0 denotes an edge from vi to vj. A nonnegative matrix A is said to be reducible if there exists a permutation matrix P such that (1)PAPT=A11A120A22,where A11 and A22 are square matrices and PT is defined by the transpose of P. Otherwise, A is called irreducible. It is called a Frobenius normal form of A if there exists a suitable permutation P such that PAPT is in block triangular of the form(2)PAPT=A11A12⋯A1p0A22⋯A2p⋮⋮⋱⋮00⋯App,where Aii is square and either irreducible or a 1×1 null matrix. The spectral radius of A, denoted by ρ(A), is defined by the largest absolute eigenvalue of A. The spectrum of A, denoted by σ(A), is defined by the set of eigenvalues of A. The dominant eigenvalue of A is called by the maximal real among the real eigenvalues of A, denoted by λmax(A).
A matrix A∈Rn×n is said to be a Metzler matrix if all of its off-diagonal elements are nonnegative. If A is a Metzler matrix, then there exist a nonnegative matrix B≥0 and some α>0 such that A=B-αI. The real dominant eigenvalue of A is defined by λmax(A)=ρ(B)-α if λmax(A)≥Re(λ) for all λ∈σ(A). There is a long stream of research dealing with M-matrices and Z-matrices instead of Metzler matrices [10]. Metzler matrices are replaced by Z-matrices; that is, -A is a Z-matrix if A is a Metzler matrix. In particular, a matrix A is called an M-matrix if any matrix A is expressed in the form A=ηI-B where B≥0 and η≥ρ(B).
Basic definitions and results of cone theory may be needed within this paper. A set K is said to be a cone if αK⊂K for all α≥0. A cone is convex if for any two points in K it contains the line segment between them. A convex cone K is solid if the interior of K is nonempty. It is pointed if K∩(-K)={0}. A closed pointed solid convex cone is called a proper cone. A cone K is said to be polyhedral if it can be expressed as the set of nonnegative combinations of a finite set of generating vectors. We adopt the notation K=Cone(C) if K coincides with the set of nonnegative combinations of the n×k matrix C.
We discuss the phase-type distribution for a random variable X≥0 in the terms of a continuous-time Markov process. A continuous-time Markov process is defined on n+1 finite state space. The row vector β gives the initial probabilities with i state probability βi. A phase-type distribution is defined as the distribution of the time to absorption in a continuous-time Markov chain (CTMC) with one absorbing state [8]. If the n+1 state is an absorbing state and all other states are transient, we define a phase-type infinitesimal generator matrix of the Markov chain, denoted by (A-,b-,β-), such that (3)A-=A-A100,b-=en+1,β-=ββn+1,where 0 refers to the column vector, row vector, or matrix with all entries equal to zero in the case without ambiguity, ei is the corresponding row vector whose ith entry is one and the others are zero, and 1 is the column vector with all entries being one. We can see that βn+1=0 if β1=1, and β1+βn+1=1 otherwise. Phase-type distributions are commonly represented by a vector-matrix tuple (β,A) that describes the transient part of the CTMC. The vector-matrix tuple (β,A) is a phase-type (Markovian) representation of a phase-type distribution if and only if β1≤1 and β≥0, A is a Metzler matrix with A1≤0 and A1≠0, and A is nonsingular.
The probability density function (PDF), cumulative distribution function (CDF), and Laplace-Stieltjes transform (LST) of the PDF, respectively, are defined by (4)fx=βexpAx-A1,F(x)=1-βexpAx1,f∗s=EexpsX=βsI-A-1-A1,where E(·) is an expectation. The ME (matrix exponential) distribution is a generalization of the phase-type distribution. A distribution function is called an ME distribution if there exists the triple (A,u,β) such that F(x)=1-βexp(Ax)u, and there is no restriction on the elements of (A,u,β).
We note that these representations of phase-type distributions are equivalent to the state space realizations of linear systems in control. Let us consider single-input, single output linear time-invariant systems of the form (5)x˙(t)=Ax(t)+Bu(t),y(t)=Cx(t).The linear system in (5) is said to be a positive linear system provided that, for any nonnegative input and nonnegative initial state, the state trajectory and the output are always nonnegative. Let a transfer function H(s)=C(sI-A)-1B be defined by the Laplace transform of the impulse response function h(t)=Cexp(At)B for t≥0 and, otherwise, h(t)=0. A triple (A,B,C) is said to be a positive realization of H(s) in a continuous-time linear positive system if and only if A is a Metzler matrix, B≥0, and C≥0. A triple (A,B,C) is denoted by a minimal realization if (A,B,C) is jointly completely controllable and completely observable.
An integral function of h(t) is defined by hI(t)=∫0th(t)dt and hI(0)=0. We can see that an augmented realization (A-,B-,C-) of hI(t) is defined by (6)A-=AB00,B-=en+1,C-=C0.
The augmented realization (A-,B-,C-) presents a state space realization of the integral function, such as hI(t)=C-exp(A-t)B-. The augmented realization (A-,B-,C-) is closely related to representation (3). We note that the positive realization of the integration of a positive system is closely related to the representation of the phase-type distribution.
3. Phase-Type Representation and the Positive Realization3.1. Some Connections
Using a generalized version of the Perron-Frobenius theorem of nonnegative matrix theory, we derived a transformation from positive realization into phase-type realization under a constraint. The Perron-Frobenius theorem asserts that a real square matrix with positive entries has a unique largest real eigenvalue and that the corresponding eigenvector has strictly positive components [10, 19]. The Perron-Frobenius results of reducible matrices are characterized by being weaker than those of irreducible matrices. The general Perron-Frobenius theorem for reducible matrices is introduced from the result in Chapter 8 in [19] as follows.
Theorem 1 (Perron-Frobenius [19]).
Let A≥0 be an n×n matrix. Then, A has a nonnegative real eigenvalue equal to its spectral radius ρ(A) and there is an eigenvector x≥0 corresponding to ρ(A).
The solvability problem of the matrix equation (I-A)y=b with constraint y≥0 was originally solved by Carlson [11]. Several generalized versions have been discussed by researchers [12, 13, 20]. We note that the excitability is closely related to the existence of the strict positive solution of (I-A)y=b (i.e., y>0) as a more general case. For a discrete-time positive system with realization (A,b,c), (A,b) is defined to be excitable if there is an n such that ∑m=0nAmb>0, and (c,A) is defined to be transparent if there is an n such that ∑m=0ncAm>0 [2, 14]. We consider some properties of excitability.
Lemma 2.
Assume that A≥0 with ρ(A)<1. Let b≥0. Then the following statements are equivalent:
There is a y>0 such that (I-A)y=b and y=limn→∞∑m=0nAmb.
(A,b) is excitable.
Proof.
(1)⇒(2): we define yn=∑m=0nAmb. Assume that y is strictly positive and limn→∞yn=y. Therefore, there is an n such that yk-y∞<ϵ for all k≥n and ϵ=y∞/2 where y∞ is defined by ∞-norm of y. Because yn>0 for sufficiently large n, it follows that (A,b) is excitable.
(2)⇒(1): because A is stable (i.e., ρ(A)<1), the sum ∑k=1nAk converges uniformly to (I-A)-1 as n→∞. We have y=limn→∞yn. Because each entry of yn is monotonically increasing with respect to n and excitable, we can see that y is strictly positive.
Consider a continuous-time system with a positive realization (A,B,C). Since A is a Metzler matrix, we can choose an α>0 such that (I+αA)≥0 and ρ(I+αA)<1. We can define the excitability and the transparency of continuous-time positive linear systems in a similar form to discrete-time ones. The pair (A,B) is excitable if there is an integer n>0 such that (7)zn=∑k=0nI+αAkB>0for some α>0. The pair (C,A) is transparent if there exists an n such that (8)wn=∑k=0nCI+αAk>0.for some α>0.
Lemma 3.
For the continuous time with a positive realization (A,B,C), assume that (A,B) is excitable, A is an asymptotically stable Metzler matrix (i.e., λmax(A)<0), and an augmented realization (A-,B-,C-) is given as (6). Then there is a positive eigenvector v of the spectral radius σ(A-)=0(i.e., A-v=0) such that v is strictly positive (i.e., v>0).
Proof.
We can choose a sufficiently large η>0 such that a positive matrix A~=A-+ηI satisfies ρ(A~)=η, and η>|x| for all x∈σ(A). By using the Perron-Frobenius Theorem 1 for the augmented nonnegative matrix A~ with the order n+1, there exists an eigenvector v≥0 corresponding to η such that A~v=ηv. If the last entry of v is zero, vn+1=0, then it induces that a vector v∗≜v1v2⋯vnT is an eigenvector of A+ηI corresponding to an eigenvalue η. It contradicts the fact that all the eigenvalues of A+ηI are less than η. Therefore, we have vn+1>0.
Set A+=I+A/η. We have A+v∗+vn+1B/η=v∗ since A~v=ηv. Thus, A+ is stable and positive. By using Lemma 2, the excitable condition (7) implies a strict positive solution v∗=vn+1(I-A+)-1B/η>0. Therefore, it is shown that v is strictly positive; that is, v>0.
We will show that a positive realization of continuous-time positive system can be transformed into a phase-type representation normalized by a positive number. Under the irreducible assumption, it was proven that the positive realization can be transformed into phase-type representation [9]. We modify the proof as a generalized version of the correspondence between positive realizations and phase-type representation. We can see that a phase-type representation is a special positive realization with excitable constraint.
Theorem 4.
Consider the continuous-time positive system with the positive realization (A,B,C) such that (A,B) is excitable, and A is an asymptotically stable and Metzler matrix. Then it is transformed into a phase-type infinitesimal generator matrix (A-+,B-+,C-+) such that (9)A-+=A+-A+100,B-+=en+1,C-+=C+0,where A+1=-B+ and C+1>0.
Proof.
First, let us define an augmented realization (A-,B-,C-) in the form of (6). Using Lemma 3, the Perron eigenvector v of A- corresponding to the Perron root ρ(A-)=0 is strictly positive. Set U-=diag(v1,…,vn,vn+1)/vn+1. Thus, we can derive that A-+=U--1A-U- with A-+1=0, B-+=U--1B-=en+1, and C-+=C-U-. Therefore, we can verify that A+1+B1=0 where A+=U-1AU, B+=U-1B, C+=CU, and U=diag(v1,…,vn)/vn+1.
An important consequence of the above theorem is that an excitable positive realization can be transformed into the form of phase-type representation. Therefore, it is remarked that he concept of positive realizations is a superset of phase-type representations.
3.2. Common Properties and Characteristics
We discuss the properties and characteristics, such as stability, irreducibility, excitability, and transparency, in positive systems and phase-type distributions. A positive system with a positive realization (A,B,C) is said to be irreducible if (A,B) is excitable and (C,A) is transparent in the terminology introduced in [2, 14]. We note that the properties and characteristics of excitability and transparency are closely related to those of the reachability and observability of positive linear systems [17, 18]. A phase-type representation (α,A) in whose graph all the state vertices are connected to the initial vertex and to the absorbing vertex is called irreducible [8]. We note that the irreducible representation is closely related to the irreducibility of the phase-type renewal process in the Markovian point process introduced in [7, 8]. Renewal processes provide simple models of point processes, which may describe an ordered set of points. We consider a renewal process with a phase-type distribution for the interrenewal intervals. In [8], the phase-type representation (α,A) for the distribution function F(x) is called irreducible if Q∗ is irreducible where an infinitesimal generator Q∗ is defined by (10)Q∗=A-A1α.We may associate a Markov process with a phase-type renewal process. A renewal function denoted by R(x) for a phase-type distribution F(x) is defined by the expected number of renewals in the interval [0,x]; that is, R(x)=E[N(x)], where N(x) denotes the number of renewals. A renewal density is defined by r(x)=dR(x)/dx. The Laplace transform of the renewal density r(x) is denoted by r∗(s), which is rewritten by (11)r∗(s)=f∗s1-f∗s=-1+1-f∗s-1.In view of control theory, the equation in (11) is equivalent to the positive feedback control. Because the state space realization of the inverse system (1-f∗(s))-1 is given by (1-f∗(s))-1=1+α(sI-Q∗)-1(-A1), we can see that its renewal density is given by r(x)=αexp(Q∗x)(-A1), which is equal to the results in [7]. For an irreducible representation (α,A), the vectors αexp(Tx) and exp(Ax)(-A1) are strictly positive [8]. We note that these results are related to the excitability and transparency of the positive linear system.
The irreducibility of a positive system can be defined in a similar manner. A positive system is irreducible if each state variable influences and is influenced by another variable [2]. It is defined by an irreducible realization for the positive system with a positive realization (A,B,C) if Q is irreducible, where Q is defined by (12)Q=A+BC.For the open-loop system (5), the associated closed loop system (positive feedback system) is given by (13)x˙t=Axt+Byt=Axt+BCxt=Qxt,where the linear state-feedback law is u(t)=Cx(t) and C is the constant feedback gain row vector. The closed loop system (13) is positive if and only if Q is a Metzler matrix. The stabilization problem of positive systems has recently been discussed in [3, 15, 16]. It is known an unstable open positive system (5) cannot be stabilized by linear state-feedback if the restriction on nonnegative control in the closed loop is imposed [15, 16].
The properties and characteristics of excitability and transparency are closely related to the reachability and observability of positive linear systems [5, 17, 18]. A reachable set R is the set of all points which the states approach from the origin by nonnegative inputs within finite time. It was shown that R=cl[cone{x∣x=exp(At)B,t≥0}], where cl(S) is a closure set of S [17]. An observable set S is the set of initial states in which the output is nonnegative for all t≥0. The observable set S can be defined by S={x∣Cexp(At)x≥0,∀t≥0}.
Theorem 5 (see [17]).
Let a transfer function H(s) be a strictly proper rational function with degree n, whose realization is given by (A,B,C). Then, H(s) has a positive realization (A∗,B+,C+) with a Metzler matrix A∗=A+-ηI if and only if there exists a generator matrix P and λ≥0 such that a polyhedral cone P=cone(P) satisfies
(A+λI)P⊂P;
R⊂P⊂S,
where P∈Rn×m, n≤m, R is a reachable set, and S is an observable set.
A positive system (A,B,C) is said to be reducible otherwise. When a representation is not irreducible, it can be simplified by discarding some states. Our next question is how to discard some unnecessary states when a representation is not irreducible. We discuss an order reduction algorithm of the asymptotical stable and unexcitable positive realization. Let a set 〈n〉={1,…,n}. For a column vector b with n entries, we define a support set(14)supp(b)=i∈n∣bi≠0.
Theorem 6.
Assume that a positive realization (A∗,b+,c+) is unexcitable and A∗ is asymptotically stable. Then there is a permutation matrix M=M1M2 such that (15)A∗M=M1M2A11A120A22,b+=M1M2e10,c+M1M2=f1f2,where A11 is an n1×n1 matrix. Furthermore, we have a reduced positive realization (A11,e1,f1) such that (16)A∗M1=M1A11,b+=M1e1,c+M1=f1.
Proof.
Because A∗ is a Metzler matrix and asymptotical stable, there is an η such that A+=A∗+ηI is a positive matrix. Define zm=∑k=0m(I+αA)kB with α=1/η. There is an integer N>0 such that supp(zm)=supp(zm+1) for all m≥N. A support set for zN is defined by Z-=supp(zN). Let n1 be the element number of Z-. We can find a permutation matrix M such that supp(MTzN)=〈n1〉. Set A1=MTAM=A11A12A21A22 where the size of A11 is n1. If A21 is a nonzero matrix, then we have supp(A1MTzN)⊈〈n1〉, but this contradicts the definition of zm. We can see that there is a permutation matrix M=M1M2 such that A∗M=M1M2A11A120A22 and supp(e1)⊂〈n1〉. We have R⊂cone(M1). Because A∗ is a Metzler matrix and c+>0, we can see cone(M1)⊂S by the definition of S. By using Theorem 5, we can derive (16).
We discussed the method to remove unnecessary states in the unexcitable case. When the transposed realization is given by (AT,CT,BT), the concept of the transparency can be interpreted by that of the excitability. A removing method of the unnecessary state in the nontransparent case is similar to that in the unexcitable case. We illustrate the previous theorems by means of an example.
Example 7.
A positive state-space realization (A,B,C) is given by (17)A=-21110-30011-10010-4,B=1010C=2413.We can see that Q=A+BC is reducible. Set α=1/5. We have z5=4.1139205.489920T. The support set is Z-={1,3}. We obtain (18)M=1000001001000001.By using Theorem 6 and removing the unexcitable part, we obtain an excitable positive realization (A11,e1,f1) such that A11=-211-1, e1=11, and f1=21T. By using Theorem 4, the phase-type representation (A,β) can be obtained as follows: (19)A=-2.0001.5000.666-1.000,β=43.
4. Discrete Phase-Type Distributions and Discrete-Time Positive Systems
A discrete phase-type (DPT) distribution is the distribution of the time until one absorbing state in a discrete-state discrete-time Markov chain (DTMC) with n transients states and one absorbing state [8, 21]. DPT distributions have received little attention in applied stochastic modeling. The main research activity has addressed continuous phase-type distributions. Let τ be the time till absorption into state n+1 in the DTMC. We say that τ is a random variable of order n and representation (A,b,β) [8]. The DPT representation has the following properties. A positive matrix A=[aij] is an (n×n) matrix grouping the transition probabilities among the transient states. A column vector b=[bi] is a positive n-dimensional column vector grouping the probabilities from any state to the absorbing state. Thus, we have ∑j=1naij=1-bi. It means that (I-A)1=b. An n vector β is defined by the initial probability vector. The augmented matrix tuple (A^,b^,β^) is denoted by the augmented phase-type representation where the one-step transition probability matrix A^ of the corresponding DTMC can be partitioned by (20)A^=Ab01,the augmented initial vector is defined by β^=β0, and b^=en+1. Its probability generating function is defined by P(z)=β(z-1I-A)-1b [8].
We can also discuss the realization between the DPH distributions and discrete-time positive systems in a similar manner. The discrete-time linear system is represented by (21)xt+1=Axt+But,yt=Cxt,where A∈R+n×n, B∈R+n, and CT∈R+n. The matrix tuple (A,B,C) denotes the positive realization (21). In the next theorem, we show that the positive realization can be transformed into a DPH representation multiplied by a positive scalar (i.e., β is not necessarily a probability vector).
Theorem 8.
Assume that a realization (A,B,C) is denoted by a positive realization satisfying (21) and (A,B) is excitable (essential reachable) and stable. Then there is a nonsingular matrix M such that the realization (A~,b~,β~), which is defined by β~=CM, A~=M-1AM, and b~=M-1B, has the properties of the DPH representation such as b~=(I-A~)1 and β≥0.
Proof.
Because A is positive, stable, and excitable, the absolute values of its eigenvalues are less than 1 and we can use Lemma 2. Thus, we obtain the fact that the entries of x=(I-A)-1B are positive. A similarity transform matrix M is defined by a diagonal matrix, M=diag(x). Compute a new realization (A~,b~,β~). We obtain the facts that (I-A~)1=b~≥0, β~≥0, and A≥0. Therefore, we can verify that the new realization satisfies the discrete phase-type distribution properties.
We can easily deploy the properties and characteristics, such as irreducibility, excitability, transparency, and order reduction, in the discrete domain in a similar manner as in the continuous case. We omit the detailed exploration for the discrete case in this paper.
5. Conclusions
We considered the relation between the positive realization and the phase-type representation in continuous time and discrete time, respectively. Using the Perron-Frobenius theorem, it was shown that a phase-type representation is a special case with excitable constraint of the positive realization. We discussed their common properties and characteristics, such as irreducibility, excitability, transparency, stabilization, and order reduction. The connection between the phase-type renewal process and the feedback control of positive system was discussed. A lot of open problems related to positive system still remain and should be addressed in future research. The communities of control and probability theory can work together on solving the remaining same open problems.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by research fund of Chungnam National University.
BenvenutiL.FarinaL.AndersonB. D. O.Filtering through combination of positive filters199946121431144010.1109/81.8095452-s2.0-0033338653FarinaL.RinaldiS.2000Wiley10.1002/9781118033029MR1784150KaczorekT.SajewskiL.20141Springer10.1007/978-3-319-04834-5MR3185228KimK.A construction method for positive realizations with an order bound201261775976510.1016/j.sysconle.2012.04.0092-s2.0-84861220467KimK.A constructive positive realization with sparse matrices for a continuous-time positive linear system20132013987814610.1155/2013/878146CumaniA.On the canonical representation of homogeneous markov processes modelling failure-time distributions198222358360210.1016/0026-2714(82)90033-62-s2.0-0020419336LatoucheG.RamaswamiV.199943SIAMNeutsM. F.1994New York, NY, USADover PublicationsMR1313503CommaultC.MocanuS.Phase-type distributions and representations: some results and open problems for system theory20037665665802-s2.0-014219933910.1080/0020717031000114986MR1979678BermanA.PlemmonsR. J.1994Society for Industrial and Applied Mathematics10.1137/1.9781611971262MR1298430CarlsonD.A note on M-matrix equations196311410271033HershkowitzD.SchneiderH.Solutions of Z-matrix equations1988106253810.1016/0024-3795(88)90020-1MR951825SchneiderH.The influence of the marked reduced graph of a nonnegative matrix on the Jordan form and on related properties: a survey19868416118910.1016/0024-3795(86)90313-7MuratoriS.RinaldiS.Excitability, stability, and sign of equilibria in positive linear systems1991161596310.1016/0167-6911(91)90030-i2-s2.0-0025995207De LeenheerP.AeyelsD.Stabilization of positive linear systems200144425927110.1016/s0167-6911(01)00146-32-s2.0-0035889361JamesG.RumchevV.Stability of positive linear discrete-time systems200553118OhtaY.MaedaH.KodamaS.Reachability, observability, and realizability of continuous-time positive systems198422217118010.1137/0322013MR732422ZBL0539.930052-s2.0-0021394486ValcherM. E.Reachability properties of continuous-time positive systems20095471586159010.1109/tac.2009.20155562-s2.0-67949083442MR2535754MeyerC.2000Society for Industrial and Applied Mathematics10.1137/1.9780898719512MR1777382VictoryJ.On nonnegative solutions of matrix equations19856340641210.1137/0606042MR791170BobbioA.HorvathA.ScarpaM.TelekM.Acyclic discrete phase type distributions: properties and a parameter estimation algorithm200354113210.1016/s0166-5316(03)00044-02-s2.0-0043123281