- Research
- Open access
- Published:
A new method for high-order boundary value problems
Boundary Value Problems volume 2021, Article number: 48 (2021)
Abstract
This paper presents a numerical algorithm for solving high-order BVPs. We introduce the construction method of multiscale orthonormal basis in \(W^{m}_{2}[0,1]\). Based on the orthonormal basis, the numerical solution of the boundary value problem is obtained by finding the ε-approximate solution. In addition, the convergence order, stability, and time complexity of the method are discussed theoretically. At last, several numerical experiments show the feasibility of the proposed method.
1 Introduction
High-order BVPs are important mathematical models in the field of electro-magnetics, fluid mechanics, and material science. Many problems in the theory of elastic stability can be handled by BVPs [1]. It is difficult to find the analytic solutions of high-order BVPs because of the complexity of the systems, many numerical algorithms for high-order BVPs have been proposed in recent years. The multistage integration method is an important method to solve the numerical solution of high-order models by reducing the order gradually [2–5]. Ref. [6–8] discuss the existence of solutions to higher-order differential equations. Cao [9] solved a class of high-order fractional ordinary differential equations by the quadratic interpolation function method. The collocation method proposed by [10] and the orthonormal Bernstein polynomials method proposed by Mirzaee [11, 12] can solve high-order linear complex differential equations effectively. Mirzaee et al. [13–22] proposed a variety of numerical algorithms for solving high-order integro-differential equations. Many scholars have also proposed many methods in the field of numerical solution of high-order partial differential equations [23–25]. Reproducing kernel space is an important Banach space which has been used in the field of numerical analysis. The reproducing kernel methods are used in the numerical solutions of high-order models, singular BVPs, and interface problems [26–31].
In this paper, we construct a set of multiscale orthonormal bases based on the idea of wavelet in the reproducing kernel space. This set of bases is orthonormal, which can improve the computational efficiency. For the numerical solution of differential equations, many literature works, such as [32, 33], use the idea of ε-approximate solution. ε-approximate solution provided the stability of the algorithm, good order of convergence in the calculation method. In this article, we construct the multiscale method for the following high-order boundary value problems (BVPs):
where \(B_{i}\) (\(i=1,2,\ldots , m\)) are bounded linear functionals on \(W_{2}^{m}[0,1]\), \(p_{i}(x)\) (\(i=1,2,\ldots , m\)) have certain smoothness.
The paper is organized as follows: In Sect. 2, we construct a set of multiscale orthonormal bases in the reproducing kernel space \(W^{m}_{2}[0,1]\). In Sect. 3, we introduce a method to obtain the numerical solution of BVPs by finding an ε-approximate solution, and verify the existence of the ε-approximate solution. In Sect. 4, the convergence, stability, and complexity of this method are discussed. In Sect. 5, we report the numerical result obtained by the present method and compare this method with other previous methods.
2 Multiscale orthonormal basis
In this section, the reproducing kernel space is defined and a set of multiscale orthonormal bases is constructed. This knowledge is very useful in the following article.
Definition 2.1
The reproducing kernel space \(W^{m}_{2}[0,1]=\{u| u^{(m-1)}\in C[0,1], u^{(m)} \in L^{2}[0,1] \}\), and the inner product of \(W^{m}_{2}\) is
Definition 2.2
The reproducing kernel space
Clearly, \(W^{m}_{2,0}[0,1]\) is the closed subspace of \(W^{m}_{2}[0,1]\). In Ref. [32], we set up a multiscale orthonormal basis in \(W^{1}_{2,0}\):
where \(i=1,2,\ldots \) , \(k=0,1,2,\ldots ,2^{i-1}-1\). The graph of \(\phi _{i,k}\) is shown in Fig. 1.
In order to solve Eq. (1.1), this paper constructs a set of orthonormal bases in \(W_{2}^{m}[0,1]\) by \(\{\phi _{i,k}\}_{i=1,k=0}^{\infty ,2^{i-1}-1}\). First, let us construct the basis functions in \(W_{2,0}^{m}[0,1]\). Note
Theorem 2.1
\(\{J^{m-1}_{0}\phi _{1,0}(x),J^{m-1}_{0}\phi _{2,0}(x),J^{m-1}_{0} \phi _{2,1}(x),\ldots ,J^{m-1}_{0}\phi _{i,k}(x),\ldots \}\) is the multiscale orthonormal basis of \(W^{m}_{2,0}[0,1]\).
Proof
We just prove the orthogonality and completeness.
First, orthonormality. Obviously,
orthogonality is true.
Second, completeness. If \(\langle u, J^{m-1}_{0}\phi _{i,k} \rangle _{W^{m}_{2,0}}=0\), then \(u\equiv 0\), which means the basis is complete. In fact,
According to Def. 2.2, \(u\equiv 0\). □
Next, we construct the orthonormal basis in \(W^{m}_{2}[0,1]\). There are \(m+1\) more conditions in space \(W^{m}_{2,0}[0,1]\) than in space \(W^{m}_{2}[0,1]\). If the basis in \(W^{m}_{2}[0,1]\) is constructed from the basis in \(W^{m}_{2,0}[0,1]\), we need to look for \(m+1\) functions \(g_{k}(x)\in W^{m}_{2}[0,1]\), \(k=0,1,2,\ldots ,m\), such that
It is clear that \(g_{1}(x)=1\), \(g_{2}(x)=x\) in \(W^{m}_{2}[0,1]\) satisfy Eq. (2.2)–Eq. (2.4). Let \(g_{k}(x)=ax^{k} \in W^{m}_{2}[0,1] \), (\(k=2,\ldots ,m\)). By the definition of inner product and Eq. (2.2)–Eq. (2.4), we can obtain \(a=\frac{1}{k!}\).
Theorem 2.2
is the multiscale orthonormal basis of \(W^{m}_{2}[0,1]\).
Proof
According to Th. 2.1 and Eq. (2.2)–Eq. (2.4), it is clear that
So \(\{\rho _{j}(x)\}_{j=1}^{\infty }\) is orthogonal.
Next, we just need to prove completeness. That is, if \(\langle u, \rho _{j} \rangle _{W^{m}_{2}}=0\), then \(u\equiv 0\). In fact,
3 ε-approximate solution of high-order BVPs
In this section, we give the ε-approximation of Eq. (1.1) and get the numerical solution of BVPs by finding the ε-approximate solution of Eq. (1.1).
Put
where \(L:W^{m}_{2}[0,1]\to L^{2}[0,1]\).
Theorem 3.1
\(L:W^{m}_{2}[0,1]\to L^{2}[0,1]\) is a bounded linear operator.
Proof
Because \(W^{m}_{2}[0,1]\) is a reproducing kernel space,
By Eq. (3.1), there exist positive constants \(M_{k}\) such that
Therefore
From Eq. (3.2) and Eq. (3.3), it follows that
where M is a positive constant. □
Then Eq. (1.1) is equivalent to the following equation:
Zhang [32] proposed the ε-approximate theory of second-order differential equations, now we define the ε-approximate solution of Eq. (3.4) based on the idea.
Definition 3.1
\(\forall \varepsilon >0\), \(\exists N>0\), when \(n>N\), if \(\|Lu_{n}-f\|_{L^{2}}^{2}+\sum_{i=1}^{m}(B_{i}u_{n}-\alpha _{i})^{2}< \varepsilon ^{2}\), \(u_{n}\) is called ε-approximate solution of Eq. (3.4).
Lemma 3.1
\(\forall \varepsilon >0\), \(\exists N>0\), when \(n>N\),
is the ε-approximate solution of Eq. (3.4), where \(c_{k}^{*}\) satisfies
Put J is a quadratic form about \(c=(c_{1},\ldots ,c_{n})\),
\(c_{k}^{*}\) is the minimum point of J. In fact, in order to find the minimum value of J, that is,
Because of
then
Let
Then Eq. (3.7) changes to
According to [32], the unique solution of Eq. (3.8) is the minimum point of J.
4 Theoretical analysis
In this section, the properties of the algorithm, such as uniform convergence, stability, and complexity, are introduced.
4.1 Convergence analysis
Theorem 4.1
Assume that u is the exact solution of Eq. (1.1), \(u_{n}\) is the ε-approximation of Eq. (1.1). If \(u^{(m+1)}\) is bounded on [0,1], then \(|u-u_{n}|^{2}\leq 2^{-2n}C\), where C is a constant.
Proof
Assume that
satisfies Eq. (3.4), where \(c_{j}=\langle u,\frac{x^{j}}{j!} \rangle _{W^{m}_{2}}\), \(c_{i,k}=\langle u,J^{2}_{0}\phi _{ik} \rangle _{W^{m}_{2}}\). Clearly, \({\lim_{n \to +\infty }}v_{n}=u\).
Due to Lemma 3.1, it can get
That is,
We can obtain \((c_{i,k})^{2} \leq (\frac{1}{2} )^{3i}C_{1}\). In fact,
According to Holder’s inequality,
Because \(u^{(m+1)}\) is bound, so \(|\int _{\frac{k}{2^{i-1}}}^{\frac{k+1}{2^{i-1}}} (u^{(m+1)})^{2}\,dx| \leq \frac{1}{2^{i-1}}M\). Then
So then
where M is a constant.
where C is a constant, \(K(x,y)\) is the reproducing kernel of \(W_{2}^{m}\). □
From Theorem 4.1, \(u_{n}\) uniformly converges to u.
4.2 Stability analysis
It is well known that if A is a reversible symmetrical matrix, then the condition number of A is
where \(\lambda _{1}\) and \(\lambda _{n}\) are the maximum and minimum eigenvalues of A respectively.
Obviously, A of Eq. (3.8) is an invertible symmetric matrix. Therefore, in order to prove the stability of the algorithm, we can first prove the boundedness of the eigenvalues.
Lemma 4.1
Suppose \(\lambda \mathbf{x=Ax}\), \(\|\mathbf{x}\|=1\), where \(\mathbf{x}=(x_{1},\ldots ,x_{n})^{T}\) is the related eigenvector of λ, then
Proof
According to \(\lambda \mathbf{x=Ax}\),
Multiply both sides of (4.2) by \(x_{i}\) (\(i = 1,2,\ldots ,n\)) and add up to get
So
In addition, \(\lambda \ge \|\sum_{i=1}^{n}x_{i}L\rho _{i}\|_{L^{2}}\).
Let \(u=\sum_{i=1}^{n}x_{i}\rho _{i}\), then \(\lambda \ge \|Lu\|_{L^{2}}^{2}\). Without loss of generality, put \(\|u\|_{W^{m}_{2,0}}=1\). According to the inverse operator theorem [21], \(\|Lu\|_{L^{2}}\ge C^{2}\|u\|^{2}_{W^{m}_{2}}\).
So \(\lambda \ge \|Lu\|_{L^{2}} \ge C^{2}\|u\|^{2}_{W^{m}_{2}}=C^{2}\).
To sum up
□
From Lemma 4.1, we get
That is, the presented method is stable.
4.3 Complexity analysis
Complexity analysis includes time complexity and space complexity. But ultimately, it is the time efficiency of the algorithm that matters. As long as the algorithm does not take up storage space that is unacceptable to the computer. So this part analyzes the time complexity.
Theorem 4.2
The time complexity of the algorithm is \(O(n^{3})\).
Proof
There are four steps to calculate ε-the approximate solution \(u_{n}(x)\) of (3.1).
First, the calculation of matrix An in Eq. (3.8). The matrix An is
Set the number of multiplication required to compute \(\langle L\rho _{k},L\rho _{j} \rangle \) and \(B_{l}\rho _{k} B_{l}\rho _{j}\) as \(C_{1}\), \(C_{2}\) respectively, \(C_{1}\), \(C_{2}\) are constant. So each term of An is evaluated \(C_{1}+mC_{2}\) times. Since An is a symmetric matrix, we only need to consider the calculation of the main diagonal and above elements. The first row of An is evaluated n times, the second row \(n-1\) times, and so on, so that the total number of multiplications required in calculation of An is
Second, the calculation of vector bn in Eq. (3.8). The vector bn is
Let the number of multiplication of \(\langle L\rho _{k},f \rangle \) and \(B_{l}\rho _{j}\alpha _{l}\) be \(C_{3}\), \(C_{4}\), \(C_{3}\), \(C_{4}\) are constant. So each term of bn is evaluated \(C_{3}+mC_{4}\) times. So the total number of multiplications required in calculation of bn is
Third, solve Eq. (3.8). We solve the system by Gaussian elimination. From the mathematical knowledge, Gaussian elimination requires operations
Forth, calculation \(u_{n}\). When calculating the \(u_{n}\), the total number of multiplications required is n.
To sum up, the total number of multiplication is
□
5 Numerical experiments
In this section, we give several numerical experiments to verify the effectiveness of the proposed algorithm. We denote by \(u_{n}(x)\) the approximation to the exact solution \(u(x)\) obtained by the numerical schemes in the present work, and we measure the errors in the following sense:
where n is the number of bases. C.R. represents the convergence order. All numerical experiments are computed by Mathematica 9.0.
Example 5.1
Ref. [34] mentioned that in order to get the shear deformation of sandwich beams, consider the following third-order BVP:
where the physical constants are \(k=5\) and \(r=1\). The function \(u(x)\) shows the shear deformation of sandwich beams. The analytic solution of this problem is
The numerical results are given in Table 1.
Example 5.2
Our second example is for fourth-order BVP [31]
Table 2 and Fig. 2 list the error comparison between the multilevel augmentation method [31] and our method for this BVP.
Example 5.3
Our third example is a fifth-order BVP
The exact solution of this example is \(u(x)=e^{x}-x-1\). Figure 3 shows the error of Example 5.3.
6 Conclusion
In summary, this study used a set of multiscale orthonormal bases to find the ε-approximate solutions of higher-order BVPs. This paper not only demonstrates the convergence and stability in theory, but also demonstrates the feasibility of the method through numerical experiments. Through theoretical analysis and numerical experiments, this method can be extended to solve general linear models, such as linear integral equations, differential equations, and fractional differential equations.
Availability of data and materials
The authors declare that all data and materials in this paper are available and veritable.
References
Timoshenko, S.P., Gere, J.M.: Theory of Elastic Stability. McGrawHill-Kogakusha, Tokyo (1961)
Jackiewicz, Z., Mittelmann, H.: Construction of IMEX DIMSIMs of high order and stage order. Appl. Numer. Math. 121, 234–248 (2017)
Ling, X., Gao, N., Hu, A.: Dynamic analysis of a planetary gear system with multiple nonlinear parameters. J. Comput. Appl. Math. 327, 325–340 (2017)
Moradi, A., Sharififi, M., Abdi, A.: Transformed implicit-explicit second derivative diagonally implicit multi stage integration methods with strong stability preserving explicit part. Appl. Numer. Math. 156, 14–31 (2020)
Lei, X., Li, J.: Transversal effects of high order numerical schemes for compressible fluid flows. Appl. Math. Mech. 3, 343–354 (2019)
Lü, X., Cui, M.: Existence and numerical method for nonlinear third-order boundary value problem in the reproducing kernel space. Bound. Value Probl. 2010, 459754 (2010)
Ji, Y., Guo, Y., Yao, Y.: Positive solutions for higher order differential equations with integral boundary conditions. Bound. Value Probl. 2015, 214 (2015)
Xu, L., Chen, H.: Existence and multiplicity of solutions for fourth-order elliptic equations of Kirchhoff type via genus theory. Bound. Value Probl. 2014, 212 (2014)
Cao, J., Wang, Z., Xu, C.: A high-order scheme for fractional ordinary differential equations with the Caputo–Fabrizio derivative. Commun. Appl. Math. Comput. 2(2), 179–199 (2020)
Toutounian, F., Tohidi, E., Shateyi, S.: A collocation method based on the Bernoulli operational matrix for solving high-order linear complex differential equations in a rectangular domain. Abstr. Appl. Anal. 4, 215–222 (2013)
Mirzaee, F., Samadyar, N., Alipour, S.: Numerical solution of high order linear complex differential equations via complex operational matrix method. SeMA J. 76, 1–13 (2019)
Mirzaee, F., et al.: Parameters estimation of HIV infection model of CD4+ T-cells by applying orthonormal Bernstein collocation method. Int. J. Biomath. 11(2), 1850020 (2018)
Raslan, K., et al.: Numerical solution of high-order linear integro differential equations with variable coefficients using two proposed schemes for rational Chebyshev functions. New Trends Math. Sci. 4(3), 22–35 (2016)
Yzba, S., Yildirim, A., et al.: A collocation approach for solving high-order linear Fredholm–Volterra integro-differential equations. Math. Comput. Model. 55, 547–563 (2012)
Mirzaee, F., Rafei, Z.: The block by block method for the numerical solution of the nonlinear two-dimensional Volterra integral equations. J. King Saud Univ., Sci. 23, 191–195 (2011)
Mirzaee, F., Hoseini, S.F.: Solving systems of linear Fredholm integro-differential equations with Fibonacci polynomials. Ain Shams Eng. J. 5, 271–283 (2014)
Mirzaee, F., Hoseini, S.F.: A new collocation approach for solving systems of high-order linear Volterra integro-differential equations with variable coefficients. Appl. Math. Comput. 311, 272–282 (2017)
Mirzaee, F., Bimesl, S.: An efficient numerical approach for solving systems of high-order linear Volterra integral equations. Sci. Iran. 21(6), 2250–2263 (2014)
Mirzaee, F., Bimesl, S., Tohidi, E.: A numerical framework for solving high-order pantograph delay Volterra integro-differential equations. Kuwait J. Sci. 43(1), 69–83 (2016)
Mirzaee, F., Bimesl, S.: Numerical solutions of systems of high-order Fredholm integro-differential equations using Euler polynomials. Appl. Math. Model. 39, 6767–6779 (2015)
Mirzaee, F., Bimesl, S.: Application of Euler matrix method for solving linear and a class of nonlinear Fredholm integro-differential equations. Mediterr. J. Math. 11, 999–1018 (2014)
Samadyar, N., Orthonormal, M.F.: Bernoulli polynomials collocation approach for solving stochastic Itô–Volterra integral equations of Abel type. Int. J. Numer. Model. 2019, e2688 (2019)
Ren, J., Shi, D., Vong, S.: High accuracy error estimates of a Galerkin finite element method for nonlinear time fractional diffusion equation. Numer. Methods Partial Differ. Equ. 36, 284–301 (2020)
Li, Z., Liang, Z., Yan, Y.: High-order numerical methods for solving time fractional partial differential equations. J. Sci. Comput. 71, 785–803 (2017)
Mirzaee, F., Bimesl, S.: A new approach to numerical solution of second-order linear hyperbolic partial differential equations arising from physics and engineering. Results Phys. 3, 241–247 (2013)
Mei, L.: A novel method for nonlinear impulsive differential equations in broken reproducing kernel space. Acta Math. Sci. 40, 723C733 (2020)
Zhao, Z., Lin, Y., Niu, J.: Convergence order of the reproducing kernel method for solving boundary value problems. Math. Model. Anal. 21(4), 466–477 (2016)
Tirmizi, I.A., Twizell, E.H., Islam, S.U.: A numerical method for third-order non-linear boundary-value problems in engineering. Int. J. Comput. Math. 82, 103–109 (2005)
Li, X., Wu, B.: Reproducing kernel method for singular multipoint boundary value problems. Math. Sci. 6(1), 1–5 (2016)
Li, X., Wu, B.: A new kernel functions based approach for solving 1-D interface problems. Appl. Math. Comput. 380, 125276 (2020)
Chen, Z., Wu, B., Xu, Y.: Multilevel augmentation methods for differential equations. Adv. Comput. Math. 24, 213–238 (2006)
Zhang, Y., Sun, H., Jia, Y., Lin, Y.: An algorithm of the boundary value problem based on multiscale orthogonal compact base. Appl. Math. Lett. 101, 106044 (2020)
Zheng, Y., Lin, Y., Shen, Y.: A new multiscale algorithm for solving second order boundary value problems. Appl. Numer. Math. 156, 528–541 (2020)
Haque, M., Baluch, M.H., Mohsen, M.F.N.: Solution of multiple point, nonlinear boundary value problems by method of weighted residuals. Int. J. Comput. Math. 19, 69–84 (1986)
Acknowledgements
Not applicable.
Authors’ information
Zhuhai Campus, Beijing Institute of Technology, Zhuhai, Guangdong, 519085, China
Funding
This work has been supported by the Basic and Applied Basic Research Project of Zhuhai City (ZH22017003200026PWC) and the Characteristic Innovative Scientific Research Project of Guangdong Province (2019KTSCX217).
Author information
Authors and Affiliations
Contributions
This is to declare that all authors have contributed equally and significantly to the contents of the paper. All authors have read and agreed to the published version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhang, Y., Mei, L. & Lin, Y. A new method for high-order boundary value problems. Bound Value Probl 2021, 48 (2021). https://doi.org/10.1186/s13661-021-01527-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-021-01527-4