Главная » Просмотр файлов » Богачев - Практикум на ЭВМ. Методы решения линейных систем и нахождения собственных значений

Богачев - Практикум на ЭВМ. Методы решения линейных систем и нахождения собственных значений (947493), страница 17

Файл №947493 Богачев - Практикум на ЭВМ. Методы решения линейных систем и нахождения собственных значений (Богачев - Практикум на ЭВМ) 17 страницаБогачев - Практикум на ЭВМ. Методы решения линейных систем и нахождения собственных значений (947493) страница 172013-09-15СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Компоненты 1+1,...,и /с-го столбца матрицы А1й), равные компонентам вектора (/а1 )/е,, уже вычислены в (10). Столбец й вычисляется не по '1й — Ц 1ч — й) общим формулам (7) для сокращения количества арифметических операций и уменыпения вычислительной погрептности. 3. Поскольку в формуле (19) (а1 ); й+1 „й+1 „— — П(хй)(а," ); й+, „й+, „П(хй), (20) (й) . (й-Ц К.Ю.Богачев Точные методы решения линейных систем "315. ПРИВЕДЕНИЕ МАТРИЦЫ К ПОЧТИ ТРЕУГОЛЬНОМУ ВИДУ 77 ~15.

ПРИВЕДЕНИЕ МАТРИЦЫ К ПОЧТИ ТРЕУГОЛЬНОМУ ВИДУ 78 то в силу леммы 1 на вычисление подматрицы (20) матрицы Айй требуется 2(п — й)г + 0(и — 1с) мультипликативных и столько же адцитивных операций. Итак, на й-ом шаге алгоритма требуется выполнить п — й+ 1+ 2(п — й)г+ 0(и — й) = 2(п — й)г+0(п — й) мультипликативных операций, п — И+1+2(п — 1с)а+ 0(п — й) = 21п — й) + 0(п — й) аддитивных операций и 2 операции извлечения корня.

Следовательно, всего для проведения алгоритма требуется выполнить ~~;(2(п — й)~+0(н — й)) = 2ип — 1)(п — 2)(2п — 3)/6)+0(п ) = ~~й+0(пг) (и — + оо) й=1 мультипликативных операций, столько же аддитивных операций и 2(п — 2) операций извлечения корня (которые по трудоемкости по порядку можно сравнить с операциями деления). Таким образом, на приведение матрицы к почти треугольному виду унитарным подобием методом отражений требуется гпг + 0(нг) (п — > оо) мультипликативных операций и столько же аддитивных операций.

Заметим, что зто количество операций в два с половиной раза меньше, чем требуется для приведения произвольной матрицы к почти треугольному виду унитарным подобием методом отражений и совпадает количеством операций, необходимым для решения линейной системы методом отражений. Теорема 2. Всякал невырожденнол самосопряженнал матрица А может быть представлена в виде А = Я Я Я~, где матрица Я вЂ” рнитарнол, а матрица  — трехдиагональная. Доказательство Совпадает с доказательством теоремы 1. Хранение матриц Ч и В в памяти осуществляется одним из способов, изложенных при обсуждении алгоритма построения ЯЛ-разложения для матрицы А методом отражений.

Трудоемкость алгоритма построения описанного вьппе разложения складывается из количества арифметических операций, необходимых для проведения самого алгоритма, и количества арифметических операций, необходимых для построения матрицы Я. Подробные выкладки были проведены при обсуждении алгоритма построения ЯВ-разложения методом отражений. К.Ю.Богачев Точные методы решения линейных систем Глава 11. МЕТОДЫ НАХОЖДЕНИЯ СОБСТВЕННЫХ ЗНАЧЕНИЙ ~ 1. ТОт1НЫ1' И ИТЕРАЦИОННЫЕ М1' ТОДЫ Определение. Метод решения линейной системы называется точным, если при отсутствии округлений точное решение системы находится этим методом за конечное число арифметических операций (например, для метода Гаусса это 2пз + Р(пг) ) На реальной вычислительной магпине точный метод дает некоторое приближение к точному решению системы. Мера близости оценена в з 1.3.

