Analiza matematyczna 2/Wykład 6: Ciągłość funkcji wielu zmiennych. Pochodne cząstkowe. Gradient

Funkcje wielu zmiennych. Ciągłość. Pochodne cząstkowe

Przypominamy przykłady funkcji wielu zmiennych, które znamy z życia codziennego. Do badania przebiegu zmienności funkcji, badania ciągłości, wyznaczania ekstremów stosujemy analizę przebiegu poziomic, a następnie wprowadzamy pochodne kierunkowe i cząstkowe.

Przykłady funkcji wielu zmiennych

Prognoza pogody

Prognoza pogody

Z funkcjami dwóch zmiennych rzeczywistych spotykamy się na co dzień. Śledząc prognozę pogody po wieczornym wydaniu wiadomości w telewizji (w prasie, w internecie) sprawdzamy, jaka temperatura jest przewidywana na najbliższą noc, na kolejny poranek, popołudnie, w dniach następnych. Temperatura podawana jest przeważnie liczbowo dla kilku regionów naszego kraju albo też - w dokładniejszej formie - na mapie z zaznaczonymi izotermami, tj. liniami, które łączą punkty o takiej samej temperaturze.

Osoby podatne na zmiany ciśnienia atmosferycznego z niepokojem śledzą informacje o spodziewanym załamaniu pogody i wahaniach ciśnienia. Przypomnijmy, że linie łączące punkty o takim samym ciśnieniu atmosferycznym nazywamy izobarami.

Zagęszczenie izobar nad danym obszarem oznacza dużą prędkość wiatru w terenie: im izobary są gęstsze, tym prędkość wiatru większa. Pamiętamy, że wiatr wieje od obszaru o wyższym ciśnieniu do obszaru o niższym ciśnieniu.

Kierunek wiatru także nie jest przypadkowy: odpowiada temu kierunkowi, w którym ciśnienie spada najszybciej, co na mapie odpowiada kierunkowi, w którym izobary najbardziej zagęszczają się.

Ze względu na czytelność map z prognozą pogody, obszary zawarte między kolejnymi poziomicami koloruje się zgodnie z umową tak, że obszary, nad którymi panuje niskie ciśnienie, bądź niska temparatura, oznacza się kolorem fioletowym, ciemno niebieskim, niebieskim. Kolory jasno zielony, zielony, jasno żółty, rezerwuje się do oznaczania obszarów o przeciętnym ciśnieniu czy temperaturze, natomiast obszary o najwyższych wartościach koloruje się na żółto, pomarańczowo, czerwono. Do umowy tej przywykliśmy. Tak bowiem pokolorowana jest mapa fizyczna (mapa hipsometryczna), np. ta przedstawiająca nasz kraj.

Prognoza pogody

Prognoza pogody

Mapa hipsometryczna

Gdybyśmy powędrowali palcem po mapie z południa na północ Polski, zaczynając od Tatr, które po polskiej stronie sięgają prawie 2500 metrów nad poziom morza, wystartowalibyśmy z obszaru pokolorowanego na brązowo, intensywnie czerwono, pomarańczowo. Kierując się do Krakowa i dalej Wyżyną Krakowsko-Częstochowską, przemierzalibyśmy obszar pokolorowany na żółto. Obszar nizinny w centralnej i północnej części naszego kraju zaznaczono na zielono, z wyjątkiem pasm wzgórz na północy, np. na Kaszubach, które zaznaczono na żółto. Jeśli spojrzymy trochę na prawo od ujścia Wisły, między Tczewem a Elblągiem, zauważymy obszar ciemnozielony, którym pokolorowano obszar depresji, tj. obszar położony poniżej poziomu morza. W końcu docieramy do brzegu Bałtyku, którego poziom stanowi umowny punkt odniesienia wysokości obszaru nad poziom morza. Pamiętamy, że głębokość dna morza na mapie również została zaznaczona różnymi kolorami: od białego (którym zaznaczono płytkie obszary tuż przy brzegu i mielizny), przez niebieski, aż po ciemnoniebieski, którym zaznaczono głębsze obszary.

Wycinek mapy Tatr

Wycinek mapy Tatr

Linie na mapie, łączące punkty o samej wysokości nad poziom morza, nazywamy poziomicami.

Wędrując po górach, w zależności od upodobania, wybieramy szlak, który krótszą, ale bardziej stromą drogą doprowadzi nas do celu, bądź też szlak mniej stromy, łagodny. Każdy, kto wędrował choć raz po górach z mapą w ręku wie, że im gęściej szlak poprzecinany jest kolejnymi poziomicami, tym jest bardziej stromy i wymaga większego wysiłku fizycznego. Szlak, który przebiega między dwiema poziomicami, prawie żadnej nie przecina, jest zdecydowanie łagodniejszy, bez stromych podejść, nie wymaga wysiłku.

Na ogół szlaki turystyczne w górach omijają obszary, gdzie poziomice przebiegają bardzo gęsto, bądź wręcz urywają się. Nic dziwnego: tak na mapie zaznaczono strome zbocza i urwiska.

Zauważmy, że poziomice odpowiadające różnym wysokościom są krzywymi rozłącznymi. Na mapie, która przeważnie przedstawia pewien prostokątny (w przybliżeniu) obszar terenu, krzywe te są zamknięte lub nie. Nasze doświadczenie podpowiada nam, że wewnątrz obszaru ograniczonego poziomicą, która jest linią zamkniętą, zawsze da się wskazać punkt położony najwyżej (np. szczyt wzniesienia) lub najniżej (np. dno doliny).

W ramach Analizy matematycznej I poznaliśmy twierdzenie, które opisuje taką sytuację: funkcja ciągła na zbiorze zwartym osiąga swoje kresy. Jest to twierdzenie Weierstrassa, które pozostaje prawdziwe nie tylko w przypadku funkcji jednej zmiennej.

Mapa fizyczna danego obszaru, mapa rozkładu ciśnienia, mapa rozkładu temperatury to przykłady graficznej reprezentacji (wykresu) funkcji dwóch zmiennych rzeczywistych (długości i szerokości geograficznej) o wartościach w zbiorze liczb rzeczywistych, bowiem wysokość punktu nad poziom morza, wartość ciśnienia atmosferycznego, temperatura to wielkości liczbowe.

Granica i ciągłość funkcji wielu zmiennych

Większość pojęć i twierdzeń, którymi będziemy posługiwać się w tym module, poznaliśmy już przy okazji omawiania własności funkcji ciągłych w przestrzeniach metrycznych. Przypomnijmy parę z nich.

Niech $(X,d)$, $(Y,\rho )$ będą przestrzeniami metrycznymi.

Będziemy zajmowali się badaniem funkcji

$f:X\mapsto Y$.

Większość praktycznych przykładów zastosowania teorii będzie dotyczyć funkcji określonych na zbiorze ${\mathbb{ℝ}}^n$, $n=2,3,\dots$, z metryką $d(x,y)=\Vert x-y\Vert$ zadaną przez pewną ustaloną normę $\Vert \cdot \Vert$ w ${\mathbb{ℝ}}^n$, np.

