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

Z Studia Informatyczne
Przejdź do nawigacjiPrzejdź do wyszukiwania
Gracja (dyskusja | edycje)
 
(Nie pokazano 73 wersji utworzonych przez 8 użytkowników)
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> \displaystyle  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 11: Linia 11:


==Długość krzywej==
==Długość krzywej==
{{definicja|15.1.||
 
Niech <math> \displaystyle  -\infty<a<b<+\infty.</math>
{{definicja|15.1.|definicja_15_1|
Niech <math>-\infty<a<b<+\infty</math>.
'''''Krzywą''''' nazywamy
'''''Krzywą''''' nazywamy
zbiór punktów
zbiór punktów


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


gdzie <math> \displaystyle  \displaystyle\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><math> \displaystyle 
<center>
K=K(\varphi,\psi):\
<math>
\left\{
K=K(\varphi,\psi):
\begin{array} {l}
\left\{ \begin{array} {l}
x=\varphi(t)\\
x=\varphi(t)\\
y=\psi(t)
y=\psi(t)
\end{array}  
\end{array}  
\right.
\right.
\qquad t\in[a,b].
\qquad t\in[a,b]
</math></center>
</math>
</center>


Powyższe równanie nazywamy też
Powyższe równanie nazywamy też '''''równaniem parametrycznym krzywej'''''.<br>}}
'''''równaniem parametrycznym krzywej'''''.<br>}}


{| border="0" align="center" cellspacing="10"
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:375px;">
|<div class="thumb"><div style="width:375px;">
<flash>file=AM1.M15.W.R01.swf|width=375|height=375</flash>
<flash>file=AM1.M15.W.R01.swf|width=375|height=375</flash>
<div.thumbcaption>AM1.M15.W.R01</div>
<div.thumbcaption>Krzywa</div>
</div></div>
</div></div>
|<div class="thumb"><div style="width:375px;">
|<div class="thumb"><div style="width:375px;">
<flash>file=Am1.M15.W.R02.swf|width=375|height=375</flash>
<flash>file=Am1.M15.W.R02.swf|width=375|height=375</flash>
<div.thumbcaption>Am1.M15.W.R02.swf</div>
<div.thumbcaption>Parametryczny opis okręgu</div>
</div></div>
</div></div>
|}
|}


<div class="thumb tright"><div style="width:260px;height:240;">
[[File:AM1.M15.W.R03.mp4|253x253px|thumb|right|Krzywa z punktem wielokrotnym (potrójnym)]]
<flashwrap>file=AM1.M15.W.R03.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R03</div>
</div></div>


<span id="przyklad_15_2">{{przyklad|15.2.||
<span id="przyklad_15_2">{{przyklad|15.2.||


Zapiszmy parametryczne równanie okręgu o
Zapiszmy parametryczne równanie okręgu o
promieniu <math> \displaystyle  R>0</math> w <math> \displaystyle  \displaystyle\mathbb{R}^2.</math>
promieniu <math>R>0</math> w <math>\mathbb{R}^2</math>
Jeśli jako parametr <math> \displaystyle  t</math> przyjmiemy kąt jaki tworzy promień
Jeśli jako parametr <math>t</math> przyjmiemy kąt jaki tworzy promień
poprowadzony do punktu <math> \displaystyle  \displaystyle (x,y)</math> na okręgu, to łatwo widzimy
poprowadzony do punktu <math>(x,y)</math> na okręgu, to widzimy
(patrz rysunek), że
(patrz rysunek), że
<math> \displaystyle  x=\cos t</math> i <math> \displaystyle  y=\sin t.</math>
<math>x=\cos t</math> i <math>y=\sin t</math>
Zatem następująca krzywa:
Zatem następująca krzywa:


<math> \displaystyle 
<math>
K:\
K:
\left\{
\left\{ \begin{array} {l}
\begin{array} {l}
x=R\cos t\\
x=R\cos t\\
y=R\sin t
y=R\sin t
Linia 77: Linia 76:
}}</span>
}}</span>


{{definicja|15.3.||
{{definicja|15.3.|definicja_15_3|


Mówimy, że punkt <math> \displaystyle  \displaystyle (x,y)\in K</math> jest
Mówimy, że punkt <math>(x,y)\in K</math> jest
'''''punktem wielokrotnym''''' krzywej <math> \displaystyle  K,</math>
'''''punktem wielokrotnym''''' krzywej <math>K</math>,
jeśli
jeśli


<br>
<center>
<center>
<math> \displaystyle  \exists t_1,t_2\in(a,b):\
<math>\exists t_1,t_2\in(a,b):
t_1\ne t_2\quad\land\quad
t_1\ne t_2\quad\land\quad
(x,y)=\big(\varphi(t_1),\psi(t_1)\big)=\big(\varphi(t_2),\psi(t_2)\big).
(x,y)=\big(\varphi(t_1),\psi(t_1)\big)=\big(\varphi(t_2),\psi(t_2)\big)
</math>
</math>
</center>
</center>
 
<br>
Krzywą <math> \displaystyle  K</math> nazywamy '''''zwyczajną''''',
Krzywą <math>K</math> nazywamy '''''zwyczajną''''',
jeśli nie zawiera punktów wielokrotnych, to znaczy
jeśli nie zawiera punktów wielokrotnych, to znaczy


<br>
<center>
<center>
<math>\begin{array}{ll} \displaystyle &  \bigg[
<math>\begin{array}{ll}\bigg[
\big(\varphi(t_1),\psi(t_1)\big)=\big(\varphi(t_2),\psi(t_2)\big),\ t_1\le t_2
\big(\varphi(t_1),\psi(t_1)\big)=\big(\varphi(t_2),\psi(t_2)\big),\ t_1\le t_2
\bigg]\\\\
\bigg]\\\\
Linia 105: Linia 106:
</center>}}
</center>}}


<div class="thumb tright"><div style="width:250px;">
[[File:AM1.M15.W.R04.svg|250x250px|thumb|right|Krzywe zwyczajne]]
<flash>file=AM1.M15.W.R04.swf|width=250|height=250</flash>
<div.thumbcaption>AM1.M15.W.R04</div>
</div></div>


{{definicja|15.4.||
{{definicja|15.4.||
Linia 115: Linia 113:


<center>
<center>
<math> \displaystyle  a
<math>a
\ =\
=
t_0
t_0
\ <\
<
t_1
t_1
\ <\
<
\ldots
\ldots
\ <\
<
t_n
t_n
\ =\
=
b
b
</math>
</math>
</center>
</center><br>


będzie podziałem przedziału <math> \displaystyle  \displaystyle [a,b].</math>
będzie podziałem przedziału <math>[a,b]</math>.
Łamaną <math> \displaystyle  p</math> łączącą punkty:
Łamaną <math>p</math> łączącą punkty:


<center><math> \displaystyle  \big(\varphi(t_0),\psi(t_0)\big),
<br>
\ \ldots,\
<center>
<math>\big(\varphi(t_0),\psi(t_0)\big),
\ \ldots,
\big(\varphi(t_n),\psi(t_n)\big)
\big(\varphi(t_n),\psi(t_n)\big)
</math></center>
</math>
</center>
<br>


nazywamy
nazywamy
'''''łamaną wpisaną w krzywą <math> \displaystyle  K</math>'''''.
'''''łamaną wpisaną w krzywą <math>K</math>'''''.
Przez <math> \displaystyle  l(p)</math> oznaczamy '''''długość'''''
Przez <math>l(p)</math> oznaczamy '''''długość'''''
łamanej <math> \displaystyle  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> \displaystyle  K</math> nazywamy liczbę:
Długością krzywej <math>K</math> nazywamy liczbę:


<center><math> \displaystyle  l(K)
<center>
\ =\
<math>l(K)
\sup_p l(p),
=
</math></center>
\sup_p l(p)
</math>
</center>


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




{| border="0" align="center" cellspacing="10"
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R05.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R05.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R05</div>
<div.thumbcaption>Łamana wpisana w krzywą</div>
</div></div>
</div></div>
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R06.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R06.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R06</div>
<div.thumbcaption>Łamana wpisana w krzywą</div>
</div></div>
</div></div>
|}
|}
Linia 169: Linia 172:
{{definicja|15.6.||
{{definicja|15.6.||


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


Niech <math> \displaystyle  \displaystyle\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math> \displaystyle  C^1</math>
Niech <math>\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math>C^1</math>
oraz niech <math> \displaystyle  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> \displaystyle  K</math> jest prostowalna.
Wówczas krzywa <math>K</math> jest prostowalna.
}}</span>
}}</span>


<div class="thumb tright"><div style="width:260px;height:240;">
[[File:AM1.M15.W.R07.mp4|253x253px|thumb|right|Łamana wpisana w krzywą]]
<flashwrap>file=AM1.M15.W.R07.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R07</div>
</div></div>


{{dowod|twierdzenia 15.7.||
{{dowod|15.7. [nadobowiązkowy]||
(Dowód nadobowiązkowy.)<br>
Niech <math>p</math> będzie dowolną łamaną wpisaną w
Niech <math> \displaystyle  p</math> będzie dowolną łamaną wpisaną w
krzywą <math>K</math>,
krzywą <math> \displaystyle  K,</math>
to znaczy istnieje podział
to znaczy istnieje podział


<center><math> \displaystyle  a
<center>
\ =\
<math>a
=
t_0
t_0
\ <\
<
t_1
t_1
\ <\
<
\ldots
\ldots
\ <\
<
t_n
t_n
\ =\
=
b
b
</math></center>
</math>
</center>


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


<center><math> \displaystyle 
<center>
\left\{
<math>
\begin{array} {l}
\left\{ \begin{array} {l}
x_i\ =\ \varphi(t_i)\\
x_i= \varphi(t_i)\\
y_i\ =\ \psi(t_i)
y_i= \psi(t_i)
\end{array}  
\end{array}  
\right.
\right.
\qquad i\in\left\{0,\ldots,n\right\}.
\qquad i\in\left \{ 0,\ldots,n\right \}
</math></center>
</math>
</center>


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


<center><math> \displaystyle  l(p)
<center>
\ =\
<math>l(p)
=
\sum_{i=1}^n
\sum_{i=1}^n
\sqrt{\big(x_i-x_{i-1}\big)^2+\big(y_i-y_{i-1}\big)^2}.
\sqrt{\big(x_i-x_{i-1}\big)^2+\big(y_i-y_{i-1}\big)^2}
</math></center>
</math>
</center>


Ponieważ <math> \displaystyle  \displaystyle\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><math> \displaystyle  x_i-x_{i-1}
<center>
\ =\
<math>x_i-x_{i-1}
=
\varphi(t_i)-\varphi(t_{i-1})
\varphi(t_i)-\varphi(t_{i-1})
\ =\
=
\varphi'(\tau_i)\left(t_i-t_{i-1}\right),</math></center>
\varphi'(\tau_i)\left(t_i-t_{i-1}\right)</math>,
<center><math> \displaystyle  y_i-y_{i-1}
</center>
\ =\
<center>
<math>y_i-y_{i-1}
=
\psi(t_i)-\psi(t_{i-1})
\psi(t_i)-\psi(t_{i-1})
\ =\
=
\psi'(\tau^*_i)\left(t_i-t_{i-1}\right),</math></center>
\psi'(\tau^*_i)\left(t_i-t_{i-1}\right)</math>,
</center>
gdzie
gdzie
<center><math>\tau_i\in \left(t_{i-1},t_i\right),\quad i=1,\ldots n,\quad \tau^*_i\in \left(t_{i-1},t_i\right),\quad i=1,\ldots n.</math></center>
<center>
<math>\tau_i\in \left(t_{i-1},t_i\right),\quad i=1,\ldots n,\quad \tau^*_i\in \left(t_{i-1},t_i\right),\quad i=1,\ldots n</math>.
</center>


Zatem
Zatem


<center><math> \displaystyle  l(p)
<center>
\ =\
<math>l(p)
=
\sum_{i=1}^n
\sum_{i=1}^n
\sqrt{\varphi'(\tau_i)^2+\psi'(\tau^*_i)^2}\cdot\left(t_i-t_{i-1}\right).
\sqrt{\varphi'(\tau_i)^2+\psi'(\tau^*_i)^2}\cdot\left(t_i-t_{i-1}\right)
</math></center>
</math>
</center>


Ponieważ <math> \displaystyle  \displaystyle\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> \displaystyle  \displaystyle [a,b]</math> jest zwarty,
i przedział <math>[a,b]</math> jest zwarty,
więc funkcje <math> \displaystyle  \displaystyle\varphi',\psi'</math> są ograniczone.<br>
więc funkcje <math>\varphi',\psi'</math> są ograniczone.<br>
Definiujemy
Definiujemy


<center><math> \displaystyle  M \ =\
<center>
<math>M =
\sup_{t\in[a,b]}\varphi'(t),
\sup_{t\in[a,b]}\varphi'(t),
\qquad
\qquad
M^*
M^*
\ =\
=
\sup_{t\in[a,b]}\psi'(t),
\sup_{t\in[a,b]}\psi'(t)
</math></center>
</math>
</center>


<center><math> \displaystyle  m \ =\
<center>
\displaystyle\inf_{t\in[a,b]}\varphi'(t),
<math>m =
\inf_{t\in[a,b]}\varphi'(t),
\qquad
\qquad
m^* \ =\
m^* =
\displaystyle\inf_{t\in[a,b]}\psi'(t).
\inf_{t\in[a,b]}\psi'(t)</math>
</math></center>
</center>


Zatem
Zatem


<center><math> \displaystyle  \sqrt{m^2+{m^*}^2}\cdot(b-a)
<center>
\ \le\
<math>\sqrt{m^2+{m^*}^2}\cdot(b-a)
\le
l(p)
l(p)
\ \le\
\le
\sqrt{M^2+{M^*}^2}\cdot(b-a).
\sqrt{M^2+{M^*}^2}\cdot(b-a)
</math></center>
</math>
</center>


Ponieważ powyższa nierówność zachodzi dla dowolnej
Ponieważ powyższa nierówność zachodzi dla dowolnej
łamanej <math> \displaystyle  p</math> wpisanej w krzywą <math> \displaystyle  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><math> \displaystyle  \sqrt{m^2+{m^*}^2}\cdot(b-a)
<center>
\ \le\
<math>\sqrt{m^2+{m^*}^2}\cdot(b-a)
\le
l(K)
l(K)
\ \le\
\le
\sqrt{M^2+{M^*}^2}\cdot(b-a),
\sqrt{M^2+{M^*}^2}\cdot(b-a)
</math></center>
</math>
</center>


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


Linia 298: 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> \displaystyle  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>)
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> \displaystyle  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> \displaystyle  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> \displaystyle  C^1,</math>
Wszystkie wypowiadane tu twierdzenia dla krzywych klasy <math>C^1</math>
stosują się także do krzywych kawałkami klasy <math> \displaystyle  C^1.</math>
stosują się także do krzywych kawałkami klasy <math>C^1</math>.
}}
}}


<div class="thumb tright"><div style="width:260px;height:240;">
[[File:AM1.M15.W.R08.mp4|253x253px|thumb|right|Krzywa <math>K(t)</math>]]
<flashwrap>file=AM1.M15.W.R08.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R08</div>
</div></div>


{{definicja|15.9.||
{{definicja|15.9.||


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


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


oraz
oraz


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


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


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


Niech <math> \displaystyle  \displaystyle\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math> \displaystyle  C^1</math>
Niech <math>\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math>C^1</math>
oraz niech <math> \displaystyle  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><math> \displaystyle  s'(t)
<center>
\ =\
<math>s'(t)
=
\sqrt{\varphi'(t)^2+\psi'(t)^2}
\sqrt{\varphi'(t)^2+\psi'(t)^2}
\qquad\forall\  t\in[a,b].
\qquad\forall\  t\in[a,b]
</math></center>
</math>
</center>


}}
}}


{{dowod|twierdzenia 15.10.||
{{dowod|15.10. [nadobowiązkowy]||
(Dowód nadobowiązkowy.)<br>
Niech <math>t_0,t_0+h\in[a,b]</math>.
Niech <math> \displaystyle  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><math> \displaystyle  \sqrt{m_h^2+{m_h^*}^2}\cdot h
<center>
\ \le\
<math>\sqrt{m_h^2+{m_h^*}^2}\cdot h
\le
s(t_0+h)-s(t_0)
s(t_0+h)-s(t_0)
\ \le\
\le
\sqrt{M_h^2+{M_h^*}^2}\cdot h,
\sqrt{M_h^2+{M_h^*}^2}\cdot h
</math></center>
</math>
</center>


gdzie
gdzie


<center><math> \displaystyle  M_h
<center>
\ =\
<math>M_h
=
\sup_{t\in[t_0,t_0+h]}\varphi'(t),
\sup_{t\in[t_0,t_0+h]}\varphi'(t),
\qquad
\qquad
M_h^*
M_h^*
\ =\
=
\sup_{t\in[t_0,t_0+h]}\psi'(t),
\sup_{t\in[t_0,t_0+h]}\psi'(t)
</math></center>
</math>
</center>


<center><math> \displaystyle  m_h
<center>
\ =\
<math>m_h
=
\inf_{t\in[t_0,t_0+h]}\varphi'(t),
\inf_{t\in[t_0,t_0+h]}\varphi'(t),
\qquad
\qquad
m_h^*
m_h^*
\ =\
=
\displaystyle\inf_{t\in[t_0,t_0+h]}\psi'(t).
\inf_{t\in[t_0,t_0+h]}\psi'(t)
</math></center>
</math>
</center>


Dzielimy wszystkie strony
Dzielimy wszystkie strony
powyższego oszacowania
powyższego oszacowania
przez <math> \displaystyle  h,</math> dostając:
przez <math>h</math>, dostając:


<center><math> \displaystyle  \sqrt{m_h^2+{m_h^*}^2}
<center>
\ \le\
<math>\sqrt{m_h^2+{m_h^*}^2}
\le
\frac{s(t_0+h)-s(t_0)}{h}
\frac{s(t_0+h)-s(t_0)}{h}
\ \le\
\le
\sqrt{M_h^2+{M_h^*}^2}.
\sqrt{M_h^2+{M_h^*}^2}
</math></center>
</math>
</center>


Ponieważ funkcje <math> \displaystyle  \displaystyle\varphi'</math> i <math> \displaystyle  \displaystyle\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><math> \displaystyle  \aligned
<center>
<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),\\
M_h^* &  \xrightarrow[h\rightarrow 0]{} & \psi'(t_0),\\
M_h^* &  \xrightarrow[h\rightarrow 0]{} & \psi'(t_0),\\
m_h^* &  \xrightarrow[h\rightarrow 0]{} & \psi'(t_0).
m_h^* &  \xrightarrow[h\rightarrow 0]{} & \psi'(t_0).
\endaligned</math></center>
\end{align}</math>
</center>


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


<center><math> \displaystyle  s'(t_0)
<center>
\ =\
<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}
\ =\
=
\sqrt{\varphi'(t_0)^2+\psi'(t_0)^2}.
\sqrt{\varphi'(t_0)^2+\psi'(t_0)^2}
</math></center>
</math>
</center>


}}
}}


<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> \displaystyle  \displaystyle\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math> \displaystyle  C^1</math>
Niech <math>\varphi,\psi\colon [a,b]\longrightarrow\mathbb{R}</math> będą klasy <math>C^1</math>
oraz niech <math> \displaystyle  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><math> \displaystyle  l(K)
<center>
\ =\
<math>l(K)
\int\limits_a^b\sqrt{\varphi'(\tau)^2+\psi'(\tau)^2}\,d\tau.
=
</math></center>
\int\limits_a^b\sqrt{\varphi'(\tau)^2+\psi'(\tau)^2}\,d\tau</math>
</center>


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


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


}}</span>
}}</span>


{{dowod|twierdzenia 15.11.||
{{dowod|15.11.||


<center><math> \displaystyle  l(K)
<center>
\ =\
<math>l(K)
=
s(b)
s(b)
\ =\
=
s(b)-\underbrace{s(a)}\limits_{=0}
s(b)-\underbrace{s(a)}\limits_{=0}
\ =\
=
\int\limits_a^b s'(\tau)\,d\tau
\int\limits_a^b s'(\tau)\,d\tau
\ =\
=
\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</math>
</math></center>
</center>


W drugim przypadku krzywą zadaną przez funkcję <math> \displaystyle  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><math> \displaystyle 
<center>
<math>
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 461: Linia 499:
\qquad
\qquad
t\in[a,b]
t\in[a,b]
</math></center>
</math>
</center>


i od razu otrzymujemy drugi ze wzorów.
i od razu otrzymujemy drugi ze wzorów.
Linia 471: Linia 510:
krzywej zadanej w postaci biegunowej:
krzywej zadanej w postaci biegunowej:


<center><math> \displaystyle  r
<center>
\ =\
<math>r
=
g(\vartheta)
g(\vartheta)
\qquad
\qquad
\vartheta\in[\alpha,\beta].
\vartheta\in[\alpha,\beta]
</math></center>}}</span>
</math>
</center>}}</span>


Przedstawmy tę krzywą
Przedstawmy tę krzywą
w postaci parametrycznej:
w postaci parametrycznej:


<center><math> \displaystyle 
<center>
<math>
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
x\ = \ r\cos\vartheta\ =\ g(\vartheta)\cos\vartheta\\
x\ = \ r\cos\vartheta= g(\vartheta)\cos\vartheta\\
y\ = \ r\sin\vartheta\ =\ g(\vartheta)\sin\vartheta.
y\ = \ r\sin\vartheta= g(\vartheta)\sin\vartheta.
\end{array}  
\end{array}  
\right.
\right . </math>
</math></center>
</center>


Liczymy
Liczymy


<center><math>\begin{array}{lll}\displaystyle
<center>
<math>\begin{array}{lll}
\quad
\quad
x'(\vartheta)^2+y'(\vartheta)^2& = &
x'(\vartheta)^2+y'(\vartheta)^2& = &
Linia 508: Linia 551:
g'(\vartheta)^2
g'(\vartheta)^2
+g(\vartheta)^2,
+g(\vartheta)^2,
\end{array}</math></center>
\end{array}</math>
</center>


Zatem
Zatem


<center><math> \displaystyle  l(K)
<center>
\ =\
<math>l(K)
\displaystyle\int\limits_{\alpha}^{\beta}
=
\sqrt{g(\vartheta)^2+g'(\vartheta)^2}\,d\vartheta.
\int\limits_{\alpha}^{\beta}
</math></center>
\sqrt{g(\vartheta)^2+g'(\vartheta)^2}\,d\vartheta</math>
</center>


{| border="0" align="center" cellspacing="10"
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R09.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R09.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R09</div>
<div.thumbcaption>Krzywa we współrzędnych biegunowych</div>
</div></div>
</div></div>
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R10.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R10.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R10</div>
<div.thumbcaption>Cykloida</div>
</div></div>
</div></div>
|}
|}
Linia 533: Linia 578:


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


[[File:AM1.M15.W.R11.svg|250x250px|thumb|right|Cykloida]]
{{przyklad|15.14.||
{{przyklad|15.14.||


Wyprowadzić wzór parametryczny cykloidy.
Wyprowadzić wzór parametryczny cykloidy.
}}
<div class="mw-collapsible mw-made=collapsible mw-collapsed"><span class="mw-collapsible-toogle mw-collapsible-toogle-default style="font-variant:small-caps">Rozwiązanie </span><div class="mw-collapsible-content" style="display:none">
 
<div class="thumb tright"><div style="width:250px;">
<flash>file=AM1.M15.W.R11.swf|width=250|height=250</flash>
<div.thumbcaption>AM1.M15.W.R11</div></div>
</div>


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


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


<center><math> \displaystyle  x
<center>
\ = \
<math>x
\ =
OF
OF
\ =\
=
ON-FN
ON-FN
\ =\
=
\widehat{NM}-MG
\widehat{NM}-MG
\ =\
=
at-a\sin t
at-a\sin t</math>,
</math></center><br>
</center><br>


<center><math> \displaystyle  y
<center>
\ = \
<math>y
\ =  
FM
FM
\ =\
=
NG
NG
\ =\
=
ND-GD
ND-GD
\ =\
=
a-a\cos t.
a-a\cos t</math>
</math></center>
</center>


Zatem
Zatem


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


{{przyklad|15.15.||
{{przyklad|15.15.||
Linia 601: Linia 639:
Obliczyć długość łuku cykloidy:
Obliczyć długość łuku cykloidy:


<center><math> \displaystyle 
<center>
<math>
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 609: Linia 648:
\right.
\right.
\qquad
\qquad
t\in [0,2\pi].
t\in [0,2\pi]</math>
</math></center>
</center>  
 
<br>
}}
<center>
 
<math>\begin{array}{lll} \sqrt{x'(t)^2+y'(t)^2}
<div class="mw-collapsible mw-made=collapsible mw-collapsed"><span class="mw-collapsible-toogle mw-collapsible-toogle-default style="font-variant:small-caps">Rozwiązanie </span><div class="mw-collapsible-content" style="display:none"> 
 
<center><math>\begin{array}{lll} \displaystyle  \sqrt{x'(t)^2+y'(t)^2}
&=&
&=&
\sqrt{a^2(1-\cos t)^2+a^2\sin^2 t}
\sqrt{a^2(1-\cos t)^2+a^2\sin^2 t}
\ =\
=
\sqrt{a^2-2a^2\cos t+a^2\cos^2 t+a^2\sin^2 t}\\
\sqrt{a^2-2a^2\cos t+a^2\cos^2 t+a^2\sin^2 t}\\
& =&\displaystyle
& =&
\sqrt{2a^2(1-\cos t)}
\sqrt{2a^2(1-\cos t)}
\ =\
=
\sqrt{4a^2\sin^2\frac{t}{2}}
\sqrt{4a^2\sin^2\frac{t}{2}}
\ =\
=
2a\sin\frac{t}{2}.\end{array}
2a\bigg|\sin\frac{t}{2}\bigg|.\end{array}
</math></center>
</math>
</center>


Zatem
Zatem


<center><math> \displaystyle  l(K)
<center>
\ =\
<math>\begin{array}{lll}l(K)
\displaystyle\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
\ =\
=
2a\displaystyle\int\limits_0^{2\pi}\sin\frac{t}{2}\,dt
2a\int\limits_0^{2\pi}\bigg|\sin\frac{t}{2}\,\bigg|dt\\
\ =\
&=&
2a\int\limits_0^{2\pi}\sin\frac{t}{2}\,dt
=
-4a\cos\frac{t}{2}\bigg|_0^{2\pi}
-4a\cos\frac{t}{2}\bigg|_0^{2\pi}
\ =\
=8a.
8a.
\end{array}
</math></center>
</math>
 
</center>}}
</div></div>


{{przyklad|15.16.||
{{przyklad|15.16.||
Linia 649: Linia 687:
Obliczyć długość łuku asteroidy:
Obliczyć długość łuku asteroidy:


<center><math> \displaystyle  x^{\frac{2}{3}}+y^{\frac{2}{3}}
<center>
\ =\
<math>x^{\frac{2}{3}}+y^{\frac{2}{3}}
a^{\frac{2}{3}}
=
</math></center>
a^{\frac{2}{3}}</math>
</center>


}}
<div class="mw-collapsible mw-made=collapsible mw-collapsed"><span class="mw-collapsible-toogle mw-collapsible-toogle-default style="font-variant:small-caps">Rozwiązanie </span><div class="mw-collapsible-content" style="display:none"> 
Równanie parametryczne asteroidy, to:
Równanie parametryczne asteroidy, to:


{| border="0" align="center" cellspacing="10"
<center>
|<div class="thumb"><div style="width:260px;height:240;">
<math>
<flashwrap>file=AM1.M15.W.R12.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R12</div>
</div></div>
|<div class="thumb"><div style="width:260px;height:240;">
<flashwrap>file=AM1.M15.W.R13.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R13</div>
</div></div>
|<div class="thumb"><div style="width:260px;height:240;">
<flashwrap>file=AM1.M15.W.R14.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R14</div>
</div></div>
|}
 
 
 
<center><math> \displaystyle 
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 684: Linia 704:
\right.
\right.
\qquad
\qquad
t\in [0,2\pi].
t\in [0,2\pi]</math>
</math></center>
</center>


Liczymy
Liczymy


<center><math> \displaystyle  \sqrt{x'(t)^2+y'(t)^2}
<center>
\ =\
<math>\sqrt{x'(t)^2+y'(t)^2}
=
3a\sin t\cos t
3a\sin t\cos t
\qquad\forall\  t\in[0,2\pi].
\qquad\forall\  t\in[0,2\pi]</math>
</math></center>
</center>


Zatem
Zatem


<center><math> \displaystyle  l(K)
<center>
\ =\
<math>l(K)
4\displaystyle\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
6a.
=
</math></center>
6a</math>
 
</center>}}
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:253px;">
<flash>file=AM1.M15.W.R12.swf|width=272|height=272</flash>
<div.thumbcaption>Asteroida</div>
</div></div>
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R13.swf|size=small</flashwrap>
<div.thumbcaption>Asteroida</div>
</div></div>
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R14.swf|size=small</flashwrap>
<div.thumbcaption>Asteroida</div>
</div></div>
</div></div>
|}


==Całka krzywoliniowa==
==Całka krzywoliniowa==


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


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


Przypuśćmy, że określona jest funkcja ciągła
Przypuśćmy, że określona jest funkcja ciągła
<math> \displaystyle  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> \displaystyle  M</math>
to znaczy funkcja, która każdemu punktowi <math>M</math>
krzywej <math> \displaystyle  K</math> przyporządkowuje pewną wartość rzeczywistą
krzywej <math>K</math> przyporządkowuje pewną wartość rzeczywistą
<math> \displaystyle  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> \displaystyle  f</math> po krzywej <math> \displaystyle  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ę]]]]
Całkę tę wprowadza sie analogicznie jak całkę Riemanna na
Całkę tę wprowadza się analogicznie jak całkę Riemanna na
odcinku. Pominiemy to jednak w tym miejscu podając jedynie wzór
odcinku. Pominiemy to jednak w tym miejscu, podając jedynie wzór
końcowy na obliczanie takiej całki:
końcowy na obliczanie takiej całki:


<center><math> \displaystyle  \displaystyle\int\limits_K f(x,y)\,ds
<center>
\ =\
<math>\int\limits_K f(x,y)\,ds
\displaystyle\int\limits_a^b
=
f\big(\varphi(t),\psi(t)\big)\sqrt{\varphi'(t)^2+\psi'(t)^2}\,dt
\int\limits_a^b
</math></center>
f\big(\varphi(t),\psi(t)\big)\sqrt{\varphi'(t)^2+\psi'(t)^2}\,dt</math>
</center>


Tę całkę stosuje się w fizyce na przykład do obliczania masy i
Tę całkę stosuje się w fizyce na przykład do obliczania masy i
Linia 737: Linia 773:
długością są pomijalne).
długością są pomijalne).


Jeśli  mamy daną krzywą (pręt) <math> \displaystyle  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> \displaystyle  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> \displaystyle  \displaystyle\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><math> \displaystyle  m
<center>
\ =\
<math>m
\displaystyle\int\limits_K
=
\varrho(x,y)\,ds.
\int\limits_K
</math></center>
\varrho(x,y)\,ds</math>
</center>


Współrzędne środka ciężkości pręta
Współrzędne środka ciężkości pręta
<math> \displaystyle  \displaystyle (x_0,y_0)</math> możemy policzyć ze wzorów
<math>(x_0,y_0)</math> możemy policzyć ze wzorów


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


{{przyklad|15.17.||
{{przyklad|15.17.||


Obliczyć masę pręta półkolistego
Obliczyć masę pręta półkolistego
<math> \displaystyle  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> \displaystyle  \displaystyle\varrho(x,y)=y^2.</math>
o gęstości <math>\varrho(x,y)=y^2</math>.
}}


<div class="mw-collapsible mw-made=collapsible mw-collapsed"><span class="mw-collapsible-toogle mw-collapsible-toogle-default style="font-variant:small-caps">Rozwiązanie </span><div class="mw-collapsible-content" style="display:none"> 
Masa krzywej o gęstości <math>\varrho</math> dana jest wzorem
Masa krzywej o gęstości <math> \displaystyle  \displaystyle\varrho</math> dana jest wzorem


<center><math> \displaystyle  m
<center>
\ =\
<math>m
\displaystyle\int\limits_K \varrho(x,y)\,ds
=
\ =\
\int\limits_K \varrho(x,y)\,ds
\displaystyle\int\limits_K y^2\,ds.
=
</math></center>
\int\limits_K y^2\,ds</math>
</center>


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


<center><math> \displaystyle 
<center>
K=K(\varphi,\psi):\
<math>
K=K(\varphi,\psi):
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 783: Linia 822:
\end{array}  
\end{array}  
\right.
\right.
\qquad t\in[0,\pi],
\qquad t\in[0,\pi]</math>,
</math></center>
</center>


mamy
mamy


<center><math> \displaystyle  m
<center>
\ =\
<math>\begin{array}{lll} m
\displaystyle\int\limits_0^{\pi}R^2\sin^2 t\sqrt{(-R\sin t)^2+(R\cos t)^2}\,dt
&=&
\ =\
\int\limits_0^{\pi}R^2\sin^2 t\sqrt{(-R\sin t)^2+(R\cos t)^2}\,dt
R^3\displaystyle\int\limits_0^{\pi}\sin^2 t\,dt
=
\ =\
R^3\int\limits_0^{\pi}\sin^2 t\,dt\\
&=&
R^3\bigg[\frac{t}{2}-\frac{1}{4}\sin 2t\bigg]_0^{\pi}
R^3\bigg[\frac{t}{2}-\frac{1}{4}\sin 2t\bigg]_0^{\pi}
\ =\
=
\frac{R^3\pi}{2}.
\frac{R^3\pi}{2}.\end{array}
</math></center>
</math>
</center>


'''Odpowiedź:'''
'''Odpowiedź:'''
Masa pręta wynosi <math> \displaystyle  \displaystyle \frac{R^3\pi}{2}</math>
Masa pręta wynosi <math>\frac{R^3\pi}{2}</math>}}.
</div></div>


