Dynamic analysis and optimal control of rumor propagation models considering different education levels and hesitation mechanisms

1.   Introduction

Partial differential equations are closely related to many physical problems in real life. For example, the ninth-order linear or non-linear boundary value problems are related to the laminar viscous flow in a semi-porous channel or the hydro-magnetic stability, and the telegraph equations are related to the vibrations within objects or the propagation of waves. There are also many different types of partial differential equations in engineering and other applied sciences. Zhang et al. ^[1] discussed cubic spline solutions of ninth-order linear and non-linear boundary value problems using a cubic B-spline. Shah et al. ^[2] proposed a new and efficient operational matrix method for solving time-fractional telegraph equations with Dirichlet boundary conditions. Nisar et al. ^[3] proposed a hybrid mesh free framework based on Pad$\acute{e}$ approximation in order to solve the numerical solutions of nonlinear partial differential equations. There are also many different types of partial differential equations in chemistry, engineering, and other applied sciences. There have been many successful conclusions about these partial differential equations.

Complex partial differential equations of analytic functions also have a wide range of applications. Bi-analytic functions, the generalizations of analytic functions, have important applications in elasticity. In 1961, Sander ^[4] studied the properties of pairs of functions $\{u(x, y), v(x, y)\}$ with binary real variables $(x, y)$, which satisfy the system of partial differential equations:

for the real constant $k (k\neq -1)$, where $\theta(x, y)$ and $\omega(x, y)$ are continuously differentiable functions of $x$ and $y$. Sander provided the definition of bi-analytic functions of type $k$, which are of great significance for studying some physical problems for $k > 0$, and extended some properties of analytic functions to bi-analytic functions.

In 1965, Lin and Wu ^[5] introduced the function class that is more extensive than Sander's function class, i.e., bi-analytic functions of the type $(\lambda, k)$ which are defined by the system of equations:

where $\theta(x, y)$ and $\omega(x, y)$ are continuously differentiable functions of $x$ and $y$, and $\lambda, k$ are real constants with $\lambda\neq0, 1, k^2$, and $0 < k < 1$. The complex form of the system is

in which $\varphi(z) = k\vartheta-i\lambda\omega$ is analytic and is called the associate function of $f(z) = u+iv$. In ^[5], the general expression and the properties of bi-analytic functions of the type $(\lambda, k)$ were researched in detail.

Hua et al. ^[6] introduced a mechanical interpretation for bi-analytic functions and promoted the corresponding function theory. Thereafter, bi-analytic functions aroused widespread attention from many scholars ^[7,8,9]. In 1994, Kumar ^[10] discussed a broader class of functions, i.e., bi-polyanalytic functions, and investigated several Riemann-Hilbert problems for systems of $n$-order partial differential equations applying polyanalytic functions ^[11] and bi-polyanalytic functions on the unit disk. In 2005, Kumar and Prakash ^[12] investigated Dirichlet problems for the Poisson equation and some boundary value problems for bi-polyanalytic functions on the unit disk. They obtained the explicit representations of the solutions and the corresponding solvable conditions. In 2006, Begehr and Kumar ^[13] discussed some complex partial differential equations of higher order. Some boundary value problems for bi-polyanalytic functions were solved on different conditions on the unit disk.

In recent years, some other boundary value problems for bi-analytic functions were solved ^{[14,15,16,17]}. With the gradual improvements of the theory for bi-analytic functions and polyanalytic functions ^{[18,19,20,21]} in the complex plane, some scholars attempted to generalize the relevant achievements to spaces of several complex variables ^[22,23].

In this paper, based on the work of the former researchers, we study a class of Schwarz problems for polyanalytic functions on the bicylinder. Then, from the perspective of series and applying the particular solution to the Schwarz problem for polyanalytic functions, we discuss a Dirichlet problem for bi-polyanalytic functions on the bicylinder.

In the following, let the bicylinder $D^2 = D_1\times D_2 = \{(z_1, z_2):|z_1| < 1, \; |z_2| < 1\},$ and let $\partial_0D^2$ denote the characteristic boundary of $D^2$. Let $C(G)$ represent the set of continuous functions within $G$.

2.   The Schwarz problem

To get the main results, we need to discuss the following Schwarz problem.

Theorem 2.1. Let $g_{\mu\nu}\in C(\partial_0D^2;\mathbb{R})$ for $1\leq \mu, \nu\leq m-1$ ($m\geq2$), and let

where

and

in which

and

Then, $\widetilde{\phi}(z)$ satisfies

Proof: 1) Let

in which

is analytic on $D^2$, and

In the following, we will show that $\widetilde{\phi}(z)$ satisfies

By (2.3), we get that

where

Therefore, for $1\leq k_1, k_2\leq m-1$,

in which

Let

For $1\leq k_1, k_2\leq m-1$, (2.6) follows that

