Analiza matematyczna 1/Wykład 15: Krzywe i bryły obrotowe: Różnice pomiędzy wersjami

Z Studia Informatyczne
Przejdź do nawigacjiPrzejdź do wyszukiwania
m Zastępowanie tekstu – „ </math>” na „</math>”
m Zastępowanie tekstu – „<math> ” na „<math>”
Linia 4: Linia 4:
zwyczajnej.
zwyczajnej.
Definiujemy długość krzywej i krzywą prostowalną.
Definiujemy długość krzywej i krzywą prostowalną.
Dowodzimy, że krzywa zwyczajna klasy <math> C^1</math> jest prostowalna.
Dowodzimy, że krzywa zwyczajna klasy <math>C^1</math> jest prostowalna.
Wyprowadzamy wzór na długość krzywej i liczymy długości
Wyprowadzamy wzór na długość krzywej i liczymy długości
cykloidy i asteroidy.
cykloidy i asteroidy.
Linia 12: Linia 12:
==Długość krzywej==
==Długość krzywej==
{{definicja|15.1.|definicja_15_1|
{{definicja|15.1.|definicja_15_1|
Niech <math> -\infty<a<b<+\infty</math>.
Niech <math>-\infty<a<b<+\infty</math>.
'''''Krzywą''''' nazywamy
'''''Krzywą''''' nazywamy
zbiór punktów
zbiór punktów


<center>
<center>
<math> K
<math>K
=
=
\big\{(x,y)\in\mathbb{R}^2:\ x=\varphi(t),\ y=\psi(t),\ t\in[a,b]\big\},
\big\{(x,y)\in\mathbb{R}^2:\ x=\varphi(t),\ y=\psi(t),\ t\in[a,b]\big\},
Linia 23: Linia 23:
</center>
</center>


gdzie <math> \varphi,\psi\colon[a,b]\longrightarrow\mathbb{R}</math> są dwiema funkcjami
gdzie <math>\varphi,\psi\colon[a,b]\longrightarrow\mathbb{R}</math> są dwiema funkcjami
ciągłymi. Piszemy:
ciągłymi. Piszemy:


<center>
<center>
<math>  
<math>
K=K(\varphi,\psi):
K=K(\varphi,\psi):
\left\{ \begin{array} {l}
\left\{ \begin{array} {l}
Linia 131: Linia 131:
<br>
<br>
<center>
<center>
<math> \big(\varphi(t_0),\psi(t_0)\big),
<math>\big(\varphi(t_0),\psi(t_0)\big),
\ \ldots,
\ \ldots,
\big(\varphi(t_n),\psi(t_n)\big)
\big(\varphi(t_n),\psi(t_n)\big)
Linia 139: Linia 139:


nazywamy
nazywamy
'''''łamaną wpisaną w krzywą <math> K</math>'''''.
'''''łamaną wpisaną w krzywą <math>K</math>'''''.
Przez <math> l(p)</math> oznaczamy '''''długość'''''
Przez <math>l(p)</math> oznaczamy '''''długość'''''
łamanej <math> p</math> (to znaczy sumę długości odcinków
łamanej <math>p</math> (to znaczy sumę długości odcinków
wchodzących w skład łamanej).}}
wchodzących w skład łamanej).}}


{{definicja|15.5.||
{{definicja|15.5.||


Długością krzywej <math> K</math> nazywamy liczbę:
Długością krzywej <math>K</math> nazywamy liczbę:


<center>
<center>
<math> l(K)
<math>l(K)
=
=
\sup_p l(p),
\sup_p l(p),
Linia 155: Linia 155:
</center>
</center>


gdzie supremum jest brane po wszystkich łamanych wpisanych w <math> K</math>.}}
gdzie supremum jest brane po wszystkich łamanych wpisanych w <math>K</math>.}}




Linia 171: Linia 171:
{{definicja|15.6.||
{{definicja|15.6.||


Jeśli <math> l(K)<+\infty</math>, to mówimy, że krzywa <math> K</math> jest
Jeśli <math>l(K)<+\infty</math>, to mówimy, że krzywa <math>K</math> jest
'''''prostowalna'''''.
'''''prostowalna'''''.
}}
}}
Linia 177: Linia 177:
<span id="twierdzenie_15_7">{{twierdzenie|15.7.||
<span id="twierdzenie_15_7">{{twierdzenie|15.7.||


Niech <math> \varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math> C^1</math>
Niech <math>\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math>C^1</math>
oraz niech <math> K=K(\varphi,\psi)</math> będzie krzywą
oraz niech <math>K=K(\varphi,\psi)</math> będzie krzywą
zwyczajną.<br>
zwyczajną.<br>
Wówczas krzywa <math> K</math> jest prostowalna.
Wówczas krzywa <math>K</math> jest prostowalna.
}}</span>
}}</span>


Linia 186: Linia 186:


{{dowod|15.7. [nadobowiązkowy]||
{{dowod|15.7. [nadobowiązkowy]||
Niech <math> p</math> będzie dowolną łamaną wpisaną w
Niech <math>p</math> będzie dowolną łamaną wpisaną w
krzywą <math> K</math>,
krzywą <math>K</math>,
to znaczy istnieje podział
to znaczy istnieje podział


<center>
<center>
<math> a
<math>a
=
=
t_0
t_0
Linia 205: Linia 205:
</center>
</center>


taki, że <math> p</math> jest łamaną o wierzchołkach
taki, że <math>p</math> jest łamaną o wierzchołkach
<math>(x_i,y_i)</math> dla <math> i=0,\ldots,n</math>, gdzie
<math>(x_i,y_i)</math> dla <math>i=0,\ldots,n</math>, gdzie


<center>
<center>
<math>  
<math>
\left\{ \begin{array} {l}
\left\{ \begin{array} {l}
x_i= \varphi(t_i)\\
x_i= \varphi(t_i)\\
Linia 219: Linia 219:
</center>
</center>


Długość łamanej <math> p</math> wyraża się wzorem:
Długość łamanej <math>p</math> wyraża się wzorem:


<center>
<center>
<math> l(p)
<math>l(p)
=
=
\sum_{i=1}^n
\sum_{i=1}^n
Linia 229: Linia 229:
</center>
</center>


Ponieważ <math> \varphi,\psi\in C^1\big([a,b];\mathbb{R}\big)</math>,
Ponieważ <math>\varphi,\psi\in C^1\big([a,b];\mathbb{R}\big)</math>,
więc z twierdzenia o wartości średniej
więc z twierdzenia o wartości średniej
(patrz [[Analiza matematyczna 1/Wykład 9: Pochodna funkcji jednej zmiennej#twierdzenie_9_37|twierdzenie 9.37.]]) mamy
(patrz [[Analiza matematyczna 1/Wykład 9: Pochodna funkcji jednej zmiennej#twierdzenie_9_37|twierdzenie 9.37.]]) mamy


<center>
<center>
<math> x_i-x_{i-1}
<math>x_i-x_{i-1}
=
=
\varphi(t_i)-\varphi(t_{i-1})
\varphi(t_i)-\varphi(t_{i-1})
Linia 241: Linia 241:
</center>
</center>
<center>
<center>
<math> y_i-y_{i-1}
<math>y_i-y_{i-1}
=
=
\psi(t_i)-\psi(t_{i-1})
\psi(t_i)-\psi(t_{i-1})
Linia 255: Linia 255:


<center>
<center>
<math> l(p)
<math>l(p)
=
=
\sum_{i=1}^n
\sum_{i=1}^n
Linia 262: Linia 262:
</center>
</center>


Ponieważ <math> \varphi',\psi'\in C\big([a,b];\mathbb{R}\big)</math>
Ponieważ <math>\varphi',\psi'\in C\big([a,b];\mathbb{R}\big)</math>
i przedział <math>[a,b]</math> jest zwarty,
i przedział <math>[a,b]</math> jest zwarty,
więc funkcje <math>\varphi',\psi'</math> są ograniczone.<br>
więc funkcje <math>\varphi',\psi'</math> są ograniczone.<br>
Linia 268: Linia 268:


<center>
<center>
<math> M =
<math>M =
\sup_{t\in[a,b]}\varphi'(t),
\sup_{t\in[a,b]}\varphi'(t),
\qquad
\qquad
Linia 278: Linia 278:


<center>
<center>
<math> m =
<math>m =
\inf_{t\in[a,b]}\varphi'(t),
\inf_{t\in[a,b]}\varphi'(t),
\qquad
\qquad
Linia 289: Linia 289:


<center>
<center>
<math> \sqrt{m^2+{m^*}^2}\cdot(b-a)
<math>\sqrt{m^2+{m^*}^2}\cdot(b-a)
\le
\le
l(p)
l(p)
Linia 298: Linia 298:


Ponieważ powyższa nierówność zachodzi dla dowolnej
Ponieważ powyższa nierówność zachodzi dla dowolnej
łamanej <math> p</math> wpisanej w krzywą <math> K</math>,
łamanej <math>p</math> wpisanej w krzywą <math>K</math>,
więc przechodząc do supremum po wszystkich takich
więc przechodząc do supremum po wszystkich takich
łamanych, dostajemy
łamanych, dostajemy


<center>
<center>
<math> \sqrt{m^2+{m^*}^2}\cdot(b-a)
<math>\sqrt{m^2+{m^*}^2}\cdot(b-a)
\le
\le
l(K)
l(K)
Linia 311: Linia 311:
</center>
</center>


a zatem krzywa <math> K</math> jest prostowalna.
a zatem krzywa <math>K</math> jest prostowalna.
}}
}}


Linia 317: Linia 317:


W powyższym twierdzeniu (jak i twierdzeniach
W powyższym twierdzeniu (jak i twierdzeniach
następnych) zakładamy, że krzywa jest klasy <math> C^1</math>.  
następnych) zakładamy, że krzywa jest klasy <math>C^1</math>.  
(to znaczy <math>\varphi,\psi</math>, są klasy <math>C^1</math>)
(to znaczy <math>\varphi,\psi</math>, są klasy <math>C^1</math>)
W zastosowaniach okaże się jednak, że często będziemy mieli do
W zastosowaniach okaże się jednak, że często będziemy mieli do
czynienia z krzywymi, które są tylko ciągłe,
czynienia z krzywymi, które są tylko ciągłe,
zwyczajne oraz
zwyczajne oraz
"kawałkami" klasy <math> C^1</math>, to znaczy krzywą można otrzymać jako
"kawałkami" klasy <math>C^1</math>, to znaczy krzywą można otrzymać jako
"sklejenie" kilku krzywych klasy <math> C^1</math>
"sklejenie" kilku krzywych klasy <math>C^1</math>
(przy sklejaniu początek następnej krzywej jest końcem
(przy sklejaniu początek następnej krzywej jest końcem
poprzedniej).
poprzedniej).
Wszystkie wypowiadane tu twierdzenia dla krzywych klasy <math> C^1</math>
Wszystkie wypowiadane tu twierdzenia dla krzywych klasy <math>C^1</math>
stosują się także do krzywych kawałkami klasy <math> C^1</math>.
stosują się także do krzywych kawałkami klasy <math>C^1</math>.
}}
}}


Linia 334: Linia 334:
{{definicja|15.9.||
{{definicja|15.9.||


Niech <math> K=K(\varphi,\psi)</math> będzie krzywą.
Niech <math>K=K(\varphi,\psi)</math> będzie krzywą.
Zdefiniujmy:
Zdefiniujmy:


<center>
<center>
<math> K(t)
<math>K(t)
\ \ \stackrel{df}{=}
\ \ \stackrel{df}{=}
\bigg\{\big(\varphi(\tau),\psi(\tau)\big):\ \tau\in[a,t]\bigg\},
\bigg\{\big(\varphi(\tau),\psi(\tau)\big):\ \tau\in[a,t]\bigg\},
Linia 347: Linia 347:


<center>
<center>
<math>  
<math>
s(t) \ \ \stackrel{df}{=}
s(t) \ \ \stackrel{df}{=}
l\big(K(t)\big)\quad</math> (długośćkrzywejK(t)) <math> .
l\big(K(t)\big)\quad</math> (długośćkrzywejK(t)) <math> .
</math>
</math>
</center>
</center>


W szczególności <math> s(b)=l(K)</math>.
W szczególności <math>s(b)=l(K)</math>.
}}
}}


{{twierdzenie|15.10.||
{{twierdzenie|15.10.||


Niech <math> \varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math> C^1</math>
Niech <math>\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math>C^1</math>
oraz niech <math> K=K(\varphi,\psi)</math> będzie krzywą
oraz niech <math>K=K(\varphi,\psi)</math> będzie krzywą
zwyczajną.<br>
zwyczajną.<br>
Wówczas
Wówczas


<center>
<center>
<math> s'(t)
<math>s'(t)
=
=
\sqrt{\varphi'(t)^2+\psi'(t)^2}
\sqrt{\varphi'(t)^2+\psi'(t)^2}
Linia 374: Linia 374:


{{dowod|15.10. [nadobowiązkowy]||
{{dowod|15.10. [nadobowiązkowy]||
Niech <math> t_0,t_0+h\in[a,b]</math>.
Niech <math>t_0,t_0+h\in[a,b]</math>.
Analogicznie do ostatniego oszacowania
Analogicznie do ostatniego oszacowania
w dowodzie [[#twierdzenie_15_7|twierdzenia 15.7.]] dostajemy:
w dowodzie [[#twierdzenie_15_7|twierdzenia 15.7.]] dostajemy:


<center>
<center>
<math> \sqrt{m_h^2+{m_h^*}^2}\cdot h
<math>\sqrt{m_h^2+{m_h^*}^2}\cdot h
\le
\le
s(t_0+h)-s(t_0)
s(t_0+h)-s(t_0)
Linia 390: Linia 390:


<center>
<center>
<math> M_h
<math>M_h
=
=
\sup_{t\in[t_0,t_0+h]}\varphi'(t),
\sup_{t\in[t_0,t_0+h]}\varphi'(t),
Linia 401: Linia 401:


<center>
<center>
<math> m_h
<math>m_h
=
=
\inf_{t\in[t_0,t_0+h]}\varphi'(t),
\inf_{t\in[t_0,t_0+h]}\varphi'(t),
Linia 413: Linia 413:
Dzielimy wszystkie strony
Dzielimy wszystkie strony
powyższego oszacowania
powyższego oszacowania
przez <math> h</math>, dostając:
przez <math>h</math>, dostając:


<center>
<center>
<math> \sqrt{m_h^2+{m_h^*}^2}
<math>\sqrt{m_h^2+{m_h^*}^2}
\le
\le
\frac{s(t_0+h)-s(t_0)}{h}
\frac{s(t_0+h)-s(t_0)}{h}
Linia 424: Linia 424:
</center>
</center>


Ponieważ funkcje <math> \varphi'</math> i <math> \psi'</math> są ciągłe,
Ponieważ funkcje <math>\varphi'</math> i <math>\psi'</math> są ciągłe,
więc dostajemy
więc dostajemy


<center>
<center>
<math> \begin{align}  
<math>\begin{align}  
M_h  &  \xrightarrow[h\rightarrow 0]{} & \varphi'(t_0),\\
M_h  &  \xrightarrow[h\rightarrow 0]{} & \varphi'(t_0),\\
m_h  &  \xrightarrow[h\rightarrow 0]{} & \varphi'(t_0),\\
m_h  &  \xrightarrow[h\rightarrow 0]{} & \varphi'(t_0),\\
Linia 439: Linia 439:


<center>
<center>
<math> s'(t_0)
<math>s'(t_0)
=
=
\lim_{h\rightarrow 0}\frac{s(t_0+h)-s(t_0)}{h}
\lim_{h\rightarrow 0}\frac{s(t_0+h)-s(t_0)}{h}
Linia 450: Linia 450:


<span id="twierdzenie_15_11">{{twierdzenie|15.11. [O długości krzywej]||
<span id="twierdzenie_15_11">{{twierdzenie|15.11. [O długości krzywej]||
Niech <math> \varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math> C^1</math>
Niech <math>\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math>C^1</math>
oraz niech <math> K=K(\varphi,\psi)</math> będzie krzywą zwyczajną.
oraz niech <math>K=K(\varphi,\psi)</math> będzie krzywą zwyczajną.
Wówczas długość krzywej wyraża się wzorem
Wówczas długość krzywej wyraża się wzorem


<center>
<center>
<math> l(K)
<math>l(K)
=
=
\int\limits_a^b\sqrt{\varphi'(\tau)^2+\psi'(\tau)^2}\,d\tau.
\int\limits_a^b\sqrt{\varphi'(\tau)^2+\psi'(\tau)^2}\,d\tau.
Linia 462: Linia 462:


W szczególności, jeśli krzywa zadana jest wykresem funkcji
W szczególności, jeśli krzywa zadana jest wykresem funkcji
<math> y=f(x)</math>, dla <math> x\in[a,b]</math>,
<math>y=f(x)</math>, dla <math>x\in[a,b]</math>,
to
to


<center>
<center>
<math> l(K)
<math>l(K)
=
=
\int\limits_a^b\sqrt{1+f'(t)^2}\,dt.
\int\limits_a^b\sqrt{1+f'(t)^2}\,dt.
Linia 477: Linia 477:


<center>
<center>
<math> l(K)
<math>l(K)
=
=
s(b)
s(b)
Linia 489: Linia 489:
</center>
</center>


W drugim przypadku krzywą zadaną przez funkcję <math> f</math>
W drugim przypadku krzywą zadaną przez funkcję <math>f</math>
możemy zapisać w postaci parametrycznej
możemy zapisać w postaci parametrycznej


<center>
<center>
<math>  
<math>
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 514: Linia 514:


<center>
<center>
<math> r
<math>r
=
=
g(\vartheta)
g(\vartheta)
Linia 526: Linia 526:


<center>
<center>
<math>  
<math>
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 561: Linia 561:


<center>
<center>
<math> l(K)
<math>l(K)
=
=
\int\limits_{\alpha}^{\beta}
\int\limits_{\alpha}^{\beta}
Linia 583: Linia 583:


'''''Cykloidą''''' nazywamy krzywą kreśloną
'''''Cykloidą''''' nazywamy krzywą kreśloną
przez ustalony punkt <math> 0</math> na okręgu toczącym się po prostej <math> l</math>.}}
przez ustalony punkt <math>0</math> na okręgu toczącym się po prostej <math>l</math>.}}


[[File:AM1.M15.W.R11.svg|250x250px|thumb|right|Cykloida]]
[[File:AM1.M15.W.R11.svg|250x250px|thumb|right|Cykloida]]
Linia 591: Linia 591:


Oznaczenia:<br>
Oznaczenia:<br>
<math> a</math> - promień okręgu;<br>
<math>a</math> - promień okręgu;<br>
<math> O</math> - początkowy punkt styczności okręgu i prostej
<math>O</math> - początkowy punkt styczności okręgu i prostej
<math> l</math>;<br>
<math>l</math>;<br>
<math> N</math> - nowy punkt styczności;<br>
<math>N</math> - nowy punkt styczności;<br>
<math> M</math> - nowe położenie punktu <math> O</math>;<br>
<math>M</math> - nowe położenie punktu <math>O</math>;<br>
<math>t=\sphericalangle NDM</math> - parametr określający
<math>t=\sphericalangle NDM</math> - parametr określający
położenie punktu <math> M</math>.
położenie punktu <math>M</math>.


Liczymy współrzędne punktu <math> M(x,y)</math>:
Liczymy współrzędne punktu <math>M(x,y)</math>:


<center>
<center>
<math> x
<math>x
\ =
\ =
OF
OF
Linia 615: Linia 615:


<center>
<center>
<math> y
<math>y
\ =  
\ =  
FM
FM
Linia 630: Linia 630:


<center>
<center>
<math>  
<math>
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 639: Linia 639:
\qquad
\qquad
t\in [0,2\pi]
t\in [0,2\pi]
\quad(</math> lub <math> \ t\in\mathbb{R}).
\quad(</math> lub <math> \ t\in\mathbb{R}).
</math>
</math>
</center>}}
</center>}}
Linia 648: Linia 648:


<center>
<center>
<math>  
<math>
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 678: Linia 678:


<center>
<center>
<math> \begin{array}{lll}l(K)
<math>\begin{array}{lll}l(K)
&=&\int\limits_0^{2\pi}
&=&\int\limits_0^{2\pi}
\sqrt{x'(t)^2+y'(t)^2}\,dt
\sqrt{x'(t)^2+y'(t)^2}\,dt
Linia 697: Linia 697:


<center>
<center>
<math> x^{\frac{2}{3}}+y^{\frac{2}{3}}
<math>x^{\frac{2}{3}}+y^{\frac{2}{3}}
=
=
a^{\frac{2}{3}}.
a^{\frac{2}{3}}.
Linia 706: Linia 706:


<center>
<center>
<math>  
<math>
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 721: Linia 721:


<center>
<center>
<math> \sqrt{x'(t)^2+y'(t)^2}
<math>\sqrt{x'(t)^2+y'(t)^2}
=
=
3a\sin t\cos t
3a\sin t\cos t
Linia 731: Linia 731:


<center>
<center>
<math> l(K)
<math>l(K)
=
=
4\int\limits_0^{\frac{\pi}{2}}3a\sin t\cos t\,dt
4\int\limits_0^{\frac{\pi}{2}}3a\sin t\cos t\,dt
Linia 755: Linia 755:
==Całka krzywoliniowa==
==Całka krzywoliniowa==


Niech <math> K</math> będzie krzywą klasy <math> C^1</math>:
Niech <math>K</math> będzie krzywą klasy <math>C^1</math>:


<center>
<center>
<math> K
<math>K
=
=
\big\{(x,y)\in\mathbb{R}^2:\ x=\varphi(t),\ y=\psi(t),\ t\in[a,b]\big\},
\big\{(x,y)\in\mathbb{R}^2:\ x=\varphi(t),\ y=\psi(t),\ t\in[a,b]\big\},
Linia 765: Linia 765:


Przypuśćmy, że określona jest funkcja ciągła
Przypuśćmy, że określona jest funkcja ciągła
<math> f\colon K\ni M\longmapsto f(M)\in\mathbb{R}</math>,
<math>f\colon K\ni M\longmapsto f(M)\in\mathbb{R}</math>,
to znaczy funkcja, która każdemu punktowi <math> M</math>
to znaczy funkcja, która każdemu punktowi <math>M</math>
krzywej <math> K</math> przyporządkowuje pewną wartość rzeczywistą
krzywej <math>K</math> przyporządkowuje pewną wartość rzeczywistą
<math> f(M)</math>.
<math>f(M)</math>.
Okazuje się, że dla takich funkcji możemy także zdefiniować
Okazuje się, że dla takich funkcji możemy także zdefiniować
całkę oznaczoną, to znaczy całkę z funkcji <math> f</math> po krzywej <math> K</math>.
całkę oznaczoną, to znaczy całkę z funkcji <math>f</math> po krzywej <math>K</math>.
[[grafika:Riemann.jpg|thumb|right||Georg Friedrich Bernhard Riemann (1826-1866)<br>[[Biografia Riemann|Zobacz biografię]]]]
[[grafika:Riemann.jpg|thumb|right||Georg Friedrich Bernhard Riemann (1826-1866)<br>[[Biografia Riemann|Zobacz biografię]]]]
Całkę tę wprowadza się analogicznie jak całkę Riemanna na
Całkę tę wprowadza się analogicznie jak całkę Riemanna na
Linia 777: Linia 777:


<center>
<center>
<math> \int\limits_K f(x,y)\,ds
<math>\int\limits_K f(x,y)\,ds
=
=
\int\limits_a^b
\int\limits_a^b
Linia 788: Linia 788:
długością są pomijalne).
długością są pomijalne).


Jeśli  mamy daną krzywą (pręt) <math> K</math> zadaną jak wyżej,
Jeśli  mamy daną krzywą (pręt) <math>K</math> zadaną jak wyżej,
o gęstości w każdym jej punkcie <math> M(x,y)</math> danej funkcją
o gęstości w każdym jej punkcie <math>M(x,y)</math> danej funkcją
ciągłą  <math> \varrho(M)</math>,
ciągłą  <math>\varrho(M)</math>,
to masa tego pręta wyraża się wzorem
to masa tego pręta wyraża się wzorem


<center>
<center>
<math> m
<math>m
=
=
\int\limits_K
\int\limits_K
Linia 805: Linia 805:


<center>
<center>
<math> \begin{align} x_0 & = & \frac{1}{m}\int\limits_K x\cdot \varrho(x,y)\,ds,\\
<math>\begin{align} x_0 & = & \frac{1}{m}\int\limits_K x\cdot \varrho(x,y)\,ds,\\
x_0 & = & \frac{1}{m}\int\limits_K y\cdot \varrho(x,y)\,ds.
x_0 & = & \frac{1}{m}\int\limits_K y\cdot \varrho(x,y)\,ds.
\end{align}</math>
\end{align}</math>
Linia 813: Linia 813:


Obliczyć masę pręta półkolistego
Obliczyć masę pręta półkolistego
<math> K=\{(x,y)\in\mathbb{R}^2:\ x^2+y^2=R^2,\ y\ge 0\}</math>
<math>K=\{(x,y)\in\mathbb{R}^2:\ x^2+y^2=R^2,\ y\ge 0\}</math>
o gęstości <math> \varrho(x,y)=y^2</math>.
o gęstości <math>\varrho(x,y)=y^2</math>.


Masa krzywej o gęstości <math> \varrho</math> dana jest wzorem
Masa krzywej o gęstości <math>\varrho</math> dana jest wzorem


<center>
<center>
<math> m
<math>m
=
=
\int\limits_K \varrho(x,y)\,ds
\int\limits_K \varrho(x,y)\,ds
Linia 831: Linia 831:


<center>
<center>
<math>  
<math>
K=K(\varphi,\psi):
K=K(\varphi,\psi):
\left\{
\left\{
Linia 863: Linia 863:
{{przyklad|15.18.||
{{przyklad|15.18.||


Obliczyć masę i współrzędne środka ciężkości odcinka <math> K</math> łączącego
Obliczyć masę i współrzędne środka ciężkości odcinka <math>K</math> łączącego
punkt <math>(0,0)</math> z punktem <math>(1,1)</math> o gęstości wprost
punkt <math>(0,0)</math> z punktem <math>(1,1)</math> o gęstości wprost
proporcjonalnej
proporcjonalnej
do odległości punktu od środka układu i równej <math> \sqrt{2}</math> w
do odległości punktu od środka układu i równej <math>\sqrt{2}</math> w
punkcie <math>(1,1)</math>.
punkcie <math>(1,1)</math>.
    
    
Skoro gęstość  <math> \varrho</math> jest proporcjonalna
Skoro gęstość  <math>\varrho</math> jest proporcjonalna
do odległości punktu od środka układu i wynosi <math> \sqrt{2}</math> w
do odległości punktu od środka układu i wynosi <math>\sqrt{2}</math> w
punkcie <math>(1,1)</math>, to
punkcie <math>(1,1)</math>, to


<center>
<center>
<math> \varrho(x,t)
<math>\varrho(x,t)
=
=
c\sqrt{x^2+y^2}
c\sqrt{x^2+y^2}
\quad</math> oraz <math> \quad
\quad</math> oraz <math> \quad
\varrho(1,1)=c\sqrt{2}=\sqrt{2},
\varrho(1,1)=c\sqrt{2}=\sqrt{2},
</math>
</math>
</center>
</center>


stąd <math> c=1</math>.
stąd <math>c=1</math>.
Parametryzacją odcinka jest na przykład
Parametryzacją odcinka jest na przykład


<center>
<center>
<math>  
<math>
K=K(\varphi,\psi):
K=K(\varphi,\psi):
\left\{
\left\{
Linia 901: Linia 901:


<center>
<center>
<math> m
<math>m
=
=
\int\limits_K\sqrt{x^2+y^2}\,ds
\int\limits_K\sqrt{x^2+y^2}\,ds
Linia 919: Linia 919:


<center>
<center>
<math> x_0
<math>x_0
=
=
\frac{1}{m}\int\limits_K x\cdot \varrho(x,y)\,ds
\frac{1}{m}\int\limits_K x\cdot \varrho(x,y)\,ds
Linia 933: Linia 933:
</center>
</center>


Z symetrii zadania wynika, że <math> y_0=\frac{2}{3}</math>.}}
Z symetrii zadania wynika, że <math>y_0=\frac{2}{3}</math>.}}


==Pole powierzchni i objętość bryły obrotowej==
==Pole powierzchni i objętość bryły obrotowej==


W dalszej części wykładu będziemy zakładać, że krzywe są klasy
W dalszej części wykładu będziemy zakładać, że krzywe są klasy
<math> C^1</math>. Podamy liczne wzory na obliczanie pól powierzchni i
<math>C^1</math>. Podamy liczne wzory na obliczanie pól powierzchni i
objętości brył obrotowych, w większości pozostawiając je bez
objętości brył obrotowych, w większości pozostawiając je bez
dowodów (podając natomiast pewne ich uzasadnienia).
dowodów (podając natomiast pewne ich uzasadnienia).
Linia 954: Linia 954:


<center>
<center>
<math> y=f_1(x)
<math>y=f_1(x)
\quad</math> i <math> \quad
\quad</math> i <math> \quad
y=f_2(x)
y=f_2(x)
\quad x\in[a,b],
\quad x\in[a,b],
Linia 964: Linia 964:


<center>
<center>
<math> |P|
<math>|P|
=
=
\int\limits_a^b\big[f_1(x)-f_2(x)\big]\,dx
\int\limits_a^b\big[f_1(x)-f_2(x)\big]\,dx
Linia 981: Linia 981:


<center>
<center>
<math> K:
<math>K:
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 989: Linia 989:
\right.,
\right.,
\qquad
\qquad
</math> dla <math> \ t\in[\alpha,\beta],
</math> dla <math> \ t\in[\alpha,\beta],
</math>
</math>
</center>
</center>
Linia 996: Linia 996:


<center>
<center>
<math> |P|
<math>|P|
=
=
\int\limits_{\alpha}^{\beta}\psi(t)\varphi'(t)\,dt.
\int\limits_{\alpha}^{\beta}\psi(t)\varphi'(t)\,dt.
Linia 1013: Linia 1013:
<span id="twierdzenie_15_21">{{twierdzenie|15.21.||
<span id="twierdzenie_15_21">{{twierdzenie|15.21.||


Jeśli obszar jest ograniczony odcinkami <math> OA</math> i <math> OB</math>
Jeśli obszar jest ograniczony odcinkami <math>OA</math> i <math>OB</math>
(gdzie <math> O=(0,0)</math>)
(gdzie <math>O=(0,0)</math>)
oraz krzywą
oraz krzywą
<math> AB</math> daną w postaci biegunowej
<math>AB</math> daną w postaci biegunowej


<center>
<center>
<math> r
<math>r
=
=
g(\vartheta),
g(\vartheta),
Linia 1030: Linia 1030:


<center>
<center>
<math> |P|
<math>|P|
=
=
\frac{1}{2}\int\limits_{\vartheta_1}^{\vartheta_2}\big[g(\vartheta)\big]^2\,d\vartheta.
\frac{1}{2}\int\limits_{\vartheta_1}^{\vartheta_2}\big[g(\vartheta)\big]^2\,d\vartheta.
Linia 1039: Linia 1039:
Obszar dzielimy na trójkąty krzywoliniowe jak na rysunku
Obszar dzielimy na trójkąty krzywoliniowe jak na rysunku


Oznaczając przez <math> P_{ABC}</math> pole trójkąta krzywoliniowego, mamy
Oznaczając przez <math>P_{ABC}</math> pole trójkąta krzywoliniowego, mamy


<center>
<center>
<math> P_{ABC} \approx
<math>P_{ABC} \approx
\frac{1}{2}\cdot g(\vartheta)\cdot g(\vartheta)\cdot\sin\Delta\vartheta
\frac{1}{2}\cdot g(\vartheta)\cdot g(\vartheta)\cdot\sin\Delta\vartheta
\approx
\approx
Linia 1049: Linia 1049:
</center>
</center>


(dla małych kątów <math> \Delta\vartheta</math> zachodzi
(dla małych kątów <math>\Delta\vartheta</math> zachodzi
<math> \Delta\sin\vartheta\approx\Delta\vartheta</math>).
<math>\Delta\sin\vartheta\approx\Delta\vartheta</math>).
Sumując pola trójkątów (analogicznie jak
Sumując pola trójkątów (analogicznie jak
sumy całkowe w całce Riemanna;
sumy całkowe w całce Riemanna;
Linia 1064: Linia 1064:


<center>
<center>
<math> K:\ y=f(x),
<math>K:\ y=f(x),
\quad</math> dla <math> \ x\in[a,b]
\quad</math> dla <math> \ x\in[a,b]
</math>
</math>
</center>
</center>


wokół osi <math> Ox</math>:
wokół osi <math>Ox</math>:


<center>
<center>
<math> |P|
<math>|P|
=
=
2\pi
2\pi
Linia 1087: Linia 1087:


<center>
<center>
<math> K:
<math>K:
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 1094: Linia 1094:
\end{array}  
\end{array}  
\right.
\right.
\quad</math> dla <math> \ t\in[\alpha,\beta]
\quad</math> dla <math> \ t\in[\alpha,\beta]
</math>
</math>
</center>
</center>


wokół osi <math> Ox</math>:
wokół osi <math>Ox</math>:


<center>
<center>
<math> |P|
<math>|P|
=
=
2\pi
2\pi
Linia 1139: Linia 1139:


<center>
<center>
<math> K:\ y=f(x),
<math>K:\ y=f(x),
\quad</math> dla <math> \ x\in[a,b]
\quad</math> dla <math> \ x\in[a,b]
</math>
</math>
</center>
</center>


wokół osi <math> Ox</math>:
wokół osi <math>Ox</math>:


<center>
<center>
<math> |V_x|
<math>|V_x|
=
=
\pi
\pi
Linia 1159: Linia 1159:


<center>
<center>
<math> P:
<math>P:
a
a
=
=
Linia 1177: Linia 1177:
to znaczy na bryły powstałe przez
to znaczy na bryły powstałe przez
obrót obszaru pod wykresem funkcji
obrót obszaru pod wykresem funkcji
<math> y=f(x)</math> dla <math> x\in[x_{i-1},x_i]</math>.
<math>y=f(x)</math> dla <math>x\in[x_{i-1},x_i]</math>.
Objętość takiego "plasterka" jest w przybliżeniu równa
Objętość takiego "plasterka" jest w przybliżeniu równa
objętości walca o promieniu podstawy <math> f(x_i)</math> i wysokości
objętości walca o promieniu podstawy <math>f(x_i)</math> i wysokości
<math> \Delta x_i=x_i-x_{i-1}</math>, czyli
<math>\Delta x_i=x_i-x_{i-1}</math>, czyli
<math> \pi f(x_i)^2\Delta x_i</math>.
<math>\pi f(x_i)^2\Delta x_i</math>.
Sumując objętości "plasterków", otrzymujemy
Sumując objętości "plasterków", otrzymujemy
sumę całkową jak w całce Riemanna
sumę całkową jak w całce Riemanna
Linia 1191: Linia 1191:


<center>
<center>
<math> K:
<math>K:
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 1198: Linia 1198:
\end{array}  
\end{array}  
\right.
\right.
\quad</math> dla <math> \ t\in[\alpha,\beta]
\quad</math> dla <math> \ t\in[\alpha,\beta]
</math>
</math>
</center>
</center>


wokół osi <math> Ox</math>:
wokół osi <math>Ox</math>:


<center>
<center>
<math> |V_x|
<math>|V_x|
=
=
\pi
\pi
Linia 1242: Linia 1242:


<center>
<center>
<math> K:\ y=f(x)
<math>K:\ y=f(x)
\quad</math> dla <math>   x\in[a,b]
\quad</math> dla <math> x\in[a,b]
</math>
</math>
</center>
</center>


wokół osi <math> Oy</math>:
wokół osi <math>Oy</math>:


<center>
<center>
<math> |V_y|
<math>|V_y|
=
=
2\pi
2\pi
Linia 1262: Linia 1262:


<center>
<center>
<math> P:
<math>P:
a
a
=
=
Linia 1280: Linia 1280:
powstałe przez
powstałe przez
obrót obszaru pod wykresem funkcji
obrót obszaru pod wykresem funkcji
<math> y=f(x)</math> dla <math> x\in[x_{i-1},x_i]</math> wokół osi <math> Oy</math>.
<math>y=f(x)</math> dla <math>x\in[x_{i-1},x_i]</math> wokół osi <math>Oy</math>.
Objętość takiego "cylindra" jest w przybliżeniu równa
Objętość takiego "cylindra" jest w przybliżeniu równa
<math> 2\pi x_i f(x_i)-2\pi x_{i-1}f(x_i)=2\pi\Delta x_i f(x_i)</math>.
<math>2\pi x_i f(x_i)-2\pi x_{i-1}f(x_i)=2\pi\Delta x_i f(x_i)</math>.
Sumując objętości "cylindrów", otrzymujemy
Sumując objętości "cylindrów", otrzymujemy
sumę całkową jak w całce Riemanna
sumę całkową jak w całce Riemanna
i przechodząc do granicy, dostajemy wzór
i przechodząc do granicy, dostajemy wzór
na <math> |V_y|</math>.<br>
na <math>|V_y|</math>.<br>
<br>
<br>
'''(2)'''
'''(2)'''
Linia 1293: Linia 1293:


<center>
<center>
<math> K:
<math>K:
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 1300: Linia 1300:
\end{array}  
\end{array}  
\right.
\right.
\quad</math> dla <math> \ t\in[\alpha,\beta]
\quad</math> dla <math> \ t\in[\alpha,\beta]
</math>
</math>
</center>
</center>


wokół osi <math> Oy</math>:
wokół osi <math>Oy</math>:


<center>
<center>
<math> |V_y|
<math>|V_y|
=
=
2\pi
2\pi
Linia 1336: Linia 1336:


<center>
<center>
<math> x^2+(y-a)^2
<math>x^2+(y-a)^2
\le
\le
r^2
r^2
Linia 1344: Linia 1344:
</center>
</center>


wokół osi <math> Ox</math>.
wokół osi <math>Ox</math>.


<center>
<center>
Linia 1372: Linia 1372:


<center>
<center>
<math> \begin{align}  
<math>\begin{align}  
(\bigstar)\quad
(\bigstar)\quad
I
I
Linia 1386: Linia 1386:


<center>
<center>
<math> I_1
<math>I_1
=
=
\arcsin\frac{x}{|r|}+c.
\arcsin\frac{x}{|r|}+c.
Linia 1392: Linia 1392:
</center>
</center>


Teraz liczymy całkę <math> I</math> inaczej:
Teraz liczymy całkę <math>I</math> inaczej:


<center>
<center>
Linia 1410: Linia 1410:


<center>
<center>
<math> r^2I_1-I_2
<math>r^2I_1-I_2
=
=
x\sqrt{r^2-x^2}+I_2,
x\sqrt{r^2-x^2}+I_2,
Linia 1419: Linia 1419:


<center>
<center>
<math> 2I_2
<math>2I_2
=
=
r^2I_1-x\sqrt{r^2-x^2}
r^2I_1-x\sqrt{r^2-x^2}
Linia 1431: Linia 1431:


<center>
<center>
<math> I_2
<math>I_2
=
=
\frac{r^2}{2}\arcsin\frac{x}{r}
\frac{r^2}{2}\arcsin\frac{x}{r}
Linia 1441: Linia 1441:


<center>
<center>
<math> \begin{align}  
<math>\begin{align}  
I
I
& = &
& = &

Wersja z 10:35, 5 wrz 2023

Krzywe i bryły obrotowe

W tym wykładzie wprowadzamy pojęcie krzywej i krzywej zwyczajnej. Definiujemy długość krzywej i krzywą prostowalną. Dowodzimy, że krzywa zwyczajna klasy C1 jest prostowalna. Wyprowadzamy wzór na długość krzywej i liczymy długości cykloidy i asteroidy. W drugiej części wykładu podajemy wzory na pola powierzchni i objętości brył obrotowych.

Długość krzywej

Definicja 15.1.

Niech <a<b<+. Krzywą nazywamy zbiór punktów

K={(x,y)2: x=φ(t), y=ψ(t), t[a,b]},

gdzie φ,ψ:[a,b] są dwiema funkcjami ciągłymi. Piszemy:

K=K(φ,ψ):{x=φ(t)y=ψ(t)t[a,b].

Powyższe równanie nazywamy też równaniem parametrycznym krzywej.

<flash>file=AM1.M15.W.R01.swf|width=375|height=375</flash>

<div.thumbcaption>Krzywa

<flash>file=Am1.M15.W.R02.swf|width=375|height=375</flash>

<div.thumbcaption>Parametryczny opis okręgu
Plik:AM1.M15.W.R03.mp4
Krzywa z punktem wielokrotnym (potrójnym)

Przykład 15.2.

Zapiszmy parametryczne równanie okręgu o promieniu R>0 w 2 Jeśli jako parametr t przyjmiemy kąt jaki tworzy promień poprowadzony do punktu (x,y) na okręgu, to widzimy (patrz rysunek), że x=cost i y=sint Zatem następująca krzywa:

K:{x=Rcosty=Rsintt[0,2π] opisuje okrąg.

Definicja 15.3.

Mówimy, że punkt (x,y)K jest punktem wielokrotnym krzywej K, jeśli


t1,t2(a,b):t1t2(x,y)=(φ(t1),ψ(t1))=(φ(t2),ψ(t2)).


Krzywą K nazywamy zwyczajną, jeśli nie zawiera punktów wielokrotnych, to znaczy


[(φ(t1),ψ(t1))=(φ(t2),ψ(t2)), t1t2][(t1=t2)(t1=a  t2=b)].

Plik:AM1.M15.W.R04.svg
Krzywe zwyczajne

Definicja 15.4.

Niech

a=t0<t1<<tn=b


będzie podziałem przedziału [a,b]. Łamaną p łączącą punkty:


(φ(t0),ψ(t0)), ,(φ(tn),ψ(tn))


nazywamy łamaną wpisaną w krzywą K. Przez l(p) oznaczamy długość łamanej p (to znaczy sumę długości odcinków

wchodzących w skład łamanej).

Definicja 15.5.

Długością krzywej K nazywamy liczbę:

l(K)=suppl(p),

gdzie supremum jest brane po wszystkich łamanych wpisanych w K.


<flashwrap>file=AM1.M15.W.R05.swf|size=small</flashwrap>

<div.thumbcaption>Łamana wpisana w krzywą

<flashwrap>file=AM1.M15.W.R06.swf|size=small</flashwrap>

<div.thumbcaption>Łamana wpisana w krzywą

Definicja 15.6.

Jeśli l(K)<+, to mówimy, że krzywa K jest prostowalna.

Twierdzenie 15.7.

Niech φ,ψ:[a,b] będą klasy C1 oraz niech K=K(φ,ψ) będzie krzywą zwyczajną.
Wówczas krzywa K jest prostowalna.

Plik:AM1.M15.W.R07.mp4
Łamana wpisana w krzywą

Dowód 15.7. [nadobowiązkowy]

Niech p będzie dowolną łamaną wpisaną w krzywą K, to znaczy istnieje podział

a=t0<t1<<tn=b

taki, że p jest łamaną o wierzchołkach (xi,yi) dla i=0,,n, gdzie

{xi=φ(ti)yi=ψ(ti)i{0,,n}.

Długość łamanej p wyraża się wzorem:

l(p)=i=1n(xixi1)2+(yiyi1)2.

Ponieważ φ,ψC1([a,b];), więc z twierdzenia o wartości średniej (patrz twierdzenie 9.37.) mamy

xixi1=φ(ti)φ(ti1)=φ(τi)(titi1),

yiyi1=ψ(ti)ψ(ti1)=ψ(τi*)(titi1),

gdzie

τi(ti1,ti),i=1,n,τi*(ti1,ti),i=1,n.

Zatem

l(p)=i=1nφ(τi)2+ψ(τi*)2(titi1).

Ponieważ φ,ψC([a,b];) i przedział [a,b] jest zwarty, więc funkcje φ,ψ są ograniczone.
Definiujemy

M=supt[a,b]φ(t),M*=supt[a,b]ψ(t),

m=inft[a,b]φ(t),m*=inft[a,b]ψ(t).

Zatem

m2+m*2(ba)l(p)M2+M*2(ba).

Ponieważ powyższa nierówność zachodzi dla dowolnej łamanej p wpisanej w krzywą K, więc przechodząc do supremum po wszystkich takich łamanych, dostajemy

m2+m*2(ba)l(K)M2+M*2(ba),

a zatem krzywa K jest prostowalna.

Uwaga 15.8.

W powyższym twierdzeniu (jak i twierdzeniach następnych) zakładamy, że krzywa jest klasy C1. (to znaczy φ,ψ, są klasy C1) W zastosowaniach okaże się jednak, że często będziemy mieli do czynienia z krzywymi, które są tylko ciągłe, zwyczajne oraz "kawałkami" klasy C1, to znaczy krzywą można otrzymać jako "sklejenie" kilku krzywych klasy C1 (przy sklejaniu początek następnej krzywej jest końcem poprzedniej). Wszystkie wypowiadane tu twierdzenia dla krzywych klasy C1 stosują się także do krzywych kawałkami klasy C1.

Plik:AM1.M15.W.R08.mp4
Krzywa K(t)

Definicja 15.9.

Niech K=K(φ,ψ) będzie krzywą. Zdefiniujmy:

K(t)  =df{(φ(τ),ψ(τ)): τ[a,t]},

oraz

s(t)  =dfl(K(t)) (długośćkrzywejK(t)) .

W szczególności s(b)=l(K).

Twierdzenie 15.10.

Niech φ,ψ:[a,b] będą klasy C1 oraz niech K=K(φ,ψ) będzie krzywą zwyczajną.
Wówczas

s(t)=φ(t)2+ψ(t)2 t[a,b].

Dowód 15.10. [nadobowiązkowy]

Niech t0,t0+h[a,b]. Analogicznie do ostatniego oszacowania w dowodzie twierdzenia 15.7. dostajemy:

mh2+mh*2hs(t0+h)s(t0)Mh2+Mh*2h,

gdzie

Mh=supt[t0,t0+h]φ(t),Mh*=supt[t0,t0+h]ψ(t),

mh=inft[t0,t0+h]φ(t),mh*=inft[t0,t0+h]ψ(t).

Dzielimy wszystkie strony powyższego oszacowania przez h, dostając:

mh2+mh*2s(t0+h)s(t0)hMh2+Mh*2.

Ponieważ funkcje φ i ψ są ciągłe, więc dostajemy

Mhh0φ(t0),mhh0φ(t0),Mh*h0ψ(t0),mh*h0ψ(t0).

Z powyższych oszacowań oraz z twierdzenia o trzech ciągach dostajemy:

s(t0)=limh0s(t0+h)s(t0)h=φ(t0)2+ψ(t0)2.

Twierdzenie 15.11. [O długości krzywej]

Niech φ,ψ:[a,b] będą klasy C1 oraz niech K=K(φ,ψ) będzie krzywą zwyczajną. Wówczas długość krzywej wyraża się wzorem

l(K)=abφ(τ)2+ψ(τ)2dτ.

W szczególności, jeśli krzywa zadana jest wykresem funkcji y=f(x), dla x[a,b], to

l(K)=ab1+f(t)2dt.

Dowód 15.11.

l(K)=s(b)=s(b)s(a)=0=abs(τ)dτ=abφ(τ)2+ψ(τ)2dτ.

W drugim przypadku krzywą zadaną przez funkcję f możemy zapisać w postaci parametrycznej

{x(t)=ty(t)=f(t),t[a,b]

i od razu otrzymujemy drugi ze wzorów.

Przykład 15.12.

Wyprowadzić wzór na długość krzywej zadanej w postaci biegunowej:

r=g(ϑ)ϑ[α,β].

Przedstawmy tę krzywą w postaci parametrycznej:

{x = rcosϑ=g(ϑ)cosϑy = rsinϑ=g(ϑ)sinϑ.

Liczymy

x(ϑ)2+y(ϑ)2=[g(ϑ)cosϑg(ϑ)sin(ϑ)]2+[g(ϑ)sinϑ+g(ϑ)cos(ϑ)]2=g(ϑ)2cos2ϑ2g(ϑ)g(ϑ)sinϑcosϑ+g(ϑ)2sin2ϑ+g(ϑ)2sin2ϑ+2g(ϑ)g(ϑ)sinϑcosϑ+g(ϑ)2cos2ϑ=g(ϑ)2+g(ϑ)2,

Zatem

l(K)=αβg(ϑ)2+g(ϑ)2dϑ.

<flashwrap>file=AM1.M15.W.R09.swf|size=small</flashwrap>

<div.thumbcaption>Krzywa we współrzędnych biegunowych

<flashwrap>file=AM1.M15.W.R10.swf|size=small</flashwrap>

<div.thumbcaption>Cykloida


Definicja 15.13.

Cykloidą nazywamy krzywą kreśloną

przez ustalony punkt 0 na okręgu toczącym się po prostej l.
Plik:AM1.M15.W.R11.svg
Cykloida

Przykład 15.14.

Wyprowadzić wzór parametryczny cykloidy.

Oznaczenia:
a - promień okręgu;
O - początkowy punkt styczności okręgu i prostej l;
N - nowy punkt styczności;
M - nowe położenie punktu O;
t=NDM - parametr określający położenie punktu M.

Liczymy współrzędne punktu M(x,y):

x =OF=ONFN=NM^MG=atasint,


y =FM=NG=NDGD=aacost.

Zatem

{x=a(tsint)y=a(1cost)t[0,2π]( lub  t).

Przykład 15.15.

Obliczyć długość łuku cykloidy:

{x=a(tsint)y=a(1cost)t[0,2π].


x(t)2+y(t)2=a2(1cost)2+a2sin2t=a22a2cost+a2cos2t+a2sin2t=2a2(1cost)=4a2sin2t2=2a|sint2|.

Zatem

l(K)=02πx(t)2+y(t)2dt=2a02π|sint2|dt=2a02πsint2dt=4acost2|02π=8a.

Przykład 15.16.

Obliczyć długość łuku asteroidy:

x23+y23=a23.

Równanie parametryczne asteroidy, to:

{x=acos3ty=asin3tt[0,2π].

Liczymy

x(t)2+y(t)2=3asintcost t[0,2π].

Zatem

l(K)=40π23asintcostdt=6a.

<flash>file=AM1.M15.W.R12.swf|width=272|height=272</flash>

<div.thumbcaption>Asteroida

<flashwrap>file=AM1.M15.W.R13.swf|size=small</flashwrap>

<div.thumbcaption>Asteroida

<flashwrap>file=AM1.M15.W.R14.swf|size=small</flashwrap>

<div.thumbcaption>Asteroida

Całka krzywoliniowa

Niech K będzie krzywą klasy C1:

K={(x,y)2: x=φ(t), y=ψ(t), t[a,b]},

Przypuśćmy, że określona jest funkcja ciągła f:KMf(M), to znaczy funkcja, która każdemu punktowi M krzywej K przyporządkowuje pewną wartość rzeczywistą f(M). Okazuje się, że dla takich funkcji możemy także zdefiniować całkę oznaczoną, to znaczy całkę z funkcji f po krzywej K.

Georg Friedrich Bernhard Riemann (1826-1866)
Zobacz biografię

Całkę tę wprowadza się analogicznie jak całkę Riemanna na odcinku. Pominiemy to jednak w tym miejscu, podając jedynie wzór końcowy na obliczanie takiej całki:

Kf(x,y)ds=abf(φ(t),ψ(t))φ(t)2+ψ(t)2dt.

Tę całkę stosuje się w fizyce na przykład do obliczania masy i środka ciężkości krzywej (pręta, którego wszystkie wymiary poza długością są pomijalne).

Jeśli mamy daną krzywą (pręt) K zadaną jak wyżej, o gęstości w każdym jej punkcie M(x,y) danej funkcją ciągłą ϱ(M), to masa tego pręta wyraża się wzorem

m=Kϱ(x,y)ds.

Współrzędne środka ciężkości pręta (x0,y0) możemy policzyć ze wzorów

x0=1mKxϱ(x,y)ds,x0=1mKyϱ(x,y)ds.

Przykład 15.17.

Obliczyć masę pręta półkolistego K={(x,y)2: x2+y2=R2, y0} o gęstości ϱ(x,y)=y2.

Masa krzywej o gęstości ϱ dana jest wzorem

m=Kϱ(x,y)ds=Ky2ds.

Stosując wzór na całkę krzywoliniową oraz korzystając z parametryzacji półokręgu:

K=K(φ,ψ):{x=φ(t)=Rcosty=ψ(t)=Rsintt[0,π],

mamy

m=0πR2sin2t(Rsint)2+(Rcost)2dt=R30πsin2tdt=R3[t214sin2t]0π=R3π2.

Odpowiedź:

Masa pręta wynosi R3π2

.

Przykład 15.18.

Obliczyć masę i współrzędne środka ciężkości odcinka K łączącego punkt (0,0) z punktem (1,1) o gęstości wprost proporcjonalnej do odległości punktu od środka układu i równej 2 w punkcie (1,1).

Skoro gęstość ϱ jest proporcjonalna do odległości punktu od środka układu i wynosi 2 w punkcie (1,1), to

ϱ(x,t)=cx2+y2 oraz ϱ(1,1)=c2=2,

stąd c=1. Parametryzacją odcinka jest na przykład

K=K(φ,ψ):{x=φ(t)=ty=ψ(t)=tt[0,1],

zatem masa wynosi

m=Kx2+y2ds=01t2+t22dt=201tdt=t2|01=1.

Pierwszą współrzędną środka ciężkości liczymy ze wzoru

x0=1mKxϱ(x,y)ds=01t2t22dt=201t2dt=23t3|01=23.

Z symetrii zadania wynika, że y0=23.

Pole powierzchni i objętość bryły obrotowej

W dalszej części wykładu będziemy zakładać, że krzywe są klasy C1. Podamy liczne wzory na obliczanie pól powierzchni i objętości brył obrotowych, w większości pozostawiając je bez dowodów (podając natomiast pewne ich uzasadnienia).

Plik:AM1.M15.W.R15.svg
Pole między wykresami funkcji

Z poprzedniego wykładu znamy już związek całki Riemanna z polem obszaru ograniczonego wykresami funkcji. Dla porządku przypomnijmy ten związek.

Uwaga 15.19.

Jeśli trapez krzywoliniowy jest ograniczony z góry i z dołu krzywymi:

y=f1(x) i y=f2(x)x[a,b],

to pole tego trapezu wynosi:

|P|=ab[f1(x)f2(x)]dx

Uzasadnienie: Wzór ten wynika bezpośrednio z geometrycznej interpretacji całki oznaczonej.

Twierdzenie 15.20.

Pole obszaru pod wykresem krzywej zadanej w postaci parametrycznej

K:{x=φ(t)y=ψ(t), dla  t[α,β],

wynosi

|P|=αβψ(t)φ(t)dt.

Uzasadnienie: Wzór ten jest konsekwencją wzoru z uwagi 15.19. i twierdzenia o całkowaniu przez podstawienie.

Plik:AM1.M15.W.R16.mp4
Pole obszaru ograniczonego półprostymi i krzywą zadaną we współrzędnych biegunowych
Plik:AM1.M15.W.R17.svg
Trójkąt krzywoliniowy

Twierdzenie 15.21.

Jeśli obszar jest ograniczony odcinkami OA i OB (gdzie O=(0,0)) oraz krzywą AB daną w postaci biegunowej

r=g(ϑ),ϑ[ϑ1,ϑ2],

to pole tego obszaru wynosi:

|P|=12ϑ1ϑ2[g(ϑ)]2dϑ.

Uzasadnienie: Obszar dzielimy na trójkąty krzywoliniowe jak na rysunku

Oznaczając przez PABC pole trójkąta krzywoliniowego, mamy

PABC12g(ϑ)g(ϑ)sinΔϑ12g(ϑ)2Δϑ

(dla małych kątów Δϑ zachodzi ΔsinϑΔϑ). Sumując pola trójkątów (analogicznie jak sumy całkowe w całce Riemanna; patrz definicja 14.4.) i przechodząc do granicy, dostajemy powyższy wzór.

Twierdzenie 15.22.

(1) Pole powierzchni powstałej z obrotu krzywej

K: y=f(x), dla  x[a,b]

wokół osi Ox:

|P|=2πab[f(x)]1+f(x)2dx.

Wzór ten pozostawiamy bez uzasadnienia.

(2) Pole powierzchni powstałej z obrotu krzywej

K:{x=φ(t)y=ψ(t) dla  t[α,β]

wokół osi Ox:

|P|=2παβ[ψ(t)]φ(t)2+ψ(t)2dt.

<flashwrap>file=AM1.M15.W.R18.swf|size=small</flashwrap>

<div.thumbcaption>Powierzchnia powstała przez obrót krzywej dookoła osi Ox

<flashwrap>file=AM1.M15.W.R21.swf|size=small</flashwrap>

<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi Ox

<flashwrap>file=AM1.M15.W.R19.swf|size=small</flashwrap>

<div.thumbcaption>Powierzchnia powstała przez obrót krzywej dookoła osi Ox

<flashwrap>file=AM1.M15.W.R20.swf|size=small</flashwrap>

<div.thumbcaption>Powierzchnia powstała przez obrót krzywej dookoła osi Ox

Twierdzenie 15.23.

(1) Objętość bryły powstałej z obrotu obszaru "pod krzywą"

K: y=f(x), dla  x[a,b]

wokół osi Ox:

|Vx|=πabf(x)2dx.

Uzasadnienie: Weźmy podział odcinka [a,b]:

P:a=x0<x1<<xn=b

oraz podzielmy bryłę na "plasterki", to znaczy na bryły powstałe przez obrót obszaru pod wykresem funkcji y=f(x) dla x[xi1,xi]. Objętość takiego "plasterka" jest w przybliżeniu równa objętości walca o promieniu podstawy f(xi) i wysokości Δxi=xixi1, czyli πf(xi)2Δxi. Sumując objętości "plasterków", otrzymujemy sumę całkową jak w całce Riemanna i przechodząc do granicy, dostajemy powyższy wzór.

(2) Objętość bryły powstałej z obrotu obszaru "pod krzywą"

K:{x=φ(t)y=ψ(t) dla  t[α,β]

wokół osi Ox:

|Vx|=παβψ(t)2φ(t)dt.

Uzasadnienie: Wzór powyższy jest konsekwencją poprzedniego wzoru oraz twierdzenia o całkowaniu przez podstawienie.

<flashwrap>file=AM1.M15.W.R22.swf|size=small</flashwrap>

<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi Ox

<flashwrap>file=AM1.M15.W.R23.swf|size=small</flashwrap>

<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi Ox

<flashwrap>file=AM1.M15.W.R24.swf|size=small</flashwrap>

<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi Oy

Twierdzenie 15.24.

(1) Objętość bryły powstałej z obrotu obszaru "pod krzywą"

K: y=f(x) dla x[a,b]

wokół osi Oy:

|Vy|=2πabxf(x)dx.

Uzasadnienie: Weźmy podział odcinka [a,b]:

P:a=x0<x1<<xn=b

oraz podzielmy bryłę na "cylindry" powstałe przez obrót obszaru pod wykresem funkcji y=f(x) dla x[xi1,xi] wokół osi Oy. Objętość takiego "cylindra" jest w przybliżeniu równa 2πxif(xi)2πxi1f(xi)=2πΔxif(xi). Sumując objętości "cylindrów", otrzymujemy sumę całkową jak w całce Riemanna i przechodząc do granicy, dostajemy wzór na |Vy|.

(2) Objętość bryły powstałej z obrotu obszaru "pod krzywą"

K:{x=φ(t)y=ψ(t) dla  t[α,β]

wokół osi Oy:

|Vy|=2παβφ(t)ψ(t)φ(t)dt.

<flashwrap>file=AM1.M15.W.R25.swf|size=small</flashwrap>

<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi Ox

<flashwrap>file=AM1.M15.W.R26.swf|size=small</flashwrap>

<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi Ox
Plik:AM1.M15.W.R27.mp4
Torus

Przykład 15.25.

Obliczyć objętość torusa powstałego przez obrót koła

x2+(ya)2r2(0<r<a)

wokół osi Ox.

|Vx|=πrr[(a+r2x2)2(ar2x2)2]dx=4πarrr2x2dx=()4πa[r22arcsinxr+x2r2x2]rr=4πa[r22π2+r22π2]=4πar2π2=2π2ar2,

gdzie wykorzystano następującą całkę:

()I=r2x2dx=r2x2r2x2dx=r2dxr2x2I1x2dxr2x2I2.

I1=arcsinx|r|+c.

Teraz liczymy całkę I inaczej:

I=r2x2dx części=xr2x2x2x2r2x2dx=xr2x2+x2r2x2dx=I2=xr2x2+I2.

Porównując to z (), otrzymujemy:

r2I1I2=xr2x2+I2,

stąd

2I2=r2I1xr2x2=r2arcsinxrxr2x2,

zatem

I2=r22arcsinxrx2r2x2.

Wstawiając do (), otrzymujemy:

I=r2arcsinxr12r2arcsinxr+12xr2x2+c=r22arcsinxr+x2r2x2+c.