Дослідження методів чисельного інтегрування (стр. 1 из 2)
Міністерство освіти і науки України
Вінницький національний технічний університет
Інститут автоматики, електроніки та комп’ютерних систем управління
Факультет АКСУ
Кафедра АІВТ
з дисципліни
Дослідження методів чисельного інтегрування
Перевірив : к.т.н., доцент _________ “___” _______2008 Кабачій В.В.
Виконав: ст. Гр. 1АМ-06 _________ “___” _______2008 Ігнатенко В.О.
Вінниця 2008
Зміст
Анотація........................................................................................................... 3
Вступ................................................................................................................ 4
1 Теоретичні відомості.................................................................................... 6
2 Алгоритм Методу......................................................................................... 7
3 Загальні відомості та функціональне призначення................................... 11
4 Аналіз результатів...................................................................................... 13
5 Інструкція користувачеві............................................................................ 14
Висновки........................................................................................................ 15
Література...................................................................................................... 16
Додаток А. Блок-схема програми................................................................. 17
Додаток Б. Лістинг програми....................................................................... 18
Анотація
В даній курсовій роботі проведено дослідження методу чисельного інтегрування. Дослідження проводилося за допомогою Методу Гауса, при обчисленні інтегралу третього, четвертого та п’ятого порядків.
Дана програма розв’язує інтеграл методом Гауса, знаходить похибку і виводить результати порівнюючи їх із розв’язком, отриманим у прикладній програмі MathCAD 2001 Professional.
Вступ
Використання сучасних персональних комп’ютерів охоплює майже всі сфери людської діяльності і поступово підпорядковує собі всі інформаційні технології. В останній час головний аспект застосування комп’ютерів зсувається з галузей, де комп’ютер був самостійним інструментом – обчислення при наукових дослідженнях та проектуванні, зберігання та обробка статистичної інформації тощо, у бік галузей, де комп’ютер розглядається як складова частина більш масштабних систем – систем автоматичного, та автоматизованого управління, інформаційно – вимірювальних систем, систем мультимедії тощо. Сьогодні практично всі фахівці у цих галузях повинні вільно володіти комп’ютерними технологіями.
За декілька десятиліть світової комп’ютерної індустрії було створено безліч різноманітних мов програмування, проте переважна їх більшість або не дуже вдало копійована, або створена для деякого вузького спеціалізованого застосування. Внаслідок своєї еволюції отримали визнання невелика кількість мов програмування, серед яких і добре відома мова Сі. Створюються системи програмування С, до складу якої входять бібліотеки, з широким набором різноманітних функцій та інтегровані середовища розробки (IDE – Integrated Development Environment).
Обчислювальна математика заснована на чисельних методах, придатних до застосування при розрахунках на ЕОМ. Сучасні ЕОМ дозволили дослідникам значно підвищити ефективність математичного моделювання складних задач науки і техніки. На сьогодні методи дослідження проникають практично в усі сфери людської діяльності, а математичні моделі стають засобами пізнання.
Значення математичних моделей неперервно зростає у зв'язку з тенденціями до оптимізації технічних пристроїв і технологічних схем планування експерименту. Реалізація моделей на ЕОМ здійснюється за допомогою різноманітних методів обчислювальної математики, яка неперервно вдосконалюється [3].
1 Теоретичні відомості
Визначений інтеграл – це площа, обмежена графіком функції y = f(x), віссю Ох і прямими f(a) і f(b). Якщо він перетинає вісь Ох, то інтеграл чисельно рівний алгебраїчній сумі площ, що знаходяться по кожну сторону вісі Ох.
Обчислення визначених інтегралів має дуже широке застосування. Так наприклад їх обрахунок необхідний в задачах, що пов’язані з оцінкою якості, аналізом інформаційно-вимірювальної техніки, радіоелектроніки.
В основу інтегрування покладено наближене обчислення площини під кривою, яка описується підінтегральною функцією інтеграла (1.1):
I =
f(x)dx (1.1)Для обчислення визначеного інтеграла від функції f(х) в тому випадку, коли можна знайти відповідний невизначений інтеграл F(x), є формула Ньютона-Лейбніца:
I =
f(x)dx = F(x) 
= F(b) – F(a), (1.2)Загальний підхід до розв’язування цієї задачі такий : визначений інтеграл I являє собою площину, обмежену кривою f(х), віссю x та прямими x = a, x =b, відрізок від a до b розбивають на множину менших відрізків, знаходять наближено площу кожної площини Si, яку отримують за таким розбиванням, значення інтеграла I знаходять як суму площ площин Sі, тобто
I =
Si.При цьому використовують два способи розбивання початкового відрізка на менші:
а) розбивання відрізка проводиться раніше, до того ж завжди відрізок вибирають рівними (метод прямокутників, трапецій, Сімпсона);
б) місцезнаходження та довжина відрізків визначаються аналізом, до того ж спочатку ставиться за мету досягти найбільшої точності з заданим числом відрізків, а потім відповідно з цим визначають їхні межі (методи Гаусса, Ньютона - Котеса, Чебишева) [6].
2 Алгоритм Методу
Формула Гаусса називається формулою найвищої алгебраїчної точності. Для формули розрахунку найвища точність може бути досягнута для поліномів степеня (2n-1), які визначаються 2n постійними ti і Ai (і=1,2,...,n).
Суть методу полягає у визначенні коефіцієнтів Ai і абсцис точок ti. Для знаходження цих постійних розглянемо виконання формули розрахунку для функцій вигляду (2.1):
f(t)=tk, k=0,1,…,2n-1. (2.1)
Враховуючи (2.2)
tkdt= 
, (2.2)отримаємо систему рівнянь (2.3) [4] :
Ai=2; Aiti=0; Aiti2=1; (2.3) Aiti2n-2= 
; Aiti2n-1=0;Ця система нелінійна, і її звичайне розв'язання пов'язане із значними обчислювальними труднощами. Але якщо використовувати систему для поліномів вигляду (2.4):
f(t)=tkPn(t), k = 0,1,…,n-1, (2.4)
де Pn(t) - поліном Лежандра, тоді її можна звести до лінійної відносно коефіцієнтів Ai з заданими точками ti. Оскільки степені поліномів в співвідношенні не перевищують 2n-1, повинна виконуватися система і дана формула приймає вигляд (2.5) :
tkPn(t)dt= AitikPn(ti) (2.5)В результаті властивості ортогональності ліва частина виразу дорівнює 0, тоді формула буде (2.6) :
AitikPn(ti)=0, (2.6) що завжди забезпечується при будь-яких значеннях Ai в точках ti, які відповідають кореням відповідних поліномів Лежандра.Підставляючи ці значення ti в систему і враховуючи перші n рівнянь, можна визначити коефіцієнти Ai.
Формула розрахунку, де ti - нулі полінома Лагранжа Pn(t), а Ai, i=1,2,...,n визначаються із системи, називається формулою Гаусса.
Значення Ai, ti для різних n наведені в довідниках.
Для довільного Інтервалу (а,b) формула для методу Гаусса приймає вигляд (2.7) :
I=
Aif(xi), (2.7)де xi обчислюється за формулою (2.8) :
xi=
+ ti. (2.8)Оцінка похибки формули Гаусса з n вузлами визначається із співвідношення (2.9) :
(2.9)де M2n- максимальне значення 2n похідної на ділянці (а,b). Враховуючи наведені вище формули розробимо алгоритм методу (рис 1) [4].
Рисунок 1 - Алгоритм методу Гауса
3 Загальні відомості та функціональні призначення
Дана програма обчислює інтеграл виду
I=
,методом Гауса третього, четвертого та п’ятого порядків. При зменшенні кроку інтегрування похибка обчислювань зменшиться.
Програма досить зручна та проста у користуванні. Результати обчислень виводяться на екран монітора разом із розв’язком, отриманим в математичному пакеті MathCAD 2001 Professional та похибкою обчислень .
Програма призначена для обчислення тільки одного інтегралу, що значно зменшує сферу її використання. Але змінивши програмний код можна досягти обчислення і кількох інтегралів.
Дана програма досить швидко проводить обрахунок і виводить відповідь безпосередньо на екран монітору.
Програма складена на Borland C++ 5.02 і для обчислення інтегралу методом Гауса третього, четвертого та п’ятого порядків потребує такі наступні системні параметри: