Rozdział 2. Różniczkowanie odwzorowań
2.19. Średnia uogólniona
Wynika stąd, że
feλ00(a)(X) = −c n
n−2
kXk2< 0, X ∈ (TaM )∗, a więc ef ma w punkcie a maksimum warunkowe, czyli
x1· · · xn 6c n
n
=x1+ · · · + xn n
n
, x1, . . . , xn> 0.
2.19. Średnia uogólniona
Ustalmy n ∈ N2 oraz liczby a1, . . . , an, t1, . . . , tn > 0 takie, że a1, . . . , an są parami różne oraz t1+ · · · + tn = 1. Niech a := (a1, . . . , an), t := (t1, . . . , tn). Rozważmy funkcję
R∗3 x7−→Sa,t
t1ax1+ · · · + tnaxn1/x
∈ R>0.
Uwaga: Formalnie można zrezygnować z założenia, że a1, . . . , an są parami różne, ponieważ (poprzez zmianę n i t1, . . . , tn) zawsze można doprowadzić do sytuacji, w której bądź n = 1, bądź a1, . . . , an są parami różne.
Lemat 2.19.1. (a) lim
x→+∞Sa,t(x) = max{a1, . . . , an} =: Sa,t(+∞).
(b) lim
x→0Sa,t(x) = at11· · · atnn=: Sa,t(0).
(c) lim
x→−∞Sa,t(x) = min{a1, . . . , an} =: Sa,t(−∞).
Dowód . (a) i (c) są elementarne (Ćwiczenie).
(b) Zauważmy, że Sa,t = exp(ϕ), gdzie ϕ(x) := 1xln(t1ax1 + · · · + tnaxn). Granica limx→0ϕ(x) jest symbolem nieoznaczonym typu 00. Stosujemy regułę d’Hôpitala 74
x→0limϕ(x)H= lim
x→0
t1ax1ln a1+ · · · + tnaxnln an
t1ax1+ · · · + tnaxn = t1ln a1+ · · · + tnln an. Definicja 2.19.2. Funkcję Sa,t : R −→ R>0 nazywamy średnią uogólnioną liczb a1, . . . , an z wagami t1, . . . , tn.
Obserwacja 2.19.3. W przypadku, gdy t1 = · · · = tn := 1/n, niech Sa := Sa,t. Wtedy dostajemy klasyczne średnie:
• Sa(2) =
ra21+ · · · + a2n
n = średnia kwadratowa;
• Sa(1) = a1+ · · · + an
n = średnia arytmetyczna;
• Sa(0) = √n
a1· · · an= średnia geometryczna;
• Sa(−1) = n
1
a1 + · · · +a1
n
= średnia harmoniczna.
Uwaga: Definiując powyższe średnie rezygnujemy z założenia, że aj 6= ak dla j 6= k.
Lemat 2.19.4. Funkcja Sa,t: R −→ R jest analityczna.
Dowód . Przypomnijmy najpierw, że dla c > 0 funcja R 3 x 7−→ cx = exp(x ln c) jest analityczna (zob. Wniosek 2.4.5). Wynika stąd łatwo, że funkcja R 3 x7−→ ln(tψ 1ax1 + · · · + tnaxn) jest analityczna.
Wiemy, że ψ(0) = 0. Wynika stąd natychmiast, że funkcja ϕ jest analityczna (ϕ(0) := ψ0(0)). Istotnie, jeżeli ψ(x) = P∞
k=1ckxk, |x| < r, to ϕ(x) = P∞
k=1ckxk−1, |x| < r. Teraz, korzystając ponownie z
Wniosku 2.4.5, dostajemy analityczność funkcji Sa,t.
Propozycja 2.19.5. Funkcja Sa,t jest ściśle rosnąca.
74 Guillaume de L’Hôpital (1661–1704).
2. Różniczkowanie odwzorowań
Dowód . Wystarczy udowodnić ścisłą monotoniczność funkcji ϕ. Mamy ϕ0(x) = −1
x2ln(t1ax1+ · · · + tnaxn) + 1 x
t1ax1ln a1+ · · · + tnaxnln an t1ax1+ · · · + tnaxn . Ustalmy x ∈ R i połóżmy uj= uj(x) := axj, j = 1, . . . , n. Widzimy, że wystarczy pokazać, że
f (u1, . . . , un) := − ln(t1u1+ · · · + tnun) +t1u1ln u1+ · · · + tnunln un
t1u1+ · · · + tnun > 0, u1, . . . , un > 0, przy czym równość zachodzi wyłącznie dla u1= · · · = un.
Jest oczywiste, że f (u, . . . , u) = 0. Ustalmy liczbę p > 0 i rozważmy zbiór M := {(u1, . . . , un) ∈ Rn>0: t1u1+ · · · + tnun= p}.
Zauważmy, że dla (u1, . . . , un) ∈ M mamy
f (u1, . . . , un) = h(u1, . . . , un) := − ln p + 1
p(t1u1ln u1+ · · · + tnunln un).
Wobec zależności limδ→0+δ ln δ = 0, funkcja h przedłuża się do funkcji ciągłej eh na zbiorze zwartym M = {(u1, . . . , un) ∈ Rn+ : t1u1+ · · · + tnun = p}. Niech u0 = (u01, . . . , u0n) ∈ M będzie punktem, w którym eh przyjmuje globalne minimum na M . Wiemy, że eh(u0) 6 0. Rozważmy dwa przypadki:
(a) u0 ∈ M : Stosujemy metodę mnożników Lagrange’a (do funkcji h z warunkiem g(u1, . . . , un) :=
t1u1+ · · · + tnun− p = 0). Warunek Lagrange’a
∂h
∂uj(u0) − λ∂g
∂uj(u0) = 0, j = 1, . . . , n, daje
1
p(ln u0j+ 1) − λ = 0, j = 1, . . . , n, a stąd u01= · · · = u0n, co oznacza, że h(u0) = 0.
(b) u0 ∈ M \ M : Możemy założyć, że dla pewnego 1 6 s 6 n − 1 mamy u01 > 0, . . . , u0s > 0, u0s+1 = · · · = u0n = 0. Zauważmy, że eh(u1, . . . , us, 0, . . . , 0) = − ln p + 1p(t1u1ln u1+ · · · + tsusln us), (u1, . . . , us) ∈ M0 := {(u1, . . . , us) ∈ Rs>0: t1u1+ · · · + tsus= p}. Stosując rozumowanie z (a) względem zmiennych u1, . . . , us, wnioskujemy, że musi być u01 = · · · = u0s = t p
1+···+ts. Wynika stąd, że eh(u0) =
− ln p +1p(t1+ · · · + ts)t p
1+···+tslnt p
1+···+ts = − ln(t1+ · · · + ts) > 0 — sprzeczność.
Oznacza to, że globalne minimum funkcji h na M jest równe 0. Wobec rozumowania w (a), widzimy również, że jeżeli u0∈ M jest taki, że h(u0) = 0, to u01= · · · = u0n. Propozycja 2.19.6. Istnieją liczby x± ∈ R takie, że funkcja Sa,t jest wypukła w przedziale (−∞, x−) i wklęsła w (x+, +∞).
Dowód . Wiadomo, że Sa,t00 = Sa,t((ϕ0)2+ ϕ00). Obliczmy ϕ00(x):
ϕ00(x) = 2
x3ln(t1ax1+ · · · + tnaxn) − 2 x2
t1ax1ln a1+ · · · + tnaxnln an t1ax1+ · · · + tnaxn +1
x
t1ax1ln2a1+ · · · + tnaxnln2an
t1ax1+ · · · + tnaxn −1 x
t1ax1ln a1+ · · · + tnaxnln an
t1ax1+ · · · + tnaxn
2
. Rozumowanie przy x −→ −∞: Wystarczy pokazać, że lim
x→−∞x3ϕ00(x) < 0. Możemy założyć, że a1< · · · < an. Podstawmy uj= uj(x) := axj i vj = vj(x) = uj/u1= (aj/a1)x. Zauważmy, że vj(x) −→ 0 gdy x −→ −∞ dla j = 2, . . . , n. Dla ξ = (ξ1, . . . , ξn) ∈ Rn>0 niech
M (ξ) : = t1ξ1+ · · · + tnξn,
P (ξ) : = t1ξ1ln ξ1+ · · · + tnξnln ξn, Q(ξ) : = t1ξ1ln2ξ1+ · · · + tnξnln2ξn.
2.19. Średnia uogólniona
Zauważmy, że M (v) −→ t1, P (v) −→ 0, Q(v) −→ 0 gdy x −→ −∞. Mamy x3ϕ00(x) = 2 ln M (u) − 2P (u)
M (u)+ Q(u)
M (u)−P (u) M (u)
2
= 2 ln M (v) − 2 ln u1− 2t1v1(ln v1+ ln u1) + · · · + tnvn(ln vn+ ln u1) M (v)
+t1v1(ln v1+ ln u1)2+ · · · + tnvn(ln vn+ ln u1)2
M (v) −t1v1(ln v1+ ln u1) + · · · + tnvn(ln vn+ ln u1) M (v)
2
= 2 ln M (v) + 2 ln u1− 2P (v)
M (v)− 2 ln u1+ Q(v)
M (v)+ 2P (v)
M (v)ln u1+ ln2u1−P (v)
M (v)+ ln u1
2
= 2 ln M (v) − 2P (v)
M (v)+ Q(v)
M (v)−P (v) M (v)
2
x→−∞−→ 2 ln t1< 0.
Rozumowanie przy x −→ +∞: Wystarczy pokazać, że lim
x→+∞x4((ϕ0(x))2+ ϕ00(x)) = −∞. Możemy założyć, że a1> · · · > an. Tak, jak poprzednio definiujemy uj, vj, M , P , Q. Zauważmy, że vj(x) −→ 0 gdy x −→ +∞ dla j = 2, . . . , n, oraz M (v) −→ t1, P (v) −→ 0, Q(v) −→ 0 gdy x −→ +∞. Mamy
x4((ϕ0(x))2+ ϕ00(x)) =
− ln M (v) + P (v) M (v)
2
+ x
2 ln M (v) − 2P (v)
M (v)+ Q(v)
M (v)−P (v) M (v)
2
x→+∞−→ ln2t1+ (+∞) · (2 ln t1) = −∞. Ćwiczenie* 2.19.7. Zbadać wypukłość funkcji Sa,t. Odnotujmy, że problem wypukłości funkcji Sa,tjest wysoce nietrywialny. Dla przykładu, korzystając z pomocy komputera można łatwo, sprawdzić, że dla n = 4, a1 := 0.1635, a2:= 4.7965, a3 := 9.3668, a4 := 1.7856, t1 := 0.0455, t2 := 0.1430, t3:= 0.0007, t4:= 0, 8108, mamy co najmniej 5 punktów przegięcia.
Całka Riemanna
3.1. Całka Riemanna na kostce Definicja 3.1.1. Kostką (domkniętą) nazywamy dowolny zbiór postaci
P := [a1, b1] × · · · × [an, bn] ⊂ Rn, gdzie aj< bj, j = 1, . . . , n. 1
Objętością kostki P nazywamy liczbę
|P | := (b1− a1) · · · (bn− an). 2
Podziałem kostki P nazywamy dowolną skończoną rodzinę kostek π = {P1, . . . , Pm}
taką, że P = P1∪ · · · ∪ Pm oraz int Pj∩ int Pk= ∅ dla j 6= k 3 . Podziałem prostym kostki P nazywamy podział postaci
{Pk1,...,kn: 1 6 kj 6 mj, 1 6 j 6 n}, gdzie
Pk1,...,kn= [x1,k1−1, x1,k1] × · · · × [xn,kn−1, xn,kn], zaś (xj,0, . . . , xj,mj) jest podziałem [aj, bj], j = 1, . . . , n 4 .
Niech π1= {P1, . . . , Pm}, π2= {Q1, . . . , Qr} będą podziałami kostki P . Powiemy, że π2jest wpisany w π1 lub też, że π2 jest zagęszczeniem π1, jeżeli dla dowolnego j ∈ {1, . . . , m} rodzina {Qk : Qk ⊂ Pj} jest podziałem Pj 5 . W tej sytuacji piszemy π24 π1.
Średnicą podziału π = {P1, . . . , Pm} nazywamy liczbę
diam π := max{diam P1, . . . , diam Pm}.
Niech (πk)∞k=1będzie ciągiem podziałów kostki P . Powiemy, że jest to normalny ciąg podziałów, jeżeli diam πk −→ 0.
Obserwacja 3.1.2. (a) Relacja4 jest przechodnia.
(b) Dla dowolnych podziałów π1, π2kostki P istnieje podział prosty π taki, że π4 πj, j = 1, 2 6 . (c) Dla dowolnego podziału π = {P1, . . . , Pm} mamy
|P | = |P1| + · · · + |Pm|.
1 Rozważamy tylko kostki niezdegenerowane.
2 Zauważmy, że |P | > 0 oraz, że |x0+ P | = |P |, x0∈ Rn.
3 Dla n = 1 podział kostki P = [a, b] możemy utożsamiać z ciągiem π = (x0, . . . , xm), gdzie a = x0< x1 < · · · <
xm= b.
4 Podział ten składa się z m1. . . mnkostek. Typowym przypadkiem jest sytuacja, gdy xj,k:= aj+ (k/mj)(bj− aj) (tzn. dzielimy krawędzie na równe części).
5 Innymi słowy, jeżeli int Pj∩ int Qk6= ∅, to Qk⊂ Pj.
6 π jest wspólnym zagęszczeniem podziałów π1 i π2. 117
3. Całka Riemanna Istotnie, najpierw sprawdzamy przypadek podziału prostego:
m1
Dla dowolnego podziału znajdujemy najpierw wpisany podział prosty π0 = {Q1, . . . , Qr},
Definicja 3.1.3. Niech f : P −→ R będzie dowolną funkcją ograniczoną na kostce P . Zdefiniujmy:
m(f, P ) := inf f (P ), M (f, P ) := sup f (P ).
Liczbę L(f, π) nazywamy sumą aproksymacyjną dolną dla funkcji f przy podziale π. Analogicznie, U (f, π) nazywamy sumą aproksymacyjną górną 7 . Zauważmy, że
m(f, P )|P | 6 L(f, π) 6 U (f, π) 6 M (f, P )|P |.
gdzie supremum i infimum bierzemy po wszystkich podziałach kostki P . LiczbęR
∗Pf nazywamy całką dolną z funkcji f po kostce P . Analogicznie, liczbęR∗
Pf nazywamy całką górną.
Powiemy, że funkcja f jest całkowalna w sensie Riemanna na kostce P 8
(f ∈ R(P )), jeżeli R
∗Pf =R∗
Pf . Wtedy wspólną wartość tych całek oznaczamy przez R
Pf i nazywamy całką Riemanna z funkcji f po kostce P .
Niech teraz ϕ : P −→ C będzie funkcją ograniczoną. Powiemy, że ϕ jest całkowalna w sensie Rie-manna na P (ϕ ∈ R(P, C)), jeżeli Re ϕ, Im ϕ ∈ R(P ). Kładziemy wtedy RPϕ := R
PRe ϕ + iR
PIm ϕ i nazywamy tę liczbę całką Riemanna z funkcji ϕ po kostce P .
Oczywiście każda funkcja stała c ∈ C jest całkowalna w sensie Riemanna iR
Pc = c|P |.
Poniżej przedstawimy szereg (mniej lub bardziej elementarnych) własności całki Riemanna; f, g : P −→ R, ϕ, ψ : P −→ C oznaczają funkcje ograniczone, π, π1, π2 — podziały kostki P .
7 Czasami mówi się o sumach Darboux.
8 Bernhard Riemann (1826–1866).
9 Przypomnijmy, że χAoznacza funkcję charakterystyczną zbioru A ⊂ P . Funkcję f nazywamy funkcją Dirichleta w kostce P .
10 Peter Dirichlet (1805–1859).
3.1. Całka Riemanna na kostce Własność 3.1.5. Dla dowolnego c ∈ R mamy:
L(f + c, π) = L(f, π) + c|P | oraz U (f + c, π) = U (f, π) + c|P |.
Wynik przenosi się natychmiast na funkcje ograniczone ϕ : P −→ C i c ∈ C:
ϕ ∈R(P, C) ⇐⇒ ϕ + c ∈ R(P, C) oraz Z
Teraz wystarczy jako π wziąć wspólne zagęszczenie π1 i π2 i skorzystać w Własności 3.1.8.
(⇐=): Ustalmy ε i niech π będzie podziałem takim jak w warunku. Wtedy Z ∗
P
f − Z
∗Pf 6 U (f, π) − L(f, π) 6 ε,
co, wobec dowolności ε, daje całkowalność.
11 Wystarczy wykorzystać poprzednie nierówności dla π i πj, gdzie π jest wspólnym zagęszczeniem podziałów π1
i π2.
3. Całka Riemanna
Definicja 3.1.10. Niech π = {P1, . . . , Pm} będzie podziałem kostki P i niech ξ = {ξ1, . . . , ξm}, ξj ∈ Pj, j = 1, . . . , m. Niech ϕ : P −→ C. Sumę
M (ϕ, π, ξ) :=
m
X
j=1
ϕ(ξj)|Pj|
nazywamy sumą aproksymacyjną pośrednią dla funkcji ϕ przy podziale π i punktach pośrednich ξ 12 . (a) Jest rzeczą widoczną, iż M (αϕ + βψ, π, ξ) = αM (ϕ, π, ξ) + βM (ψ, π, ξ), α, β ∈ C.
(b) W szczególności, M (ϕ, π, ξ) = M (Re ϕ, π, ξ) + iM (Im ϕ, π, ξ).
(c) Ponadto, |M (ϕ, π, ξ)|6 M (|ϕ|, π, ξ).
(d) Dla f : P −→ R mamy: L(f, π) 6 M (f, π, ξ) 6 U (f, π).
Własność 3.1.11. Dla każdego normalnego ciągu podziałów (πk)∞k=1 mamy:
L(f, πk) −→
Z
∗P
f, U (f, πk) −→
Z ∗ P
f.
Dowód . Ograniczymy się do sum górnych. Można założyć, że f > 0 (zastępując f przez f + c). Weźmy ε > 0 i niech π = {Q1, . . . , Qm} będzie podziałem takim, że U (f, π) −R∗
Pf 6 ε. Niech πk = {Pk,1, . . . , Pk,mk}, k > 1.
Wtedy
U (f, πk) = X
j∈{1,...,mk}
∃i∈{1,...,m}:Pk,j⊂Qi
M (f, Pk,j)|Pk,j| + X
j∈{1,...,mk}
∀i∈{1,...,m}:Pk,j6⊂Qi
M (f, Pk,j)|Pk,j| 6 U(f, π) + M(f, P )ηk,
gdzie
ηk := X
j∈{1,...,mk}
∀i∈{1,...,m}: Pk,j6⊂Qi
|Pk,j|, k > 1.
Zauważmy, że ηk −→ 0. Istotnie, jeżeli Pk,j 6⊂ Qi, i = 1, . . . , m, to Pk,j musi przecinać którąś ze ścian którejś z kostek Q1,. . . , Qm. Stąd ηk 6 c diam πk, k > 1, gdzie c > 0 jest pewną stałą (Ćwiczenie).
Ostatecznie
Z ∗ P
f 6 lim inf
k→+∞U (f, πk) 6 lim sup
k→+∞
U (f, πk) 6 Z ∗
P
f + ε,
co, wobec dowolności ε > 0, kończy dowód.
Własność 3.1.12. Następujące warunki są równoważne:
(i) ϕ ∈R(P, C);
(ii) istnieje c ∈ C takie, że dla dowolnego normalnego ciągu podziałów (πk)∞k=1 oraz dla dowolnego wyboru punktów pośrednich (ξk)∞k=1 13
mamy M (ϕ, πk, ξk) −→ c 14 ;
(ii’) istnieje c ∈ C takie, że dla dowolnego ε > 0 istnieje δ > 0 taka, że dla dowolnego podziału π o średnicy 6 δ oraz dla dowolnego wyboru punktów pośrednich ξ mamy |M (ϕ, π, ξ) − c| 6 ε;
(iii) istnieje c ∈ C oraz normalny ciągu podziałów (πk)∞k=1 takie, że dla dowolnego wyboru punktów pośrednich (ξk)∞k=1, mamy M (ϕ, πk, ξk) −→ c 15 .
12 Czasami mówimy o sumie Cauchy’ego–Riemanna.
13 Tzn. ξkjest zbiorem punktów pośrednich dla πk.
14
Jest rzeczą oczywistą, że jeżeli mamy dwa ciągi zbieżne liczb rzeczywistych (ak)∞k=1 i (bk)∞k=1 oraz wiemy, że ciąg a1, b1, a2, b2, . . . jest również zbieżny, to lim
k→+∞ak = lim
k→+∞bk. Z tej trywialnej uwagi wynika, że warunek (ii) jest równoważny następującemu warunkowi:
dla dowolnego normalnego ciągu podziałów (πk)∞k=1 oraz dla dowolnego wyboru punktów pośrednich (ξk)∞k=1, ciąg (M (ϕ, πk, ξk))∞k=1 jest zbieżny do granicy skończonej.
15 Podobnie jak poprzednio, warunek (iii) jest równoważny następującemu warunkowi: istnieje normalny ciąg po-działów (πk)∞k=1taki, że dla dowolnego wyboru punktów pośrednich (ξk)∞k=1, ciąg (M (ϕ, πk, ξk))∞k=1jest zbieżny do granicy skończonej.
3.1. Całka Riemanna na kostce
Dowód . Wobec Własności 3.1.5 oraz (b), możemy założyć, że ϕ = f : P −→ R.
(i) =⇒ (ii) wynika z Własności 3.1.11 i nierówności (d).
Równoważność (ii) ⇐⇒ (ii’) jest elementarna (Ćwiczenie).
Implikacja (ii) =⇒ (iii) jest trywialna.
Dla dowodu implikacji (iii) =⇒ (i) wystarczy zauważyć, że istnieją ciągi punktów pośrednich (ξk0)∞k=1, (ξk00)∞k=1 takie, że
M (f, πk, ξk0) − L(f, πk) 6 1
k, U (f, πk) − M (f, πk, ξk00) 6 1
k, k > 1. 16 Wtedy, na podstawie Własności 3.1.11,
M (f, πk, ξk0) −→
Z
∗P
f, M (f, πk, ξ00k) −→
Z ∗ P
f, a zatemR
∗Pf =R∗
Pf , czyli f ∈R(P ).
Własność 3.1.13. R(P, C) jest zespoloną przestrzenią wektorową, a operator R(P, C) 3 ϕ 7−→Z
Pϕ ∈ C jest C–liniowy.
Dowód . Dowód wynika natychmiast z (a) oraz Własności 3.1.12.
Własność 3.1.14. Jeżeli ψ ∈R(P, C) oraz
|ϕ(x0) − ϕ(x00)| 6 |ψ(x0) − ψ(x00)|, x0, x00∈ P, to ϕ ∈R(P, C).
Dowód . Ponieważ |x + iy| 6 |x| + |y|, zatem
U (Re ϕ, π) − L(Re ϕ, π) 6 U (Re ψ, π) − L(Re ψ, π) + U (Im ψ, π) − L(Im ψ, π), co, wobec Własności 3.1.9, daje całkowalność Re ϕ. Podobnie postępujemy z Im ϕ.
Własność 3.1.15. Jeżeli ϕ ∈R(P, C), to |ϕ| ∈ R(P ) oraz
Z
P
ϕ 6
Z
P
|ϕ|.
Dowód . Całkowalność wynika z Własności 3.1.14. Wzór wynika z (c). Własność 3.1.16. Jeżeli ϕ, ψ ∈R(P, C), to ϕ · ψ ∈ R(P, C).
Dowód . Wystarczy wykazać, że f · g ∈ R(P ) dla f, g ∈ R(P ) (Ćwiczenie). Ustalmy f, g ∈ R(P ) i niech
|f |, |g| 6 C. Mamy:
|f (x0)g(x0) − f (x00)g(x00)| 6 C(|f (x0) − f (x00)| + |g(x0) − g(x00)|), x0, x00∈ P.
Teraz możemy użyć rozumowania takiego, jak w dowodzie Własności 3.1.14. Własność 3.1.17. Operator
R(P, C) × R(P, C) 3 (ϕ, ψ)7−→Φ Z
P
ϕψ ∈ C
jest półtoraliniowy, Φ(ψ, ϕ) = Φ(ϕ, ψ) oraz Φ(ϕ, ϕ) > 0. W szczególności, na podstawie nierówności Schwarza, mamy
Z
P
ϕψ
2
6Z
P
|ϕ|2Z
P
|ψ|2
, ϕ, ψ ∈R(P, C).
Ponadto, funkcjaR(P, C) 3 ϕ 7−→ (RP|ϕ|2)1/2 jest seminormą 17 . Własność 3.1.18. C(P, C) ⊂ R(P, C).
16 Np. ξ0k,jwybieramy tak, by f (ξk,j0 ) − m(f, Pk,j) 6 1/(k|P |).
17 Tzn. spełnia wszystkie warunki normy z wyjątkiem tego, że może się zerować na niezerowych wektorach.
3. Całka Riemanna
Dowód . Wystarczy udowodnić, że C(P ) ⊂ R(P ). Ustalmy ε > 0. Wobec jednostajnej ciągłości funkcji f na P , istnieje δ > 0 takie, że |f (x0) − f (x00)| 6 ε o ile kx0− x00k 6 δ. Niech π będzie podziałem P o średnicy6 δ. Wtedy U (f, π) − L(f, π) 6 ε|P |. Na mocy Własności 3.1.9 dowód jest zakończony. Własność 3.1.19 (Twierdzenie o wartości średniej). Dla dowolnej funkcji f ∈ C(P ) istnieje ξ ∈ P taki, że
f (ξ) = 1
|P | Z
P
f.
Dowód . Wobec nierówności
min f (P ) 6 1
|P | Z
P
f 6 max f (P ),
wystarczy skorzystać z zasady Darboux.
Definicja 3.1.20. Mówimy, że zbiór A ⊂ Rn ma objętość zero (|A| = 0), jeżeli dla dowolnego ε > 0 istnieje skończona rodzina kostek P1, . . . , Pmtaka, że A ⊂ P1∪ · · · ∪ Pmi |P1| + · · · + |Pm| 6 ε.
Obserwacja 3.1.21. (a) Każdy zbiór o objętości zero jest ograniczony.
(b) Każdy zbiór skończony ma objętość zero.
(c) Podzbiór zbioru o objętości zero ma objętość zero.
(d) Suma skończonej liczby zbiorów o objętości zero ma objętość zero.
(e) Jeżeli Rn 3 aν−→ a0∈ Rn, to zbiór {aν : ν > 0} ma objętość zero.
(f) Jeżeli A ⊂ Rn ma objętość zero, to dla dowolnego zbioru ograniczonego B ⊂ Rm zbiór A × B ma objętość zero. W szczególności, dla dowolnej kostki ograniczonej Q ⊂ Rn 18
mamy |∂Q| = 0.
(g) Jeżeli |A| = 0, to |A| = 0.
Propozycja 3.1.22. (a) Niech Q ⊂ Rd będzie kostką domkniętą i niech ϕ : Q −→ Rn−d będzie dowol-nym odwzorowaniem ciągłym, 16 d 6 n − 1. Wtedy wykres A := {(t, ϕ(t)) : t ∈ Q} ma objętość zero.
(b) Niech Q ⊂ Rd będzie kostką domkniętą i niech p : Q −→ Rn będzie odwzorowaniem spełniającym warunek Höldera z wykładnikiem α. Wówczas:
(b1) Jeżeli αn > d, to zbiór A := p(Q) ma objętość zero. 19
(b2) Jeżeli αn > d, to zbiór A := p(Z) ma objętość zero dla dowolnego zbioru Z ⊂ Q o objętości zero.
(b3) Jeżeli p spełnia warunek Lipschitza (np. p ∈ C1(Ω, Rn), gdzie Ω jest otoczeniem Q), to:
• jeżeli n > d, to |p(Q)| = 0,
• jeżeli n> d, to |p(Z)| = 0 dla dowolnego zbioru Z ⊂ Q o objętości zero. 20
(c) Jeżeli M jest d–wymiarową podrozmaitością w Rn klasy C1, 06 d 6 n − 1, to każdy zwarty podzbiór A ⊂ M ma objętość zero 21 .
Dowód . (a) Ustalmy ε > 0. Wobec jednostajnej ciągłości odwzorowania ϕ istnieje δ > 0 taka, że max
j=1,...,n−d|ϕj(t0) − ϕj(t00)| 6 ε o ile kt0− t00k 6 δ.
Niech π = {Q1, . . . , Qm} będzie podziałem kostki Q takim, że diam π 6 δ. Ustalmy tj ∈ Qj i niech Pj := Qj× (ϕ(tj) + [−ε, ε]n−d).
Oczywiście A ⊂ P1∪ · · · ∪ Pmoraz
|P1| + · · · + |Pm| = |Q1|(2ε)n−d+ · · · + |Qm|(2ε)n−d= |Q|(2ε)n−d.
18 Q jest (niezdegenerowaną) kostką ograniczoną, jeżeli Q jest kostką domkniętą i int Q ⊂ Q.
19 Warto w tym miejscu przypomnieć o istnieniu krzywej γ : [0, 1] −→ R2 takiej, że γ([0, 1]) = [0, 1] × [0, 1]. Z (b1) wynika, że γ nie może spełniać warunku Höldera z wykładnikiem α > 1/2.
20 Jeżeli n < d, to biorąc jako p projekcję pr
Rd−n, zaś jako Z zbiór postaci A × B ⊂ Rn× Rd−n, gdzie |A| = 0, a B jest ograniczony, widzimy, że p(Z) = B nie musi być zbiorem o objętości zero nawet, gdy p jest liniowe.
21 Dla przykładu, |Sn−1| = 0.
3.1. Całka Riemanna na kostce (b) Przypuśćmy, że
j=1,...,nmax |pj(t0) − pj(t00)| 6 Ckt0− t00kα, t0, t00∈ Q.
Ustalmy ε > 0. Przypadki (b1) i (b2) rozpatrzymy jednocześnie. Niech X := Q (odp. X := Z) i niech κ := 2|Q| (odp. κ > 0, gdzie κ jest dowolne).
Niech π = {Q1, . . . , Qm} będzie układem kostek takim, że
• Qj = uj+ [−η, η]d, gdzie uj∈ Rd, j = 1, . . . , m, η ∈ (0, 1],
• X ⊂ Q1∪ · · · ∪ Qm, Qj∩ Q 6= ∅, j = 1, . . . , m,
• m(2η)d = |Q1| + · · · + |Qm| 6 κ.
Niech δ := C(2√
dη)α. Ustalmy tj ∈ Qj∩ Q i niech
Pj := p(tj) + [−δ, δ]n, j = 1, . . . , m.
Oczywiście A ⊂ P1∪ · · · ∪ Pmoraz
|P1| + · · · + |Pm| = m(2C(2√
dη)α)n6 2(α+1)n−dκ(Cdα/2)nηαn−d.
Gdy αn > d, to dla małych η > 0 ostatnia liczba będzie dowolnie mała (odp. gdy αn = d, to dobierając κ > 0 stosownie małe, możemy uczynić tę ostatnią liczbę dowolnie małą).
(b3) wynika z (b1) i (b2).
(c) wynika z (b).
Własność 3.1.23. Jeżeli zbiór
NP(ϕ) := {a ∈ P : ϕ nie jest ciągła w punkcie a}
ma objętość zero, to ϕ ∈R(P, C) (por. Własność 3.1.18).
Dowód . Wystarczy rozważyć przypadek, gdy ϕ = f : P −→ R. Ustalmy ε > 0 i niech P1, . . . , Pm ⊂ P będą kostkami takimi, że NP(f ) ⊂ intPP1∪ · · · ∪ intPPm 22 , |P1| + · · · + |Pm| 6 ε. Niech K :=
P \
m
S
j=1
intPPj. Zbiór K jest zwarty i f jest ciągła na K. Zatem jest jednostajnie ciągła. Niech δ > 0 będzie takie, że |f (x0) − f (x00)| 6 ε dla x0, x00∈ K, kx0− x00k 6 δ. Rozważmy podział π = {Q1, . . . , Qr} kostki P taki, że diam π6 δ, Q1, . . . , Qs⊂ K, Qs+1, . . . , Qr⊂
m
S
j=1
Pj oraz dla dowolnego j ∈ {1, . . . , m}
rodzina {Qk: Qk ⊂ Pj} stanowi podział Pj. Wtedy
U (f, π) − L(f, π) 6 ε(|Q1| + · · · + |Qs|) + (M (f, P ) − m(f, P ))(|Qs+1| + · · · + |Qr|)
6 ε(|P | + M (f, P ) − m(f, P )). 23
Własność 3.1.24. Niech π = {P1, . . . , Pm} będzie ustalonym podziałem. Wtedy ϕ ∈R(P, C) ⇐⇒ ϕ|Pj ∈ R(Pj, C), j = 1, . . . , m. Ponadto,
Z
P
ϕ = Z
P1
ϕ + · · · + Z
Pm
ϕ.
Dowód . Wystarczy rozważyć przypadek, gdy ϕ = f : P −→ R. Niech π0 będzie podziałem P powstałym przez „zsumowanie” podziałów π1, . . . , πm, gdzie πj jest podziałem Pj, j = 1, . . . , m. Wtedy
U (f, π0) − L(f, π0) =
m
X
j=1
U (f |Pj, πj) − L(f |Pj, πj) , co daje równoważność całkowalności na P i na każdym Pj, j = 1, . . . , m.
Ponadto,
m
P
j=1
L(f |Pj, πj) = L(f, π0) 6 R
Pf , skąd wynika, że
m
P
j=1
R
Pjf 6 R
Pf . Rozumowanie na
sumach górnych daje przeciwną nierówność.
22
intPA oznacza wnętrze w topologii indukowanej.
23 Nierówność |Qs+1| + · · · + |Qr| 6 |P1| + · · · + |Pm| wynika z tego, że kostki Qkdzielą kostki Pj.
3. Całka Riemanna
Ponadto, na podstawie Własności 3.1.24, mamy:
jest relacją równoważnościową wR(P, C); całka Riemanna jest dobrze określonym operatorem liniowym R(P, C)/∼ −→ C.
Wynika stąd całkowalność funkcji f . Zbieżność całek wynika bezpośrednio z nierówności
jest iloczynem skalarnym. Niestety przestrzeń C(P, C) z tym iloczynem skalarnym nie jest przestrzenią Hilberta. Istotnie, niech P := [0, 1] ⊂ R,
24 Por. dowód Własności 3.1.23.
3.1. Całka Riemanna na kostce f0:= χ(1/2,1]. WtedyR
P(fν− f0)2−→ 0 (Ćwiczenie). W szczególności, (fν)∞ν=1jest ciągiem Cauchy’ego w C(P, C), który nie ma granicy 25 .
Własność 3.1.30. Jeżeli 06 f ∈ R(P ) iR
Pf = 0, to zbiór Zf := {x ∈ P : f (x) > 0} jest przeliczalną sumą zbiorów o objętości zero.
Dowód . Wystarczy pokazać, że dla dowolnego c > 0 zbiór A := {f > c} ma objętość zero. Weźmy ε > 0 i podział π = {P1, . . . , Pm} taki, że U (f, π) 6 ε. Wtedy
ε > U (f, π) > c X
j∈{1,...,m}:
A∩Pj6=∅
|Pj|.
Wynika stąd, że A można pokryć kostkami o łącznej objętości6 ε/c. Przykład 3.1.31. Zbiór Zf := {f > 0} we Własności 3.1.30 może nie mieć objętości zero. Dla przykładu, niech P ∩ Qn = {r1, r2, . . . } i niech f : P −→ [0, 1],
f (x) :=
(0, x /∈ Qn
1
j, x = rj . Oczywiście, zbiór Zf = P ∩ Qn nie ma objętości zero, ale R
Pf = 0. Istotnie, dla dowolnego k ∈ N rozważmy podział π = {P1, . . . , Pm} (m > k) taki, że rj ∈ intPPj, j = 1, . . . , k, |P1| + · · · + |Pk| 6 k1. Wtedy
U (f, π) =
k
X
j=1
M (f, Pj)|Pj| +
m
X
j=k+1
M (f, Pj)|Pj| 6
k
X
j=1
|Pj| + 1 k + 1
m
X
j=k+1
|Pj| 6 1 k+ 1
k + 1|P |.
Definicja 3.1.32. Mówimy, że zbiór ograniczony A ⊂ Rn jest mierzalny w sensie Jordana, jeżeli χA∈ R(P ) dla pewnej kostki P ⊃ A. Liczbę |A| := RPχA nazywamy miarą Jordana (objętością) zbioru A
26.
Własność 3.1.33. Jeżeli A jest zbiorem ograniczonym takim, że |∂A| = 0 w sensie Definicji 3.1.20, to A jest mierzalny w sensie Jordana. Ponadto, jeżeli |A| = 0 w sensie Definicji 3.1.20, to |A| = 0 w sensie Definicji 3.1.32.
Dowód . Wystarczy skorzystać z Własności 3.1.23.
Obserwacja 3.1.34. Sumy aproksymacyjne pośrednie można zdefiniować dla dowolnego odwzorowania ϕ : P −→ F , gdzie F jest przestrzenią unormowaną. Pozwala to przenieść pojęcie całki Riemanna:
odwzorowanie ϕ : P −→ F nazywamy całkowalnym w sensie Riemanna (ϕ ∈ R(P, F )), jeżeli istnieje c ∈ F takie, że dla dowolnego normalnego ciągu podziałów (πk)∞k=1oraz dla dowolnego wyboru punktów pośrednich mamy: M (ϕ, πk, ξk) −→ c. Element c nazywamy wtedy całką Riemanna odwzorowania ϕ po kostce P i oznaczamy R
Pϕ. Własność 3.1.12 gwarantuje zgodność definicji dla F = C. Pojawia się tu pewna subtelność: nowa definicja całki obejmuje formalnie funkcje nieograniczone, a Własność 3.1.12 dotyczy tylko funkcji ograniczonych. Dla usunięcia tego problemu wystarczy zauważyć, że jeżeli funkcja ϕ : P −→ F jest całkowalna w nowym sensie, to musi być ograniczona. Istotnie, przypuśćmy np. że supPkϕk = +∞ i niech P 3 aν −→ a0 ∈ P będzie ciągiem takim, że kϕ(aν)k −→ +∞. Weźmy normalny ciąg podziałów (πk)∞k=1, πk = {Pk,1, . . . , Pk,mk}, taki, że a0 ∈ intPPk,1, k > 1. Ustalmy k.
Wybierzmy w sposób dowolny punkty pośrednie ξk,2, . . . , ξk,mk. Ponieważ supPk,1kϕk = +∞, to zawsze znajdziemy punkt ξk,1taki, że kM (ϕ, πk, ξk)k > k. Mamy więc kM (ϕ, πk, ξk)k −→ +∞, czyli ϕ nie może być całkowalna w nowym sensie.
25 Przypuśćmy, że granica taka istnieje i jest nią funkcja g0 ∈ C(P, C). Mamy R
P|f0 − g0|2 = 0. Stąd, wobec Obserwacji 3.1.29, g0= 0 na [0,12) i g0= 1 na (12, 1]; sprzeczność.
26 Łatwo widać (por. Własność 3.1.24), że jeżeli χA∈R(P ) dla pewnej kostki P ⊃ A, to χA∈R(Q) dla dowolnej kostki Q ⊃ A. Ponadto,R
PχA=R
QχA.
3. Całka Riemanna
3.2. Całka Riemanna na zbiorze regularnym
Niech A ⊂ Rn będzie dowolnym zbiorem ograniczonym takim, że |∂A| = 0 — takie zbiory bę-dziemy roboczo nazywać regularnymi. Naszym celem jest zbudowanie teorii całki Riemanna R
Af . Dla prostoty ograniczymy się do funkcji ograniczonych f : A −→ R. Przeniesienie teorii na funkcje zespolone pozostawiamy jako Ćwiczenie.
Ćwiczenie 3.2.1. (a) Jeżeli |A| = 0, to A jest regularny.
(b) Każda kostka jest regularna.
(c) Jeżeli zbiór A jest regularny, to zbiór int A jest regularny; implikacja przeciwna nie jest prawdziwa.
(d) Jeżeli zbiór A jest regularny, to zbiór A jest regularny; implikacja przeciwna nie jest prawdziwa.
(e) Jeżeli zbiory A, B są regularne, to zbiory A ∪ B, A ∩ B, A \ B są regularne.
(f) Jeżeli zbiory A ⊂ Rn, B ⊂ Rmsą regularne, to zbiór A×B jest regularny (por. Obserwacja 3.1.21(f)).
Definicja 3.2.2. Niech A ⊂ P , gdzie P jest kostką. Dla f : A −→ R niech f0:=
(f, na A 0, na P \ A.
Powiemy, że f jest całkowalna w sensie Riemanna na A (f ∈R(A)), jeżeli f0∈R(P ). Kładziemy wtedy R
Af :=R
Pf0. Bez trudu widać, że taka definicja jest poprawna, tzn. nie zależy od wyboru kostki P ⊃ A.
Poniżej, przedstawimy listę podstawowych własności tak zdefiniowanej całki Riemanna. Dowody poszczególnych własności polegają w większości przypadków na wykorzystaniu własności całki Riemanna na kostce i dlatego też podamy tylko szkice dowodów.
Obserwacja 3.2.3. (a) Jeżeli |A| = 0, to każda funkcja ograniczona f : A −→ R jest całkowalna oraz R
Af = 0 (por. Ćwiczenie 3.2.1(a)).
Istotnie, zbiór {x ∈ P : f0(x) 6= 0} ⊂ A i wystarczy skorzystać z Własności 3.1.25.
(b) Jeżeli zbiór NA(f ) := {a ∈ A : f nie jest ciągła w a} ma objętość zero, to f ∈ R(A) (jest to uogólnienie Własności 3.1.23).
Istotnie, wtedy NP(f0) = {a ∈ P : f0 nie jest ciągła w punkcie a} ⊂ NA(f ) ∪ ∂A i możemy skorzystać z Własności 3.1.23.
W szczególności, BC(A) ⊂R(A).
(c) 1 ∈R(A) oraz RA1 = |A|, gdzie |A| oznacza miarę Jordana (objętość) zbioru A w sensie Definicji 3.1.32.
(d) R(A) jest algebrą, operator R(A) 3 f 7−→ RAf ∈ R jest liniowy i monotoniczny. W szczególności, jeżeli A ⊂ B, gdzie B jest także regularny, to |A|6 |B|.
(e) Jeżeli f ∈R(A), to |f| ∈ R(A) oraz | RAf | 6R
A|f |.
(f) OperatorR(A) × R(A) 3 (f, g) 7−→ RAf g ∈ R jest dwuliniowy, symetryczny i nieujemny. W szcze-gólności, zachodzi nierówność Schwarza:
Z
A
f g2
6Z
A
f2Z
A
g2
, f, g ∈R(A).
(g) Jeżeli zbiór {a ∈ A : f (a) 6= g(a)} ma objętość zero, to f ∈R(A) wtedy i tylko wtedy, gdy g ∈ R(A).
Ponadto,R
Af =R
Ag (jest to uogólnienie Własności 3.1.25).
(h) Jeżeli B ⊂ A jest zbiorem regularnym, to dla funkcji f : A −→ R mamy: f |B ∈R(B) ⇐⇒ f · χB ∈ R(A).
(i) Jeżeli A = A1∪ A2, gdzie A1, A2są zbiorami regularnymi i |A1∩ A2| = 0, to f ∈R(A) wtedy i tylko wtedy, gdy f |Aj ∈R(Aj), j = 1, 2. PonadtoR
Af =R
A1f +R
A2f .
Istotnie, jeżeli f ∈ R(A), to f · χAj ∈R(A) 27 , a zatem f |Aj ∈R(Aj), j = 1, 2. Odwrotnie, jeżeli f |Aj ∈R(Aj), to fj := f · χAj ∈R(A), j = 1, 2, a zatem f1+ f2∈R(A). Teraz wystarczy tylko zauważyć, że {f1+ f26= f } ⊂ A1∩ A2i skorzystać z (g).
W szczególności (por. Ćwiczenie 3.2.1):
• jeżeli B ⊂ A jest zbiorem regularnym i f ∈R(A), to f|B∈R(B);
• dla dowolnej funkcji ograniczonej f : A −→ R mamy: f ∈ R(A) ⇐⇒ f ∈ R(A) orazR
Af =R
Af ;
27 Bo χAj∈R(A) i R(A) jest algebrą (por. (d)).
3.3. Własności całki Riemanna
• dla dowolnej funkcji ograniczonej f : A −→ R mamy: f ∈ R(A) ⇐⇒ f ∈ R(int A) oraz R
Af =R
int Af .
Propozycja 3.2.4. Niech A ⊂ Rn−1 będzie zbiorem regularnym, niech f ∈ R(A) i niech B oznacza wykres funkcji f , tzn. zbiór
B := {(x, f (x)) : x ∈ A} ⊂ Rn. Wtedy |B| = 0.
Dowód . Niech A ⊂ P , gdzie P jest kostką, i niech f0 oznacza standardowe przedłużenie funkcji f . Zauważmy, że |B| = 0 ⇐⇒ |B0| = 0 (Ćwiczenie). Wynika stąd w szczególności, że możemy założyć, że A = P jest kostką.
Dla ε > 0, niech π = {P1, . . . , Pm} będzie podziałem kostki P takim, że U (f, π) − L(f, π) 6 ε.
Połóżmy, Qj:= Pj×[m(f, Pj), M (f, Pj)], j = 1, . . . , m. Wtedy B ⊂ Q1∪· · ·∪Qmoraz |Q1|+· · ·+|Qm| =
U (f, π) − L(f, π) 6 ε.
Twierdzenie 3.2.5 (Twierdzenie o całkach iterowanych). Niech A ⊂ Rn, B ⊂ Rm będą zbiorami re-gularnymi i niech f ∈ R(A × B) (por. Ćwiczenie 3.2.1(f)). Załóżmy, że f(x, ·) ∈ R(B) dla dowolnego x ∈ A 28 . Wtedy funkcja
A 3 x 7−→
Z
Bf (x, y)dy ∈ R jest całkowalna na A oraz
Z
A×B
f (x, y)dxdy = Z
A
Z
B
f (x, y)dy
dx. 29 W szczególności, powyższy wzór stosuje się dla funkcji klasy BC(A × B).
Dowód . W oczywisty sposób dowód sprowadza się do przypadku, gdy A = P i B = Q są kostkami.
Ustalmy dowolne podziały π0= {P1, . . . , Pr} kostki P i π00= {Q1, . . . , Qs} kostki Q. Wtedy rodzina π := {Pj× Qk: j = 1, . . . , r, k = 1, . . . , s}
jest podziałem kostki P × Q. Wybierzmy dowolne punkty pośrednie ξ dla podziału π0. Wtedy m(f, Pj× Qk)|Qk| 6
Z
Qk
f (ξj, y)dy 6 M (f, Pj× Qk)|Qk|, 30
skąd po pomnożeniu przez |Pj|, zsumowaniu najpierw względem k, a potem względem j, mamy:
L(f, π) 6
r
X
j=1
Z
Q
f (ξj, y)dy|Pj| 6 U(f, π).
Pozostaje rozważyć normalne ciągi podziałów i zastosować Własność 3.1.12. 3.3. Własności całki Riemanna
Propozycja 3.3.1. Niech A ⊂ Rn, B ⊂ Rm będą zbiorami regularnymi i niech f ∈ R(A), g ∈ R(B).
Wtedy f ⊗ g ∈R(A × B) 31 oraz Z
A×B
f (x)g(y)dxdy =Z
A
f (x)dxZ
B
g(y)dy .
Dowód . Jedynym problemem jest całkowalność f ⊗ g na A × B. Możemy założyć, że A = P , B = Q są kostkami. Ustalmy dowolne podziały π0= {P1, . . . , Pr} kostki P i π00= {Q1, . . . , Qs} kostki Q i niech π będzie jak w dowodzie Twierdzenia 3.2.5. Załóżmy, że |f |, |g|6 c. Wtedy
U (f ⊗ g, π) − L(f ⊗ g, π) 6 c
(U (f, π0) − L(f, π0))|Q| + (U (g, π00) − L(g, π00))|P |
i dalej możemy rozumować standardowo.
28 Odnotujmy, że nie wynika to z całkowalności f na A × B. Dla przykładu: A = B = [0, 1], f (0, ·) = χQ∩[0,1], f (x, ·) := 0 dla x ∈ (0, 1].
29 Stosujemy tu wygodny tradycyjny zapis całki Riemanna z uwidocznieniem zmiennych.
30 Z całkowalności funkcji f (ξj, ·) na Q wynika jej całkowalność na Qk.
31 (f ⊗ g)(x, y) := f (x)g(y)
3. Całka Riemanna
Propozycja 3.3.2 (Całka jako miara objętości). Niech Q będzie kostką w Rn−1 i niech α, β ∈ C(Q), α 6 β. Połóżmy
A := {(x, y) ∈ Q × R : α(x) 6 y 6 β(x)}.
Wtedy A jest regularnym zbiorem zwartym oraz dla dowolnej funkcji f ∈ C(A) funkcja Q 3 x 7−→
Dowód . Regularność zbioru wynika z Propozycji 3.1.22(a). Wybierzmy c > 0 takie, że −c 6 α 6 β 6 c.
Niech P := Q × [−c, c]. Wtedy na podstawie twierdzenia o całkach iterowanych mamy:
Z
Twierdzenie 3.3.3 (Twierdzenie o funkcjach danych całką). Niech Ω będzie zbiorem otwartym w Rm, niech P ⊂ Rn będzie kostką i niech f : Ω × P −→ R będzie odwzorowaniem takim, że dla pewnego dla x ∈ B(a, δ) i t ∈ P . Stąd, na podstawie twierdzenia o przyrostach skończonych, otrzymujemy:
Twierdzenie 3.3.4 (Twierdzenie o funkcjach danych całką niewłaściwą). Niech Ω będzie zbiorem otwar-tym w Rm, niech −∞ < a < b6 +∞ i niech f : Ω × [a, b) −→ R będzie odwzorowaniem takim, że dla
3.3. Własności całki Riemanna
• dla dowolnego |α|6 k istnieje odwzorowanie ciągłe gα: [a, b) −→ R takie, że |Dα(ft)(x)| 6 gα(t) dla x ∈ Ω i t ∈ [a, b) oraz takie, że całka niewłaściwaRb
agα(t)dt jest zbieżna.
Wtedy odwzorowanie
ϕ(x) :=
Z b a
f (x, t)dt, x ∈ Ω, jest klasy Ck(Ω) oraz
Dαϕ(x) = Z b
a
Dα(ft)(x)dt, x ∈ Ω, |α| 6 k. 34
Dowód . Ustalmy ciąg a < bν < b, bν % b i niech ϕν(x) :=
Z bν a
f (x, t)dt, x ∈ Ω.
Na podstawie poprzedniego twierdzenia wnioskujemy, że ϕν ∈ Ck(Ω) oraz Dαϕν(x) =
Z bν a
Dα(ft)(x)dt, x ∈ Ω, |α| 6 k.
Zauważmy, że ϕν−→ ϕ jednostajnie na Ω. Istotnie,
|ϕν(x) − ϕ(x)| 6 Z b
bν
g0(t)dt −→ 0.
Wynika stąd, że ϕ jest ciągła. Z tych samych powodów, dla dowolnego |α| 6 k ciąg (Dαϕν)∞ν=1 jest zbieżny jednostajnie. Teraz wystarczy już tylko wykorzystać twierdzenie o różniczkowaniu ciągu wyraz
po wyrazie.
Obserwacja 3.3.5. Przedział [a, b) można zastąpić przedziałem (a, b] lub też przedziałem (a, b).
Przykład 3.3.6 (Funkcja Γ Eulera 35 ). Niech Γ (x) :=
Z ∞ 0
tx−1e−tdt, x > 0.
Wtedy
• Γ ∈ C∞(R>0),
• Γ (1) = 1, xΓ (x) = Γ (x + 1), x > 0.
W szczególności, Γ (n + 1) = n!, n ∈ Z+. Dowód . Niech f (x, t) := tx−1e−t, x, t > 0. Wtedy
Γ (x) = Z 1
0
f (x, t)dt + Z ∞
1
f (x, t)dt, x > 0.
Zauważmy, że ∂x∂kfk(x, t) = tx−1(ln t)ke−t. Wobec poprzedniego twierdzenia, dla dowodu, że Γ jest klasy C∞ wystarczy pokazać, że dla dowolnego przedziału [α, β] ⊂ (0, +∞) i dla dowolnego k ∈ Z+ istnieją funkcje ciągłe gk : (0, 1] −→ R, hk: [1, +∞) −→ R takie, że
• |tx−1(ln t)ke−t| 6 gk(t), x ∈ [α, β], t ∈ (0, 1],
• |tx−1(ln t)ke−t| 6 hk(t), x ∈ [α, β], t ∈ [1, ∞),
• całkiR1
0 gk(t)dt iR∞
1 hk(t)dt są zbieżne.
Przyjmujemy:
gk(t) := tα−1| ln t|ke−t, 0 < t 6 1; całkowalność gk wynika np. stąd, że dla dowolnego 0 < ε < α funkcja tα−1−εjest całkowalna 36
oraz tα−1| ln t|ke−t6 const(ε)tα−1−ε.
hk(t) := N !tβ−1+k−N > tβ−1(ln t)ke−t, t> 1, gdzie N > β + k; całkowalność hkjest oczywista 37 . Wzory występujące w tezie uzyskujemy przez proste całkowanie przez części.
34 Zauważmy, że nasze założenia gwarantują zbieżność wszystkich występujących w tezie całek niewłaściwych.
35 Leonhard Euler (1707–1783).
36 R1
0 tc−1dt =1c dla c > 0.
37 R∞
1 tc−1dt = −1c dla c < 0.
3. Całka Riemanna
Twierdzenie 3.3.7 (Twierdzenie o zmianie zmiennych w całce Riemanna). Niech Φ : U −→ V będzie dyfeomorfizmem klasy C1 zbiorów otwartych U, V ⊂ Rn i niech A ⊂⊂ U będzie zbiorem regularnym (ze względu na całkę Riemanna). Wtedy:
• zbiór Φ(A) jest regularny,
• dla dowolnej funkcji f ∈R(Φ(A)) funkcja (f ◦Φ)·|Φ0| jest całkowalna na Φ(A), gdzie |Φ0| oznacza moduł wyznacznika macierzy Jacobiego odwzorowania Φ 38 ,
• R
Φ(A)f =R
A(f ◦ Φ) · |Φ0|, f ∈R(Φ(A)).
Dowód zostanie podany później (Wniosek 4.12.10). Teraz tylko zauważmy, że regularność Φ(A) wy-nika z Propozycji 3.1.22(b) oraz że twierdzenie jest prawdziwe dla translacji.
Ćwiczenie 3.3.8. Wyznaczyć det Φ0, sprawdzić czy Φ|U jest dyfeomorfizmem oraz wyznaczyć V := Φ(U ) dla następujących transformacji i obszarów U (a, b, c > 0 oznaczają stałe):
• (współrzędne biegunowe) Φ : R2−→ R2, Φ(r, ϕ) := (ar cos ϕ, br sin ϕ),
U := R>0× (0, 2π);
• (współrzędne walcowe) Φ : R3−→ R3, Φ(r, ϕ, z) := (ar cos ϕ, br sin ϕ, cz), U := R>0× (0, 2π) × R;
• (współrzędne sferyczne) Φ : R3−→ R3, Φ(r, ϕ, θ) := (ar cos ϕ cos θ, br sin ϕ cos θ, cr sin θ), U := R>0× (0, 2π) × (−π2,π2);
• (współrzędne sferyczne w Rn) Φ : Rn−→ Rn,
Φ1(r, ω) :=r cos ω1cos ω2cos ω3. . . cos ωn−2cos ωn−1, Φ2(r, ω) :=r sin ω1cos ω2cos ω3. . . cos ωn−2cos ωn−1, Φ3(r, ω) :=r sin ω2cos ω3. . . cos ωn−2cos ωn−1,
. . .
Φn−1(r, ω) :=r sin ωn−2cos ωn−1, Φn(r, ω) :=r sin ωn−1,
U := R>0× (0, 2π) × (−π2,π2)n−2. Ćwiczenie 3.3.9. Niech a1, . . . , an> 0,
E :=
n
(x1, . . . , xn) ∈ Rn: x21
a21 + · · · +x2n a2n 6 1o
. Obliczyć |E|.
3.4. Całki krzywoliniowe. Wzór Greena
Uwaga 3.4.1. Przypomnijmy, że długością krzywej γ : [a, b] −→ X w przestrzeni metrycznej (X, %) nazywamy liczbę L%(γ) ∈ [0, +∞] daną wzorem:
L%(γ) := sup
π
{S%(γ, π)}, gdzie
S%(γ, π) :=
m
X
j=1
%(γ(tj−1), γ(tj)),
a supremum jest brane po wszystkich podziałach π = (t0, . . . , tm) przedziału [a, b] 39 . Krzywą γ nazywamy %–prostowalną, gdy L%(γ) < +∞. Wiemy, że:
38 Tzn. |Φ0|(x) := | det[∂Φ∂xj
k(x)]j,k=1,...,n|.
39 Zauważmy, że jeżeli π04 π, to S%(γ, π0) > S%(γ, π).
3.4. Całki krzywoliniowe. Wzór Greena
• L%(γ) nie zależy od parametryzacji γ (i dlatego możemy zawsze ograniczyć się do krzywych γ : [0, 1] −→ X) 40 ,
• L%( γ) = L%(γ), tzn. długość krzywej nie zależy od orientacji 41 ,
• L%(γ1⊕ γ2) = L%(γ1) + L%(γ2) 42 .
Propozycja 3.4.2. Dla każdego normalnego ciągu podziałów (πk)∞k=1 przedziału [0, 1] mamy:
S%(γ, πk) −→ L%(γ).
Dowód . Wystarczy pokazać, że dla dowolnej liczby ` < L%(γ) istnieje δ > 0 takie, że S%(γ, π) > ` o ile diam π 6 δ. Niech π0 = (t00, . . . , t0m
0) będzie podziałem przedziału [0, 1] takim, że S%(γ, π0) > `.
Dla 0 < δ < min{t0j − t0j−1 : j = 1, . . . , m0} weźmy dowolny podział π = (t0, . . . , tm) odcinka [0, 1]
o średnicy 6 δ. Punkty t0, . . . , tm dzielimy na m0–grup, zaliczając do i-tej grupy te punkty tj, dla których t0i−16 tj< t0i, i = 1, . . . , m0− 1, a do grupy m0 — pozostałą „końcówkę” punktów. Punkty i-tej grupy numerujemy kolejno tni, . . . , tni+1−1, gdzie n1 = 0, nm0+1 = m + 1. Wobec nierówności trójkąta mamy:
S%(γ, π) − S%(γ, π0) >
m0
X
i=2
%(γ(tni−1), γ(tni)) − %(γ(tni−1), γ(t0i−1)) − %(γ(t0i−1), γ(tni))
> −2(m0− 1)ωγ(δ),
gdzie ωγ(δ) oznacza moduł ciągłości odwzorowania γ zdefiniowany następująco ωγ(δ) := sup{%(γ(t0), γ(t00)) : |t0− t00| 6 δ}.
Teraz pozostaje już tylko wykorzystać jednostajną ciągłość odwzorowania γ 43 , z której wynika, że dla
dostatecznie małych δ mamy S%(γ, π) > `.
W dalszym ciągu będziemy zainteresowani krzywymi γ : [0, 1] −→ Rn (i sytuacją, gdy na Rn rozwa-żamy odległość euklidesową).
Powiemy, że krzywa γ : [0, 1] −→ Rn jest kawałkami klasy Ck (γ ∈ bCk([0, 1], Rn)) jeżeli istnieje podział (t0, . . . , tm) odcinka [0, 1] taki, że γ|[tj−1,tj] ∈ Ck([tj−1, tj]), 1 6 j 6 m 44 . Droga to krzywa kawałkami klasy C1.
Propozycja 3.4.3. Dowolna droga γ : [0, 1] −→ Rn jest prostowalna 45 oraz L(γ) =
Z 1 0
kγ0(t)kdt. 46
Dowód . Możemy założyć, że γ jest klasy C1. Dla podziału π = (t0, . . . , tm) niech ξ := (t0, . . . , tm−1).
Wtedy, korzystając z twierdzenia o przyrostach skończonych, dostajemy:
|S(γ, π) − M (kγ0k, π, ξ)| =
m
X
j=1
kγ(tj) − γ(tj−1)k −
m
X
j=1
kγ0(tj−1)k(tj− tj−1) 6
m
X
j=1
kγ(tj) − γ(tj−1) − γ0(tj−1)(tj− tj−1)k
6
m
X
j=1
max{kγ0(t) − γ0(tj−1)k : t ∈ [tj−1, tj]}
(tj− tj−1) 6 ωγ0(diam π). 47
40 Przypomnijmy, że zmiana parametryzacji krzywej γ : [a, b] −→ X polega na zastąpieniu jej przez krzywą γ ◦ σ : [c, d] −→ X, gdzie σ : [c, d] −→ [a, b] jest pewną bijekcją rosnącą.
41 Przypomnijmy, że dla krzywej γ : [0, 1] −→ X kładziemy: ( γ)(t) := γ(1 − t).
42 Przypomnijmy, że dla krzywych γj: [0, 1] −→ X, j = 1, 2, takich, że γ1(1) = γ2(0) definiujemy: (γ1⊕ γ2)(t) :=
γ1(2t) dla 0 6 t 612 i (γ1⊕ γ2)(t) := γ2(2t − 1) dla 12 6 t 6 1.
43 Co w języku modułów ciągłości oznacza, że dla dowolnego ε > 0 istnieje δ > 0 taka, że ωγ(δ) 6 ε.
44 W szczególności, dla dowolnego t ∈ [0, 1], można mówić o pochodnych jednostronnych γ−(j)(t), γ+(j)(t), j = 1, . . . , k.
45 Względem odległości euklidesowej.
46 Jeżeli % jest metryką euklidesową, to będziemy pomijać wskaźnik%. Zauważmy, że całka po prawej stronie istnieje.
3. Całka Riemanna
Pozostaje wykorzystać Propozycję 3.4.2, Własność 3.1.12 i skorzystać z jednostajnej ciągłości
odwzoro-wania γ0.
Przykład 3.4.4. (a) Niech γ : [a, b] −→ R3, γ(t) := (r cos t, r sin t, st) 48 . Wtedy L(γ) =p
r2+ s2(b − a).
(b) Niech γ : [a, b] −→ R2, γ(t) := (A cos t, B sin t). Wtedy L(γ) =
Z b a
pA2sin2t + B2cos2tdt = ?
Niech γ : [a, b] −→ Rn będzie krzywą i niech f : γ∗−→ R będzie dowolną funkcją ograniczoną 49 . Dla podziału π = (t0, . . . , tm) przedziału [a, b] i dla punktów pośrednich ξ = (ξ1, . . . , ξm) niech
Niech γ : [a, b] −→ Rn będzie krzywą i niech f : γ∗−→ R będzie dowolną funkcją ograniczoną 49 . Dla podziału π = (t0, . . . , tm) przedziału [a, b] i dla punktów pośrednich ξ = (ξ1, . . . , ξm) niech