Analiza matematyczna 2/Wykład 11: Twierdzenie Fubiniego. Twierdzenie o zmianie zmiennych: Różnice pomiędzy wersjami

Z Studia Informatyczne
Przejdź do nawigacjiPrzejdź do wyszukiwania
m (Zastępowanie tekstu - "<div class="thumb t(.*)"><div style="width:(.*);"> <flash>file=(.*)\.swf\|width=(.*)\|height=(.*)<\/flash> <div\.thumbcaption>(.*)<\/div> <\/div><\/div>" na "$4x$5px|thumb|$1|$6")
 
(Nie pokazano 13 wersji utworzonych przez 4 użytkowników)
Linia 61: Linia 61:
  
 
<center><math>\displaystyle \iint\limits_{[0,1]\times[0,1]}xy\ dxdy
 
<center><math>\displaystyle \iint\limits_{[0,1]\times[0,1]}xy\ dxdy
\ =\
+
=
 
\displaystyle\int\limits_0^1\left(\displaystyle\int\limits_0^1 xy dx\right) dy
 
\displaystyle\int\limits_0^1\left(\displaystyle\int\limits_0^1 xy dx\right) dy
\ =\
+
=
 
\displaystyle\int\limits_0^1\left(\displaystyle\int\limits_0^1 xy dy\right) dx.
 
\displaystyle\int\limits_0^1\left(\displaystyle\int\limits_0^1 xy dy\right) dx.
 
</math></center>
 
</math></center>
Linia 81: Linia 81:
  
 
<center><math>\displaystyle \displaystyle\int\limits_{A\times B}f(x,y)\ dxdy=\displaystyle\int\limits_A\left(\displaystyle\int\limits_B f(x,y) dy\right)dx
 
<center><math>\displaystyle \displaystyle\int\limits_{A\times B}f(x,y)\ dxdy=\displaystyle\int\limits_A\left(\displaystyle\int\limits_B f(x,y) dy\right)dx
\ =\
+
=
 
\displaystyle\int\limits_B\left(\displaystyle\int\limits_A f(x,y) dx\right)dy.
 
\displaystyle\int\limits_B\left(\displaystyle\int\limits_A f(x,y) dx\right)dy.
 
</math></center>
 
</math></center>
Linia 87: Linia 87:
 
}}
 
}}
  
<div class="thumb tleft"><div style="width:375px;">
+
[[File:AM2.M11.W.r01.svg|375x375px|thumb|left|Rysunek do twierdzenia Fubiniego]]
<flash>file=AM2.M11.W.R01.swf|width=375|height=375</flash>
 
<div.thumbcaption>Rysunek do twierdzenia Fubiniego</div>
 
</div></div>
 
 
{{uwaga|11.2.||
 
{{uwaga|11.2.||
  
Linia 101: Linia 98:
 
<center>
 
<center>
 
<math>\displaystyle \displaystyle\int\limits_{A\times B}f(x,y)\ dxdy
 
<math>\displaystyle \displaystyle\int\limits_{A\times B}f(x,y)\ dxdy
\ =\
+
=
 
\displaystyle\int\limits_A\left(\displaystyle\int\limits_B
 
\displaystyle\int\limits_A\left(\displaystyle\int\limits_B
 
f(x,y) dy\right) dx=\displaystyle\int\limits_B\left(\displaystyle\int\limits_A f(x,y) dx\right)dy.
 
f(x,y) dy\right) dx=\displaystyle\int\limits_B\left(\displaystyle\int\limits_A f(x,y) dx\right)dy.
Linia 139: Linia 136:
  
 
<center><math>\displaystyle \displaystyle\int\limits_{[a,b]\times [c,d]}f(x,y)\ dxdy
 
<center><math>\displaystyle \displaystyle\int\limits_{[a,b]\times [c,d]}f(x,y)\ dxdy
\ =\
+
=
 
\displaystyle\int\limits_a^b\left(\displaystyle\int\limits_c^d f(x,y)dy\right) dx
 
\displaystyle\int\limits_a^b\left(\displaystyle\int\limits_c^d f(x,y)dy\right) dx
\ =\
+
=
 
\displaystyle\int\limits_c^d\left(\displaystyle\int\limits_a^b f(x,y) dx\right)dy.
 
\displaystyle\int\limits_c^d\left(\displaystyle\int\limits_a^b f(x,y) dx\right)dy.
 
</math></center>
 
</math></center>
Linia 154: Linia 151:
 
<center>
 
<center>
 
<math>\displaystyle \displaystyle\int\limits_{[a,b]\times [c,d]}f(x,y)\ dxdy
 
<math>\displaystyle \displaystyle\int\limits_{[a,b]\times [c,d]}f(x,y)\ dxdy
\ =\
+
=
 
\displaystyle\int\limits_a^b\left(\displaystyle\int\limits_c^d  f(x,y) dy\right) dx.
 
\displaystyle\int\limits_a^b\left(\displaystyle\int\limits_c^d  f(x,y) dy\right) dx.
 
</math>
 
</math>
Linia 166: Linia 163:
 
<center>
 
<center>
 
<math>\displaystyle d_2((x_1,y_1),(x_2,y_2))
 
<math>\displaystyle d_2((x_1,y_1),(x_2,y_2))
\ =\
+
=
 
\sqrt{(x_2-x_1)^2+(y_2-y_1)^2}.
 
\sqrt{(x_2-x_1)^2+(y_2-y_1)^2}.
 
</math>
 
</math>
Linia 192: Linia 189:
  
 
<center><math>\displaystyle g(x)
 
<center><math>\displaystyle g(x)
\ :=\
+
\ :=
 
\displaystyle\int\limits_c^d  f(x,y) dy
 
\displaystyle\int\limits_c^d  f(x,y) dy
 
</math></center>
 
</math></center>
Linia 262: Linia 259:
 
odcinków, czyli:
 
odcinków, czyli:
  
<center><math>\displaystyle \aligned
+
<center><math>\displaystyle \begin{align}
 
a_i^n
 
a_i^n
&= a+\frac{i}{n}(b-a), \quad b_i^n=a_i^n+\frac{1}{n}(b-a),\
+
&= a+\frac{i}{n}(b-a), \quad b_i^n=a_i^n+\frac{1}{n}(b-a),
 
i=0,1,\ldots,n-1,\\
 
i=0,1,\ldots,n-1,\\
 
c_j^n
 
c_j^n
 
&=
 
&=
c+\frac{j}{n}(d-c), \quad d_j^n=c_j^n+\frac{1}{n}(d-c),\
+
c+\frac{j}{n}(d-c), \quad d_j^n=c_j^n+\frac{1}{n}(d-c),
 
j=0,1,\ldots,n-1,
 
j=0,1,\ldots,n-1,
\endaligned</math></center>
+
\end{align}</math></center>
  
 
a następnie biorąc iloczyn kartezjański tych
 
a następnie biorąc iloczyn kartezjański tych
Linia 276: Linia 273:
  
 
<center><math>\displaystyle K_{ij}^n
 
<center><math>\displaystyle K_{ij}^n
\ :=\
+
\ :=
 
[a_i^n,b_i^n]\times[c_j^n,d_j^n],\ \ i,j=0,\ldots,n-1.
 
[a_i^n,b_i^n]\times[c_j^n,d_j^n],\ \ i,j=0,\ldots,n-1.
 
</math></center>
 
</math></center>
Linia 291: Linia 288:
  
 
<center><math>\displaystyle S_n
 
<center><math>\displaystyle S_n
\ :=\
+
\ :=
 
\sum_{i,j=0}^{n-1}f(p_{ij}^n)(b_i^n-a_i^n)(d_j^n-c_j^n).
 
\sum_{i,j=0}^{n-1}f(p_{ij}^n)(b_i^n-a_i^n)(d_j^n-c_j^n).
 
</math></center>
 
</math></center>
Linia 298: Linia 295:
  
 
<center><math>\displaystyle \lim_{n\to\infty}S_n
 
<center><math>\displaystyle \lim_{n\to\infty}S_n
\ =\
+
=
 
\displaystyle\int\limits_{[a,b]\times [c,d]}f(x,y)dxdy.
 
\displaystyle\int\limits_{[a,b]\times [c,d]}f(x,y)dxdy.
 
</math></center>
 
</math></center>
Linia 312: Linia 309:
  
 
<center><math>\displaystyle \bigg|\displaystyle\int\limits_a^b\left(\displaystyle\int\limits_c^d  f(x,y)dy\right)dx-S_n\bigg|
 
<center><math>\displaystyle \bigg|\displaystyle\int\limits_a^b\left(\displaystyle\int\limits_c^d  f(x,y)dy\right)dx-S_n\bigg|
\ <\
+
<
 
\varepsilon.
 
\varepsilon.
 
</math></center>
 
</math></center>
Linia 355: Linia 352:
  
 
<center><math>\displaystyle \bigg|\displaystyle\int\limits_c^df(x,y)dy-\sum_{j=0}^{n-1}f(p_{ij}^n)(d_j^n-c_j^n)\bigg|
 
<center><math>\displaystyle \bigg|\displaystyle\int\limits_c^df(x,y)dy-\sum_{j=0}^{n-1}f(p_{ij}^n)(d_j^n-c_j^n)\bigg|
\ <\
+
<
 
\varepsilon'(d-c).
 
\varepsilon'(d-c).
 
</math></center>
 
</math></center>
Linia 408: Linia 405:
 
zastosować twierdzenie Fubiniego. Otrzymamy
 
zastosować twierdzenie Fubiniego. Otrzymamy
  
<center><math>\displaystyle \aligned
+
<center><math>\displaystyle \begin{align}
 
\iint\limits_K xy-y^2 dxdy
 
\iint\limits_K xy-y^2 dxdy
 
&=
 
&=
 
\displaystyle\int\limits_1^2dx\displaystyle\int\limits_3^4(xy-y^2)dy
 
\displaystyle\int\limits_1^2dx\displaystyle\int\limits_3^4(xy-y^2)dy
\ =\
+
=
 
\displaystyle\int\limits_1^2\left((x\frac{y^2}{2}-\frac{y^3}{3})\bigg|_3^4\right)dx\\
 
\displaystyle\int\limits_1^2\left((x\frac{y^2}{2}-\frac{y^3}{3})\bigg|_3^4\right)dx\\
 
&=\displaystyle\int\limits_1^2\left(\frac{7x}{2}-\frac{37}{3}\right)dx
 
&=\displaystyle\int\limits_1^2\left(\frac{7x}{2}-\frac{37}{3}\right)dx
\ =\
+
=
 
\left(\frac{7x^2}{4}-\frac{37x}{3}\right)\bigg|_1^2=-\frac{85}{12}.
 
\left(\frac{7x^2}{4}-\frac{37x}{3}\right)\bigg|_1^2=-\frac{85}{12}.
\endaligned</math></center>}}
+
\end{align}</math></center>}}
 
<br>
 
<br>
 
<center>
 
<center>
Linia 425: Linia 422:
 
     <param name="coloring" value="maple">
 
     <param name="coloring" value="maple">
 
     <param name="model" value="images/b/bd/Am2.m11.w.r02.mgs.zip">
 
     <param name="model" value="images/b/bd/Am2.m11.w.r02.mgs.zip">
     <param name="scale" value="1.0 1.0 1.0">
+
     <param name="scale" value="1.0 1.0 0.2">
 
<param name="shading" value="0.2">
 
<param name="shading" value="0.2">
 
</applet>
 
</applet>
Linia 445: Linia 442:
  
 
<center><math>\displaystyle D
 
<center><math>\displaystyle D
\ :=\
+
\ :=
 
\{(x,y)\in\mathbb{R}^2 : a\leq x\leq b, h_1(x)\leq y\leq h_2(x)\}
 
