Fast computation of optimal damping parameters for linear vibrational systems

by   N. Jakovcevic Stor, et al.

We formulate the quadratic eigenvalue problem underlying the mathematical model of a linear vibrational system as an eigenvalue problem of a diagonal-plus-low-rank matrix A. The eigenvector matrix of A has a Cauchy-like structure. Optimal viscosities are those for which trace(X) is minimal, where X is the solution of the Lyapunov equation AX+XA^*=GG^*. Here G is a low-rank matrix which depends on the eigenfrequencies that need to be damped. After initial eigenvalue decomposition of linearized problem which requires O(n^3) operations, our algorithm computes optimal viscosities for each choice of external dampers in O(n^2) operations, provided that the number of dampers is small. Hence, the subsequent optimization is order of magnitude faster than in the standard approach which solves Lyapunov equation in each step, thus requiring O(n^3) operations. Our algorithm is based on O(n^2) eigensolver for complex symmetric diagonal-plus-rank-one matrices and fast O(n^2) multiplication of linked Cauchy-like matrices.


page 1

page 2

page 3

page 4


Fast optimization of viscosities for frequency-weighted damping of second-order systems

We consider frequency-weighted damping optimization for vibrating system...

A fast iterative algorithm for near-diagonal eigenvalue problems

We introduce a novel iterative eigenvalue algorithm for near-diagonal ma...

Optimal Diagonal Preconditioning: Theory and Practice

Preconditioning has been a staple technique in optimization and machine ...

SuperDC: Stable superfast divide-and-conquer eigenvalue decomposition

For dense Hermitian matrices with small off-diagonal (numerical) ranks a...

Rank-1 matrix differential equations for structured eigenvalue optimization

A new approach to solving eigenvalue optimization problems for large str...

Fast Coherent Point Drift

Nonrigid point set registration is widely applied in the tasks of comput...

Please sign up or login with your details

Forgot password? Click here to reset