Программные продукты и системы

Теория нечетких множеств широко применяется для решения разнообразных задач. Однако случайные изменения функций принадлежности до сих пор не рассмотрены, хотя некоторые из них имеют стохастический характер. Их можно определить как функции принадлежности, значения которых являются случайными величинами.

Прогнозирование лавинной опасности с помощью этой теории описано, например, в работах [1, 2]. Вычисляются функции принадлежности заданной ситуации к различным ее уровням. В то же время для конкретной даты в известном районе величины этих функций меняются случайным образом для различных лет, что обусловлено стохастическим изменением метеоусловий и сейсмической обстановки в разные годы. Кроме того, реально доступной является только малая выборка, так как время наблюдения за лавинными очагами, как правило, ограничено.

Таким образом, информация о случайной изменчивости указанных функций очень важна для обеспечения безопасности работ в горах, как и для оценки других опасных или нежелательных ситуаций, связанных с природными, социально-эко­номическими или производственными процессами. Поэтому разработка соответствующего ПО представляет значительный научно-практический интерес.

При проектировании, в частности, различных сооружений и дорог в горной местности требуется оценить возможность появления лавинной опасности различного уровня. За известный период наблюдений лавины могут не возникнуть, но функции принадлежности ситуаций к лавиноопасным вычисляемы. Тогда, зная математическое ожидание и среднее квадратичное отклонение функций принадлежности, можно оценить риск с помощью правила трех сигм [3] или правила шести сигм [4, 5]. Однако все это требует применения стохастических функций принадлежности.

 Кроме того, если плотность вероятности функции принадлежности ситуации к какому-либо опасному или нежелательному состоянию известна, можно выполнять и более сложные вероятностные расчеты. Например, на каждой реализации метода Монте-Карло вычислить отношение нагрузки от снеговоздушной волны к допускаемой (тоже стохастической) нагрузке на сооружение, затем построить соответствующую плотность вероятности и определить вероятность разрушения.

Классические методы восстановления зависимостей эффективны только при использовании выборок достаточно большого объема [6]. Если выборка невелика, лучшие результаты получаются при применении метода структурной минимизации риска [6, 7], идея которого состоит в следующем. Если на допустимом множестве решений задать структуру, то наряду с минимизацией эмпирического риска внутри ее элементов появляется дополнительная возможность оптимизации качества оценки по элементам структуры, что позволяет найти решение, дающее более глубокий гарантируемый минимум среднему риску, чем решение, доставляющее минимум эмпирическому риску на всем допустимом множестве [6].

Тем не менее при достаточно малой выборке метод структурной минимизации риска обычно выбирает наиболее простую зависимость. Действительно, верхняя граница среднего риска при восстановлении плотности вероятности оценивается по формуле

, (1)

где N – объем выборки; k – степень аппроксимирующего полинома; [y] – вектор значений эмпирической функции распределения; [F(α)] – вектор значений аппроксимирующего полинома в экспериментальных точках; αi – коэффициенты аппроксимирующего полинома; R – обратная ковариационная матрица вектора [y]; 1–η – вероятность, с которой справедлива оценка (1).

При k=1, N=10, η=0,95 знаменатель (3) равен 0,274, а при k=2, N=10, η=0,95 он составляет 0,184, то есть его уменьшение происходит очень быстро.

Кроме того, на развитие опасных социально-экономических, природных и техногенных процессов обычно влияет большое количество случайных факторов. Так, согласно [1], возникновение снежных лавин зависит от суммы осадков, их интенсивности, температуры воздуха и т.д. Поэтому здесь вполне может быть использована центральная предельная теорема, согласно которой достаточно большая сумма сравнительно малых случайных величин распределена приближенно по нормальному закону [8].

Случайная величина, распределенная таким образом, изменяется от –∞ до ∞. Поскольку область изменения функций принадлежности обычно представляет собой интервал [0; 1], плотность вероятности отыскивается в виде

,                                           (2)

где а – математическое ожидание случайной величины; σ – среднее квадратичное отклонение случайной величины; b – коэффициент, который всегда больше единицы.

Математическое ожидание и среднее квадратичное отклонение вычисляются (см. [3]) по формулам:

,                                                                (3)

.                                                 (4)

Величина коэффициента b должна быть такой, чтобы выполнялось условие

