Как найти длину окружности эллипса - Исправление недочетов и поиск решений вместе с Examum.ru

Перейти к содержанию

Длина эллипса (овала)

На чтение 1 мин

Эллипс – это множество точек плоскости, для которых сумма расстояний до фокусов эллипса постоянна и больше расстояния между фокусами.

a, b – полуоси эллипса

О – центр эллипса

Длина эллипса (L) равна произведению суммы его полуосей (a, b) на число π:

Вам также может понравиться

Дуга – это часть окружности, отсекаемая хордой.

0156

Окружность – геометрическое место точек, равноудаленных

0124

Окружность – геометрическое место точек, равноудаленных

0125

Многоугольник – это геометрическая фигура, которая

0141

Шестиугольник – это многоугольник, который имеет шесть углов.

0442

Пятиугольник – это многоугольник, который имеет пять углов.

0349

Трапеция – это четырехугольник, у которого параллельна

088

Квадрат – это параллелограмм, у которого все углы и

091

Источник

При помощи нашего калькулятора вы легко сможете узнать периметр эллипса.

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

P= 2π

√

(a²+b²)/2
формула для вычисления P – периметра по полуосям a, b

P= 2π

√

(A²+B²)/8
формула для вычисления P – периметра осям A, B

Где P – периметр, a и b – полуоси, A и B – оси эллипса, π – число Пи которое всегда примерно равно 3,14.

Источник

Калькулятор периметра эллипса

Введите длину большой и малой полуосей эллипса, укажите точность расчета и нажмите «Посчитать». Калькулятор выполнит расчет периметра эллипса (расчет приблизительный).

Калькулятор

Эллипсом называется множество всех точек плоскости, для каждой из которых сумма расстояний до двух данных точек F1 и F2 этой плоскости есть величина постоянная, больше расстояния между F1 и F2.

Точки F1 и F2 называют фокусами эллипса, а расстояние между ними — фокусным расстоянием.

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

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

Точка пересечения большой и малой осей эллипса называется его центром.

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

Формулу периметра эллипса нельзя выразить при помощи простейших функций.

Что такое эллипс: формула длины окружности эллипса

В астрономии, когда рассматривают движение космических тел по орбитам, часто применяют понятие «эллипс», поскольку их траектории характеризуются именно этой кривой. Рассмотрим в статье вопрос, что представляет собой отмеченная фигура, а также приведем формулу длины эллипса.

Что такое эллипс?

Согласно математическому определению, эллипс — это замкнутая кривая, для которой сумма расстояний от любой ее точки до двух других определенных точек, лежащих на главной оси, и носящих название фокусов, является постоянной величиной. Ниже приведен рисунок, который поясняет это определение.

Вам будет интересно: Возраст зрелости у ребенка. Признаки, психология, акселерация

На рисунке сумма расстояний PF’ и PF равна 2 * a, то есть PF’ + PF = 2 * a, где F’ и F — фокусы эллипса, «a» — длина его большой полуоси. Отрезок BB’ называется малой полуосью, а расстояние CB = CB’ = b — длина малой полуоси. Здесь точка C определяет центр фигуры.

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

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

Формула длины окружности эллипса

Хотя рассматриваемая фигура является достаточно простой, длину ее окружности точно можно определить, если вычислить так называемые эллиптические интегралы второго рода. Однако, индусский математик-самоучка Рамануджан еще в начале XX века предложил достаточно простую формулу длины эллипса, которая приближается к результату отмеченных интегралов снизу. То есть рассчитанное по ней значение рассматриваемой величины будет немного меньше, чем реальная длина. Эта формула имеет вид: P ≈ pi * [3 * (a+b) — √((3 * a + b) * (a + 3 * b))], где pi = 3,14 — число пи.

Например, пусть длины двух полуосей эллипса будут равны a = 10 см и b = 8 см, тогда его длина P = 56,7 см.

Каждый может проверить, что если a = b = R, то есть рассматривается обычная окружность, тогда формула Рамануджана сводится к виду P = 2 * pi * R.

Отметим, что в школьных учебниках часто приводится другая формула: P = pi * (a + b). Она является более простой, но и менее точной. Так, если ее применить для рассмотренного случая, то получим значение P = 56,5 см.

Периметр эллипса

Эллипс — замкнутая кривая на плоскости, которая может быть получена как пересечение плоскости и кругового цилиндра или как ортогональная проекция окружности на плоскость. Окружность является частным случаем эллипса.

