หน้าแรก คลังความรู้ แคลคูลัสหลายตัวแปร อินทิกรัลหลายชั้น
Multivariable & Vector Calculus ดูครบทุกส่วน · 13 หน้า 4 · ฟังก์ชันเวกเตอร์และแคลคูลัสเวกเตอร์ 5 · การประยุกต์ทางวิศวกรรม
// บทที่ 24.5 · แคลคูลัสหลายตัวแปร · Double & Triple Integrals
อินทิกรัลสองชั้นและสามชั้น (Double & Triple Integrals): พิกัดขั้ว ทรงกระบอก ทรงกลม และ Jacobian
อินทิกรัลตัวแปรเดียวรวม “พื้นที่ใต้เส้นโค้ง” ส่วนอินทิกรัลสองชั้น (double integral) รวม “ปริมาตรใต้ผิว” และอินทิกรัลสามชั้น (triple integral) รวมปริมาณใด ๆ ทั่วเนื้อของวัตถุ ไม่ว่าจะเป็นมวล ประจุ พลังงาน หรือโมเมนต์ความเฉื่อยของชิ้นส่วนเครื่องจักร หน้านี้สอนทั้งวิธีคำนวณทีละชั้น การเลือกบริเวณ การสลับลำดับ การเปลี่ยนไปใช้พิกัดขั้ว ทรงกระบอก ทรงกลมพร้อมตัวคูณ Jacobian และการประยุกต์หามวลกับจุดศูนย์กลางมวล พร้อมแล็บที่เห็นทั้งผลรวมรีมันน์และค่าที่แท้จริง
สรุปใน 30 วินาที
∬ R f d A = \iint_R f\,dA = ∬ R f d A = ปริมาตรใต้ผิว (เมื่อ f ≥ 0 f \ge 0 f ≥ 0 ) · คำนวณทีละชั้น (Fubini): ในสุดก่อน ถือตัวนอกเป็นค่าคงที่
บริเวณ type I a ≤ x ≤ b , g 1 ( x ) ≤ y ≤ g 2 ( x ) a \le x \le b,\ g_1(x) \le y \le g_2(x) a ≤ x ≤ b , g 1 ( x ) ≤ y ≤ g 2 ( x ) · type II c ≤ y ≤ d , h 1 ( y ) ≤ x ≤ h 2 ( y ) c \le y \le d,\ h_1(y) \le x \le h_2(y) c ≤ y ≤ d , h 1 ( y ) ≤ x ≤ h 2 ( y ) · อินทิเกรตไม่ออกให้สลับลำดับ
พิกัดขั้ว: d A = r d r d θ dA = r\,dr\,d\theta d A = r d r d θ · ทรงกระบอก: d V = r d z d r d θ dV = r\,dz\,dr\,d\theta d V = r d z d r d θ · ทรงกลม: d V = ρ 2 sin φ d ρ d φ d θ dV = \rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta d V = ρ 2 sin φ d ρ d φ d θ · ตัวคูณคือ ∣ J ∣ |J| ∣ J ∣
ประยุกต์: มวล M = ∬ ρ d A M = \iint\rho\,dA M = ∬ ρ d A · x ˉ = 1 M ∬ x ρ d A \bar x = \tfrac1M\iint x\rho\,dA x ˉ = M 1 ∬ x ρ d A · I x = ∬ y 2 ρ d A I_x = \iint y^2\rho\,dA I x = ∬ y 2 ρ d A
ระดับ พื้นฐาน ปี 1 · แคลคูลัส 1กลาง แคลคูลัส 3 · คณิตศาสตร์วิศวกรรมขั้นสูง พิสูจน์ · ทฤษฎี
ค้นหาสูตรในหน้านี้
โหมดฝึกทำ: ซ่อนวิธีทำ
อินทิกรัลสองชั้นคืออะไร: ผลรวมรีมันน์ในสองมิติ พื้นฐาน
ตัวแปรเดียว: แบ่งช่วงเป็นชิ้นแคบ Δ x \Delta x Δ x คูณความสูง f f f แล้วบวก สองตัวแปร: แบ่งบริเวณ R R R เป็นช่องเล็กพื้นที่ Δ A \Delta A Δ A แต่ละช่องมี “แท่ง” สูง f ( x i ∗ , y j ∗ ) f(x_i^*, y_j^*) f ( x i ∗ , y j ∗ ) ปริมาตรรวมของแท่งทั้งหมดคือผลรวมรีมันน์ (Riemann sum) เมื่อช่องเล็กลงจนไม่สิ้นสุด ผลรวมเข้าใกล้อินทิกรัลสองชั้น ในแล็บด้านล่างเลื่อน n n n จาก 1 ถึง 40 แล้วดูแท่งปริมาตรกับค่า “ที่แท้จริง” วิ่งเข้าหากัน
M1 Double integral (Riemann sum) อินทิกรัลสองชั้น (ผลรวมรีมันน์)
∬ R f ( x , y ) d A = lim m , n → ∞ ∑ i = 1 m ∑ j = 1 n f ( x i ∗ , y j ∗ ) Δ A , Δ A = Δ x Δ y \iint_R f(x, y)\,dA = \lim_{m,n\to\infty}\sum_{i=1}^{m}\sum_{j=1}^{n} f(x_i^{*}, y_j^{*})\,\Delta A,\qquad \Delta A = \Delta x\,\Delta y ∬ R f ( x , y ) d A = m , n → ∞ lim i = 1 ∑ m j = 1 ∑ n f ( x i ∗ , y j ∗ ) Δ A , Δ A = Δ x Δ y
f ≥ 0 f \ge 0 f ≥ 0 : ปริมาตรใต้ผิว z = f ( x , y ) z = f(x, y) z = f ( x , y ) เหนือ R R R · จุด ( x i ∗ , y j ∗ ) (x_i^*, y_j^*) ( x i ∗ , y j ∗ ) เลือกกึ่งกลางช่อง = midpoint rule · f = 1 f = 1 f = 1 ให้พื้นที่ A ( R ) A(R) A ( R )
M2 Fubini's theorem (rectangle) ทฤษฎีบทฟูบินี (สี่เหลี่ยมผืนผ้า)
∬ [ a , b ] × [ c , d ] f d A = ∫ a b ( ∫ c d f ( x , y ) d y ) d x = ∫ c d ( ∫ a b f ( x , y ) d x ) d y \iint_{[a,b]\times[c,d]} f\,dA = \int_a^b\!\left(\int_c^d f(x, y)\,dy\right)dx = \int_c^d\!\left(\int_a^b f(x, y)\,dx\right)dy ∬ [ a , b ] × [ c , d ] f d A = ∫ a b ( ∫ c d f ( x , y ) d y ) d x = ∫ c d ( ∫ a b f ( x , y ) d x ) d y
เมื่อ f f f ต่อเนื่อง สองลำดับให้ค่าเท่ากัน · ชั้นในอินทิเกรตก่อน โดยถือตัวแปรนอกเป็นค่าคงที่ · ถ้า f = g ( x ) h ( y ) f = g(x)h(y) f = g ( x ) h ( y ) แยกเป็นผลคูณ ∫ g d x ⋅ ∫ h d y \int g\,dx\cdot\int h\,dy ∫ g d x ⋅ ∫ h d y
M3 Average value ค่าเฉลี่ยของฟังก์ชัน
f avg = 1 A ( R ) ∬ R f d A , A ( R ) = ∬ R d A f_{\text{avg}} = \frac{1}{A(R)}\iint_R f\,dA,\qquad A(R) = \iint_R dA f avg = A ( R ) 1 ∬ R f d A , A ( R ) = ∬ R d A
เชิงเส้นและบวกกันได้ตามบริเวณ · ค่าเฉลี่ยทางวิศวกรรม เช่น อุณหภูมิเฉลี่ยของแผ่น ความสูงเฉลี่ยของพื้นที่ (ใช้หาระดับสมดุลในบทโยธา )
ตัวอย่าง M-1 · Riemann sum กับค่าที่แท้จริง หา ∬ R ( x 2 + y 2 ) d A \displaystyle\iint_R (x^2 + y^2)\,dA ∬ R ( x 2 + y 2 ) d A เมื่อ R = [ 0 , 1 ] × [ 0 , 2 ] R = [0, 1]\times[0, 2] R = [ 0 , 1 ] × [ 0 , 2 ] และเทียบกับผลรวมรีมันน์แบบ 2×2 ช่อง (กึ่งกลาง) เปิดในแล็บ (n = 2)
ชั้นใน (ตาม y y y , ถือ x x x คงที่): ∫ 0 2 ( x 2 + y 2 ) d y = [ x 2 y + y 3 3 ] 0 2 = 2 x 2 + 8 3 \displaystyle\int_0^2 (x^2 + y^2)\,dy = \left[x^2 y + \dfrac{y^3}{3}\right]_0^2 = 2x^2 + \dfrac83 ∫ 0 2 ( x 2 + y 2 ) d y = [ x 2 y + 3 y 3 ] 0 2 = 2 x 2 + 3 8
ชั้นนอก: ∫ 0 1 ( 2 x 2 + 8 3 ) d x = 2 3 + 8 3 = 10 3 ≈ 3.3333 \displaystyle\int_0^1\left(2x^2 + \dfrac83\right)dx = \dfrac{2}{3} + \dfrac83 = \dfrac{10}{3} \approx 3.3333 ∫ 0 1 ( 2 x 2 + 3 8 ) d x = 3 2 + 3 8 = 3 10 ≈ 3.3333
Riemann sum 2×2: Δ x = 0.5 , Δ y = 1 , Δ A = 0.5 \Delta x = 0.5,\ \Delta y = 1,\ \Delta A = 0.5 Δ x = 0.5 , Δ y = 1 , Δ A = 0.5 · จุดกึ่งกลาง x ∈ { 0.25 , 0.75 } x \in \{0.25, 0.75\} x ∈ { 0.25 , 0.75 } , y ∈ { 0.5 , 1.5 } y \in \{0.5, 1.5\} y ∈ { 0.5 , 1.5 }
ค่า f f f : 0.3125 , 0.8125 , 2.3125 , 2.8125 0.3125,\ 0.8125,\ 2.3125,\ 2.8125 0.3125 , 0.8125 , 2.3125 , 2.8125 · ผลรวม = 6.25 = 6.25 = 6.25 · × 0.5 = 3.125 \times\,0.5 = 3.125 × 0.5 = 3.125 (ต่ำกว่าจริง 6.25%) · เพิ่ม n n n เป็น 6 ได้ 3.3102 3.3102 3.3102 (ต่ำกว่า 0.7%)
∬ R ( x 2 + y 2 ) d A = 10 3 \iint_R (x^2 + y^2)\,dA = \dfrac{10}{3} ∬ R ( x 2 + y 2 ) d A = 3 10 · midpoint 2×2 ได้ 3.125
บริเวณที่ไม่ใช่สี่เหลี่ยม: type I, type II และการสลับลำดับ พื้นฐาน
บริเวณในงานจริงแทบไม่เป็นสี่เหลี่ยมผืนผ้า ให้วาดบริเวณก่อน แล้วเลือกว่าจะ “กวาด” ด้วยแถบแนวตั้งหรือแนวนอน: type I ขอบบนและล่างเป็นฟังก์ชันของ x x x (แถบตั้ง) · type II ขอบซ้ายและขวาเป็นฟังก์ชันของ y y y (แถบนอน) ขอบเขตของชั้นในได้เป็นฟังก์ชันของตัวแปรนอก ได้ แต่ขอบเขตของชั้นนอกต้องเป็นค่าคงที่เสมอ
x y type I: แถบแนวตั้ง g₁(x) ≤ y ≤ g₂(x) y = x² y = x x y type II: แถบแนวนอน h₁(y) ≤ x ≤ h₂(y) x = √y x = y ∫₀¹ ∫ f dy dx (y: x² → x) ∫₀¹ ∫ f dx dy (x: y → √y) บริเวณ type I (แถบแนวตั้ง: y y y ระหว่าง g 1 ( x ) g_1(x) g 1 ( x ) กับ g 2 ( x ) g_2(x) g 2 ( x ) ) และ type II (แถบแนวนอน: x x x ระหว่าง h 1 ( y ) h_1(y) h 1 ( y ) กับ h 2 ( y ) h_2(y) h 2 ( y ) ) ของบริเวณเดียวกัน — การสลับลำดับการอินทิเกรตคือการเปลี่ยนวิธีแบ่ง
M4 Type I region บริเวณ type I
R = { a ≤ x ≤ b , g 1 ( x ) ≤ y ≤ g 2 ( x ) } : ∬ R f d A = ∫ a b ∫ g 1 ( x ) g 2 ( x ) f ( x , y ) d y d x R = \{a \le x \le b,\ g_1(x) \le y \le g_2(x)\}:\qquad \iint_R f\,dA = \int_a^b\!\int_{g_1(x)}^{g_2(x)} f(x, y)\,dy\,dx R = { a ≤ x ≤ b , g 1 ( x ) ≤ y ≤ g 2 ( x )} : ∬ R f d A = ∫ a b ∫ g 1 ( x ) g 2 ( x ) f ( x , y ) d y d x
แถบตั้ง: ลากเส้นตั้งผ่านบริเวณ เส้นเข้าที่ g 1 g_1 g 1 ออกที่ g 2 g_2 g 2 · ขอบชั้นนอก a , b a, b a , b เป็นค่าคงที่
M5 Type II region บริเวณ type II
R = { c ≤ y ≤ d , h 1 ( y ) ≤ x ≤ h 2 ( y ) } : ∬ R f d A = ∫ c d ∫ h 1 ( y ) h 2 ( y ) f ( x , y ) d x d y R = \{c \le y \le d,\ h_1(y) \le x \le h_2(y)\}:\qquad \iint_R f\,dA = \int_c^d\!\int_{h_1(y)}^{h_2(y)} f(x, y)\,dx\,dy R = { c ≤ y ≤ d , h 1 ( y ) ≤ x ≤ h 2 ( y )} : ∬ R f d A = ∫ c d ∫ h 1 ( y ) h 2 ( y ) f ( x , y ) d x d y
แถบนอน: ลากเส้นนอนผ่านบริเวณ เส้นเข้าที่ h 1 h_1 h 1 ออกที่ h 2 h_2 h 2 · เปลี่ยนลำดับ = เปลี่ยนจาก type I เป็น type II (วาดรูปแล้วอ่านขอบใหม่)
ตัวอย่าง M-2 · บริเวณสามเหลี่ยม หา ∬ R x y d A \iint_R xy\,dA ∬ R x y d A เมื่อ R R R เป็นสามเหลี่ยมที่มีจุดยอด ( 0 , 0 ) , ( 2 , 0 ) , ( 0 , 2 ) (0,0),\ (2,0),\ (0,2) ( 0 , 0 ) , ( 2 , 0 ) , ( 0 , 2 ) เปิดในแล็บ
วาดบริเวณ: เส้นตรงจาก ( 2 , 0 ) (2,0) ( 2 , 0 ) ถึง ( 0 , 2 ) (0,2) ( 0 , 2 ) คือ y = 2 − x y = 2 - x y = 2 − x · type I: 0 ≤ x ≤ 2 , 0 ≤ y ≤ 2 − x 0 \le x \le 2,\ 0 \le y \le 2 - x 0 ≤ x ≤ 2 , 0 ≤ y ≤ 2 − x
ชั้นใน: ∫ 0 2 − x x y d y = x ⋅ ( 2 − x ) 2 2 \displaystyle\int_0^{2-x} xy\,dy = x\cdot\dfrac{(2-x)^2}{2} ∫ 0 2 − x x y d y = x ⋅ 2 ( 2 − x ) 2
ชั้นนอก: 1 2 ∫ 0 2 x ( 4 − 4 x + x 2 ) d x = 1 2 [ 2 x 2 − 4 x 3 3 + x 4 4 ] 0 2 = 1 2 ( 8 − 32 3 + 4 ) = 1 2 ⋅ 4 3 = 2 3 \displaystyle\dfrac12\int_0^2 x(4 - 4x + x^2)\,dx = \dfrac12\left[2x^2 - \dfrac{4x^3}{3} + \dfrac{x^4}{4}\right]_0^2 = \dfrac12\left(8 - \dfrac{32}{3} + 4\right) = \dfrac12\cdot\dfrac43 = \dfrac23 2 1 ∫ 0 2 x ( 4 − 4 x + x 2 ) d x = 2 1 [ 2 x 2 − 3 4 x 3 + 4 x 4 ] 0 2 = 2 1 ( 8 − 3 32 + 4 ) = 2 1 ⋅ 3 4 = 3 2
ตรวจด้วยอีกลำดับ (type II: 0 ≤ y ≤ 2 , 0 ≤ x ≤ 2 − y 0 \le y \le 2,\ 0 \le x \le 2 - y 0 ≤ y ≤ 2 , 0 ≤ x ≤ 2 − y ): โดยสมมาตรของ x y xy x y ได้ 2 3 \tfrac23 3 2 เท่ากัน ✓
∬ R x y d A = 2 3 \iint_R xy\,dA = \dfrac23 ∬ R x y d A = 3 2
ตัวอย่าง M-3 · สลับลำดับเมื่อชั้นในอินทิเกรตไม่ได้ คำนวณ ∫ 0 1 ∫ x 1 e y 2 d y d x \displaystyle\int_0^1\!\int_x^1 e^{y^2}\,dy\,dx ∫ 0 1 ∫ x 1 e y 2 d y d x
∫ e y 2 d y \int e^{y^2}\,dy ∫ e y 2 d y ไม่มีรูปปิดจากฟังก์ชันพื้นฐาน จึงอินทิเกรตตาม y y y ก่อนไม่ได้
อ่านบริเวณ: 0 ≤ x ≤ 1 0 \le x \le 1 0 ≤ x ≤ 1 และ x ≤ y ≤ 1 x \le y \le 1 x ≤ y ≤ 1 คือสามเหลี่ยมใต้เส้น y = 1 y = 1 y = 1 เหนือเส้น y = x y = x y = x (จุดยอด ( 0 , 0 ) , ( 0 , 1 ) , ( 1 , 1 ) (0,0),(0,1),(1,1) ( 0 , 0 ) , ( 0 , 1 ) , ( 1 , 1 ) )
เปลี่ยนเป็น type II: 0 ≤ y ≤ 1 0 \le y \le 1 0 ≤ y ≤ 1 และ 0 ≤ x ≤ y 0 \le x \le y 0 ≤ x ≤ y · ∫ 0 1 ∫ 0 y e y 2 d x d y \displaystyle\int_0^1\!\int_0^y e^{y^2}\,dx\,dy ∫ 0 1 ∫ 0 y e y 2 d x d y
ชั้นใน (ตาม x x x , e y 2 e^{y^2} e y 2 คงที่): y e y 2 y\,e^{y^2} y e y 2 · ชั้นนอก: ∫ 0 1 y e y 2 d y = [ 1 2 e y 2 ] 0 1 = e − 1 2 ≈ 0.8591 \displaystyle\int_0^1 y\,e^{y^2}\,dy = \left[\tfrac12 e^{y^2}\right]_0^1 = \dfrac{e - 1}{2} \approx 0.8591 ∫ 0 1 y e y 2 d y = [ 2 1 e y 2 ] 0 1 = 2 e − 1 ≈ 0.8591
e − 1 2 ≈ 0.8591 \dfrac{e - 1}{2} \approx 0.8591 2 e − 1 ≈ 0.8591
ตัวอย่าง M-4 · พื้นที่ระหว่างเส้นโค้ง หาพื้นที่ที่ล้อมด้วย y = x 2 y = x^2 y = x 2 และ y = x y = x y = x เปิดในแล็บ
จุดตัด: x 2 = x x^2 = x x 2 = x → x = 0 , 1 x = 0, 1 x = 0 , 1 · บน 0 < x < 1 0 < x < 1 0 < x < 1 เส้นตรง y = x y = x y = x อยู่เหนือพาราโบลา
A = ∫ 0 1 ∫ x 2 x d y d x = ∫ 0 1 ( x − x 2 ) d x = 1 2 − 1 3 = 1 6 A = \displaystyle\int_0^1\!\int_{x^2}^{x} dy\,dx = \int_0^1 (x - x^2)\,dx = \dfrac12 - \dfrac13 = \dfrac16 A = ∫ 0 1 ∫ x 2 x d y d x = ∫ 0 1 ( x − x 2 ) d x = 2 1 − 3 1 = 6 1
A = 1 6 A = \dfrac16 A = 6 1 (นี่คืออินทิกรัลสองชั้นของ f = 1 f = 1 f = 1 · ตรงกับอินทิกรัลตัวแปรเดียว ∫ ( บน − ล่าง ) d x \int (\text{บน} - \text{ล่าง})\,dx ∫ ( บน − ล่าง ) d x )
พิกัดขั้ว (Polar Coordinates) ในอินทิกรัลสองชั้น กลาง
บริเวณวงกลม ครึ่งวงกลม เซกเตอร์ และวงแหวนเขียนในพิกัดฉากยากมาก แต่ในพิกัดขั้วเรียบง่าย: 0 ≤ r ≤ R , 0 ≤ θ ≤ 2 π 0 \le r \le R,\ 0 \le \theta \le 2\pi 0 ≤ r ≤ R , 0 ≤ θ ≤ 2 π กฎคือแทน x = r cos θ , y = r sin θ x = r\cos\theta,\ y = r\sin\theta x = r cos θ , y = r sin θ แล้วเปลี่ยน d A = r d r d θ dA = r\,dr\,d\theta d A = r d r d θ ตัวคูณ r r r มาจากการที่ “ช่องขั้ว” กว้าง Δ r \Delta r Δ r และยาวเป็นส่วนโค้ง r Δ θ r\,\Delta\theta r Δ θ ยิ่งไกลจุดกำเนิด ช่องเดียวกันใน ( r , θ ) (r,\theta) ( r , θ ) ยิ่งใหญ่
x y θ Δr r Δθ θ + Δθ θ ความกว้าง Δr ความยาวของส่วนโค้ง r Δθ ΔA ≈ r Δr Δθ ⇒ dA = r dr dθ ตัวคูณ r = จาโคเบียนของ (r, θ) ช่องในพิกัดขั้ว (polar rectangle): กว้าง Δ r \Delta r Δ r ยาว r Δ θ r\,\Delta\theta r Δ θ จึงมีพื้นที่ Δ A ≈ r Δ r Δ θ \Delta A \approx r\,\Delta r\,\Delta\theta Δ A ≈ r Δ r Δ θ ซึ่งเป็นที่มาของตัวคูณ r r r
M6 Double integral in polar coordinates อินทิกรัลสองชั้นในพิกัดขั้ว
∬ R f ( x , y ) d A = ∫ α β ∫ r 1 ( θ ) r 2 ( θ ) f ( r cos θ , r sin θ ) r ⏟ Jacobian d r d θ \iint_R f(x, y)\,dA = \int_{\alpha}^{\beta}\!\int_{r_1(\theta)}^{r_2(\theta)} f(r\cos\theta,\ r\sin\theta)\;\underbrace{r}_{\text{Jacobian}}\;dr\,d\theta ∬ R f ( x , y ) d A = ∫ α β ∫ r 1 ( θ ) r 2 ( θ ) f ( r cos θ , r sin θ ) Jacobian r d r d θ
ใช้เมื่อบริเวณหรือ f f f มี x 2 + y 2 x^2 + y^2 x 2 + y 2 · r 2 = x 2 + y 2 r^2 = x^2 + y^2 r 2 = x 2 + y 2 · พื้นที่ A = 1 2 ∫ ( r 2 2 − r 1 2 ) d θ A = \tfrac12\int (r_2^2 - r_1^2)\,d\theta A = 2 1 ∫ ( r 2 2 − r 1 2 ) d θ · ระวังขอบเขต θ \theta θ ของบริเวณ (ครึ่งบน 0 ≤ θ ≤ π 0 \le \theta \le \pi 0 ≤ θ ≤ π , ควอดแรนต์แรก 0 ≤ θ ≤ π / 2 0 \le \theta \le \pi/2 0 ≤ θ ≤ π /2 )
ลองในแล็บ: โหมดขั้วใช้ r r r ระหว่าง r 1 ( θ ) r_1(\theta) r 1 ( θ ) ถึง r 2 ( θ ) r_2(\theta) r 2 ( θ ) และ θ \theta θ ระหว่างสองค่า เลือกพรีเซ็ต “พื้นที่วงกลม” “ปริมาตรพาราโบลอยด์” หรือ “ใบกลีบ r = 2 sin 2 θ r = 2\sin 2\theta r = 2 sin 2 θ ” แล้วดูแผนผังช่องที่เป็นรูปพัด
ตัวอย่าง M-5 · ปริมาตรพาราโบลอยด์ หาปริมาตรใต้ผิว z = 4 − x 2 − y 2 z = 4 - x^2 - y^2 z = 4 − x 2 − y 2 เหนือระนาบ x y xy x y เปิดในแล็บ
ผิวตัดระนาบ z = 0 z = 0 z = 0 ที่ x 2 + y 2 = 4 x^2 + y^2 = 4 x 2 + y 2 = 4 · บริเวณคือแผ่นวงกลมรัศมี 2: 0 ≤ r ≤ 2 , 0 ≤ θ ≤ 2 π 0 \le r \le 2,\ 0 \le \theta \le 2\pi 0 ≤ r ≤ 2 , 0 ≤ θ ≤ 2 π
แทน x 2 + y 2 = r 2 x^2 + y^2 = r^2 x 2 + y 2 = r 2 : V = ∫ 0 2 π ∫ 0 2 ( 4 − r 2 ) r d r d θ V = \displaystyle\int_0^{2\pi}\!\int_0^2 (4 - r^2)\,r\,dr\,d\theta V = ∫ 0 2 π ∫ 0 2 ( 4 − r 2 ) r d r d θ
ชั้นใน: ∫ 0 2 ( 4 r − r 3 ) d r = [ 2 r 2 − r 4 4 ] 0 2 = 8 − 4 = 4 \displaystyle\int_0^2 (4r - r^3)\,dr = \left[2r^2 - \dfrac{r^4}{4}\right]_0^2 = 8 - 4 = 4 ∫ 0 2 ( 4 r − r 3 ) d r = [ 2 r 2 − 4 r 4 ] 0 2 = 8 − 4 = 4
ชั้นนอก: ∫ 0 2 π 4 d θ = 8 π ≈ 25.13 \displaystyle\int_0^{2\pi} 4\,d\theta = 8\pi \approx 25.13 ∫ 0 2 π 4 d θ = 8 π ≈ 25.13
V = 8 π ≈ 25.13 V = 8\pi \approx 25.13 V = 8 π ≈ 25.13 หน่วยปริมาตร (ตรงกับค่าที่แล็บคำนวณ)
ตัวอย่าง M-6 · อินทิกรัลเกาส์เซียน (ไฟฟ้า) แสดงว่า ∫ − ∞ ∞ e − x 2 d x = π \displaystyle\int_{-\infty}^{\infty} e^{-x^2}dx = \sqrt{\pi} ∫ − ∞ ∞ e − x 2 d x = π ด้วยอินทิกรัลสองชั้น
ให้ I = ∫ − ∞ ∞ e − x 2 d x I = \int_{-\infty}^{\infty} e^{-x^2}dx I = ∫ − ∞ ∞ e − x 2 d x · I 2 = ∫ e − x 2 d x ∫ e − y 2 d y = ∬ R 2 e − ( x 2 + y 2 ) d A I^2 = \int e^{-x^2}dx\int e^{-y^2}dy = \iint_{\mathbb R^2} e^{-(x^2+y^2)}\,dA I 2 = ∫ e − x 2 d x ∫ e − y 2 d y = ∬ R 2 e − ( x 2 + y 2 ) d A
พิกัดขั้ว: ทั่วระนาบคือ 0 ≤ r < ∞ , 0 ≤ θ ≤ 2 π 0 \le r < \infty,\ 0 \le \theta \le 2\pi 0 ≤ r < ∞ , 0 ≤ θ ≤ 2 π · I 2 = ∫ 0 2 π ∫ 0 ∞ e − r 2 r d r d θ I^2 = \displaystyle\int_0^{2\pi}\!\int_0^{\infty} e^{-r^2}\,r\,dr\,d\theta I 2 = ∫ 0 2 π ∫ 0 ∞ e − r 2 r d r d θ
∫ 0 ∞ r e − r 2 d r = [ − 1 2 e − r 2 ] 0 ∞ = 1 2 \displaystyle\int_0^\infty r e^{-r^2}dr = \left[-\tfrac12 e^{-r^2}\right]_0^\infty = \tfrac12 ∫ 0 ∞ r e − r 2 d r = [ − 2 1 e − r 2 ] 0 ∞ = 2 1 · I 2 = 2 π ⋅ 1 2 = π I^2 = 2\pi\cdot\tfrac12 = \pi I 2 = 2 π ⋅ 2 1 = π → I = π I = \sqrt{\pi} I = π
ผลนี้คือตัวทำให้เป็นปกติ (normalizing constant) ของการแจกแจงเกาส์เซียน 1 σ 2 π e − x 2 / 2 σ 2 \dfrac{1}{\sigma\sqrt{2\pi}}e^{-x^2/2\sigma^2} σ 2 π 1 e − x 2 /2 σ 2 ที่ใช้อธิบายสัญญาณรบกวนในระบบสื่อสารและการวัด
∫ − ∞ ∞ e − x 2 d x = π \int_{-\infty}^{\infty} e^{-x^2}dx = \sqrt{\pi} ∫ − ∞ ∞ e − x 2 d x = π (ใช้ตัวคูณ r r r จึงอินทิเกรตได้ — ในพิกัดฉากไม่มีรูปปิด)
ตัวอย่าง M-7 · พื้นที่กลีบกุหลาบ หาพื้นที่หนึ่งกลีบของ r = 2 sin 2 θ r = 2\sin 2\theta r = 2 sin 2 θ (0 ≤ θ ≤ π / 2 0 \le \theta \le \pi/2 0 ≤ θ ≤ π /2 ) เปิดในแล็บ
A = ∫ 0 π / 2 ∫ 0 2 sin 2 θ r d r d θ = ∫ 0 π / 2 ( 2 sin 2 θ ) 2 2 d θ = 2 ∫ 0 π / 2 sin 2 2 θ d θ A = \displaystyle\int_0^{\pi/2}\!\int_0^{2\sin 2\theta} r\,dr\,d\theta = \int_0^{\pi/2}\dfrac{(2\sin2\theta)^2}{2}\,d\theta = 2\int_0^{\pi/2}\sin^2 2\theta\,d\theta A = ∫ 0 π /2 ∫ 0 2 s i n 2 θ r d r d θ = ∫ 0 π /2 2 ( 2 sin 2 θ ) 2 d θ = 2 ∫ 0 π /2 sin 2 2 θ d θ
sin 2 2 θ = 1 2 ( 1 − cos 4 θ ) \sin^2 2\theta = \tfrac12(1 - \cos 4\theta) sin 2 2 θ = 2 1 ( 1 − cos 4 θ ) → 2 ⋅ 1 2 [ θ − sin 4 θ 4 ] 0 π / 2 = π 2 2\cdot\tfrac12\left[\theta - \tfrac{\sin4\theta}{4}\right]_0^{\pi/2} = \dfrac{\pi}{2} 2 ⋅ 2 1 [ θ − 4 s i n 4 θ ] 0 π /2 = 2 π
A = π 2 ≈ 1.571 A = \dfrac{\pi}{2} \approx 1.571 A = 2 π ≈ 1.571
อินทิกรัลสามชั้น: พิกัดฉาก ทรงกระบอก และทรงกลม กลาง
เมื่อเพิ่มมิติ z z z ได้ ∭ E f d V \iiint_E f\,dV ∭ E f d V ที่รวม f f f ทั่วบริเวณตัน E E E ถ้า f = 1 f = 1 f = 1 ได้ปริมาตร ถ้า f = ρ f = \rho f = ρ (ความหนาแน่น) ได้มวล เลือกพิกัดตามรูปทรง: ฉาก (กล่อง ปริซึม เตตระฮีดรอน) · ทรงกระบอก (ท่อ ถัง กรวย พาราโบลอยด์) · ทรงกลม (ลูกบอล เปลือก ซีกทรงกลม)
ทรงกระบอก (cylindrical) x y z P (r, θ, z) r θ z x = r cos θ, y = r sin θ, dV = r dr dθ dz ทรงกลม (spherical) x y z P (ρ, φ, θ) ρ φ θ x = ρ sin φ cos θ, dV = ρ² sin φ dρ dφ dθ พิกัดทรงกระบอก (cylindrical coordinates) ( r , θ , z ) (r, \theta, z) ( r , θ , z ) และพิกัดทรงกลม (spherical coordinates) ( ρ , φ , θ ) (\rho, \varphi, \theta) ( ρ , φ , θ ) พร้อมความสัมพันธ์กับ ( x , y , z ) (x, y, z) ( x , y , z )
M7 Triple integral (rectangular) อินทิกรัลสามชั้น (พิกัดฉาก)
∭ E f d V = ∫ a b ∫ g 1 ( x ) g 2 ( x ) ∫ u 1 ( x , y ) u 2 ( x , y ) f ( x , y , z ) d z d y d x , V = ∭ E d V , M = ∭ E ρ d V \iiint_E f\,dV = \int_a^b\!\int_{g_1(x)}^{g_2(x)}\!\int_{u_1(x,y)}^{u_2(x,y)} f(x, y, z)\,dz\,dy\,dx,\qquad V = \iiint_E dV,\quad M = \iiint_E \rho\,dV ∭ E f d V = ∫ a b ∫ g 1 ( x ) g 2 ( x ) ∫ u 1 ( x , y ) u 2 ( x , y ) f ( x , y , z ) d z d y d x , V = ∭ E d V , M = ∭ E ρ d V
ขอบชั้นในเป็นฟังก์ชันของตัวแปรนอกได้ · วาดเงาของ E E E บนระนาบ x y xy x y เพื่อหาขอบ x x x และ y y y แล้วระบุผิวบน-ล่างของ z z z
M8 Cylindrical coordinates พิกัดทรงกระบอก
x = r cos θ , y = r sin θ , z = z ; r 2 = x 2 + y 2 ; d V = r d z d r d θ x = r\cos\theta,\ \ y = r\sin\theta,\ \ z = z;\qquad r^2 = x^2 + y^2;\qquad dV = r\,dz\,dr\,d\theta x = r cos θ , y = r sin θ , z = z ; r 2 = x 2 + y 2 ; d V = r d z d r d θ
เหมาะกับรูปที่สมมาตรรอบแกน z z z · ∭ f d V = ∭ f ( r cos θ , r sin θ , z ) r d z d r d θ \iiint f\,dV = \iiint f(r\cos\theta, r\sin\theta, z)\,r\,dz\,dr\,d\theta ∭ f d V = ∭ f ( r cos θ , r sin θ , z ) r d z d r d θ
M9 Spherical coordinates พิกัดทรงกลม
x = ρ sin φ cos θ , y = ρ sin φ sin θ , z = ρ cos φ ; ρ 2 = x 2 + y 2 + z 2 ; d V = ρ 2 sin φ d ρ d φ d θ x = \rho\sin\varphi\cos\theta,\ \ y = \rho\sin\varphi\sin\theta,\ \ z = \rho\cos\varphi;\quad \rho^2 = x^2 + y^2 + z^2;\quad dV = \rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta x = ρ sin φ cos θ , y = ρ sin φ sin θ , z = ρ cos φ ; ρ 2 = x 2 + y 2 + z 2 ; d V = ρ 2 sin φ d ρ d φ d θ
ρ ≥ 0 \rho \ge 0 ρ ≥ 0 ระยะจากจุดกำเนิด · 0 ≤ φ ≤ π 0 \le \varphi \le \pi 0 ≤ φ ≤ π มุมจากแกน + z +z + z · 0 ≤ θ ≤ 2 π 0 \le \theta \le 2\pi 0 ≤ θ ≤ 2 π มุมราบ · ระวังตำราบางเล่นสลับ θ \theta θ กับ φ \varphi φ (ฟิสิกส์ใช้แบบนี้ ตำราคณิตศาสตร์บางเล่มสลับ)
แล็บสามชั้นรับ f ( x , y , z ) f(x, y, z) f ( x , y , z ) เสมอในพิกัดฉาก โปรแกรมแทนพิกัดใหม่และคูณจาโคเบียนให้ จุดสีในภาพแสดงตัวอย่างจุดภายในรูปทรง เลือกพรีเซ็ต แล้วลองแก้ขอบเขตเอง
ตัวอย่าง M-8 · เตตระฮีดรอน (พิกัดฉาก) หาปริมาตรของรูปตันที่ล้อมด้วยระนาบพิกัดสามระนาบกับ x + y + z = 1 x + y + z = 1 x + y + z = 1 เปิดในแล็บ
เงาบนระนาบ x y xy x y : สามเหลี่ยม 0 ≤ x ≤ 1 , 0 ≤ y ≤ 1 − x 0 \le x \le 1,\ 0 \le y \le 1 - x 0 ≤ x ≤ 1 , 0 ≤ y ≤ 1 − x · z z z จาก 0 ถึง 1 − x − y 1 - x - y 1 − x − y
V = ∫ 0 1 ∫ 0 1 − x ∫ 0 1 − x − y d z d y d x = ∫ 0 1 ∫ 0 1 − x ( 1 − x − y ) d y d x = ∫ 0 1 ( 1 − x ) 2 2 d x = 1 6 V = \displaystyle\int_0^1\!\int_0^{1-x}\!\int_0^{1-x-y} dz\,dy\,dx = \int_0^1\!\int_0^{1-x}(1 - x - y)\,dy\,dx = \int_0^1 \dfrac{(1-x)^2}{2}\,dx = \dfrac16 V = ∫ 0 1 ∫ 0 1 − x ∫ 0 1 − x − y d z d y d x = ∫ 0 1 ∫ 0 1 − x ( 1 − x − y ) d y d x = ∫ 0 1 2 ( 1 − x ) 2 d x = 6 1
V = 1 6 V = \dfrac16 V = 6 1 (ตรงกับสูตรปริมาตรพีระมิด 1 3 ⋅ ฐาน 1 2 ⋅ สูง 1 \tfrac13 \cdot \text{ฐาน}\tfrac12 \cdot \text{สูง}1 3 1 ⋅ ฐาน 2 1 ⋅ สูง 1 )
ตัวอย่าง M-9 · กรวย (พิกัดทรงกระบอก) หาปริมาตรรูปตันที่อยู่เหนือกรวย z = r z = r z = r และใต้ระนาบ z = 2 z = 2 z = 2 เปิดในแล็บ
กรวยตัดระนาบ z = 2 z = 2 z = 2 ที่ r = 2 r = 2 r = 2 · ขอบเขต: 0 ≤ θ ≤ 2 π , 0 ≤ r ≤ 2 , r ≤ z ≤ 2 0 \le \theta \le 2\pi,\ 0 \le r \le 2,\ r \le z \le 2 0 ≤ θ ≤ 2 π , 0 ≤ r ≤ 2 , r ≤ z ≤ 2
V = ∫ 0 2 π ∫ 0 2 ∫ r 2 r d z d r d θ = ∫ 0 2 π ∫ 0 2 r ( 2 − r ) d r d θ V = \displaystyle\int_0^{2\pi}\!\int_0^2\!\int_r^2 r\,dz\,dr\,d\theta = \int_0^{2\pi}\!\int_0^2 r(2 - r)\,dr\,d\theta V = ∫ 0 2 π ∫ 0 2 ∫ r 2 r d z d r d θ = ∫ 0 2 π ∫ 0 2 r ( 2 − r ) d r d θ
ชั้นใน: [ r 2 − r 3 3 ] 0 2 = 4 − 8 3 = 4 3 \left[r^2 - \dfrac{r^3}{3}\right]_0^2 = 4 - \dfrac83 = \dfrac43 [ r 2 − 3 r 3 ] 0 2 = 4 − 3 8 = 3 4 · ชั้นนอก: 2 π ⋅ 4 3 = 8 π 3 2\pi\cdot\dfrac43 = \dfrac{8\pi}{3} 2 π ⋅ 3 4 = 3 8 π
ตรวจด้วยสูตร: ทรงกระบอกรัศมี 2 สูง 2 มีปริมาตร 8 π 8\pi 8 π และกรวยมี 1 3 \tfrac13 3 1 ของ 8 π 8\pi 8 π = 8 π 3 \tfrac{8\pi}{3} 3 8 π → รูปตันที่ถามคือ “กรวยกลับหัว” ปริมาตร 8 π 3 \tfrac{8\pi}{3} 3 8 π ✓
V = 8 π 3 ≈ 8.378 V = \dfrac{8\pi}{3} \approx 8.378 V = 3 8 π ≈ 8.378
ตัวอย่าง M-10 · ทรงกลมและซีกทรงกลม (พิกัดทรงกลม) (ก) ปริมาตรทรงกลมรัศมี R R R (ข) จุดศูนย์กลางมวลของซีกทรงกลมตันรัศมี 1 ความหนาแน่นสม่ำเสมอ เปิดในแล็บ (ซีกทรงกลม)
(ก) V = ∫ 0 2 π ∫ 0 π ∫ 0 R ρ 2 sin φ d ρ d φ d θ = 2 π ⋅ 2 ⋅ R 3 3 = 4 π R 3 3 V = \displaystyle\int_0^{2\pi}\!\int_0^{\pi}\!\int_0^{R}\rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta = 2\pi\cdot 2\cdot\dfrac{R^3}{3} = \dfrac{4\pi R^3}{3} V = ∫ 0 2 π ∫ 0 π ∫ 0 R ρ 2 sin φ d ρ d φ d θ = 2 π ⋅ 2 ⋅ 3 R 3 = 3 4 π R 3 (R = 2 ได้ 32 π 3 = 33.51 \tfrac{32\pi}{3} = 33.51 3 32 π = 33.51 )
(ข) ซีกบน: 0 ≤ φ ≤ π / 2 0 \le \varphi \le \pi/2 0 ≤ φ ≤ π /2 · ปริมาตร = 2 π 3 = \tfrac{2\pi}{3} = 3 2 π
∭ z d V = ∫ 0 2 π ∫ 0 π / 2 ∫ 0 1 ( ρ cos φ ) ρ 2 sin φ d ρ d φ d θ = 2 π ⋅ 1 4 ⋅ 1 2 = π 4 \iiint z\,dV = \int_0^{2\pi}\!\int_0^{\pi/2}\!\int_0^1 (\rho\cos\varphi)\,\rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta = 2\pi\cdot\dfrac14\cdot\dfrac12 = \dfrac{\pi}{4} ∭ z d V = ∫ 0 2 π ∫ 0 π /2 ∫ 0 1 ( ρ cos φ ) ρ 2 sin φ d ρ d φ d θ = 2 π ⋅ 4 1 ⋅ 2 1 = 4 π
z ˉ = π / 4 2 π / 3 = 3 8 \bar z = \dfrac{\pi/4}{2\pi/3} = \dfrac38 z ˉ = 2 π /3 π /4 = 8 3 และโดยสมมาตร x ˉ = y ˉ = 0 \bar x = \bar y = 0 x ˉ = y ˉ = 0
(ก) V = 4 3 π R 3 V = \tfrac43\pi R^3 V = 3 4 π R 3 · (ข) ( x ˉ , y ˉ , z ˉ ) = ( 0 , 0 , 3 8 R ) (\bar x, \bar y, \bar z) = (0, 0, \tfrac38 R) ( x ˉ , y ˉ , z ˉ ) = ( 0 , 0 , 8 3 R ) สำหรับซีกทรงกลมรัศมี R R R (ตรงกับตารางในหนังสือกลศาสตร์วิศวกรรม)
การเปลี่ยนตัวแปรและ Jacobian ขั้นสูง
พิกัดขั้วเป็นกรณีเฉพาะของการเปลี่ยนตัวแปรทั่วไป x = x ( u , v ) , y = y ( u , v ) x = x(u, v),\ y = y(u, v) x = x ( u , v ) , y = y ( u , v ) ซึ่งแปลงช่องสี่เหลี่ยมเล็กบนระนาบ u v uv uv ให้เป็นสี่เหลี่ยมด้านขนานโค้งบนระนาบ x y xy x y ขนาดพื้นที่ใหม่เท่ากับ ∣ J ∣ |J| ∣ J ∣ คูณพื้นที่เดิม โดย J J J คือดีเทอร์มิแนนต์ของเมทริกซ์จาโคเบียน (จากบทอนุพันธ์ย่อย ) ใช้แล็บด้านล่างลากแถบ ( u 0 , v 0 ) (u_0, v_0) ( u 0 , v 0 ) และขนาดช่อง แล้วเทียบ “พื้นที่ภาพจริง” กับ ∣ det J ∣ Δ u Δ v |\det J|\Delta u\Delta v ∣ det J ∣Δ u Δ v
M10 Change of variables (2D) เปลี่ยนตัวแปร (สองมิติ)
∬ R f ( x , y ) d x d y = ∬ S f ( x ( u , v ) , y ( u , v ) ) ∣ ∂ ( x , y ) ∂ ( u , v ) ∣ d u d v , ∂ ( x , y ) ∂ ( u , v ) = ∣ x u x v y u y v ∣ \iint_R f(x, y)\,dx\,dy = \iint_S f\big(x(u,v),\,y(u,v)\big)\left|\frac{\partial(x, y)}{\partial(u, v)}\right|du\,dv,\qquad \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix} x_u & x_v\\ y_u & y_v\end{vmatrix} ∬ R f ( x , y ) d x d y = ∬ S f ( x ( u , v ) , y ( u , v ) ) ∂ ( u , v ) ∂ ( x , y ) d u d v , ∂ ( u , v ) ∂ ( x , y ) = x u y u x v y v
การแปลงต้องหนึ่งต่อหนึ่ง (อาจยกเว้นขอบ) · พิกัดขั้ว J = r J = r J = r · ∣ J ∣ |J| ∣ J ∣ ใช้ค่าสัมบูรณ์เสมอ (พื้นที่ไม่ติดลบ) · เลือกตัวแปรใหม่ให้บริเวณ S S S กลายเป็นสี่เหลี่ยมผืนผ้า
M11 Jacobians of standard coordinates จาโคเบียนของพิกัดมาตรฐาน
ขั้ว: ∣ ∂ ( x , y ) ∂ ( r , θ ) ∣ = r ทรงกระบอก: r ทรงกลม: ∣ ∂ ( x , y , z ) ∂ ( ρ , φ , θ ) ∣ = ρ 2 sin φ \text{ขั้ว: } \left|\tfrac{\partial(x,y)}{\partial(r,\theta)}\right| = r\qquad \text{ทรงกระบอก: } r\qquad \text{ทรงกลม: } \left|\tfrac{\partial(x,y,z)}{\partial(\rho,\varphi,\theta)}\right| = \rho^2\sin\varphi ขั้ว : ∂ ( r , θ ) ∂ ( x , y ) = r ทรงกระบอก : r ทรงกลม : ∂ ( ρ , φ , θ ) ∂ ( x , y , z ) = ρ 2 sin φ
ตัวคูณทั้งสามคือที่มาของ r d r d θ , r d z d r d θ , ρ 2 sin φ d ρ d φ d θ r\,dr\,d\theta,\ r\,dz\,dr\,d\theta,\ \rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta r d r d θ , r d z d r d θ , ρ 2 sin φ d ρ d φ d θ · สามตัวแปร J J J เป็นดีเทอร์มิแนนต์ 3×3
u v ระนาบ (u, v) Δu Δv ΔA′ x y ระนาบ (x, y) x = x(u, v), y = y(u, v) r_u Δu r_v Δv ΔA ≈ |det J| Δu Δv, J = ∂(x, y)/∂(u, v) จาโคเบียน (Jacobian): ช่องเล็ก Δ u × Δ v \Delta u\times\Delta v Δ u × Δ v ถูกแปลงเป็นสี่เหลี่ยมด้านขนานที่ขึงด้วย ( x u , y u ) Δ u (x_u, y_u)\Delta u ( x u , y u ) Δ u และ ( x v , y v ) Δ v (x_v, y_v)\Delta v ( x v , y v ) Δ v พื้นที่ใหม่คือ ∣ det J ∣ Δ u Δ v |\det J|\,\Delta u\,\Delta v ∣ det J ∣ Δ u Δ v
ตัวอย่าง M-11 · พื้นที่วงรี หาพื้นที่ภายในวงรี x 2 a 2 + y 2 b 2 ≤ 1 \dfrac{x^2}{a^2} + \dfrac{y^2}{b^2} \le 1 a 2 x 2 + b 2 y 2 ≤ 1
เปลี่ยนตัวแปร x = a u , y = b v x = au,\ y = bv x = a u , y = b v · วงรีกลายเป็นแผ่นวงกลมหนึ่งหน่วย u 2 + v 2 ≤ 1 u^2 + v^2 \le 1 u 2 + v 2 ≤ 1
J = ∣ a 0 0 b ∣ = a b J = \begin{vmatrix} a & 0\\ 0 & b\end{vmatrix} = ab J = a 0 0 b = ab
A = ∬ d x d y = ∬ u 2 + v 2 ≤ 1 a b d u d v = a b ⋅ π ⋅ 1 2 = π a b A = \iint dx\,dy = \iint_{u^2+v^2\le1} ab\,du\,dv = ab\cdot\pi\cdot1^2 = \pi ab A = ∬ d x d y = ∬ u 2 + v 2 ≤ 1 ab d u d v = ab ⋅ π ⋅ 1 2 = πab
ตัวอย่างเชิงตัวเลข: a = 2 , b = 1 a = 2,\ b = 1 a = 2 , b = 1 ได้ A = 2 π = 6.283 A = 2\pi = 6.283 A = 2 π = 6.283 (แล็บ Jacobian โหมด “วงรี” ให้ det J = 2 u \det J = 2u det J = 2 u เมื่อรวม x = 2 u cos v x = 2u\cos v x = 2 u cos v ในพิกัดขั้ววงรี)
A = π a b A = \pi ab A = πab
มวล จุดศูนย์กลางมวล และโมเมนต์ความเฉื่อย กลาง
แผ่นบางที่ความหนาแน่นผิว ρ ( x , y ) \rho(x, y) ρ ( x , y ) (มวลต่อพื้นที่) มีมวล M = ∬ ρ d A M = \iint\rho\,dA M = ∬ ρ d A จุดศูนย์กลางมวลคือค่าเฉลี่ยถ่วงน้ำหนักของพิกัด และโมเมนต์ความเฉื่อยวัดความ “ต้านการหมุน” ซึ่งใช้ในสมการ τ = I α \tau = I\alpha τ = I α ของเครื่องกล และ I I I ของหน้าตัดคาน (ใส่ ρ = 1 \rho = 1 ρ = 1 ) ในงานโยธา เปิดสวิตช์ “ถือว่า f f f = ความหนาแน่น” ในแล็บเพื่อให้โปรแกรมคำนวณทั้งหมดพร้อมกัน
M12 Mass and center of mass มวลและจุดศูนย์กลางมวล
M = ∬ R ρ d A , x ˉ = 1 M ∬ R x ρ d A , y ˉ = 1 M ∬ R y ρ d A M = \iint_R \rho\,dA,\qquad \bar x = \frac{1}{M}\iint_R x\,\rho\,dA,\qquad \bar y = \frac{1}{M}\iint_R y\,\rho\,dA M = ∬ R ρ d A , x ˉ = M 1 ∬ R x ρ d A , y ˉ = M 1 ∬ R y ρ d A
ρ = 1 \rho = 1 ρ = 1 : จุดศูนย์ถ่วงทางเรขาคณิต (centroid) · สามมิติ z ˉ = 1 M ∭ z ρ d V \bar z = \tfrac1M\iiint z\rho\,dV z ˉ = M 1 ∭ z ρ d V · บริเวณสมมาตรให้ x ˉ \bar x x ˉ หรือ y ˉ \bar y y ˉ อยู่ที่แกนสมมาตร
M13 Moments of inertia โมเมนต์ความเฉื่อย
I x = ∬ y 2 ρ d A , I y = ∬ x 2 ρ d A , I 0 = I x + I y = ∬ ( x 2 + y 2 ) ρ d A , I z = ∭ ( x 2 + y 2 ) ρ d V I_x = \iint y^2\rho\,dA,\qquad I_y = \iint x^2\rho\,dA,\qquad I_0 = I_x + I_y = \iint (x^2 + y^2)\rho\,dA,\qquad I_z = \iiint (x^2 + y^2)\rho\,dV I x = ∬ y 2 ρ d A , I y = ∬ x 2 ρ d A , I 0 = I x + I y = ∬ ( x 2 + y 2 ) ρ d A , I z = ∭ ( x 2 + y 2 ) ρ d V
I x I_x I x ถ่วงด้วยระยะจากแกน x x x (จึงเป็น y 2 y^2 y 2 ) · I 0 I_0 I 0 = โมเมนต์เชิงขั้ว · ทฤษฎีบทแกนขนาน I = I c + M d 2 I = I_c + Md^2 I = I c + M d 2 · รัศมีไจเรชัน k = I / M k = \sqrt{I/M} k = I / M
ตัวอย่าง M-12 · แผ่นครึ่งวงกลม แผ่นครึ่งวงกลมรัศมี 1 (y ≥ 0 y \ge 0 y ≥ 0 ) ความหนาแน่นสม่ำเสมอ ρ = 1 \rho = 1 ρ = 1 หา M , y ˉ , I x , I y M,\ \bar y,\ I_x,\ I_y M , y ˉ , I x , I y เปิดในแล็บ
ขอบเขตขั้ว: 0 ≤ r ≤ 1 , 0 ≤ θ ≤ π 0 \le r \le 1,\ 0 \le \theta \le \pi 0 ≤ r ≤ 1 , 0 ≤ θ ≤ π · M = ∫ 0 π ∫ 0 1 r d r d θ = π 2 M = \displaystyle\int_0^\pi\!\int_0^1 r\,dr\,d\theta = \dfrac{\pi}{2} M = ∫ 0 π ∫ 0 1 r d r d θ = 2 π
∬ y d A = ∫ 0 π sin θ d θ ∫ 0 1 r 2 d r = 2 ⋅ 1 3 = 2 3 \iint y\,dA = \displaystyle\int_0^\pi\sin\theta\,d\theta\int_0^1 r^2\,dr = 2\cdot\dfrac13 = \dfrac23 ∬ y d A = ∫ 0 π sin θ d θ ∫ 0 1 r 2 d r = 2 ⋅ 3 1 = 3 2 → y ˉ = 2 / 3 π / 2 = 4 3 π ≈ 0.4244 \bar y = \dfrac{2/3}{\pi/2} = \dfrac{4}{3\pi} \approx 0.4244 y ˉ = π /2 2/3 = 3 π 4 ≈ 0.4244 · x ˉ = 0 \bar x = 0 x ˉ = 0 (สมมาตร)
I x = ∬ y 2 d A = ∫ 0 π sin 2 θ d θ ∫ 0 1 r 3 d r = π 2 ⋅ 1 4 = π 8 I_x = \iint y^2\,dA = \displaystyle\int_0^\pi\sin^2\theta\,d\theta\int_0^1 r^3\,dr = \dfrac{\pi}{2}\cdot\dfrac14 = \dfrac{\pi}{8} I x = ∬ y 2 d A = ∫ 0 π sin 2 θ d θ ∫ 0 1 r 3 d r = 2 π ⋅ 4 1 = 8 π · I y = ∬ x 2 d A = π 8 I_y = \iint x^2\,dA = \dfrac{\pi}{8} I y = ∬ x 2 d A = 8 π
I 0 = π 4 I_0 = \dfrac{\pi}{4} I 0 = 4 π และเมื่อ ρ \rho ρ เป็นมวลต่อพื้นที่ M = π 2 M = \dfrac{\pi}{2} M = 2 π ได้ I 0 = 1 2 M R 2 I_0 = \tfrac12 M R^2 I 0 = 2 1 M R 2 เหมือนจานกลมเต็มดวง (ซีกก็มี I 0 I_0 I 0 อัตราส่วนเดียวกัน)
M = π 2 M = \dfrac{\pi}{2} M = 2 π · y ˉ = 4 3 π \bar y = \dfrac{4}{3\pi} y ˉ = 3 π 4 · I x = I y = π 8 I_x = I_y = \dfrac{\pi}{8} I x = I y = 8 π · I 0 = π 4 I_0 = \dfrac{\pi}{4} I 0 = 4 π
ตัวอย่าง M-13 · โมเมนต์ความเฉื่อยของกระบอกและทรงกลม (เครื่องกล) แสดงว่าทรงกระบอกตันรัศมี R R R สูง h h h มี I z = 1 2 M R 2 I_z = \tfrac12 MR^2 I z = 2 1 M R 2 และทรงกลมตันมี I = 2 5 M R 2 I = \tfrac25 MR^2 I = 5 2 M R 2 (ความหนาแน่น ρ 0 \rho_0 ρ 0 สม่ำเสมอ)
ทรงกระบอก (พิกัดทรงกระบอก): I z = ∫ 0 2 π ∫ 0 R ∫ 0 h ρ 0 r 2 ⋅ r d z d r d θ = ρ 0 ⋅ 2 π ⋅ R 4 4 ⋅ h = π ρ 0 h R 4 2 I_z = \displaystyle\int_0^{2\pi}\!\int_0^R\!\int_0^h \rho_0\,r^2\cdot r\,dz\,dr\,d\theta = \rho_0\cdot2\pi\cdot\dfrac{R^4}{4}\cdot h = \dfrac{\pi\rho_0 hR^4}{2} I z = ∫ 0 2 π ∫ 0 R ∫ 0 h ρ 0 r 2 ⋅ r d z d r d θ = ρ 0 ⋅ 2 π ⋅ 4 R 4 ⋅ h = 2 π ρ 0 h R 4
มวล M = ρ 0 π R 2 h M = \rho_0\pi R^2h M = ρ 0 π R 2 h ดังนั้น I z = 1 2 M R 2 I_z = \dfrac{1}{2}MR^2 I z = 2 1 M R 2 ✓
ทรงกลม (พิกัดทรงกลม): ( x 2 + y 2 ) = ρ 2 sin 2 φ (x^2 + y^2) = \rho^2\sin^2\varphi ( x 2 + y 2 ) = ρ 2 sin 2 φ · I = ρ 0 ∫ 0 2 π ∫ 0 π ∫ 0 R ρ 2 sin 2 φ ⋅ ρ 2 sin φ d ρ d φ d θ = ρ 0 ⋅ R 5 5 ⋅ 4 3 ⋅ 2 π = 8 π ρ 0 R 5 15 I = \rho_0\displaystyle\int_0^{2\pi}\!\int_0^\pi\!\int_0^R \rho^2\sin^2\varphi\cdot\rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta = \rho_0\cdot\dfrac{R^5}{5}\cdot\dfrac43\cdot2\pi = \dfrac{8\pi\rho_0R^5}{15} I = ρ 0 ∫ 0 2 π ∫ 0 π ∫ 0 R ρ 2 sin 2 φ ⋅ ρ 2 sin φ d ρ d φ d θ = ρ 0 ⋅ 5 R 5 ⋅ 3 4 ⋅ 2 π = 15 8 π ρ 0 R 5
มวล M = 4 π ρ 0 R 3 3 M = \dfrac{4\pi\rho_0R^3}{3} M = 3 4 π ρ 0 R 3 → I = 8 π / 15 4 π / 3 M R 2 = 2 5 M R 2 I = \dfrac{8\pi/15}{4\pi/3}MR^2 = \dfrac25MR^2 I = 4 π /3 8 π /15 M R 2 = 5 2 M R 2 ✓
กระบอกตัน I z = 1 2 M R 2 I_z = \tfrac12MR^2 I z = 2 1 M R 2 · ทรงกลมตัน I = 2 5 M R 2 I = \tfrac25MR^2 I = 5 2 M R 2 (ใช้ ∫ 0 π sin 3 φ d φ = 4 3 \int_0^\pi\sin^3\varphi\,d\varphi = \tfrac43 ∫ 0 π sin 3 φ d φ = 3 4 )
ตัวอย่าง M-14 · ประจุรวมบนแผ่นกลม (ไฟฟ้า) แผ่นตัวนำกลมรัศมี a = 5 a = 5 a = 5 cm มีความหนาแน่นประจุผิว σ ( r ) = σ 0 e − r / b \sigma(r) = \sigma_0 e^{-r/b} σ ( r ) = σ 0 e − r / b เมื่อ σ 0 = 20 \sigma_0 = 20 σ 0 = 20 nC/m² และ b = 2 b = 2 b = 2 cm จงหาประจุรวม
Q = ∬ σ d A = ∫ 0 2 π ∫ 0 a σ 0 e − r / b r d r d θ = 2 π σ 0 ∫ 0 a r e − r / b d r Q = \displaystyle\iint\sigma\,dA = \int_0^{2\pi}\!\int_0^{a}\sigma_0 e^{-r/b}\,r\,dr\,d\theta = 2\pi\sigma_0\int_0^a r e^{-r/b}dr Q = ∬ σ d A = ∫ 0 2 π ∫ 0 a σ 0 e − r / b r d r d θ = 2 π σ 0 ∫ 0 a r e − r / b d r
อินทิเกรตทีละส่วน: ∫ 0 a r e − r / b d r = b 2 [ 1 − ( 1 + a b ) e − a / b ] \displaystyle\int_0^a r e^{-r/b}dr = b^2\left[1 - \left(1 + \dfrac{a}{b}\right)e^{-a/b}\right] ∫ 0 a r e − r / b d r = b 2 [ 1 − ( 1 + b a ) e − a / b ]
แทน a / b = 2.5 a/b = 2.5 a / b = 2.5 : e − 2.5 = 0.08208 e^{-2.5} = 0.08208 e − 2.5 = 0.08208 · วงเล็บ = 1 − 3.5 ( 0.08208 ) = 0.7127 = 1 - 3.5(0.08208) = 0.7127 = 1 − 3.5 ( 0.08208 ) = 0.7127
Q = 2 π ( 20 × 10 − 9 ) ( 0.02 ) 2 ( 0.7127 ) = 3.58 × 10 − 11 Q = 2\pi(20\times10^{-9})(0.02)^2(0.7127) = 3.58\times10^{-11} Q = 2 π ( 20 × 1 0 − 9 ) ( 0.02 ) 2 ( 0.7127 ) = 3.58 × 1 0 − 11 C
Q ≈ 35.8 Q \approx 35.8 Q ≈ 35.8 pC (ถ้าประจุสม่ำเสมอที่ σ 0 \sigma_0 σ 0 ทั้งแผ่นจะได้ 157 pC ซึ่งมากกว่า 4.4 เท่า เพราะประจุจริงกระจุกใกล้ศูนย์กลาง)
นำไปใช้ในงานวิศวกรรม พื้นฐาน
สาขา อินทิกรัลหลายชั้นที่ใช้ ผลลัพธ์
ไฟฟ้า Q = ∬ σ d A Q = \iint\sigma\,dA Q = ∬ σ d A , Q = ∭ ρ v d V Q = \iiint\rho_v\,dV Q = ∭ ρ v d V ประจุรวมบนแผ่น/ในปริมาตร (ตัวอย่าง M-14)
ไฟฟ้า W = 1 2 ε 0 ∭ E 2 d V W = \tfrac12\varepsilon_0\iiint E^2\,dV W = 2 1 ε 0 ∭ E 2 d V พลังงานที่สะสมในสนามไฟฟ้า (บทไฟฟ้า )
เครื่องกล I = ∭ r ⊥ 2 ρ d V I = \iiint r_\perp^2\rho\,dV I = ∭ r ⊥ 2 ρ d V โมเมนต์ความเฉื่อยของล้อ เพลา ใบพัด (ตัวอย่าง M-13)
เครื่องกล r ˉ = 1 M ∭ r ρ d V \bar{\mathbf r} = \tfrac1M\iiint\mathbf r\rho\,dV r ˉ = M 1 ∭ r ρ d V จุดศูนย์กลางมวลของชิ้นส่วน (บทเครื่องกล )
โยธา V = ∬ ( h − z 0 ) d A V = \iint (h - z_0)\,dA V = ∬ ( h − z 0 ) d A ปริมาตรดินขุด-ถม (บทโยธา )
โยธา A = ∬ d A , I x = ∬ y 2 d A A = \iint dA,\ I_x = \iint y^2\,dA A = ∬ d A , I x = ∬ y 2 d A พื้นที่และโมเมนต์ความเฉื่อยของหน้าตัดคาน (แล็บคาน )
ทุกสาขา ∬ f d A / A \iint f\,dA / A ∬ f d A / A ค่าเฉลี่ยบนพื้นที่ เช่นอุณหภูมิหรือความเข้มรังสีเฉลี่ย
ข้อผิดพลาดที่พบบ่อย พื้นฐาน
ไม่วาดบริเวณ ก่อนเขียนขอบเขต · ส่วนใหญ่ที่ผิดเพราะอ่านขอบจากสมการโดยไม่เห็นรูป
ขอบชั้นนอกเป็นตัวแปร · ขอบชั้นนอกต้องเป็นค่าคงที่ (ผลลัพธ์สุดท้ายต้องเป็นตัวเลข)
ลืมตัวคูณ r r r ในพิกัดขั้วและทรงกระบอก และ ρ 2 sin φ \rho^2\sin\varphi ρ 2 sin φ ในพิกัดทรงกลม
สลับลำดับโดยไม่เปลี่ยนขอบ · เปลี่ยนลำดับ d y d x dy\,dx d y d x เป็น d x d y dx\,dy d x d y ต้องอ่านขอบจากรูปใหม่ทุกตัว
ใช้ φ \varphi φ กับ θ \theta θ สลับกัน ระหว่างตำราต่างเล่ม · ตรวจนิยามก่อน (ในเว็บนี้ φ \varphi φ = มุมจากแกน z z z , θ \theta θ = มุมราบ)
ได้ปริมาตรติดลบ · ตรวจว่า “บน − ล่าง” ถูกด้าน; พื้นที่และปริมาตรไม่เป็นลบ
ไปต่อที่ฟังก์ชันเวกเตอร์ ซึ่งอธิบายการเคลื่อนที่ของวัตถุในอวกาศด้วย r ( t ) \mathbf r(t) r ( t ) และเตรียมไปสู่อินทิกรัลตามเส้นและบนผิว
คำถามที่พบบ่อย อินทิกรัลสองชั้น (double integral) คืออะไร คำนวณอย่างไร? คือผลรวมของค่า f ( x , y ) f(x,y) f ( x , y ) คูณพื้นที่ช่องเล็ก Δ A \Delta A Δ A ทั่วบริเวณ R R R เมื่อ f ≥ 0 f \ge 0 f ≥ 0 จะเป็นปริมาตรใต้ผิว z = f ( x , y ) z = f(x, y) z = f ( x , y ) เหนือ R R R คำนวณด้วยทฤษฎีบทฟูบินี (Fubini) โดยทำทีละชั้น: อินทิเกรตตัวแปรในสุดก่อนโดยถือตัวแปรนอกเป็นค่าคงที่ แล้วอินทิเกรตผลลัพธ์ตามตัวแปรนอก เช่น ∫ 0 1 ∫ 0 2 ( x 2 + y 2 ) d y d x = ∫ 0 1 ( 2 x 2 + 8 3 ) d x = 10 3 \int_0^1\!\int_0^2 (x^2+y^2)\,dy\,dx = \int_0^1 (2x^2 + \tfrac83)\,dx = \tfrac{10}{3} ∫ 0 1 ∫ 0 2 ( x 2 + y 2 ) d y d x = ∫ 0 1 ( 2 x 2 + 3 8 ) d x = 3 10
เมื่อไรควรเปลี่ยนเป็นพิกัดขั้ว (polar) ในอินทิกรัลสองชั้น? เมื่อบริเวณเป็นวงกลม ครึ่งวงกลม วงแหวน หรือเซกเตอร์ หรือเมื่อ integrand มีกลุ่ม x 2 + y 2 x^2 + y^2 x 2 + y 2 แทน x = r cos θ , y = r sin θ x = r\cos\theta,\ y = r\sin\theta x = r cos θ , y = r sin θ และเปลี่ยน d A dA d A เป็น r d r d θ r\,dr\,d\theta r d r d θ อย่าลืมตัวคูณ r r r เด็ดขาด เพราะช่องเล็กในพิกัดขั้วกว้าง d r dr d r และยาว r d θ r\,d\theta r d θ ตัวอย่าง ∬ e − ( x 2 + y 2 ) d A \iint e^{-(x^2+y^2)}\,dA ∬ e − ( x 2 + y 2 ) d A ทั่วระนาบ = π \pi π คำนวณได้ง่ายในพิกัดขั้ว แต่แทบทำไม่ได้ในพิกัดฉาก
พิกัดทรงกระบอกกับพิกัดทรงกลมเลือกใช้อย่างไร? พิกัดทรงกระบอก ( r , θ , z ) (r, \theta, z) ( r , θ , z ) เหมาะกับรูปที่มีแกนสมมาตร เช่นทรงกระบอก กรวย พาราโบลอยด์ และมี d V = r d z d r d θ dV = r\,dz\,dr\,d\theta d V = r d z d r d θ ส่วนพิกัดทรงกลม ( ρ , φ , θ ) (\rho, \varphi, \theta) ( ρ , φ , θ ) เหมาะกับรูปที่สมมาตรรอบจุด เช่นทรงกลม ซีกทรงกลม เปลือกทรงกลม และมี d V = ρ 2 sin φ d ρ d φ d θ dV = \rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta d V = ρ 2 sin φ d ρ d φ d θ ข้อสังเกต: ถ้าสมการขอบเขตมี x 2 + y 2 + z 2 x^2 + y^2 + z^2 x 2 + y 2 + z 2 ให้นึกถึงทรงกลมก่อน ถ้ามี x 2 + y 2 x^2 + y^2 x 2 + y 2 อย่างเดียวให้นึกถึงทรงกระบอก
สลับลำดับการอินทิเกรต (change order of integration) ทำอย่างไร? วาดบริเวณก่อน แล้วอธิบายบริเวณเดิมด้วยลำดับใหม่ เช่น ∫ 0 1 ∫ x 1 e y 2 d y d x \int_0^1\!\int_x^1 e^{y^2}\,dy\,dx ∫ 0 1 ∫ x 1 e y 2 d y d x อินทิเกรต e y 2 e^{y^2} e y 2 ตาม y y y ไม่ได้ แต่บริเวณคือสามเหลี่ยม 0 ≤ x ≤ y ≤ 1 0 \le x \le y \le 1 0 ≤ x ≤ y ≤ 1 สลับเป็น ∫ 0 1 ∫ 0 y e y 2 d x d y = ∫ 0 1 y e y 2 d y = e − 1 2 \int_0^1\!\int_0^y e^{y^2}\,dx\,dy = \int_0^1 y e^{y^2}\,dy = \tfrac{e-1}{2} ∫ 0 1 ∫ 0 y e y 2 d x d y = ∫ 0 1 y e y 2 d y = 2 e − 1 ใช้สลับเมื่อลำดับเดิมอินทิเกรตไม่ได้หรือขอบเขตซับซ้อน
จาโคเบียน (Jacobian) คืออะไร ทำไมต้องคูณเวลาเปลี่ยนตัวแปร? จาโคเบียนคือดีเทอร์มิแนนต์ของเมทริกซ์อนุพันธ์ย่อยของการแปลง ∂ ( x , y ) / ∂ ( u , v ) \partial(x,y)/\partial(u,v) ∂ ( x , y ) / ∂ ( u , v ) มันบอกว่าช่องเล็กบนระนาบ u v uv uv ถูกขยายหรือหดเป็นกี่เท่าเมื่อแปลงไปเป็นระนาบ x y xy x y จึงต้องคูณ ∣ J ∣ |J| ∣ J ∣ เพื่อรักษาพื้นที่ที่ถูกต้อง ในพิกัดขั้ว J = r J = r J = r ในพิกัดทรงกระบอก J = r J = r J = r ในพิกัดทรงกลม J = ρ 2 sin φ J = \rho^2\sin\varphi J = ρ 2 sin φ สำหรับวงรี x = a u , y = b v x = au,\ y = bv x = a u , y = b v ได้ J = a b J = ab J = ab จึงได้พื้นที่ π a b \pi ab πab
มด
อาจารย์มด · ธีรวุฒิ แสวงบุญ อาจารย์สาขาวิชาวิศวกรรมไฟฟ้า มหาวิทยาลัยราชภัฏธนบุรี สมุทรปราการ · ดูโปรไฟล์