Matematyka i fizyka od podstaw
PLEN

Matematyka › Część I. Wektory i analiza jednej zmiennej › 1.6 gotowa · PDF do druku

Zanim zaczniesz 1.2 Ciągi, granice i ciągłość1.5 Całki

Nowe pojęcia w tej lekcji (4)

szereg Tayloraszereg geometrycznyszereg liczbowy i jego zbieżnośćszereg potęgowy, promień zbieżności

Matematyka 1.6

Szeregi liczbowe i potęgowe

Czy można dodać nieskończenie wiele liczb? Czasem wychodzi liczba skończona, jak w paradoksie Zenona, a czasem nie, choć dodawane liczby maleją do zera. Ta lekcja uczy rozpoznawać, kiedy suma ma sens: od szeregu geometrycznego, przez zaskakująco rozbieżny szereg harmoniczny (i stos klocków wystający dowolnie daleko poza stół), do szeregów potęgowych, które zamieniają funkcje w nieskończone wielomiany. Na końcu policzysz z nich ee, π\pi i ln⁡2\ln2.

Czas
ok. 3 h teorii + 3 h zadań
Zadania
27 zadań w treści, z podpowiedziami i odpowiedziami
Źródła
Feynmana wykłady z fizyki, t. I, 8-2, 15-8; K. A. Stroud, Engineering Mathematics

0 Prolog: piłka, która odbija się nieskończenie wiele razy

Piłka spada z wysokości 2 m, a po każdym odbiciu wznosi się na 60% poprzedniej wysokości. W idealnym modelu odbija się nieskończenie wiele razy. Czy w takim razie skacze bez końca? To ten sam problem, co u Zenona: nieskończenie wiele kroków w rozumowaniu nie musi oznaczać nieskończenie długiego czasu.FLP I-8-2

Pierwszy spadek trwa t0=2h0/g≈0,639t_0 = \sqrt{2h_0/g} \approx 0{,}639 s. Wysokość kolejnego odbicia to 0,60{,}6 poprzedniej, a czas spadku z wysokości hh to 2h/g\sqrt{2h/g} (z h=12gt2h = \tfrac12gt^2, lekcja 0.1 Liczby, potęgi, logarytmy i funkcje elementarne), więc każdy następny lot w górę i w dół trwa 0,6≈0,775\sqrt{0{,}6} \approx 0{,}775 długości poprzedniego (pierwszy taki lot trwa 2t0q2t_0q). Całkowity czas to

T=t0+2t0q+2t0q2+2t0q3+…,q=0,6.T = t_0 + 2t_0q + 2t_0q^2 + 2t_0q^3 + \dots,\qquad q = \sqrt{0{,}6}.

To nieskończona suma. W tej lekcji nauczysz się ją policzyć i sprawdzić, czy w ogóle ma sens.

Pytanie na start. Po jakim czasie piłka przestaje skakać? Odpowiedź w przykładzie 1.1.

1 Szereg i jego suma

Szereg to nieskończona suma a1+a2+a3+⋯=∑n=1∞ana_1 + a_2 + a_3 + \dots = \sum_{n=1}^{\infty} a_n. Sensu nadaje mu granica. Dodajemy kolejno skończenie wiele wyrazów i dostajemy sumy częściowe

SN=∑n=1Nan=a1+a2+⋯+aN.S_N = \sum_{n=1}^{N} a_n = a_1 + a_2 + \dots + a_N .

Jeśli ciąg SNS_N ma skończoną granicę SS (lekcja 1.2 Ciągi, granice i ciągłość), szereg jest zbieżny, a SS to jego suma. W przeciwnym razie szereg jest rozbieżny.

1.1 Szereg geometryczny

Najważniejszy szereg ma stały iloraz kolejnych wyrazów: a+aq+aq2+…a + aq + aq^2 + \dots. Sumę częściową znamy z lekcji 0.2 Równania, nierówności, wielomiany i ułamki proste: SN=a1−qN1−qS_N = a\dfrac{1 - q^N}{1 - q} dla q≠1q \neq 1. Dla ∣q∣<1|q| < 1 mamy qN→0q^N \to 0, więc

Szereg geometryczny

∑n=0∞aqn=a+aq+aq2+⋯=a1−qdla ∣q∣<1.\sum_{n=0}^{\infty} aq^n = a + aq + aq^2 + \dots = \frac{a}{1 - q}\qquad\text{dla } |q| < 1 .

Dla ∣q∣≥1|q| \ge 1 (i a≠0a \neq 0) szereg jest rozbieżny.

Zenon: 100+10+1+⋯=1001−0,1=111,1‾100 + 10 + 1 + \dots = \dfrac{100}{1 - 0{,}1} = 111{,}\overline1 m. Ułamki okresowe to też szeregi geometryczne: 0,9‾=910+9100+⋯=9/101−1/10=10{,}\overline{9} = \tfrac{9}{10} + \tfrac{9}{100} + \dots = \dfrac{9/10}{1 - 1/10} = 1. Liczba 0,999…0{,}999\ldots to dokładnie 1, a nie „prawie 1”.

Symulacja 1 Sumy częściowe SNS_N kilku szeregów. Szereg geometryczny szybko się stabilizuje. Harmoniczny rośnie tak wolno, że na wykresie wydaje się zbieżny, a jednak nie ma granicy. Naprzemienny zbliża się do swojej sumy, skacząc raz powyżej, raz poniżej.

Przykład 1.1 Ile trwa skakanie piłki?

