Студопедия

КАТЕГОРИИ:


Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748)

Вступ. Загальні відомості про чисельні інтегрування




Лекція 14-15.Чисельне інтегрування функцій

Приклад роботи програми

Для демонстрації програми візьмемо, диференційне рівняння f =x+y і початкові значення: х0=1, х=2, у0=0 і n=10.

Вирішимо дане рівняння трьома методами: Рунге-Кутта, Адамса, Крилова. Розв’язання можна побачити при вході в четвертий пункт меню.

При розв’язанні цього рівняння на екрані виведуться результати трьома методами.

 

    Метод
х Адамса Крилова Рунге-Кута
  1,0 0.110530417 0.110780417 0.110341666
  1,1 0.134999983 0.243803443 0.242805142
  1,2 0.310937496 0.401416230 0.399716997
  1,3 0.504006507 0.586213602 0.583648486
  1,4 0.728881703 0.801064221 0.797441291
  1.5 0.989889204 1.049139457 1.044235940
  1.6 1.288596224 1.333945346 1.327503276
  1.7 1.629020812 1.649357884 1.651079157
  1.8 2.015833021 2.029662071 2.019202866
  1.9 2.453864896 2.449595121 2.436559545
  2.0 2.948468326 2.924394106 2.908327041

 

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

 

Список літератури:

1. Авторський колектив: В.М Глушков, И.Н. Молчанов, А.А. Стогний, Б.Н. Брусникин, Г.И. Визнюк, Е.Г.Геец, В.Ф. Гук, А.А. Дородницына й ін. Програмне забезпечення ЕОМ «СВІТ-1» і «МИР-2» т.2, 1976.

 

2. В.Е. Краскевич, К.Х. Зеленський, В.Н. Грегко Чисельні методи в інженерних дослідженнях. –київ: Вища школа, 1986.

 

3. Б.П. Димидович, И.А.Марон, Шувалова Э.З. Чисельні методи аналізу.


Визначений інтеграл з межами інтегрування a і b можна трактувати як площу фігури обмеженої ординатами в точках a і b, віссю абсциси x і графіком підінтегральної функції f(x) (мал.1).

a

S = ò f(x)dx

b

Y

f(x)

 

S

x

a b

Мал. 1.

Визначений інтеграл, у якому відома його первісна F(x), обчислюється по формулі Ньютона – Лейбніца:

S = F(b) – F(a).

Тому досить обчислити значення функції F(x) в точках a і b.

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

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

Використані на практиці методи чисельного інтегрування можна згрупувати у залежності від методу апроксимації підінтегральної функції.

Такі методи, як метод Ньютона - Котеса засновані на поліноміальній апроксимації. Методи цього класу відрізняються один від одного ступенем полінома, від якого, у свою чергу, залежить кількість вузлів, де необхідно обчислити функції f(x). Алгоритми методів прості і легко піддаються програмній реалізації.

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

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

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

Незалежно від обраного методу в процесі чисельного інтегрування необхідно обчислити наближене значення інтегралу і оцінити похибку R. Похибка буде меншою при збільшенні кількості відрізків, на які розбивається інтервал інтегрування [ a, b ], за рахунок більш точної апроксимації функції, однак при цьому буде зростати похибка за рахунок підсумовування часткових інтегралів, і остання похибка починаючи з визначеного значення m стає переважною. Ця обставина повинна застерегти від вибору досить великого числа m і привести до необхідності розробки способу оцінки похибки R вибору способу оцінки похибки R обраного методу.




Поделиться с друзьями:


Дата добавления: 2014-01-07; Просмотров: 515; Нарушение авторских прав?; Мы поможем в написании вашей работы!


Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет



studopedia.su - Студопедия (2013 - 2024) год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! Последнее добавление




Генерация страницы за: 0.009 сек.