Пример: Эффекты TOL и метода определенного интегрирования
Параметр допуска
Посмотрите, как допуск сходимости системных переменных (TOL) отражается на результатах вычисления определенных интегралов. TOL можно задать на вкладке Расчет (Calculation) в группе Параметры документа (Worksheet Settings) или в самом документе.
1. Вычислите следующий интеграл.
Для вычисления ответа используется значение по умолчанию TOL:
2. Сделайте допуск свободнее и вычислите интеграл заново.
3. Сделайте допуск строже и вычислите интеграл заново.
Прерывные функции
Прерывные функции могут оказаться нестабильными при некоторых значениях интегралов, если у них большие амплитуды и резкие прерывности. Перед интегрированием необходимо определить интервал интегрирования, который содержит основную часть площади. Можно также поэкспериментировать с TOL.
1. Используйте функцию Heaviside StepΦ для определения прерывистой пилообразной функции.
2. Определите переменную допуска сходимости.
Минимальное значение TOL (10-15) слишком мало для крайне прерывистых функций, и алгоритм может вернуть неточные оценки.
3. Постройте график пилообразной функции f(x) и соответствующей интегральной функции Int(x). Масштабируйте интегральную функцию по коэффициенту 4.
Интегральная функция имеет пик в окрестностях x = 15. Это может произойти при интегрировании прерывных функций. Уменьшение TOL ниже 10-10 только ухудшит результат.
4. Чтобы получить точный ответ, разбейте интеграл на части, соответствующие прерывности функции.
5. Постройте график функции f(x) и вновь определенной функции Int2(x).
Пик в окрестностях x = 15 исчезает.
Ограничения интеграции
Одним из ограничений численного интегрирования является то, что очень узкие импульсы в функции с обычно нулевым значением часто интегрируются до нуля. Обычно если подынтегральное выражение равно нулю на более чем 95 % области интегрирования, то алгоритм может не вычислить подынтегральное выражение в одной из ненулевых точек.
1. Определите и постройте график узкой пульсации шириной 0.05 в сигнале с нулевым значением шириной 1.0.
2. Вычислите числовой интеграл функции.
Результат должен равняться площади пульсации 0.05x1.0 или 0.05.
3. Устраните эту проблему путем интеграции по небольшой области, содержащей ненулевую часть подынтегрального выражения.