Tawfiq LNM* and Al-Abrahemee KMM
Department of Mathematics, College of Education for Pure Science/Ibn Al-Haitham, University of Baghdad, Baghdad, Iraq
Received Date: February 04, 2014; Accepted Date: March 25, 2014; Published Date: April 14, 2014
Citation: Tawfiq LNM, Al-Abrahemee KMM (2014) Design Neural Network to Solve Singular Perturbation Problems. J Appl Computat Math 3:160 doi: 10.4172/2168-9679.1000160
Copyright: © 2014 Tawfiq LNM, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Applied & Computational Mathematics
The aim of this paper is to design neural network to present a method to solve Singular perturbation problems (SPP) by using network having one hidden layer with 5 hidden units (neurons) and one linear output unit, the sigmoid activation of each hidden units is tansigmoid. The neural network trained by the back propagation with different algorithms such as quasi-Newton, Levenberg-Marquardt, and Bayesian Regulation. Finally the results of numerical experiments are compared with the exact solution in illustrative examples to confirm the accuracy and efficiency of the presented scheme.
Singularly perturbed problems; Neural network; Quasi- Newton; Levenberg-Marquardt, Bayesian regulation
Singularly perturbed problems (SPP) are common in applied sciences and engineering. They often occur in, for example, fluid dynamics, quantum mechanics, chemical reactions, electrical networks, etc. A well known fact is that the solution of such problems has a multi scale character, i.e., there are thin transition layers where the solution varies very rapidly, while away from the layers the solution behaves regularly and varies slowly. For a detailed discussion on the analytical and numerical treatment of such problems one may refer to the Malley , Doolan et al. , Roos et al. , and Miller et al. . Numerically, the presence of the perturbation parameter leads to difficulties when classical numerical techniques are used to solve such problems, this is due to the presence of the boundary layers in these problems, see for example [1,5]. Even in the case when only the approximate solution is required.
Many methods have been developed so far solving Singularly perturbed boundary value problems (SPBVP) , nowadays there is a new way of computing denominated artificial intelligence which through different methods is capable of managing the imprecision's and uncertainties that appear when trying to solve problems related to the real world, offering strong solution and of easy implementation. One of those techniques is known as Artificial Neural Networks (ANN). Inspired, in their origin, in the functioning of the human brain, and entitled with some intelligence. These are the combination of a great amount of elements of process–artificial neurons interconnected that operating in a parallel way get to solve problems related to aspects of classification. The construction of any given ANN we can identify, depending on the location in the network, three kind of computational neurons: input, output and hidden.
The term "perturbation problem" is generally used in mathematics when one deals with the following situation: There is a family of problems depending on a small parameter ϵ> 0,, which we denote by Pϵ, when ϵ=0, we have the reduced problem Po.
We want to study the relationship between the solution of Pϵ and the solution of Po under appropriate assumptions. The perturbation problem (PP), may consist of an ordinary differential equation, or a system of differential equations, dong with some given conditions which illustrate the problem. Thus, the general form of the 2nd order singularly perturbed problems (SPP):
where f are n-dimensional vector functions, x is a scalar variable in a given interval.
A perturbation problem (1) is called SPP if ϵ→ 0, the solution yϵ(x) converges to yo(x) only in some interval of x, but not throughout the entire interval, thus giving rise to an "boundary layers" phenomena at both end-Points .
Artificial neural network is a simplified mathematical model of the human brain. It can be implemented by both electric elements and computer software. It is a parallel distributed processor with large numbers of connections; it is an information processing system that has certain performance characters in common with biological neural networks .
The arriving signals, called inputs, multiplied by the connection weights (adjusted) are first summed (combined) and then passed through a transfer function to produce the output for that neuron. The transfer (activation) function acts on the weighted sum of the neuron’s inputs and the most commonly used transfer function is the sigmoid function (tansig.) .
There are two main connection formulas (types): feedback (recurrent) and feed forward connections. Feedback is one type of connection where the output of one layer routes back to the input of a previous layer, or to the same layer. Feed forward neural network (FFNN) does not have a connection back from the output to the input neurons . There are many different training algorithms, but the most often used training algorithm is the back propagation (BP) rule. A NN is trained to map a set of input data by iterative adjustment of the weights. Information from inputs is fed forward through the network to optimize the weights between neurons. Optimization of the weights is made by backward propagation of the error during training phase. The ANN reads the input and output values in the training data set and changes the value of the weighted links to reduce the difference between the predicted and target (observed) values. The error in prediction is minimized across many training cycles (iteration or ePoch) until network reaches specified level of accuracy. A complete round of forward backward passes and weight adjustments using all input output pairs in the data set is called an ePoch or iteration. In order to perform a supervised training we need a way of evaluating the ANN output error between the actual and the expected outputs. A Popular measure is the mean squared error (MSE) or root mean squared error (RMSE) .
In an ANN expressions structure, architecture or toPology, express the way in which computational neurons are organized in the network .
Particularly, these terms are focused in the description of how the nodes are connected and in how the information is transmitted through the network. As it has been mentioned, the distribution of computational in the following:
Number of levels or layers
Neurons in the neural network is done forming levels or layers of a determined number of nodes each one. As there are input, output and hidden neurons, we can talk about an input layer, an output layer and single layer or multilayer hidden layers. By the peculiarity of the behavior of the input nodes some authors consider just two kinds of layers in the ANN, the hidden and the output.
Depending on the links between the elements of the different layers. The ANN can be classified as: totally connected, when all the outputs from a level get to all and each one of the nodes in the following level, if some of the links in the network are lost, then we say that the network is partially connected.
Another classification of the ANN is obtained by considering the direction of the flow of the through the layers, when any output of the neurons is input of neurons of the same level or preceding levels, the network is described as feed forward. In counter Position if there is at least one connected exit as entrance of neurons of previous levels or of the same level, including themselves, the network is denominated of feedback.
This section illustrates how our approach can be used to the approximation solution of the singular perturbation problem of the form:
where x ϵ D, D ⊂ R denoted the domain and y(x) is the solution to be computed.
If yt (x, p) denoted a trial solution with adjustable parameters p, the problem is transformed to a discredited from:
In our proPosed approach, the trial solution yt employs an ANN and the parameters p corresPond to the weights and biases of the neural architecture.
We choose a form for the trial function (x) such that achieved by writing as a sum of two terms:
where N(x, p) is ANN with parameters p and n input units fed with the input vector x. The term (x) contains no adjustable parameters and satisfies the BCs. The second term G is constructed so as not to contribute to the BCs, since (x) satisfy them. This term can be formed by using ANN whose weights and biases are to be adjusted in order to deal with the minimization problem.
In this section, we describe solution of SPP using ANN. To illustrate the method, we will consider the 2nd order SPP:
where x ∈ [a, b] and the BC: y(a)=A, y(b)=B; a trial solution can be written as:
where N(x, p) is the output of the ANN with one input unit for x and weights p.
The error quantity to be minimized is given by:
where the xi ∈ [a, b]. Since
It is straight forward to compute the gradient of the error with respect to the parameters p.
In this section, we rePort numerical result, using amulti-layer ANN having one hidden layer with 5 hidden units (neurons) and one linear output unit. The sigmoid activation of each hidden unit is tansig, the analytic solution (x) was known in advance. Therefore, we test the accuracy of obtained solutions computing the mean square error (MSE).
Consider the following 2nd order singular perturbation problem:
BC′s (Dirishlit case): y(0)=-1, y(1)=3.9995, and the analytic solution :
c1 = 2.9995
according to the equation (6 ) the trial neural form of the solution is taken to be:
The ANN trained using a grid of ten equidistant Points in 0, 1. Figure 1 displays the analytic and neural solutions with different training algorithm. The neural results with different types of training algorithm such as: Levenberg–Marquardt (trainlm), quasi–Newton (trainbfg), Bayesian Regulation (trainbr) introduced in Table 1 and its errors gave in Table 2, Table 3 gives the performance of the train with ePoch and time and Table 4 gives the weight and bias of the designer network.
|input||Analytic solution||Out of suggested ANN yt(x) for different training algorithm||Numerical solution|
Table 1: Analytic and Neural solution of Example 1.
|The error E(x) = | yt(x)-ya(x) | where yt(x) computed by the following training algorithm|
Table 2: Accuracy of solutions for Example 1.
|Train Function||Performance of train||EPoch||Time||Msereg.|
Table 3: The performance of the train with ePoch and time.
|Weights and bias for trainlm|
|Weights and bias for trainbfg|
|Weights and bias for trainbr|
Table 4: Initial weight and bias of the network for different training algorithm.
Consider the following 2nd order singular perturbed problem:
with BC: y(0)=0, y(1)=0 and x∈ (0, 1).
The analytic solution is :
The ANN trained using a grid of ten equidistant Points in 0, 1. Figure 2 display the analytic and neural solutions with different training algorithm. The neural results with different types of training algorithm such as: Levenberg–Marquardt (trainlm), quasi–Newton (trainbfg), Bayesian Regulation (trainbr) introduced in Table 5 and its errors gave in Tables 6 and 7 gives the performance of the train with ePoch and time and Table 8 gives the weight and bias of the designer network.
|In-put||Analytic solution||Out of suggested FFNN yt(x) for different training algorithm||Numerical Patching method|
Table 5: Analytic and Neural solution of Example 2.
|The error E(x) thit(x) -y,(x)||1 where yt(x) computed by the follming training algorithm|
Table 6: Accuracy of solutions for Example 2.
|Msereg.||Time||EPoch||Performance of train||Train Function|
Table 7: The performance of the train with ePoch and time.
|Weights and bias for trainhn|
|Weights and bias for trainhn|
|Weights and bias for trainhn|
Table 8: Initial weight and bias of the network for different training algorithm.
In this paper, we design neural network to solve singular perturbed problem. A fast and efficient algorithm (LM) for ANN with one hidden layer has been presented and tested on two examples. Through the comparison with exact solutions show that the ANN method has good accuracy and efficiency and results obtained using the ANN method is with low error.