$\begin{array}{rl}\Vert x{\Vert }_p& =(|x_1{|}^p+|x_2{|}^p+\dots +|x_n{|}^p{)}^{\frac{1}{p}},\text{ dla }1\le p<\mathrm{\infty }\\ & \text{ w szczególności }\\ \Vert x{\Vert }_1& =|x_1|+|x_2|+\dots +|x_n|\\ \Vert x{\Vert }_2& =\sqrt{|x_1{|}^2+|x_2{|}^2+\dots +|x_n{|}^2}\\ & \text{ bądź też }\\ \Vert x{\Vert }_{\mathrm{\infty }}& =\max \{|x_1|,|x_2|,\dots ,|x_n|\}.\end{array}$

Zbiorem wartości funkcji $f$ najczęściej będzie zbiór liczb rzeczywistych $\mathbb{ℝ}$ z metryką zadaną przez wartość bezwzględną, tj. $\rho (a,b)=|a-b|$.

Definicja 6.1.

Mówimy, że $g\in Y$ jest granicą funkcji $f:X\mapsto Y$ w punkcie $x$ będącym punktem skupienia dziedziny funkcji $f$, jeśli

$\mathrm{\forall }ϵ>0\mathrm{\exists }\delta >0:\mathrm{\forall }y:0<d(x,y)<\delta ⟹\rho (g,f(y))<ϵ$.

Definicja 6.2.

Mówimy, że funkcja $f:X\mapsto Y$ jest ciągła w punkcie x, jeśli

$\mathrm{\forall }ϵ>0\mathrm{\exists }\delta >0:\mathrm{\forall }y:d(x,y)<\delta ⟹\rho (f(x),f(y))<ϵ$.

Pamiętamy również, że zachodzi następujące

Twierdzenie 6.3.

Niech $X,Y$ będą przestrzeniami metrycznymi i niech $f:X\mapsto Y$ będzie funkcją. Wówczas następujące warunki są równoważne:

1) funkcja $f$ jest ciągła w punkcie $a\in X$,

2) istnieje granica $\underset{x\to a}{\lim }f(x)$ i jest równa wartości funkcji $f(a)$.

Niech $X$, $Y$, $Z$ będą przestrzeniami metrycznymi.

Twierdzenie 6.4.

Złożenie $g\circ f:X\mapsto Z$ funkcji ciągłych $f:X\mapsto Y$ i $g:Y\mapsto Z$ jest funkcją ciągłą.

Twierdzenie 6.5.

Jeśli $f:X\mapsto \mathbb{ℝ}$ oraz $g:X\mapsto \mathbb{ℝ}$ są funkcjami ciągłymi, to suma $f+g$ oraz iloczyn $f\cdot g$ są funkcjami ciągłymi. Ponadto odwrotność $\frac{1}{g}:Z\ni x\mapsto \frac{1}{g(x)}\in \mathbb{ℝ}$ oraz iloraz $\frac{f}{g}:Z\ni x\mapsto \frac{f(x)}{g(x)}\in \mathbb{ℝ}$

są funkcjami ciągłymi na zbiorze $Z:=X\setminus \{x\in X:g(x)=0\}$.

Przypomnijmy jeszcze następujące twierdzenie Weierstrassa o osiąganiu kresów przez funkcję ciągłą na zbiorze zwartym.

Twierdzenie 6.6.

Jeśli $f:X\mapsto \mathbb{ℝ}$ jest funkcją ciągłą określoną na przestrzeni zwartej $X$, to istnieją punkty $a,b\in X$, w których funkcja $f$ osiąga kresy: kres dolny $\inf \{f(x),x\in X\}=f(a)$ i kres górny $\sup \{f(x),x\in X\}=f(b)$.

Rozważmy przykład, który pokazuje, że zachowanie funkcji wielu zmiennych może wykraczać poza naszą intuicję ukształtowaną podczas badania funkcji jednej zmiennej.

Przykład 6.7.

Funkcja $f(x,y)=\frac{xy}{x^2+y^2}$ określona jest we wszystkich punktach płaszczyzny ${\mathbb{ℝ}}^2$ z wyjątkiem punktu $(0,0)$. Wyraźmy ją we współrzędnych biegunowych

$\mathrm{\Phi }:(r,\phi )\mapsto \{\begin{array}{r}x(r,\phi )=r\cos \phi \\ y(r,\phi )=r\sin \phi \end{array}$

W punktach leżących poza początkiem układu współrzędnych, tj. gdy $r>0$, otrzymamy:

$(f\circ \mathrm{\Phi })(r,\phi )=\frac{r^2\cos \phi \sin \phi }{r^2({\cos }^2\phi +{\sin }^2\phi )}=\frac{1}{2}\sin 2\phi$.

Zauważmy, że zbiorem wartości tej funkcji jest przedział $[-\frac{1}{2},\frac{1}{2}]$. Ponadto funkcja $(r,\phi )\mapsto \frac{1}{2}\sin 2\phi$ nie zależy od zmiennej $r$. Oznacza to, że zacieśnienie funkcji $f$ do którejkolwiek półprostej danej równaniem $\phi ={\phi }_0$ (tj. półprostej, która tworzy z dodatnią półosią osi rzędnych ustalony kąt ${\phi }_0$) jest funkcją o stałej wartości $\frac{1}{2}\sin 2{\phi }_0$, niezależnej od odległości $r$ punktu od początku układu współrzędnych. Każde z tych zacieśnień do prostej $\{\phi ={\phi }_0\}\cup \{\phi ={\phi }_0+\pi \}$ ma granicę przy $r\to 0$ równą $\frac{1}{2}\sin 2{\phi }_0$. Jednak wartość ta zależy od wyboru kąta ${\phi }_0$, stąd nie istnieje granica funkcji $(x,y)\mapsto f(x,y)$, gdy $(x,y)\to (0,0)$. Zauważmy, że gdy rozważymy osobno funkcje jednej zmiennej, ustalając drugą zmienną, tzn. $y$ lub odpowiednio $x$:

$\begin{array}{rl}{f}_y& =f(\cdot ,y):\mathbb{ℝ}\ni x\mapsto \frac{xy}{x^2+y^2}\\ f_x& =(x,\cdot ):\mathbb{ℝ}\ni y\mapsto \frac{xy}{x^2+y^2},\end{array}$

to zarówno $\underset{x\to 0}{\lim }{f}_y(x)=0$, jak też $\underset{y\to 0}{\lim }{f}_x(y)=0$, a więc w szczególności istnieją granice iterowane

$\begin{array}{r}\underset{x\to 0}{\lim }(\underset{y\to 0}{\lim }f(x,y))=0,\\ \underset{y\to 0}{\lim }(\underset{x\to 0}{\lim }f(x,y))=0\end{array}$ i są równe.

<applet code="JavaviewModApplet.class" archive="images/a/a0/Javaview.jar,images/b/be/JavaviewModApplet.jar" width="450" height="400">

   <param name="colors" value="-1:#7e0b9f -0.80:#ce3bf8 -0.60:#2c4ae5 -0.40:#2c85e5 -0.20:#2ecca9 -0.05:#2ecc5b 0.05:#2ecc5b 0.20:#97cc2e 0.40:#edff27 0.60:#ffba27 0.80:#ff6e27 1:#d42525">
   <param name="coloring" value="maple">
   <param name="model" value="images/7/74/Am2m05.0080.mgs.zip">
  <param name="scale" value="1.0 1.0 1.0">

<param name="shading" value="0.2"> </applet>

<div.thumbcaption>Wykres funkcji $f(x,y)=\frac{xy}{x^2+y^2}$

Przykład pokazuje więc, że

Wniosek 6.8.

Z istnienia granic iterowanych

$\begin{array}{r}\underset{x\to a}{\lim }(\underset{y\to b}{\lim }f(x,y))\\ \underset{y\to b}{\lim }(\underset{x\to a}{\lim }f(x,y))\end{array}$

i równości tych granic nie wynika istnienie granicy funkcji $f$ w punkcie $(a,b)$.

Prawdziwa natomiast jest implikacja:

Uwaga 6.9.

Jeśli funkcja $f:\mathbb{ℝ}\times \mathbb{ℝ}\mapsto \mathbb{ℝ}$ ma granicę w punkcie $(a,b)$, to istnieją obie granice iterowane

$\begin{array}{r}\underset{x\to a}{\lim }(\underset{y\to b}{\lim }f(x,y))\\ \underset{y\to b}{\lim }(\underset{x\to a}{\lim }f(x,y))\end{array}$

i są równe granicy funkcji $f$ w punkcie $(a,b)$.

Uwaga ta stanowi warunek konieczny istnienia granicy $\underset{(x,y)\to (a,b)}{\lim }f(x,y)$. Jeśli bowiem nie istnieje któraś z granic iterowanych, bądź nie są one równe, to funkcja $f$ nie ma granicy w punkcie $(a,b)$. Podkreślmy jeszcze raz fakt, że istnienie i równość obu granic iterowanych nie gwarantuje istnienia granicy funkcji.

Poziomice

Niech $f:X\mapsto \mathbb{ℝ}$ będzie funkcją określoną na przestrzeni metrycznej $X$ o wartościach rzeczywistych.

Definicja 6.10.

Poziomicą funkcji $f$

odpowiadającą wartości $a\in \mathbb{ℝ}$ nazywamy zbiór $\{f=a\}=\{x\in X:f(x)=a\}$, czyli przeciwobraz zbioru jednopunktowego $\{a\}$ przez funkcję $f$.

Często pobieżna nawet analiza przebiegu poziomic pozwala ustalić, czy dana funkcja osiąga ekstrema (zob. definicja ekstremum w module 9 Analizy matematycznej 1), czy też nie.

Przykład 6.11.

Niech $f(x,y)=x^2+y^2-4$.
wykres
Poziomica $\{f=a\}=\{(x,y):x^2+y^2-4=a\}$ jest okręgiem o środku w punkcie $(0,0)$ i promieniu $\sqrt{4+a}$, gdy $a>-4$. Poziomica $\{f=-4\}$ składa się tylko z jednego punktu $(0,0)$, natomiast jeśli $a<-4$, to poziomica $\{f=a\}$ jest zbiorem pustym. Funkcja $f$ osiąga minimum globalne w punkcie $(0,0)$ równe $f(0,0)=-4$.

Przykład 6.12.

Niech $f(x,y)=x^2-y^2$.
wykres

Poziomica zerowa $\{f=0\}=\{(x,y):x^2-y^2=0\}=\{x=y\}\cup \{x=-y\}$ jest sumą dwóch prostych: $x=y$ i $x=-y$. Jeśli $a\ne 0$ poziomica $\{f=a\}=\{x^2-y^2=a\}$ jest hiperbolą o asymptotach $x=y$ i $x=-y$. Przebieg poziomic pozwala ustalić, że funkcja $f$ w żadnym punkcie płaszczyzny nie osiąga ekstremum, bowiem w dowolnie małym otoczeniu każdego punktu $(x,y)$ potrafimy z łatwością wskazać punkty, w których funkcja przyjmuje zarówno wartości mniejsze jak i większe od wartości funkcji $f$ w punkcie $(x,y)$.

Przykład 6.13.

Niech $f(x,y)=|x|+|y|$.
wykres
Funkcja $f$ jest normą w ${\mathbb{ℝ}}^2$, przyjmuje więc wyłącznie wartości nieujemne, stąd $\{f=a\}=\mathrm{\varnothing }$, gdy $a<0$. Poziomica zerowa $\{f=0\}=\{(0,0)\}$ składa się tylko z jednego punktu. Gdy $a>0$, poziomica $\{f=a\}=\{|x|+|y|=a\}$ jest kwadratem o wierzchołkach $(a,0)$, $(0,a)$, $(-a,0)$, $(0,-a)$. Funkcja $f$ osiąga minimum globalne w punkcie $(0,0)$, gdyż $f(x,y)>0$ w dowolnym punkcie $(x,y)\ne (0,0)$. Podobnie jak w poprzednim przykładzie przebieg poziomic pozwala ustalić, że funkcja $f$ w żadnym punkcie płaszczyzny poza punktem $(0,0)$ nie osiąga ekstremum, bowiem w dowolnie małym otoczeniu każdego punktu $(x,y)\ne (0,0)$ potrafimy z łatwością wskazać punkty, w których funkcja przyjmuje zarówno wartości mniejsze jak i większe od wartości funkcji $f$ w punkcie $(x,y)$.

Przykład 6.14.

Niech $f(x,y)=|x{|}^{\frac{2}{3}}+|y{|}^{\frac{2}{3}}$.
wykres
Funkcja $f$ przyjmuje wyłącznie wartości nieujemne, stąd $\{f=a\}=\mathrm{\varnothing }$, gdy $a<0$. Poziomica zerowa $\{f=0\}=\{(0,0)\}$ składa się tylko z jednego punktu. Gdy $a>0$, poziomica $\{f=a\}=\{|x{|}^{\frac{2}{3}}+|y{|}^{\frac{2}{3}}=a\}$ jest krzywą zawartą w kwadracie o wierzchołkach $(\sqrt{a^3},0)$, $(0,\sqrt{a^3})$, $(-\sqrt{a^3},0)$, $(0,-\sqrt{a^3})$. Krzywą tę nazywamy asteroidą. Zauważmy, że funkcja osiąga minimum globalne w punkcie $(0,0)$, gdyż $f(x,y)>0$, w dowolnym punkcie $(x,y)\ne (0,0)$. Podobnie jak w poprzednich przykładach przebieg poziomic pozwala stwierdzić, że jest to jedyne ekstremum tej funkcji na płaszczyźnie ${\mathbb{ℝ}}^2$.

Przykład 6.15.

Niech $f(x,y)=xy(1-x-y)$.
Poziomicą zerową $\{f=0\}$ tej funkcji jest suma trzech prostych: $x=0$, $y=0$ oraz $x+y=1$. Zauważmy, że w dowolnie małym otoczeniu któregokolwiek punktu leżącego na tej poziomicy znajdziemy punkty, w których funkcja osiąga zarówno wartości mniejsze jak i większe od zera. Stąd w żadnym punkcie zbioru $\{f=0\}$ funkcja $f$ nie osiąga ekstremum. Zwróćmy jednak uwagę, że wnętrze trójkąta o wierzchołkach $(0,0)$, $(1,0)$, $(0,1)$ zawarte jest w zbiorze $\{f<0\}$ tych punktów, w których funkcja przyjmuje wartości ujemne. W sumie mnogościowej z brzegiem trójkąta o podanych wierzchołkach zbiór ten jest zwarty. Z twierdzenia Weierstrassa o osiąganiu kresów przez funkcję ciągłą na zbiorze zwartym wynika, że we wnętrzu tego trójkąta funkcja $f$ osiąga minimum. Dalsza analiza przebiegu poziomic we wnętrzu trójkąta nie jest efektywnym narzędziem prowadzącym do precyzyjne określenia, w którym punkcie minimum to jest osiągane, ponieważ kreślenie poziomic $\{f=a\}$, gdy $a\ne 0$, nie jest prostym zadaniem. Warto jednak zauważyć, że żaden punkt leżący na poziomicy zerowej funkcji $f$ nie może być jej punktem ekstremalnym, gdyż w otoczeniu któregokolwiek punktu z tej poziomicy funkcja przyjmuje zarówno wartości dodatnie jak i ujemne.

<applet code="JavaviewModApplet.class" archive="images/a/a0/Javaview.jar,images/b/be/JavaviewModApplet.jar" width="450" height="400">

  <param name="colors" value="-1:#7e0b9f -0.80:#ce3bf8 -0.60:#2c4ae5 -0.40:#2c85e5 -0.20:#2ecca9 -0.05:#2ecc5b 0.05:#2ecc5b 0.20:#97cc2e 0.40:#edff27 0.60:#ffba27 0.80:#ff6e27 1:#d42525">
   <param name="coloring" value="maple">
   <param name="model" value="images/6/65/Am2m05.0130.mgs.zip">
   <param name="scale" value="1.0 1.0 0.1">

<param name="shading" value="0.2"> </applet>

<div.thumbcaption>Wykres funkcji $f(x,y)=xy(1-x-y)$

Kartezjusz (1596-1650)
Zobacz biografię

Przykład 6.16.

Niech $f(x,y)=x^3+y^3-3xy$.
wykres

Poziomicą zerową $\{f=0\}$ tej funkcji jest nieograniczona krzywa, którą nazywamy liściem Kartezjusza. Krzywa ta ma asymptotę o równaniu $x+y+1=0$. W pierwszej ćwiartce $\{(x,y)\in {\mathbb{ℝ}}^2:x\ge 0,y\ge 0\}$ tworzy petlę ograniczającą obszar, we wnętrzu którego funkcja $f$ przyjmuje wartości ujemne. Podobnie jak w poprzednim przykładzie, z twierdzenia Weierstrassa o osiąganiu kresów przez funkcję ciągłą na zbiorze zwartym wynika, że we wnętrzu pętli liścia Kartezjusza funkcja $f$ osiąga minimum. Dalsza analiza przebiegu poziomic nie prowadzi efektywnie do precyzyjne określenia, w którym punkcie minimum to jest osiągane. Warto jednak zauważyć, że w żadnym z punktów liścia Kartezjusza funkcja $f$ nie może osiągać ekstremum, gdyż w dowolnie małym otoczeniu któregokolwiek punktu tej krzywej funkcja $f$ osiąga wartości dodatnie jak i ujemne.

Przykład 6.17.

Niech $f(x,y)=(x^2+y^2{)}^2-2(x^2-y^2)$.
wykres

Jakob Bernoulli (1654-1705)
Zobacz biografię

Poziomicą zerową $\{f=0\}$ tej funkcji jest krzywa, zwana lemniskatą Bernoullego. Przebieg lemniskaty $\{f=0\}$ najwygodniej zbadać we współrzędnych biegunowych:

$\begin{array}{rl}(x^2+y^2{)}^2& =2(x^2-y^2)\\ (r^2{\cos }^2\phi +r^2{\sin }^2\phi {)}^2& =2(r^2{\cos }^2\phi -r^2{\sin }^2\phi )\\ r^4& =2r^2\cos 2\phi \\ r=0& \text{ lub }r=\sqrt{2\cos 2\phi },\end{array}$

przy czym wyrażenie pod znakiem pierwiastka jest nieujemne dla $\phi \in [-\frac{\pi }{4},\frac{\pi }{4}]\cup [\frac{3\pi }{4},\frac{5\pi }{4}]$. Lemniskata Bernoullego jest więc zawarta w części wspólnej koła o promieniu $\sqrt{2}$ i dwóch obszarów wyciętych z płaszczyzny półprostymi tworzącymi z osią rzędnych kąty $-\frac{\pi }{4},\frac{\pi }{4},\frac{3\pi }{4},-\frac{3\pi }{4}$. Zauważmy, że we wnętrzu obszaru ograniczonego lemniskatą Bernoullego funkcja $f$ osiąga wartości ujemne. Na zewnątrz zaś - dodatnie. Podobnie jak w obu poprzednich przykładach twierdzenie Weierstrassa gwarantuje istnienie minimum lokalnego w obszarze $\{(x,y):f(x,y)\le 0\}$ ograniczonym lemniskatą Bernoullego. Z kolei w dowolnie małym otoczeniu każdego punktu poziomicy zerowej $\{f=0\}$ funkcja przyjmuje zarówno wartości dodatnie jak i ujemne. Funkcja $f$ nie osiąga więc ekstremum w żadnym punkcie należącym do swojej poziomicy zerowej.

Trzy ostatnie przykłady każą nam szukać doskonalszego narzędzia do precyzyjnego lokalizownia punktów ekstremalnych funkcji wielu zmiennych, gdy analiza przebiegu poziomic jest niewystarczająca. Tym narzędziem są

Pochodna kierunkowa i pochodne cząstkowe

Niech $A\subset X$ będzie otwartym podzbiorem przestrzeni unormowanej $X$. Niech $v\ne 0,v\in X$ będzie ustalonym niezerowym wektorem tej przestrzeni.

Definicja 6.18.

Mówimy, że funkcja $f:A\mapsto \mathbb{ℝ}$ ma w punkcie $a$ pochodną kierunkową w kierunku wektora $v$, jeśli

istnieje granica ilorazu różnicowego: $\underset{h\to 0}{\lim }\frac{f(a+hv)-f(a)}{h}$.

Granicę tę oznaczamy symbolem ${\partial }_{v}f(a)$ i nazywamy pochodną kierunkową funkcji $f$ w kierunku wektora $v$ w punkcie $a$.

Zwróćmy uwagę, że zbiór $\{a+tv,t\in \mathbb{ℝ}\}$ jest prostą przechodzącą przez punkt $a$ równoległą do wektora $v$. Stąd pochodna ${\partial }_{v}f(a)$ jest w istocie pochodną w punkcie $t=0$ funkcji jednej zmiennej rzeczywistej $t\mapsto f(a+tv)$, czyli restrykcji funkcji $f$ do podzbioru otwartego $A\cap \{a+tv,t\in \mathbb{ℝ}\}$ rozważanej prostej $\{a+tv,t\in \mathbb{ℝ}\}$. Wobec tego możemy powtórzyć jednowymiarowy warunek konieczny istnienia ekstremum w punkcie, w którym istnieje pochodna kierunkowa funkcji (zob. moduł 9, Analiza matematyczna I).

Twierdzenie 6.19.

Niech $A\subset X$ będzie otwartym podzbiorem przestrzeni unormowanej $X$ i niech $v\in X$, $v\ne 0$. Jeśli funkcja $f:A\mapsto \mathbb{ℝ}$ osiąga ekstremum w punkcie $a\in A$ i istnieje pochodna kierunkowa ${\partial }_{v}f(a)$, to

pochodna ta zeruje się.

Dowód 6.19.

Jeśli funkcja $A\ni x\mapsto f(x)\in \mathbb{ℝ}$ osiąga maksimum (odpowiednio: minimum) w punkcie $a$, to funkcja jednej zmiennej $t\mapsto f(a+tv)$ osiąga maksimum (odpowiednio: minimum) w punkcie $t=0$. Z warunku koniecznego istnienia ekstremum funkcji jednej zmiennej wynika, że pochodna (o ile istnieje) funkcji

$t\mapsto f(a+tv)$ zeruje się w punkcie $t=0$. Stąd ${\partial }_{v}f(a)=0$

O ile w przypadku funkcji określonej na otwartym przedziale prostej $\mathbb{ℝ}$ sytuacja jest oczywista (na prostej mamy tylko jeden kierunek), o tyle już w przypadku funkcji dwóch zmiennych (na płaszczyźnie można wskazać nieskończenie wiele kierunków!) powstaje pytanie o liczbę kierunków, które należy ustalić, aby rozwiązać praktyczny problem lokalizacji punktów ekstremalnych danej funkcji. Warto zauważyć, że w przypadku funkcji określonej na $n$ wymiarowej przestrzeni unormowanej $X$ nie ma potrzeby rozważać pochodnych kierunkowych w kierunku wektorów liniowo zależnych, a więc większej niż wynosi wymiar przestrzeni. Wyróżnijmy wobec tego pochodne kierunkowe w kierunku wektorów bazowych.

Niech $X={\mathbb{ℝ}}^n$ i niech $e_1=(1,0,0,\dots ,0)$, $e_2=(0,1,0,\dots ,0)$, ..., $e_n=(0,0,0,\dots ,1)$ będzie bazą kanoniczną tej przestrzeni. Niech $A$ będzie otwartym podzbiorem przestrzeni ${\mathbb{ℝ}}^n$.

Definicja 6.20.

Pochodne kierunkowe (o ile istnieją) ${\partial }_{e_1}f(a)$, ${\partial }_{e_2}f(a)$, ..., ${\partial }_{e_n}f(a)$ funkcji $f:A\mapsto \mathbb{ℝ}$ w kierunku wektorów bazy $\{e_1,e_2,\dots ,e_n\}$ nazywamy pochodnymi cząstkowymi funkcji $f$ w punkcie $a$. Pochodną cząstkową funkcji $(x_1,x_2,\dots ,x_n)\mapsto f(x_1,x_2,\dots ,x_n)\in \mathbb{ℝ}$ w kierunku wektora

$e_i$ oznaczamy tradycyjnie symbolem: $\frac{\partial f}{\partial x_i}(a),\frac{\partial }{\partial x_i}f(a),f_{x_i}(a)\text{ lub }f{\text{'}}_{x_i}(a)$.

W przypadku, gdy nie numerujemy współrzędnych argumentu funkcji $(x,y,z)\mapsto f(x,y,z)$ pochodne cząstkowe oznaczamy symbolami

$\frac{\partial f}{\partial x}(a),\frac{\partial f}{\partial y}(a),\frac{\partial f}{\partial z}(a)$.

Przeformułujmy warunek konieczny istnienia ekstremum funkcji określonej na zbiorze otwartym $A\subset {\mathbb{ℝ}}^n$.

Twierdzenie 6.21.

Jeśli funkcja $f:A\mapsto \mathbb{ℝ}$ osiąga ekstremum w punkcie $a\in A$, w którym istnieją pochodne cząstkowe $\frac{\partial }{\partial x_k}f(a)$, $k\in \{1,2,\dots ,n\}$, to pochodne te zerują się w tym punkcie, tj.

$\mathrm{\forall }k\in \{1,2,\dots ,n\}:\frac{\partial }{\partial x_k}f(a)=0$.

Zwróćmy uwagę, że twierdzenie podaje jedynie warunek konieczny istnienia ekstremum. Punkt $a$, który spełnia układ równań:

$\{\begin{array}{rl}\frac{\partial f}{\partial x_1}(a)& =0\\ \frac{\partial f}{\partial x_2}(a)& =0\\ & ⋮\\ \frac{\partial f}{\partial x_n}(a)& =0\end{array}$

nie musi być punktem ekstremalnym funkcji $f$.

Wróćmy do przykładów, w których stwierdziliśmy potrzebę znalezienia dokładniejszego narzędzia do lokalizacji ekstremów.

Przykład 6.22.

Z przebiegu poziomicy zerowej funkcji $f(x,y)=xy(1-x-y)$ wywnioskowaliśmy - w oparciu o twierdzenie Weierstrassa o osiąganiu kresów przez funkcję ciągłą na zbiorze zwartym - że funkcja ta osiąga minimum w pewnym punkcie wewnątrz trójkąta o wierzchołkach $(0,0)$, $(1,0)$, $(0,1)$. Rozwiązując układ dwóch równań

$\{\begin{array}{r}\frac{\partial f}{\partial x}=0\\ \frac{\partial f}{\partial y}=0\end{array}⟺\{\begin{array}{r}y-2xy-y^2=0\\ x-x^2-2xy=0\end{array}$.

otrzymujemy układ

$\{\begin{array}{r}y=0\text{ lub }1-2x-y=0\\ x=0\text{ lub }1-x-2y=0\end{array}$,

który spełniają współrzędne czterech punktów $P_1=(0,0)$, $P_2=(1,0)$, $P_3=(0,1)$, $P_4=(\frac{1}{3},\frac{1}{3})$. Jedynym punktem z wnętrza wskazanego trójkąta jest punkt $P_4$, w którym funkcja $f$ osiąga minimum równe $f(P_4)=\frac{1}{27}$. Pozostałe punkty $P_1$, $P_2$, $P_3$ leżą na poziomicy zerowej funkcji $f$, która - jak już sprawdziliśmy - nie może zawierać żadnego punktu ekstremalnego funkcji $f$ (zob. przykład 6.15.).

Przykład 6.23.

Z przebiegu poziomicy zerowej funkcji $f(x,y)=x^3+y^3-3xy$ wywnioskowaliśmy - w oparciu o twierdzenie Weierstrassa o osiąganiu kresów przez funkcję ciągłą na zbiorze zwartym - że funkcja ta osiąga minimum w pewnym punkcie wewnątrz pętli liścia Kartezjusza. Rozwiązując układ dwóch równań

$\{\begin{array}{r}\frac{\partial f}{\partial x}=0\\ \frac{\partial f}{\partial y}=0\end{array}⟺\{\begin{array}{r}3x^2-3y=0\\ 3y^2-3x=0\end{array}$.

otrzymujemy układ

$\{\begin{array}{rl}y=0& \text{ lub }y=1\\ x& =y^2\end{array}$,

który spełniają współrzędne dwóch punktów $P_1=(0,0)$, $P_2=(1,1)$. Jedynym punktem z wnętrza obszaru ograniczonego przez pętlę liścia Kartezjusza jest punkt $P_2$, w którym funkcja $f$ osiąga minimum równe $f(P_2)=-1$. Punkt $P_1$ leży na poziomicy zerowej funkcji $f$, która - jak już sprawdziliśmy - nie może zawierać żadnego punktu ekstremalnego funkcji $f$ (zob. przykład 6.16.).

Przykład 6.24.

Podobnie jak w obu poprzednich przykładach z przebiegu poziomicy zerowej funkcji $f(x,y)=(x^2+y^2{)}^2-2(x^2-y^2)$ wywnioskowaliśmy - w oparciu o twierdzenie Weierstrassa o osiąganiu kresów przez funkcję ciągłą na zbiorze zwartym - że funkcja ta osiąga minimum w pewnym punkcie wewnątrz obszaru ograniczonego lemniskatą Bernoullego. Rozwiązując układ dwóch równań

$\{\begin{array}{r}\frac{\partial f}{\partial x}=0\\ \frac{\partial f}{\partial y}=0\end{array}⟺\{\begin{array}{r}2(x^2+y^2)2x-4x=0\\ 2(x^2+y^2)2y+4y=0\end{array}$.

otrzymujemy układ

$\{\begin{array}{r}x=0\text{ lub }{x}^2+y^2-1=0\\ y=0\text{ lub }{x}^2+y^2+1=0\end{array}$,

który spełniają współrzędne trzech punktów $P_1=(0,0)$, $P_2=(-1,0)$, $P_3=(1,0)$. We wnętrzu obszaru ograniczonego lemniskatą Bernoullego leżą punkty $P_2$ i $P_3$, w których funkcja $f$ osiąga minima równe $f(P_2)=f(P_3)=-1$. Punkt $P_1$ leży na poziomicy zerowej funkcji $f$, która - jak już sprawdziliśmy - nie może zawierać żadnego punktu ekstremalnego funkcji $f$ (zobacz przykład 6.17.).

Pochodne cząstkowe wyższych rzędów

Rozważmy funkcję $\frac{\partial f}{\partial x_i}$, która punktowi $x\in U$ przyporządkowuje pochodną cząstkową funkcji $f$ po zmiennej $x_i$ w punkcie $a$, czyli funkcję

$\frac{\partial f}{\partial x_i}:U\ni a\mapsto \frac{\partial f}{\partial x_i}(a)\in \mathbb{ℝ}$.

Definicja 6.25.

Jeśli w punkcie $a\in U$ istnieje pochodna cząstkowa funkcji $\frac{\partial f}{\partial x_i}$ po zmiennej $x_j$, to mówimy, że funkcja $f$ ma pochodną cząstkową rzędu drugiego po zmiennych $x_i$ oraz $x_j$. Pochodną tę oznaczamy symbolem $\frac{\partial }{\partial x_j}\frac{\partial }{\partial x_i}f(a)$, bądź krótko $\frac{{\partial }^2}{\partial x_j\partial x_i}f(a)$ lub $\frac{{\partial }^{2}f(a)}{\partial x_j\partial x_i}$. Gdy $i=j$

piszemy $\frac{{\partial }^{2}f(a)}{\partial x_i^2}$ zamiast $\frac{{\partial }^{2}f(a)}{\partial x_i\partial x_i}$.

Uwaga 6.26.

Jeśli $f:{\mathbb{ℝ}}^n\ni (x,y,z,\dots ,t)\mapsto f(x,y,z,\dots ,t)\in \mathbb{ℝ}$ jest funkcją $n$ zmiennych, to często zamiast pisać

$\frac{{\partial }^{2}f(a)}{\partial x^2},\frac{{\partial }^{2}f(a)}{\partial x\partial y},\frac{{\partial }^{2}f(a)}{\partial x\partial z},\dots$

piszemy

$f_{xx}(a),f_{xy}(a),f_{xz}(a),\dots$,

bądź

$f{\text{'}}_{xx}(a),f{\text{'}}_{xy}(a),f{\text{'}}_{xz}(a),\dots$

Powstaje naturalne pytanie, czy zachodzi równość między pochodnymi $\frac{\partial }{\partial x_j}\frac{\partial }{\partial x_i}f(a)$ oraz $\frac{\partial }{\partial x_i}\frac{\partial }{\partial x_j}f(a)$, jeśli obie istnieją.

Zanim sformułujemy twierdzenie, które stanowi pozytywną odpowiedź na pytanie, rozważmy następujący

Przykład 6.27.

Funkcja $f(x,y)=\{\begin{array}{rlr}& \frac{xy(x^2-y^2)}{x^2+y^2},& \text{ gdy }(x,y)\ne (0,0)\\ & 0,& \text{ gdy }(x,y)=(0,0)\end{array}$. ma w punkcie $(0,0)$ obie pochodne cząstkowe mieszane $\frac{\partial }{\partial x}\frac{\partial }{\partial y}f(0,0)$ oraz $\frac{\partial }{\partial y}\frac{\partial }{\partial x}f(0,0)$, lecz są one różne. A mianowice $\frac{\partial }{\partial x}\frac{\partial }{\partial y}f(0,0)=1$, podczas gdy $\frac{\partial }{\partial y}\frac{\partial }{\partial x}f(0,0)=-1$.

Okazuje się jednak, że wystarczy przyjąć naturalne założenie o ciągłości pochodnych cząstkowych mieszanych $\frac{\partial }{\partial x}\frac{\partial }{\partial y}f$ oraz $\frac{\partial }{\partial y}\frac{\partial }{\partial x}f$ w otoczeniu punktu $a$, aby mieć gwarancję ich równości w danym punkcie.

Uwaga 6.28.

Jeśli $f:{\mathbb{ℝ}}^n\supset U\ni x\mapsto f(x)\in \mathbb{ℝ}$ jest funkcją, która w punkcie $a\in U$ ma ciągłe

pochodne cząstkowe $\frac{\partial }{\partial x_j}\frac{\partial }{\partial x_i}f$ oraz $\frac{\partial }{\partial x_i}\frac{\partial }{\partial x_j}f$, to w punkcie $a$ są one równe, tj. $\frac{\partial }{\partial x_j}\frac{\partial }{\partial x_i}f(a)=\frac{\partial }{\partial x_i}\frac{\partial }{\partial x_j}f(a)$

Dowód uwagi pomijamy (można go znaleźć np. w podręczniku Ryszarda Rudnickiego Wykłady z analizy matematycznej, Wydawnictwo Naukowe PWN, Warszawa 2001).

W podobny sposób definiujemy pochodne cząstkowe wyższych rzędów. Wprowadźmy wygodne oznaczenie pochodnych cząstkowych za pomocą wielowskaźników $\alpha =({\alpha }_1,{\alpha }_2,\dots ,{\alpha }_n)\in {\mathbb{\mathbb{N}}}_0^n$. Niech $f:{\mathbb{ℝ}}^n\supset U\mapsto \mathbb{ℝ}$ będzie funkcją określoną na zbiorze otwartym $U$.

Oznaczmy symbolem $\frac{{\partial }^{{\alpha }_i}}{\partial x_i^{{\alpha }_i}}$ operację, która funkcji $f$ przypisuje pochodną cząstkową rzędu ${\alpha }_i$ po zmiennej $x_i$, o ile ta pochodna istnieje.

Definicja 6.29.

Załóżmy, że istnieją kolejno pochodne cząstkowe

$\frac{{\partial }^{{\alpha }_n}}{\partial x_n^{{\alpha }_n}}(\dots \frac{{\partial }^{{\alpha }_2}}{\partial x_2^{{\alpha }_2}}(\frac{{\partial }^{{\alpha }_1}}{\partial x_1^{{\alpha }_1}}f)\dots )(a)$

i nie zależą od kolejności różniczkowania. Mówimy wówczas, że funkcja $f$ ma pochodną

cząstkową $\frac{{\partial }^{|\alpha |}f(a)}{\partial x^{\alpha }}:=\frac{{\partial }^{{\alpha }_n}}{\partial x_n^{{\alpha }_n}}(\dots \frac{{\partial }^{{\alpha }_2}}{\partial x_2^{{\alpha }_2}}(\frac{{\partial }^{{\alpha }_1}}{\partial x_i^{{\alpha }_1}}f)\dots )(a)$

rzędu $|\alpha |={\alpha }_1+{\alpha }_2+\dots +{\alpha }_n$ w punkcie $a$. Pochodną tę notujemy też często symbolem $D^{\alpha }f(a)$.

Pochodne cząstkowe w fizyce. Elementy teorii pola

Niech $f:D\mapsto \mathbb{ℝ}$ będzie funkcją określoną na pewnym zbiorze otwartym $D\subset {\mathbb{ℝ}}^n$. Załóżmy, że w pewnym punkcie $a\in D$ istnieją pochodne cząstkowe $\frac{\partial f}{\partial x_1}(a),\frac{\partial f}{\partial x_2}(a),\dots ,\frac{\partial f}{\partial x_n}(a)$.

Definicja 6.30.

Wektor $\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}f(a)=(\frac{\partial f}{\partial x_1}(a),\frac{\partial f}{\partial x_2}(a),\dots ,\frac{\partial f}{\partial x_n}(a))\in {\mathbb{ℝ}}^n$ nazywamy gradientem funkcji $f$ w punkcie $a$. Wektor ten oznaczamy też często symbolem nabla: $\mathrm{\nabla }f(a)$. Punkt $a$, w którym wyznaczamy gradient funkcji $f$, zapisujemy czasem w formie indeksu dolnego: $\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}{}_{a}f$, ${\mathrm{\nabla }}_{a}f$.

Uwaga 6.31.

Jeśli funkcje $f,g:{\mathbb{ℝ}}^n\supset D\mapsto \mathbb{ℝ}$ mają w punkcie $a\in D$ pochodne cząstkowe $\frac{\partial f}{\partial x_i}(a)$, $\frac{\partial g}{\partial x_i}(a)$, $i=1,2,\dots ,n$, to

a) $\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}(f+g)(a)=\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}f(a)+\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}g(a)$,

b) $\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}(fg)(a)=g(a)\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}f(a)+f(a)\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}g(a)$.

Dowód 6.31.

Korzystając z twierdzenia o pochodnej sumy i iloczynu funkcji $f,g$, wyznaczamy kolejne współrzędne wektorów $\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}(f+g)(a)$ oraz $\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}(fg)(a)$:

$\frac{\partial }{\partial x_i}(f+g)(a)=\frac{\partial }{\partial x_i}f(a)+\frac{\partial }{\partial x_i}g(a)$

oraz

$\frac{\partial }{\partial x_i}(fg)(a)=g(a)\frac{\partial }{\partial x_i}f(a)+f(a)\frac{\partial }{\partial x_i}g(a)$,

gdy $i=1,2,\dots ,n$. Stąd otrzymujemy równości a) oraz b).

W ramach następnego modułu wykażemy, że

Uwaga 6.32.

Pochodna kierunkowa w kierunku gradientu jest największa.

W fizyce funkcję $f:{\mathbb{ℝ}}^3\mapsto \mathbb{ℝ}$ o wartościach liczbowych nazywa się funkcją skalarną, natomiast funkcję $F:{\mathbb{ℝ}}^3\mapsto {\mathbb{ℝ}}^3$ nazywa się polem (wektorowym). Przykładem funkcji skalarnych są np. temperatura, potencjał pola grawitacyjnego. Przykładem pola, które znamy ze szkoły, jest pole grawitacyjne. Jeśli w początku układu współrzędnych w przestrzeni ${\mathbb{ℝ}}^3$ znajduje się punkt materialny o masie $M$, to - zgodnie z prawem powszechnego ciążenia Newtona - na dowolny inny punkt materialny położony w punkcie $\overrightarrow{r}=(x,y,z)$ o masie $m$ działa siła $F=(F_x,F_y,F_z)$, której składowe wynoszą:

$\begin{array}{rl}{F}_x(\overrightarrow{r})& =-k\frac{x}{r^3},\\ F_y(\overrightarrow{r})& =-k\frac{y}{r^3},\\ F_z(\overrightarrow{r})& =-k\frac{z}{r^3},\end{array}$

gdzie $k=GmM$ jest iloczynem mas obu punktów materialnych i stałej grawitacji

$G=6,67259...\cdot 10^{-11}N\cdot m^2\cdot k{g}^{-2}$,

natomiast $r=\Vert \overrightarrow{r}{\Vert }_2=\Vert (x,y,z){\Vert }_2=\sqrt{x^2+y^2+z^2}$ jest odległością obu punktów. Zwróćmy uwagę, że

$F(\overrightarrow{r})=-\frac{k}{r^3}\overrightarrow{r}$,

stąd

$\Vert F(\overrightarrow{r}){\Vert }_2=\frac{k}{r^3}\Vert \overrightarrow{r}{\Vert }_2=\frac{k}{r^3}r=\frac{k}{r^2}$

siła ta jest odwrotnie proporcjonalna do kwadratu odległości obu punktów materialnych.

Definicja 6.33.

Pole wektorowe $F:{\mathbb{ℝ}}^3\supset D\mapsto {\mathbb{ℝ}}^3$ nazywamy polem potencjalnym, jeśli istnieje funkcja skalarna $U:D\mapsto \mathbb{ℝ}$ taka, że $\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}U(a)=F(a)$ w dowolnym punkcie $a$ zbioru otwartego $D\subset {\mathbb{ℝ}}^3$.

Funkcję $U$ nazywamy wówczas potencjałem pola wektorowego $F$.

Uwaga 6.34.

Pole grawitacyjne $F(\overrightarrow{r})=-\frac{k}{r^3}\overrightarrow{r}$ jest polem potencjalnym. Potencjałem tego pola jest funkcja skalarna $U(\overrightarrow{r})=\frac{k}{r}$, gdzie (jak powyżej) $\overrightarrow{r}=(x,y,z)$

oraz $r=\Vert \overrightarrow{r}{\Vert }_2=\sqrt{x^2+y^2+z^2}$.

Dowód 6.34.

Policzmy pochodne cząstkowe funkcji $U(x,y,z)=\frac{k}{\sqrt{x^2+y^2+z^2}}=\frac{k}{r}$ określonej w zbiorze otwartym $D={\mathbb{ℝ}}^3\setminus \{0\}$, czyli wszędzie w przestrzeni ${\mathbb{ℝ}}^3$ poza początkiem układu współrzędnych. Mamy

$\begin{array}{rl}\frac{\partial }{\partial x}U(\overrightarrow{r})& =\frac{\partial }{\partial x}\frac{k}{r}=-\frac{k}{r^2}\cdot \frac{\partial r}{\partial x}=-\frac{k}{r^2}\cdot \frac{2x}{2r}=-\frac{k}{r^3}x\\ \frac{\partial }{\partial y}U(\overrightarrow{r})& =\frac{\partial }{\partial y}\frac{k}{r}=-\frac{k}{r^2}\cdot \frac{\partial r}{\partial y}=-\frac{k}{r^2}\cdot \frac{2y}{2r}=-\frac{k}{r^3}y\\ \frac{\partial }{\partial z}U(\overrightarrow{r})& =\frac{\partial }{\partial z}\frac{k}{r}=-\frac{k}{r^2}\cdot \frac{\partial r}{\partial z}=-\frac{k}{r^2}\cdot \frac{2z}{2r}=-\frac{k}{r^3}z,\end{array}$

czyli

$\begin{array}{rl}\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}U(\overrightarrow{r})& =\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}U(x,y,z)=(-\frac{k}{r^3}x,-\frac{k}{r^3}y,-\frac{k}{r^3}z)\\ & =-\frac{k}{r^3}(x,y,z)=-\frac{k}{r^3}\overrightarrow{r}\\ & =F(\overrightarrow{r}).\end{array}$

Definicja 6.35.

Dywergencją pola wektorowego $F=(F_x,F_y,F_z):{\mathbb{ℝ}}^3\supset D\mapsto {\mathbb{ℝ}}^3$ w punkcie $a\in D$ nazywamy liczbę

$\mathrm{d}\mathrm{i}\mathrm{v}F(a)=\frac{\partial F_x}{\partial x}(a)+\frac{\partial F_y}{\partial y}(a)+\frac{\partial F_z}{\partial z}(a)$,

o ile istnieją pochodne cząstkowe $\frac{\partial F_x}{\partial x}(a),\frac{\partial F_y}{\partial y}(a),\frac{\partial F_z}{\partial z}(a)$. Jeśli w dowolnym punkcie $a\in D$ dywergencja $\mathrm{d}\mathrm{i}\mathrm{v}F(a)=0$, to pole wektorowe $F$ nazywamy polem bezźródłowym.

Uwaga 6.36.

Pole grawitacyjne $F(\overrightarrow{r})=-\frac{k}{r^3}\overrightarrow{r}$ jest polem bezźródłowym w ${\mathbb{ℝ}}^3\setminus \{0\}$.

Dowód 6.36.

W dowolnym punkcie $\overrightarrow{r}=(x,y,z)\ne 0$ mamy

