VPF::[С#, Алгоритм] Методы численного интегрирования

Zhotick

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

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

Загрузка ...

Новичок

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

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

Здравствуйте, помогите пожалуйста.
Я пытаюсь реализовать метод численного интегрирования методом "Выходящих прямоугольников", вот теоретическая база: http://ru.wikipedia.org/wiki/%D0%A2%D1%80%....86.D0.B8.D0.B9
Где ошибка? Может кто поможет?
Код ниже:

Код


using System;
class Интеграл1 
{
struct Integral
{
double a, b;
int m;

public Integral(double ina, double inb, int inm)
{
a = ina < inb ? ina : inb;
b = ina < inb ? inb : ina;
m = inm;
}
public override string ToString()
{
string s;
s = String.Format("Точное значение= {0:f3}\nПриближённое значение={1:f3}",
ИнтТрапеции, ИнтЛейбниц);
double Delta = (ИнтТрапеции - ИнтЛейбниц) / ИнтЛейбниц * 100;
Delta = Math.Abs(Delta);
s += String.Format("\nПогрешность = {0:f3} %",Delta);
return s;
}
public double fx(double x)
{
return Math.Sin(x) - x;
}
public double Fx(double x)
{
return -Math.Cos(x) - 0.5*x*x;
}
public double p
{ get { return Math.Pow((b - a), 3)*Math.Sqrt(3) / (48 * m * m); } }


public double h
{
get
{
return (b-a)/m;
}
}
public double ИнтЛейбниц
{
get
{
return Fx(b) - Fx(a);
}
}
public double ИнтТрапеции
{
get
{
double sum = (fx(a)+fx(b));
for (double i = a + h; i < b; i += h) sum += fx(i);
sum *= h;
sum += p;
return sum;
}
}
}
static void Main()
{
Integral obj=new Integral( 3.1415926535897932384626433832795,4*3.1415926535897932384626433832795/3,800);
Console.WriteLine(obj.ToString() );

}
}

Спасибо.

walter66

Дата 10.6.2009, 08:00 (ссылка)

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

Загрузка ...

Новичок

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

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

Милочка, ну кто же так цикл организует? :
for (double i = a + h; i < b; i += h) sum += fx(i);

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

int i,n=(b-a)/h;
for(i=0;i<=n;i++)
sum+=fx(a+i*h);

Остальное не смотрел.

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