משתמש:אור שחף/133 - הרצאה/3.4.11: הבדלים בין גרסאות בדף
מ (←סוג א) |
|||
שורה 32: | שורה 32: | ||
* ולפי סימפסון: {{left|<math>\begin{array}{l}\frac h3\left(f(x_0)+4\sum_{k=1}^{2}f(x_{2k-1})+2\sum_{k=1}^{1}f(x_{2k})+f(x_4)\right)\\=\frac1{12}\left(1+4\frac45+2\frac23+4\frac47+\frac12\right)\\\approx\underline{0.693}253968\end{array}</math>}} נחשב את סדר הגודל של הטעות בקירוב סימפסון: {{left|<math>f'(x)=-x^{-2}\implies f''(x)=2x^{-3}\implies f^{(3)}(x)=-6x^{-4}\implies f^{(4)}(x)=24x^{-5}</math>}} ולכן <math>M=\max_{x\in[1,2]}\left|24x^{-5}\right|=24</math> והטעות R בקירוב מקיימת <math>|R|\le\frac{Mh^4}{180}(2-1)=\frac1{1920}<5.21\cdot10^{-4}</math> | * ולפי סימפסון: {{left|<math>\begin{array}{l}\frac h3\left(f(x_0)+4\sum_{k=1}^{2}f(x_{2k-1})+2\sum_{k=1}^{1}f(x_{2k})+f(x_4)\right)\\=\frac1{12}\left(1+4\frac45+2\frac23+4\frac47+\frac12\right)\\\approx\underline{0.693}253968\end{array}</math>}} נחשב את סדר הגודל של הטעות בקירוב סימפסון: {{left|<math>f'(x)=-x^{-2}\implies f''(x)=2x^{-3}\implies f^{(3)}(x)=-6x^{-4}\implies f^{(4)}(x)=24x^{-5}</math>}} ולכן <math>M=\max_{x\in[1,2]}\left|24x^{-5}\right|=24</math> והטעות R בקירוב מקיימת <math>|R|\le\frac{Mh^4}{180}(2-1)=\frac1{1920}<5.21\cdot10^{-4}</math> | ||
=אינטגרל לא אמיתי | =אינטגרל לא אמיתי, סוג I= | ||
עד עתה הגדרנו אינטגרלים מסויימים רק עבור פונקציות חסומות בקטעים סופיים. אם הפונקציה לא חסומה ו/או הקטע לא חסום עדיין ניתן להגדיר "אינטגרל לא אמיתי". | עד עתה הגדרנו אינטגרלים מסויימים רק עבור פונקציות חסומות בקטעים סופיים. אם הפונקציה לא חסומה ו/או הקטע לא חסום עדיין ניתן להגדיר "אינטגרל לא אמיתי" (improper integral). | ||
אינטגרלים של קטעים אינסופיים מהסוג <math>\int\limits_{-\infty}^b f,\ \int\limits_a^\infty f,\ \int\limits_{-\infty}^\infty f</math>. | אינטגרלים של קטעים אינסופיים מהסוג <math>\int\limits_{-\infty}^b f,\ \int\limits_a^\infty f,\ \int\limits_{-\infty}^\infty f</math>. | ||
שורה 47: | שורה 47: | ||
עבור f מוגדרת בכל <math>\mathbb R</math> נאמר שהיא אינטגרבילית מקומית אם היא אינטגרבילית מקומית בכל קטע סופי, ואם כן נגדיר <math>\int\limits_{-\infty}^\infty f=\int\limits_{-\infty}^a f+\int\limits_a^\infty f</math> עבור <math>a\in\mathbb R</math> כרצוננו עבורו שני האינטגרלים באגף ימין מתכנסים. | עבור f מוגדרת בכל <math>\mathbb R</math> נאמר שהיא אינטגרבילית מקומית אם היא אינטגרבילית מקומית בכל קטע סופי, ואם כן נגדיר <math>\int\limits_{-\infty}^\infty f=\int\limits_{-\infty}^a f+\int\limits_a^\infty f</math> עבור <math>a\in\mathbb R</math> כרצוננו עבורו שני האינטגרלים באגף ימין מתכנסים. | ||
==דוגמאות== | |||
# <math>\int\limits_1^\infty\frac{\mathrm dx}{x^2}</math>. נחשב: {{left|<math>\begin{align}\int&=\lim_{R\to\infty}\int\limits_1^R\frac{\mathrm dx}{x^2}\\&=\lim_{R\to\infty}\left[-\frac1x\right]_{x=1}^R\\&=\lim_{R\to\infty}\left(-\frac1R+1\right)\\&=1\end{align}</math>}} ניתן גם לכתוב בקיצור: <math>\int\limits_1^\infty\frac{\mathrm dx}{x^2}=\left[-\frac1x\right]_{x=1}^\infty=0-(-1)=1</math>. | # <math>\int\limits_1^\infty\frac{\mathrm dx}{x^2}</math>. נחשב: {{left|<math>\begin{align}\int&=\lim_{R\to\infty}\int\limits_1^R\frac{\mathrm dx}{x^2}\\&=\lim_{R\to\infty}\left[-\frac1x\right]_{x=1}^R\\&=\lim_{R\to\infty}\left(-\frac1R+1\right)\\&=1\end{align}</math>}} ניתן גם לכתוב בקיצור: <math>\int\limits_1^\infty\frac{\mathrm dx}{x^2}=\left[-\frac1x\right]_{x=1}^\infty=0-(-1)=1</math>. | ||
# <math>\int\limits_1^\infty\frac{\mathrm dx}x=[\ln|x|]_{x=1}^\infty=\ln(\infty)-0=\infty</math>, כלומר האינטגרל מתבדר (אך מתכנס במובן הרחב). | # <math>\int\limits_1^\infty\frac{\mathrm dx}x=[\ln|x|]_{x=1}^\infty=\ln(\infty)-0=\infty</math>, כלומר האינטגרל מתבדר (אך מתכנס במובן הרחב). |
גרסה מ־11:19, 2 במאי 2011
מבוא לאינטגרציה נומרית (המשך)
בהרצאה הקודמת הצגנו את כלל סימפסון לקירוב האינטגרל המסויים והראנו שהטעות בחישוב בקטע [math]\displaystyle{ [a,b] }[/math] חסומה ע"י [math]\displaystyle{ \max_{x\in[a,b]}\left|f^{(4)}(x)\right|\frac{h^4(b-a)}{45} }[/math] כאשר h המרחק בין שתי נקודות סמוכות בחלוקה שבחרנו. ניתן גישה אחרת למציאת הטעות, שהיא יותר קצרה ונותנת ערך יותר קטן לחסם של הטעות, אבל היא פחות אינטואיטיבית:
ע"י נירמול מספיק לחשב בקירוב [math]\displaystyle{ \int\limits_{-h}^h f\approx\frac h3(f(-h)+4f(0)+f(h)) }[/math] כאשר f גזירה 4 פעמים בסביבת 0. נגדיר פונקציה חדשה [math]\displaystyle{ G(h):=\int\limits_{-h}^h f-\frac h3(f(-h)+4f(0)+f(h)) }[/math]. G מוגדרת ורציפה בסביבה של 0 וגם (לפי הצבה) [math]\displaystyle{ G(0)=0 }[/math]. עבור F קדומה ל-f מתקיים (לפי המשפט היסודי)
לכן [math]\displaystyle{ \lim_{h\to0}G'(h)=f(0)+f(0)-\frac13(f(0)+4f(0)+f(0)-0=2f(0)-\frac63f(0)=0 }[/math]. ע"פ הלמה השנייה בהרצאה הקודמת [math]\displaystyle{ G'(0) }[/math] קיים ושווה ל-0. נגזור שוב את G ונקבל [math]\displaystyle{ \frac{\mathrm d^2}{\mathrm dh^2}G(h)=\frac13(f'(h)-f'(-h))-\frac h3(f''(-h)+f''(h)) }[/math]. מכאן ש-[math]\displaystyle{ \lim_{h\to0}G''(h)=\frac13(f'(0)-f'(0))-0=0=G''(0) }[/math]. נמשיך לגזור פעמיים נוספות ונקבל [math]\displaystyle{ G'''(0)=G^{(4)}(0)=0 }[/math] וגם [math]\displaystyle{ G^{(4)}(h)=-\frac13(-f'''(-h)+f'''(h))-\frac h3(f^{(4)}(-h)+f^{(4)}(h)) }[/math]. עתה:
[math]\displaystyle{ }[/math] | [math]\displaystyle{ \frac{G(h)-G(0)}{h^5-0^5} }[/math] | [math]\displaystyle{ = }[/math] | [math]\displaystyle{ \frac{G(h)}{h^5} }[/math] | [math]\displaystyle{ }[/math] | ||
לפי משפט קושי קיים [math]\displaystyle{ h_1\in(0,h) }[/math] עבורו: | [math]\displaystyle{ }[/math] | [math]\displaystyle{ \frac{G'(h_1)}{5h_1^4} }[/math] | [math]\displaystyle{ = }[/math] | [math]\displaystyle{ }[/math] | [math]\displaystyle{ }[/math] | |
[math]\displaystyle{ }[/math] | [math]\displaystyle{ \frac{G'(h_1)-G'(0)}{5h_1^4-5\cdot0^4} }[/math] | [math]\displaystyle{ = }[/math] | [math]\displaystyle{ }[/math] | [math]\displaystyle{ }[/math] | ||
קיים [math]\displaystyle{ h_2\in(0,h_1) }[/math] עבורו: | [math]\displaystyle{ }[/math] | [math]\displaystyle{ \frac{G''(h_2)}{20h_2^3} }[/math] | [math]\displaystyle{ = }[/math] | [math]\displaystyle{ }[/math] | [math]\displaystyle{ }[/math] | |
קיים [math]\displaystyle{ h_3\in(0,h_2) }[/math] עבורו: | [math]\displaystyle{ }[/math] | [math]\displaystyle{ \frac{G'''(h_3)}{60h_3^2} }[/math] | [math]\displaystyle{ = }[/math] | [math]\displaystyle{ }[/math] | [math]\displaystyle{ }[/math] | |
קיים [math]\displaystyle{ h_4\in(0,h_3) }[/math] עבורו: | [math]\displaystyle{ }[/math] | [math]\displaystyle{ \frac{G^{(4)}(h_4)}{120h_4} }[/math] | [math]\displaystyle{ = }[/math] | [math]\displaystyle{ }[/math] | [math]\displaystyle{ }[/math] |
כעת נגדיר [math]\displaystyle{ M:=\max_{x\in[a,b]}\left|f^{(4)}(x)\right| }[/math]. לפי משפט לגראנז' קיים [math]\displaystyle{ c\in(-h,h) }[/math] כך ש-[math]\displaystyle{ \frac{f^{(3)}(h)-f^{(3)}(-h)}{2h}=f^{(4)}(c) }[/math] ולכן [math]\displaystyle{ \left|f^{(3)}(h)-f^{(3)}(-h)\right|\le2hM }[/math]. מכל זה נובע
עתה [math]\displaystyle{ \left|\frac{G(h)}{h^5}\right|=\left|\frac{G^{(4)}(h_4)}{120h_4}\right|\le\frac1{120h_4}\frac43Mh_4=\frac M{90} }[/math] וקיבלנו ש-[math]\displaystyle{ |G(h)|\le\frac{Mh^5}{90} }[/math], כלומר הטעות בכל קטע מהסוג [math]\displaystyle{ [x_{k-1},x_{k+1}] }[/math] חסומה ע"י [math]\displaystyle{ \frac{Mh^5}{90} }[/math]. ב-[math]\displaystyle{ [a,b] }[/math] יש [math]\displaystyle{ \frac n2=\frac{(b-a)}{2h} }[/math] ולפיכך הטעות חסומה ע"י [math]\displaystyle{ \frac{b-a}{2h}\cdot\frac{Mh^5}{90}=\frac{Mh^4}{180}(b-a) }[/math]. [math]\displaystyle{ \blacksquare }[/math]
דוגמה
נקרב [math]\displaystyle{ \int\limits_1^2\frac{\mathrm dx}x=\ln(2)\approx0.69314718 }[/math]. נבחר [math]\displaystyle{ h=\frac14 }[/math]. נציב:
- הקירוב לפי סכום רימן הוא [math]\displaystyle{ \sum_{k=0}^4f(x_k)h=\frac14\left(1+\frac45+\frac23+\frac47+\frac12\right)\approx\underline{0.6}34523809 }[/math].
- כעת נעשה קירוב בשיטת הטרפזים: [math]\displaystyle{ \frac{f(x_0)+f(x_4)}2h+h\sum_{k=1}^3f(x_k)=\frac18\left(1+\frac12\right)+\frac14\left(\frac45+\frac23+\frac47\right)\approx\underline{0.69}7023792 }[/math]
- ולפי סימפסון: [math]\displaystyle{ \begin{array}{l}\frac h3\left(f(x_0)+4\sum_{k=1}^{2}f(x_{2k-1})+2\sum_{k=1}^{1}f(x_{2k})+f(x_4)\right)\\=\frac1{12}\left(1+4\frac45+2\frac23+4\frac47+\frac12\right)\\\approx\underline{0.693}253968\end{array} }[/math]נחשב את סדר הגודל של הטעות בקירוב סימפסון:[math]\displaystyle{ f'(x)=-x^{-2}\implies f''(x)=2x^{-3}\implies f^{(3)}(x)=-6x^{-4}\implies f^{(4)}(x)=24x^{-5} }[/math]ולכן [math]\displaystyle{ M=\max_{x\in[1,2]}\left|24x^{-5}\right|=24 }[/math] והטעות R בקירוב מקיימת [math]\displaystyle{ |R|\le\frac{Mh^4}{180}(2-1)=\frac1{1920}\lt 5.21\cdot10^{-4} }[/math]
אינטגרל לא אמיתי, סוג I
עד עתה הגדרנו אינטגרלים מסויימים רק עבור פונקציות חסומות בקטעים סופיים. אם הפונקציה לא חסומה ו/או הקטע לא חסום עדיין ניתן להגדיר "אינטגרל לא אמיתי" (improper integral).
אינטגרלים של קטעים אינסופיים מהסוג [math]\displaystyle{ \int\limits_{-\infty}^b f,\ \int\limits_a^\infty f,\ \int\limits_{-\infty}^\infty f }[/math].
הגדרה: תהי f פונקציה מוגדרת בקטע מהסוג [math]\displaystyle{ [a,\infty) }[/math]. נאמר ש-f אינטגרבילית מקומית (locally integrable) בקטע זה אם לכל [math]\displaystyle{ b\gt a }[/math] f אינטגרבילית בקטע [math]\displaystyle{ [a,b] }[/math].
למשל, אם f רציפה למקוטעין אז היא אינטגרבילית מקומית.
הגדרה: תהי f מוגדרת ואינטגרבילית מקומית ב-[math]\displaystyle{ [a,\infty) }[/math] נגדיר [math]\displaystyle{ \int\limits_a^\infty f:=\lim_{R\to\infty}\int\limits_a^R f }[/math]. אם הגבול קיים נאמר שהאינטגרל מתכנס, אחרת הוא מתבדר.
אינטגביליות מקומית מוגדרת באופן דומה עבור קטע מהצורה [math]\displaystyle{ (-\infty,b] }[/math] ואם f אינטגרבילית מקומית שם נגדיר [math]\displaystyle{ \int\limits_{-\infty}^b f:=\lim_{R\to-\infty}\int\limits_R^b f }[/math].
עבור f מוגדרת בכל [math]\displaystyle{ \mathbb R }[/math] נאמר שהיא אינטגרבילית מקומית אם היא אינטגרבילית מקומית בכל קטע סופי, ואם כן נגדיר [math]\displaystyle{ \int\limits_{-\infty}^\infty f=\int\limits_{-\infty}^a f+\int\limits_a^\infty f }[/math] עבור [math]\displaystyle{ a\in\mathbb R }[/math] כרצוננו עבורו שני האינטגרלים באגף ימין מתכנסים.
דוגמאות
- [math]\displaystyle{ \int\limits_1^\infty\frac{\mathrm dx}{x^2} }[/math]. נחשב: [math]\displaystyle{ \begin{align}\int&=\lim_{R\to\infty}\int\limits_1^R\frac{\mathrm dx}{x^2}\\&=\lim_{R\to\infty}\left[-\frac1x\right]_{x=1}^R\\&=\lim_{R\to\infty}\left(-\frac1R+1\right)\\&=1\end{align} }[/math]ניתן גם לכתוב בקיצור: [math]\displaystyle{ \int\limits_1^\infty\frac{\mathrm dx}{x^2}=\left[-\frac1x\right]_{x=1}^\infty=0-(-1)=1 }[/math].
- [math]\displaystyle{ \int\limits_1^\infty\frac{\mathrm dx}x=[\ln|x|]_{x=1}^\infty=\ln(\infty)-0=\infty }[/math], כלומר האינטגרל מתבדר (אך מתכנס במובן הרחב).
- שאלה מארה"ב מלפני הרבה שנים: חצוצרה אינסופית תתקבל מסיבוב הגרף [math]\displaystyle{ y=\frac1x }[/math] סביב ציר ה-x ב-[math]\displaystyle{ [1,\infty) }[/math]. איך צובעים אותו מבפנים?
פתרון: לכאורה זה בלתי אפשרי, כי שטח הפנים של החצוצרה הוא [math]\displaystyle{ \pi\int\limits_1^\infty\frac2x\sqrt{1+\left(-\frac1{x^2}\right)^2}\mathrm dx\gt \pi\int\limits_1^\infty\frac2x\mathrm dx=\infty }[/math], כלומר אין מספיק צבע בעולם. אך מכיוון שכמות הצבע נמדדת ביחידות נפח ולא שטח, ומכיוון שהנפח בתוך החצוצרה הוא [math]\displaystyle{ \pi\int\limits_1^\infty\frac{\mathrm dx}{x^2}=\pi }[/math], יספיקו לנו [math]\displaystyle{ \pi }[/math] יחידות מעוקבות של צבע ואפילו ישאר לנו עודף.
שאלה: האם התכנסות האינטגרל [math]\displaystyle{ \int\limits_1^\infty f }[/math] גוררת ש-[math]\displaystyle{ \lim_{x\to\infty} f(x)=0 }[/math] (בדומה לטורים)?
תשובה: לא. נגדיר פונקציה f שהגרף שלה הוא
אזי
כלומר האינטגרל מתכנס, אבל [math]\displaystyle{ \lim_{x\to\infty} f(x) }[/math] לא קיים.