You are viewing a javascript disabled version of the site. Please enable Javascript for this site to function properly.
Go to headerGo to navigationGo to searchGo to contentsGo to footer
In content section. Select this link to jump to navigation

Quantitative reachability analysis of generalized possibilistic decision processes

Abstract

The verification of reachability properties of fuzzy systems is usually based on the fuzzy Kripke structure or possibilistic Kripke structure. However, fuzzy Kripke structure or possibilistic Kripke structure is not enough to describe nondeterministic and concurrent fuzzy systems in real life. In this paper, firstly, we propose the generalized possibilistic decision process as the model of nondeterministic and concurrent fuzzy systems, and deduce the possibilities of sets of paths of the generalized possibilistic decision process relying on defining of schedulers. Then, we give fuzzy matrices calculation methods of the maximal possibilities and the minimal possibilities of eventual reachability, always reachability, constrained reachability, repeated reachability and persistent reachability. Finally, we propose a model checking approach to convert the verification of safety property into the analysis of reachabilities.

1Introduction

Reachability is one of the central problems in model checking [1, 2], program analysis [3] and verification [3], which is about whether a system in one state can reach other states [4]. Reachability is widely applied in urban transportation planning [5], human geography [6], regional economics and computer science [7]. In computer science, the use of reachability decision algorithm can avoid the loop detection of useless state space, save the memory occupancy of the algorithm, and improve the efficiency of the algorithm. The use of reachability optimization algorithm can reduce the complexity of the algorithm [4].

Reachability problems based on classical model checking were proposed in [8]. Classical model checking is to verify the qualitative characteristics of systems. Qualitative reachability aims at obtaining exact values of certain events. However, in real life, there are some randomness, uncertainties and inconsistencies that classical model checking is unable to express those information, such as a 90 percent probability of system crashing during operation [12]. To embed uncertain information into the model, quantitative model checking methods are proposed, including probabilistic model checking [4], possibilistic model checking [17], multi-valued model checking [10, 11] and fuzzy model checking [2, 12, 18]. Quantitative reachability is to find the maximum probability [15] or the minimum probability [17].

This study extends the existing approach of the generalized possibilistic model checking to solve quantitative reachability problems. Li et al. [17, 18, 21, 22] propose possibility model checking and generalized possibility model checking in combination with measure theory, providing a solution to the existing problems of fuzzy model checking. In [17], the reachability problem based on the possibility measure is given, and the reachability problem is studied by using a possibilistic Kripke structure as a model. In the existing generalized possibilistic model checking, the generalized possibilistic Kripke structure is the formal model of the representation system. However, in real life, a generalized possibilistic Kripke structure cannot describe nondeterministic and concurrent fuzzy systems. For example, suppose three experts want to formulate treatment schemes for a new bacterial infection. Because different experts have different understandings of the disease, each expert has a different scheme. We use a generalized possibilistic Kripke structure K = (S, P, I, AP, L) to model the patient’s treatment process respectively. It describes three treatment processes given by three experts(α, β and γ) in Fig. 1(a), (b), (c). However, the generalized possibilistic Kripke structure can only represent the treatment process of a scheme, cannot describe the joint consultation of three experts. To solve this problem, we propose a generalized possibilistic decision process to describe nondeterministic and concurrent fuzzy systems.

Fig. 1

The treatment processes by three experts.

The treatment processes by three experts.

In this paper, we mainly study the maximum possibility problem and the minimum possibility problem of eventual reachability, always reachability, constrained reachability, repeated reachability and persistent reachability under a generalized possibilistic decision process. First, we propose a generalized possibilistic decision process as the model of fuzzy systems, and deduce the possibilities of sets of paths of the generalized possibilistic decision process relying on defining of schedulers. Then, the fuzzy matrices calculation methods of reachability under a generalized possibilistic decision process are given, and the results show that the compositional operations of the fuzzy matrix is polynomial time, which is better than the exponential level of existing algorithms. Finally, we propose a model checking approach to convert the verification of safety property into the analysis of reachabilities.

The rest of this paper is organized as follows. Section 2 gives some preliminary knowledge. In Section 3, we define generalized possibilistic decision processes. In Section 4, we present the fuzzy matrix calculation methods of eventual reachability, always reachability, constrained reachability, repeated reachability and persistent reachability. In Section 5, we use good prefixes to analyze the possibilistic regular safety property. Section 6 is the conclusion. The proofs of some theorems of this paper can be found in the Appendix.

2Preliminaries

To model and verify fuzzy systems, we provide some necessary knowledge including the fuzzy set, fuzzy set operation, fuzzy matrix operation, closure and others.

Definition 1. [23, 26, 28] Let X be a universal set. A fuzzy set A of X is an function which associates each element in X a value in the interval [0, 1], i.e., A : X ⟶ [0, 1]. For x ∈ X, A (x) is the membership of x in the fuzzy set A.

We use F(X) to represent all fuzzy sets in X, i.e. F(X)={AA:X[0,1]} .

Definition 2. [26, 28] Let A,BF(X) , we use A ∪ B, A ∩ B, Ac to represent the union, intersection and complement of A and B. The definition is as follows.

(A ∪ B) (x) = A (x) ∨B (x) = max {A (x) , B (x)},

(A ∩ B) (x) = A (x) ∧ B (x) = min {A (x) , B (x)},

Ac (x) =1 - A (x).

