Scaling, Self-Similarity, and Systems of Fractional OrderView this Special Issue
Hybrid Bernstein Block-Pulse Functions Method for Second Kind Integral Equations with Convergence Analysis
We introduce a new combination of Bernstein polynomials (BPs) and Block-Pulse functions (BPFs) on the interval [0, 1]. These functions are suitable for finding an approximate solution of the second kind integral equation. We call this method Hybrid Bernstein Block-Pulse Functions Method (HBBPFM). This method is very simple such that an integral equation is reduced to a system of linear equations. On the other hand, convergence analysis for this method is discussed. The method is computationally very simple and attractive so that numerical examples illustrate the efficiency and accuracy of this method.
In recent years, many different basic functions have been used for solving integral equations, such as Block-Pulse functions [1, 2], Triangular functions , Haar functions , Hybrid Legendre and Block-Pulse functions , Hybrid Chebyshev and Block-Pulse functions [6, 7], Hybrid Taylor and Block-Pulse functions , and Hybrid Fourier and Block-Pulse functions .
Bernstein polynomials have been applied in various fields of mathematics. For example, some researchers applied the Bernstein polynomials for solving high order differential equations , some classes of integral equations , partial differential equations, and optimal control problems . Also, we introduced new operational matrices of fractional derivative and integral operators by Bernstein polynomials and then used them for solving fractional differential equations [15–17], system of fractional differential equations , and fractional optimal control problems [19, 20].
In this work, we combine the Bernstein polynomials (BPs) and Block-Pulse functions (BPFs) on the interval . Then, we use these bases for finding an approximate solution of the second kind integral equation. We call this method Hybrid Bernstein Block-Pulse Functions Method (HBBPFM). In this method the integral equation is reduced to a system of linear equations. Also, we discuss the convergence analysis for this method. Furthermore, we compare the accuracy of obtained results of BPFs, BPs, and HBBPFM by some examples.
The rest of this paper is as follows. In Section 2, HBBPFs are introduced; therefore we approximate functions by using HBBPFs and also we discuss best approximation and convergence analysis in Section 3. Then we apply HBBPF method to find an approximate solution for the second kind integral equations and we survey error analysis for proposed method in Section 4. Also, we apply the proposed method on some examples. We observe that the accuracy and efficiency of this method are more than the near methods. Finally, Section 6 concludes our work in this paper.
2. Hybrid of Bernstein and Block-Pulse Functions
In this section, we recall some definitions and properties of Bernstein polynomials and Block-Pulse functions.
Lemma 1 (see ). The Bernstein polynomials (BPs) of th-degree are defined on the interval as follows: where Then in Hilbert space is a complete basis. Therefore, any polynomial of degree can be expanded in terms of linear combination of .
Lemma 2. Let a set of Block-Pulse functions (BPFs) , be on the interval such that. Then, the following properties for these functions satisfy the following:(i) disjointness,(ii) orthogonality,(iii) completeness.
Proof. The disjointness property can be clearly obtained from the definition of Block-Pulse functions as follows:
The other property is orthogonality. It is clear that where and is the Kroneker delta.
The third property is completeness. For every , when approaches the infinity, Parseval's identity holds: where .
Definition 3 (Hybrid Bernstein Block-Pulse Functions (HBBPFs)). , , , have three arguments; and are the order of BPFs and BPs, respectively, and is the normalized time. HBBPFs are defined on the interval as follows:
In the next section, we deal with the problem of approximation of these functions.
3. Approximation of Functions by Using HBBPFs and Convergence Analysis
Theorem 4. Suppose that the function is times continuously differentiable, and . Then is the best approximation out of with the following inner product: where and . Also, one can obtain the following inequality: where .
Proof. We prove that is the best approximation for out of . We can prove that is a convex subset of a real inner product space (see ). Therefore, for any , is its best approximation in if and only if it satisfies where the inner product is defined by . Then for any , its best approximation is unique. Also, we know that is a convex and closed finite-dimensional subset of an inner product space . Then for any , there is a unique element such that . Therefore, there exist the unique coefficients , such that On the other hand, we can consider that is a basis for polynomials space of degree . Therefore we define . Hence, from Taylor expansion we have where . Since is the best approximation out of , and we assume that , therefore, we have Then by taking square roots, the proof is complete.
The previous theorem shows that the error vanishes as .
Corollary 5. One can write , such that one defines that is a matrix and is said dual matrix of , and one can obtain
Proof. We know therefore, the proof is complete.
Corollary 6. A function may be expanded as follows: If the infinite series in (16) is truncated, then we have where Therefore we can get Then where where by using (7), is defined as follows: We can also approximate the function as follows: where is an matrix that we can obtain as follows:
Theorem 7. Let the function be times continuously differentiable; then we have where .
Proof. By using Theorem 4 we get where and . Therefore by taking square roots, the proof is complete.
The above theorem shows that the approximation error vanishes as .
4. HBBPFs for the Second Kind Integral Equations and Error Analysis
In this section, we are dealing with the following Fredholm equations of the second kind: where , , and is an unknown function.
Let us approximate , , and by (18) and (25) as follows: By substituting (29) in (28) we obtain Therefore we have the following linear system: that by solving this linear system we can obtain the vector .
Theorem 8. Suppose that is exact solution of (28) and is approximate solution by HBBPFs for and is perturbation function that depends only on . Let . Then as .
Proof. Suppose is the error function of approximate solution to the exact solution . Therefore we get
By taking absolute value and using Holder inequality we get
Now, by taking norm we obtain
Finally, from Theorem 7 we can write
Therefore, we can show that as .
5. Numerical Examples
In this section we discuss the implementation of the new method and investigate its accuracy by applying it to different examples. In the following examples, we suppose that , , and are approximate solutions by BPFs, BPs, and HBBPFM for the exact solution , respectively.
Example 1. Consider the following integral equation: We know that the exact solution is . The obtained results of BPFs, BPs, and HBBPFs are reported in Tables 1 and 2 and are plotted in Figures 1 and 2. We compare the obtained results and observe that HBBPFM is very effective and accuracy of approximate solutions in this method is more than methods of BPFs and BPs.
Example 2. Consider the following integral equation: with exact solution . We obtain the computational by BPFs, BPs, and HBBPFM with , , and , ; then we compare them together. The results are reported in Tables 3 and 4 and are plotted in Figures 3 and 4. Similar to the previous example, we see that the method HBBPFM is very effective and accuracy of solution in this method is more than methods of BPFs and BPs.
In this paper, HBBPFs are used to solve second kind integral equations we call this method with HBBPFM. This method converts second kind integral equations to systems of linear equations whose answers are coefficient of HBBPFs expansion of the solution of second kind integral equations. Also, by using several lemmas and theorems, we have discussed convergence analysis of the proposed method. Numerical examples show the efficiency and accuracy of the method. Moreover we see that accuracy of solutions in HBBPFM is more satisfactory than the methods of BPFs and BPs.
Conflict of Interests
The authors declare that there is no conflict of interests in this paper.
E. Babolian and Z. Masouri, “Direct method to solve Volterra integral equation of the first kind using operational matrix with block-pulse functions,” Journal of Computational and Applied Mathematics, vol. 220, no. 1-2, pp. 51–57, 2008.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
K. Maleknejad and B. Rahimi, “Modification of block pulse functions and their application to solve numerically Volterra integral equation of the first kind,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 6, pp. 2469–2477, 2011.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
H. R. Marzban, H. R. Tabrizidooz, and M. Razzaghi, “A composite collocation method for the nonlinear mixed Volterra-Fredholm-Hammerstein integral equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 3, pp. 1186–1194, 2011.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
X. T. Wang and Y. M. Li, “Numerical solutions of integrodifferential systems by hybrid of general block-pulse functions and the second Chebyshev polynomials,” Applied Mathematics and Computation, vol. 209, no. 2, pp. 266–272, 2009.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
G. Prasada Rao, Piecewise Constant Orthogonal Functions and their Application to Systems and Control, Springer, Berlin, Germany, 1983.View at: MathSciNet
B. M. Mohan and K. B. Datta, Orthogonal Functions in Systems and Control, 1995.View at: MathSciNet
M. Alipour and D. Rostamy, “Solving nonlinear fractional differential equations by Bernstein polynomials operational matrices,” The Journal of Mathematics and Computer Science, vol. 5, no. 3, pp. 185–196, 2012.View at: Google Scholar
D. Rostamy, M. Alipour, H. Jafari, and D. Baleanu, “Solving multi-term orders fractional differential equations by operational matrices of BPs with convergence analysis,” Romanian Reports in Physics, vol. 65, no. 2, pp. 334–349, 2013.View at: Google Scholar
D. Baleanu, M. Alipour, and H. Jafari, “The Bernstein operational matrices for solving the fractional quadratic Riccati differential equations with the Riemann-Liouville derivative,” Abstract and Applied Analysis, vol. 2013, Article ID 461970, 7 pages, 2013.View at: Publisher Site | Google Scholar | MathSciNet
M. Alipour and D. Rostamy, “BPs operational matrices for solving time varying fractional optimal control problems,” The Journal of Mathematics and Computer Science, vol. 6, pp. 292–304, 2013.View at: Google Scholar
K. Atkinson and W. Han, Theoretical Numerical Analysis: A Functional Analysis Framework, Springer, 2nd edition, 2000.