Допусковое множество решений для интервальных систем уравнений со связанными коэффициентами Текст научной статьи по специальности «Математика»

Вычислительные технологии

Том 14, № 3, 2009

Допусковое множество решений для интервальных систем уравнений со связанными коэффициентами*

И. А. ШАРЛЯ, С. П. ШАРЫй Учреждение Российской академии наук Институт вычислительных технологий СО РАН, Новосибирск, Россия e-mail: sharaya@ict.nsc.ru, shary@ict.nsc.ru

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

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

Введение

Интервальной системой линейных алгебраических уравнений (ИСЛАУ) вида Ax = b с интервальной матрицей коэффициентов A G IMmxn и интервальной правой частью b G IRm называется семейство вещественных систем линейных уравнений вида Ax = b, в котором матрица коэффициентов пробегает интервальную матрицу A, а правая часть — интервальный вектор b.

Для интервальной линейной системы уравнений Ax = b можно определить различные множества решений [1-3]. Наиболее популярны объединенное, допусковое и управляемое множества решений.

В этой работе рассмотрим допусковое множество решений (ДМР). Оно обозначается через Stoi(A, b) и состоит из всех возможных векторов x, таких что для всякой матрицы коэффициентов A из интервальной матрицы A значение произведения Ax не выходит за границы интервального вектора b:

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

* Публикуется при финансовой поддержке Президентской программы "Ведущие научные школы РФ" (грант № НШ-931.2008.9).

© ИВТ СО РАН, 2009.

(1)

Интервальной системой линейных алгебраических уравнений Ах = Ь со связью & на коэффициенты системы назовем семейство уравнений Ах = Ь, в котором матрица коэффициентов пробегает множество А П &, а правая часть — интервал Ь.

Допусковым множеством решений для интервальной системы Ах = Ь со связью & на коэффициенты назовем множество

Из определений (1) и (2) очевидно, что при А П & = 0 наложение связи на коэффициенты способствует расширению допускового множества решений:

Исследованиям интервальных линейных систем уравнений со связанными параметрами посвящено много публикаций. Среди их авторов — И. Рон, Х. Янссон, З.М. Румп, А. Ноймайер, Л.В.Колев, Е.Д. Попова, Б.С. Добронец, С.П. Шарый. Целая серия работ опубликована в соавторстве Г. Алефельдом, В. Крейновичем и Г. Майером. Но до сих пор в работах, посвященных интервальным линейным системам со связанными параметрами, рассматривалось только объединенное множество решений. Данная статья посвящена не объединенному, а допусковому множеству решений для ИСЛАУ со связанными параметрами.

Структура статьи такова: в разделе 1 приведены обозначения, необходимые понятия и факты, в разделе 2 изложен и обоснован метод отыскания ДМР для ИСЛАУ с выпуклой многогранной связью на коэффициенты, в разделах 3 и 4 предложены упрощенные варианты этого метода для частных случаев связи, в разделе 5 рассмотрен числовой пример на применение упрощенных вариантов метода и для демонстрации влияния связи коэффициентов ИСЛАУ на допусковое множество решений.

1. Обозначения, необходимые понятия и факты

В этой работе, начиная с введения, мы следуем обозначениям, принятым в [4]. В частности, Ж := {[х, у] | х,у £ К, х < у} — это множество интервалов на вещественной оси, *Ж := {[х,у] | х € К и {-то}, у € К и {то}, х < у} — множество расширенных интервалов.

Особенности использования шрифтов в статье следующие. Вещественные числа, векторы и матрицы набраны обычным математическим курсивом. Например, х £ — вещественный вектор длины п. Каллиграфический шрифт применяется для множеств, например, S С Ктхп — подмножество в множестве всех прямоугольных вещественных матриц размера т на п. Жирный курсив, в соответствии с [4], используется для интервалов, интервальных векторов и интервальных матриц. Например, А £ Жтхп — интервальная матрица, имеющая т строк и п столбцов. Такая интервальная матрица А представляет собой в Ктхп множество специального типа — прямоугольный параллелепипед, размерность которого не выше тп, а ребра параллельны координатным осям.

Нижний индекс для вектора обозначает его компоненту или, с геометрической точки зрения, проекцию этого вектора на соответствующую координатную ось, например, для интервального вектора Ь £ Жт имеем Ь = (Ь\, Ь2,... , Ьт)Т. Для матрицы пара нижних

0, если А П& = 0,

, иначе.

(2)

А П& = 0 Ь) С ^о!(А П&, Ь).

индексов указывает на ее элемент, например, для Б £ Ктхп через Б^ выражен элемент, стоящий на пересечении ¿-й строки и ]-го столбца, т. е. проекция Б на -ю координатную ось. Двоеточие в качестве нижнего индекса обозначает весь диапазон значений этого индекса, например, для Б £ Ктхп через Б^ обозначается ¿-я строка матрицы Б, т. е. проекция матрицы Б на координатное подпространство ¿-й строки:

БГ: = (Бг1) Бi2, • • • , Бгп) •

Аналогично нижний индекс для множества будем использовать для проекции множества на координатное подпространство, например, для Б С Ктхп через Б^ будем выражать ортогональную проекцию множества Б на координатное подпространство ¿-й строки, т. е. множество всех возможных значений ¿-й строки для матриц из множества Б:

Si: = (Б, | Б £5} = |(Ба,...,Бт) | Б

Символ "©" будем применять для поэлементного умножения множеств, например, Б © х = (Бх | Б £ 5} для Б С Ктхп и х £ Кп — это множество всех произведений Бх, в которых вещественная матрица Б является элементом множества Б. Символ " будет, как обычно, обозначать прямое произведение, например, для множества Б С 1Кп запись Б = Бз- означает, что Б совпадает с прямым произведением своих проекций з

на координатные оси.

Нам понадобятся некоторые определения и факты из выпуклого анализа, приведенные ниже.

Множество в конечномерном векторном пространстве над полем К называется выпуклым многогранным множеством, если оно может быть получено как пересечение конечного числа замкнутых полупространств или совпадает со всем пространством. То есть выпуклое многогранное множество в Кп (или Ктхп) — это такое множество, которое может быть описано конечной системой нестрогих линейных неравенств. Пересечение конечного числа выпуклых многогранных множеств является выпуклым многогранным множеством. Ограниченное выпуклое многогранное множество называется выпуклым многогранником.

Точка выпуклого многогранного множества называется его вершиной, если в этом множестве нет отрезка, для которого она служит внутренней точкой. Множество вершин выпуклого многогранного множества конечно. Множество всех вершин выпуклого многогранника V обозначим ует! V. Всякий выпуклый многогранник является выпуклой оболочкой множества своих вершин:

V = еопу(уег1 V).

Пусть с, х £ Кп, d = £ Включение-принадлежность

стх £ d, (3)

или, что то же самое, двустороннее нестрогое линейное неравенство вида

^ < стх <

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

при с = 0, й = й — множество, ограниченное двумя параллельными гиперплоскостями;

при с = 0, й = й — гиперплоскость; при с = 0, 0 € d — все пространство; при с = 0, 0 € d — пустое множество. Множество решений конечной системы элементарных линейных включений является пересечением конечного числа гиперполос. Элементарное линейное включение вида (3), в котором интервал d вырожден (т. е. й = й), является линейным уравнением.

Наряду с элементарными нам понадобятся и матричные линейные включения, т. е. включения вида

Сх € d, где х € Кп, С € Ктхп, d € 1Кт.

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

стх € d, где с € Кп, d € *Ж,

будем называть линейным ограничением. При d € Ж — это элементарное линейное включение, при d = [й, то] или d = где й € К, — это нестрогое линейное

неравенство стх > й или стх < й соответственно.

2. Метод отыскания ДМР для ИСЛАУ

с выпуклой многогранной связью на коэффициенты

2.1. Переход к решению включения V © х С Ь, где V — выпуклый многогранник

Обратимся к отысканию допускового множества решений интервальной линейной системы уравнений Ах = Ь со связью 0 на коэффициенты системы. Прежде всего преобразуем определение этого множества, избавившись в (2) от кванторной приставки. Получим

'0, если А ПО = 0, ^(А П0, Ь) = {, _ , (4)

х € Кп | (А ПО) © х С Ь |, иначе.

Теперь рассмотрим пересечение множеств А и 0 .В качестве связи 0 мы договорились брать только выпуклые многогранные множества. Интервальная матрица А представляет собой выпуклый многогранник (ограниченное выпуклое многогранное множество) в Ктхп. Поэтому пересечение А П 0 является выпуклым многогранным множеством (как пересечение двух выпуклых многогранных множеств) и ограничено (так как А ограниченное множество). Значит, А П 0 представляет собой выпуклый многогранник. Обозначим его через V.

Таким образом, для отыскания ДМР системы Ах = Ь с выпуклой многогранной связью 0 на коэффициенты нам надо научиться решать относительно х включение

V © х С Ь, где V — выпуклый многогранник в Ктхп.

2.2. Свойства включения V © x Ç b

В этом подразделе рассмотрим несколько свойств включения V © x Ç b, которые помогут нам построить метод его решения. Первое из свойств имеет место для произвольного множества S матриц в Rmxn, а не только для выпуклого многогранника V.

Свойство 1. Пусть S Ç Rmxn, x G Rn, b G IRm. Включение

S© x Ç b (5)

и система из m включений

& (Si: © x Ç bi) (6)

г

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

