بي-سبلاين

في التحليل العددي ، تُعدّ دالة B-spline (اختصارًا لـ base spline) نوعًا من دوال spline المصممة بحيث يكون لها حد أدنى من الدعم (التداخل) لدرجة معينة ، ونعومة ، ومجموعة من نقاط الانقطاع (العُقد التي تقسم نطاقها ) ، مما يجعلها لبنة أساسية لجميع دوال spline من تلك الدرجة. تُعرَّف دالة B-spline بأنها متعددة حدود مجزأة من الدرجة .ن{\displaystyle n}، بمعنى درجة منن-1{\displaystyle n-1}تُبنى هذه الدالة من مقاطع تلتقي عند هذه العقد، حيث يعتمد استمرار الدالة ومشتقاتها على عدد مرات تكرار كل عقدة (تعددها). يمكن التعبير عن أي دالة تكعيبية من درجة معينة بشكل فريد كتركيبة خطية من دوال تكعيبية من نفس الدرجة عند نفس العقد، [ 1 ] وهي خاصية تجعلها متعددة الاستخدامات في النمذجة الرياضية. يستخدم نوع فرعي خاص، هو الدوال التكعيبية الأساسية، عقدًا متساوية البعد .

يعود مفهوم B-splines إلى القرن التاسع عشر، عندما استكشف نيكولاي لوباتشيفسكي أفكارًا مماثلة في جامعة كازان في روسيا، [ 2 ] على الرغم من أن مصطلح "B-spline" قد صاغه إسحاق جاكوب شوينبيرج [ 3 ] في عام 1967، مما يعكس دورها كدوال أساسية . [ 4 ]

تُستخدم منحنيات B-spline على نطاق واسع في مجالات مثل التصميم بمساعدة الحاسوب (CAD) ورسومات الحاسوب ، حيث تقوم بتشكيل المنحنيات والأسطح من خلال مجموعة من نقاط التحكم ، وكذلك في تحليل البيانات لمهام مثل تركيب المنحنيات والتفاضل العددي للبيانات التجريبية.

منحنى سبلاين مرسوم كمجموع مرجح لمنحنيات بي-سبلاين مع نقاط تحكم/مضلع تحكم، ومنحنيات مكونة مميزة

تعريف

دالة B-spline تربيعية أساسية ذات متجه عقدة (0، 0، 0، 1، 2، 3، 3، 3) ونقاط تحكم (0، 0، 1، 0، 0)، ومشتقتها الأولى
دالة B-spline مكعبة أساسية ذات متجه عقدة (−2، −2، −2، −2، −1، 0، 1، 2، 2، 2، 2) ونقاط تحكم (0، 0، 0، 6، 0، 0، 0)، ومشتقتها الأولى
دالة بي-سبلاين رباعية أساسية ذات متجه عقدة (0، 0، 0، 0، 0، 1، 2، 3، 4، 5، 5، 5، 5، 5) ونقاط تحكم (0، 0، 0، 0، 1، 0، 0، 0، 0)، ومشتقاتها الأولى والثانية

منحنى بي-سبلاين من الرتبةص+1{\displaystyle p+1}هي مجموعة من الدوال متعددة الحدود المجزأةبأنا،ص(ت){\displaystyle B_{i,p}(t)}درجة علميةص{\displaystyle p}في متغيرت{\displaystyle t}قيمت{\displaystyle t}تُعرف نقاط التقاء أجزاء متعددة الحدود بالعقد، ويُرمز لها بـت0،ت1،ت2،...،تم{\displaystyle t_{0},t_{1},t_{2},\ldots ,t_{m}}وتم ترتيبها بترتيب تصاعدي.

بالنسبة لتسلسل معين من العقد، يوجد، حتى عامل قياس، منحنى فريدبأنا،ص(ت){\displaystyle B_{i,p}(t)}مُرضٍ

بأنا،ص(ت)={غير صفريلو تأنات<تأنا+ص+1،0لو ت<تأناتأنا+ص+1ت.{\displaystyle B_{i,p}(t)={\begin{cases}{\text{غير صفري}}&{\text{إذا كان }}t_{i}\leq t<t_{i+p+1},\\0&{\text{إذا كان }}t<t_{i}\lor t_{i+p+1}\leq t.\end{cases}}}

إذا أضفنا القيد الإضافي الذي

أنا=0م-ص-1بأنا،ص(ت)=1{\displaystyle \sum _{i=0}^{mp-1}B_{i,p}(t)=1}

للجميعت{\displaystyle t}بين العقدتص{\displaystyle t_{p}}وتم-ص{\displaystyle t_{mp}}ثم عامل القياس لـبأنا،ص(ت){\displaystyle B_{i,p}(t)}تُصبح ثابتة. العُقد الموجودة بينها (ولا تشملها)تص{\displaystyle t_{p}}وتم-ص{\displaystyle t_{mp}}وتسمى هذه العقد بالعقد الداخلية.

يمكن إنشاء دوال B-spline باستخدام صيغة كوكس-دي بور التكرارية. نبدأ بدالات B-spline من الدرجةص=0{\displaystyle p=0}أي كثيرات الحدود الثابتة القطعية.

بأنا،0(ت):={1لو تأنات<تأنا+1،0خلاف ذلك.{\displaystyle B_{i,0}(t):={\begin{cases}1&{\text{إذا كان }}t_{i}\leq t<t_{i+1},\\0&{\text{فيما عدا ذلك}}.\end{cases}}}

كلما ارتفع(ص+1){\displaystyle (p+1)}يتم تعريف منحنيات B-spline من الدرجة n عن طريق التكرار.

بأنا،ص(ت):=ت-تأناتأنا+ص-تأنابأنا،ص-1(ت)+تأنا+ص+1-تتأنا+ص+1-تأنا+1بأنا+1،ص-1(ت).{\displaystyle B_{i,p}(t):={\dfrac {t-t_{i}}{t_{i+p}-t_{i}}}B_{i,p-1}(t)+{\dfrac {t_{i+p+1}-t}{t_{i+p+1}-t_{i+1}}}B_{i+1,p-1}(t).}

ملكيات

دالة B-spline هي عبارة عن مجموعة من الأشرطة المرنة التي يتم التحكم بها بواسطة عدد من النقاط تُسمى نقاط التحكم ، مما يُنتج منحنيات سلسة. تُستخدم هذه الدوال لإنشاء وإدارة الأشكال والأسطح المعقدة باستخدام عدد من النقاط. وتُستخدم دوال B-spline ودوال بيزير على نطاق واسع في أساليب تحسين الأشكال . [ 5 ]

منحنى بي-سبلاين من الرتبةن{\displaystyle n}هي دالة متعددة الحدود مجزأة من الدرجةن-1{\displaystyle n-1}في متغيرx{\displaystyle x}يتم تعريفها علىن{\displaystyle n}المواقعتأنا{\displaystyle t_{i}}وتسمى هذه النقاط بالعقد أو نقاط الانقطاع، والتي يجب أن تكون بترتيب غير تنازلي.تأناتأنا+1{\displaystyle t_{i}\leq t_{i+1}}لا تُساهم دالة B-spline إلا في النطاق بين العقدة الأولى والأخيرة، وتكون قيمتها صفرًا في أي مكان آخر. إذا كانت كل عقدة مفصولة بنفس المسافةح{\displaystyle h}(أينح=تأنا+1-تأنا{\displaystyle h=t_{i+1}-t_{i}}) من سابقتها، يُطلق على متجه العقدة و B-splines المقابلة اسم "موحد" (انظر B-spline الأساسي أدناه).

لكل فترة عقدة محدودة حيث تكون القيمة غير صفرية، فإن دالة B-spline هي متعددة حدود من الدرجةن-1{\displaystyle n-1}دالة B-spline هي دالة متصلة عند نقاط العقد. [ ملاحظة 1 ] عندما تكون جميع نقاط العقد التابعة لدالة B-spline متميزة، فإن مشتقاتها تكون متصلة أيضًا حتى مشتقة الدرجة t.ن-2{\displaystyle n-2}إذا كانت العقد متطابقة عند قيمة معينة لـx{\displaystyle x}يقلّ معامل استمرارية رتبة المشتقة بمقدار 1 لكل عقدة متطابقة إضافية. قد تشترك دوال B-spline في مجموعة فرعية من عقدها، لكنّ دالتين B-spline معرفتين على نفس العقد تمامًا تكونان متطابقتين. بعبارة أخرى، تُعرَّف دالة B-spline بشكل فريد من خلال عقدها.

يُفرّق المرء بين العُقد الداخلية ونقاط النهاية. تغطي العُقد الداخليةx{\displaystyle x}-المجال الذي يهتم به المرء. بما أن منحنى B-spline واحد يمتد بالفعل علىن{\displaystyle n}وبالتالي، فإن العقد الداخلية تحتاج إلى تمديدها معن-1{\displaystyle n-1}تُحدد نقاط النهاية على كل جانب لتوفير دعم كامل لأول وآخر منحنى B-spline، مما يؤثر على فترات العقد الداخلية. لا تُؤخذ قيم نقاط النهاية في الاعتبار، وعادةً ما يتم تكرار أول أو آخر عقدة داخلية.

تكمن فائدة دوال B-spline في حقيقة أن أي دالة spline من الرتبةن{\displaystyle n}يمكن التعبير عن مجموعة معينة من العقد كتركيبة خطية من دوال B-spline:

Sن،ت(x)=أناαأنابأنا،ن(x).{\displaystyle S_{n,\mathbf {t} }(x)=\sum _{i}\alpha _{i}B_{i,n}(x).}

تؤدي دوال B-spline دور الدوال الأساسية لفضاء دوال spline ، ومن هنا جاءت التسمية. وتنتج هذه الخاصية عن حقيقة أن جميع الأجزاء لها نفس خصائص الاستمرارية، ضمن نطاق دعمها الفردي، عند العقد. [ 6 ]

يمكن اشتقاق تعبيرات الأجزاء متعددة الحدود عن طريق صيغة كوكس دي بور التكرارية [ 7 ].

بأنا،0(x):={1لو تأناx<تأنا+1،0خلاف ذلك.{\displaystyle B_{i,0}(x):={\begin{cases}1&{\text{إذا كان }}t_{i}\leq x<t_{i+1},\\0&{\text{فيما عدا ذلك}}.\end{cases}}}
بأنا،ك(x):=x-تأناتأنا+ك-تأنابأنا،ك-1(x)+تأنا+ك+1-xتأنا+ك+1-تأنا+1بأنا+1،ك-1(x).{\displaystyle B_{i,k}(x):={\frac {x-t_{i}}{t_{i+k}-t_{i}}}B_{i,k-1}(x)+{\frac {t_{i+k+1}-x}{t_{i+k+1}-t_{i+1}}}B_{i+1,k-1}(x).}

إنه،بأنا،0(x){\displaystyle B_{i,0}(x)}قيمة ثابتة جزئياً، إما واحد أو صفر، تشير إلى أي عقدة تقع x (صفر إذا تكررت العقدة i ). تتكون معادلة التكرار من جزأين:

x-تأناتأنا+ك-تأنا{\displaystyle {\frac {x-t_{i}}{t_{i+k}-t_{i}}}}

تتناقص القيم من الصفر إلى الواحد عندما تتغير قيمة x منتأنا{\displaystyle t_{i}}لتأنا+ك{\displaystyle t_{i+k}}، و

تأنا+ك+1-xتأنا+ك+1-تأنا+1{\displaystyle {\frac {t_{i+k+1}-x}{t_{i+k+1}-t_{i+1}}}}

تتناقص القيم من واحد إلى صفر عندما تتغير قيمة x منتأنا+1{\displaystyle t_{i+1}}لتأنا+ك+1{\displaystyle t_{i+k+1}}تكون قيم B المقابلة صفرًا خارج تلك النطاقات المحددة. على سبيل المثال،بأنا،1(x){\displaystyle B_{i,1}(x)}هي دالة مثلثية تساوي صفرًا من الأسفلx=تأنا{\displaystyle x=t_{i}}، منحدرات تؤدي إلى واحد عندx=تأنا+1{\displaystyle x=t_{i+1}}والعودة إلى الصفر عند وما بعدهx=تأنا+2{\displaystyle x=t_{i+2}}ومع ذلك، نظرًا لأن دوال الأساس B-spline لها دعم محلي ، يتم حساب B-splines عادةً بواسطة خوارزميات لا تحتاج إلى تقييم دوال الأساس حيث تكون قيمتها صفرًا، مثل خوارزمية دي بور .

تؤدي هذه العلاقة مباشرة إلى خوارزمية BSPLV المكتوبة بلغة FORTRAN ، والتي تولد قيم B-splines من الرتبة n عند x . [ 8 ] يوضح المخطط التالي كيف أن كل جزء من الرتبة n هو عبارة عن تركيبة خطية من أجزاء B-splines من الرتبة n  1 إلى يساره.

000بأنا-2،2بأنا-1،1بأنا،0بأنا-1،2بأنا،10بأنا،200{\displaystyle {\begin{matrix}&&0\\&0&\\0&&B_{i-2,2}\\&B_{i-1,1}&\\B_{i,0}&&B_{i-1,2}\\&B_{i,1}&\\0&&B_{i,2}\\&0&\\&&0\end{matrix}}}

تطبيق صيغة التكرار مع العقد عند(0،1،2،3){\displaystyle (0,1,2,3)}يعطي المكونات القطعية لـ B-spline المنتظم من الرتبة 3

ب0=x2/2ب1=(-2x2+6x-3)/2ب2=(3-x)2/2{\displaystyle {\begin{aligned}b_{0}&=x^{2}/2\\b_{1}&=(-2x^{2}+6x-3)/2\\b_{2}&=(3-x)^{2}/2\end{aligned}}}

أين

ب0،2(x)={ب0لو 0x<1،ب1لو 1x<2،ب2لو 2x<3،0خلاف ذلك.{\displaystyle B_{0,2}(x)={\begin{cases}b_{0}&{\text{if }}0\leq x<1,\\b_{1}&{\text{if }}1\leq x<2,\\b_{2}&{\text{if }}2\leq x<3,\\0&{\text{otherwise}}.\end{cases}}}

تظهر هذه الأجزاء في الرسم التخطيطي. كما هو موضح أدناه، يتم توضيح خاصية استمرارية دالة التجزئة التربيعية ومشتقتها الأولى عند العقد الداخلية.

في x=1: ب0=ب1=0.5، دب0دx=دب1دx=1.في x=2: ب1=ب2=0.5، دب1دx=دب2دx=-1.{\displaystyle {\begin{aligned}&{\text{At }}x=1\colon \ b_{0}=b_{1}=0.5,\ {\frac {db_{0}}{dx}}={\frac {db_{1}}{dx}}=1.\\[6pt]&{\text{At }}x=2\colon \ b_{1}=b_{2}=0.5,\ {\frac {db_{1}}{dx}}={\frac {db_{2}}{dx}}=-1.\end{aligned}}}

المشتقة الثانية لدالة B-spline من الدرجة 2 تكون غير متصلة عند العقد:

د2ب0دx2=1، د2ب1دx2=-2، د2ب2دx2=1.{\displaystyle {\frac {d^{2}b_{0}}{dx^{2}}}=1,\ {\frac {d^{2}b_{1}}{dx^{2}}}=-2,\ {\frac {d^{2}b_{2}}{dx^{2}}}=1.}

تم اقتراح نسخ أسرع من خوارزمية دي بور، لكنها تعاني من استقرار أقل نسبياً. [ 9 ] [ 10 ]

منحنى بي-سبلاين الأساسي

تتميز دالة B-spline الأساسية بمسافة ثابتة h بين العقد. وتُعدّ دوال B-spline الأساسية من رتبة n معينة نسخًا مُزاحة لبعضها البعض. ويمكن الحصول عليها من التعريف الأبسط. [ 11 ]

بأنا،ن،ت(x)=x-تأناحن[0،...،ن](-تأنا)+ن-1.{\displaystyle B_{i,n,t}(x)={\frac {x-t_{i}}{h}}n[0,\dots ,n](\cdot -t_{i})_{+}^{n-1}.}

يُستخدم رمز "العنصر النائب" للإشارة إلى أن الفرق المقسم رقم n للدالة(ت-x)+ن-1{\displaystyle (t-x)_{+}^{n-1}}يتم تحديد المتغيرين t و x عن طريق تثبيت x والنظر في(ت-x)+ن-1{\displaystyle (t-x)_{+}^{n-1}}كدالة لـ t فقط.

تحتوي دالة B-spline الأساسية على عقد متباعدة بشكل منتظم، وبالتالي فإن الاستيفاء بين العقد يساوي الالتفاف مع نواة التنعيم.

مثال، إذا أردنا إدخال ثلاث قيم بين عقد B-spline (ب{\displaystyle {\textbf {b}}}يمكننا كتابة الإشارة على النحو التالي:

x=[ب1،0،0،ب2،0،0،ب3،0،0،...،بن،0،0].{\displaystyle \mathbf {x} =[\mathbf {b} _{1},0,0,\mathbf {b} _{2},0,0,\mathbf {b} _{3},0,0,\dots ,\mathbf {b} _{n},0,0].}

التفاف الإشارةx{\displaystyle \mathbf {x} }باستخدام دالة مستطيلةح=[1/3،1/3،1/3]{\displaystyle \mathbf {h} =[1/3,1/3,1/3]}يُعطي قيم B-spline المُستكمَلة من الدرجة الأولى. أما استيفاء B-spline من الدرجة الثانية فهو عبارة عن عملية التفاف مع دالة مستطيلة مرتين.y=x*ح*ح{\displaystyle \mathbf {y} =\mathbf {x} *\mathbf {h} *\mathbf {h} }; من خلال الترشيح التكراري باستخدام دالة مستطيلة، يتم الحصول على استيفاء من الدرجة الأعلى.

يمكن إجراء استيفاء B-spline سريع على نطاق عينة منتظم باستخدام ترشيح المتوسط ​​التكراري. بدلاً من ذلك، تساوي دالة المستطيل دالة sinc في مجال فورييه . لذلك، فإن استيفاء التكعيبية يساوي ضرب الإشارة في مجال فورييه في sinc 4 .

انظر توزيع إيروين-هول# حالات خاصة للتعبيرات الجبرية للخطوط الأساسية B-splines من الدرجة 1-4.

منحنى P

يشير مصطلح P-spline إلى "B-spline المعاقب". ويشير إلى استخدام تمثيل B-spline حيث يتم تحديد المعاملات جزئيًا بواسطة البيانات المراد ملاءمتها ، وجزئيًا بواسطة دالة جزاء إضافية تهدف إلى فرض السلاسة لتجنب التجاوز في الملاءمة . [ 12 ]

يمكن لتقريبات P-spline ثنائية الأبعاد ومتعددة الأبعاد للبيانات أن تستخدم ضرب تقسيم الوجوه للمصفوفات لتقليل عمليات الحساب. [ 13 ]

التعبيرات المشتقة

إن مشتقة B-spline من الدرجة k هي ببساطة دالة لـ B-splines من الدرجة k  1: [ 14 ]

دبأنا،ك(x)دx=ك(بأنا،ك-1(x)تأنا+ك-تأنا-بأنا+1،ك-1(x)تأنا+ك+1-تأنا+1).{\displaystyle {\frac {dB_{i,k}(x)}{dx}}=k\left({\frac {B_{i,k-1}(x)}{t_{i+k}-t_{i}}}-{\frac {B_{i+1,k-1}(x)}{t_{i+k+1}-t_{i+1}}}\right).}

وهذا يعني أن

ددxأناαأنابأنا،ك=أنا=ر-ك+2s-1كαأنا-αأنا-1تأنا+ك-تأنابأنا،ك-1على[تر،تs]،{\displaystyle {\frac {d}{dx}}\sum _{i}\alpha _{i}B_{i,k}=\sum _{i=r-k+2}^{s-1}k{\frac {\alpha _{i}-\alpha _{i-1}}{t_{i+k}-t_{i}}}B_{i,k-1}\quad {\text{on}}\quad [t_{r},t_{s}],}

مما يدل على وجود علاقة بسيطة بين مشتقة دالة السبلاين ودوال السبلاين من الدرجة الأقل بواحد.

لحظات دوال بي-سبلاين أحادية المتغير

يمكن استخدام دوال B-spline أحادية المتغير، أي دوال B-spline حيث تقع مواضع العقد في بُعد واحد، لتمثيل دوال كثافة الاحتمال أحادية البعد.ص(x){\displaystyle p(x)}ومن الأمثلة على ذلك المجموع المرجح لـأنا{\displaystyle i}دوال أساس B-spline من الرتبةن{\displaystyle n}، والتي يتم تطبيع كل منها إلى الوحدة (أي لا يتم تقييمها مباشرة باستخدام خوارزمية دي بور القياسية)

ص(x)=أناجأنابأنا،ن،معيار(x){\displaystyle p(x)=\sum _{i}c_{i}\cdot B_{i,n,{\textbf {norm}}}(x)}

ومع قيد ثابت التطبيعأناجأنا=1{\displaystyle \sum _{i}c_{i}=1}اللحظة الخام رقم kμك{\displaystyle \mu _{k}}من منحنى B-spline المعياريبأنا،ن،معيار{\displaystyle B_{i,n,{\textbf {norm}}}}يمكن كتابتها على شكل متوسط ​​كارلسون ديريشليRك{\displaystyle R_{k}}[ 15 ] والتي بدورها يمكن حلها بدقة عبر التكامل الكفافي والمجموع التكراري [ 16 ] كما

μك=Rك(م؛ت)=-xكبأنا،ن،معيار(x|ت1...تج)دx=Γ(ك+1)Γ(م)Γ(م+ك)دك(م،ت){\displaystyle \mu _{k}=R_{k}(\mathbf {m} ;\mathbf {t} )=\int _{-\infty }^{\infty }x^{k}\cdot B_{i,n,{\textbf {norm}}}(x\mid t_{1}\dots t_{j})\,dx={\frac {\Gamma (k+1)\Gamma (m)}{\Gamma (m+k)}}\cdot D_{k}(\mathbf {m} ,\mathbf {t} )}

مع

دك=1كu=1ك[(أنا=1جمأناتأناu)دك-u]{\displaystyle D_{k}={\frac {1}{k}}\sum \limits _{u=1}^{k}\left[\left(\sum \limits _{i=1}^{j}m_{i}\cdot {t_{i}}^{u}\right)D_{k-u}\right]}

ود0=1{\displaystyle D_{0}=1}. هنا،ت{\displaystyle \mathbf {t} }يمثل متجهًا معج{\displaystyle j}مواضع العقد وم{\displaystyle \mathbf {m} }متجه ذو تعدد عقدي متناسب. وبالتالي، يمكن حساب أي عزم لدالة كثافة الاحتمال.ص(x){\displaystyle p(x)}يتم تمثيلها بمجموع دوال أساس B-spline بدقة، دون اللجوء إلى التقنيات العددية.

العلاقة بمنحنى بيزير متعدد الأجزاء/المركب

منحنى بيزير هو أيضًا منحنى متعدد الحدود يمكن تعريفه باستخدام علاقة تكرارية من منحنيات ذات درجة أقل من نفس الفئة، ويتم ترميزه بدلالة نقاط التحكم، ولكن الفرق الرئيسي هو أن جميع الحدود في العلاقة التكرارية لجزء من منحنى بيزير لها نفس مجال التعريف (عادةً[0،1]{\displaystyle [0,1]})، بينما تختلف دعامات الحدين في تكرار B-spline (الفترات الفرعية الخارجية غير مشتركة). هذا يعني أن منحنى بيزير من الدرجةن{\displaystyle n}مقدم منمن{\displaystyle m\gg n}تتكون نقاط التحكم من حواليم/ن{\displaystyle m/n}تتميز معظم أجزاء منحنى بيزير باستقلاليتها، بينما ينتقل منحنى بيزير ذو المعاملات نفسها بسلاسة من فترة فرعية إلى أخرى. وللحصول على منحنى مماثل، يلزم فرض شرط سلاسة على الانتقالات بين الأجزاء، مما ينتج عنه نوع من منحنى بيزير (حيث تُحدد العديد من نقاط التحكم بشرط السلاسة).

منحنى بيزير متعدد القطع/المركب هو سلسلة من منحنيات بيزير متصلة ببعضها البعض على الأقلج0{\displaystyle C^{0}}الاستمرارية (تتطابق النقطة الأخيرة لمنحنى ما مع نقطة بداية المنحنى التالي). وبحسب التطبيق، توجد متطلبات إضافية للنعومة (مثلج1{\displaystyle C^{1}}أوج2{\displaystyle C^{2}}يمكن إضافة (الاستمرارية). [ 17 ]ج1{\displaystyle C^{1}}المنحنيات المتصلة لها مماسات متطابقة عند نقطة الانقطاع (حيث يلتقي المنحنيان).ج2{\displaystyle C^{2}}المنحنيات المتصلة لها نفس الانحناء عند نقطة الانقطاع. [ 18 ]

ملاءمة المنحنيات

في عملية مطابقة المنحنيات ، عادةً ما تُطابق مجموعة من نقاط البيانات مع منحنى مُعرَّف بدالة رياضية. على سبيل المثال، تستخدم الأنواع الشائعة من مطابقة المنحنيات متعددة الحدود أو مجموعة من الدوال الأسية . عندما لا يتوفر أساس نظري لاختيار دالة المطابقة، يمكن مطابقة المنحنى بدالة تكعيبية مُكوَّنة من مجموع دوال تكعيبية من النوع B، باستخدام طريقة المربعات الصغرى . [ 19 ] [ ملاحظة 2 ] وبالتالي، فإن دالة الهدف لتقليل المربعات الصغرى، بالنسبة لدالة تكعيبية من الدرجة k ، هي:

يو=الجميع x{دبليو(x)[y(x)-أناαأنابأنا،ك،ت(x)]}2،{\displaystyle U=\sum _{{\text{all}}~x}\left\{W(x)\left[y(x)-\sum _{i}\alpha _{i}B_{i,k,t}(x)\right]\right\}^{2},}

حيث W ( x ) هو وزن، و y ( x ) هي قيمة البيانات عند x . المعاملاتαأنا{\displaystyle \alpha _{i}}هذه هي المعايير المطلوب تحديدها. يمكن تثبيت قيم العقد أو التعامل معها كمعايير.

تكمن الصعوبة الرئيسية في تطبيق هذه العملية في تحديد عدد العقد المستخدمة ومواقعها. يقترح دي بور استراتيجيات متنوعة لمعالجة هذه المشكلة. على سبيل المثال، يتم تقليل المسافة بين العقد بما يتناسب مع انحناء البيانات (المشتقة الثانية). وقد نُشرت بعض التطبيقات، منها استخدام دوال B-spline لمطابقة منحنيات لورنتزية وغوسية مفردة . تم حساب دوال السبلاين المثلى من الدرجات 3 إلى 7، استنادًا إلى ترتيبات متناظرة من 5 و6 و7 عقد، وطُبقت الطريقة لتنعيم وتفاضل المنحنيات الطيفية. [ 20 ] في دراسة مقارنة، حقق الإصدار ثنائي الأبعاد من ترشيح سافيتزكي-غولاي وطريقة السبلاين نتائج أفضل من المتوسط ​​المتحرك أو ترشيح تشيبيشيف . [ 21 ]

التصميم بمساعدة الحاسوب والرسومات الحاسوبية

في تطبيقات التصميم بمساعدة الحاسوب ورسومات الحاسوب ، يُشار أحيانًا إلى منحنى السبلاين على النحو التالي:ج(ت){\displaystyle C(t)}، منحنى بارامتري لبعض المعلمات الحقيقيةت{\displaystyle t}في هذه الحالة المنحنىج(ت){\displaystyle C(t)}يمكن التعامل معها كدالتين أو ثلاث دوال إحداثية منفصلة(x(ت)،y(ت)){\displaystyle (x(t),y(t))}، أو(x(ت)،y(ت)،z(ت)){\displaystyle (x(t),y(t),z(t))}دوال الإحداثياتx(ت){\displaystyle x(t)}،y(ت){\displaystyle y(t)}وz(ت){\displaystyle z(t)}كل منها عبارة عن دوال انحناء، مع مجموعة مشتركة من قيم العقدت1،ت2،...،تن{\displaystyle t_{1},t_{2},\ldots ,t_{n}}.

بما أن دوال B-spline تشكل دوال أساسية، فإنه يمكن التعبير عن كل دالة من دوال الإحداثيات كمجموع خطي لدوال B-spline، لذلك لدينا

X(ت)=أناxأنابأنا،ن(ت)،Y(ت)=أناyأنابأنا،ن(ت)،Z(ت)=أناzأنابأنا،ن(ت).{\displaystyle {\begin{aligned}X(t)&=\sum _{i}x_{i}B_{i,n}(t),\\Y(t)&=\sum _{i}y_{i}B_{i,n}(t),\\Z(t)&=\sum _{i}z_{i}B_{i,n}(t).\end{aligned}}}

الأوزانxأنا{\displaystyle x_{i}}،yأنا{\displaystyle y_{i}}وzأنا{\displaystyle z_{i}}يمكن دمجها لتشكيل نقاطPأنا=(xأنا،yأنا،zأنا){\displaystyle P_{i}=(x_{i},y_{i},z_{i})}في الفضاء ثلاثي الأبعاد. هذه النقاطPأنا{\displaystyle P_{i}}تُعرف هذه النقاط عادةً باسم نقاط التحكم.

بالعمل عكسيًا، يُحدد تسلسل نقاط التحكم وقيم العقد وترتيب منحنى B-spline منحنىً بارامتريًا. يتميز هذا التمثيل للمنحنى بواسطة نقاط التحكم بعدة خصائص مفيدة:

  1. نقاط التحكمPأنا{\displaystyle P_{i}}قم بتعريف منحنى. إذا تم تحويل جميع نقاط التحكم معًا بطريقة ما، مثل نقلها أو تدويرها أو تغيير حجمها أو تحريكها بواسطة أي تحويل خطي ، فإن المنحنى المقابل يتم تحويله بنفس الطريقة.
  2. لأن دوال B-splines لا تساوي الصفر إلا لعدد محدود من فترات العقد، فإذا تم تحريك نقطة تحكم واحدة، فإن التغيير المقابل في المنحنى البارامتري يكون فوق نطاق المعلمات لعدد صغير من فترات العقد.
  3. لأنأنابأنا،ن(x)=1{\displaystyle \sum _{i}B_{i,n}(x)=1}وفي جميع الأوقاتبأنا،ن(x)0{\displaystyle B_{i,n}(x)\geq 0}إذاً، يبقى المنحنى داخل المربع المحيط بنقاط التحكم. كما أن المنحنى، بمعنى ما، يتبع نقاط التحكم بشكل عام.

من السمات الأقل استحساناً أن المنحنى البارامتري لا يقوم باستيفاء نقاط التحكم. عادةً لا يمر المنحنى عبر نقاط التحكم.

منحنيات بي-سبلاين المكعبة

منحنى بي-سبلاين مكعبج(ت){\displaystyle \mathbf {C} (t)}بمعامل معياريت[0،1]{\displaystyle t\in [0,1]}يتم تعريفها بأربع عقد (أي نقاط تحكم )ب0{\displaystyle {\textbf {b}}_{0}}،ب1{\displaystyle {\textbf {b}}_{1}}،ب2{\displaystyle {\textbf {b}}_{2}}، وب3{\displaystyle {\textbf {b}}_{3}}يشكل هذا متعدد حدود من الدرجة الثالثة يمكن كتابته على النحو التالي

ج(ت)=16[ت3ت2ت1][-13-313-630-30301410][ب0ب1ب2ب3]{\displaystyle \mathbf {C} (t)={\frac {1}{6}}\;{\begin{bmatrix}t^{3}&t^{2}&t&1\end{bmatrix}}{\begin{bmatrix}-1&3&-3&1\\3&-6&3&0\\-3&0&3&0\\1&4&1&0\end{bmatrix}}{\begin{bmatrix}\mathbf {b} _{0}\\\mathbf {b} _{1}\\\mathbf {b} _{2}\\\mathbf {b} _{3}\end{bmatrix}}}.

يتوافق هذا مع كثيرات حدود بي-سبلاين

ب0(ت)=16(-ت3+3ت2-3ت+1)ب1(ت)=16(3ت3-6ت2+4)ب2(ت)=16(-3ت3+3ت2+3ت+1)ب3(ت)=16ت3{\displaystyle {\begin{aligned}B_{0}(t)&={\frac {1}{6}}(-t^{3}+3t^{2}-3t+1)\\B_{1}(t)&={\frac {1}{6}}(3t^{3}-6t^{2}+4)\\B_{2}(t)&={\frac {1}{6}}(-3t^{3}+3t^{2}+3t+1)\\B_{3}(t)&={\frac {1}{6}}t^{3}\end{aligned}}}

ويمكن تقييم المنحنى على النحو التاليج(ت)=أنا=03بأنا(ت)بأنا{\displaystyle \mathbf {C} (t)=\sum _{i=0}^{3}B_{i}(t)\,\mathbf {b} _{i}}بتوسيع هذا، يمكننا كتابة الصيغة الكاملة لكثير الحدود كما يلي

ج(ت)=16((-ب0+3ب1-3ب2+ب3)ت3+(3ب0-6ب1+3ب2)ت2+(-3ب0+3ب2)ت+(ب0+4ب1+ب2)){\displaystyle \mathbf {C} (t)={\frac {1}{6}}{\biggl (}(-\mathbf {b} _{0}+3\mathbf {b} _{1}-3\mathbf {b} _{2}+\mathbf {b} _{3})t^{3}+(3\mathbf {b} _{0}-6\mathbf {b} _{1}+3\mathbf {b} _{2})t^{2}+(-3\mathbf {b} _{0}+3\mathbf {b} _{2})t+(\mathbf {b} _{0}+4\mathbf {b} _{1}+\mathbf {b} _{2}){\biggr )}}.

بما أن هذه دالة كثيرة الحدود من الدرجة الثالثة ، فيمكننا أيضًا كتابتها كمنحنى بيزير من الدرجة الثالثة مع نقاط تحكم. P0{\displaystyle {\textbf {P}}_{0}}،P1{\displaystyle {\textbf {P}}_{1}}،P2{\displaystyle {\textbf {P}}_{2}}، وP3{\displaystyle {\textbf {P}}_{3}}بحيث

P0=16(ب0+4ب1+ب2)،P1=13(2ب1+ب2)،P2=13(ب1+2ب2)،P3=16(ب1+4ب2+ب3).{\displaystyle {\begin{aligned}\mathbf {P} _{0}&={\frac {1}{6}}(\mathbf {b} _{0}+4\mathbf {b} _{1}+\mathbf {b} _{2}),\\\mathbf {P} _{1}&={\frac {1}{3}}(2\mathbf {b} _{1}+\mathbf {b} _{2}),\\\mathbf {P} _{2}&={\frac {1}{3}}(\mathbf {b} _{1}+2\mathbf {b} _{2}),\\\mathbf {P} _{3}&={\frac {1}{6}}(\mathbf {b} _{1}+4\mathbf {b} _{2}+\mathbf {b} _{3}).\end{aligned}}}

يتم تشكيل منحنى B-spline المكعب المجزأ بواسطة مجموعة من العقد، وتحدد كل أربع عقد متتالية جزءًا مكعبًا من المنحنى بالصيغة المذكورة أعلاه.

شرائح B المثبتة

إذا كان الأول والأخيرن+1{\displaystyle n+1}العقد (ن{\displaystyle n}(درجة) متجه العقدة لها نفس القيمة، وتتداخل نقاط البداية والنهاية للمنحنى مع نقاط التحكم، مما ينتج عنه B-splines مثبتة.

مثال على متجهات العقد لشرائح B-spline المكعبة المثبتة معم{\displaystyle M}نقاط التحكم هي كالتالي:

م{\displaystyle M}متجه العقدةنطاقت{\displaystyle t}
4{\displaystyle 4}[00001111]{\displaystyle {\begin{bmatrix}0&0&0&0&1&1&1&1\end{bmatrix}}}[0,1]
5{\displaystyle 5}[000012222]{\displaystyle {\begin{bmatrix}0&0&0&0&1&2&2&2&2\end{bmatrix}}}[0,2]
6{\displaystyle 6}[0000123333]{\displaystyle {\begin{bmatrix}0&0&0&0&1&2&3&3&3&3\end{bmatrix}}}[0,3]
7{\displaystyle 7}[00001234444]{\displaystyle {\begin{bmatrix}0&0&0&0&1&2&3&4&4&4&4\end{bmatrix}}}[0,4]

يمكن تقسيم منحنى B-spline المثبت إلى أجزاء B-spline مجزأة باستخدام معلمة معياريةت[0،1]{\displaystyle t\in [0,1]}، محددة بأربع نقاط تحكمبأنا{\displaystyle {\textbf {b}}_{i}}،بأنا+1{\displaystyle {\textbf {b}}_{i+1}}،بأنا+2{\displaystyle {\textbf {b}}_{i+2}}، وبأنا+3{\displaystyle {\textbf {b}}_{i+3}}ويمكن كتابتها على النحو التالي:

جأنا(ت)=[ت3ت2ت1]Rأنا[بأنابأنا+1بأنا+2بأنا+3]{\displaystyle \mathbf {C} _{i}(t)={\begin{bmatrix}t^{3}&t^{2}&t&1\end{bmatrix}}\mathbf {R} _{i}{\begin{bmatrix}\mathbf {b} _{i}\\\mathbf {b} _{i+1}\\\mathbf {b} _{i+2}\\\mathbf {b} _{i+3}\end{bmatrix}}}،

أينRأنا{\displaystyle \mathbf {R} _{i}}هي مصفوفة معاملات .

يمكننا أيضًا كتابتها كمنحنى بيزير مكعب مع نقاط تحكم Pأنا0{\displaystyle \mathbf {P_{i}} _{0}}،Pأنا1{\displaystyle \mathbf {P_{i}} _{1}}،Pأنا2{\displaystyle \mathbf {P_{i}} _{2}}، وPأنا3{\displaystyle \mathbf {P_{i}} _{3}}بحيث

[Pأنا0Pأنا1Pأنا2Pأنا3]=Rبz-1Rأنا[بأنابأنا+1بأنا+2بأنا+3]،wحهـرهـRبz=[-13-313-630-33001000]{\displaystyle {\begin{bmatrix}\mathbf {P_{i}} _{0}\\\mathbf {P_{i}} _{1}\\\mathbf {P_{i}} _{2}\\\mathbf {P_{i}} _{3}\end{bmatrix}}=\mathbf {R_{bz}} ^{-1}\mathbf {R} _{i}{\begin{bmatrix}\mathbf {b} _{i}\\\mathbf {b} _{i+1}\\\mathbf {b} _{i+2}\\\mathbf {b} _{i+3}\end{bmatrix}},\qquad where\quad \mathbf {R_{bz}} ={\begin{bmatrix}-1&3&-3&1\\3&-6&3&0\\-3&3&0&0\\1&0&0&0\end{bmatrix}}}.

متجهات العقد، مصفوفة المعاملاتRأنا{\displaystyle \mathbf {R} _{i}}ومصفوفة التحويلRبz-1Rأنا{\displaystyle \mathbf {R_{bz}} ^{-1}\mathbf {R} _{i}}بالنسبة لكل جزء من أجزاء B-spline، تكون القيم كما يلي: [ 22 ] [ 23 ]

م{\displaystyle M}أنا{\displaystyle i}متجه العقدةمصفوفة المعاملاتRأنا{\displaystyle \mathbf {R} _{i}}مصفوفة التحويلRبz-1Rأنا{\displaystyle \mathbf {R_{bz}} ^{-1}\mathbf {R} _{i}}
4{\displaystyle 4}0{\displaystyle 0}[00001111]{\displaystyle {\begin{bmatrix}0&0&0&0&1&1&1&1\end{bmatrix}}}[-13-313-630-33001000]{\displaystyle {\begin{bmatrix}-1&3&-3&1\\3&-6&3&0\\-3&3&0&0\\1&0&0&0\end{bmatrix}}}[1000010000100001]{\displaystyle {\begin{bmatrix}1&0&0&0\\0&1&0&0\\0&0&1&0\\0&0&0&1\end{bmatrix}}}
5{\displaystyle 5}0{\displaystyle 0}[00001222]{\displaystyle {\begin{bmatrix}0&0&0&0&1&2&2&2\end{bmatrix}}}[-174-1143-92320-33001000]{\displaystyle {\begin{bmatrix}-1&{\frac {7}{4}}&-1&{\frac {1}{4}}\\3&-{\frac {9}{2}}&{\frac {3}{2}}&0\\-3&3&0&0\\1&0&0&0\end{bmatrix}}}[100001000121200141214]{\displaystyle {\begin{bmatrix}1&0&0&0\\0&1&0&0\\0&{\frac {1}{2}}&{\frac {1}{2}}&0\\0&{\frac {1}{4}}&{\frac {1}{2}}&{\frac {1}{4}}\end{bmatrix}}}
1{\displaystyle 1}[-1-1-101111]{\displaystyle {\begin{bmatrix}-1&-1&-1&0&1&1&1&1\end{bmatrix}}}[-141-74134-32340-3403401412140]{\displaystyle {\begin{bmatrix}-{\frac {1}{4}}&1&-{\frac {7}{4}}&1\\{\frac {3}{4}}&-{\frac {3}{2}}&{\frac {3}{4}}&0\\-{\frac {3}{4}}&0&{\frac {3}{4}}&0\\{\frac {1}{4}}&{\frac {1}{2}}&{\frac {1}{4}}&0\end{bmatrix}}}[141214001212000100001]{\displaystyle {\begin{bmatrix}{\frac {1}{4}}&{\frac {1}{2}}&{\frac {1}{4}}&0\\0&{\frac {1}{2}}&{\frac {1}{2}}&0\\0&0&1&0\\0&0&0&1\end{bmatrix}}}
6{\displaystyle 6}0{\displaystyle 0}[00001233]{\displaystyle {\begin{bmatrix}0&0&0&0&1&2&3&3\end{bmatrix}}}[-174-1112163-92320-33001000]{\displaystyle {\begin{bmatrix}-1&{\frac {7}{4}}&-{\frac {11}{12}}&{\frac {1}{6}}\\3&-{\frac {9}{2}}&{\frac {3}{2}}&0\\-3&3&0&0\\1&0&0&0\end{bmatrix}}}[1000010001212001471216]{\displaystyle {\begin{bmatrix}1&0&0&0\\0&1&0&0\\0&{\frac {1}{2}}&{\frac {1}{2}}&0\\0&{\frac {1}{4}}&{\frac {7}{12}}&{\frac {1}{6}}\end{bmatrix}}}
1{\displaystyle 1}[-1-1-101222]{\displaystyle {\begin{bmatrix}-1&-1&-1&0&1&2&2&2\end{bmatrix}}}[-14712-7121434-54120-341412014712160]{\displaystyle {\begin{bmatrix}-{\frac {1}{4}}&{\frac {7}{12}}&-{\frac {7}{12}}&{\frac {1}{4}}\\{\frac {3}{4}}&-{\frac {5}{4}}&{\frac {1}{2}}&0\\-{\frac {3}{4}}&{\frac {1}{4}}&{\frac {1}{2}}&0\\{\frac {1}{4}}&{\frac {7}{12}}&{\frac {1}{6}}&0\end{bmatrix}}}[1471216002313001323001671214]{\displaystyle {\begin{bmatrix}{\frac {1}{4}}&{\frac {7}{12}}&{\frac {1}{6}}&0\\0&{\frac {2}{3}}&{\frac {1}{3}}&0\\0&{\frac {1}{3}}&{\frac {2}{3}}&0\\0&{\frac {1}{6}}&{\frac {7}{12}}&{\frac {1}{4}}\end{bmatrix}}}
2{\displaystyle 2}[-2-2-101111]{\displaystyle {\begin{bmatrix}-2&-2&-1&0&1&1&1&1\end{bmatrix}}}[-161112-74112-54340-12-1434016712140]{\displaystyle {\begin{bmatrix}-{\frac {1}{6}}&{\frac {11}{12}}&-{\frac {7}{4}}&1\\{\frac {1}{2}}&-{\frac {5}{4}}&{\frac {3}{4}}&0\\-{\frac {1}{2}}&-{\frac {1}{4}}&{\frac {3}{4}}&0\\{\frac {1}{6}}&{\frac {7}{12}}&{\frac {1}{4}}&0\end{bmatrix}}}[1671214001212000100001]{\displaystyle {\begin{bmatrix}{\frac {1}{6}}&{\frac {7}{12}}&{\frac {1}{4}}&0\\0&{\frac {1}{2}}&{\frac {1}{2}}&0\\0&0&1&0\\0&0&0&1\end{bmatrix}}}
7{\displaystyle \geq 7}0{\displaystyle 0}[00001234]{\displaystyle {\begin{bmatrix}0&0&0&0&1&2&3&4\end{bmatrix}}}[-174-1112163-92320-33001000]{\displaystyle {\begin{bmatrix}-1&{\frac {7}{4}}&-{\frac {11}{12}}&{\frac {1}{6}}\\3&-{\frac {9}{2}}&{\frac {3}{2}}&0\\-3&3&0&0\\1&0&0&0\end{bmatrix}}}112[120000120006600372]{\displaystyle {\frac {1}{12}}{\begin{bmatrix}12&0&0&0\\0&12&0&0\\0&6&6&0\\0&3&7&2\end{bmatrix}}}
1{\displaystyle 1}[-1-1-101234]{\displaystyle {\begin{bmatrix}-1&-1&-1&0&1&2&3&4\end{bmatrix}}}[-14712-121634-54120-341412014712160]{\displaystyle {\begin{bmatrix}-{\frac {1}{4}}&{\frac {7}{12}}&-{\frac {1}{2}}&{\frac {1}{6}}\\{\frac {3}{4}}&-{\frac {5}{4}}&{\frac {1}{2}}&0\\-{\frac {3}{4}}&{\frac {1}{4}}&{\frac {1}{2}}&0\\{\frac {1}{4}}&{\frac {7}{12}}&{\frac {1}{6}}&0\end{bmatrix}}}112[3720084004800282]{\displaystyle {\frac {1}{12}}{\begin{bmatrix}3&7&2&0\\0&8&4&0\\0&4&8&0\\0&2&8&2\end{bmatrix}}}
2أنام-6{\displaystyle 2\leq i\leq M-6}[-2-2-101233]{\displaystyle {\begin{bmatrix}-2&-2&-1&0&1&2&3&3\end{bmatrix}}}[ ملاحظة 3 ]16[-13-313-630-30301410]{\displaystyle {\frac {1}{6}}{\begin{bmatrix}-1&3&-3&1\\3&-6&3&0\\-3&0&3&0\\1&4&1&0\end{bmatrix}}}16[1410042002400141]{\displaystyle {\frac {1}{6}}{\begin{bmatrix}1&4&1&0\\0&4&2&0\\0&2&4&0\\0&1&4&1\end{bmatrix}}}
م-5{\displaystyle M-5}[-3-2-101222]{\displaystyle {\begin{bmatrix}-3&-2&-1&0&1&2&2&2\end{bmatrix}}}[-1612-7121412-1120-1201201623160]{\displaystyle {\begin{bmatrix}-{\frac {1}{6}}&{\frac {1}{2}}&-{\frac {7}{12}}&{\frac {1}{4}}\\{\frac {1}{2}}&-1&{\frac {1}{2}}&0\\-{\frac {1}{2}}&0&{\frac {1}{2}}&0\\{\frac {1}{6}}&{\frac {2}{3}}&{\frac {1}{6}}&0\end{bmatrix}}}112[2820084004800273]{\displaystyle {\frac {1}{12}}{\begin{bmatrix}2&8&2&0\\0&8&4&0\\0&4&8&0\\0&2&7&3\end{bmatrix}}}
م-4{\displaystyle M-4}[-3-2-101111]{\displaystyle {\begin{bmatrix}-3&-2&-1&0&1&1&1&1\end{bmatrix}}}[-161112-74112-54340-12-1434016712140]{\displaystyle {\begin{bmatrix}-{\frac {1}{6}}&{\frac {11}{12}}&-{\frac {7}{4}}&1\\{\frac {1}{2}}&-{\frac {5}{4}}&{\frac {3}{4}}&0\\-{\frac {1}{2}}&-{\frac {1}{4}}&{\frac {3}{4}}&0\\{\frac {1}{6}}&{\frac {7}{12}}&{\frac {1}{4}}&0\end{bmatrix}}}112[273006600012000012]{\displaystyle {\frac {1}{12}}{\begin{bmatrix}2&7&3&0\\0&6&6&0\\0&0&12&0\\0&0&0&12\end{bmatrix}}}

نيربس

منحنى NURBS منحنى متعدد الحدود معرف في إحداثيات متجانسة (أزرق) وإسقاطه على المستوى منحنى كسري (أحمر)

في التصميم بمساعدة الحاسوب ، والتصنيع بمساعدة الحاسوب ، ورسومات الحاسوب ، يُعدّ منحنى بي-سبلاين غير المنتظم العقلاني (NURBS) امتدادًا قويًا لمنحنيات بي-سبلاين. وNURBS هي في الأساس منحنيات بي-سبلاين في إحداثيات متجانسة . وكما هو الحال مع منحنيات بي-سبلاين، تُعرَّف منحنيات بي-سبلاين برتبتها، ومتجه العقد، ومجموعة من نقاط التحكم، ولكن على عكس منحنيات بي-سبلاين البسيطة، فإن لكل نقطة تحكم وزنًا. عندما يكون الوزن مساويًا لـ 1، يكون منحنى NURBS هو ببساطة منحنى بي-سبلاين، وبذلك يُعمِّم NURBS كلاً من منحنيات بي-سبلاين ومنحنيات وأسطح بيزير ، والفرق الرئيسي هو وزن نقاط التحكم الذي يجعل منحنيات NURBS "عقلانية".

من خلال تقييم NURBS عند قيم مختلفة للمعلمات، يمكن تتبع المنحنى عبر الفضاء؛ وبالمثل، من خلال تقييم سطح NURBS عند قيم مختلفة للمعلمتين، يمكن تمثيل السطح في الفضاء الديكارتي.

كما هو الحال في منحنيات B-spline، تحدد نقاط التحكم في NURBS شكل المنحنى. تُحسب كل نقطة من نقاط المنحنى بأخذ مجموع مرجح لعدد من نقاط التحكم. يختلف وزن كل نقطة تبعًا للمعامل المتحكم. بالنسبة لمنحنى من الدرجة d ، يكون تأثير أي نقطة تحكم غير صفري فقط في d + 1 فترات (امتدادات العقد) من فضاء المعاملات . ضمن هذه الفترات، يتغير الوزن وفقًا لدالة متعددة الحدود (دوال الأساس) من الدرجة d . عند حدود الفترات، تتلاشى دوال الأساس بسلاسة إلى الصفر، وتُحدد درجة هذه السلاسة بدرجة متعددة الحدود.

متجه العقد هو سلسلة من قيم المعاملات التي تحدد مكان وكيفية تأثير نقاط التحكم على منحنى NURBS. عدد العقد يساوي دائمًا عدد نقاط التحكم مضافًا إليه درجة المنحنى زائد واحد. في كل مرة تدخل فيها قيمة المعامل نطاق عقد جديد، تُفعَّل نقطة تحكم جديدة، بينما تُهمل نقطة التحكم القديمة.

يأخذ منحنى NURBS الشكل التالي: [ 24 ]

ج(u)=أنا=1كشمالأنا،ن(u)wأناPأناأنا=1كشمالأنا،ن(u)wأنا{\displaystyle C(u)={\frac {\sum _{i=1}^{k}N_{i,n}(u)w_{i}P_{i}}{\sum _{i=1}^{k}N_{i,n}(u)w_{i}}}}

تُستخدم الرموز التالية: u هو المتغير المستقل (بدلاً من x )، و k هو عدد نقاط التحكم، و N هو دالة B-spline (تُستخدم بدلاً من Bوn هي درجة متعددة الحدود، و P هي نقطة تحكم، و w هو وزن. المقام هو عامل توحيد يساوي واحدًا إذا كانت جميع الأوزان تساوي واحدًا.

من المعتاد كتابة هذا على النحو التالي

ج(u)=أنا=1كRأنا،ن(u)Pأنا{\displaystyle C(u)=\sum _{i=1}^{k}R_{i,n}(u)P_{i}}

حيث الوظائف

Rأنا،ن(u)=شمالأنا،ن(u)wأناج=1كشمالج،ن(u)wج{\displaystyle R_{i,n}(u)={\frac {N_{i,n}(u)w_{i}}{\sum _{j=1}^{k}N_{j,n}(u)w_{j}}}}

تُعرف باسم دوال الأساس الكسرية.

يتم الحصول على سطح NURBS كحاصل ضرب موتر لمنحنيين NURBS، وبالتالي باستخدام معلمتين مستقلتين u و v (مع المؤشرات i و j على التوالي): [ 25 ]

S(u،v)=أنا=1كج=1Rأنا،ج(u،v)Pأنا،ج{\displaystyle S(u,v)=\sum _{i=1}^{k}\sum _{j=1}^{\ell }R_{i,j}(u,v)P_{i,j}}

مع

Rأنا،ج(u،v)=شمالأنا،ن(u)شمالج،م(v)wأنا،جص=1كq=1شمالص،ن(u)شمالq،م(v)wص،q{\displaystyle R_{i,j}(u,v)={\frac {N_{i,n}(u)N_{j,m}(v)w_{i,j}}{\sum _{p=1}^{k}\sum _{q=1}^{\ell }N_{p,n}(u)N_{q,m}(v)w_{p,q}}}}

كدوال أساسية نسبية.

انظر أيضاً

ملحوظات

  1. بالمعنى الدقيق للكلمة، تُعرَّف منحنيات B-spline عادةً بأنها متصلة من اليسار.
  2. يقدم دي بور إجراءات FORTRAN لتركيب البيانات التجريبية بطريقة المربعات الصغرى.
  3. متجه العقدة مكافئ لـ[-3-2-101234]{\displaystyle {\begin{bmatrix}-3&-2&-1&0&1&2&3&4\end{bmatrix}}}وذلك لأن العناصر الموجودة عند طرفي متجه العقدة لا تُستخدم في حساب الدوال الأساسية. ويمكن أن تكون قيماً عشوائية.

مراجع

  1. ^ هارتموت براوتش. فولفجانج بوهم؛ ماركو بالوزني (2002). تقنيات بيزييه و بي سبلاين . الرياضيات والتصور. برلين، هايدلبرغ: سبرينغر للعلوم والإعلام التجاري. ص.  63. دوى : 10.1007/978-3-662-04919-8 . رقم ISBN 978-3-540-43761-1. OCLC 851370272 . 
  2. فارين، جي إي (2002). المنحنيات والأسطح للتصميم بمساعدة الحاسوب: دليل عملي . مورغان كوفمان. ص 119.
  3. دي بور، ص 114.
  4. غاري د. نوت (2000)، استيفاء الدوال التكعيبية . سبرينغر. ص 151.
  5. طالبي توتي، ر.؛ شجاعي فرد، م.ح.؛ يارمحمديستري، صادق (2015). "تحسين تصميم شكل الخزان الأسطواني باستخدام منحنيات بي-سبلاين". الحوسبة والموائع . 109 : 100-112 . doi : 10.1016/j.compfluid.2014.12.004 .
  6. دي بور، ص 113.
  7. دي بور، ص 131.
  8. دي بور، ص 134.
  9. لي، إي تي واي (ديسمبر 1982). "روتين حسابي مبسط لـ B-Spline". الحوسبة . 29 (4): 365-371 . doi : 10.1007/BF02246763 . S2CID 2407104 . 
  10. لي، إي تي واي (1986). "تعليقات على بعض خوارزميات بي-سبلاين". الحوسبة . 36 (3): 229-238 . doi : 10.1007/BF02240069 . S2CID 7003455 . 
  11. دي بور، ص 322.
  12. إيلرز، بي إتش سي وماركس، بي دي (1996). التنعيم المرن باستخدام دوال بي-سبلاين والعقوبات (مع التعليقات والردود). العلوم الإحصائية 11(2): 89-121.
  13. إيلرز، بول إتش سي؛ ماركس، برايان دي. (2003). "المعايرة متعددة المتغيرات مع تفاعل درجة الحرارة باستخدام انحدار الإشارة المعاقب ثنائي الأبعاد". علم القياسات الكيميائية وأنظمة المختبرات الذكية . 66 (2): 159-174 . doi : 10.1016/S0169-7439(03)00029-7 .
  14. دي بور، ص 115.
  15. كارلسون، بي سي (1991). "دوال بي-سبلاين، والدوال فوق الهندسية، ومتوسطات ديريشليه" . مجلة نظرية التقريب . 67 (3): 311-325 . doi : 10.1016/0021-9045(91)90006-V .
  16. غلوسنكامب، ت. (2018). "المعالجة الاحتمالية لعدم اليقين الناتج عن الحجم المحدود لبيانات مونت كارلو الموزونة". EPJ Plus . 133 (6): 218. arXiv : 1712.01293 . Bibcode : 2018EPJP..133..218G . doi : 10.1140/epjp/i2018-12042-x . S2CID 125665629 . )
  17. يوجين ف. شيكين؛ ألكسندر إ. بليس (14 يوليو 1995). دليل استخدام الدوال المنحنية . مطبعة سي آر سي. الصفحات 96 وما بعدها. رقم ISBN  978-0-8493-9404-1.
  18. ويرنيك، جوزي (1993). "8" . دليل المخترع: برمجة الرسومات ثلاثية الأبعاد الموجهة للكائنات باستخدام أوبن إنفينتور، الإصدار 2 ( الطبعة الأولى). بوسطن، ماساتشوستس، الولايات المتحدة الأمريكية: شركة أديسون-ويسلي لونغمان للنشر. ISBN  978-0201624953.
  19. ^ دي بور، الفصل الرابع عشر، ص. 235.
  20. غانز، بيتر؛ جيل، ج. برنارد (1984). "تنعيم وتفاضل المنحنيات الطيفية باستخدام دوال التجزئة". علم الأطياف التطبيقي . 38 (3): 370-376 . Bibcode : 1984ApSpe..38..370G . doi : 10.1366/0003702844555511 . S2CID 96229316 . 
  21. فيسيك، ماريا؛ نيل، شارون ل.؛ وارنر، إيزايا م. (1986). "ترشيح بيانات التألق ثنائية الأبعاد في المجال الزمني" . مطيافية تطبيقية . 40 (4): 542-548 . رمز Bibcode : 1986ApSpe..40..542V . doi : 10.1366/0003702864508773 . S2CID 28705788. مؤرشف من الأصل في 23 يونيو 2017. 
  22. كوهين، إيلين؛ ريزنفيلد، ريتشارد ف. (1982). "التمثيلات المصفوفية العامة لمنحنيات بيزير و بي-سبلاين". الحواسيب في الصناعة . 3 : 13-15 .
  23. "مقدمة لمنحنيات بي-سبلاين لمهندسي التنفيذ" . كيتا .
  24. بيغل وتيلر، الفصل 4، القسم 2
  25. بيغل وتيلر، الفصل 4، القسم 4

المراجع

  • كارل دي بور (1978). دليل عملي للخطوط . سبرينغر-فيرلاغ. رقم ISBN 978-3-540-90356-7.
  • بيغل، ليس؛ تيلر، واين (1997). كتاب NURBS (الطبعة الثانية  ). سبرينغر. ISBN 978-3-540-61545-3.

للمزيد من القراءة

  • ريتشارد هـ. بارتلز؛ جون س. بيتي؛ ​​برايان أ. بارسكي (1987). مقدمة في استخدام الدوال التكعيبية في رسومات الحاسوب والنمذجة الهندسية . مورغان كوفمان. ISBN 978-1-55860-400-1.
  • جان غالييه (1999). المنحنيات والأسطح في النمذجة الهندسية: النظرية والخوارزميات . مورغان كوفمان. الفصل 6. منحنيات بي-سبلاين.هذا الكتاب غير متوفر في الأسواق ويمكن الحصول عليه مجاناً من المؤلف.
  • هارتموت براوتش؛ فولفجانج بوهم؛ ماركو بالوزني (2002). تقنيات بيزييه و بي سبلاين . سبرينغر العلوم والإعلام التجاري. رقم ISBN 978-3-540-43761-1.
  • ديفيد سالومون (2006). المنحنيات والأسطح في رسومات الحاسوب . سبرينغر. الفصل 7. تقريب بي-سبلاين. ISBN 978-0-387-28452-1.
  • هوفي، تشاد (2022). صياغة وتنفيذ هندسة بيزير وهندسة بي-سبلاين باستخدام لغة بايثون. SAND2022-7702C . (153 صفحة)