Главная » Просмотр файлов » 1611688890-f641c9ec8276824e4686da772eb56520

1611688890-f641c9ec8276824e4686da772eb56520 (826652), страница 48

Файл №826652 1611688890-f641c9ec8276824e4686da772eb56520 (Шарый Курс вычислительных методов) 48 страница1611688890-f641c9ec8276824e4686da772eb56520 (826652) страница 482021-01-26СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

Прямые методы решения линейных системкоторая отличается от единичной матрицы наличием одного дополнительного ненулевого элемента ri1 на месте (i, 1). Исключение поддиагональных элементов первого столбца матрицы СЛАУ — это последовательное домножение обеих частей этой системы слева на матрицы1 r 21 0 . . .001..0.1и так далее до11 0 r31 . . .,010 0 . .. 0rn1.11..,0100..01.11.Нетрудно убедиться, что умножение матриц выписанного выше ви-2923.

Численные методы линейной алгебрыда выполняется по простому правилу: ri11 01...1..0.1     ·     rk1 ri1= rk111....010...1....10101.1.Оно также остаётся верным в случае, когда у матриц-сомножителейна несовпадающих местах в первом столбце присутствует более одного ненулевого элемента. Следовательно, обнуление поддиагональныхэлементов первого столбца и соответствующие преобразования правойчасти в методе Гаусса — это не что иное, как умножение обеих частейСЛАУ слева на матрицу1 r21E1 =  r31 . . .rn101010...1.(3.62)Аналогично, обнуление поддиагональных элементов j-го столбцаматрицы СЛАУ и соответствующие преобразования правой части мож-2933.6.

Прямые методы решения линейных системно интерпретировать как умножение системы слева на матрицуEj = 1..0.10rj+1,j...1...rnj1.(3.63)В целом метод Гаусса представляется как последовательность умножений обеих частей решаемой СЛАУ слева на матрицы Ej вида (3.63),j = 1, 2, . .

. , n − 1. При этом матрицей системы становится матрица(3.64)En−1 · · · E2 E1 A = U,которая является верхней треугольной матрицей.Коль скоро все Ej — нижние треугольные матрицы, их произведение также является нижним треугольным. Кроме того, все Ej неособенны (нижние треугольные с единицами по главной диагонали). Поэтомунеособенно и их произведение En−1 · · · E2 E1 .

Если определитьL = En−1 · · · E2 E1−1,то, как нетрудно понять, L — тоже нижняя треугольная матрица с единицами по главной диагонали. Для этой матрицы в силу (3.64) справедливо равенствоA = LU.Получается, что исходная матрица СЛАУ оказалась представленной в виде произведения нижней треугольной L и верхней треугольной U матриц.

