due to the effects of rapid globalization in recent years , large airports around the world , including narita airport in japan , are forced to increase the number of takeoffs and landings they handle . therefore , taxiing aircraft in the airport , which causes congestion , has been an issue of concern . the taxiing of aircraft is controlled by ground control ( gc ) , which manages taxiing near the terminal , and air traffic controllers , who look down upon the entire airport from the control tower . to start , as shown in figure 1(a ) , departing aircraft  pushback  on instruction from gc from the point called  spot  to passenger embarkation and debarkation as shown from ( 1 ) and ( 2 ) . pushback is an operation that pulls an aircraft that structurally can not taxi backward . in response to instructions for the route to the runway from the air traffic controllers after the pushback , from gate , which as the entrance of the ramp area , aircraft move from ( 3 ) out of the ramp area near the terminal , travel along it to the runway as shown in ( 4 ) , and take off . arriving aircraft , as shown in figure 1(b ) , land on the runway in ( 1 ) , receive instructions for taxiing route from air traffic controllers , and travel to the gate shown in ( 2 ) . then , gc issues a permit of entry for the ramp area at the right time for any state aircraft or vehicle that may be driving through the ramp area to the spot of interest shown in ( 3 ) . gc must control the aircraft ( of which there are more than 30 at the same time in an airport at peak times ) quickly and safely to allow these aircraft to taxi . for example , as the highlighted areas of figure 1(a ) have a bidirectional travel route , there is a need for appropriate control to ensure that deadlock does not occur . for ease of traffic congestion , it is necessary to direct all in concert to control the start timing of the pushback , and the entry timing of the ramp area . the task is performed heuristically , but there are currently concerns as to if the heuristics can cope with a future increase in flights . therefore , it is required to computationally generate control rules of gc , for variety in schedules of aircraft , which we call robustness of the solution in this paper . the aircraft control problem can be regarded as an optimization problem of constructing decision - making rules for a multiagent system [ 1 , 2 ] . this term refers to a system that hopes to achieve complex goals through the operation of multiple agents repeating actions based on a set of rules and observation of the environment . there are a variety of examples present in the literature other than aircraft control . in [ 3 , 4 ] , they introduced the concept of using agent - based systems for the construction of a flexible production system with high efficiency . they define two kinds of agents : the first one is to make a work plan for each device and the second one is to integrate the plans . in the paper , they design behavioral rules for the agents and measure the efficiency of the system as a whole . in , multiagent system is proposed for traffic flow maximization via traffic signal control . with each traffic signal independently determining the congestion state of its surroundings , it can continuously change the time interval of the signal accordingly and can be successfully efficient at maximizing traffic flow as a result . in , they model an agent - based system for workflow involving business loans from a bank , such as requests received from customers . this system independently models the performance of actions such as requests received by each agent . decision - making rules in any agent - based systems such as these , where decisions are not reversible in all situations , must not just be robust but in light of maintenance must also be very readable . the proposed system here for controlling aircraft is intended to be implemented by humans using rules . as a result , it is necessary to calculate highly readable rules that can be operated and understood by humans . for the above reasons , it is important to use a method of generating behavioral rules for multiagent system that takes into account both robustness and readability . from the viewpoint of the robust optimization problem , in , a robust optimization technique is proposed for known linear programming problems with variation . deterministic algorithms guarantee the robustness of calculating solutions , but a modeled problem must be linear to be solved . for many systems , application to nonlinear problems is difficult to realize , and robust optimization techniques using simulation are required . in , a robust optimization technique with monte carlo method is presented . in this method , simulations are repeatedly performed by using random numbers , and the robustness of the solution is calculated from the results . however , if the simulation is made on a large scale , computation time becomes enormous , and the accuracy of any solution calculated within a realistic time is very low . in [ 9 , 10 ] they propose methods of reducing the number of simulations by combining it with a curved surface approximation of the evaluation function through response surface methodology . however , these methods have been applied to problems with continuous solution space , and applying them to discrete multimodal problems such as rule generation is difficult . further , as described above , there is still no method that contains measures to reduce the number of simulations without surface approximation , which would be necessary for rule generation for multiagent systems . in , an approximate optimization method is utilized on problems with a tree - structure solution that mimics the evolution of organisms called genetic programming ( gp ) . this optimization technique has been studied many times , its rules taking the form of polynomial conditions . however , this method suffers from a phenomenon called bloat , where solutions become increasingly complicated ; it calculates solutions with low readability . in order to suppress bloat , in , a method is proposed for setting a new evaluation function that combines the depth of the tree structure with what the solution evaluation function should be . further , in , a method is presented of solving a multiobjective optimization problem treating the depth of the tree structure as a new evaluation function . however , these methods are discussed only in terms of improving the robustness of the rule , with no consideration towards the readability of the solution . in , rules are generated by optimizing tempo and pitch to match user preferences , which is represented as a polynomial using arithmetic operations , and so forth , over several parameters , but it has very low readability . from the point of view of software maintenance , in , they enumerate about twenty different elements that affect the readability of source code of a computer program , and in , an evaluation expression for measuring element readability is proposed . in addition , in , although a rule generation method for multiagent systems has been proposed , evaluation was not sufficient from the viewpoint of the scale of the problems to evaluate . as shown in the previous section , research on rule generation has been focused on optimization from the viewpoint of a kind of efficiency or robustness , with little focus on readability . therefore , the purpose of this study is the construction of a simulation - based rule generation method that considers both robustness and readability . in this paper , we propose a method that incorporates a method both with pado   and with operator equalization   using metrics on readability . then , we apply our method to the aircraft control problem in two settings and evaluate the results . there is more than one agent in the environment , with each agent performing observations to produce the state space based on the information obtained . an action is decided using predetermined rules based on the state space , and this is repeated several times . we assume that the state is defined in advance in accordance with the problem , that any error in the environment that may occur is known , but that there is no error when retrieving information . what actions the agents may select is defined in advance according to the type of problem and on the basis of rules . we evaluate the rules in the simulation environment multiple times with different settings . in this paper , as an evaluation value , we utilize total stopping time of all the agents ( aircraft ) , which does not take negative value and whose optimal number is equal to zero ( so - called s - type characteristics ) . then we adopt the index for measuring the robustness of the rules , as shown in(1)robustness = log1e2,e :  evaluation value in each simulation.in addition , our rules are evaluated using the readability formula described in section 2.2 . in this study , we evaluate the rules using the readability formula proposed by   ( 2)readability=11+ez , z=8.870.033v+0.40  lines1.5  entropy , v :  halstead 's volume  = programlengthlogprogramvocabulary , programlength :  number  of  words  in  the  documents , programvocabulary :  number  of  different  words  in  the  documents , lines :  number  of  lines , entropy :  hx=i=1npxilog2pxi , pxi = count(xi)j=1ncount(xj),countx :  number  of  appearances  of  the  word   x,n :  number  of  words.this formula measures the readability of source code through a variety of indicators . in formula ( 2 ) , v is an index called halstead 's volume   that represents the complexity of the document as a function of the length of the document , increasing with complexity . line represents the number of lines in the document , which we have quantified as the number of rows . entropy represents the amount of information in the document , which takes an increasing value as the probability of occurrence of each word falls . according to the readability index mentioned in the previous section , line breaks in the document help to improve readability . therefore , we propose a hierarchical two - step representation of rules , condition and action . here we then represent the rules by a tree - structure inequation or equation , to graphically represent the order of evaluation . the graph determines whether the agent starts to take a specific action or not , such as to start moving . the graph has three types of nodes : a start node , condition nodes , and action nodes , as shown in figure 2 . condition nodes are associated with equality / inequality conditions , which are represented by trees as later defined , and action nodes are associated with whether the agent to take the action ( y ) or not ( n ) . a condition node has an outdegree of two , represented by the arrows . on receipt of a state vector , the rule outputs whether to take a single action or not . the state vector is given as input , starting from the start node and arriving at the condition node . the process then proceeds to the next node , following the red arc if false or the blue arc if true . the process outputs the behavior associated with the action node when it arrives at the node , and that action is then taken by the previously defined agent . the equality / inequality conditions associated with the condition nodes of the graph are represented by trees consisting of three kinds of nodes , root nodes , nonterminal nodes , and terminal nodes , as shown in figure 3 . the root node represents a comparison operator ( < , = ,  , and > ) , the right - hand side represents a constant ( in this tree , 0 ) , and the left - hand side represents arithmetic operators , a decomposition of a polynomial in the state variables and constants . terminal nodes on the left - hand side represent the state variables or constants , while nonterminal nodes represent arithmetic operators . when evaluating a formula represented in such a tree , the state vector is input , and the result calculated by substituting the values of the state variables in the expression and using the comparison operator of the root node against 0 on the right side is output . in the case of figure 3 , it expresses an inequality equation of t(a  8)(w  8) < 0 , where a , t , and w are state variables . table 1 describes the variables appearing in formula ( 2 ) and their correspondence to both the source code and the representation of rules . here we understand that the rules are similar to the source codes in that humans need to read and understand them . in this representation , since the total combination of tree and graph forms a ruleset corresponding to the entire document in the source code , we associate the nodes in the graph with lines of source code and treat the number of nodes across all conditional trees as the number of words . an overall picture of the algorithm is shown in figure 4 . by the problem setting described in the previous section , the simulation model to be calculated is our input , and it is assumed that the simulation model at this time meets the requirements as described in the previous section . the main components of the algorithm are the following four sections.simulator;candidate rule generator;constraints;readability calculator.in the simulator , using the simulation model as an input , our method generates an evaluation solution using candidate rules also serving as input and outputs the evaluation value using an index that is defined in advance . as we repeatedly create new solution candidates , we approach the optimal solution . without attachment to generate a solution structure that is layered as described in the previous section , the solution candidate generator will be in a hierarchical structure . further , in the constraint section , every time the simulation is performed , we automatically extract and save the constraints of the rules for a feasible solution . then , in generating a new solution candidate by regenerating and referencing previous solution candidates , they must also satisfy these constraints . when calculation first begins , the graph portion of the solution candidates is generated by the graph generating unit . then , it generates a set of trees ( equal to the number of conditional nodes ) in the amount required for the graph tree generator . for each candidate solution , we evaluate the simulator , with random variables varying according to the distribution of each . the simulation is performed multiple times , and the robustness of the candidate solution calculated according to formula ( 1 ) . we believe that an infeasible solution that can not be calculated becomes feasible with repeated minor modifications of the tree through rule generation , and the evaluation value falls largely within the solution space . after the evaluation value of some graphs is calculated in this manner , we use the gathered information to produce an improved graph using a technique known as pado , first proposed by , and through a method called operator equalization . when improving robustness and readability is treated as a multiobjective optimization problem , such as in , the problem , as stated in , will fall into local optimal solutions , which undesirably lowers the quality of the solution . in this study , we prevent this by building in optimization of the graph by pado via a technique called operator equalization . this technique , first proposed in , has been suggested as a method of suppressing the excess complication of solutions through gp via the most fundamental tree optimization methods . with this method we take into account both the robustness and the complexity of each individual in the previous generation in generating the children that make up the next generation of rules . like gp , pado as used in this method is also a method that approximates the solution by mimicking natural selection ; in this case , it is akin to repeated cross survival selection and mutation by performing optimization of the graph by the concept of operator equalization . we show the algorithm in a concrete fashion below . creating a distribution of evaluation values from readability was calculated by applying formula ( 2 ) to the current generation , such as shown in the upper graph in figure 5 . to determine the distribution of the complexity of the next generation individual using this distribution , the next generation is allocated in proportion to the magnitude of evaluation values for the complexity of the previous generation . that is , assuming total population in pado as n , it is determined in the ratio , as shown in the lower table in figure 5 . after the distribution has been determined , individuals are produced by crossing the individual child candidates , to perform the calculation of the evaluation value of a decision tree that will lead to a condition node by the local search method described above and to calculate the readability by the readability index . for example , in case of figure 6 , n is set as 7 and maximal numbers of each category are obtained in advance , as 1 for the category whose readability band is between 0 and 0.2 , 2 for the band between 0.2 and 0.4 , 4 for the band between 0.4 and 0.6 , 1 for the band between 0.6 and 0.8 . the individual can be assigned to the band between 0.2 and 0.4 because it does not exceed the number of individuals allocated to the band . the graphs are first selected not to exceed the number of individuals allocated to the corresponding readability band in figure 5 . if the value of formula ( 1 ) for a certain individual is best among others and it violates the readability band when it is selected , the graph is selected by slightly changing the number of individuals allocated to each band . simulations are performed using the entire narita airport to evaluate the performance of the proposed method , by showing how the proposed one can generate feasible control rules of taxiing aircraft . the evaluation is made by comparison with the present control rule in narita airport and that with other methods . figure 7 represents the entire narita airport , with about 300 spots . as a result of discussion with the airport officials , it was found that there are two specific local areas that most required an improved control method for congestion relief , as shown in figure 7 . in area ( i ) of figure 7 , aircraft begin with a pushback as they attempt to start out . attempts to head to or from the runway are complicated , and nearly deadlock situation of the loop due to pushback operations can force the expensive act of waiting increasingly long in the air while the runway is unusable due to a waiting column of aircraft . therefore , the aim is to alleviate congestion by controlling the proper timing of the pushback in this spot . in area ( ii ) of figure 7 , as described in section 1 , the ramp area above the gate is a bidirectional traveling path . aircraft must wait at this gate point depending on the status of the ramp area , and so forth , for deadlock prevention . however , there is a possibility that waiting on the spot so long will stop the travel of the aircraft , which leads to congestion . therefore , the aim is to alleviate congestion by optimizing the wait - passing rules of the gate . the control rules of these two areas affect each other , and so it is necessary to design them simultaneously with each other . evaluating the proposed method by using a multiagent system model as described above will be with reference to the following : the specific scheme in control of gate control and pushback and the state variables used in each rule . ( i ) control of pushback . departing aircraft at the spots asks for instruction to the air traffic controllers when it becomes their departure time . air traffic controllers issue instructions to either start or wait in order starting with the aircraft that have been waiting longest . the aircraft receives an instruction either to start the pushback or to report back to the controller after a certain time . this time interval is equal to the minimum time interval in the simulation , which is set to one second in the present study . state variables. consider the following.the number of taxiing aircraft in area ( i ) : expressed as   airn.the number of arrivals coming to area ( i ) : expressed as   usesp.the number of aircraft pushing back in area ( i ) : expressed as   pushn.the number of waiting aircraft in area ( i ) : expressed as   waitn.the waiting time of the departing aircraft : expressed as   time . the number of taxiing aircraft in area ( i ) : expressed as   airn . the number of arrivals coming to area ( i ) : expressed as   usesp . the number of aircraft pushing back in area ( i ) : expressed as   pushn . the number of waiting aircraft in area ( i ) : expressed as   waitn . when the aircraft passes through the appropriate gate , it reports to the air traffic controllers . the air traffic controllers evaluate the rules and instruct the aircraft to either pass or wait . an aircraft that receives the instruction to pass enters the ramp area to pass through the gate . an aircraft that receives the instruction to wait must wait to report back to the controller after a certain time . this time interval is equal to the minimum time interval in the simulation , which is set to one second in the present study . state variables. consider the following.whether the aircraft is waiting on each spot : expressed as   si_sn  ( 1 when waiting at spot i , 0 otherwise).whether the aircraft is pushing back on each spot : expressed as   si_pb  ( 1 when with pushback at spot i , 0 otherwise).the waiting time of an aircraft on each spot : expressed as   si_st  ( waiting time at spot i).the waiting time of an aircraft on the stop line : expressed as   sus_st.the number of waiting aircraft by the stop line : expressed as   sus_sn.whether an aircraft is taxiing on the ramp area : expressed as   onramp ( 1 when taxiing on the ramp area , 0 otherwise).the kind of aircraft ( departing or arriving ) taxiing on the ramp area : expressed as   ramptype  ( 1 when the departing aircraft is taxiing , 2 when arriving aircraft is taxiing , 0 otherwise).we assume this simulation is during the busy time in narita , time slot from 13:30 to 19:30 , when the departing aircraft and arriving aircraft are equally three per 30 minutes , which means that totally six aircraft are departing or arriving within 30 minutes in the airport . we run the simulations a fixed number of times ( 50 in this paper ) , with each candidate rule system , setting n ( total population in pado ) = 20 , using the characteristic of total stopping time of all aircraft as the evaluation value for the robustness of the set of rules as a solution candidate as in formula ( 1 ) . departing and arriving times of all aircraft determined randomly assuming uniform distribution in 30 minutes . whether the aircraft is waiting on each spot : expressed as   si_sn  ( 1 when waiting at spot i , 0 otherwise ) . whether the aircraft is pushing back on each spot : expressed as   si_pb  ( 1 when with pushback at spot i , 0 otherwise ) . the waiting time of an aircraft on each spot : expressed as   si_st  ( waiting time at spot i ) . the waiting time of an aircraft on the stop line : expressed as   sus_st . the number of waiting aircraft by the stop line : expressed as   sus_sn . whether an aircraft is taxiing on the ramp area : expressed as   onramp ( 1 when taxiing on the ramp area , 0 otherwise ) . the kind of aircraft ( departing or arriving ) taxiing on the ramp area : expressed as   ramptype  ( 1 when the departing aircraft is taxiing , 2 when arriving aircraft is taxiing , 0 otherwise ) . the result of the above is shown in figure 8 , showing that it was possible to calculate the pareto optimal solution of the two problems described above . the rule at the start of the pushback region in the solution of ( i ) is to wait if it has not been over 16 seconds , and if it has been 16 seconds or more , then proceed to evaluation of condition 01 . the aircraft can start only after waiting even if empty , and then only if no aircraft is taxiing ( represented by the variable a corresponding to state variable 1 ) . in the solution of ( ii ) , since the condition node 00 is always true in any state , the rule for passing through to the gate region in area ( ii ) is to always allow entry . this solution was the most readable among solutions , and when the simulation was performed using these rules , the probability that deadlock did not occur was 99.1% . the rule at the start of the pushback region in the solution of ( i ) involves the variables   airn  ( the amount of aircraft taxiing ) ,   usesp  ( the number of arrivals coming to area ( i ) ) , and   waitn  ( the number of waiting aircraft ) . this condition can be interpreted as a rule that depends solely on the relationship between   waitn  and   airn  regardless of   usesp ;  if  airn  is greater than or equal to   waitn , conditional expression 00 is false , and if   airn  is less than   waitn , conditional expression 00 is true . the agent instructs the aircraft to wait if the condition is found to be false but proceeds to evaluation of condition 01 if true . condition 01 is a second inequality containing an additional state variable   pushn  ( the number of aircraft pushing back in area ( i ) ) but , given that   pushn  is over the natural numbers , this condition ends up equivalent to the rule that   pushn  is greater than 0 . conditional expression 01 is then to wait if there are no aircraft pushing back and start otherwise . the solution for area ( ii ) changes according to the value of the variable   s9_pb  ( 1 when with pushback at spot 9 , 0 otherwise ) . because the range of this variable is 0 - 1 , condition 00 is always true . thus this is effectively equivalent to always pausing in the same manner as the previous solution ; it is a rule that passes through the line . this solution has the second best scores in both robustness and readability among the calculated pareto solutions . the probability that these rules do not generate the deadlock was determined by simulation to be 97.9% . such rules would have been hard to produce by hand but could be obtained by using the proposed method . as for the current control rules of narita airport , an aircraft is allowed to start when it reaches area ( i ) , and it is passed through the gate if there are no aircraft ramps in area ( ii ) . however , under this rule , according to the results of our simulation , the probability that deadlock does not occur is a mere 60% . the number of takeoffs and landings at narita airport continues to grow every year , and it is expected to increase by a factor of 1.5 from 2012 to 2015 . this rule was created in 2012 , when it was not possible to assume an increase in the departure and arrival numbers and has more recently had a low executable probability . enormous cost could result under the present method , but with the results of the proposed method , the chance of avoiding deadlock is about 99% , which shows the effectiveness of the proposed method . the experiments were carried out by comparing the proposed method and the one proposed in . we perform each approach over the course of 10 hours and compare the robustness and readability of the final solutions . figure 11 shows the maximum and minimum values and the average value of the readability of the final solution by the former method and the proposed method . we calculated the result of the t - test with a significance level of 5% with respect to readability distribution , and it was found that there is significant difference in the distribution , with the proposed method improved in readability as compared to the former method . figure 12 shows the maximum and minimum values and the average value of the robustness of the final solution by the former method and the proposed method . as before , we calculated the result of the t - test with a significance level of 5% with respect to robustness distribution , and it was found that there is again an improvement with respect to robustness through the proposed method . figures 13 and 14 plot the solution of each generation of the proposed and former method , respectively . from these figures , the proposed method has shown an improvement in both robustness and readability in parallel . in comparison , although a solution with high robustness given enough time does eventually come out through the former method , an improvement of readability does not . we have thus found that our method produces excellent results as compared to the former method . various systems in the real world can be modeled as multiagent systems , and appropriately designing the behavior rules of the agents is important to optimize the performance of that system . as important as improving robustness , however , is improving readability , for the sake of maintenance and operation . as a result , the object of the present study was the construction of a rule generation method that considers both robustness and readability . in this paper , we proposed a combined method of pado proposed by   and operator equalization proposed by   and used in , but incorporating metrics on readability from   on our the rule generation method . by applying our proposed method to the problem of aircraft ground control , it was possible to obtain a set of rules that was both more robust and more readable than the existing rule set , as well as when compared with the rules generated by the former method . in the future , we hope to consider further discussion about other aspects of readability . in this paper , it was in the form of applying indicators of readability on the source code of generated rules , but a different point of view , such as ease of understanding , may become necessary when talking about humans successfully using a set of rules to act . being able to calculate the usability of rules is necessary in order to further deepen such discussions .
rule generation method is proposed for an aircraft control problem in an airport . 
 designing appropriate rules for motion coordination of taxiing aircraft in the airport is important , which is conducted by ground control . 
 however , previous studies did not consider readability of rules , which is important because it should be operated and maintained by humans . 
 therefore , in this study , using the indicator of readability , we propose a method of rule generation based on parallel algorithm discovery and orchestration ( pado ) . by applying our proposed method to the aircraft control problem 
 , the proposed algorithm can generate more readable and more robust rules and is found to be superior to previous methods .