Next Article in Journal
Spectrum Allocation Based on an Improved Gravitational Search Algorithm
Previous Article in Journal
Special Issue on Computational Intelligence and Nature-Inspired Algorithms for Real-World Data Analytics and Pattern Recognition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Evolutionary Algorithm for Designing Robust Analog Filters

1
Key Laboratory of Advanced Manufacturing Technology of Ministry of Education, Guizhou University, Guiyang 550025, China
2
School of Mechanical Engineering, Guizhou University, Guiyang 550025, China
3
Department of Computer Science and Engineering, University of South Carolina, Columbia, SC 29208, USA
*
Author to whom correspondence should be addressed.
Algorithms 2018, 11(3), 26; https://doi.org/10.3390/a11030026
Submission received: 12 November 2017 / Revised: 24 February 2018 / Accepted: 27 February 2018 / Published: 1 March 2018

Abstract

:
Designing robust circuits that withstand environmental perturbation and device degradation is critical for many applications. Traditional robust circuit design is mainly done by tuning parameters to improve system robustness. However, the topological structure of a system may set a limit on the robustness achievable through parameter tuning. This paper proposes a new evolutionary algorithm for robust design that exploits the open-ended topological search capability of genetic programming (GP) coupled with bond graph modeling. We applied our GP-based robust design (GPRD) algorithm to evolve robust lowpass and highpass analog filters. Compared with a traditional robust design approach based on a state-of-the-art real-parameter genetic algorithm (GA), our GPRD algorithm with a fitness criterion rewarding robustness, with respect to parameter perturbations, can evolve more robust filters than what was achieved through parameter tuning alone. We also find that inappropriate GA tuning may mislead the search process and that multiple-simulation and perturbed fitness evaluation methods for evolving robustness have complementary behaviors with no absolute advantage of one over the other.

Graphical Abstract

1. Introduction

Topologically open-ended computational design by genetic programming (GP) has been used for engineering design innovation, with many success stories in a variety of domains including analog circuits, digital circuits, molecular design, and mechatronic systems [1,2]. This approach uses genetic programming as an open-ended search method for functional design innovation—achieving a specified behavior without pre-specifying the design topology—and has achieved considerable success. However, in practical engineering design, there is another criterion in addition to functional specifications that should be considered during the design process. Robustness, as the ability of a system to maintain function even with changes in internal structure (including variations in parameters from nominal values) or external environment [3], is also critical to engineering design decisions. Engineering design, in reality, does not normally take into account every uncertainty or variation to which the engineered artifacts are subject to, such as manufacturing variation, degradation or non-uniformity of material properties, environmental changes, and changing operating conditions. We are interested in the following two kinds of system robustness: One is the system robustness with respect to perturbation in the parameters of the system, which has been well studied in the traditional robust design community [4] and also in evolutionary robust design [5]. The other kind is the system robustness with respect to topological perturbation—for example, accidental removal or failure of components. Reliable systems, having the least sensitivity of performance to variations or failures in the system components or environmental conditions, are highly desirable. However, there are relatively few studies that explore how a GP-based open-ended topology search may contribute to the design of robust systems such that they can withstand internal or external perturbations.
In this paper, we are interested in exploring GP algorithms for the design of robust systems such as analog filters. Specifically, we proposed a GP-based algorithm to facilitate the design of more robust systems with respect to parameter variations or uncertainty of the design variables. A set of systematic experiments are conducted to verify the following hypotheses about robust design using genetic programming: (1) that dynamic systems with high functional performance evolved by GP without considering a robustness criterion during the evolutionary process may have unacceptably low robustness with respect to parameter perturbation, (2) that the robustness of a system is strongly constrained by its topological/functional structure, and that the amount of robustness improvement available through parameter tuning is limited, and (3) that topologically open-ended synthesis by GP allows for the evolution of more robust solutions than traditional robust design approaches with parameter tuning.
Evolving robustness has a unique role in GP-based automated design. As most of such computational synthesis approaches depend on software simulators to evaluate the performance of candidate solutions, the quality of the final implementation of the evolved solutions strongly depends on how well the simulator approximates the physical behavior of the physical systems. A related issue arises when the objective function is replaced by a surrogate function (approximation) to speed up simulation [6]. One strategy to ensure effective transfer from evolved solutions to reality is to improve the robustness of the solutions by introducing noise or environmental variation during simulation or fitness evaluation, as proposed by Jakobi [7] in evolutionary robotics, for example. If the evolved solution is robust to such variations, the performance of the model evolved may sufficiently resemble (such that the model is useful) that of its physical implementation.
In natural evolution, the evolution of the functional structures of living organisms contributes significantly to their robustness with respect to the noisy environment. Functional characteristics such as modularity, redundancy, and negative feedback control are all typical features of living systems in response to survival requirements for robustness. Our hypothesis here is that control factors (design variables) as used in the robust design framework in [8] should not be limited to changing the dimensions (or sizing) and other parameters of the system, and that the system’s robustness should be considered when designing or evolving the functional or topological structure of the system. As any given function of a dynamic system can be implemented in various ways, we believe that the topological or functional design in the conceptual design phase may have a significant role in determining the ultimate robustness of the design solutions achievable through applications of traditional robust design procedures. There is already a body of research on how the structure of a system affects its functional robustness. Ferrer et al. [9] analyzed topological patterns in electric circuits and their relationship with the properties of the system behavior. Balling and Sobieszczanski-Sobieski [10,11], among others, discussed how the coupling structure of the system may affect robust parameter design. Robust design for truss topology has also been studied in [12]. However, a systematic methodology and investigation of robust design of dynamic systems based on a topologically open-ended search by genetic programming is still not available.
To examine the role of topology search in designing robust systems, design problems of two analog filters, i.e., a lowpass filter and a highpass filter, are synthesized using genetic programming. For each synthesis problem, three experiments are conducted: (a) evolutionary synthesis using GP without considering robustness criteria, (b) improvement in robustness of these evolved filters by tuning their parameters using a genetic algorithm (GA), and (c) evolution of robust filters (topological structure and parameters) using GP with a robustness criterion in the fitness function. These filter design problems are selected as they are perhaps the most popular problems in evolutionary synthesis research by either GA or GP [13,14,15].
The rest of the paper is organized as follows. Section 2 presents a short survey of applications of evolutionary algorithms in robust design. Section 3 introduces our GPRD (GP-based robust design) methodology, which exploits GP and bond graphs for automated synthesis of robust systems. Section 4 compares the experimental results of these approaches. Finally, conclusions and future research are highlighted in Section 5.

2. Related Work

Robust design [16] has been widely investigated in the engineering design community since the 1980s and remains an important topic [17]. In traditional robust design, one seeks to determine the parameter settings that produce desirable values of the mean (nominal) system performance, while at the same time minimizing the variance of the performance [16]. Many aspects of robust engineering design have been studied, such as the determination of the most robust operating points by system-wide parameter design [11], performance sensitivity distribution, and robust parameter design techniques [17]. Weng et al. [18] studied how simulation can be used to design photonic circuits that are robust to fabrication process variations. Sun et al. [19] proposed the use of norm metrics to quantify the robustness of both design parameters and performance uncertainty and then used a surrogating procedure to achieve robustness in a design space and to reduce uncertainty in performance space. Methods of representing and incorporating uncertainty in engineering design and optimization are reviewed in [4]. However, most of these approaches start from a design solution and only aim to determine its robust operating parameters with respect to different variations. The relationship between the topological or functional structure of a system and its robustness is often not treated. Especially, how robustness criteria should be incorporated into the conceptual functional design stage is not addressed. This is because the prevailing approach for engineering designs is a top-down procedure starting from the functional design to the detailed design and the robustness criterion is hard to evaluate without detailed design parameters, which are only available after the detailed design stage.
Evolutionary algorithms may provide a potential approach to address the limitations of traditional robust design. In the pioneering studies of evolutionary robotics, the Sussex group [7,20,21,22] developed a systematic methodology to evolve robust neural network controllers for robots. Husbands and Harvey [20] suggested an evolutionary approach to evolve neural network control systems for robots that are equipped with low-resolution sensors to work in uncertain environments. An incremental GA was used to evolve both the topology and connection weights of the controller. The discrepancy of simulation with reality was addressed by incorporating noise into the simulation process. They also argued that adaptive noise-tolerant networks are better building blocks than programs evolved by genetic programming for controller synthesis. Jakobi et al. [7] later showed that adding noise to simulations needs to be managed carefully. In their experiments, the correspondence between the robot behavior in simulation and that in the real world was demonstrated to be strongly correlated with the correspondence of the noise levels in the two cases. A more systematic approach to remove the gaps in transfer from simulation to reality is proposed in [23]. To ensure that the evolved robots do not exploit specifics of the simulation environments, Jakobi suggested making all implementation aspects unreliable by perturbing them randomly and in an irregular way [23]. That work suggested that it is better to resort to the evolutionary process to exploit intricate interactions and system dynamics if functional robustness is to be achieved.
Application of evolutionary algorithms to traditional parametric robust design has been attracting increasing attention in the past decade [24,25,26]. Tsutsui et al. [24] proposed a GA-based robust solution-searching scheme ( R S 3 ) to evolve robust solutions. This approach works by adding perturbation noise to the design variables before fitness evaluation. It thus does not rely on multiple re-sampling to evaluate the expected fitness, as was used by Wiesmann et al. [25]. In Wiesmann et al.’s approach, each individual is simulated t times to estimate its expected loss function (fitness) defined as 1 t · i = 1 t ( τ f ˜ ( x , δ i ) ) 2 , where δ i , i = 1 , , t are normally distributed noise with zero mean and standard deviation σ ; τ is the ideal target value. Their experiments showed that multilayer optical coating designs found by the evolutionary algorithm were substantially more robust to parameter variations than the reference design, but usually at the cost of reduced performance in undisturbed situations. This observation motivated the later work of using an evolutionary multi-objective approach to figure out the trade-off map between robustness and optimal functional performance [26,27,28]. Forouraghi [26] introduced an interval computation method to avoid artificial insertion of Gaussian noise to parameter variables in order to build tolerance against internal or external perturbations. In this approach, each design variable is represented as an interval. The robustness of a candidate solution is defined as the signal-to-noise ratio S/N, which considers both target performance and variations in face of perturbations. As a minimization criterion, they used
S / N = η s = 10 log ( i = 1 n y i 2 n )
and for maximization criterion,
S / N = η l = 10 log ( i = 1 n 1 / y i 2 n )
where n is the number of simulations, and y i is the performance of the i-th simulation. This robustness measure ensures that either improvement in the mean result or a reduction in variability will lead to greater S/N values. They then used a multi-objective GA to evolve a solution with multiple robustness criteria using a design-of-experiments (DOE) method with orthogonal arrays to evaluate the fitness. In a more explicit way of addressing the balance of performance and robustness, Ray [27] expressed the robust design problem as a three-criterion multi-objective problem, simultaneously optimizing an individual’s performance without perturbation, the mean performance of its neighbors resulting from perturbations, and the standard deviation of its neighbors’ performances. The perturbations in their experiments are simply modeled by deviations of ± 1 % beyond the variable values. They also observed a degradation in performance of the evolved robust solutions. In the most existing evolutionary robust design, multiple simulations are needed to evaluate the robustness of a single individual, which severely increases the computational burden. To address this problem, Jin et al. [28] used two methods for estimating the robustness measures of an individual—by exploiting its neighbor individuals in the current population as first proposed by Branke [29] and by using all evaluated individuals. They then evolved robust solutions with both performance and robustness as criteria in a two-criterion optimization as suggested by Das [30]. Jin’s robustness estimation approach can greatly reduce the number of function evaluations, when it is applicable. However, it is difficult to apply this method for evolving robust designs with variable structures as in topologically open-ended automated synthesis using GP. In these cases, it is difficult to define a neighborhood for a given individual, since most other individuals have different numbers and types of components.
A related research area to robust design focuses on evolving robustness in a noisy environment [31,32,33]. In these problems, the variation in the environment leads to uncertainty in the fitness function evaluation and the true fitness of a candidate solution needs to be evaluated based on sampling in multiple environments. In the evolutionary robotics area, Lee et al. [34] and Jakobi [23] evolved robust mobile robot controllers by training them in multiple trials of simulation, using genetic programming and a genetic algorithm, respectively. The active area of evolving robust systems is evolvable hardware [35,36,37]. Most of these approaches employ genetic algorithms or evolution strategies as the search procedures. Very recent work is the evolution of robust digital circuits [38,39]. In this work, Miller, Hartmann, and their collaborators examine the feasibility of evolving robust digital circuits using a type of “messy gate”. Hartmann et al. [40] investigated how evolution may exploit non-perfect digital gates to achieve fault tolerance, including tolerance to output noise and gate failure. However, the noise introduced to improve robustness is not applied to parametric values of the components, but to the analog outputs of the messy gates, and an evolution strategy (ES) is used as the open-ended topology search tool. This method is thus not as instructive as might be desired in exploring effects of alternative topologies.
In addition to robust design and optimization in noisy environments, evolutionary algorithms have also been successfully applied to other optimization problems with uncertainty such as time-varying optimization problems. For a complete overview on the general topic of evolutionary optimization under uncertainty including the major techniques and issues, readers are referred to an excellent survey of Jin and Branke [28].
In this paper, we chose analog filter design as our test problems. The pioneering work of Koza in automated analog circuit synthesis, including lowpass, highpass, and asymmetric bandpass filters, is described in [13,41]. However, they did not specifically work on evolving robust circuits. Lohn and Colombano [15] proposed a linear representation approach to evolve analog circuits in which several lowpass filters were used as test problems. In another work, by Goh and Li [42], a genetic algorithm with a linear direct encoding approach is used to evolve lowpass and bandpass analog filters. What is special in their work is that they constrained the component values to be chosen among the manufacturer’s preferred component values. Finally, we have applied GPBG approach to a few benchmark and real-world problems including the lowpass filter design problem [43] and pump design problem [2].
Much work has been done in both evolutionary robust design and analog circuit synthesis. However, few studies specifically address how GP-based topologically open-ended synthesis may provide a new method of open-ended robust design. This may enable us to move robust design forward to the conceptual/functional design stage, i.e., at the very beginning, which will augment the current practice of designing for robustness in the parametric design stage.

3. Methods: A Robust Analog Filter Design Using Evolutionary Algorithms

3.1. Analog Filter Design using Bond Graphs and GP

In this work, we used bond graph modeling [44] as the representation scheme for analaog filter design. A bond graph is a graphical representation of a physical dynamic system, which allows the conversion of the system into a state-space representation. This capability of seamlessly incorporating multiple domains makes our results widely applicable in the evolution/design of physical systems. Here, we briefly introduce the bond graph and show how GP can be used to evolve analog filters based on bond graph representation with details explained in our previous work [2,45,46].

3.1.1. Bond Graphs

The bond graph is a multi-domain modeling tool for the analysis and design of dynamic systems, especially hybrid multi-domain systems including mechanical, electrical, pneumatic, and hydraulic components [44]. One advantage of using bond graphs for open-ended design exploration is that the complex loops typical in electric circuit schematics can be transformed into tree-like structures by the bond graph’s 1-junction (serial connection) and 0-junction (parallel connection) concepts, which tend to be easier to evolve in general. Another advantage is that the multi-domain nature of bond graph modeling facilitates the evolution of mechatronic systems. Many researchers have explored the bond graph as a tool for system design [47]. Details of notation and methods of system analysis related to the bond graph representation can be found in [44]. Figure 1 illustrates a small bond graph that represents the accompanying electrical system. Figure 2 shows the complex bond graph model of a lowpass filter. A typical simple bond graph model is composed of inductors (I), resistors (R), capacitors (C), transformers (TF), gyrators (GY), 0-Junctions (J0), 1-junctions (J1), sources of effort (SE), and sources of flow (SF). In this paper, we are only concerned with linear dynamic systems, or more specifically, analog filters as a case study, represented as bond graphs, which are composed of I/R/C/SE(as input signals)/SF(as output signal access points). Due to the limitation of our bond graph simulation, we did not include transistors into our model so the system will have difficulty evolving active analog filters. They can be evolved using SPICE simulation and modeling software, which is beyond the scope of this paper.

3.1.2. Evolving Analog Filters Using Bond Graphs and GP: The GPBG Framework