Доказательство. По определению операции "©" имеем S © x = {Sx | S G S}.

Поэтому включение S © x Ç b можно переписать в виде & (Sx G b).

SeS

Для интервального вектора b всякое матричное линейное включение Sx G b представляет собой систему элементарных линейных включений & (Si:x G bi). Значит, S©x Ç b

i

— это сжатая запись для

& & (Srx G bi). (7)

ses i

Перегруппируем элементарные линейные включения в (7):

& & (Srx G bi). i ses

Каждый блок вида & (Si:x G bi) заменим на & (Si:x G bi). При этом мы просто ses Si: eSi:

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

& & (Si:x G bi). (8)

i Si: eSi:

Теперь представим выражение (8) в сжатом виде, воспользовавшись для этого еще раз определением операции " ©":

& (Si: © x Ç bi).

i

Смысл свойства 1 в том, что переход от (5) к (6) состоит в исключении заведомых повторов элементарных линейных включений.

Свойство 2. Пусть V Ç Rmxn — выпуклый многогранник, x G Rn, b G IRm. Тогда

V © x Ç b ^^ (vert V) © x Ç b. (9)

Рис. 1. К доказательству свойства 2

Доказательство (рис. 1).

Очевидно, поскольку vert V Ç V.

^=) Множество b выпукло, поэтому для всякого его подмножества взятие выпуклой оболочки не выводит из множества b:

(vert V) © x Ç b conv((vert V) © x) Ç b.

Умножение на вектор x — линейное отображение. Оно переводит всякую выпуклую комбинацию конечного числа точек в комбинацию образов этих точек с теми же коэффициентами. Поэтому

(conv(vert V)) © x Ç conv((vert V) © x).

Остается использовать тот факт, что многогранник V представляет собой выпуклую оболочку множества своих вершин: conv(vert V) = V. □

Свойство 2 позволяет перейти от включения V © x Ç b, обозначающего обычно бесконечную систему элементарных линейных включений, к включению (vert V) © x Ç b, всегда обозначающему конечную систему элементарных линейных включений. Суть этого перехода — удаление тех элементарных линейных включений, которые заведомо являются следствиями оставшихся.

Утверждение 1. Пусть V С Rmxn — выпуклый многогранник, b G IRm — интервальный вектор, x G — вещественный вектор. Тогда

V © x Ç b ^ & ((vert(Vi:)) © x Ç b*). (10)

Доказательство. Воспользовавшись свойством 1, перепишем включение V©x Ç b в виде

& (Vi: © x Ç bi).

г

Проекции Vi:, i = 1,... , m, выпуклого многогранника V сами являются выпуклыми многогранниками. Поэтому для них имеет место свойство 2:

Vi: © x Ç bi ^^ (vert(Vi:)) © x Ç bi. □

В утверждении 1, как и в свойстве 2, мы получили конечную систему элементарных линейных включений, эквивалентную включению V © х С 5. Давайте сравним эти системы. Для удобства сравнения запишем их в сходном виде, выделяя в каждой системе в г-й блок строки с правой частью 5г:

Система Из формулы Исходный вид Блочный вид

1 (9) (уех! V) © х С 5 & ( & ^х С 5г)

г V V €уей V /

2 (10) & ((уег1(Н:)) © х С 5г) & ( & то С 5г)

г г Ч'иеуег^Ю /

При таком блочном представлении, где к тому же выписаны явно все элементарные линейные включения, становится очевидным, что система 2 является подсистемой системы 1. Поясним это с помощью рис. 2. Строки коэффициентов г-го блока системы 1 соответствуют проекциям Vi: всех вершин V многогранника V. А строки коэффициентов г-го блока системы 2 соответствуют только вершинам проекции многогранника V. Для каждой вершины многогранника можно указать проецируемую в нее вершину многогранника V. При этом часть вершин многогранника V может проецироваться в точки, которые не являются вершинами многогранника Может получиться и так, что несколько вершин многогранника V проецируются в одну вершину многогранника Итак, из двух сравниваемых эквивалентных конечных систем элементарных включений система 2 предпочтительнее, поскольку она является подсистемой системы 1.

Заметим, что утверждение 1 удобно использовать при нахождении ДМР для ИСЛАУ без связей. Действительно, исключив в определении (1) кванторную приставку, получим, что е4о1(А, 5) совпадает с множеством решений включения А © х С 5. Интервальная матрица А — это выпуклый многогранник, и из утверждения 1 следует, что (А, 5) описывается системой

& (7уег1;(А:)) © х С 5г) . (11)

° | — у:, где v € уей V О — v € уег^г:)

г

Рис. 2. Проекции вершин многогранника

V

и

вершины его проекции

Множество уег^А^.) содержит не более 2П элементов, поэтому запись (11) обозначает систему элементарных линейных включений, в которой не более т-2га строк. Например, допусковое множество решений для системы

