Analiza matematyczna 1/Wykład 15: Krzywe i bryły obrotowe

Z Studia Informatyczne
Przejdź do nawigacjiPrzejdź do wyszukiwania

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.