CC BY
32
Вычислительные технологии
Том 9, № 3, 2004
ОГРАНИЧЕНО ЛИ ДОПУСТИМОЕ МНОЖЕСТВО РЕШЕНИЙ ИНТЕРВАЛЬНОЙ СИСТЕМЫ?
И. А. ШАРАЯ
Институт вычислительных технологий СО РАН, Новосибирск, Россия
e-mail: shary@ict.nsc.ru http://www.ict.nsc.ru/lab1.2/irash
A criterion is proved for the tolerable solution set to be unbounded.
Введение
В работе используются обозначения интервального анализа, предложенные в [1], и классическая интервальная арифметика (см., например, [2, 3]).
Для интервальной системы линейных алгебраических уравнений
Ax = b,
где A 6 IMmxn — интервальная матрица размерности m х n; b 6 IRm — интервальный вектор длины m; x 6 — вещественный вектор длины n, допустимое множество решений £ описывается так [4]:
£ = {x | (VA 6 A) (3 b 6 b) (Ax = b)}.
Критерий принадлежности вектора x допустимому множеству решений на языке интервальной арифметики имеет вид (см., например, [5])
x 6 £ ^^ Ax С b. (1)
Легко показать [6], что £ — выпуклое многогранное множество в Можно ли по виду матрицы A судить об ограниченности множества £? Этот вопрос поставил в переписке профессор Зенон Кульпа (http://www.ippt.gov.pl/~zkulpa/) — автор диаграмма-тического метода исследования множеств решений интервальных систем линейных алгебраических уравнений [7, 8]. Поблагодарим его за интересную постановку и попробуем ответить.
1. Критерий неограниченности
Докажем сначала два вспомогательных утверждения. Первое позволяет по специальному виду матрицы сказать, что допустимое множество решений неограничено, а второе, наоборот, показывает, что если допустимое множество решений неограничено, то матрица имеет этот специальный вид.
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2004.
В соответствии с [9] будем называть конечное множество {а3} вещественных векторов линейно зависимым, если существует соответствующее множество {с3} вещественных чисел, не все из которых равны нулю, такое, что ^ а3с3 = 0. Заметим, что множество из
з
одного вектора линейно зависимо тогда и только тогда, когда этот вектор нулевой.
Утверждение 1. Пусть допустимое множество решений 2 непусто. Если в матрице А есть линейно зависимые вещественные столбцы, то 2 неограничено.
Доказательство. В силу (1) непустота 2 означает, что существует X £ Кп такой, что АХ С Ь. Расписав произведение АХ по столбцам А :з, получим
п
:3X С Ь. (2)
з=1
Пусть 3 — множество номеров линейно зависимых вещественных столбцов А :з. Тогда (2) можно переписать в виде
п
А .зХз + 5] А :зХз С Ь, (3)
з'€3 з=1
№
а линейную зависимость вещественных столбцов выразить формулой
А :зСз =0, (4)
зез
где сз £ К и ^ |сз| > 0. зез
Домножая (4) на произвольное вещественное добавляя к (3) и пользуясь законом дистрибутивности для вещественных чисел, получим
п
VI £ К 52 А :з (Хз + ¿Сз) + ^ А :зХз С Ь. (5)
з'ез з=1
з ^
Введем вектор с = (с1, . . . , сп), дополнив множество коэффициентов линейной зависимости нулевыми для ] £ 3. Тогда (5) можно переписать в виде
п
Vг £ к ^ А :з (Хз + ¿сз) с ь,
з=1
что в матричной форме выглядит так:
V г £ к а(х + гс) с ь.
По критерию принадлежности (1) это означает, что вместе с решением X в множество 2 попадает прямая, проходящая через X и параллельная вектору с. Значит, 2 неограничено. □
Утверждение 2. Пусть допустимое множество решений 2 непусто. Если оно неогра-ничено, то в матрице А есть линейно зависимые вещественные столбцы.
Доказательство. Как упоминалось во введении, 5 — это выпуклое многогранное множество. Если 5 неограничено, значит, неограничено его пересечение с каким-нибудь ортан-том. Тогда в этом ортанте лежит выпуклое многогранное неограниченное подмножество 5, из которого можно выбрать какой-нибудь луч (X + £с), где X — начало луча, с — вектор направления, £ £ — параметр, задающий точки луча.
Так как луч (X + ¿с) целиком лежит в 5, по критерию принадлежности (1)
V£ £ Е+ А(Х + ¿с) С Ь. С другой стороны, луч (X + ¿с) целиком лежит в одном ортанте, поэтому
V; £{1,...,п} х3с, > 0
и в (6) можно раскрыть скобки по правилу дистрибутивности [2, 3]:
V£ £ Е+ АХ + А(£с) С Ь. Для произвольных интервальных векторов х, у, ъ имеет место очевидное свойство
(6)
(7)
хСу х + ъ С у
|ъ| < widу,
где | | — модуль, wid — ширина интервальных векторов. Используя это свойство в (7) для х = АХ, ъ = А(£с), у = Ь, получим, что должно выполняться неравенство
V £ £ | А(£с) | < wid Ь.
(8)
Для вещественного £ применимо правило дистрибутивности [2, 3]. Это позволяет переписать (8) в виде
V£ £ |£(Ас)| < wid Ь.
Положительное £ можно вынести за знак модуля и разделить на него обе части неравенства. Тогда
> г ^ , . , wid Ь V£ £ Е+ \{0} |Ас| <
£
Это возможно только при | Ас| = 0, что эквивалентно
Ас = 0.
Равенство (9) означает, что
п
" =0,
(9)
ЕА 3 с
3=1
т. е. линейная комбинация столбцов интервальной матрицы А с коэффициентами, соответствующими компонентам вектора с, равна нулю. Остается заметить, что ненулевые коэффициенты в этой сумме могут стоять только при вещественных столбцах матрицы. В противном случае радиус линейной комбинации
rad (Ас) = ^^ |с, |rad А
(10)
3 = 1
п
будет отличен от нуля, что противоречит (9). □
Следствием утверждений 1 и 2 является
Критерий неограниченности. Пусть допустимое множество решений £ непусто. Оно неограничено тогда и только тогда, когда в матрице A есть линейно зависимые вещественные столбцы.
Можно ли представить, как выглядит непустое неограниченное £? Давайте попробуем. Пусть c £ — произвольный вектор коэффициентов, для которого линейная комбинация столбцов интервальной матрицы A равна нулю, т. е. Ac = 0. Обозначим через L пространство всех таких векторов.
В силу (10) у вектора c могут отличаться от нуля только компоненты, соответствующие вещественным столбцам матрицы A. Это значит, что размерность пространства L равна (k' — k''), где k' — число всех вещественных столбцов, k'' — максимальное число линейно независимых вещественных столбцов матрицы A. Например, если в матрице A все вещественные столбцы нулевые, то k'' = 0 и размерность пространства L равна числу нулевых столбцов k'.
Опираясь на доказательство утверждения 1, можно сказать, что £ представляет собой объединение прямых, параллельных произвольному вектору c из L. Следовательно, допустимое множество состоит из пространств, полученных параллельными сдвигами L. Для выпуклого многогранного множества это означает, что все его грани лежат в гиперплоскостях, параллельных L.
2. Выводы
На основании критерия неограниченности можно по виду матрицы A судить об ограниченности допускового множества решений, если оно непусто. Делается это так.
Если в матрице A нет линейно зависимых вещественных столбцов, то £ ограничено; если есть такие столбцы, то неограничено. В простейших случаях можно обойтись без вычислений, например:
— если в матрице A нет вещественных компонент, то £ ограничено;
— если в матрице A каждый столбец имеет невещественную компоненту, то £ ограничено;
— если в матрице A есть нулевые столбцы, то £ неограничено;
— если в матрице A есть пропорциональные вещественные столбцы, то £ неограничено;
— если в матрице A число вещественных столбцов больше числа строк этой матрицы, то £ неограничено.
Неограниченное допустимое множество решений представляет собой выпуклое многогранное множество, все грани которого параллельны пространству таких векторов c £ для которых Ac = 0. Размерность этого пространства равна (k' — k''), где k' — число всех вещественных столбцов, k'' — максимальное число линейно независимых вещественных столбцов матрицы A.
Список литературы
[1] KearfüTT R.B., Nakaü M.T., Neumaier A. ET AL. Standardized notation in interval analysis // Reliable Computing (в печати). (http://www.mat.univie.ac.at/~neum/software/int)
[2] Алефельд Г., Херцбергер Ю. Введение в интервальные вычисления. М.: Мир, 1987.
[3] Шарая И.А. О дистрибутивности в классической интервальной арифметике // Вы-числ. технологии. 1997. Т. 2, № 1. C. 71-83.
(http://www.ict.nsc.ru/lab1.2/Irene/ct97.ps)
[4] Шарый С.П. Алгебраический подход к анализу линейных статических систем с интервальной неопределенностью // Изв. РАН. Теория и системы управления. 1997. № 3. С. 51-61. (http://www.ict.nsc.ru/lab1,2/Shary/IzvAN.ps)
[5] SHARY S.P. Solving the linear interval tolerance problem // Mathematics and Computers in Simulation. 1995. Vol. 39. P. 53-85.
(http://www.ict.nsc.ru/lab1.2/Shary/LinTol.pdf)
[6] Шайдуров В.В., Шарый С.П. Решение интервальной алгебраической задачи о допусках // Красноярск, 1988 (Препр. АН. Сиб. отд-ние. Вычисл. центр. № 5).
[7] Kulpa Z. Diagrammatic analysis of interval linear equations. Pt I // Reliable Computing. 2003. Vol. 9, N 1. P. 1-20.
[8] Kulpa Z. Diagrammatic analysis of interval linear equations. Pt II // Reliable Computing. 2003. Vol. 9, N 3. P. 205-228.
[9] ХАлмош П. Конечномерные векторные пространства. М.: Физматгиз, 1963.
Поступила в редакцию 6 февраля 2004 г., в переработанном виде — 9 марта 2004 г.
Правила для Авторов
<http://www.ict.nsc.ru/mathpub/comp-tech/>
1. Статья должна быть представлена в редакцию в двух экземплярах в одной из двух форм:
а) рукопись, отпечатанная на одной стороне листа стандартного формата A4 + файлы рукописи в формате ЖТХ (2.09, 2е) и файлы рисунков на дискете;
б) рукопись, отпечатанная на одной стороне листа стандартного формата A4 + ее электронная версия, набранная в текстовом формате Microsoft Word (RTF), и файлы рисунков на дискете.
Время прохождения издательского цикла для рукописей, представленных в первой форме, минимально, во второй — максимально.
2. Все файлы предоставляются на дискете 3.5" формата 1440 Кбайт. Возможна пересылка файлов по электронной почте jct@ict.nsc.ru в виде *.zip архива.
3. На отдельной странице прилагаются на русском и английском языках название статьи, имена авторов, аннотация (не более 300 знаков) и ключевые слова.
4. Статья должна сопровождаться разрешением на опубликование от учреждения, в котором выполнена данная работа. В сопроводительном письме необходимо указать почтовый адрес, телефоны, e-mail автора, с которым будет проводиться переписка.
5. Для каждого автора должна быть представлена (на русском и английском языках) в виде отдельного файла следующая информация:
о Фамилия, имя, отчество о место работы и должность о ученая степень и звание о год рождения о почтовый адрес
о телефоны с кодом города (дом. и служебный), факс, e-mail, URL домашней страницы о область научных интересов (краткое резюме)
6. Материалы следует направлять по адресу: редакция журнала "Вычислительные технологии", Институт вычислительных технологий СО РАН, просп. Акад. Лаврентьева, 6, 630090, Новосибирск, 90, Россия, Игорю Алексеевичу Пестунову (отв. секретарь) — тел.: +7(3832)343785, e-mail: jct@ict.nsc.ru; Галине Григорьевне Митиной (зав. РИО).
Рекомендации по оформлению статьи в L-TeX-В редакцию следует направлять исходный файл, подготовленный в формате ЖТХ (версии 2.09, 2е) в стиле (классе) jctart (допускается использование стандартного стиля article).
Соответствующие файлы стиля jctart.sty, jctart-e.sty (для статей на английском языке) и класса jctart.cls можно скачать с сайта ЖВТ: http://www.ict.nsc.ru/win/mathpub/comp-tech/. 1. Структура файла в формате LT^X 2.09: \documentstyle[12pt,twoside]{jctart}
\setlength{\textwidth}{170mm}\setlength{\textheight}{240mm}
\begin{document}
\pagestyle{myheadings}
\markboth{<^ О. Фамилия автора(ов)>}{<КРАТКОЕ НАЗВАНИЕ СТАТЬИ (ДО 40 СИМВОЛОВ)>} ^^^{<НАЗВАНИЕ CTАТЬИ>\footnote{<Ссылка на поддержку (факультативно)>.}} \author{\sc{<^ О. Фамилия первого автора>}\\
\^{<Место работы первого автора>}\\[2mm] \sc{<^ О. Фамилия второго автора>}\\ \^{<Место работы второго автора>}\\[2mm] ...} \maketitle \begin{abstract} <Текст аннотации> \end{abstract} <Текст статьи> \begin{thebibliography} <Библиография (\item-список)>
\item {\sc Иванов~И.И., Иванова~И.И.} К вопросу о вычислительных технологиях // Вычисл. технологии. 1999. Т.~11, №~11. С.~1123--1135.
\end{thebibliography} \end{document}
(В конце файла даются:
<Перевод названия статьи на английский язык (или на русский, если статья на английском)> <аннотации на английский язык (или на русский, если статья на английском)>)
2. При подготовке исходного файла в системе L-TeX- 2е соответственно используются команда выбора класса документа и дополнительные пакеты макрокоманд, т. е. в структуре файла меняется только заголовок.
\documentclass[12pt,twoside]{jctart}
\setlength{\textwidth}{170mm}\setlength{\textheight}{240mm} \usepackage{amsmath}
\begin{document}
3. Список литературы составляется по ходу упоминания работы в тексте и оформляется по образцу: Книга
Шокин Ю.И. Метод дифференциального приближения. Новосибирск: Наука, 1979. 222 с. Бренстед А. Введение в теорию выпуклых многогранников: Пер. с англ. М.: Мир, 1988. Рояк М.Э., СоловЕйчик Ю.Г., Шурина Э.П. Сеточные методы решения краевых задач математической физики: Учеб. пособие. Новосибирск: Изд-во НГТУ, 1998.
Finlayson B.A. The Method of Weighted Resuduals and Variational Principles. N.Y.: Acad. Press, 1972. Книга четырех авторов
Проблемы вычислительной математики / А.Ф. Воеводин, В.В. Остапенко, В.В. Пивоваров, С.М. Шургин. Новосибирск: Изд-во СО РАН, 1995.
Статья из продолжающегося тематического сборника
Федорова А.А., Черных Г.Г. О численном моделировании струйных течений вязкой несжимаемой жидкости // Моделирование в механике: Сб. науч. тр. / РАН. Сиб. отд-ние. Вычисл. центр. Ин-т теор. и прикл. механики. 1992. Т. 6 (23). С. 129-140. Статья из журнала
Игнатьев Н.А. Выбор минимальной конфигурации нейронных сетей // Вычисл. технологии. 2001. Т. 6, № 1. С. 23-28.
Venkatakrishnan V. Newton solution of inviscid and viscous problems // AIAA J. 1989. Vol. 27, N 7. P. 285-291.
Труды конференции
Ivanov I.I. Problems in computational techologies // Intern. Conf. Comput. Techs. Novosibirsk, 1988. P. 225-229. Препринт
Гуськов А.Е., Федотов А.М., Молородов Ю.И. Информационная система"Конференции". Новосибирск, 2003 (Препр. РАН. Сиб. отд-ние. ИВТ. № 1-03). Диссертация
Деменков А.Г. Численное моделирование турбулентных следов в однородной жидкости: Дис. ... канд. физ.-мат. наук. Новосибирск, 1997. 123 с.
4. Иллюстрации вставляются в текст статьи с помощью команд special{em:graph <имя файла рисунка>} (LTEX 2.09) и includegraphics{<имя файла рисунка>} (LTEX 2е), например:
\begin{figure}[htbp]
\hspace*{<сдвиг рисунка по горизонтали в мм>mm} \special{em:graph <fig1.bmp>} \vspace*{<высота рисунка в мм>mm} \caption{<Подрисуночная подпись.>} \end{figure}
\begin{figure}[htbp] \centering
\includegraphics{fig1.eps} \caption{<Подрисуночная подпись.>} \end{figure}
Наиболее предпочтительной формой представления иллюстраций являются файлы рисунков в векторном формате PostScript (.eps) или черно-белых растровых в форматах .pcx, .bmp, .tif с разрешением 300 dpi.
Все надписи на рисунках (обозначение осей и т.п.) должны быть выполнены в том же начертании (гарнитура "Roman"), что и в тексте статьи. Латинские символы — курсивом, из математической моды (x[k], z х 10-3, ф, P,...), цифровые обозначения на графиках — наклонно (№ кривой — 1, 2,...), единицы измерения — по-русски (кг, м,...), цифры по осям — прямо (10, 15,...).
Instructions für Authors
<http://www.ict.nsc.ru/mathpub/comp-tech/>
1. Papers may be submitted to the editorial board in one of the following forms:
a) two copies of the manuscript typed on one side of the standard A4 sheet (297x210 mm) + files of the manuscript in LTEX (2.09, 2e) format and files of the figures on a diskette;
b) two copies of the manuscript typed on one side of the standard A4 sheet (297x210 mm) + files of the manuscript in Microsoft Word (RTF) format and files of the figures on a diskette.
The publishing cycle for manuscripts, submitted in the first form, is the longest one and for manuscripts in the second form — the shortest one.
2. All files should be submitted on a 3.5"floppy disc (1440 Kbytes) or sent to jct@ict.nsc.ru as a *.zip - archive.
3. A separate page should contain a title, names of the authors, an abstract (not more than 300 characters) and keywords.
4. The paper should be accompanied by the publication permission from the organization, where the work was done. The enclosed letter should contain the postal address, phone numbers and e-mail of the corresponding author.
5. A separate file should contain the following information on each author:
o First name, second name, last name o Affiliation, position o Academic degree and title o Year of birth o Postal address
o Office and home phone numbers (including area code), fax number, e-mail address, homepage URL o Scientific interests (brief curriculum vitae)
6. All materials should be sent to the following address: Dr. Igor A. Pestunov (executive secretary), Journal of Computational Technologies, Institute of Computational Technologies SB RAS, Academician Lavrentyev Ave. 6, Novosibirsk, 630090, Russia. Phone +7(3832)343785, E-mail: jct@ict.nsc.ru; Galina G. Mitina (publishing department manager).
Recommendations on submitting paper in LTeX
The source file should be submitted in LTeX format (2.09, 2e versions) using jctart style file (class file) (standard article style (class)can also be used).
The files of appropriate style — jctart.sty, jctart-e.sty (for papers in English) and jctart.cls class file can be downloaded from JCT web site: http://www.ict.nsc.ru/win/mathpub/comp-tech/. 1. The file structure in LTeX 2.09 format: \documentstyle[12pt,twoside]{jctart}
\setlength{\textwidth}{170mm}\setlength{\textheight}{240mm}
\begin{document}
\pagestyle{myheadings}
\markboth{<Name(s) of author(s)>}{<SHORT TITLE (LESS THAN 40 CHARACTERS)>} \title{<TITLE OF PAPER>\footnote{<Reference to supporting organization (optional)>.}} \author{\sc{<Name of the first author>}\\
\it{<Affiliation of the first author>}\\[2mm] \sc{<Name of the second author>}\\
\it{<Affiliation of the second author>}\\[2mm] ...} \maketitle \begin{abstract} <Abstract> \end{abstract} <Text of paper> \begin{thebibliography} <References (\item-спнсок)>
\item {\sc Ivanov~I.I., Ivanova~I.I.} On computational technologies // Computational technologies. 1999. Vol.~11, No.~11. P.~1123--1135.
\end{thebibliography} \end{document}
2. When submitting the source file in LTeX 2e format, the documentclass command and additional packages are used, therefore only the header is changed.
\documentclass[12pt,twoside]{jctart}
\setlength{\textwidth}{170mm}\setlength{\textheight}{240mm} \usepackage{amsmath}
\begin{document}
3. A list of the references should be sorted according to the order of citations in the text and it should be written as in the following example:
Book
Finlayson B.A. The method of weighted residuals and variational principles. N.Y.: Acad. Press, 1972. Book by four authors
Problems of computational mathematics / A.F. Voevodin, V.V. Ostapenko, V.V. Pivovarov, S.M. Shurgin. Novosibirsk: SB RAS Publishing House, 1995.
Paper from continued subject transactions
Fedorova A.A., Chernykh G.G. On numerical modelling of viscous incompressible jet fluid flows // Modelling in mechamics: Scientific transactions / RAS. Siberian branch. Computing Center. Institute of Theoretical and Applied Mechanics. 1992. Vol. 6 (23). P. 129-140. Paper from journal
Venkatakrishnan V. Newton solution of inviscid and viscous problems // AIAA J. 1989. Vol. 27, N 7. P. 285-291.
Conference proceedings
Ivanov I.I. Problems in computational techologies // Intern. Conf. Comput. Techs. Novosibirsk, 1988. P. 225-229.
Dissertation
Demenkov A.G. Numerical modelling of turbulent wakes in homogeneous fluid: Dissertation for degree of candidate of physical and mathematical sciences. Novosibirsk, 1997. 123 p.
4. Figures should be included into the text using commands \special{em:graph <figure file name>} (LTEX 2.09) and \includegraphics{<figure file name>} (LTEX 2e), for example:
\begin{figure}[htbp]
\hspace*{<horizontal shift of figure in mm>mm} \special{em:graph <fig1.bmp>} \vspace*{<height of figure in mm>mm} \caption{<Figure caption.>} \end{figure}
\begin{figure}[htbp] \centering
\includegraphics{fig1.eps} \caption{<Figure caption.>} \end{figure}
The preferred presentation form for illustrations is a figure file in vector format PostScript (.eps) or black and white bitmap formats .pcx, .bmp, .tif with 300 dpi resolution.
All figure inscriptions (axes definitions, etc.) should be done by the same font as in the text of paper ("Roman"type family). Latin characters should be done in italics in mathematical mode (x[k], zx10-3, —, P,...), figures on axes — by straight font.
In papers, which are written in Russian, the units of measurement should be written in Russian.
В ближайших номерах / Forthcoming papers
Petrovskaya N.B. The Analysis of The Grid Refinement Procedure for Vector Functions
ПЕТРОВСКАЯ Н.Б. Анализ процедуры разбиения сетки для многокомпонентного решения
Воропаева О.Ф. Численные модели динамики безымпульсного турбулентного следа в устойчиво стратифицированной среде
VOROPAYEVA O.F. The numerical models of momentumless turbulent wake dynamics in a stably stratified medium
Денисенко В.В. Применение энергетических методов при построении векторных полей
DENISENKO V.V. Energy methods application for vector fields calculation
Мартюшов С.Н., Мартюшова Я.Г. Численное моделирование струйных течений методом конечного объема на основе TVD-схемы 2-го порядка точности
MartyuSHOV S.N., MARTYUSHOVA Y.G. Numerical simulation of jet flows on the basis of 2-nd order accuracy Harten's scheme
Михеев С.Е. Выпуклая квадратичная аппроксимация MlHEEV S.E. Convex Quadratic Approximation
НемировскиЙ Ю.В., Янковский А.П. Интегрирование задачи динамического упругопластического изгиба армированных стержней переменного поперечного сечения обобщенными методами Рунге — Кутты
NEMIROVSKY J.V., YANKOVSKY A.P. Integration of a problem dynamic elastic-plastic bending of reinforced cores of variable cross-section section with help generalized Runge — Kutta methods
ТАТОСОВ А.В. Схема расчета нестационарных течений газа в пневматической системе
TATOSOV A.V. The calculation scheme nonstationari gas motion in pneumatic system
Объединенный семинар "Информационно-вычислительные технологии". Аннотации докладов за весенний семестр 2004 года
Joint seminar "Computational technologies". Abstracts of presented talks. Lent term 2004.
