VPF::Помогите с функцией Лапласа!!! - Форум программистов

Kamdi

Дата 12.12.2004, 15:11 (ссылка)

(голосов: 0)

Загрузка ...

Unregistered

Уважаемые друзья! Помогите мне с поиском библиотеки в С++, которая позволяет вычислить значение функции Лапласа для нормального распределения. Как она называется?! smile

Обращаю внимание на то, что необходима именно функция, а не уравнение Лапласа...
С уважением, Kamdi.

_hunter

Дата 14.12.2004, 00:28 (ссылка)

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

Загрузка ...

Эксперт

Профиль
Группа: Участник Клуба
Сообщений: 8564
Регистрация: 24.6.2003
Где: Europe::Ukraine:: Kiev

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

вряд ли такая есть...
IMHO проще самому посчитать ( как -- в алгоритмах должны помочь )

--------------------

Tempora mutantur, et nos mutamur in illis...

Hroft

Дата 15.12.2004, 10:50 (ссылка)

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

Загрузка ...

Опытный

Профиль
Группа: Участник
Сообщений: 310
Регистрация: 20.10.2003
Где: Москва

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

Так функция Лапласа - эт же Ф(x) = integr(0,x,"exp(-t*t/2)"), или нет? Если так, то ее любым методом вычисления интеграла надо считать...
PS: какая связь у функции Лапласа с уравнением Лапласа?

Hroft

Дата 17.12.2004, 10:11 (ссылка)

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

Загрузка ...

Опытный

Профиль
Группа: Участник
Сообщений: 310
Регистрация: 20.10.2003
Где: Москва

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

Код

#include <cmath>
#include <iostream>

const double eps = 0.1e-5;

double func(double x)
{
	return std::exp(-x*x/2);
}

double simpsonIntegr(double a,double b)
{
	return ((b-a) / 6) * (func(a) + 4*func((a+b)/2) + func(b));
}

double LaplasFunction(double x,double eps = ::eps)
{
	double a = 0.0;
	double b = x;
	double tst, tst1, incr;
	unsigned int c = 2;
	tst1 = simpsonIntegr(a,b);
	do {
  tst = tst1;
  tst1 = 0;
  incr = (b - a) / c;
  for (unsigned int i=0;i<c;++i) {
  	tst1 = tst1 + simpsonIntegr(a + incr * i, a + incr * (i + 1));
  }
  c += 1;
	} while (fabs(tst - tst1) >= eps);
	return tst1;
}

int main()
{
	for (double x=-4;x<4;x+=0.5) {
  std::cout<<"Fi("<<x<<") == "<<LaplasFunction(x)<<std::endl;
	}
	return 0;
}

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Общие вопросы \| Следующая тема »