Applying the properties of the Cauchy kernels on $D^2$ and $D$, for $z\in \partial_0D^2$, (2.7) leads to

which means $\Re\partial^{\mu}_{\bar{z}_1}\partial^{\nu}_{\bar{z}_2}\widetilde{\phi}(z) = g_{\mu\nu}(z)$ for $z\in \partial_0D^2$.

In addition, by (2.7), we get that

Similarly, we have that

Therefore,

2) In the expression of $u_{\widetilde{v_1}\widetilde{v_2}}(z)$ determined by (2.4),

Moreover, we have that

in which

and

where

Similarly, we get that

in which

Plugging (2.9)–(2.11) into (2.8) gives that

In addition, as the result of

in which

we have that

which follows that

Similarly, we have that

and

in which

Therefore,

as the result of

On the other hand, due to the analyticity of the function $u_{\widetilde{v_1}\widetilde{v_2}}(z)$, it can be expressed as

Plugging (2.12), (2.14), (2.16), and (2.17) into (2.4), and considering the Eq (2.18), we get (2.2). Moreover, (2.18) and (2.3) lead to (2.1). Therefore, from the result in the first part (1), it can be concluded that $\widetilde{\phi}(z)$ determined by (2.1) satisfies

3.   The Dirichlet problem

Theorem 3.1. Let $\varphi\in C(\partial_0D^2;\mathbb{C})$, $\lambda\in\mathbb{R}\setminus\{-1, 0, 1\}$, and let $g_{\mu\nu}\in C(\partial_0D^2;\mathbb{R})$ for $1\leq \mu, \nu\leq m-1$ ($m\geq2$). Then, the problem

with the conditions

is solvable and the solution is

where $u_{m_1, m_2}^{m-v_1, m-v_2}$ is determined by (2.2) (in which $\widetilde{v_1}$ and $\widetilde{v_2}$ are replaced by $v_1$ and $v_2$, respectively), and $u_{m_1, m_2}^{0, 0}$, $b_{m_1, m_2}$ are determined by the following:

(ⅰ) for $1\leq m_1\leq m-1$ and $1\leq m_2\leq m-1$,

(ⅱ) for $m_1\geq 1$ and $m_2\geq 1$,

(ⅲ) with $u^{0, 0}_{1, 1}$ being determined by (3.1),