Тогда .                            (5)

Плотность вероятности типа (2) иллюстрируется рисунком.

После заполнения таблицы Initial Data необходимо нажать кнопку Reconstructing. Далее в таблице Results of Analysis будут напечатаны результаты вычислений.

В качестве контрольного примера можно привести восстановление плотности вероятности по случайной выборке, приведенной в столбце 2 таблицы. Эта выборка получена с помощью датчика случайных чисел программного комплекса ANSYS. Задавались математическое ожидание 0,5 и среднее квадратичное отклонение 0,05. В результате расчета было получено: а=0,507, σ=0,041, b=1. Действительно, случайные числа, распределенные по этому закону, практически всегда будут находиться в интервале [0; 1], поэтому коэффициент b должен быть равен единице.

Элементы выборки для восстановления плотности вероятности

Другим примером является восстановление плотности вероятности случайной величины х по выборке из 6 элементов, приведенных в столбце 3 таблицы. В этом случае а=0,458, σ=0,195, b=1,056.

И еще один пример – восстановление плотности вероятности случайной величины х по выборке из 6 элементов, приведенных в столбце 4 таблицы. Плотность вероятности, восстановленная по ним, имеет следующие параметры: а=0,342, σ=0,255, b=1,042.

Таким образом, разработанный комплекс программ позволяет дополнять нечеткий анализ стохастическим, улучшая тем самым точность оценки различных сложных ситуаций, обусловленных опасными или нежелательными социально-экономическими, природными и техногенными процессами. Он требует около 2,2 Мб памяти и выполняет вычисления за незначительное время.

Литература

1.     Зимин М.И. Прогнозирование лавинной опасности. СПб: Гидрометеоиздат, 2000. 16 с.

2.     Зимин М.И. Прогнозирование опасных процессов на основе бионического подхода и его использование в системах автоматизации проектирования // Естественные и технические науки. 2011. № 3. C. 405–412.

3.     Пискунов Н.С. Дифференциальное и интегральное исчисление. М.: Наука, 1976. Т. 2. 576 с.

4.     Зорин А.А. Время шести сигм // Методы менеджмента качества. 2006. № 4. С. 32–36.

5.     Nanda V., Robinson J.A., Six Sigma Software Quality Improvement. Success Stories from Leaders in the High Tech Industry, NY, McGraw Hill, 2011, 611 p.

6.     Вапник В.Н., Глазкова Т.Г., Кощеев В.А., Михальский А.И., Червоненкис А.Я. Алгоритмы и программы восстановления зависимостей. М.: Наука, 1984. 816 с.

7.     Вапник В.Н. Восстановление зависимостей по эмпирическим данным. М.: Наука, 1979. 220 с.

8.     Колемаев В.А. Теория вероятностей в примерах и задачах. М.: Гос. ун-т управления, 2001. 87 с.

References

1.     Zimin M.I., Prognozirovanie lavinnoy opasnosti [Forecasting danger of avalanches], St. Petersburg, Gydrometeoizdat, 2000.

2.     Zimin M.I., Estestvennye i tekhnicheskie nauki [Natural and engineering sciences], 2011, no. 3, pp. 405–412.

3.     Piskunov N.S., Differentsialnoe i integralnoe ischislenie [Differential and integral calculus], Vol. 2, Мoscow, Nauka, 1976.

4.     Zorin A.A., Metody menegmenta kachestva [Quality management methods], 2006, no. 4, pp. 32–36.

5.     Nanda V., Robinson J.A., Six Sigma Software Quality Improvement (Success Stories from Leaders in the High Tech Industry), NY, McGraw Hill, 2011.

6.     Vapnik V.N., Glazkova T.G., Koshcheev V.А., Mikhalsky А.I., Chervonenkis А.Ya., Algoritmy i programmy vosstanovleniya zavisimostey [Dependencies recovery algorithms and programs], Мoscow, Nauka, 1984.

7.     Vapnik V.N., Vosstanovlenie zavisimostey po empiricheskim dannym [Dependencies recovery with empirical data], Мoscow, Nauka, 1979.

8.     Kolemaev V.A., Teoriya veroyatnostey v primerakh i zadachakh [The theory of probability in examples and exercises], Мoscow, State Univ. of Management Publ., 2001.