Spline (matematika) - Spline (mathematics)

1/3 va 2/3 dagi bitta tugunlar uch kubikli polinomlarning splinasini tashkil qiladi C2 uzluksizlik. Intervalning ikkala uchidagi uchta tugun egri chiziqning so'nggi nuqtalarni interpolatsiya qilishini ta'minlaydi

Yilda matematika, a spline maxsus funktsiya belgilangan qismli tomonidan polinomlar.In interpolatsiya qilish muammolar, spline interpolatsiyasi ko'pincha afzaldir polinom interpolatsiyasi chunki u past darajadagi polinomlardan foydalanganda ham, qochish paytida ham shunga o'xshash natijalarni beradi Runge fenomeni yuqori darajalar uchun.

In Kompyuter fanlari subfields kompyuter yordamida loyihalash va kompyuter grafikasi, spline atamasi ko'proq polinomni anglatadi (parametrli) egri chiziq. Splines ushbu pastki maydonlarda mashhur egri chiziqlardir, chunki ularning tuzilishining soddaligi, baholashning qulayligi va to'g'riligi va murakkab shakllarni taxminiy qilish qobiliyatlari egri chiziq va interfaol egri dizayni.

Spline atamasi egiluvchanlikdan kelib chiqadi spline kema quruvchilar tomonidan ishlatiladigan qurilmalar va chizmachilar silliq shakllar chizish uchun.

Kirish

"Spline" atamasi ma'lumotlarning interpolatsiyasini va / yoki tekislashni talab qiladigan dasturlarda ishlatiladigan keng funktsiyalar sinfiga murojaat qilish uchun ishlatiladi. Ma'lumotlar bir o'lchovli yoki ko'p o'lchovli bo'lishi mumkin. Interpolyatsiya uchun spline funktsiyalari odatda interpolatsiya cheklovlariga rioya qilingan holda pürüzlülüğün mos o'lchovlarini minimayzerlari sifatida belgilanadi (masalan, integral kvadrat egrilik). Silliqlashtiruvchi splinlarni interpolatsiya splini umumlashtirilishi deb hisoblash mumkin, bu erda funktsiyalar kuzatilgan ma'lumotlar va pürüzlülük o'lchovi bo'yicha o'rtacha kvadratik taxminiy xatolikning og'irlashtirilgan kombinatsiyasini minimallashtirish uchun belgilanadi. Pürüzlülük o'lchovining bir qator mazmunli ta'riflari uchun spline funktsiyalari cheklangan o'lchovli bo'lib, bu ularning hisoblash va tasvirlashda foydaliligining asosiy sababi hisoblanadi. Ushbu bo'limning qolgan qismida biz butunlay bir o'lchovli, polinomial splinelarga e'tibor qaratamiz va "spline" atamasini ushbu cheklangan ma'noda ishlatamiz.

Ta'rif

Biz o'z bahsimizni cheklash bilan boshlaymiz bir o'zgaruvchan polinom ishi. Bunday holda, spline a qismli polinom funktsiyasi.Ushbu funktsiya, uni chaqiring S, qiymatlarni intervaldan oladi [a,b] va ularni xaritalar , to'plami haqiqiy raqamlar,

Biz xohlaymiz S qismlarga bo'linib belgilanishi kerak. Buni amalga oshirish uchun [a,b] tomonidan qamrab olinishi kerak k er-xotin disjoint bilan buyurtma qilingan subintervallar ichki qismlar,

Ularning har biri bo'yicha k "qismlari" ning [a,b], biz polinomni aniqlamoqchimiz, uni chaqiring Pmen.

.

Ustida men[ning subintervallaria,b], S bilan belgilanadi Pmen,

Berilgan k + 1 ochkolar tmen deyiladi tugunlar. Vektor deyiladi a tugun vektori spline uchun. Agar tugunlar intervalda teng masofada taqsimlansa [a,b] biz spline deb aytamiz bir xil, aks holda biz shunday deymiz bir xil bo'lmagan.

Agar polinom bo'laklari bo'lsa Pmen ularning har biri eng yuqori darajaga ega n, keyin spline deb aytiladi daraja (yoki of.)buyurtma n + 1).

Agar mahallasida tmen, keyin spline beof deb aytiladi silliqlik (kamida) da tmen. Ya'ni, da tmen ikki qism Pi-1 va Pmen 0 buyrug'i (funktsiya qiymati) ning hosilasi orqali tartib hosilasi orqali yuqoridagi umumiy hosila qiymatlarini baham ko'ring rmen (boshqacha qilib aytganda, ikkita qo'shni polinom bo'lagi bilan bog'lanadi silliqlikni yo'qotish ko'pi bilan n - rmen).

Vektor shplenning silliqligi bor da tmen uchun deyiladi a silliqlik vektori spline uchun.

Tugun vektori berilgan , daraja nva silliqlik vektori uchun , barcha darajadagi splinelar to'plamini ko'rib chiqish mumkin tugunli vektorga ega va silliqlik vektori . Ikkita funktsiyani qo'shish (nuqtali qo'shish) va funktsiyalarning haqiqiy ko'paytmalarini olish bilan jihozlangan ushbu to'plam haqiqiy vektor makoniga aylanadi. Bu spline bo'sh joy odatda tomonidan belgilanadi .

Ko'p polinomlarni matematik o'rganishda ikkita tugun, deylik, nima bo'ladi degan savol tmen va tmen+1, birgalikda ko'chirish oson javobga ega. Polinom bo'lagiPmen(t) yo'qoladi va parchalarPmen−1(t) va Pmen+1(t) uchun uzluksiz yo'qotishlarning yig'indisi bilan qo'shilingtmen va tmen+1.Anavi,

qayerda

Bu tugun vektorini yanada kengroq tushunishga olib keladi va har qanday nuqtada uzluksizlikni yo'qotish natijasi deb hisoblanishi mumkinbir nechta tugun o'sha nuqtada joylashgan va spline turi uning darajasi bilan to'liq tavsiflanishi mumkin n va uning kengaytirilgan tugun vektori

qayerda tmen takrorlanadi jmen vaqt uchun .

A parametrik egri oralig'ida [a,b]

a spline egri chizig'i agar ikkalasi bo'lsa X va Y spline funktsiyalari, xuddi shu darajadagi uzaytirilgan tugun vektorlari bilan bir xil darajada.

Misollar

Faraz qilaylik [a,b] [0,3] va subintervallar [0,1], [1,2] va [2,3]. Areto 2 darajali polinom bo'laklari, va [0,1] va [1,2] qismlar qiymat va birinchi hosilaga qo'shilishi kerak deylik (at t= 1) bo'lsa, [1,2] va [2,3] qismlar oddiy qiymati bo'yicha qo'shiladi (at t = 2) .Bu spline turini belgilaydi S(t) buning uchun

shu turdagi a'zolar bo'lar edi, shuningdek

Ushbu turdagi a'zolar bo'lar edi. (Izoh: polinom bo'lagi 2 bo'lsat kvadratik emas, natijada baribir kvadratik spline deyiladi. Bu shpinning darajasi uning polinom qismlarining maksimal darajasidir.) Ushbu turdagi spline uchun kengaytirilgan tugun vektori (0, 1, 2, 2, 3) bo'ladi.

