ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 20.09.2025
Просмотров: 1192
Скачиваний: 0
СОДЕРЖАНИЕ
Л.В. Маркова, е.А. Корчевская,
Глава 1 Элементы теории погрешностей п 1.1 Источники погрешностей
П 1.2 Вычисление абсолютной и относительной погрешностей
П 1.4 Вычисление погрешностей арифметических операций
П 1.5 Оценка погрешности по способу границ
Глава 2 объектно-ориентированный подход к программированию методов линейной алгебры
П 2.1 Создание матричной иерархии классов
П 2.2 Создание иерархии классов вычислительных методов алгебры
Глава 3 решение систем линейных алгебраических уравнений
П 3.1 Метод Гаусса решения систем линейных алгебраических уравнений
П 3.2 Метод Гаусса с выбором главного элемента для решения систем линейных алгебраических уравнений
П 3.3 Решение системы линейных алгебраических уравнений методом Жордана-Гаусса
П 3.4 Метод квадратного корня для решения систем линейных алгебраических уравнений
П 3.5 Вычисления определителя и нахождения обратной матрицы
П 3.6 Решение системы линейных алгебраических уравнений методом прогонки
П 3.7 Метод простых итераций решения систем линейных алгебраических уравнений
П 3.8 Метод Зейделя решения систем линейных алгебраических уравнений
П 3.9 Итерационные методы вариационного типа решения систем линейных алгебраических уравнений
Глава 4 вычисление собственных значений и собственных векторов матриц
П 4.1 Метод Данилевского для нахождения собственных значений и собственных векторов
П 4.3 qr-алгоритм для нахождения собственных значений матрицы
П 4.4 Метод Якоби для нахождения собственных значений и собственных векторов
П р и л о ж е н и я Приложение 1 Основные сведения о матрицах
.
Расширенная матрица – это матрица, полученная «склеиванием» двух матриц, как правило, для осуществления одних и тех же элементарных преобразований со строками сразу в двух матрицах. Например, для матриц
,
расширенная матрица будет выглядеть
следующим образом:
.
Суммой матриц A и B называется матрица C такая, что элемент
.
Матрица
называетсяпроизведением
матриц A
и B,
если элемент
. (1)
Причем умножение
матриц не является коммутативным, т.е.
(за исключением умножения на единичную
матрицу:
).
Поэтому для
будет справедлива формула:
. (2)
Данные формулы (1) и (2) будут называться соответственно левостороннее и правостороннее произведение матриц.
Вектор – это матрица, состоящая из одного столбца (вектор- столбец), т.е.
.
.
При умножении матрицы A на вектор X мы получим некоторый вектор
,
элементы которого вычисляются по
формуле:
.
Квадратные матрицы А и В одинакового порядка называются подобными, если существует невырожденная матрица P такого же порядка, такая что
.
Матрицей, обратной данной матрице A, называется матрица
такая, что произведениеA
на
равняется единичной матрице
.
Ортогональной называют такую квадратную матрицу A, для которой выполняется равенство
.
Невырожденной называют квадратную матрицу А, определитель которой не равен 0.
Квадратная матрица называется симметрической, если
для любыхi,
j,
т.е. ее элементы расположены симметрично
относительно главной диагонали.
Транспонированная матрица – матрица
,
полученная из исходной матрицыА
заменой строк на столбцы:

.
Признаки положительной определенности матрицы
1. Критерий Сильвестра. Чтобы матрица А была положительно определенной, необходимо и достаточно, чтобы все ее главные миноры были положительны.
2. Достаточное
условие. Диагональное преобладание,
т.е. свойство
влечет
положительную определенность матрицы.
Нормой вектора Х называется поставленное в соответствие этому вектору неотрицательное число
,
удовлетворяющее аксиомам:
Положительная определенность, т.е. для любого ненулевого вектора его норма больше нуля и равна нулю только для ноль вектора

Однородность


Существует несколько способов введения нормы вектора. Наиболее употребительными являются следующие:
первая (кубическая)
;вторая (октаэдрическая)
;третья (сферическая)
.
Нормой матрицы А называется поставленное этой матрице в соответствие неотрицательное число
такое, что
;
;
;
.
Здесь Н – линейное пространство квадратных матриц n-го порядка.
Норма матрицы, как и норма вектора, может быть определена по-разному.
первая
;
вторая
;
,
где
наибольшее собственное значение матрицы
.
Если для любой матрицы А и любого вектора Х выполняется неравенство
то говорят, что норма матрицы согласована
с данной нормой вектора.
Нормой матрицы А, подчиненной данной норме вектора, называется число
верхняя грань (т.е. максимальное число)
множество норм такого вида.
Собственным значением квадратной матрицы А n-го порядка
называется такое
число λ,
при котором для некоторого ненулевого
вектора
имеет место
равенство
АХ= λХ.
Любой ненулевой вектор X, удовлетворяющий этому равенству, называется собственным вектором матрицы А, соответствующим собственному значению λ.
Матрица
называется характеристической матрицей
данной матрицыА.
Уравнение
называется
характеристическим уравнением матрицы
А,
а полином
характеристическим
полиномом.
Совокупность всех собственных значений λ1, λ2, …, λn матрицы А называется спектром этой матрицы.
Спектральным радиусом (А) матрицы Аназывается максимум из модулей собственных значений этой матрицы.
Собственные значения треугольной матрицы равны ее диагональным элементам.

.
Приложение 2
Функции MathCad
Для контроля вычислений можно использовать стандартные функции пакета MathCad:
eigenvals (A) – возвращает вектор из собственных значений матрицы A.
eigenvec (A, λ) – нормированный собственный вектор матрицы A, соответсвующий ее собственному значению λ.
eigenvecs (A) – возвращает матрицу, чьими столбцами являются собственные векторы матрицы A. Порядок расположения собственных векторов соответсвует порядку собственных значений, возвращаемых функцией eigenvals.
max(A), min(A) – вычисление максимального и минимального элементов матрицы или вектора.
norm1(A) – возвращает первую норму матрицы или вектора.
norm2(A) – возвращает вторую норму матрицы или вектора.
norme(A) – возвращает третью норму матрицы или вектора.
Для решения системы уравнений необходимо выполнить следующее:
• Задать начальное приближение для всех неизвестных, входящих в систему уравнений, так как Mathcad решает систему с помощью итерационных методов.
• Оформить блок решения системы с помощью ключевого слова Given. Оно указывает Mathcad, что далее следует система уравнений.
• Записать систему уравнений. Для печати символа = необходимо использовать комбинацию клавиш [Ctrl]+[ =] .
• Для решения системы вызвать функцию Find или Minerr.
Find(x1,x2,…,xn), Minerr(x1,x2,…,xn) – возвращает решение системы уравнений в виде вектора. В качестве аргументов указываются имена искомых неизвестных. Число аргументов должно быть равно числу неизвестных.
Функция Minerr использует тот же самый алгоритм, что и функция Find. Различие состоит в следующем: если в результате поиска решения не может быть получено дальнейшее уточнение текущего приближения к решению, Minerr возвращает это приближение, а функция Find возвращает в этом случае сообщение об ошибке “решение не найдено”.
Л и т е р а т у р а
Архангельский А.Я. С++ Builder 6: справоч. пособие. – М.: Бином-Пресс, 2002. Кн. 1. Язык С++. – 564 с.
Архангельский А.Я. С++ Builder 6: справоч. пособие. – М.: Бином-Пресс, 2002. Кн. 2. Классы и компоненты. – 528 с.
Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. М.: Бином, 2004. – 636 с.
Бахвалов Н.С. Численные методы в задачах и упражнениях: учеб. пособие. М.: Высшая школа (Высшая математика), 2000. 190 с.
Икрамов Х.Д. Нессиметричная проблема собственных значений. Численные методы. М.: Наука, 1991. – 240 с.
Исаков В.Н. Элементы численных методов. – М.: Изд. Центр «Академия», 2003.
Калиткин Н.Н. Численные методы: учеб. пособие. – М.: Наука, 1978. – 512 с.
Киреев В.И., Пантелеев А.В. Численные методы в примерах и задачах: учеб. пособие. М.: «Высшая школа», 2006. – 480 с.
Костомаров Д.П., Фаворский А.П. Программирование и численные методы. – М.: МГУ, 2004.
Крылов В.И., Бобков В.В., Монастырный П.И. Вычислительные методы высшей математики: учеб. пособие. – Минск: Вышэйшая школа, 1972. – 585 с.
Маркова Л.В., Адаменко Н.Д., Казанцева О.Г., Корчевская Е.А. Формирование профессиональных компетенций у студентов специальности «Прикладная математика» // Вестн. Витебск. гос. ун-та. 2012. № 1(67). – С. 116121.
Маркова Л.В., Корчевская Е.А. Обучение вычислительной математике. Современные аспекты // Инновационные технологии обучения физико-математическим дисциплинам: материалы междунар. науч.-практ. интернет-конф., посвященной 60-летию доктора физико-математических наук, профессора Н.Т. Воробьева, Витебск, 21–22 июня 2011 г. Витебск, 2011. – С. 128–129.
Маркова Л.В., Корчевская Е.А., Красоткина А.Н. Объектная реализация методов вычислительной алгебры // Вестн. Витебск. гос. ун-та. 2013. № 2(74). – С. 1822.
Марчук П.И. Методы вычислительной математики. М.: Наука, 1989. – 608 с.
Павловская Т.А. С/С++. Программирование на языке высокого уровня: учебник для вузов. СПб.: Питер, 2003. – 461 с.
Рыжиков Ю.И. Вычислительные методы. – СПб.: БХВ-Петербург, 2007.
Ращиков В.И., Рошаль А.С. Численные методы решения физических задач. – СПб.: «Лань», 2005.
Самарский А.А. Введение в численные методы. М.: Наука, 1983. – 272 с.
Самарский А.А., Гулин А.В. Численные методы. М.: Наука, 1989. – 432 с.
Сборник задач по методам вычислений: учеб. пособие для студ. учреждений, обеспечивающих получение высш. образования по физико-математическим спец. / под ред. П.И. Монастырного. Минск: Изд. центр БГУ, 2007. 376 с.
Семенов В.А., Тарлапан О.А. Объектно-ориентированный подход к программированию прямых методов линейной алгебры // Вопросы кибернетики. Приложения системного программирования / под ред. В.П. Иванникова. – М.: Науч. совет по комплексной проблеме «Кибернетика» РАН, 1996. Вып. 2. – С. 147170.
Турчак Л.И. Основы численных методов: учеб. пособие для студ. вузов. 2-е изд., перераб. и доп. М.: Физматлит, 2002. 300 с.
Устинов С.М., Зимницкий В.А. Вычислительная математика. – СПб.: БХВ-Петербург, 2009.
Формалев В.Ф., Ревезников Д.Л. Численные методы. М.: ФИЗМАТЛИТ, 2004. – 400 с.