Dla piłki z prologu (h0=2h_0 = 2 m, każde odbicie na 60% wysokości) oblicz całkowitą drogę i całkowity czas skakania.

  1. Droga (r=0,6r = 0{,}6; w lekcji 0.2 Równania, nierówności, wielomiany i ułamki proste po pięciu wzlotach wyszło 7,533 m): spadek h0h_0, potem każde odbicie to lot w górę i w dół: h0+2h0r+2h0r2+⋯=h0+2h0r1−r=2+2,40,4=8h_0 + 2h_0r + 2h_0r^2 + \dots = h_0 + \dfrac{2h_0r}{1 - r} = 2 + \dfrac{2{,}4}{0{,}4} = 8 m.
  2. Czas: T=t0+2t0q1−qT = t_0 + \dfrac{2t_0q}{1 - q} z q=0,6≈0,7746q = \sqrt{0{,}6} \approx 0{,}7746 i t0=4/9,81≈0,6386t_0 = \sqrt{4/9{,}81} \approx 0{,}6386 s.
  3. T≈0,6386+2⋅0,6386⋅0,77460,2254≈5,03T \approx 0{,}6386 + \dfrac{2\cdot0{,}6386\cdot0{,}7746}{0{,}2254} \approx 5{,}03 s.

Wynik: piłka przebywa 8 m i przestaje skakać po ok. 5,0 s, choć w modelu odbija się nieskończenie wiele razy.

Sprawdzenie: po 20 odbiciach do końca zostaje tylko 0,027 s, więc ogon szeregu naprawdę znika ✓. Prawdziwa piłka przestaje odbijać się wcześniej, bo przy małych prędkościach model stałego współczynnika odbicia przestaje działać.

Zadania · część 1

5 zadań · ok. 25 min

  1. Oblicz sumy: 1+13+19+…1 + \tfrac13 + \tfrac19 + \dots;     4−2+1−12+…\;\;4 - 2 + 1 - \tfrac12 + \dots;     ∑n=1∞(25)n\;\;\sum_{n=1}^{\infty}\left(\tfrac{2}{5}\right)^n.

    Odpowiedź

    32\tfrac32;   41+1/2=83\;\tfrac{4}{1 + 1/2} = \tfrac83;   2/51−2/5=23\;\tfrac{2/5}{1 - 2/5} = \tfrac23 (pierwszy wyraz to 25\tfrac25).

  2. Zamień na ułamek zwykły: 0,3‾0{,}\overline{3},   0,12‾\;0{,}\overline{12},   2,15‾\;2{,}1\overline{5}.

    Odpowiedź

    13\tfrac13;   1299=433\;\tfrac{12}{99} = \tfrac{4}{33};   2,1+0,051−0,1=2,1+118=9745\;2{,}1 + \tfrac{0{,}05}{1 - 0{,}1} = 2{,}1 + \tfrac{1}{18} = \tfrac{97}{45}.

  3. Dla jakich xx zbieżny jest szereg 1+2x+4x2+8x3+…1 + 2x + 4x^2 + 8x^3 + \dots? Ile wynosi jego suma?

    Odpowiedź

    Iloraz q=2xq = 2x, więc ∣x∣<12|x| < \tfrac12; suma 11−2x\tfrac{1}{1 - 2x}.

  4. Lek podaje się co 24 godziny w dawce 100 mg, a w ciągu doby organizm usuwa 75% leku. Ile leku jest w organizmie tuż po bardzo wielu dawkach?

    Wskazówka

    Tuż po NN-tej dawce: 100+100⋅0,25+⋯+100⋅0,25N−1100 + 100\cdot0{,}25 + \dots + 100\cdot0{,}25^{N-1}.

    Odpowiedź

    1001−0,25≈133\dfrac{100}{1 - 0{,}25} \approx 133 mg. Poziom ustala się szybko: po trzech dawkach 131 mg.

  5. Wyjaśnij przez sumy częściowe, dlaczego szereg 1−1+1−1+…1 - 1 + 1 - 1 + \dots nie ma sumy, choć „grupując” wyrazy można dostać (1−1)+(1−1)+⋯=0(1 - 1) + (1 - 1) + \dots = 0 albo 1−(1−1)−⋯=11 - (1 - 1) - \dots = 1.

    Odpowiedź

    Sumy częściowe to 1,0,1,0,…1, 0, 1, 0, \dots, a ten ciąg nie ma granicy. Grupowanie wybiera tylko podciąg sum częściowych, więc może dać dowolną z dwóch wartości. Dla szeregów rozbieżnych zwykłe reguły arytmetyki nie działają.

2 Szereg harmoniczny i stos klocków

Warunek konieczny zbieżności jest oczywisty: jeśli szereg jest zbieżny, to an→0a_n \to 0. Inaczej sumy częściowe wciąż przyrastałyby o coś niemałego. Ale ten warunek nie wystarcza. Pokazuje to szereg harmoniczny:

1+12+13+14+⋯=∞.1 + \frac12 + \frac13 + \frac14 + \dots = \infty .

Wyrazy dążą do zera, a suma jest nieskończona. Dowód z XIV wieku (Mikołaj z Oresme) grupuje wyrazy:

1+12+13+14⏟>2⋅14=12+15+⋯+18⏟>4⋅18=12+19+⋯+116⏟>12+…1 + \frac12 + \underbrace{\frac13 + \frac14}_{> 2\cdot\frac14 = \frac12} + \underbrace{\frac15 + \dots + \frac18}_{> 4\cdot\frac18 = \frac12} + \underbrace{\frac19 + \dots + \frac1{16}}_{> \frac12} + \dots

Każda grupa daje więcej niż 12\tfrac12, a grup jest nieskończenie wiele. Szereg rośnie jednak bardzo wolno: suma częściowa HNH_N jest bliska ln⁡N+0,577\ln N + 0{,}577. Żeby przekroczyć 10, trzeba 12 367 wyrazów, a do 100 potrzeba ok. 104310^{43}.