[1, 2] [3, 4] 5 [6,7]

х =

[0, 6] [-1, 0]

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

1x1 + 3х2 € [0, 6], 2x1 + 3x2 € [0, 6], 1х1 + 4х2 € [0, 6], 2х1 + 4х2 € [0, 6], 5х1 + 6х2 € [-1, 0], ^5х1 + 7х2 € [-1, 0].

2.3. Метод отыскания ДМР

Вернемся от рассмотрения свойств включения V © х С Ь к поиску допускового множества решений 24о1(А П Ь) для ИСЛАУ Ах = Ь с выпуклой многогранной связью $ на коэффициенты системы. Утверждение 1, доказанное в подразделе 2.2, завершает следующую цепочку рассуждений.

1. Сначала в подразделе 2.1, преобразовав определение (2), мы получили равенство (4). Оно означает, что если А П $ пусто, то ^о1(А П $, Ь) тоже пусто, а если А П $ непусто, то Е^01(А П Ь) совпадает с множеством решений включения (А П $) © х С Ь относительно неизвестного х.

2. Далее в подразделе 2.1 мы отметили, что для выпуклого многогранного множества $ пересечение А П $ является выпуклым многогранником.

3. И наконец, на основании утверждения 1 стало возможным заменить включение

(А П $) © х С Ь, где А П $ — выпуклый многогранник, конечной системой элементарных линейных включений, сжатая форма которой

& ((уеЛ((А П £)*)) © х С Ьг).

Эта цепочка рассуждений

1) доказывает, что допусковое множество решений для ИСЛАУ Ах = Ь с выпуклой многогранной связью $ на коэффициенты системы с геометрической точки зрения представляет собой пересечение конечного числа гиперполос, потому что оно либо пусто, либо описывается конечной системой элементарных линейных включений;

2) позволяет предложить следующий метод отыскания множества !Е4о1(А П 0, Ь):

(12)

Мы не будем конкретизировать этап II метода (12) из-за обилия соответствующих ему постановок задач (среди которых — получение оптимальных двойственных описаний и отыскание оценок множествами различной формы и с разными требованиями к расположению по отношению к оцениваемому множеству), а также из-за того, что эти задачи уже достаточно известны. Договоримся считать, что мы нашли множество П 0, Ь), если нам удалось выполнить этап I.

Для выполнения этапа I метода (12) можно предложить следующий порядок действий:

Этап ¡метода (12) — универсальный вариант

1. Найти описание выпуклого многогранного множества 0 в виде конечной системы линейных неравенств.

2. Включить в эту систему неравенства, выражающие принадлежность интервальной матрице А.

3. Найти из полученной системы множество вершин многогранника А П 0.

4. Если полученное множество пусто, то Е^сДА П 0, Ь) = 0. Если множество вершин многогранника А П 0 непусто, спроецировать это множество на координатные подпространства отдельных строк и удалить из проекций лишние точки.

Этот путь достаточно универсальный, но трудоемкий.

Учитывая конкретный вид начальных данных, этот путь можно сократить. В разделах 3 и 4 мы рассмотрим два специальных случая связи 0 и предложим для этих случаев более простые, по сравнению с универсальным, способы выполнения этапа I метода (12).

3. Случай 1

3.1. Особенности выпуклой многогранной связи 0

Рассмотрим случай, когда выпуклая многогранная связь 0 удовлетворяет следующим трем условиям.

Метод отыскания допускового множества решений для интервальной линейной системы уравнений Ах = Ь с выпуклой многогранной связью 0 на коэффициенты системы

Этап I. Найти множества уегЦ(А П 0)г:), г = 1,... , т.

Если хотя бы одно из них пустое, то Е^^А П 0, Ь) = 0. Если все множества непусты, перейти к этапу II. Этап II. Дать описание множества решений

конечной системы элементарных линейных включений

& ((уеЛ((А П 0)*)) © х с Ьг) в нужном для использования виде.

1. Множество 0 задано в виде конечной системы линейных ограничений:

О € 0 ^ & (v С^Оц € ¿И , (13)

г,]

где 5 € N Скц € К, dk € *Ж.

2. Каждый элемент матрицы О может входить с ненулевым коэффициентом только в одно из линейных ограничений системы (13), т.е.

Сщ = 0 (Ук €{1,...,5}Ш) (Скг] =0).

3. Никакие два элемента из одной строки матрицы О не могут входить с ненулевыми коэффициентами в одно и то же линейное ограничение системы (13), т.е.

Скгг = 0 (V? € {1, . . . , п} \ {Г}) (Скг] = 0).

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

О] + О]г = 0, г = 1,..., п — 1, ? = г + 1,..., п. 3.2. Влияние особенностей связи на этап I метода (12)

