Интеллектуализация обработки информации при использовании дискретных методов исследования динамических систем

INTELLIGENT MULTI-AGENT SYSTEMS MODELING FOR AUDITING QUALITY EDUCATION TECHNICAL UNIVERSITY

A.P. Ragovsky

The report focuses on solving the problem of modeling intelligent multi-agent system to ensure audit quality education in technical colleges. The article discusses the architecture and principles of intellectual multi-agent system. Architecture of an intelligent multi-agent system is based on the modeling of the main component with the

agent-oriented approach. Model of intelligent agents is built on the principles of the organization of the extended network, where the ratio of communication plays a central role. Mechanism mediated communication agents, as well as the proposed model allows agents to share the basic process of searching for directly making and management processes. Since it entered into the knowledge base model of the application domain can be subject to both structural and qualitative changes, support setting up a system for the subject area and the types of tasks.

Key words: intelligent application systems, multiagent systems, knowledge engineering, quality education.

УДК 519.240

ИНТЕЛЛЕКТУАЛИЗАЦИЯ ОБРАБОТКИ ИНФОРМАЦИИ ПРИ ИСПОЛЬЗОВАНИИ ДИСКРЕТНЫХ МЕТОДОВ ИССЛЕДОВАНИЯ

ДИНАМИЧЕСКИХ СИСТЕМ

И.Н. Статников, Г.И. Фирсов

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

Ключевые слова: планирование вычислительных экспериментов, ЛП-сетки, имитационное моделирование.

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

вило, выведенное Полем Фитсом [1, 2]: время достижения цели обратно пропорционально ее размеру и дистанции до нее. Если объем исходной области поиска обозначить через Б, а объем области, содержащей предпочтительные решения, как S, то число вычислительных экспериментов может быть определено по формуле: N = а + Ь 1оя2[СО/+1], где а и Ь - некоторые константы. Перефразируя известное выражение, можно сказать, что как могущество человека прирастает могуществом создаваемой им техники, так и эффективность используемой техники определяется во многом интеллектуальными возможностями человека, т.е. мыслительными способностями человека освоить получаемую информацию (т-tellctus - разум; мыслительные способности). Это положение полностью подтверждается историей развития вычислительной техники. Как только появились ЭВМ (пусть еще и не-

совершенные в современном смысле), стало возможным использование методов Монте-Карло, которые мы и рассматриваем как представителей идеологии «слепого» перебора [3]. Привлекала как их независимость от вида широкого класса исследуемых функций, так и простота алгоритмической реализации на ЭВМ (хотя всегда есть проблемы качества используемых датчиков случайных или квазислучайных (псевдослучайных) чисел). Но возникало противоречие: колоссальные объемы числовой информации и возможности интерпретации исследователем таких объемов. Статистическая обработка числовой информации, добываемой при использовании дискретных методов, представляет собой, безусловно, пример интеллектуализации процесса анализа этой информации. Однако стандартный набор статистических характеристик (средние, среднеквадратичные отклонения и т.д.) все же оставлял без ответа очень важный вопрос: а что между двумя последовательными испытаниями функции? Поэтому постоянно присутствует потребность в дополнении существующих стандартных пакетов математической обработки методами, повышающими степень их интеллектуализации.

Полагаем, что одним из таких методов является ПЛП-поиск, основная идея которого заключается в синтезе двух подходов: дискретного случайного поиска и его особой организации в пространстве варьируемых параметров исследуемой функции [4]. Благодаря реализации в ПЛП-поиске этой идеи удается статистическим путем (и визуально; графически) оценить в среднем чувствительность исследуемой функции к изменениям параметров, в нее входящих со всеми вытекающими из этого обстоятельства возможностями при решении задач исследования (и проектирования) динамических систем.

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

или об окончании решения

векторов а , рассчитываемых по ЛПТ-сеткам [5], и которая оказывается возможной благодаря тому, что весь вычислительный эксперимент проводится сериями. В ПЛП-поиске на сегодняшний день можно варьировать одновременно значения до 51-го параметров (J = 51). Для рандомизации (случайного смешения уровней варьируемых параметров aijh) дискретного обзора G(a) могут быть

использованы многие существующие таблицы равномерно распределенных по вероятности целых чисел. В целях экономии памяти ЭВМ в ПЛП-поиске алгоритм рандомизации построен на использовании датчика псевдослучайных чисел q (0 < q <1) из [5]. Рандомизация состоит в том, что для каждой h-й серии экспериментов (h=1,..., H( i, j )), где H( i, j) - объем выборки из элементов ФуА для одного критерия, вычисляется свой вектор случайный номеров строк J = (j 1 h, j2 h Jph ) в таблице направляющих числителей (ТНЧ) по формуле:

j ph = [R х q] + 1,

а значения ау в h-й серии рассчитываются

с помощью линейного преобразования

а..,=а,+ хДа., где Да.=а„-а,, а *,а, -

iJh j ihJph J J J J J J

соответственно верхние и нижние границы области G(d); p = 1, ..., J; R - любое целое число (в ПЛП-поиске R = 51); J - фиксированный номер варьируемого параметра; i = 1,., M (J) - номер уровня J -го параметра в h-й серии; M (J ) - число уровней, на которое разбивается J -й параметр; в общем случае jph Ф J (в чем и состоит одна из целей рандомизации). Можно показать с помощью критерия Романовского [6], что числа jph, вырабатываемые по формуле (1), оказываются совокупностью равномерно распределенных по вероятности целых чисел. Обратим внимание, что M ( J ) и есть количество экспериментов, реализуемых в одной серии. И если M (J) = M = const и H (i, J) = H = const,

то в этом случае параметры N0, M и H связаны простым соотношением:

N0 = M х Я,

где N0 - общее число вычислительных экспериментов (ВЭ), при этом длина выборки из ф й в точности равна H . Но в общем

случае, когда M (j) = var, то и H (i,j ) = var, и тогда формула (2) для одного критерия примет такой вид:

м

N 0 = £ H (i, j). i=1

Для проведения однофакторного дисперсионного анализа [7] по всем параметрам для каждого критерия производится сортировка результатов вычислений, полученных при вычислении в точках матрицы планируемых экспериментов (МПЭ). В результате сортировки для одного критерия будет получено J матриц, состоящих из элементов ф , а для K критериев будет

получено J х K матриц, состоящих из элементов ф , где k - номер критерия. Этот

анализ позволяет принять (или отвергнуть) с требуемой вероятностью P = 1 -а, где а -заданный уровень значимости, следующую нулевую гипотезу: средние значения ф к

не существенно (случайно) отличаются от общего среднего значения k -го критерия Фм. Если принят положительный ответ (гипотеза принята), то допускается на следующем этапе решения задачи несущественно влияющий параметр а. не варьировать, а зафиксировать одно из его значений, например, а. = ау для такого i, где ф

имеет наилучшее значение в смысле искомого экстремума.

Иначе говоря, инструментальная основа ПЛП-поиска, реализованного в среде MAT-LAB [8, 9], состоит в следующем [10]. В пространстве варьируемых (исследуемых) параметров а}- (j=1,...,J) строится матрица планируемых экспериментов (МПЭ) размером (N0 *J), где N0 - число строк этой матрицы (число вычислительных экспериментов), а J - число столбцов. При этом весь численный эксперимент проводится сериями, в каждой из которых реализуется алгоритм ЛП-поиска [5]. При этом каждая серия вычислительных экспериментов рассчитывается по совокупности J строк из таблицы направ-

ляющих числителей {V jm}, где j=1,...,51,

m=1,...,20 [5]. Номера строк для каждой серии берутся из совокупности равномерно распределенных по вероятности целых чисел, вырабатываемой в ПЛП-поиске с помощью того же датчика псевдослучайных чисел qu (0 < qu < 1), что используется и при выработке случайных значений агу, где i=1,..., N0. При реализации ЛП-поиска в каждой серии каждый варьируемый параметр в диапазоне своего изменения разбивается на M сечений. В ПЛП-поиске предусмотрены варианты, когда M=2l, M=2l+1, M= 21 (l=2,3,4,...). В свою очередь также реализуются для каждого из трех перечисленных случаев варианты при M=const и M=var.

В сформированной указанным способом МПЭ производятся вычисления по каждому критерию в каждой строке матрицы. Далее определенным выше образом производится сортировка результатов вычислений, привязанная к МПЭ, и появляются возможности использования оценок математической статистики для определения статистического влияния каждого варьируемого параметра на каждый вычислявшийся критерий.

