Optimal control of dengue fever model with a logistically growing human population

1.   Introduction

Throughout this paper, the complex $m \times n$ matrix space is denoted by ${\mathbb{C}}^{m \times n}$. The conjugate transpose, the Moore-Penrose inverse, the range space and the null space of a complex matrix $A \in {\mathbb{C}}^{m \times n}$ are denoted by $A^\mathrm{H}$, $A^+$, $\mathcal{R}(A)$ and $\mathcal{N}(A)$, respectively. $I_n$ denotes the $n\times n$ identity matrix. $P_{\mathcal{L}}$ stands for the orthogonal projector on the subspace $\mathcal{L}\subset \mathbb{C}^n.$ Furthermore, for a matrix $A \in{\mathbb{C}}^{m \times n},$ $E_A$ and $F_A$ stand for two orthogonal projectors: $E_A = I_m-AA^+ = P_{\mathcal{N}(A^\mathrm{H})}$, $F_A = I_n-A^+A = P_{\mathcal{N}(A)}.$

A number of papers have been published for solving linear matrix equations. For example, Chen et al. ^[1] proposed LSQR iterative method to solve common symmetric solutions of matrix equations $AXB = E$ and $CXD = F$. Zak and Toutounian ^[2] studied the matrix equation of $AXB = C$ with nonsymmetric coefficient matrices by using nested splitting conjugate gradient (NSCG) iteration method. By applying a Hermitian and skew-Hermitian splitting (HSS) iteration method, Wang et al. ^[3] computed the solution of the matrix equation $AXB = C$. Tian et al. ^[4] obtained the solution of the matrix equation $AXB = C$ by applying the Jacobi and Gauss-Seidel-type iteration methods. Liu et al. ^[5] solved the matrix equation $AXB = C$ by employing stationary splitting iterative methods. In addition, some scholars studied matrix equations by direct methods. Yuan and Dai ^[6] obtained generalized reflexive solutions of the matrix equation $AXB = D$ and the optimal approximation solution by using the generalized singular value decomposition. Zhang et al. ^[7] provided the explicit expression of the minimal norm least squares Hermitian solution of the complex matrix equation $AXB+CXD = E$ by using the structure of the real representations of complex matrices and the Moore-Penrose inverse. By means of the definitions of the rank and inertias of matrices, Song and Yu ^[8] obtained the existence conditions and expressions of the nonnegative (positive) definite and the Re-nonnegative (Re-positive) definite solutions to the matrix equations $AXA^{\mathrm{H}} = C$ and $BXB^{\mathrm{H}} = D$.

In this paper, we will focus on the restricted solutions to the following well-known linear matrix equation

where $A \in {\mathbb{C}}^{m \times n}, B \in {\mathbb{C}}^{n \times p}$ and $C \in {\mathbb{C}}^{m \times p}.$ We observe that the structured matrix, such as Hermitian matrix, skew-Hermitian matrix, Re-nonnegative definite matrix and Re-positive definite matrix, is of important applications in structural dynamics, numerical analysis, stability and robust stability analysis of control theory and so on ^{[9,10,11,12,13,14]}. Conditions for the existence of Hermitian solutions to Eq $(1)$ were studied in ^[15,16,17]. A solvability criterion for the existence of Re-nonnegative definite solutions of Eq $(1)$ by using generalized inner inverses was investigated by Cvetković-Ilić ^[19]. Recently, a direct method for solving Eq $(1)$ by using the generalized inverses of matrices and orthogonal projectors was proposed by Yuan and Zuo ^[21]. In addition, the Re-nonnegative definite and Re-positive definite solutions to some special cases of Eq $(1)$ were considered by Wu ^[22], Wu and Cain ^[23] and Groß ^[24]. In ^[25], necessary and sufficient conditions for the existence of common Re-nonnegative definite and Re-positive definite solutions to the matrix equations $AX = C, XB = D$ were discussed by virtue of the extremal ranks of matrix polynomials.

