VPF::Какое условие добавить в метод трапеций, для обх..

ЧерныйМальчик

Дата 26.12.2008, 23:52 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 26.12.2008

Репутация: нет
Всего: нет

Ищу интеграл функции (a1*sin(a2*x)*sin(a3*x)+a4*sin(a5*x)*cos(a6*x))/ (b1*sin(b2*x)*cos(b3*x)) методом трапеций. В некоторых х у функции резко увеличивается значение и кол-во таких экстремумов увеличивается с увелечением количества точек интегрирования. А эти экстремумы сильно влияют на итоговое значение интеграла. Т.е. при 100 и 200 точках интегрирования получаются разные значения интеграла.

Как можно обходить эти "скачки" ?

Заранее спасибо кто что хорошего подскажет.

Цитата

h:=(x2-x1)/(n-1);
xv[1]:=x1; //нижняя граница интегрирования
null:=false;
prom:=0;//результат интеграла накапливает
h1:=h;

for i:=1 to n do
begin
fx[i]:=F(xv[i], a1, a2, a3, a4, a5, a6, b1, b2, b3);
if i<>1 then
prom:=prom+(fx[i-1]+ fx[i])*h1/2;

xv[i+1]:= xv[i]+h;
end;

end; //for i

hkdkest

Дата 27.12.2008, 23:45 (ссылка)

(нет голосов)

Загрузка ...

Опытный

Профиль
Группа: Участник
Сообщений: 300
Регистрация: 30.11.2008

Репутация: нет
Всего: 1

вот функция вроде норм работает : Вычисление определенного интеграла методом трапеций с заданной точностью

THandle

Дата 28.12.2008, 00:05 (ссылка)

(нет голосов)

Загрузка ...

Хранитель Клуба
Group Icon

Награды: 1

Профиль
Группа: Админ
Сообщений: 3639
Регистрация: 31.7.2007
Где: Moscow, Dubai

Репутация: 50
Всего: 372

hkdkest, вообще этот материал взят из drkb, который я бы посоветовал скачать ТС. Там много чего полезного есть.

Это сообщение отредактировал(а) THandle - 28.12.2008, 00:32

ЧерныйМальчик

Дата 28.12.2008, 11:02 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 2
Регистрация: 26.12.2008

Репутация: нет
Всего: нет

Спасибо.

Это сообщение отредактировал(а) ЧерныйМальчик - 28.12.2008, 22:10

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Delphi: Для новичков \| Следующая тема »