Приближённая формула для нахождения периметра эллипса:

, где a — большая полуось, где b — малая полуось эллипса.

Максимальная погрешность этой формулы

0,63% при эксцентриситете эллипса

0,988 (соотношение осей

1/6,5). Погрешность всегда положительная.

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

Онлайн калькуляторы

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

Актуальная информация

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

источники:

http://1ku.ru/obrazovanie/10686-chto-takoe-ellips-formula-dliny-okruzhnosti-ellipsa/

http://calculatorium.ru/math/ellipse-perimeter

Источник

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

Калькулятор

Давайте вначале посчитаем, оценим.

Если есть вопросы, предложения по калькулятору или заметили ошибку, буду очень рад обратной связиx

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

Эллипс:

a:
b:

Параметры дуги:

Get a better browser, bro…

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

Длина дуги, как сумма хорд

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

Рис.1. Сумма хорд

Формула такова:

(1) Сумма хорд

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

В коде выглядит так:

//*****************************************************************

// Нахождение длины дуги эллипса, заданного прямоугольником ARect

// По аналогии с GDI+, углы в градусах, по часовой. Хордами

// ∑ √ (Δx² + Δy²)

//*****************************************************************

function CalcArcEllipseLen1(ARect : TxRect; // эллипс

startAngle : Extended; // старт.угол

sweepAngle : Extended; // угол дуги

deltaAngle : Extended=0.01 // дискрета

): Extended;

var tmp : Extended; // наращиваемый угол

pnt : TxPoint; // новые координаты по tmp

val : TxPoint; // координаты в предыдущей итерации

dw : Extended; // убывающий параметр цикла

l : extended; // длина текущей хорды

begin

// переводим градусы в радианы

startAngle := startAngle * pi/180;

sweepAngle := sweepAngle * pi/180;

deltaAngle := deltaAngle * pi/180;

// инициализируем

tmp := startAngle;

dw := sweepAngle;

result := 0;

// стартовая координата

val := CalcEllipsePointCoord(ARect,tmp);

repeat

// уменьшаем параметр цикла

dw := dw — deltaAngle;

// определяем значение угла

if dw < 0 then tmp := tmp + (deltaAngle+dw)

else tmp := tmp + deltaAngle;

// считаем новую точку на эллипсе

pnt := CalcEllipsePointCoord(ARect,tmp);

// длина хорды

l := sqrt(sqr(val.X—pnt.x)+sqr(val.Y—pnt.y));

// сумма, интеграл

result := result+l;

val := pnt;

until dw < 0;

end;

Проверим с планетарным размахом. По последним данным Международной Службы Вращения Земли (IERS — 1996) большая экваториальная полуось Земли равна 6 378 136,49 м, а полярная малая — 6 356 751,75 м.

Посчитаем периметр меридиана каким-нибудь онлайн-калькулятором, получаем 40 007 859.543 (некоторые могут дать другое число, т.к. используют приближенные формулы для вычисления периметра).

Представленная выше функция за 109 милисекунд выдала результат 40 007 996.265 при дельте 0.001. Это нельзя назвать точным результатом.

Длина дуги, как интеграл

Длиной некоторой дуги называется тот предел, к которому стремится длина вписанной ломаной, когда длина наибольшего ее звена стремится к нулю:

(2) Предел суммы хорд при максимальной длине хорды, стремящейся к нулю

(3) Длина t-го звена (хорды) вписанной ломаной

Таким образом, длина дуги эллипса может быть описана интегральным уравнением:

(4) Интегрально-дифференциальное уравнение дуги эллипса

Используя параметрическое уравнение эллипса, приходим к уравнению:

Где t1 и t2 – параметры для начала и конца дуги параметрического уравнения эллипса. Параметром является некий угол к оси абсцисс. Что такое и как найти параметр для угла эллипса подробно изложено тут.

Зная, что (cos t)’ = — sin t, (sin t)’ = cos t (подробный вывод приведен тут и тут), получаем следующую формулу:

(5) Длина дуги эллипса

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

// Получить разницу между углами A2-A1

function CalcSweepAngle (A1,A2 : Extended;

ANormalize : boolean = true) : Extended;

begin

result := A2—A1;

if (result < 0) and (ANormalize) then

result := 2*pi + result;

end;

//*****************************************************************

// Нахождение длины дуги эллипса, заданного через a и b

