Главная » Просмотр файлов » Дж. Деммель - Вычислительная линейная алгебра

Дж. Деммель - Вычислительная линейная алгебра (1156793), страница 43

Файл №1156793 Дж. Деммель - Вычислительная линейная алгебра (Дж. Деммель - Вычислительная линейная алгебра) 43 страницаДж. Деммель - Вычислительная линейная алгебра (1156793) страница 432019-09-18СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 43)

Одним из таких множеств является малая окрестность матрицы О 1 О О 1 О Ь О О вЂ” 5 О О О 1 О Число Ь получается умножением машинного эпсилона на константу, равную нескольким тысячам. Чтобы устранить подобные случаи несходимости, в современные реализации алгоритма был добавлен еще один «особый сдвиг». Однако построение стратегии сдвигов, которая гарантировала бы быструю сходимость для всех матриц, остается все еще открытой проблемой. 4.5, Друтие типы несимметричных спектральных задач 4.5.1. Регулярные матричные пучки н каноническая форма Вейерштр асса В стандартной проблеме собственных значений разыскиваются числа и, для которых матрица А — г1 вырожденна; эти числа называются собственными значениями.

Существует несколько важных обобщений этого понятия. Определение 4.6. Пусть А и  — матрицы размера т х и. Матрица А — ЛВ называется матричным пучком, или просто пучком. Здесь Л вЂ” параметр, а не конкретное число. Определение 4Л. Если А и  — квадратные матрицы и бег(А — ЛВ) не равен нулю тождественно, то пучок А — ЛВ называют регулярным. В противном случае, пучок н зывагтся сингулярным. Для регулярного пучка А — ЛВ многочлен р(Л) = «)е»(А — ЛВ) называется характеристическим многочленом, а собственныг значения пучка А — ЛВ определяются как: 186 Глава 4. Несимметричная проблема собственных значений (1) корни многочлена р(Л) = О и, если г1еб(р) < п, (2) сю (кратности, равной п — Лей(р)). Пример 4.12. Пусть А — ЛВ= 1 — Л О Тогда р(Л) = г(ег(А — ЛВ) = (1 — 2Л) (1 — ОЛ) . (Π— Л) = (2Л вЂ” 1)Л, поэтому собственными значениями являются —, О и со.

1 О Матричные пучки возникают естественным образом во многих математических моделях физических систем; ниже мы приведем соответствующие примеры. В следующем предложении указывается связь между собственными значениями регулярного пучка А — ЛВ и собственными значениями некоторой матрицы. Предложение 4.6.

Пусть пучок А — ЛВ регуллрен. Если матрица В невырожденна, то все собственные значения пучка конечны и совпадают с собственными значениами матрицы АВ г или матприцы В гА. Если  — вырожденная матрица, то А — ЛВ имеет собственное значение со кратности п — гапЦВ). Если матприца А невырожденна, то собственные значения пучка А — ЛВ суть числа, обратные к собственным значениям матрицы А 'В или матрицы ВА г. При этпом нулевому собстпвенному значению матприцы А гВ соответствует бесконечное собственное значение пучка А — ЛВ. Доказательство. Если Л' — собственное значение пучка и матрица В невырожденна, то О = г(ег(А — Л'В) = г1еь(АВ ' — Л'1) = г1ес(В ~А — Л'1), поэтому Л' является собственным значением и для матриц АВ ' и В 'А. Для вырожденной В положим р(Л) = г(ег(А — ЛВ) и подставим сюда сингулярное разложение В = УЕИт матрицы В.

Имеем р(Л) = деГ(А — ЛПЕ(тт) = бег(П(ПтАР— ЛЕ)Ит) = ~г(ег(ПтАР' — ЛЕ) Поскольку гапк(В) = гапк(Е), то на диагонали матрицы ПтА)г — ЛЕ параметр Л появляется только гела(В) раз. Поэтому степень многочлена г1е1(ПтАИ вЂ” ЛЕ) равна гапк(В). Если А — невырожденная матрица, то равенство г(еь(А — ЛВ) = О эквивалентно равенству г1еь(1 — ЛА тВ) = О или равенству т1еС(1 — ЛВА г) = О. Эти последние равенства могут выполняться лишь, если Л ,-Е О и -„' есть собственное значение матриц А гВ и ВА '.

П Определение 4.8. Пусть Л' — конечное собственное значение регулярного пучка А — ЛВ. Ненулевой вектор х, такой, что (А — Л'В)х = О или, что эквивалентно, Ах = Л'Вх, называется правым собственным вектором. Длл собственного значения Л' = оо правым собственным вектором называетсл вектпор х ф О, тпакой, чтпо Вх = О. Левый собственный вектор пучка А — ЛВ определяется как правый собственный вектор пучка А* — ЛВ'. Пример 4.13. Рассмотрим пучок А — ЛВ из примера 4.12. Так как А и В— диагональные матрицы, то правые и левые собственные векторы — это просто столбцы единичной матрицы. О 187 4.5. Другие тяпы несимметричных спектральных зааач Пример 4.14. Рассмотрим связанную систему материальных точек из примера 4.1. С этой системой естественным образом связаны два матричных пучка.

Во-первых, задачу на собственные значения — М ' — М 'К можно записать в виде =Л Эта постановка задачи может быть более удачной, если матрица М очень плохо обусловлена, поскольку в этом случае трудно вычислить хорошие приближения к матрицам М 'В и М 'К. Во-вторых, часто рассматривают случай, когда В = 0 (нет затухания). Исходное дифференциальное уравнение в этом случае превращается в Мх(г) + Кх(1) = О. Если искать его решения вида х,(1) = е~*'»х,(0), то получим Лзе""Мх;(0) + е~нКхг(0) = О, или ЛгМх,(0) + Кх,(0) = О. ДРУгими словамн, — Лз есть собственное значение, а х;(0) — правый собственный вектор пучка К вЂ” ЛМ.

Поскольку, по предположению, матрица М невырожденна, — Лз и х,(0) являются собственным значением и собственным вектором и для матрицы М !К. О Бесконечные собственные значения тоже возникают в практических задачах вполне естественным образом. Например, ниже в этом разделе мы покажем, как бесконечные собственные значения связаны с импульсным откликом системы, описываемой обыкновенными дифференциальными уравнениями с линейными ограничениями, иначе говоря, дифференции«ьно-алгебраическими уравнениями (41].

См. также вопрос 4.16, где обсуждаются приложения матричных пучков к вычислительной геометрии и машинной графике. Вспомним, что и теория, и алгоритмы для матричной проблемы собственных значений предусматривали приведение исходной матрицы А с помощью подобия 5 'АЯ к некоторой форме, более «простой», чем А.

Приводимое ниже определение указывает, как обобщить понятие подобия на случай матричных пучков. После этого будут обобщены для пучков и формы Жордана и Шура. Определение 4.9. Пусть Рь и Рл — произвольные невирожденные матри- цы. Роворят, что пучки А — ЛВ и РсАРя — ЛРьВРя эквивалентны. Предложение 4.7. Эквивалснтнмс регулярные пучки А — ЛВ и Р»АРя— ЛРьВРя имеюта одни и те лсе собственные значенил.

Вектор х тогда и только тогда является правым собственным вектором пучка А — ЛВ, когда Р 'х является правым собствснны.м вектором пучка РвАРя — ЛРьВРя. Вектор у тогда и только тогда является левым собственным вектором пучка А — ЛВ, когда (Р~) 1у является левым собственнькм вектором пучка РсАРя — ЛРьВРя.

Доказательство. Равенство с(ег(А — ЛВ) = 0 эквивалентно равенству г(ес(Рь(А — ЛВ)Ря) = О. Равенство (А — ЛВ)х = 0 эквивалентно равенству Рг. (А — ЛВ)РяР 'х = О. 188 Глава 4. Несимметричная проблема собственных значений Равенство (А — ЛВ)'у=О эквивалентно равенству Рл(А — ЛВ)'Р;(Р,*) 'у=О.ГС В следующей теореме каноническая форма Жордана обобщается на случай регулярных матричных пучков. Теорема 4.10 (каноническая форма Вейерштрасса). Для регулярного пучка А — ЛВ найдутся невьсрожденные матрицы Рс. и Рн, такие, что Рь(А — ЛВ)Рй = гйа8(дп,(Лс) Л1пс1 ° Упс(Лп,) Л1п, «п«««п, Дгп«) где У„,.

(Л,) — жорданов блок порядка пс с собственным значением Л,: Л, 1 .Упс(Лс) = 1 Л, а Усспм — гжорданов блок для собственного значения оо кра«пности тс»м 1 Л = 1, — Л,У,. (0). Л 1 Доказательство теоремы можно найти в [110]. Приложения форм Жордаиа и Вейерштрасса к дифференциальным уравнениям Рассмотрим задачу Коши для линейного дифференциального уравнения х(г) = Ах(1) + У(1), х(0) = хо. Ке решение можно выписать в явном виде: х(1) = елсхо + у е"!с «>у(т)дт. Коли мы умеем приводить А к форме Жордана: с лс, А = ЯУЯ ', то, производя в дифференциальном уравнении замену переменного у(с) = Я 'х(1), получим у(с) = Уу(г) + Б с)(С), у(0) = уо = Я 'хо. Решение этой задачи имеет вид у(с) = емуо + ) ед' '>$ сУ(т)сУт.

Имеется явная форз с-« мула, позволяющая найти езс нли любую другую функцию от жордановой матрицы 1. (Эту формулу не следует использовать в машинных вычислениях! В вопросе 4.4 обсуждается более надежный подход к вычислению функций от матрицы.) Предположим, что функция У задана своим рядом Тейлора У(г) = ~ с о —,.; —, а .У вЂ” жорданов блок:,У = Л1+ М; в матрице М первая нэддиагональ состоит из единиц, а прочие элементы равны нулю. Имеем 10> (О) (ЛУ+ 1Ч) с с! с=о У<сс(0) с«с — г (')г-гяс ф««г и с! с=о 1=0 189 4.5. Другие тилы несимметричных спектральных задач /О)(0) /4 «.)»' 'нс Рч«у г! У=е ~=1 = ~ а* ~ Г',"' (*) у=с Л 1 пхн " 'Л /ОО(Л) З з=е /(1) =/ 1 Л у(Л) ур(Л) У" (л): %++л11 (4.6) /'(Л) /(Л) Итак, 1(з) — верхняя треугольная матрица, на у-й налдиагонали которой стоит число /®(Л)/у!.

Для решения более общей задачи Вх = Ах + /(1), где пучок А — ЛВ регулярен, используем форму Вейерштрасса. Предположим, что пучок Рь(А— ЛВ)Рн находится в форме Вейерштрасса, и перепишем уравнение в виде РьВРнР,'х = РьАРнР,'х + Рь/(1). Положим Рл'х = у и Рь/Я = д(1). Теперь задача разлагается в сумму подзадач: 1 „ В последнем равенстве использовано то обстоятельство, что № =О для у)п — 1. Заметим, что при у' < п матрица № содержит единицы на у-й налдиагонали и нули во всех прочих позициях.

Наконец, заметим, что ~, . —,, (,'./ Л' т есть разложение Тейлора для /®(Л)/Я. Таким образом, 100 Глава 4. Несимметричная проблема собственных значений Каждая подзадача у = .7„, (Л,)у+д(т) =,Уу+у(т) есть стандартное линейное дифференциальное уравнение, решение которого имеет вид ,г у(Ю) = у(0)ел+ / езр ~~д(т)дт, о Решение системы д (0)у = у + д(т) достигается обратной подстановкой, исходя из последнего уравнения. Запишем систему более подробно: 0 1 д1 у1 у1 ды Уы Из т-го (т.е. последнего) уравнения имеем 0 = у + д, или у = — д Уравнение с номером г' есть у,, = у;+д,, откуда у; = уг+, — д; и, как следствие, Следовательно, решение зависит от производных функции д, а не выражается интегралом от д, как для обычного дифференциального уравнения. Поэтому для непрерывной, но не дифференцируемой функции д мы можем получить разрывное решение. Такую ситуацию иногда характеризуют термином импульсный отклик.

Характеристики

Тип файла
PDF-файл
Размер
40,83 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
А знаете ли Вы, что из года в год задания практически не меняются? Математика, преподаваемая в учебных заведениях, никак не менялась минимум 30 лет. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6390
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее