Skip to main content

On a singular system of fractional nabla difference equations with boundary conditions

Abstract

In this article, we study a boundary value problem of a class of linear singular systems of fractional nabla difference equations whose coefficients are constant matrices. By taking into consideration the cases that the matrices are square with the leading coefficient matrix singular, square with an identically zero matrix pencil and non-square, we provide necessary and sufficient conditions for the existence and uniqueness of solutions. More analytically, we study the conditions under which the boundary value problem has a unique solution, infinite solutions and no solutions. Furthermore, we provide a formula for the case of the unique solution. Finally, numerical examples are given to justify our theory.

1 Introduction

Difference equations of fractional order have recently proven to be valuable tools in the modeling of many phenomena in various fields of science and engineering. Indeed, we can find numerous applications in viscoelasticity, electrochemistry, control, porous media, electromagnetism, and so forth [1–7]. There has been a significant development in the study of fractional differential/difference equations and inclusions in recent years; see the monographs of Baleanu et al. [1], Kaczorek [4], Klamka et al. [8], Malinowska et al. [5], Podlubny [7], and the survey by Agarwal et al. [9]. For some recent contributions on fractional differential/difference equations, see [1, 4, 5, 8–27] and the references therein. In this article we provide an introductory study for a boundary value problem of a class of singular fractional nabla discrete time systems. If we define N α by N α ={α,α+1,α+2,…}, α integer, and n such that 0<n<1 or 1<n<2, then the nabla fractional operator in the case of Riemann-Liouville fractional difference of n th order for any Y k : N a → R m × 1 is defined by, see [5, 10–12, 23],

∇ α n Y k = 1 Γ ( − n ) ∑ j = α k ( k − j + 1 ) − n − 1 ¯ Y j ,

where the raising power function is defined by

k α ¯ = Γ ( k + α ) Γ ( k ) .

The following problem will then be considered. The singular fractional discrete time systems of the form

F ∇ 0 n Y k =G Y k ,k=1,2,…,N
(1)

with known boundary conditions

A 1 Y 0 = B 1 , A 2 Y N = B 2 ,
(2)

where F,G∈M(r×m;F) (i.e., the algebra of matrices with elements in the field ℱ) with Y k , V k ∈M(m×1;F), A 1 ∈M( r 1 ×m;F), A 2 ∈M( r 2 ×m;F), B 1 ∈M( r 1 ×1;F) and B 2 ∈M( r 2 ×1;F). For the sake of simplicity, we set M m =M(m×m;F) and M r m =M(r×m;F). The matrices F and G can be non-square (when r≠m) or square (r=m) and F singular (detF=0). The main purpose will be to provide necessary and sufficient conditions for the existence and uniqueness of solutions of the above boundary value problem, i.e., to study the conditions under which the system has unique, infinite and no solutions and to provide a formula for the case of the unique solution (if it exists). Many authors use matrix pencil theory to study linear discrete time systems with constant matrices; see, for instance, [28–43]. A matrix pencil is a family of matrices sF−G, parametrized by a complex number s, see [39, 41, 44, 45]. When G is square and F= I m , where I m is the identity matrix, the zeros of the function det(sF−G) are the eigenvalues of G. Consequently, the problem of finding the nontrivial solutions of the equation

sFX=GX

is called the generalized eigenvalue problem. Although the generalized eigenvalue problem looks like a simple generalization of the usual eigenvalue problem, it exhibits some important differences. In the first place, it is possible for F, G to be non-square matrices. Moreover, even with F, G square it is possible (in the case F is singular) for det(sF−G) to be identically zero, independent of s. Finally, even if we assume F, G square matrices with a non-zero pencil, it is possible (when F is singular) for the problem to have infinite eigenvalues. To see this, write the generalized eigenvalue problem in the reciprocal form

FX= s − 1 GX.

If F is singular with a null vector X, then FX= 0 m , 1 , so that X is an eigenvector of the reciprocal problem corresponding to eigenvalue s − 1 =0; i.e., s=∞.

Definition 1.1 Given F,G∈ M r m and an arbitrary s∈F, the matrix pencil sF−G is called:

  1. 1.

    Regular when r=m and det(sF−G)≠0.

  2. 2.

    Singular when r≠m or r=m and det(sF−G)≡0.

The paper is organized as follows. In Section 2, we study the existence of solutions of the system (1) when its pencil is regular. In Section 3 we study the case of the system (1) with a singular pencil, and Section 3 contains numerical examples.

2 Regular case

In this section, we consider the case of the system (1) with a regular pencil. The class of sF−G is characterized by a uniquely defined element, known as complex Weierstrass canonical form, s F w − Q w , see [39, 41, 44, 45], specified by the complete set of invariants of sF−G. This is the set of elementary divisors (e.d.) obtained by factorizing the invariant polynomials into powers of homogeneous polynomials irreducible over the field ℱ. In the case where sF−G is regular, we have e.d. of the following type:

  • e.d. of the type ( s − a j ) p j are called finite elementary divisors (f.e.d.), where a j is a finite eigenvalue of algebraic multiplicity p j ;

  • e.d. of the type s ˆ q = 1 s q are called infinite elementary divisors (i.e.d.), where q is the algebraic multiplicity of the infinite eigenvalues.

We assume that ∑ i = 1 ν p j =p and p+q=m.

Definition 2.1 Let B 1 , B 2 ,…, B l be elements of M m . The direct sum of them denoted by B 1 ⊕ B 2 ⊕…⊕ B l is the blockdiag[ B 1 B 2 ⋯ B l ].

From the regularity of sF−G, there exist nonsingular matrices P,Q∈ M m such that

PFQ= F w = I p ⊕ H q
(3)

and

PGQ= G w = J p ⊕ I q .
(4)

The complex Weierstrass form s F w − Q w of the regular pencil sF−G is defined by

s F w − Q w :=s I p − J p ⊕s H q − I q ,

where the first normal Jordan-type element is uniquely defined by the set of the finite eigenvalues of sF−G and has the form

s I p − J p :=s I p 1 − J p 1 ( a 1 )⊕…⊕s I p ν − J p ν ( a ν ).

The second uniquely defined block s H q − I q corresponds to the infinite eigenvalues of sF−G and has the form

s H q − I q :=s H q 1 − I q 1 ⊕…⊕s H q σ − I q σ .

The matrix H q is a nilpotent element of M q with index q ∗ =max{ q j :j=1,2,…,σ}, where

H q q ∗ = 0 q , q

and I p j , J p j ( a j ), H q j are defined as

I p j = [ 1 0 ⋯ 0 0 0 1 ⋯ 0 0 ⋮ ⋮ ⋱ ⋮ ⋮ 0 0 ⋯ 0 1 ] ∈ M p j , J p j ( a j ) = [ a j 1 ⋯ 0 0 0 a j ⋯ 0 0 ⋮ ⋮ ⋱ ⋮ ⋮ 0 0 ⋯ a j 1 0 0 ⋯ 0 a j ] ∈ M p j , H q j = [ 0 1 ⋯ 0 0 0 0 ⋯ 0 0 ⋮ ⋮ ⋱ ⋮ ⋮ 0 0 ⋯ 0 1 0 0 ⋯ 0 0 ] ∈ M q j .

For algorithms about the computations of the Jordan matrices, see [39, 41, 44, 45].

Definition 2.2 If for the system (1) with boundary conditions (2) there exists at least one solution, the boundary value problem (1)-(2) is said to be consistent.

For the regular matrix pencil of the system (1), there exist nonsingular matrices P,Q∈ M m as applied in (3), (4). Let

Q=[ Q p Q q ],
(5)

where Q p ∈ M m p is a matrix with columns p linear independent (generalized) eigenvectors of the p finite eigenvalues of sF−G, and Q q ∈ M m q is a matrix with columns q linear independent (generalized) eigenvectors of the q infinite eigenvalues of sF−G.

Lemma 2.1 Consider the system (1) with a regular pencil. Then the system (1) is divided into two subsystems:

∇ 0 n Z k p = J p Z k p

and

H q ∇ 0 n Z k q = Z k q .

Proof Consider the transformation

Y k =Q Z k
(6)

and by substituting (6) into (1), we obtain

F ∇ 0 n Q Z k =GQ Z k

or, equivalently,

FQ ∇ 0 n Z k =GQ Z k .

Whereby multiplying by P, we arrive at

F w ∇ 0 n Z k = G w Z k .

Moreover, let

Z k = [ Z k p Z k q ] ,

where Z k p ∈ M p 1 , Z k q ∈ M q 1 , and by using (3) and (4), we obtain

[ I p 0 p , q 0 q , p H q ] ∇ α n [ Z k p Z k q ] = [ J p 0 p , q 0 q , p I q ] [ Z k p Z k q ] .

From the above expressions, we arrive easily at the subsystems

∇ 0 n Z k p = J p Z k p
(7)

and

H q ∇ 0 n Z k q = Z k q .
(8)

The proof is completed. □

Definition 2.3 With F n , n ( J p ( k + n ) n ¯ ) we denote the discrete Mittag-Leffler function with two parameters defined by

F n , n ( J p ( k + n ) n ¯ ) = ∑ i = 0 ∞ J p i ( k + n ) i n ¯ Γ ( ( i + 1 ) n ) .
(9)

See [10–12, 23, 46].

Proposition 2.1 The subsystem (7) has the solution

Z k p = ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p ) Z 0 p
(10)

if and only if

∥ J p ∥<1,
(11)

where ∥⋅∥ is an induced matrix norm and F n , n ( J p ( k + n ) n ¯ ) is the discrete Mittag-Leffler function with two parameters as defined by Definition  2.3.

Proof From [10–12, 23, 46] the solution of (7) can be calculated and given by the formula

Z k p = ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p ) Z 0 p

or, equivalently, by

Z k p = ( k + 1 ) n − 1 ¯ ( ∑ i = 0 ∞ J p i ( k + n ) i n ¯ Γ ( ( i + 1 ) n ) ) ( I p − J p ) Z 0 p .

The existence and uniqueness of the above solution depends on the convergence of the matrix power series

∑ i = 0 ∞ J p i ( k + n ) i n ¯ Γ ( ( i + 1 ) n )

or, equivalently, if and only if

lim i → ∞ ∥ J p ( i + 1 ) ( k + n ) ( i + 1 ) n ¯ Γ ( ( i + 2 ) n ) ∥ ∥ J p i ( k + n ) i n ¯ Γ ( ( i + 1 ) n ) ∥ <1

or, equivalently,

∥ J p ∥< lim i → ∞ ( k + n ) i n ¯ Γ ( ( i + 1 ) n ) ( k + n ) ( i + 1 ) n ¯ Γ ( ( i + 2 ) n ) .

By using the property

Γ(z+1)=zΓ(z),

we get

∥ J p ∥< lim i → ∞ ( ( k − 1 ) + ( i + 1 ) n ) ⋯ ( 1 + ( i + 1 ) n ) ( ( i + 1 ) n ) ( ( k − 1 ) + ( i + 2 ) n ) ⋯ ( 1 + ( i + 2 ) n ) ( ( i + 2 ) n )

or, equivalently,

∥ J p ∥<1.

The proof is completed. □

Proposition 2.2 The subsystem (8) has the unique solution

Z k q = 0 q , 1 .
(12)

Proof Let q ∗ be the index of the nilpotent matrix H q , i.e., H q q ∗ = 0 q , q . Then if we obtain the following equations:

H q ∇ 0 n Z k q = Z k q , H q 2 ∇ 0 2 n Z k q = H q ∇ 0 n Z k q , H q 3 ∇ 0 3 n Z k q = H q 2 ∇ 0 2 n Z k q , H q 4 ∇ 0 4 n Z k q = H q 3 ∇ 0 3 n Z k q , ⋮ H q q ∗ − 1 ∇ 0 ( q ∗ − 1 ) n Z k q = H q q ∗ − 2 ∇ 0 ( q ∗ − 2 ) n Z k q , H q q ∗ ∇ 0 q ∗ n Z k q = H q q ∗ − 1 ∇ 0 ( q ∗ − 1 ) n Z k q

by taking the sum of the above equations and using the fact that H q q ∗ = 0 q , q , we arrive easily at the solution (12). The proof is completed. □

Theorem 2.1 Consider the system (1) with a regular pencil and boundary conditions of type (2). Then the boundary value problem (1)-(2) is consistent if and only if:

  1. 1.

    The pencil sF−G has p distinct eigenvalues and all lie within the open disk

    |s|<1;
  2. 2.
    [ B 1 B 2 ] ∈colspan [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] .
    (13)

