Банк рефератов содержит более 364 тысяч рефератов, курсовых и дипломных работ, шпаргалок и докладов по различным дисциплинам: истории, психологии, экономике, менеджменту, философии, праву, экологии. А также изложения, сочинения по литературе, отчеты по практике, топики по английскому.
Полнотекстовый поиск
Всего работ:
364150
Теги названий
Разделы
Авиация и космонавтика (304)
Административное право (123)
Арбитражный процесс (23)
Архитектура (113)
Астрология (4)
Астрономия (4814)
Банковское дело (5227)
Безопасность жизнедеятельности (2616)
Биографии (3423)
Биология (4214)
Биология и химия (1518)
Биржевое дело (68)
Ботаника и сельское хоз-во (2836)
Бухгалтерский учет и аудит (8269)
Валютные отношения (50)
Ветеринария (50)
Военная кафедра (762)
ГДЗ (2)
География (5275)
Геодезия (30)
Геология (1222)
Геополитика (43)
Государство и право (20403)
Гражданское право и процесс (465)
Делопроизводство (19)
Деньги и кредит (108)
ЕГЭ (173)
Естествознание (96)
Журналистика (899)
ЗНО (54)
Зоология (34)
Издательское дело и полиграфия (476)
Инвестиции (106)
Иностранный язык (62792)
Информатика (3562)
Информатика, программирование (6444)
Исторические личности (2165)
История (21320)
История техники (766)
Кибернетика (64)
Коммуникации и связь (3145)
Компьютерные науки (60)
Косметология (17)
Краеведение и этнография (588)
Краткое содержание произведений (1000)
Криминалистика (106)
Криминология (48)
Криптология (3)
Кулинария (1167)
Культура и искусство (8485)
Культурология (537)
Литература : зарубежная (2044)
Литература и русский язык (11657)
Логика (532)
Логистика (21)
Маркетинг (7985)
Математика (3721)
Медицина, здоровье (10549)
Медицинские науки (88)
Международное публичное право (58)
Международное частное право (36)
Международные отношения (2257)
Менеджмент (12491)
Металлургия (91)
Москвоведение (797)
Музыка (1338)
Муниципальное право (24)
Налоги, налогообложение (214)
Наука и техника (1141)
Начертательная геометрия (3)
Оккультизм и уфология (8)
Остальные рефераты (21697)
Педагогика (7850)
Политология (3801)
Право (682)
Право, юриспруденция (2881)
Предпринимательство (475)
Прикладные науки (1)
Промышленность, производство (7100)
Психология (8694)
психология, педагогика (4121)
Радиоэлектроника (443)
Реклама (952)
Религия и мифология (2967)
Риторика (23)
Сексология (748)
Социология (4876)
Статистика (95)
Страхование (107)
Строительные науки (7)
Строительство (2004)
Схемотехника (15)
Таможенная система (663)
Теория государства и права (240)
Теория организации (39)
Теплотехника (25)
Технология (624)
Товароведение (16)
Транспорт (2652)
Трудовое право (136)
Туризм (90)
Уголовное право и процесс (406)
Управление (95)
Управленческие науки (24)
Физика (3463)
Физкультура и спорт (4482)
Философия (7216)
Финансовые науки (4592)
Финансы (5386)
Фотография (3)
Химия (2244)
Хозяйственное право (23)
Цифровые устройства (29)
Экологическое право (35)
Экология (4517)
Экономика (20645)
Экономико-математическое моделирование (666)
Экономическая география (119)
Экономическая теория (2573)
Этика (889)
Юриспруденция (288)
Языковедение (148)
Языкознание, филология (1140)

Курсовая работа: Рішення систем нелінійних рівнянь. Метод ітерацій. Метод Ньютона–Канторовича

Название: Рішення систем нелінійних рівнянь. Метод ітерацій. Метод Ньютона–Канторовича
Раздел: Рефераты по экономико-математическому моделированию
Тип: курсовая работа Добавлен 07:08:21 24 апреля 2010 Похожие работы
Просмотров: 2682 Комментариев: 2 Оценило: 2 человек Средний балл: 3 Оценка: неизвестно     Скачать

ОБЛАСНИЙ КОМУНАЛЬНИЙ ВИЩИЙ НАВЧАЛЬНИЙ ЗАКЛАД "ІНСТИТУТ ПІДПРИЄМНИЦТВА "СТРАТЕГІЯ"

КАФЕДРА ЕКОНОМІЧНОЇ КІБЕРНЕТИКИ

Курсова робота

З дисципліни: "Обчислювальні методи"

На тему: "Рішення систем нелінійних рівнянь. Метод ітерацій. Метод Ньютона - Канторовича."

Студента Іощенка І.Г.

группа С-05-51

Керівник Андрейшина Н.Б.

Філімоненко М.І.

м. Жовті Води 2007

Зміст

Вступ

1. Рішення систем нелінійних рівнянь

1.1 Метод ітерацій

1.1.1 Приклад рішення системи нелінійних рівнянь методом ітерацій

1.2 Метод найшвидшого спуску

1.2.1 Приклад рішення системи нелінійних рівнянь методом спуска

1.3 Метод Ньютона-Канторовича

Вступ

При рішенні систем нелінійних і трансцендентних рівнянь дуже складно знайти точне рішення, тому точним рішення рівняння не є. Задача пошуку кореня системи рівняння може вважатися практично вирішеною, якщо ми зуміємо визначити корінь з потрібним ступенем точності і вказати межі можливої погрішності. Умови збіжності метода Ньютона для системи досліджувалися Виллерсом, Стениним, Канторовичем.

У наш час рішення систем нелінійних рівнянь досить актуальна тема, адже її можна застосовувати на практиці для рішення кола задач. Прикладом цього є задачі, які виникають у геодезії.

Цілю моєї курсової роботи є опис методів рішення систем нелінійних рівнянь, а також продемонструвати на практиці рішення системи рівнянь методом Ньютона - Канторовича та написання програми до цього методу.

1. Рішення систем нелінійних рівнянь

Задачі, які виникають при математичній обробці результатів вимірювання, як правило, зводяться до рішення нелінійних систем алгебраїчних або трансцендентних рівнянь:

або у векторній формі

F (X) = 0.

Як і у випадку одного рівняння, рішення нелінійних систем рівнянь поділяється на два етапи:

знаходження приблизного рішення системи;

уточнення приблизного рішення.

Для знаходження приблизного значення коренів системи рівнянь не існує загальних методів. Завжди кожна нелінійна система повинна розглядатися як спеціальна задача.

Для уточнення коренів розробленні загальні методи. Найбільш розповсюдженні в нинішній час є метод ітерацій, метод спуска, метод Ньютона та деякі їх модифікації.

1.1 Метод ітерацій

Нехай дана система нелінійних рівнянь спеціального виду


(1)

де функції , ,... ., дійсно визначенні та непереривні на деякій області ізольованого рішення цієї системи.

Розглядаючи вектори і (x) = (1 (x), 2 (x), …. .,n (x)), систему (1) можна записати у виді:

x = (x) (2)

Наприклад, для рішення системи двох нелінійних рівнянь з двома невідомими

потрібно перейти до рівностей:

Нехай вибрано початкове приближення (,), тоді


і k+1 приближення буде розраховуватися за формулами

Відомо, що процес ітерації зводиться до рішення системи, якщо усі числа матриці

по модулю менше одиниці. Більш простою вимогою, використовуваною на практиці, є наступне: сума модулів частних похідних по кожному стовбці матриці повинна бути менша одиниці

У випадку використання методу ітерацій до системи n рівнянь, k+1 ітерація буде будуватися по формулам


Тоді вимога сходження матиме вигляд:

Слід відмітити, що ця вимога виповняється для дуже малого числа функцій, і тому метод ітерації дуже рідко використовується на практиці, не дивлячись на його простоту.

1.1.1 Приклад рішення системи нелінійних рівнянь методом ітерацій

Рішить систему рівнянь

Ця система еквівалентна системі рівнянь:


Виберемо початкові приближення та провіримо умови

сходження процесу. Часні похідні мають вигляд

Маємо

Звідси слідує, що процес сходиться. Розрахунки на правому приближенні дають:

x (1 ) = 1+0.85=1.85

y (1 ) =0.842-1.32=-0.478

x (2 ) =0.888+0.85=1.738

y (2 ) =0.961-1.32=0.359

x (3 ) =0.936+0.85=1.786

y (3 ) =0.986-1.32=0.334

x (4 ) =0.945+0.85=1.795

y ( 4) =0.977-1.32=0.343

x (5) =0.9408+0.85=1.7908

y (5) =0.9750-1.32= - 0.3450

x (6 ) = 0.9411+0.85=1.7911

y (6 ) = 0.9759-1.32=0.3441

x ( 7) = 0.9414+0.85=1.7914

y ( 7) = 0.9758-1.32=-0.3442.

1.2 Метод найшвидшого спуску

Нехай маємо систему рівнянь:

або в матричному вигляді:

де

Допустимо, що функція дійсно непереривна та непреривно диференційована в загальній області визначення. Розглянемо функцію

Тоді рішення даної системи зводиться до мінімізації цієї функції.

Для мінімізації по методу спуску вибирається початковий вектор Х0 , а потім шукається напрямлення спуска до рішення , таке щоб

для векторів Х (1) виду . Тут - скалярна величина, постійна для даної ітерації і знаходить величину шагу за напрямом .

Методи спуску розрізняються в залежності від вибору напрямлення спуска. Одним із найкращих направлень є напрямлення градієнта

Функція Ф (Х (і)) задається в n-мірному просторі сімейства гіперповерхонь і градієнт вирішує напрям найшвидшого спуска. Тому саме воно використовується у методі найшвидшого спуска для мінімізації функції.

Другою проблемою в методах найшвидшого спуску є вибір величини шагу , на який потрібно про двинутися вздовж напряму зменшення функції.

Спробуємо вибрати оптимальний шаг для - ітерації методу найшвидшого спуска і побудувати вектор

для якого функція приймає менше значення, чим . Розкладемо функцію


в ряд Тейлора та обмежившись членами другого порядку меншості получимо

(3)

Тоді значення , для якого функція прийме мінімальне значення, визначається із умови Про диференціювавши рівняння (3) по і враховуючи, що получимо

(4)

Оскільки в методі найшвидшого спуску компоненти градієнта мають вигляд

то формула (4) після підстановки цих рівнянь перейде до вигляду


(5)

Формула (5) дуже складна оскільки потребує рахування других часних похідних.

На практиці завжди використовується наступний варіант знаходження .

Нехай значення Ф (Х) змінюється вздовж напрямку градієнта . Розглянемо точку пересікання кривої та касатільної в точці з осю .

Вона буде розраховуватися наступним чином:

. (6)

Як бачимо, в цьому випадку рахується просто, але сходження метода може бути дуже повільно. Тому інколи на практиці використовують наступну модифікацію.

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


(7)

Практика показує, що хоча цей варіант більш громіздкий, так як у порівнянні з формулою (5) доводиться додатково рахувати два значення функції , але метод сходиться набагато швидше.

Інколи характер Ф (Х) такий, що аналітичне рівняння для частних похідних має надто складний вигляд і рахувати їх надто складно.

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

Порядок рахування в методі найшвидшого спуска наступний:

знаходиться аналітичне рівняння для градієнта ;

вибирають початкове приближення вектора невідомих ;

вираховують координати градієнта в точці ;

вираховують шаг по градієнту по формулам (6) або (7);

вираховують уточнений вектор невідомих .

Далі процес повторюється з пункту 3 до сходження.

1.2.1 Приклад рішення системи нелінійних рівнянь методом спуска

Методом найшвидшого спуска приблизно розрахувати корені системи

розміщенні в області початку координат.

Маємо:

Тут та

Підставляємо нульове приближення, будемо мати:

та

по формулам получимо перше приближення

Аналогічно находимо друге приближення . Маємо:


.

1.3 Метод Ньютона-Канторовича

Метод Ньютона-Канторовича, придатний для проведення розрахунків в Excel. Як і в методі Ньтона для нелінійних рівнянь для знаходження кореня системи нелінійних рівнянь необхідно спочатку якимсь чином знайти початкове наближення до цього кореня (тобто вектор

),

а потім вже використовуються ітераційні формули методу проводиться його уточнення до досягнення заданої точності. Виклад методу (і його використання) зручніше проводити в матричній формі запису. При цьому, окрім векторів, , и (. (i - номер ітерації,, i ³ 0) ) використовується також матриця A (розмірності n ´ n), що складається з приватних похідних по всіх компонентах вектора :

:

Розглянемо ці методи для випадку n=2, тобто коли рівнянь в системі два і невідомих теж дві. В цьому випадку

, та .

Ідея методу полягає в розкладанні вектор-функції в ряд Тейлора в околиці початкового наближення із збереженням тільки доданків першого ступеня. Позначимо найдене (якимсь чином) початкове приближення до шуканого кореня через . Тоді можна приблизно записати


, (8)

На основі формула (8) будується ітараційна формула. А саме, вибирається так, щоб .

Тоді (у загальному вигляді) ітераційна формула матиме вигляд

(9)

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

(10)

У координатному вигляді формула (10) представляє систему з двох рівнянь щодо двох невідомих xi+1 и yi+1.

У матричному вигляді рішення її матиме вигляд

допоміжний вектор-стовпець z, що містить n елементів.


(11)

Ітераційна формула методу в матричному записі має наступний вигляд

zj = - A-1 (xj) × F (xj)

xj +1 = xj + zj , (12)

тут j - номер ітерації, - початкове наближення шуканого кореня. Процес ітерацій завершується, якщо всі елементи останнього вектора z по абсолютній величині стануть менше заданої точності (кажучи точніше, коли норма вектора z стане менше заданої точності).

Обчислення даним методом зручно проводити в Excel з використанням функцій матричної алгебри. Результати розрахунків представляються у вигляді таблиці.

Для випадку n=2 система рівнянь найчастіше має такий вигляд:

Як змінна х1 тут виступає змінна х , а як змінна х2 - змінна y. Матриця А, вектори F і z в цьому випадку приймуть вигляд:

А = , F = , z = ,

Порядок рішення системи нелінійних рівнянь методом Ньютона-Канторовича полягає в послідовному виконанні наступних дій:

Знайти початкове (нульове) наближення х0 шуканого кореня заданої системи рівнянь. Для випадку n=2 це можна зробити графічним методом, побудувавши графіки кожної з функцій і приблизно визначивши координати точок перетинів графіків. В цьому випадку вектор початкового наближення може мати вигляд ;

Привести задану систему до вигляду (1), перенести все з правої частини рівняння в ліву;

Записати в аналітичному вигляді матрицю А, використовуючи формулу (8);

Приймемо j=0;

Підставимо значення хj в аналітичні вирази для матриці А і вектора F;

Знайдемо зворотну матрицю А-1 ;

По формулах (12) знайдемо вектор zj і вектор хj+1 ;;;

Знайдемо норму вектор zj ;

Якщо норма вектора zj більше заданої точності обчислення (норма більша за ε) - наростимо значення j на одиницю і повернемося до пункту 5 цього переліку;

За знайдене рішення приймемо останнього набутого значення вектора х .

Оценить/Добавить комментарий
Имя
Оценка
Комментарии:
Где скачать еще рефератов? Здесь: letsdoit777.blogspot.com
Евгений07:34:02 19 марта 2016
Кто еще хочет зарабатывать от 9000 рублей в день "Чистых Денег"? Узнайте как: business1777.blogspot.com ! Cпециально для студентов!
15:58:38 25 ноября 2015

Работы, похожие на Курсовая работа: Рішення систем нелінійних рівнянь. Метод ітерацій. Метод Ньютона–Канторовича
Економічна модель оптимізації закупівель та поставок кондитерських ...
Міністерство освіти і науки України Національний гірничий університет Кафедра Економічної кібернетики та інформаційних технологій ПОЯСНЮВАЛЬНА ЗАПИСКА ...
Зокрема, якщо в околиці стаціонарної точки х0 функція W(x) двічі безупинно дифференцируема і її матриця Гессе позитивно визначена в х0, то функція W(x) буде опуклої, а х0- точкою ...
Викладемо ідею методу для випадку функції двох перемінних f (x, y).
Раздел: Рефераты по экономико-математическому моделированию
Тип: дипломная работа Просмотров: 2204 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать
Интегралы, дифуры, матрицы
Інтегральне числення Невизначений інтеграл 1. Поняття первісної Означення: Функція F(x) називається первісною для ф-ії f(x) на проміжку І, якщо на ...
В тривимірному випадку для ф-ії u=f(x;y;z), коли дуга кривої L задана параметричними рівняннями x=x(t), y=y(t), z=z(t), a $ t $ b. Формула має вигляд:
Вектор з координатами , який характеризує напрям максимального зростання ф-ії z=f(x;y) в точці P0=(x0;y0)
Раздел: Топики по английскому языку
Тип: шпаргалка Просмотров: 2015 Комментариев: 7 Похожие работы
Оценило: 2 человек Средний балл: 2 Оценка: неизвестно     Скачать
Числові методи
МІНІСТЕРСТВО ОСВІТИ УКРАЇНИ ЧЕРНІВЕЦЬКИЙ ДЕРЖАВНИЙ УНІВЕРСИТЕТ ІМ. Ю. ФЕДЬКОВИЧА КОНТРОЛЬНА РОБОТА з дисципліни " Числові методи " Варіант 16. Виконав ...
- вектор невідомих, - вектор вільних членів, - обернена матриця до матриці з коєфіцієнтів при невідомих.
Y- вектор значень функції f(x,y) в точках x[0],x[1],...x[n].
Раздел: Рефераты по информатике, программированию
Тип: контрольная работа Просмотров: 143 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать
Дослідження процесів масопереносу при фільтрації підземних вод
дипломний проекту на тему: " Дослідження процесів масопереносу при фільтрації підземних вод " Івано-Франківськ 2008 р. Зміст Вступ 1.Моделювання й ...
Швидкість поширення розчинного у фільтраційному потоці речовини U(x,y,t) у цьому випадку дорівнює дійсній швидкості руху підземних вод V(x,y,t) , яка зв'язана зі швидкістю ...
де D - коефіцієнт конвективної дифузії в м/сут, c й N - концентрація речовин, що дифундують, у г/л або кг/м відповідно в рідкій і твердій фазах; vx(x, y, t) і vy(x, y, t ...
Раздел: Рефераты по экологии
Тип: дипломная работа Просмотров: 281 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать
Численные методы
... РІВНЯНЬ. Розглянемо чисельні методи розв"язання систем лінійних алгебраїчних рівнянь Ax=f T (1) де A - матриця m*m, x = ( x1, x2 , ... ,xm ) - шуканий ...
де A - матриця m*m, x = ( x1, x2 , ... ,xm ) - шуканий вектор,
Матриця системи (5) має вигляд:
Раздел: Рефераты по математике
Тип: реферат Просмотров: 495 Комментариев: 2 Похожие работы
Оценило: 1 человек Средний балл: 3 Оценка: неизвестно     Скачать
Методи розв"язування одновимірних та багатовимірних нелінійних ...
Міністерство освіти і науки України Полтавський національний технічний університет імені Юрія Кондратюка Факультет інформаційно-телекомунікаційних ...
3. Розв"язання задачі мінімізації за допомогою методу Ньютона і методу найшвидшого спуску
Отже матриця Гессе матиме вигляд:
Раздел: Рефераты по математике
Тип: курсовая работа Просмотров: 2123 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать
Регресійний аналіз інтервальних даних
МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ ДНІПРОПЕТРОВСЬКИЙ НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ ДИПЛОМНА РОБОТА Регресійний аналіз інтервальних даних Виконала ...
В якості локальних змінних оголошенні kol_strok- кількість строк матриці matrica_X що вводиться з клавіатури користувачем, kol_stolbcov- кількість стовпчиків даної матриці, X ...
Цей поріг, де починається стабілізація, не залежить від вектора похибки і приблизно дорівнює 50.
Раздел: Рефераты по математике
Тип: дипломная работа Просмотров: 215 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать

Все работы, похожие на Курсовая работа: Рішення систем нелінійних рівнянь. Метод ітерацій. Метод Ньютона–Канторовича (3145)

Назад
Меню
Главная
Рефераты
Благодарности
Опрос
Станете ли вы заказывать работу за деньги, если не найдете ее в Интернете?

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



Результаты(150725)
Комментарии (1839)
Copyright © 2005-2016 BestReferat.ru bestreferat@mail.ru       реклама на сайте

Рейтинг@Mail.ru