ПАРАДОКСЫ ВЫЧИСЛИТЕЛЬНОЙ ЛИНЕЙНОЙ АЛГЕБРЫ И СПЕКТРАЛЬНЫЕ ПОРТРЕТЫ МАТРИЦ Годунов С.К. Эти слайды были приготовлены в феврале 2011 г. для конференции, проведенной. - презентация

1 ПАРАДОКСЫ ВЫЧИСЛИТЕЛЬНОЙ ЛИНЕЙНОЙ АЛГЕБРЫ И СПЕКТРАЛЬНЫЕ ПОРТРЕТЫ МАТРИЦ Годунов С.К. Эти слайды были приготовлены в феврале 2011 г. для конференции, проведенной Президиумом РАН к 100- летию со дня рождения М.В. Келдыша

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

3 В основу вычислительной линейной алгебры естественно положить Постулат: Только такие числовые функцииот матрицы можно вычислять, для которых справедливо неравенство в котором- известная функция При этом условии, знаяи точность можно дать гарантированную оценку точности для вычисленной Пример вычислимой функции - число обусловленности матрицы где если Хорошо известны алгоритмы решения системы линейных уравнений, при выполнении которых одновременно с решением вычисляется 3

4 Изложение понятия о решении системы уравнений Обычно начинается с введения определителя Реальное вычисление определителя приводит к серьёзным проблемам: матрицам Точные неравенства т.е. для любого N существуют матрицы ПРИМЕР: 4

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

6 M - оценка амплитуды L - характерное время (декремент затухания) Типичное поведение затухающих решений При в оценке решения Можно ли это считать устойчивостью? Если A – NxN матрица то И.М. Гельфанд, Г.Е. Шилов, 1958 г. 6

7 Теорема Островского (о непрерывной зависимости ) Если все элементы матрицыи матрицы подчинены неравенствам то для каждого найдется такое что В нашем случае Пример теореме Островского не противоречит. Формальная непрерывность имеет место. 7

8 Определение-спектра принадлежит -спектру, если Спектральный портрет матрицы A 8

9 9

10 Спектральные портреты симплектических матриц Рассмотрим симплектическую матрицу вида: Матрицы С, S, P имеют следующую структуру: 10 Изучим поведение спектральных портретов при изменении параметра t

11 Спектральные портреты симплектических матриц 11

12 Еще один поучительный пример (к вопросу о расчёте собственных значений матриц) 12

13 Эксперимент: Эксперимент: Собственные числа матрицы С найденные с использованием пакетов MATLAB, MAPLE, SCILAB и библиотеки IMSL 13 MATLABMAPLESCILABIMSL

14 В действительности Точные значения : ВСЕ собственные значения вычисленные при помощи пакетов являются точными точками спектра матрицы С, при Этот ε-спектр покрывает круг 14

15 Решение систем линейных уравнений с матрицей MATLAB SCILAB Решения получены с помощью коммерческого MATLAB и свободно распростроняемого SCILAB (НГУ, ИМ СО РАН)

16 Резюме проведенного обсуждения Стоит ли заниматься расчетом ??? Нет гарантии, что их можно вычислить с приемлемой точностью. (речь идет о несимметричных матрицах ) ВОПРОС: Зачем в приложениях интересуются ??? ОТВЕТ: Часто требуется убедится, что или, что на прямой нет Предлагается решать более общий вопрос: Есть ли на той или иной кривой ? Если кривая не проходит через то всюду на этой кривой Для гладкой кривой конечной длины при этом 16

17 Удобно критерий отсутствия на кривой формулировать как Для кривых конечной длины предполагается, что Важное неравенство Критерий дихотомии спектра кривой Дихотомия спектра 17

18 a Спектральные зоны – полосы содержащие точки спектра Одномерный спектральный портрет 18 - числовая функция от матрицы критерий дихотомии спектрапрямой

19 Одномерный радиальный спектральный портрет критерий дихотомии спектраокружностью 19

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

21 1) исследование «устойчивости» решений дифференциальных уравнений Вопрос: для всех ли решений справедливо утверждение Критерий устойчивости: ? Универсальная оценка Н -- матрица Ляпунова – решения матричного уравнения Исследование устойчивости (по Ляпунову) 21 Дихотомия спектрапрямой имеет место, если существует матрица Грина при

22 Дихотомия прямой 22 Дихотомия спектрапрямой имеет место, если существует матрица Грина при Критерий дихотомии

23 Сходится ли итерационный процесс к решению системы ? Критерий сходимости: Н -- матрица решения дискретного матричного уравнения Ляпунова 23

24 Алгоритм анализа радиальной дихотомии спектра находимиз систем: после чего вычисляем Если то 24

25 Если то 25 - критерий дихотомии спектраокружностью Дискретное уравнение Ляпунова (обобщение): Оценки:

26 Одномерный радиальный спектральный портрет Мы показали как рассчитать и следовательно как нарисовать этот спектральный портрет Портреты дихотомии прямыми рассчитываются аналогично

27 APPLICATION OF NEW MATHEMATICAL TOOL ONE-DIMENTIONAL SPECTRAL PORTRAITS OF MATRIX TO THE PROBLEM OF AEROELASTICITY VIBRATION Godunov S.KNovosibirsk Kurzin V.B.Novosibirsk Bunkov V.G.Jukovskii Sadkane M.Brest (France) Из доклада, прочитанного на конференции по аэроупругости (Москва, октябрь 2006) 27

28 The simple flatter model Without the aerodynamic effect: Modeling of aerodynamic effects (v is the flow velocity) 28

29 29

30

31 The same example V 31

32 Упорядоченная последовательность букв Рассмотрим 6 букв алфавита:а б и п р т Рассмотрим большую (периодическую) последовательность букв: …ритатипбратарбатпиратритатипбратарбат… В этой последовательности: за буквой а следует 1 раз за период буква р и 4 раза буква т, за буквой б 1 раз следует буква а и 1 раз буква р… Таблица вероятностей следования букв

33 Этой последовательности соответствует матрица: Можно рассмотреть 32 буквы алфавита и любые длинные тексты, написанные с их помощью. Например, произведения разных писателей. Каждому произведению аналогичным способом сопоставляется 32х32 матрица. Можно ли идентифицировать писателя по спектральному портрету такой матрицы ? Упорядоченная последовательность букв

34 Характерные двумерные спектральные портреты писателей

35 Характерный одномерный спектральный портрет Б. Акунина

36 Характерные одномерные спектральные портреты писателей

37 Нормы проекторов ||P|| для спектральных точек одномерных спектральных портретов писателей

38 Литература

39 1.Уилкинсон Д. Х. Алгебраическая проблема собственных значений. М.: Наука, с. 2.Trefethen L.N. Psendospectra of matrices. D.F.Griffiths and G.A.Wotson, editors Numerical Analysis, P Longman Scientific and Technical. 3.Kostin V.I. Rassakov. On convergence of the power orthogonal method of spectrum Computing. Trans. Inst. Math. Sib. Branch Acad. Sci 6, P Годунов С.К., Булгаков А.. Круговая дихотомия матричного спектра. Сиб. Мат. журнал, Т. 29, 5. С Годунов С.К., Кирилюк О., Костин В. Спектральные портреты матриц. (Препринт Ин-та математики СО АН СССР, 3 ), с. 6.Годунов С.К.. Гарантированная точность в несимметричных задачах линейной алгебры. Дополнение к книге: Малышев А.Н. Введение в вычислительную линейную алгебру. Новосибирск, С Godunov S. Spectral portraits of matrices and criteria of spectrum dichotomy Computer Arithmetic and Enclosure Methods proc. 3-rd Intern. Sympos Oldenburg. Amsterdam P Godunov S.K., Sadkane M. Elliptic dichotomy of a matrix spectrum. Linear Algebra and its Appl Vol P Годунов С.К. Современные аспекты линейной алгебры. Научная книга. Новосибирск, Английский перевод Providence A.M.S p. (Transl. Math. Monogr. Vol. 169). 10.Годунов С.К. Лекции по современным аспектам линейной алгебры. Научная книга. Новосибирск, с. 11.Godunov S.K., Kursin V.B., Bunkov V.G., Sadkane M. Application of a new Mathematical tool «One- Dimensionsl spectral portraits of matrices» to the Problem of aeroelasticity vibrations of turbine blade cascades. Turbomachines: Aeroelasticity, aeroacoustics, and unsteady aerodynamics. Aeroelastic analysis of Bladerows, Moscow. TORUS PRESS Ltd., C Годунов С.К., Курзин В.Б., Буньков В.Г., Sadkane M. Применение нового математического аппарата «одномерные спектральные портреты матриц» к решению проблемы аэроупругих колебаний решеток лопастей // Ученые записки ЦАГИ, Т.XL, 6. С

40 Спасибо за внимание !