Szereg harmoniczny ma piękne zastosowanie fizyczne. Układamy na krawędzi stołu stos jednakowych klocków o długości LL, każdy przesunięty względem niższego. Stos stoi, dopóki dla każdego klocka (i dla stołu) środek masy wszystkich klocków leżących wyżej nie wychodzi poza jego krawędź. Przy największym wysięgu leży dokładnie na krawędzi: to równowaga graniczna, w praktyce zostawiamy mały zapas. Jeśli na każdym poziomie leży jeden klocek, największy wysięg daje ułożenie, w którym najwyższy klocek wystaje o L2\tfrac{L}{2} poza drugi, drugi o L4\tfrac{L}{4} poza trzeci, trzeci o L6\tfrac{L}{6} i tak dalej. Całkowity wysięg nn klocków to

dn=L2(1+12+13+⋯+1n)=L2Hn.d_n = \frac L2\left(1 + \frac12 + \frac13 + \dots + \frac1n\right) = \frac L2 H_n .

Już 4 klocki wystają poza stół o więcej niż długość klocka (d4=1,04 Ld_4 = 1{,}04\,L). A ponieważ szereg harmoniczny jest rozbieżny, odpowiednio wysoki stos może wystawać dowolnie daleko. Tylko że na wysięg 3L3L potrzeba już 227 klocków. (Przy kilku klockach na jednym poziomie, z przeciwwagami, stos da się wysunąć jeszcze dalej.)

Symulacja 2 Stos klocków na krawędzi stołu (pionowa linia). Każdy klocek wysunięty o L2k\tfrac{L}{2k} względem niższego. Zwiększaj liczbę klocków: wysięg rośnie bez ograniczeń, ale coraz wolniej, jak L2ln⁡n\tfrac L2\ln n.

Zadania · część 2

4 zadania · ok. 25 min

  1. Który szereg na pewno jest rozbieżny na podstawie warunku koniecznego: ∑nn+1\sum\dfrac{n}{n + 1},   ∑1n2\;\sum\dfrac1{n^2},   ∑cos⁡1n\;\sum\cos\dfrac1n?

    Odpowiedź

    Pierwszy (an→1a_n \to 1) i trzeci (an→1a_n \to 1). Dla drugiego an→0a_n \to 0, co nic nie przesądza (jest zbieżny, część 3).

  2. Ile wystają poza stół stosy 2 i 3 klocków o długości 20 cm?

    Odpowiedź

    10⋅1,5=1510\cdot1{,}5 = 15 cm;   10⋅1,833≈18,3\;10\cdot1{,}833 \approx 18{,}3 cm.

  3. Uzasadnij wysięg L2k\tfrac{L}{2k}: środek masy kk górnych klocków musi leżeć dokładnie na krawędzi klocka k+1k + 1. Sprawdź to dla k=1k = 1 i k=2k = 2.

    Wskazówka

    Mierz położenia środków klocków od krawędzi klocka poniżej. Środek masy kk jednakowych klocków to średnia położeń ich środków.

    Odpowiedź

    k=1k = 1: środek górnego klocka musi być nad krawędzią drugiego, więc górny wystaje o L2\tfrac L2. k=2k = 2: wspólny środek masy dwóch górnych ma być nad krawędzią trzeciego. Jeśli drugi wystaje o dd, jego środek jest w d−L2d - \tfrac L2, a środek pierwszego w dd. Średnia d−L4=0d - \tfrac L4 = 0 daje d=L4d = \tfrac L4 ✓.

  4. Grupując wyrazy jak Oresme, pokaż, że H2m≥1+m2H_{2^m} \ge 1 + \tfrac m2 (ostro dla m≥2m \ge 2). Ile wyrazów wystarczy, żeby suma przekroczyła 6?

    Odpowiedź

    Po 2m2^m wyrazach mamy 1, 12\tfrac12 i m−1m - 1 grup, każda większa niż 12\tfrac12: razem więcej niż 1+m21 + \tfrac m2. Dla m=10m = 10: H1024>6H_{1024} > 6. (Naprawdę wystarczy 227 wyrazów; oszacowanie jest zgrubne, ale wystarcza do dowodu rozbieżności.)

3 Kryteria zbieżności

Zwykle nie umiemy policzyć sumy, ale umiemy rozstrzygnąć, czy istnieje. Dla szeregów o wyrazach dodatnich służą do tego trzy narzędzia.

3.1 Porównanie

Jeśli 0≤an≤bn0 \le a_n \le b_n i ∑bn\sum b_n jest zbieżny, to ∑an\sum a_n też. Jeśli ∑an\sum a_n jest rozbieżny, to ∑bn\sum b_n też. Szereg mniejszy od zbieżnego jest zbieżny, większy od rozbieżnego rozbieżny. Porównujemy zwykle z szeregiem geometrycznym albo z szeregami ∑1np\sum\tfrac{1}{n^p} z części 3.2.

3.2 Kryterium całkowe i szeregi ∑1np\sum\frac1{n^p}

Suma ∑f(n)\sum f(n) dla dodatniej, malejącej funkcji ff to pole słupków o szerokości 1. Słupek nad [n; n+1][n;\ n + 1] ma wysokość f(n)f(n), więc leży nad wykresem, a przesunięty o 1 w lewo leży pod nim: ∫1N+1f(x) dx≤SN≤f(1)+∫1Nf(x) dx\int_1^{N+1} f(x)\,\mathrm dx \le S_N \le f(1) + \int_1^{N} f(x)\,\mathrm dx. Dlatego szereg i całka niewłaściwa ∫1∞f(x) dx\int_1^\infty f(x)\,\mathrm dx są jednocześnie zbieżne albo rozbieżne (lekcja 1.5 Całki), o ile ff jest dodatnia i malejąca na [1; ∞)[1;\ \infty). Stąd:

Szeregi pp

∑n=1∞1np  jest zbiez˙ny dla p>1 i rozbiez˙ny dla p≤1.\sum_{n=1}^{\infty}\frac{1}{n^p}\ \text{ jest zbieżny dla } p > 1 \text{ i rozbieżny dla } p \le 1 .

Harmoniczny to przypadek graniczny p=1p = 1, odpowiadający całce ∫dxx=ln⁡x→∞\int\tfrac{\mathrm dx}{x} = \ln x \to \infty. Euler pokazał, że ∑1n2=π26≈1,6449\sum\tfrac{1}{n^2} = \tfrac{\pi^2}{6} \approx 1{,}6449.

3.3 Kryterium d'Alemberta (ilorazowe)

Policz granicę ilorazu kolejnych wyrazów q=lim⁡n→∞an+1anq = \lim\limits_{n\to\infty}\dfrac{a_{n+1}}{a_n}. Dla q<1q < 1 szereg jest zbieżny (od pewnego miejsca maleje szybciej niż geometryczny), dla q>1q > 1 (także q=∞q = \infty) rozbieżny, a dla q=1q = 1 kryterium nic nie mówi. To najwygodniejsze kryterium dla szeregów z silnią i potęgami.

Przykład 3.1 Trzy szeregi, trzy kryteria

Zbadaj zbieżność: (a) ∑n2n\sum\dfrac{n}{2^n}; (b) ∑1n2+n\sum\dfrac{1}{n^2 + n}; (c) ∑1n\sum\dfrac{1}{\sqrt n}.

  1. (a) d'Alembert: an+1an=n+12n→12<1\dfrac{a_{n+1}}{a_n} = \dfrac{n + 1}{2n} \to \dfrac12 < 1. Zbieżny (suma wynosi 2).
  2. (b) Porównanie: 1n2+n<1n2\dfrac{1}{n^2 + n} < \dfrac{1}{n^2}, a ∑1n2\sum\tfrac1{n^2} jest zbieżny. Zbieżny. Suma: 1n(n+1)=1n−1n+1\tfrac{1}{n(n + 1)} = \tfrac1n - \tfrac1{n + 1} (ułamki proste, lekcja 0.2 Równania, nierówności, wielomiany i ułamki proste), więc SN=1−1N+1→1S_N = 1 - \tfrac1{N + 1} \to 1.
  3. (c) Szereg pp z p=12≤1p = \tfrac12 \le 1. Rozbieżny.

Wynik: (a) i (b) zbieżne, z sumami 2 i 1; (c) rozbieżny.

Sprawdzenie: dla (c) d'Alembert daje n/(n+1)→1\sqrt{n/(n + 1)} \to 1, czyli nic nie rozstrzyga; dlatego potrzebne było inne kryterium ✓. Suma (b) „składa się jak teleskop”: wyrazy środkowe się znoszą.

Zadania · część 3