Eng oddiy spline 0 darajaga ega. U a deb ham nomlanadi qadam funktsiyasi.Quyidagi eng oddiy spline 1 darajaga ega. U shuningdek a deb ham nomlanadi chiziqli spline. Tekislikdagi yopiq chiziqli spline (ya'ni birinchi tugun va oxirgisi bir xil) shunchaki a ko'pburchak.

Umumiy spline - bu tabiiy kubik spline uzluksizligi bilan 3 daraja C2. "Tabiiy" so'zi shpin polinomlarining ikkinchi hosilalari interpolatsiya oralig'ining so'nggi nuqtalarida nolga tengligini anglatadi.

Bu splinni intervaldan tashqarida to'g'ri chiziq bo'lishga majbur qiladi, shu bilan birga uning silliqligini buzmaydi.

Tabiiy kubik splini hisoblash algoritmi

Kubik splinelar shaklga ega .
Berilgan koordinatalar to'plami to'plamini topishni xohlaymiz splinelar uchun

Ular quyidagilarni qondirishi kerak:

  • .

Keling, bitta kubik splini aniqlaymiz 5 karra sifatida qayerda va ilgari ko'rsatilgan shakldagi koeffitsientlarga mos keladi va ga teng

Tabiiy kubik splini hisoblash algoritmi:
Kirish: koordinatalar to'plami , bilan
Chiqish: tarkibiga o'rnatilgan spline-larni o'rnating n 5-gilzalar.

  1. Yangi qator yarating a hajmi n + 1 va uchun o'rnatilgan
  2. Yangi massivlar yarating b va d har bir o'lcham n.
  3. Yangi qator yarating h hajmi n va uchun o'rnatilgan
  4. Yangi qator yarating a hajmi n va uchun o'rnatilgan .
  5. Yangi massivlar yarating v, l, mva z har bir o'lcham .
  6. O'rnatish
  7. Uchun
    1. O'rnatish .
    2. O'rnatish .
    3. O'rnatish .
  8. O'rnatish
  9. Uchun
    1. O'rnatish
    2. O'rnatish
    3. O'rnatish
  10. Yangi Splines to'plamini yarating va uni output_set deb nomlang. Uni to'ldiring n splinelar S.
  11. Uchun
    1. O'rnatish Smen,a = amen
    2. O'rnatish Smen,b = bmen
    3. O'rnatish Smen,v = vmen
    4. O'rnatish Smen,d = dmen
    5. O'rnatish Smen,x = xmen
  12. Chiqish_set

Izohlar

Buning ma'nosi nimadan ko'proq ekanligini so'rashi mumkin n tugunli vektorda bir nechta tugun bor, chunki bu kabi davomiylikka olib keladi

bu juda ko'p sonli joyda. An'anaga ko'ra, har qanday bunday vaziyat ikkita qo'shni polinom bo'lagi orasidagi oddiy uzilishni ko'rsatadi. Bu degani, agar tugun bo'lsa tmen dan ko'proq ko'rinadi n + 1 marta kengaytirilgan tugun vektorida, uning barcha holatlari (n + 1) th spline xarakterini o'zgartirmasdan o'chirilishi mumkin, chunki barcha ko'paytmalar n + 1, n + 2, n + 3 va boshqalar bir xil ma'noga ega. Odatda har qanday spline turini belgilaydigan har qanday tugun vektori shu tarzda bekor qilingan deb taxmin qilinadi.

Klassik spline turi n raqamli tahlilda ishlatiladigan uzluksizlikka ega

bu shuni anglatadiki, har ikkala qo'shni polinom bo'lagi o'z qiymatiga ko'ra va birinchi bo'lib uchrashadi n - har bir tugunda 1 ta hosilalar. Modellashtirilgan matematik spline tekis spline kubik (n = 3), ikki marta doimiy ravishda farqlanadigan (C2), tabiiy spline, bu so'nggi nuqtalarda qo'shimcha shartlar qo'yilgan ushbu klassik turdagi spline a va b.

Grafika, masalan, rasm chizish dasturlarida juda ko'p ishlatiladigan yana bir spline turi Adobe Illustrator dan Adobe tizimlari, kubik bo'laklarga ega, lekin ko'pi bilan davomiylikka ega

Ushbu spline turi ham ishlatiladi PostScript shuningdek, ba'zi bir kompyuter tipografik shriftlarining ta'rifida.

Yuqori darajadagi grafikalar va animatsiya uchun mo'ljallangan ko'plab kompyuterli dizayn tizimlari, masalan, kengaytirilgan tugun vektorlaridan foydalanadi Mayya dan Taxalluslar.Kompyuter yordamida loyihalash tizimlarida ko'pincha a deb nomlanuvchi splinening kengaytirilgan tushunchasi ishlatiladi Bir xil bo'lmagan ratsional B-spline (NURBS).

Agar funktsiya yoki jismoniy ob'ektdan olingan namunalar mavjud bo'lsa, spline interpolatsiyasi bu ma'lumotlarga yaqinlashadigan spline yaratishga yondashuv.

A uchun umumiy ifoda C2 Kubik shpani interpolatsiya qilish

Uchun umumiy ifoda menth C2 bir nuqtada interpolyatsion kubik splini x tabiiy holat bilan formuladan foydalanib topish mumkin

qayerda

  • da ikkinchi lotin qiymatlari mentugun.
  • funktsiya qiymatlari mentugun.

Vakillar va ismlar

Berilgan interval uchun [a,b] va shu oraliqda berilgan kengaytirilgan tugun vektori, daraja splini n shakl vektor maydoni. Qisqacha aytganda, bu ma'lum bir turdagi har qanday ikkita spline qo'shilsa, ushbu turdagi spline hosil bo'ladi va berilgan turdagi splinani har qanday doimiyga ko'paytirilsa, ushbu turdagi spline hosil bo'ladi. The o'lchov kengaytirilgan tugun vektoridan ma'lum turdagi barcha splinelarni o'z ichiga olgan bo'shliqni hisoblash mumkin:

O'lchov daraja yig'indisiga va ko'paytmalarga teng

Agar spline turiga qo'shimcha chiziqli shartlar qo'yilgan bo'lsa, unda hosil bo'lgan spline pastki bo'shliqda yotadi. Masalan, barcha tabiiy kubik splinallarining maydoni barcha kublarning fazosining pastki fazosidir C2 splinelar.

Spline adabiyoti maxsus spline turlari uchun nomlar bilan to'ldirilgan bo'lib, ular quyidagi nomlar bilan bog'langan:

  • Splini ko'rsatish uchun qilingan tanlovlar, masalan:
  • Kengaytirilgan tugun vektorini shakllantirishda qilingan tanlovlar, masalan:
    • uchun bitta tugun yordamida Cn-1 uzluksizlik va bu tugunlarning oralig'ini bir tekisda [a,b] (bizga berish bir xil splinelar)
    • masofani cheklashsiz tugunlardan foydalanish (bizga berish bir xil bo'lmagan splinelar)
  • Spline-ga o'rnatilgan har qanday maxsus shartlar, masalan:
    • at nol sonli hosilalarni bajarish a va b (bizga berish tabiiy splinelar)
    • berilgan ma'lumotlarning splinada bo'lishini talab qilish (bizga berish splinlarni interpolatsiya qilish)

Ko'pincha yuqorida ko'rsatilgan ikkita yoki undan ko'p narsalarni qondiradigan spline turi uchun maxsus nom tanlangan. Masalan, Hermit spline - bu alohida polinom qismlarini har birini ifodalash uchun Hermit polinomlari yordamida ifodalangan spline. Ular ko'pincha bilan ishlatiladi n = 3; ya'ni Hermit kubiklari. Ushbu daraja sifatida ular qo'shimcha ravishda faqat tangens-uzluksiz tanlanishi mumkin (C1); bu barcha ichki tugunlarning ikki baravarligini anglatadi. Bunday splinelarni berilgan ma'lumotlarga moslashtirish uchun bir necha usullar ixtiro qilingan; ya'ni ularni interpolyatsion splinelarga aylantirish va buni har ikkala polinom bo'laklari uchrashadigan (bizni beradigan) mantiqiy teginish qiymatlarini baholash orqali amalga oshirish. Kardinal splinelar, Catmull-Rom splines va Kochanek-Bartels splines, ishlatilgan usulga qarab).

Namoyishlarning har biri uchun spline qiymatlari talabga binoan ishlab chiqarilishi uchun ba'zi baholash vositalarini topish kerak. Har bir alohida polinom bo'lagini ifodalaydigan vakillar uchun Pmen(t) daraja uchun ba'zi asoslar bo'yicha n polinomlar, bu kontseptual jihatdan sodda:

  • Argumentning berilgan qiymati uchun t, u yotadigan intervalni toping
  • Ushbu interval uchun tanlangan polinom asosini qidiring
  • Har bir asosli polinomning qiymatini toping t:
  • Ushbu intervalda spline beradigan asosli polinomlarning chiziqli birikmasining koeffitsientlarini qidiring v0, ..., vk-2
  • Spline qiymatini olish uchun asosli polinom qiymatlarining chiziqli kombinatsiyasini qo'shing t:

Biroq, baholash va yig'ish bosqichlari ko'pincha aqlli usullar bilan birlashtiriladi. Masalan, Bernshteyn polinomlari polinomlar uchun asos bo'lib, ularni chiziqli kombinatsiyalarda maxsus takrorlanish munosabatlari yordamida samarali baholash mumkin. Bu mohiyat De Kastelxau algoritmi, qaysi xususiyatlari Bézier egri chiziqlari va Bézier splines.

Spline-ni asosiy spline-larning chiziqli kombinatsiyasi sifatida belgilaydigan vakillik uchun yanada murakkabroq narsa kerak. The de Boor algoritmi baholashning samarali usuli hisoblanadi B-splinalar.

Tarix

Kompyuterlardan foydalanishdan oldin raqamli hisob-kitoblar qo'l bilan bajarilgan. Kabi qismlarga bo'linib belgilangan funktsiyalar bo'lsa ham belgi funktsiyasi yoki qadam funktsiyasi ishlatilgan, odatda polinomlarga ustunlik berilgan, chunki ular bilan ishlash osonroq edi. Kompyuterlarning paydo bo'lishi bilan splinelar muhim ahamiyat kasb etdi. Dastlab ular interpolatsiyadagi polinomlarni almashtirish sifatida, so'ngra kompyuter grafikalarida silliq va egiluvchan shakllarni qurish vositasi sifatida ishlatilgan.

Splinlarga birinchi matematik ma'lumot 1946 yilgacha bo'lgan qog'oz deb qabul qilingan Shoenberg, bu "spline" so'zining silliq, qismli polinomik yaqinlashish bilan bog'liq holda ishlatilishining birinchi joyidir. Biroq, g'oyalar o'zlarining samolyotlari va kemasozlik sanoatida ildiz otgan. (Bartels va boshq., 1987) ga kirish so'zida, Robin Forrest tasvirlaydi "balandlik ", davomida Britaniya aviatsiya sanoatida qo'llaniladigan texnika Ikkinchi jahon urushi yupqa yog'och chiziqlardan o'tib, samolyotlar uchun shablonlarni yaratishsplinelar ") loftning katta qavatida yotqizilgan nuqtalar orqali kema-korpus dizaynidan olingan usul. Bir necha yillar davomida kema dizayni amaliyotida kichkinagina dizayni uchun modellar ishlatilgan. Muvaffaqiyatli dizayn keyinchalik grafik qog'ozga tushirilgan va uchastkaning asosiy nuqtalari kattaroq kattalikdagi grafik qog'ozga qayta tushirildi.Yupqa yog'och chiziqlar asosiy nuqtalarni silliq egri chiziqlar bilan interpolyatsiyasini ta'minladi, chiziqlar diskret nuqtalarda (Forrest tomonidan "o'rdak" deb nomlangan) ushlab turilar edi. Shoenberg "itlar" yoki "kalamushlar" dan foydalangan) va shu nuqtalar orasida minimal kuchlanish energiyasi shakllari mavjud bo'lar edi. Forrestning so'zlariga ko'ra, ushbu jarayon uchun matematik model uchun mumkin bo'lgan turtki bu butun samolyot uchun juda muhim dizayn qismlarining yo'qolishi bo'lishi mumkin. loftni dushman bomba bilan urish kerakmi? Bu "konusli lofting" ni vujudga keltirdi, u erda konus kesimlari yordamida o'rdaklar orasidagi egri chiziqning o'rnini modellashtirish mumkin edi. Konik lofting o'rnini biz 1960 yillarning boshlarida spline deb ataydigan narsalar bilan almashtirdik. tomonidan nashr etilgan J. C. Fergyuson da Boeing va (birozdan keyin) tomonidan M.A.Sabin da Britaniya aviatsiya korporatsiyasi.

"Spline" so'zi dastlab an bo'lgan Sharqiy Angliya dialekt so'z.

Avtomobil korpuslarini modellashtirish uchun splinlardan foydalanish bir nechta mustaqil boshlanishlarga o'xshaydi. Kredit nomidan talab qilinadi de Casteljau da Citroen, Per Bezier da Renault va Birxof, Garabedian va de Boor da General Motors (qarang: Birkhoff va de Boor, 1965), barchasi 1960 yillarning boshlarida yoki 1950 yillarning oxirlarida sodir bo'lgan ishlar uchun. 1959 yilda de Kastelxauning hech bo'lmaganda bittasi nashr etilgan, ammo keng tarqalmagan. De Burning ishi General Motors 60-yillarning boshlarida bir qator hujjatlar, shu jumladan ba'zi bir fundamental ishlarni nashr etishga olib keldi B-splinalar.

Pratt & Whitney Aircraft aviakompaniyasida ham ishlar olib borildi, u erda mualliflarning ikkitasi (Ahlberg va boshq., 1967) - birinchi splinallarni uzunlikdagi muolajasi ishlagan va Devid Teylor model havzasi, Feodor Teygeymer tomonidan. Ish General Motors (Birkhoff, 1990) va (Young, 1997) da yaxshi yoritilgan. Devis (1997) ushbu materialning bir qismini qisqacha bayon qiladi.

Adabiyotlar

  • Fergyuson, Jeyms C, Ko'p o'zgaruvchan egri interpolatsiya, J. ACM, vol. 11, yo'q. 2, 221-228 betlar, 1964 yil aprel.
  • Ahlberg, Nilson va Uolsh, Spline nazariyasi va ularning qo'llanilishi, 1967.
  • Birkhoff, Suyuqlik dinamikasi, reaktorni hisoblash va sirtni namoyish qilish, Stiv Nash (tahr.), Ilmiy hisoblash tarixi, 1990.
  • Bartels, Bitti va Barski, Kompyuter grafikasi va geometrik modellashtirishda foydalanish uchun splinallarga kirish, 1987.
  • Birxof va de Bur, Parcha-parcha polinom interpolyatsiyasi va yaqinlashuvi, H. L. Garabedian (tahr.), Proc. General Motors simpoziumi 1964 yil, 164-190 betlar. Elsevier, Nyu-York va Amsterdam, 1965 yil.
  • Devis, B-splinelar va geometrik dizayn, SIAM News, jild 29, yo'q. 5, 1997 yil.
  • Epperson, Splinlar tarixi, NA Digest, jild 98, yo'q. 26, 1998 yil.
  • Stoer & Bulirsch, Raqamli tahlilga kirish. Springer-Verlag. p. 93-106. ISBN  0387904204
  • Schoenberg, Analitik funktsiyalar bo'yicha teng masofada joylashgan ma'lumotlarni yaqinlashtirish muammosiga qo'shgan hissasi, Kvart. Qo'llash. Matematik., jild 4, 45–99 va 112-141 betlar, 1946 y.
  • Yosh, Garrett Birxof va amaliy matematik, AMS xabarnomalari, jild 44, yo'q. 11, 1446–1449-betlar, 1997 y.
  • Chapra, Kanal, "Muhandislar uchun raqamli usullar" 5-nashr.

Tashqi havolalar

Nazariya

Excel funktsiyasi

Onlayn yordam dasturlari

Kompyuter kodi