Калькулятор длины дуги по правилу средних точек

Калькулятор длины дуги по правилу средних точек приближает длину кривой: он делит интеграл длины дуги на прямоугольники и берёт значение в центре каждого из них. Введите функцию, задайте интервал, выберите число подынтервалов, и калькулятор покажет приближение, полную таблицу средних точек, оценку погрешности и живую диаграмму прямоугольников. Выполните эти шаги, чтобы получить правильный ответ.

1. Выберите режим расчёта. Используйте «Явная функция» для кривой, записанной как y = f(x), например y = x² или y = sin(x). Используйте «Параметрическая кривая» для кривой, записанной как x = x(t) и y = y(t), например окружности или эллипса, заданного параметром.
2. Введите функцию. Введите функцию в поле ввода, используя ^ для степеней, функцию квадратного корня для корней, а sin(), cos(), exp() и ln() для стандартных функций. Калькулятор находит производную автоматически, или вы можете ввести производную сами, если уже знаете её.
3. Задайте интервал. Введите нижний предел a и верхний предел b. Интервал [a, b] это участок кривой, длину которого нужно найти. Если ввести b меньше a, калькулятор всё равно вернёт положительную длину, потому что длина дуги никогда не бывает отрицательной.
4. Выберите число подынтервалов n. Используйте ползунок или введите значение. Предустановки включают 1, 2, 4, 5, 8, 10, 20, 50, 100, 500, 1000 и 5000. Чем больше n, тем уже прямоугольники и тем ближе приближение. Начните с n = 4, чтобы ясно видеть ход решения, затем увеличивайте n и следите, как ответ устанавливается.
5. Выберите точность отображения. Выберите 2, 4, 6, 8 или 10 знаков после запятой. Калькулятор хранит полную точность внутри и округляет только отображаемый результат, поэтому смена точности никогда не меняет сам расчёт.
6. Нажмите «Вычислить длину дуги». Главный результат это приближение L_n по правилу средних точек, чётко помеченное как приближение. Ниже вы получите ширину подынтервала, первую и последнюю средние точки, первое и последнее значения подынтегральной функции, сумму значений подынтегральной функции, длину хорды между конечными точками, высокоточное эталонное значение длины, а также оценки абсолютной и относительной погрешности.
7. Изучите диаграммы. На первой диаграмме показана исходная кривая с выделенной дугой и отмеченными средними точками. На второй диаграмме показана подынтегральная функция длины дуги с прямоугольниками средних точек под ней. Общая площадь этих прямоугольников и есть приближение, поэтому две диаграммы нельзя путать: прямоугольники принадлежат подынтегральной функции, а не исходной кривой.
8. Прочитайте таблицу расчёта. Каждый подынтервал получает строку со своими границами, средней точкой, производной в средней точке, значением подынтегральной функции и вкладом прямоугольника. В последней строке показана сумма, равная главному результату.
9. Сравните сходимость и попробуйте экстраполяцию Ричардсона. Панель сходимости показывает приближения для n, 2n и 4n рядом, вместе с разностями между ними. Опция Ричардсона объединяет два приближения в уточнённую оценку. Обе возможности объяснены в разделах руководства ниже.
10. Скопируйте, экспортируйте или сбросьте. Используйте «Копировать результат» для ответа, экспортируйте таблицу в CSV для своей работы или нажмите «Сбросить», чтобы очистить панель. Кнопки предустановок мгновенно загружают проверенные примеры: параболу, прямую, горизонтальную прямую, четверть окружности, кубическую кривую и синусоиду.

Попробуйте пример. y = x² на [0, 2] при n = 4 даёт L₄ = 4.6267234873 единиц при эталонном значении 4.6467837624 единиц. y = 3x + 1 на [0, 4] при n = 4 даёт ровно 12.6491106407 единиц. y = 2 на [1, 5] при n = 4 даёт ровно 4.0000000000 единиц. Параметрическая четверть окружности радиуса 5 даёт 7.8539816340 единиц. y = x³ на [0, 1] при n = 4 даёт 1.5330289642 единиц, а y = sin(x) на [0, π] при n = 8 даёт 3.8201978627 единиц.

Важные советы. Всегда сначала дифференцируйте функцию: в подынтегральной функции используется f'(x), а не f(x). Берите значения в средних точках, а не в концах, иначе вы вычислите другое правило. Для тригонометрических функций держите углы в радианах. Не округляйте средние точки и значения подынтегральной функции в середине расчёта, потому что эти малые округления суммируются по многим прямоугольникам. Если функция не определена в средней точке, калькулятор сообщит об ошибке, а не пропустит значение молча.

Численное интегрирование это искусство оценивать определённый интеграл, когда точную первообразную найти трудно или когда известны только отдельные значения функции. Правило средних точек один из самых простых и надёжных методов численного интегрирования. Оно оценивает площадь под кривой, накрывая её прямоугольниками и складывая их площади.

Сумма Римана делит интервал [a, b] на n равных подынтервалов, каждый шириной Δx = (b - a) / n, и строит на каждом подынтервале один прямоугольник. Высота каждого прямоугольника равна значению функции в одной точке выборки внутри этого подынтервала. Суммы по левым концам берут значение на левом краю, суммы по правым концам на правом краю, а правило средних точек берёт значение в центре каждого подынтервала.

Выборка в центре это ключевая идея. На коротком подынтервале кривая, которая растёт с одной стороны от средней точки, примерно так же падает с другой стороны, поэтому завышение на одной половине примерно компенсирует занижение на другой. Именно эта компенсация делает правило средних точек точнее правил левых и правых концов при том же числе прямоугольников, а его погрешность убывает как квадрат ширины подынтервала.

Для функции g(x) на [a, b] составное правило средних точек записывается так:

Здесь x_i* это средняя точка i-го подынтервала, а сумма берётся по всем n подынтервалам. Каждое слагаемое Δx × g(x_i*) это площадь одного прямоугольника, а сумма площадей всех прямоугольников это приближение интеграла.

