# Numerical Algorithm Based on Haar-Sinc Collocation Method for Solving the Hyperbolic PDEs

Authors: A. Pirkhedri; H. H. S. Javadi; H. R. Navidi
Journal: The Scientific World Journal
 (2014)

---

## Abstract

The present study investigates the Haar-Sinc collocation method for the solution of the hyperbolic partial telegraph equations. The advantages of this technique are that not only is the convergence rate of Sinc approximation exponential but the computational speed also is high due to the use of the Haar operational matrices. This technique is used to convert the problem to the solution of linear algebraic equations via expanding the required approximation based on the elements of Sinc functions in space and Haar functions in time with unknown coefficients. To analyze the efficiency, precision, and performance of the proposed method, we presented four examples through which our claim was confirmed.


---

## Body

## 1. Introduction

In recent years, the field of hyperbolic partial differential equations has attracted the attention of scientists in several areas and has also been used to solve many actual problems modeled in mathematical physics, such as the vibrations of structures (e.g., buildings, beams, and machines), fluid mechanics, and atomic physics [1].The telegraph equation has typically been used for transmission and propagation of electrical signals [2], wave propagation model [3], random walk theory [4], and so forth. Let us consider an infinitesimal piece of a telegraph wire which consists of resistor of resistance R, a coil of inductance L, a resistor of conductance G, and a capacitor of capacitance C. The telegraph equation is concerned with the changes in voltage between the ends of the piece of the wire which can be formulated as follows [5]:
(1)∂2y∂t2(x,t)+GC+RL∂y∂t(x,t)+GCRLy(x,t)=1LC∂2y∂x2(x,t)+f(x,t),(x,t)∈[0,1]×[0,1],
with initial and boundary conditions:
(2)y(x,0)=g0(x),∂y∂t(x,0)=g1(x),(3)y(0,t)=f0(t),y(1,t)=f1(t),
where y(x,t) denotes the voltage at position x and time t and f(x,t) is the external source term. Suppose ρ=G/C, σ=R/L and η2=1/LC; thus, we have
(4)∂2y∂t2(x,t)+ρ+σ∂y∂t(x,t)+ρσy(x,t)=η2∂2y∂x2(x,t)+f(x,t).
Because of the lack of appropriate mathematical methods, most of the analytical solutions for linear and nonlinear partial differential equations are challenging to acquire. Therefore, approximation and numerical techniques such as Adomian decomposition method, local radial basis collocation method, variational iteration, homotopy perturbation, and Laplace transform method have been applied [6–14]. The numerical methods for solving the second-order telegraph equation are well developed. Lakestani and Saray [1] solved this problem through expanding the obtained approximate solution as the elements of interpolating scaling functions. Furthermore, the numerical approximation based on differential transform method (DTM) was considered to solve telegraph equation [5]. Using DTM, it is possible to find the exact solution or a closed approximate solution for an equation. The fully discrete local discontinuous Galerkin finite element method based on a finite difference scheme in time was introduced by [15]. Chen et al. [16] used the method of separation of variables for the analytical solution of the nonhomogeneous telegraph equation under three types of nonhomogeneous boundary conditions.We are inspired to have an algorithm which not only is appropriate for long-time calculations but also reflects the global behavior of exact solutions. The novelty of the present paper is that we investigate the behaviour of the combination of two different groups of orthogonal functions from two different intervals. The combination of piecewise orthogonal Haar functions defined on[0,1) interval, with continuous orthogonal Sinc functions defined on (-∞,+∞) interval. It is essential to be said that [17–24] have previously used the Haar and Sinc functions separately for solving optimal control problems and some nonlinear ordinary differential equations. For example, Karimi and Lohmann [17] applied the Haar functions for modeling and robust control of bounce and pitch vibration for the enginebody vibration structure. Also, the Haar wavelet method has been investigated for optimal control of time-varying state-delayed [18], linear singularly perturbed systems [19], and second-order linear systems [20]. In addition, based on the properties of orthogonal Sinc functions, it is apparent that the convergence rate of approximation is exponential [21]. By using this property, the authors of [22–24] studied the Sinc collocation method for solving nonlinear singular equations like Thomas-Fermi, Lane-Emden, and Blasius equations.Another advantage of proposed method is that it transforms the problem into a system of algebraic equations, so that the computation becomes simple and computer oriented. In the new proposed algorithm, we extend the solution of the problem to the sum of basis functions and take good advantage of the orthogonality of Haar and Sinc functions to build a set of equations for the coefficients of the solution. Although the test model is a very simple one, the proposed method is also applicable to many other problems such as fractional and two-dimensional nonlinear PDEs.The organization of the rest of the paper is as follows. In Section2, we present a brief introduction to the essential definitions of the Haar and Sinc functions from which are derived some tools for developing our method. In Section 3, the convergence rate analysis of the Haar and Sinc functions is given. In Section 4, we apply the method of Haar-Sinc collocation for solving the model equation. In Section 5, the proposed method is used in some types of telegraph equations and it is compared with the current analytic solutions revealed in different published works within the literature. The conclusion is presented in the final section.


