Полная версия этой страницы: численное интегрирование
Форум разработчиков электроники ELECTRONIX.ru > Сайт и форум > В помощь начинающему > Программирование
ищу исходники взятия интеграла методом прямоугольников и сплайнов. нужны исходники на С. Наверняка такое в институтах на лабах есть.
Рекомендую книгу "numerical recipes in с".
http://astronu.jinr.ru/wiki/upload/d/d6/Nu...lRecipesinC.pdf
Можно найти и исходники, если из текста брать неудобно.
http://astronu.jinr.ru/wiki/upload/d/d6/Nu...lRecipesinC.pdf
Можно найти и исходники, если из текста брать неудобно.
Цитата(TigerSHARC @ Aug 5 2012, 13:00) 
ищу исходники взятия интеграла методом прямоугольников и сплайнов. нужны исходники на С.
Код
double Integral(double * x, double * y, int num){
double r = 0;
for (int i = 1; i < num; i++) r += y[i] * (x[i] - x[i-1]);
// for (int i = 1; i < num-1; i++){r += 0.5L * y[i] * (x[i+1]-x[i-1]);} r += 0.5L * (y[0] * (x[1]-x[0]) + y[num-1] * (x [num-1]-x[num-2]));
// for (int i = 1; i < num; i++) r += 0.5L * (y[i] + y[i-1]) * (x[i] - x[i-1]);
return r;
}
double r = 0;
for (int i = 1; i < num; i++) r += y[i] * (x[i] - x[i-1]);
// for (int i = 1; i < num-1; i++){r += 0.5L * y[i] * (x[i+1]-x[i-1]);} r += 0.5L * (y[0] * (x[1]-x[0]) + y[num-1] * (x [num-1]-x[num-2]));
// for (int i = 1; i < num; i++) r += 0.5L * (y[i] + y[i-1]) * (x[i] - x[i-1]);
return r;
}
есть некие различия для неравномерной сетки по х.
еще есть промежуточный вариант симпсона, где не через две точки проводится прямая, как в коде выше, а через три - парабола, ну или через 4 кубический полином. Радость по сравнению со сплайнами, что не надо по массиву несколько раз бегать, для построения честных кубических сплайнов, можно считать на ходу.
за исходниками сюда: alglib
Для просмотра полной версии этой страницы, пожалуйста, пройдите по ссылке.