Все описанные выше методы являются точными. Определение. Метод решения линейной системы называется итерационным, если он состоит в вычислении последовательности (хь), сходящейся к точному решению: хь — ~ х при й -+ сс. Итерационный метод за конечное число арифметических операций дает только некоторое приближение хь, к точному решению. Теория итерационных методов будет изложена в курсе "Численные методы". Определение. Метод нахождения собственных значений называется итерационным, если он состоит в вычислении последовательности (Ль), сходящейся к точному собственному значению: Ль — > Л при й -+ со.

Итерационный метод за конечное число арифметических операций дает только некоторое приближение Ль, к точному собственному значению. Теорема 1. (Без доказательства.) Не может существовать точного метода нахожденил всех собственных значений произвольной матрицы А Е М„при и ) 5. Другими словами, за конечное число арифметических операций нельзя найти все собственные значения произвольной матрицы А Е М„при и ) 5. К.Ю.Богачев Методы нахождения собственных значений ~2. ЛОКАЛИЗАЦИЯ СОБСТВЕННЫХ ЗНАЧЕНИЙ 80 ~ 2. ЛОКАЛИЗАЦИЯ СОБСТВЕННЫХ ЗНАЧЕНИЙ Для матрицы А Е М„обозначим В;'(А) = ~ ~а; ~, т'=1 хн ст(А) = 1Ль..., Л„) — множество собственных значений. Теорема 1 (Гершгорин).

Для всякой матрицы А Е М„справедливо соотношениет о(А) с Ц 1г е С": )г — аа( < В';(А)) . (1) Кроме тпого, если обьединение й, 1 < й < и из зтих кругов есть связная область, не пересекающаяся с остальными п — 1с кругами, то в ней находится ровно й собственнтях значений матрицтя А.

Доказательство. Пусть Л вЂ” собственное значение матрицы А и х (хм...,х„) ф Π— соответствующий собственный вектор: Ах = Лх. Обозначим ~хр~ = тпах; т „~х;~ ф О. Так как Ах = Лх, то Лхр — — (Лх)р — — (Ах)р — — ~ , 'а хт хр(Л вЂ” арр) = ~ ар х. 1=1 1Фр Следовательно, ~хрОЛ вЂ” арр~ < т ~арДхД < ~хр~;т ~арт~ у=т Ыр 1Фр Таким образом, собственное значение Л е о(А) принадлежит объединению кругов в (1).

Докажем второе утверждение теоремы. Положим А = П + В, где П = йаи(ан,...,а„„) Е М„, В = А — П; А, = П+ еВ, е > О. Тогда Н;'(А,) = Я';(еВ) = еЛ,'(А). Без ограничения общности мы можем считать что первые Ж К.Ю.Богачев Методы нахождения собственных значений ~2. ЛОКАЛИЗАЦИЯ СОБСТВЕННЫХ ЗНАЧЕНИЙ 81 кругов в (1) образуют связную область, не пересекаюшуюся с остальными и — й кругами.

Обозначим Сь(А) = [ ) «э е С": !я — ан! < Н,'(А)), и=1 С„ь(А) = Д 1Я Е С": !г — аи! < Л,'(А)), причем по условию Сь(А) П С д(А) = И. Рассмотрим введенные множества для матрицы А,: (2) Сь(А,) = Д (г Е С": !э — ап! < е~к(А)), в=1 С„ь(А,) = О (г Е С": !я — ан! < еВ!(А)).

