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

УДК 51-7;519.6;519.8

ДИСКРЕТНАЯ МОДЕЛЬ СТРУКТУРНОГО РАЗРУШЕНИЯ СЛОЖНЫХ СИСТЕМ1

A.A. Кочкаров

(1)Институт прикладной математики им. М.В. Келдыша РАН, г. Москва;

(2)Карачаево-Черкесская государственная технологическая академия, г. Черкесск

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

ВВЕДЕНИЕ

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

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

Ряд аварий в электроэнергетических системах в крупных городах России (Москва, 2005 г.), Европы (Лондон, 2003 и 2006 гг.; Париж, 2006 г.) США и Канады (Детройт, Нью-Йорк, Кливленд, Оттава,

1 Работа выполнена при поддержке РФФИ (проект № 07-01-00618) и РГНФ (проект № 05-03-03188).

Торонто, 2003 г.), показал, что развитие чрезвычайных ситуаций в коммуникационных системах с сетевой структурой происходит по «принципу домино» (в случае электроэнергетических систем — это веерные отключения). Один вышедший из строя объект (элемент системы) сильно повышает вероятность аварии на остальных, что приводит к возникновению лавины аварий. О последствиях таких аварий красноречиво свидетельствует следующий факт.

14 августа 2003 г. в ряде крупнейших городов восточного побережья США и Канады — Детройте, Нью-Йорке, Кливленде, Оттаве, Торонто и др. — произошло девятисекундное отключение электроснабжения, приведшее к веерным отключениям электроэнергии на площади более 24 тыс. кв. км и получившее название «Блэкаут-2003». Причина — перегрузка на энергетическом каскаде Ниагара — Мохок (на границе США и Канады). Авария затронула более 50 млн чел. в восьми штатах США и провинции Онтарио Канады, привела к остановке более 100 электростанций, в том числе 22-х атомных реакторов. Ликвидация последствий аварии заняла более 30-ти часов. Сумма нанесенного США финансового ущерба составила не менее 6 млрд. долл. США.

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

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

Описанная в настоящей работе математическая модель структурного разрушения сложной системы представляет собой дополнение к модели распространения внешних воздействий по структуре системы [3]. Эти модели в схеме развития чрезвычайных ситуаций (инициирование ЧС ^ развитие ЧС ^ выход ЧС за пределы системы) описывают первый и второй этапы.

1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ

СТРУКТУРНОГО РАЗРУШЕНИЯ

Изменения, происходящие в структуре сложной системы, могут быть описаны простейшими теоретико-графовыми операциями [4]: стягиванием ребра, удалением (добавлением) ребра, удалением (добавлением) вершины. Изменения структуры системы могут быть разовыми, а могут быть постоянными (периодическими, регулярными). Для последнего случая разумно ввести понятие структурной динамики — изменение структуры системы с течением времени. Несомненно, для описания структурной динамики лучше всего подходит аппарат теории графов [4].

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

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

Обозначим через G = (V, E) граф, соответствующий структуре исследуемой системы, V — множество вершин, соответствующее элементам системы, а Е — множество ребер, соответствующее

связям между элементами системы. Каждой вершине V е V припишем веса ^(у) и w (V), отражающие текущую загрузку и предельную загрузку элемента системы. В случае, когда текущая загрузка w(v) элемента системы достигает предельного значения w (у), то элементы системы выходят из строя. А проходящие через него потоки перераспределяются по «соседним» элементам системы. Выход из строя элемента системы в теоретико-графовой терминологии соответствует удалению из графа системы вершины с инцидентными ей ребрами. А перераспределение весов в тривиальном случае соответствует равному разделению веса

w (у) удаленной вершины по вершинам, смежным с удаляемой.

При выходе из строя одного или нескольких элементов системы возможны несколько сценариев дальнейшего развития событий. Один из них, если система функционирует в предельном состоянии, т. е. загрузка элементов близка к предельному значению, то возможен «быстрый» переход системы в критическое состояние. Структурное разрушение, вообще говоря, процесс динамический. Не нарушая общности, будем считать, что wt(v) — текущая загрузка вершины V е V в момент времени t = 0, 1, 2,..., Т, ... Если через V* = {у*} с V,

j = 1 2, ..., V |, обозначить множество вершин, вышедших из строя в момент времени t, т. е. те, у