{{przyklad|15.18.||
{{przyklad|15.18.||


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


<div class="mw-collapsible mw-made=collapsible mw-collapsed"><span class="mw-collapsible-toogle mw-collapsible-toogle-default style="font-variant:small-caps">Rozwiązanie </span><div class="mw-collapsible-content" style="display:none"> 
<center>
Skora gęstość  <math> \displaystyle  \displaystyle\varrho</math> jest proporcjonalna
<math>\varrho(x,t)
do odległości punktu od środka układu i wynosi <math> \displaystyle  \displaystyle\sqrt{2}</math> w
=
punkcie <math> \displaystyle  \displaystyle (1,1),</math> to
 
<center><math> \displaystyle  \varrho(x,t)
\ =\
c\sqrt{x^2+y^2}
c\sqrt{x^2+y^2}
\quad </math> oraz <math> \displaystyle  \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> \displaystyle  c=1.</math>
stąd <math>c=1</math>.
Parametryzacją odcinka jest na przykład
Parametryzacją odcinka jest na przykład


<center><math> \displaystyle 
<center>
K=K(\varphi,\psi):\
<math>
K=K(\varphi,\psi):
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 835: Linia 875:
\end{array}  
\end{array}  
\right.
\right.
\qquad t\in[0,1],
\qquad t\in[0,1]</math>,
</math></center>
</center>


zatem masa wynosi
zatem masa wynosi


<center><math> \displaystyle  m
<center>
\ =\
<math>m
\displaystyle\int\limits_K\sqrt{x^2+y^2}\,ds
=
\ =\
\int\limits_K\sqrt{x^2+y^2}\,ds
\displaystyle\int\limits_0^1\sqrt{t^2+t^2}\sqrt{2}\,dt
=
\ =\
\int\limits_0^1\sqrt{t^2+t^2}\sqrt{2}\,dt
2\displaystyle\int\limits_0^1t\,dt
=
\ =\
2\int\limits_0^1t\,dt
=
t^2\bigg|_0^1
t^2\bigg|_0^1
\ =\
=
1.
1</math>
</math></center>
</center>


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


<center><math> \displaystyle  x_0
<center>
\ =\
<math>x_0
\frac{1}{m}\displaystyle\int\limits_K x\cdot \varrho(x,y)\,ds
=
\ =\
\frac{1}{m}\int\limits_K x\cdot \varrho(x,y)\,ds
\displaystyle\int\limits_0^1 t\cdot\sqrt{2t^2}\sqrt{2}\,dt
=
\ =\
\int\limits_0^1 t\cdot\sqrt{2t^2}\sqrt{2}\,dt
2\displaystyle\int\limits_0^1 t^2\,dt
=
\ =\
2\int\limits_0^1 t^2\,dt
=
\frac{2}{3}t^3\bigg|_0^1
\frac{2}{3}t^3\bigg|_0^1
\ =\
=
\frac{2}{3}.
\frac{2}{3}</math>
</math></center>
</center>


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


==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> \displaystyle  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).


<div class="thumb tright"><div style="width:250px;">
[[File:AM1.M15.W.R15.svg|250x250px|thumb|right|Pole między wykresami funkcji]]
<flash>file=AM1.M15.W.R15.swf|width=250|height=250</flash>
<div.thumbcaption>AM1.M15.W.R15</div>
</div></div>


Z poprzedniego wykładu znamy już związek całki Riemanna z
Z poprzedniego wykładu znamy już związek całki Riemanna z
Linia 893: Linia 931:
krzywymi:
krzywymi:


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


to pole tego trapezu wynosi:
to pole tego trapezu wynosi:


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


'''Uzasadnienie:'''
'''Uzasadnienie:'''
Linia 916: Linia 957:
parametrycznej
parametrycznej


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


wynosi
wynosi


<center><math> \displaystyle  |P|
<center>
\ =\
<math>|P|
\displaystyle\int\limits_{\alpha}^{\beta}\psi(t)\varphi'(t)\,dt.
=
</math></center>
\int\limits_{\alpha}^{\beta}\psi(t)\varphi'(t)\,dt</math>
</center>


'''Uzasadnienie:'''
'''Uzasadnienie:'''
Linia 939: Linia 982:
}}</span>
}}</span>