Тогда для всех е Е [0,1] С,(А,) с С„(А,) = С,(А), С„,(А.) с С„ „(А,) = С„,(А). (8) Л;(А,) ~ Сь(Ая) [)Сч ь(А.). (4) По доказанному выгпе Сь(А,) П С„ь(А,) = 9 (б) для всех е е [О, 1[, причем Л;(А~) = Л;(А) = Л,, Л;(Ае) = аи. Следовательно, Л;(Ае) = аи е Сд(Ае) для 1 < з < й Л;(Ае) = ан Е С„ь(Ае) для й+ 1 < з < п.

(б) Методы нахождения собственных значений К.Ю.Богачев причем множество С~(А,) = Сь(А) связно по условию. В силу (2), (3) Сь(А,) П С ь(А,) = 9 для всех е е [0,1!. Обозначим через Л;(А,) з-е собственное значение матрицы А„ н(А,) = 1Л1(А,),..., Л„(А,)) - спектр матрицы А,.

По построению Л;(А1) = Л;(А) = Л;, Л;(Ае) = ан. Собственное значение Л;(А,) матрицы А, является корнем характеристического многочлена этой матрицы. Корни многочлена являются непрерывными функциями его коэффициентов. Следовательно, собственные значения Л;(А,) матрицы А, являются непрерывными функциями элементов этой матрицы, которые, в свою очередь, являются непрерывными фунциями параметра е. Таким образом, как композиции непрерывных функций, собственные значения Л;(А,) матрицы А, являются непрерывными фунциями параметра е. По доказанному первому утверждению теоремы н(А,) С Сь(А,) () С„ь(А,) для всех е Е [О, 1[, т.е.

~3. ОШИБКИ ПРИ НАХОЖДЕНИИ СОБСТВЕННЫХ ЗНАЧЕНИЙ 82 Поскольку Лс(А,) — непрерывная функция и принимает все промежуточные зна- чения, то из (4), (5), (6) вытекает, что Л;(А,) Е Сь(А,) для 1 < с < й Л;(А,) = аа Е С„с»(А,) для й+ 1 < з < н для всех е Е (О, Ц. При е = 1 мы получаем второе утвеждение теоремы. 8 3.

ОШИБКИ ПРИ НАХО:~КДЕНИИ СОБСТВЕННЫХ ЗНАЧЕНИЙ Пусть находятся собственные значения матрицы А Е М„Пусть 6 алгоритм нахождения собственных значений, д(А) = 6(А) — полученный этим алгоритмом спектр матрицы А. Этот спектр не совпадает с истинным значением о(А) спектра матрицы А, поскольку в силу теоремы 1.1 найти все собственные значения матрицы А за конечное число арифметических операций невозможно. Обозначим через А + Е матрицу, спектром которой является набор д(А), т.е. д(А) = о.(А+ Е).

Оценим погрешность при нахождении собственных значений (т.е. д(А) — сг(А) ) через величину матрицы .Е. Лемма 1. Пусть Р = йаи(Л»,..., Л„) — диагональная матрица, Е = (ес ) е М„. Тогда длл всякого Л Е о (Р+ Е) собственного значения магирицы .Р+ Е, существует Л = Л; Е о (Р) — собственное значение лсатрицьс .Р, такое, что )Л вЂ” Л! < )/Е!) (где )) (! — максимальнол строчнал норма матрицьс, см. стр.

9). Доказательство. Применяя теорему Гершгорина к матрице Р+Е, находим, что существует такое с, 1 < с < н, что Л ~ г ~ С": (г — Л; — еа~ < Б,'(Е) = с ')ес,( с=» сФ» Так как (г — Л; — ен( > )г — Лс( — (есс), то Л е г е С": )г — Лс( < ~ (ес) 1=» Но ~~, )ес,! < шах ~~, (ес ( = !)Е(! Поэтому Л Е «г Е С": )г — Лс! < )/Е!)„) что доказывает утверждение леммы. К.Ю.Богачев Методы нахождения собственных значений ~4.

СТЕПЕННОЙ МЕТОД Теорема 1. Пусть А Е М„диагонализируемая матрица, тп.е. су»цестпвуют невь»розтсденная матрица С и диагональная матрица Л = »1»афЛ»,..., Л„) такие, что А = СЛС '. Пустпь Е = (е; ) Е М„. Тогда для всякого Л Е о(А+ Е) — собственного значения матриць» А+ Е, су»цестпвует Л = Л; е о(А) собственное значение матрицы А, п»акое, что (Л вЂ” Л( < к (С)()Е)! где к (С) — число обусловленности матрицы С относитпельно максимальной строчной нормы матрицы (см. стпр. 13).

Доказательство. Матрицы А+Е и С '(А+Е)С = Л+С 'ЕС подобны и потому имеют одинаковые собственные значения. В силу леммы 1 для всякого Л е в(А + Е) = »т(Л + С 'ЕС) существует Л = Л; е о(Л) = о(А) такое, что )Л вЂ” Л) < ()С 'ЕС!) < !)С ')( )(Е!) ))Сй = к (С))(Е!) Теорема доказана.

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

Тип файла
DJVU-файл
Размер
497,22 Kb
Тип материала
Учебное заведение
Неизвестно

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

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