# A Coupled Pseudospectral-Differential Quadrature Method for a Class of Hyperbolic Telegraph Equations

Authors: Fangzong Wang; Yong Wang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Pseudospectral methods and differential quadrature methods are two kinds of important meshless methods, both of which have been widely used in scientific and engineering calculation. The Lagrange interpolation polynomials are used as the trial function of the two methods, and the same distribution of grid points is used. This paper points out that the differential quadrature method is a special form of the pseudospectral method. On the basis of the above, a coupled pseudospectral-differential quadrature method (PSDQM) is proposed to solve a class of hyperbolic telegraph equations. Theoretical analysis and numerical tests show that the new method has spectral precision convergence in spatial domain and has A-stability in time domain. And it is suitable for solving multidimensional telegraph equations.


---

## Body

## 1. Introduction

The spectral methods [1] are a class of new numerical methods proposed by Gottlieb and Orszag for solving partial differential equations. As a kind of meshless methods, it can be essentially expressed as using special orthogonal polynomials (such as Chebyshev and Legendre polynomials) as test functions to discretize partial differential equations. Specifically, it is to define a set of orthogonal functions as basis functions on the calculation domain and use basis functions to approximate the desired variable as its spectral approximation form. Its notable feature is that only a few discrete points are needed to get a higher approximation accuracy [2]. Spectral methods can be divided into the Galerkin method and the spectral collocation method (SCM) or the pseudospectral method (PSM) according to the different approximation schemes [1]. For the latter, the differentiation matrix [3] is used to compute derivatives at the collocation points, which obtains great convenience for the SCM. So the differentiation matrix plays an important role in the SCM for numerical calculation [3].The differential quadrature methods (DQM) [4] also are a class of meshless methods for solving differential or partial differential equations, which were proposed by Bellman and Casti in the early 1970s as an extension of the idea of conventional integral quadrature. The differential quadrature method was proved to be essentially equivalent to the general collocation method [5]. The basic idea of DQM is to use the weighted linear sum of the function values of all discrete points in a coordinate direction over the whole domain to represent the derivative of the function at that discrete point. The differential quadrature method is simple in mathematical principle and its implementation does not depend on the variational principle [5]. In addition, it is also considered to be easy in programming and achieve higher computational efficiency. In order to solve the weighting coefficients or the weighting coefficients matrix of differential quadrature method, the researchers have deduced the implicit matrix expression and the explicit recursive algebraic formulation, respectively, using the general polynomial and the Lagrange interpolation polynomial as the trial function [6–8]. In [7], it was pointed out that the weighting coefficients of DQM depend only on the choice of trial function and the distribution of grid points and have nothing to do with the actual research problems.In this study, a coupled pseudospectral-differential quadrature method is presented for solving a class of hyperbolic telegraph equations. The hyperbolic telegraph equation is first reduced to a set of ordinary differential equations by spatial discretization based on the spectral collocation method. Then thes-stages-order time domain differential quadrature method is employed to solve the resulting initial value problem. Finally, the resulting algebraic equations can be solved by either direct or iterative methods. As a comparison, the method in [9] is also used for solving the telegraph equation.The rest of the paper is arranged as follows. In Section2, the spectral collocation method or pseudospectral method is briefly introduced. In Section 3, the space domain differential quadrature method is firstly introduced. Then, we point out the equivalence relationship of pseudospectral method and differential quadrature method in the spatial domain. In Section 4, we deduce how to use the proposed coupled pseudospectral-differential quadrature methods to solve one-dimensional and two-dimensional hyperbolic telegraph equations in detail. In Section 5, three specific numerical examples are computed, which are given to verify the computational accuracy and numerical stability of the proposed method. Conclusions are then given in Section 6.


## 2. Pseudospectral Method