In this paper, necessary and sufficient conditions for the existence of Hermitian (skew-Hermitian), Re-nonnegative (Re-positive) definite, and Re-nonnegative (Re-positive) definite least-rank solutions to Eq $(1)$ are deduced by using the Moore-Penrose inverse of matrices, and the explicit representations of the general solutions are given when the solvability conditions are satisfied. Compared with the approaches proposed in ^{[18,19,20,21]}, the coefficient matrices of Eq $(1)$ have no any constraints and the method in this paper is straightforward and easy to implement.

2.   Preliminaries

Definition 1. A matrix $A \in {\mathbb{C}}^{n \times n}$ is said to be Re-nonnegative definite $($Re-nnd$)$ if $H(A): = \frac{1}{2}(A+A^\mathrm{H})$ is Hermitian nonnegative definite $($i.e., $H(A)\geq 0)$, and $A$ is said to be Re-positive definite $($Re-pd$)$ if $H(A)$ is Hermitian positive definite $($i.e., $H(A) > 0)$. The set of all Re-nnd $($Re-pd$)$ matrices in ${\mathbb{C}}^{n \times n}$ is denoted by $\mathbb{RND}^{n \times n}$ $(\mathbb{RPD}^{n \times n})$.

Lemma 1. ^[26] Let $A \in {\mathbb{C}}^{m \times n}, B \in {\mathbb{C}}^{n \times p}$ and $C \in {\mathbb{C}}^{m \times p}.$ Then the matrix equation $AXB = C$ is solvable if and only if $AA^+CB^+B = C.$ In this case, the general solution can be written in the following parametric form

where $L_1, L_2 \in {\mathbb{C}}^{n \times n}$ are arbitrary matrices.

Lemma 2. ^[27,28] Let $B_1 \in {\mathbb{C}}^{l \times q}, D_1 \in {\mathbb{C}}^{l \times l}.$ Then the matrix equation

has a solution $Y \in {\mathbb{C}}^{l \times q}$ if and only if

In which case, the general solution is

where $V \in {\mathbb{C}}^{l \times q}$ is an arbitrary matrix.

Lemma 3. ^[29] Let $A \in {\mathbb{C}}^{m \times n}$ and $B \in {\mathbb{C}}^{m \times p}$. Then the Moore-Penrose inverse of the matrix $[A, B]$ is

where $C = (I-AA^{+})B$ and $T = A^{+}B(I-C^{+}C).$

Lemma 4. ^[30] Suppose that a Hermitian matrix $M$ is partitioned as

where $M_{11}, \ M_{22}$ are square. Then

$(\mathrm{i}).$ $M$ is Hermitian nonnegative definite if and only if

In which case, $M$ can be expressed as

where $H_1$ is an arbitrary matrix and $H_2$ is an arbitrary Hermitian nonnegative definite matrix.

$(\mathrm{ii}).$ $M$ is Hermitian positive definite if and only if

In the case, $M$ can be expressed as

where $H_3$ is an arbitrary Hermitian positive definite matrix.

Lemma 5. ^[31] Let

Then the general least-rank solution to Eq.$(1)$ can be written as

where $R_1 \in {\mathbb{C}}^{(p+n) \times n}, R_2 \in {\mathbb{C}}^{n \times (m+n)}$ are arbitrary matrices.

3.   The Hermitian and skew-Hermitian solutions of Eq (1)

Theorem 1. Eq $(1)$ has a Hermitian solution $X$ if and only if

where $\mathcal{T} = \mathcal{R}(A^\mathrm{H})\cap \mathcal{R}(B).$ In which case, the general Hermitian solution of Eq $(1)$ is

where

and $V_1, V_2 \in {\mathbb{C}}^{n \times n}$ are arbitrary matrices.

Proof. By Lemma $1$, if the first condition of $(2)$ holds, then the general solution of Eq $(1)$ is given by $(3)$. Now, we will find $L_1$ and $L_2$ such that $AXB = C$ has a Hermitian solution, that is,

Clearly, Eq $(6)$ can be equivalently written as

where $A_1 = [F_A, -E_B], \ X_1 = [L_1^\mathrm{H}, L_2], \ D_1 = A^+CB^+-(A^+CB^+)^\mathrm{H}.$