$\begin{array}{lll}\frac{\partial F_x}{\partial x}(\overrightarrow{r})& =& \frac{\partial }{\partial x}(-\frac{k}{r^3}x)=-k(\frac{\partial x}{\partial x}\frac{1}{r^3}+x\frac{\partial }{\partial x}\frac{1}{r^3})\\ & =& -k(\frac{1}{r^3}+x\cdot \frac{(-3)}{r^4}\frac{\partial r}{\partial x})=-k(\frac{1}{r^3}-\frac{3x^2}{r^5})\end{array}$ i podobnie $\frac{\partial F_y}{\partial y}(\overrightarrow{r})=-k(\frac{1}{r^3}-\frac{3y^2}{r^5})\text{ oraz }\frac{\partial F_z}{\partial z}(\overrightarrow{r})=-k(\frac{1}{r^3}-\frac{3z^2}{r^5})$.

Stąd

$\begin{array}{rl}\mathrm{d}\mathrm{i}\mathrm{v}F(\overrightarrow{r})& =\frac{\partial F_x}{\partial x}(\overrightarrow{r})+\frac{\partial F_y}{\partial y}(\overrightarrow{r})+\frac{\partial F_z}{\partial z}(\overrightarrow{r})\\ & =-k(\frac{1}{r^3}-\frac{3x^2}{r^5})-k(\frac{1}{r^3}-\frac{3y^2}{r^5})-k(\frac{1}{r^3}-\frac{3z^2}{r^5})\\ & =-k(\frac{3}{r^3}-\frac{3(x^2+y^2+z^2)}{r^5})=-k(\frac{3}{r^3}-\frac{3r^2}{r^5})=0.\end{array}$

