Computational aspects of the numerical solution of SDEs

In the last 30 to 40 years, many researchers have combined to build the knowledge base of theory and solution techniques that can be applied to the case of differential equations which include the effects of noise. This class of ``noisy'' differential equations is now known as stochastic differential equations (SDEs). Markov diffusion processes are included within the field of SDEs through the drift and diffusion components of the Itô form of an SDE. When these drift and diffusion components are moderately smooth functions, then the processes' transition probability densities satisfy the Fokker-Planck-Kolmogorov (FPK) equation -- an ordinary partial differential equation (PDE). Thus there is a mathematical inter-relationship that allows solutions of SDEs to be determined from the solution of a noise free differential equation which has been extensively studied since the 1920s. The main numerical solution technique employed to solve the FPK equation is the classical Finite Element Method (FEM). The FEM is of particular importance to engineers when used to solve FPK systems that describe noisy oscillators. The FEM is a powerful tool but is limited in that it is cumbersome when applied to multidimensional systems and can lead to large and complex matrix systems with their inherent solution and storage problems. I show in this thesis that the stochastic Taylor series (TS) based time discretisation approach to the solution of SDEs is an efficient and accurate technique that provides transition and steady state solutions to the associated FPK equation. The TS approach to the solution of SDEs has certain advantages over the classical techniques. These advantages include their ability to effectively tackle stiff systems, their simplicity of derivation and their ease of implementation and re-use. Unlike the FEM approach, which is difficult to apply in even only two dimensions, the simplicity of the TS approach is independant of the dimension of the system under investigation. Their main disadvantage, that of requiring a large number of simulations and the associated CPU requirements, is countered by their underlying structure which makes them perfectly suited for use on the now prevalent parallel or distributed processing systems. In summary, l will compare the TS solution of SDEs to the solution of the associated FPK equations using the classical FEM technique. One, two and three dimensional FPK systems that describe noisy oscillators have been chosen for the analysis. As higher dimensional FPK systems are rarely mentioned in the literature, the TS approach will be extended to essentially infinite dimensional systems through the solution of stochastic PDEs. In making these comparisons, the advantages of modern computing tools such as computer algebra systems and simulation software, when used as an adjunct to the solution of SDEs or their associated FPK equations, are demonstrated.

Refined Jensen-based inequality approach to stability analysis of time-delay systems

© 2015 The Institution of Engineering and Technology. In this study, the authors derive some new refined Jensen-based inequalities, which encompass both the Jensen inequality and its most recent improvement based on the Wirtinger integral inequality. The potential capability of this approach is demonstrated through applications to stability analysis of time-delay systems. More precisely, by using the newly derived inequalities, they establish new stability criteria for two classes of time-delay systems, namely discrete and distributed constant delays systems and interval time-varying delay systems. The resulting stability conditions are derived in terms of linear matrix inequalities, which can be efficiently solved by various convex optimisation algorithms. Numerical examples are given to show the effectiveness and least conservativeness of the results obtained in this study.

Demographic demand systems with application to equivalence scales estimation and inequality analysis : the Australian evidence

This paper proposes and applies an alternative demographic procedure for extending a demand system to allow for the effect of household size and composition changes, along with price changes, on expenditure allocation. The demographic procedure is applied to two recent demand functional forms to obtain their estimable demographic extensions. The estimation on pooled time series of Australian Household Expenditure Surveys yields sensible and robust estimates of the equivalence scale, and of its variation with relative prices. Further evidence on the usefulness of this procedure is provided by using it to evaluate the nature and magnitude of the inequality bias of relative price changes in Australia over a period from the late 1980s to the early part of the new millennium.

Stochastic stability of nonlinear discrete-time Markovian jump systems with time-varying delay and partially unknown transition rates

This paper is concerned with stochastic stability of a class of nonlinear discrete-time Markovian jump systems with interval time-varying delay and partially unknown transition probabilities. A new weighted summation inequality is first derived. We then employ the newly derived inequality to establish delay-dependent conditions which guarantee the stochastic stability of the system. These conditions are derived in terms of tractable matrix inequalities which can be computationally solved by various convex optimized algorithms. Numerical examples are provided to illustrate the effectiveness of the obtained results.

Static output feedback frequency stabilization of time-delay power systems with coordinated electric vehicles state of charge control

In this paper, for the first time, electric vehicles are used for both the primary and secondary frequency controls to support power plants to rapidly suppress fluctuations in the system frequency due to load disturbances. Via networked control and wide-area communication infrastructures, multiple interval time-varying delays exist in the communication channels between the control center, power plant, and an aggregation of electric vehicles. By coordinating batteries’ state of charge control, the behaviors of the vehicle owners and the uncertainties imposed by the changes of the batteries’ state of charge are taken intoconsideration. A power system model incorporating multiple time-varying delays and uncertainties is first proposed. Then, a robust static output feedback frequency controller is designed to guarantee the resulting closed-loop system stable with an H∞ attenuation level. By utilizing a novel integral inequality, namely refined-Jensen inequality, and an improved reciprocally convex combination, the design conditions are formulated in terms of tractable linear matrix inequalities which can be efficiently solved by various computational tools. The effectiveness of the proposed control scheme is verified by extensive simulations.

Convex hulls as an hypothesis language bias

