Tensor Computation: A New Framework for High-Dimensional Problems in EDA

Many critical EDA problems suffer from the curse of dimensionality, i.e. the very fast-scaling computational burden produced by large number of parameters and/or unknown variables. This phenomenon may be caused by multiple spatial or temporal factors (e.g. 3-D field solvers discretizations and multi-rate circuit simulation), nonlinearity of devices and circuits, large number of design or optimization parameters (e.g. full-chip routing/placement and circuit sizing), or extensive process variations (e.g. variability/reliability analysis and design for manufacturability). The computational challenges generated by such high dimensional problems are generally hard to handle efficiently with traditional EDA core algorithms that are based on matrix and vector computation. This paper presents "tensor computation" as an alternative general framework for the development of efficient EDA algorithms and tools. A tensor is a high-dimensional generalization of a matrix and a vector, and is a natural choice for both storing and solving efficiently high-dimensional EDA problems. This paper gives a basic tutorial on tensors, demonstrates some recent examples of EDA applications (e.g., nonlinear circuit modeling and high-dimensional uncertainty quantification), and suggests further open EDA problems where the use of tensor computation could be of advantage.Comment: 14 figures. Accepted by IEEE Trans. CAD of Integrated Circuits and System

A Perceptual Comparison of “Black Box” Modeling Algorithms for Nonlinear Audio Systems

Nonlinear systems identification is a widespread topic of interest, particularly within the audio industry, as these techniques are employed to synthesize black box models of nonlinear audio effects. Given the myriad approaches to black box modeling, questions arise as to whether an “optimal” approach exists, or one that achieves valid subjective results as a model with minimal computational expense. This thesis uses ABX listening tests to compare black box models of three hardware audio effects using two popular nonlinear implementations, along with two proposed modified implementations. Models were constructed in the Hammerstein form using sine sweeps and a novel measurement technique for the filters and nonlinearities, respectively. Testing revolved around null hypotheses assuming no change in model identification regardless of the device modeled, implementation used, or program material of the model stimulus. Results provide clear evidence of an effect on all of these accounts, and support a full rejection of the null hypotheses. Outcomes demonstrate a preferable implementation out of the algorithms tested, and suggest the removal of certain implementations as valid approaches altogether

Stochastic response determination and spectral identification of complex dynamic structural systems

Uncertainty propagation in engineering mechanics and dynamics is a highly challenging problem that requires development of analytical/numerical techniques for determining the stochastic response of complex engineering systems. In this regard, although Monte Carlo simulation (MCS) has been the most versatile technique for addressing the above problem, it can become computationally daunting when faced with high-dimensional systems or with computing very low probability events. Thus, there is a demand for pursuing more computationally efficient methodologies. Further, most structural systems are likely to exhibit nonlinear and time-varying behavior when subjected to extreme events such as severe earthquake, wind and sea wave excitations. In such cases, a reliable identification approach is behavior and for assessing its reliability. Current work addresses two research themes in the field of stochastic engineering dynamics related to the above challenges. In the first part of the dissertation, the recently developedWiener Path Integral (WPI) technique for determining the joint response probability density function (PDF) of nonlinear systems subject to Gaussian white noise excitation is generalized herein to account for non-white, non-Gaussian, and non-stationary excitation processes. Specifically, modeling the excitation process as the output of a filter equation with Gaussian white noise as its input, it is possible to define an augmented response vector process to be considered in the WPI solution technique. A significant advantage relates to the fact that the technique is still applicable even for arbitrary excitation power spectrum forms. In such cases, it is shown that the use of a filter approximation facilitates the implementation of the WPI technique in a straightforward manner, without compromising its accuracy necessarily. Further, in addition to dynamical systems subject to stochastic excitation, the technique can also account for a special class of engineering mechanics problems where the media properties are modeled as non-Gaussian and non-homogeneous stochastic fields. Several numerical examples pertaining to both single- and multi-degree-of freedom systems are considered, including a marine structural system exposed to flow-induced non-white excitation, as well as a beam with a non-Gaussian and non-homogeneous Young’s modulus. Comparisons with MCS data demonstrate the accuracy of the technique. In the second part of the dissertation, a novel multiple-input/single-output (MISO) system identification technique is developed for parameter identification of nonlinear time-variant multi-degree-of-freedom oscillators with fractional derivative terms subject to incomplete non-stationary data. The technique utilizes a representation of the nonlinear restoring forces as a set of parallel linear subsystems. In this regard, the oscillator is transformed into an equivalent MISO system in the wavelet domain. Next, a recently developed L1-norm minimization procedure based on compressive sampling theory is applied for determining the wavelet coefficients of the available incomplete non-stationary input-output (excitation-response) data. Finally, these wavelet coefficients are utilized to determine appropriately defined time- and frequency-dependent wavelet based frequency response functions and related oscillator parameters. A nonlinear time-variant system with fractional derivative elements is used as a numerical example to demonstrate the reliability of the technique even in cases of noise corrupted and incomplete data

SCEE 2008 book of abstracts : the 7th International Conference on Scientific Computing in Electrical Engineering (SCEE 2008), September 28 – October 3, 2008, Helsinki University of Technology, Espoo, Finland

This report contains abstracts of presentations given at the SCEE 2008 conference.reviewe

Low Complexity Regularization of Linear Inverse Problems

Inverse problems and regularization theory is a central theme in contemporary signal processing, where the goal is to reconstruct an unknown signal from partial indirect, and possibly noisy, measurements of it. A now standard method for recovering the unknown signal is to solve a convex optimization problem that enforces some prior knowledge about its structure. This has proved efficient in many problems routinely encountered in imaging sciences, statistics and machine learning. This chapter delivers a review of recent advances in the field where the regularization prior promotes solutions conforming to some notion of simplicity/low-complexity. These priors encompass as popular examples sparsity and group sparsity (to capture the compressibility of natural signals and images), total variation and analysis sparsity (to promote piecewise regularity), and low-rank (as natural extension of sparsity to matrix-valued data). Our aim is to provide a unified treatment of all these regularizations under a single umbrella, namely the theory of partial smoothness. This framework is very general and accommodates all low-complexity regularizers just mentioned, as well as many others. Partial smoothness turns out to be the canonical way to encode low-dimensional models that can be linear spaces or more general smooth manifolds. This review is intended to serve as a one stop shop toward the understanding of the theoretical properties of the so-regularized solutions. It covers a large spectrum including: (i) recovery guarantees and stability to noise, both in terms of $\ell^2$-stability and model (manifold) identification; (ii) sensitivity analysis to perturbations of the parameters involved (in particular the observations), with applications to unbiased risk estimation ; (iii) convergence properties of the forward-backward proximal splitting scheme, that is particularly well suited to solve the corresponding large-scale regularized optimization problem

On the Influence of Piecewise Defined Contact Geometries on Friction Dampers

Diese Arbeit betrachtet Dämpfer, die sich nicht auf eine Schwingungsreduktionsstrategie beschränken, sondern mehrere kombinieren, um optimale Ergebnisse zu erzielen. Die Möglichkeiten herkömmlicher Reibungsdämpfer werden durch stetige, stückweise definierten Kontaktgeometrien erweitert. Dies führt zu Reibungsdämpfern, die ihr Verhalten je nach Amplitude der Schwingungen ändern. Der passive, abgestimmte Keildämpfer wird entworfen und untersucht. Dieser Dämpfer bringt Dämpfung bei hohen Schwingungsamplituden in System ein und nutzt Tilgung bei niedrigen Schwingungsamplituden aus. Es werden numerische und analytische Untersuchungen durchgeführt. Um das qualitative Verhalten des Dämpfers zu validieren, wird ein Dämpferprototyp konstruiert und erprobt. Zudem wurde auch eine aktive Variante des abgestimmten Keildämpfers betrachtet. Es werden zwei Regelstrategien entworfen, die adaptive Mehrmodellregelung und die langsame, frequenzbasierte Regelung. Diese werden mit einer State-of-the-Art- Regelungsstrategie in transienten, quasistationären und Anwendungsszenarien verglichen. Die Untersuchungen zum passiven, abgestimmten Keildämpfer zeigen, dass Dämpfung und Tilgung entkoppelt werden. Eine Optimierung der Dämpferparameter ergibt im Frequenzgang eine Reduktion der Maximalamplitude von 87.47% unter Beibehaltung der Tilgung. Die Experimente validieren den Entkopplungseffekt sowie den qualitativen Einfluss der Parameter. Die aktiven Systeme erreichen mit Amplitudenabsenkungen von 91.11% das beste Ergebnis

