Multiresolution Time-Domain Scheme for Terminal Response of Two-Conductor Transmission Lines

This paper derives a multiresolution time-domain (MRTD) scheme for the two-conductor lossless transmission line equations based on Daubechies’ scaling functions. And a method is proposed to generate the scheme at the terminal and near the terminal of the lines. The stability and numerical dispersion of this scheme are studied, and the proposed scheme shows a better dispersion property than the conventional FDTD method. Then the MRTD scheme is extended to the two-conductor lossy transmission line equations. The MRTD scheme is implemented with different basis functions for both lossless and lossy transmission lines. Numerical results show that theMRTD schemes which use the scaling functions with high vanishingmoment obtainmore accurate results.


Introduction
The multiresolution time-domain (MRTD) scheme proposed by [1] provides an efficient algorithm for electromagnetic field computation and shows excellent capability to approximate exact solution with low sampling rates.However, the Battle-Lemarie wavelet function used in [1] is not compact supported, which means the iterative equations contain infinite terms.We must cut off the iterative equations in the actual computation and this may introduce truncation errors.So different wavelet bases, which are compact supported with some numbers of vanishing moments, have been used to improve this method [2][3][4][5].This makes a great development for MRTD schemes.As a kind of numerical method, the MRTD schemes show great advantages in numerical dispersion properties [6][7][8][9]; meanwhile, these schemes need a more rigorous stable condition than the conventional FDTD method [10].For containing more terms in the iterative equations, the terminal conditions or absorbing boundary conditions are more complicated to process in MRTD schemes; this disadvantage has limited the application of the MRTD scheme.To overcome this limitation, some works on the perfect match layer have been made [11][12][13]; however, other terminal conditions also need to be analyzed specifically.For the transmission lines equations, the resistive terminal conditions could be equivalent as a general Thevenin circuit; this paper will solve this kind of terminal condition in the MRTD scheme.
Since the appearance of the telegraph equations, studies on transmission lines have had a considerable development.Several equivalent forms of transmission line theory have been proposed to describe the influence of the incident electromagnetic field to the transmission lines [14][15][16].In [17], the classical theory of the transmission line has been summarized and the theory on the high frequency radiation effects to the transmission lines is introduced.In the monograph [18], the multiconductor transmission lines (MTL) theory has been comprehensively studied in detail.For the twoconductor lossless transmission lines, there are several methods, which contain the series solution, the SPICE solution, the time-domain to frequency-domain (TDFD) transformation method, and the FDTD method [18].However, the MRTD scheme has not been used to calculate the terminal response of transmission lines.In this paper we will derive a MRTD scheme for this problem.
In this paper, we focus on the calculation of the terminal response of two-conductor transmission lines equations by using MRTD scheme.In Section 2, the MRTD scheme is derived based on Daubechies' scaling functions for the two-conductor lossless transmission line equations, and, for the resistive terminations, the iterative equations for the terminal voltages are derived, a method is proposed to update the iterative equations which contain some terms whose indices exceed the index range in the MRTD scheme, and then the stability and the numerical dispersion are studied.In Section 3, the MRTD scheme is extended to the two-conductor lossy transmission line.In Section 4, the numerical results are presented on the terminal response of both lossless and lossy transmission lines using the MRTD scheme and compared to the FDTD method at different space discretization numbers and different Courant numbers.