Pseudospectral method, also known as the collocation points method [1, 2, 9], is an efficient and meshless numerical method for solving partial differential equation or equations. Without loss of generality, consider the differential equation u(x) of a one-dimensional function:(1)Lux=fx,x∈c,d,where L is a linear differential operator and f(x) is a given smooth function of space; u(x) is a variable that is unknown. x is one point in the interval [c,d]. The core idea of the spectral method and pseudospectral method is consistent; namely, the variable u(x) can be interpolated and approximated by the truncated series of basis function ϕn(x); that is,(2)ux≅uNx=∑n=0Nαnϕnx,where N is the number of divided intervals; αn is the spectral coefficients. uN(x) is the approximate value of u(x).To measure the approximate effect of (2), define the approximate residual value r(x,α) as follows: (3)rx,α=LuN-fx.We can see from the above analysis that the essential part of spectral method is to find spectral coefficients  αn, which make the approximate residual values r(x,α) as small as possible.In general, basis functions include discrete Fourier series, Chebyshev polynomials, Lagrange interpolation basis function, and so on [3]. The collocation points include the following four typical types (take standardized interval [-1,1], for example), which are defined as follows [4]:(1)
Uniform grid points(4)xi=2iN-1,i∈0,N.(2)
Chebyshev-Gauss-Lobatto grid points(5)xi=cos⁡iπN,i∈0,N.(3)
Gauss-Legendre grid points(6)x0=-1;xi,i∈1,N-1are  the  zero  points  ofdN-1x2-1N-1dxN-1=0;xN=1.(4)
Gauss-Lobatto grid points(7)xi=121-cos⁡iπN,i∈0,N,where xi∈0,1,i∈0,N. Therefore, the Gauss-Lobatto grid points on the interval -1,1 can be obtained by a linear transformation.Using the collocation points method to obtain the approximate value ofu(x), each collocation point xii∈0,N satisfies (3); that is,(8)rxi,αi=LuNxi-fxi=0.Solving the algebraic equations about spectral coefficientsαi, the approximate solution uN of the variable u(x) is obtained.To avoid solving directly (8), consider the PSM or collocation point s method with the collocation points x0,…,xN and the corresponding approximate value u0,…,uN. Denote Luxi=fi,i∈0,N; the p,p≥1 order differential matrix DNp is defined in the corresponding collocation points under the linear differential operator L; that is, (9)F=DNpV,where F=f0,…,fNT; V=u0,…,uNT.From formula (9), it can be seen that using the differentiation matrix DNp brings great convenience for numerical calculation. The solution procedure of the differentiation matrices of the pseudospectral method is given in [3] in detail. It is necessary to point out that the differentiation matrix formation of the pseudospectral method is related to the interpolation basis function and the collocation points, and it has nothing to do with the specific problem [3]. In other words, if basis function and collocation points are given, the unique differentiation matrices can be formed. The p-order differential matrix DNp of pseudospectral method with the Lagrange interpolation basis function and the space discrete points using the Chebyshev-Gauss-Lobatto grid points is presented as follows [3, 10]:(10)DNp=d00p⋯d0Np⋮⋱⋮dN0p⋯dNNp∈RN+1×N+1,where, in the case of p=1, and(11)d001=2N2+16;dNN1=-2N2+16;dii1=-xi21-xi2,i=1,2,…,N-1;dij1=γiγj-1i+jxi-xj,i≠j,i,j=0,1,…,N-1;γi=2i=0,N1i=2,…,N-1,whenN=1,2,3, the first-order differential matrices are presented as follows:(12)D11=12-1212-12;D21=32-212120-12-122-32;D31=196-443-121-13-113-13113-112-434-196.The second-order differentiation matrixDN2 of the PSM can be obtained directly by the following formula [3, 10]:(13)DN2=DN1DN1.When N=1,2,3, the second-order differentiation matrices are given in the formula (13):(14)D12=0000;D22=1-211-211-21;D22=163-283203-83103-16383-23-2383-163103-83203-283163.For convenience, the above method is named Chebyshev pseudospectral method in next sections.


## 3. Differential Quadrature Method

Differential quadrature method is also a kind of meshless numerical method [11]. Without loss of generality, suppose one-dimensional function g(x) is sufficiently smooth in the whole interval a,b; and the interval a,b is divided into M sections; that is,(15)Lgx=∑i=0MAixgxi,where L means linear differential operator; Ai(x)i∈(0,M) is called the interpolation basis function; a=x0≤xi≤xM=b, xi is the (i+1)th discrete point.IfLg(x)=g(p)(x), (15) can be rewritten as the following formula:(16)dpgxdxpx=xi=∑j=0Mwijpgxj,i∈0,M,where wijp is the p-order derivative weighting coefficient s of DQM.In DQM, using general polynomial as test functions, it is difficult to solve the weighting coefficients whenM increases to a certain value [8]. Using Lagrange interpolation polynomial as test functions, the explicit expression of the weighting coefficients is shown as follows (take standardized interval [-1,1], for example) [6, 8]:(17)qij=Jxixi-xjJxj,i≠j,i,j=0,1,…,M,qii=-∑i=0,i≠jMqij,i=j,i=0,1,…,M,p=1,where T(x)=∏k=0M(x-xk), J(x)=T′(x).(18)wijp=pqijwiip-1-wijp-1xi-xj,i≠j,i,j∈0,M,wiip=-∑i=0,i≠jMwijp,i=j,i=0,M,p≥2.Shu has proved that the weighting coefficients are equivalent either using general polynomial or Lagrange interpolation polynomial as test functions.In general, considering the grid pointsσi,i∈(0,M), in the regularized interval -1,1, obviously, the transformation between xi and σi can be achieved by the following linear transformation:(19)xi=b-a2σi+b+a2.Similarly, the coefficientsξij,(i,j=0,1,…,M) on the general interval a,b can be obtained by using the following transformation formula [12]:(20)ξij=2b-apψij,where ψij,(i,j=0,1,…,M) is the p-order derivative weighting coefficient of the regularized region -1,1.The weighting coefficients matrixWMp of the p-order derivative of g(x) is defined as follows:(21)WMp=w00p⋯w0Mp⋮⋱⋮wM0p⋯wMMp∈RM+1×M+1.DefineLg(x)=g(p)(x) and Lg(xi)=hi,i∈(0,M); (15) can be rewritten as the following matrix form:(22)H=WMpG,where H=h0,h1,…,hMT; G=g(x0),g(x1),…,g(xM)T.Whenp=1, M=1,2,3, the first-order weighting coefficients matrices with the grid points using the Chebyshev points in the regularized region -1,1 are shown as follows:(23)W11=12-1212-12;W21=32-212120-12-122-32;W31=196-443-121-13-113-13113-112-434-196.Similarly, the second-order derivative weighting coefficients matrices can expressed as follows:(24)W12=0000;W22=1-211-211-21;W32=163-283203-83103-16383-23-2383-163103-83203-283163.Obviously, the above DQM is mainly for the numerical computation in the spatial domain. The basic content of time domain differential quadrature method has been analyzed in detail in [5, 7]. Reference [7] indicates that classicals-stages-order time domain DQM can be transformed into equivalent implicit Runge-Kutta method. From the definition, test functions, interpolation points, and differentiation matrix or weighting matrix of the DQM and PSM, in the spatial domain, and the differential matrices of the PSM are equivalent to the weighting coefficients matrices of the DQM as long as the test functions and the interpolation points are consistent with each other. In other words, the PSM and the DQM are also equivalent in the spatial domain.Inspired by this, a new coupled numerical method is proposed, which is called pseudospectral-differential quadrature methods (PSDQM). In this paper, PSDQM is used to solve the hyperbolic telegraph equations. And the effectiveness of PSDQM is illustrated by comparing with the method used in [9].In the following, the four-stage time domain differential quadrature method using uniform grid points will be given as an example; that is,(25)W~41=-1036-213-83083-132-61031-16312-16253;W0=-W~41e=-113-131T,where e is a column vector of dimension four, all components of which are one.


## 4. Solving Hyperbolic Telegraph Equations By PSDQM

