Numerical solution of integral equations system of the second kind by Block–Pulse functions
Introduction
In recent years, many different basic functions have been used to estimate the solution of integral equations. In this paper we use BPF (Block–Pulse functions) as a simple base for solving a system of integral equations. This set of functions was first introduced to electrical engineers by Harmuth. Then several researchers (Gopalsami and Deekshatulu, 1997 [8]; Chen and Tsay, 1977 [9]; Sannuti, 1977 [10]) discussed the Block–Pulse functions and their operational matrix [1], [2].
Section snippets
Definition of BPF
An m-set of BPF is defined as follows:with t ∈ [0, T) and . Fig. 1, Fig. 2, Fig. 3, Fig. 4 illustrate Φ1(t) to Φ4(t).
Now we explain the properties of BPF.
- (i)
Disjointness:We clearly havet ∈ [0, T),i, j = 1, 2, …, m. This property is obtained from definition of BPF.
- (ii)
Orthogonality:We havet ∈ [0, T), i, j = 1, 2, …, m. This property is obtained from the disjointness property.
- (iii)
Completeness:
For every f ∈ L2,{Φ} is complete if
Linear integral equations system
Consider following integral equations system (see [5], [6], [7]):
Our problem is to determine Black–Pulse coefficient of yj(x), j = 1, 2, …, n in the interval x ∈ [α, β) from the known functions fi(x), i = 1, 2, …, n and the kernels kij(x, t) i, j = 1, 2, …, n. Usually we set α = 0 to facilities the use of Black–Pulse functions. In case α ≠ 0 we set where T = mh.
We approximate fi, yj, kij by its BPF as follows:
Numerical examples
Example 1 Consider the integral equations systemwith exact solution y1(t) = t, y2(t) = t2 results with m = 16 and m = 32 shown in Table 1. Example 2 Consider the linear integral equations systemwith exact solution y1(t) = et and y2(t) = e−t. The computational results for m = 16 and m = 32 together with the exact solution are given in Table 2.
Conclusion
Block–Pulse functions which had previously been used in control theory have been developed in this paper to solve a linear Fredholm integral equations system of the second kind by means of references [6], [7]. As examples indicate, in order to increase the accuracy of the numerical results, it is necessary to increase m. Other orthogonal functions can also be used for solving integral equations system of the second kind.
References (10)
- et al.
Solving second kind integral equations by Galerkin methods with hybrid Legendre and Block–Pulse functions
Applied Mathematics and Computation
(2003) - et al.
Piecewise Constant Orthogonal Functions and their Application to System and Control
(1983) - B.M. Mohan, K.B. Datta, Orthogonal Function in Systems and Control,...
- Chin-Fan Chen, Chi, Hung Itsia, Design of piecewise constant gains for optimal control via Walsh functions, IEEE...
- M. Razzaghi, J. Nazarzadeh, Walsh functions, in: Wiley Encyclopedia of Electrical and Electronics Engineering, vol. 23,...
Cited by (86)
Bernoulli wavelet method for numerical solution of linear system of Fredholm integral equation of the second kind
2023, Alexandria Engineering JournalApproximation methods for system of linear Fredholm integral equations of second kind
2021, Applied Mathematics and ComputationA new hybrid orthonormal Bernstein and improved block-pulse functions method for solving mathematical physics and engineering problems
2020, Alexandria Engineering JournalNumerical solution of fractional Mathieu equations by using block-pulse wavelets
2019, Journal of Ocean Engineering and ScienceFractional-order Bernoulli functions and their applications in solving fractional Fredholem–Volterra integro-differential equations
2017, Applied Numerical Mathematics