Гіперкульовий сегмент

Матеріал з testwiki
Перейти до навігації Перейти до пошуку

Шаблон:Недостатньо джерел

Приклад тримірного кульового сегмента (прозора частина). Решта частина кулі є другим кульовим сегментом.

Гіперкульовий сегмент - вибрана частина гіперкулі, за перерізу гіперплощиною, або сегмент гіперсферичної поверхні. Термін поєднує залежності пов'язані з алгебраїчним узагальненням рівнянь сфери та площини для довільної вимірності[1].

Відповідно не є геометричним об'єктом, але містить загальні властивості цілої групи геометричних об'єктів. Якщо геометричне тіло не має самостійної назви, чи вказується зв'язок із узагальненням, для його позначення додають значення вимірності відкидаючи префікс гіпер.

Гіперсферичний сегмент - гіперповерхня сферичної частини кульового сегмента[2].

Алгебраїчно фігура задається трьома параметрами, часто два з яких є геометричними а третій є вимірністю простору.

Кутові параметри

Просторова вимірність кутових параметрів складна, але їх можна цілісно виділити як незалежні у просторі параметрів.

cos(θ)=rhr=(1(ar)2)12=t

a=h(2rh)=r1t2=r|sin(θ)|

В двомірному випадку є простий зв'язок з об'ємом сектора та його конічною частиною Vdome=VsectorVcone, Vcone(r,h,n)=rhnSbase.

n=2,Vdome=r2θa(rh)

У тримірному, простий зв'язок є з супровідними тригонометричними функціями. Для більших вимірностей простору такий зв'язок менш очевидний, через те що подобою кута в радіанах є узагальнений стерадіан, який пов'язаний із плоским аналогом більш складно.

Лінійні параметри

  • Радіус гіперсферичного сегмента r - радіус гіперсфери якій належить його поверхня.
  • Основа гіперкульового сегмента — куля попередньої вимірності з радіусом a=h(2rh), що належить площині перерізу.
  • Висота гіперкульового сегмента h — найменша відстань між дотичними до нього паралельними гіперплощинами, дорівнює довжині сегменту прямої що нормально проходить через центр основи[3].

Об'єм

Вирази через радіус гіперсферичного сегмента r та висоту h, для евклідового простору вимірності n:

За обмеження аргумента виразу по висоті, отримання значень в необхідному діапазоні досягається заміною hs2rh та Vdome(r,h,n)Vball(r,n)Vdome(r,hs,n).

З використанням спеціальних функцій

гіпергеометричної функції 2F1, чи регуляризованої неповної бета-функції[4] Ix(a,b):


Vdome(r,h,n)=12Vball(r,n)Isin2θ(n+12,12)=rVball(a,n1)(2F1(12,n12;n+12;(ar)2)(rhr))

де rh, Vball(r,k)=rkπk/2Γ(k/2+1), Γ (гамма-функція), sin2θ=(ar)2=h(2rh)r2, a=h(2rh)


або без обмежень по висоті:

Vdome(r,h,n)=Vball(r,n)(12rhrΓ(1+n2)πΓ(n+12)2F1(12,1n2;32;(rhr)2))

За наявності алгебраїчних функцій у розрахунках, вирази корисні можливістю точно оцінювати залежності, в тому числі за інтегрування чи диференціювання. Розрахунково, відстань у декілька ланок операторів від бібліотечних функцій з точністю останнього біта, кінцевий вираз робить якісною апроксимацією. Але попри те інколи достатні більш прості наближення чи менші необхідні можливості.

Через кінцеву ітерацію

за k=nmod(n,2)2 кроків, F(m,t)=t1(1x2)m2dx:

Алгебраїчна форма
обмеження кроки вираз параметри
n1 кінцевий Vdome(r,h,n)=rVball(r,n1)F(n1,t) p=mod(n1,2), t=cos(θ)=rhr
i1 k F(i+2,t)=1i+3((i+2)F(i,t)t(1t2)i+22) F(p,t)={1tp=0arccos(t)p=1
i0 kp Vball(r,i+2)=2πr2i+2Vball(r,i) Vball(r,p)={1p=02rp=1

"Гонча" серед загальних виразів для невеликих значень вимірності, мінімальні вимоги до оперативної пам'яті та елементарні оператори, але попри те що має характер схильний до представлення добутком, тобто заявку на точність, містить в тілі ітерації сумування яке суттєво її зменшує, тому за "великих дистанцій", розрізнює тільки "головну доріжку" - залежність від вимірності. Оскільки вираз для ітерації нескладно перевести в диференціальний вигляд, є корисним "на швидкості" тобто за аналізу поведінки похідних.

За алгебраїчного використання недоліків пов'язаних із розрахунковою помилкою немає. Допускає незначне прискорення за узгодження початкових значень вимірності.

Приклад створеного виразу для частини залежної від висоти, у шестивимірному просторі n=6,k=3,p=1:

F(n1,t)=16(t(1t2)52+54(t(1t2)32+32(t(1t2)12+arccos(t)))).

Через кінцеві суми

У замкненій формі. Cумування для чисельних розрахунків, через зменшення помилки починається з малих доданків, що подвоює порівняно з прямим порядком кількість операцій, та потребує додаткової оперативної пам'яті через використання масивів. Або за оптимізації по розміру, зводить складність обрахунку до квадратичної.

Vdome(r,h,n)={πk1r2k(arccos(t)k!2k1t2(2k1)!!i=0k1(1)iαit2i+1)p=0,k1πkr2k+1(2k(2k+1)!!1k!i=0k(ki)(1)it2i+12i+1)p=1,k0

n1,k=np2,p=mod(n,2),t=cos(θ)=rhr,αi>0,(ki)={k(k1)(k2)(ki+1)i!i11i=0 - біноміальний коефіцієнт.

  • Для парних вимірностей значення коефіцієнтів αi можна отримати ітерацією αi+1=(ki+1)(2i+2)αi2i+3, α0=1(2k1)!!(2k)!!. Для оберненого сумування αi1=(ki)(2i+1)αi2i, перший коефіцієнт має простий вигляд αk1=12k. Зв'язок слідує із умови рівності похідних розрахункового та інтегрального виразів для об'єму (1t2)k=(2k1)!!(2k)!!tS(t)+(1t2)dS(t)dt, S(t)=i=0k1(1)iαit2i+1. Вираз вираз можна узагальнити для об'єму нульової вимірності n=0,p=0,k=0,S(t)=0,Vdome(r,t,n)=arccos(t)π. Мінімальна складність через отримання коефіцієнтів ітерацією - подвійна.
  • Для непарних вимірностей, за швидкого оцінкового розрахунку чи аналітичних виразів доречно використовувати прямий порядок. Для оберненого порядку аргумент суми інваріантний відносно заміни граничних значень індекса i=0k=i=k0.

Метод є компромісом між ітераційним шляхом та спеціальними функціями, через близькість до параметрів дозволяє робити аналітичні висновки в областях де спеціальні функції малодосліджені, та має меншу порівняно із ітераційним методом складність кінцевого алгебраїчного виразу.

Приклад створеного виразу, у шестивимірному просторі n=6,k=3,p=0:

Vdome(r,t,n)=π2r6(arccos(t)681t215(1116t1324t3+16t5)).

Рекурентний зв'язок

- в явному вигляді, походить з ітераційного розрахунку:

Vdome(r,t,i+2)=2πr2i+2(Vdome(r,t,i)rVball(r,i1)(i+1)t(1t2)i+12)

де i1,t=cos(θ)=rhr,Vball(r,0)=1

- між сусідніми вимірностями, пірамідне інтегрування через основу:

Vdome(r,h,i)=02aVdome(r(x),h(x),i1)dx

де i1,rh,r(x)=r1(axr)2,h(x)=r(x)(rh),a=h(2rh)

Інтегральна форма

За алгебраїчного використання, коли вираз є проміжним, чи для ілюстрації залежності від параметрів:

  • інтегрування через основу, циліндричне

Vdome(r,h,n)=r(n1)Vball(1,n1)0axn2((1(xr)2)12(1(ar)2)12)dx

де a=r1t2, t=cos(θ)=rhr=(1(ar)2)12.

  • інтегрування через висоту, пірамідне

Vdome(r,h,n)=rhrVball(a(x),n1)dx=rVball(r,n1)t1(1x2)n12dx

де a(x)=r2x2, t=cos(θ)=rhr.

  • інтегрування через кут

Vdome(r,h,n)=rVball(r,n1)0θsinnϕdϕ

де cos(θ)=rhr.

Асимптотична поведінка, зв'язок із площею

Оцінки нормованого об'єму сегмента: (1t2)n/22πn та 12tn2π.

За нескінченної вимірності, поведінка нормованого об'єму сферичного сегмента подібна до сходинкової функції Гевісайда H(x)=1+sgn(x)2, за сталих висоти та радіуса сферичної поверхні:

limninfVdome(r,h,n)Vball(r,n)={00<h<r,1r<h<2r, за тотожності з півкулею h=r завжди 1/2.

Характер наближення для малих висот - швидший за показниковий, у точках максимальної кривини збіжність значно менша - ступенева.

Для малих значень висоти, верхня t>0, або нижня t<0 межа:

VdomeVball=1sgn(t)2+(1t2)n/2sgn(t)2πn, sgn(t)=t|t|,t=cos(θ)=rhr,n>2,|t|>0,7.

Яка за незначних змін, може застосовуватись для більшого діапазону висот в об'єднаній оцінці.

Точки максимальної кривини нормованого об'єму, верхня оцінка: t=±ln(n)n, t=cos(θ), n>10.

Оцінка досить груба для малих вимірностей, тому обмежена як апроксимація та є радше теоретичною. Наприклад з її допомогою можна обґрунтовувати подібність поведінки нормованого об'єму сферичного сегмента до сходинкової функції. Оскільки за великої вимірності кривина показує наявність "точок зламу" біля тотожності з півкулею 1R2πnln(n)(1+2π)3/2, які знаходяться в нулі VdomeVball12πnln(n)1/2 або одиниці. Центр кривини, як і більш точні позицію та радіус краще отримувати безпосередньо, межеве значення похідної показує принципово не нульове або не нескінченне значення: ddtVdomeVball12π, тобто знаходження в "області зламу".

Безпосередньо позиція може бути визначена через значущу частину похідної кривини по косинусу:

f=(π(t2(n2)1)(1t2)n1(1+t2(2n1))(Γ(1+n/2)Γ((1+n)/2))2)=0,n>2

За використання граничних переходів Γ(1+n/2)Γ((1+n)/2)n2,(1x2)nexp(nx2), отримуються вирази що також походять з erf апроксимації. Де взявши похідну по вимірності fn від виразу для позиції f(n,nt2)=(n(1+2nt2)2πexp(nt2)(nt21))=0, та замінивши значення з експонентами із первинного, можна отримати зручне диференційне рівняння для допоміжного виразу nt2. Вибір допоміжного виразу пов'язаний із поведінкою другої похідної, позиція її максимума для апроксимації nt21=0 або (n2)t21=0 для точного виразу. Розв'язком межевого диференційного рівняння є логарифмічна функція.

Поведінка біля тотожності з півкулею t=0, оцінка у формі апроксимації:

- Для верхньої межі t>0 можна користуватись обмеженням ступеня вираза розрахунку об'єма через суми, що тотожньо розкладу вираза з гіпергеометричною функцією:

Vdome(r,h,n)Vball(r,n)hr(12Γ(n2+1)πΓ(n+12)(tn16t3))n1hr(12n2π(tn6t3))

- Для нижньої t>0 лінійної оцінки достатньо відкинути доданок із кубічною залежністю, тоді її крайні межі: 1π2n<hr<1+π2n або |t|<π2n,n1.

Апроксимація, поведінку біля тотожності з півкулею за великих вимірностей добре описує наближення через стандартний нормальний розподіл[5], яке можна виразити через функцію помилок erf,:

Апроксимація нормованого об'єму, erf, та її оцінка exp(t2n/2)t2πn.

Vdome(r,t,n)=12Vball(r,n)(1erf(tn/2)), erf(z)=2π0zexp(x2)dx,t=cos(θ)=rhr,n>10

з апроксимації походить проста верхня межа t>0:

VdomeVball=1sgn(t)2+exp(t2n/2)t2πn

що добре працює біля нуля |t|<ln(n)n, але, як і саме наближення, значно завищує відносні значення за віддалення.

Об'єднана оцінка, вираз для малих висот може бути значно розширений по діапазону до меж лінійної оцінки біля тотожності з півкулею, зміною знаменника подібно до межі отриманої з erf апроксимації:

VdomeVball=1sgn(t)2+(1t2)n/2t2πn, |t|>π2n,n>2

Обидві оцінки нескладно отримати з інтегральної форми через висоту, перша є наслідком спрощення до функції обмеженого росту множника (1+x)=2 в інтегральному виразі (1x2)=(1x)(1+x)2(1x) та відновленням форми після інтегрування, друга використовує визначення експоненти (1x2)n=(1nx2n)nexp(nx2). Але їх очевидне поєднання, через подібність у висновках, стосується областей де не працюють припущення що з ними пов'язані, тому доведення справедливості дії об'єднаної оцінки потребує додаткових кроків.

Асимптотична поведінка площі, для сферичної частини сегмента Scup подібна до об'єму Vdome, що можна побачити з рекурентного зв'язку. Тому заміною відношень об'ємів у виразах на відношення площ Vdome/VballScup/Ssphere отримуються необхідні вирази.

Відмінності є в оцінці біля тотожності з півкулею для частини що пов'язана з невеликими значеннями вимірності, де за переходу потрібно провести додаткову заміну nn2, за тримірного та двомірного випадків можна користуватись точними чи похідними з них виразами. Алгебраїчно це пов'язано зі зміною поведінки гіпергеометричної функції, яка за додатних значень другого аргументу має вигляд нескінченного ряду замість кінцевих сум.

Площа

У просторі вимірності n, поверхня може бути представлена поєднанням n мірного сферичного сегмента та поверхні основи, площа останньої Sbase дорівнює об'єму відповідної кулі:

Sdome=Scup+Sbase, де Sbase(r,h,n)=Vball(a,n1), Vball(a,k)=akπk/2Γ(k/2+1), a=h(2rh)=rsin(θ)=r1t2,t=cos(θ).

За обмеження аргумента виразу по висоті, отримання значень в необхідному діапазоні досягається заміною hs2rh та Scup(r,h,n)Ssphere(r,n)Scup(r,hs,n).

Інколи можна зустріти позначення площі поверхні з прив'язкою до лінійної вимірності, але це є неоднозначним, оскільки не вироджена до площини поверхня має додаткову параметричну інформацію та в просторі параметрів має вимірність геометричного простору якому вона належить. Перехід до залежності від лінійних параметрів є граничним, і тільки у спеціальних випадках простим.

З використанням спеціальних функцій

гіпергеометричної функції 2F1, чи регуляризованої неповної бета-функції[4] Ix(a,b):


Scup(r,h,n)=12Ssphere(r,n)Isin2θ(n12,12)=Vball(a,n1)2F1(12,n12;n+12;(ar)2)


де rh, sin2θ=(ar)2=h(2rh)r2, Ssphere(r,k)=krVball(r,k)=2πk/2rk1Γ(k/2) - площа гіперсфери.

або без обмежень по висоті:

Scup(r,h,n)=Ssphere(r,n)(12rhrΓ(n2)πΓ(n12)2F1(12,3n2;32;(rhr)2))

Через кінцеву ітерацію

За k=nmod(n,2)2 кроків, F(m,t)=t1(1x2)m2dx:

Алгебраїчна форма
обмеження кроки вираз параметри
n2 кінцевий Scup(r,h,n)=(n1)Vball(r,n1)F(n3,t) p=mod(n1,2), t=cos(θ)=rhr
i1 k1 F(i+2,t)=1i+3((i+2)F(i,t)t(1t2)i+22) F(p,t)={1tp=0arccos(t)p=1
i0 kp Vball(r,i+2)=2πr2i+2Vball(r,i) Vball(r,p)={1p=02rp=1

Ітераційні частини є такими самими як і для об'єму, відмінності: частково у кількості кроків, мінімальній вимірності та кінцевому виразі.

Коли додатково необхідно розраховувати площу основи, більш оптимальним є розрахунок одиничної кулі, з використанням її об'єму як множника з радіусами сферичного сегмента та основи Vball(x,m)=Vball(1,m)xm. В будь якому разі не варто очікувати високу продуктивність на великих вимірностях, де зведення до ступеню у бібліотечних функціях реалізовано змінним бітовим зсувом а факторіали обраховуються за формулою Стірлінга, кінцеве ж значення просто узгоджується налаштованим градієнтним спуском чи методом Ньютона. Перевага простоти у "дешевизні останнього значущого біта" який можна отримувати прямо з виразу без додаткових узгоджень, але "ціною" звуження області значень вхідних параметрів.

Через кінцеві суми

У замкненій формі. Має такі ж зауваження щодо порядку сумування як і для об'єму.

Scup(r,h,n)={2πkr2k+1(arccos(t)k!2k1t2(2k1)!!i=0k1(1)iαit2i+1)p=0,k12πk+1r2(k+1)(2k(2k+1)!!1k!i=0k(ki)(1)it2i+12i+1)p=1,k0

n3,k=np22,p=mod(n,2),t=cos(θ)=rhrαi>0,, (ki)={k(k1)(k2)(ki+1)i!i11i=0 - біноміальний коефіцієнт.

  • Для парних вимірностей значення коефіцієнтів αi можна отримати ітерацією αi+1=(ki+1)(2i+2)αi2i+3, α0=1(2k1)!!(2k)!!. Для оберненого сумування αi1=(ki)(2i+1)αi2i, перший коефіцієнт має простий вигляд αk1=12k. Зв'язок слідує із умови рівності похідних розрахункового та інтегрального виразів для площі (1t2)k=(2k1)!!(2k)!!tS(t,k)+(1t2)dS(t,k)dt, S(t,k)=i=0k1(1)iαit2i+1. Через складну поведінку оператора суми для від'ємних граничних значень індекса, в загальному виразі обмежена мінімальна вимірність, але за потреби її можна доповнити явно n=2,p=0,k=0,S(t,k)=0,Scup(r,t,n)=2rarccos(t), наприклад через додаткову вхідну перевірку. Мінімальна складність через отримання коефіцієнтів ітерацією - подвійна.
  • Для не парних вимірностей, з оберненим порядком сумування, її аргумент інваріантний відносно заміни граничних значень індекса i=0k=i=k0.

Рекурентний зв'язок

- в явному вигляді, походить з ітераційного розрахунку:

Scup(r,t,i+2)=2πr2i(Scup(r,t,i)t(1t2)i12Vball(r,i1))=2πr2i(Scup(r,t,i)tSbase(r,t,i))

де i2,t=cos(θ)=rhr, Sbase(r,t,n)=Vball(a,n1), a=r1t2.

Для площі основи вираз походить з рекурентного зв'язку об'ємів кулі:

Sbase(r,t,i+2)=2πr2i+1Sbase(r,t,i)

Інтегральна форма

  • інтегрування через основу, циліндричне

Scup(r,h,n)=(n1)Vball(1,n1)0axn2(1(xr)2)12dx

де rh, a=h(2rh)=r1t2, t=cos(θ)=rhr=(1(ar)2)12.

  • інтегрування через висоту

Scup(r,h,n)=(n1)Vball(r,n1)t1(1x2)n32dx

де t=cos(θ)=rhr.

  • інтегрування через кут

Scup(r,h,n)=(n1)Vball(r,n1)0θsinn2ϕdϕ

де cos(θ)=rhr.

Зв'язок із об'ємом

походить з порівняння ітераційних форм:

Vdome=rn(Scupt(1t2)n12Vball(r,n1))=rn(ScuptSbase)=rn(Sdome(1+t)Sbase)

де n2, Scup=SdomeSbase, Sbase=Vball(a,n1), a=h(2rh)=r1t2, t=cos(θ)=rhr

за використання рекурентного зв'язку:

Scup(r,h,n+2)=2πrVdome(r,h,n)

де n0,Vdome(r,t,0)=arccos(t)π, рекурентність підтверджує те що двомірність для поверхні є межевою, тобто менші вимірності пов'язані із другою групою аксіом, алгеброю, де об'єкт лише частково пов'язаний із геометричним. Доповнення значення для нуль вимірності розглянуто в розділі розрахунку об'єму через суми.

Асимптотична поведінка сферичної частини сегмента подібна об'єму, так як для сфери та кулі справедливий такий самий рекурентний зв'язок Ssphere(r,n+2)=2πrVball(r,n), за великих значень вимірності n1 їх відношення тотожні Scup/SsphereVdome/Vball, а для малих значень достатньо зробити заміну nn2 в об'ємній частині.

Зв'язок із іншими фігурами

Об'єм гіперкульового сектора можна отримати через площу сферичного сегмента[4]:

Vsector=Scup(r,h,n)Ssphere(r,n)Vball(r,n), Vball(r,n)=rnπn/2Γ(n/2+1), Ssphere(r,n)=2πn/2rn1Γ(n/2)=nrVball(r,n)

Так як об'єм сектора складається з об'ємів конічної частини та сегмента, Vdome+Vcone=Vsector, Vcone=rhnSbase, r>h, можна більш цілісно записати вираз для зв'язку об'єму гіперкульового сегмента та площі його сферичної частини :

Vdome+Vcone=Vdome+rhnSbase=rnScup=Vsector

Узгодженість математичного терміну округлення

Якщо представляти відрізок послідовністю одномірних куль, то округлення веде себе подібно послідовності нескінченно вимірних, через асимптотичну поведінку нормованого об'єму сегмента, зв'язок із похибкою пов'язаний із апроксимацію через функцію помилок. Тобто термін округлення математично узгоджений з поведінкою нормованого об'єму гіперкульового сегмента у гіперпросторі.

Розвиток термінології як і наукових напрямків пов'язаний з академічністю, де доказовість це оператор що застосовується до множини аксіом, вибір останніх пов'язаний із аксіоматичним методом, що за спрощеного тлумачення стосовно термінології, вимагає найбільшого узгодження серед наукових та других напрямків видів діяльності.

Див. також

Примітки

Шаблон:ReflistШаблон:Math-stub

  1. Загальну форму доречно використовувати у підготовці виразів для конкретних вимірностей, аналітичних процедур пов'язаних з оптимізацією коли вираз є проміжним, чи у випадках коли неможливо передбачити вимірність переходу. Оптимізовані для конкретних вимірностей, значень параметрів, вирази мають не меншу точність але гарантовано вищу швидкість розрахунку.
  2. Визначення об'єму через граничні поверхні інколи приводить до асоціації сегмента із поверхнею сферичної частини, але таке позначення менш узгоджене. Якщо вважати відрізок площею із нульовим об'ємом, січна площина має другу вимірність, або потребує зв'язку із більшою вимірністю, та тому є радше аналогією додатності порівняно з об'ємним розглядом. В одновимірному випадку довжину, об'єм, проміжка обмежує однаково вимірні поверхні - точки. Відповідно у гіперпросторі сегмент віднесений до властивостей поверхні є залежним від об'ємного об'єкта терміном.
  3. Через зв'язок терміну із фізикою, цей параметр відповідає габаритному розміру - один з розмірів описаного гіперпрямокутника з мінімальною площею, за мінімальної потенціальної енергії тіла. Математично умова має складний вигляд, тому оптимізацію обирається проста просторова варіація позиції, очевидно що за спрощення втрачається загальність, тому зустрічається неоднозначність у визначеннях. У геометрії переходи у ізотропному просторі однаково зважені, тому висотами також інколи називають частини прямих що мають подібну побудову. У транспортній сфері висоту пов'язують із оптимізацією доступу та зберігання.
  4. 4,0 4,1 4,2 Шаблон:Статья
  5. Шаблон:Cite journal