Обозначим как 1пх множество всевозможных пар индексов (г?) для элементов матрицы О:

1пх = {(г?) | г = 1,..., т, ? = 1,..., п}.

Через 1пхк обозначим множество тех пар индексов (г?), для которых О] входит в к-е ограничение с ненулевым коэффициентом:

1пхк = {(г?) € 1пх | Сг]к = 0}.

Пересечение множеств А и 0 обозначим буквой Б.

В силу первого условия на множество 0 матрицы Б из множества Б описываются системой линейных ограничений

& Скг] Бг] € ^

к=1--9 \(г])&пхк

& (Бг] € Аг]).

_ (г])е1пх

Пользуясь вторым условием на множество 0, разобьем эту систему на блоки так, чтобы разные блоки не имели общих переменных:

Скг] Бг] € dk ;

& { (г])е1пхк

к=1'""9 {Бг] € Аг], (г?) € 1пхк; (14)

& (Бг] € Аг]).

(г]) е 1пх 1пхк

к

Здесь отдельный блок соответствует каждому индексу к. Кроме того, отдельный блок, состоящий из одного ограничения £ А^, соответствует каждой паре индексов (гу), для которой = 0). При таком разбиении очевидно, что значения переменных

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

(Vi)

S

Определим множество из единственного блока системы (14), содержащего переменную .

Если (гу) £ 1пх^, то единственный блок, содержащий переменную , имеет вид

CW Sir £ dk,

(ir)einxfc

£ Air, (1r) G Inxk . Перепишем (15), выделяя переменную Sj:

£ ( dk —

У^ Ckir Sir

(ir) e inxfc \{(j)}

Sir G Air, (Ir) G Inxk \{(ij)},

Sij £ *

Учитывая специальный вид ограничения (16с), имеем

(15)

(16a)

(16b) (16c)

(17)

где — гу-я проекция множества решений системы включений (16а) и (16Ь).

Ограничения (16а) и (16Ь) позволяют мыслить множество как множество значений многозначной функции

(1г) е 1пхй \{(у)}

на интервале . При таком взгляде на становится очевидно, что

(гг)етпхк \{(г^)}

ij

и

dk - Ckir Sir

SlreAlr, \ (ir) e Inxk \{(ij)}

(ir) e Inxfc \{(ij)}

kij *

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

S.

i? 1 ^k ^ ^ CkirAir Ckij *

(ir) e inxk \{(ij)}

dk -

Ckir Ai

(18)

(Этот же результат для множества можно получить более длинным путем, выполнив добросовестно все преобразования метода Фурье для исключения переменных Б¡г, (/г) € 1пхк\{(г;)}, из системы включений (16а) и (16Ь).)

Таким образом, для (г?) € 1пхк множество определяется из (17) и (18). Если

же (г^') € 1пх \У 1пхк, то единственный блок, содержащий переменную Б^, имеет вид к

Б^ € А^, и потому = А^. В целом,

АуП dfc - CkirЛг / Ckij , если (ij) G Inxfc,

Sij = { VV (ir)eInxk \{(ij)} // / (19)

Aj, если (ij) G Inx \ U Inxk .

k

Множество Б^ пусто или является интервалом.

Итак, мы показали, что для выпуклой многогранной связи G, соответствующей случаю 1, множество (АП^)г-. является прямым произведением множеств Б^, вычисляемых по правилу (19). Поэтому этап I метода (12) для случая 1 можно выполнить следующим образом:

Этап I метода (12) для случая 1

1. Найти множества Sij, (ij) G Inx, по правилу (19).

2. Если хотя бы одно из этих множеств пусто, то Stoi(A П G, b) = 0. Если все множества Sij, (ij) G Inx, отличны от пустого, определить vert ((А П G)i:) для каждого индекса i

как множество вершин интервального вектора Sij и перейти к этапу II метода (12). j

(20)

3.3. Пример постановки задачи

В качестве примера постановки задачи, соответствующей случаю 1, рассмотрим задачу об отыскании допускового множества решений для интервальной модели межотраслевого баланса с условием рентабельности отраслей. Аналогичный пример, но только для объединенного множества решений, был решен И. Роном в [5]. Возьмем уравнение Леонтьева межотраслевого баланса в виде

(I - д)х = у, (21)

в котором п — число секторов (отраслей) модели, I — единичная матрица размера п, Q € Мпхп — матрица коэффициентов прямых затрат, х € Кп — вектор объемов производства, у € Кп — вектор объемов конечного продукта. Заметим, что число 1 — ^^ Qij, равное сумме элементов j-го столбца матрицы I — Q, представляет собой добавленную стоимость на единицу производства в j-м секторе экономики. Относительно уравнения (21) рассмотрим следующую задачу.

Задача

Дано:

1) для значений коэффициентов прямых затрат Qгj известны только границы, т. е. дана такая интервальная матрица Q £ Жгахга, что Q £ ф;

2) для компонент уг вектора объемов конечного продукта тоже указаны только границы, т. е. задан интервальный вектор у £ Жга;

3) заданы границы добавленной стоимости на единицу производства в каждом секторе экономики в виде интервального вектора d £ Жга.

