- •В.В. Чуркин численные методы
- •Содержание
- •Нелинейных уравнений Краткие сведения
- •1 Метод деления пополам (метод дихотомии, метод бисекций)
- •2 Метод хорд
- •3 Метод касательных (метод Ньютона)
- •4 Метод секущих
- •5 Метод итераций
- •5 Комбинированные методы решений нелинейных уравнений
- •Решение нелинейных уравнений в системе Mathcad
- •Пример построения графика функции и решения нелинейного уравнения
- •Лабораторная работа 1
- •Задание
- •Контрольные вопросы
- •Краткие сведения
- •Интерполирование в системе Mathcad
- •Лабораторная работа 2
- •Содержание отчета
- •Контрольные вопросы
- •Краткие сведения Алгебра матриц
- •Алгоритмы формирования матриц
- •Методы разложения матриц
- •Методы обращения матриц
- •Операции с векторами и матрицами в системе Mathcad
- •Лабораторная работа 3
- •Задание
- •Содержание отчета
- •Контрольные вопросы
- •Методы решений систем линейных алгебраических уравнений (слау) Краткие сведения
- •Прямые методы решений слау Метод Гаусса
- •Метод ортогонализации строк
- •Метод решения системы с ленточными матрицами
- •Метод Холецкого
- •Метод квадратного корня Пусть требуется решить слау с симметрической положительно определенной матрицей Матрица приводится к виду где
- •Метод прогонки
- •Метод вращений
- •Итерационные методы решений слау
- •Метод релаксации
- •Вычисление матричных выражений
- •Пример решения слау в системе Mathcad
- •Лабораторная работа 4
- •Контрольные вопросы
- •Лабораторная работа 5
- •Контрольные вопросы
- •Краткие сведения
- •Выполнение аппроксимации (регрессии) в системе Mathcad
- •Пример проведения регрессий – линейной и линейной общего вида
- •Лабораторная работа 6
- •Задание
- •Содержание отчета
- •Контрольные вопросы
- •Краткие сведения
- •Вычисление первообразных и интегралов в системе Mathcad
- •Лабораторная работа 7
- •Задание
- •Варианты вычисляемых интегралов и методов (формул) вычислений представлены в таблицах 1 и 2
- •Контрольные вопросы
- •Краткие сведения
- •Метод интегрирования оду с помощью ряда Тейлора
- •Метод Эйлера
- •Метод Рунге-Кутта третьего порядка рк3
- •Ошибки методов
- •Интегрирование систем оду и оду высших порядков
- •Методы прогноза и коррекции
- •Первый вариант метода Адамса
- •Второй вариант метода Адамса
- •Метод на основе методов Милна и Адамса-Башфорта
- •Метод Хемминга
- •Интегрирование систем оду в системе Mathcad
- •Пример 2
- •Контрольные вопросы
Содержание отчета
Задание
Формулы метода
Результаты выполнения задания в Mathcad’е .
Блок-схема алгоритма и таблица идентификаторов
Исходный код
Результаты выполнения работы в виде значений коэффициентов интерполирующего полинома и графика полученного полинома с нанесенными на график исходными данными
Библиографический список.
Контрольные вопросы
Что означают понятия интерполяции и экстраполяции?
Как решают задачу интерполяции и экстраполяции?
Каков результат решения задачи интерполяции и экстраполяции?
Какие методы интерполяции и почему применяют при небольшом числе узлов интерполяции?
Какие методы интерполяции и почему применяют при большом числе узлов интерполяции?
Сравните интерполяции по Лагранжу и сплайнами.
Какие недостатки имеют кусочно-линейная интерполяция и экстраполяция?
Расскажите о достоинствах и недостатках тригонометрической интерполяции.
Что нужно сделать с исходными данными, чтобы получить на графике существенное отличие сплайновой интерполяции от кусочно-линейной?
Как изменяется качество интерполяции для перечисленных методов при увеличении числа узлов интерполяции?
Сравните методы интерполяции по объему вычислений, т.е. по затратам машинного времени.
АЛГЕБРА, ФОРМИРОВАНИЕ,
РАЗЛОЖЕНИЕ И ОБРАЩЕНИЕ МАТРИЦ
Краткие сведения Алгебра матриц
Система чисел (действительных или комплексных)расположенных в прямоугольной таблице изnстрок и mстолбцов
,
называется матрицей (числовой)Числа - элементы матрицыПользуются также сокращенными записями и причем говорятчто матрица имеет размер
Если то матрицу называютпрямоугольнойа в случае-квадратнойМатрица размераназываетсявектором-строкойа размера-вектором-столбцомЧисло (скаляр) можно считать матрицей размераКвадратная матрицаимеющая ненулевые элементы только на главной диагоналиназываетсядиагональнойЕсли в диагональной матрице ненулевые элементы равны 1то матрицу называютединичнойи обозначают через
С квадратной матрицей связан определитель(число)определяемое по правилугде сумма распространена на всевозможные перестановкиэлементовиследовательносодержитслагаемыхпричемесли перестановка четнаяиесли перестановка нечетная.
Рангданной матрицы размераесть такое число, что по крайней мере один определитель -го порядкаполучаемый при удалении некоторых строк и (или) столбцовотличен от нуляа все определители-го порядка равны нулю. Дляневырожденнойматрицы.
След(spur) квадратной матрицыравен сумме элементов на главной диагонали:
Для матриц произвольного типа и размера используют три нормы:
(m-норма);
(l-норма);
(k-норма).
Матрицы равныесли они имеют один и тот же размер и равны их соответствующие элементыСумма и разность матриц определена для матриц одинакового типаПравило умножения матрицэлементматрицы-произведения представляет собой сумму произведений элементовi-й строки первой матрицы на элементыj- го столбца второй матрицыСледовательноесли размер первой матрицыа второй -перемножение возможно приа матрица-произведение имеет размер
Определитель матрицы-произведения равен произведению определителей матриц-сомножителей
Транспонированнаяматрицаполучается из исходнойзаменой строк столбцами
Матрица называется симметрической (симметричной)если она совпадает со своей транспонированнойт.е. если. В случаематрицу называюткососимметрической.
Обратнойпо отношению к даннойназывается матрицакотораябудучи умноженной как справатак и слева на данную матрицудает единичную матрицуНахождение обратной матрицы для данной называется обращением данной матрицы
При вычислении матрицыобратной матрице, точность оценивается с помощью модифицированнойk-нормы матрицы ошибок, равной,
Квадратная матрица называется неособеннойесли определитель ее не равен нулюВ противном случае матрица называетсяособеннойОбратную матрицу имеет всякая неособенная матрица
Присоединенную(илисоюзную) матрицу получают транспонированием матрицысоставленной из алгебраических дополнений элементов исходной матрицыРазделив присоединенную матрицу на определитель исходной матрицыполучают матрицуобратную для исходной
С помощью горизонтальных и вертикальных перегородокидущих вдоль всей матрицыпоследнюю разбивают на матрицы низших порядков (клетки или блоки)Матрицуразбитую на клеткиназываютклеточнойилиблочной
Разбиение матрицы на клетки может быть осуществлено различными способамиЕсли все клетки являются квадратными матрицами и главные диагонали всех клеток находятся на главной диагонали исходной матрицыа вне клеток стоят нулито такую клеточную матрицу называютквазидиагональной
Окаймленныематрицы представляют собой другой важный частный случай клеточных матрицИз исходной матрицыокаймленную матрицу матрицуполучают следующим образом , где
- матрица порядкаn-1; - матрица-столбец; - матрица-строка и- число
Клеточные матрицы одинакового типа и с одинаковым разбиением называют конформными. Удобство клеточных матриц состоит в томчто действия над ними совершаются формально по тем же правиламчто и над обыкновенными матрицами.
Матрицавсе элементы которой равны нулюназываетсянулевойи обозначается через.
Матрицы общего виданазываютплотнымиилизаполненными
Практические задачинапример анализ протяженной электрической цепи с локальными связями между элементамиприводят кразреженнымматрицамв которых нулевых элементов значительно большечем ненулевых. Пример разреженной матрицы –трехдиагональнаяматрицавсе ненулевые элементы которой расположены на главной и двух соседних с ней диагоналях.
Матрица называется ленточнойс полушириной лентыравной,если для . Все ненулевые элементы расположены наближайших к главной диагоналях матрицы; числоназывают шириной ленты. В случае ленточная матрица является разреженной.
Если в квадратной матрице элементыстоящие выше (ниже) главной диагоналиравны нулюто такую матрицу называютнижней(верхней)треугольнойматрицейДиагональная матрица является частным случаем как верхнейтак и нижней треугольной матрицыОпределитель треугольной матрицы равен произведению ее диагональных элементов.
Всякую квадратную матрицу имеющую отличные от нуля главные диагональные миноры ….; можно представить в виде произведения двух треугольных матриц (верхней и нижней)причем это разложение будет единственнымесли зафиксировать диагональные элементы одной из матриц (напримерпринять их равными 1).
Если главные диагональные миноры симметрической матрицы положительныто она называетсяположительно определенной.
Действительнаяматрицаявляетсяортогональнойесли ее транспонированная матрицасовпадает с обратнойт.е.или . Ортогональная матрица имеет следующие свойства.
Строки (столбцы) попарно ортогональныт.е.
приипри.
2. Сумма квадратов элементов каждой строки (столбца) равна 1.
3. Определитель равен 1.
4. Транспонированная и обратная матрицы ортогональной матрицы также являются ортогональными матрицами.
Две матрицы называются эквивалентнымиесли одна получается из другой с помощью элементарных преобразований: перестановка двух строк или столбцов; умножение всех элементов какой-либо строки (столбца) на одно и то же число; прибавление к элементам какой-либо строки (столбца) соответствующих элементов другой строки (столбца)умноженных на одно и то же число.
Указание для изменения порядка следования строк (столбцов) матрицы задают в виде вектора транспозиции. Например, если для строк матрицы размера 5*6 задан вектор транспозиции , то матрица не будет изменена, а в случаепоменяются местами строки: первая – с пятой, третья – со второй.