By Lemma $2$, Eq $(7)$ has a solution $X_1$ if and only if

The first condition of $(8)$ is obviously satisfied. And note that

Thus, the second condition of $(8)$ is equivalent to

where $\mathcal{T} = \mathcal{R}(A^\mathrm{H})\cap \mathcal{R}(B),$ which is the second condition of $(2)$. In which case, the general solution of Eq $(7)$ is

where $V = [V_1, V_2]$ is an arbitrary matrix. By Lemma $3$, we have

where $C_1 = -(I-F_AF_A^+)E_B$ and $T_1 = -F_{A}^+E_{B}(I-C_1^+C_1).$ Upon substituting $(10)$ into $(9)$, we can get $(4)$ and $(5)$.

Corollary 1. Eq $(1)$ has a skew-Hermitian solution $X$ if and only if

where $\mathcal{T} = \mathcal{R}(A^\mathrm{H})\cap \mathcal{R}(B).$ In which case, the general skew-Hermitian solution of Eq $(1)$ is

where

and $V_3, \ V_4 \in {\mathbb{C}}^{n \times n}$ are arbitrary matrices.

4.   The Re-nnd and Re-pd solutions of Eq (1)

Theorem 2. Let $A \in {\mathbb{C}}^{m \times n}, \ B \in {\mathbb{C}}^{n \times p}$, $C \in {\mathbb{C}}^{m \times p}$ and $\mathcal{T} = \mathcal{R}(A^\mathrm{H})\cap \mathcal{R}(B).$ Assume that the spectral decomposition of $P_\mathcal{T}$ is

where $U = [U_1, U_2]\in \mathbb{C}^{n \times n}$ is a unitary matrix and $s = \mathit{\mbox{dim}}(\mathcal{T}).$ Then

$(a).$ Eq $(1)$ has a Re-nnd solution if and only if

In which case, the general Re-nnd solution of $(1)$ is

where

$V_5, V_6 \in {\mathbb{C}}^{n \times n}, H_1 \in {\mathbb{C}}^{s \times (n-s)}$ are arbitrary matrices, and $H_2 \in {\mathbb{C}}^{(n-s) \times (n-s)}$ is an arbitrary Hermitian nonnegative definite matrix.

$(b).$ Eq $(1)$ has a Re-pd solution if and only if

In which case, the general Re-pd solution of $(1)$ is

where

$J_1, J_2, D_3, C_3, T_3, P_3, Q_3, W_3, Z_3$ and $K_{11}$ are given by $(18)$ and $(19),$ $K_{12} \in {\mathbb{C}}^{s \times (n-s)}$ is an arbitrary matrix and $H_3 \in {\mathbb{C}}^{(n-s) \times (n-s)}$ is an arbitrary Hermitian positive definite matrix.

Proof. By Lemma $1$, if the first condition of $(16)$ holds, then the general solution of Eq $(1)$ is given by $(17)$. Now, we will find $J_1$ and $J_2$ such that $AXB = C$ has a Re-nnd (Re-pd) solution, that is, we will choose suitable matrices $J_1$ and $J_2$ such that

Clearly, Eq $(22)$ can be equivalently written as

where $A_3 = [F_A, E_B], \ X_3 = [J_1^\mathrm{H}, J_2], \ D_3 = K-A^+CB^+-(A^+CB^+)^\mathrm{H}.$

By Lemma $2$, Eq $(23)$ has a solution $X_1$ if and only if

The first condition of $(24)$ is obviously satisfied. And note that

Then the second condition of $(24)$ is equivalent to

where $\mathcal{T} = \mathcal{R}(A^\mathrm{H})\cap \mathcal{R}(B).$ By $(15),$ we can obtain

Let

where $U = [U_1, U_2]$. By $(25),$ we can obtain

it is easily known that $K\geq 0$ $(K > 0)$ if and only if $U^\mathrm{H}KU\geq 0$ $(U^\mathrm{H}KU > 0).$ And $X$ is Re-nnd (Re-pd) if and only if $K\geq 0$ $(K > 0).$ Thus, by $(26)$ and $(27)$, we can get