MRTD Scheme for Two-Conductor Lossless Transmission Lines
2.1.MRTD Formulation.In this section, the MRTD scheme is applied to the following scalar transmission lines equations for two-conductor lossless lines [18]: where  and  are the per-unit-length inductance and capacitance, respectively.Based on the method outlined in [1], the voltage and current can be expanded as follows: where    and  +1/2 +1/2 are the coefficients for the voltages and currents in terms of scaling functions, respectively.The indices  and  are the discrete spatial and temporal indices related to space and time coordinates via  = Δ and  = Δ, where Δ and Δ represent the spatial and temporal discretization intervals in  and  direction.The function ℎ  () is defined as with the rectangular pulse function The function   () is defined as where () represents Daubechies' scaling function.Figure 1 shows Daubechies' scaling function with two vanishing moments.
For deriving the MRTD scheme for (1a) and (1b), we need the following integrals: where  ,  represents the Kronecker symbol.Consider where   denotes the effective support size of the basis functions.The coefficients () are called connection coefficients and can be calculated by (9).Taking Daubechies' scaling functions as the basis functions, Table 1 shows () for 0 ≤  ≤   − 1, which are zeros for  >   − 1, and for  < 0 it can be obtained by the symmetry relation where φ() represents the Fourier transform of ().
Following the theory in [3] and making use of ( 10), ( 5) is modified to In spite of the support of the scaling functions [20], singlepoint sampling of the total voltages and currents can be taken at integer points with negligible error.Taking voltage at spatial point Δ and at time Δ, we obtain where  is the Dirac delta function.Equation ( 12) means the voltage value at each integer point is equal to the coefficient.The current values have the same character at each half integer point.Therefore, we will use    and  +1/2 +1/2 directly to represent the voltage at the point (Δ, Δ) and the current at the point (( + 1/2)Δ, ( + 1/2)Δ) in this paper.
The modified   () in (11) also satisfy integrals (7) and (8).Applying the Galerkin technique to (1a) and (1b), we can obtain the following iterative equations for the voltages and currents:

Terminal Iterative Equations for Resistive Load in MRTD
Scheme.We will consider the terminal conditions for the two-conductor lossless transmission lines equations in this section.Equations (1a) and (1b) are homogeneous linear equations; we need to add the terminal conditions to obtain the unique solution.
Considering the two-conductor lines shown in Figure 2, we assume the length of the total line is  and the resistive loads are   and   .The line is divided uniformly into NDZ segments with the space interval Δ and the total solution time is divided into NDT steps with the uniform time interval Δ.Similar to the conventional FDTD, we will calculate the interlace voltages,   0 ,   1 , . . .,   NDZ , and currents,  +1/2 1/2 ,  +1/2 3/2 , . . .,  +1/2 NDZ−1/2 , in both space domain For the resistive terminations, we note the voltage at the source ( = 0) as   () and the current at the source as   (), the external voltage at the load ( = ) as   (), and the current at the load as   ().The discrete voltages and currents at the source are denoted as    ≡   (Δ) and    ≡   (Δ), and the discrete voltages and currents at the load are denoted as    ≡   (Δ), and    ≡   (Δ), then the terminal characterizations could be written in terms of a generalized Thevenin equivalent as Equations ( 14a) and (14b) denote the discretization terminal conditions for the case of resistive terminations, so we need to introduce these conditions to the iterative equations (13a) and (13b) to obtain the numerical solution.
Notice that, in the iterative equations (13a) and (13b), not only the iterative equations of the terminal voltages  +1 0 and  +1 NDZ should be derived and the iterative equations of voltages and currents "near" the terminals also need to be updated.The voltages and currents "near" the terminals we mean are the voltages  +1  and   NDZ− for  = 1, 2, . . .,   − 1,and the currents  +1/2 +1/2 and  +1/2 NDZ−+1/2 for  = 0, 1, . . .,   − 2. All of these voltages and currents contain some terms that exceed the index range in iterative equations (13a) and (13b).Figure 3 shows the discretization of the terminal voltages and the voltages and currents near the terminal.
We will derive the MRTD scheme at the terminal firstly.For updating the iterative equations for the terminal voltages, we need to decompose iterative equations ((13a), (13b)).Since the coefficients () satisfy the following relation [4] substituting ( 15) into (13a), we can obtain Considering the corresponding terms with , we can decompose (13a) as [21] for  = 0, 1, . . .,   − 1. Equation (13b) could make the analogous decomposition.We could view the MRTD scheme for two-conductor transmission lines as the weighted mean of the conventional FDTD method with spatial discretization step (2 + 1)Δ for  = 0, 1, . . .,   − 1, and the weighting coefficient for each term is (2+1)().Besides, for the MRTD scheme whose coefficients () satisfy relationship (15), the analogous decomposition could be made.This relationship between the MRTD scheme and the conventional FDTD method is useful for us to update the iterative equations.Taking  +1 0 as an example to derive the iterative equations at the terminal, Following steps of ( 16) and ( 17), we can decompose (18) as Here, we could view each equation in (19a), (19b), and (19c) as a central difference scheme, (19a) is the central difference scheme related to points  = −Δ/2 and  = Δ/2, (19b) is the central difference scheme related to points  = −3Δ/2 and  = 3Δ/2, and so on, but the terms , whose subscripts exceed the index range, make (19a), (19b), and (19c) out of work.So we need to make some update for the iterative equations.Using the forward difference scheme to replace the central difference scheme, we change the difference points in (19a) to become  = 0 and  = Δ/2 and change the difference points in (19b) to become  = 0 and  = 3Δ/2; the others follow the same step.Keeping the weighting coefficient unchanged in each equation, we can obtain . . .
where the terminal current  +1/2  = (   +  +1  )/2 and    can be derived from (14a) Summing up all the equations in (20a), (20b), and (20c), Substituting ( 21) into (22), we can obtain the iterative equation at the source With the same steps, we can obtain the iterative equation at the load After deriving the iterative equations at the terminal, we will put forward a truncation method to update iterative equations which contain some terms whose indices exceed the index range in the MRTD scheme.
Noticing the first  terms in (25a), (25b), (25c), (25d), and (25e), there is no term exceeding the index range in each equation.Meanwhile, the equations which contain the exceeding indices terms are all appearing in the rest of   −  terms.As mentioned before, we can view the MRTD scheme as the weighted mean of the conventional FDTD method, but (25a), (25b), (25c), (25d), and (25e) show that the last   −  equations are unavailable for forming the iterative equations in MRTD scheme.To solve this problem, we make truncation here.We update the iterative equation of  +1  by using the summation of the first  terms in (25a), (25b), (25c), (25d), and (25e) that means we use the weighted mean of the first  to approximate the summation of all the   terms.
For the iterative equation of the current, there is a little difference from the voltage's.As shown in Figure 3, the interlace currents appear at the half integer points which means all the currents are located at the interior points.So we only need to modify the currents near the terminals.Following the same steps of the derivation of voltages iterative equations near the terminal, we could obtain the current iterative equations near the terminals.

