Model Order Reduction Using Genetic Algorithm
Model Order Reduction Using Genetic Algorithm
Model Order Reduction Using Genetic Algorithm
Ahmed Adel
Alexandria University
Alexandria, Egypt
ahmedadel1993@student.alx.edu.eg
Khaled Salah
Mentor Graphics
Cairo, Egypt
Khaled_Mohamed@mentor.com
ABSTRACT
Model order reduction has been one of the most challenging
topics in the past years. Conventional mathematical methods have
been used to obtain a reduced order model of high order complex
models. In this paper, genetic algorithm (GA) which is one of the
artificial intelligence algorithms is used to approximate high-order
transfer functions (TFs) as lower-order TFs. Genetic algorithm is
considered as one of the evolutionary techniques which are used
in optimization problems. In this approach, genetic algorithm is
applied to model order reduction putting in consideration
improving accuracy and preserving the properties of the original
model which are two important issues for improving the
performance of simulation and computation and maintaining the
behavior of the original complex models being reduced. The
proposed technique could be used in EDA tools.
Keywords
Model reduction, crossover, selection, mutation,
algorithm, optimization, MSE, transfer function.
genetic
1. INTRODUCTION
In nature, biological structures that are more successful in
grappling with their environment survive and reproduce at a
higher rate[2]. Genetic algorithm is an artificial intelligence
algorithm based on natural selection which has been used in
optimization problems in the past years. Chromosomes are the
basic elements of genetics. Chromosomes of parents combine
with each other producing new individuals with new
chromosomes which have modified properties different from
properties of parents chromosomes. Successive generations of
individuals create a great probability for improvement over time.
Genetic algorithms were first developed by John Holland, his
colleagues, and his students at the University of Michigan for the
purpose of creating adaptive algorithms of natural systems and
designing system software based on them [1]. Holland introduced
the concept of genetic algorithm as a principle of Charles
Darwinian theory of evolution to natural biology [2]. Genetic
algorithm was proven to be a powerful technique in optimization
problems and it proves its efficiency mostly in problems requiring
efficient and effective search. The main reason for being powerful
in this area is although the algorithm is computationally simple
1.2 Contribution
In our work, we introduce a genetic algorithm model. The
proposed model is fully generic and able to be optimized. The
proposed model reduces the order of complex models of high
order transfer function to any desired order model. The model
targets high accuracy and preserving the basic properties of the
original model. The results of the proposed model simulations is
compared to powerful conventional mathematical methods like:
Pade Approximation which is based on moments matching, Routh
approximation, and Truncated balanced method in terms of
accuracy of step response, time required for reducing the model,
and frequency response. The proposed model proved efficiency in
2. Background
Genetic algorithm is based on iterative procedures of search for
an optimal solution for a problem which have multiple local
minima or maxima. The algorithm passes through steps of
recombination including crossover, and mutation then selection
which increase the probability of finding the most optimum
solution, which is the reduced model with the least error
compared to the original transfer model. The error is compared
with the original transfer function in terms of fitness function.
Before applying the genetic operators, a method of encoding
should be chosen to represent the data either in float form which is
the raw form of the data or binary representation or any other
representation.
The crossover operator is a process of creating new individuals
by selecting two or more parents and passing them through
crossover procedures producing one or two or more individuals.
Unlike real life there is no obligation to be abided by nature rules,
so the new individual can have more than two parents. There is
more than one method for crossover process like simple crossover
which includes exchange of genes between the two chromosomes
according to a specified crossover rate.
Arithmetic crossover occurs by choosing two or more
individuals randomly from the current generation and multiplying
them by one random in the case of the Inevitability of the
presence of a certain defined domain and more than one random if
there is no physical need for a defined search domain.
The second process of genetic operators is mutation which is a
process that occurs to prevent falling of all solutions of the
population into a local optimum of the problem. The third genetic
operator is selection which is the process of choosing a certain
number of individuals for the next generation to be recombined
generating new individuals aiming to find the most optimum
solution. There is more than one technique to select individuals
for the next generation.
The Elitism selection is simply selecting the fittest individuals
from the current population for the next population. This method
guarantees a high probability of getting closer to the most
optimum solution due to passing of the fittest chromosomes to the
crossover operator producing fitter individuals. Another technique
is Roulette wheel selection; in this kind of selection the
parenthood probability is directly proportional to the fitness of the
(1)
3. Proposed Method
Consider we have a linear system of nth order of q inputs and r
outputs described using state space in time domain as follows:
=
( ) +
( )=
( )
( )
(2)
(3)
( )
( )
(4)
( )
( )
(5)
Reduced TF
Iterations
MSE
GA MOR
Desired order
MODEL
Gen. Population
Time elapsed
step/freq. response
Apply arithmetic
cros sover using 2 random
numbers
4. Implementation
The model is implemented using genetic operators: crossover
and selection. The mutation operator wasnt used as it proves no
extra improvement for the model however it adds more processing
time. After generating an initial population from a selected
domain, crossover operator is applied. Arithmetic crossover was
used as an operator by generating two random numbers which add
extra randomness and extend the domain of search of the genetic
algorithm.
There was no physical need for a fixed domain of search for the
coefficients of the transfer function. Selection operator was based
on elitism where the fittest individuals in the current population
(with the least mean square error) are chosen and passed to the
next generation. This method guarantees that the best individuals
are passed to the next generation and at the same time it provides
a high degree of randomness unlike the other methods. Other
methods like roulette wheel selection and methods based on
giving weight to the individuals have a great disadvantage where
individuals with higher weight may only pass to the next
population which decrease randomness and decrease the
probability of generating individuals of different properties. The
chosen fitness function is the mean square error. The proposed
methodology is shown in Figure. 2.
A pseudo code for the implemented model that describes the
procedures of implementation is given as follows:
Listing 1
Pseudo code of GA
Input: problem, Population size, desired order.
Output: best_solution
Create initial population
Current population=initial population
For (i=1 to iterations)
For (j=1: to population size)
Choose 2 random parents (current population)
Do crossover producing 2 individuals
End for
Add best solution to best solutions set
Choose the best individuals of new generation
No
Iterations=max
iterations?
Yes
Compare the fittes t of al l
generations and choos e
the most o ptimum
solution
5. Testing
The simulations are carried on Matlab 2012b software, core i7
processor, 6 GB Ram memory. The reduced order model is
compared to the original one in terms of mean square error.
The testing procedures includes a plot of step response for the
reduced and the original TF, frequency response using Bode plot,
mean square error, and the time elapsed.
The time taken for most of the original transfer functions to be
reduced to lower order is nearly the same in all transfer functions
but, sometimes some functions need more time than the others.
The testing procedures was carried on more than 30 different
transfer functions. The transfer functions have different forms of
step responses and frequency responses where the step response
of some of them has an exponential form and other responses
were sinusoidal. The proposed model showed a great efficiency in
terms of accuracy and preserving the original models properties.
The first transfer function is a 10th order transfer function
expressed as follows:
.
G1(s)=
.
.
.
.
.
(6)
.
.
.
.
(7)
( )=
(9)
Step Response
1.2
original model
1
reduced model
0.8
Step Response
2.5
original model
Amplitude
0.6
reduced model
2
0.4
0.2
Amplitude
1.5
-0.2
-0.4
10
0.5
15
Time (seconds)
Time (seconds)
-20
-40
Magnitude (dB)
Magnitude(dB)
Bode Diagram
original model
reduced model
-60
90
-40
-60
0
45
Phase(deg)
-20
original model
Phase (deg)
0
-45
-90
-135
-2
10
-45
-135
10
10
-180
Frequency (rad/s)
-1
R1(s)
10
Iterations
MSE (%)
100
Time
0.07
10th
23 sec
G2(s)=
.
.
.
.
.
.
.
.
10
10
10
10
Frequency (rad/s)
10
Metric
reduced model
-90
(8)
Population
Iterations
MSE
(%)
Time
R2(s)
10
100
0.18
23 sec
The first comparison was made between the results of the reduced
order model of G1(s) and the mathematical methods. The step
responses are shown in Figure. 7 while the frequency responses
are shown in figure 8. The results of the proposed model and the
mathematical methods are compared in Table. 3.
Bode Diagram
Magnitude (dB)
20
0
original
-20
balanced
Pade
Routh
-40
reduced
-60
540
360
Phase (deg)
180
0
-180
-360
-2
10
-1
10
10
10
10
10
Frequency (rad/s)
Metric
GA
Pade
Routh
Balanced
0.07
1760
5.28
0.05
0.77
34
100
0.63
GA
Pade
Routh
Balanced
Freq. response
phase MSE
200
5*104
1.492*
103
1.2*105
0.07
1760
5.28
0.05
Time
23 sec
0.1 sec
0.1 sec
0.1 sec
0.50
0.33
23.8
0.19
Freq. response
phase MSE
34
291.86
1*103
3.7*103
Time
23 sec
0.1 sec
0.1 sec
0.1 sec
Step Response
2.5
Step Response
1.2
2
Amplitude
0.8
Amplitude
0.6
1.5
0.4
original transf er f unction
0.2
Blanced truncated
0.5
blanced truncated
Routh approximation
-0.2
reduced model
-0.4
Pade approximation
Routh approximation
Original model
10
15
Time (seconds)
Reduced model
0
Time (seconds)
Electrical,
Computer,
Energetic,
Electronic
Communication Engineering Vol:6, No:9, 2012.
Bode Diagram
20
M
agnitude(dB
)
[4]
[5]
-20
original
-40
balance
-60
Pade
Routh
-80
reduced
-100
45
0
P
hase(deg)
and
-45
-90
-135
-180
0
10
10
Frequency
10
(rad/s)
6. Conclusion
Model order reduction using genetic algorithm proved to be a
promising solution that can reduce complex high order models.
The proposed methodology showed efficiency in terms of
accuracy and preserving the basic properties of the original
complex model over conventional mathematical methods
increasing the opportunity for enhancing simulations and
increasing the accuracy of modeling.
REFERENCES
[1] J. Holland (1975), Adaptation in Natural and Artificial
Systems, University of Michigan Press.
[2] John R. Koza, Genetic Programming On the Programming of
Computers by Means of Natural Selection, 1992
Massachusetts Institute of Technology.
[3] S. Panda, J. S. Yadav, N. P. Patidar and C. Ardil,
Evolutionary Techniques for Model Order Reduction of
Large Scale Linear Systems, World Academy of Science,
Engineering and Technology International Journal of