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

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

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

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

Добавлен: 30.11.2023

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

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

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


b

h

n

2
    

  h      




IT f

a

xdx



i1

fi1 fi

2 f0

2 f1

2 f2

...

2 fn1

fn . (3.6)

Для удобства вычислений формулу (3.6) записывают следующим обра-

зом:


b
 


f0 fn

n1

IT f

xdx h

2

fi .

ai1 

Вид представленной формулы позволяет сделать вывод, что она может быть сформирована также исходя из других соображений, так как получае- мый с ее применением результат является средним арифметическим резуль- татов использования формул левых (3.3) и правых (3.4) прямоугольников.

Поэтому

I IL IRT 2

и указанное утверждение справедливо для случаев

равных и неравных шагов интегрирования.

Замечание. Несмотря на то, что в методе трапеций для аппроксимации подынтегральной функции используются полиномы первой степени, по сравнению с методами прямоугольников, которые используют полиномы ну- левой степени, точность метода трапеций оказывается ниже точности метода средних прямоугольников. Более высокая точность метода средних прямо- угольников объясняется «удачным» выбором узловых точек для вычисления значений подынтегральной функции.
      1. Метод Симпсона (метод парабол)


Аппроксимацию функции

fx

по методу трапеций можно

интерпретировать как замену исходной функции

fx

некоторой кусочно-


линейной функцией. Ошибка метода в данном случае определяется

«грубостью» предложенного способа аппроксимации функции. Естественно

допустить, что если исходную функцию

fx

аппроксимировать на



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

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

функция

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.


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

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