Definition 3. [26, 28] Let R and S be two fuzzy matrices with m rows and n columns, i.e., R = (rijm×n, S = (sijm×n. We introduce some set operators.

R = S if and only if rij = sij for all i, j.

R ⊆ S if and only if rij ≤ sij for all i, j.

R ∪ S = (rijsijm×n.

R ∩ S = (rij ∧ sijm×n.

Rc = (1 - rijm×n.

Definition 4. [26, 28] Let R be a fuzzy matrix with m rows and n columns, and S be a fuzzy matrix with n rows and l columns, i.e., R = (rijm×n and S = (sijn×l. The composition operation of R and S is R ∘ S = (tijm×l, where tij=k=1n(rikskj) , (i = 1, 2, . . . , m, j = 1, 2, . . . , l). For fuzzy matrices R, S, T, the composition operation has some operation laws.

(R ∘ S) ∘ T = R ∘ (S ∘ T);

(R ∪ S) ∘ T = (R ∘ T) ∪ (S ∘ T).

Let X be a universal set. For the fuzzy matrix R = (R (s, t)) s,tX, we use R+ to denote its transitive closure. When X is finite, X has ∣X∣ elements, then R+ = R ∪ R2 ∪ . . . ∪ RX, where Rk+1 = Rk ∘ R for any positive integer number k. The Kleene closure R* = R0 ∪ R+, for each 1≤ s, t ≤ ∣ S ∣, R0(s,t)={1s=t0st . Possibility measure theory is a kind of uncertainty theory, which mainly deals with incomplete information and uncertain information. In addition, the possibility measure does not require additivity, which is more applicable to deal with practical application systems.

Definition 5. [21] Let X be a nonempty set,and Ω be a set composed of some subsets of X elements. We call Ω a σ- algebra, which is closed to countable and take complement set operations. The possibility measure on σ- algebra Ω is a mapping POS : Ω → [0, 1], which satisfies the following conditions:

  • (1) POS (∅) =0;

  • (2) POS (X) =1;

  • (3) If Ei ∈ Ω, i ∈ I, then POS (⋃ iIEi) = ⋁ iIPOS (Ei)

If only conditions(1) and (3), then POS is called a generalized possibility measure. If POS is a generalized possibility measure on the power set 2X,for any A ⊆ X, there is POS (A) = ⋁ aAPOS ({ a }).

Definition 6. [18] A Generalized possibilistic Kripke structure(GPKS) is a tuple K = (S, P, I, AP, L), where

  • (1) S is a countable, nonempty states set;

  • (2) P : S × S → [0, 1] is the possibility transition distribution, for any states s, there is a state t, such that P (s, t) >0;

  • (3) I : S → [0, 1] is the initial distribution of possibility and there exists a state s, such that I (s) >0;

  • (4) AP is a set of atomic propositions;

  • (5) L : S × AP → [0, 1] is a label function, L (s, a) represents the true value of atomic proposition a in state s.

For a GPKS K, its path is defined as an infinite sequence of states π = s0s1s2 ⋯ ∈ s ω, for any i so that P (si, si+1) > 0. Let Paths (s) and Pathsfin (s) represent the set of all infinite and finite paths starting from the state s in K. Paths (K) represents the set of all infinite paths in K, Pathsfin (K) represents the set of all finite paths in K, such as πˆ=s0s1sn .

3Generalized Possibilistic Decision Processes

In this section, first, we give the notion of the generalized possibilistic decision process. Then, the definition of the scheduler is proposed to solve the nondeterministic the generalized possibilistic decision process. Finally, we solve the generalized possibility measure problems.

Nondeterminism is absent in GPKS. The generalized possibilistic decision process can be viewed as a variant of GPKS that permits both possibility and nondeterministic choices.

Definition 7. A Generalized possibilistic decision process(GPDP) is a tuple M = (S, Act, P, I, AP, L), where

  • (1) S is a countable, nonempty set of states;

  • (2) Act is a set of actions;

  • (3) P : S × Act × S ⟶ [0, 1] is a function, called possibilistic transition distribution function. For all states s ∈ S and actions α ∈ Act, there exits a state t ∈ S such that P (s, α, t) >0;

  • (4) I : S ⟶ [0, 1] is a function, there exits states s such that I (s) >0;

  • (5) AP is a set of atomic propositions;

  • (6) L : S × AP ⟶ [0, 1] is a possibilistic labeling function, which can be viewed as function mapping a state s to the fuzzy set of atomic proposition, L (s, a) denotes the possibility or truth value of atomic proposition a which is hold in state s.

Furthermore, if the set S, Act and AP are finite sets, then M is a finite GPDP. In this paper, we always assume that GPDPs are finite. A GPDP has a unique initial distribution I. For all states s, t ∈ S and actions α ∈ Act, P (s, α, t) denotes the possibility from state s under action α to state t. An action α is enabled in state s if and only if P (s, α, t) >0. Let Act (s) denote the set of enabled actions in state s. For any states s ∈ S, it is required that Act (s)≠ ∅. Each state t for which P (s, α, t) >0 is called an α-successor of s. The set of direct α-successors of s is defined as:

Post (s, α) = {t ∈ S ∣ P (s, α, t) >0}.

The set of α-predecessors of s is defined by:

Pref (t) = {(s, α) ∈ S × Act (s) ∣ P (s, α, t) >0}.

We also use the P (s, α, T) to denote the possibility from the state s under the action α to the set T of states, that is, P(s,α,T)=tTP(s,α,t) .

Paths in GPDP M are defined as infinite alternating sequences π = s0 α0s1 α1s2 ⋯ ∈ (S × Act ω such that P (si, αi, si+1) >0 for all i ∈ I. Paths (s) denotes the set of all infinite paths in M that start in state s. Similarly, Pathsfin (s) denotes the set of all finite path fragment πˆ=s0α0s1α1s2αn-1sn such that s0 = s. Paths (M) and Pathsfin (M) denote the set of all infinite paths and finite paths in M respectively. The trace of the infinite path fragment π = s0s1⋯ is defined as trace(π) = L (s0L (s1)⋯. The trace of the finite path fragment πˆ = s0s1 ⋯ sn is defined as trace( πˆ ) = L (s0L (s1) ⋯ L (sn).The trace of all infinite paths starting from state s is defined as Traces (s) = trace (Paths (s)).

Example 1. Let us consider the joint consult of three experts in Fig. 1(a),(b),(c). GPDP in Fig. 2, where states are represented by ovals and transitions by labeled edges, and state names are depicted outside the ovals, and labeling functions of the states are depicted inside the ovals.

Fig. 2

Treatment process by three experts’ consultation.

Treatment process by three experts’ consultation.

S = {s0, s1, s2}, Act = {α, β, γ}, AP = {P, G, E}. P indicates that the patient is in “bad” health, G indicates that the patient is in “normal” health, E indicates that the patient is in “good” health.

For state s0, the labeling functions are L (s0, P) =0.85, L (s0, G) =0.3, L (s0, E) =0.2. L (s0, P) =0.85 indicates that the degree of “bad” health in state s0 is 0.85. L (s0, G) =0.3 indicates that the degree of “normal” health in state s0 is 0.3. L (s0, E) =0.2 indicates that the degree of “good” health in state s0 is 0.2.

Act (s0) = {α, β, γ}, α, β, γ indicate three different treatment indicates. P (s0, α, s1) =0.8 indicates that the possibility of the patient’s health condition changing from state s0 to state s1 is 0.8 after the expert treats the patient with the α scheme. P (s0, β, s1) =0.6 indicates that the possibility of the patient’s health condition changing from state s0 to state s1 is 0.6 after the expert treats the patient with the β scheme. P (s0, γ, s1) =0.2 indicates that the possibility of the patient’s health condition changing from state s0 to state s1 is 0.2 after the expert treats the patient with the γ scheme.

Post (s0, α) = {s1, s2},Post (s0, α) indicates all α successor states of state s0.

Pref (s0) = {(s0, α) , (s0, β) , (s0, γ) , (s1, α) , (s1, β) ,

(s1, γ) , (s2, α) , (s2, β) , (s2, γ)}, Pref (s0) indicates all predecessor states of state s0.

For the GPKS K, 2Paths (K) is the algebra that is generated by {Cyl(πˆ)πˆPathsfin(K)} on Paths (K), but the GPDPs are not augmented with a unique possibility measure. Instead, deducing possibilities of sets of paths of a GPDP rely on the resolution of nondeterminism. This resolution is performed by a scheduler. A scheduler chooses in any state s one of the actions set A ⊆ Act (s). It does not impose any constraint on the possibilistic choice that is resolved once α ∈ A has been chosen.

Definition 8. Let M = (S, Act, P, I, AP, L) be a GPDP, a scheduler for M is a function Adv : S → 2Act so that Adv (s) ⊆ Act (s) for all s ∈ S.

Let PathAdv (s) and PathAdvfin(s) denote the set of paths and finite paths from state s under the decision of the scheduler Adv. Let PathAdv (M) and PathAdvfin(M) denote the set of paths and finite paths in the M under the decision of the scheduler Adv.

Given a GPDP M = (S, Act, P, I, AP, L), α ∈ Act, possibility distribution function P : S × α × S ⟶ [0, 1] can be represented by a fuzzy matrix. For convenience, the fuzzy matrix is written as P α, so that P α = (P (s, α, t)) s,tS, called the fuzzy transition matrix of M corresponding to scheduler α. Using fuzzy matrix Pmax=i=0nPαi denotes the maximal possibility transition matrix, so that

(Pmax(s,t))s,tS=(αAct(s)P(s,α,t))s,tS.

Using fuzzy matrix Pmin=i=0nPαi denotes the minimal possibility transition matrix, so that

(Pmin(s,t))s,tS=(αAct(s)P(s,α,t))s,tS.

The GPKS Kmax = (S, Pmax, I, AP, L) and GPKS Kmin = (S, Pmin, I, AP, L) can be constructed from matrix Pmax and matrix Pmin respectively.

Example 2. As shown in the Example 1, the order of s0s1s2 is used to give the fuzzy matrices Pmax, Pmin.

Pmax=(0.80.80.30.50.90.90.30.81),Pmin=(0.30.20.10.20.50.40.10.30.7),

Figure 3(a) shows the result of GPKS with respect to Pmax, Fig. 3(b) shows the result of GPKS with respect to Pmin.

Fig. 3

GPKS with respect to Pmax and Pmin.

GPKS with respect to Pmax and Pmin.

Although GPDP can induce the GPKS under the consideration of some schedulers, reasoning about quantitative reachabilities requires a formalization of the possiblities for sets of paths. This formalization is based on possibility measure theory, in particular possibility spaces and generalized possibility measure theory.

Definition 9. Given a GPDP M = (S, Act, P, IAP, L), let πˆ=s0α0s1α1s2sn-1αn-1snPathsAdvfin(M) and Adv be the max or min scheduler, then the cylinder set of the finite path πˆ is defined as

(1)
Cyl(πˆ)={πPathsAdv(M)πˆPref(π)},
where Pref(πAdv)={πPathsAdvfin(M)π is a finite prefix of πAdv }.

The cylinder set spanned by the finite adv-paths πˆ consists of all infinite adv-paths that start with πˆ . The cylinder sets serve as basis events of the measurable space 2PathsAdv (M) associated with M. From classical concepts of possibility theory, it follows that there exists a unique possibility measure GPo on the measurable space 2PathsAdv (M) associated with M are given by Definition 10.

Definition 10. Let M = (S, Act, P, I, AP, L) be a finite GPDP, the function GPo : PathsAdv  (M) → [0, 1] is defined:

(2)
GPo(π)=I(s0)i0PAdv(si,αi,si+1)
where π = s0 α0s1 α1s2 ⋯ ∈ PathsAdv  (M). Furthermore, we define GPo (E) = ∨ {GPo (π) ∣ π ∈ E}, for any E ⊆ PathsAdv  (M), then we have the function
GPo:2PathsAdv(M)[0,1]
is called generalized possibility measure over Ω = 2PathsAdv(M).

Let M = (S, Act, P, I, AP, L) be a GPDP, s ∈ S, αi ∈ Act, i ≥ 0, rAdv : S → [0, 1] is defined as following, which denotes the maximal possibility measure of all Adv-paths from state s in M:

(3)
rAdv(s)={i0PAdv(si,αi,si+1)s1=s,siS,αiAct}.

The role of the function rAdv is to help us to calculate the possibility of Adv-paths in M. The Theorem 1 gives a fuzzy matrix calculation for rAdv.

Theorem 1. Let M = (S, Act, P, I, AP, L) be a finite GPDP, for any s ∈ S, then we have

(4)
rAdv(s)=tS(PAdv+(s,t)PAdv+(t,t)).

In particular, PAdv is normal iff rAdv (s) =1 for any states s.

Proof. See the Appendix. ■

In the matrix notation, we have

(5)
rAdv=PAdv+DAdv,
where DAdv=(PAdv+(t,t))tS .

The computational complexity of rAdv (s) mainly depends on the time of computational possibilistic transition closure. In [30], they gave an optimal algorithm to calculate DAdv, then we get the time complexity is O (n2 log n), where n =∣ S ∣.

Example 3. According to the GPDP in Example 1, for any s ∈ S,rmax and rmin corresponding to the maximal scheduler and the minimal scheduler are given,

Dmax=(Pmax+(t,t))tS=(0.80.91),

Dmin=(Pmin+(t,t))tS=(0.30.50.7),

Pmax+=(0.80.80.80.50.90.90.50.81),

Pmin+=(0.30.20.10.20.50.40.20.30.7),

rmax=Pmax+Dmax=(0.80.91),

rmin=Pmin+Dmin=(0.30.50.7).

rmax (s0) = 0.8 indicates that under the maximal possibilistic scheduler. That is, the maximal possibility of all paths from state s0 is 0.8. rmin (s0) = 0.3 indicates that under the minimal possibilistic scheduler. That is, the maximal possibility of all paths from state s0 is 0.3.

Based on Theorem 1, we can get Theorem 2, which can convert the calculation of generalized possibility measure of infinite path into the calculation of possibility of finite path.

Theorem 2. Let M = (S, Act, P, I, AP, L) be a finite GPDP, then the generalized possibility measure of cylinder set πˆ=s0α0s1α1αn-1snPathsAdvfin(M) is:

GPo(Cyl(s0α0s1α1αn-1sn))=I(s0)i=0n-1PAdv(si,αi,si+1)rAdv(sn),
where GPo (Cyl (s0)) = I (s0) ∧ rAdv (s0).

Proof. See the Appendix. ■

Example 4. According to the GPDP of Example 1, under the maximal possibilistic measure and the minimal possibilistic measure, the generalized possibility measure of cylinder set of the corresponding finite paths πˆ=s0α0s1 are respectively,

GPomax  (Cyl (s0 α0s1) ) = I (s0) ∧ P (s0, α0, s1) ∧ rmax (s1) =0.8

GPomin (Cyl (s0 α0s1)) = I (s0) ∧ P (s0, α0, s1) ∧ rmin (s1) =0.5 .

4Reachability possibility

A typical task for the quantitative analysis of GPDPs is to compute the minimum or the maximum possibility for some reachabilities under consideration of the min or the max scheduler. This corresponds to the worst-case or best-case analysis possibility of GPDPs. Let M = (S, Act, P, I, AP, L) be a finite GPDP and B be a fuzzy set of target states. The fuzzy set B may represent a set of certain bad states that shall be visited with the minimum possibility, or dually, a set of good states that shall be visited with the maximum possibility. Some special reachabilities, such as eventual reachability, always reachability, constrained reachability, repeated reachability and persistent reachability are considered in this section.

4.1Eventual reachability possibility

The event of eventual reachability is denoted ⋄B. We use the B : S ⟶ [0, 1] to represent the fuzzy set of states. For the given GPDP M, π = s0 α0s1 … ∈ PathsAdv (M), ∣∣ ⋄ B (π) ∣∣ = ⋁ i≥0B (si). In the following, the quantitative analysis of eventual reachability is reduced to the computational range for all strategies of the minimum possibility or the maximum possibility of reaching a certain fuzzy set B of states.

Theorem 3. Let M = (S, Act, P, I, AP, L) be GPDP, max corresponds to the maximum possibility scheduler and min corresponds to the minimum possibility scheduler, then we have,

(6)
GPomax(B)=(GPomax(sB))sS=Pmax*DBrmax
(7)
GPomin(B)=(GPomin(sB))sS=Pmin*DBrmin.

Proof. See the Appendix. ■

Example 5. Based on Example 1, consider the generalized possibility of eventual reachability event ⋄B. Given the fuzzy state set B = (L (s, E)) sS, only the reachability of the event is considered at this point, and the label function in GPDP is defined as L (s) = s.

Then the maximum possibility and minimum possibility of the event ײ7B are respectively:

GPomax(7B)=Pmax*DBrmax=(0.80.80.30.50.90.90.30.81)*(0.20000.70000.9)(0.80.91)=(0.80.90.9),
GPomin(7B)=Pmin*DBrmin=(0.20.30.10.20.50.40.10.30.7)*(0.20000.70000.9)(0.30.50.7)=(0.20.50.7).

GPo (◊ ≤7B) ∥ max (s0) =0.8, ∥GPo (◊ ≤7B) ∥ min (s0) =0.2, describes the maximum possibility of the patient’s final state of health with the maximum possibility of “good” being 0.8 and the minimum possibility of “good” being 0.2 after 7 days treatment from state s0.

4.2Always reachability possibility

The event of always reachability is denoted B. We use the B : S ⟶ [0, 1] to represent the fuzzy set of states. For the given GPDPs M, π = s0 α0s1 … ∈ PathsAdv (M), B (π) = ⋀ i≥0B (si). Under the maximum possibility scheduler and the minimum possibility scheduler, the methods of computing GPomax (s ⊨  B) and GPomin (s ⊨  B) are given.

GPomax(sB)=πPathsmax(s)(GPomax(π)B(π))=(GPomax(π)i0B(si)).

GPomin(sB)=(GPomin(π)B(π))=(GPomin(π)i0B(si)).
From these results, we can get Theorem 4.

Theorem 4. Let M = (S, Act, P, I, AP, L) be GPDPs, max corresponds to the maximum possibility scheduler and min corresponds to the minimum possibility scheduler, then we have

(8)
GPomax(B)=(GPomax(sB))sS=υ.fBmax(Z),

(9)
GPomin(B)=(GPomin(sB))sS=υ.fBmin(Z),
where fBmax (Z) = B ∧ Pmax ∘ DZ ∘ rmax, fBmin (Z) = B ∧ Pmin ∘ DZ ∘ rmin. υ . fBmax (Z) denotes the maximal fixed point of the operator fBmax (Z), υ . fBmin (Z) denotes the minimal fixed point of the operator fBmin (Z). We have given the solution to the maximal fixed point in [18].

Example 6. Based on Example 1, considering the generalized possibility of always reachability event ≤7B, given the fuzzy state set B = (L (s, E)) sS, the maximum possibility and the minimum possibility of all states in M satisfied event ≤7B are respectively:

f(z1)=B(Pmax*Dzrmax)=(0.20.70.9)((0.80.80.30.50.90.90.30.81)*(100010001)(0.80.91))=(0.20.70.9),

f(z2)=B(Pmax*Dz1rmax)=(0.20.70.9)((0.80.80.30.50.90.90.30.81)*(0.20000.70000.9)(0.80.91))=(0.20.70.9).

fz1 = fz2, Therefore

GPomax(7B)=(0.20.70.9).

Calculated by the same method

GPomin(7B)=(0.20.50.7).

GPo ( E) ∥ max (s0) =0.2, ∥GPo ( E) ∥ min (s0) =0.2, it shows that from the state s0, it is unlikely that the patient’s health always reach a “good” state after 7 days treatment.

4.3Constrained reachability possibility

Given the GPDP M, and the fuzzy set of states B, C : S ⟶ [0, 1], consider the event of reaching B via a finite path fragment which ends in a fuzzy state s ∈ B, and visits only states in set of fuzzy states C prior to reaching fuzzy states s. This event is denoted by C ⊔ B. For n ≥ 0, the event C ⊔ nB has the same meaning as C ⊔ B, and it is required to reach B(via fuzzy state C) within n steps. Formally, C ⊔ nB is the union of the basic cylinders spanned by path fragments s0 α0s1 α1 … sm so that m ≤ n with possibility C (si) for all 0 ≤ i ≤ m with possibility B (sk).

Theorem 5. Let M = (S, Act, P, I, AP, L) be a GPDP, max corresponds to the maximum possibility scheduler and min corresponds to the minimum possibility scheduler, then we have

(10)
GPomax(CnB)=GPomax(CB)=(DCPmax)*DBrmax,

(11)
GPomin(CnB)=GPomin(CB)=(DCPmin)*DBrmin.

Proof. See the Appendix. ■

Example 7. Next, let’s continue to consider the generalized possibility of constrained reachability event C ⊔ ≤7B. Given the fuzzy state set C = (L (s, P)) sS, B = (L (s, E)) sS, the maximum possibility and the minimum possibility of all states in M satisfying event C ⊔ ≤7B are respectively:

GPomax(C7B)=(DCPmax)*DBrmax=((0.850000.40000.1)(0.80.80.30.50.90.90.30.81))*(0.20000.70000.9)(0.80.91)=(0.70.70.9),
GPomin(C7B)=(DCPmin)*DBrmin=((0.850000.40000.1)(0.20.30.10.20.50.40.10.30.7))*(0.20000.70000.9)(0.30.50.7)=(0.20.50.7).

GPo (C ⊔ ≤7B) ∥ max (s0) =0.7, ∥GPo (C ⊔ ≤7B) ∥ min (s0) =0.2, which means that the patient’s health state is “bad” at state s0. The expert adopts three treatment schemes after 7 days of treatment, the maximum possibility that a patient’s health states turn to a “good” state is 0.7 and the minimum possibility is 0.2.

4.4Repeated reachability possibility

Let B : S ⟶ [0, 1] be a set of fuzzy states in GPDPs. For the event of repeated reachability, we use the  ⋄ B to represent it. The set of all paths that visit B infinitely is often measurable. Let π = s0 α0s1 α1 … ∈ PathsAdv (s0), B(π)=i0jiB(sj) , we consider to calculate GPoAdv (s ⊨  ⋄ B) under a scheduler Adv.

Theorem 6. Let M = (S, Act, P, I, AP, L) be a GPDP, B : S ⟶ [0, 1] is a fuzzy set of states, max corresponds to the maximum possibility scheduler and min corresponds to the minimum possibility scheduler. Then we have,

(12)
GPomax(B)=Pmax+diag(Pmax+(t,t))tSB,

(13)
GPomin(B)=Pmin+diag(Pmin+(t,t))tSB.

Proof. See the Appendix. ■

Example 8. Next, let’s continue to consider the generalized possibility of the repeated reachability event  ⋄ B. Given the fuzzy state set B = (L (s, P)) sS, the maximum possibility and the minimum possibility of all states in M satisfying event  ⋄ B are respectively:

GPomax(B)=Pmax+diag(Pmax+(t,t))tSB=(0.80.80.30.50.90.90.30.81)+(0.80000.90001)(0.850.40.1)=(0.80.40.1),
GPomin(B)=Pmin+diag(Pmin+(t,t))tSB=(0.20.30.10.20.50.40.10.30.7)+(0.30000.500007)(0.850.40.1)=(0.30.40.1).

When the state is s0, ∥GPo ( ⋄ B) ∥ max (s0) =0.8, ∥GPo ( ⋄ B) ∥ min (s0) =0.3, indicates that the patients in this state are more likely to relapse. When the state is s2, ∥GPo ( ⋄ B) ∥ max (s2) =0.1, ∥GPo ( ⋄ B) ∥ min (s0) =0.1, indicates that once the patient’s disease in this state is cured, and it basically will not recur.

4.5Persistent reachability possibility

Let us consider persistent reachability properties events of the form ⋄  B. Let π = s0 α0s1 α1 … ∈ PathsAdv (s0) and B : S ⟶ [0, 1] be a fuzzy set of states in GPDPs. Thus, B(π)=i0jiB(sj) . Let us calculate GPoAdv (s ⊨ ⋄  B) under the maximum possibility scheduler and the minimum possibility scheduler.

Theorem 7. Let M = (S, Act, P, I, AP, L) be a GPDP, B : S ⟶ [0, 1] is the fuzzy set of states, max corresponds to the maximum possibility scheduler and min corresponds to the minimum possibility scheduler.

(14)
GPomax(B)=Pmax*rDBPmax,
(15)
GPomin(B)=Pmin*rDBPmin.

Proof. See the Appendix. ■

Example 9. Next, continue to consider the generalized possibility of persistent reachability event ⋄  B. Given the fuzzy state set B = (L (s, E)) sS, the maximum possibility and the minimum possibility of all states in M satisfied event ⋄  B are respectively,

DBPmax=(0.20000.70000.9)(0.80.80.30.50.90.90.30.81)=(0.20.20.20.50.70.70.30.80.9),
GPomax (⋄  B) =Pmax*rDBPmax

=(0.80.80.30.50.90.90.30.81)*(0.20.70.9)=(0.20.70.9),

DBPmin=(0.20000.70000.9)(0.20.30.10.20.50.40.10.30.7)=(0.20.20.10.20.50.40.10.30.7),
GPomin(B)=Pmin*rDBPmax

=(0.20.30.10.20.50.40.10.30.7)*(0.20.50.7)=(0.20.50.7), where DB={L(s,E)s=t0otherwise .

The possibility that the patient’s health continue to be “good” after the expert adopts the combination of three schemes. When the state is s0, ∥GPo (⋄  B) ∥ max (s0) =0.2, ∥GPo (⋄  B) ∥ min (s0) =0.2, which indicates that the patient is in a bad state of health and is already in a state of illness. For state s2, ∥GPo (⋄  B) ∥ max (s2) =0.9, ∥GPo (⋄  B) ∥ min (s2) =0.7, which indicates that the patient’s state health has been in a “good” state.

5Reachability in the safety properties

The analysis of safety properties and the techniques for checking safety properties are closely related to reachability. Safety properties are often characterized as “nothing bad should happen”. In classic model checking, safety properties are defined as linear property which does not include a bad prefix, and analyses the eventual reachability and repeated reachability. Since it is difficult to define the notion of a bad prefix in fuzzy logic or possibility logic, Li [21] uses the good prefixes to define the fuzzy safety property, and computes the always reachability possibility and persistent reachability possibility. In the following section, we use the good prefixes to analyze the possibilistic of regular safety property.

For a GPDP M = (S, Act, P, I, AP, L), we assume the alphabet Σ = lAP for some finite subset l ⊆ [0, 1] in the following.

Definition 11.c4:def1 For the safety property Psafe, we define a possibilistic language Gref (Psafe) : Σ* ⟶ [0, 1] as Gpref (Psafe) (θ) = ⋁  {Psafe (θσ) ∣ σ ∈ Σ ω}. For all θ ∈ Σ*, Gpref (Psafe) is called the good prefixes of Psafe.

Psafe is called possibilistic regular safety property if Psafe (σ) = ⋀  {Gpref (Psafe) (θ) ∣ θ ∈ Gpref (σ) }, for all σ ∈ Σ ω, Pref (σ) =  {θ ∈ Σ* ∣ σ = θσ′, σ′ ∈ Σ ω} is called the set of prefixes of σ.

We call Psafe a generalized possibilistic regular safety property if Psafe is a generalized possibilistic safety property and Gpref (Psafe) is a fuzzy regular language over Σ.

Definition 12. [22] A fuzzy finite automata is a five tuple A = (Q, Σ, δ, Q0, F), where

  • (1) Q is a finite nonempty set of states;

  • (2) Σ a finite nonempty set of input symbols;

  • (3) δ : Q × Σ × Q ⟶ [0, 1] is a fuzzy transition function. For any p, q ∈ Q, a ∈ Σ, δ (p, a, q) denotes the possibility of state p reaching state Q under the action of input letter a;

  • (4) Q0 : Q ⟶ [0, 1] represent the fuzzy initial state of fuzzy automata A,for q ∈ Q,Q0 (q) denotes q is the possibility of the initial state;

  • (5) F : Q ⟶ [0, 1] represent the fuzzy acceptance state of fuzzy automata A,F (q) denotes q is the possibility of the acceptance state;

For a GPDP M = (S, Act, P, I, AP, L) and a fuzzy finite automata A = (Q, Σ, δ, Q0, F), the product of M and A is defined as follows.

Definition 13. Let M = (S, P, Act, I, AP, L) be a GPDP, A = (Q, Σ, δ, Q0, F) be a fuzzy finite automata, M ⊗ A = (S × Q, P′, Act′, I′, AP′, L′), where P′ (< s, q > , α, < s′, q′ >) = P (s, α, s′) ∧ δ (q, L (s′) , q′); Act′ (< s, q >) = Act (s); I′ (< s, q >) = I (s) ∧ ⋁ q0QQ0 (q0) ∧ δ (q0, L (s) , q); for all <s, q > ∈ S × Q, AP′ = S × Q; L′ (< s, q >) = < s, q >.

For each path π = s0 α0s1 α1s2 α2 … ∈ PathsAdv (M) in M, TraceAdv (π) = L (s0L (s1L (s2)…, the fuzzy automata A has a unique run q0q1q2… and in M ⊗ A there exists π+ =< s0, q1 > α0 < s1, q2 > α1 … corresponding to them.Similarly, the path in M ⊗ A with state <s0, δ (q0, L (s))> also corresponds to the path in M and the run in A.

Theorem 8.Let M = (S, Act, P, I, AP, L) be a GPDP, Psafe is the generalized possibilistic regular safety property accepted by a deterministic fuzzy finite automata A = (Q, Σ, δ, Q0, F), and Adv is the set of all strategies begin state s. We have:

(16)
GPoAdvM(sPsafe)=GPoAdvMA(<s,qs>B),
where qs = δ (q0, L (s)), B = S× F = ∑sS,qQF (q)/< s, q >, which means that B (< s, q >) = F (q), for all <s, q > ∈ S × Q.

Proof. See the Appendix. ■

Let GPoAdv ( B) = GPoAdv (s ⊨  BsS, then

GPoAdv ( B) can be solved by the maximum fixed point, where the maximum fixed point of the operator is fB (Z) = B ∧ PAdv ∘ DZ ∘ rAdv. For the scheduler Adv taking the maximum possibility scheduler max and the minimum possibility scheduler min, the maximum possibility and the minimum possibility that the state s satisfies the generalized possibilistic regular safety property P are respectively:

(17)
GPomax(B)=(GPomax(sB))sS=υZ.fBmax(Z),
(18)
GPomin(B)=(GPomin(sB))sS=υZ.fBmin(Z),
where
fBmax(Z)=BPmaxDZrmax,fBmin(Z)=BPminDZrmin.

Example 10. We use Example 1 as a sample. To discuss the generalized possibility measure of generalized possibilistic regular safety property in the alphabet Σ = lAP, in which

l={0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1}.
Generalized possibilistic regular safety property Psafe = {A0, A1, … ∈ Σ ω ∣ ∀ i ≥ 0, Ai (E) >0} a describes the patient’s health states is good (E). Because the safety property Psafe has a good prefix Gref (Psafe) = {A0, A1, …, An ∈ Σ* ∣ n ≥ 0, ∀ i ≥ 0, Ai (E) >0}, and for any infinite words δ ∈ Σ ω, if w ∈ Pref (δ), w ∈ Gref (Psafe), then δ ∈ Psafe. Gref (Psafe) can be accepted by the definite finite fuzzy automata shown in the Fig. 4. Therefore, Psafe is generalized possibilistic regular safety property, where the letter E in the automata A denotes the atomic proposition for any A ∈ Σ such that A (E) >0 is an atomic proposition.

Fig. 4

Finite automata A corresponding to Gref (Psafe).

Finite automata A corresponding to Gref (Psafe).

Given a GPDP M and generalized possibilistic regular safety property Psafe = {A0, A1, … ∈ Σ ω ∣ ∀ i ≥ 0, Ai (E) >0}, Adv is the scheduler defined in M. The solution processes of GPomax (s0 ⊨ Psafe) and GPomin (s0 ⊨ Psafe) corresponding to Adv is the maximum scheduler and the minimum scheduler are given.

The product M ⊗ A of a GPDP M and Psafe good prefix Gref (Psafe) corresponding to a finite automata A is shown in Fig. 5. According to formulas 17 and 18, there B = S × {q1} = {1, 1, 1} is

Fig. 5

Product GPDPs M ⊗ A.

Product GPDPs M ⊗ A.

f(z1)=B(Pmax*Dzrmax)=(111)((0.80.80.30.50.90.90.30.81)*(100010001)(0.80.91))=(0.80.91),
f(z2)=B(Pmax*Dz1rmax)=(111)((0.80.80.30.50.90.90.30.81)*(0.80000.90001)(0.80.91))=(0.80.91),

fz1 = fz2, Therefore

GPomax(B)=(0.80.91).

Calculated by the same method

GPomin(B)=(0.30.50.7).

GPomax (s0 ⊨ Psafe) =0.8,GPomin (s0 ⊨ Psafe) =0.3 are obtained, which shows that the maximum possibility and minimum possibility of generalized possibilistic regular safety property Psafe in GPDPs M are 0.8 and 0.3.

6Conclusion

In this paper, firstly, we propose GPDPs as the models of nondeterministic and concurrent fuzzy systems. Then, we give fuzzy matrices calculation methods of the maximal possibilities and the minimal possibilities of reachabilities. Finally, we propose a model checking approach to convert the verification of safety property into the analysis of reachabilities. We have given the method of optimization algorithm for reachability. In the future, we will investigate the optimization algorithm for reachability to solve the verification of other properties, such as liveness and fairness in fuzzy systems.

Acknowledgments

The work is partially supported by National Natural Science Foundation of China(Grant No. 61962001), Graduate Innovation Project of North Minzu University(Grant No. YCX22184 and No. YCX22196).

Appendix

Proof. [Proof of Theorem 1]

Since S and Act is finite, given the scheduler Adv, the possibility transition matrix PAdv is also finite. Obviously, the fuzzy meet operator ∧ does not generate new element,so the set {i0PAdv(si,αi,si+1)siS,αiAct} is finite.Since S and Act is finite, for any states s ∈ S under consideration of the scheduler Adv, there exits t ∈ S and i < j so that si = sj = t.

In this case, PAdv (s, α0, s1)∧ PAdv (s1, α1, s2) ∧ ⋯

PAdv (s, α0, s1) ∧ ⋯ ∧ PAdv (si-1, αi-1, t) ∧ (PAdv (t, αi, si+1)

∧⋯ PAdv (sj-1, αj, t)) ∧ ⋯

PAdv (s, α0, s1) ∧ ⋯ ∧ PAdv (si-1, αi-1, t) ∧ (PAdv (t, αi, si+1)

∧ ⋯ ∧ PAdv (sj-1, αj, t))

PAdv+(s,Act(s),t)PAdv+(t,Act(t),t) .

Hence, rAdv(s)PAdv+(s,Act(s),t)PAdv+(t,Act(t),t) .

Conversely, for any t ∈ S and schedulers Adv, by the definition of PAdv+ , there exists ss1 ⋯ si = t ∈ S and si+1 ⋯ sj such that

PAdv+(s,t)=PAdv(s,s1)PAdv(si-1,t)

PAdv+(t,t)=PAdv(t,s1)PAdv(sj,t) .

Let πAdv = s0 α0s1 α1 ⋯ si-1 αi-1i (si+1 ⋯ sj  αjt ω,then PAdv+(s,t)PAdv+(t,t)=PAdv(s,s1)PAdv(s1,s2).

Hence, PAdv+(s,t)PAdv+(t,t)rAdv(s) , {PAdv+(s,t)PAdv+(t,t)tS}rAdv(s).

Therefore, rAdv(s)={PAdv+(s,t)PAdv+(t,t)tS} .

Proof. [Proof of Theorem 2]

From Cyl(πˆ)={πPathsAdv(M)πˆPref(π)} , we have Cyl (s0 α0s1 α1 ⋯ αn-1sn) = ⋃ {πAdv ∈ PathsAdv (M) ∣ s0 α0 ⋯ αn-1sn ∈ Pref (πAdv)} .

Therefore, GPo  (Cyl (s0 α0 ⋯ αn-1sn) )

= ⋁ {GPo (πAdv) ∣ s0 α0 ⋯ αn-1sn ∈ Pref (πAdv)}

={I(s0)i0PAdv(si,αi,si+1)s0α0s1αn-1snPref(πAdv)}={I(s0)i=0n-1PAdv(si,αi,si+1)}{jnPAdv(sj,αj,sj+1)}

=I(s0)i=0n-1PAdv(si,αi,si+1)}rAdv(sn).

Proof. [Proof of Theorem 3]

GPomax (s ⊨ ⋄ B)

= ⋁ πPathsmax(s)  (GPomax (π) ∧ ∥ ⋄ B ∥ (π) )

=(i=0αiActP(si,αi,si+1)j=0B(sj))

=i=0(Pmax(s,s1)Pmax(si-1,si)B(si))j=iPmax(sj,sj+1)

=i=0(Pmax(s,s1)Pmax(si-1,si)B(si)rmax(si))

=i=0(PmaxiDBrmax)(s)

=((i=0Pmaxi)DBrmax)(s)

=(Pmax*DBrmax)(s)

where Pmaxis the maximum possibility transition matrix under the max scheduler, Kleene closure of Pmax is Pmax* , DB is the diagonal matrix diag (B (s)) sS, rmax=Pmax+D , Pmax+=PmaxPmax* , D=(Pmax+(t,t))tS .

GPomin (s ⊨ ⋄ B)

= ⋁ πPathsmin(s)  (GPomin (π) ∧ ∥ ⋄ B ∥ (π) )

=(i=0αiActP(si,αi,si+1)j=0B(sj))

=i=0(Pmin(s,s1)Pmin(si-1,si)B(si))j=iPmin(sj,sj+1)

=i=0(Pmin(s,s1)Pmin(si-1,si)B(si)rmin(si))

=i=0(PminiDBrmin)(s)

=((i=0Pmini)DBrmin)(s)

=(Pmin*DBrmin)(s),

where Pmin is the minimum possibility transition matrix under the min scheduler.

Proof. [Proof of Theorem 5]

GPomax (s ⊨ C ⊔ nB) = ⋁ πPathsmax(s)  (GPomax (π) ∧ ∥ C ⊔ nB ∥ (π) )

= ⋁  (⋁  α0Act(s0)P (s0, α0, s1) ∧ ⋁  α1Act(s1)P (s1, α1, s2)

(j0nB(sj)ijC(si)))

=(Pmax(s0,s1)Pmax(s1,s2)(j0nB(sj)ijC(si)))

=(B(s0)rmax(s0))(j0nC(sj)ijPmax(si-1,si)C(si)Pmax(sj-1,sj)B(sj)rmax(sj))

=(i=0n(DCPmax)iDBrmax)(s).

GPomin (s ⊨ C ⊔ nB)

= ⋁ πPathsmin(s)  (GPomin (π) ∧ ∥ C ⊔ nB ∥ (π) )

=(α0Act(s0)P(s0,α0,s1)α1Act(s1)P(s1,α1,s2)(j0nB(sj)ijC(si)))

=(Pmin(s0,s1)Pmin(s1,s2)(j0nB(sj)ijC(si)))

=(B(s0)rmin(s0))(j0nC(sj)ijPmin(si-1,si)C(si)Pmin(sj-1,sj)B(sj)rmin(sj))

=(i=0n(DCPmin)iDBrmin)(s).

Proof. [Proof of Theorem 6]

GPoAdv (s ⊨  ⋄ B) = ⋁ πPathsAdv(s)  (GPo (π) ∧ ∥  ⋄ B ∥ (π) ). Let π = s0 α0s1 α1 … ∈ PathsAdv (s0), inf (π) denotes the set of states that occur infinitely many times on the path π, then ∥  ⋄ B ∥ (π) ≤ ⋁ tinf(π)B (t). Furthermore, for any t ∈ inf (π), π ⊨  ⋄ t, we can get GPoAdv (π) ≤ GPoAdv  ({π ∈ Pathsfin (s) ∣ π ⊨  ⋄ t} ) , thus GPoAdv (π) ∧ ∥  ⋄ B ∥ (π) ≤ ⋁ tinf(π) > (B (t) ∧ GPoAdv (s ⊨  ⋄ t) ) ≤ ⋁ tS  (B (t) ∧ GPoAdv (s ⊨  ⋄ t) ) . Therefore, GPoAdv (s ⊨  ⋄ B) ≤ ⋁ tS  (B (t) ∧ GPoAdv (s ⊨  ⋄ t) ) .

Conversely, for any state t ∈ S, and any path π ∈ PathsAdv (s) that satisfies the event  ⋄ t,we have B (s) ≤ ∥  ⋄ B ∥ (π). It follows that ⋁tS  (B (t) ∧ GPoAdv (s ⊨  ⋄ t) ) ≤ ⋁ tS  (B (t) ∧ GPoAdv (s ⊨  ⋄ t) ), therefore, ⋁tS  (B (t) ∧ GPoAdv (s ⊨  ⋄ t) ) ≤ GPoAdv (s ⊨  ⋄ B).

In conclusion, GPoAdv (s⊨  ⋄ B) = ⋁ tSB (t) ∧

GPoAdv (s ⊨  ⋄ t) .

We can get the method to calculate GPoAdv (s ⊨  ⋄ t) from these results, that is GPoAdv(st)=PAdv+(s,t)PAdv+(t,t) . Then we obtain GPoAdv(sB)=tSB(t)PAdv+(s,t)PAdv+(t,t).

When the Adv is the maximum possibility scheduler and the minimum possibility scheduler, we can get Theorem 6.

Proof. [Proof of Theorem 7]

GPomax (s ⊨ ⋄  B)

= ⋁ πPathsmax(s)  (GPomax (π) ∧ ∥ ⋄  B ∥ (π) )

=(GPomax(π)i=0j=iB(sj))

=i=0(α0ActP(s,α0,s1)α1ActP(s1,α1,s2)αi-1ActP(si-1,αi-1,si)B(si)αiActP(si,αi,si+1)

B (si+1) ∧ ⋁  αi+1ActP (si+1, αi+1, si+2) ∧ B (si+2) ∧ … )

=i=0(Pmax(s,s1)Pmax(s1,s2)Pmax(si-1si)B(si)Pmax(si,si+1)B(si+1)Pmax(si+1,si+2)

B (si+2) ∧ … )

=i=0(Pmax(s,s1)Pmax(s1,s2)Pmax(si-1si)

∧ (Dmax ∘ Pmax) (si, si+1) ∧ (Dmax ∘ Pmax) (si+1, si+2)

∧ … )

=i=0(Pmax(s,s1)Pmax(s1,s2)Pmax(si-1si)

rDBPmax (si) )

=i=0(siS(Pmax(s,s1)Pmax(s1,s2))

Pmax (si-1si) ∧ rDBPmax (si) )

=i=0(PmaxirDBPmax(s))=(i=0Pmaxi)rDBPmax(s)

=Pmax*rDBPmax(s),

GPomin (s ⊨ ⋄  B)

= ⋁ πPathsmin(s)  (GPomin (π) ∧ ∥ ⋄  B ∥ (π) )

=(GPomin(π)i=0j=iB(sj))

=i=0(α0ActP(s,α0,s1)α1ActP(s1,α1,s2)αi-1ActP(si-1,αi-1,si)B(si)αiActP(si,αi,si+1)B(si+1)αi+1ActP(si+1,αi+1,si+2)B(si+2))

=i=0(Pmin(s,s1)Pmin(s1,s2)Pmin(si-1si)B(si)Pmin(si,si+1)B(si+1)Pmin(si+1,si+2)B(si+2))

=i=0(Pmin(s,s1)Pmin(s1,s2)Pmin(si-1si)(DminPmin)(si,si+1)(DminPmin)(si+1,si+2))

=i=0(Pmin(s,s1)Pmin(s1,s2)Pmin(si-1si)rDBPmin(si))=i=0(siS(Pmin(s,s1)Pmin(s1,s2))Pmin(si-1si)rDBPmin(si))

=i=0(PminirDBPmin(s))=(i=0Pmini)rDBPmin(s)=Pmin*rDBPmin(s).

Proof. [Proof of Theorem 8]

GPoAdvM(sPsafe)

= ⋁ πPathsAdv(s)  (GPoAdv (π) ∧ Psafe (TraceAdv (π)) )

= ⋁  (GPoAdv (π) ∧ ⋀ j≥0  {L (A) (θ) ∣ θ ∈ Pref (TraceAdv (π)) } )

= ⋁  (GPoAdv (π) ∧ ⋀ j≥0  {F (qj) ∣ qjδ* (q0, L (sL (s1) … L (sj)) } ) ⋁  (GPoAdv (π) ∧ ⋀ i≥0F (qi) ) .

For π = 0s1 α1 … ∈ PathsAdv (s), we define the run sequence q0q1q2… of the deterministic fuzzy finite automata by qi+1 = δ (qi, L (si)). Vice versa, for the same run sequence q0q1q2… of the deterministic fuzzy finite automata, we have

GPoAdvMA(<a,qs>B)

= ⋁ π+PathsAdv(<s,qs>)  (GPoAdv (π+) ∧ ⋀ i≥0B (π+ [i]) )

= ⋁ πPathsAdv(s)  (GPoAdv (π) ∧ ⋀ i≥0F (qi) ) ,

References

[1] 

Mahmood I. A Verification Framework for ComponentBased Modeling and Simulation Putting the pieces together[J], arXiv preprint arXiv:2301.03088, 2023.

[2] 

Pan H. , Li Y. , Cao Y. et al. Model checking fuzzy computation tree logic[J]. Fuzzy sets and Systems 262: ((2015) ), 60–77.

[3] 

Laroussinie F. et al. Principles of Model Checking[J]. The Computer Journal 53: (5) ((2010) ), 615–616.

[4] 

Hansen H.A. , Schneider G. and Steffen M. , Reachability analysis of complex planar hybrid systems[J]. Science of Computer Programming 78: (12) ((2013) ), 2511–2536.

[5] 

Herbert S.L. , Bansal S. , Ghosh S. et al. Reachability-based safety guarantees using efficient initializations[C], 2019 IEEE 58th Conference on Decision and Control (CDC), IEEE, 2019, 4810–4816.

[6] 

Wang J. , Deng Y. and Xu G. , Reachability analysis of real-time systems using time Petri nets. IEEE Trans Syst Man Cybern B Cybern 30: (5) ((2000) ), 725–736. doi: 10.1109/3477.875448. PMID: 18252405

[7] 

Bouyer P. , Markey N. , Randour M. et al. Reachability in networks of register protocols under stochastic schedulers[J], arXiv preprint arXiv:1602.05928, 2016.

[8] 

Xie D. , Xiong W. , Bu L. et al. Deriving unbounded reachability proof of linear hybrid automata during bounded checking procedure[J]. IEEE Transactions on Computers 66: (3) ((2016) ), 416–430.

[9] 

Ganguli S. , Iyer C.V.K. , Pandey V. Reachability Embeddings: Scalable Self-Supervised Representation Learning from MobilityTrajectories for Multimodal Geospatial Computer Vision[C], 2022 23rd IEEE International Conference on Mobile Data Management (MDM), IEEE, 2022, 44–53.

[10] 

Song Y.H. , Shahid S. and Chung E.S. , Differences in multi-model ensembles of CMIP5 and CMIP6 projections for future droughts in South Korea[J]. International Journal of Climatology 42: (5) ((2022) ), 2688–2716.

[11] 

Li Y. , Droste M. and Lei L. , Model checking of linear-time properties in multi-valued systems[J]. Information Sciences 377: ((2017) ), 51–74.

[12] 

Ma Z. , Li Z. , Li W. et al. Model Checking Fuzzy Computation Tree Logic Based on Fuzzy Decision Processes with Cost[J]. Entropy 24: (9) ((2022) ), 1183.

[13] 

Lee N. , Kim Y. , Kim M. et al. Directed Model Checking for Fast Abstract Reachability Analysis[J]. IEEE Access 9: ((2021) ), 158738–158750.

[14] 

Wang Y. , Yuan Y. , Zhang W. et al. Reachability-Driven Influence Maximization in Time-dependent Road-social Networks[C], 2022 IEEE 38th International Conference on Data Engineering (ICDE), IEEE, 2022, 367–379.

[15] 

Wu D. and Koutsoukos X. , Reachability analysis of uncertain systems using bounded-parameter Markov decision processes[J]. Artificial Intelligence 172: (8-9) ((2008) ), 945–954.

[16] 

Wu D. and Koutsoukos X. , Reachability analysis of uncertain systems using bounded-parameter Markov decision processes[J]. Artificial Intelligence 172: (8-9) ((2008) ), 945–954.

[17] 

Li Y. and Li L. , Model checking of linear-time properties based on possibility measure[J]. IEEE Transactions on Fuzzy Systems 21: (5) ((2012) ), 842–854.

[18] 

Li Y. and Wei J. , Possibilistic fuzzy linear temporal logic and itsmodel checking[J]. IEEE Transactions on Fuzzy Systems 29: (7) ((2020) ), 1899–1913.

[19] 

Jiang J. , Zhang P. and Ma Z. , The μ-calculus model-checkingalgorithm for generalized possibilistic decision process[J]. Applied Sciences 10: (7) ((2020) ), 2594.

[20] 

Am A. , Mp B. , Sr C. et al. Model-checking graded computation-tree logic with finite path semantics[J]. Theoretical Computer Science 806: ((2020) ), 577–586.

[21] 

Li Y. , Quantitative model checking of linear-time properties based on generalized possibility measures[J]. Fuzzy sets and Systems 320: ((2017) ), 17–39.

[22] 

Li Y. and Ma Z. , Quantitative computation tree logic model checking based on generalized possibility measures[J]. IEEE Transactions on Fuzzy Systems 23: (6) ((2015) ), 2034–2047.

[23] 

Ying S. and Ying M. , Reachability analysis of quantum Markov decision processes[J], . Information and Computation 263: ((2018) ), 31–51.

[24] 

Huang S. and Cleaveland R. , A tableau construction for finite linear-time temporal logic[J]. Journal of Logical and Algebraic Methods in Programming 125: ((2022) ), 100743.

[25] 

Quatmann T. , Junges S. , Katoen J.P. Markov automata with multiple objectives[C], International Conference on Computer Aided Verification, Springer, Cham, 2017, 140–159.

[26] 

Shojafar M. , Javanmardi S. , Abolfazli S. et al. FUGE: A joint meta-heuristic approach to cloud job scheduling algorithm using fuzzy theory and a genetic method[J]. Cluster Computing 18: (2) ((2015) ), 829–844.

[27] 

Wei X. and Li Y. , Fuzzy alternating automata over distributive lattices[J]. Information Sciences 425: ((2018) ), 34–47.

[28] 

Madhumangal Pal Fuzzy Matrices with Fuzzy Rows and Columns, 1 Jan. 2016, 561–573.

[29] 

Wang Q. , Li R. , Ying M. Equivalence Checking of Sequential Quantum Circuits[J], IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, 2021. doi: 10.1109/TCAD.2021.3117506.

[30] 

Garmendia L. , Salvador , Gonzĺćlez del R. , Campo , Lĺöpez et al. An algorithm to compute the transitive closure, a transitive approximation and a transitive opening of a fuzzy proximity[J]. Mathware and Soft Computing 16: (2) ((2009) ), 175–191.

[31] 

Hartmanns A. , Junges S. , Katoen J.P. et al. Multi-cost bounded reachability in MDP[C], International Conference on Tools and Algorithms for the Construction and Analysis of Systems, Springer, Cham, 2018, 320–339.

[32] 

Chondrogiannis T. , Nascimento M.A. , Bouros P. Relative reachability analysis as a tool for urban mobility planning: Position paper[C], Proceedings of the 12th ACM SIGSPATIAL International Workshop on Computational Transportation Science, 2019, 1–4.