Furthermore, when the boundary value problem (1)-(2) is consistent, it has a unique solution if and only if:

  1. 1.
    p≤ r 1 + r 2 ;
    (14)
  2. 2.
    rank [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] =p.
    (15)

In this case the unique solution is then given by

Y k = Q p ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p )C,
(16)

where C is the unique solution of the algebraic system

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] C= [ B 1 B 2 ] .
(17)

Proof By applying the transformation (6) into the system (1), we get the systems (7), (8) with solutions (10), (12) respectively. Note that from Proposition 2.1 the solution (10) exists if and only if

∥ J p ∥<1,

where J p is the Jordan matrix related to the p finite eigenvalues of the pencil sF−G, which is equivalent to the fact that the finite eigenvalues of the pencil must be distinct and all lie within the unit disk |s|<1. Based on these results, the solution of (1) can be written as

Y k =Q Z k =[ Q p Q q ] [ Z k p Z k q ]

or, equivalently,

Y k = Q p Z k p + Q q Z k q ,

or, equivalently, by using (10), (12)

Y k = Q p ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p ) Z 0 p .

The initial value Z 0 p of the subsystem (7) is not known and can be replaced by a constant vector C∈ M p 1

Y k = Q p ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p )C.

The above solution exists if and only if

A 1 Y 0 = B 1 , A 2 Y N = B 2

or, equivalently,

A 1 Q p C= B 1 , A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p )C= B 2 ,

or, equivalently,

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] C= [ B 1 B 2 ] .

For the above algebraic system, there exists at least one solution if and only if

[ B 1 B 2 ] ∈colspan [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] .

The algebraic system (17) contains r 1 + r 2 equations and p unknowns. Hence the solution is unique if and only if

p≤ r 1 + r 2

and

rank [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] =p,

where C is then the unique solution of (17). This can be proved as follows. If we assume that the algebraic system has two solutions C 1 and C 2 , then

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] C 1 = [ B 1 B 2 ]

and

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] C 2 = [ B 1 B 2 ]

or, equivalently,

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] ( C 1 − C 2 )= 0 p , 1 .

But the matrix [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] is left invertible since it is assumed to have p linear independent columns and r 1 + r 2 ≥p and hence

C 1 = C 2 .

The unique solution is then given from (16). The proof is completed. □

3 Singular case

In this section, we consider the case of the system (1) with a singular pencil. The class of sF−G in this case is characterized by a uniquely defined element, s F K − Q K , known as the complex Kronecker canonical form, see [39, 41, 44, 45], specified by the complete set of invariants of the singular pencil sF−G. This is the set of the elementary divisors (e.d.) and the minimal indices (m.i.). Unlike the case of the regular pencils, where the pencil is characterized only from the e.d., the characterization of a singular matrix pencil apart from the set of the determinantal divisors requires the definition of additional sets of invariants, the minimal indices. The distinguishing feature of a singular pencil sF−G is that either r≠m or r=m and det(sF−G)≡0. Let N r , N l be the right and the left null space of a matrix respectively. Then the equations

(sF−G)U(s)= 0 r , 1 ,
(18)
V T (s)(sF−G)= 0 1 , m ,
(19)

where ( ) T is the transpose tensor, have solutions in U(s), V(s), which are vectors in the rational vector spaces N r (sF−G) and N l (sF−G) respectively. The binary vectors U(s) and V T (s) express dependence relationships among the columns or rows of sF−G respectively. Note that U(s)∈ M m 1 and V(s)∈ M r 1 are polynomial vectors. Let d=dim N r (sF−G) and t=dim N l (sF−G). It is known, see [39, 41, 44, 45], that N r (sF−G) and N l (sF−G) as rational vector spaces are spanned by minimal polynomial bases of minimal degrees

ϵ 1 = ϵ 2 =⋯= ϵ g =0< ϵ g + 1 ≤⋯≤ ϵ d
(20)

and

ζ 1 = ζ 2 =⋯= ζ h =0< ζ h + 1 ≤⋯≤ ζ t
(21)

respectively. The set of minimal indices ϵ 1 , ϵ 2 ,…, ϵ d and ζ 1 , ζ 2 ,…, ζ t are known as column minimal indices (c.m.i.) and row minimal indices (r.m.i.) of sF−G respectively. To sum up, in the case of a singular pencil, we have invariants of the following type:

  • finite elementary divisors of the type ( s − a j ) p j ;

  • infinite elementary divisors of the type s ˆ q = 1 s q ;

  • column minimal indices of the type ϵ 1 = ϵ 2 =⋯= ϵ g =0< ϵ g + 1 ≤⋯≤ ϵ d ;

  • row minimal indices of the type ζ 1 = ζ 2 =⋯= ζ h =0< ζ h + 1 ≤⋯≤ ζ t .

The Kronecker canonical form, see [39, 41, 44, 45], is defined by

s F K − G K :=s I p − J p ⊕s H q − I q ⊕s F ϵ − G ϵ ⊕s F ζ − G ζ ⊕ 0 h , g ,
(22)

where s I p − J p , s H q − I q are defined as in Section 2. The matrices F ϵ , G ϵ , F ζ and G ζ are defined by

F ϵ =blockdiag{ L ϵ g + 1 , L ϵ g + 2 ,…, L ϵ d },
(23)

where L ϵ =[ I ϵ ⋮ 0 ϵ , 1 ] for ϵ= ϵ g + 1 ,…, ϵ d

G ϵ =blockdiag{ L ¯ ϵ g + 1 , L ¯ ϵ g + 2 ,…, L ¯ ϵ d },
(24)

where L ¯ ϵ =[ 0 ϵ , 1 ⋮ I ϵ ] for ϵ= ϵ g + 1 ,…, ϵ d . The matrices F ζ , G ζ are defined as

F ζ =blockdiag{ L ζ h + 1 , L ζ h + 2 ,…, L ζ t },
(25)

where L ζ = [ I ζ 0 1 , ζ ] for ζ= ζ h + 1 ,…, ζ t

G ζ =blockdiag{ L ¯ ζ h + 1 , L ¯ ζ h + 2 ,…, L ¯ ζ t },
(26)