Classification learning is dominated by systems which induce large numbers of small axis-orthogonal decision surfaces which biases such systems towards particular hypothesis types. However, there is reason to believe that many domains have underlying concepts which do not involve axis orthogonal surfaces. Further, the multiplicity of small decision regions mitigates against any holistic appreciation of the theories produced by these systems, notwithstanding the fact that many of the small regions are individually comprehensible. We propose the use of less strongly biased hypothesis languages which might be expected to model' concepts using a number of structures close to the number of actual structures in the domain. An instantiation of such a language, a convex hull based classifier, CHI, has been implemented to investigate modeling concepts as a small number of large geometric structures in n-dimensional space. A comparison of the number of regions induced is made against other well-known systems on a representative selection of largely or wholly continuous valued machine learning tasks. The convex hull system is shown to produce a number of induced regions about an order of magnitude less than well-known systems and very close to the number of actual concepts. This representation, as convex hulls, allows the possibility of extraction of higher level mathematical descriptions of the induced concepts, using the techniques of computational geometry.

State and input simultaneous estimation for a class of nonlinear systems

This paper addresses the problem of estimating simultaneously the state and input of a class of nonlinear systems. Here, the systems nonlinear part comprises a Lipschitz nonlinear function with respect to the state and input, and a state-dependent unknown function including additive disturbance as well as uncertain/nonlinear/time-varying terms. Upon satisfying some conditions, the observer design problem can be solved via a Riccati inequality or a LMI-based technique with asymptotic estimation guaranteed. A numerical example is included for illustration.

State and input simultaneous estimation for a class of time-delay systems with uncertainties

This brief addresses the problem of estimation of both the states and the unknown inputs of a class of systems that are subject to a time-varying delay in their state variables, to an unknown input, and also to an additive uncertain, nonlinear disturbance. Conditions are derived for the solvability of the design matrices of a reduced-order observer for state and input estimation, and for the stability of its dynamics. To improve computational efficiency, a delay-dependent asymptotic stability condition is then developed using the linear matrix inequality formulation. A design procedure is proposed and illustrated by a numerical example.

Robust strictly positive real synthesis of polynomial segments for discrete time systems

By using the result of robust strictly positive real synthesis of polynomial segments for continuous time systems, it is proved that, for any two n-th order polynomials a(z) and b(z), the Schur stability of their convex combination is necessary and sufficient for the existence of an n-th order polynomial c(z) such that c(z)/a(z) and c(z)/b(z) are both strictly positive real. We also provide the construction method of c(z). Illustrative examples are provided to show the effectiveness of this method.

Convex hulls in concept induction

Classification learning is dominated by systems which induce large numbers of small axis-orthogonal decision surfaces. This strongly biases such systems towards particular hypothesis types but there is reason believe that many domains have underlying concepts which do not involve axis orthogonal surfaces. Further, the multiplicity of small decision regions mitigates against any holistic appreciation of the theories produced by these systems, notwithstanding the fact that many of the small regions are individually comprehensible. This thesis investigates modeling concepts as large geometric structures in n-dimensional space. Convex hulls are a superset of the set of axis orthogonal hyperrectangles into which axis orthogonal systems partition the instance space. In consequence, there is reason to believe that convex hulls might provide a more flexible and general learning bias than axis orthogonal regions. The formation of convex hulls around a group of points of the same class is shown to be a usable generalisation and is more general than generalisations produced by axis-orthogonal based classifiers, without constructive induction, like decision trees, decision lists and rules. The use of a small number of large hulls as a concept representation is shown to provide classification performance which can be better than that of classifiers which use a large number of small fragmentary regions for each concept. A convex hull based classifier, CH1, has been implemented and tested. CH1 can handle categorical and continuous data. Algorithms for two basic generalisation operations on hulls, inflation and facet deletion, are presented. The two operations are shown to improve the accuracy of the classifier and provide moderate classification accuracy over a representative selection of typical, largely or wholly continuous valued machine learning tasks. The classifier exhibits superior performance to well-known axis-orthogonal-based classifiers when presented with domains where the underlying decision surfaces are not axis parallel. The strengths and weaknesses of the system are identified. One particular advantage is the ability of the system to model domains with approximately the same number of structures as there are underlying concepts. This leads to the possibility of extraction of higher level mathematical descriptions of the induced concepts, using the techniques of computational geometry, which is not possible from a multiplicity of small regions.

Improvement on delay dependent absolute stability of Lurie control systems with multiple time delays

Absolute stability of Lurie control systems with multiple time-delays is studied in this paper. By using extended Lyapunov functionals, we avoid the use of the stability assumption on the main operator and derive improved stability criteria, which are strictly less conservative than the criteria in [2,3].

Parameter-dependent H∞ control for time varying delay polytopic systems

This paper addresses the robust stabilization and H_∞ control problem for a class of linear polytopic systems with continuously distributed delays. The control objective is to design a robust H_∞ controller that satisfies some exponential stability constraints on the closed-loop poles. Using improved parameter-dependent Lyapunov Krasovskii functionals, new delay-dependent conditions for the robust H_∞ control are established in terms of linear matrix inequalities.

Partial state estimation for linear systems with output and input time delays.

This paper deals with the problem of partial state observer design for linear systems that are subject to time delays in the measured output as well as the control input. By choosing a set of appropriate augmented Lyapunov-Krasovskii functionals with a triple-integral term and using the information of both the delayed output and input, a novel approach to design a minimal-order observer is proposed to guarantee that the observer error is ε-convergent with an exponential rate. Existence conditions of such an observer are derived in terms of matrix inequalities for the cases with time delays in both the output and input and with output delay only. Constructive design algorithms are introduced. Numerical examples are provided to illustrate the design procedure, practicality and effectiveness of the proposed observer.