Добрый день.
Передо мной встала задача вычисления значений функции определяемой рекурсивно через интегрирования самой себя порядком ниже.
В упрощеном виде всё выглядит примерно так:
Функции
и
- обычные, не рекурсивные. Индекс
у
призван показать что для разных
она вычисляется по-разному.
Интегралы не берутся, значения нахожу численно, квадратурными формулами.
У меня есть с чем сравнивать полученные результаты (некоторые таблицы значений данной фунции представлены в литературе, также я посчитал её методом Монте-Карло) и результаты выглядят корректно.
Однако у меня вы получается оценить погрешность вычисления.
Предположим, для определенности, что мы интегрируем методом средних прямоугольников. Оценить погрешность при вычислении лишь самого верхнего в рекурсии интеграла будет неверным - ведь значение функции в каждом из узлов на которые разбиваем интервал интегрирования так же получено с погрешностью.
Оценим погрешность уровня
как
где
а
, где
- количество узлов.
Но значение функции в каждом из
узлов получено с погрешностью
, что дает суммарную погрешность от вычисления функции в узлах
И т.д.
Каким же образом верно оценить погрешность?
Перемножать все эти погрешности, видимо, будет неверным, сделав это получим огромное число в
Насколько обосновано в таком случае пользоваться методом двойного счета?
Подсчитав интегралы с
и
узлами на
каждом уровне имеем ли мы право считать относительную погрешность как
?
Посоветуйте что можно предпринять в данной задаче.
Спасибо.