Haar wavelet approach to linear stiff systems☆
Introduction
When the solution of a system contains components which change at significantly different rates for given changes in the independence variable the system is said to be “stiff” [5], [12]. In practical solution of stiff problems, the choice of the solution steps is critical. Large steps will lose some fast changing properties of the system, while small steps will introduce too many round-off errors and cause numerical instability.
In [2], Carroll presents an exponentially fitted scheme for solving the stiff systems of initial-value problems. A long-standing work by Pavlov and Rodionova [13], [14] retains nonlinear terms of the original equations. However, these methods cannot reduce computational effort. In this paper, we solve the linear stiff systems via Haar wavelets by taking advantage of some nice properties of Haar wavelets.
Section snippets
Preliminary works
The orthogonal set of Haar wavelets hi(t) is a group of square waves with magnitude of ±1 in some intervals and zeros elsewhere [3], [4], [6], [7], [8], [9], [10], [11]. In general, Namely, each Haar wavelet contains one and just one square wave, and is zero slsewhere. Just these zeros make Haar wavelets to be local and very useful in solving stiff systems.
Any function y(t), which is square integrable in the interval [0,1], can
Solution of linear stiff system
Consider a linear stiff system of the following form: In the normalized interval with τ=m·t Now expressing , , and in single-term Haar series as and the following recursive relationship is obtained with P=1/2, where . Substituting , into Eq. (13), we obtain nonlinear algebraic equations and they may be
Numerical examples
Consider linear stiff systems in the following examples. Example 1 This example is due to Alt [1] The exact solution is
Example 2 We take the stiff system of Rosenbrock and Storey [17] The exact solution is
Example 3 Consider the following linear system [12]
Additional comparisons
We would like to point out that this is a common property for many transform approaches to linear systems including Haar, Walsh, block pulse, Fourier and many others; and all result in a generalized Lyapunov equation. The kernel of the complete solution lies on how to solve this generalized Lyapunov equation. Compare the proposed STHW approach with the known approaches as follows.
References (17)
A-stable one-step methods with step-size control for stiff systems of ordinary differential equations
J. Comput. Appl. Math.
(1978)A matricial exponentially fitted scheme for the numerical solution of stiff initial-value problems
Comput. Math. Appl.
(1993)- et al.
Haar wavelet method for solving lumped and distributed-parameter systems
IEE Proc. Pt. D Control Theory Appl.
(1997) State analysis of linear time delayed systems via Haar wavelets
Math. Comput. Simul.
(1997)- et al.
State analysis of time-varying singular nonlinear systems via Haar wavelets
Math. Comput. Simul.
(1999) - et al.
State analysis of time-varying singular bilinear systems via Haar wavelets
Math. Comput. Simul.
(2000) - et al.
Wavelet approach to optimising dynamic systems
IEE Proc. Pt. D Control Theory Appl.
(1999) - E. Hairer, G. Wanner, Solving ordinary differential equations II, stiff and differential-algebraic problems,...
Cited by (0)
- ☆
To be submitted to Mathematics and Computers in Simulation. This research was supported in part by the National Science Council of Taiwan under Grant NSC.