\{(x,y)\in\mathbb{R}^2 : a\leq x\leq b, h_1(x)\leq y\leq h_2(x)\}
 
</math></center>
 
</math></center>
Linia 461: Linia 458:
  
 
<center><math>\displaystyle D
 
<center><math>\displaystyle D
\ =\
+
=
 
\{(x,y,z)\in\mathbb{R}^3 : (x,y)\in A,                g_1(x,y)
 
\{(x,y,z)\in\mathbb{R}^3 : (x,y)\in A,                g_1(x,y)
\ \leq\
+
\ \leq
 
z\leq g_2(x,y)\}.
 
z\leq g_2(x,y)\}.
 
</math></center>
 
</math></center>
Linia 481: Linia 478:
  
 
{| border="0" align="center" cellspacing="10"
 
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:375px;">
+
|[[File:AM2.M11.W.R03.svg|375x375px|thumb|center|Zbiór normalny względem osi <math>Ox</math>]]
<flash>file=AM2.M11.W.R03.swf|width=375|height=375</flash>
+
|[[File:AM2.M11.W.R04.svg|375x375px|thumb|center|Zbiór normalny względem osi <math>Oy</math>]]
<div.thumbcaption>Zbiór normalny względem osi <math>Ox</math></div>
 
</div></div>
 
|<div class="thumb"><div style="width:375px;">
 
<flash>file=AM2.M11.W.R04.swf|width=375|height=375</flash>
 
<div.thumbcaption>Zbiór normalny względem osi <math>Oy</math></div>
 
</div></div>
 
 
|}
 
|}
<div class="thumb tright"><div style="width:375px;">
+
[[File:AM2.M11.W.R05.svg|375x375px|thumb|right|Zbiór regularny, który nie jest zbiorem normalnym]]
<flash>file=AM2.M11.W.R05.swf|width=375|height=375</flash>
 
<div.thumbcaption>Zbiór regularny, który nie jest zbiorem normalnym</div>
 
</div></div>
 
  
 
Definicje normalności i regularności można oczywiście
 
Definicje normalności i regularności można oczywiście
Linia 506: Linia 494:
 
<br><center>
 
<br><center>
 
<math>\displaystyle A
 
<math>\displaystyle A
\ :=\
+
\ :=
 
\{(x,y)\in\mathbb{R}^2 : a\leq x\leq b, h_1(x)\leq y\leq h_2(x)\},
 
\{(x,y)\in\mathbb{R}^2 : a\leq x\leq b, h_1(x)\leq y\leq h_2(x)\},
 
</math>
 
</math>
Linia 528: Linia 516:
 
<br><center>
 
<br><center>
 
<math>\displaystyle \iint\limits_Af(x,y)dxdy
 
<math>\displaystyle \iint\limits_Af(x,y)dxdy
\ =\
+
=
 
\displaystyle\int\limits_a^bdx\displaystyle\int\limits_{h_1(x)}^{h_2(x)}f(x,y)dy.
 
\displaystyle\int\limits_a^bdx\displaystyle\int\limits_{h_1(x)}^{h_2(x)}f(x,y)dy.
 
</math>
 
</math>
Linia 537: Linia 525:
  
 
<center><math>\displaystyle \iiint\limits_D f(x,y,z)dxdydz
 
<center><math>\displaystyle \iiint\limits_D f(x,y,z)dxdydz
\ =\
+
=
 
\displaystyle\int\limits_a^bdx\displaystyle\int\limits_{h_1(x)}^{h_2(x)}dy\displaystyle\int\limits_{g_1(x,y)}^{g_2(x,y)}f(x,y,z)dz.
 
\displaystyle\int\limits_a^bdx\displaystyle\int\limits_{h_1(x)}^{h_2(x)}dy\displaystyle\int\limits_{g_1(x,y)}^{g_2(x,y)}f(x,y,z)dz.
 
</math></center>
 
</math></center>
Linia 549: Linia 537:
  
 
Możemy teraz policzyć następującą całkę.
 
Możemy teraz policzyć następującą całkę.
<div class="thumb tright"><div style="width:375px;">
+
[[File:AM2.M11.W.R06.svg|375x375px|thumb|right|Trójkąt <math>T</math>]]
<flash>file=AM2.M11.W.R06.swf|width=375|height=375</flash>
 
<div.thumbcaption>Trójkąt <math>T</math></div>
 