Отметим, что в ПЛП-поиске представляется возможным варьировать до 51 параметра и разбивать варьируемые параметры на число сечений до 1024, хотя реально употребляются от 4 до 256 сечений.

В результате применения ПЛП-поиска можно с требуемой вероятностью получить одновременно следующие результаты (или некоторые отдельно):

1) определить существенные (релевантные) параметры а (m < J + U) в смысле их m

влияния на значения каждого критерия Ф (а) ; иначе говоря, статистическим путем

К

оценить в среднем изменения производных

дФ^ /да^ при а^ е ; здесь

£ = 1, J + и, 3 - число варьируемых параметров, и - число параметров управления динамической системой; к = 1, К , где К - число предъявленных критериев качества к системе;

2) выделить области концентрации О^ (а) наилучших решений по каждому кри-

терию Ф^ (а) при заданной метрике р(Ф^ (а),Ф^), где Ф+ - экстремальное значение к-го критерия качества исследуемой (проектируемой) динамической системы, заранее известное или определяемое по ходу проведения вычислительных экспериментов; иначе говоря, отыскать область О^ (а), удовлетворяющую одновременно двум условиям: с вероятностью Р > Рз найти множество п ^ N

таких точек ат (т = 1,..., п), что в каждой из

этих точек выполнится неравенство

А

и при этом (n/N) > 1 - 5, здесь X е G^(X) , 0 <S << 1, 0 < s^ << 1,

Рз - заданная (требуемая) вероятность;

3) построение на основе определенных существенных параметров и выделенных подобластей регрессионных зависимостей;

4) выделение в ^-мерном пространстве критериев множества с высокой концентрацией точек Парето (или, если возможно, построение поверхности Парето); а в случае задания какой-либо схемы компромисса -выделение подобласти G0 (х) с G(x), содержащую наибольшую концентрацию требуемых компромиссных решений.

Метод ПЛП-поиска легко реализуется в среде MATLAB [8].

Подведем некоторые итоги. Доказывая эффективность ПЛП-поиска по сравнению со «слепым» способом поиска экстремумов, не предлагается отбрасывать последний. Однако очевидно, что ПЛП-поиск не только может помогать скорейшему поиску экстремумов, если это требуется, но и дает информацию о влиятельности варьируемых параметров и, что также важно, способствует возникновению вопросов, связанных с исследованием пространства варьируемых параметров (в частности, результаты ПЛП-поиска могут помочь выбрать эффективную схему компромисса). Итак, можно утверждать, что если принято решение исследовать сначала сформулированную задачу дискретным способом (что полезно даже в том случае, когда удается получить аналитические зависимости, но очень сложные), то ПЛП-поиск представляется весьма эффективным методом

компьютерных технологии в смысле ранее упомянутой их второИ составляющей. Кроме того, все алгоритмы ПЛП-поиска легко программируются, допускают диалоговый режим работы. То, что результаты применения ПЛП-поиска носят вероятностный характер, позволяет нам утверждать, перефразируя суждение известного французского философа ХХ в. Габриэля Марселя, высказанное в его диалоге с философом Полем Рикером [9], о его методе исследования природы и социальных систем, что ПЛП-поиск - это метод честных вопросов, а не фальшивых ответов.

Литература

1. Fitts P.M. The information capacity of the human motor system in controlling the amplitude of movement // Journal of Experimental Psychology. 1954. V. 47. №. 6. P. 381-391.

2. Зуев А.С. Графические интерфейсы как средства управления работой информационных систем // Информационные модели экономики. М.: МГАПИ, 2006. С. 80-84.

3. Метод статистических испытаний (метод Монте-Карло) / Бусленко Н.П. [и др.]; под ред. Ю.А. Шрейдера. М.: Физматгиз, 1962.

4. Статников И.Н., Андреенков Е.В. ПЛП-поиск -эвристический метод решения задач математического программирования. М.: ИИЦ МГУДТ, 2006.

5. Соболь И.М. Многомерные квадратурные формулы и функции Хаара. М.: Наука: ГРФМЛ, 1969.

6. Митропольский А.К. Техника статистических вычислений. М.: Наука: ГРФМЛ, 1971.

7. Шеффе Г. Дисперсионный анализ / пер. с англ. М.: Наука: ГРФМЛ, 1980.