Найти

все такие векторы объемов производства х, при которых для всех матриц коэффициентов прямых затрат Q из ф, удовлетворяющих условиям 1 — г Qгj £ dj, ] = 1,...,п, вектор объемов конечного продукта у не выйдет за границы интервального вектора у.

(22)

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

Покажем, что связь G соответствует условиям случая 1.

1. Связь G описывается конечной системой линейных ограничений:

(е с

С ££ ^ , & > Сг, £ dj

(При согласовании с (13) надо считать, что Сщ = 1 при к = ] и 0 — во всех остальных случаях. Затем надо заменить к на ] и исключить ненужные индексы.)

2. Каждый элемент Сг, матрицы С входит только в одно ограничение, а именно в ограничение с номером ].

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

Итак, все требования к связи выполнены, и для решения задачи можно воспользоваться методом (12) в упрощенном варианте, соответствующем случаю 1.

Этап I метода (12)

1. С учетом конкретной связи правило (19) для определения множеств , где Б = А П£, примет вид

= Аг, П I dj ^^ Apj I

V р=г )

Поскольку А = I—ф, это правило можно переписать в терминах начальных данных:

(1 — фгг) П ^ + Е , если г =

V Р=г )

( —) П — 1 + Е , иначе.

\ Р=г /

2. Если какое-то из множеств , (г;) € 1пх, пусто, то множество решений задачи (22) тоже пусто.

Если все множества , (г;) € 1пх, непустые, то множество решений задачи (22) описывается системой

& ( (уег!(Бг:)) © х с у Л , где Б = 0 , г = 1,..., п. (23)

Этап II метода (12)

Дать описание множества решений конечной системы элементарных линейных включений (23) в нужном для использования виде. На этом этапе можно учесть типичное для задачи (22) требование неотрицательности вектора объемов производства, добавив к системе (23) неравенство х > 0.

4. Случай 2

4.1. Особенности выпуклой многогранной связи $

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

1. Множество $ задано в параметрической форме

О(р), где р — вектор параметров,

причем каждый элемент О^ (р) матрицы О(р) пропорционален одной из компонент вектора параметров:

Сгз (р) = Рк(г,з), % € К \ {0}, рк(^) € {р1, . . . ,ря}.

2. Никакие два элемента матрицы О(р), пропорциональные одной компоненте вектора параметров, не лежат в одной строке, т. е. для всех к = 1,... , ^ имеет место следование ( ) ( )

(г1) € 1пхк (V; €{1,...,п}\{/}) ((г;?) € 1пхк),

где 1пхк := {(г;) | Оу(р) = сз-рк} — множество пар индексов всех тех элементов матрицы О(р), которые пропорциональны параметру рк.

Очевидно, что в рассматриваемом случае множество $ — линейное подпространство в Етхп

В виде, соответствующем случаю 2, можно представить, в частности, следующие множества матриц: симметричные матрицы, кососимметричные матрицы, циклические матрицы, матрицы Ганкеля, матрицы Гурвица, матрицы Тёплица. Например, для множества кососимметричных матриц одна из возможных параметризаций, соответствующих случаю 2, имеет вид

О(р) =

р1 р2 рз . . рп

- р2 рп+1 рп+2 . . р2п-1

-рз -рп+2 р2п . . рзп-з

\-рп

-р2га-1 -р3п-3

\

р1^. . , рп(п+1)/2 € К

рп(п+1)/2/

4.2. Влияние особенностей связи на этап I метода (12)

Если связь G удовлетворяет первому условию случая 2, то множество S := AflG состоит из всех таких матриц S, что

& (Sij = Cjpfc, (ij) G Inxfc), fc=l,...,q (24)

_ Sjj G Aj, i = 1,..., m, j = 1,..., n.

Заметим, что в (24) не известны не только значения переменных Sj при i = 1,... , m и j = 1,... , n, но и значения параметров p. при k = 1,... , q. Перепишем (24) в виде системы из q блоков:

& |Sij = QjPfc, (ij) G Inxk, (25)

k=1...... \Sjj G Ajj, (ij) G Inxfc .

В (25) никакие два блока не имеют общих переменных, поэтому значения переменных Sj, (ij) G Inxfc, из k-го блока не зависят от значений переменных из других блоков.

В силу второго требования к множеству G элементы одной строки матрицы S не могут входить в один блок. Поэтому для каждого i переменные Sj, j = 1,... , n, будут

независимы между собой, т.е. (Vi) (Si: = Sj).

j

Найдем множество Sj из единственного блока системы (25), содержащего переменную Sj. Если (ij) G Inxfc, это будет блок

Sir = Qrpfc, (/r) G Inxfc, Sir G Air, (1r) G Inxfc .

После очевидных эквивалентных преобразований он примет вид

Sir = QrPfc, (/r) G Inxfc, CirPfc G Air, (1r) G Inxfc .

Равенства Sir = Qrpk, (/r) G Inxfc, при pk G R описывают в пространстве переменных Sir, (/r) G Inxfc, прямую, проходящую через начало координат. Так как все константы Cir, (/r) G Inxfc, отличны от нуля, эта прямая не перпендикулярна никакой из координатных осей. Каждое включение qrpk G Air при qr = 0 ограничивает множество значений параметра pk отрезком Air/qr. Учитывая действие всех таких ограничений, получим, что множество значений параметра pk равно Р| Air/qr. Это множество либо пусто,

(ir)einxfc

либо является интервалом. Итак,

Sir = CirPk, (/r) G Inxk, Pk G f| Air/Cir.

(ir)einxfc

Поэтому множество значений переменной Sj, удовлетворяющих этой системе, можно найти по правилу

Sjj = Cjj © I Pi Air/Cir I для (ij) G Inxk . (26)

\(ir)einxfc J

Заметим, что если группа элементов матрицы С, пропорциональных параметру , состоит только из одного элемента (у), то множество Б^ совпадает с А^:

¡ПХЛ = {(у)} вгз _ Ау. (27)

Мы показали, что этап I метода (12) можно выполнить следующим образом:

Этап I метода (12) для случая 2

1. Найти множества Sj, i = 1,... , m, j = 1,... , n, по правилу (26).

2. Если хотя бы одно из них пустое, то £toi(A П G, b) = 0.

Если все множества Sj, i = 1,... , m, j = 1,...,n, отличны от пустого, определить vert ((A П G)i:) для каждого индекса i как множество вершин интервального вектора Sj и перейти к этапу II метода (12). j

(28)

5. Числовой пример

Рассмотрим числовой пример, чтобы пояснить, как пользоваться методом (12) в упрощенных вариантах.

Пусть интервальная матрица А и интервальный вектор Ь имеют вид

А _ ло, 1] [-5,-т Ь _ л-1,1] А _ \[0, 2] [1, 2] )> Ь _ \[—2, 2]