// Сумма от корня суммы квадратов производных, через sin/cos по ∂t

// ∑ [√(∂x²/∂t² + ∂y²/∂t²)] ∂t

//*****************************************************************

function CalcArcEllipseLen2(

a,b : extended; // полуоси

startAngle : Extended; // старт.угол

sweepAngle : Extended; // угол откл.

deltaAngle : Extended=0.01 // дискрета

): Extended;

var sn,cs: Extended;// синус, косинус

tmp : Extended; // наращиваемый угол

bet : Extended; // старый параметр

dw : Extended; // убывающий параметр цикла

t : Extended; // новый параметр

l : extended; // значение функции в t

begin

// переводим градусы в радианы

startAngle := startAngle * pi/180;

sweepAngle := sweepAngle * pi/180;

deltaAngle := deltaAngle * pi/180;

// инициализируем

tmp := startAngle;

dw := sweepAngle;

result := 0;

// параметр для начального угла

SinCos(tmp,sn,cs);

bet := ArcTan2(a*sn, b*cs);

repeat

// уменьшаем параметр цикла

dw := dw — deltaAngle;

// определяем значение угла

if dw < 0 then

tmp := tmp + (deltaAngle+dw)

else tmp := tmp + deltaAngle;

// находим параметр для угла

SinCos(tmp,sn,cs);

t := ArcTan2(a*sn, b*cs);

// разница между параметрами dt=t-bet

bet := CalcSweepAngle(bet,t);

// эллиптический интеграл II рода

SinCos(t,sn,cs);

l := sqrt ( sqr(a)*sqr(sn) + sqr(b)*sqr(cs))*bet;

// запоминаем параметр

bet := t;

// сумма, интеграл

result := result+l;

until dw <= 0;

end;

Результат при dt = 0.001 равен 40 007 859,543 за 109 милисекунд. Отлично!

Длина дуги через эксцентриситет

Это еще не окончательный вид уравнения. В ряде интересных параметров для эллипса есть такая числовая характеристика, показывающая степень отклонения эллипса от окружности, как эксцентриситет. Формула для эллипса:

(6) Эксцентриситет эллипса через полуоси

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

Выразим b² = a²(1 — e²), подставим в формулу (5), помним, что sin²t + cos²t = 1 (справочник [1]), убираем a² за знак корня, и, как постоянную величину, за знак интеграла тоже, получаем:

(7) Окончательный вид уравнения для длины дуги эллипса

Пишем функцию, ожидаем более шустрого выполнения:

//*****************************************************************

// Нахождение длины дуги эллипса, заданного через a и b

// По аналогии с GDI+, углы в градусах, по часовой. Считаем E(t,e)

// a ∑ [√(1-e²cos²t)] ∂t

//*****************************************************************

function CalcArcEllipseLen3(

a,b : extended; // полуоси

startAngle : Extended; // старт.угол

sweepAngle : Extended; // угол откл.

deltaAngle : Extended=0.01 // дискрета

): Extended;

var sn,cs: Extended;// синус, косинус

tmp : Extended; // наращиваемый угол

bet : Extended; // старый параметр

dw : Extended; // убывающий параметр цикла

ex : Extended; // эксцентриситет

t : Extended; // новый параметр

l : extended; // значение функции в t

begin

// переводим градусы в радианы

startAngle := startAngle * pi/180;

sweepAngle := sweepAngle * pi/180;

deltaAngle := deltaAngle * pi/180;

// считаем эксцентриситет

ex := 1 —sqr(b)/sqr(a);

// инициализируем

tmp := startAngle;

dw := sweepAngle;

result := 0;

// параметр для начального угла

SinCos(tmp,sn,cs);

bet := ArcTan2(a*sn, b*cs);

repeat

// уменьшаем параметр цикла

dw := dw — deltaAngle;

// определяем значение угла

if dw < 0 then

tmp := tmp + (deltaAngle+dw)

else tmp := tmp + deltaAngle;

// находим параметр для угла

SinCos(tmp,sn,cs);

t := ArcTan2(a*sn, b*cs);

// разница между параметрами dt=t-bet

bet := CalcSweepAngle(bet,t);

// эллиптический интеграл II рода

l := sqrt (1—ex*sqr(cos(t)))*bet;

// запоминаем параметр

bet := t;

// сумма, интеграл

result := result+l;

until dw <= 0;

result := a*result;

end;

