37836

РЕШЕНИЕ СИСТЕМ НЕЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ МЕТОДОМ НЬЮТОНА

Лабораторная работа

Физика

Метод Ньютона Многие прикладные задачи радиофизики и электроники требуют решения систем нелинейных алгебраических уравнений СНАУ или в векторной форме 2. Для численного решения таких систем используются итерационные методы. Построение k1го приближения в этой схеме осуществляется посредством решения линейной системы 2.3 при этом вектор поправки находится путем решения системы линейных алгебраических уравнений 2.

Русский

2013-09-25

247 KB

25 чел.

Лабораторная работа № 2

РЕШЕНИЕ  СИСТЕМ  НЕЛИНЕЙНЫХ  АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ  МЕТОДОМ  НЬЮТОНА

ЦЕЛЬ РАБОТЫ: изучить и программно реализовать на языке высокого уровня метод Ньютона, исследовать его точность и эффективность на тестовых задачах.

Метод Ньютона

Многие прикладные задачи радиофизики и электроники требуют решения  систем нелинейных алгебраических  уравнений (СНАУ)

или в векторной форме

, (2.1)

где – вектор-столбец переменных, – вектор-столбец функций, – n-мерное векторное пространство.

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

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

,

при этом – [nxn]-неособенная матрица, задающая итерационный процесс, – числовой параметр. Построение (k+1)-го приближения в этой схеме осуществляется посредством решения линейной системы

, (2.2)

где

.

Если  для всех , здесь– [nxn]-единичная матрица, то итерационный метод называют явным, так как в этом случае  Метод является стационарным, когда  и  не зависят от номера итерации, и нестационарным в противном случае.

Качество итерационных методов оценивают по скорости сходимости, определяя ее как степень уменьшения нормы вектора погрешности при выполнении одного итерационного шага:

,

где – коэффициент сжатия, – порядок метода. Если , то итерационный метод имеет линейную сходимость, при – квадратичную сходимость.

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

Предположим, что известно k-е приближение  к точному решению  системы (2.1). Следующее (k+1)-е приближение в методе Ньютона вычисляется как

,  (2.3)

при этом вектор поправки  находится путем решения системы линейных алгебраических уравнений

, (2.4)

где  – [nxn]-матрица Якоби, определяемая следующим образом:

.

Из сравнения соотношений (2.2) и (2.4) следует, что метод Ньютона является одношаговым, неявным, нестационарным итерационным правилом.

На каждом шаге итерационного ньютоновского процесса необходимо вычислить вектор невязки , матрицу Якоби , решить систему линейных алгебраических уравнений (2.4) относительно вектора-поправки , определить новое приближение  по уточняю-щей формуле (2.3).

Критерием завершения итерационного процесса является одновременное выполнение условий:

     и   , (2.5)

где

,

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

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

Описание алгоритма

Алгоритм решения систем нелинейных алгебраических уравнений методом Ньютона реализуется следующим образом:

Алгоритм 2.1

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

.

  1.  Вычислить матрицу Якоби:

.

  1.  Решить систему линейных алгебраических уравнений

.

  1.  Уточнить решение:

.

  1.  Вычислить по формулам (2.5) и вывести на экран текущие значения  и , текущий номер итерации.
  2.  Проверить критерий (2.5) завершения итерационного процесса. Если этот критерий выполняется, то выйти из программы.
  3.  Проверить условие . Если это условие имеет место, то выйти из итерационного процесса с сообщением .

10. Положить  и перейти к п. 3.

Задание

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

Содержание электронного  отчета

  1.  Текст  программы.
  2.  Задача, результаты ее решения, характеристики , представленные таблично:

1

2

Таблица 2.1

Система уравнений

Начальное приближение

1

(1;   1)

2

(0.5;   0.2)

3

(-1.5;   1.5)

(-1;   1)

4

(1;   0)

5

(1;   1)

(2;   1.5)

(-3;   -1.5)

6

(1;   1)

7

(1;   -1)

(-1;   1)

8

(3;   2)

(3;   -2)

9

(1;   1)

10

(1;   1)

11

(1;   1); (-1;   -1)

12

(1;   0)

13

(1;   1);  (1;   -1)

14.

(1;   1)

Продолжение табл. 2.1

Система уравнений

Начальное приближение

15

(1.2;   1.3)

16

(0;   1)

17

(1;   1)

18

(1;   1)

(-1.;   1)

19

(1;   1)

(-1;   -1)

20

(0;   0;   0)

21

(1;   1;   1)

22

(1;   2.2;   2)

23

(1;   1;   1)

24

(1;   1;   1;   1)

(10; 10; 10; 10)

(100;  100;

100;  100)

14


 

А также другие работы, которые могут Вас заинтересовать

64055. Организация рационального использования земель сельского поселения 1.21 MB
  Организация рационального использования земель и их охраны является важнейшим фактором и главными условием развития всех отраслей экономики страны и определяющей частью земельной политики государства. Отечественная практика и мировой опыт показывают что землеустройство является основным механизмом осуществления проводимой земельной политики...
64056. Значение гигиены полости рта в профилактике заболеваний пародонта 1.39 MB
  Я выбрал эту тему потому что она актуально в наше время современные эпидемиологические данные свидетельствуют не только о значительной распространенности патологии пародонта у детей и взрослых но и о влиянии на частоту заболевания зубных отложений гигиены полости рта...
64057. Проект частной мастерской по ремонту кузовов легковых автомобилей в условиях с.Ребриха 1.26 MB
  Зубило используют для прямолинейной рубки деталей из тонких листов для обрезания фасок на толстых металлических листах при отсутствии шлифовальной машинки а также для рубки металлических прутков и т. Рубку листа осуществляют зубилом установленным под углом к обрабатываемой поверхности...
64058. Разработка технологического процесса сборки и сварки подкрановой балки 773.5 KB
  Устройство поясов из нескольких листов, сваренных по продольным кромкам, менее предпочтительно ввиду необходимости сварки протяжённых швов, сложности обеспечения плотного прилегания поясных листов друг к другу, неравномерности распределения напряжений.
64059. РАЗРАБОТКА ПРОИЗВОДСТВЕННОЙ И ТЕХНОЛОГИЧЕСКОЙ ПРОГРАММ ДЛЯ РЕСТОРАНА КЛАССА «ЛЮКС» НА 100 ПОСАДОЧНЫХ МЕСТ (ШВЕДСКИЙ СТОЛ) 633 KB
  Основными типами являются: Городские рестораны и кафе Располагаются в городах предполагают разнообразный набор блюд специализируются в предоставлении обедов ужинов работают в определенные часы и имеют многочисленную клиентуру; Рестораны и кафе самообслуживания...
64061. Создание благоприятного социально-психологического климата на уроках английского языка 643 KB
  В работе рассматривается вопрос использования различных подходов к созданию благоприятного психологического климата на уроках в совокупности с современными методиками обучения. В литературе имеется немало исследований посвященных проблеме формирования благоприятного социально-психологического климата...
64062. Фонетические, лексические и грамматические особенности говора пос. Береславка 378 KB
  Таким образом наблюдается смешение и параллельное существование черт Северорусского и Южнорусского наречий а также Среднерусских говоров. Одним из важнейших определителей отнесённости говора к определённому наречию является такая характеристика как тип предударного вокализма.
64063. Примениние холодной электродуговой сварки чугуна в строительстве 348 KB
  Всё изделие за исключением зоны прилегающей к сварочной ванне остаётся практически холодным на всём протяжении процесса сварки. Металл наплавляется узкими короткими валиками отдельными участками распределёнными по поверхности сварки с тем чтобы избежать...