Это представление называют треугольным разложением матрицы или LU-разложением.14 Соответственно, преобразованияматрицы A в прямом ходе метода Гаусса (3.59) можно трактовать какеё разложение на нижний треугольный L и верхний треугольный Uмножители.14 От английских слов lower (нижний) и upper (верхний). Нередко для обозначения этого же понятия можно встретить кальки с иностранных терминов — «LUфакторизация» и «LU-декомпозиция».2943. Численные методы линейной алгебрыОтметим, что если LU-разложение матрицы A исходной СЛАУ ужедано, то система Ax = b может быть переписана в равносильной формеL (U x) = b.Тогда её решение сводится к решению двух треугольных систем линейных алгебраических уравнений(Ly = b,(3.65)Ux = yс помощью прямой и обратной подстановок соответственно. В действительности LU-разложение, получаемое с помощью версии метода Гаусса (3.59)–(3.60), таково, что в нижней треугольной матрице L по диагонали стоят все единицы.

Тогда при реализации метода Гаусса на компьютере для экономии машинной памяти можно хранить треугольныесомножители L и U на месте A, так как диагональ L имеет фиксированный вид.3.6гМетод Гаусса с выборомведущего элементаИ в прямом, и в обратном ходе метода Гаусса встречаются операции деления, которые не выполнимы в случае, когда делитель равеннулю. Тогда не может быть выполнен и метод Гаусса в целом. Этотраздел посвящен тому, как модифицировать метод Гаусса, чтобы онбыл применим для решения любых СЛАУ с неособенными матрицами.Ведущим элементом в методе Гаусса называют элемент матрицырешаемой системы, на который выполняется деление при исключенииподдиагональных элементов очередного столбца.15 В алгоритме, описанном в §3.6б, ведущим всюду берётся фиксированный диагональныйэлемент ajj , вне зависимости от его значения, но желательно модифицировать метод Гаусса так, чтобы ведущий элемент, по возможности,всегда был отличен от нуля.

С другой стороны, при решении конкретных СЛАУ, даже в случае ajj 6= 0, по соображениям устойчивости алгоритма более предпочтительным иногда может оказаться выбор другогоэлемента в качестве ведущего.15 Иногдав русской математической литературе его назыают главным элементом.3.6. Прямые методы решения линейных систем295Отметим, что любое изменение порядка уравнений в системе приводит к равносильной системе уравнений, хотя при этом в матрицеСЛАУ переставляются строки и она существенно меняется.

Воспользуемся этим наблюдением для организации успешного выполнения метода Гаусса.Назовём активной подматрицей j-го шага прямого хода методаГаусса подматрицу исходной матрицы СЛАУ, образованную последними n−j+1 строками и столбцами. Именно эта подматрица подвергаетсяпреобразованиям на j-ом шаге прямого хода, тогда как первые j − 1строк и столбцов остаются уже неизменными.(j −1j −1(0активнаяподматрица0j-ый столбецРис.

3.13. Структура матрицы СЛАУ перед началомj-го шага прямого хода метода Гаусса.Частичным выбором ведущего элемента на j-ом шаге прямого хода метода Гаусса называют его выбор, как максимального по модулюэлемента из всех элементов j-го столбца, лежащих не выше диагонали.Соответственно, частичный выбор ведущего элемента сопровождаетсянеобходимой перестановкой строк матрицы и компонент правой части(т. е. уравнений СЛАУ), когда этот максимальный по модулю элементстановится диагональным.

Следует пояснить, что именно максимальным по модулю, а не просто ненулевым, ведущий элемент выбираетсядля того, чтобы обеспечить наибольшую численную устойчивость алгоритма в условиях вычислений с конечной точностью.2963. Численные методы линейной алгебрыПредложение 3.6.1 Метод Гаусса с частичным выбором ведущегоэлемента всегда выполним для систем линейных алгебраических уравнений с неособенными квадратными матрицами.Доказательство. Преобразования прямого хода метода Гаусса сохраняют свойство определителя матрицы системы быть неравным нулю.Перед началом j-го шага метода Гаусса эта матрица имеет блочнотреугольный вид, изображённый на Рис.

3.13, и поэтому её определитель равен произведению определителей ведущей подматрицы порядка(j −1) и активной подматрицы порядка n−j +1. Следовательно, активная подматрица имеет ненулевой определитель, т. е. в первом её столбцеобязан найтись хотя бы один ненулевой элемент. Максимальный по модулю из этих ненулевых элементов также ненулевой, и его мы делаемведущим. Как следствие, прямой ход метода Гаусса выполним.Обратный ход также не встречает деления на нуль, поскольку полученная в прямом ходе верхняя треугольная матрица неособенна, т. е.все её диагональные элементы должны быть ненулевыми.Каково матричное представление метода Гаусса с выбором ведущегоэлемента? Введём элементарные матрицы перестановок1...0···1← i-ая строка1......P =(3.66)...1← j-ая строка1···0...1(называемые также матрицами транспозиции), которые получаютсяиз единичной матрицы перестановкой двух её строк (или столбцов) сномерами i и j.

Умножение на такую матрицу слева приводит к перестановке i-ой и j-ой строк, а умножение справа — к перестановке i-го иj-го столбцов. Тогда для метода Гаусса с частичным выбором ведущегоэлемента справедливо следующее матричное представление(En−1 Pn−1 ) · · · (E1 P1 )A = U,3.6. Прямые методы решения линейных систем297где Ej — матрицы преобразований вида (3.63), введённые в предыдущем разделе, а P1 , P2 , .

. . , Pn−1 — элементарные матрицы перестановок(3.66), при помощи которых выполняется необходимая перестановкастрок на 1-м, 2-м, . . . , (n − 1)-м шагах прямого хода метода Гаусса.Несмотря на то, что метод Гаусса с частичным выбором ведущего элемента теоретически всегда спасает положение, на практике длянекоторых «плохих» СЛАУ он всё-таки может работать недостаточноустойчиво. Это происходит в случае, когда ведущий элемент оказывается малым, и потому коэффициенты rij из прямого метода Гаусса(3.59) получаются большими по абсолютной величине. По этой причинедля обеспечения желаемой устойчивости вычислительного процесса пометоду Гаусса иногда имеет смысл выбирать ведущий элемент болеетщательно, чем это делается при описанном выше частичном выборе.Заметим, что ещё одним простым способом равносильного преобразования системы уравнений является перенумерация переменных.

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

Метод Гаусса с полным выбором выдущегоэлемента имеет следующее матричное представление(En−1 P̌n−1 ) · · · (E1 P̌1 )AP̂1 · · · P̂n−1 = U,где P̌i — элементарные матрицы перестановок, при помощи которыхвыполняется перестановка строк, P̂j — элементарные матрицы перестановок, с помощью которых выполняется перестановка столбцов насоответствующих шагах прямого хода метода Гаусса.Напомним, что матрицей перестановок называется матрица, получающаяся из единичной матрицы перестановкой произвольного числаеё строк (или столбцов). Матрица перестановок может быть представлена как произведение нескольких элементарных матриц перестановоквида (3.66) (см. подробности, к примеру, в [7]).Теорема 3.6.1 Для неособенной матрицы A существуют матрицыперестановок P̌ и P̂ , такие чтоP̌ AP̂ = LU,2983.

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

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

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

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