The problem of automated synthesis of bond graphs involves two basic searches: the search for a good topology and the search for good parameters for each topology, in order to be able to evaluate its performance. Building upon Koza’s work on automated synthesis of electronic circuits [41], we created a developmental GP system for open-ended synthesis of mechatronic systems represented as bond graphs [2]. It includes the following major components: (1) an embryo bond graph with modifiable sites at which further topological operations can be applied to grow the embryo into a functional system, (2) a GP function set, composed of a set of topology manipulation and other primitive instructions which will be assembled into a GP tree by the evolutionary process (the execution of this GP program leads to topological and parametric manipulation of the developing embryo bond graph), and (3) a fitness function to evaluate the performance of candidate solutions.
Choosing a good function set for bond graph synthesis is not easy. In our earliest work [43], a basic GP function set was used for evolutionary synthesis of analog filters. In that approach, the GP functions for topological operation included {Insert_J0/J1, Add_C/I/R, and Replace_C/I/R}, which allowed evolution of a large variety of bond graph topologies. The shortcoming of this approach is that it tended to evolve redundant and sometimes causally ill-posed bond graphs [48]. Later, we used a causally well-posed modular GP function set to evolve more concise bond graphs with much less redundancy [45]. However, that encoding had a strong bias toward a chain-type topology and thus may have limited the scope of topology search [46]. In this paper, we used the following hybrid function set approach to reduce redundancy while enjoying the flexibility of topological exploration.
F = {Insert_J0E, Insert_J1E, Add_C/I/R, EndNode, EndBond, ERC}.
Details of this GP function set for bond-graph-based analog filter synthesis are available in [46,49]. Figure 3 shows a GP tree that develops an embryo bond graph into a complete bond graph solution. Our comparison experiments [46] showed that this function set is more effective on both an eigenvalue and an analog filter test problem, so the new set was used in this paper.
In this study, two types of analog filters including lowpass and highpass filters were synthesized. In such GPBG-based filter design problems [43], a bond-graph-represented analog filter composed of capacitors, resistors, and inductors is evolved such that the magnitude of its frequency response approximates a specified filter frequency response specification. An embryo bond graph and its equivalent circuit are illustrated in Figure 1 and are used in all experiments. Note that the 0-junction is the initial modifiable site, where further topological developments can proceed as instructed by a GP program tree. The voltage at this 0-junction is the output signal.
Instead of using the sophisticated SPICE simulation program as used in the existing analog filter synthesis [13,50], the frequency response of a bond graph can be simulated in a more convenient way: first, a state equation of a bond graph is derived automatically from the model, which generates the A, B, C, and D matrices [44]. These state space models can then be simulated on Linux PCs using C++ simulation code generated from Matlab Compiler 3.0.
The detailed specifications of the three analog filter synthesis problems are as follows:
  • The lowpass filter synthesis problem is extracted from [13], in which the frequency response performance of a candidate filter is defined as the weighted sum of deviations from an ideal frequency response magnitude over 101 points:
    F k o z a ( t ) = i = 0 100 [ W ( d ( f i ) , f i ) × d ( f i ) ]
    where f i is the sampling frequency, d ( x ) is the absolute deviation of candidate frequency response from target response at frequency x, and W(x,y) is the weight function specifying the penalty level for a given frequency response at a specified frequency range. The sampling points range from 1 Hz to 100 KHz, evenly distributed on a logarithmic scale. If the deviation from the ideal magnitude is less than 0.03 V, the weight W is 1. If the deviation is more than 0.03 V, the weight W is 10. The pass band is [1, 1K] Hz, and the stop band is [2K, 10K] Hz. A “don’t care” band between 1 and 2 KHz neglects any deviation from the target response there.
  • The highpass filter synthesis problem has a similar configuration to the lowpass filter except for the complementary definitions of the pass and stop bands. The pass band is now defined as [2K, 10K] Hz, while the stop band is [1, 1K] Hz.
To evolve an analog filter without considering robustness, the fitness function is defined as follows. First, we calculate the raw fitness defined as the average absolute deviation between the frequency response magnitude of the candidate solution and the target frequency response over all 101 sampling frequencies:
f r a w = 1 101 · F k o z a ( t ) = 1 101 · i = 0 100 [ W ( d ( f i ) , f i ) × d ( f i ) ] .
Note that this f r a w definition differs from Koza’s raw fitness definition in Equation (1) by a multiplier equal to the number of sampling frequencies. We use the average deviation rather than sum of deviations to remove the influence of the number of sampling points.
The raw fitness of a candidate solution is then normalized as follows:
f n o r m = N O R M N O R M + f r a w .
where NORM is a normalization term aiming to adjust the f n o r m into the range of [0, 1]. This transformation also serves to convert the minimization of deviation from the target frequency response into the maximizing of fitness process used by our GP system. Since tournament selection is used as the selection operator, the normalization term can be arbitrary positive numbers. Here, the tournament selection operator is used due to its capability of achieving a good balance between exploiting and exploration and is widely used in GP [51]. For the lowpass and highpass filter problems here, NORM is set to 10, which gives a fitness range of [0, 1].
In this paper, a standard strongly-typed multi-population generational GP enhanced with the above features is used to evolve analog filters represented as bond graphs. The running parameters are specified in Section 4.

3.2. Evolving Robust Analog Filters With Respect to Parameter Variation: The Unified Approach

In this paper, we examine how to evolve robust analog filters with respect to parameter variations. There are many ways to achieve this by combining GP for an open-ended topological search and GA or differential evolution (DE) for fine-tuning the evolved solutions from GP. In this section, first the common characteristics of the various methods for evolving robustness are presented, followed by details of each individual approach.
The typical approach for evolving robust designs [32] is to use multiple Monte Carlo samplings with different environmental or system configurations (e.g., perturbation of parameter values of the system) to calculate a worst-case or an average fitness for a given candidate solution, as shown in Equation (3). This robust-by-multiple-simulation (RMS) method is used in [25]. Another method is to simply add a perturbation to the design variable before evaluation. This perturbation, however, is not incorporated into the genome, making it different from a “normal” parameter mutation operator or Larmarckian-style evolution algorithms. This robust-by-perturbed-evaluation (RPE) method is used in [24] and is suggested to be more efficient by Jin et al. [28]. Both methods are tested in this work. For the RPE method, no special fitness function is needed. One only needs to perturb the parameters before normal fitness evaluation. For the RMS multi-simulation method, the raw fitness for a design solution with robustness criterion is defined as the sum of a number SPI of 101-point deviations from the target frequency response curve:
f r o b u s t r a w = k = 1 S P I f r a w k
where SPI is the number of Monte Carlo sampling evaluations for each individual, and f r a w k is the raw fitness of the k-th sampled evaluation with a different Monte Carlo perturbation of the parameters as defined in Equation (2). With this raw robustness fitness, we then calculate the final fitness according to Equation (3).
In the experiments reported below, the perturbation of the component values during evolution is implemented by adding to each component parameter Gaussian noise N ( μ , σ ) with a mean μ of 0 and the standard deviation σ set to 20% of the parameter value. This perturbation model has been widely used by previous researchers [24,25]; while it may not be appropriate for all manufacturing processes, it is sufficient for our purposes as an approximation to the real component value degradation model in some situations. When more sophisticated models are used in the future to capture all parameter variations due to process, temperature, and voltage changes and simulate them in our system design evolution, the evolved systems can be more robust to these variations in reality. We also need to consider the parasitic effects existing in the evolved circuits. One difficulty with this definition is that if the original parameter value is zero, then no perturbation will be generated. Although this is rare in evolutionary experiments, we have incorporated a special step to check whether the component value is zero. If that is the case, we set the standard deviation for the perturbation to 1.

3.3. Genetic Algorithms for Robust Analog Filter Design: The Traditional Robust Design

In this section, we follow the practice of traditional robust design methodology to improve system robustness by tuning its parameters. We first use the improved GPBG approach as described in Section 3.1.2 to evolve a high-performance analog filter, without incorporating any robustness criterion in the fitness function. We shall then apply to the evolved filter a state-of-the-art real parameter genetic algorithm—the G3PCX-GA proposed by Deb [52]—to tune its parameters so as to improve its robustness with respect to parameter perturbation, while its functional structure remains unchanged. Since parameter search for dynamic systems is a complex highly multi-modal search problem, we choose G3PCX-GA to approximate the highest robustness achievable through parameter tuning within a given budget of functional evaluations. When the multi-simulation method (RMS) to evolve robustness is used, we call this hybrid GP/GA method the GPGARMS method. When the robust-by-perturbed-evaluation (RPE) method is used to evolve robustness, we call the hybrid GP/GA method the GPGARPE method, which, however, is not evaluated in this study.
In the minimization by G3PCX-GA, we used the robust raw fitness defined in Equation (4) as the final fitness of an individual. We believe this fitness measure is better than the average of normalized final fitness of each sampling evaluation for its lower distortion of the optimization objective values.

3.4. GPRD: Robust Analog Filter Design Using GP

In natural evolution, many structural patterns, such as modularity, emerge naturally as responses to a robustness requirement for survival. In a similar way, open-ended evolutionary synthesis by GP also provides such a capability for the evolution of robust systems. Unlike traditional engineering design, in this framework, topology search and parameter tuning are conducted simultaneously, so there is no distinct boundary between the conceptual design stage and the detailed design stage.
Two GP-based approaches for synthesizing robust designs are tested in this paper. The first is the GP with robustness-by-multi-simulation (GPRMS), which uses multiple simulation to estimate the robustness fitness of a candidate solution. The second is GP with robustness-by-perturbed-evaluation (GPRPE), which temporarily adds perturbation noise to the design variables before evaluating its fitness.
To be consistent with the experiments in Section 3.1.2, we use the generational GP, which is also employed by the experiments of Section 3.1.2. The final fitness of an individual, calculated from the sampling fitnesses, is the same as that defined in Equation (2), where f i k is defined as F i t n e s s n o r m in Equation (3).

3.5. Evaluation Criteria

For all evolved filters, we used the following criteria to evaluate their robustness and functional performance. For each candidate filter, we run SPI = 10,000 simulations, each with different perturbations of its component parameters. This SPI number is chosen so that the statistical values from a set of simulations remain stable. The perturbation is modeled as Gaussian noise N ( μ , σ ) with the mean μ set to 0 and the standard deviation σ set to 10–50% of parameter values, in steps of 10%.
There are two ways to calculate the robustness of a candidate design. One approach is the definition of robustness fitness as used by Forouraghi [26] and defined in Equation (4). This Type-I robustness f r o b u s t I definition uses the average (or sum) of the total 101-point deviation from the target frequency response curve for SPI simulations as the final robustness. The benefit of this robustness definition is that it incorporates requirements on both performance and its variation: increasing mean performance and reducing performance variation both lead to increased robustness. This measure is especially good for evolving robustness during which the absolute performance is also important in addition to the stability of functions. The disadvantage of this definition is that it cannot distinguish the robustness if a specified level of performance is satisfied by all solutions being compared. For such a case, it is more desirable to use the standard deviation of the 101-point deviations from the target frequency response curve for all 10,000 simulations as the robustness criterion. We call this robustness definition Type-II robustness. Note that this is not the other form of robustness (to component failure) but rather another definition of robustness to parameter variation. The ideal case is to use a multi-objective method to separate performance and its variation completely, as suggested by Das [30] and Jin et al. [28]. We defined this Type-II robustness criterion as in Equation (5). In this paper, both criteria will be evaluated.
f r o b u s t I I = k = 1 S P I ( f r a w k f ¯ r a w ) 2 S P I
where f ¯ r a w is the average 101-point-deviation from the target frequency response curve; f r a w is the 101-point-deviation of a single simulation from target frequency response curve; SPI is the number of simulations used to evaluate the fitness of an individual.

4. Experiments and Results

In this section, the following four different methods for evolving robust or non-robust lowpass and highpass filters are presented:
  • standard genetic programming (GP) without considering robustness requirements;
  • genetic programming with robustness-by-multi-simulation (GPRMS);
  • genetic programming with robustness-by-perturbed-evaluation (GPRPE);
  • hybrid GP/GA robust design method (GPGARMS).
For all these methods, a fixed number of function evaluations is allocated to ensure fairness of comparison. Here, one function evaluation means one simulation of a design for measuring its performance, which is different from one fitness evaluation for an individual. In the multi-simulation robust design method GPRMS, multiple simulations (SPI = 10) are used to evaluate the fitness of a single design. For the lowpass and highpass filter design problem, the computation budget is 1,000,000 function evaluations. Note that, for methods that use multiple simulations to estimate the robustness fitness, each simulation is counted as one function evaluation. In addition, for the hybrid GPGARMS method, we allocate 500,000 function evaluations for GP evolution and the remaining 500,000 for GA evolution for robustness.
We used only one parameter perturbation model in the evolutionary synthesis experiments—Gaussian noise N ( μ , σ ) with a mean μ of 0 and the standard deviation σ set to 20% of the parameter value. However, the evaluation of the evolved filters are subject to multiple perturbation magnitudes with extensive simulation.
To assess the statistical significance of the performance discrepancy between these methods, for each target filter type and each synthesis method, 15 runs were conducted. The size of these experiments was determined by the computing resources available. However, since we found that the results are quite stable across multiple runs, the level of replication appears to be sufficient for the purposes of this paper.
All experiments described below used the same embryo bond graph shown in Figure 1. The component values of source resistor R s and load resistor R l o a d are both 1 Ω for lowpass and highpass filter synthesis.
In the following sections, we first describe separately the experimental configuration of each method, the best evolved bond graph models of each filter type, and the frequency response of the best solutions from each method. These results provide some general ideas regarding how robustness is evolved with respect to the parameter perturbations. Then, in Section 4.4, a statistical comparison of the performances of five algorithms for evolving robust filters are presented. These analyses will be used to support the hypothesis outlined in the introduction section.

4.1. Evolving Analog Filters Using GP without Considering Robustness in the Fitness Function

In this experiment, 15 analog lowpass and highpass filters are evolved using standard GP without incorporating a robustness criterion in the fitness function Equation (3). The following common running parameters (Table 1) are used throughout all GP experiments in this paper:
Note that the maximum tree size here is 8 rather than 17, used in many other works. This parameter was selected by considering the reduced tree sizes due to our simplified parameter representation method and available computational resources. The maximum tree depth of 8 allows for the synthesis of analog filters with up to 100 components, which is sufficient for our purpose.
We select the evolved filter with the highest performance to test its noise tolerance over the degradation or variation in the component parameters with different perturbation magnitudes. As described above, the evaluation of robustness with respect to parameter perturbation is conducted by running 10,000 simulations of the configurations of the Gaussian parameter perturbations.
Figure 2 and Figure 4 show the topology of the evolved lowpass and highpass filters with the highest performance out of 15 runs. The frequency responses of these filters are summarized in Figure 5 along with the results of evolved robust filters. The evolved best lowpass and highpass filters have 39 and 27 components, respectively. The lowpass and highpass filters approximate the ideal frequency response closely, with the sum of deviations over 101 points being only 6.43 and 0.32, respectively.
Figure 5a,b show the frequency responses of the four best lowpass and four highpass filters evolved using GP, GPGARMS, GPRMS, and GPRPE. The performances of the evolved filters in an ideal condition without parametric perturbation were compared. Here, the “best” means the highest fitness of each approach out of 15 runs. Only standard GP has no consideration of robustness. The other three methods considered robustness in their fitness function definition.
From Figure 5, first, one can see that the frequency responses of the filters evolved by GP, GPRMS, and GPRPE are very similar. That means introducing robustness requirements does not necessarily decrease performance with nominal parameters significantly. The only exception is GPGARMS, which distinguishes itself from the others with much worse nominal performance (larger deviations). After checking the standard deviation of the sum of deviations of 10,000 runs with 20%-level parametric perturbation, it appears that the nominal performance is sacrified by the GA for better robustness. One can also observe that the filters with highest performance for both filter types are evolved using GPRPE. This can be explained by the actual number of genetic topological and parametric operations. Here, standard GP only ran 500,000 fitness/function evaluations during which genetic operations were conducted (we did not run 1,000,000 function evaluations for GP because standard GP always converges to local optima long before running out of 500,000 evaluations). GPRMS actually only ran 100,000 fitness evaluations (with 1,000,000 function evaluations) during which genetic operations were conducted, because 10 simulations were used for each individual during fitness evaluation, while genetic information was kept unchanged. The GPRPE, however, ran 1,000,000 fitness (function) evaluations during which genetic operations were conducted.

4.2. Evolving Robust Analog Filters Using Genetic Algorithms: The Classical Robust Design

In this experiment, the G3PCX-GA [52] is used to improve the robustness of the best filters evolved using standard GP through parameter tuning while keeping functional structure unchanged. As we can see from Figure 2 and Figure 4, these two evolved filters are very complex, with 39 and 27 parameters to search. As the objective function of this optimization is highly multi-modal, this is hard optimization even for G3PCX-GA, as the experiment demonstrates. The running parameters for this experiment are summarized in Table 2.
Figure 6a,b show the frequency response distributions of the filters and their perturbed variations evolved by standard GP and fine-tuned by G3PCX-GA. These perturbed filters are obtained by adding to the component values of the original solution a Gaussian parameter perturbation with a magnitude subject to Gaussian distribution N ( 0 , 0.2 ) times the original component values. The GA-based fine tuning of the filter improves its performance, which is indicated by its robustness fitness, as defined in Equation (5): the filter originally evolved by GP has a f r o b u s t I I of 7.21, while the filter then improved by GA has a f r o b u s t I I of only 1.15. More systematic comparison of their performance and robustness is described in Section 4.4.
An interesting observation about G3PCX-GA’s result is that the increase in Type II robustness is accompanied by a decrease in Type I robustness, which is impacted by both performance and function stability. This phenomenon can be easily seen in Figure 7. One explanation of this severe degradation of function is that GP, where robustness is not considered, has evolved a high-fitness filter that is, however, highly sensitive to parameter perturbations. In this case, parameter tuning by G3PCX-GA can easily disrupt the highly coupled parameter relationship of these components’ values and thus decrease the functional performance and improve robustness. Another possible factor is that this parameter optimization problem with 65 variables is too difficult for G3PCX-GA to optimize well within the allocated computational effort. We also found that the best fitness of the population stagnates after about 200,000 evaluations, which may be because the best fitness has approximated or reached the limit of the best possible fitness with the given topological structure. The population may also stagnate due to premature convergence.

4.3. Evolving Robust Analog Filters Using Genetic Programming: Open-Ended Topology Innovation for Robust Design

