5,924 research outputs found
Computation and verification of Lyapunov functions
Lyapunov functions are an important tool to determine the basin of attraction of equilibria in Dynamical Systems through their sublevel sets. Recently, several numerical construction methods for Lyapunov functions have been proposed, among them the RBF (Radial Basis Function) and CPA (Continuous Piecewise Affine) methods. While the first method lacks a verification that the constructed function is a valid Lyapunov function, the second method is rigorous, but computationally much more demanding. In this paper, we propose a combination of these two methods, using their respective strengths: we use the RBF method to compute a potential Lyapunov function. Then we interpolate this function by a CPA function. Checking a finite number of inequalities, we are able to verify that this interpolation is a Lyapunov function. Moreover, sublevel sets are arbitrarily close to the basin of attraction. We show that this combined method always succeeds in computing and verifying a Lyapunov function, as well as in determining arbitrary compact subsets of the basin of attraction. The method is applied to two examples
Inverse heat conduction problems by using particular solutions
Based on the method of fundamental solutions, we develop in this paper a new computational method to solve two-dimensional transient heat conduction inverse problems. The main idea is to use particular solutions as radial basis functions (PSRBF) for approximation of the solutions to the inverse heat conduction problems. The heat conduction equations are first analyzed in the Laplace transformed domain and the Durbin inversion method is then used to determine the solutions in the time domain. Least-square and singular value decomposition (SVD) techniques are adopted to solve the ill-conditioned linear system of algebraic equations obtained from the proposed PSRBF method. To demonstrate the effectiveness and simplicity of this approach, several numerical examples are given with satisfactory accuracy and stability.Peer reviewe
Designing Illumination Lenses and Mirrors by the Numerical Solution of Monge-Amp\`ere Equations
We consider the inverse refractor and the inverse reflector problem. The task
is to design a free-form lens or a free-form mirror that, when illuminated by a
point light source, produces a given illumination pattern on a target. Both
problems can be modeled by strongly nonlinear second-order partial differential
equations of Monge-Amp\`ere type. In [Math. Models Methods Appl. Sci. 25
(2015), pp. 803--837, DOI: 10.1142/S0218202515500190] the authors have proposed
a B-spline collocation method which has been applied to the inverse reflector
problem. Now this approach is extended to the inverse refractor problem. We
explain in depth the collocation method and how to handle boundary conditions
and constraints. The paper concludes with numerical results of refracting and
reflecting optical surfaces and their verification via ray tracing.Comment: 16 pages, 6 figures, 2 tables; Keywords: Inverse refractor problem,
inverse reflector problem, elliptic Monge-Amp\`ere equation, B-spline
collocation method, Picard-type iteration; OCIS: 000.4430, 080.1753,
080.4225, 080.4228, 080.4298, 100.3190. Minor revision: two typos have been
corrected and copyright note has been adde
- …