которых wt(v.) 1 w (у.), а через £(( () = {у/} — окру-

t 1 1 / У

жение вершины у/ (или множество вершин, смежных с вершиной у/), |^(( )| = degVj, 1 = 1, 2, ...,

|^(( )|, то процесс структурного разрушения формально описывается следующим образом.

В момент времени t = 0 необходимо произвести проверку по всем вершинам у е V, и сформировать

множество Р1 из вершин, для которых справедливо Wo(V/) 1 w (). Во все последующие моменты времени t = 1, 2,..., Т, ... следует воспользоваться правилом

^ + ( ) = ) + Е;W ( ^

1 = 1, 2,..., |$( )|, j = 1, 2,..., |й|.

Если wt + 1(у/) 1 w(у/), то вершина у\ удаляется из графа G и добавляется в множество V* + 1.

Коэффициент Е; — параметр распределения загрузки. Он может зависеть от различных факторов, в простейшем случае он равномерно распределяет предельную загрузку удаляемой вершины по соседним, т. е. для каждой вершины у- он вычисляется как Е; = 1/deg у/. Структурное разрушение при

таком значении параметра распределения загрузки будем называть равномерным.

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

2. ХАРАКТЕРИСТИКИ И ОСОБЕННОСТИ СТРУКТУРНОГО РАЗРУШЕНИЯ СЛОЖНОЙ СИСТЕМЫ

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

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

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

Для исследования процесса структурного разрушения систем предлагается использовать следующие критерии отказа.

Критерий полного разрушения ст0(к). Система считается вышедшей из строя, если в системе выйдут из строя все элементы (будут удалены все вершины графа — структуры системы). Критерий полного разрушения ст0(к) зависит от одного параметра к — числа удаленных вершин в начальный момент времени структурного разрушения.

Критерий связности ст1(к). Система считается вышедшей из строя, если нарушена связность ее структуры при удалении вершин. Критерий связности <г1(к) зависит от одного параметра к — числа удаленных вершин в начальный момент времени структурного разрушения.

Компонентный критерий ст2(к, т). Система считается вышедшей из строя, если число компонентов в структуре системы при ее разрушении станет равным (или больше) заданного числа т. Компонентный критерий ст2(к, т) выхода системы из строя зависит от двух параметров: к — числа удаленных вершин в начальный момент времени структурного разрушения и (т — 1) — максимально допустимого числа компонентов структуры при ее разрушении.

Диаметральный критерий <г3(к, В). Система считается вышедшей из строя, если диаметр хотя бы одного из компонентов структуры системы в процессе разрушения окажется меньше заданного значения В. Диаметральный критерий <г3(к, В) выхода системы из строя зависит от двух параметров: к — числа удаленных вершин в начальный момент времени структурного разрушения и В — минимально допустимого диаметра компонентов структуры при ее разрушении.

Множество Ф^) элементов, вышедших из строя (удаленных из структуры) в момент времени t = 1, будем называть множеством эпицентров структурного разрушения. В критериях ст0(к), ст1(к), ст2(к, т) и ст3(к, В) число к соответствует числу эпицентров структурного разрушения системы.

В настоящей работе раскрыты некоторые аспекты равномерного структурного разрушению ациклических графов (цепей и деревьев) с равными значениями начальных загрузок w0(v) и равными значениями предельных загрузок w (у) для всех их вершин.

3. СТРУКТУРНОЕ РАЗРУШЕНИЕ ГРАФОВ-ЦЕПЕЙ

3.1. Структурное разрушение графов-цепей

по критерию связности

Всякий связный ациклический граф называется деревом. Частным случаем дерева G является граф-цепь С = (Ус, Ес). Множество вершин VC графа-цепи С состоит из двух висячих вершин — концов цепи со степенями, равными единице, и внутренних вершин со степенями, равными двум.

Рассмотрим граф-цепь С = (VC, Ес), |VC| = п, с равными для всех его вершин у е Vс весами w0(v),

W (у).

Ввиду того, что всякий граф-цепь утратит связность при удалении хотя бы одной невисячей (внутренней) вершины, справедлива

Лемма 1. Всякий граф-цепь С = (Vс, Ес), | ^с| = п, будет разрушен по критерию <з1(к), где 1 < к < п — 2, при удалении хотя бы одной невисячей вершины. ♦

3.2. Структурное разрушение графов-цепей

по компонентному критерию

Лемма 2. Всякий граф-цепь С = (Vс, Ес), | ^с| = п, п 1 3, будет разрушен по критерию ст2(к, т), где

2 < т т (п + 1)/2 при нечетном п и 2 < т < п/2 при четном п, если число попарно несмежных внутренних вершин-эпицентров равно к = т — 1.4

Доказательство. Удаление одной невисячей вершины графа-цепи С = (Vс, Ес) непременно приведет к распадению его на два компонента С' = (VC ,, Ес) и С" = (^ „, Ес„). «Простейшим» компонентом в таком случае может быть граф-вершина. Всякий граф-цепь С = (VC, Ес), |VC| = п, п 1 3, имеет (п — 1)/2 попарно несмежных внутренних вершин, если п — нечетное, и п/2 — 1, если п — четное. Удаление из графа-цепи С = (Ус, Ес) всех попарно несмежных внутренних вершин приведет к появлению (п + 1)/2 компонентов в первом случае и п/2 — во втором, причем каждый компонент будет представлять собой граф-вершину. ♦

3.3. Структурное разрушение графов-цепей по диаметральному критерию

Пусть H = (W, Q) есть n-вершинный связный граф [4]. Длина кратчайшей цепи, соединяющей пару вершин w, v е W, называется расстоянием между вершинами w и v [4] и обозначается через p(w, v). Заметим, что введенное таким образом расстояние удовлетворяет известным аксиомам евклидовой метрики.

Для фиксированной вершины w е W величина

sw = max p(w, v) называется эксцентриситетом

v е W

вершины w е W. Максимальный среди всех эксцентриситетов вершин графа H = (W, Q) называется диаметром графа H и обозначается через d(H),

т. е. d(H) = max s(w). Если пара вершин u, w е W

v е W

соединяется кратчайшей цепью длины p(u, w) = = d(H), то эта цепь называется диаметральной. Вершина w называется периферийной, если s(w) = d(H). Радиус графа H обозначается через r(H) и вычисляется по формуле r(H) = min s(w). Вершина w

w е W

называется центральной, если s(w) = r(H). Центром графа H называется множество центральных вершин.

Лемма 3. Всякий граф-цепь C = (VC, EC), | VC| = n, n l 3, будет разрушен по критерию ст3(1, r(C)) при удалении центральной вершины (т. е. когда эпицентром является центральная вершина). Причем диаметры появившихся в результате структурного разрушения компонентов будут равны r(C) — 1 и d(C) - r (C) - 1>

Доказательство. Для графа-цепи C = = (VC, EC), |VC| = n, расстояние между двумя его висячими вершинами vj и v2 совпадает с диаметром и длиной графа-цепи p(vp v2) = d(C) = n — 1, а сами вершины vj и v2 являются периферийными.

Всякий граф-цепь C = (VC, EC), | VC| = n, имеет центр, состоящий из двух вершин, если число вершин цепи | VC| = n — четное, и состоящий из одной вершины, если число вершин | VC| = n — нечетное. В первом случае r (C) = n/2, во втором r (C) = = (n —1)/2. Удаление центральной вершины графа-цепи C = (VC, EC) приведет к распадению его

на две цепи C1 = (Vc , Ec ) и C2 = (, Ec ) с

соответствующими диаметрами d(CJ) = r (C) — 1 и d(CJ) = d(C) — r (C) — 1. Очевидно, что при нечетном n диаметры компонентов Cj = (, Ec ) и

C2 = ( (, Ec2 ) будут равны: d(C1) = d(C2), а при n четном d(CJ) > d(C2)>

Для эксцентриситета всякой внутренней вершины у е Vс графа-цепи С = (Vс, Ес) справедливо неравенство г (С) < е(у) < С(С). Так же, как и в случае с удалением из графа-цепи центральной вершины, удаление любой внутренней вершины приведет к распадению графа-цепи С = (Vс, Ес) на

два компонента С1 = (( , Ес ) и С2 = ( , Ес ),

с соответствующими диаметрами С(С1) = е(у) — 1 и С(С2) = С(С) — е(у) — 1. Поэтому из леммы 3 очевидным образом вытекает

Лемма 4. Всякий граф-цепь С = (Ус, Ес), | ^| = п, п 1 3, будет разрушен по критерию ст3(1, е(у)) при удалении внутренней вершины у е Vс (т. е. когда эпицентром является внутренняя вершина у е Vс). Причем диаметры появившихся в результате структурного разрушения компонентов будут равны е(у) — 1 и С(С) — е(у) — 1.4

3.4. Структурное разрушение графов-цепей по критерию полного разрушения

Рассмотрим вопрос о разности начальной (текущей) w0(v) и предельной загрузки w (у) элементов графа-цепи С = (VC, Ес).

Предположим, что эпицентром на графе-цепи С = (Vс, Ес) является одна из его внутренних вершин у * е Vс. Пусть также w (у *) — w0(v *) > w (у *)/2, тогда после удаления эпицентра у * е Vс текущая загрузка смежных ему вершин у, у" е Vс станет равной

Wl (у') = Wl (у") = w0(v') + w(у*)/2 =

= W0(v *) + w (у *)/2 < w (у *).

Поэтому процесс структурного разрушения прекратится, а его длительность будет Тсг = 1.

В другом случае, когда w (у *) — w0(v *) < w (у *)/2, длительность процесса структурного разрушения Тсг > 1, если система не выйдет из строя при достижении установленного критерия разрушения.

Иная ситуация складывается, когда эпицентром является одна из висячих вершин у1, у2 е Vс

графа-цепи С = (Ус, Ес). В таком случае w (у1) —

— w0(v1) < w (у1), поэтому для смежной эпицентру вершине у"' е Vс будет справедливо выражение

w1(v"') = w0(v"') + w (у1) > w (у"). А это значит, что вершина у'" е Vс выйдет из строя (будет удалена) в следующий момент времени t = 2. Причем вершина у"' окажется висячей для цепи С/у1, что при-

ведет к удалению смежной ей вершине в момент времени t = 3. Этот процесс будет продолжаться до тех пор, пока не будут удалены все вершины графа-цепи С = (Vс, Ес), т. е. будет достигнут критерий разрушения ст0(1). Из проделанных рассуждений вытекает

Лемма 5. Всякий граф-цепь С = (Vс, Ес), | Vс| = п, п 1 3, будет разрушен по критерию ст0(1) при удалении одной из висячих вершин за время Тсг = п. 4

4. СТРУКТУРНОЕ РАЗРУШЕНИЕ ДЕРЕВЬЕВ

4.1. Структурное разрушение деревьев

по критерию связности

У деревьев Т = ^т, Ет), как и у графов-цепей С = ( Ус, Ес), центр может состоять либо из одной, либо из двух вершин. Если диаметральная цепь дерева Т = (Vт, Ет) имеет четную длину, т. е. диаметр с (Т) четный, то центр дерева состоит из одной вершины и из двух в противном случае, т. е. когда диаметр С(Т) нечетный. У всякого дерева Т = (VT, Ет) не менее чем две висячие вершины. Все остальные, как в случае с графами-цепями, будем называть внутренними.

Поскольку всякое дерево утратит связность при удалении хотя бы одной невисячей (внутренней) вершины, то справедлива

Лемма 6. Всякое дерево Т = (Кр Ет), | Vт| = п, будет разрушено по критерию ст1(к), где 1 < к < п — пт, пт — число висячих вершин, при удалении хотя бы одной внутренней вершины за время Тсг = 1.4

4.2. Структурное разрушение деревьев по компонентному критерию

Лемма 7. Всякое дерево Т = (^ Ет), | VT| = п, будет разрушено по критерию ст2(1, т) при удалении одной внутренней вершины у е VT за время Тсг = 1, причем т = deg(v).♦

Доказательство. Удаление из дерева Т = = Ет) хотя бы одной внутренней вершины приведет к распадению его на компоненты, причем число компонентов будет зависеть от степени удаляемой вершины. ♦

Лемма 8. Всякое дерево Т = (Vт, Ет), | VT| = п, будет разрушено по критерию ст2(к, т) при удалении к попарно несмежных внутренних вершин у. е VT за

k

время Тсг = 1, причем т = ^ (deg(vг.) — 1) + 1.4

I = 1

Доказательство. Удаление из дерева Т = = (Vт, Ет) всех к внутренних вершин у. е VT про-

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

После удаления первой вершины у1 є Ут дерево распадется на deg(v1) компонентов. Поскольку эпицентры являются попарно несмежными и невисячими вершинами дерева Т = (Ут, Ет), то вершина у2 є Ут, принадлежащая какому-то из полученных при удалении вершины у1 компоненту также не будет являться для своего компонента висячей вершиной. Поэтому при удалении вершины у2 компонент, которому она принадлежит, распадется на deg(v2) компонентов. А общее число компонентов, на которые распадется само дерево Т = (Ут, Ет), после удаления вершин у1 и у2 станет равным

deg(v1) + deg(v2) - 1 =

= ^(^) - 1) + ^^) - 1) + 1.

И далее, каждое удаление одной из вершин V. є Ут, і = 3, 4,..., к, будет увеличивать число компонентов на deg(vг.) — 1. А это значит, что при одновременном удалении всех эпицентров, соответствующих условиям леммы, дерево Т = (Ут, Ет) k

распадется на т = ^ (deg(vг.) — 1) + 1 компонен-

і = 1

тов. Длительность процесса структурного разрушения составит Т = 1. ♦

ЗАКЛЮЧЕНИЕ

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

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

ЛИТЕРАТУРА

1. Управление риском / В.А. Владимиров, В.В. Кульба, Г.Г. Малинецкий и др. — М.: Наука, 2000.

2. Архипова Н.И., Кульба В.В. Управление в чрезвычайных ситуациях. — М.: РГГУ, 1998.

3. Кочкаров А.А., Малинецкий Г.Г. Управление безопасностью и стойкостью сложных систем в условиях внешних воздействий // Проблемы управления. — 2005. — № 5. — С. 70—76.

4. Лекции по теории графов / В.А. Емеличев, О.И. Мельников, В.И. Сарванов, Р.И. Тышкевич. — М.: Наука, 1990.

5. Форд Л., Фалкерсон Д. Потоки в сетях. — М.: Мир, 1966.

6. Тоффоли Т., Марголус Н. Машины клеточных автоматов. — М.: Мир, 1991.

в (495) 250-79-71, e-mail: Azret_Kochkarov@mail.ru

Статья представлена к публикации членом редколлегии

Г.Г. Малинецким.

книги

Алексеев В.Е. Графы и алгоритмы. Структуры данных. Модели вычислений. — М.: Бином, 2006. — 318 с.

Иванов Г.Е. Слабо выпуклые множества и функции. — М.: ФИЗМАТЛИТ, 2006. — 351 с.

Камаев В.А. Технологии программирования. — М.: Высшая школа, 2005. — 359 с.

Квасов Б.И. Методы изогеометрической аппроксимации сплайнами. — М.: ФИЗМАТЛИТ, 2006. — 360 с.

Климонтович Ю.Л. Турбулентное движение и структура хаоса. — М.: URSS, 2007. — 323 с.

Корис Р. Справочник инженера-схемотехника. — М.: Техносфера, 2006. — 607 с.

Левин Б.Р. Вероятностные модели и методы в системах связи и управления. — М.: Радио и связь, 1985. — 312 с.

Ориентация и навигация подвижных объектов. — М.: ФИЗМАТЛИТ, 2006. — 422 с.

Пантелеев А.В. Спектральный метод анализа нелинейных стохастических систем управления. - М.: Вузовская книга, 2006. — 392 с. Разработка моделей и методов повышения эффективности управления некоммерческими организациями. — М.: ИПУ, 2007. — 75 с.

Шестаков А.А. Обобщённый прямой метод Ляпунова для систем с распределёнными параметрами. — М.: URSS, 2007. — 318 с. Юрко В.А. Введение в теорию обратных спектральных задач. — М.: ФИЗМАТЛИТ, 2007. — 384 с.

Васильев А.Н. Научные вычисления в Microsoft Excel. — М.: Диалектика, 2004. — 512 с.

Гусев В.С. Яндекс: эффективный поиск. Краткое руководство. — М.: Диалектика, 2007. — 224 с.

Йордон Э. Объектно ориентированный анализ и проектирование систем. — М.: Лори, 2007. — 264 с.