4,232 research outputs found
A new approach for solving nonlinear Thomas-Fermi equation based on fractional order of rational Bessel functions
In this paper, the fractional order of rational Bessel functions collocation
method (FRBC) to solve Thomas-Fermi equation which is defined in the
semi-infinite domain and has singularity at and its boundary condition
occurs at infinity, have been introduced. We solve the problem on semi-infinite
domain without any domain truncation or transformation of the domain of the
problem to a finite domain. This approach at first, obtains a sequence of
linear differential equations by using the quasilinearization method (QLM),
then at each iteration solves it by FRBC method. To illustrate the reliability
of this work, we compare the numerical results of the present method with some
well-known results in other to show that the new method is accurate, efficient
and applicable
Comparing Solution Methods for Dynamic Equilibrium Economies
This paper compares solution methods for dynamic equilibrium economies. We compute and simulate the stochastic neoclassical growth model with leisure choice using Undetermined Coefficients in levels and in logs, Finite Elements, Chebyshev Polynomials, Second and Fifth Order Perturbations and Value Function Iteration for several calibrations. We document the performance of the methods in terms of computing time, implementation complexity and accuracy and we present some conclusions about our preferred approaches based on the reported evidence.Dynamic Equilibrium Economies, Computational Methods, Linear and Nonlinear Solution Methods
A Parameterized multi-step Newton method for solving systems of nonlinear equations
We construct a novel multi-step iterative method for solving systems of nonlinear equations by introducing a parameter. to generalize the multi-step Newton method while keeping its order of convergence and computational cost. By an appropriate selection of theta, the new method can both have faster convergence and have larger radius of convergence. The new iterative method only requires one Jacobian inversion per iteration, and therefore, can be efficiently implemented using Krylov subspace methods. The new method can be used to solve nonlinear systems of partial differential equations, such as complex generalized Zakharov systems of partial differential equations, by transforming them into systems of nonlinear equations by discretizing approaches in both spatial and temporal independent variables such as, for instance, the Chebyshev pseudo-spectral discretizing method. Quite extensive tests show that the new method can have significantly faster convergence and significantly larger radius of convergence than the multi-step Newton method.Peer ReviewedPostprint (author's final draft
Comparing solution methods for dynamic equilibrium economies
This paper compares solution methods for dynamic equilibrium economies. The authors compute and simulate the stochastic neoclassical growth model with leisure choice using Undetermined Coefficients in levels and in logs, Finite Elements, Chebyshev Polynomials, Second and Fifth Order Perturbations and Value Function Iteration for several calibrations. The authors document the performance of the methods in terms of computing time, implementation complexity and accuracy and they present some conclusions about their preferred approaches based on the reported evidence.
A Probabilistic Approach to Robust Optimal Experiment Design with Chance Constraints
Accurate estimation of parameters is paramount in developing high-fidelity
models for complex dynamical systems. Model-based optimal experiment design
(OED) approaches enable systematic design of dynamic experiments to generate
input-output data sets with high information content for parameter estimation.
Standard OED approaches however face two challenges: (i) experiment design
under incomplete system information due to unknown true parameters, which
usually requires many iterations of OED; (ii) incapability of systematically
accounting for the inherent uncertainties of complex systems, which can lead to
diminished effectiveness of the designed optimal excitation signal as well as
violation of system constraints. This paper presents a robust OED approach for
nonlinear systems with arbitrarily-shaped time-invariant probabilistic
uncertainties. Polynomial chaos is used for efficient uncertainty propagation.
The distinct feature of the robust OED approach is the inclusion of chance
constraints to ensure constraint satisfaction in a stochastic setting. The
presented approach is demonstrated by optimal experimental design for the
JAK-STAT5 signaling pathway that regulates various cellular processes in a
biological cell.Comment: Submitted to ADCHEM 201
- …