Search results for "linear equations"
showing 10 items of 64 documents
Beyond Viner: Smoothed Cost Curves and Co-Detetermination of Output and Production Capacity
2017
We address two problems of traditional cost functions: the discontinuity caused by the production capacity (the marginal cost abruptly becomes infinite when production capacity is reached) and the production capacity is artificially exogenous. So, we introduce a smoothed form of marginal cost function. It progressively tends to infinity when it approaches to the production capacity. Then, we prove that it is perfectly possible to determine directly output, fixed costs and production capacity simultaneously, even if this could lead to a system of equations that is not so elementary to solve because it includes a Lambert function. We also show that the smoothed cost function prevails in both …
Field estimation in wireless sensor networks using distributed kriging
2012
In this paper, we tackle the problem of spatial interpolation for distributed estimation in Wireless Sensor Networks by using a geostatistical technique called kriging. We present a novel Distributed Iterative Kriging Algorithm (DIKA) which is composed of two main phases. First, the spatial dependence of the field is exploited by calculating semivariograms in an iterative way. Second, the kriging system of equations is solved by an initial set of nodes in a distributed manner, providing some initial interpolation weights to each node. In our algorithm, the estimation accuracy can be improved by iteratively adding new nodes and updating appropriately the weights, which leads to a reduction i…
Delay-Range-Dependent Linear Matrix Inequality Approach to Quantized H∞ Control of Linear Systems with Network-Induced Delays and Norm-Bounded Uncert…
2010
This paper deals with a convex optimization approach to the problem of robust network-based H∞ control for linear systems connected over a common digital communication network with static quantizers. Both the polytopic and the norm-bounded uncertainties are taken into consideration separately. First, the effect of both the output quantization levels and the network conditions under static quantizers is investigated. Second, by introducing a descriptor technique, using a Lyapunov—Krasovskii functional and a suitable change of variables, new required sufficient conditions are established in terms of delay-range-dependent linear matrix inequalities for the existence of the desired network-bas…
A fast 3D dual boundary element method based on hierarchical matrices
2008
AbstractIn this paper a fast solver for three-dimensional BEM and DBEM is developed. The technique is based on the use of hierarchical matrices for the representation of the collocation matrix and uses a preconditioned GMRES for the solution of the algebraic system of equations. The preconditioner is built exploiting the hierarchical arithmetic and taking full advantage of the hierarchical format. Special algorithms are developed to deal with crack problems within the context of DBEM. The structure of DBEM matrices has been efficiently exploited and it has been demonstrated that, since the cracks form only small parts of the whole structure, the use of hierarchical matrices can be particula…
On Discovering Low Order Models in Biochemical Reaction Kinetics
2007
We develop a method by which a large number of differential equations representing biochemical reaction kinetics may be represented by a smaller number of differential equations. The basis of our technique is a conjecture that the high dimension equations of biochemical kinetics, which involve reaction terms of specific forms, are actually implementing a low dimension system whose behavior requires right hand sides that can not be biochemically implemented. For systems that satisfy this conjecture, we develop a simple approximation scheme based on multilinear algebra that extracts the low dimensional system from simulations of the high dimension system. We demonstrate this technique on a st…
Nonresistive dissipative magnetohydrodynamics from the Boltzmann equation in the 14-moment approximation
2018
We derive the equations of motion of relativistic, non-resistive, second-order dissipative magnetohydrodynamics from the Boltzmann equation using the method of moments. We assume the fluid to be composed of a single type of point-like particles with vanishing dipole moment or spin, so that the fluid has vanishing magnetization and polarization. In a first approximation, we assume the fluid to be non-resistive, which allows to express the electric field in terms of the magnetic field. We derive equations of motion for the irreducible moments of the deviation of the single-particle distribution function from local thermodynamical equilibrium. We analyze the Navier-Stokes limit of these equati…
Efficient numerical methods for pricing American options under stochastic volatility
2007
Five numerical methods for pricing American put options under Heston's stochastic volatility model are described and compared. The option prices are obtained as the solution of a two-dimensional parabolic partial differential inequality. A finite difference discretization on nonuniform grids leading to linear complementarity problems with M-matrices is proposed. The projected SOR, a projected multigrid method, an operator splitting method, a penalty method, and a componentwise splitting method are considered. The last one is a direct method while all other methods are iterative. The resulting systems of linear equations in the operator splitting method and in the penalty method are solved u…
A fast hierarchical dual boundary element method for three-dimensional elastodynamic crack problems
2010
In this work a fast solver for large-scale three-dimensional elastodynamic crack problems is presented, implemented, and tested. The dual boundary element method in the Laplace transform domain is used for the accurate dynamic analysis of cracked bodies. The fast solution procedure is based on the use of hierarchical matrices for the representation of the collocation matrix for each computed value of the Laplace parameter. An ACA (adaptive cross approximation) algorithm is used for the population of the low rank blocks and its performance at varying Laplace parameters is investigated. A preconditioned GMRES is used for the solution of the resulting algebraic system of equations. The precond…
Guaranteed lower bounds for cost functionals of time-periodic parabolic optimization problems
2019
In this paper, a new technique is shown for deriving computable, guaranteed lower bounds of functional type (minorants) for two different cost functionals subject to a parabolic time-periodic boundary value problem. Together with previous results on upper bounds (majorants) for one of the cost functionals, both minorants and majorants lead to two-sided estimates of functional type for the optimal control problem. Both upper and lower bounds are derived for the second new cost functional subject to the same parabolic PDE-constraints, but where the target is a desired gradient. The time-periodic optimal control problems are discretized by the multiharmonic finite element method leading to lar…
Scheduled Relaxation Jacobi method: improvements and applications
2016
Elliptic partial differential equations (ePDEs) appear in a wide variety of areas of mathematics, physics and engineering. Typically, ePDEs must be solved numerically, which sets an ever growing demand for efficient and highly parallel algorithms to tackle their computational solution. The Scheduled Relaxation Jacobi (SRJ) is a promising class of methods, atypical for combining simplicity and efficiency, that has been recently introduced for solving linear Poisson-like ePDEs. The SRJ methodology relies on computing the appropriate parameters of a multilevel approach with the goal of minimizing the number of iterations needed to cut down the residuals below specified tolerances. The efficien…