8. Статников И.Н., Фирсов Г.И. ПЛП-поиск и его реализация в среде MATLAB // Проектирование инженерных и научных приложений в среде MATLAB. М.: ИПУ РАН, 2004. С. 398-411.

9. Марсель Г. Трагическая мудрость философии. Избранные работы. М.: Издательство гуманитарной литературы, 1995.

References

1. Fitts P.M. The information capacity of the human motor system in controlling the amplitude of movement // Journal of Experimental Psychology. 1954. V. 47. № 6. P. 381-391.

2. Zuev A.S. Graficheskie interfeysy kak sredstva upravleniya rabotoy informatsionnykh sistem // Informatsionnye modeli ekonomiki. M.: MGAPI, 2006. S. 80-84.

3. Metod statisticheskikh ispytaniy (metod Monte-Karlo) / Buslenko N.P. [i dr.]; pod red. Yu.A. Shreydera. M.: Fizmatgiz, 1962.

4. Statnikov I.N., Andreenkov E.V. PLP-poisk -evristicheskiy metod resheniya zadach matemati-cheskogo programmirovaniya. M.: IITs MGUDT, 2006.

5. Sobol' I.M. Mnogomernye kvadraturnye formuly i funktsii Khaara. M.: Nauka: GRFML, 1969.

6. Mitropol'skiy A.K Tekhnika statisticheskikh vy-chisleniy. M.: Nauka: GRFML, 1971.

7. Sheffe G. Dispersionnyy analiz / per. s angl. M.: Nauka: GRFML, 1980.

8. Statnikov I.N., Firsov G.I. PLP-poisk i ego reali-zatsiya v srede MATLAB // Proektirovanie in-zhenernykh i nauchnykh prilozheniy v srede MATLAB. M.: IPU RAN, 2004. S. 398-411.

9. Marsel' G. Tragicheskaya mudrost' filosofii. Iz-brannye raboty. M.: Izdatel'stvo gumanitarnoy li-teratury, 1995.

INTELLECTUALIZATION INFORMATION PROCESSING USING DIGITAL METHODS OF DYNAMIC SYSTEMS

I.N. Statnikov, G.I. Firsov

The application for solving multi-criteria method of synthesis of dynamic systems PLP search, which not only allows on the basis of the simulation model to carry out experiments show the parameter space in a given range of changes, but also as a result of the special nature of the planning of randomized experiments to apply quantitative statistical assessments of the varied parameters and their paired combinations on the analyzed properties of the dynamical system. It can be argued that if the decision is to explore the problem first formulated by a discrete manner (which is useful even in the case when it is possible to obtain analytical expressions, but very difficult), the PLP-Search is a very effective method of computer technology. In addition, all PLP-search algorithms are easy to program, allow dialog mode.

Key words: scheduling of computational experiments, the PL-grid, imitation simulation.

УДК 539.4

МОДЕЛИРОВАНИЕ ПРОЦЕССОВ РАЗРУШЕНИЯ КРЕМНИЯ ПРИ ДЕЙСТВИИ ВЫСОКИХ ЛОКАЛЬНЫХ НАПРЯЖЕНИЙ

А.И. Тюрин, И.А. Шуварин, А.М. Купряшкин, М.О. Воробьев

Цель работы - определение условий трещинообразования монокристаллического бездислокационного Si в зависимости от типа используемого индентора (Виккерса, Берковича и Бирбаума) и степени локального нагружения, при котором выполняются критерии определения параметров трещиностойкости. В работе представлены результаты моделирования процессов разрушения монокристаллического кремния методом динамического микро- и наноиндентирования. При использовании различных типов инден-торов смоделировано поведение исследуемого материала при различных видах локального механического взаимодействия. Установлены интервалы критических нагрузок и глубин отпечатка, при которых для каждого типа индентора выполняются критерии корректного расчета параметров трещиностойкости; определены числовые значения величин, характеризующих трещиностойкость материала. Показано, что внутри выявленных интервалов значения коэффициентов вязкости разрушения исследуемого монокристаллического кремния остаются постоянными (Kc=0,95±0,05 МПа-м1/2) и не зависят от типа применяемого индентора.

Ключевые слова: вязкость разрушения, трещинообразование, индентирование, нагружение.

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

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

Для корректного применения расчетных моделей локального разрушения в зоне укола