Terminal Iterative Equations for Inductive
Resistance in MRTD Scheme.Since we have discussed the resistive load in Section 2.2, we will consider a more complicated terminal load which consisted of a resistance and inductance shown in Figure 4.The resistance and the inductance at the load are noted as   and   , respectively.Keeping the source as a resistive terminal and changing the load including inductance, we note the voltage at the source ( = 0) as   () and the current at the source as   (), the external voltage at the load ( = ) as   (), and the current at the load as   ().The terminal conditions could be written as follows: Expanding  0 (),  0 (),   (), and   () as (2a) and (2b) and sampling them at the time discreting point ( + 1/2)Δ, the terminal conditions become It could be seen from Section 2.2 that the change of the terminal condition only affects the iterative equations at the terminals, so we just need to consider two iterative equations.Since we keep the source as a resistive termination, the iterative equation at the source should be the same as (23).Actually, if we substitute (33a) into (22), we could obtain (23).So the only iterative equation we should derive is located at the load.If we set   = 0 in (33b), the terminal condition at the load will have the analogous form with the source, which means the load with a resistance and inductance degenerates to be a resistance.For those   ̸ = 0, transforming (33b) as ) . ( Following the steps we get the iterative equation of  +1 0 in Section 2.2; we can obtain the iterative equation at the load where for  = 2, 3, . . ., NDT.

Stability Analysis.
For the purpose of stability analysis, (13a) and (13b) can be rewritten as Following the procedures in [22], the finite-difference approximations of the time derivations on the left side of (37a) and (37b) can be written as an eigenvalue problem In order to avoid instability during normal time stepping, the imaginary part of  must satisfy As we consider the lossless two-conductor transmission lines in (1a) and (1b), the transient values of voltages and currents distributed in space can be Fourier transformed with respect to -coordinates to provide a spectrum of sinusoidal modes.Assuming an eigenmode of the spectral domain with   , the voltages and currents can be written as where   represents the voltage at point Δ,  +1/2 represents the current at point ( + 1/2)Δ, and   0 and   0 are the amplitudes of the voltages and currents.Substituting (40a) and (40b) into (38a) and (38b) and (37a) and (37b), we obtain In (41),  is pure imaginary and where V = 1/ √  is the velocity of the wave along with the lines.Numerical stability is maintained for every spatial mode only when the range of eigenvalues given by (42) is contained entirely within the stable range of time-difference eigenvalues given by (39).Since both ranges are symmetrical around zero, it is adequate to set the upper bound of (42) to be smaller or equal to (39); we can obtain the stability condition Noting the Courant number the maximum values of  required by a stable algorithm, which are listed in Table 1 as  max , can be calculated from the connection coefficients.

Dispersion Analysis.
To calculate the numerical dispersion, substituting a time harmonic trial solution into (37a) and (37b), we can obtain where  is the wave angular frequency and   is the numerical wave number.
Using the number of cells per wavelength   =  REAL /Δ and the wave number   = (2)/ NUM , we obtain the dispersion relationship where  =  REAL / NUM is the ratio between the theoretical and numerical wavelength.The value of  could be computed by Newton iterative method and using the formula |( − 1) × 360| calculates the phase error in degrees.Figure 5 shows the phase errors of different MRTD schemes versus the samples per wavelength.Compared with the conventional FDTD method at the same sampling numbers, the MRTD schemes show better numerical dispersion properties than FDTD method, which means the discretization error of the MRTD schemes is smaller than the FDTD method.

MRTD Scheme for Two-Conductor Lossy Transmission Lines
In this section, we will extend the MRTD scheme to the twoconductor lossy transmission lines.

MRTD
where , , , and  are the per-unit-length resistance, inductance, conductance, and capacitance, respectively.Compared to the lossless transmission lines, the lossy case must consider the resistance losses along the lines and the losses in the medium.However, the steps to obtain the MRTD scheme are the same.Similar to the lossless case, we firstly extend the voltage and current with Daubechies' scaling functions and the rectangular function and use    and  +1/2 +1/2 representing the voltage at the point (Δ, Δ) and the current at the point ((+1/2)Δ, (+1/2)Δ), respectively, and Δ and Δ represent the spatial and temporal discretization intervals.By applying the Galerkin technique, we can obtain the iterative equations for the lossy transmission lines equations where () is the connection coefficient and  1 ,  2 ,  1 , and  2 are constants The differences between (48a) and (48b) and (13a) and (13b) are the coefficients of terms in the iterative equations that are caused by the unit-per-length resistance  and inductance .If we set  = 0 and  = 0, the coefficients  1 ,  2 and  1 ,  2 are all equal to 1, and the iterative equations (48a) and (48b) will degenerate to iterative equations (13a) and (13b).
For the lossy case, we should also modify the iterative equations at the terminal and near the terminal of the lines.Considering the two-conductor lines shown in Figure 2, we assume the length of the total line is  and the resistive loads are   and   .The line is divided uniformly into NDZ segments with the space interval Δ and the total solution time is divided into NDT steps with the uniform time interval Δ.Following the same steps in Section 2.2, we could obtain the modified iterative equations.
For the voltage iterative equation at the source For the voltage iterative equation at the load For the voltage iterative equations near the source for  = 1, 2, . . .,   − 1.

Stability Analysis.
To study the stability of the MRTD scheme for lossy case, we need to make some changes.Rewrite (48a) and (48b) as follows: The finite-difference approximations of the time derivations on the left side of (61a) and (61b) are different from the left side of (37a) and (37b), but we can also write them as an eigenvalue problem In order to avoid instability during normal time stepping, the imaginary part of  also must satisfy Then following the steps in Section 2.4, using the Fourier transform, expand the transient values of voltages and currents distributed in space.We can also obtain the same stability condition as (43) However, the dispersion analysis for the iterative equations of the lossy transmission lines is quite different from the lossless case.Since +1/2 on the left side of (61a) and (61b), if we substitute a time harmonic into the iterative equations into (61a) and (61b), there are some other terms that could affect the ratio between the theoretical and numerical wavelength except samples per wavelength, so we could not obtain a brief dispersion relation as (45).However, the numerical results in Section 4.2 show that the MRTD method could obtain a more accurate result than the FDTD method at the same space interval and time interval.

Numerical Result
4.1.Lossless Transmission Lines.In this section, the twoconductor lossless transmission lines as shown in Figure 2 are considered to calculate the terminal voltages.The length of the lines is  = 400 m and the per-unit-length capacitance and inductance are  = 100 pF/m and  = 0.25 H/m, respectively.This corresponds to RG58U coaxial cable [18].The characteristic impedance of the line is   = √/ = 50 Ω and the velocity of propagation is V = 1/ √  = 200 m/s.The terminal load is 100 Ω and the source resistance is 50 Ω.
We will use a trapezoidal pulse as shown in Figure 6 as the source voltage, whose initial value is 0 V and amplitude is  = 30 V, total time is 10 s, rise time is   = 1 s, and fall time is   = 1 s, respectively.The bandwidth of the pulse is approximate to BW = 1/  , so the segment length Δ should be electrically short at this frequency requiring Δ ≤ (1/10)(V/ max ) ( max = 1/  ) [18].As we take the rise time of the trapezoidal pulse to be 1 s, the approximate bandwidth of the pulse is 1 MHz and the maximum section length of Δ should be less than 20 m, noting Δ max = 20.Since we divide the line into NDZ segments uniformly, NDZ should be greater than /Δ max = 20.The time step will be calculated by (44) with different Courant number .
Under the conditions of the space discretization number NDZ = 20 and the Courant number  = 0.5, we calculate the terminal voltages using the MRTD schemes with  2 ,  3 , and  4 wavelets' scaling functions.We also calculate the terminal voltages by the FDTD method under the same conditions.Figure 7 shows the numerical results by different methods and   -MRTD represents the MRTD scheme using Daubechies' scaling functions with  vanishing moment as basis functions, where  = 2, 3, 4.
For the time-dependent discrete terminal voltages, the relative error is defined as follows [23]: Here, () represents the numerical results of FDTD and MRTD at each time discretization point and x() is the numerical result of series solution at each time discretization point, which can be regarded as the exact solution [18].Taking x() as the exact result, we calculate the relative errors of FDTD method and MRTD scheme.Table 2 shows the relative errors and runtime for different schemes.Since the iterative equations in MRTD schemes contain more terms than the conventional FDTD method, the MRTD schemes expend more runtime.When we use  2 -MRTD, the numerical result shows a larger relative error.The reason is that the vanishing moment of the  2 wavelet's scaling function is not high enough.For the wavelets' scaling functions whose vanishing moment is high enough, like  3 wavelet and  4 wavelet, the numerical results show smaller relative errors.
Figure 8 shows the relative errors versus the space discretization numbers.The space interval will decrease with the increase of the space discretization number.We can see from the figures that the relative errors for  2 -MRTD increase with the increase of NDZ.And, for  3 -MRTD and  4 -MRTD, the relative errors decrease with the increase of NDZ and show a little smaller relative error than the conventional FDTD.It means we can get a more accurate result at same space interval and time interval by  3 -MRTD and  4 -MRTD.
Figure 9 describes the relative errors versus the Courant numbers.The time interval will decrease with the decrease of the Courant number.The results show that, with the decrease of the Courant numbers, the relative errors are almost unchanged for the conventional FDTD method, while the relative errors are quite different with different Courant numbers for all the three MRTD schemes.That means the choice of the Courant number significantly affects the relative errors of the MRTD schemes, and we could choose the proper Courant number to optimize the MRTD schemes.
To validate the stability of the MRTD schemes, we increase the initial value of the pulse to 20 V and keep other parameters unchanged.Figure 10 shows the relative errors versus the pace discretization numbers.It can be seen from the figures that the  2 -MRTD scheme shows a larger relative error and the relative errors of  3 -MRTD and  4 -MRTD schemes are smaller than the conventional FDTD method.
The numerical results for the lossless transmission lines also show that the  2 -MRTD does not perform better than the FDTD method; meanwhile,  3 -MRTD and  4 -MRTD schemes show better quality in accuracy and stability.The reason for this phenomenon is that the scaling function of the  2 wavelet does not have enough high vanishing moment.When we use Daubechies' scaling functions to expand the voltages and currents in the two-conductor transmission lines equations, the vanishing moment decides the accuracy of the approximation.The scaling functions with high vanishing moments could approximate voltages and currents more accurately; however, the scaling functions with low  vanishing moments like  2 wavelet's may introduce a larger error.So we can see from Figure 5 that  2 -MRTD has a better numerical dispersion property than the FDTD method, but it gets larger relative errors in numerical computation.For the scaling functions which have high enough vanishing moments like  3 wavelet and  4 wavelet, the numerical results show smaller relative errors and are in agreement with the dispersion analysis.However, a high vanishing moment for scaling function may increase the computation complexity in MRTD schemes.So the vanishing moments of the wavelet's scaling functions have a great effect on the accuracy of the MRTD scheme; it is necessary to choose a scaling function with proper vanishing moment when we use the MRTD scheme for the numerical computation.

Lossy Transmission Lines.
In this section, we will consider the lossy two-conductor transmission lines shown in Figure 11.Two conductors of rectangular cross section of width  = 20 m and thickness  = 10 m are separated by  = 20 m and placed on one side of a silicon substrate (  = 12) of thickness ℎ = 100 m; the total line length is  = 20cm.The near end is a source with a   = 50 Ω resistance and the far end is a load with a   = 50 Ω resistance and   = 0.5 H inductance in series.The perunit-length inductance and capacitance were computed as  = 0.805969 H/m and  = 88.2488pF/m.This gives a velocity of V = 1.18573 × 10 8 m/s and a one-way time delay of   = 1.68672 ns, which gives an effective dielectric constant of (   = 6.4) and characteristic impedance of   = 95.566Ω.The per-unit-length dc resistance is computed as  = 1/() = 86.207Ω/m [18].Dielectric loss is not included in these calculations, which means the per-unitlength conductance is  = 0.
The source is a ramp function as shown in Figure 12, the initial value of the ramp function is 0 V, and the amplitude is  = 1 V with a rise time of   = 0.5 ns.The total computing time is 20 ns.The bandwidth of the source is approximate to BW = 1/  = 2 GHz.The space discretization step for the MRTD was chosen to be /10, so the space discretization  We calculate the near end voltage and the far end voltage of the lossy PCB with space discretization number NDZ = 40 and the Courant number  = 0.5.Figure 13 shows the computing results.Since  2 -MRTD may introduce a larger error in the computation as shown in Section 4.1, we use  3 -MRTD and  4 -MRTD to compute the terminal voltages.
Here, we choose the time-domain to frequency-domain transformation method (TDFD), which is a straightforward adaptation of a common analysis technique for lumped, linear circuits and systems [18], to validate the computing results of MRTD schemes and FDTD method.Table 3 shows the relative errors and the runtime of the MRTD schemes and FDTD method.It can be seen that MRTD schemes spend much time to obtain more accurate results.
Figure 14 describes the relative errors versus the space discretization numbers.For both  3 -MRTD and  4 -MRTD, the MRTD schemes show a little smaller relative error than the conventional FDTD.That means the MRTD schemes could obtain a more accurate solution under the same time interval and space interval.This is because  3 -MRTD and  4 -MRTD have better dispersion property than the FDTD method and the scaling functions of  3 wavelet and  4 wavelet have high enough vanishing moments.
Figure 15 shows the relative errors versus the Courant numbers.With the decrease of the Courant number, the relative errors of the FDTD method increase, because the time-domain is oversampled for the FDTD method.However, the MRTD schemes perform decreasing relative errors.That means the MRTD schemes could obtain a more accurate result with a smaller time interval.And even when the FDTD method is oversampled in the time-domain, the MRTD schemes perform well.

Conclusion
In this paper, we derived the MRTD scheme for the twoconductor transmission lines and studied the stability and the numerical dispersion of this scheme.By viewing the MRTD schemes as the weighted mean of the conventional FDTD method at different space interval, we derived the iterative equations for the terminal voltages when the terminals are pure resistive, and a method is proposed to update the iterative equations which contain some terms whose indices exceed the index range in the MRTD scheme.Using the same method, we derived the iterative equation for the inductive load.Then we extended the MRTD scheme to the lossy transmission lines.Using different wavelets' scaling functions as basis functions, the MRTD schemes are implemented for both lossless case and lossy case and the numerical results are compared to the conventional FDTD method.The numerical results show the MRTD schemes need more runtime to obtain more accurate results.And the vanishing moment of the wavelet's scaling functions will significantly affect the quality of the MRTD scheme; using a scaling function with a proper vanishing moment as basis function in MRTD scheme could obtain a more accurate result.

Figure 2 :
Figure 2: A two-conductor line in time-domain.

Figure 4 :
Figure 4: A two-conductor line with an inductive resistance.

Figure 5 :
Figure 5: Phase error (in degrees) for different MRTD schemes.The Courant numbers are  = 0.5 for (a) and  = 0.1 for (b).

Figure 6 :
Figure 6: Trapezoidal pulse with the rise time   and fall time   .

Figure 7 :
Figure 7: Terminal voltage calculated by different numerical methods.The space discretization number and the Courant number are NDZ = 20 and  = 0.5.

Figure 8 :Figure 9 :
Figure 8: Relative error for different discretization numbers.The Courant numbers are  = 0.5 for (a) and  = 0.1 for (b).

Figure 10 :
Figure 10: Relative error for different space discretization numbers for the initial value of the trapezoidal pulse to be 20 V. The Courant numbers are  = 0.5 for (a) and  = 0.1 for (b).

Figure 11 :
Figure 11: A lossy printed circuit board: (a) line dimensions and terminations and (b) cross-sectional dimensions.

Figure 12 :
Figure 12: Representation of the source voltage waveform.

Figure 13 :
Figure 13: Terminal voltages for lossy transmission lines: (a) for the near end and (b) for the far end.The space discretization number is NDZ = 40 and the Courant number is  = 0.5.

Figure 14 :Figure 15 :
Figure 14: Relative error for the terminal voltages with different space discretization numbers: (a) for the near end and (b) for the far end.The Courant number is  = 0.5.

Table 1 :
Connection coefficients () and the first-order moments  1 of Daubechies' scaling functions.

Table 3 :
Relative errors and runtime of different methods. 0 represents the near end voltage and   represents the far end voltage (NDZ = 40,  = 0.5).