ЛАБОРАТОРНАЯ РАБОТА ’5
Работа добавлена на сайт samzan.net:
ЛАБОРАТОРНАЯ РАБОТА №5.
ВЫЧИСЛЕНИЕ ИНТЕГРАЛОВ МЕТОДАМИ ПРЯМОУГОЛЬНИКОВ, ТРАПЕЦИЙ, СИПСОНА.
I. ЦЕЛЬ РАБОТЫ.
1. Изучение основных определений и положений теории численного интегрирования.
2. Изучение основных квадратурных формул численного интегрирования.
3. Разработка численного алгоритма и программ для вычисления на ЭВМ интегралов по квадратурным формулам прямоугольников, трапеций, Симпсона. Оценка погрешности этих формул по правилу Рунге.
II. КРАТКИЕ ТЕОРЕТИЧЕСКИЕ ПОЛОЖЕНИЯ.
1. Основные определения. Определенным интегралом от функции f(x) на отрезке [a,b] называется предел интегральной суммы:
(2.1)
где, n - количество элементарных отрезков [xi-xi-1], i=1,...,n; на которые разбивается отрезок интегрирования [a,b], xi=(xi-xi-1) - длина i-ого отрезка, i - точка на отрезке [xi-1,xi].
Когда функция f(x) задана аналитически в виде формулы и интеграл удается свести к табличному, то интеграл (2.1) вычисляется с помощью таблиц неопределенных интегралов и формулы Ньютона-Лейбница, например:
(2.2)
где F(x) - первообразная, т.е. F(x)=f(x).
Однако на практике обычно интеграл (2.2) не сводится известными приемами к табличному интегралу, даже тогда, когда под интегральная функция задана аналитически, не говоря уже о тех случаях, когда значения под интегральной f(x) заданы в виде таблицы. В этом случае используют численные методы
2. Основные квадратурные формулы. Для вычисления определенных интегралов используется приближенное соотношение:
(2.3)
которое называется квадратурной формулой с узлами i и весами qi.
В формуле (2.3) интеграл приближенно заменяется конечной суммой, члены которой представляют произведение значений функций в некоторых узлах на некоторую величину. Наиболее часто используются следующие квадратурные формулы:
а) формула прямоугольников:
(2.4)
где
Для всего отрезка [a,b] имеем:
(2.5)
Погрешность формулы (2.5), полученная с помощью ряда Тейлора равна:
. (2.6)
- максимальное значение второй производной на отрезке [a,b].
б) формула трапеции:
(2.7)
где fi=f(xi), i=1,...n.
Для всего отрезка имеем:
(2.8)
при этом погрешность равна:
.
в) формула Симпсона (формула парабол):
(2.9)
Для всего отрезка:
(2.10)
где
Погрешность формулы Симпсона равна:
(2.11)
- Метод ячеек для вычисления кратных интегралов. Пусть требуется вычислить двукратный интеграл в области G(axb, cyd):
С помощью узлов xi (i=0,1,..n) и yj (j=0,1,...,m) и прямых, проходящих через эти узлы: x=xi и y=yj, разобьем область G на (nm) прямоугольных ячеек, имеющих площадь:
Выбираем в этой ячейке центральную точку:
Будем считать, что интеграл для каждой ячейке приближенно равен:
(2.12)
Суммируя по всем ячейкам имеем:
(2.13)
при этом погрешность, когда все ячейки имеют одинаковую площадь будет равна
; (2.14)
где S - площадь области G, m и n - количество узлов по координатам x,y; , - максимальное значение вторых частных производных по соответствующим координатам.
4. Правило Рунге практической оценки погрешности и уточнению по Ричардсону. Пусть I - точное значение интеграла, Ih - значение интеграла, вычисленное по квадратурной формуле с шагом h, а Ih/2- значение того же интеграла, вычисленное для шага h/2.
Можем записать:
(2.15)
где c - константа.
Величина chk- называется главной частью погрешности квадратурной формулы с порядком точности k по шагу h. Остальная часть погрешности обозначена как 0(hk+1) и имеет порядок k+1.
Вычитая из первого уравнения (2.15) второе получаем соотношение, которое с точностью порядка 0(hk+1) позволяет вычислить значение главной части погрешности:
(2.16)
Данная формула называется практической оценкой погрешности по правилу Рунге:
Подставляя (2.16) в первую формулу (2.15) получаем формулу для уточнения значение интеграла по Ричардсону:
(2.17)
Для формул прямоугольников, трапеций и ячеек имеем k=2, для формул Симпсона - k=4.
III. ЗАДАНИЕ.
1.Написать соотношения для приближенного вычисления интеграла для функции, взятой из таблицы заданий с использованием заданной квадратурной формулы.
2.Определить величину шага исходя из заданной точности.
3.Для вычисления интеграла применить уточнение по Ричардсону.
4.Написать программу и рассчитать на ЭВМ интеграл от заданной функции.
IV. ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ.
Задание. Вычислить интеграл методом прямоугольников с точностью не ниже 10-4.
1. Для приближенного вычисления интеграла от под интегральной функции f(x)=sin(sin(x)) используем квадратурную формулу прямоугольников (2.5):
(4.1)
2. Определяем число узлов интегрирования. Для этого с помощью соотношения (2.6) вначале выбираем промежуточный шаг:
(4.2)
Далее оцениваем величину
(4.3)
Знаем все значения: , , , поэтому согласно (4.2) имеем
(4.4)
Определяем число узлов для шага hp=3.5*10-2:
(4.5)
3. Определяем уточненное значение шага для выбранного числа узлов N=30:
h=(b-a)/N=1/30. (4.6)
4. Для уточнения квадратурной формулы (4.1) используем метод Ричардсона. Согласно (2.17) имеем:
(4.7)
где Ih/2 - значение интеграла, вычисленное по формуле (4.1) для шага h/2.
- Пример текста программ на Mathcad и на Delphy (в консольном режиме) для приближенного вычисления интеграла по формуле прямоугольников.
program lab5;
{Вычисление интегралов по формуле прямоугольников}
{a -нижний, b- верхний пределы интегрирования}
{r-точность вычисления}
{fp - максимальное значение второй производной функции}
{Ih - значение интеграла с шагом h}
{Ih2 - значение интеграла c шагом h/2}
{Ihh - уточненное значение интеграла по Ричардсону}
var a,b,yp,fp,r,Ih,Ih2,Ihh,hp,h,y,Iz,s : real;
var x : array[0..1000] of real;
var i,j,n : integer;
begin
writeln(‘Введите значение нижнего предела интегрирования a’);
readln (a);
writeln(‘Введите значение верхнего предела интегрирования b’);
readln (b);
writeln(‘Введите максимальное значение второй производной fp’);
readln (fp);
writeln(‘Введите значение погрешности r’);
readln (r);
hp:=sqrt(24.*r/((b-a)*fp)); вычисление по (4.4)
n:=round((b-a)/hp)+1; вычисление по (4.5)
h:=(b-a)/n; вычисление по (4.6)
for j:=1 to 2 do begin
if j=1 then hp:=h else hp:=h/2;
n:=round((b-a)/hp);
hp:=(b-a)/n;
s:=0;
x[0]:=a;
for i:=1 to n do begin
x[i]:=x[0]+i*hp;
y:=(x[i]+x[i-1])/2; вычисление i по (4.1)
s:=s+sin(sin(y)); вычисление суммы в (4.1)
end;
Iz:=hp*s; вычисление по (4.1)
if j=1 then Ih:=Iz else Ih2:=Iz;
end;
Ihh:=(4*Ih2-Ih)/3.; вычисление по (4.7)
writeln(‘Шаг h=’,h,‘ Знач. интеграла Ih=’,Ih);
writeln(‘Шаг h/2=’,hp,’ Знач. интеграла Ih2=’,Ih2);
writeln(Уточн. по Ричардсону знач. интеграла Ihh=’,Ihh);
end.
6. Заполняем таблицу:
|
Ih |
0,430635 |
|
Ih/2 |
0,430614 |
|
Ihh |
0,430606 |
V. СОДЕРЖАНИЕ ОТЧЕТА.
1. Название лабораторной работы.
2. Индивидуальное задание.
3. Теоретическая часть.
4. Текст программы.
5. Результаты расчета.
Пункты 1-4 должны быть оформлены до начала лабораторной работы.
VI. КОНТРОЛЬНЫЕ ВОПРОСЫ.
1. Понятие определенного интеграла.
2. Определение квадратурной формулы.
3. Обусловленность задачи численного интегрирования.
4. Формулы прямоугольников, трапеций, Симпсона.
5. Погрешности основных квадратурных формул.
6. Формула численного интегрирования с помощью сплайнов.
7. Метод ячеек.
8. Погрешность метода ячеек.
9. Правило Рунге.
10. Уточнение по Ричардсону.
11. Метод Монте-Карло.
12. Погрешность метода Монте-Карло.
VII. ТАБЛИЦА ИНДИВИДУАЛЬНЫХ ЗАДАНИЙ.
|
№ |
Определенный интеграл |
Метод |
Точность метода |
|
1 |
трапеций |
10-2 |
|
|
2 |
прямоугольников |
10-2 |
|
|
3 |
трапеций |
10-2 |
|
|
4 |
прямоугольников |
10-2 |
|
|
5 |
трапеций |
10-2 |
|
|
6 |
прямоугольников |
10-2 |
|
|
7 |
ячеек |
10-2 |
|
|
8 |
Симпсона |
10-4 |
|
|
9 |
прямоугольников |
10-4 |
|
|
10 |
Симпсона |
10-2 |
|
|
11 |
ячеек |
10-2 |
|
|
12 |
ячеек |
10-2 |
|
|
13 |
трапеций |
10-2 |
|
|
14 |
прямоугольников |
10-2 |
|
|
15 |
трапеций |
10-2 |
|
|
16 |
Симпсона |
10-4 |
|
|
17 |
прямоугольников |
10-3 |
|
|
18 |
ячеек |
10-2 |
|
|
19 |
ячеек |
10-2 |
|
|
20 |
трапеций |
10-4 |
2. тема предусматривающая осуществление для большей чисти электрооборудования двух видов ремонта- текущего и
3. Лабораторная работа 1 Тема- Изучение возможностей операционной системы Windows Цель работы
4. Реферат- Тоталитарный режим
5. го по 1903 pp. Надвиробництво продукції погіршення її збуту зумовили скорочення або повну зупинку діяльності б
6. Опять Изавель умышляет восхитить виноградник Навуфеев и умышляет изгнать в горы святого Илию
7. Зарождение криптографии.html
8. Банкротство юридического лица как способ его ликвидации
9. При испытании измерительной системы дифманометр вторичный прибор в нормальных условиях эксплуатации при
10. Информационная система сети салонов красоты
11. Финансовый рынок
12. Реферат- Административное принуждение, его формы
13. Понятие состава преступления- совокупность предусмотренных уголовным законом объективных и субъективны
14. Характеристика промышленности
15. специальную В основу периодизации истории медицины положена- гражданская периодизация статисти
16. Биология белого медведя.html
17. Первобытнообщинный строй на территории Северного Кавказа в т
18. ЛАБОРАТОРНА РОБОТА 2 ВИМІРЮВАННЯ ТЕМПЕРАТУР І ПЕРЕБІЛЬШЕНЬ ТЕМПЕРАТУР Мета- засвоїти методи ви
19. 1 Відрахування на пенсійне забезпечення
20. тематической экспедицией Главтюменьгеологии и утверждены ГКЗ СССР в 1985 году