Large-scale estimation of dominant poles of a transfer function by an interpolatory framework

DOI10.1137/21M1434349zbMATH Open1501.65018arXiv2107.08286OpenAlexW3186034741MaRDI QIDQ5097600FDOQ5097600

Authors: Emre Mengi 

Publication date: 25 August 2022

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

Abstract: We focus on the dominant poles of the transfer function of a descriptor system. The transfer function typically exhibits large norm at and near the imaginary parts of the dominant poles. Consequently, the dominant poles provide information about the points on the imaginary axis where the ${mathcalL}_{i}nfty$ norm of the system is attained, and they are also sometimes useful to obtain crude reduced-order models. For a large-scale descriptor system, we introduce a subspace framework to estimate a prescribed number of dominant poles. At every iteration, the large-scale system is projected into a small system, whose dominant poles can be computed at ease. Then the projection spaces are expanded so that the projected system after subspace expansion interpolates the large-scale system at the computed dominant poles. We prove an at-least-quadratic-convergence result for the framework, and provide numerical results confirming this. On real benchmark examples, the proposed framework appears to be more reliable than SAMDP [IEEE Trans. Power Syst. 21, 1471-1483, 2006], one of the widely used algorithms due to Rommes and Martins for the estimation of the dominant poles.

Full work available at URL: https://arxiv.org/abs/2107.08286

Recommendations

• Computing Transfer Function Dominant Poles of Large-Scale Second-Order Dynamical Systems
• scientific article
• Convergence of the Dominant Pole Algorithm and Rayleigh Quotient Iteration
• Eigenvalue methods for calculating dominant poles of a transfer function and their applications in small-signal stability
• Model reduction of descriptor systems by interpolatory projection methods

zbMATH Keywords

Hermite interpolationmodel order reductiondescriptor systemlarge scale problemdominant pole

Mathematics Subject Classification ID

Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Transformation and reduction of functional-differential equations and systems, normal forms (34K17) Large-scale systems (93A15) Linear systems in control theory (93C05)

Cites Work

• The AAA algorithm for rational approximation
• Title not available (Why is that?)
• Approximation of Large-Scale Dynamical Systems
• Port-Hamiltonian Systems Theory: An Introductory Overview
• Interpolatory model reduction of large-scale dynamical systems
• Stability radius for structured perturbations and the algebraic Riccati equation
• \(H_\infty\)-norm-based optimization for the identification of gray-box LTI state-space model parameters
• $\mathcal{H}_2$ Model Reduction for Large-Scale Linear Dynamical Systems
• Interpolatory projection methods for structure-preserving model reduction
• Model reduction of descriptor systems by interpolatory projection methods
• Model Reduction of MIMO Systems via Tangential Interpolation
• A regularity result for the singular values of a transfer matrix and a quadratically convergent algorithm for computing its \(L_{\infty}\)-norm
• A fast algorithm to compute the \(H_{\infty}\)-norm of a transfer function matrix
• Fast approximation of the \(H_\infty\) norm via optimization over spectral value sets
• Title not available (Why is that?)
• Convergence of the Dominant Pole Algorithm and Rayleigh Quotient Iteration
• Large-scale computation of \(\mathcal{L}_\infty\)-norms by a greedy subspace method
• Approximation of stability radii for large-scale dissipative Hamiltonian systems
• Stability radii for linear Hamiltonian systems with dissipation under structure-preserving perturbations
• A structured pseudospectral method for \(\mathcal {H}_{\infty}\)-norm computation of large-scale descriptor systems
• Calculating the \(H_{\infty}\)-norm using the implicit determinant method
• Hybrid expansion-contraction: a robust scaleable method for approximating the \(H_{\infty}\) norm
• Faster and more accurate computation of the \(\mathcal{H}_\infty\) norm via optimization
• A Subspace Framework for ${\mathcal H}_\infty$-Norm Minimization
• Derivative Interpolating Subspace Frameworks for Nonlinear Eigenvalue Problems
• Title not available (Why is that?)

Cited In (3)

• Estimation of input pulse locations from the response of an all-pole transfer system using tapered rank reduction
• Computing Transfer Function Dominant Poles of Large-Scale Second-Order Dynamical Systems
• Title not available (Why is that?)

This page was built for publication: Large-scale estimation of dominant poles of a transfer function by an interpolatory framework

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5097600)