Главная страница
Навигация по странице:

  • Исследуемая система: Дополнение

  • Итерация Значение x Значение y Значения первой функции

  • Точность Количество итераций Погрешность x Погрешность y

  • Решение методом релаксации

  • Курсовая работа Стрелковский С.П. 493110200003. Защищена с оценкой руководитель зудов Р. И подпись, дата инициалы, фамилия пояснительная записка к курсовой работе численное решение системы нелинейных уравнений По дисциплине Численные методы


    Скачать 103.87 Kb.
    НазваниеЗащищена с оценкой руководитель зудов Р. И подпись, дата инициалы, фамилия пояснительная записка к курсовой работе численное решение системы нелинейных уравнений По дисциплине Численные методы
    Дата25.05.2022
    Размер103.87 Kb.
    Формат файлаdocx
    Имя файлаКурсовая работа Стрелковский С.П. 493110200003.docx
    ТипРешение
    #549352

    Санкт-Петербургский политехнический университет Петра Великого

    Институт электроники и телекоммуникаций

    Высшая школа прикладной физики и космических технологий

    КУРСОВАЯ РАБОТА

    ЗАЩИЩЕНА С ОЦЕНКОЙ ______
    РУКОВОДИТЕЛЬ _____________________ Зудов Р.И

    подпись, дата инициалы, фамилия
    ПОЯСНИТЕЛЬНАЯ ЗАПИСКА

    К КУРСОВОЙ РАБОТЕ

    Численное решение системы

    нелинейных уравнений


    По дисциплине «Численные методы»

    РАБОТУ ВЫПОЛНИЛ

    Студент гр. № _____________________

    подпись, дата инициалы, фамилия


    Санкт-Петербург, 2022

    Оглавление

    Введение…………………………………………………………………………………………….…….3

    Исследуемая система………………………………………………………………………….……3
    Аналитическое решение……………………………………………………………………….….4
    Решение методом Ньютона………………………………………………………………….….5
    Решение методом релаксации….………………………………………………………….…7
    Выводы……………………………………………………………………………………………...……10

    Введение
    В этой курсовой работе исследуются различные методы решения систем нелинейных уравнений – трансцендентных уравнений в виде . Требуется сравнить эффективности методов по критерию числа потребовавшихся вычислений в зависимости от заданной точности.
    Исследуемые методы:


    1. Ньютона

    2. Релаксации (при оптимальном значении параметра релаксации и отличном от него)


    Исследуемая система:


    Дополнение



    Вся работа проводится в вычислительной среде MATLAB.

    Вычисления проводятся для трех различных точностей.

    Критерием останова итерационного процесса считается достижение разницей численных решений i-той и (i-1)-й итераций заданной точности.

    Аналитическое решение



    Чтобы понять верность решений, посчитанных методами Ньютона и релаксации, необходимо изначально самостоятельно решить заданную систему. Для этого воспользуемся онлайн-калькулятором.

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



    Как можно заметить из графика –

    Решение всего одно. Можно

    приступать к поиску решения

    системы на калькуляторе.

    К данным значениям должна

    сходиться система при ее решении

    от них будет считаться погрешность

    при окончании работы цикла

    итераций.

    Решение методом Ньютона
    Из различных источников можно слышать, что метод Ньютона – один из самых эффективных методов поиска численного решения уравнений или их систем. Мы подтвердим это или опровергнем в данном пункте работы.

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

    Чтобы численно решить уравнение (или подобную систему уравнений выпуклых функций двух и более переменных), его необходимо привести к такой эквивалентной записи, что итерационная процедура вычисления, последовательность которой по теореме Банаха о неподвижной точке стремится к решению, будет иметь вид , где – матрица Якоби, а i – номер итерации. Отсюда очевидно, что W, то есть det(W), должна быть отличной от нуля. Последнее суждение – условие сходимости.

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



    После описания метода, приступаем к вычислениям. Результаты приведем в виде трех таблиц. Первая и вторая – показательные. Они приводят пример работы алгоритма при одной из точностей, включают в себя номер выполненной итерации и значения полученные после исполнения этой итерации. Третья – сводная. Она демонстрирует результаты работы при различных точностях, включает в себя заданную точность, количество исполненных итераций и погрешность значений по итогу работы цикла.
    Условие сходимости выполнено ( ).


    Итерация

    Значение x

    Значение y

    Значения первой функции

    Значения второй функции

    1

    0.5816494931458931

    0.6640182638099812

    -0.00042419

    1.5137e-05

    2

    0.5819261487318413

    0.6645944447270368

    -3.1974e-08

    2.7135e-08

    3

    0.5819261517547013

    0.6645944783626472

    0

    -1.1102e-16

    Таблица 1. Работа метода Ньютона при точности 10-4.


    Итерация

    Значение x

    Значение y

    Значения первой функции

    Значения второй функции

    1

    0.5816494931458931

    0.6640182638099812

    -0.00042419

    1.5137e-05

    2

    0.5819261487318413

    0.6645944447270368

    -3.1974e-08

    2.7135e-08

    3

    0.5819261517547013

    0.6645944783626472

    0

    -1.1102e-16

    Таблица 2. Работа метода Ньютона при точности 10-8.


    Точность

    Количество итераций

    Погрешность x

    Погрешность y

    10-4

    3

    2.615175470133124e-05

    5.521637352789988e-06

    10-8

    3

    2.615175470133124e-05

    5.521637352789988e-06

    10-12

    4

    2.615175470144226e-05

    5.521637352789988e-06

    10-16

    4

    2.615175470144226e-05

    5.521637352789988e-06

    Таблица 3. Работа метода Ньютона при различных точностях.
    На практике мы увидели, что метод Ньютона работает быстро и точно, достигая значения с погрешностью, близкой к нулю за 3–4 единицы итераций.

    Решение методом релаксации



    О методе релаксации можно сказать, что он является модификацией метода простой итерации и чем-то схож с методом Якоби. Сам метод практически повторяет таковой для линейных. Его суть заключается в дополнительном смещении после вычисления компонент (n+1)-й итерации по формуле метода Зейделя, а итерационная процедура вычисления выглядит так:

    Здесь – параметр релаксации. Необходимым условием сходимости является неравенство , а достаточным использование оптимального параметра, находящегося по формуле:

    Где a и b это граничные точки окрестности решения. По ним оптимальное значение параметра составляет:


    Значение отличного от оптимального параметра возьмем равным:

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


    Итерация

    Значение x

    Значение y

    Значения первой функции

    Значения второй функции

    1

    0.5668742647106437

    0.6543800975352182

    -0.0020363

    -0.020019

    2

    0.5737071902376472

    0.6582883758007742

    -0.0018169

    -0.010214

    3

    0.5771933837387153

    0.660823828383443

    -0.0011787

    -0.0057447

    4

    0.5791541664185149

    0.6623599659645627

    -0.00071415

    -0.0033393

    5

    0.5802939222751571

    0.6632736795841919

    -0.00042478

    -0.0019614

    6

    0.5809633781738257

    0.6638143405794922

    -0.00025135

    -0.0011559

    7

    0.5813579120373132

    0.6641337819064133

    -0.00014851

    -0.00068199

    8

    0.581590690107891

    0.664322438338797

    -8.7706e-05

    -0.00040255

    9

    0.5817280901290038

    0.6644338420412322

    -5.1791e-05

    -0.00023766

    10

    0.5818092070174642

    0.6644996249574008

    -3.0582e-05

    -0.00014032

    Таблица 3. Работа метода релаксации с оптимальным параметром при точности 10-4.



    Итерация

    Значение x

    Значение y

    Значения первой функции

    Значения второй функции

    1

    0.5668742647106437

    0.6543800975352182

    -0.0020363

    -0.020019

    2

    0.5771933837387153

    0.660823828383443

    -0.0011787

    -0.0057447

    3

    0.5802939222751571

    0.6632736795841919

    -0.00042478

    -0.0019614

    4

    0.5813579120373132

    0.6641337819064133

    -0.00014851

    -0.00068199

    5

    0.5817280901290038

    0.6644338420412322

    -5.1791e-05

    -0.00023766

    6

    0.5818571002312557

    0.6645384688569449

    -1.8058e-05

    -8.285e-05

    7

    0.5819020762116377

    0.6645749493908435

    -6.2965e-06

    -2.8886e-05

    8

    0.5819260276872681

    0.6645943777230602

    -3.2448e-08

    -1.4886e-07

    9

    0.5819261511153384

    0.664594477844016

    -1.3754e-09

    -6.3098e-09

    10

    0.5819261516230654

    0.6645944782558681

    -3.4427e-11

    -1.5794e-10

    Таблица 3. Работа метода релаксации с оптимальным параметром при точности 10-10.


    Точность

    Количество итераций

    Погрешность x

    Погрешность y

    10-6

    19

    2.513123303982123e-05

    6.349451770470971e-06

    10-10

    36

    2.615162306540686e-05

    5.52174413182005e-06

    10-16

    65

    2.615175470133124e-05

    5.52163735290101e-06

    Таблица 4. Работа метода релаксации с оптимальным параметром при различных точностях.

    Итерация

    Значение x

    Значение y

    Значения первой функции

    Значения второй функции

    1

    0.5551821463806376

    0.6500409964024187

    -0.00015456

    -0.039114

    2

    0.5625730860363963

    0.651279837427896

    -0.0028346

    -0.025557

    3

    0.5592821479387742

    0.6505279516194495

    -0.0018358

    -0.031396

    4

    0.5652520239994032

    0.6521742791997948

    -0.0033721

    -0.021084

    5

    0.5721952178117653

    0.6559230486308542

    -0.0033626

    -0.010902

    6

    0.5777100440412889

    0.6604807228410198

    -0.0018039

    -0.0043727

    7

    0.5789367858631125

    0.6616457100469777

    -0.0013094

    -0.0030691

    8

    0.5800278708729323

    0.6627096044742893

    -0.00084302

    -0.0019369

    9

    0.5807180776224974

    0.6633917013737656

    -0.00053939

    -0.0012295

    10

    0.5811566339441454

    0.6638274521608554

    -0.00034432

    -0.00078233

    Таблица 5. Работа метода релаксации с отличном от оптимального параметром при точности 10-4


    Итерация

    Значение x

    Значение y

    Значения первой функции

    Значения второй функции

    1

    0.5625730860363963

    0.651279837427896

    -0.0028346

    -0.025557

    2

    0.5751896115474781

    0.6582389339063137

    -0.0026719

    -0.0072002

    3

    0.5807180776224974

    0.6633917013737656

    -0.00053939

    -0.0012295

    4

    0.5816136278193151

    0.6642827250376045

    -0.00014002

    -0.0003175

    5

    0.5819260335221983

    0.6645943603848202

    -5.2993e-08

    -1.2008e-07

    6

    0.5819261428758604

    0.6645944695029307

    -3.9796e-09

    -9.0177e-09

    7

    0.5819261472357298

    0.6645944738534091

    -2.0255e-09

    -4.5896e-09

    8

    0.5819261494547271

    0.664594476067627

    -1.0309e-09

    -2.3359e-09

    9

    0.5819261507087289

    0.6645944773189276

    -4.6882e-10

    -1.0623e-09

    10

    0.581926150919583

    0.6645944775293275

    -3.7431e-10

    -8.4818e-10

    Таблица 5. Работа метода релаксации с отличном от оптимального параметром при точности 10-10



    Точность

    Количество итераций

    Погрешность x

    Погрешность y

    10-6

    71

    1.853928636319502e-05

    1.311766356926203e-05

    10-10

    152

    2.61509195830234e-05

    5.522470672425506e-06

    10-16

    277

    2.615175470066511e-05

    5.521637353678166e-06

    Таблица 6. Работа метода релаксации с отличном от оптимального параметром при различных точностях.

    Вывод



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

    Таким образом при решении данной системы метод Ньютона показал себя значительно лучше и является менее трудозатратным.

    Список литературы

    1. Кудряшова Т. Ю. Численные методы: лекции СПБПУ, 2022.

    2. Решение систем нелинейных уравнений в MatLab: Статья - https://codetown.ru/matlab/reshenie-sistem-nelinejnyh-uravnenij/

    3. Релаксационные методы: Статья - http://www.machinelearning.ru/wiki/index.php?title=%D0%A0%D0%B5%D0%BB%D0%B0%D0%BA%D1%81%D0%B0%D1%86%D0%B8%D0%BE%D0%BD%D0%BD%D1%8B%D0%B5_%D0%BC%D0%B5%D1%82%D0%BE%D0%B4%D1%8B




    написать администратору сайта