## 2. Basic Definitions

### 2.1. Haar Functions

The orthogonal set of Haar functions is a number of square waves with magnitude of±1 in some intervals and zeros elsewhere [25]. The Haar functions are defined on the interval C=[0,b) by
(5)hu(t)=1,l1≤t<l1/2,-1,l1/2≤t<l0,0,otherwise,
where,
(6)lq=n-q2mb,q=0,12,1.
The value of u is defined by two parameters m and n as
(7)u=2m+n-1,m=0,1,2,…,M,n=1,2,3,…,2m.
The integer m indicates the level of the wavelet and the maximal level of resolution is the integer M. Also, h0(t) is defined for m=n=0 and is given by
(8)h0(t)=1,0≤t<b.
We can expand any function y(t)∈L2(0,1) in first k terms of Haar functions as
(9)y(t)≃∑u=0k-1‍xuhu(t)=XTϕk(t),
where
(10)k=2s+1,s=0,1,2,….
The Haar functions coefficient vector X and Haar functions vector ϕk(t) are defined as
(11)X=x0,x1,…,xk-1T,(12)ϕk(t)=h0t,h1t,…,hk-1(t)T.
The matrix Φk×k can be expressed as
(13)Φk×k=ϕk12k,ϕk32k,…,ϕk2k-12k.
Furthermore, the integration of the ϕk(t) defined in (12) is given by
(14)∫0t‍ϕk(s)ds≃Jϕk(t),
where J=Jk×k is the k×k operational matrix for integration and is given in [26] as
(15)Jk×k=b2k2kJ(k/2)×(k/2)-Φ(k/2)×(k/2)Φ(k/2)×(k/2)-10,
where Φ1=[1], J1=[1/2].


### 2.2. Sinc Functions

The Sinc function is defined on the whole real lineI=(-∞,∞) by
(16)Sinc(x)=sin(πx)πx,x≠0,1,x=0.
For each integer m and the mesh size h, the Sinc basis functions are defined on R by [22]
(17)Smh,x≡Sincx-mhh=sin⁡π/hx-mhπ/hx-mh,x≠mh,1,x=mh.
The Sinc functions form an interpolatory set of functions; that is,
(18)Sm(h,jh)=δmj=1,j=m,0,j≠m.
If a function f(x) is defined on the real axis, then for h>0 the series,
(19)C(f,h)(x)=∑m=-∞∞‍fmhSincx-mhh,
is called the Whittaker cardinal expansion of f whenever this series converges. The properties of the Whittaker cardinal expansion have been extensively studied in [27]. These properties are derived in the infinite strip DS of the complex ω-plane, where for d>0,
(20)DS={ω=t+is:s<d≤π2}.
Approximations can be constructed for infinite, semiinfinite, and finite intervals. To construct approximations on the interval (0,1) which is used in this paper, the eye-shaped domain in thez-plane,
(21)DE={z=x+iy:argz1-z<d≤π2},
is mapped conformally onto the infinite strip DS via
(22)ω=ψ(z)=ln⁡z1-z.
The basis functions on (0,1) are taken to be the composite translated Sinc functions:
(23)Sm(x)≡S(m,h)∘ψ(x)=Sincψx-mhh,
where S(m,h)∘ψ(x) is defined by S(m,h)(ψ(x)). The inverse map of ω=ψ(z) is
(24)z=ψ-1(ω)=eω1+eω.
Thus we may define the inverse images of the real line and of the evenly spaced nodes {mh}m=-∞m=+∞ as
(25)Γ={ψ-1(u)∈DE:-∞<u<+∞}=(0,1),(26)xm=ψ-1mh=emh1+emh,m=0,±1,±2,….
Also, thenth derivative of the function f at some points xm can be approximated [24]:
(27)δm,j(0)=S(m,h)∘ψ(x)x=xj=1,m=j,0,m≠j,(28)δm,j(1)=ddψS(m,h)∘ψ(x)x=xj=1h0,m=j,-1j-mj-m,m≠j,(29)δm,j(2)=d2dψ2S(m,h)∘ψ(x)x=xj=1h2-π23,m=j,-2-1j-mj-m2,m≠j.