Definicja 6.37.

Rotacją pola wektorowego $F=(F_x,F_y,F_z):{\mathbb{ℝ}}^3\supset D\mapsto {\mathbb{ℝ}}^3$ w punkcie $a\in D$

nazywamy wektor $\mathrm{r}\mathrm{o}\mathrm{t}F(a)=(\frac{\partial F_z}{\partial y}(a)-\frac{\partial F_y}{\partial z}(a),\frac{\partial F_x}{\partial z}(a)-\frac{\partial F_z}{\partial x}(a),\frac{\partial F_y}{\partial x}(a)-\frac{\partial F_x}{\partial y}(a))$.

Wektor ten oznaczamy też czasem symbolem $\mathrm{\nabla }\times F(a)$. Jeśli w każdym punkcie $a\in D$ rotacja $\mathrm{r}\mathrm{o}\mathrm{t}F(a)=0$, to pole wektorowe $F$ nazywamy bezwirowym.

Uwaga 6.38.

Pole grawitacyjne $F(\overrightarrow{r})=-\frac{k}{r^3}\overrightarrow{r}$ jest polem bezwirowym w ${\mathbb{ℝ}}^3\setminus \{0\}$.

Dowód 6.38.

W dowolnym punkcie $\overrightarrow{r}=(x,y,z)\ne 0$ mamy

