REINFORCEMENT learning (RL) [1] is a useful ap-proach to learning an optimal policy from sample behaviors of a controlled system with inherent stochasticity, e.g., a Markov decision process (MDP) [2], when the probabilities associated with the controlled system are unknown a priori. In RL, we use a reward function that assigns a reward to each transition in the behaviors and evaluate a control policy by the return, namely an expected (discounted) sum of the rewards along the behaviors. One of the recent trends is to apply RL to synthesis problems under linear temporal logic (LTL) constraints. LTL is a formal language with rich expressivity and thus suitable for describing complex missions assigned to a controlled system [3], [4]. It is known that any LTL formula can be converted into an -automaton with a B¨uchi or a Rabin acceptance condition [3]. In many studies on LTL synthesis problems using RL, reward functions are formed systematically from automata corresponding to the LTL specification. This direction was first investigated by Sadigh et al. [5], where they defined a reward function based on the acceptance condition of a deterministic Rabin automaton [3] that accepts all words satisfying the LTL
This work was partially supported by JST-ERATO HASUO Project Grant Number JPMJER1603, Japan, and JST-Mirai Program Grant Number JP-MJMI18B4, Japan. The work of A. Sakakibara was supported by the Grant-in-Aid for Japan Society for the Promotion of Science Research Fellow under Grant JP19J13487.
The authors are with the Graduate School of Engineering Science, Osaka University, Toyonaka 560-8531, Japan (e-mail: r-oura, sakakibara@hopf.sys.es.osaka-u.ac.jp; ushio@sys.es.osaka-u.ac.jp).
constraint. Reward functions defined on specification automata were also proposed for a deep reinforcement learning method [6] and for an extension of LTL in collaboration with a control barrier function [7].
Recently, a limit-deterministic B¨uchi automaton (LDBA) or a generalized one (LDGBA) is paid much attention to as an -automaton corresponding to the LTL specification [8]. The RL-based approaches to the synthesis of a control policy using LDBAs or LDGBAs have been proposed in [9]– [12]. An (LD)GBA has multiple accepting sets and accepts behaviors visiting all accepting sets infinitely often. One possible approach to generalized B¨uchi acceptance conditions is to convert a GBA into a corresponding BA, which has a single accepting set. The conversion, however, fixes the order of visits to accepting sets of the GBA [3] and causes the sparsity of the reward, which is a critical issue in RL-based controller synthesis. Another approach to RL-based synthesis for generalized B¨uchi conditions is the accepting frontier function introduced in [9], [10], based on which the reward function is defined. However, the function is memoryless, that is, it does not provide information of accepting sets that have been visited, which is important to improve learning performance.
In this letter, we propose a novel method to augment an LDGBA converted from a given LTL formula. Then, we define a reward function based on the acceptance condition of the product MDP of the augmented LDGBA and the controlled system represented as the MDP. The rest of the letter is organized as follows. Section II reviews an MDP, LTL, and automata. Section III proposes a novel RL-based method for the synthesis of a control policy. Section IV presents a numerical example to show the effectiveness of our proposed method.
Notations: is the set of nonnegative integers. is the set of nonnegative real numbers. For sets A and B, AB denotes their concatenation, i.e., . denotes the infinite concatenation of the set A and denotes the finite one, i.e., and , respectively. is the empty string.
A. Markov Decision Process
Definition 1: A (labeled) Markov decision process (MDP) is a tuple , where S is a finite set of states, A is a finite set of actions, is a transition probability function, is the initial state, AP is a finite set of atomic propositions, and is a labeling function that assigns a set of atomic propositions to each transition. Let S s.t. . Note that for any state and action .
In the MDP M, an infinite path starting from a state is defined as a sequence such that for any . A finite path is a finite sequence in . In addition, we sometimes represent as to emphasize that starts from . For a path , we define the corresponding labeled path . resp., is defined as the set of infi-nite (resp., finite) paths starting from in the MDP M. For each finite path denotes its last state.
Definition 2: A policy on an MDP M is defined as a mapping . A policy is a positional policy if for any and any , it holds that ; and there exists one action such that if , and for any with .
Let (resp., ) be the set of infi- nite (resp., finite) paths starting from in the MDP M under a policy . The behavior of the MDP M under a policy is defined on a probability space .
A Markov chain induced by the MDP M with a positional policy is a tuple , where for and such that . The state set of can be represented as a disjoint union of a set of transient states and closed irreducible sets of recurrent states with , i.e., [13]. In the following, we say a “recurrent class” instead of a “closed irreducible set of recurrent states” for simplicity.
In the MDP M, we define a reward function . The function returns the immediate reward received after the agent performs an action a at a state s and reaches a next state as a result.
Definition 3: For a policy on an MDP M, any state , and a reward function R, we define the expected discounted reward as
where denotes the expected value given that the agent follows the policy from the state s and is a discount factor. The function is often referred to as a state-value function under the policy . For any stateaction pair , we define an action-value function under the policy as follows.
Definition 4: For any state , a policy is optimal if
where is the set of positional policies over the state set S.
B. Linear Temporal Logic and Automata
In our proposed method, we use linear temporal logic (LTL) formulas to describe various constraints or properties and to systematically assign corresponding rewards. LTL formulas are constructed from a set of atomic propositions, Boolean operators, and temporal operators. We use the standard notations for the Boolean operators: (true), (negation), and (conjunction). LTL formulas over a set of atomic propositions AP are defined as
where , and are LTL formulas. Additional Boolean operators are defined as , and . The operators X and U are called “next” and “until”, respectively. Using the operator U, we define two temporal operators: (1) eventually, Fand (2) always, G.
Definition 5: For an LTL formula and an infinite path of an MDP M with , the satisfaction relation is recursively defined as follows.
where be the i-th suffix .
The next operator X requires that is satisfied by the next state suffix of . The until operator U requires that holds true until becomes true over the path . In the following, we write for simplicity without referring to M.
For any policy , the probability of all paths starting from on the MDP M that satisfy an LTL formula under the policy , or the satisfaction probability under is defined as
We say that an LTL formula is positively satisfied by a positional policy if
Any LTL formula can be converted into various -automata, namely finite state machines that recognize all infinite words satisfying . We review a generalized B¨uchi automaton at the beginning, and then introduce a limit-deterministic generalized B¨uchi automaton [10].
Definition 6: A transition-based generalized B¨uchi automaton (tGBA) is a tuple , where X is a finite set of states, is the initial state, is an input alphabet including is a set of transitions, and is an acceptance condition, where for each is a set of accepting transitions and called an accepting set. We refer to a tGBA with one accepting set as a tBA.
An infinite sequence is called an infinite run if for any . An infinite word is accepted by if and only if there exists an infinite run starting from such that for each , where inf(r) is the set of transitions that occur infinitely often in the run r.
Definition 7: A transition-based limit-deterministic generalized B¨uchi automaton (tLDGBA) is a tGBA B = such that X is partitioned into two disjoint sets and such that
• F, ..., n},
• |{(x, σ, x, x,
• |{(x, σ, x, x=0, ,
• x, ε, xδ, x .
An -transition enables the tLDGBA to change its state with no input and reflects a single “guess” from to . Note that by the construction in [8], the transitions in each part are deterministic except for -transitions from to . It is known that, for any LTL formula , there exists a tLDGBA that accepts all words satisfying [8]. We refer to a tLDGBA with one accepting set as a tLDBA. In particular, we represent a tLDGBA recognizing an LTL formula as , whose input alphabet is given by .
We introduce an automaton augmented with binary-valued vectors. The automaton can ensure that transitions in each accepting set occur infinitely often.
Let be a set of binary-valued vectors, and let 1 and 0 be the n-dimentional vectors with all elements 1 and 0, respectively. In order to augment a tLDGBA , we introduce three functions , and as follows. For any , where if and otherwise. For any if v = 1 and reset(v) = v otherwise. For any , where for any .
Each vector v is called a memory vector and represents which accepting sets have been visited. The function visitf returns a binary vector whose i-th element is 1 if and only if a transition in the accepting set occurs. The function reset returns the zero vector 0 if at least one transition in each accepting set has occurred after the latest reset. Otherwise, it returns the input vector without change.
Definition 8: For a tLDGBA , its augmented automaton is a tLDGBA , where , , is defined as
Clearly from the construction, we have for any . Thus, is a run of starting from .
We now show that for each . Since for each , we have
for any , where for each . By the construction of , therefore, there are infinitely many indices with . Let be arbitrary nonnegative integers such that , and for any . Then,
where is the j-th element of . Hence, we have for each , which implies .
: Consider any . Then, there exists a run of such that and for each , i.e.,
For the run , we construct a sequence r = such that and for any . It is clear that r is a run of B starting from . It holds by Eq. (1) that for each , which implies .
For example, shown in Figs. 1 and 2 are a tLDGBA and its augmented automaton, respectively, associated with the following LTL formula.
The acceptance condition F of the tLDGBA is given by F = , where and . Practically, states in a strongly connected component that contains no accepting transitions can be merged as shown in Fig. 2.
Fig. 1. The tLDGBA recognizing the LTL formula GF. Red arcs are accepting transitions that are numbered in accordance with the accepting sets they belong to.
Fig. 2. The augmented automaton for the tLDGBA in Fig. 1 recognizing the LTL formula GF, where the initial state is Red arcs are accepting transitions that are numbered in accordance with the accepting sets they belong to. All states corresponding to are merged into
We modify the standard definition of a product MDP to deal with -transitions in the augmented automaton.
Definition 9: Given an augmented tLDGBA and an MDP M, a tuple is a product MDP, where is the finite set of states; is the finite set of actions such that s.t. , where is the action for the -transition to the state is the initial state; is the transition probability function defined as
where sx, v)) and s, v; δ, a, sP , a) > 0} is
the set of transitions; and is the acceptance condition, where for each .
Definition 10: The reward function is defined as
, a, s
where is a positive value. Remark 1: When constructing a tBA from a tGBA, the order of visits to accepting sets of the tGBA is fixed. Consequently, the rewards based on the acceptance condition of the tBA tends to be sparse. The sparsity is critical in RL-based controller synthesis problems. The augmentation of the tGBA relaxes the sparsity since the augmented tGBA has all of the order of visits to all accepting sets of the original tGBA. For the acceptance condition F of the tGBA, the size of the state space of the augmented tGBA is about
the tBA constructed from the tGBA. However, the number of accepting transitions to all transitions in the augmented tGBA is much greater than the tBA. Therefore, our proposed method is expected to be better than using the tLDBA in the sense of sample efficiency. In the following, we show that a positional policy positively satisfying on the product MDP is synthesized by using the reward function R in Definition 10, which is based on the acceptance condition of .
with a positional policy , let be the set of states in , where is the set of transient states and is the recurrent class for each , and let be the set of all recurrent states in . Let be the set of transitions in a recurrent class , namely , and let : be the transition probability function under .
the Markov chain , satisfies one of the following conditions.
where is the probability that the transition reoccurs after it occurs in k time steps. Eq. (3) means that each transition in occurs infinitely often. However, the memory vector v is never reset in by the assumption. This directly contradicts Eq. (3).
Lemma 1 implies that, for an LTL formula , if a policy does not satisfy , then the agent obtains no reward in recurrent classes; otherwise there exists at least one recurrent class where the agent obtains rewards infinitely often.
Theorem 1: For a product MDP of an MDP M and an augmented tLDGBA corresponding to a given LTL formula and the reward function given by Definition 10, if there exists a positional policy positively satisfying on , then there exists a discount factor such that any algorithm that maximizes the expected discounted reward with will find a positional policy positively satisfying on .
Proof: Suppose that is an optimal policy but does not satisfy the LTL formula . Then, for any recurrent class in the Markov chain and any accepting set of the product MDP holds by Lemma 1. Thus, the agent under the policy can obtain rewards only in the set of transient states. We consider the best scenario in the assumption. Let be the probability of going to a state in k time steps after leaving the state s, and let be the set of states in recurrent classes that can be transitioned from states in by one action. For the initial state in the set of transient states, it holds that
where the action a is selected by . By the property of the transient states, for any state in , there exists a bounded positive value m such that [13]. Therefore, there exists a bounded positive value such that . Let be a positional policy satisfying . We consider the following two cases.
1) Assume that the initial state is in a recurrent class for some . For any accepting set , holds by the definition of . The expected discounted reward for is given by
where the action a is selected by . Since is in , there exists a positive number [13]. We consider the worst
whereas all states in are positive recurrent because [14]. Obviously, holds for any by the Chapman-Kolmogorov equation [13]. Furthermore, we have by the property of irreducibility and positive recurrence [15]. Hence, there exists such that for any and we have
Therefore, for any and any , there exists such that implies
2) Assume that the initial state is in the set of transient states . holds by the definition of . For a recurrent class such that for each accepting set , there exist a number , a state in , and a subset of transient states such that , for , and by the property of transient states. Hence, it holds that for the state . Thus, by ignoring rewards in , we have
In this section, we apply the proposed method to a path planning problem of a robot in an environment consisting
Fig. 3. The environment consisting of eight rooms and one corridor. Red arcs are the transitions that we want to occur infinitely often, while blue arcs are the transitions that we never want to occur. is the initial state.
of eight rooms and one corridor as shown in Fig. 3. The state is the initial state and the action space is speci-fied with A(s) = {Right, Left, Up, Down} for any state and , where means attempting to go to the state for . The robot moves in the intended direction with probability 0.9 and it stays in the same state with probability 0.1 if it is in the state . In the states other than , it moves in the intended direction with probability 0.9 and it moves in the opposite direction to that it intended to go with probability 0.1. If the robot tries to go to outside the environment, it stays in the same state. The labeling function is as follows.
In the example, the robot tries to take two transitions that we want to occur infinitely often, represented by arcs labeled by {a} and {b}, while avoiding unsafe transitions represented by the arcs labeled by {c}. This is formally specified by the LTL formula given by Eq. (2). The LTL formula requires the robot to keep on entering the two rooms and from the corridor regardless of the order of entries, while avoiding entering the four rooms , and . The tLDGBA and its augmented automaton are shown in Figs. 1 and 2, respectively.
Through the above scenario, we compare our approach with 1) a case where we first convert the tLDGBA into a tLDBA, for which the augmentation makes no change, and thus a reward function in Definition 10 is based on a single accepting set; and 2) the method using a reward function based on the accepting frontier function [9], [10]. For the three methods, we use Qlearning1 with an epsilon-greedy policy. The epsilon-greedy parameter is given by , where is the number of visits to state within t time steps [16], so that it will gradually reduce to 0 to learn an optimal policy asymptotically. We set the positive reward and the discount factor . The learning rate varies in accordance with the Robbins-Monro condition. We train the agent in 10000 iterations and 1000 episodes for 100 learning sessions.
Fig. 4. The mean of average reward in each episode for 100 learning sessions obtained from our proposed method (left) and the method using tLDBA (right). They are plotted per 50 episodes and the green areas represent the range of standard deviations.
Fig. 5. The optimal policy obtained from our proposed method (left) and the method in [9], [10] (right).
Results
1) Fig. 4 shows the average rewards obtained by our proposed method and the case using a tLDBA converted from , respectively. Both methods eventually acquire an optimal policy satisfying . As shown in Fig. 4, however, our proposed method converges faster. This is because the order of entrances to the rooms and is determined according to the tLDBA. Moreover, the number of transitions with a positive reward in is larger than that in .
2) We use the accepting frontier function [9], [10] for the tLDGBA . Initializing a set of transitions F with the set of the all accepting transitions in , the function receives the transition that occurs and the set F. If is in F, then Acc removes the accepting sets containing from F. For the product MDP of the MDP M and the tLDGBA , the reward function is based on the removed sets of . Then, we synthesize a positional policy on the product MDP derived from the tLDGBA .
Fig. 5 shows the optimal policies obtained by our proposed method and the method in [9], [10], respectively. The policy obtained by the method with the accepting frontier function fails to satisfy the LTL specification2 because it is impossible that the transitions labeled with {a} and {b} occur from infinitely often by any positional policy on the product MDP with shown in Fig. 1. More specifically, as shown in Fig. 6, the agent cannot take each accepting transition colored with red by any deterministic stationary action selection at the product state . In our proposed method, the augmented tLDGBA includes the information of the (path-dependent) domain of the accepting frontier function explicitly as memory vectors, which enables us to synthesize a positional policy satisfying on the product MDP.
Fig. 6. A part of the product MDP of the MDP shown in Fig. 3 and the tLDGBA shown in Fig. 1, where each transitions is labeled with either an action or the transition probability.
The letter proposed a novel RL-based method for the synthesis of a control policy for an LTL specification using an augmented tLDGBA. The proposed method improved the learning performance compared to existing methods. It is future work to investigate a method that maximizes the satisfaction probability.
[1] R. S. Sutton and A. G. Barto, Reinforcement Learning: An Introduction, 2nd ed. MIT press, 2018.
[2] M. L. Puterman, Markov Decision Processes: Discrete Stochastic Dynamic Programming. John Wiley & Sons, 1994.
[3] C. Baier and J.-P. Katoen, Principles of Model Checking. MIT Press, 2008.
[4] C. Belta, B. Yordanov, and E. Aydin Gol, Formal Methods for DiscreteTime Dynamical Systems. Springer International Publishing, 2017.
[5] D. Sadigh, E. S. Kim, S. Coogan, S. S. Sastry, and S. A. Seshia, “A learning based approach to control synthesis of Markov decision processes for linear temporal logic specifications,” in 53rd IEEE Conference on Decision and Control, 2014, pp. 1091–1096.
[6] Q. Gao, D. Hajinezhad, Y. Zhang, Y. Kantaros, and M. M. Zavlanos, “Reduced variance deep reinforcement learning with temporal logic specifications,” in Proceedings of the 10th ACM/IEEE International Conference on Cyber-Physical Systems, 2019, pp. 237–248.
[7] X. Li, Z. Serlin, G. Yang, and C. Belta, “A formal methods approach to interpretable reinforcement learning for robotic planning,” Science Robotics, vol. 4, no. 37, pp. 1–16, 2019.
[8] S. Sickert, J. Esparza, S. Jaax, and J. Kˇret´ınsk`y, “Limit-deterministic B¨uchi automata for linear temporal logic,” in International Conference on Computer Aided Verification. Springer, 2016, pp. 312–332.
[9] M. Hasanbeig, A. Abate, and D. Kroening, “Logically-constrained reinforcement learning,” arXiv preprint arXiv:1801.08099, 2018.
[10] M. Hasanbeig, Y. Kantaros, A. Abate, D. Kroening, G. J. Pappas, and I. Lee, “Reinforcement learning for temporal logic control synthesis with probabilistic satisfaction guarantees,” arXiv preprint arXiv:1909.05304, 2019.
[11] E. M. Hahn, M. Perez, S. Schewe, F. Somenzi, A. Trivedi, and D. Wojtczak, “Omega-regular objectives in model-free reinforcement learning,” in International Conference on Tools and Algorithms for the Construction and Analysis of Systems. Springer, 2019, pp. 395–412.
[12] A. K. Bozkurt, Y. Wang, M. M. Zavlanos, and M. Pajic, “Control synthesis from linear temporal logic specifications using model-free reinforcement learning,” arXiv preprint arXiv:1909.07299, 2019.
[13] R. Durrett, Essentials of Stochastic Processes, 2nd ed. Springer, 2012.
[14] L. Breuer, “Introduction to stochastic processes,” [Online]. Available: https://www.kent.ac.uk/smsas/personal/lb209/files/sp07.pdf.
[15] S. M. Ross, Stochastic Processes, 2nd ed. Wiley New York, 1996.
[16] S. Singh, T. Jaakkola, M. L. Littman, and C. Szepesv´ari, “Convergence results for single-step on-policy reinforcement-learning algorithms,” Machine learning, vol. 38, no. 3, pp. 287–308, 2000.