פולינום ברנשטיין

פולינומי ברנשטיין לקירוב עקומה
בתחום האנליזה הנומרית, פולינום ברנשטיין, הקרוי על שם ממציאו, סרגיי ברנשטיין, הוא פולינום בתצורת ברנשטיין, כלומר הוא מהווה צירוף ליניארי של פולינומי הבסיס של ברנשטיין.

דרך יציבה מבחינה נומרית להעריך פולינומים בתצורת ברנשטיין היא בעזרת אלגוריתם דה-קסטלז'ו.

פולינומים בתצורת ברנשטיין היו בשימוש לראשונה על ידי ברנשטיין, בהוכחה קונסטרוקטיבית למשפט הקירוב של סטון–ויירשטראס. עם כניסת הגרפיקה הממוחשבת, הפכו פולינומי ברנשטיין (על הקטע החסום <math>[1,0]</math>) לשימושיים וחשובים ביצירת עקומות בזייה.

== הגדרה ==
n + 1 פולינומי הבסיס של ברנשטיין מדרגה n מוגדרים כך:
: <math>b_{\nu,n}(x) = {n \choose \nu} x^{\nu} \left( 1 - x \right)^{n - \nu}, \quad \nu = 0, \ldots, n.</math>
כאשר <math>{n \choose \nu}</math> הוא המקדם הבינומי.

פולינומי הבסיס של ברנשטיין מדרגה n יוצרים בסיס למרחב הווקטורי <math>\Pi_n</math> של פולינומים בדרגה של לכל היותר n.

צירוף ליניארי של פולינומי הבסיס של ברנשטיין,
:<math>B_n(x) = \sum_{\nu=0}^{n} \beta_{\nu} b_{\nu,n}(x)</math>
נקרא פולינום ברנשטיין או פולינום בתצורת ברנשטיין מדרגה n. מקדמי <math>\beta_\nu</math> נקראים מקדמי ברנשטיין או מקדמי בזייר.

== דוגמה ==
מספר פולינומי הבסיס הראשונים של ברנשטיין הם:
: <math>
\begin{align}
b_{0,0}(x) & = 1, \\
b_{0,1}(x) & = 1 - x, & b_{1,1}(x) & = x \\
b_{0,2}(x) & = (1 - x)^2, & b_{1,2}(x) & = 2x(1 - x), & b_{2,2}(x) & = x^2 \\
b_{0,3}(x) & = (1 - x)^3, & b_{1,3}(x) & = 3x(1 - x)^2, & b_{2,3}(x) & = 3x^2(1 - x), & b_{3,3}(x) & = x^3 \\
b_{0,4}(x) & = (1 - x)^4, & b_{1,4}(x) & = 4x(1 - x)^3, & b_{2,4}(x) & = 6x^2(1 - x)^2, & b_{3,4}(x) & = 4x^3(1 - x), & b_{4,4}(x) & = x^4
\end{align}
</math>

== תכונות ==
לפולינומי הבסיס של ברנשטיין התכונות הבאות:
* <math>b_{\nu, n}(x) = 0</math>, אם <math>\nu < 0</math> או <math>\nu > n</math>.
* <math>b_{\nu, n}(0) = \delta_{\nu, 0}</math> ו <math>b_{\nu, n}(1) = \delta_{\nu, n}</math> כאשר <math>\delta</math> היא הדלתא של קרונקר.
* <math>b_{\nu, n}(x)</math> יש שורש עם כפילות <math>\nu</math> בנקודה <math>x = 0</math> (הערה: אם <math>\nu = 0</math>, אין שורש ב-0).
* <math>b_{\nu, n}(x)</math> יש שורש עם כפילות <math>\left( n - \nu \right)</math> בנקודה <math>x = 1</math> (הערה: אם <math>\nu = n</math>, אין שורש ב-1).
* <math>b_{\nu, n}(x) \ge 0</math> עבור <math>x \in [0,\ 1]</math>.
* <math>b_{\nu, n}\left( 1 - x \right) = b_{n - \nu, n}(x)</math>.
* הנגזרת יכולה להיכתב כצירוף של שני פולינומים מדרגה נמוכה יותר:
*: <math>b'_{\nu, n}(x) = n \left( b_{\nu - 1, n - 1}(x) - b_{\nu, n - 1}(x) \right).</math>
* האינטגרל קבוע עבור <math>n</math> נתון:
*: <math>\int_{0}^{1}b_{\nu, n}(x)dx = \frac{1}{n+1} ; \forall \nu = 0,1 \dots n</math>
* אם <math>n \ne 0</math>, אז <math>b_{\nu, n}(x)</math> בעל מקסימום ייחודי מקומי באינטרוול <math>[0,\ 1]</math> ב <math>x = \frac{\nu}{n}</math>. המקסימום הנ"ל בעל הערך:
*: <math>\nu^\nu n^{-n} \left( n - \nu \right)^{n - \nu} {n \choose \nu}.</math>
* פולינומי הבסיס של ברנשטיין מדרגה <math>n</math> יוצרים חלוקת יחידה:
*: <math>\sum_{\nu = 0}^n b_{\nu, n}(x) = \sum_{\nu = 0}^n {n \choose \nu} x^\nu \left( 1 - x \right)^{n - \nu} = \left(x + \left( 1 - x \right) \right)^n = 1.</math>
* אם לוקחים את האיבר הראשון של <math>(x+y)^n</math> כאשר<math>y = 1-x</math>, ניתן להראות כי
*: <math>\sum_{\nu=0}^{n}\nu b_{\nu, n}(x) = nx</math>
* האיבר השני <math>(x+y)^n</math> כאשר <math>y = 1-x</math> משמש להראות כי 
*: <math>\sum_{\nu=1}^{n}\nu(\nu-1) b_{\nu, n}(x) = n(n-1)x^2</math>
* A פולינום ברנשטיין ניתן להיכתב כקומבינציה ליניארית של פולינומים מדרגה גבוהה יותר:
*: <math>b_{\nu, n - 1}(x) = \frac{n - \nu}{n} b_{\nu, n}(x) + \frac{\nu + 1}{n} b_{\nu + 1, n}(x).</math>

