Точность численного интегрирования в Matlab

Я пытаюсь интегрировать аналитическую функцию (смесь функций sqrt и trig) в области прямоугольника. У этого нет никакого особенного пункта в области и, кажется, прекрасный кандидат, чтобы использовать dblquad, Мой вопрос состоит в том, как оценить точность числового значения, которое Matlab предоставил мне. Не зная точного значения интеграции, как мы можем обосновать значащие цифры? Когда вам нужно дать значение с определенными цифрами точности, вы должны быть в состоянии оправдать. Возможно ли достичь этого, учитывая, что значение рассчитывается с помощью Matlab?

2 ответа

Если вы не установите это иначе, dblquad использует порог допуска по умолчанию (10-6 в последних выпусках) для абсолютной квадратурной ошибки. Аппроксимация интеграла будет в пределах погрешности, не превышающей указанный допуск.

Вы могли бы взглянуть на исходный код для dblquadгде-то он будет использовать определенное количество "шагов". Я полагаю, вы могли бы создать новый m-файл с важными битами, которые заставят работать интеграл, и поиграть с количеством шагов, пока компьютер не займет много времени и не изменит результат. Лично я использую пользовательское правило Симпсонов для числовой интеграции и просто меняю N (количество шагов) на какое-то большое число.

Другие вопросы по тегам