First, we can consider the following one-dimensional hyperbolic telegraph equations:(26)∂2χx,t∂t2+k0∂χx,t∂t=a02∂2χx,t∂x2+b0χx,t,where k0>0, a0>0, and b0<0 are constants related to the transmission line parameters; variable χ(x,t) denotes voltage or current on transmission line; its domain of definition is (x,t)∈0,τ×0,T; with initial condition (27)χx,0=g1x,x∈0,τ;χtx,0=g2x,x∈0,τ.With Dirichlet boundary condition (28)χ0,t=m1t,t∈0,T;χτ,t=m2t,t∈0,T.By adopting transformationu~(x,t)=ek0t/2χ(x,t) in [13], the new one-dimensional telegraph equation is obtained by transformation; that is,(29)∂2u~x,t∂t2=a02∂2u~x,t∂x2+b0+14k02u~x,t.With new initial condition (30)u~x,0=g1x,x∈0,τ;u~tx,0=g2x+12k0g1x,x∈0,τ.With new Dirichlet boundary condition (31)u~0,t=ek0t/2m1t,t∈0,T;u~τ,t=ek0t/2m2t,t∈0,T.Using the PSM with Lagrange interpolation basis function and Chebyshev grid points (or uniform grid points) on the interval[-1,1] used as collocation points xi, the second-order partial derivatives u~xx can be expressed as follows:(32)u~xxx~i,t=∑j=0Nd~ij2u~x~j,t,i=0,1,…,N,where d~ij2=(2/τ)2dij2; x~i=τ/2xi+τ/2,i∈(0,N).Formula (32) is substituted into (29), and the following formula can be obtained: (33)d2U~dt2=2τ2a02DN2+b0+14k02IN+1U~,where U~=u~x~0,t,u~x~1,t,…,u~x~N-1,t,u~x~N,tT; IN+1 is the identity matrix of dimension N+1.LetV~=dU~/dt; formula (33) can be reduced to the first-order equations:(34)dV~dt=SU~,dU~dt=V~.If we defineY as follows:(35)Y=V~U~,obviously, formula (34) can be rewritten as the first-order linear ordinary differential equations: (36)dYdt=0SIN+10Y,where S=(2/τ)2a02DN2+(b0+(1/4)k02)IN+1.Using classic time domain DQM based on uniform grid points to solve first-order ordinary differential equation (36), the specific calculation process is given in [5] and there is no repeat here. So far, the basic idea of using PSDQM for solving one-dimensional telegraph equations is summarized as follows. First of all, by using the pseudospectral method to discretize the spatial variables and obtain a linear ordinary differential equations in time domain, then using thes-stages-order time domain DQM based on uniform grid points to solve the obtained first-order linear equations, and, finally, using the inverse transform χ^x~,t=e-k0t/2u~x~,t, the numerical solution of (26) on collocation points can be obtained.On the basis of solving the one-dimensional telegraph equations, PSDQM is extended to solve the two-dimensional hyperbolic telegraph equations [13]; that is,(37)∂2χ∂t2+k0∂χ∂t=a02∂2χ∂x2+∂2χ∂y2+b0χ,χx,y,0=g~1x,y,x,y∈Ω,χtx,y,0=g~2x,y,x,y∈Ω,χx,y,t=0,0≤t≤T,x,y∈∂Ω,where χ=χx,y,t; (x,y)∈Ω=0,τ×0,τ; ∂Ω is a closed curve of a enclosing region named Ω; the meaning of other similar variables is shown in (26).The transformationu^=ek0t/2χ in [13] is introduced; the new two-dimensional telegraph equations are obtained as follows:(38)∂2u^∂t2=a02∂2u^∂x2+∂2u^∂y2+b0+14k02u^,u^x,y,0=g~1x,y,x,y∈Ω,u^tx,y,0=g~2x,y+12k0g~1x,y,x,y∈Ω,u^x,y,t=0,0≤t≤T,x,y∈∂Ω.Using the Chebyshev pseudospectral method again, the second-order partial derivativesu^xx and u^yy are expressed as follows:(39)u^xxx^i,y^k,t=∑j=0N^d^ij2u^x^j,y^k,t,i=0,1,…,N^,(40)u^yyx^i,y^k,t=∑j=0M^d^kj2u^x^i,y^j,t,k=0,1,…,M^,where d^ij(2)=(2/τ)2dij(2); d^kj(2)=(2/τ)2dkj(2); x^i=τ/2xi+τ/2,i∈(0,N^); y^k=τ/2yk+τ/2,k∈(0,M^).On the basis of formulas (34)~(36) and (39)~(40), (38) can be rewritten as follows: (41)dzdt=IN^+1⊗Bz,where z=z0T,z2T,…,zN^-1T,zN^TT; zi=V^i,U^i; the symbol ⊗ denotes the tensor product; IN^+1 is the identity matrix of dimension (N^+1); U^i=u^x^i,y^1,t,u^x^i,y^2,t,…,u^x^i,y^M^,t, i∈(0,N^); V^i=u^˙x^i,y^1,t,u^˙x^i,y^2,t,…,u^˙x^i,y^M^,t, i∈(0,N^); B∈R2(M^+1)×2(M^+1), and its specific expressions can expressed as follows: (42)B=0S^IM^+10,where S^=(2/τ)2a02DM^2+(b0+1/4k02)IM^+1; IM^+1 is the identity matrix of dimension (M^+1).Similarly, classics-stages-order DQM based on uniform grid points is used to solve equation (41), and the solution procedure is exactly the same as the one-dimensional telegraph equation. Finally, the inverse transformation χ-x^,y^,t=e-k0t/2u^x^,y^,t is employed to obtain the numerical solution of telegraph equation (37).