Правило средних точек применяется далеко за пределами домашних заданий по математическому анализу. Инженеры используют его для оценки интегралов в расчётах прочности и обработке сигналов, физики в численных моделях, а компьютерная графика для измерения длин кривых и площадей поверхностей. Сочетание простоты и точности второго порядка делает его стандартным первым выбором всегда, когда у интеграла нет удобного замкнутого выражения.

Почему средние точки лучше концов. Возьмите g(x) = x на [0, 1] при n = 1. Сумма по левым концам даёт 0, сумма по правым концам даёт 1, а сумма по средним точкам даёт 0.5, что и есть точный ответ. Средняя точка отражает среднее поведение функции на подынтервале, а конец видит только одну его сторону. Это преимущество сохраняется для любой гладкой функции, и именно поэтому правило средних точек это сумма Римана по умолчанию в большинстве численных работ.

Длина дуги это расстояние, измеренное вдоль кривой, а не расстояние по прямой между её концами. Дорога, которая вьётся через холмы, длиннее прямой, проведённой между её концами на карте. Изогнутая проволока длиннее промежутка, который она перекрывает. Длина дуги измеряет сам извилистый путь.

Чтобы найти эту длину, представьте, что кривая разрезана на множество крошечных кусочков. Каждый крошечный кусочек так короток, что он почти прямой, поэтому его длина близка к прямой, соединяющей его концы. По теореме Пифагора крошечный кусочек, который смещается на Δx по горизонтали и на Δy по вертикали, имеет длину, близкую к √(Δx² + Δy²). Сложение всех крошечных кусочков даёт общую длину, а когда кусочки стремятся к нулю, сумма превращается в интеграл.

Для кривой y = f(x) крошечный горизонтальный шаг dx даёт вертикальное изменение dy = f'(x) dx, поэтому крошечный элемент длины равен:

Сложение всех элементов длины от a до b даёт интеграл длины дуги:

Эта формула объясняет, почему здесь нужна производная. Крутая кривая быстро меняет высоту, поэтому f'(x) велика, подынтегральная функция √(1 + [f'(x)]²) велика, и каждая единица горизонтального расстояния даёт больше одной единицы длины кривой. У плоской кривой f'(x) близка к нулю, подынтегральная функция близка к 1, а длина дуги близка к горизонтальному расстоянию b - a.

Длина дуги всегда не меньше хорды. Прямая между концами это кратчайший путь между ними, поэтому кривая никогда не может быть короче её. Для параболы y = x² на [0, 2] хорда от (0, 0) до (2, 4) равна √(20) = 4.4721359550 единиц, а истинная длина дуги равна 4.6467837624 единиц, то есть длиннее на 0.1746478074 единиц. Для синусоиды y = sin(x) на [0, π] хорда равна π = 3.1415926536 единиц, а дуга равна 3.8201977890 единиц, то есть длиннее на 0.6786051354 единиц. Проверка, что дуга длиннее хорды, это быстрая проверка разумности любого ответа о длине дуги.

У большинства интегралов длины дуги нет элементарной первообразной. Интеграл параболы ∫ √(1 + 4x²) dx можно взять подстановкой, а интеграл синуса ∫ √(1 + cos²x) dx вообще нельзя выразить элементарными функциями. Именно здесь правило средних точек занимает своё место: оно оценивает интеграл численно, без первообразной.

Расчёт длины дуги по правилу средних точек начинается с точного интеграла длины дуги, а интеграл заменяется суммой площадей прямоугольников. Определим подынтегральную функцию длины дуги:

Тогда точная длина дуги равна L = ∫_a^b g(x) dx, и правило средних точек приближает этот интеграл. Разделите [a, b] на n равных подынтервалов. Ширина подынтервала равна:

Средняя точка i-го подынтервала, при счёте i от 0, равна:

Вычислите подынтегральную функцию в каждой средней точке, умножьте каждое значение на Δx и сложите результаты:

Каждый символ имеет точное значение:

Конкретная проверка. Для y = x² на [0, 2] при n = 4 ширина равна Δx = (2 - 0) / 4 = 0.5, средние точки равны 0.25, 0.75, 1.25 и 1.75, а значения подынтегральной функции равны 1.1180339887, 1.8027756377, 2.6925824036 и 3.6400549446. Их сумма равна 9.2534469747, а умножение на 0.5 даёт L₄ = 4.6267234873 единиц. Каждое число в этой цепочке показано на панели пошагового решения калькулятора и в его таблице расчёта.

Для параметрической кривой x = x(t), y = y(t) та же идея применяется к интервалу параметра [t₁, t₂], где скорость v(t) = √([x'(t)]² + [y'(t)]²) играет роль подынтегральной функции. Параметрическая формула рассмотрена в отдельном разделе ниже.

Расчёт по правилу средних точек состоит из десяти шагов, и калькулятор выполняет их все автоматически. Если один раз пройти их вручную, пользоваться инструментом станет гораздо проще. В разобранном примере ниже используется y = x² на [0, 2] при n = 4, тот же пример, который калькулятор загружает как первую предустановку.

1. Определите функцию. f(x) = x².
2. Найдите её первую производную. f'(x) = 2x.
3. Определите интервал интегрирования. a = 0, b = 2.
4. Выберите число подынтервалов. n = 4.
5. Вычислите ширину подынтервала. Δx = (2 - 0) / 4 = 0.5.
6. Определите среднюю точку каждого подынтервала. x₀* = 0.25, x₁* = 0.75, x₂* = 1.25, x₃* = 1.75. Каждая средняя точка находится ровно посередине между границами своего подынтервала: 0.25 это середина между 0 и 0.5, 0.75 середина между 0.5 и 1, и так далее.
7. Вычислите подынтегральную функцию длины дуги в каждой средней точке. g(x) = √(1 + 4x²), поэтому g(0.25) = 1.1180339887, g(0.75) = 1.8027756377, g(1.25) = 2.6925824036 и g(1.75) = 3.6400549446.
8. Умножьте каждое значение на ширину подынтервала. Четыре вклада прямоугольников равны 0.5 × 1.1180339887 = 0.5590169944, 0.5 × 1.8027756377 = 0.9013878189, 0.5 × 2.6925824036 = 1.3462912018 и 0.5 × 3.6400549446 = 1.8200274723.
9. Сложите все вклады. 0.5590169944 + 0.9013878189 + 1.3462912018 + 1.8200274723 = 4.6267234873.
10. Истолкуйте приближённую длину дуги. L₄ = 4.6267234873 единиц. Высокоточное эталонное значение равно 4.6467837624 единиц, поэтому это приближение короче на 0.0200602751 единиц, а относительная погрешность составляет 0.4317%.