In the following experiment, we try to evolve robust analog filters that have a higher tolerance of the variation in component values and have graceful performance degradation. Two GP algorithms for evolving robust solutions are used: GPRMS and GPRPE. The configurations for this experiment are the same as those used in Section 4.1.
Figure 8 shows the best lowpass and highpass filters evolved using GPRMS with robustness requirements. Compared with the filters evolved using standard GP in Figure 2 and Figure 4, these two filters use fewer components, while the functional performance remains similar. The robustness of these two filters is compared to that of filters evolved using GP and GPGARMS in Figure 6 for the highpass filter problem. As can be seen, both GP methods with robustness requirements (GPRMS and GPRPE) evolved filters with higher robustness, as shown by their much less dispersed frequency response distributions. Statistical evidence is presented in Section 4.4. Compared to GPGARMS, one distinguishing characteristic of GPRMS and GPRPE is that they evolved more robust filters without significant degradation of functional performance. This observation can also be made in Figure 7, which shows that GPGARMS has the lowest performance variation (standard deviation) across all perturbation levels but also has the highest average-101-point deviation from the target frequency response curve.
From Figure 7a, however, we find that filters evolved by both GPRMS and GPRPE, in a trade-off for higher robustness, have an average 101-point deviation higher than that evolved by standard GP for all perturbation levels less than 0.3. However, when the perturbation level is greater than 0.3, GPRMS achieves an average-101-point deviation that is lower than that achieved by GP, and GPRPE’s performance remains worse. This means that GP with the multiple-simulation method, compared with the perturbed fitness evaluation method, evolved designs with a higher tolerance to parametric perturbation. From Figure 7b, we find that GPGARMS shows a clear advantage in terms of f r o b u s t I I , which is, however, achieved by sacrificing functional performance. In other words, the performance is more stable with regard to its average frequency response. Nevertheless, the total deviation from the target response becomes larger.
To further compare the characteristics of the evaluated algorithms, we did an extensive analysis of the performances of all filters evolved by the five algorithms, and this analysis is summarized in Table 3. For each algorithm, 15 runs were performed to evolve 15 lowpass filters. Each filter was then subjected to 10,000 simulations with random perturbations of their component parameters. We calculated the average 101-point deviation ( f r o b u s t I ) and its standard deviation ( f r o b u s t I I ) from the target frequency response, for each filter. The upper half of the table shows the average of f r o b u s t I and its standard error over 15 runs for each algorithm under different perturbation environment. The lower half shows the average of f r o b u s t I I and its standard error over 15 runs for each algorithm under different perturbation environments. The perturbation magnitudes are represented by the magnitudes of the Gaussian perturbation as a percentage of component parameter values.
In the top half of the table, we can see that the GPRPE method evolved filters with the highest average performance or f r o b u s t I . These filters are not only best at a perturbation level of 0.2, which is applied during the evolution process, but also best across all other perturbation levels tested in the “evaluation” phase. This implies that ideal analog filters have continuous performance landscapes with satisfactory degradation properties. What is surprising here is that not only are the f r o b u s t I values of the filters evolved by GPRPE the best but also their standard error values are lowest at all perturbation levels. This means that the GPRPE has the maximum number of fitness evaluations; thus, the genetic operations algorithm evolves high-performance filters with the most reliability.
Interestingly, there is a complementary relationship between the performances of GPRPE and GPRMS. Examining the lower half of Table 3, the GPRMS turns out to be the best algorithm to evolve robust filters with the highest robustness in terms of evaluation criterion f r o b u s t I I , with minimum variation in the frequency response, and this advantage holds true for all perturbation levels. Comparing GPRPE and GPRMS, it appears that GPRMS with its multi-simulation method can evolve more stable solutions due to its superior estimation of the true robustness property of a candidate solution. GPRPE, instead, suffers from its robustness estimation error but enjoys a higher amount of genetic operations. These results show that there is no simple conclusion about which method (of GPRPE and GPRMS) is better than the other. The choice depends on the design objectives.
Another interesting observation is the performance of GPGARMS. For highpass filters, we showed that GPGARMS evolved filters with the highest type II robustness f r o b u s t I I . However, it is not the case for the lowpass filter problem. The GPGARMS shows the worst results in terms of both Type I and Type II robustness measures. How could GA produce worse results than standard GP? We have two possible explanations. First, the fitness function for GPGARMS only uses 10 simulations to estimate the robustness. This may result in large estimation error, thus misleading the evolution. Since the robustness scores in Table 3 are calculated by 10,000 perturbed simulations, results produced by a misled GPGARMS could be even worse than its starting solutions. Another projection is that the parameters of lowpass filters are located in a much more rugged landscape, so a misled GPGARMS tends to drift away quickly from the highly coupled parameters of the starting solution evolved by GP. These explanations are under investigation and will be reported elsewhere.
From these analyses, we find that, compared to evolutionary design where robustness is not considered, allowing topological variation during the evolution of robust systems can yield more robustness and a superior structure with respect to perturbations of the component parameters of the filters, through topological innovation. To further consolidate our observations, in the following section, we report a series of statistical tests run to verify the significance of the differences seen between these algorithms.

4.4. Statistical Results of the Three Methods for Evolving Robust Filters

4.4.1. GA-RMS Improves the Robustness of Standard GP Results

For the highpass filter problem, we did a t-test to compare the robustness of the evolved solutions by GPGARMS and standard GP in terms of variation (Equation (5)) at the 0.2 perturbation level. A significance level of p 0.001 is achieved, strongly indicating that GPGARMS improved the robustness of the filters evolved by standard GP. However, as shown in Figure 7, this improvement is at the cost of degraded performance. It is impressive that, in this highpass problem, GPGARMS achieves the highest robustness in terms of standard deviation for all 5 perturbation levels. However, as shown in Table 3, the GA could also potentially generate statistically significantly worse solutions or ruin the results evolved by standard GP. For the lowpass filter problem, we found that there is no significant difference in Type-I robustness between GPGARMS and GP. GA also decreased the Type-II robustness in this problem.

4.4.2. Topological Innovation Using GP-Evolved Filters with Higher Robustness than Parametric Robust Design Using GA

For the highpass problem, we applied the t-test to compare GPGARMS and GPRMS. The 95-percent confidence interval for differences in the means of Type-I robustness f r o b u s t I ranges from 51.617 to 39.841 , showing that GPGARMS degraded Type-I robustness. The difference in the mean values of the two groups is greater than would be expected by chance ( p 0.001 ). For GPGARMS and GPRPE, the difference is even more significant. The 95-percent confidence interval for difference of means ranges from 57.482 to 44.846 with p 0.001 . In terms of f r o b u s t I I , the 95-percent confidence interval for difference of means between GPRMS and GPGARMS ranges from 9.920 to 2.538 , while there is no significant difference between GPRPE and GPGARMS.
For the lowpass problem, the 95-percent confidence interval for the difference in the means of Type-I robustness f r o b u s t I between GPRPE and GPGARMS ranges from 51.617 to 39.841 with p 0.001 . There is no significant difference between GPRMS and GPGARMS in terms of f r o b u s t I . In terms of f r o b u s t I I , the 95 percent confidence interval for the difference in the means between GPRMS and GPGARMS ranges from 9.920 to 2.538 , while there is no significant difference between GPRPE and GPGARMS.

4.4.3. GP with Robustness Requirements Constrains Bloating

One critical issue in GP-based evolutionary synthesis is the bloating problem, or the growing sizes of the solutions produced during the evolution process. To evaluate how a robustness requirement influences this phenomenon, we compared the number of components of the filters evolved by methods discussed in the beginning of this section. To date, we have only made 10 runs for each algorithm. The results are shown in Table 4. The standard GP runs with 100,000 and 500,000 function evaluations with no robustness requirements are compared with those of two GP methods with robustness requirements. Note that, for the GPRMS method, as each fitness evaluation of an individual costs 10 function evaluations, the actual number of genetic operations is only one-tenth of the number of such operations performed in the 1,000,000 function evaluations of GPRPE.
First, we found that the average sizes of the design solutions evolved by GP (in 100,000 function evaluations) and GPRSM (in 1,000,000) are similar. As both methods used the same number of genetic operations, it seems that the sizes of the evolved filters are determined more by the actual number of genetic operations rather than the number of fitness evaluations. Since the number of evolved solutions by GP and GPRSM are too small to satisfy the normality condition, a Mann–Whitney rank sum test was used to show that any difference in these two groups of sizes is not statistically significant, at the level of p = 0.880. Second, we compared the solution sizes of GP (100,000) with GP (500,000) using the same Mann–Whitney rank sum test. A statistically significant difference exists at the level of p < 0.001 . This demonstrates that the solution sizes indeed grow with an increasing number of fitness evaluations. A third Mann–Whitney rank sum test was applied to show the size difference between the solutions evolved using GP (500,000) and GPRPE (1,000,000). There was no statistically significant evidence that shows that these two size groups are different (p = 0.127). This means that the robustness requirement in the GPRPE method may have constrained the bloating process such that even doubling the number of fitness evaluations does not lead to larger filter sizes. For highpass filter problems, we reached similar conclusions. Of course, a larger number of replicates may later show these differences to be statistically significant.

4.4.4. Comparison with Other Approaches for Robust Filter Design

There are several previous works focusing on evolving (robust) passive analog filters. Ando and Iba [50] suggested a simple linear genome method to evolve lowpass and bandpass filters with a small number (<50) of components. Particularly related to this work, they used 5%, 10%, and 20% disturbances to the component values to evaluate the robustness of the evolved circuits. Although they did not evolve robust circuits in their work, they suggested that a consideration of manufacturing errors (tolerances) of component values is needed and that in evolvable hardware, the component errors intrinsic in physical components can be absorbed during the evolutionary process. Dib and El-Asir [53] proposed using a symbiotic organisms search (SOS) algorithm to determine the values of the passive components (resistors and capacitors) used in active filters. They compared SOS performance with those obtained using other optimization methods, such as particle swarm optimization, seeker optimization algorithms, and differential evolution, and found that SOS works best. However, they used a fixed topology of the active filters and did not consider the robustness of the evolved filters. In another, similar work [54], analytical characterization of a 2nd order Gm-C Sallen-Key LP filter was considered, as was the corresponding symbolic expressions for their sensitivity to component value variations. When the sensitivity is used as a goal minimization function in the optimization process, parasitic capacitances are canceled and robust designs are obtained. Recently, Lovay [55] proposed using multi-objective optimization, by minimizing the sensitivity to variations in the parameters and the design errors, to determine the values of the passive components (resistors and capacitors) of robust active filters. In summary, these related works all focused on parameter search instead of open-ended topological exploration to evolve robust analog filters.

