Abdella, Kenzu

In this work, various aspects of neural networks, pre-trained with denoising autoencoders (DAE) are explored. To saturate neurons more quickly for feature learning in DAE, an activation function that offers higher gradients is introduced. Moreover, the introduction of sparsity functions applied to the hidden layer representations is studied. More importantly, a technique that swaps the activation functions of fully trained DAE to logistic functions is studied, networks trained using this technique are reffered to as SPAF-networks. For evaluation, the popular MNIST dataset as well as all \(3\) sub-datasets of the Chars74k dataset are used for classification purposes. The SPAF-network is also analyzed for the features it learns with a logistic, ReLU and a custom activation function. Lastly future roadmap is proposed for enhancements to the SPAF-network.

Author Keywords: Artificial Neural Network, AutoEncoder, Machine Learning, Neural Networks, SPAF network, Unsupervised Learning

Comprehending the dynamics of infectious diseases is very important in formulating public health policies to tackling their prevalence. Mathematical epidemiology (ME) has played a very vital role in achieving the above. Nevertheless, classical mathematical epidemiological models do not explicitly model the behavioural responses of individuals in the presence of prevalence of these diseases. Economic epidemiology (EE) as a field has stepped in to fill this gap by integrating economic and mathematical concepts within one framework. This thesis investigated two issues in this area. The methods employed are the standard linear analysis of stability of dynamical systems and numerical simulation. Below are the investigations and the findings of this thesis:

Firstly, an investigation into the stability properties of the equilibria of EE

models is carried out. We investigated the stability properties of modified EE systems studied by Aadland et al. [6] by introducing a parametric quadratic utility function into the model, thus making it possible to model the maximum number of contacts made by rational individuals to be determined by a parameter. This parameter in particular influences the level of utility of rational individuals. We have shown that if rational individuals have a range of possible contacts to choose from, with the maximum of the number of contacts allowable for these individuals being dependent on a parameter, the variation in this parameter tends to affect the stability properties of the system. We also showed that under the assumption of permanent recovery for

disease coupled with individuals observing or not observing their immunity, death

and birth rates can affect the stability of the system. These parameters also have

effect on the dynamics of the EE SIS system.

Secondly, an EE model of syphilis infectivity among &ldquo men who have sex with men &rdquo (MSM) in detention centres is developed in an attempt at looking at the effect of behavioural responses on the disease dynamics among MSM. This was done by explicitly incorporating the interplay of the biology of the disease and the behaviour of the inmates. We investigated the stability properties of the system under rational expectations where we showed that: (1) Behavioural responses to the prevalence of

the disease affect the stability of the system. Therefore, public health policies have the tendency of putting the system on indeterminate paths if rational MSM have complete knowledge of the laws governing the motion of the disease states as well as a complete understanding on how others behave in the system when faced with risk-benefit trade-offs. (2) The prevalence of the disease in the long run is influenced by incentives that drive the utility of the MSM inmates. (3) The interplay between the dynamics of the biology of the disease and the behavioural responses of rational MSM tends to put the system at equilibrium quickly as compared to its counterpart (that is when the system is solely dependent on the biology of the disease) when subjected to small perturbation.

Author Keywords: economic and mathematical epidemiology models, explosive path, indeterminate-path stability, numerical solution, health gap, saddle-path stability, syphilis,

In this thesis, we explore the application of the Sinc-Collocation method to an oceanography model. The model of interest describes a wind-driven current with depth-dependent eddy viscosity and is formulated in two different systems; a complex-velocity system and a real-value coupled system. In general, the Sinc-based methods excel over other traditional numerical methods due to their exponentially decaying errors, rapid convergence and handling problems in the presence of singularities at end-points. In addition, the Sinc-Collocation approach that we utilize exploits first derivative interpolation, whose integration is less sensitive to numerical errors. We present several model problems to demonstrate the accuracy, and stability of the method. We compare the approximate solutions determined by the Sinc-Collocation technique with exact solutions and also with those obtained by the Sinc-Galerkin approach in earlier studies. Our findings indicate that the method we utilized outperforms those used in past studies.

Author Keywords: Boundary Value Problems, Eddy Viscosity, Oceanography, Sinc Numerical Methods, Wind-Driven Currents

In this thesis, a new sinc-collocation method based upon derivative interpolation is developed for solving linear and nonlinear boundary value problems involving differential as well as integro-differential equations. The sinc-collocation method is chosen for its ease of implementation, exponential convergence of error, and ability to handle to singularities in the BVP. We present a unique method of treating boundary conditions and introduce the concept of the stretch factor into the conformal mappings of domains. The result is a method that achieves great accuracy while reducing computational cost. In most cases, the results from the method greatly exceed the published results of comparable methods in both accuracy and efficiency. The method is tested on the Blasius problem, the Lane-Emden problem and generalised to cover Fredholm-Volterra integro-differential problems. The results show that the sinc-collocation method with derivative interpolation is a viable and preferable method for solving nonlinear BVPs.

Author Keywords: Blasius, Boundary Value Problem, Exponential convergence, Integro-differential, Nonlinear, Sinc

The time-dependent Gross-Pitaevskii Equation, describing the movement of parti-

cles in quantum mechanics, may not be solved analytically due to its inherent non-

linearity. Hence numerical methods are of importance to approximate the solution.