Те же десять шагов работают для любой функции. Для y = x³ на [0, 1] при n = 4 производная равна f'(x) = 3x², ширина равна Δx = 0.25, средние точки равны 0.125, 0.375, 0.625 и 0.875, значения подынтегральной функции равны 1.0010980300, 1.0853471867, 1.5405489332 и 2.5051217067, их сумма равна 6.1321158566, а умножение на 0.25 даёт L₄ = 1.5330289642 единиц при эталонном значении 1.5478656547 единиц.

Второй разобранный пример с тригонометрической функцией. Для y = sin(x) на [0, π] при n = 8 производная равна f'(x) = cos(x), ширина равна Δx = π / 8 = 0.3926990817, а восемь средних точек это π/16, 3π/16, 5π/16, 7π/16, 9π/16, 11π/16, 13π/16 и 15π/16. Значения подынтегральной функции по порядку равны 1.4006926023, 1.3005159423, 1.1439660326, 1.0188524102, 1.0188524102, 1.1439660326, 1.3005159423 и 1.4006926023. Их сумма равна 9.7280539749, а умножение на 0.3926990817 даёт L₈ = 3.8201978627 единиц при эталонном значении 3.8201977890 единиц, то есть абсолютная погрешность составляет всего 0.0000000736 единиц.

Каждый прямоугольник средних точек это небольшая наглядная картина одного кусочка интеграла. У него есть три части, и у каждой части есть ясный смысл.

Для параболы y = x² на [0, 2] при n = 4 все четыре прямоугольника имеют ширину 0.5. Их высоты равны 1.1180339887, 1.8027756377, 2.6925824036 и 3.6400549446, а их площади равны 0.5590169944, 0.9013878189, 1.3462912018 и 1.8200274723. Общая площадь равна 4.6267234873 единиц, что в точности равно приближению L₄ по правилу средних точек.

Почему сумма площадей приближает интеграл. Интеграл ∫_a^b g(x) dx это точная площадь под кривой подынтегральной функции. Каждый прямоугольник заменяет изогнутую полоску над одним подынтервалом плоской полоской той же ширины, используя высоту в средней точке как представительную высоту. Там, где подынтегральная функция изгибается вверх, плоская вершина срезает немного площади с одной стороны и добавляет немного с другой. Поскольку значение берётся в центре, эти две погрешности почти компенсируют друг друга, и площадь прямоугольника близка к истинной площади полоски. Сложение прямоугольников даёт оценку всей площади, а значит, и длины дуги.