where L ¯ ζ = [ 0 1 , ζ I ζ ] for ζ= ζ h + 1 ,…, ζ t .

For algorithms about the computations of these matrices, see [39, 41, 44, 45].

Following the above given analysis, there exist nonsingular matrices P, Q with P∈ M r , Q∈ M m such that

P F Q = F K , P G Q = G K .
(27)

Let

Q=[ Q p Q q Q ϵ Q ζ Q g ],
(28)

where Q p ∈ M m p , Q q ∈ M m q , Q ϵ ∈ M m ( d − g ) , Q ζ ∈ M m ( t − h ) and Q g ∈ M m g .

Lemma 3.1 The system (1) is divided into five subsystems:

∇ 0 n Z k p = J p Z k p ,
(29)

the subsystem

H q ∇ 0 n Z k q = Z k q ,
(30)

the subsystem

F ϵ ∇ 0 n Z k ϵ = G ϵ Z k ϵ ,
(31)

the subsystem

F ζ ∇ 0 n Z k ζ = G ζ Z k ζ ,
(32)

and the subsystem

0 h , g ⋅ ∇ 0 n Z k + 1 g = 0 h , g ⋅ Z k g .
(33)

Proof Consider the transformation

Y k =Q Z k .
(34)

Substituting the previous expression into (1), we obtain

FQ ∇ 0 n Z k =GQ Z k .

Whereby multiplying by P and using (27), we arrive at

F K ∇ 0 n Z k = G K Z k .
(35)

Moreover, let

Z k = [ Z k p Z k q Z k ϵ Z k ζ Z k g ] ,

where Z p k ∈ M p 1 , Z q k ∈ M q 1 , Z ϵ k ∈ M ( d − g ) 1 , Z ζ k ∈ M ( t − h ) 1 and Z g k ∈ M g 1 . Taking into account the above expressions, we arrive easily at the subsystems (29), (30), (31), (32), and (33). The proof is completed. □

Solving the system (1) is equivalent to solving subsystems (29), (30), (31), (32) and (33). The solutions of the systems (29), (30) are given by (10) and (12) respectively; see Propositions 2.1 and 2.2.

Proposition 3.1 The subsystem (31) has infinite solutions and can be taken arbitrarily

Z k ϵ = C k , 1 .
(36)

Proof If we set

Z k ϵ = [ Z k ϵ g + 1 Z k ϵ g + 2 ⋮ Z k ϵ d ]

by using (23), (24), the system (31) can be written as

blockdiag{ L ϵ g + 1 ,…, L ϵ d } [ ∇ 0 n Z k ϵ g + 1 ∇ 0 n Z k ϵ g + 2 ⋮ ∇ 0 n Z k ϵ d ] =blockdiag{ L ¯ ϵ g + 1 ,…, L ¯ ϵ d } [ Z k ϵ g + 1 Z k ϵ g + 2 ⋮ Z k ϵ d ] .
(37)

Then, for the non-zero blocks, a typical equation from (37) can be written as

L ϵ i ∇ 0 n Z k ϵ i = L ¯ ϵ i Z k ϵ i ,i=g+1,g+2,…,d
(38)

or, equivalently,

[ I ϵ i ⋮ 0 ϵ i , 1 ] ∇ 0 n Z k ϵ i =[ 0 ϵ i , 1 ⋮ I ϵ i ] Z k ϵ i ,

or, equivalently,

[ 1 0 ⋯ 0 0 0 1 ⋯ 0 0 ⋮ ⋮ ⋯ ⋮ ⋮ 0 0 ⋯ 1 0 ] [ ∇ 0 n z k ϵ i , 1 ∇ 0 n z k ϵ i , 2 ⋮ ∇ 0 n z k ϵ i , ϵ i ∇ 0 n z k ϵ i , ϵ i + 1 ] = [ 0 1 ⋯ 0 0 0 0 ⋯ 0 0 ⋮ ⋮ ⋯ ⋮ ⋮ 0 0 ⋯ 0 1 ] [ z k ϵ i , 1 z k ϵ i , 2 ⋮ z k ϵ i , ϵ i z k ϵ i , ϵ i + 1 ] ,

or, equivalently,

∇ 0 n z k ϵ i , 1 = z k ϵ i , 2 + v k ϵ i , 1 , ∇ 0 n z k ϵ i , 2 = z k ϵ i , 3 + v k ϵ i , 2 , ⋮ ∇ 0 n z k ϵ i , ϵ i = ∇ 0 n z k ϵ i , ϵ i + 1 .
(39)

The system (39) is a regular-type system of difference equations with ϵ i equations and ϵ i +1 unknowns. It is clear from the above analysis that in every one of the d−g subsystems one of the coordinates of the solution has to be arbitrary by assigned total. The solution of the system can be assigned arbitrarily

Z k ϵ = C k , 1 .

The proof is completed. □

Proposition 3.2 The solution of the system (32) is unique and is the zero solution

Z k ζ = 0 t − h , 1 .
(40)

Proof From (25), (26) the subsystem (32) can be written as

blockdiag{ L ζ h + 1 ,…, L ζ t } [ ∇ 0 n Z k ζ h + 1 ∇ 0 n Z k ζ h + 2 ⋮ ∇ 0 n Z k ζ t ] =blockdiag{ L ¯ ζ h + 1 ,…, L ¯ ζ t } [ Z k ζ h + 1 Z k ζ h + 2 ⋮ Z k ζ t ] .
(41)

Then for the non-zero blocks, a typical equation from (41) can be written as

L ζ j ∇ 0 n Z k ζ j = L ¯ ζ j Z k ζ j ,j=h+1,h+2,…,t
(42)

or, equivalently,

[ I ζ j ⋯ 0 1 , ζ j ] ∇ 0 n Z k ζ j = [ 0 1 , ζ j ⋯ I ζ j ] Z k ζ j ,

or, equivalently,

[ 1 0 ⋯ 0 0 1 ⋯ 0 ⋮ ⋮ ⋯ ⋮ 0 0 ⋯ 1 0 0 ⋯ 0 ] [ ∇ 0 n z k ζ j , 1 ∇ 0 n z k ζ j , 2 ⋮ ∇ 0 n z k ζ j , ζ j ] = [ 0 0 ⋯ 0 1 0 ⋯ 0 ⋮ ⋮ ⋯ ⋮ 0 0 ⋯ 0 0 0 ⋯ 1 ] [ z k ζ j , 1 z k ζ j , 2 ⋮ z k ζ j , ζ j ] ,

or, equivalently,

∇ 0 n z k ζ j , 1 = 0 , ∇ 0 n z k ζ j , 2 = z k ζ j , 1 , ⋮ ∇ 0 n z k ζ j , ζ j − 1 = z k ζ j , ζ j − 2 , ∇ 0 n z k ζ j , ζ j = z k ζ j , ζ j − 1 , 0 = z k ζ j , ζ j .
(43)

We have a system of ζ j +1 difference equations and ζ j unknowns. Starting from the last equation, we get the solutions

z k ζ j , ζ j = 0 , z k ζ j , ζ j − 1 = 0 , z k ζ j , ζ j − 2 = 0 , ⋮ z k ζ j , 1 = 0 ,

which means that the solution of the system (32) is unique and is the zero solution. The proof is completed. □

Proposition 3.3 The subsystem (33) has an infinite number of solutions that can be taken arbitrarily

Z k g = C k , 2 .
(44)

Proof It is easy to observe that the subsystem

0 h , g ⋅ ∇ 0 n Z k + 1 g = 0 h , g ⋅ Z k g

does not provide any non-zero equations. Hence all its solutions can be taken arbitrarily. The proof is completed. □

We can now state the following theorem.

Theorem 3.1 Consider the system (1) with a singular pencil and known boundary conditions of type (2). Then the boundary value problem (1)-(2) is consistent if and only if:

  1. 1.
    ∥ J p ∥<1;
  2. 2.

    the column minimal indices are zero, i.e.,

    dim N r (sF−G)=0;
    (45)
  3. 3.
    [ B 1 B 2 ] ∈colspan [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] .
    (46)

Furthermore, when the boundary value problem (1)-(2) is consistent, it has a unique solution if and only if

  1. 1.
    p≤ r 1 + r 2 ;
    (47)
  2. 2.
    rank [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] =p.
    (48)

In this case the unique solution is given by the formula

Y k = Q p ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p )C,
(49)

where C is the unique solution of the algebraic system

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] C= [ B 1 B 2 ] .
(50)

In any other case the system has infinite solutions.

Proof First we consider that the system has non-zero column minimal indices and non-zero row minimal indices. By using the transformation (34), the solutions of the subsystems (29), (30), (31), (32) and (33) are given by (10), (12), (36), (40) and (44) respectively. Note that from Proposition 2.1 the solution (10) exists if and only if

∥ J p ∥<1.

Furthermore, if

Z k = [ Z k p Z k q Z k ϵ Z k ζ Z k g ] = [ ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p ) Z 0 p 0 q , 1 C k , 1 0 t − h , 1 C k , 2 ] .

Since Z 0 p is unknown, it can be replaced with the unknown vector C. Then

Y k =Q Z k =[ Q p Q q Q ϵ Q ζ Q g ] [ ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p ) C 0 q , 1 C k , 1 0 t − h , 1 C k , 2 ]

or, equivalently,

Y k = Q p ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p )C+ Q ϵ C k , 1 + Q g C k , 2 .

Since C k , 1 and C k , 2 can be taken arbitrarily, it is clear that the general singular discrete time system for every suitable defined boundary condition has an infinite number of solutions. It is clear that the existence of the column minimal indices is the reason that the systems (31) and consequently (33) exist. These systems as shown in Propositions 3.1 and 3.3 have always infinite solutions. Thus a necessary condition for the system to have a unique solution is not to have any column minimal indices which are equal to

dim N r (sF−G)=0.

In this case the Kronecker canonical form of the pencil sF−G has the following form:

s F K − G K :=s I p − J p ⊕s H q − I q ⊕s F ζ − G ζ .
(51)

Then the system (1) is divided into three subsystems (29), (30), (32) with solutions (10), (12), (40) respectively. Thus

Y k =Q Z k =[ Q p Q q Q ζ ] [ ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p ) C 0 q , 1 0 t − h , 1 ]

or, equivalently,

Y k = Q p ( k + 1 ) n − 1 ¯ F n , n ( J p ( k + n ) n ¯ ) ( I p − J p )C.

The solution exists if and only if

A 1 Y 0 = B 1 , A 2 Y N = B 2

or, equivalently,

A 1 Q p C= B 1 , A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p )C= B 2 ,

or, equivalently,

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] C= [ B 1 B 2 ] .

For the above algebraic system, there exists at least one solution if and only if

[ B 1 B 2 ] ∈colspan [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] .

The algebraic system (50) contains r 1 + r 2 equations and p unknowns. Hence the solution is unique if and only if

p≤ r 1 + r 2

and

rank [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] =p,

where C is then the unique solution of (50). The uniqueness of C can be proved as follows. If we assume that the algebraic system has two solutions C 1 and C 2 , then

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] C 1 = [ B 1 B 2 ]

and

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] C 2 = [ B 1 B 2 ]

or, equivalently,

[ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] ( C 1 − C 2 )= 0 p , 1 .

But the matrix [ A 1 Q p A 2 Q p ( N + 1 ) n − 1 ¯ F n , n ( J p ( N + n ) n ¯ ) ( I p − J p ) ] is left invertible since it is assumed to have p linear independent columns and r 1 + r 2 ≥p and hence

C 1 = C 2 .

The unique solution is then given from (49). The proof is completed. □

4 Numerical examples

Example 1

Assume the system (1) for k=0,1,2,3 and n= 3 2 . Let

F= [ 0 1 0 0 0 0 0 0 0 − 1 1 − 1 0 1 − 1 1 − 1 0 2 − 1 0 0 0 0 0 ]
(52)

and

G= [ 0 0 1 − 1 1 0 0 0 0 − 1 4 1 2 − 1 2 0 − 1 2 1 2 1 2 − 1 2 1 − 1 2 3 4 0 0 0 1 − 1 ] .
(53)

Then det(sF−G)=(s− 1 2 )s(s− 1 4 ) and the pencil is regular. We assume the boundary conditions (2) with

A 1 = [ 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 1 ] , A 2 = [ 1 0 0 0 1 ] , B 1 = [ 0 0 36 36 ]