## 2.1. Haar Functions

The orthogonal set of Haar functions is a number of square waves with magnitude of±1 in some intervals and zeros elsewhere [25]. The Haar functions are defined on the interval C=[0,b) by
(5)hu(t)=1,l1≤t<l1/2,-1,l1/2≤t<l0,0,otherwise,
where,
(6)lq=n-q2mb,q=0,12,1.
The value of u is defined by two parameters m and n as
(7)u=2m+n-1,m=0,1,2,…,M,n=1,2,3,…,2m.
The integer m indicates the level of the wavelet and the maximal level of resolution is the integer M. Also, h0(t) is defined for m=n=0 and is given by
(8)h0(t)=1,0≤t<b.
We can expand any function y(t)∈L2(0,1) in first k terms of Haar functions as
(9)y(t)≃∑u=0k-1‍xuhu(t)=XTϕk(t),
where
(10)k=2s+1,s=0,1,2,….
The Haar functions coefficient vector X and Haar functions vector ϕk(t) are defined as
(11)X=x0,x1,…,xk-1T,(12)ϕk(t)=h0t,h1t,…,hk-1(t)T.
The matrix Φk×k can be expressed as
(13)Φk×k=ϕk12k,ϕk32k,…,ϕk2k-12k.
Furthermore, the integration of the ϕk(t) defined in (12) is given by
(14)∫0t‍ϕk(s)ds≃Jϕk(t),
where J=Jk×k is the k×k operational matrix for integration and is given in [26] as
(15)Jk×k=b2k2kJ(k/2)×(k/2)-Φ(k/2)×(k/2)Φ(k/2)×(k/2)-10,
where Φ1=[1], J1=[1/2].


## 2.2. Sinc Functions

The Sinc function is defined on the whole real lineI=(-∞,∞) by
(16)Sinc(x)=sin(πx)πx,x≠0,1,x=0.
For each integer m and the mesh size h, the Sinc basis functions are defined on R by [22]
(17)Smh,x≡Sincx-mhh=sin⁡π/hx-mhπ/hx-mh,x≠mh,1,x=mh.
The Sinc functions form an interpolatory set of functions; that is,
(18)Sm(h,jh)=δmj=1,j=m,0,j≠m.
If a function f(x) is defined on the real axis, then for h>0 the series,
(19)C(f,h)(x)=∑m=-∞∞‍fmhSincx-mhh,
is called the Whittaker cardinal expansion of f whenever this series converges. The properties of the Whittaker cardinal expansion have been extensively studied in [27]. These properties are derived in the infinite strip DS of the complex ω-plane, where for d>0,
(20)DS={ω=t+is:s<d≤π2}.
Approximations can be constructed for infinite, semiinfinite, and finite intervals. To construct approximations on the interval (0,1) which is used in this paper, the eye-shaped domain in thez-plane,
(21)DE={z=x+iy:argz1-z<d≤π2},
is mapped conformally onto the infinite strip DS via
(22)ω=ψ(z)=ln⁡z1-z.
The basis functions on (0,1) are taken to be the composite translated Sinc functions:
(23)Sm(x)≡S(m,h)∘ψ(x)=Sincψx-mhh,
where S(m,h)∘ψ(x) is defined by S(m,h)(ψ(x)). The inverse map of ω=ψ(z) is
(24)z=ψ-1(ω)=eω1+eω.
Thus we may define the inverse images of the real line and of the evenly spaced nodes {mh}m=-∞m=+∞ as
(25)Γ={ψ-1(u)∈DE:-∞<u<+∞}=(0,1),(26)xm=ψ-1mh=emh1+emh,m=0,±1,±2,….
Also, thenth derivative of the function f at some points xm can be approximated [24]:
(27)δm,j(0)=S(m,h)∘ψ(x)x=xj=1,m=j,0,m≠j,(28)δm,j(1)=ddψS(m,h)∘ψ(x)x=xj=1h0,m=j,-1j-mj-m,m≠j,(29)δm,j(2)=d2dψ2S(m,h)∘ψ(x)x=xj=1h2-π23,m=j,-2-1j-mj-m2,m≠j.


