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 – „\displaystyle ” na „”
m Zastępowanie tekstu – „\displaystyle” na „”
Linia 23: Linia 23:
</center>
</center>


gdzie <math> \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:


Linia 177: Linia 177:
<span id="twierdzenie_15_7">{{twierdzenie|15.7.||
<span id="twierdzenie_15_7">{{twierdzenie|15.7.||


Niech <math> \displaystyle\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>
Linia 229: Linia 229:
</center>
</center>


Ponieważ <math> \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
Linia 262: Linia 262:
</center>
</center>


Ponieważ <math> \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>[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 279: Linia 279:
<center>
<center>
<math> m =
<math> m =
\displaystyle\inf_{t\in[a,b]}\varphi'(t),
\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>
Linia 358: Linia 358:
{{twierdzenie|15.10.||
{{twierdzenie|15.10.||


Niech <math> \displaystyle\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>
Linia 407: Linia 407:
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>
</math>
</center>
</center>
Linia 424: Linia 424:
</center>
</center>


Ponieważ funkcje <math> \displaystyle\varphi'</math> i <math> \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


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> \displaystyle\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
Linia 539: Linia 539:


<center>
<center>
<math>\begin{array}{lll}\displaystyle
<math>\begin{array}{lll}
\quad
\quad
x'(\vartheta)^2+y'(\vartheta)^2& = &
x'(\vartheta)^2+y'(\vartheta)^2& = &
Linia 563: Linia 563:
<math> l(K)
<math> l(K)
=
=
\displaystyle\int\limits_{\alpha}^{\beta}
\int\limits_{\alpha}^{\beta}
\sqrt{g(\vartheta)^2+g'(\vartheta)^2}\,d\vartheta.
\sqrt{g(\vartheta)^2+g'(\vartheta)^2}\,d\vartheta.
</math>
</math>
Linia 679: Linia 679:
<center>
<center>
<math> \begin{array}{lll}l(K)
<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}\bigg|\sin\frac{t}{2}\,\bigg|dt\\
2a\int\limits_0^{2\pi}\bigg|\sin\frac{t}{2}\,\bigg|dt\\
&=&
&=&
2a\displaystyle\int\limits_0^{2\pi}\sin\frac{t}{2}\,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}
Linia 733: Linia 733:
<math> l(K)
<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.
6a.
Linia 777: Linia 777:


<center>
<center>
<math> \displaystyle\int\limits_K f(x,y)\,ds
<math> \int\limits_K f(x,y)\,ds
=
=
\displaystyle\int\limits_a^b
\int\limits_a^b
f\big(\varphi(t),\psi(t)\big)\sqrt{\varphi'(t)^2+\psi'(t)^2}\,dt.
f\big(\varphi(t),\psi(t)\big)\sqrt{\varphi'(t)^2+\psi'(t)^2}\,dt.
</math>
</math>
Linia 790: Linia 790:
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> \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


Linia 796: Linia 796:
<math> m
<math> m
=
=
\displaystyle\int\limits_K
\int\limits_K
\varrho(x,y)\,ds.
\varrho(x,y)\,ds.
</math>
</math>
Linia 805: Linia 805:


<center>
<center>
<math> \begin{align} x_0 & = & \frac{1}{m}\displaystyle\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}\displaystyle\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>
</center>
</center>
Linia 814: Linia 814:
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> \displaystyle\varrho(x,y)=y^2.</math>
o gęstości <math> \varrho(x,y)=y^2.</math>


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


<center>
<center>
<math> m
<math> m
=
=
\displaystyle\int\limits_K \varrho(x,y)\,ds
\int\limits_K \varrho(x,y)\,ds
=
=
\displaystyle\int\limits_K y^2\,ds.
\int\limits_K y^2\,ds.
</math>
</math>
</center>
</center>
Linia 848: Linia 848:
<math>\begin{array}{lll} m
<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}
Linia 866: Linia 866:
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> \displaystyle\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> \displaystyle\varrho</math> jest proporcjonalna
Skoro gęstość  <math> \varrho</math> jest proporcjonalna
do odległości punktu od środka układu i wynosi <math> \displaystyle\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


Linia 903: Linia 903:
<math> m
<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
Linia 921: Linia 921:
<math> x_0
<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
Linia 966: Linia 966:
<math> |P|
<math> |P|
=
=
\displaystyle\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
</math>
</math>
</center>
</center>
Linia 998: Linia 998:
<math> |P|
<math> |P|
=
=
\displaystyle\int\limits_{\alpha}^{\beta}\psi(t)\varphi'(t)\,dt.
\int\limits_{\alpha}^{\beta}\psi(t)\varphi'(t)\,dt.
</math>
</math>
</center>
</center>
Linia 1032: Linia 1032:
<math> |P|
<math> |P|
=
=
\frac{1}{2}\displaystyle\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.
</math>
</math>
</center>
</center>
Linia 1049: Linia 1049:
</center>
</center>


(dla małych kątów <math> \displaystyle\Delta\vartheta</math> zachodzi
(dla małych kątów <math> \Delta\vartheta</math> zachodzi
<math> \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;
Linia 1075: Linia 1075:
=
=
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.
Linia 1104: Linia 1104:
=
=
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.
Linia 1150: Linia 1150:
=
=
\pi
\pi
\displaystyle\int\limits_a^b
\int\limits_a^b
f(x)^2\,dx.
f(x)^2\,dx.
</math>
</math>
Linia 1180: Linia 1180:
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> \displaystyle\Delta x_i=x_i-x_{i-1},</math> czyli
<math> \Delta x_i=x_i-x_{i-1},</math> czyli
<math> \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
Linia 1208: Linia 1208:
=
=
\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>
Linia 1253: Linia 1253:
=
=
2\pi
2\pi
\displaystyle\int\limits_a^b
\int\limits_a^b
x\,f(x)\,dx.
x\,f(x)\,dx.
</math>
</math>
Linia 1310: Linia 1310:
=
=
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>
Linia 1350: Linia 1350:
|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{(\bigstar)}{=} &
& \stackrel{(\bigstar)}{=} &
Linia 1395: Linia 1395:


<center>
<center>
<math>\begin{array}{lll} \displaystyle
<math>\begin{array}{lll}  
I& = &\int\sqrt{r^2-x^2}\,dx
I& = &\int\sqrt{r^2-x^2}\,dx
\ \begin{array}{c}\text{części}\\=\end{array} x\sqrt{r^2-x^2}
\ \begin{array}{c}\text{części}\\=\end{array} x\sqrt{r^2-x^2}

Wersja z 10:20, 28 sie 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.