Результат 40 007 859.543, за чуть меньшее время, 94 милисекунды.

Длина дуги через эксцентриситет с подготовкой

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

Рис.2. Параметры дуги эллипса.

То есть, вместо того, чтобы идти из точки A в точку B, мы заранее считаем параметры t1 и t2. И в цикле больше нахождением параметров не занимаемся, а только считаем очередное приращение.

Берем код последней реализации и улучшаем:

//*****************************************************************

// Нахождение длины дуги эллипса, заданного через a и b

// По аналогии с GDI+, углы в градусах, по часовой. Считаем E(t,e)

// a ∑ [√(1-e²cos²t)] ∂t

// Но параметры и угол отклонения параметра считаем заранее

//*****************************************************************

function CalcArcEllipseLen4(

a,b : extended; // полуоси

startAngle : Extended; // старт.угол

sweepAngle : Extended; // угол откл.

deltaAngle : Extended=0.01 // дискрета

): Extended;

var sn,cs: Extended;// синус, косинус

tmp : Extended; // наращиваемый угол

bet : Extended; // старый параметр

dw : Extended; // убывающий параметр цикла

ex : Extended; // эксцентриситет

l : extended; // значение функции в t

begin

// переводим градусы в радианы

startAngle := startAngle * pi/180;

sweepAngle := sweepAngle * pi/180;

deltaAngle := deltaAngle * pi/180;

// считаем эксцентриситет

ex := 1 —sqr(b)/sqr(a);

// инициализируем

tmp := startAngle;

result := 0;

// параметр для начального угла

SinCos(tmp,sn,cs);

bet := ArcTan2(a*sn, b*cs);

// считаем параметр для цикла

SinCos(tmp+sweepAngle,sn,cs);

dw := ArcTan2(a*sn, b*cs);

dw := CalcSweepAngle(bet, dw);

// для корректной работы цикла, запоминаем в бет старт.параметр

tmp := bet;

// если sweepAngle = 360, результат будет равен 0

// если sweepAngle > 360, результат будет от оставшейся разницы

if sweepAngle >= 2*pi then

dw := 2*pi + dw;

repeat

// уменьшаем параметр цикла

dw := dw — deltaAngle;

bet := tmp;

// определяем значение угла

if dw < 0 then tmp := tmp + (deltaAngle+dw)

else tmp := tmp + deltaAngle;

// разница между параметрами dt=tmp-bet

bet := CalcSweepAngle(bet,tmp);

// эллиптический интеграл II рода

l := sqrt (1—ex*sqr(cos(tmp)))*bet;

// сумма, интеграл

result := result+l;

until dw <= 0;

result := a*result;

end;

Результат — 40 007 859,543 за 30 милисекунд! Мысль явно здравая.

Длина дуги и хорошо забытые хорды

Но вернемся к сумме хорд. Что там-то не так? Казалось бы, все просто, понятно и должно работать, но результат, мягко говоря, не точен.

Работая с вещественными типами, результат зависит от «мелкости» приращения, от разрядности используемых типов. Накапливаемые погрешности, неточное представление дробного числа и прочие прелести чисел с плавающей запятой.

Тип Extended – самый точный из существующих в Delphi (и, как следствие, самый медленный). Перепишем нахождение длины дуги хордами, используя этот тип, и без вызовов внешних функций нахождения координат на эллипсе.

//*****************************************************************

// Нахождение длины дуги эллипса, заданного через a и b

// По аналогии с GDI+, углы в градусах, по часовой. Считаем хордами

// Абсолютно то же самое, что (1), но типы и точность изменены

// ∑ √ (Δx² + Δy²)

//*****************************************************************

function CalcArcEllipseLen5(

a,b : extended; // полуоси

startAngle : Extended; // старт.угол

sweepAngle : Extended; // угол откл.

deltaAngle : Extended=0.01 // дискрета

): Extended;

var sn,cs: Extended;// синус, косинус

tmp : Extended; // наращиваемый угол

dw : Extended; // убывающий параметр цикла

l : extended; // длина текущей дуги

vx,vy : extended;

px,py : Extended;

begin

// переводим градусы в радианы

startAngle := startAngle * pi/180;

sweepAngle := sweepAngle * pi/180;

deltaAngle := deltaAngle * pi/180;

// инициализируем

tmp := startAngle;

dw := sweepAngle;

result := 0;

// находим параметр (некий угол) для уравнения

SinCos(tmp,sn,cs);

SinCos (ArcTan2(a*sn, b*cs), sn, cs);

vX := a * cs;

vY := b * sn;

repeat

// уменьшаем параметр цикла

dw := dw — deltaAngle;

// определяем значение угла

if dw < 0 then tmp := tmp + (deltaAngle+dw)

else tmp := tmp + deltaAngle;

// считаем новую точку на эллипсе

SinCos(tmp,sn,cs);

SinCos (ArcTan2(a*sn,b*cs),sn,cs);

pX := a * cs;

pY := b * sn;

// длина хорды

l := sqrt(sqr(pX—vx)+sqr(pY—vy));

// сумма, интеграл

result := result+l;

vx := px;

vy := py;

until dw < 0;

end;

Результат 40 007 859.542 за 94 милисекунды. Разница в одну тысячную, и такое же время. Весьма неплохо!

Практика

Скачать исходник + исполняемый файл

Качаем, запускаем. Видим такое окно:

Рис.3. Пример: Вид при запуске

При запуске в полях полуосей находятся «земные» значения. Эллипс сильно смахивает на окружность. Можно сразу нажать «Расчет», увидеть значения полярного меридиана, и начинать уже эксперименты.

Для начала проверим корректность получаемых результатов. Я попытался найти какие-нибудь онлайн калькуляторы для подсчета длины дуги эллипса, в которых можно было бы указать стартовый угол и произвольное отклонение, либо просто пару углов, определяющих дугу. Но все калькуляторы считают периметр эллипса. А с произвольной дугой как-то… не встретились.

Картинка из вики:

Длина дуги эллипса (s) в зависимости от его параметра (θ) Автор: Illustr — собственная работа, CC BY-SA 4.0, Ссылка

Заскриним несколько моментов и посчитаем корректность:

Рис.4. Скрины расчетов

Введем a=2, b=0.5, стартовый угол равен 0, отклонение θ=5.642.

Рис.5. Подсчет длины дуги для a=2, b=0.5, стартовый угол=0, отклонение θ=5.642

Видим результат у всех один 8.055. Т.е правы все.

Аналогично поступаем с остальными тремя скринами:

отклонение θ=1.154. У нас получился результат 1.333. На скрине видим результат s=1.334. Как так? Давайте увеличим «Дельту» в 10 раз, т.е. вместо 0.001, сделаем 0,01. У всех интегральных 2) 3) 4) результат станет 1.334. В то время, как у 1) и 5), т.е. примитивно-неинтегрально-хордовых останется 1.333.

Какой результат более истинный? Смотрите сами. Уменьшение дельты, т.е. угла для подсчетов, ведет к более точному результату. На 0.001 интегральные функции выдали результат как хордовые. При более грубой дельте, интегральные чуть изменились, а хордовые верны своему результату.

Сделаем дельту очень мелкой, равной 0.00001. Результат у всех, кроме первой, тот же, 1.333.

Лично я начинаю верить 5-ой формуле.

отклонение θ= 206. У нас получился результат 4.322. На скрине результат s= 4.322. Дельта 0.001. Все отлично.
отклонение θ= 4,488. У нас получился результат 5,989. На скрине результат s= 5,989. Дельта 0.001. Все отлично.

Вывод : Формулы 2-5 работают как надо. Фаворит 5. За меньшее время, т.е. при более «грубой» дельте, находит правильный результат.

Проверим на ужасных эллипсах. Т.е. где эксцентриситет очень близок к единице.

Рис.6. Проверка на очень «сплющенном» эллипсе

Можем заметить следующее: при 0.00001 функции 2 и 3 дали результат, близкий к результату функции 4, полученный при дельте 0.001. При дельте 0.00001 функция 4 дала результат, близкий к результату функции 5. Сама же функция 5 слабо колеблется в показаниях, что при дельте в 0.001, что при 0.00001.

Аналогичную ситуацию можно пронаблюдать при сильно вытянутом эллипсе:

Рис.7. Проверка на очень «вытянутом» эллипсе

Таким образом, имеет смысл использовать функции 4 и 5. Одну, как представительницу интегрального сословия, самую быструю и более точную из них. Другую, как представительницу очевидного и простого метода, работающую, между тем, лучше своих интегральных коллег при минимальных ресурсных затратах.

Небольшая инструкция

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

Левая кнопка мыши задает конечную точку дуги, т.е. угол отклонения.