## 3. Convergence Rate Analysis

### 3.1. Haar Functions

Theorem 1.
Assume thaty(t)∈L2(R) with the bounded first derivative on (0,1); then, the error norm at Mth level satisfies the following inequality:
(30)eM(t)≤K7C2(-3/2)M,
where K,C are some real constants [28, 29].Proof.
The error atMth level may be defined as
(31)eMt=yt-yMt=∑i=2M+1+1∞cihi(t),eM(t)2=∫-∞∞‍∑i=2M+1+1∞cihi(t),∑q=2M+1+1∞cqhq(t)eMt2=∑i=2M+1+1∞∑q=2M+1+1∞‍cicq∫-∞∞‍hi(t)hq(t)dt,eM(t)2≤∑i=2M+1+1∞ci2,
where yM(t)=∑i=12M+1cihi(t). But |ci|≤C2-3i/2max⁡|y′(l)|, where C=∫01‍|th2(t)|dt and l∈(k2-m,(k+1)2-m). Then
(32)eM(t)2≤∑i=2M+1+1∞KC22-3i,eM(t)2≤KC2172-3M,eM(t)≤K7C2-(3/2)M,
where y′(t)≤K∀t∈(0,1) and K is positive constant.


### 3.2. Sinc Functions

The following theorem for which the proof can be found in [21] shows that the convergence rate of Sinc approximation is exponential.Definition 2.
LetH2(DE) be the class of functions f which are analytic in DE (the eye-shaped domain defined in (21)) satisfy
(33)∫ψ-1(p+T)fzdz⟶0,x⟶±∞,
where T={iq:|q|<d≤π/2} and the function f satisfies the following equation on the boundary of DE:
(34)N(f)=∫∂Dfzdz<∞.Theorem 3.
Assume thatfψ′∈H2(DE); then, for all z∈(0,1),
(35)Ef,hz=f(z)-∑k=-∞∞f(kh)S(k,h)∘ψ(z)≤N(fψ′)2πdsinh⁡(πd/h)≤2N(fψ′)πde-πd/h.
Moreover, if |fz|≤Ce-α|ψz|,z∈Γ for some positive constants C, α if the selection h=πd/αN≤2πd/ln⁡(2), then
(36)f(z)-∑k=-NNf(kh)S(k,h)∘ψ(z)≤KNe(-πdαN),
where K depends only on f, d, and α.


## 3.1. Haar Functions

Theorem 1.
Assume thaty(t)∈L2(R) with the bounded first derivative on (0,1); then, the error norm at Mth level satisfies the following inequality:
(30)eM(t)≤K7C2(-3/2)M,
where K,C are some real constants [28, 29].Proof.
The error atMth level may be defined as
(31)eMt=yt-yMt=∑i=2M+1+1∞cihi(t),eM(t)2=∫-∞∞‍∑i=2M+1+1∞cihi(t),∑q=2M+1+1∞cqhq(t)eMt2=∑i=2M+1+1∞∑q=2M+1+1∞‍cicq∫-∞∞‍hi(t)hq(t)dt,eM(t)2≤∑i=2M+1+1∞ci2,
where yM(t)=∑i=12M+1cihi(t). But |ci|≤C2-3i/2max⁡|y′(l)|, where C=∫01‍|th2(t)|dt and l∈(k2-m,(k+1)2-m). Then
(32)eM(t)2≤∑i=2M+1+1∞KC22-3i,eM(t)2≤KC2172-3M,eM(t)≤K7C2-(3/2)M,
where y′(t)≤K∀t∈(0,1) and K is positive constant.


## 3.2. Sinc Functions

The following theorem for which the proof can be found in [21] shows that the convergence rate of Sinc approximation is exponential.Definition 2.
LetH2(DE) be the class of functions f which are analytic in DE (the eye-shaped domain defined in (21)) satisfy
(33)∫ψ-1(p+T)fzdz⟶0,x⟶±∞,
where T={iq:|q|<d≤π/2} and the function f satisfies the following equation on the boundary of DE:
(34)N(f)=∫∂Dfzdz<∞.Theorem 3.
Assume thatfψ′∈H2(DE); then, for all z∈(0,1),
(35)Ef,hz=f(z)-∑k=-∞∞f(kh)S(k,h)∘ψ(z)≤N(fψ′)2πdsinh⁡(πd/h)≤2N(fψ′)πde-πd/h.
Moreover, if |fz|≤Ce-α|ψz|,z∈Γ for some positive constants C, α if the selection h=πd/αN≤2πd/ln⁡(2), then
(36)f(z)-∑k=-NNf(kh)S(k,h)∘ψ(z)≤KNe(-πdαN),
where K depends only on f, d, and α.