and

B 2 =24.

The three finite eigenvalues (p=3) of the pencil are 1 2 , 0, 1 4 , and the Jordan matrix J p has the form

J p = [ 1 2 0 0 0 0 0 0 0 1 4 ] .

It is easy to observe that

∥ J p ∥<1.

By calculating the eigenvectors of the finite eigenvalues, we get the matrix Q p

Q p = [ 1 1 0 0 1 0 0 0 0 0 0 1 0 0 1 ] .

Moreover,

A 1 Q p = [ 1 1 0 0 1 0 0 0 1 0 0 1 ] .
(54)

From (9) we get the Mittag-Leffler function

F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) = ∑ i = 0 ∞ J p i ( 3 + 3 2 ) i 3 2 ¯ Γ ( ( i + 1 ) 3 2 )

or, equivalently,

F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) = ∑ i = 0 ∞ J p i Γ ( 3 + 3 2 ( i + 1 ) ) Γ ( 3 + 3 2 ) Γ ( ( i + 1 ) 3 2 ) ,

or, equivalently,

F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) = ∑ i = 0 ∞ J p i ( 2 + 3 2 ( i + 1 ) ) ( 1 + 3 2 ( i + 1 ) ) ( 3 2 ( i + 1 ) ) Γ ( 3 2 ( i + 1 ) ) Γ ( 3 + 3 2 ) Γ ( ( i + 1 ) 3 2 ) ,

or, equivalently,

F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) = 1 Γ ( 3 + 3 2 ) ∑ i = 0 ∞ J p i ( 2 + 3 2 ( i + 1 ) ) ( 1 + 3 2 ( i + 1 ) ) ( 3 2 ( i + 1 ) )

and since ∥ J p ∥<1, by using the sum ∑ i = 0 ∞ x= ( 1 − x ) − 1 for ∥x∥<1, we calculate the sum of the matrix power series ∑ i = 0 ∞ J p i (2+ 3 2 (i+1))(1+ 3 2 (i+1))( 3 2 (i+1)), and we get

F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) = 1 Γ ( 3 + 3 2 ) ( − 1 9 J p 2 + 20 9 J p + 35 9 I p ) ( I p − J p ) − 4 .

And since

A 2 Q p =[ 1 1 1 ],

by using the above expression

A 2 Q p ( 4 ) 1 2 ¯ F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) ( I p − J p )= 1 36 [ 1 1 1 ] [ 358 0 0 0 35 0 0 0 24 ]

or, equivalently,

A 2 Q p ( 4 ) 1 2 ¯ F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) ( I p − J p )= 1 36 [ 358 35 24 ],
(55)

it is easy to observe that the conditions (13), (14) and (15) are satisfied. Thus from Theorem 2.1 the unique solution of the boundary value problem (1)-(2) is given by

Y k = Q p ( k + 1 ) 3 2 − 1 ¯ F 3 2 , 3 2 ( J p ( k + 3 2 ) 3 2 ¯ ) ( I p − J p )C

or, equivalently, by

Y k = 1 Γ ( k + 1 ) ∑ i = 0 ∞ Γ ( k + 3 2 ( i + 1 ) ) Γ ( 3 2 ( i + 1 ) ) [ 1 2 i + 1 0 0 0 0 0 0 0 0 0 0 3 4 i + 1 0 0 3 4 i + 1 ] C,

where C is the unique solution of the algebraic system

[ A 1 Q p A 2 Q p ( 4 ) 1 2 ¯ F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) ( I p − J p ) ] C= [ B 1 B 2 ]

or, equivalently,

C= [ 0 0 36 ] ,

and thus the unique solution of the boundary value problem is

Y k = 27 Γ ( k + 1 ) ∑ i = 0 ∞ Γ ( k + 3 2 ( i + 1 ) ) Γ ( 3 2 ( i + 1 ) ) [ 0 0 0 1 4 i 1 4 i ] .
(56)

Example 2

We assume the system (1) as in Example 1 but with different boundary conditions. Let

A 1 = [ 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 1 ] , A 2 = [ 1 0 0 0 1 ] , B 1 = [ 0 0 0 36 ]

and

B 2 =24.

It is easy to observe that

[ B 1 B 2 ] ∉colspan [ A 1 Q p A 2 Q p ( 4 ) 1 2 ¯ F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) ( I p − J p ) ]

since

[ 0 0 0 36 24 ] ∉colspan [ 36 36 0 0 36 0 0 0 36 0 0 36 358 35 24 ]

and thus from Theorem 2.1, and since (13) does not hold, the boundary value problem is not consistent.

Example 3

Consider the system (1) and let

F= [ 1 1 1 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 0 1 1 0 1 0 0 0 0 1 1 1 ] ,G= [ 1 2 2 1 2 0 2 2 0 2 1 2 2 2 3 0 2 3 1 3 0 0 0 0 0 1 0 1 0 0 ] .

Since the matrices F, G are non-square, the matrix pencil sF−G is singular and has invariants such as the finite elementary divisors s−2, s−1, an infinite elementary divisor of degree 1 and the row minimal indices ζ 1 =0, ζ 2 =1. Since the Jordan matrix has the form

J p = [ 2 0 0 1 ]

with

∥ J p ∥>1

for every induced matrix norm, from Theorem 3.1 the boundary value problem (1)-(2) is non-consistent.

Example 4

Consider the system (1) for k=0,1,2,3 and n= 3 2 . Let

F= [ 1 0 0 1 0 0 ]
(57)

and

G= [ 3 4 0 0 0 0 1 ] .
(58)

Since the matrices F, G are non-square, the matrix pencil sF−G is singular and has invariants such as a finite elementary divisor s− 3 4 and the row minimal indices ζ 1 =0, ζ 2 =1. We assume the boundary conditions (2) with

A 1 =[ 1 1 ], A 2 = [ 0 1 0 2 ] , B 1 =2

and

B 2 = [ 0 0 ] .

The Jordan matrix is J p = 3 4 with ∥ J p ∥<1 for every induced matrix norm. By calculating the matrix Q p , we get

Q p = [ 1 0 ] .

Moreover,

A 1 Q p =1,
(59)

and since

A 2 Q p = [ 0 0 ] ,

we get

A 2 Q p ( 4 ) 1 2 ¯ F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) ( I p − J p )= [ 0 0 ] .
(60)