$\frac{\partial }{\partial y}{F}_z(\overrightarrow{r})=\frac{\partial }{\partial y}(-k\frac{1}{r^3}z)=-kz\frac{\partial }{\partial y}\left(\frac{1}{r^3}\right)=-kz\frac{(-3)}{r^4}\frac{\partial r}{\partial y}=3kz\frac{y}{r^5}=zy\frac{3k}{r^5}$

oraz podobnie

$\frac{\partial }{\partial z}{F}_y(\overrightarrow{r})=yz\frac{3k}{r^5}$

Pierwsza współrzędna wektora rotacji jest więc równa zeru, gdyż

$\frac{\partial F_z}{\partial y}(\overrightarrow{r})-\frac{\partial F_y}{\partial z}(\overrightarrow{r})=zy\frac{3k}{r^5}-yz\frac{3k}{r^5}=0$.

W ten sam sposób sprawdzamy, że również druga i trzecia współrzędna wektora rotacji zerują się:

$\begin{array}{rl}\frac{\partial F_x}{\partial z}(\overrightarrow{r})-\frac{\partial F_z}{\partial x}(\overrightarrow{r})& =xz\frac{3k}{r^5}-zx\frac{3k}{r^5}=0\\ \frac{\partial F_y}{\partial x}(\overrightarrow{r})-\frac{\partial F_x}{\partial y}(\overrightarrow{r})& =yx\frac{3k}{r^5}-xy\frac{3k}{r^5}=0.\end{array}$

Stąd $\mathrm{r}\mathrm{o}\mathrm{t}F(\overrightarrow{r})=0$, dla $\overrightarrow{r}\ne 0$.