Банк рефератов содержит более 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)

Курсовая работа: Дослідження чисельних методів вирішення нелінійних рівнянь

Название: Дослідження чисельних методів вирішення нелінійних рівнянь
Раздел: Рефераты по информатике, программированию
Тип: курсовая работа Добавлен 15:24:05 01 декабря 2009 Похожие работы
Просмотров: 197 Комментариев: 2 Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать

Міністерство освіти і науки України

Вінницький національний технічний університет

Інститут АЕКСУ

Кафедра АІВТ

Дослідження чисельних методів вирішення нелінійних рівнянь

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

з дисципліни:

“Обчислювальні методи та застосування ЕОМ”

Виконав: ст. гр. 3АВ-06_________ Пересунько А.О.

Перевірив: професор, д.т.н. каф.АІВТ___Квєтний Р.Н.

2008


Анотація

В даній курсовій роботі розглянуто наближені методи розв’язку нелінійних рівнянь для методів Ньютона та хорд, складено блок-схеми та написано програму, за допомогою якої розв’язується задане рівняння. Проведено аналіз як самого рівняння і методів його розв’язання, так і результатів обрахунку.


Зміст

Вступ

1.Короткі теоретичні відомості

2.Аналіз заданого рівняння

3.Алгоритми методів

3.1.Вибір інструментальних засобів

3.2.Вхідні та вихідні дані

3.3. Структура програми

3.4. Інструкція користувачеві

4. Аналіз результатів розрахунку

5. Висновки

Література

Додатки

Додаток А. Алгоритм методів

Додаток Б. Блок-схема програми

Додаток В. Лістінг програми


Вступ

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

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

Сучасний світ неможливо уявити без використання комп’ютерних технологій. Зараз комп’ютер використовується у багатьох сферах людського життя. Зараз обчислення залишаються одним із основних видів застосування ЕОМ. Хоча комп’ютер дуже швидко виконує прості арифметичні дії, без спеціальних програм він не в змозі проводити складні обчислення. Тому постає задача алгоритмізувати поставлене завдання, тобто перевести його в зрозумілу для ЕОМ форму.


1.Короткі теоретичні відомості

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

Метод половинного ділення

В цьому методі спочатку обчислюється значення функції в точках що розташовані через рівні інтервали на осі х. Коли f(xn ) if(xn +1 ) мають протилежні знаки, знаходять , f(xcp ). Якщо знак f(xcp ) збігається зі знаком f(xn ), то надалі замість хn використовується хср . Якщо ж f(xcp ) має знак, протилежний f(xn ), тобто збігається зі знаком f(xn +1 ), то на хср замінюється xn +1 . За умову припинення ітераційного процесу доцільно брати умову | xn +1 – xn | < e, де e - задана похибка. Похибка розвязку Δ через n ітерацій знаходиться в межах

Δ< . (1.1)

Метод має малу швидкість збіжності, оскільки інтервал, де знаходиться корінь, з кожним кроком зменшується не більше ніж в два рази.

Метод простої ітерації

Цей метод можна використовувати лише якщо доведена збіжність ітераційного алгоритму. В цьому методі процес розв’язання потрібно починати з пошуку інтервалу збіжності. Умовою збіжності є те що максимальне значення І-ї похідної правої частини рівняння Х=g(x) (1)(до такого вигляду потрібно привести вихідне рівняння f(x)=0 ) повинна бути менша за 1. Якщо умова не виконується, то алгоритм не збіжний. Коли в інтервалі збіжності немає коренів, треба застосовувати інші методи або приходити до рівняння (1) через інші способи.

Похибка ж методу на n – ій ітерації обчислюється так:

Δ< (1.2)

В даній курсовій роботі розглядаються два методи розв’язку нелінійних рівнянь – це метод Ньютона та метод хорд, тому розглянемо їх більш детально. Метод хибного положення (хорд)

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

. (1.3)

Після цього визначають f(xn +1 ) і порівнюють його з f(xn ). Надалі користуються xn +1 замість того значення, з яким воно збіглося за знаком. Якщо | xn +1 – xn | < e, то вся процедура повторюється спочатку.

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

Похибка розв'язку оцінюється за формулою:

, (1.4)

де М1 , m1 – відповідно, найбільше та найменше значення модуля першої похідної на відрізку.

Рис 1.1.Метод хорд

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

Метод Ньютона полягає в побудові дотичної до графіка функції в обраній точці. Наступне наближення знаходиться як точка перетину дотичної з віссю ОХ.В основі цього методу лежить розкладання функції в ряд Тейлора:

(1.5)

Члени що містять h у другому і більших степенях відкидаються і в результаті отримується наближена формула для оцінки Хn +1 :

хn +1 = хn , (1.6)

але оскільки цей метод є наближеним, то логічно буде якщо для нього задавати певну похибку і тоді наближене значення кореня буде визначатися з виконання наступної умови: < Δ, де дельта певна задана похибка. Швидкість збіжності цього алгоритму значною мірою залежить від вірного вибору початкової точки. Коли в процесі обчислень кут нахилу дотичної f’(x)перетворюється на нуль, застосування цього методу ускладнюється. Можна також показати, що у випадку дуже великих значень f ’’(x) чи кратних коренів метод Ньютона стає неефективним.

Початкове наближення слід вибирати з умови:

(1.7)

Грубо оцінити похибку для методу можна так:

, (1.8)

де М2 – найбільше за модулем значення другої похідної на інтервалі

[xn , xn +1 ].

Рис. 1.2. Метод Ньютона


2.Аналіз заданого рівняння

В даній курсовій роботі необхідно розв’язати нелінійне рівняння 5-го порядку, яке відповідно матиме п’ять коренів. Для того, щоб розв’язати це рівняння методами Ньютона та січних, необхідно визначити початкове приблизне наближення, це можна зробити за допомогою графіка цього рівняння (Рис. 2.1), побудувавши його за допомогою математичного пакета Mathcad 2001 Professional.

Рис. 2.1

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

Отже, потрібно брати такі приблизні значення початкових наближень: -1,8; 0,1; 1,8. Дійсно ці значення задовольняють необхідну умову . При значенні х0 =0,1 (f’(x0 ))2 =99,9, a добуток f’’(x0 )f(x0 )=5*10-10 , що є меншим за значення першої похідної, піднесеної до квадрату. При значеннях х0 = -1,8 ; х0 =1,8 також виконуються дані співвідношення.

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

Щоб знайти комплексні корені нелінійного рівняння розкладемо задане рівняння на найпростіші множники. В результаті отримаємо такий многочлен:, де 1,75; 1,8; 0,1 – дійсні корені рівняння. Прирівнявши перший множник до нуля , знайдемо приблизні значення комплексних коренів рівняння, яких має бути два:

;дискримінантD=0.0025-12.61= -12.6<0

;


3. Алгоритм методів

На основі теоретичного матеріалу, розглянутого в попередньому розділі роботи, було розроблено алгоритми методів.

Алгоритм розвязку нелінійного рівняння методом Ньютона за допомогою ЕОМ є досить простим і полягає в тому, що спочатку задається дане вихідне рівняння, його похідна, а також допустима похибка. Потім використовуючи вищеописану ітераційну формулу знаходять ряд значень х:

хn +1 = хn -,

де хn +1 – значення х на наступній ітерації, а хn – значення х на попередній ітерації. Ця операція повторюється до тих пір, поки не виконається умова <Δ, тобто різниця значень наступної ітерації і попередньої менше за задану похибку.

Алгоритм розвязку цього ж рівняння за методом хорд полягає в тому, щовизначаються значення функції до зміни знаку при переході від до поки кінці інтервалів xn +1 , хn не будуть мати різні знаки.

Після цього визначають f(x*) і порівнюють його з f(xn ). Надалі користуються xn +1 замість того значення, з яким воно збіглося за знаком.

Знаходимо ряд значень х до тих пір, поки не виконається умова< Δ де Δ - задана допустима похибка.

Блок-схеми даних методів наведені в додатку А.


3.1 Вибір інструментальних засобів

Для вирішення цієї задачі було обрано середовище програмування С, так як воно має ряд вагомих переваг перед іншими середовищами і мовами програмування. Зокрема такими перевагами є те, що:

- не вимагає великих затрат як апаратної частини комп’ютера так і програмної

- дозволяє досить просто реалізовувати поставлені задачі

- є дуже візуальним і наглядним що робить його зручним інструментом в користуванні

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

3.2 Вхідні та вихідні дані

Для даних методів розв’язку нелінійного рівняння вхідними даними є початкове рівняння; похідна від нього; початкові наближення х0 (1,75; -1,8; 0,1) і допустима похибка Δ, яка вводиться з клавіатури .

Вихідними даними є знайдені корені х, які задовольняють умову:

< Δ та кількість кроків для отримання розвязку із заданою похибкою.

3.3 Структура програми

Програму можна умовно розділити на чотири такі частини:

1. блок опису вхідних та вихідних даних

2. введення початкових даних

3. процедури розвязку рівняння методом хорд та Ньютона

4. виведення результатів

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

Для виведення отриманих значень х використовується функціяres.

Лістінг програми наведений в додатку В.

3.4 Інструкція користувачеві

Для завантаження програми необхідно запустити програмний файл KURSOVA.EXE. При цьому з’явиться вікно (Рис. 3.4)

Рис. 3.4.

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


4. Аналіз результатів розрахунку

Після проведення розрахунку по знаходженню коренів нелінійного рівняння за методами Ньютона та хорд отримано такі результати: рівняння має п’ять коренів, а саме три дійсних і два комплексних.

За методом хорд: х1 =-1,801176 (5 ітерацій); х2 =0,00001(3 ітерації); х3= 1,748379 (5 ітерацій).

За методом Ньютона: х1 =-1,802453 (3 ітерації); х2 =0,00001(2 ітерації); х3= 1,752384 (3 ітерації).

Комплексні корені : ,

Порівнявши отримані результати з наступними результатами,

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


Висновки

В даній курсовій роботі було проаналізовано розв’язок нелінійних рівнянь методами Ньютона та хорд. В результаті роботи було досліджено існуючі методи для розв’язання таких рівнянь, а більш детально розглянуті вищезгадані два методи. Для цих методів було складено блок-схему, а також написано програму на мові програмування С++. В результаті роботи за допомогою складеної програми було отримано певні корені заданого рівняння і порівняно їх зі значеннями коренів цього ж рівняння, але розв’язаного за допомогою спеціалізованого математичного програмного пакету Mathcad. Також було доведено, що метод Ньютона має значно вищу швидкість збіжності і для знаходження коренів потрібно значно менше ітерацій.


Література

1. Квєтний Р. Н. Методи комп’ютерних обчислень: Навчальний посібник. – Вінниця.: ВДТУ, 2001.–С. 35.

2. Вержбицький В. М. Основы численных методов.– М.: Высшая школа, 2002. –С. 43.

3. Волков Е. А. Численные методы. – М.: Наука, 1982. –С. 102.

4. Лященко М.Я., Головань М.С. Чисельні методи : Підручник. –К.: Либідь, 1996. –С. 144.

5. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. – М.: Наука, 1987. –С. 83.


Додаток А

(Алгоритм методів)

Блок-схемаметоду хорд

Блок-схемаметоду Ньютона


Додаток Б

Блок-схема програми


Додаток В

Лістінг програми

#include <math.h>

#include <stdio.h>

#include <conio.h>

double e; // похибка

// f(x)=x^5-10x+1 - функция

double f(double x)

{

return pow(x,5) - 10*x + 1;

}

// f'(x)=5x^4-10 - перша похідна f(x) для метода Ньютона

double f_der(double x)

{

return 5*pow(x,4) - 10;

}

int chords(double xa, double xb, double &x, int &iter)

// метод Хорд

{

double xlast;

x = 0;

if (f(xa)*f(xb) >= 0)

return 0;

iter = 0;

do

{

xlast = x;

x = xb - f(xb)*(xb-xa)/(f(xb)-f(xa));

if (f(x)*f(xa) > 0)

xa = x;

else

xb = x;

iter++;

}

while (fabs(x-xlast)>e);

return 1;

}

double Nuton(double x, int &iter)

// метод Ньютона

{

double xlast;

iter = 0;

do

{

xlast = x;

x = x - f(x)/f_der(x);

iter++;

}

while (fabs(x-xlast)>e);

return x;

}

int main()

// основна програма

{

double res;

int iter;

clrscr();

printf(" Enter the accuracy : ");

scanf("%f", &e);

printf("\n-------------------------Chord's method--------------\n");

for (int i = -20; i < 20; i++)

if (chords(i, i+1, res, iter))

printf("\n x= %f (%d iterations)\n", res, iter);

printf("\n-------------- Nuton's method----------------------\n");

res = Nuton(-2, iter);

printf("\n x= %f (%d iterations)\n", res, iter);

res = Nuton(0, iter);

printf("\n x= %f (%d iterations)\n", res, iter);

res = Nuton(2, iter);

printf("\n x= %f (%d iterations)\n", res, iter);

getch();

return 0;

}

Оценить/Добавить комментарий
Имя
Оценка
Комментарии:
Где скачать еще рефератов? Здесь: letsdoit777.blogspot.com
Евгений06:46:49 19 марта 2016
Кто еще хочет зарабатывать от 9000 рублей в день "Чистых Денег"? Узнайте как: business1777.blogspot.com ! Cпециально для студентов!
19:52:36 28 ноября 2015

Работы, похожие на Курсовая работа: Дослідження чисельних методів вирішення нелінійних рівнянь

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

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



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

Рейтинг@Mail.ru