Файл: Учебнометодическое пособие для выполнения лабораторных работ Волгоград, 2019 удк 519. 6(075. 8) Ббк в19я7 Печатается по решению редакционноиздательского совета Волгоградского государственного университета.docx

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 04.12.2023

Просмотров: 165

Скачиваний: 3

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

СОДЕРЖАНИЕ

1. ЦЕЛЬ И ЗАДАЧИ ОСВОЕНИЯ ДИСЦИПЛИНЫ

2. ЛАБОРАТОРНЫЕ РАБОТЫ

2.1 Лабораторная работа 1. «Теория погрешностей и машинная арифметика»

2.2 Лабораторная работа 2. « Решение нелинейных уравнений»

2) Расчетные формулы методов решения нелинейного уравнения .

2.3 Лабораторная работа 3. « Интерполирование функций»

2.4 Лабораторная работа 4.« Решение систем линейных уравнений» Цель: изучение основных численных методов решения систем линейных уравнений; разработка численного алгоритма и решения на ЭВМ систем линейных уравнений методами прогонки и итераций. Порядок выполнения работы Теоретическая часть Основные определения Система уравнений вида: или в сокращенной записи:называется линейной алгебраической системой из n уравнений с n-неизвестными xi (i=1,...,n). В матричной форме она записывается следующим образом:г де A - квадратная матрица, В и Х - векторы столбцы вида:1.2 Прямые методы1.Метод Крамера.Метод Крамера основан на использовании определителей в решении систем линейных уравнений. Это значительно ускоряет процесс решения. Метод Крамера может быть использован в решении системы стольких линейных уравнений, сколько в каждом уравнении неизвестных. Если определитель системы не равен нулю, то метод Крамера может быть использован в решении, если же равен нулю, то не может. Кроме того, метод Крамера может быть использован в решении систем линейных уравнений, имеющих единственное решение.Решение системы имеет вид:xi* = DAi / DA, i = 1, n, 1 2. Метод Гаусса. Этот метод основан на приведении методом исключения системы к треугольному виду (прямой ход ), а затем решение этой системы начиная с xn и т.д. (обратный ход). Если система сразу сводится к диагональному виду, то такой метод называется методом Жордана-Гаусса. Для уменьшения погрешности округления при сведении матрицы А к треугольному виду выбирается максимальный элемент в столбце или в строке и с помощью перестановок он делается главным (схема частичного выбора). Если главный элемент выбирается во всей матрице, то схема носит название глобального выбора. Алгоритм решения системы из n уравнений методом Гаусса с выбором главного элемента по столбцам выглядит следующим образом. Алгоритм прямого хода: Шаг 1. Примем k=1 Шаг 2. Выбираем рабочую строку.Если akk ≠ 0, то k-ая строка – рабочая.Если нет, меняем k-ю строку на m-ю (n≥m>k), в которой amk ≠ 0, . Если такой строки нет, система вырожденная, решение прекратить. Ш аг 3. Для строк i=k+1, k+2, …, n вычисляются новые значения коэффициентов: и новые правые части Шаг 4. Увеличиваем k = k + 1. Если k = n, прямой ход завершен, иначе алгоритм повторяется со второго шага. Получаем верхнюю треугольную матрицу А:А лгоритм обратного хода: Шаг 1. ВычислимШ аг 2. Вычислим:3.Метод прогонки. Данный метод применяется для решения трех диагональных систем. Метод состоит из двух этапов прямой прогонки - и обратной прогонки.Прямая прогонка: величину xi выразим через xi+1 с помощью коэффициентов Ai,Bi : xi= Ai xi+1+ Bi. Из первого уравнения находим значения A1 и B1. Подставляя x1=A1·x2+B1 во второе уравнение системы получаем: a2(A1x2+B1)+b2x2+c2x3=d2. Выражаем x2 и находим A2 и B2,т.е. зная A1 и B1 по этой формуле мы можем вычислить A2 и B2. Аналогично подставляя значение xi-1=Ai-1xi+Bi-1 в i уравнение имеем: ai(Ai-1xi+Bi-1)+bixi+cixi+1=di, i=1,2,...n. После n шагов получим значения An и Bn. Так как cn=0, то An=0. Следовательно, имеем: xn=Bn. Обратная прогонка состоит в последовательных вычислениях значений xn-1, xn-2 и т.д. до x1. 1.3 Итерационные методы:а) метод простой интерполяции (Метод Якоби);б) метод Гаусса-Зейделя Эти методы используются обычно при решении уравнений большого порядка, поскольку при итерационном процессе не накапливается ошибка округления. Задается некоторое приближенное решение x(0), затем производится цикл вычислений ( итераций ) и вычисляется новое приближение x(1). Процесс продолжается до получения решения с заданной точностью. Формула итерационного процесса:Условия завершения итерационного процесса d£e: Условие сходимости итерационного процесса:2. Контрольные вопросы Напишите матричную форму записи системы линейных уравнений. Перечислите основные прямые методы решения системы линейных уравнений. Напишите формулы Крамера. В чем суть метода Гаусса? В чем суть метода прогонки? Перечислите основные итерационные методы решения системы линейных уравнений. В чем суть итерационного метода? Напишите формулу итерационного процесса и условие завершения итерационного процесса. Практические задания 1.Решить систему линейных алгебраических уравнений, коэффициенты которой приведены в таблице заданий методами Крамера, Гаусса, прогонки, итерационным методом. Предварительно систему привести к трех диагональному виду. 2. Показать, что используемый метод имеет единственное решение в случае использования прямого метода или сходится в случае итерационного метода. 3. Написать программу и решить на ЭВМ с помощью этих методов систему уравнений и сравнить результаты. 4. Результаты занести в таблицу: Метод x1 x2 x3 x4 Крамера Гаусса Прогонки Гаусса-Зейделя ПРИЛОЖЕНИЕ 4ВАРИАНТЫ ДЛЯ ВЫПОЛНЕНИЯ ИНДИВИДУАЛЬНЫХ ЗАДАНИЙ