Самая частая ошибка с диаграммой. Прямоугольники нужно рисовать под кривой подынтегральной функции g(x) = √(1 + [f'(x)]²), а не под исходной кривой f(x). Исходная кривая и подынтегральная функция это разные графики с разными высотами: для y = x² на [0, 2] исходная кривая идёт от 0 до 4, а подынтегральная функция от 1 до √(17) = 4.1231056256. Прямоугольники, нарисованные под исходной кривой, в сумме дали бы площадь под f(x), которая не имеет отношения к длине дуги. Калькулятор держит два графика отдельно, а диаграмма выше показывает прямоугольники только под подынтегральной функцией.

Как читать прямоугольник. Наведение курсора на прямоугольник в калькуляторе или его выбор показывает номер подынтервала, его левую и правую границы, среднюю точку, производную в средней точке, значение подынтегральной функции и вклад прямоугольника. Именно эта детализация делает таблицу полезной для домашних заданий: каждое число в итоговом ответе можно проследить до одного прямоугольника.

Число подынтервалов n это самый важный элемент управления калькулятором, потому что именно оно решает, насколько мелкими будут прямоугольники. Малое n означает несколько широких прямоугольников. Большое n означает много узких прямоугольников, которые близко следуют за кривой подынтегральной функции.

Что происходит при малом n. При n = 1 весь интервал это один прямоугольник, а ответ равен Δx × g(средняя точка). Для параболы y = x² на [0, 2] при n = 1 получается 2 × √(5) = 4.4721359550 единиц, что короче истинных 4.6467837624 единиц на 0.1746478074 единиц. Один прямоугольник не может следовать за подынтегральной функцией, пока она растёт от 1 до 4.1231056256, поэтому погрешность велика. При n = 2 ответ улучшается до 4.5764912225 единиц. При n = 4 он равен 4.6267234873 единиц. При n = 8 он равен 4.6417306581 единиц. При n = 16 он равен 4.6455205322 единиц, а при n = 32 он равен 4.6464679593 единиц.

Как убывает погрешность. Для гладкой подынтегральной функции погрешность правила средних точек пропорциональна 1 / n². Удвоение n сокращает погрешность примерно в четыре раза. Числа для параболы ясно это показывают: погрешность падает с 0.0200602751 единиц при n = 4 до 0.0050531044 единиц при n = 8, затем до 0.0012632302 единиц при n = 16, затем до 0.0003158031 единиц при n = 32. Каждое удвоение делит погрешность примерно на четыре. Такое поведение 1 / n² называется сходимостью второго порядка, и именно поэтому правило средних точек так эффективно: небольшое увеличение n даёт большое улучшение точности.

Почему улучшение не всегда монотонно. Погрешность обычно убывает с ростом n, но она не обязана убывать на каждом шаге для каждой функции. Синусоида показывает мягкий вариант этого: при n = 2 получается 3.8476494905 единиц, что превышает эталонное значение 3.8201977890 единиц, а при n = 4 получается 3.8204517876 единиц, что гораздо ближе. Приближение перешло от значения выше истинного к значению чуть выше истинного, а погрешность в целом всё равно уменьшилась. Для функций с резкими изгибами, колебаниями или почти особенностями погрешность может колебаться при изменении n, хотя долгосрочная тенденция направлена вниз. Никогда не считайте, что большее n всегда даёт лучший ответ на одном шаге: вместо этого проверяйте таблицу сходимости.

Цена большего n. Больше подынтервалов означает больше вычислений: больше значений производной, больше значений подынтегральной функции и больше строк таблицы. При n = 5000 калькулятор выполняет 5000 вычислений, что быстро для компьютера, но даёт таблицу, слишком длинную для построчного чтения, поэтому калькулятор разбивает её на страницы. Практическое правило: увеличивайте n, пока ответ не перестанет меняться при нужной вам точности. Если при n = 100 и n = 200 ответы совпадают до шести знаков после запятой, то n = 100 достаточно, а n = 5000 добавляет только работу.

Выбор n на практике. Используйте n = 4 или n = 8 во время обучения, потому что таблица остаётся достаточно короткой для проверки вручную. Используйте n = 50 или n = 100 для домашних заданий, где нужны четыре или пять верных знаков после запятой. Используйте n = 1000 и более только тогда, когда нужна очень высокая точность или когда подынтегральная функция быстро меняется на интервале. Панель сходимости показывает n, 2n и 4n вместе, поэтому вы видите, как ответ устанавливается, без догадок.

Правило средних точек это один из трёх классических методов численного интегрирования, и каждый из них берёт значения подынтегральной функции по-разному. Знание их различий помогает выбрать правильный инструмент и понять, почему их ответы не совпадают.

Правило средних точек. Берёт значение в центре каждого подынтервала и строит прямоугольники. Его погрешность имеет порядок O(1 / n²), тот же порядок, что и у правила трапеций, но правило средних точек обычно примерно вдвое точнее правила трапеций при том же n, потому что выборка в центре сокращает главную погрешность. Ему нужно всего n вычислений, по одному на подынтервал, и оно работает при любом положительном n.

Правило трапеций. Берёт значения на концах и соединяет их прямыми, образуя трапеции. У него тот же порядок погрешности O(1 / n²), что и у правила средних точек, но его погрешность обычно имеет противоположный знак: для выпуклой подынтегральной функции правило трапеций завышает, а правило средних точек занижает, или наоборот. Усреднение этих двух даёт комбинацию трапеций и средних точек, которая лежит в основе правила Симпсона. Попробуйте калькулятор правила трапеций, чтобы сравнить два метода на одной функции.

Правило Симпсона. Строит параболу через два конца и среднюю точку каждой пары подынтервалов, для чего нужно чётное число подынтервалов. Его погрешность имеет порядок O(1 / n⁴), поэтому для гладких функций оно сходится гораздо быстрее двух других. Цена этого ограничение на чётность n и большее число вычислений на подынтервал. Попробуйте калькулятор правила Симпсона, когда нужна высокая точность при малом числе подынтервалов.

Сравнение чисел. Для параболы y = x² на [0, 2] при n = 4 правило средних точек даёт 4.6267234873 единиц, правило трапеций даёт 4.6870560064 единиц, а правило Симпсона даёт 4.6502010785 единиц при эталонном значении 4.6467837624 единиц. Погрешность средних точек равна 0.0200602751 единиц, погрешность трапеций 0.0402722440 единиц, а погрешность Симпсона 0.0034173160 единиц. Здесь выигрывает Симпсон, как обычно бывает для гладких функций, но правило средних точек проще, не требует чётного n и является естественным выбором, когда доступны только значения в средних точках.

Какое правило использовать? Используйте правило средних точек, когда нужна простая надёжная оценка, когда вы берёте значения функции в равномерно расположенных внутренних точках или когда вы изучаете суммы Римана. Используйте правило трапеций, когда значения на концах это естественные данные, например измерения, взятые на концах каждого интервала. Используйте правило Симпсона, когда функция гладкая и нужна высокая точность при малом n. Для длины дуги правило средних точек имеет практическое преимущество: оно никогда не вычисляет подынтегральную функцию на концах, где функция может быть неопределена или иметь особенность, поэтому оно справляется с кривыми, с которых два других правила не могут начать.

Каждое приближение по правилу средних точек отличается от истинной длины дуги на некоторую величину, и понимание этой погрешности превращает число в заслуживающий доверия ответ. Есть три разных способа говорить о погрешности, и их нельзя путать.

Абсолютная погрешность. Абсолютная погрешность это величина разрыва между приближением и истинным значением: |L - L_n|. Для параболы при n = 4 абсолютная погрешность равна |4.6467837624 - 4.6267234873| = 0.0200602751 единиц. Для синусоиды при n = 8 она равна |3.8201977890 - 3.8201978627| = 0.0000000736 единиц. Абсолютная погрешность измеряется в тех же единицах, что и сама длина.

Относительная погрешность. Относительная погрешность делит абсолютную погрешность на величину истинного значения и обычно выражается в процентах: |L - L_n| / |L| × 100. Для параболы при n = 4 относительная погрешность равна 0.4317%. Для кубической кривой при n = 4 она равна 0.9585%. Для синусоиды при n = 8 она равна 0.00000193%. Относительная погрешность это честный способ сравнивать точность задач разного масштаба: погрешность 0.02 единиц мало значит для длины 100 единиц и много значит для длины 0.1 единиц.

Индикатор погрешности на основе измельчения. В реальной работе истинное значение неизвестно, поэтому абсолютную погрешность нельзя вычислить напрямую. Вместо этого сравнивают два приближения: разность |M_2n - M_n| показывает, насколько ответ сдвинулся при измельчении разбиения. Для параболы |M₈ - M₄| = |4.6417306581 - 4.6267234873| = 0.0150071707 единиц. Эта разность это практический предупреждающий сигнал: если она велика, ответ ещё движется, и n следует увеличить. Это индикатор, а не гарантированная оценка: истинная погрешность может быть больше или меньше разности измельчения, и обычно она составляет около трети от неё для правила средних точек, потому что погрешность убывает в четыре раза при каждом удвоении.

Теоретическая оценка погрешности. Для достаточно гладкой подынтегральной функции g(x) погрешность составного правила средних точек удовлетворяет неравенству:

Для этой оценки нужен максимум второй производной подынтегральной функции на интервале, который часто трудно найти точно. Для подынтегральной функции параболы g(x) = √(1 + 4x²) вторая производная равна g''(x) = 4 / (1 + 4x²)^3/2, а её наибольшее значение на [0, 2] равно 4 при x = 0. Поэтому оценка при n = 4 равна (8 / (24 × 16)) × 4 = 0.0833333333 единиц, а фактическая погрешность 0.0200602751 единиц спокойно укладывается в неё. При n = 8 оценка равна 0.0208333333 единиц. Оценка честна, но неточна: здесь она завышает реальную погрешность в четыре раза, что типично.

Когда оценка не применяется. Оценка погрешности предполагает, что подынтегральная функция гладкая и её вторая производная ограничена. Функция с изломом, скачком или особенностью нарушает это предположение, и оценка либо не существует, либо бесполезно велика. Калькулятор не выдумывает оценку, когда не может её подтвердить. В таком случае он показывает разность измельчения и объясняет её пределы, и это честный ответ.

Почему гладкость имеет значение. Поведение погрешности 1 / n² зависит от того, что подынтегральная функция достаточно гладкая для работы компенсации в средних точках. Гладкая подынтегральная функция, такая как √(1 + 4x²), сходится чисто. Подынтегральная функция с острым углом сходится медленнее, а подынтегральная функция с особенностью может сходиться очень медленно или не сходиться вовсе. Если удвоение n почти не меняет ответ, функция может быть гладкой и уже точной. Если удвоение n сильно меняет ответ, либо функция грубая, либо n всё ещё слишком мало.

Экстраполяция Ричардсона это простой приём, который превращает два приближения по правилу средних точек в одну гораздо лучшую оценку. Она использует известный характер погрешности правила средних точек, чтобы сократить главный член погрешности.

Для гладкой подынтегральной функции приближение по правилу средних точек при n подынтервалах имеет погрешность, которая ведёт себя как C / n², где C это постоянная, зависящая от функции, но не от n. Удвоение числа подынтервалов до 2n делит погрешность на четыре. Это даёт два уравнения с двумя неизвестными, истинным значением L и постоянной C:

Решение относительно L исключает C и даёт экстраполированную оценку Ричардсона:

Разобранный пример. Для параболы y = x² на [0, 2] M₄ = 4.6267234873 единиц, а M₈ = 4.6417306581 единиц. Оценка Ричардсона равна L_R = (4 × 4.6417306581 - 4.6267234873) / 3 = 4.6467330483 единиц. Истинное значение равно 4.6467837624 единиц, поэтому погрешность экстраполяции составляет всего 0.0000507141 единиц, что гораздо меньше исходной погрешности M₄, равной 0.0200602751 единиц. Экстраполяция дала точность гораздо большего n, используя только вычисления при n = 4 и n = 8.

Разобранный пример с синусоидой. Для y = sin(x) на [0, π] M₈ = 3.8201978627 единиц, а M₁₆ = 3.8201977890 единиц. Оценка Ричардсона равна L_R = (4 × 3.8201977890 - 3.8201978627) / 3 = 3.8201977645 единиц при эталонном значении 3.8201977890 единиц, то есть погрешность составляет всего 0.0000000245 единиц.

Разобранный пример с кубической кривой. Для y = x³ на [0, 1] M₄ = 1.5330289642 единиц, а M₈ = 1.5441590533 единиц. Оценка Ричардсона равна L_R = (4 × 1.5441590533 - 1.5330289642) / 3 = 1.5478690830 единиц при эталонном значении 1.5478656547 единиц, то есть погрешность равна 0.0000034283 единиц. Исходная погрешность M₄ составляла 0.0148366905 единиц, поэтому экстраполяция сократила её примерно в 4327.7 раза.

Чем экстраполяция Ричардсона не является. Экстраполированное значение это уточнённая оценка, а не результат правила средних точек и не гарантированно точный ответ. Она предполагает, что погрешность действительно ведёт себя как C / n², что верно для гладких подынтегральных функций и неверно для грубых. Если у подынтегральной функции есть излом или особенность, экстраполяция может оказаться хуже исходного приближения, поэтому калькулятор всегда сначала показывает исходный результат правила средних точек, а экстраполированное значение помечает отдельно. Экстраполяция Ричардсона это также идея, лежащая в основе интегрирования Ромберга, которое повторяет тот же приём на нескольких уровнях, чтобы достичь очень высокой точности.

Некоторые кривые нельзя записать как y = f(x). Окружность, например, даёт два значения y для большинства значений x, поэтому ни одна функция её не описывает. Параметрические кривые решают это: обе координаты задаются как функции параметра t: x = x(t) и y = y(t). Когда t пробегает от t₁ до t₂, точка (x(t), y(t)) описывает кривую.

Для параметрической кривой крошечный элемент длины получается из движения по обеим координатам: ds = √(dx² + dy²). Запись dx = x'(t) dt и dy = y'(t) dt даёт параметрический интеграл длины дуги:

Величина под квадратным корнем это скорость v(t) = √([x'(t)]² + [y'(t)]²), то есть скорость, с которой точка движется вдоль кривой при изменении t. Правило средних точек применяется к интервалу параметра точно так же, как к интервалу x: разделите [t₁, t₂] на n равных частей шириной Δt = (t₂ - t₁) / n, возьмите скорость в средних значениях параметра t_i* = t₁ + (i + ½)Δt и сложите вклады:

Разобранный пример: четверть окружности. Возьмите x(t) = 5 cos(t) и y(t) = 5 sin(t) на [0, π/2] при n = 4. Производные равны x'(t) = -5 sin(t) и y'(t) = 5 cos(t). Скорость равна v(t) = √(25 sin²(t) + 25 cos²(t)) = √(25) = 5, то есть постоянна, потому что sin²(t) + cos²(t) = 1. Ширина параметра равна Δt = (π/2 - 0) / 4 = π/8 = 0.3926990817, а средние значения параметра равны π/16, 3π/16, 5π/16 и 7π/16. Каждое значение скорости равно 5, поэтому каждый вклад равен 0.3926990817 × 5 = 1.9634954085, а общая сумма равна 4 × 1.9634954085 = 7.8539816340 единиц, что в точности равно 5π/2. Правило средних точек здесь точно при любом n, потому что скорость постоянна.

Почему параметрическая форма имеет значение. Многие реальные кривые естественно параметрические: окружности, эллипсы, спирали, кривые Безье в программах дизайна и траектории движущихся объектов в физике и робототехнике. Явная форма y = f(x) не может описать кривую, которая поворачивает обратно на себя, а параметрическая форма может. Калькулятор обрабатывает обе формы одним и тем же механизмом средних точек, той же таблицей и тем же анализом сходимости. Подробнее о параметрических кривых смотрите в калькуляторе длины дуги параметрической кривой.

Предупреждение о параметре. Параметр t это не то же самое, что длина дуги, если скорость не равна 1. Равные шаги по t дают равные шаги вдоль кривой только при постоянной скорости. Для общей параметрической кривой правило средних точек всё равно работает правильно, потому что оно интегрирует скорость, но прямоугольники равны по t, а не по расстоянию вдоль кривой. Это нормально: формула учитывает это через v(t).

Правило средних точек полезно везде, где интеграл нужно оценить численно, а длина дуги это только одна из его задач. Его простота, точность второго порядка и терпимость к проблемам на концах делают его рабочей лошадкой во многих областях.

Численное исчисление и обучение. Правило средних точек часто становится первым методом численного интегрирования, с которым знакомятся студенты, потому что прямоугольники легко рисовать, а идею средней точки легко понять. Оно учит суммам Римана, сходимости и анализу погрешности в одном пакете. Таблица и диаграммы этого калькулятора построены именно для этой цели: каждое число в ответе можно проследить до одного прямоугольника.

Инженерное дело. Инженеры оценивают длины изогнутых деталей, кабелей, труб и поворотов дорог по математическим моделям. Когда кривая задана многочленом или тригонометрической функцией, интеграл длины дуги часто не имеет замкнутой формы, поэтому численный метод это единственный вариант. Правило средних точек даёт быстрый ответ с ясной оценкой погрешности, а увеличение n уточняет его настолько, насколько нужно проекту.

Физика. Физики интегрируют вдоль путей, чтобы вычислить работу, пройденное расстояние и величины полей. Частица, движущаяся по изогнутому пути, проходит расстояние, равное длине дуги её траектории, а правило средних точек оценивает это расстояние по данным о скорости. Тот же метод оценивает интегралы в термодинамике, оптике и квантовой механике везде, где аналитическая первообразная недостижима.

Компьютерная графика и САПР. Кривые на экране рисуются как множество коротких прямых отрезков, а общая длина кривой оценивается сложением этих отрезков, что является формой численного интегрирования. Контуры шрифтов, траектории анимации и траектории инструментов САПР зависят от точного измерения длины кривой. Правило средних точек даёт обоснованный способ измерять эти длины вместо догадок по картинке.

Измерение кривых и робототехника. Рука робота, перемещающая захват по дуге окружности, станок с ЧПУ, вырезающий изогнутый край, или геодезист, измеряющий изогнутую границу, все нуждаются в длине кривой, известной только через её уравнение или через отдельные точки. Правило средних точек превращает эти значения в длину с оценкой погрешности, которая говорит, насколько ответу можно доверять.

Научное моделирование и вычислительная геометрия. Модели в биологии, финансах и климатологии интегрируют функции, известные только в отдельных точках. Правило средних точек это естественный выбор, когда данные это значения, взятые в центрах интервалов, что часто встречается в показаниях датчиков и результатах моделирования. В вычислительной геометрии оно измеряет периметры изогнутых фигур и длины сплайн-кривых.

Когда его не использовать. Правило средних точек это плохой выбор, когда у подынтегральной функции есть особенность внутри интервала, когда функция разрывна или когда очень высокая точность нужна при малом числе вычислений: правило Симпсона или адаптивный метод справятся там лучше. Это также неверный инструмент, когда сам интеграл не является вопросом, например когда замкнутая форма длины дуги уже известна. Четверть окружности радиуса 5 имеет дугу ровно 5π/2 = 7.8539816340 единиц, и никакой численный метод не улучшит этот результат.

Большинство неверных ответов по правилу средних точек происходят из короткого списка ошибок. Проверка каждой из них занимает секунды и спасает от большинства ошибок.

1. Использование концов вместо средних точек. Выборка на левом или правом краю каждого подынтервала вычисляет левую или правую сумму Римана, а не правило средних точек. Для параболы y = x² на [0, 2] при n = 4 сумма по левым концам равна 3.75 единиц, а сумма по правым концам 5.75 единиц, тогда как правило средних точек даёт 4.6267234873 единиц. Только значение в средней точке является ответом правила средних точек.
2. Забыли продифференцировать функцию. Подынтегральная функция равна g(x) = √(1 + [f'(x)]²), а не √(1 + [f(x)]²). Для y = x² использование f(x) вместо f'(x) даёт интеграл от √(1 + x⁴), совершенно другое число. Всегда сначала дифференцируйте.
3. Неверная подынтегральная функция длины дуги. Подынтегральная функция должна быть √(1 + [f'(x)]²). Если убрать 1, забыть квадрат или взять квадратный корень не от той величины, ответы будут неверными. Для прямой с наклоном 3 подынтегральная функция это постоянная √(10) = 3.1622776602, а не 3 и не √(1 + 3).
4. Неправильное деление интервала. Ширина равна Δx = (b - a) / n, где сверху стоит полная длина интервала. Деление на n + 1 или использование b / n вместо (b - a) / n сдвигает каждую среднюю точку и портит ответ. Для [0, 2] при n = 4 ширина равна 0.5, а не 0.4 и не 0.6666666667.
5. Неверное число средних точек. Средних точек ровно n, по одной на подынтервал, а не n + 1. Концы не являются средними точками. Подсчёт n + 1 точек выборки означает, что один подынтервал посчитан дважды или конец взят по ошибке.
6. Путаница исходной кривой с графиком подынтегральной функции. Прямоугольники средних точек находятся под подынтегральной функцией g(x), а не под исходной кривой f(x). Если нарисовать их под f(x) и сложить их площади, получится площадь под f(x), а не длина дуги. Два графика это разные кривые, и их нужно показывать отдельно.
7. Округление промежуточных значений. Округление средних точек, значений производной или значений подынтегральной функции до сложения вносит малые погрешности, которые накапливаются по многим прямоугольникам. Сохраняйте полную точность на протяжении всего расчёта и округляйте только итоговый отображаемый результат. Калькулятор делает это автоматически.
8. Убеждение, что каждое приближение точно. Правило средних точек даёт приближение, а не истинное значение, за исключением особых случаев, таких как прямые и параметрические кривые с постоянной скоростью. Для параболы L₄ = 4.6267234873 единиц близко к истинным 4.6467837624 единиц, но не равно им. Всегда сообщайте результат как приближение и приводите рядом оценку погрешности.
9. Путаница оценочной погрешности с гарантированной оценкой погрешности. Разность измельчения |M_2n - M_n| это индикатор того, насколько ответ ещё движется, а не гарантия того, что истинная погрешность меньше этой разности. Теоретическая оценка ((b - a)³ / (24n²)) × max|g''(x)| это гарантия, но только когда известен максимум второй производной и подынтегральная функция гладкая.
10. Игнорирование ограничений области определения функции. Такая функция, как √(x), не определена при отрицательных x, а ln(x) не определён при x ≤ 0. Если средняя точка попадает вне области определения, значение подынтегральной функции не существует, и приближение бессмысленно. Калькулятор в таком случае сообщает об ошибке, а не выдумывает значение.
11. Неверные параметрические производные. Для параметрической кривой нужны обе производные x'(t) и y'(t), а скорость равна √([x'(t)]² + [y'(t)]²). Использование только одной из них или дифференцирование не по той переменной даёт неверную скорость и неверную длину. Для четверти окружности x'(t) = -5 sin(t) и y'(t) = 5 cos(t), и обе нужны, чтобы получить постоянную скорость 5.
12. Отношение к негладкой или особой функции так, как будто стандартные предположения о погрешности применяются автоматически. Поведение погрешности 1 / n² и теоретическая оценка предполагают гладкую подынтегральную функцию. Функция с изломом, скачком или особенностью внутри интервала может сходиться гораздо медленнее, а экстраполяция Ричардсона может не сработать. Проверьте таблицу сходимости: если удвоение n не сокращает изменение примерно в четыре раза, стандартные предположения для этой функции не выполняются.

Итоговая проверка, которая ловит большинство из них. Сверьте результат с длиной хорды и эталонным значением. Приближение должно быть не короче хорды, оно должно стремиться к эталонному значению при росте n, сумма строк таблицы должна равняться главному результату, а для прямой или кривой с постоянной скоростью оно должно равняться точному значению. Если какая-то из этих проверок не проходит, причина в одной из двенадцати ошибок выше.

Эти шесть примеров охватывают основные режимы калькулятора и типы функций. Каждый из них вычислен тем же механизмом средних точек, который питает инструмент, и каждый результат независимо проверен по высокоточному эталону адаптивного метода Симпсона.

Пример 1: Парабола, y = x² на [0, 2], n = 4. Производная равна f'(x) = 2x, ширина равна Δx = 0.5, а средние точки равны 0.25, 0.75, 1.25 и 1.75. Значения подынтегральной функции равны 1.1180339887, 1.8027756377, 2.6925824036 и 3.6400549446, их сумма равна 9.2534469747, а L₄ = 0.5 × 9.2534469747 = 4.6267234873 единиц. Эталонное значение равно 4.6467837624 единиц, поэтому абсолютная погрешность равна 0.0200602751 единиц, а относительная погрешность 0.4317%. Удвоение до n = 8 даёт 4.6417306581 единиц, а оценка Ричардсона из этих двух равна 4.6467330483 единиц.

Пример 2: Прямая, y = 3x + 1 на [0, 4], n = 4. Производная равна f'(x) = 3, поэтому подынтегральная функция это постоянная g(x) = √(10) = 3.1622776602. Ширина равна Δx = 1, средние точки равны 0.5, 1.5, 2.5 и 3.5, а каждое значение подынтегральной функции равно 3.1622776602. Сумма равна 4 × 3.1622776602 = 12.6491106407, а L₄ = 1 × 12.6491106407 = 12.6491106407 единиц, что в точности равно 4√(10). Правило средних точек здесь точно при любом n, потому что подынтегральная функция постоянна и каждый прямоугольник имеет правильную высоту.

Пример 3: Горизонтальная прямая, y = 2 на [1, 5], n = 4. Производная равна f'(x) = 0, поэтому подынтегральная функция это постоянная g(x) = √(1 + 0) = 1. Ширина равна Δx = 1, средние точки равны 1.5, 2.5, 3.5 и 4.5, а каждый вклад равен 1 × 1 = 1. Итого L₄ = 4.0000000000 единиц, что в точности равно горизонтальному расстоянию 5 - 1 = 4. У горизонтальной прямой нет вертикального изменения, поэтому её длина дуги просто равна длине интервала.

Пример 4: Параметрическая четверть окружности, x(t) = 5 cos(t), y(t) = 5 sin(t) на [0, π/2], n = 4. Производные равны x'(t) = -5 sin(t) и y'(t) = 5 cos(t), а скорость это постоянная v(t) = 5. Ширина параметра равна Δt = π/8 = 0.3926990817, средние значения параметра равны π/16, 3π/16, 5π/16 и 7π/16, а каждый вклад равен 0.3926990817 × 5 = 1.9634954085. Итого L₄ = 7.8539816340 единиц, ровно 5π/2. Полная окружность составила бы 10π = 31.4159265359 единиц, а полуокружность 5π = 15.7079632679 единиц.

Пример 5: Кубическая кривая, y = x³ на [0, 1], n = 4. Производная равна f'(x) = 3x², ширина равна Δx = 0.25, а средние точки равны 0.125, 0.375, 0.625 и 0.875. Значения подынтегральной функции равны 1.0010980300, 1.0853471867, 1.5405489332 и 2.5051217067, их сумма равна 6.1321158566, а L₄ = 0.25 × 6.1321158566 = 1.5330289642 единиц. Эталонное значение равно 1.5478656547 единиц, поэтому абсолютная погрешность равна 0.0148366905 единиц, а относительная погрешность 0.9585%. Удвоение до n = 8 даёт 1.5441590533 единиц, а оценка Ричардсона равна 1.5478690830 единиц.

Пример 6: Тригонометрическая функция, y = sin(x) на [0, π], n = 8. Производная равна f'(x) = cos(x), ширина равна Δx = π/8 = 0.3926990817, а восемь средних точек это π/16, 3π/16, 5π/16, 7π/16, 9π/16, 11π/16, 13π/16 и 15π/16. Значения подынтегральной функции равны 1.4006926023, 1.3005159423, 1.1439660326, 1.0188524102, 1.0188524102, 1.1439660326, 1.3005159423 и 1.4006926023. Их сумма равна 9.7280539749, а L₈ = 0.3926990817 × 9.7280539749 = 3.8201978627 единиц. Эталонное значение равно 3.8201977890 единиц, поэтому абсолютная погрешность составляет всего 0.0000000736 единиц, а относительная погрешность 0.00000193%. У этого интеграла нет элементарной первообразной, поэтому правило средних точек и эталон это практичный способ получить ответ.

Нажмите любую кнопку предустановки в калькуляторе, чтобы загрузить один из этих примеров и увидеть полное пошаговое решение, живые диаграммы, таблицу расчёта и сравнение сходимости для этих входных данных.

Длина дуги по правилу средних точек встречается везде, где изогнутое расстояние нужно измерить по формуле, а не рулеткой. Один и тот же расчёт служит студентам, инженерам и дизайнерам.

Дороги, треки и кабели. Изогнутая дорога, поворот беговой дорожки или кабель, провисший между двумя опорами, все имеют длины, которые следуют кривой. Когда кривая моделируется функцией, правило средних точек оценивает её длину по модели, а оценка погрешности говорит, насколько ответу можно доверять. Кабель, смоделированный как y = x² на [0, 2], требует 4.6467837624 единиц кабеля, а не 4.4721359550 единиц прямого промежутка.

Производство и листовой металл. Изогнутая панель, согнутая труба или завёрнутый край начинаются как плоский материал известной длины. Правило средних точек переводит изогнутый проект в плоскую длину для резки, а таблица сходимости подтверждает длину с точностью, нужной мастерской. Четверть окружности радиуса 5 требует 7.8539816340 единиц материала вдоль кривой.

Архитектура. Арки, купола и изогнутые лестницы измеряются по их дугам. Полукруглая арка радиуса 5 требует 5π = 15.7079632679 единиц отделки вдоль кривой. Для не круглых арок, смоделированных параболами или синусоидами, правило средних точек даёт длину отделки напрямую из уравнения проекта.

Физика и движение. Расстояние, которое частица проходит по изогнутому пути, равно длине дуги её траектории. Груз маятника, снаряд или спутник на эллиптической орбите движутся по кривым, длины которых важны для энергии, времени и расчётов топлива. Правило средних точек оценивает эти расстояния по уравнениям движения.

Компьютерная графика, анимация и игры. Объекты, движущиеся по изогнутым путям, движения камеры и анимированные персонажи следуют сплайнам и параметрическим кривым. Расстояние вдоль пути решает скорость, время и расстановку кадров. Правило средних точек измеряет эти длины путей по уравнениям кривых, а параметрический режим обрабатывает кривые, которые поворачивают обратно.

Медицина и биология. Кровеносные сосуды, позвоночник и изогнутые поверхности органов измеряются по сканам и моделям. Длина дуги вдоль центральной линии сосуда помогает планировать стенты и катетеры, а правило средних точек превращает смоделированную центральную линию в длину с оценкой погрешности.

Связанные калькуляторы на этом сайте. Для точных длин дуг названных кривых попробуйте калькулятор длины дуги параболы, калькулятор длины дуги эллипса, калькулятор длины дуги полуокружности, калькулятор длины дуги окружности и калькулятор длины дуги сектора. Для пошагового точного интегрирования попробуйте калькулятор длины дуги с шагами, калькулятор длины дуги параметрической кривой, калькулятор длины дуги полярной кривой и калькулятор длины дуги в 3D. Для других численных методов попробуйте калькулятор правила трапеций и калькулятор правила Симпсона.

Эти таблицы собирают поведение сходимости двух основных разобранных примеров в одном месте, чтобы вы могли увидеть, как быстро правило средних точек устанавливается и как экстраполяция Ричардсона улучшает ответ. Все значения вычислены механизмом средних точек, а эталон взят из адаптивного интегрирования Симпсона.

Парабола y = x² на [0, 2], эталонное значение 4.6467837624 единиц.

Синусоида y = sin(x) на [0, π], эталонное значение 3.8201977890 единиц.

Как читать таблицы. Погрешность параболы падает примерно в четыре раза при каждом удвоении, с 0.0200602751 при n = 4 до 0.0050531044 при n = 8 и до 0.0012632302 при n = 16, что и есть ожидаемое поведение 1 / n². Синусоида сходится ещё быстрее, потому что её подынтегральная функция почти плоская возле середины интервала: к n = 8 погрешность уже равна 0.0000000736 единиц, а к n = 16 приближение совпадает с эталоном до всех десяти показанных знаков после запятой. Строки Ричардсона показывают ту же идею за один шаг: объединение двух приближений превосходит каждое из них по отдельности, с погрешностью 0.0000507141 единиц для параболы и 0.0000000245 единиц для синусоиды. Используйте эти таблицы как эталон: если ваш собственный расчёт для одного из этих примеров не совпадает с таблицей, проверьте производную, средние точки и подынтегральную функцию, прежде чем доверять ответу.