equivalently,

which are the second conditions of $(16)$ and $(20)$. In which case, by Lemma 4,

where $H_1 \in {\mathbb{C}}^{s \times (n-s)}$ is an arbitrary matrix, $H_2 \in {\mathbb{C}}^{(n-s) \times (n-s)}$ is an arbitrary Hermitian nonnegative definite matrix and $H_3 \in {\mathbb{C}}^{(n-s) \times (n-s)}$ is an arbitrary Hermitian positive definite matrix. And the general solution of Eq $(23)$ is

where $V = [V_5, V_6]$ is an arbitrary matrix. By Lemma $3$, we have

where $C_3 = (I-F_AF_A^+)E_B$ and $T_3 = F_{A}^+E_{B}(I-C_3^+C_3).$ Upon substituting $(29)$ into $(28)$, we can get $(18)$ and $(19)$.

5.   The Re-nnd and Re-pd least-rank solutions of Eq (1)

Theorem 3. Let $A \in {\mathbb{C}}^{m \times n}, B \in {\mathbb{C}}^{n \times p},$ $C \in {\mathbb{C}}^{m \times p}$ and $\mathcal{\tilde{T}} = \mathcal{N}(N_1^\mathrm{H})\cap \mathcal{N}(S_1).$ Assume that the spectral decomposition of $P_\mathcal{\tilde{T}}$ is

where $\tilde{U} = [\tilde{U}_1, \tilde{U}_2]\in \mathbb{C}^{n \times n}$ is a unitary matrix and $k = \mathit{\mbox{dim}}(\mathcal{\tilde{T}})$, and $M, N, S, N_1, S_1$ are given by Lemma $5$. Then

$(a).$ Eq $(1)$ has a Re-nnd least-rank solution if and only if

In which case, the general Re-nnd least-rank solution of Eq $(1)$ is

where

$V_7 \in {\mathbb{C}}^{n \times (p+n)},$ $V_8 \in {\mathbb{C}}^{n \times (m+n)},$ $\tilde{H}_1 \in {\mathbb{C}}^{k \times (n-k)}$ are arbitrary matrices, and $\tilde{H}_2 \in {\mathbb{C}}^{(n-k) \times (n-k)}$ is an arbitrary Hermitian nonnegative definite matrix.

$(b).$ Eq $(1)$ has a Re-pd least-rank solution if and only if

In which case, the general Re-nnd least-rank solution of Eq $(1)$ is

where

$R_1, R_2, D_4, C_4, T_4, P_4, Q_4, W_4, Z_4$ and $\tilde{K}_{11}$ are given by $(33)$ and $(34),$ $\tilde{K}_{12} \in {\mathbb{C}}^{k \times (n-k)}$ is an arbitrary matrix and $\tilde{H}_3 \in {\mathbb{C}}^{(n-k) \times (n-k)}$ is an arbitrary Hermitian positive definite matrix.

Proof. By Lemmas $1$ and $5$, if the first condition of $(31)$ holds, then the least-rank solution of Eq $(1)$ is given by $(32)$. Now, we will find $R_1$ and $R_2$ such that $AXB = C$ has a Re-nnd (Re-pd) least-rank solution, that is, we will choose suitable matrices $R_1$ and $R_2$ such that

Clearly, Eq $(37)$ can be equivalently written as

where $A_4 = [N_1, S_1^\mathrm{H}], \ X_4 = [R_1^\mathrm{H}, R_2], \ D_4 = \tilde{K}+NM^+S+(NM^+S)^\mathrm{H}.$

By Lemma $2$, Eq $(38)$ has a solution $X_4$ if and only if

The first condition of $(39)$ is obviously satisfied. And note that

Thus, the second condition of $(39)$ is equivalent to

where $\mathcal{\tilde{T}} = \mathcal{N}(N_1^\mathrm{H})\cap \mathcal{N}(S_1).$ By $(30),$ we can obtain

Let

