78,233 research outputs found
Optimal designs for rational function regression
We consider optimal non-sequential designs for a large class of (linear and
nonlinear) regression models involving polynomials and rational functions with
heteroscedastic noise also given by a polynomial or rational weight function.
The proposed method treats D-, E-, A-, and -optimal designs in a
unified manner, and generates a polynomial whose zeros are the support points
of the optimal approximate design, generalizing a number of previously known
results of the same flavor. The method is based on a mathematical optimization
model that can incorporate various criteria of optimality and can be solved
efficiently by well established numerical optimization methods. In contrast to
previous optimization-based methods proposed for similar design problems, it
also has theoretical guarantee of its algorithmic efficiency; in fact, the
running times of all numerical examples considered in the paper are negligible.
The stability of the method is demonstrated in an example involving high degree
polynomials. After discussing linear models, applications for finding locally
optimal designs for nonlinear regression models involving rational functions
are presented, then extensions to robust regression designs, and trigonometric
regression are shown. As a corollary, an upper bound on the size of the support
set of the minimally-supported optimal designs is also found. The method is of
considerable practical importance, with the potential for instance to impact
design software development. Further study of the optimality conditions of the
main optimization model might also yield new theoretical insights.Comment: 25 pages. Previous version updated with more details in the theory
and additional example
Optimal stability polynomials for numerical integration of initial value problems
We consider the problem of finding optimally stable polynomial approximations
to the exponential for application to one-step integration of initial value
ordinary and partial differential equations. The objective is to find the
largest stable step size and corresponding method for a given problem when the
spectrum of the initial value problem is known. The problem is expressed in
terms of a general least deviation feasibility problem. Its solution is
obtained by a new fast, accurate, and robust algorithm based on convex
optimization techniques. Global convergence of the algorithm is proven in the
case that the order of approximation is one and in the case that the spectrum
encloses a starlike region. Examples demonstrate the effectiveness of the
proposed algorithm even when these conditions are not satisfied
On the interpretation and identification of dynamic Takagi-Sugenofuzzy models
Dynamic Takagi-Sugeno fuzzy models are not always easy to interpret, in particular when they are identified from experimental data. It is shown that there exists a close relationship between dynamic Takagi-Sugeno fuzzy models and dynamic linearization when using affine local model structures, which suggests that a solution to the multiobjective identification problem exists. However, it is also shown that the affine local model structure is a highly sensitive parametrization when applied in transient operating regimes. Due to the multiobjective nature of the identification problem studied here, special considerations must be made during model structure selection, experiment design, and identification in order to meet both objectives. Some guidelines for experiment design are suggested and some robust nonlinear identification algorithms are studied. These include constrained and regularized identification and locally weighted identification. Their usefulness in the present context is illustrated by examples
Rational minimax approximation via adaptive barycentric representations
Computing rational minimax approximations can be very challenging when there
are singularities on or near the interval of approximation - precisely the case
where rational functions outperform polynomials by a landslide. We show that
far more robust algorithms than previously available can be developed by making
use of rational barycentric representations whose support points are chosen in
an adaptive fashion as the approximant is computed. Three variants of this
barycentric strategy are all shown to be powerful: (1) a classical Remez
algorithm, (2) a "AAA-Lawson" method of iteratively reweighted least-squares,
and (3) a differential correction algorithm. Our preferred combination,
implemented in the Chebfun MINIMAX code, is to use (2) in an initial phase and
then switch to (1) for generically quadratic convergence. By such methods we
can calculate approximations up to type (80, 80) of on in
standard 16-digit floating point arithmetic, a problem for which Varga, Ruttan,
and Carpenter required 200-digit extended precision.Comment: 29 pages, 11 figure
- âŠ