Zero stability


Zero-stability, also known as D-stability in honor of Germund Dahlquist,[1] refers to the stability of a numerical scheme applied to the simple initial value problem .

A linear multistep method is zero-stable if all roots of the characteristic equation that arises on applying the method to have magnitude less than or equal to unity, and that all roots with unit magnitude are simple.[2] This is called the root condition[3] and means that the parasitic solutions of the recurrence relation will not grow exponentially.

Example

The following third-order method has the highest order possible for any explicit two-step method[2] for solving : . If identically, this gives a linear recurrence relation with characteristic equation . The roots of this equation are and and so the general solution to the recurrence relation is . Rounding errors in the computation of would mean a nonzero (though small) value of so that eventually the parasitic solution would dominate. Therefore, this method is not zero-stable.

References

  1. Dahlquist, Germund (1956). "Convergence and stability in the numerical integration of ordinary differential equations". Mathematica Scandinavica. 4 (4): 33–53. doi:10.7146/math.scand.a-10454. JSTOR 24490010. Retrieved 19 July 2022.
  2. Hairer, Ernst; Nørsett, Syvert; Wanner, Gerhard (1987). Solving Ordinary Differential Equations I. Berlin: Springer-Verlag. pp. 326–328.
  3. Butcher, John C (1987). The Numerical Analysis of Ordinary Differential Equations. Wiley. p. 11.
This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.