Требуется представить графически множество таких вещественных векторов х £ Кга, для которых при всех кососимметричных матрицах А из А значение Ах лежит в интервале Ь.

5.1. Первый способ решения

5.1.1. Выбор описания связи

Множество G кососимметричных матриц размером 2 х 2 представим в виде

С _ с12) ^^ С12 + С21 _ 0 уС21 С22)

Такое описание связи G удовлетворяет условиям случая 1, при этом д _ 1, _ 0, С112 _ С121 _ 1, С111 _ С122 _ 0, 1пх1 _ {(12), (21)}. Для связи такого типа этап I метода решения (12) можно выполнить в упрощенном варианте (20).

5.1.2. Этап I в виде (20)

1. Найдем множества , г _ 1, 2, ] _ 1, 2, по правилу (19). Получим 5ц _ Ац _[0,1],

Б12 _ А12 П (0 - А21) _ [-5,-1] П (-[0, 2])_[-5,-1] П [-2,0] _ [-2,-1],

521 _ А21 П (0 - А12) _ [0, 2] П (-[-5, -1]) _ [0, 2] П [1, 5] _ [1, 2],

522 _ А22 _[1, 2].

2. Все множества Sj, i =1, 2, j = 1, 2, непустые, поэтому

(A П G) 1: = (Sil S12) = ([0,1] [-2,-1]) , (A П G)2: = (S21 S22) = ([1, 2] [1,2]) .

Отсюда

vert((A П G) 1:) = {(0,-2), (1,-2), (0,-1), (1,-1)}, vert((A nG)2:) = {(1,1), (2,1), (1, 2), (2, 2)}.

5.1.3. Этап II

На этом этапе нам надо представить графически множество решений конечной системы элементарных линейных включений

-2X2 е [-1,1],

x1 - 2x2 е [-1,1],

-X2 е [-1,1],

x1 - x2 е [-1,1],

X1 + X2 е [-2, 2],

2X1 + X2 е [-2, 2],

X1 + 2X2 е [-2, 2],

2x1 + 2x2 е [-2, 2].

Рис. 3. Множество решений числового примера (заштриховано)

Первую строку системы разделим на -2, последнюю строку — на 2, третью строку (как следствие первой строки системы) и пятую строку (как следствие последней строки) удалим. Получим систему элементарных линейных включений:

1 1"

Х2 е - 2, 2

х1 - 2х2 е [-1,1]

х1 - х2 е [-1,1]

2x1 + Х2 е [-2, 2]

Х1 + 2x2 е [-2, 2]

Х1 + Х2 е [-1,1]

(I)

(II)

(III)

(IV)

(V)

(VI)

множество решений которой представлено на рис. 3. В рассмотренном нами числовом примере множество решений — это выпуклый шестиугольник с вершинами (0,0.5), (0.5,0.5), (1, 0), (0, -0.5), (-0.5, -0.5), (-1, 0).

5.2. Второй способ решения 5.2.1. Выбор описания связи

Множество G кососимметричных матриц размером 2 х 2 представим в виде

Такое описание связи G удовлетворяет условиям случая 2, поэтому можно применить упрощенный вариант (28) этапа I метода (12).

5.2.2. Этап I в виде (28)

1. Поиск множеств , г =1, 2, ] = 1, 2.

Множества 1пх1 = {(11)} и 1пх3 = {(2 2)} содержат по одному элементу. Поэтому на основании (27) получаем Бц = Ац = [0,1], Б22 = А22 = [1, 2].

Множества Б12 и Б21 определим по правилу (26), опираясь на то, что 1пх2 = {(12), (21)}, С12 = 1, С21 = -1:

¿12 = С12 © ((А12/С12) п (А21/С21)) = 1 © ([-5, -1] П [-2, 0]) = [-2, -1],

¿21 = С21 © ((А12/С12) П (А21 /С21)) = -1 © ([-5, -1] П [-2, 0]) = [1, 2].

2. Как в первом способе решения.

5.2.3. Этап II

Как в первом способе решения.

5.3. Влияние связи

Во введении было замечено, что при А П О = 0 наложение связи на коэффициенты способствует расширению допускового множества решений:

А ПО = 0 Вы (А, Ь) С ВЫ(А ПО, Ь).

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

На рис. 4 представлены множества в^А, Ь) и в4ы(А П О, Ь) для случая, когда

. /[0,1] [-5,-1]\ . /[-1,1]\ ,, А = 1 [о 2] [1 2] ) , Ь = 1 [ 2 2П , 0 — множество кососимметричных матриц.

Множество в4о1(А П О , Ь) взято из рассмотренной числовой задачи. Множество 54о1(А, Ь) можно определить, в соответствии с (11), из системы элементарных линейных включений:

-5x2 е [-1,1],

XI — 5х2 е [-1,1],

-Х2 е [-1,1],

XI — Х2 е [-1,1],

Х2 е [-2, 2],

2X1 + Х2 е [-2, 2],

2X2 е [-2, 2],

2х1 + 2х2 е [-2, 2].

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

Другой способ построения множества £toi(A, b) — воспользоваться программой, расположенной по адресу http://www.nsc.ru/interval/Programming/AEsolset.ps и предназначенной для визуализации множеств решений интервальной системы Ax = b с матрицей A размером 2 х 2.

Множество Stoi(A, b) представляет собой выпуклый шестиугольник с вершинами (0, 0.2), (0.8, 0.2), (1, 0), (0, -0.2), (-0.8, -0.2), (-1,0). На рис. 4 он заштрихован для сопоставления с содержащим его шестиугольником, соответствующим множеству Stoi(A П G, b).

Список литературы

[1] Шарый С.П. Алгебраический подход к анализу линейных статических систем с интервальной неопределенностью // Известия РАН. Теория и системы управления. 1997. № 3. С. 51-61. http://www.nsc.ru/interval/shary/Papers/IzvAN.pdf

[2] Шарый С.П. Конечномерный интервальный анализ: электронная книга. http://www.nsc.ru/interval/Library/InteBooks/SharyBook.pdf

[3] Shary S.P. A new technique in the systems analysis under interval uncertainty and ambiguity // Reliable Computing. 2002. Vol. 8, N 5. P. 321-419.

http://www.nsc.ru/interval/shary/Papers/ANewTech.pdf

[4] Kearfott R.B., Nakao M.T., Neumaier A., Rump S., Shary S., van Hentenryok P. Standardized notation in interval analysis // Интервальный анализ: Труды XIII Байкальской международной школы-семинара «Методы оптимизации и их приложения». Т. 4. Иркутск: ИСЭМ СО РАН. 2005. С. 106-113. http://www.nsc.ru/interval/INotation.pdf,

http://www.nsc.ru/interval/Conferences/Baikal-2005/IntervalAnalysis.pdf, а также http://www.mat.univie.ac.at/~neum/interval.html

[5] Rohn J. Interval linear systems with prescribed column sums // Linear algebra and its applications. 1981. Vol. 39. P. 143-148.

Поступила в редакцию 9 февраля 2009 г.