(ⅳ) for $\{\substack{m_1\geq m\\m_2\geq m}$ or $\{\substack{1\leq m_1\leq m-1\\m_2\geq m}$ or $\{\substack{m_1\geq m\\1\leq m_2\leq m-1}$,

(ⅴ) for $m_1, \; m_2\geq 1$,

in which $u_{(m_1+1), (m_2+1)}^{0, 0}$, $u_{(m_1+1), 1}^{0, 0}$, and $u_{1, (m_2+1)}^{0, 0}$ are determined by (3.6).

Proof: 1) By Theorem 2.1,

where $u_{m_1, m_2}^{m-v_1, m-v_2}$ is determined by (2.2) (in which $\widetilde{v_1}$ and $\widetilde{v_2}$ are replaced by $v_1$ and $v_2$, respectively), and $u_0(z)$ is analytic on $D^2$. Thus, $u_0(z)$ can be represented as

in which $u_{m_1, m_2}^{0, 0}$ is to be determined.

Let

where

Thus, $\partial_{\bar{z}_1}\partial_{\bar{z}_2}\phi_1(z) = \phi(z)$. Let

and let

Then, we get that $\partial_{z_1}\partial_{z_2}\widetilde{u_0}(z) = u_0(z)$ and $\partial_{z_1}\partial_{z_2}\widetilde{u}_{v_1, v_2}(z) = u_{v_1, v_2}(z).$ Let

which follows that

Therefore,

which means that

is a special solution to

So, the solution of the problem is

where $\psi(z) = \sum_{m_1, m_2 = 0}^{+\infty}b_{m_1, m_2}z_1^{m_1}z_2^{m_2}$ is analytic on $D^2$, and $u_{m_1, m_2}^{0, 0}$ and $b_{m_1, m_2}$ are to be determined.

2) In this part, we seek the expressions of $u_{m_1, m_2}^{0, 0}$ and $b_{m_1, m_2}$.

For $z\in\partial_0D^2$, let $z_1 = e^{it_1}$ and $z_2 = e^{it_2}$ ($t_1, t_2\in[0, 2\pi]$). Then, we get that

(ⅰ) In the case of $0\leq r_1, r_2\leq m-2$, multiplying both sides of the Eq (3.11) by $e^{it_1(m-r_1)}e^{it_2(m-r_2)}$, and then integrating with respect to $t_1, t_2\in[0, 2\pi]$ yields that

and the last equation is due to

Therefore,

which leads to (3.1) for $1\leq m_1, \; m_2\leq m-1$.

(ⅱ) Multiplying both sides of the Eq (3.11) by $e^{it_1}e^{it_2}$, and then integrating with respect to $t_1, t_2\in[0, 2\pi]$ yields that

which leads to (3.2).

Additionally, for $r_2\geq 1$, multiplying both sides of the Eq (3.11) by $e^{it_1}e^{i(1-r_2)t_2}$, and then integrating with respect to $t_1, t_2\in[0, 2\pi]$ yields that

and the last equation is due to

Therefore, for $r_2\geq 1$,

For $r_1\geq 1$, multiplying both sides of the equation (3.11) by $e^{i(1-r_1)t_1}e^{it_2}$ and then integrating with respect to $t_1, t_2\in[0, 2\pi]$, similar to (3.12), yields that

(ⅲ) In addition, integrating both sides of the Eq (3.11) with respect to $t_1, t_2\in[0, 2\pi]$ yields that

which follows (3.5), where $u^{0, 0}_{1, 1}$ is determined by (3.1).

(ⅳ) In the case of $r_1, r_2\geq m+1$, multiplying both sides of the Eq (3.11) by $e^{ir_1t_1}e^{ir_2t_2}$, and then integrating with respect to $t_1, t_2\in[0, 2\pi]$ yields that

Therefore,

for $r_1, r_2\geq m$.

Similarly, in the case of $r_1\geq m+1$ and $0\leq r_2\leq m-2$, multiplying both sides of the Eq (3.11) by $e^{ir_1t_1}e^{it_2(m-r_2)}$, and then integrating with respect to $t_1, t_2\in[0, 2\pi]$ yields that $u_{r_1, r_2}^{0, 0}$ ($r_1\geq m, \; 1\leq r_2\leq m-1$) has the same representation as (3.14). In the case of $0\leq r_1\leq m-2$ and $r_2\geq m+1$, multiplying both sides of the Eq (3.11) by $e^{it_1(m-r_1)}e^{ir_2t_2}$, and then integrating with respect to $t_1, t_2\in[0, 2\pi]$ yields that $u_{r_1, r_2}^{0, 0}$ ($1\leq r_1\leq m-1, \; r_2\geq m$) has the same representation as (3.14).

(ⅴ) Similarly, for $r_1, r_2\geq 1$, multiplying both sides of the Eq (3.11) by $e^{-ir_1t_1}e^{-ir_2t_2}$, and then integrating with respect to $t_1, t_2\in[0, 2\pi]$, we can get the expression of $b_{r_1r_2}$ which leads to (3.7). For $r_1\geq 1$, multiplying both sides of the Eq (3.11) by $e^{-ir_1t_1}$, and then integrating with respect to $t_1, t_2\in[0, 2\pi]$ yields the expression of $b_{r_10}$ which leads to (3.8). For $r_2\geq 1$, we can get $b_{0r_2}$ which leads to (3.9).

Remark 3.2. The results obtained in this article extend the existing conclusions about ployanalytic functions and bi-ployanalytic functions. On this basis, we can study other partial differential equation problems. For example, it would be interesting to discuss whether bi-polyanalytic or even ployanalytical function solutions exist for some nonlocal integrable partial differential equations (see, e.g., ^[24]), which need to explore the solutions to the corresponding Riemann-Hilbert problems. Additionally, the non-existence of solutions to Cauchy problems on the real line for first-order nonlocal differential equations (see, e.g., ^[25]) indicates that we can attempt to discuss the generalizations of analytical solutions for partial differential equation problems.

4.   Conclusions

With the help of the series expansion of polyanalytic functions, and applying the properties of Cauchy kernels on the bicylinder and the unit disk, we first discuss a class of Schwarz problems with the conditions concerning the real and imaginary parts of high-order partial differentiation for polyanalytic functions on the bicylinder. On this basis, we investigate a type of boundary value problem for bi-polyanalytic functions with Dirichlet boundary conditions on the bicylinder. From the perspective of series, we obtain the specific representation of the solution to the Dirichlet problem. The method used in this article, with the help of series expansion, is different from the previous methods for solving boundary value problems. It is a very effective method and can be used to solve other types of problems regarding complex partial differential equations of bi-polyanalytic functions in high-dimensional complex spaces. The conclusions of this article also lay a necessary foundation for further research on polyanalytic and bi-polyanalytic functions.

Author contributions

Yanyan Cui: conceptualization, project administration, writing original, writing-review and editing; Chaojun Wang: investigation, writing original, writing-review and editing. All authors have read and approved the final version of the manuscript for publication.

Acknowledgments

The authors are grateful to the anonymous referees for their valuable comments and suggestions which improved the quality of this article. This work was supported by the NSF of China (No. 11601543), the NSF of Henan Province (Nos. 222300420397 and 242300421394), and the Science and Technology Research Projects of Henan Provincial Education Department (No. 19B110016).

Conflict of interest

The authors declare no conflict of interest.