## 4. Haar-Sinc Collocation Method

A discrete approximation to the(∂2y/∂t2)(x,t) can be expanded into 2n+1 Sinc functions and k Haar functions as
(37)∂2yn,k∂t2(x,t)=∑i=-nn∑j=0k-1cijSi(x)hj(t).Lemma 4.
Letxm be Sinc collocation points, given in (26). Then the following relations hold:
(38)∂yn,k∂t(xm,t)=g1(xm)+∑j=0k-1∑l=0k-1cmlJljhj(t),yn,kxm,t=g0(xm)+g1(xm)t+∑j=0k-1∑l=0k-1cmlJlj2hj(t),∂2yn,k∂x2xm,t=∂2g0(x)+g1(x)t∂x2xm+∑i=-nn∑j=0k-1∑l=0n-1wim(2)cilJlj2hj(t),
where
(39)Jlj=Jk×kl,j,wim(2)=ψ′′(xm)δi,m(1)+ψ′xm2δi,m(2).Proof.
Employing (2), (14), (27), and (37) we have
(40)∂yn,k∂txm,t=g1(xm)+∫0t∑i=-nn∑j=0k-1cijSi(xm)hj(t)dt=g1(xm)+∫0t∑i=-nn∑j=0k-1ci,jδi,m(0)hj(t)dt=g1(xm)+∫0t∑j=0k-1cmjhj(t)dt=g1(xm)+∑j=0k-1∑l=0k-1cmlJljhj(t).
Also, using (40) we get
(41)yn,kxm,t=g0(xm)+∫0t∂y∂t(x,t)dtxm=g0(xm)+∫0tg1(xm)dt+∫0t∑j=0k-1∑l=0k-1cmlJljhj(t)dt=g0(xm)+g1xmt+∑j=0k-1∑l=0k-1cmlJlj2hj(t).
In addition, using (28), (29), and (41) we have
(42)∂2yn,k∂x2xm,t=∂2g0(x)+g1(x)t∂x2xm︷A+∑i=-nn∑j=0k-1∑l=0k-1cilJlj2d2Sixdx2x=xmhj(t)=A+∑i=-nn∑j=0k-1∑l=0k-1cilJlj2ψ′′xddϕSixhhhhhhhhhhhhhhhhhihh+ψ′x2ddϕ2Si(x)x=xmhj(t)=A+∑i=-nn∑j=0k-1∑l=0k-1cilJlj2ψ′xm2ψ′′xmδi,m(1)hhhhhhhhhhhhhhhhhhh+ψ′xm2δi,k(2)hj(t),
which completes our proof.The residualRn,k(x,t) for (4) can be written as
(43)Rn,kx,t=∂2yn,k∂t2x,t+ρ+σ∂yn,k∂t(x,t)+ρσyn,k(x,t)-η2∂2yn,k∂x2(x,t)-f(x,t).
The equations for obtaining the (2n+1)k coefficients {cij}i=-n⋯nj=0⋯k-1 arise from equalizing Rn,k(x,t) to zero at 2n+1 Sinc points and k Haar points are defined by
(44)xm=ψ-1(mh)=emh1+emh,m=0,±1,±2,…,±n,ti=2i-12k,i=1,2,…,k.
By substitution collocation points in Rn,k(x,t) and equalizing to zero we have
(45)Rn,k(xm,ti)=0,m=-n,…,n,i=1,2,…,k.
Equation (45) gives (2n+1)k linear algebraic equations which can be solved for the unknown coefficients {cij}i=-n⋯nj=0⋯k by using the Newton’s method. Consequently, y(x,t) given in (4) can be calculated.


## 5. Illustrative Examples

