Please use this identifier to cite or link to this item: http://dx.doi.org/10.14279/depositonce-14675
For citation please use:
Main Title: Computation of Stability Radii for Large-Scale Dissipative Hamiltonian Systems
Author(s): Aliyev, Nicat
Mehrmann, Volker
Mengi, Emre
Type: Research Paper
URI: https://depositonce.tu-berlin.de/handle/11303/15902
http://dx.doi.org/10.14279/depositonce-14675
License: http://rightsstatements.org/vocab/InC/1.0/
Abstract: A linear time-invariant dissipative Hamiltonian (DH) system x = (J-R)Q x, with a skew-Hermitian J, an Hermitian positive semi-definite R, and an Hermitian positive definite Q, is always Lyapunov stable and under weak further conditions even asymptotically stable. In various applications there is uncertainty on the system matrices J, R, Q, and it is desirable to know whether the system remains asymptotically stable uniformly against all possible uncertainties within a given perturbation set. Such robust stability considerations motivate the concept of stability radius for DH systems, i.e., what is the maximal perturbation permissible to the coefficients J, R, Q, while preserving the asymptotic stability. We consider two stability radii, the unstructured one where J, R, Q are subject to unstructured perturbation, and the structured one where the perturbations preserve the DH structure. We employ characterizations for these radii that have been derived recently in SIAM J. Matrix Anal. Appl., 37, pp. 1625-1654, 2016 and propose new algorithms to compute these stability radii for large scale problems by tailoring subspace frameworks thatear rate in theory. At every iteration, they first solve a reduced problem and then expand the subspaces in order to attain certain Hermite interpolation properties between the full andreduced problems. The reduced problems are solved by means of the adaptations of existing level-set algorithms for H∞-norm computation in the unstructured case, while, for the structured radii, we benefit from algorithms that approximate the objective eigenvalue function with a piece-wise quadratic global underestimator. The performance of the new approaches is illustrated with several examples including a system that arises from a finite-element modeling of an industrial disk brake.
Subject(s): linear time-invariant dissipative Hamiltonian system
port-Hamiltonian system
robust stability
stability radius
eigenvalue optimization
subspace projection
structure preserving subspace framework
hermite interpolation
Issue Date: 26-Sep-2018
Date Available: 17-Dec-2021
Language Code: en
DDC Class: 510 Mathematik
MSC 2000: 65F15 Eigenvalues, eigenvectors
93D09 Robust stability
Series: Preprint-Reihe des Instituts für Mathematik, Technische Universität Berlin
Series Number: 2018, 09
ISSN: 2197-8085
TU Affiliation(s): Fak. 2 Mathematik und Naturwissenschaften » Inst. Mathematik
Appears in Collections:Technische Universität Berlin » Publications

Files in This Item:
AliMM18_ppt.pdf
Format: Adobe PDF | Size: 1.04 MB
DownloadShow Preview
Thumbnail

Item Export Bar

Items in DepositOnce are protected by copyright, with all rights reserved, unless otherwise indicated.