Fuzzification of convex bodies in Rn / Фазификация выпуклых тел в 𝑹𝒏

Оганян Виктор Кароевич доктор физико-математических наук профессор, зав. кафедрой теории вероятностей и математической статистики, Ереванский государственный университет 0025, Армения, г. Ереван, ул. Алекса Манукяна, 1 Ohanyan Viktor Karoevich Doctor of Physics and Mathematics Professor, the head of the department of Probability Theory and Mathematical Statistics, Yerevan State University 0025, Armeniya, g. Erevan, ul. Aleksa Manukyana, 1 victoohanyan@ysu.am Бардахчян Вардан Геворкович кандидат физико-математических наук ассистент, кафедра Актуарной Математики и финансовой Математики, Ереванский государственный университет 0025, Армения, г. Ереван, ул. Алекса Манукяна, 1 Bardakhchyan Vardan Gevorkovich PhD in Physics and Mathematics Assistant, Department of Actuarial Mathematics and Financial Mathematics, Yerevan State University 0025, Armeniya, g. Erevan, ul. Aleksa Manukyana, 1 vardanbardakchyan@gmail.com Симонян Арсен Рафикович кандидат физико-математических наук доцент, кафедра прикладной математики и информатики, Сочинский государственный университет 354008, Россия, Краснодарский край, г. Сочи, ул. Политехническая, 7 Simonyan Arsen Rafikovich PhD in Physics and Mathematics Associate Professor, Department of Applied Mathematics and Computer Science, Sochi State University 354008, Russia, Krasnodarskii krai, g. Sochi, ul. Politekhnicheskaya, 7 oppm@mail.ru

Улитина Елена Ивановна кандидат физико-математических наук доцент, кафедра прикладной математики и информатики, Сочинский государственный университет 354008, Россия, Краснодарский край, г. Сочи, ул. Политехническая, 7 Ulitina Elena Ivanovna PhD in Physics and Mathematics Associate Professor, Department of Applied Mathematics and Computer Science, Sochi State University 354008, Russia, Krasnodarskii krai, g. Sochi, ul. Politekhnicheskaya, 7 ulitina@rambler.ru

The problem of identifying of convex bodies is widely known problem in integral geometry and several tools have been introduced to deal with it (starting from support function determination, and integral representation ending with random line length distribution and more complex structures like V-lines) ^[1-4].

There is also a general way of identifying the convex body through its covariogram, determination of which requires other sort of measurements as well. One of the techniques suggests taking line segments’ length distributions (for all directions) and recover from that covariogram of the body (for several types of convex bodies exact formula of reversion are known ^[5-6]).

Methodology

Matheron’s theorem provides differential tool connecting covariogram and length distribution of a parallel line of given direction. However, in most cases, measurements will contain errors. This makes precise identification of body flawed. Instead of dealing with error bounds (in covariogram) and finding upper values for it, we deal with other form of error incorporation, fuzzy sets. By doing so, we try to describe all possible convex bodies which we could be recovered for given.

So we deal with the inverse problem in some sense. We define a body and ask what result we will get for length distribution for the given fuzzy body.

Through this procedure, we generalize Matheron’s theorem for this case. We extensively use fuzzy statistics and fuzzy random variables to extend convex bodies and length distribution functions to a fuzzy case. Use several properties of fuzzy numbers and fuzzy calculus techniques (mainly Aumann integration).

First, we generalize a fuzzy distribution function for our purposes. Then define fuzzy convex bodies as a natural generalization of the convex body to the fuzzy case.

Next, we derive some useful properties of fuzzy convex bodies and fuzzy covariogram. And derive Matheron’s theorem generalization for fuzzy case.

Results

1. General fuzzy distribution

Suppose we have fuzzy random variable, i.e. random variable which takes on fuzzy numbers as values.

More formally we have , with

(1)

α-levels. As fuzzy numbers these α-levels must obey inclusion principle, i.e. for all ω-s for .

Generally, fuzzy densities are given by precise definition (for example in Viertl’s book ^{[7, p. 55]}; ^[16]). However, no case of discrete distribution is discussed. To define fuzzy distribution, we should analyze the idea of taking values smaller than or equal to given x.

