Файл: Методические указания по выполнению курсовой работы по дисциплине (модулю) Введение в математический анализ.docx

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.11.2023

Просмотров: 141

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.


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

которой на частичных отрезках интегрирования заменяется исходная

функция

fx, выбрана парабола.


Пояснение. Более высокая точность вычисления интегралов обеспечивается при использовании параболической аппроксимации (полиномом второй

степени)

y ax2 bx c

по трем соседним точкам. Для нахождения

коэффициентов a, bи c полинома, проходящего через точки x0 , y0 ,

x1 , y1 , x2 , y2 , требуется найти решение следующей системы линейных уравнений

y0 ax2 bx

  • c,

0 0


1
y1 ax2bx1 c,

y2 ax2 bx c,

2 2

относительно неизвестных a, b, c.

Разделим отрезок интегрирования a,b

на четное число nравных от-

резков с шагами h. На каждом отрезке xi1 , xi1

длиной

2h, содержащем


три узла xi1, xi, xi1 , заменим подынтегральную функцию полиномом вто- рой степени ix ax2 bx c. Пример представлен на рис. 5.



Рис. 5. Пример замены функции

y fx параболой.


Для рассматриваемого примера значения коэффициентов a, bи cмо- гут быть вычислены следующим образом:

ah2 bh c

fi1 ,



cf,

i



ah2bhc

fi1 ,


откуда c

fi,

2ah2 2 fi

fi1

fi1 , a

fi1 2 fi

2h2

fi1 , а b

fi1

2h

fi1 .

В результате рассматриваемый полином второй степени примет вид

fx i

x

fi

fi1

2h

fi1 x

fi1 2 fi

2h2

fi1 x2 .

Интегрируя приведенное выражение на отрезке xi1 , xi1 , получим


xi1

xi1

f f

f 2 f f 2 h

ixdx

fii1 i1 xi1 i i1 x

dx

fi1 4 fi fi1 .

(3.7)

xi1

xi1 2h

2h2 3

Приближенное значение интеграла на исходном отрезке интегрирова-

ния a,b

получим суммированием частичных интегралов (3.7) по всем от-

резкам xi1 , xi1 :


b
  

hn21

 

I f

a

xdx



3 i0

f2i

4 f2i1

f2i2

(3.8)

h f

3 0

4 f1

2 f2

4 f3

2 f4

... 2 f
n2

  • 4 f

n1

fn.


Полученное соотношение называется формулой Симпсона или

формулой парабол, в соответствии с которой искомый определенный


интеграл вычисляется как суммарная площадь параболических сегментов на всех частичных отрезках интегрирования xi1 , xi1 .

Если подынтегральную функцию

fx

заменять полиномом второй

степени на отрезках xi1, xi,

i 1, n, с привлечением дополнительных точек

x xi xi1 ,

i 2

i 1, n

середин данных отрезков, то число отрезков разбиения

nможет быть любым (не обязательно четным), а формула (4.8) будет иметь вид:

b hn





f xdx 6 fi1 4 fi

fi

a i1

h f 4 f 2 f 4 f


... 2 f

4 f f

(3.9)


6 0 1 1 2

n1 n n


Напомним, что

fi

fxi,

fi1

fxi1 ,

fi

fxi, i 1, n.

Замечание 1. Формула (3.8) может быть использована для вычисления интегралов от функций, заданных как аналитическим выражением, так и таблично, тогда как формула (3.9) применима только в тех случаях, когда подынтегральная функция задана аналитически.

Замечание 2. Рассмотренные формулы, используемые для приближенного вычисления интегралов, называются квадратурными формулами. Нетрудно заметить, что все рассмотренные выше формулы имеют следующую структуру:

b n

fxdx Aifi,

a i0

где

fi значения подынтегральной функции в узловых точках

xi, а

Ai

весовые коэффициенты.
    1. 1   ...   4   5   6   7   8   9   10   11   ...   15