Приближенное вычисление определенного интеграла при помощи квадратурной формулы Чебышева
Курсовой проект - Математика и статистика
Другие курсовые по предмету Математика и статистика
ен сумме площадей прямоугольных трапеций, где основание трапеции какая-либо малая величина (точность), и сумма площадей прямоугольников, где основание прямоугольника какая-либо малая величина (точность), а высота определяется по точке пересечения верхнего основания прямоугольника, которое график функции должен пересекать в середине. Соответственно получаем формулы площадей
для метода трапеций:
,
для метода средних прямоугольников:
.
1.4. Общая формула Симпсона (параболическая формула)
Пусть n=2m есть четное число и yi=f(xi) (i=0,1,2...n) - значения функции y=f(x) для равноотстоящих точек а=x0,x1, ... ,xn=b с шагом
Применив формулу Симпсона к каждому удвоенному промежутку [x0,x2], [x2,x4] ... [x2m-2,x2m] длины 2h и введя обозначения
1=y1+y2+ ... +y2m-1
2=y2+y4+ ... +y2m
получим обобщенную формулу Симпсона:
Остаточный член формулы Симпсона в общем виде:
где k I (x2к-2,x2к)
1.5. Квадратурная формула Чебышева
Рассмотрим квадратурную формулу вида:
функцию f(x) будем исать в виде когда f(x) многочлен вида f(x)=ao+a1x+...+anxn . Проинтегрировав, преобразовав и подставив значения многочлена в узлах
f(x1)=a0+a1x1+a2x12+a3x13+...+anx1n
f(x2)=a0+a1x2+a2x22+a3x23+...+anx2n
f(x3)=a0+a1x3+a2x32+a3x33+...+anx3n
. . . . . . . . . . . . . . . .
f(xn)=a0+a1xn+a2xn2+a3xn3+...+anxnn
получим формулу Чебышева.
Значения х1,х2,..,хn для различных n приведены в таблице 3.
Таблица 3 Значения х1,х2,..,хn для различных n.
n I ti n i ti 21;2 0,577350 61;6 0,86624731;3 0,7071072;5 0,4225192 03;4 0,26663541;4 0,794654 71;7 0,8838622;3 0,1875922;6 0,52965751;5 0,8324983;5 0,3219122;4 0,3745414 03 02. Решение контрольного примера
где a=0 ; b= ; при n=5;
f(x) = sin(x);
i xi yi1 0,131489 0,1311182 0,490985 0,4714943 0,785 0,7068254 0,5090150,4873175 0,8685110,763367
x1= /4+/4*t1=/4+/4(-0,832498)=0,131489
x2= /4+/4*t2=/4+/4(-0,374341)=0,490985
x3= /4+/4*t3=/4+/4*0=0,785
x4=1- x2=1-0,490985 = 0,509015
x5=1- x1=1-0,131489=0,868511
y1=sin(x1) = sin(0,131489)=0,131118
y2=sin(x2) = sin(0,490985)=0,471494
y3=sin(x3) = sin(0,785)=0,706825
y4=sin(x4) = sin(0,509015)=0,487317
y5=sin(x5) = sin(0,868511)=0,763367
I = /10(0,131118+0,471494+0,706825+0,487317+0,763367) =
=/10*2,560121=0,8038779.
3. Описание программы Integral. pas. Алгоритм.
Процедура VVOD - заполняет массив, содержащий в себе аргументы xi
Процедура FORM - используя массив, содержащий аргументы xi заполняет массив yi
Процедура CHEB - используя массивы xi и yi, высчитывает по квадратурной формуле Чебышева приближенное значение интеграла.
Процедура TABL - это подпрограмма, осуществляющая вывод таблицы узлов (аргумент - функция)
При запуске программы нужно ввести границы интегрирования.
После ввода границ интегрирования используется процедура VVOD, а затем высчитывается и выводиться на экран шаг табулирования функции h.
После этого используем процедуры FORM и CHEB .
Получив результат, выводим таблицу ( процедура TABL ) и интеграл.
4. Заключение и выводы.
Таким образом очевидно, что при вычислении определенных интегралов с помощью квадратурных формул, а в частности по формуле Чебышева не дает нам точного значения, а только приближенное.
Чтобы максимально приблизиться к достоверному значению интеграла нужно уметь правильно выбрать метод и формулу, по которой будет вестись расчет. Так же очень важно то, какой будет взят шаг интегрирования.
Хотя численные методы и не дают очень точного значения интеграла, но они очень важны, так как не всегда можно решить задачу интегрирования аналитическим способом.
5.Список литературы:
1. Ракитин Т.А., Первушин В.А. “Практическое руководство по численным методам с приложением программ на языке Basic“
2. Крылов В.И. “Приближенные вычисления интегралов“ - М. : Физмат.
3. Демидович и Марон “Основы вычислительной математики“
4. Копченова и Марон “Вычислительная математика в примерах и задачах”
5. Вольвачев А.Н., Крисевич В.С. Программирование на языке Паскаль для ПЭВМ ЕС. Минск.: 1989 г.
6. Зуев Е.А. Язык программирования Turbo Pascal. М.1992 г.
7. Скляров В.А. Знакомьтесь: Паскаль. М. 1988 г.
6. Листинг программы.
Программа написана на языке Tubro Pascal 7.0 для MS-DOS. Ниже приведен ее листинг:
program integral;
uses crt;
const n=5;
k=-0.832498;
l=-0.374541;
z=0.0;
type aa=array[1..n] of real;
var x,y:aa;
a,b,h,ich:real;
{ заполнение х-сов в массив х[5] }
procedure vvod(var a,b:real;var c:aa);
var i:integer;
t:aa;
Begin
t[1]:=k;
t[2]:=l;
t[3]:=z;
t[4]:=l;
t[5]:=k;
for i:=1 to n-1 do
c[i]:=((b+a)/2+(b-a)/2*t[i]);
for i:=n-1 to n do
c[i]:=1 - c[n+1-i];
end;
{ заполнение y-ков в массиве у[5] }
procedure form(var x:aa; var y:aa);
var i:integer;
Begin
for i:=1 to n do
y[i]:=sin(x[i]); {функция}
end;
{ процедура для расчета интеграла по квадратурной
формуле Чебышева }
procedure cheb(var y:aa;var ich:real);
var i:integer;
Begin
ich:=0;
for i:=1 to n do
ich:=ich+y[i]*h;
end;
{ процедура вывода таблицы}
procedure tabl;
var i:integer;
Begin
writeln( ___________________________________ );
writeln(| i | t | x | y |);
writeln( ___________________________________ );
writeln(| 1 |,k:9:6,