Suppose we have two possible outcomes i.e. , and we have their probabilities .

Suppose we take and . Now we consider some x.

It is obvious that if then . It is also obvious that if then .

Now let’s consider the case where x is in between two fuzzy number.

So if we have then

However, if x belongs to the interval of α-levels we get the following:

For simplicity let’s suppose that , so that they have no intersection (or strictly ).

And suppose . Here we either can have that outcome corresponding either took value less than or equal to x, or value greater than x. So in this case

So in general will not be fuzzy number for each x.

The general case for arbitrary Ω would be.

This is general formula, for fuzzy distributions.

When we are speaking about fuzzy densities we require some interval for values of cdf-s.

The general framework is as follows.

(3)

with finite integral of both and . And such that there is some true density function in the interval with highest belief.

And fuzzy distributions for continuous case will be defined

(4)

So this obviously become fuzzy function in sense that for each x it will be fuzzy number.

This structure is somewhat intuitive, while the first structure is purely constructive.

For our case we can unintendedly encounter non continuous probabilities. So we should use general formulation.

2. Fuzzy convex bodies

Now let’s skip to our main problem. First of all, let’s define a fuzzy case of convex body . We will do it in the following way.

(5)

where is some fuzzy number with support in and with α-cuts. Here D is understood as fuzzy number with support in , whose membership function is D for any level α. Next

(6)

where minus is understood as Hukuhara difference. So is some fuzzy number for which it is possible. must be satisfied.

Now the final part is

(7)

and is defined as a fuzzy concept having as its α-levels ^[8].

Proposition 1. One can note that each level is a convex set, from the fact that Minkowski sum of two convex bodies is again convex. Putting it another way, if , then , for .

Proof. For cases a=0, and a=1, the proposition is obvious. Next, we give two simple arguments.

and are convex bodies, and so is convex. Moreover

(8)

and

(9)

Summing up, we get

(10)

Proposition 2. For , .

Proof: We must show, that each element of is also contained in .

Let’s first note that and . From which we have

(11)

is a convex set that must be summed with to get D. is a convex set that must be summed with to get D. As , . So we have

(12)

From which obviously if , then . So

(13)

and as G is convex. So .

Remark 1. So defined satisfy all criteria of fuzzy number ^[9], except that of elements of are convex bodies, not numbers.

3. Fuzzy covariogram and its properties.

Now we can define covariogram for (for covariogram in classical case ^[10]; ^[17]). Though as a function of this fuzzy argument it has not developed approaches for further investigation, we try to show that it is a fuzzy number, and will try to recover its membership function.

Let be the (n-1)-dimensional unit sphere with center at the origin.

Definition1. We will name fuzzy covariogram of the fuzzy concept, , with following α-levels

(14)

where and , with being n-dimensional Lebesgue measure.

Remark2. So fuzzy covariogram has α-levels, and this are collection of numbers.

Proposition 3.For , we have .

Proof. From Proposition 2 we have , so if , than , and

(15)

Proposition 4. For each , is a fuzzy number.

Proof. By the inclusion property for (proposition 2), we have, if

(16)

Next for completeness we should proof the continuity of membership function, which can be done in the following way. First, we will show, that the membership function satisfies the intermediate value property. And from boundedness we can imply that it is continuous.

To do it take the two layers. and , with . Let’s look at the smallest values and . From the inclusion property we have. . Whenever the equality is hold the intermediate value property is immediate. Suppose now .

Let’s take and , with and . This is guaranteed by the properties of . Now look at the convex combination of this two: . For each there is .

Obviously, there is such that for each , and

(17)

This proves that intermediate value property is satisfied. And from bounded of membership function, and monotonicity of its left and right bounds of each α-levels we can say that it is continuous. The same can be done with maximal values of covariogram. And for the remaining part, the function is constant.

Remark 3. In the above proof we extensively used the property of convex bodies, for dealing with Lebesgue measure of Minkowski sum of convex bodies, and to interchange the sign of intersection and addition. In general this is forbidden, as the sum may be even non-measurable ^[11].

Next, if we define covariogram as a function of the magnitude of translation (i.e. ), then we have the following. Here the translation is a vector (t,u), where defines direction, and t determines magnitude.