== קירוב פונקציות רציפות ==
נניח ש-ƒ פונקציה רציפה בקטע [0,1]. נבחן את פולינום ברנשטיין להלן:
: <math>B_n(f)(x) = \sum_{\nu = 0}^n f\left( \frac{\nu}{n} \right) b_{\nu,n}(x).</math>
ניתן להראות כי:
: <math>\lim_{n \to \infty}{ B_n(f)(x) } = f(x) \,</math>
מתכנס במידה שווה בקטע [0,1]. זו הצהרה חזקה יותר מאשר ההנחה שהגבול קיים בנפרד בכל ערך של x; זו תהיה התכנסות נקודתית במקום התכנסות במידה שווה. באופן ספציפי, המילים "במידה שווה" מציינות:
: <math>\lim_{n \to \infty} \sup \left\{\, \left| f(x) - B_n(f)(x) \right| \,:\, 0 \leq x \leq 1 \,\right\} = 0.</math>
פולינומי ברנשטיין לכן מאפשרים דרך אחת להוכיח את משפט הקירוב של ויירשטראס שכל פונקציה רציפה ממשית על הקטע הממשי [a,b] ניתנת להערכה במידה שווה על ידי פונקציות פולינום מעל R.

טענה כללית יותר על פונקציה שגזירה ברציפות k פעמים היא:
: <math>{\left\| B_n(f)^{(k)} \right\|}_\infty \le \frac{ (n)_k }{ n^k } \left\| f^{(k)} \right\|_\infty \text{ and } \left\| f^{(k)}- B_n(f)^{(k)} \right\|_\infty \to 0</math>
כאשר בנוסף
: <math>\frac{ (n)_k }{ n^k } = \left( 1 - \frac{0}{n} \right) \left( 1 - \frac{1}{n} \right) \cdots \left( 1 - \frac{k - 1}{n} \right)</math>
הוא ערך עצמי של B<sub>n</sub>; והפונקציה העצמית המקבילה היא פולינום מדרגה k.

=== הוכחה ===
נניח ש-K משתנה מקרי המפוזר כמספר ההצלחות מתוך n ניסויי ברנולי בלתי תלויים, עם הסתברות x של הצלחה בכל ניסוי; במילים אחרות, k מתפלג בינומית עם פרמטרים n ו-x. מכאן שערך התוחלת <math>E\left[K / n\right] = x</math>.

בעזרת הגרסה החלשה של חוק המספרים הגדולים מתורת ההסתברות,
: <math>\lim_{n \to \infty}{ P\left( \left| \frac{K}{n} - x \right|>\delta \right) } = 0</math>
לכל <math>\delta>0</math>. יתרה מכך, הקשר מתקיים במידה שווה ב-x, כפי שניתן להראות מההוכחה בעזרת אי-שוויון צ'בישב, אם לוקחים בחשבון שהשונות של <math>K/n</math>, השווה ל-<math>x(1-x)/n</math>, חסומה מלמעלה על ידי <math>1/4n</math> ללא תלות ב-x.

מכיוון ש-ƒ, הרציפה על תחום סגור וחסום, חייבת להיות רציפה במידה שווה על התחום, ניתן לנסח את הטענה כי
: <math>\lim_{n \to \infty}{ P\left( \left| f\left( \frac{K}{n} \right) - f\left( x \right) \right| > \varepsilon \right) } = 0</math>
במידה שווה ב-x. אם לוקחים בחשבון ƒ חסום (בקטע נתון), ניתן לצפות ל:
: <math>\lim_{n \to \infty}{ E\left( \left| f\left( \frac{K}{n} \right) - f\left( x \right) \right| \right) } = 0</math>
במידה שווה ב-x. כאן אפשר לראות כי התוחלת מתחלקת לשני חלקים. בחלק אחד ההפרש אינו עולה על ε; חלק זה יכול לתרום יותר מאשר ε. בעל החלק השני ההפרש עולה על ε, אך אינו עולה על 2M, כאשר M מהווה חסם עליון ל-|(ƒ(x)|; חלק זה אינו יכול לתרום יותר מ-2M פעמים את ההסתברות הקטנה כי ההפרש גדול מ-ε.

לסיכום, ראינו כי הערך המוחלט של ההפרש בין התוחלות לא עולה אף פעם על התוחלת של הערך המוחלט של ההפרש, וכן כי <math>E\left[f\left(K / n\right)\right]</math> הוא בדיוק פולינום ברנשטיין <math>B_n(f)(x)</math>.

ראו לדוגמה.

== ראו גם ==
* עקום בזייה
* צורת לגראנז'

==קישורים חיצוניים==

* 

== הערות שוליים ==

קטגוריה:אנליזה נומרית
קטגוריה:פולינומים