where $\tilde{U} = [\tilde{U}_1, \tilde{U}_2]$. By $(40),$ we can obtain

it is easily known that $\tilde{K}\geq 0$ $(\tilde{K} > 0)$ if and only if $\tilde{U}^\mathrm{H}\tilde{K}\tilde{U}\geq 0$ $(\tilde{U}^\mathrm{H}\tilde{K}\tilde{U} > 0).$ And $X$ is Re-nnd (Re-pd) least-rank solution if and only if $\tilde{K}\geq 0$ $(\tilde{K} > 0).$ Thus, by $(41)$ and $(42)$, we can get

equivalently,

which are the second conditions of $(31)$ and $(35)$. In which case, by Lemma 4,

where $\tilde{H}_1 \in {\mathbb{C}}^{k \times (n-k)}$ is an arbitrary matrix, $\tilde{H}_2 \in {\mathbb{C}}^{(n-k) \times (n-k)}$ is an arbitrary Hermitian nonnegative definite matrix and $\tilde{H}_3 \in {\mathbb{C}}^{(n-k) \times (n-k)}$ is an arbitrary Hermitian positive definite matrix. And the general solution of Eq $(38)$ is

where $V = [V_7, V_8]$ is an arbitrary matrix. By Lemma $3$, we have

where $C_4 = (I-N_1N_1^+)S_1^\mathrm{H}$ and $T_4 = N_1^+S_1^\mathrm{H}(I-C_4^+C_4).$ Upon substituting $(44)$ into $(43)$, we can get $(33)$ and $(34)$.

6.   Numerical examples

The following example comes from ^[9].

Example 1. Consider a $7$-DOF system modelled analytically with the first three measured modal data given by

and the corrected symmetric mass matrix $M$ and symmetric stiffness matrix $K$ should satisfy the orthogonality conditions, that is,

Since

which means that the conditions of (2) are satisfied. Choose $V_1 = 0, V_2 = 0.$ Then, by the equation of (3), we can obtain a corrected mass matrix given by

and

which implies that $M$ is a symmetric solution of $X^\top MX = I_3$.

Since

which means that the conditions of (2) are satisfied. Choose $V_1 = 0, V_2 = 0.$ Then, by the equation of (3), we obtain a corrected stiffness matrix given by

and

which implies that $K$ is a symmetric solution of $X^\top KX = \Lambda$.

Example 2. Given matrices

Since $\|AA^+CB^+B-C\| = 9.3907 \times 10^{-12},$ and the eigenvalues of $K_{11}$ are

which means that the conditions of (16) are satisfied. Choose $V, H_1$ and $H_2$ as

Then, by the equation of (17), we get a solution of Eq (1):

with corresponding residual

Furthermore, it can be computed that the eigenvalues of $X+X^\mathrm{H}$ are

which implies that $X$ is a Re-nnd solution of Eq (1).

Example 3. Let the matrices $A$ and $B$ be the same as those in Example 2 and the matrix $C$ be given by

Since $\|AA^+CB^+B-C\| = 1.7944 \times 10^{-11},$ and the eigenvalues of $K_{11}$ are

which means that the conditions of (20) are satisfied. If select $V, K_{12}$ and $H_3$ as

Then, by the equation of (21), we can achieve a Re-pd solution of Eq (1):

and the eigenvalues of $X+X^\mathrm{H}$ are

Furthermore, it can be computed that

7.   Concluding remarks

In this paper, we mainly consider some special solutions of Eq (1). By imposing some constraints on the expression $X = A^+CB^++F_AL_1+L_2E_B,$ we succeed in obtaining a set of necessary and sufficient conditions for the existence of the Hermitian, skew-Hermitian, Re-nonnegative definite, Re-positive definite, Re-nonnegative definite least-rank and Re-positive definite least-rank solutions of Eq (1), respectively. Moreover, we give the explicit expressions for these special solutions, when the consistent conditions are satisfied.

Acknowledgments

The authors wish to give special thanks to the editor and the anonymous reviewers for their helpful comments and suggestions which have improved the presentation of the paper.

Conflict of interest

The authors declare no conflict of interest.