5. Conclusions

We proposed a GPRD, a novel algorithm for designing robust dynamic systems based on genetic programming and bond graph modeling, and applied it to the synthesis of analog filters. Comparative studies of evolving robustness by GAs (parameter tuning) and by GP (simultaneous topology and parameter search) showed that traditional parameter-tuning-based robust design puts severe limits on the possible robustness achievable through parameter-tuning. We also found that our GPRD robust design algorithm enabled us to find analog filters with much higher robustness (Type-I, Type-II, or both) with respect to the variations in its parameters compared to the parameter-tuning-type evolutionary algorithms (G3PCX-GA) for robust design of fixed functional structures.
There are several issues related to GP-based robust design. One is the computational complexity, since GP is notorious for its high computational demand. Unfortunately, the multiple samplings used to evaluate the robustness of candidate solutions add additional computational load to the evolutionary robust design process. Improved Monte Carlo sampling strategies are needed. The second issue is the need to develop algorithms that evolve robust systems with respect to topology perturbation or component failures. Another related question is how the perturbation pattern affects structures. For example, how does the evolution system respond differentially to component removal as an absolute number or as a percentage of components. Trade-offs between improving robustness by parameter search and by topological innovation will be studied and methods to control it will also be developed. We also observed that, while the proposed methodology is effective for the case study on analog filter design, it may require more refinement; for example, topology constraints may need to be added to the evolved topology to achieve a practical, realistic design.

Acknowledgments

This research was partially supported by the National Natural Science Foundation of China under Grant Nos. 51741101 and 51475097 and Science and Technology Foundation of Guizhou Province under Grant Nos. JZ[2014]2004, Talents[2015]4011 and [2016]5013.

Author Contributions