2.6 Лабораторная работа 6. «Методы численного интегрирования функций»

Рекомендуемая литература


  • шаг таблицы: h = хi+1 - xi (i = 0, 1, 2, ..., n) = const;

Первая интерполяционная формула Ньютона для интерполирования вперед:



П ри n=1 – линейная интерполяция:P1 (x) = у0 + t у0

При n=2 – квадратичная интерполяция:

В торая интерполяционная формула Ньютона для интерполирования назад:

График интерполяционного полинома у = j (х) проходит через заданные точки, т. е., значения полинома j (х) и данной функции у = f (х) совпадают в узлах х = хi (i = 0, 1, .., . n). Если функция f (х) сама является полиномом степени n, то имеет место тождественное равенство f (х) = j (х). В общем случае в точках, отличных от узлов интерполяции

R(x) = f (х) - j (х). Эта разность и есть погрешность интерполяции, и называется остаточным членом интерполяционной формулы.

  • О статочный член интерполяционного полинома Лагранжа:



- производная n+1-го порядка.

  • Остаточный член первой интерполяционной формулы Ньютона:





  • Остаточный член второй интерполяционной формулы Ньютона:



1.3. Экстраполяция

Как много значат для современного общества информационные технологии. С трудом можно представить крупные организации, предприятия, фирмы без применения компьютеров, Интернета, программных и аппаратных средств, так упрощающих повседневную, рутинную работу.

Между тем работа с информацией, представляющей для ее владельцев ценность, всегда сопровождаются угрозами нарушения ее целостности, расширением доступности, нарушением конфиденциальности. Вследствие реализации этих угроз рядовой пользователь домашнего ПК потеряет не много, в то время как организации, банки, крупные фирмы могут понести серьезный материальный ущерб.

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

Существует масса вариантов как все-таки сберечь «информационные активы» от потерь. Одним из таких способов является прогнозирование вторжений в информационную систему. Прогноз объединяет в себе опыт прошлых лет и материалов для построения стратегии защиты информации в бедующем. Данный материал является основой планирования защиты информации. Процесс прогнозирования может быть выражен в словесной, математической, графической или другой форме. Прогноз позволяет построить возможную модель проведения злоумышленником его атак. Возможность реализации данной модели, а значит и всего прогноза, есть величина вероятностная. На сегодняшний день существует около 150 методов прогнозирования.

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

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


  • четкое определение задачи и выдвижение гипотез о возможном развитии прогнозируемого объекта;

  • выбор системы параметров, определение единиц измерения, относящихся к каждому параметру;

  • сбор и систематизация данных (обязательно проводится проверка на однородность обрабатываемых данных);

  • выявление данных в ходе экстраполяции тенденции или симптомы изменения изучаемых величин.

Эффективность экстраполяции заключается в том, что она позволяет не только составить прогноз, но и своевременно зафиксировать объективно намечающиеся изменения, лежащие в зачатках назревающей тенденции.


К основным и наиболее простым методам экстраполяции можно отнести:

  • экстраполяция на основе среднего значения временного ряда (прогнозируемый уровень равен среднему значению уровней в прошлом.);

  • экстраполяция по скользящей и экспоненциальной средней (применяется для краткосрочного прогнозирования).

  • экстраполяция на основе среднего темпа (получение экстраполяции функции в результате использования среднего значения изменения темпа).

Следует выделить следующие типы функций экстраполяции:

  • однофакторные прогнозирующие функции;

  • многофакторные прогнозирующие функции.

К однофакторным прогнозирующим функциям относят функции, в которых прогнозируемый показатель зависит только от одного факториального признака. Зачастую в качестве факториального признака принимают время. В качестве фактора-аргумента в однофакторной прогнозирующей функции можно использовать не только время, но и другие факторы, если известна их количественная оценка на перспективу.

Общий вид формулы не всегда известен, следовательно, вид эмпирической формулы может быть произвольным. Предпочтение отдается обычно наиболее простым формулам, обладающим достаточной точностью.

Но каждый прогнозируемый доказательству t (t=1,2…,n)можно рассматривать как функцию от нескольких параметров:

  • в виде линейной многофакторной модели;

  • в виде нелинейной многофакторной модели (степенного типа), которая путем логарифмирования преобразуется в линейную.

Более сложные виды нелинейных многофакторных моделей редко используются в практике прогнозирования и планирования.

Коэффициенты в моделях определяются с помощью метода наименьших квадратов из системы нормальных уравнений, представляющих собой частные производные равные нулю.

1.4.Метод наименьших квадратов

Пусть для исходных данных xi, fi, i=1,…,N, выбран вид эмпирической зависимости:                 с неизвестными коэффициентами . Запишем сумму квадратов отклонений между вычисленными по эмпирической формуле и заданными опытными данными:

Параметры        будем находить из условия минимума функции   .

Известно, что в точке минимума все частные производные


от     по        равны нулю:

В качестве эмпирической функции рассмотрим полином:

Формула для определения суммы квадратов отклонений примет вид:



Вычислим производные:



Приравнивая эти выражения к нулю и собирая коэффициенты при неизвестных , получим следующую систему линейных уравнений:



Это нормальная система уравнений. Находим коэффициенты .

Если m=1, т.е.    , система нормальных уравнений примет вид:



При m=2 имеем:



Как правило, выбирают несколько эмпирических зависимостей. По МНК находят коэффициенты этих зависимостей и среди них находят наилучшую по минимальной сумме отклонений.

1.5. Прогнозирование

Прогнозирование проводилось в два этапа. На первом этапе был спрогнозирована статистика спама на ноябрь и декабрь 2013 года. Затем была рассмотрена статистика вложений вредоносных программ в письма, на основании которой был составлен прогноз вредоносного вложений на ноябрь и декабрь 2013 года. В том и другом случае были задействованы однофакторные функции, а именно: линейная, квадратичная, гиперболическая. По результатам анализа квадратичная зависимость дала более точные результаты. Оптимальное количество факторов-аргументов можно установить с помощью так называемого метода исключений.

В качестве примера, иллюстрирующего метод прогнозирования с помощью экстраполяции, можно рассмотреть прогнозирование появления спама в электронных почтовых ящиках пользователей Рунета. Проблема спама как для создателей антивирусов, так и для обычных пользователей стоит довольно остро. Письма с красивой рекламой, заманчивыми предложениями, акциями, в том числе и предложениями о легком «счастье» от мошенников, по большей своей части не только «захламляют» почту пользователя, но и угрожают опасными вредоносными вложениями, которые могут с легкостью нарушить работу системы. Спектр вредоносного вложения богат и может похвастаться разнообразием, которое непрерывно дает пищу для размышлений разработчикам антивирусных программ и систем защиты информации. По результатам статистики за период с мая 2013г. по октябрь 2013г., предоставленной Лабораторией Касперского, были проведены расчеты экстраполяции функции, представленной в виде набора значений, в котором определяющим параметром было выбрано время.

2. Контрольные вопросы

  1. Каковы цель и постановка задачи интерполирования?

  2. В чем заключается и как решается задача интерполирования?

  3. Как построить интерполяционный многочлен по 1-й формуле Ньютона?

  4. Как построить интерполяционный многочлен по 2-й

формуле Ньютона?

  1. Как определяется остаточный член?


3. Практические задания

  1. Составить таблицу значений функции, приведенной в приложении 3, разбив заданный интервал на 10 отрезков с h=const.

  2. Составить таблицу конечных разностей.

  3. По первой и второй интерполяционным формулам Ньютона определить значение функции в точках

х = х0±0.25h;

х = х0 ±0.5h;

х = xl0 ±0.25h;

х = xl0 ±0.5h;

  1. Построить эмпирическую формулу для заданной функции.

  2. Выполнить программную реализацию построения таблицы конечных разностей и эмпирической формулы, а так же расчет значений функций в точках из п. 3.

  3. Основные расчетные формулы свести в таблицу:

    Диапазон ячеек

    Формула







  4. Результаты вычислений свести в таблицу вида :

x

y(x)

q

Pn(x)

Погрешность