It is wellknown that the stability analysis of step-by-step numerical methods for differential equations often reduces to the analysis of linear difference equations with variable coefficients. This class of difference equations leads to a family F of matrices depending on some parameters and the behaviour of the solutions depends on the convergence properties of the products of the matrices of F. To date, the techniques mainly used in the literature are confined to the search for a suitable norm and for conditions on the parameters such that the matrices of F are contractive in that norm. In general, the resulting conditions are more restrictive than necessary. An alternative and more effective approach is based on the concept of joint spectral radius of the family F, r(F). It is known that all the products of matrices of F asymptotically vanish if and only if r (F) < 1. The aim of this chapter is that to discuss the main theoretical and computational aspects involved in the analysis of the joint spectral radius and in applying this tool to the stability analysis of the discretizations of differential equations as well as to other stability problems. In particular, in the last section, we present some recent heuristic techniques for the search of optimal products in finite families, which constitute a fundamental step in the algorithms which we discuss. The material we present in the final section is part of an original research which is in progress and is still unpublished.

Stability of linear problems: joint spectral radius of sets of matrices

ZENNARO, MARINO
2014

Abstract

It is wellknown that the stability analysis of step-by-step numerical methods for differential equations often reduces to the analysis of linear difference equations with variable coefficients. This class of difference equations leads to a family F of matrices depending on some parameters and the behaviour of the solutions depends on the convergence properties of the products of the matrices of F. To date, the techniques mainly used in the literature are confined to the search for a suitable norm and for conditions on the parameters such that the matrices of F are contractive in that norm. In general, the resulting conditions are more restrictive than necessary. An alternative and more effective approach is based on the concept of joint spectral radius of the family F, r(F). It is known that all the products of matrices of F asymptotically vanish if and only if r (F) < 1. The aim of this chapter is that to discuss the main theoretical and computational aspects involved in the analysis of the joint spectral radius and in applying this tool to the stability analysis of the discretizations of differential equations as well as to other stability problems. In particular, in the last section, we present some recent heuristic techniques for the search of optimal products in finite families, which constitute a fundamental step in the algorithms which we discuss. The material we present in the final section is part of an original research which is in progress and is still unpublished.
File in questo prodotto:
File Dimensione Formato  
Lect. Notes Math. 2082(2014), 265-313.pdf

non disponibili

Descrizione: Articolo principale
Tipologia: Documento in Versione Editoriale
Licenza: Digital Rights Management non definito
Dimensione 1.62 MB
Formato Adobe PDF
1.62 MB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/11368/2707485
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 10
  • ???jsp.display-item.citation.isi??? 7
social impact