Stochastic dynamics and wavelets techniques for system response analysis and diagnostics: Diverse applications in structural and biomedical engineering

In the first part of the dissertation, a novel stochastic averaging technique based on a Hilbert transform definition of the oscillator response displacement amplitude is developed. In comparison to standard stochastic averaging, the requirement of “a priori” determination of an equivalent natural frequency is bypassed, yielding flexibility in the ensuing analysis and potentially higher accuracy. Further, the herein proposed Hilbert transform based stochastic averaging is adapted for determining the time-dependent survival probability and first-passage time probability density function of stochastically excited nonlinear oscillators, even endowed with fractional derivative terms. To this aim, a Galerkin scheme is utilized to solve approximately the backward Kolmogorov partial differential equation governing the survival probability of the oscillator response. Next, the potential of the stochastic averaging technique to be used in conjunction with performance-based engineering design applications is demonstrated by proposing a stochastic version of the widely used incremental dynamic analysis (IDA). Specifically, modeling the excitation as a non-stationary stochastic process possessing an evolutionary power spectrum (EPS), an approximate closed-form expression is derived for the parameterized oscillator response amplitude probability density function (PDF). In this regard, IDA surfaces are determined providing the conditional PDF of the engineering demand parameter (EDP) for a given intensity measure (IM) value. In contrast to the computationally expensive Monte Carlo simulation, the methodology developed herein determines the IDA surfaces at minimal computational cost. In the second part of the dissertation, a novel multiple-input/single-output (MISO) system identification technique is developed for parameter identification of nonlinear and time-variant oscillators with fractional derivative terms subject to incomplete non-stationary data. The technique utilizes a representation of the nonlinear restoring forces as a set of parallel linear sub-systems. Next, a recently developed L1-norm minimization procedure based on compressive sensing theory is applied for determining the wavelet coefficients of the available incomplete non-stationary input-output (excitation-response) data. Several numerical examples are considered for assessing the reliability of the technique, even in the presence of incomplete and corrupted data. These include a 2-DOF time-variant Duffing oscillator endowed with fractional derivative terms, as well as a 2-DOF system subject to flow-induced forces where the non-stationary sea state possesses a recently proposed evolutionary version of the JONSWAP spectrum. In the third part of this dissertation, a joint time-frequency analysis technique based on generalized harmonic wavelets (GHWs) is developed for dynamic cerebral autoregulation (DCA) performance quantification. DCA is the continuous counter-regulation of the cerebral blood flow by the active response of cerebral blood vessels to the spontaneous or induced blood pressure fluctuations. Specifically, various metrics of the phase shift and magnitude of appropriately defined GHW-based transfer functions are determined based on data points over the joint time-frequency domain. The potential of these metrics to be used as a diagnostics tool for indicating healthy versus impaired DCA function is assessed by considering both healthy individuals and patients with unilateral carotid artery stenosis. Next, another application in biomedical engineering is pursued related to the Pulse Wave Imaging (PWI) technique. This relies on ultrasonic signals for capturing the propagation of pressure pulses along the carotid artery, and eventually for prognosis of focal vascular diseases (e.g., atherosclerosis and abdominal aortic aneurysm). However, to obtain a high spatio-temporal resolution the data are acquired at a high rate, in the order of kilohertz, yielding large datasets. To address this challenge, an efficient data compression technique is developed based on the multiresolution wavelet decomposition scheme, which exploits the high correlation of adjacent RF-frames generated by the PWI technique. Further, a sparse matrix decomposition is proposed as an efficient way to identify the boundaries of the arterial wall in the PWI technique

The concept of nonlinear modes applied to friction-damped systems

[no abstract