This study develops a discrete scheme in time and space to simulate the solution

defined in a finite domain by using the Crank-Nicolson difference method and Sinc

Collocation Methods (SCM), respectively. In theory and practice, the time discretiz-

ing system decays errors in the second-order of accuracy, and SCMs are decaying

errors exponentially. A new SCM with a unique boundary treatment is proposed

and compared with the original SCM and other similar numerical techniques in time

costs and numerical errors. As a result, the new SCM decays errors faster than the

original one. Also, to attain the same accuracy, the new SCM interpolates fewer

nodes than the original SCM, which saves computational costs. The new SCM is

capable of approximating partial differential equations under different boundary con-

ditions, which can be extensively applied in fitting theory.

Author Keywords: Crank-Nicolson difference method, Gross-Pitaevskii Equation, Sinc-Collocation methods

In this thesis, we study modelling with non-linear ordinary differential equations, and the existence of positive solutions for Boundary Value Problems (BVPs). These problems have wide applications in many areas. The focus is on the extensions of previous work done on non-linear second-order differential equations with boundary conditions involving first-order derivative. The contribution of this thesis has four folds. First, using a fixed point theorem on order intervals, the existence of a positive solution on an interval for a non-local boundary value problem is obtained. Second, considering a different boundary value problem that consists of the first-order derivative in the non-linear term, an increasing solution is obtained by applying the Krasnoselskii-Guo fixed point theorem. Third, the existence of two solutions, one solution and no solution for a BVP is proved by using fixed point index and iteration methods. Last, the results of Green's function unify some methods in studying the existence of positive solutions for BVPs of nonlinear differential equations. Examples are presented to illustrate the applications of our results.

Author Keywords: Banach Space, Boundary Value Problems, Differential Equations, Fixed Point, Norm, Positive Solutions

With wide applications in many fields such as engineering, physics, chemistry, biology and social sciences, semi-linear equations have attracted great interests of researchers from various areas. In the study of existence of solutions for such class of equations, a general and commonly applied method is the compression cone method for fixed-point index. The main idea is to construct a cone in an ordered Banach space based on the linear part so that the nonlinear part can be examined in a relatively smaller region.

In this thesis, a new class of cone is proposed as a generalization to previous work. The construction of the cone is based on properties of both the linear and nonlinear part of the equation. As a result, the method is shown to be more adaptable in applications. We prove new results for both semi-linear integral equations and algebraic systems.

Applications are illustrated by examples. Limitations of such new method are also discussed.

Keywords: Algebraic systems; compression cone method; differential equations; existence of solutions; fixed point index; integral equations; semi-linear equations.

Author Keywords: algebraic systems, differential equations, existence of solutions, fixed point index, integral equations, semi-linear equations

This thesis is intended to support dependent-on-crops farmers to hedge the price risks of their crops. Firstly, we applied one-factor model, which incorporated a deterministic function and a stochastic process, to predict the future prices of crops (soybean). A discrete form was employed for one-month-ahead prediction. For general prediction, de-trending and de-cyclicality were used to remove the deterministic function. Three candidate stochastic differential equations (SDEs) were chosen to simulate the stochastic process; they are mean-reverting Ornstein-Uhlenbeck (OU) process, OU process with zero mean, and Brownian motion with a drift. Least squares methods and maximum likelihood were used to estimate the parameters. Results indicated that one-factor model worked well for soybean prices. Meanwhile, we provided a two-factor model as an alternative model and it also performed well in this case. In the second main part, a zero-cost option package was introduced and we theoretically analyzed the process of hedging. In the last part, option premiums obtained based on one-factor model could be compared to those obtained from Black-Scholes model, thus we could see the differences and similarities which suggested that the deterministic function especially the cyclicality played an essential role for the soybean price, thus the one-factor model in this case was more suitable than Black-Scholes model for the underlying asset.

Author Keywords: Brownian motion, Least Squares Method, Maximum Likelihood Method, One-factor Model, Option Pricing, Ornstein-Uhlenbeck Process

Artificial Neural Networks (ANN) have been applied to many areas of research. These techniques use a series of object attributes and can be trained to recognize different classes of objects. The Self-Organizing Map (SOM) is an unsupervised machine learning technique which has been shown to be successful in the mapping of high-dimensional data into a 2D representation referred to as a map. These maps are easier to interpret and aid in the classification of data. In this work, the existing algorithms for the SOM have been extended to generate 3D maps. The higher dimensionality of the map provides for more information to be made available to the interpretation of classifications. The effectiveness of the implementation was verified using three separate standard datasets. Results from these investigations supported the expectation that a 3D SOM would result in a more effective classifier.

The 3D SOM algorithm was then applied to an analysis of galaxy morphology classifications. It is postulated that the morphology of a galaxy relates directly to how it will evolve over time. In this work, the Spectral Energy Distribution (SED) will be used as a source for galaxy attributes. The SED data was extracted from the NASA Extragalactic Database (NED). The data was grouped into sample sets of matching frequencies and the 3D SOM application was applied as a morphological classifier. It was shown that the SOMs created were effective as an unsupervised machine learning technique to classify galaxies based solely on their SED. Morphological predictions for a number of galaxies were shown to be in agreement with classifications obtained from new observations in NED.

Author Keywords: Galaxy Morphology, Multi-wavelength, parallel, Self-Organizing Maps