Jianjun Hu and Shaobo Li conceived and designed the experiments; Jianjun Hu developed the algorithm. Jianjun Hu and Wang Zou performed the experiments; Jianjun Hu and Wang Zou analyzed the data; Jianjun Hu, Shaobo Li, and Wang Zou wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Koza, J.R.; Keane, M.A.; Streeter, M.J.; Mydlowec, W.; Yu, J.; Lanza, G. Genetic Programming IV: Routine Human-Competitive Machine Intelligence; Kluwer Academic Publishers: Alphen aan den Rijn, The Netherlands, 2003. [Google Scholar]
  2. Seo, K.; Fan, Z.; Hu, J.; Goodman, E.D.; Rosenberg, R.C. Toward an Automated Design Method for Multi-Domain Dynamic Systems Using Bond. Mechatronics 2003, 13, 851–885. [Google Scholar] [CrossRef]
  3. Carlson, J.M.; Doyle, J. Complexity and robustness. Proc. Natl. Acad. Sci. USA 2002, 99, 2538–2545. [Google Scholar] [CrossRef] [PubMed]
  4. Du, X.; Chen, W. Towards a Better Understanding of Modeling Feasibility Robustness in Engineering Design. ASME J. Mech. Des. 2000, 122, 385–394. [Google Scholar] [CrossRef]
  5. Jin, Y.; Branke, J. Evolutionary optimization in uncertain environments—A survey. IEEE Trans. Evolut. Comput. 2005, 9, 303–317. [Google Scholar] [CrossRef]
  6. Jin, Y. A Framework for Evolutionary Optimization with Approximate Fitness Functions. IEEE Trans. Evolut. Comput. 2002, 6, 481–494. [Google Scholar]
  7. Jakobi, N.; Husbands, P.; Harvey, I. Noise and the Reality Gap: The Use of Simulation in Evolutionary Robotics. In Proceedings of the Third European Conference on Artificial Life (ECAL95), Granada, Spain, 4–6 June 1995; Springer: Berlin/Heidelberg, Germany, 1995; pp. 704–720. [Google Scholar]
  8. Chen, W.; Allen, J.K.; Mistree, F.; Tsui, K.-L. A Procedure for Robust Design: Minimizing Variations Caused by Noise Factors and Control Factors. ASME J. Mech. Des. 1996, 118, 478–485. [Google Scholar] [CrossRef]
  9. Cancho, R.F.; Janssen, C.; Sole, R.V. Topology of technology graphs: Small world patterns in electronic circuits. Phys. Rev. E 2001, 64, 046119. [Google Scholar] [CrossRef] [PubMed]
  10. Balling, R.J.S.S. Optimization of Coupled Systems: A Critical Review of Approaches. AIAA J. 1996, 34, 6–17. [Google Scholar] [CrossRef]
  11. Suri, R.; Otto, K. Manufacturing System Robustness Through Integrated Modeling. J. Mech. Des. 2001, 123, 630–636. [Google Scholar] [CrossRef]
  12. Ben-Tal, A.; Nemirovski, A. Robust Truss Topology Design via Semidefinite Programming. SIAM J. Optim. 1997, 7, 991–1016. [Google Scholar] [CrossRef]
  13. Koza, J.R.; Bennett, F.H., III; Andre, D.; Keane, M.A.; Dunlap, F. Automated Synthesis of Analog Electrical Circuits by Means of Genetic Programming. IEEE Trans. Evolut. Comput. 1997, 1, 109–128. [Google Scholar] [CrossRef]
  14. Zebulum, R.S.; Pacheco, M.A.; Vellasco, M. Comparison of Different Evolutionary Methodologies Applied to Electronic Filter Design. In Proceedings of the 1998 IEEE World Congress on Computational Intelligence, Anchorage, AK, USA, 4–9 May 1998; pp. 434–439. [Google Scholar]
  15. Lohn, J.; Colombano, S. A circuit representation technique for automated circuit design. IEEE Trans. Evolut. Comput. 1999, 3, 205–219. [Google Scholar] [CrossRef]
  16. Taguchi, G. Taguchi on Robust Technology Development: Bringing; American Society Of Mechanical Engineers: New York, NY, USA, 1993. [Google Scholar]
  17. Zhu, J. Performance Distribution Analysis and Robust Design. J. Mech. Des. 2001, 123, 11–17. [Google Scholar] [CrossRef]
  18. Weng, T.-W. Stochastic simulation and robust design optimization of integrated photonic filters. Nanophotonics 2017, 6, 299–308. [Google Scholar] [CrossRef]
  19. Sun, J.; Xiao, L.; Tian, J.; Zhou, H.; Roveda, J. Surrogating circuit design solutions with robustness metrics. Integr. VLSI J. 2016, 52, 1–9. [Google Scholar] [CrossRef]
  20. Husbands, P.; Harvey, I. Evolution versus design. Controlling autonomous robots. In Integrating Perception, Planning and Action, Proceedings of the Third Annual Conferences on Artificial Intelligence, Perth, Australia, 8–10 July 1992; IEEE Press: Los Alamitos, CA, USA, 1992. [Google Scholar]
  21. Husband, P.; Harvey, I.; Cliff, D. An evolutionary approach to situated AI. In Prospects for Artificial Intelligence, Proceedings of the Ninth bi-annual Conference of the Society for the Study of Artificial Intelligence and the Simulation of Behavior, Birmingham, UK, 29 March–2 April 1993; IO Press: Amsterdam, The Netherlands, 1993; pp. 61–70. [Google Scholar]
  22. Jakobi, N. Evolving Motion-Tracking Behavior for a Panning Camera Head. In From Animals to Animats 5, Proceedings of the 5th International Conference on Simulation of Adaptive Behavior; MIT Press: Cambridge, MA, USA, 1997. [Google Scholar]
  23. Jakobi, N. The Minimal Simulation Approach to Evolutionary Robotics. In Evolutionary Robotics—From Intelligent Robots to Artificial Life (ER’98); Gomi, T., Ed.; AAI Books: Ottawa, ON, Canada, 1998. [Google Scholar]
  24. Tsutsui, S.; Ghosh, A. Genetic algorithms with a robust solution searching scheme. IEEE Trans. Evolut. Comput. 1997, 1, 201–208. [Google Scholar] [CrossRef]
  25. Wiesmann, D.; Hammel, U.; Bäck, T. Robust design of multilayer optical coatings by means of evolutionary algorithms. IEEE Trans. Evolut. Comput. 1998, 2, 162–167. [Google Scholar] [CrossRef]
  26. Forouraghi, B. A genetic algorithm for multiobjective robust design. Appl. Intell. 2000, 12, 151–161. [Google Scholar] [CrossRef]
  27. Ray, T. Constrained robust optimal design using a multi-objective evolutionary algorithm. In Proceedings of the Congress on Evolutionary Computation, Honolulu, HI, USA, 12–17 May 2002; IEEE Press: Los Alamitos, CA, USA, 2002; pp. 419–424. [Google Scholar]
  28. Jin, Y.; Sendhoff, B. Trade-off between optimality and robustness: An evolutionary multi-objective approach. In Proceedings of the Second International Conference on Evolutionary Multi-criterion Optimization, Faro, Portugal, 8–11 April 2003; Springer: Berlin/Heidelberg, Germany, 2003; pp. 237–251. [Google Scholar]
  29. Branke, J. Creating robust solutions by means of an evolutionary algorithms. In Proceedings of the Parallel Problem Solving from Nature Number 1498 in LNCS; Eiben, A., Bäck, T., Schoenauer, M., Schwefel, H.P., Eds.; Springer: Berlin/Heidelberg, Germany, 1998; pp. 119–128. [Google Scholar]
  30. Das, I. Robustness Optimization for constrained nonlinear programming problems. Eng. Optim. 2000, 32, 585–618. [Google Scholar] [CrossRef]
  31. Hammel, U.; Bäck, T. Evolution Strategies on Noisy Functions. How to Improve Convergence Properties. In Parallel Problem Solving from Nature; Davidor, Y., Manner, R., Schwefel, H.P., Eds.; Springer: Heidelberg, Germany, 1994; pp. 159–168. [Google Scholar]
  32. Branke, J. Evolutionary Optimization in Dynamic Environments; Kluwer: Alphen aan den Rijn, The Netherlands, 2001. [Google Scholar]
  33. Branke, J.; Schmeck, H. Designing evolutionary algorithms for dynamic optimization problems. In Theory and Application of Evolutionary Computation: Recent Trends; Tsutsui, S., Ghosh, A., Eds.; Springer: Heidelberg, Germany, 2002. [Google Scholar]
  34. Lee, W.P.; Hallam, J.; Lund, H.H. Applying Genetic Programming to Evolve Behavior Primitives and Arbitrators for Mobile Robots. In Proceedings of the IEEE 4th International Conference on Evolutionary Computation, Indianapolis, IN, USA, 13–16 April 1997; IEEE Press: Los Alamitos, CA, USA, 1997; Volume 1. [Google Scholar]
  35. Thompson, A. On the automatic design of robust electronics through artificial evolution. In Proceedings of the International Conference on Evolvable Systems, Lausanne, Switzerland, 23–25 September 1998; Springer: Berlin/Heidelberg, Germany, 1998; pp. 13–24. [Google Scholar]
  36. Thompson, A.; Layzell, P. Evolution of robustness in an electronics design. In Proceedings of the Evolvable Systems, from Biology to Hardware (ICES 2000), Edinburgh, UK, 17–19 April 2000; Springer: Berlin/Heidelberg, Germany, 2000; pp. 218–228. [Google Scholar]
  37. Tyrrell, A.; Hollingworth, G.; Smith, S. Evolutionary strategies and intrinsic fault tolerance. In Proceedings of the Third NASA/DoD Workshop on Evolvable Hardware (EH 2001), Long Beach, CA, USA, 12–14 July 2001; Keymeulen, D., Stoica, A., Lohn, J., Zebulum, R.S., Eds.; IEEE Computer Society: Long Beach, CA, USA, 2001; pp. 98–106. [Google Scholar]
  38. Miller, J.F.; Hartmann, M. Evolving messy gates for fault tolerance: Some preliminary findings. In Proceedings of the Third NASA/DoD Workshop on Evolvable Hardware (EH 2001), Long Beach, CA, USA, 12–14 July 2001; Keymeulen, D., Stoica, A., Lohn, J., Zebulum, R.S., Eds.; IEEE Computer Society: Long Beach, CA, USA, 2001; pp. 116–123. [Google Scholar]
  39. Hartmann, M.; Haddow, P. Evolution of fault-tolerant and noise-robust digital designs. In IEE Proceedings of Computers and Digital Techniques; IEEE Press: Los Alamitos, CA, USA, 2004; pp. 287–294. [Google Scholar]
  40. Hartmann, M.; Eskelund, F.; Haddow, P.C.; Miller, J.F. Evolving Fault Tolerance On An Unreliable Technology Platform. In Proceedings of the Genetic and Evolutionary Computation Conference (GECCO 2002), New York, NY, USA, 9–13 July 2002; Langdon, W.B., Cantú-Paz, E., Mathias, K., Roy, R., Davis, D., Poli, R., Balakrishnan, K., Honavar, V., Rudolph, G., Wegener, J., et al., Eds.; Morgan Kaufmann Publishers: New York, NY, USA, 2002; pp. 171–177. [Google Scholar]
  41. Koza, J.R.; Andre, D.; Bennett, F.H., III; Keane, M. Genetic Programming 3: Darwinian Invention and Problem Solving; Morgan Kaufman: Burlington, MA, USA, 1999. [Google Scholar]
  42. Goh, C.; Li, Y. GA Automated Design and Synthesis of Analog Circuits with Practical Constraints. In Proceedings of the 2001 Congress on Evolutionary Computation 2001, Seoul, Korea, 27–30 May 2001; IEEE Press: Los Alamitos, CA, USA, 2001; pp. 170–177. [Google Scholar]
  43. Fan, Z.; Hu, J.; Seo, K.; Goodman, E.D.; Rosenberg, R.C.; Zhang, B. Bond Graph Representation and GP for Automated Analog Filter Design. In Proceedings of the 2001 Genetic and Evolutionary Computation Conference Late Breaking Papers, San Francisco, CA, USA, 7–11 July 2001; pp. 81–86. [Google Scholar]
  44. Karnopp, D.; Margolis, D.L.; Rosenberg, R.C. System Dynamics: Modeling and Simulation of Mechatronic Systems, 3rd ed.; John Wiley & Sons, Inc.: New York, NY, USA, 2000. [Google Scholar]
  45. Hu, J.; Goodman, E.; Rosenberg, R. Topological search in automated mechatronic system synthesis using bond graphs and genetic programming. In Proceedings of the American Control Conference (ACC 2004), Boston, MA, USA, 30 June–2 July 2004. [Google Scholar]
  46. Hu, J.; Goodman, E.; Seo, K.; Fan, Z.; Rosenberg, R. The Hierarchical Fair Competition (HFC) Framework for Sustainable Evolutionary Algorithms. Evolut. Comput. 2005, 13, 241–247. [Google Scholar] [CrossRef] [PubMed]
  47. Tay, E.; Flowers, W.; Barrus, J. Automated Generation and Analysis of. Res. Eng. Des. 1998, 10, 15–29. [Google Scholar] [CrossRef]
  48. Seo, K.; Fan, Z.; Hu, J.; Goodman, E.D.; Rosenberg, R.C. Dense and Switched Modular Primitives for Bond Graph Model Design. In Proceedings of the Genetic and Evolutionary Computation (GECCO-2003), Chicago, IL, USA, 12–16 July 2003; Cantú-Paz, E., Foster, J.A., Deb, K., Davis, D., Roy, R., O’Reilly, U.M., Beyer, H.G., Standish, R., Kendall, G., Wilson, S., Eds.; Springer: Chicago, IL, USA, 2003; Volume 2724, pp. 1764–1775. [Google Scholar]
  49. Hu, J.; Li, S.; Goodman, E. Evolutionary Robust Design of Analog Filters Using Genetic Programming. In Evolutionary Computation in Dynamic and Uncertain Environments; Springer: Berlin/Heidelberg, Germany, 2005; pp. 479–496. [Google Scholar]
  50. Ando, S.; Iba, H. Linear Genome Methodology for Analog Circuit Design; Technical Report; University of Tokyo: Tokyo, Japan, 2000. [Google Scholar]
  51. Fang, Y.; Li, J. A Review of Tournament Selection in Genetic Programming. In Advances in Computation and Intelligence; Cai, Z., Hu, C., Kang, Z., Liu, Y., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 181–192. [Google Scholar]
  52. Deb, K.; Anand, A.; Joshi, D. A Computationally Efficient Evolutionary Algorithm for Real-Parameter Optimization. Evolut. Comput. 2002, 10, 345–369. [Google Scholar] [CrossRef] [PubMed]
  53. Dib, N.; El-Asir, B. Optimal design of analog active filters using symbiotic organisms search. Int. J. Numer. Model. Electron. Netw. Devices Fields 2018. [Google Scholar] [CrossRef]
  54. Leitão, P.M.V.; Fino, H. Robust Optimization-Based High Frequency Gm-C Filter Design. In Technological Innovation for Value Creation; Camarinha-Matos, L.M., Shahamatnia, E., Nunes, G., Eds.; Springer: Berlin/Heidelberg, Germany, 2012; pp. 465–474. [Google Scholar]
  55. Lovay, M.; Peretti, G.; Romero, E. Application of genetic algorithms in the design of robust active filters. In Proceedings of the 2015 Argentine School of Micro-Nanoelectronics, Technology and Applications (EAMTA), Villa Maria, Argentina, 30–31 July 2015; pp. 1–6. [Google Scholar]