Proposition 5. defined by α-levels

(18)

Is the fuzzy function of t, i.e. the function of , which take fuzzy values.

Remark 4. The later proposition is crucial, as it gives the opportunity to use Hukuhara differentiation and Aumann integral concepts.

Proposition 6. is not t-Hukuhara-differentiable a.s. in .

The idea behind is that is getting “thinner”. Or in other word is 0 starting from some t. Suppose, then is not definable in Hukuhara sense.

So we should use Aumann integration concept instead.

First, we should define distribution functions for our case. The difference from the crisp(classical) case is that the projection space is different for different elements of . So first we should define .

(19)

where .

We know that this is homeomorphic to ^[12]

And we assign a measure to each element of , .

Though we can combine them into a fuzzy concept, it is rather easier to leave it like this.

To use the definition of (1) for the fuzzy random variable, we should have the same universe of discourse, or in other words we should have the same Ω. However, Ω is different for different shapes. To bypass this problem, we first note that Ω has not to be the same for all α-levels. So we can deal with different Ω–s on different levels, but we should exactly fix Ω for any α-level. So we will consider the following distribution functions.

(20)

Where is length of intersection of a line parallel to u and intersecting or at some point ω.

The crucial theorem here is the following

Theorem 1.

Proof: To show this, we should literally show that any value of has corresponding with for which .

Not to embark our paper with much text, let’s just give the idea of proof. For any value of and , the maximal value is fixed by which corresponds to

(21)

Note that this set is always non-empty. This is guaranteed by the convexity of all sets considered. So we can use the convex hull of the union of sets. Moreover, this set contains only one element. The next bound is also fixed. So any value in-between can be obtained by considering a convex combination of these two bounds. No matter which convex body is considered, while in a given direction and for given magnitude it has the desired value of distribution function.

So, when we fix the Ω we can easily speak about .

And next theorem which comes along as generalization of already known result

Theorem 2. . (compare with the classical case in ^[3]). Here the integral used is Aumann integral.

Proof. This is immediate consequence of the previous theorem and Matheron theorem about covariogram. We have

(22)

where in the last line integration in Aumann sense is used, and on the previous line in Riemann-Lebesgue sense.

Discussion

Measurement error could be incorporated into model trough different error addition. An example can be error described random variable from given distribution (or distribution family). However, this approach can have several issues. First of all, this could be misleading, as the reconstruction of covariogram will require reversion. More simply suppose that you have some error field added to the given convex body. From that point, you should get rid of the error part to get the body itself. However, this will again lead to many possible alternatives of convex bodies, which you could get, as you don’t know the error exactly. So the covariogram will itself represent the random field. In that situation, other approaches like support function representation can be more fruitful [13 – 15].

Next problem is concerned with possible unboundedness of error term. If we have Gaussian field the possible infinite of (very big error can occur), which is not the case of measurement error. So measurement error when dealing with convex bodies, are more likely to be bounded ^[16].

For further we are going to try to get a generalization of other formulas, connecting direction dependent distribution of length, covariogram and distribution of length of line segments ^[12].

Conclusion

Fuzzy convex bodies are the natural generalization of convex bodies when we want to incorporate measurement error into our model and derive all possible convex bodies that could give rise of the given direction dependent length distribution. From this point of view, fuzzy convex bodies can be seen as a collection of convex bodies (though at each level it contains not only convex bodies).

We introduced the natural generalization of fuzzy distribution function based on the ordinary distribution function. Introduced the concept of a fuzzy convex body. Discussed some of its properties making it similar up to domain with fuzzy numbers.

We introduced fuzzy covariogram based on fuzzy convex bodies.

Next, we construct a generalized fuzzy length distribution based on fuzzy distributions introduced previously. And lastly, we give equation connecting it with fuzzy covariogram, thus generalizing Matheron’s theorem to the fuzzy case.

Acknowledgements

The present research was partially supported by the RA MES State Committee of ‎Science and Russian Federation Foundation of Innovation Support in the frame ‎of the joint research project 18RF-019 and 18-51-05010 Арм-а accordingly.‎