By using (59), (60), it is easy to observe that the conditions (45), (46), (47) and (48) are satisfied and thus from Theorem 3.1 the unique solution of the boundary value problem (1)-(2) is given by

Y k = Q p ( k + 1 ) 1 2 ¯ F 3 2 , 3 2 ( J p ( k + 3 2 ) 3 2 ¯ ) ( I p − J p )C

or, equivalently, by

Y k = 1 Γ ( k + 1 ) ∑ i = 0 ∞ Γ ( k + 3 2 ( i + 1 ) ) Γ ( 3 2 ( i + 1 ) ) [ 3 i 4 i + 1 0 ] C,

where C is the unique solution of the algebraic system

[ A 1 Q p A 2 Q p ( 4 ) 1 2 ¯ F 3 2 , 3 2 ( J p ( 3 + 3 2 ) 3 2 ¯ ) ( I p − J p ) ] C= [ B 1 B 2 ]

or, equivalently,

C=2,

and thus the unique solution of the boundary value problem is

Y k = 2 Γ ( k + 1 ) ∑ i = 0 ∞ Γ ( k + 3 2 ( i + 1 ) ) Γ ( 3 2 ( i + 1 ) ) [ 3 i 4 i + 1 0 ] .
(61)

5 Conclusions

In this article, we study the boundary value problem of a class of a singular system of fractional nabla difference equations whose coefficients are constant matrices. By taking into consideration the cases that the matrices are square with the leading coefficient singular, square with an identically zero matrix pencil and non-square, we study the conditions under which the boundary value problem has unique, infinite and no solutions. Furthermore, we provide a formula for the case of the unique solution. As a further extension of this article, one can study the stability, the behavior under perturbation and possible applications in economics and engineering of singular matrix difference/differential equations of fractional order. For all this, there is already some research in progress.

References

  1. Baleanu D, Diethelm K, Scalas E: Fractional Calculus: Models and Numerical Methods. World Scientific, Singapore; 2012.

    MATH  Google Scholar 

  2. Glockle WG, Nonnenmacher TF: A fractional calculus approach to self-similar protein dynamics. Biophys. J. 1995, 68(1):46–53. 10.1016/S0006-3495(95)80157-8

    Article  Google Scholar 

  3. Hilfe R (Ed): Applications of Fractional Calculus in Physics. World Scientific, River Edge; 2000.

    Google Scholar 

  4. Kaczorek T 411. In Selected Problems of Fractional Systems Theory. Springer, Berlin; 2011.

    Chapter  Google Scholar 

  5. Malinowska AB, Torres DFM: Introduction to the Fractional Calculus of Variations. Imperial College Press, London; 2012. ISBN:978–1-84816–966–1/hbk

    Book  MATH  Google Scholar 

  6. Metzler R, Schick W, Kilian HG, Nonnenmacher TF: Relaxation in filled polymers: a fractional calculus approach. J. Chem. Phys. 1995, 103(16):7180–7186. 10.1063/1.470346

    Article  Google Scholar 

  7. Podlubny I Mathematics in Science and Engineering. In Fractional Differential Equations. Academic Press, San Diego; 1999.

    Google Scholar 

  8. Klamka J: Controllability and minimum energy control problem of fractional discrete-time systems. In New Trends in Nanotechnology and Fractional Calculus. Springer, New York; 2010:503–509.

    Chapter  Google Scholar 

  9. Agarwal RP, Benchohra M, Hamani S: A survey on existence results for boundary value problems of nonlinear fractional differential equations and inclusions. Acta Appl. Math. 2010, 109(3):973–1033. 10.1007/s10440-008-9356-6

    Article  MATH  MathSciNet  Google Scholar 

  10. Ahrendt, K, Castle, L, Holm, M, Yochman, K: Laplace transforms for the nabla-difference operator and a fractional variation of parameters formula. Commun. Appl. Anal. (2011)

    MATH  Google Scholar 

  11. Atici FM, Eloe PW: Linear systems of fractional nabla difference equations. Rocky Mt. J. Math. 2011, 41(2):353–370. 10.1216/RMJ-2011-41-2-353

    Article  MATH  MathSciNet  Google Scholar 

  12. Atici FM, Eloe PW: Initial value problems in discrete fractional calculus. Proc. Am. Math. Soc. 2009, 137(3):981–989.

    Article  MATH  MathSciNet  Google Scholar 

  13. Baleanu D, Mustafa OG, Agarwal RP: Asymptotically linear solutions for some linear fractional differential equations. Abstr. Appl. Anal. 2010., 2010: Article ID 865139

    Google Scholar 

  14. Baleanu D, Mustafa OG: On the global existence of solutions to a class of fractional differential equations. Comput. Math. Appl. 2010, 59(5):1835–1841. 10.1016/j.camwa.2009.08.028

    Article  MATH  MathSciNet  Google Scholar 

  15. Baleanu D, Babakhani A: Employing of some basic theory for class of fractional differential equations. Adv. Differ. Equ. 2011., 2011: Article ID 296353

    Google Scholar 

  16. Bastos NRO, Ferreira RAC, Torres DFM: Necessary optimality conditions for fractional difference problems of the calculus of variations. Discrete Contin. Dyn. Syst. 2011, 29(2):417–437.

    Article  MATH  MathSciNet  Google Scholar 

  17. Bastos NRO, Ferreira RAC, Torres DFM: Discrete-time fractional variational problems. Signal Process. 2011, 91(3):513–524. 10.1016/j.sigpro.2010.05.001

    Article  MATH  Google Scholar 

  18. Bastos NRO, Mozyrska D, Torres DFM: Fractional derivatives and integrals on time scales via the inverse generalized Laplace transform. Int. J. Math. Comput. 2011, 11: 1–9.

    MathSciNet  Google Scholar 

  19. Debbouche A: Fractional nonlocal impulsive quasilinear multi-delay integro-differential systems. Adv. Differ. Equ. 2011., 2011: Article ID 5

    Google Scholar 

  20. Debbouche A, Baleanu D: Controllability of fractional evolution nonlocal impulsive quasilinear delay integro-differential systems. Comput. Math. Appl. 2011, 62(3):1442–1450. 10.1016/j.camwa.2011.03.075

    Article  MATH  MathSciNet  Google Scholar 

  21. Debbouche A, Baleanu D, Agarwal RP: Nonlocal nonlinear integrodifferential equations of fractional orders. Adv. Differ. Equ. 2012., 2012: Article ID 78

    Google Scholar 

  22. Ferreira RAC, Torres DFM: Fractional h -difference equations arising from the calculus of variations. Appl. Anal. Discrete Math. 2011, 5(1):110–121. 10.2298/AADM110131002F

    Article  MATH  MathSciNet  Google Scholar 

  23. Hein J, McCarthy Z, Gaswick N, McKain B, Speer K: Laplace transforms for the nabla-difference operator. Panam. Math. J. 2011, 21(3):79–97.

    MATH  MathSciNet  Google Scholar 

  24. Kaczorek T: Positive stable realizations of fractional continuous-time linear systems. Int. J. Appl. Math. Comput. Sci. 2011, 21(4):697–702.

    Article  MATH  MathSciNet  Google Scholar 

  25. Kaczorek T: Application of the Drazin inverse to the analysis of descriptor fractional discrete-time linear systems with regular pencils. Int. J. Appl. Math. Comput. Sci. 2013, 23(1):29–33.

    Article  MATH  MathSciNet  Google Scholar 

  26. Klamka J: Controllability of dynamical systems. Mat. Stosow. 2008, 50(9):57–75.

    MathSciNet  MATH  Google Scholar 

  27. Klamka J, Wyrwał J: Controllability of second-order infinite-dimensional systems. Syst. Control Lett. 2008, 57(5):386–391. 10.1016/j.sysconle.2007.10.002

    Article  MATH  MathSciNet  Google Scholar 

  28. Dai L Lecture Notes in Control and Information Sciences. Singular Control Systems 1988. (edited by M Thoma and A Wyner)

    Google Scholar 

  29. Dassios IK: On non homogeneous linear generalized linear discrete time systems. Circuits Syst. Signal Process. 2012, 31(5):1699–1712. 10.1007/s00034-012-9400-7

    Article  MATH  MathSciNet  Google Scholar 

  30. Dassios IK, Kalogeropoulos G: On a non-homogeneous singular linear discrete time system with a singular matrix pencil. Circuits Syst. Signal Process. 2013. doi:10.1007/s00034–012–9541–8

    Google Scholar 

  31. Dassios I: On solutions and algebraic duality of generalized linear discrete time systems. Discrete Math. Appl. 2012, 22(5–6):665–682.

    Article  MATH  MathSciNet  Google Scholar 

  32. Dassios I: On stability and state feedback stabilization of singular linear matrix difference equations. Adv. Differ. Equ. 2012., 2012: Article ID 75

    Google Scholar 

  33. Dassios I: On robust stability of autonomous singular linear matrix difference equations. Appl. Math. Comput. 2012, 218(12):6912–6920. 10.1016/j.amc.2011.12.067

    Article  MATH  MathSciNet  Google Scholar 

  34. Dassios I: On a boundary value problem of a class of generalized linear discrete time systems. Adv. Differ. Equ. 2011., 2011: Article ID 51

    Google Scholar 

  35. Grispos E, Giotopoulos S, Kalogeropoulos G: On generalised linear discrete-time regular delay systems. J. Inst. Math. Comput. Sci., Math. Ser. 2000, 13(2):179–187.

    MATH  MathSciNet  Google Scholar 

  36. Grispos E, Kalogeropoulos G, Mitrouli M: On generalised linear discrete-time singular delay systems. In Proceedings of the 5th Hellenic-European Conference on Computer Mathematics and Its Applications (HERCMA 2001). Edited by: Lipitakis EA. LEA, Athens; 2002:484–486. 2 volumes

    Google Scholar 

  37. Grispos E, Kalogeropoulos G, Stratis I: On generalised linear discrete-time singular delay systems. J. Math. Anal. Appl. 2000, 245(2):430–446. 10.1006/jmaa.2000.6761

    Article  MATH  MathSciNet  Google Scholar 

  38. Grispos E: Singular generalised autonomous linear differential systems. Bull. Greek Math. Soc. 1992, 34: 25–43.

    MATH  MathSciNet  Google Scholar 

  39. Kalogeropoulos, GI: Matrix pencils and linear systems. PhD thesis, City University, London (1985)

    Google Scholar 

  40. Kalogeropoulos G, Stratis IG: On generalized linear regular delay systems. J. Math. Anal. Appl. 1999, 237(2):505–514. 10.1006/jmaa.1999.6458

    Article  MATH  MathSciNet  Google Scholar 

  41. Karcanias N, Kalogeropoulos G: Geometric theory and feedback invariants of generalized linear systems: a matrix pencil approach. Circuits Syst. Signal Process. 1989, 8(3):375–397. 10.1007/BF01598421

    Article  MATH  MathSciNet  Google Scholar 

  42. Rugh WJ: Linear System Theory. Prentice Hall International, London; 1996.

    MATH  Google Scholar 

  43. Sandefur JT: Discrete Dynamical Systems. Academic Press, San Diego; 1990.

    MATH  Google Scholar 

  44. Gantmacher RF: The Theory of Matrices. Vols. I, II. Chelsea, New York; 1959.

    MATH  Google Scholar 

  45. Mitrouli M, Kalogeropoulos G: A compound matrix algorithm for the computation of the Smith form of a polynomial matrix. Numer. Algorithms 1994, 7(2–4):145–159.

    Article  MATH  MathSciNet  Google Scholar 

  46. Nagai A: Discrete Mittag-Leffler function and its applications. Publ. Res. Inst. Math. Sci. 2003, 1302: 1–20.

    MathSciNet  Google Scholar 

Download references

Acknowledgements

We would like to express our sincere gratitude to Professor GI Kalogeropoulos for his helpful and fruitful discussions that clearly improved this article. Moreover, we are very grateful to the anonymous referees for their valuable suggestions that improved the article.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ioannis K Dassios.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

IKD wrote the first draft of the manuscript and DB correct it and prepared the final version of it. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Dassios, I.K., Baleanu, D.I. On a singular system of fractional nabla difference equations with boundary conditions. Bound Value Probl 2013, 148 (2013). https://doi.org/10.1186/1687-2770-2013-148

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-148

Keywords