Figure 1. A bond graph and its equivalent electrical circuit. The dotted boxes in the left bond graph indicate modifiable sites at which further topological manipulations can be applied (to be explained in the next section).
Figure 1. A bond graph and its equivalent electrical circuit. The dotted boxes in the left bond graph indicate modifiable sites at which further topological manipulations can be applied (to be explained in the next section).
Algorithms 11 00026 g001
Figure 2. The bond graph structure of a lowpass filter evolved by the GPBG method with 500,000 function evaluations. This filter has 39 components without counting the embryo components. (Component sizing values are omitted in the figure for simplicity).
Figure 2. The bond graph structure of a lowpass filter evolved by the GPBG method with 500,000 function evaluations. This filter has 39 components without counting the embryo components. (Component sizing values are omitted in the figure for simplicity).
Algorithms 11 00026 g002
Figure 3. An example of a GP tree, composed of topology operators applied to an embryo, generating a bond graph after depth-first execution (numeric ERC nodes are omitted). Note that the 010 and 001 are the flag bits showing the presence or absence of attached C/I/R components.
Figure 3. An example of a GP tree, composed of topology operators applied to an embryo, generating a bond graph after depth-first execution (numeric ERC nodes are omitted). Note that the 010 and 001 are the flag bits showing the presence or absence of attached C/I/R components.
Algorithms 11 00026 g003
Figure 4. The topology of the best highpass analog filters evolved with standard GP with 500,000 function evaluations without considering a robustness requirement in the fitness function(parameters are omitted). This filter has 27 C/I/R components without counting the original embryo components. The best evolved lowpass filter is shown in Figure 2. These topologies are the results of a simplification procedure that removes redundancy in the original evolved bond graphs while their functional behaviors are maintained.
Figure 4. The topology of the best highpass analog filters evolved with standard GP with 500,000 function evaluations without considering a robustness requirement in the fitness function(parameters are omitted). This filter has 27 C/I/R components without counting the original embryo components. The best evolved lowpass filter is shown in Figure 2. These topologies are the results of a simplification procedure that removes redundancy in the original evolved bond graphs while their functional behaviors are maintained.
Algorithms 11 00026 g004
Figure 5. Frequency responses of best lowpass and highpass filters evolved using GP, GPGARMS, GPRMS, and GPRPE with 20% Gaussian perturbation of their components: (a) Frequency response distribution of the best four lowpass filters. The total sums of 101-point deviations from the target response curve are 6.43 (GP), 33.03 (GPGARMS), 9.61 (GPRMS), and 4.13 (GPRPE). (b) Frequency response distribution of the best four highpass filters. The total deviations from the target response curve are 0.32 (GP), 42.69 (GPGARMS), 2.53 (GPRMS), and 0.19 (GPRPE).
Figure 5. Frequency responses of best lowpass and highpass filters evolved using GP, GPGARMS, GPRMS, and GPRPE with 20% Gaussian perturbation of their components: (a) Frequency response distribution of the best four lowpass filters. The total sums of 101-point deviations from the target response curve are 6.43 (GP), 33.03 (GPGARMS), 9.61 (GPRMS), and 4.13 (GPRPE). (b) Frequency response distribution of the best four highpass filters. The total deviations from the target response curve are 0.32 (GP), 42.69 (GPGARMS), 2.53 (GPRMS), and 0.19 (GPRPE).
Algorithms 11 00026 g005
Figure 6. Frequency responses of the best lowpass and highpass filters evolved using GP, GPGARMS, GPRMS, and GPRPE in the perturbed-parametric environment: (a) Frequency response distribution of the best filter evolved using standard GP. (b) Frequency response distribution of the best filter evolved using standard GPGARMS. (c) Frequency response distribution of the best filter evolved using standard GPRMS. (d) Frequency response distribution of the best filter evolved using standard GPRPE.
Figure 6. Frequency responses of the best lowpass and highpass filters evolved using GP, GPGARMS, GPRMS, and GPRPE in the perturbed-parametric environment: (a) Frequency response distribution of the best filter evolved using standard GP. (b) Frequency response distribution of the best filter evolved using standard GPGARMS. (c) Frequency response distribution of the best filter evolved using standard GPRMS. (d) Frequency response distribution of the best filter evolved using standard GPRPE.
Algorithms 11 00026 g006
Figure 7. Robustness of highpass filters evolved by GP, GPGARMS, GPRMS, and GPRPE. For (a,b), lower values correspond to higher robustness. (a) The average Type-I robustness defined in Section 3.5. (b) The average Type-II robustness defined in Equation (5).
Figure 7. Robustness of highpass filters evolved by GP, GPGARMS, GPRMS, and GPRPE. For (a,b), lower values correspond to higher robustness. (a) The average Type-I robustness defined in Section 3.5. (b) The average Type-II robustness defined in Equation (5).
Algorithms 11 00026 g007
Figure 8. Evolved robust lowpass and highpass filters with complexity much lower than the filters evolved without considering robustness in the fitness in Figure 2 and Figure 4. (a) Topology of the most robust lowpass filter with only 13 evolved components using GPRMS. (b) Topology of the most robust lowpass filter with only 16 evolved components using GPRMS. (c) Equivalent analog circuit of the most robust lowpass filter in (b).
Figure 8. Evolved robust lowpass and highpass filters with complexity much lower than the filters evolved without considering robustness in the fitness in Figure 2 and Figure 4. (a) Topology of the most robust lowpass filter with only 13 evolved components using GPRMS. (b) Topology of the most robust lowpass filter with only 16 evolved components using GPRMS. (c) Equivalent analog circuit of the most robust lowpass filter in (b).
Algorithms 11 00026 g008
Table 1. Experimental parameters for analog filter synthesis without robustness criterion.
Table 1. Experimental parameters for analog filter synthesis without robustness criterion.
Total Population Size: 2000 (400/400/400/400/400)Number of Subpopulations: 5
Migration interval: 5 generationsMigration size: 30 individuals
Max tree depth: 8Crossover probability: 0.7
InitTreeDepth: 3–5Standard mutation probability: 0.1
Flag bit mutation rate: 0.1Swapping-tree mutation rate: 0.1
Tournament size: 7Parametric mutation probability: 0.5
Max evaluations: 1,000,000Flag mutation probability: 0.3
Pool size of elite individuals: 20Elite pool update frequency: 5 generations
Table 2. Experimental parameters for robust design by G3PCX-GA.
Table 2. Experimental parameters for robust design by G3PCX-GA.
Total Population Size: 200Max Evaluations: 500,000
Number of parents in crossover: 3Family size: 2
σ ζ : 0.1 σ η : 0.1
SPI: 10Perturbation noise percentage: 20%
Table 3. The performance and robustness of evolved lowpass filters using GP, GPGARMS, GPRMS, and GPRPE.
Table 3. The performance and robustness of evolved lowpass filters using GP, GPGARMS, GPRMS, and GPRPE.
Perturbation Magnitude0.10.20.30.40.5
MeanStderrMeanStderrMeanStderrMeanStderrMeanStderr
f r o b u s t I GP36.575.5540.045.0045.334.5654.684.2561.914.06
GPGARMS65.166.6050.086.0564.333.6673.853.8479.815.28
GPRMS48.833.9951.363.6756.093.4663.913.5272.713.72
GPRPE14.082.6519.642.2227.251.8237.511.8148.982.21
f r o b u s t I I GP2.440.405.590.7914.371.6130.202.8342.783.76
GPGARMS5.421.8510.811.7218.182.1233.571.8554.054.45
GPRMS1.350.344.580.5413.581.0528.461.9540.882.64
GPRPE3.920.468.130.7820.853.7530.892.3643.112.92
Table 4. The sizes of filters evolved using four algorithms with different numbers of fitness evaluations. Each algorithm was run 10 times. * A total of 1,000,000 function evaluations are required to do 100,000 actual fitness evaluations.
Table 4. The sizes of filters evolved using four algorithms with different numbers of fitness evaluations. Each algorithm was run 10 times. * A total of 1,000,000 function evaluations are required to do 100,000 actual fitness evaluations.
Algorithm (Fitness Evaluations)12345678910MeanStdev
GP (100,000)34263524142024252646278.9100
GP (500,000)4647374757404045487951.750015.3200
GPRMS (100,000) *2815362128144326312026.20009.1100
GPRPE (1,000,000)6736463637572963273643.400014

Share and Cite

MDPI and ACS Style

Li, S.; Zou, W.; Hu, J. A Novel Evolutionary Algorithm for Designing Robust Analog Filters. Algorithms 2018, 11, 26. https://doi.org/10.3390/a11030026

AMA Style

Li S, Zou W, Hu J. A Novel Evolutionary Algorithm for Designing Robust Analog Filters. Algorithms. 2018; 11(3):26. https://doi.org/10.3390/a11030026

Chicago/Turabian Style

Li, Shaobo, Wang Zou, and Jianjun Hu. 2018. "A Novel Evolutionary Algorithm for Designing Robust Analog Filters" Algorithms 11, no. 3: 26. https://doi.org/10.3390/a11030026

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop