Численное интегрирование разными методами

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Содержание:

Задание

Исходные данные

1. Вычисление интеграла аналитически, методом средних прямоугольников, методом трапеций, методом Симпсона

1.1 Аналитически

1.2 Метод средних прямоугольников

1.2.1 Метод средних прямоугольников при n=1

1.2.2 Метод средних прямоугольников при n=2

1.3 Метод трапеций

1.3.1 Метод трапеций при n=1

1.3.2 Метод трапеций при n=2

1.4 Метод Симпсона

1.4.1 Метод Симпсона при n =1

1.4.2 Метод Симпсона при n=2

2. Вычисление интеграла методом Гаусса

2.1 Одноточечная схема метода Гаусса

2.2 Двухточечная схема метода Гаусса

2.3 Трехточечная схема метода Гаусса

3. Сравнительный анализ точности полученных результатов

4. Вычисление интеграла

4.1 Аналитически

4.2 Метод Гаусса

4.2.1Одноточечная схема

4.2.2 Двухточечная схема

Вывод

Задание:

1. Вычислить интеграл аналитически, а затем численно методами средних прямоугольников, трапеций, Симпсона, принимая n=1 и n=2 (n-количество разбиений отрезка интегрирования)

2. Вычислить этот же интеграл методом Гаусса по одноточечной, двух точечной и трехточечной схемам интегрирования

3. Сделать сравнительный анализ точности полученных результатов

4. Вычислить интеграл аналитически, а затем численно методом Гаусса по одноточечной и двухточечной схемам. Сравнить результаты

Исходные данные:

a = 0

b =

f(x) = x cos(x)

f(x,y) = 2x² + y²



1. Вычисление интерграла аналитически, методом средних прямоугольников, методом трапеций, методом Симпсона.

1.1 Аналитически

Для вычисления данного интеграла используем формулу интегрирования по частям:

.

Пусть:

u = x

dv = cos (x) dx

Тогда:

du = dx

v = sin(x)

Следовательно:

Откуда получаем:

1.2 Метод средних прямоугольников

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

По методу средних прямоугольников на интервале [x_i,x_i+h]имеем

1.2.1 Метод средних прямоугольников для n=1

При n=1:

x_i=a

x_i+h=b

=

h=(b-a)/n

Исходя из этого получаем уравнение:

= ((b-a)/n) * ( cos())

Подставляя в полученное уравнение исходные данные получаем:

=(( - 0)/1)*((0+(( - 0)/1)/2)*cos(0+(( - 0)/1)/2)) = 0

= 0

Погрешность:

R==-2,584

1.2.2 Метод средних прямоугольников для n=2

В данном случае исходный отрезок интегрирования разбивается на 2 отрезка интегрирования: (0;) и (;).

В связи с этим:

X₁=0

X₂=

=

h=(b-a)/n

По полученным данным получаем уравнение:

= ((b-a)/2)*(((x₁+(b-a)/4)*(cos(x₁+(b-a)/4))+((x₂+(b-a)/4)*

*(cos(x₂+(b-a)/4))

Подставляем в полученное уравнение начальные данные и получаем:

= ((-0)/2)*(((0+( -0)/4)*(cos(0+( -0)/4))+(( +( -

0)/4)* *(cos( + ( -0)/4)) = -1,745

= -1,745

Погрешность:

R==-0,636



1.3 Метод трапеций

Метод трапеций -- метод численного интегрирования функции одной переменной, заключающийся в замене на каждом элементарном отрезке подинтегральной функции на многочлен первой степени, то есть линейную функцию.

Если отрезок (x_i;x_i+h) является элементарным и не подвергается дальнейшему разбиению, значение интеграла можно найти по формуле

=

1.3.1 Метод трапеций при n=1

На отрезке (a;b) заменяем f(x) полином первой степени.

В этом случае

X_i=a

X_i+h=b

h=(b-a)/n

и уравнение имеет вид:

= ((b-a)/n)*[(a cos(a))+(b cos(b))]/2

Подставим исходные данные и получим:

= ((-0)/1)*[(0 * cos(0))+( *cos(]/2 = -4,935

= -4,935

Погрешность:

R₀=(-h³/12)*f'''(x_i)= 8,117

1.3.2 Метод трапеций при n=2

Отрезке (a;b) разбиваем на 2 отрезка (a;c) и (c;b) и замеим f(x) полином первой степени на обоих отрезках.

В этом случае

X₁=a

X₂ =c = (a+b)/2

X₃ =b

h=(b-a)/n

и уравнение имеет вид:

= ((b-a)/n)*(([(a cos(a))+(c cos(c))]/2)+([( c cos(c))+(b

cos(b))]/2))

Подставим исходные данные и получим:

= ((-0)/1)*(([(0 * cos(0))+( *cos(]/2)+ ([( *

cos( ))+( * cos())]/2)) = -2,467

= -2,467

Погрешность:

R₀=(-h³/12)*f'''(x_i)=0,646

численный интегрирование симпсон гаусс



1.4 Метод Симпсона

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

Формула Симпсона:

= h[f(x₀)+4f(x₁)+f(x₂)]/3

1.4.1 Метод Симпсона при n =1

Подынтегральную функцию f(x) заменим полином второй степени P₂(x) - параболой, проходящей через равноотстоящие точки x₀,x₁,x₂.

В данном случае

X₀=a

X₁=c = (a+b)/2

X₂ = b

H = (b-a)/2n

Полученное уравнение имеет вид:

= (b-a)/2n [(a cos(a))+(c cos(c))+(b cos (b))]/3

Подставляем исходные данные и получаем:

= (-0)/2*1 [(0 cos(0))+( cos( ))+( cos ( ))]/3 = -

1,645

=-1,645



1.4.2 Метод Симпсона при n=2

В данном случае исходный отрезок [a;b] разбиваем на 2: [a;c] и [c;b].

Подынтегральную функцию f(x) заменим полином второй степени P₂(x) - параболой, проходящей через равноотстоящие точки x₀,x₁,x₂ - на первом отрезке и x₂,x₃,x₄ - на втором отрезке.

В данном случае

X₀=a

X₁=d = (a+c)/2

X₂ = с= (a+b)/2

X₃ = e = (c+b)/2

X₄ =b

h=(b-a)/2n

Полученное уравнение имеет вид:

= (b-a)/2n( ([(a cos(a))+(d cos(d))+(c cos (c))]/3)+([(c

cos(c)) + 4(e cos(e)) + (b cos(b))]/3))

Подставляем исходные данные и получаем:

= (-0)/2*1(( [(0 cos(0))+( cos( )) + ( cos (

))]/3) + ([( cos ( + 4( cos( +( cos ( )) = -1,986

=-1,986



2. Вычисление интеграла методом Гаусса

Метод Гаусса -- метод численного интегрирования, позволяющий повысить алгебраический порядок точности методов на основе интерполяционных формул путём специального выбора узлов интегрирования без увеличения числа используемых значений подынтегральной функции. Метод Гаусса позволяет достичь максимальной для данного числа узлов интегрирования алгебраической точности.

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

Формула:

=

,

Где t- координаты узлов

C_k - весовые коэффициенты

2.1 Одноточечная схема метода Гаусса.

При наличии 1-го узла

t=0

c_k = 2



Формула имеет вид:

=

Подставим в данную формулу исходные данные и получим:

= =0

=0

2.2 Двухточечная схема метода Гаусса

При наличии 2-х узлов

t₁=-1/

t₂ = 1/

c₁ = 1

c₂ = 1

Формула имеет вид:

=

Подставим в данную формулу исходные данные и получим:

=

= -2,244

=-2,244

2.3 Трехточечная схема метода Гаусса

При наличии 3-х узлов

t₁= -

t₂ = 0

t₃ =

c₁ = 5/9

c₂ = 8/9

c₃ =5/9

Формула имеет вид:

=

Подставим в данную формулу исходные данные и получим:

=

= -1,992

=-1,992

3. Сравнительный анализ точности полученных результатов

Метод	Точное значение интеграла =	Погрешность
Аналитически	-2	-
Средних прямоугольников, n=1	0	-2,584
Средних прямоугольников, n=2	-1,745	-0,636
Трапеций, n=1	-4,935	8,117
Трапеций, n=2	-2,467	0,646
Симпсона, n=1	-1,645	-0,355
Симпсона, n=2	-1,986	0,014
Гаусса, n=1	0	-2
Гаусса, n=2	-2,244	0,244
Гаусса, n=3	-1,992	-0,008



4. Вычисление интеграла

4.1 Аналитически

Вычислим внутренний интеграл

i=.

Интеграл суммы равен сумме интегралов, следовательно:

=+ = 2*+x = (2/3+) =

20,639 + 3,14

Подставим полученное значение во внешний интеграл и вычислим его.

= + =20,639y +

3,14 = 64,807 + 32, 404 = 97,211

4.2 Метод Гаусса.

Двойной интеграл вычисляется методом Гаусса аналогично одномерному случаю.

4.2.1Одноточечная схема.

При наличии 1-го узла

=0

c_i,j = 2

Формула имеет вид:

=

Подставим в данную формулу исходные данные и получим:

=

= 46,472

4.2.2 Двухточечная схема

При наличии 2-х узлов

t₁=-1/

t₂ = 1/

c₁ = 1

c₂ = 1

=

= 31



Вывод:

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

Размещено на Allbest.ru