5 zadań · ok. 35 min

  1. Zbadaj zbieżność: ∑1n3\sum\dfrac{1}{n^3};   ∑1n0,9\;\sum\dfrac{1}{n^{0{,}9}};   ∑3nn!\;\sum\dfrac{3^n}{n!}.

    Odpowiedź

    Zbieżny (p=3p = 3); rozbieżny (p=0,9p = 0{,}9); zbieżny (d'Alembert: 3n+1→0\tfrac{3}{n + 1} \to 0).

  2. Zbadaj kryterium d'Alemberta zbieżność szeregu ∑n!10n\sum\dfrac{n!}{10^n}.

    Odpowiedź

    n+110→∞\tfrac{n + 1}{10} \to \infty: rozbieżny. Silnia w końcu pokonuje każdą funkcję wykładniczą. Zresztą an→∞a_n \to \infty, więc wystarczy warunek konieczny.

  3. Zbadaj ∑2n+1n3+5\sum\dfrac{2n + 1}{n^3 + 5} przez porównanie.

    Wskazówka

    Dla dużych nn wyraz zachowuje się jak 2nn3=2n2\tfrac{2n}{n^3} = \tfrac{2}{n^2}.

    Odpowiedź

    2n+1n3+5≤3nn3=3n2\dfrac{2n + 1}{n^3 + 5} \le \dfrac{3n}{n^3} = \dfrac{3}{n^2} dla n≥1n \ge 1, a ∑3n2\sum\tfrac3{n^2} jest zbieżny. Zbieżny.

  4. Oblicz sumę teleskopową ∑n=1∞2n(n+2)\sum_{n=1}^{\infty}\dfrac{2}{n(n + 2)}.

    Odpowiedź

    2n(n+2)=1n−1n+2\tfrac{2}{n(n + 2)} = \tfrac1n - \tfrac1{n + 2}. Zostają 1+121 + \tfrac12, więc suma 32\tfrac32.

  5. Pokaż kryterium całkowym, że ∑1nln⁡n\sum\dfrac{1}{n\ln n} (od n=2n = 2) jest rozbieżny, choć jego wyrazy maleją szybciej niż 1n\tfrac1n.

    Odpowiedź

    ∫2Xdxxln⁡x=ln⁡ln⁡X−ln⁡ln⁡2→∞\int_2^X\dfrac{\mathrm dx}{x\ln x} = \ln\ln X - \ln\ln2 \to \infty (podstawienie u=ln⁡xu = \ln x, lekcja 1.5 Całki). Rozbieżny, ale niewyobrażalnie wolno: żeby suma przekroczyła 10, liczba wyrazów musiałaby mieć tysiące cyfr.

4 Szeregi naprzemienne

Gdy znaki wyrazów się zmieniają, szereg zbiega łatwiej, bo wyrazy częściowo się znoszą. Kryterium Leibniza: jeśli bnb_n maleje do zera, to szereg b1−b2+b3−b4+…b_1 - b_2 + b_3 - b_4 + \dots jest zbieżny. Co więcej, błąd sumy częściowej jest mniejszy od pierwszego pominiętego wyrazu: ∣S−SN∣<bN+1|S - S_N| < b_{N+1}. Tę regułę wykorzystaliśmy już w lekcji 1.4 Rachunek różniczkowy w praktyce.

Przykład: naprzemienny szereg harmoniczny jest zbieżny, choć harmoniczny nie jest:

1−12+13−14+⋯=ln⁡2≈0,6931.1 - \frac12 + \frac13 - \frac14 + \dots = \ln2 \approx 0{,}6931 .

Takie szeregi, zbieżne tylko dzięki zmianie znaków, nazywamy warunkowo zbieżnymi. Są delikatne: Riemann pokazał, że przestawiając ich wyrazy, można dostać dowolną sumę. Szeregi zbieżne także po wzięciu wartości bezwzględnych (bezwzględnie zbieżne) można przestawiać bez obaw.

Zadania · część 4

3 zadania · ok. 20 min

  1. Ile wyrazów szeregu 1−12+13−…1 - \tfrac12 + \tfrac13 - \dots wystarczy według oszacowania Leibniza, żeby oszacować ln⁡2\ln2 z błędem mniejszym niż 0,01?

    Odpowiedź

    1N+1<0,01\tfrac{1}{N + 1} < 0{,}01, czyli N=100N = 100. Zbieżność jest wolna. Naprawdę wystarczy 50 wyrazów (błąd 0,0099), bo oszacowanie jest z zapasem. Z 10 wyrazów wychodzi 0,6456, błąd 0,048.

  2. Które szeregi są bezwzględnie, a które warunkowo zbieżne: ∑(−1)nn2\sum\dfrac{(-1)^n}{n^2};   ∑(−1)nn\;\sum\dfrac{(-1)^n}{\sqrt n}?

    Odpowiedź

    Pierwszy bezwzględnie (∑1n2\sum\tfrac1{n^2} zbieżny). Drugi warunkowo: Leibniz działa, ale ∑1n\sum\tfrac1{\sqrt n} jest rozbieżny.

  3. Oszacuj ∑n=1∞(−1)n+1n3\sum_{n=1}^\infty\dfrac{(-1)^{n+1}}{n^3} z dokładnością do 0,001.

    Odpowiedź

    1(N+1)3<0,001\tfrac{1}{(N + 1)^3} < 0{,}001 dla N+1=11N + 1 = 11, więc wystarczy 10 wyrazów: ≈0,9011\approx 0{,}9011 (dokładnie 0,901 540{,}901\,54).

5 Szeregi potęgowe i szereg Taylora

Szereg potęgowy to „nieskończony wielomian” ∑n=0∞cnxn=c0+c1x+c2x2+…\sum_{n=0}^\infty c_nx^n = c_0 + c_1x + c_2x^2 + \dots (albo z potęgami (x−a)(x - a)). Dla każdego xx to zwykły szereg liczbowy, który może być zbieżny albo nie. Okazuje się, że zbiór dobrych xx to zawsze przedział: istnieje promień zbieżności RR taki, że szereg jest zbieżny dla ∣x∣<R|x| < R i rozbieżny dla ∣x∣>R|x| > R. Na brzegu ∣x∣=R|x| = R trzeba sprawdzić osobno. Promień zwykle daje kryterium d'Alemberta, o ile granica istnieje (dla sin⁡x\sin x co drugi współczynnik jest zerem, więc porównujemy kolejne niezerowe wyrazy): ∣cn+1xn+1cnxn∣<1\left|\dfrac{c_{n+1}x^{n+1}}{c_nx^n}\right| < 1.

5.1 Szereg Taylora

W lekcji 1.4 Rachunek różniczkowy w praktyce funkcję przybliżaliśmy wielomianem Taylora z resztą RnR_n. Jeśli dla danego xx reszta Rn→0R_n \to 0 przy n→∞n \to \infty, wielomiany zbiegają do funkcji i dostajemy szereg Taylora:

f(x)=∑n=0∞f(n)(a)n!(x−a)n.f(x) = \sum_{n=0}^\infty \frac{f^{(n)}(a)}{n!}(x - a)^n .

Szeregi do zapamiętania (wokół x=0x = 0)

ex=∑n=0∞xnn!=1+x+x22+x36+…(R=∞)e^x = \sum_{n=0}^\infty\frac{x^n}{n!} = 1 + x + \frac{x^2}{2} + \frac{x^3}{6} + \dots\qquad (R = \infty) sin⁡x=x−x33!+x55!−…,cos⁡x=1−x22!+x44!−…(R=∞)\sin x = x - \frac{x^3}{3!} + \frac{x^5}{5!} - \dots,\qquad \cos x = 1 - \frac{x^2}{2!} + \frac{x^4}{4!} - \dots\qquad (R = \infty) 11−x=1+x+x2+…,ln⁡(1+x)=x−x22+x33−…(R=1)\frac{1}{1 - x} = 1 + x + x^2 + \dots,\qquad \ln(1 + x) = x - \frac{x^2}{2} + \frac{x^3}{3} - \dots\qquad (R = 1) (1+x)α=1+αx+α(α−1)2x2+…(R=1, gdy α∉N)(1 + x)^\alpha = 1 + \alpha x + \frac{\alpha(\alpha - 1)}{2}x^2 + \dots\qquad (R = 1,\ \text{gdy } \alpha \notin \mathbb N)

Wewnątrz przedziału zbieżności, czyli dla ∣x∣<R|x| < R, szeregi potęgowe można różniczkować i całkować wyraz po wyrazie. To potężne narzędzie. Całkując szereg geometryczny 11+t=1−t+t2−…\tfrac{1}{1 + t} = 1 - t + t^2 - \dots od 0 do xx, dostajesz szereg logarytmu. Całkując 11+t2=1−t2+t4−…\tfrac{1}{1 + t^2} = 1 - t^2 + t^4 - \dots, dostajesz arctan⁡x=x−x33+x55−…\arctan x = x - \tfrac{x^3}{3} + \tfrac{x^5}{5} - \dots, a dla x=1x = 1 wzór Leibniza π4=1−13+15−…\tfrac\pi4 = 1 - \tfrac13 + \tfrac15 - \dots Punkt x=1x = 1 leży na brzegu. Szereg jest tam zbieżny z kryterium Leibniza, a to, że jego suma to nadal arctan⁡1\arctan 1, wymaga osobnego argumentu (twierdzenia Abela). Tak samo dla ln⁡2\ln2 z szeregu logarytmu.

Promień zbieżności bywa zaskakujący. Funkcja 11+x2\tfrac{1}{1 + x^2} jest gładka na całej prostej, a jej szereg działa tylko dla ∣x∣<1|x| < 1. Przyczyna leży poza osią rzeczywistą: mianownik zeruje się w liczbach zespolonych x=±ix = \pm i, odległych od zera o 1 (lekcja 1.7 Liczby zespolone i fazory).

Symulacja 3 Sumy częściowe szeregu potęgowego (cienkie linie) i funkcja (gruba). Zwiększaj liczbę wyrazów NN. Wewnątrz przedziału zbieżności (zacieniony) sumy przylegają do funkcji coraz lepiej. Na zewnątrz rozjeżdżają się tym bardziej, im większe NN.

Przykład 5.1 Liczby ee i π\pi z szeregów

Oblicz ee z szeregu ∑1n!\sum\tfrac{1}{n!} z dokładnością 10−510^{-5} i oceń, jak szybko zbiega szereg Leibniza dla π\pi.

  1. e=1+1+12+16+124+…e = 1 + 1 + \tfrac12 + \tfrac16 + \tfrac1{24} + \dots Po 9 wyrazach (do 18!\tfrac1{8!}): 2,718 2792{,}718\,279. Ogon szacujemy bez znajomości ee: 19!(1+110+1100+… )=19!⋅109≈3,1⋅10−6\tfrac1{9!}\left(1 + \tfrac1{10} + \tfrac1{100} + \dots\right) = \tfrac{1}{9!}\cdot\tfrac{10}{9} \approx 3{,}1\cdot10^{-6}. Wyrazy maleją szybciej niż geometrycznie, więc zbieżność jest błyskawiczna.
  2. Leibniz: π=4(1−13+15−… )\pi = 4\left(1 - \tfrac13 + \tfrac15 - \dots\right). Błąd po NN wyrazach jest mniejszy od 42N+1\tfrac{4}{2N + 1}. Po 1000 wyrazach: 3,140 593{,}140\,59, błąd 0,001.
  3. Szybciej: wzór Machina π4=4arctan⁡15−arctan⁡1239\tfrac\pi4 = 4\arctan\tfrac15 - \arctan\tfrac1{239}. Tu xx jest małe, więc wyrazy maleją jak (1/5)2n(1/5)^{2n}. Pięć wyrazów każdego szeregu daje 3,141 592 683{,}141\,592\,68, czyli 8 poprawnych cyfr.

Wynik: e≈2,71828e \approx 2{,}71828 z 9 wyrazów; Leibniz potrzebuje tysiąca wyrazów na 3 cyfry, Machin pięciu na 8.

Sprawdzenie: e=2,718 281 8…e = 2{,}718\,281\,8\ldots, π=3,141 592 65…\pi = 3{,}141\,592\,65\ldots ✓. Morał: o szybkości decyduje nie liczba wyrazów, ale to, jak szybko maleją, czyli jak daleko od brzegu przedziału zbieżności jesteśmy.

Zadania · część 5

10 zadań · ok. 60 min

  1. Wyznacz promień zbieżności: ∑xnn\sum\dfrac{x^n}{n};   ∑n! xn\;\sum n!\,x^n;   ∑xnn!\;\sum\dfrac{x^n}{n!}.

    Odpowiedź

    R=1R = 1;   R=0\;R = 0 (zbieżny tylko dla x=0x = 0);   R=∞\;R = \infty.

  2. Oblicz sin⁡0,5\sin0{,}5 z dwóch i trzech pierwszych wyrazów szeregu. Porównaj z kalkulatorem.

    Odpowiedź

    0,5−0,020 83=0,479 170{,}5 - 0{,}020\,83 = 0{,}479\,17;   + 0,000 26=0,479 43\;+\,0{,}000\,26 = 0{,}479\,43. Kalkulator: 0,479 4260{,}479\,426.

  3. Napisz pierwsze cztery wyrazy szeregu e−x2e^{-x^2} (podstaw −x2-x^2 w miejsce xx w szeregu exe^x).

    Odpowiedź

    1−x2+x42−x661 - x^2 + \tfrac{x^4}{2} - \tfrac{x^6}{6}.

  4. W jakim przedziale zbieżny jest ∑(x−2)n3n\sum\dfrac{(x - 2)^n}{3^n}? Sprawdź też brzegi.

    Odpowiedź

    To szereg geometryczny z q=x−23q = \tfrac{x - 2}{3}: ∣x−2∣<3|x - 2| < 3, czyli x∈(−1; 5)x \in (-1;\ 5). Na brzegach q=±1q = \pm1, rozbieżny.

  5. Różniczkując wyraz po wyrazie szereg exe^x, sprawdź, że (ex)′=ex(e^x)' = e^x. Zrób to samo dla sin⁡x\sin x.

    Odpowiedź

    (1+x+x22+x36+… )′=0+1+x+x22+…(1 + x + \tfrac{x^2}{2} + \tfrac{x^3}{6} + \dots)' = 0 + 1 + x + \tfrac{x^2}{2} + \dots ✓.   (x−x36+x5120−… )′=1−x22+x424−⋯=cos⁡x\;(x - \tfrac{x^3}{6} + \tfrac{x^5}{120} - \dots)' = 1 - \tfrac{x^2}{2} + \tfrac{x^4}{24} - \dots = \cos x ✓.

  6. Oblicz 17=4(1+116)1/2\sqrt{17} = 4\left(1 + \tfrac{1}{16}\right)^{1/2} z dwóch i trzech wyrazów szeregu dwumiennego.

    Odpowiedź

    4⋅1,031 25=4,1254\cdot1{,}031\,25 = 4{,}125;   4(1+132−12048)=4,123 05\;4\left(1 + \tfrac1{32} - \tfrac1{2048}\right) = 4{,}123\,05. Dokładnie 4,123 114{,}123\,11.

  7. Oblicz sumę ∑n=1∞n2n\sum_{n=1}^\infty\dfrac{n}{2^n}, różniczkując szereg geometryczny ∑xn=11−x\sum x^n = \tfrac{1}{1 - x}.

    Wskazówka

    ∑nxn−1=1(1−x)2\sum nx^{n-1} = \tfrac{1}{(1 - x)^2}; pomnóż przez xx i podstaw x=12x = \tfrac12.

    Odpowiedź

    ∑nxn=x(1−x)2\sum nx^n = \tfrac{x}{(1 - x)^2}; dla x=12x = \tfrac12: 1/21/4=2\tfrac{1/2}{1/4} = 2.

  8. Feynman rozwija masę relatywistyczną m=m0(1−v2/c2)−1/2m = m_0(1 - v^2/c^2)^{-1/2} w szereg (lekcja 1.4 Rachunek różniczkowy w praktyce). Dla jakich prędkości ten szereg w ogóle jest zbieżny?

    Odpowiedź

    To szereg dwumienny (1+x)−1/2(1 + x)^{-1/2} z x=−v2/c2x = -v^2/c^2, o promieniu 1. Zbieżny dla v2<c2v^2 < c^2, czyli dla każdej prędkości mniejszej od prędkości światła. Szybko zbieżny tylko dla v≪cv \ll c.FLP I-15-8

  9. Pokaż, że funkcja 11+x2\tfrac{1}{1 + x^2} ma szereg o promieniu 1, choć jest określona i gładka dla każdego rzeczywistego xx. Co dzieje się z sumami częściowymi dla x=2x = 2?

    Odpowiedź

    Szereg geometryczny z q=−x2q = -x^2: zbieżny dla x2<1x^2 < 1. Dla x=2x = 2 wyrazy (−4)n(-4)^n rosną, a sumy częściowe skaczą coraz mocniej: 1,−3,13,−51,…1, -3, 13, -51, \dots, choć funkcja ma tam wartość 0,2.

  10. Całkując wyraz po wyrazie, znajdź szereg dla ∫0xe−t2 dt\displaystyle\int_0^{x}e^{-t^2}\,\mathrm dt i oblicz tę całkę dla x=1x = 1 z dokładnością do 0,001. (Funkcji pierwotnej e−t2e^{-t^2} nie da się wyrazić przez funkcje elementarne; to podstawa rozkładu normalnego, lekcja 5.1 Prawdopodobieństwo i rozkłady.)

    Odpowiedź

    x−x33+x510−x742+x9216−…x - \tfrac{x^3}{3} + \tfrac{x^5}{10} - \tfrac{x^7}{42} + \tfrac{x^9}{216} - \dots Dla x=1x = 1: 1−0,3333+0,1−0,0238+0,0046−0,0008≈0,7471 - 0{,}3333 + 0{,}1 - 0{,}0238 + 0{,}0046 - 0{,}0008 \approx 0{,}747 (dokładnie 0,746 820{,}746\,82). Szereg naprzemienny, błąd mniejszy od pierwszego pominiętego wyrazu.

6 Błędne intuicje

„Jeśli wyrazy dążą do zera, szereg jest zbieżny.”
To warunek konieczny, a nie wystarczający. Szereg harmoniczny ma an→0a_n \to 0 i jest rozbieżny.
„0,999…0{,}999\ldots jest trochę mniejsze od 1.”
To suma szeregu geometrycznego równa dokładnie 1. Między tymi liczbami nie ma żadnej innej.
„Nieskończona suma liczb dodatnich musi być nieskończona.”
12+14+18+⋯=1\tfrac12 + \tfrac14 + \tfrac18 + \dots = 1. Liczy się, jak szybko wyrazy maleją.
„Kolejność dodawania nie ma znaczenia.”
Dla szeregów warunkowo zbieżnych ma: przestawiając wyrazy, można dostać dowolną sumę (twierdzenie Riemanna).
„Funkcja gładka wszędzie ma szereg Taylora zbieżny wszędzie.”
11+x2\tfrac{1}{1 + x^2} jest gładka na całej prostej, a jej szereg działa tylko dla ∣x∣<1|x| < 1 (symulacja 3).
„Więcej wyrazów zawsze daje lepsze przybliżenie.”
Wewnątrz przedziału zbieżności błąd w końcu dąży do zera, ale nie musi maleć od pierwszego wyrazu: dla e−xe^{-x} przy x=5x = 5 sumy częściowe to 1, −4, 8,5, …1,\ -4,\ 8{,}5,\ \dots Poza przedziałem dodatkowe wyrazy nie pomagają, a sumy częściowe się rozjeżdżają.

7 Sprawdź się

1 Ile wynosi 1+12+14+18+…1 + \frac12 + \frac14 + \frac18 + \dots?

2.

2 Czy szereg harmoniczny jest zbieżny?

Nie, choć rośnie bardzo wolno (jak ln⁡N\ln N).

3 Dla jakich pp jest zbieżny ∑1np\sum\frac1{n^p}?

Dla p>1p > 1.

4 Co mówi kryterium d'Alemberta, gdy q=1q = 1?

Nic; trzeba użyć innego kryterium.

5 Jaki jest promień zbieżności szeregu exe^x?

Nieskończony.

6 Ile wynosi 1−12+13−14+…1 - \frac12 + \frac13 - \frac14 + \dots?

ln⁡2≈0,693\ln2 \approx 0{,}693.

7 Jak daleko może wystawać stos klocków poza stół?

Dowolnie daleko, bo szereg harmoniczny jest rozbieżny; ale wysięg rośnie tylko jak L2ln⁡n\tfrac L2\ln n.

8 Karta wzorów i co dalej

Szeregi liczbowe

∑an=lim⁡SN\sum a_n = \lim S_N,   SN=a1+⋯+aN\;S_N = a_1 + \dots + a_N

geometryczny: ∑n≥0aqn=a1−q\sum_{n\ge0}aq^n = \tfrac{a}{1 - q} dla ∣q∣<1|q| < 1

warunek konieczny: zbieżny ⇒an→0\Rightarrow a_n \to 0

harmoniczny rozbieżny, HN≈ln⁡N+0,577H_N \approx \ln N + 0{,}577

∑1np\sum\tfrac1{n^p}: zbieżny ⇔p>1\Leftrightarrow p > 1;   ∑1n2=π26\;\sum\tfrac1{n^2} = \tfrac{\pi^2}{6}

porównanie; całkowe; d'Alembert: q=lim⁡an+1an<1q = \lim\tfrac{a_{n+1}}{a_n} < 1 zbieżny, >1> 1 rozbieżny

Leibniz: bn↘0⇒∑(−1)nbnb_n \searrow 0 \Rightarrow \sum(-1)^nb_n zbieżny, błąd <bN+1< b_{N+1}

Szeregi potęgowe

promień RR: zbieżny dla ∣x∣<R|x| < R; brzegi osobno

ex=∑xnn!e^x = \sum\tfrac{x^n}{n!};   sin⁡x=∑(−1)nx2n+1(2n+1)!\;\sin x = \sum\tfrac{(-1)^nx^{2n+1}}{(2n + 1)!};   cos⁡x=∑(−1)nx2n(2n)!\;\cos x = \sum\tfrac{(-1)^nx^{2n}}{(2n)!}

11−x=∑xn\tfrac{1}{1 - x} = \sum x^n;   ln⁡(1+x)=∑n≥1(−1)n+1xnn\;\ln(1 + x) = \sum_{n\ge1}\tfrac{(-1)^{n+1}x^n}{n};   arctan⁡x=∑(−1)nx2n+12n+1\;\arctan x = \sum\tfrac{(-1)^nx^{2n+1}}{2n + 1} (R=1R = 1)

(1+x)α=1+αx+α(α−1)2x2+…(1 + x)^\alpha = 1 + \alpha x + \tfrac{\alpha(\alpha - 1)}{2}x^2 + \dots (R=1R = 1; dla α∈N\alpha \in \mathbb N to wielomian)

dla ∣x∣<R|x| < R: różniczkowanie i całkowanie wyraz po wyrazie

Do przeczytania u Feynmana

  • Obowiązkowo: 8-2 (Zenon: nieskończenie wiele kroków w skończonym czasie).
  • Dla chętnych: koniec 15-8 i początek 15-9 (rozwinięcie masy relatywistycznej i 12mv2\tfrac12mv^2).

Do przeczytania w podręcznikach matematyki

  • K. A. Stroud, Engineering Mathematics: rozdziały o szeregach, kryteriach zbieżności i szeregach potęgowych.
  • Dowolny podręcznik analizy dla pierwszego roku: dowody kryteriów i twierdzenie Riemanna o przestawianiu wyrazów.

Szereg exe^x z x=iθx = i\theta da wzór Eulera eiθ=cos⁡θ+isin⁡θe^{i\theta} = \cos\theta + i\sin\theta (1.7 Liczby zespolone i fazory). Szeregi funkcji zamiast potęg, czyli sumy sinusów i cosinusów, to szeregi Fouriera (4.2 Szeregi Fouriera). A szereg z silnią w mianowniku wróci przy rozkładzie Poissona (5.1 Prawdopodobieństwo i rozkłady).

Korzystają z tej lekcji 1.7 Liczby zespolone i fazory4.2 Szeregi Fouriera5.1 Prawdopodobieństwo i rozkładyAutomatyka 3.4 Sterowanie cyfrowe: próbkowanie, transformata Z i regulator w kodzie