In this section, we apply the proposed method for solving (4) and show the efficiency of the method with the numerical results of some examples. In all examples we choose h=π/2n.Example 1.
Consider the linear telegraph equation [1]:
(46)∂2y∂x2-4e-2tsinh⁡(x)=∂2y∂t2+4∂y∂t+y,
with the following initial and boundary conditions:
(47)y(x,0)=sinh⁡(x),yt(x,0)=-2sinh⁡(x),y(0,t)=0,y(1,t)=e-2tsinh⁡(1).
The exact solution to this problem is
(48)y(x,t)=e-2tsinh⁡(x).
Table 1 shows the absolute error function |yexactx,t-yn,kx,t| obtained by the present method with n=3 and different values of k.
Also, Figure1 displays the convergence rate of our method with k=4,8,16,32 for t=0.5. It is seen from the Figure 1 that for each fixed point (x,t) the absolute errors get smaller and smaller as k increases. Furthermore, we can see that the presented method provides accurate results even by using n=3.
The maximum absolute errors fork=8 and different values of n (Sinc collocation points) are shown graphically in Figure 2 for t=0.5. In Figure 2, we observe that the values of maximum absolute error decay exponentially as expected from Theorem 3.Table 1
The values of absolute error withn=3 for Example 1.
x
t
|
y
exact
(
x
,
t
)
-
y
n
,
k
(
x
,
t
)
|
k
=
8
k
=
16
k
=
32
0.3
0.25
2.72
e
-
3
1.31
e
-
3
6.56
e
-
4
0.50
3.97
e
-
3
2.01
e
-
3
1.04
e
-
3
0.75
4.65
e
-
3
2.34
e
-
3
1.19
e
-
3
0.6
0.25
3.72
e
-
3
1.75
e
-
3
8.41
e
-
4
0.50
5.32
e
-
3
2.57
e
-
3
1.21
e
-
3
0.75
6.08
e
-
3
2.81
e
-
3
1.18
e
-
3
0.9
0.25
1.51
e
-
3
6.84
e
-
4
2.89
e
-
4
0.50
1.84
e
-
3
6.57
e
-
4
7.14
e
-
5
0.75
1.61
e
-
3
1.81
e
-
4
5.27
e
-
5Figure 1
Plot of the absolute errors for different values ofk with t=0.5 and n=3 for Example 1.Figure 2
Plot of the maximum absolute errors for different values ofn with t=0.5 and k=8 for Example 1.Example 2.
Consider the linear telegraph equation [5]:
(49)∂2y∂x2=∂2y∂t2+4∂y∂t+4y,
with the following initial and boundary conditions:
(50)y(x,0)=ex,yt(x,0)=-ex,y(0,t)=e-t,y(1,t)=e1-t.
The exact solution is given by
(51)y(x,t)=ex-t.
Table 2 shows the absolute error values using the proposed method with n=3 and k=8,16,32.
Figure3 displays the values of maximum absolute error with k=4,8,16,32 for t=0.5,0.75. This figure demonstrates the validity and applicability of the present technique for this problem. Also, we observe that the convergence rate of our method for t=0.75 is lower than t=0.5.
Figure4 displays the convergence rate of our method for different values of n for t=0.75 and k=8. We see from Figure 4 that our method is in good agreement with the actual rate of convergence and the values of maximum absolute error decay exponentially.Table 2
The values of absolute error withn=3 and k=8,16,32 for Example 2.
x
t
=
0.5
t
=
0.75
k
=
8
k
=
16
k
=
32
k
=
8
k
=
16
k
=
32
0.1
1.22
e
-
3
4.16
e
-
4
2.12
e
-
5
1.17
e
-
3
3.83
e
-
4
1.32
e
-
4
0.2
2.85
e
-
3
1.30
e
-
3
5.44
e
-
4
3.38
e
-
3
1.38
e
-
3
3.98
e
-
4
0.3
4.22
e
-
3
2.08
e
-
3
1.03
e
-
3
5.35
e
-
3
2.60
e
-
3
1.23
e
-
3
0.4
5.06
e
-
3
2.52
e
-
3
1.28
e
-
3
6.50
e
-
3
3.24
e
-
3
1.62
e
-
3
0.5
5.38
e
-
3
2.65
e
-
3
1.32
e
-
3
6.89
e
-
3
3.39
e
-
3
1.65
e
-
3
0.6
5.25
e
-
3
2.55
e
-
3
1.24
e
-
3
6.64
e
-
3
3.17
e
-
3
1.46
e
-
3
0.7
4.64
e
-
3
2.22
e
-
3
1.04
e
-
4
5.77
e
-
3
2.64
e
-
3
1.09
e
-
3
0.8
3.47
e
-
3
1.59
e
-
3
6.72
e
-
4
4.13
e
-
3
1.69
e
-
3
4.83
e
-
4
0.9
1.66
e
-
3
6.06
e
-
4
8.98
e
-
5
1.68
e
-
3
3.76
e
-
4
3.02
e
-
4Figure 3
Plot of the maximum absolute errors for different values ofk with t=0.5,0.75 and n=3 for Example 2.Figure 4
Plot of the maximum absolute errors for different values ofn with t=0.75 and k=8 for Example 2.Example 3.
Consider the nonlinear telegraph equation [5]:
(52)∂2y∂x2=∂2y∂t2+2∂y∂t+y2-e2x-4t+ex-2t,
with the following initial and boundary conditions:
(53)yx,0=ex,ytx,0=-2ex,y(0,t)=e-2t,y(1,t)=e1-2t.
The exact solution is given by
(54)y(x,t)=ex-2t.
Figures 5, 6, 7, and 8 show the absolute error function |yexactx,t-yn,kx,t| obtained by the present method with n=3 and k=4,8,16 and 32. We can see clearly that better accuracy can be achieved by increasing the number of Haar collocation points and using the arbitrary precision ability of Mathematica software, we are able to establish more accurate results.Figure 5
Plot of the absolute errore(x,t)=|yx,t-yn,kx,t| with n=3,k=4 for Example 3.Figure 6
Plot of the absolute errore(x,t)=|yx,t-yn,kx,t| with n=3,k=8 for Example 3.Figure 7
Plot of the absolute errore(x,t)=|yx,t-yn,kx,t| with n=3,k=16 for Example 3.Figure 8
Plot of the absolute errore(x,t)=|yx,t-yn,kx,t| with n=3,k=32 for Example 3.Example 4.
Consider the numerical computation of the nonlinear telegraph equations [31]:
(55)∂2y∂x2=∂2y∂t2+2∂y∂t-y3+y.
The exact solution is given by
(56)y(x,t)=12+12tanhx8+3t8+5.
We extract the initial and boundary conditions from the exact solution. Table 3 shows the comparison of the L∞ error between approximations obtained by the radial basis functions (IMQ, TSP) [30] and the Haar-Sinc collocation method with n=3 and k=8,16,32. The results obtained in the table show that the L∞ error between the numerical and the exact solution can be reduced by increasing the values of k.Table 3
Results for the value ofL∞ error obtained by Haar-Sinc collocation method with n=3 and k=8,16,32 for Example 4.
t
L
∞
- error [30]
L
∞
- error (present method)
TPS
IMQ
k
=
8
k
=
16
k
=
32
0.5
1.014
e
-
5
3.599
e
-
6
1.725
e
-
6
6.598
e
-
7
7.91
e
-
8
1.0
1.668
e
-
5
6.375
e
-
7
4.237
e
-
6
3.440
e
-
6
3.93
e
-
7
1.5
1.087
e
-
5
1.156
e
-
6
7.089
e
-
6
6.578
e
-
6
6.31
e
-
7
2.0
3.633
e
-
6
6.750
e
-
7
7.248
e
-
6
6.890
e
-
6
7.13
e
-
7
3.0
2.159
e
-
6
4.154
e
-
7
7.873
e
-
6
7.869
e
-
6
7.18
e
-
7


## 6. Conclusion

A numerical method for solving the telegraph equations based on the combination of two orthogonal Haar and Sinc functions was proposed. Also, [32, 33] have previously applied the combination of the Sinc functions with the other different groups of orthogonal Legendre and Chebyshev polynomials for solving similar PDEs. It is worth mentioning that since the convergence rate of the Sinc approximation is exponential, we expand the problem in space with the Sinc basis functions. Furthermore, according to the initial conditions of the problem and the Haar operational matrix for integration, the approximation in time was expanded via the elements of Haar functions. This method can be applied to solve similar problems in physics and provides a powerful alternative for physicians to investigate such types of nonlinear PDEs. The effectiveness of the method was examined via comparing the obtained results with the exact solutions. Based on the numerical results, it is obvious that it can more entirely simulate the global property of the exact solution and can provide more information about the structures of the problem. Also, the absolute error may be decreased if we take more collocation points.



---

*Source: 340752-2014-11-16.xml*