## 5. Numerical Examples

In this section, Lagrange interpolation polynomial is used as the test function of PSDQM. Two collocation points, Chebyshev-Gauss-Lobatto grid points and uniform grid points, are selected for PSM. In order to illustrate the effectiveness of the PSDQM, the method used in [9] is compared with the PSDQM in this paper. And the time steps of both methods are h. Three numerical examples are given to test the effectiveness of PSDQM and compare the numerical results of PSDQM with the methods in [9].Example 1.
Solving one-dimensional telegraph equation in formula (26), let a0=2, b0=-3π2, k0=8π, and τ=1, and the initial conditions and boundary conditions are as follows:(43)g1x=sinπx;g2x=-πsin⁡πx;m1t=m2t=0.The exact solution of (26) is χ(x,t)=e-πtsin⁡πx.
In order to facilitate the analysis of the computational accuracy and convergence of the PSDQM and the method in [9], the relative errors defined in [14] are employed to evaluate the numerical results. The formula for relative error κ for one-dimensional telegraph equations is given as follows:(44)κ=∑i=0N∑j=0lχx~i,tj-χ^x~i,tj2∑i=0N∑j=0lχ2x~i,tj1/2.
Accordingly, the formula for relative errorκ~ for two-dimensional telegraph equations is as follows [14]: (45)κ~=∑i=0N^∑kM^∑j=0lχx^i,y^k,tj-χ-x^i,y^k,tj2∑i=0N^∑k=0M^∑j=0lχx^i,y^k,tj21/2.
The selected computed parameters areN1=10, h=0.01, and l=T/h. The results are shown in Table 1, which compares the relative errors of the PSDQM (s=4) in this article with the method in [9].
It can be seen that the calculation accuracy and numerical stability of PSDQM using Chebyshev-Gauss-Lobatto grid points are better than those of uniform grid points, and the relative error of the PSDQM (s=4) with Chebyshev points is nearly two orders of magnitude smaller than the method in [9].
To further verify the computational accuracy of PSDQM and the numerical stability in time direction, the number of collocation pointsN1 in space is varied and the time-step h=0.001 of DQM in the time domain remains unchanged. Figure 1 shows that the calculation accuracy of the PSDQM with Chebyshev grid points is used in spatial domain at different N1 when T=2, which can be seen that the relative error is convergent in form of spectral accuracy.
Table2 shows that, in time direction, PSDQM has unconditional stability, namely, A-stability [7, 15]. Therefore, PSDQM is better than UDDQM [16] in time domain.Table 1
Relative errors of Example1.
T
PSDQM
Method in [9]
Chebyshev points
Uniform points
0.5
4.2546
×
1
0
-
8
3.5798
×
1
0
-
8
1.0160
×
1
0
-
6
1.0
5.1883
×
1
0
-
8
5.2462
×
1
0
-
8
1.2390
×
1
0
-
6
2.0
5.3156
×
1
0
-
8
6.2491
×
1
0
-
7
1.2694
×
1
0
-
6
3.0
5.3165
×
1
0
-
8
1.1185
×
1
0
-
5
1.2696
×
1
0
-
6Table 2
Relative errors of Example1 with different N1 and T.
N
1
T
Chebyshev points
6
4
1.6729
×
1
0
-
5
6
6
1.6729
×
1
0
-
5
6
8
1.6729
×
1
0
-
5
6
10
1.6729
×
1
0
-
5
12
4
3.4613
×
1
0
-
12
12
6
3.4613
×
1
0
-
12
12
8
3.4613
×
1
0
-
12
12
10
3.4613
×
1
0
-
12Figure 1
Relative errors of Example1 with different N1 at T=2.Example 2.
Solving one-dimensional telegraph equation in formula (26), let a0=1/π, b0=-4, k0=4, and τ=2, and the initial conditions and boundary conditions are as follows:(46)g1x=cosh⁡πx;g2x=-coshπx;m1t=e-t;m2t=e-tcosh⁡2π.The analytic solution of (26) is χ(x,t)=e-tcoshπx.
Select calculation parameters:s=4, N2=18, h=0.1, and l=T/h. The calculation results at this case are shown in Tables 3 and 4. Obviously, Tables 3 and 4 further verify the conclusions given in Example 1.Table 3
Relative errors of Example2.
T
PSDQM
Method in [9]
Chebyshev points
Uniform points
0.5
8.2086
×
1
0
-
9
8.2208
×
1
0
-
9
1.9458
×
1
0
-
7
1.0
1.3742
×
1
0
-
8
1.3872
×
1
0
-
8
3.2409
×
1
0
-
7
2.0
2.4822
×
1
0
-
8
2.8261
×
1
0
-
8
4.6166
×
1
0
-
7
3.0
1.2167
×
1
0
-
8
1.3834
×
1
0
-
7
5.4131
×
1
0
-
7Table 4
Relative errors of Example2 with different N2 and T.
N
2
T
Chebyshev points
16
24
3.6662
×
1
0
-
5
16
26
3.6662
×
1
0
-
5
16
28
3.6662
×
1
0
-
5
16
30
3.6662
×
1
0
-
5
18
24
8.5851
×
1
0
-
6
18
26
8.5851
×
1
0
-
6
18
28
8.5851
×
1
0
-
6
18
30
8.5851
×
1
0
-
6Example 3.
Solving the two-dimensional telegraph equation in formula (37), let a0=2, b0=-3π2, k0=8π, and τ=1. The initial conditions are as follows:(47)g~1x,y=sin⁡πxsin⁡πy;g~2x,y=-πsin⁡πxsin⁡πy.In this case, the analytic solution of (37) is χx,y,t=e-πtsinπxsin⁡πy.
Select the computational parameters:N3=M3=8, h=0.01, and l=T/h. The corresponding calculation results are shown in Table 5. As can be seen from Table 5, in the two-dimensional case, the PSDQM (s=4) using Chebyshev points still has higher computational accuracy than the methods in [9].
In order to compare the computational accuracy of PSDQM with differents, N3=M3=12, the time-step h=0.1 of DQM in time domain remains unchanged. At this situation, the results are shown in Figure 2, which shows that, with the increase of s in DQM, because of the spectral accuracy in space, the accuracy of PSDQM increases with the increase of grid points’ number in time domain. And this also shows the advantage of pseudospectral method.
Table6 shows that PSDQM is unconditionally stable in the time direction and in the two-dimensional case.Table 5
Relative errors of Example3.
T
PSDQM
Method in [9]
Chebyshev points
Uniform points
0.5
8.0566
×
1
0
-
8
1.3426
×
1
0
-
6
1.0031
×
1
0
-
6
1.0
8.2824
×
1
0
-
8
1.5568
×
1
0
-
6
1.2218
×
1
0
-
6
2.0
8.3131
×
1
0
-
8
1.5885
×
1
0
-
6
1.2516
×
1
0
-
6
3.0
8.3133
×
1
0
-
8
1.5887
×
1
0
-
6
1.2518
×
1
0
-
6Table 6
Relative errors of Example3 with different N3 and T.
N
3
T
Chebyshev points
12
23
1.3953
×
1
0
-
10
12
26
1.3953
×
1
0
-
10
12
29
1.3953
×
1
0
-
10
12
31
1.3953
×
1
0
-
10
18
23
1.4380
×
1
0
-
10
18
26
1.4380
×
1
0
-
10
18
29
1.4380
×
1
0
-
10
18
31
1.4380
×
1
0
-
10Figure 2
Relative errors of Example3 with different s at T=3.


## 6. Conclusion

Based on the comparison and analysis to the basic definition, interpolation points, and the test function, differentiation matrix or weighting coefficient matrix of pseudospectral method and differential quadrature method, a coupled numerical method, namely, pseudospectral-differential quadrature method (PSDQM), is proposed for solving the multidimensional telegraph equations. And this paper pointed out the equivalence of the two meshfree methods in the space domain. By analyzing three specific numerical examples, the numerical results show that PSDQM has spectral accuracy convergence in space domain and high accuracy and A-stability in time domain. In general, the coupled method inherits the advantages of the pseudospectral method and differential quadrature method, and the two methods can be mutually beneficial. In all, the PSDQM is suitable for solving multidimensional telegraph equations.



---

*Source: 9013826-2017-12-20.xml*