<div class="thumb tleft"><div style="width:260px;height:240;">
[[File:AM1.M15.W.R16.mp4|253x253px|thumb|left|Pole obszaru ograniczonego półprostymi i krzywą zadaną we współrzędnych biegunowych]]
<flashwrap>file=AM1.M15.W.R16.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W. R16</div>
</div></div>


<div class="thumb tright"><div style="width:250px;">
[[File:AM1.M15.W.R17.svg|250x250px|thumb|right|Trójkąt krzywoliniowy]]
<flash>file=AM1.M15.W.R17.swf|width=250|height=250</flash>
<div.thumbcaption>AM1.M15.W.R17</div>
</div></div>


<span id="twierdzenie_15_21">{{twierdzenie|15.21.||
<span id="twierdzenie_15_21">{{twierdzenie|15.21.||


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


<center><math> \displaystyle  r
<center>
\ =\
<math>r
=
g(\vartheta),
g(\vartheta),
\quad
\quad
\vartheta\in[\vartheta_1,\vartheta_2],
\vartheta\in[\vartheta_1,\vartheta_2]</math>,
</math></center>
</center>


to pole tego obszaru wynosi:
to pole tego obszaru wynosi:


<center><math> \displaystyle  |P|
<center>
\ =\
<math>|P|
\frac{1}{2}\displaystyle\int\limits_{\vartheta_1}^{\vartheta_2}\big[g(\vartheta)\big]^2\,d\vartheta.
=
</math></center>
\frac{1}{2}\int\limits_{\vartheta_1}^{\vartheta_2}\big[g(\vartheta)\big]^2\,d\vartheta</math>
</center>


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


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


<center><math> \displaystyle  P_{ABC} \approx
<center>
<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
\frac{1}{2} g(\vartheta)^2\Delta\vartheta
\frac{1}{2} g(\vartheta)^2\Delta\vartheta
</math></center>
</math>
</center>


(dla małych kątów <math> \displaystyle  \displaystyle\Delta\vartheta</math> zachodzi
(dla małych kątów <math>\Delta\vartheta</math> zachodzi
<math> \displaystyle  \displaystyle\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;
patrz [[Analiza matematyczna 1/Wykład 14: Całka Riemanna funkcji jednej zmiennej#definicja_14_4|definicja 14.4.]])
patrz [[Analiza matematyczna 1/Wykład 14: Całka Riemanna funkcji jednej zmiennej#definicja_14_4|definicja 14.4.]])
i przechodząc do granicy
i przechodząc do granicy,
dostajemy powyższy wzór.
dostajemy powyższy wzór.
}}</span>
}}</span>
<div class="thumb tright"><div style="width:260px;height:240;">
<flashwrap>file=AM1.M15.W.R18.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R18</div>
</div></div>


<span id="twierdzenie_15_22">{{twierdzenie|15.22.||
<span id="twierdzenie_15_22">{{twierdzenie|15.22.||
Linia 1000: Linia 1036:
Pole powierzchni powstałej z obrotu krzywej
Pole powierzchni powstałej z obrotu krzywej


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


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


<center><math> \displaystyle  |P|
<center>
\ =\
<math>|P|
=
2\pi
2\pi
\displaystyle\int\limits_a^b
\int\limits_a^b
\big[f(x)\big]
\big[f(x)\big]
\sqrt{1+f'(x)^2}\,dx.
\sqrt{1+f'(x)^2}\,dx</math>
</math></center>
</center>


Wzór ten pozostawiamy bez uzasadnienia.<br>
Wzór ten pozostawiamy bez uzasadnienia.<br>
Linia 1019: Linia 1058:
Pole powierzchni powstałej z obrotu krzywej
Pole powierzchni powstałej z obrotu krzywej


<center><math> \displaystyle  K:\
<center>
<math>K:
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 1025: Linia 1065:
y=\psi(t)
y=\psi(t)
\end{array}  
\end{array}  
\right.,
\right.
\quad </math> dla <math> \displaystyle  \ t\in[\alpha,\beta]
\quad</math> dla <math>\ t\in[\alpha,\beta]
</math></center>
</math>
</center>


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


<center><math> \displaystyle  |P|
<center>
\ =\
<math>|P|
=
2\pi
2\pi
\displaystyle\int\limits_{\alpha}^{\beta}
\int\limits_{\alpha}^{\beta}
\big[\psi(t)\big]
\big[\psi(t)\big]
\sqrt{\varphi'(t)^2+\psi'(t)^2}\,dt.
\sqrt{\varphi'(t)^2+\psi'(t)^2}\,dt</math>
</math></center>
</center>


}}</span>
}}</span>
{| border="0" align="center" cellspacing="10"
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R18.swf|size=small</flashwrap>
<div.thumbcaption>Powierzchnia powstała przez obrót krzywej dookoła osi <math>Ox</math></div>
</div></div>
|<div class="thumb"><div style="width:253px;">
<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 <math>Ox</math></div>
</div></div>
|}
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R19.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R19.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R19</div>
<div.thumbcaption>Powierzchnia powstała przez obrót krzywej dookoła osi <math>Ox</math></div>
</div></div>
</div></div>
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R20.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R20.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R20</div>
<div.thumbcaption>Powierzchnia powstała przez obrót krzywej dookoła osi <math>Ox</math></div>
</div></div>
</div></div>
|}
|}
<div class="thumb tright"><div style="width:260px;height:240;">
<flashwrap>file=AM1.M15.W.R21.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R21</div>
</div></div>


<span id="twierdzenie_15_23">{{twierdzenie|15.23.||
<span id="twierdzenie_15_23">{{twierdzenie|15.23.||
Linia 1063: Linia 1109:
obszaru "pod krzywą"
obszaru "pod krzywą"


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


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


<center><math> \displaystyle  |V_x|
<center>
\ =\
<math>|V_x|
=
\pi
\pi
\displaystyle\int\limits_a^b
\int\limits_a^b
f(x)^2\,dx.
f(x)^2\,dx</math>
</math></center>
</center>


'''Uzasadnienie:'''
'''Uzasadnienie:'''
Weźmy podział odcinka <math> \displaystyle  \displaystyle [a,b]</math>:
Weźmy podział odcinka <math>[a,b]</math>:


<center><math> \displaystyle  P:\
<center>
<math>P:
a
a
\ =\
=
x_0
x_0
\ <\
<
x_1
x_1
\ <\
<
\ldots
\ldots
\ <\
<
x_n
x_n
\ =\
=
b
b
</math></center>
</math>
</center>


oraz podzielmy bryłę na "plasterki",
oraz podzielmy bryłę na "plasterki",
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> \displaystyle  y=f(x)</math> dla <math> \displaystyle  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> \displaystyle  f(x_i)</math> i wysokości
objętości walca o promieniu podstawy <math>f(x_i)</math> i wysokości
<math> \displaystyle  \displaystyle\Delta x_i=x_i-x_{i-1},</math> czyli
<math>\Delta x_i=x_i-x_{i-1}</math>, czyli
<math> \displaystyle  \displaystyle\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
i przechodząc do granicy dostajemy powyższy wzór.<br>
i przechodząc do granicy, dostajemy powyższy wzór.<br>
<br>
<br>
'''(2)'''
'''(2)'''
Linia 1109: Linia 1160:
obszaru "pod krzywą"
obszaru "pod krzywą"


<center><math> \displaystyle  K:\
<center>
<math>K:
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 1115: Linia 1167:
y=\psi(t)
y=\psi(t)
\end{array}  
\end{array}  
\right.,
\right.
\quad </math> dla <math> \displaystyle  \ t\in[\alpha,\beta]
\quad</math> dla <math>\ t\in[\alpha,\beta]
</math></center>
</math>
</center>


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


<center><math> \displaystyle  |V_x|
<center>
\ =\
<math>|V_x|
=
\pi
\pi
\displaystyle\int\limits_{\alpha}^{\beta}
\int\limits_{\alpha}^{\beta}
\psi(t)^2\varphi'(t)\,dt.
\psi(t)^2\varphi'(t)\,dt</math>
</math></center>
</center>


'''Uzasadnienie:'''
'''Uzasadnienie:'''
Linia 1134: Linia 1188:


{| border="0" align="center" cellspacing="10"
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R22.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R22.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R22</div>
<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi <math>Ox</math></div>
</div></div>
</div></div>
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R23.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R23.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R23</div>
<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi <math>Ox</math></div>
</div></div>
</div></div>
|}
|}


<div class="thumb tright"><div style="width:260px;height:240;">
<div class="thumb tright"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R24.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R24.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R24</div>
<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi <math>Oy</math></div>
 
</div></div>
</div></div>


Linia 1155: Linia 1210:
obszaru "pod krzywą"
obszaru "pod krzywą"


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


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


<center><math> \displaystyle  |V_y|
<center>
\ =\
<math>|V_y|
=
2\pi
2\pi
\displaystyle\int\limits_a^b
\int\limits_a^b
x\,f(x)\,dx.
x\,f(x)\,dx</math>
</math></center>
</center>


'''Uzasadnienie:'''
'''Uzasadnienie:'''
Weźmy podział odcinka <math> \displaystyle  \displaystyle [a,b]</math>:
Weźmy podział odcinka <math>[a,b]</math>:


<center><math> \displaystyle  P:\
<center>
<math>P:
a
a
\ =\
=
x_0
x_0
\ <\
<
x_1
x_1
\ <\
<
\ldots
\ldots
\ <\
<
x_n
x_n
\ =\
=
b
b
</math></center>
</math>
</center>


oraz podzielmy bryłę na "cylindry"
oraz podzielmy bryłę na "cylindry"
powstałe przez
powstałe przez
obrót obszaru pod wykresem funkcji
obrót obszaru pod wykresem funkcji
<math> \displaystyle  y=f(x)</math> dla <math> \displaystyle  x\in[x_{i-1},x_i]</math> wokół osi <math> \displaystyle  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> \displaystyle  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> \displaystyle  |V_y|.</math><br>
na <math>|V_y|</math>.<br>
<br>
<br>
'''(2)'''
'''(2)'''
Linia 1200: Linia 1260:
obszaru "pod krzywą"
obszaru "pod krzywą"


<center><math> \displaystyle  K:\
<center>
<math>K:
\left\{
\left\{
\begin{array} {l}
\begin{array} {l}
Linia 1206: Linia 1267:
y=\psi(t)
y=\psi(t)
\end{array}  
\end{array}  
\right.,
\right.
\quad </math> dla <math> \displaystyle  \ t\in[\alpha,\beta]
\quad</math> dla <math>\ t\in[\alpha,\beta]
</math></center>
</math>
</center>


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


<center><math> \displaystyle  |V_y|
<center>
\ =\
<math>|V_y|
=
2\pi
2\pi
\displaystyle\int\limits_{\alpha}^{\beta}
\int\limits_{\alpha}^{\beta}
\varphi(t)\psi(t)\varphi'(t)\,dt.
\varphi(t)\psi(t)\varphi'(t)\,dt</math>
</math></center>
</center>


}}</span>
}}</span>


{| border="0" align="center" cellspacing="10"
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R25.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R25.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R25</div>
<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi <math>Ox</math></div>
</div></div>
</div></div>
|<div class="thumb"><div style="width:260px;height:240;">
|<div class="thumb"><div style="width:253px;">
<flashwrap>file=AM1.M15.W.R26.swf|size=small</flashwrap>
<flashwrap>file=AM1.M15.W.R26.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R26</div>
<div.thumbcaption>Bryła powstała przez obrót obszaru "pod" wykresem krzywej wokół osi <math>Ox</math></div>
</div></div>
</div></div>
|}
|}


<div class="thumb tright"><div style="width:260px;height:240;">
[[File:AM1.M15.W.R27.mp4|253x253px|thumb|right|Torus]]
<flashwrap>file=AM1.M15.W.R27.swf|size=small</flashwrap>
<div.thumbcaption>AM1.M15.W.R27</div>
</div></div>


{{przyklad|15.25.||
{{przyklad|15.25.||
Linia 1242: Linia 1302:
koła
koła


<center><math> \displaystyle  x^2+(y-a)^2
<center>
\ \le\
<math>x^2+(y-a)^2
\le
r^2
r^2
\qquad
\qquad
(0<r<a)
(0<r<a)
</math></center>
</math>
</center>


wokół osi <math> \displaystyle  Ox.</math>}}
wokół osi <math>Ox</math>.


<div class="mw-collapsible mw-made=collapsible mw-collapsed"><span class="mw-collapsible-toogle mw-collapsible-toogle-default style="font-variant:small-caps">Rozwiązanie </span><div class="mw-collapsible-content" style="display:none">  
<center>
 
<math>\begin{array}{lll}
<center><math>\begin{array}{lll} \displaystyle
|V_x|
|V_x|
& = &
& = &
\pi\displaystyle\int\limits_{-r}^r
\pi\int\limits_{-r}^r
\bigg[\big(a+\sqrt{r^2-x^2}\big)^2
\bigg[\big(a+\sqrt{r^2-x^2}\big)^2
-\big(a-\sqrt{r^2-x^2}\big)^2\bigg]\,dx
-\big(a-\sqrt{r^2-x^2}\big)^2\bigg]\,dx
\ =\
=
4\pi a\displaystyle\int\limits_{-r}^r
4\pi a\int\limits_{-r}^r
\sqrt{r^2-x^2}\,dx\\
\sqrt{r^2-x^2}\,dx\\
& \stackrel{(\star)}{=} &
& \stackrel{(\bigstar)}{=} &
4\pi a
4\pi a
\bigg[\frac{r^2}{2}\arcsin\frac{x}{r}+\frac{x}{2}\sqrt{r^2-x^2}\bigg]_{-r}^r
\bigg[\frac{r^2}{2}\arcsin\frac{x}{r}+\frac{x}{2}\sqrt{r^2-x^2}\bigg]_{-r}^r
=\
=
4\pi a
4\pi a
\bigg[\frac{r^2}{2}\cdot\frac{\pi}{2}+\frac{r^2}{2}\cdot\frac{\pi}{2}\bigg]
\bigg[\frac{r^2}{2}\cdot\frac{\pi}{2}+\frac{r^2}{2}\cdot\frac{\pi}{2}\bigg]
\ =\
=
4\pi a\frac{r^2\pi}{2}
4\pi a\frac{r^2\pi}{2}
\ =\
=
2\pi^2 ar^2,
2\pi^2 ar^2,
\end{array}</math></center>
\end{array}</math>
</center>


gdzie wykorzystano następującą całkę:<br>
gdzie wykorzystano następującą całkę:<br>


<center><math> \displaystyle  \aligned
<center>
(\star)\quad
<math>\begin{align}
(\bigstar)\quad
I
I
& = &
& = &
\int\sqrt{r^2-x^2}\,dx
\int\sqrt{r^2-x^2}\,dx
\ =\
=
\int\frac{r^2-x^2}{\sqrt{r^2-x^2}}\,dx
\int\frac{r^2-x^2}{\sqrt{r^2-x^2}}\,dx
\ =\
=
r^2\underbrace{\int\frac{dx}{\sqrt{r^2-x^2}}}_{I_1}
r^2\underbrace{\int\frac{dx}{\sqrt{r^2-x^2}}}_{I_1}
-\underbrace{\int\frac{x^2 dx}{\sqrt{r^2-x^2}}}_{I_2}.
-\underbrace{\int\frac{x^2 dx}{\sqrt{r^2-x^2}}}_{I_2}.
\endaligned</math></center>
\end{align}</math>
</center>


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


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


<center><math>\begin{array}{lll} \displaystyle
<center>
<math>\begin{array}{lll}  
I& = &\int\sqrt{r^2-x^2}\,dx
I& = &\int\sqrt{r^2-x^2}\,dx
\ \begin{array}{c}\textrm{części}\\=\end{array} x\sqrt{r^2-x^2}
\ \begin{array}{c}\text{części}\\=\end{array} x\sqrt{r^2-x^2}
-\int x\frac{-2x}{2\sqrt{r^2-x^2}}dx\\
-\int x\frac{-2x}{2\sqrt{r^2-x^2}}dx\\
& =&
& =&
Linia 1303: Linia 1369:
+\underbrace{\int\frac{x^2}{\sqrt{r^2-x^2}}\,dx}_{=I_2}=
+\underbrace{\int\frac{x^2}{\sqrt{r^2-x^2}}\,dx}_{=I_2}=
x\sqrt{r^2-x^2}+I_2.
x\sqrt{r^2-x^2}+I_2.
\end{array}</math></center>
\end{array}</math>
</center>


Porównując to z <math> \displaystyle  \displaystyle (\star),</math>
Porównując to z <math>(\bigstar)</math>,
otrzymujemy:
otrzymujemy:


<center><math> \displaystyle  r^2I_1-I_2
<center>
\ =\
<math>r^2I_1-I_2
x\sqrt{r^2-x^2}+I_2,
=
</math></center>
x\sqrt{r^2-x^2}+I_2</math>,
</center>


stąd
stąd


<center><math> \displaystyle  2I_2
<center>
\ =\
<math>2I_2
=
r^2I_1-x\sqrt{r^2-x^2}
r^2I_1-x\sqrt{r^2-x^2}
\ =\
=
r^2\arcsin\frac{x}{r}
r^2\arcsin\frac{x}{r}
-x\sqrt{r^2-x^2},
-x\sqrt{r^2-x^2}</math>,
</math></center>
</center>


zatem
zatem


<center><math> \displaystyle  I_2
<center>
\ =\
<math>I_2
=
\frac{r^2}{2}\arcsin\frac{x}{r}
\frac{r^2}{2}\arcsin\frac{x}{r}
-\frac{x}{2}\sqrt{r^2-x^2}.
-\frac{x}{2}\sqrt{r^2-x^2}</math>
</math></center>
</center>


Wstawiając do <math> \displaystyle  \displaystyle (\star),</math> otrzymujemy:
Wstawiając do <math>(\bigstar)</math>, otrzymujemy:


<center><math> \displaystyle  \aligned
<center>
<math>\begin{align}
I
I
& = &
& = &
Linia 1339: Linia 1410:
-\frac{1}{2}r^2\arcsin\frac{x}{r}
-\frac{1}{2}r^2\arcsin\frac{x}{r}
+\frac{1}{2}x\sqrt{r^2-x^2}+c
+\frac{1}{2}x\sqrt{r^2-x^2}+c
\ =\
=
\frac{r^2}{2}\arcsin\frac{x}{r}
\frac{r^2}{2}\arcsin\frac{x}{r}
+\frac{x}{2}\sqrt{r^2-x^2}+c.
+\frac{x}{2}\sqrt{r^2-x^2}+c.
\endaligned</math></center>
\end{align}</math>
 
</center>}}
</div></div>

Aktualna wersja na dzień 08:25, 24 lip 2024

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.