4,160 research outputs found
Stochastic macromodeling of nonlinear systems via polynomial chaos expansion and transfer function trajectories
A novel approach is presented to perform stochastic variability analysis of nonlinear systems. The versatility of the method makes it suitable for the analysis of complex nonlinear electronic systems. The proposed technique is a variation-aware extension of the Transfer Function Trajectory method by means of the Polynomial Chaos expansion. The accuracy with respect to the classical Monte Carlo analysis is verified by means of a relevant numerical example showing a simulation speedup of 1777 X
Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions
Low-rank matrix approximations, such as the truncated singular value decomposition and the rank-revealing QR decomposition, play a central role in data analysis and scientific computing. This work surveys and extends recent research which demonstrates that randomization offers a powerful tool for performing low-rank matrix approximation. These techniques exploit modern computational architectures more fully than classical methods and open the possibility of dealing with truly massive data sets. This paper presents a modular framework for constructing randomized algorithms that compute partial matrix decompositions. These methods use random sampling to identify a subspace that captures most of the action of a matrix. The input matrix is then compressed—either explicitly or
implicitly—to this subspace, and the reduced matrix is manipulated deterministically to obtain the desired low-rank factorization. In many cases, this approach beats its classical competitors in terms of accuracy, robustness, and/or speed. These claims are supported by extensive numerical experiments and a detailed error analysis. The specific benefits of randomized techniques depend on the computational environment. Consider the model problem of finding the k dominant components of the singular value decomposition of an m × n matrix. (i) For a dense input matrix, randomized algorithms require O(mn log(k))
floating-point operations (flops) in contrast to O(mnk) for classical algorithms. (ii) For a sparse input matrix, the flop count matches classical Krylov subspace methods, but the randomized approach is more robust and can easily be reorganized to exploit multiprocessor architectures. (iii) For a matrix that is too large to fit in fast memory, the randomized techniques require only a constant number of passes over the data, as opposed to O(k) passes for classical algorithms. In fact, it is sometimes possible to perform matrix approximation with a single pass over the data
Stochastic Behavior Analysis of the Gaussian Kernel Least-Mean-Square Algorithm
The kernel least-mean-square (KLMS) algorithm is a popular algorithm in nonlinear adaptive filtering due to its
simplicity and robustness. In kernel adaptive filters, the statistics of the input to the linear filter depends on the parameters of the kernel employed. Moreover, practical implementations require a finite nonlinearity model order. A Gaussian KLMS has two design parameters, the step size and the Gaussian kernel bandwidth. Thus, its design requires analytical models for the algorithm behavior as a function of these two parameters. This paper studies the steady-state behavior and the transient behavior of the
Gaussian KLMS algorithm for Gaussian inputs and a finite order nonlinearity model. In particular, we derive recursive expressions for the mean-weight-error vector and the mean-square-error. The model predictions show excellent agreement with Monte Carlo simulations in transient and steady state. This allows the explicit analytical determination of stability limits, and gives opportunity
to choose the algorithm parameters a priori in order to achieve prescribed convergence speed and quality of the estimate. Design examples are presented which validate the theoretical analysis and illustrates its application
Finding structure with randomness: Probabilistic algorithms for constructing approximate matrix decompositions
Low-rank matrix approximations, such as the truncated singular value
decomposition and the rank-revealing QR decomposition, play a central role in
data analysis and scientific computing. This work surveys and extends recent
research which demonstrates that randomization offers a powerful tool for
performing low-rank matrix approximation. These techniques exploit modern
computational architectures more fully than classical methods and open the
possibility of dealing with truly massive data sets.
This paper presents a modular framework for constructing randomized
algorithms that compute partial matrix decompositions. These methods use random
sampling to identify a subspace that captures most of the action of a matrix.
The input matrix is then compressed---either explicitly or implicitly---to this
subspace, and the reduced matrix is manipulated deterministically to obtain the
desired low-rank factorization. In many cases, this approach beats its
classical competitors in terms of accuracy, speed, and robustness. These claims
are supported by extensive numerical experiments and a detailed error analysis
Measure transformation and efficient quadrature in reduced-dimensional stochastic modeling of coupled problems
Coupled problems with various combinations of multiple physics, scales, and
domains are found in numerous areas of science and engineering. A key challenge
in the formulation and implementation of corresponding coupled numerical models
is to facilitate the communication of information across physics, scale, and
domain interfaces, as well as between the iterations of solvers used for
response computations. In a probabilistic context, any information that is to
be communicated between subproblems or iterations should be characterized by an
appropriate probabilistic representation. Although the number of sources of
uncertainty can be expected to be large in most coupled problems, our
contention is that exchanged probabilistic information often resides in a
considerably lower dimensional space than the sources themselves. In this work,
we thus use a dimension-reduction technique for obtaining the representation of
the exchanged information. The main subject of this work is the investigation
of a measure-transformation technique that allows implementations to exploit
this dimension reduction to achieve computational gains. The effectiveness of
the proposed dimension-reduction and measure-transformation methodology is
demonstrated through a multiphysics problem relevant to nuclear engineering
- …