</div></div>
 
 
{{przyklad|11.8.||
 
{{przyklad|11.8.||
  
Linia 590: Linia 575:
 
&=&
 
&=&
 
\bigg(\frac{1}{10}x^5-\frac{1}{6}x^3\bigg)\bigg|_1^2+\bigg(\frac{-3}{10}x^5+\frac{3}{2}x^4-\frac{3}{2}x^3\bigg)\bigg|_2^3
 
\bigg(\frac{1}{10}x^5-\frac{1}{6}x^3\bigg)\bigg|_1^2+\bigg(\frac{-3}{10}x^5+\frac{3}{2}x^4-\frac{3}{2}x^3\bigg)\bigg|_2^3
\ =\
+
=
 
\frac{57}{10}+\frac{29}{15}
 
\frac{57}{10}+\frac{29}{15}
\ =\
+
=
 
\frac{229}{30}.
 
\frac{229}{30}.
 
\end{array}</math>
 
\end{array}</math>
Linia 599: Linia 584:
 
<center>
 
<center>
 
<div class="thumb"><div style="width:450px;">
 
<div class="thumb"><div style="width:450px;">
<applet code="applet.class" archive="images/a/a0/Javaview.jar,images/e/e1/Applet.jar" width="450" height="400">
+
<applet code="JavaviewModApplet.class" archive="images/a/a0/Javaview.jar,images/b/be/JavaviewModApplet.jar" width="450" height="400">
 
     <param name="colors" value="-1:#7e0b9f -0.80:#ce3bf8 -0.60:#2c4ae5 -0.40:#2c85e5 -0.20:#2ecca9 -0.05:#2ecc5b 0.05:#2ecc5b 0.20:#97cc2e 0.40:#edff27 0.60:#ffba27 0.80:#ff6e27 1:#d42525">
 
     <param name="colors" value="-1:#7e0b9f -0.80:#ce3bf8 -0.60:#2c4ae5 -0.40:#2c85e5 -0.20:#2ecca9 -0.05:#2ecc5b 0.05:#2ecc5b 0.20:#97cc2e 0.40:#edff27 0.60:#ffba27 0.80:#ff6e27 1:#d42525">
 
     <param name="coloring" value="maple">
 
     <param name="coloring" value="maple">
 
     <param name="model" value="images/d/d6/Am2.m11.w.r07.mgs.zip">
 
     <param name="model" value="images/d/d6/Am2.m11.w.r07.mgs.zip">
 
     <param name="scale" value="1.0 1.0 0.1">
 
     <param name="scale" value="1.0 1.0 0.1">
<param name="shading" value="0.2">
+
    <param name="shading" value="0.2">
 +
    <param name="animation" value="stop">
 
</applet>
 
</applet>
 
<div.thumbcaption>Wykres funkcji <math>f(x,y)=x^2y</math> nad <math>T</math></div>
 
<div.thumbcaption>Wykres funkcji <math>f(x,y)=x^2y</math> nad <math>T</math></div>
Linia 629: Linia 615:
  
 
<center> Jac <math>\displaystyle  _x\varphi
 
<center> Jac <math>\displaystyle  _x\varphi
\ =\
+
=
 
\left[
 
\left[
 
\begin{array} {ccc}\displaystyle  
 
\begin{array} {ccc}\displaystyle  
Linia 655: Linia 641:
  
 
<center><math>\displaystyle \displaystyle\int\limits_Df(y)dy_1\ldots dy_n
 
<center><math>\displaystyle \displaystyle\int\limits_Df(y)dy_1\ldots dy_n
\ =\
+
=
 
\displaystyle\int\limits_Bf(\varphi(x))|\det  </math> Jac <math>\displaystyle  _x\varphi|dx_1\ldots dx_n.
 
\displaystyle\int\limits_Bf(\varphi(x))|\det  </math> Jac <math>\displaystyle  _x\varphi|dx_1\ldots dx_n.
 
</math></center>
 
</math></center>
Linia 667: Linia 653:
  
 
<center><math>\displaystyle \displaystyle\int\limits_Df(y)dy
 
<center><math>\displaystyle \displaystyle\int\limits_Df(y)dy
\ =\
+
=
 
\displaystyle\int\limits_Bf(\varphi(x))\varphi'(x)dx.
 
\displaystyle\int\limits_Bf(\varphi(x))\varphi'(x)dx.
 
</math></center>
 
</math></center>
Linia 684: Linia 670:
  
 
==Zmiana zmiennych na dwuwymiarowe współrzędne biegunowe==
 
==Zmiana zmiennych na dwuwymiarowe współrzędne biegunowe==
<div class="thumb tright"><div style="width:375px;">
+
[[File:AM2.M11.W.R08.svg|375x375px|thumb|right|Współrzędne biegunowe]]
<flash>file=AM2.M11.W.R08.swf|width=375|height=375</flash>
 
<div.thumbcaption>Współrzędne biegunowe</div>
 
</div></div>
 
  
 
Niech zbiorem <math>\displaystyle D</math> będzie <math>\displaystyle  \mathbb{R}^2\setminus \{(x,0) :x\geq
 
Niech zbiorem <math>\displaystyle D</math> będzie <math>\displaystyle  \mathbb{R}^2\setminus \{(x,0) :x\geq
Linia 702: Linia 685:
  
 
<center>
 
<center>
<math>\displaystyle x\ =\ r\cos\alpha,\qquad
+
<math>\displaystyle x= r\cos\alpha,\qquad
y\ =\ r\sin\alpha.
+
y= r\sin\alpha.
 
</math>
 
</math>
 
</center>
 
</center>
Linia 751: Linia 734:
  
 
<center><math>\displaystyle \iint\limits_Dx^2+y^2  dxdy=\iint\limits_B(r^2)rdrd\alpha
 
<center><math>\displaystyle \iint\limits_Dx^2+y^2  dxdy=\iint\limits_B(r^2)rdrd\alpha
\ =\
+
=
 
\displaystyle\int\limits_0^{2\pi}d\alpha\displaystyle\int\limits_0^R r^3 dr
 
\displaystyle\int\limits_0^{2\pi}d\alpha\displaystyle\int\limits_0^R r^3 dr
\ =\
+
=
 
\displaystyle\int\limits_0^{2 \pi}
 
\displaystyle\int\limits_0^{2 \pi}
\frac{R^4}{4} d\alpha=2\pi R^4,
+
\frac{R^4}{4} d\alpha=\frac{\pi}{2} R^4,
 
</math></center>
 
</math></center>
  
Linia 772: Linia 755:
 
    
 
    
 
Stosujemy taką samą zmianę zmiennych. Tym razem <math>\displaystyle r</math> zmienia się
 
Stosujemy taką samą zmianę zmiennych. Tym razem <math>\displaystyle r</math> zmienia się
także od <math>\displaystyle 0</math> do <math>\displaystyle R,</math> natomiast <math>\displaystyle \displaystyle\alpha\</math> zmienia się od <math>\displaystyle\frac{\pi}{2}</math> do <math>\displaystyle\pi</math>.
+
także od <math>\displaystyle 0</math> do <math>\displaystyle R,</math> natomiast <math>\displaystyle \displaystyle\alpha</math> zmienia się od <math>\displaystyle\frac{\pi}{2}</math> do <math>\displaystyle\pi</math>.
 
Tak więc <math>\displaystyle B=[0,R]\times \bigg[\frac{\pi}{2}, \pi\bigg]:</math>
 
Tak więc <math>\displaystyle B=[0,R]\times \bigg[\frac{\pi}{2}, \pi\bigg]:</math>
  
Linia 782: Linia 765:
 
drd\alpha=\displaystyle\int\limits_{\frac{\pi}{2}}^{\pi}d\alpha\displaystyle\int\limits_0^R r^2\cos\alpha
 
drd\alpha=\displaystyle\int\limits_{\frac{\pi}{2}}^{\pi}d\alpha\displaystyle\int\limits_0^R r^2\cos\alpha
 
dr\\
 
dr\\
\ &=&\
+
\ &=&
 
\displaystyle\int\limits_{\frac{\pi}{2}}^{\pi}\frac{R^3}{3}\cos\alpha  d\alpha
 
\displaystyle\int\limits_{\frac{\pi}{2}}^{\pi}\frac{R^3}{3}\cos\alpha  d\alpha
 
=\displaystyle  
 
=\displaystyle  
 
\frac{R^3}{3}(-\sin\alpha)\bigg|_{\frac{\pi}{2}}^{\pi}
 
\frac{R^3}{3}(-\sin\alpha)\bigg|_{\frac{\pi}{2}}^{\pi}
\ =\
+
=
 
\frac{R^3}{3}.\end{array}
 
\frac{R^3}{3}.\end{array}
 
</math></center>}}
 
</math></center>}}
Linia 792: Linia 775:
  
 
{| border="0" align="center" cellspacing="10"
 
{| border="0" align="center" cellspacing="10"
|<div class="thumb"><div style="width:375px;">
+
|[[File:AM2.M11.W.R09.svg|375x375px|thumb|center|Zbiór <math>D</math>]]
<flash>file=AM2.M11.W.R09.swf|width=375|height=375</flash>
+
|[[File:AM2.M11.W.R10.svg|375x375px|thumb|center|Zbiór <math>B</math>]]
<div.thumbcaption>Zbiór <math>D</math></div>
 
</div></div>
 
|<div class="thumb"><div style="width:375px;">
 
<flash>file=AM2.M11.W.R10.swf|width=375|height=375</flash>
 
<div.thumbcaption>Zbiór <math>B</math></div>
 
</div></div>
 
 
|}
 
|}
 +
 
==Zmiana zmiennych na trójwymiarowe współrzędne sferyczne==
 
==Zmiana zmiennych na trójwymiarowe współrzędne sferyczne==
<div class="thumb tright"><div style="width:375px;">
+
[[File:AM2.M11.W.R11.svg|375x375px|thumb|right|Współrzędne sferyczne]]
<flash>file=AM2.M11.W.R11.swf|width=375|height=375</flash>
 
<div.thumbcaption>Współrzędne sferyczne</div>
 
</div></div>
 
 
Podobnie do współrzędnych biegunowych w <math>\displaystyle \displaystyle\mathbb{R}^2</math>  
 
Podobnie do współrzędnych biegunowych w <math>\displaystyle \displaystyle\mathbb{R}^2</math>  
 
definiujemy współrzędne sferyczne w<math>\displaystyle \displaystyle\mathbb{R}^3.</math> Mamy:
 
definiujemy współrzędne sferyczne w<math>\displaystyle \displaystyle\mathbb{R}^3.</math> Mamy:
Linia 853: Linia 828:
  
 
<center>
 
<center>
<math>\displaystyle \aligned
+
<math>\displaystyle \begin{align}
 
\iiint\limits_D z^2  dxdydz&=\iiint\limits_{B} r^3 \sin\beta \cos\beta
 
\iiint\limits_D z^2  dxdydz&=\iiint\limits_{B} r^3 \sin\beta \cos\beta
 
d\alpha d\beta dr = \displaystyle\int\limits_0^{2\pi}d\alpha \displaystyle\int\limits_0^{\frac{\pi}{2}}
 
d\alpha d\beta dr = \displaystyle\int\limits_0^{2\pi}d\alpha \displaystyle\int\limits_0^{\frac{\pi}{2}}
Linia 859: Linia 834:
 
&= \frac{R^4}{4} \displaystyle\int\limits_0^{2\pi}d\alpha \displaystyle\int\limits_0^{\frac{\pi}{2}}
 
&= \frac{R^4}{4} \displaystyle\int\limits_0^{2\pi}d\alpha \displaystyle\int\limits_0^{\frac{\pi}{2}}
 
\sin\beta \cos\beta  d\beta
 
\sin\beta \cos\beta  d\beta
\ =\
+
=
 
\frac{R^4}{4} \displaystyle\int\limits_0^{2\pi}\left(\frac{1}{2}\sin^2\beta\bigg|_0^{\frac{\pi}{2}}\right)
 
\frac{R^4}{4} \displaystyle\int\limits_0^{2\pi}\left(\frac{1}{2}\sin^2\beta\bigg|_0^{\frac{\pi}{2}}\right)
 
d\alpha
 
d\alpha
\ =\
+
=
 
\frac{R^4}{4}\pi.
 
\frac{R^4}{4}\pi.
\endaligned</math>
+
\end{align}</math>
 
</center>}}
 
</center>}}
  
 
==Zmiana zmiennych na trójwymiarowe współrzędne walcowe==
 
==Zmiana zmiennych na trójwymiarowe współrzędne walcowe==
<div class="thumb tright"><div style="width:375px;">
+
[[File:AM2.M11.W.R12.svg|375x375px|thumb|right|AM2.M11.W.R12]]
<flash>file=AM2.M11.W.R12.swf|width=375|height=375</flash>
 
<div.thumbcaption>AM2.M11.W.R12</div>
 
</div></div>
 
  
 
Ta zmiana zmiennych jest w zasadzie zmianą na współrzędne
 
Ta zmiana zmiennych jest w zasadzie zmianą na współrzędne
Linia 909: Linia 881:
  
 
<center>
 
<center>
<math>\displaystyle \aligned
+
<math>\displaystyle \begin{align}
 
\iiint\limits_D z  dxdydz
 
\iiint\limits_D z  dxdydz
 
&=
 
&=
 
\iiint\limits_B rz d\alpha dr dz
 
\iiint\limits_B rz d\alpha dr dz
\ =\
+
=
 
\displaystyle\int\limits_0^{2\pi}d\alpha\displaystyle\int\limits_0^R dr\displaystyle\int\limits_0^H rz dz\\
 
\displaystyle\int\limits_0^{2\pi}d\alpha\displaystyle\int\limits_0^R dr\displaystyle\int\limits_0^H rz dz\\
\ &=\
+
\ &=
 
