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

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

Название: Дослідження методів чисельного інтегрування
Раздел: Рефераты по информатике, программированию
Тип: курсовая работа Добавлен 09:32:59 03 февраля 2010 Похожие работы
Просмотров: 87 Комментариев: 2 Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать

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

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

Інститут автоматики, електроніки та комп’ютерних систем управління

Факультет АКСУ

Кафедра АІВТ

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

з дисципліни

«Моделювання на ЕОМ»

Дослідження методів чисельного інтегрування

Перевірив : к.т.н., доцент _________ “___” _______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 і для обчислення інтегралу методом Гауса третього, четвертого та п’ятого порядків потребує такі наступні системні параметри:

- Процесор типу Pentium-2;

-256 Мb ОЗУ;

-Операційні системи MS-Windows 98/95/ХР

-Video пам’ять 32 Мб.

Вхідними даними для програми є:

а) Межі інтегруваня: a = 0, b = 1;

б) Крок інтегрування: h = 0,1; 0,2; 0,5;

в) Підінтегральна функція :

f(x) = ;


Вихідними результатами є:

Наближене значення інтегралу отримане за допомогою Mat-hcad2001;

Результат інтегралу

I=

на проміжку [0;1]

Відносна похибка при обрахунках .

Результат виконання програми рисунок 2 буде мати такий вигляд:

Рисунок 2 – Результат виконання програми

Розв’яжемо задачу в пакеті MathCAD 2001 Professional:

Значення інтегралу І=2,680695

4 Аналіз результатів

Для початку аналізу запустимо програму. Послідовно задаємо крок обчислень: 0.1; 0.2; 0.5. Отримані результати порівняємо з розвязком 2,6806954 отриманим за допомогою пакету MathCAD 2001 Professional. Результати запишемо в таблицю:

Таблиця 1

Порядок, похибка Крок обрахунку Значення MathCAD
0,1 0,2 0,5
Значення інтегралу 3-го порядку 2,7305584726 2,6816429227 2,6864903393 2,680695457
Відносна похибка, % 0,049862943 0,0353440217 0,2161708502
Значення інтегралу 4-го порядку 2,7303869675 2,6806954572 2,6806954723 2,680695457
Відносна похибка, % 0,049691510 0,0000000030 0,0000005638
Значення інтегралу 5-го порядку 2,7303870222 2,6806955109 2,6806955116 2,680695457
Відносна похибка, % 0,049691565 0,0000020051 0,0000020303

Неточності обрахунків пов’язані з недоліками наближених методів, із так званою методичною похибкою. Але незважаючи на ці недоліки формулу Гауса називають формулою найвищої алгебраїчної точності [4].


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

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

Програма починає працювати після підключення модуля введення-виведення iostream та модуля math для підключення математичних функцій. Далі вводяться дані для інтегрування та відповідні функції.

Робота програми не вимагає великих системних ресурсів. Щоб забезпечити нормальне функціонування програми, необхідно впевнитись в існуванні стандартних бібліотек мови С ++.

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


Висновок

В цій курсовій роботі було наведено дослідження вирішення визначеного інтегралу виду

I=,

методом Гауса третього, четвертого, п’ятого порядку.

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

Отже при обрахунку даного інтегралу було отримані наступні результати:

2,73055; 2,68164; 2,68649 – для значення інтегралу 3-ого порядку кроком 0.1, 0.2, 0.3 відповідно;

2,73038; 2,6806954; 2,2606955 - для значення інтегралу 4-ого порядку кроком 0.1, 0.2, 0.3 відповідно;

2,73039; 2,680695511; 2,680695512 - для значення інтегралу 5-ого порядку кроком 0.1, 0.2, 0.3 відповідно;

Щоб впевнитись у вірності роботи програми для перевірки результатів обчислень, було використано математичний пакет MatchCAD 2001. Похибка обрахунку була оцінена як різниця між точним значенням, що отримане у Mathcard та тим, що отримане в результати роботи програми див (табл. 1).

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


Література

1. Л.М. Круподьорова, А.М. Пєтух. Технологія програмування мовою Сі. - 183 с.

2. В.С. Проценко, П.Й. Чапенко, А.В.Ставровський. Техніка програмування мовою Сі. - 212 с.

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

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

5. В.М. Дубовий, Р.Н.Квєтний. Програмування комп’ютеризованих систем управління і автоматики. - В.: ВДТУ, 1997.- 208 c.

6. Н.В. Богомолов. Практические занятия по матиматике.- Киев: Вища школа, 1979.- 472 с.

7. В.Т.Маликов, Р.Н.Кветный. Вычислительные методы и применение ЭВМ. – К.: Вища школа, 1989.- 213 с.


Додаток А

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


Додаток Б

Лістинг

#include <іоstream.h>

#include <stdio.h>

#include <math.h>

#include <conio.h>

#include <stdlib.h>

double ti[4][6]={

{-0.677459667,0.0,0.677459667},{-0.86113631,-0.33998104,0.33998104,0.86113631}, // 4 порядку

{0},

{-0.93246951,-0.66120939,-0.238619119,0.238619119,0.66120939,0.93246951 }

};

double ai[4][6]={

{(double)5/9,(double)8/9,(double)5/9},

{0.34785484,0.65214516,0.65214516,0.34785484},

{0},

{0.17132450,0.36076158,0.46791394,0.46791394,0.36076158,0.17132450}

};

double real=(2.68069545714219);

double f (double x)

{

return exp(-2*x*x+1.5);

}

double gaus(double a, double b, double h, int por)

{

int i;

double sum=0,sum1,aa,bb,x,xi;

for (x=a; x<b; x+=h)

{

aa=x; bb=x+h;

sum1=0;

for (i=0; i<por; i++)

{

xi=(aa+bb)/2+(bb-aa)*ti[por-3][i]/2;

sum1+=ai[por-3][i]*f(xi);

}

sum1*=(bb-aa)/2;

sum+=sum1;

}

return sum;

}

void compare(double _real,double _given,double h,double por,char *method)

{

double ebs;

ebs=fabs(_real-_given)/_real*100;

printf("Result recieved by method %s with step %f\n",method,por);

printf("\tI=%.10f\n",_given);

printf("\tEbs=%.10f(%%)\n",ebs);

printf("\n\n++++++++++++++++++++++++++++++++++++++++++++++++++++++++\n\n\n");

}

void main()

{

double I;

double h;

double a=0,b=1;

printf("\t Nablugene znachennja Math Cad 2001 = 2.68069545714219\n\n");

printf("\tEnter step: ");

scanf("%lf",&h);

I=gaus(a,b,h,3);

compare(real,I,h,3,"Gaus");

I=gaus(a,b,h,4);

compare(real,I,h,4,"Gaus");

I=gaus(a,b,h,6);

compare(real,I,h,5,"Gaus");

getch ();

}

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

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

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

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



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

Рейтинг@Mail.ru