\frac{H^2}{2}\displaystyle\int\limits_0^{2\pi}d\alpha\displaystyle\int\limits_0^R r dr
 
\frac{H^2}{2}\displaystyle\int\limits_0^{2\pi}d\alpha\displaystyle\int\limits_0^R r dr
 
=
 
=
 
\frac{H^2}{2}\frac{R^2}{2}\displaystyle\int\limits_0^{2\pi}d\alpha
 
\frac{H^2}{2}\frac{R^2}{2}\displaystyle\int\limits_0^{2\pi}d\alpha
\ =\
+
=
 
\pi\frac{H^2R^2}{2}.
 
\pi\frac{H^2R^2}{2}.
\endaligned</math>
+
\end{align}</math>
 
</center>}}
 
</center>}}
  
 
Ciekawsze przykłady policzymy na ćwiczeniach.
 
Ciekawsze przykłady policzymy na ćwiczeniach.

Aktualna wersja na dzień 11:09, 3 paź 2021

Twierdzenie Fubiniego. Twierdzenie o zmianie zmiennych

W tym wykładzie prezentujemy twierdzenie Fubiniego (z dowodem tylko dla kostki w ) oraz twierdzenie o zmianie zmiennych w całce. Podajemy przykłady zmiany zmiennych w na współrzędne biegunowe oraz w na współrzędne walcowe i sferyczne.

Twierdzenie Fubiniego

Ten wykład poświęcony jest dwóm najważniejszym twierdzeniom dotyczącym całek wielokrotnych. Twierdzenie Fubiniego pozwala liczyć całki wielokrotne (podwójne, potrójne, itd.) po odpowiednich obszarach za pomocą kolejnego liczenia pewnych całek pojedynczych w odpowiednich granicach. Drugim z twierdzeń jest twierdzenie o zmianie zmiennych w całce, odpowiednik twierdzenia o całkowaniu przez podstawienie dla całki jednej zmiennej, także bardzo ważne dla obliczania całek.

Na ćwiczeniach do poprzedniego wykładu policzyliśmy z definicji gdzie

Policzmy teraz traktując jako stałą. Dostaniemy oczywiście

Następnie policzmy czyli całkę "z tego" co otrzymaliśmy wyżej. Dostaniemy

Policzyliśmy zatem

Jeśli policzymy "w drugą stronę", czyli najpierw całkę względem a potem względem to dostaniemy

następnie

zatem także

Otrzymaliśmy zatem następujące równości:

W takim razie możemy zapytać: czy może takie równości zachodzą zawsze? Okazuje się, że (przy rozsądnych założeniach) faktycznie tak jest - mówi o tym Twierdzenie Fubiniego.

Twierdzenie 11.1. [Twierdzenie Fubiniego]

Niech będzie kostką w a kostką w Zmienne w oznaczmy przez a w przez Weźmy funkcję Załóżmy, że dla każdego ustalonego funkcja jest całkowalna w sensie Riemanna na oraz że dla każdego ustalonego funkcja jest całkowalna w sensie Riemanna na Wtedy

Rysunek do twierdzenia Fubiniego
Uwaga 11.2.

(1) W szczególności, gdy funkcja jest ciągła na to obie funkcje i są całkowalne i zachodzą powyższe równości, czyli

(2) Nietrudno zauważyć, że w twierdzeniu Fubiniego zamiast kostek i możemy wziąć dowolne zbiory J-mierzalne - bo i tak całkowanie po dowolnych zbiorach J-mierzalnych sprowadziliśmy do całkowania po kostkach (patrz poprzedni wykład).
(3) Całki i nazywamy całkami iterowanymi.

Dowód 11.2.

Dowód twierdzenia Fubiniego przedstawimy tylko dla przypadku, gdy i są kostkami w (czyli jest kostką (prostokątem) w W tym przypadku twierdzenie i dowód łatwo zilustrować rysunkiem (patrz obok). Idea dowodu dla kostek wyżej wymiarowych jest dokładnie taka sama. Dla dodatkowego uproszczenia dowodu założymy, że funkcja jest ciągła. A zatem wypiszmy:

End of proof.gif

Twierdzenie 11.3. [Twierdzenie Fubiniego dla funkcji ciągłej na prostokącie]

Niech będzie kostką w Niech będzie funkcją ciągłą. Wówczas istnieją całki iterowane i oraz zachodzą równości

Euklides (365-300 p.n.e.)
Zobacz biografię

Dowód 11.3. [nadobowiązkowy]

Wykażemy istnienie całki i równość

Istnienia drugiej z całek iterowanych i drugiej równości dowodzi się analogicznie. Niech oznacza metrykę euklidesową w czyli

Krok I. Istnienie całki
I.1. Zauważmy, że dla dowolnego istnieje takie, że

dla z kostki Faktycznie, skoro funkcja jest ciągła, a zbiór jest zwarty, to funkcja jest jednostajnie ciągła (patrz Analiza matematyczna 1 uwaga 2.39). To dokładnie oznacza, że spełniona jest powyższa implikacja.
I.2. Wykażemy, że funkcja

jest funkcją ciągłą.

Ponieważ jest funkcją ciągłą, istnieje dla dowolnego Aby wykazać, że jest funkcją ciągłą, weźmy dowolne Szukamy takiego, że spełnione jest wynikanie:

Weźmy teraz Do tego dobierzmy tak jak w punkcie I.1. Mamy zatem w szczególności:

czyli, podstawiając do wzoru na otrzymujemy

Całkując te nierówności stronami (korzystamy z monotoniczności całki Riemanna dla funkcji jednej zmiennej), otrzymujemy:

czyli

zatem

przy co dowodzi ciągłości funkcji
I.3. Zauważmy, że skoro jest funkcją ciągłą na to istnieje a to dowodzi istnienia całki

Krok II. Równość
II.1. Z części I dowodu i z założeń twierdzenia wynika, że całki po obu stronach równości istnieją. Wystarczy zatem znaleźć granicę sum całkowych dla pewnego normalnego ciągu podziałów.
II.2. Zdefiniujmy normalny ciąg podziałów , dzieląc każdy z odcinków i na równych odcinków, czyli:

a następnie biorąc iloczyn kartezjański tych odcinków:

Kostkami podziału są więc kostki Objętość takiej kostki to
II.3. Weźmy teraz dla każdego podziału ciąg punktów pośrednich, czyli

Utwórzmy sumę całkową:

Skoro istnieje całka podwójna, to

Wystarczy zatem pokazać, że granicą ciągu jest też
II.4. Pokażemy, że

Musimy zatem pokazać, że dla ustalonego istnieje takie, że dla mamy

Ustalmy zatem Weźmy Do tego dobierzmy tak jak w punkcie I.1. dowodu. Dobierzmy takie, by W takim razie, jeśli dla mamy to a zatem (z I.1.),

czyli

Całkując te nierówności względem po przedziale dostaniemy (dla ustalonego ):

czyli

Weźmy teraz sumę powyższych nierówności dla (i dla ). Dostaniemy:

Tak więc

Całkując tę nierówność po przedziałach a następnie sumując wszystkie całki dla dostaniemy

a zatem, po zsumowaniu

co należało dowieść.

End of proof.gif
Uwaga 11.4. [Zapis całek iterowanych]

Całki iterowane, na przykład będziemy, w celu uniknięcia pisania dużej ilości nawiasów, zapisywali tak:

podobnie, zamiast napiszemy

Przykład 11.5.

Policzyć całkę

gdzie

Nasza funkcja jest ciągła na prostokącie, zatem możemy zastosować twierdzenie Fubiniego. Otrzymamy


<applet code="JavaviewModApplet.class" archive="images/a/a0/Javaview.jar,images/b/be/JavaviewModApplet.jar" width="450" height="400">

   <param name="colors" value="-1:#7e0b9f -0.80:#ce3bf8 -0.60:#2c4ae5 -0.40:#2c85e5 -0.20:#2ecca9 -0.05:#2ecc5b 0.05:#2ecc5b 0.20:#97cc2e 0.40:#edff27 0.60:#ffba27 0.80:#ff6e27 1:#d42525">
   <param name="coloring" value="maple">
   <param name="model" value="images/b/bd/Am2.m11.w.r02.mgs.zip">
   <param name="scale" value="1.0 1.0 0.2">

<param name="shading" value="0.2"> </applet>

<div.thumbcaption>Am2.m11.w.r02


Najczęściej spotykanymi obszarami, po których będziemy chcieli całkować, nie są jednak kostki, tylko tak zwane zbiory normalne. Zdefiniujmy:

Definicja 11.6.

(1) Niech będzie odcinkiem w niech i będą funkcjami ciągłymi na takimi, że Wtedy zbiór

nazywamy zbiorem normalnym względem osi
(2) Analogicznie definiujemy zbiór normalny względem osi
(3) Zbiór zawarty w jest normalny względem współrzędnej jeśli istnieje pewien zbiór normalny zawarty w płaszczyźnie oraz istnieją dwie funkcje takie, że oraz

(4) Analogicznie definiujemy zbiór normalny względem pozostałych współrzędnych.
(5) Zbiorem normalnym będziemy nazywać zbiór normalny względem jakiejś współrzędnej. Zbiorem regularnym będziemy nazywać zbiór, który można podzielić na sumę zbiorów regularnych o rozłącznych wnętrzach.


Zbiór normalny względem osi
Zbiór normalny względem osi
Zbiór regularny, który nie jest zbiorem normalnym

Definicje normalności i regularności można oczywiście uogólnić na więcej wymiarów, ale nie będziemy tego robić.

Jak już wspomnieliśmy, w praktyce najczęściej będziemy chcieli całkować funkcje po zbiorach normalnych. Wypiszmy więc, jak w przypadku takich zbiorów wygląda twierdzenie Fubiniego.

Niech zatem będzie zbiorem normalnym w zadanym jako



gdzie są jak w definicji. Niech będzie zbiorem normalnym w danym jako



gdzie są jak w definicji. Mamy:

Twierdzenie 11.7. [Twierdzenie Fubiniego dla zbiorów normalnych w i )]

(1) Jeśli jest funkcją ciągłą, to



(2) Jeśli jest funkcją ciągłą, to

Dowód tej wersji Twierdzenia Fubiniego można dostać jako wniosek z ogólnej wersji twierdzenia (dowodząc, że zbiory regularne są J-mierzalne) albo można udowodnić to twierdzenie bezpośrednio, nieco modyfikując dowód twierdzenia 11.3.

Możemy teraz policzyć następującą całkę.

Trójkąt

Przykład 11.8.

Policzyć całkę

gdzie jest trójkątem ograniczonym prostymi:

Zauważmy, że zbiór jest normalny względem osi Ponieważ jednak funkcja ograniczająca ten zbiór od dołu jest sklejeniem dwóch funkcji ( oraz ), to wygodniej będzie podzielić na dwa zbiory normalne (o rozłącznych wnętrzach). Pierwszy z tych zbiorów to trójkąt ograniczony prostymi: a drugi to trójkąt ograniczony prostymi: jest więc zbiorem regularnym. Z twierdzenia Fubiniego mamy:


<applet code="JavaviewModApplet.class" archive="images/a/a0/Javaview.jar,images/b/be/JavaviewModApplet.jar" width="450" height="400">

   <param name="colors" value="-1:#7e0b9f -0.80:#ce3bf8 -0.60:#2c4ae5 -0.40:#2c85e5 -0.20:#2ecca9 -0.05:#2ecc5b 0.05:#2ecc5b 0.20:#97cc2e 0.40:#edff27 0.60:#ffba27 0.80:#ff6e27 1:#d42525">
   <param name="coloring" value="maple">
   <param name="model" value="images/d/d6/Am2.m11.w.r07.mgs.zip">
   <param name="scale" value="1.0 1.0 0.1">
   <param name="shading" value="0.2">
   <param name="animation" value="stop">

</applet>

<div.thumbcaption>Wykres funkcji nad


Twierdzenie o zmianie zmiennych

Jeszcze jedno twierdzenie bardzo nam się przyda do liczenia całek wielowymiarowych. Jest to uogólnienie na więcej wymiarów znanego już z teorii całki jednej zmiennej twierdzenia o całkowaniu przez podstawienie. W przypadku wielowymiarowym nosi ono nazwę twierdzenia o zmianie zmiennych.

Załóżmy, że mamy zbiory J-mierzalne i w oraz odwzorowanie które jest -dyfeomorfizmem (to znaczy, że jest bijekcją klasy i odwzorowanie odwrotne do też jest tej klasy). Dla odwzorowania możemy wypisać macierz Jacobiego, czyli macierz pochodnych cząstkowych (w punkcie ):

Jac

Wyznacznik tej macierzy (w punkcie ) nazywamy jakobianem w punkcie . Gdy jest dyfeomorfizmem, to Jac .

Współrzędne w zbiorze oznaczmy przez

Twierdzenie o zmianie zmiennych brzmi następująco.

Twierdzenie 11.9. [Twierdzenie o zmianie zmiennych]

Przy oznaczeniach i założeniach jak wyżej, niech będzie funkcją ciągłą. Wtedy

Jac
Uwaga 11.10.

Zauważmy, że dla dostajemy znane twierdzenie o całkowaniu przez podstawienie:

Dowód twierdzenia 11.9. pomijamy. Przedstawimy natomiast kilka użytecznych przykładów.

Uwaga 11.11.

W powyższym twierdzeniu nie trzeba zakładać, że odwzorowanie jest dyfeomorfizmem na całym zbiorze wystarczy założyć, że istnieje podzbiór taki, że oraz jest dyfeomorfizmem.

Zmiana zmiennych na dwuwymiarowe współrzędne biegunowe

Współrzędne biegunowe

Niech zbiorem będzie Określamy odwzorowanie prowadzące ze zbioru następująco:

gdzie najczęściej zapisujemy jako


Tak więc a zatem jest odległością punktu od początku układu współrzędnych. Kąt jest kątem, jaki tworzy wektor o początku w i końcu w z dodatnią częścią osi

Licząc jakobian tej zmiany zmiennych dostajemy Jac (trzeba policzyć pochodne cząstkowe i po i a następnie wyznacznik macierzy Jacobiego). Tak więc tu jakobian jest zawsze dodatni.

Tę zmianę zmiennych stosujemy najczęściej, gdy obszarem całkowania (zbiorem ) jest koło, pierścień lub ich wycinek. Jak wtedy wygląda zbiór obrazują przykłady poniżej.

W dalszych rozważaniach najczęściej nie będziemy rozróżniać pomiędzy i (lub i ) gdzie () i, choć nie jest to w pełni poprawne, będziemy pisać o zmianie zmiennych z do a nie z do ignorując fakt, że zmiana zmiennych może nie być dyfeomorfizmem na jakimś zbiorze miary zero.

Przykład 11.12.

Policzyć całkę

gdzie jest kołem o promieniu i środku w punkcie zatem

Skoro to promień zmienia się w przedziale a kąt zmienia się w całym zakresie

Tak więc czyli mamy

gdzie pierwsza równość zachodzi na podstawie twierdzenia o zmianie

zmiennych, a druga na podstawie twierdzenia Fubiniego.

Przykład 11.13.

Policzyć całkę

gdzie jest ćwiartką koła o promieniu i środku w punkcie leżącą w drugiej ćwiartce płaszczyzny.

Stosujemy taką samą zmianę zmiennych. Tym razem zmienia się także od do natomiast zmienia się od do . Tak więc


Zbiór
Zbiór

Zmiana zmiennych na trójwymiarowe współrzędne sferyczne

Współrzędne sferyczne

Podobnie do współrzędnych biegunowych w definiujemy współrzędne sferyczne w Mamy:

gdzie

Teraz jest odległością punktu od początku układu współrzędnych, jest kątem, jaki tworzy wektor z dodatnią częścią osi a jest kątem, jaki tworzy wektor z dodatnią częścią osi

Jakobian tej zmiany zmiennych wynosi , a zatem jest dodatni, bo

Przykład 11.14.

Policz całkę

gdzie jest górną połową kuli o środku w i promieniu

Kula opisana jest nierównością w takim razie zmienia się w przedziale Górną połowę kuli zadaje nierówność zatem musi być czyli a zatem Na nie mamy żadnych dodatkowych warunków, więc Zatem Tak więc

Zmiana zmiennych na trójwymiarowe współrzędne walcowe

AM2.M11.W.R12

Ta zmiana zmiennych jest w zasadzie zmianą na współrzędne biegunowe w Opisana jest wzorami:

gdzie Jakobian tej zmiany zmiennych wynosi

Przykład 11.15.

Policzyć całkę

gdzie jest walcem o podstawie i o wysokości

Skoro to na kąt nie mamy dodatkowych warunków, natomiast skoro wysokość walca wynosi to Tak więc

Ciekawsze przykłady policzymy na ćwiczeniach.