學習目標與模型假設
將體密度或面密度轉化為質量和力矩,計算加權質心,並區分標量曲面面積與向量通量。以下模型都是假設的數學例子。密度非負且可積;質心公式還要求總質量為正。
對具有體密度 ρ \rho ρ 的立體 E E E ,質量為 M = ∭ E ρ d V M=\iiint_E\rho\,dV M = ∭ E ρ d V 。質心的坐標為 x ‾ = M − 1 ∭ E x ρ d V \overline x=M^{-1}\iiint_E x\rho\,dV x = M − 1 ∭ E x ρ d V ,其他坐標類似。繞 z z z 軸的轉動慣量為 I z = ∭ E ( x 2 + y 2 ) ρ d V I_z=\iiint_E(x^2+y^2)\rho\,dV I z = ∭ E ( x 2 + y 2 ) ρ d V ,因為 x 2 + y 2 x^2+y^2 x 2 + y 2 是到該軸的垂直距離的平方。若密度單位為質量/體積,則這些量的單位分別為質量、長度以及質量乘長度的平方。
對平面薄片,將 d V dV d V 換為 d A dA d A 並使用面積密度;對曲面,使用面密度乘 d S dS d S 。帶符號權重可以定義代數矩,但不自動構成物理上的非負質量密度;總權重為零時,不能用作質心的分母。
1. 加權三角形的中心不一定是均勻質心
設 L > 0 L>0 L > 0 、D = { x ≥ 0 , y ≥ 0 , x + y ≤ L } D=\{x\ge0,y\ge0,x+y\le L\} D = { x ≥ 0 , y ≥ 0 , x + y ≤ L } ,薄片密度為 ρ ( x , y ) = κ ( x + y ) \rho(x,y)=\kappa(x+y) ρ ( x , y ) = κ ( x + y ) ,其中 κ > 0 \kappa>0 κ > 0 。若坐標有長度單位,κ \kappa κ 的單位為質量/長度的三次方,從而 ρ \rho ρ 是面積密度。由連續簡單區域積分,
M = κ ∫ 0 L ∫ 0 L − x ( x + y ) d y d x = κ L 3 3 . M=\kappa\int_0^L\int_0^{L-x}(x+y)\,dy\,dx=\frac{\kappa L^3}{3}. M = κ ∫ 0 L ∫ 0 L − x ( x + y ) d y d x = 3 κ L 3 .
x x x 坐標的一階矩為
∬ D x ρ ( x , y ) d A = κ ∫ 0 L [ x 2 ( L − x ) + x ( L − x ) 2 2 ] d x = κ L 4 8 . \iint_Dx\rho(x,y)\,dA
=\kappa\int_0^L\left[x^2(L-x)+\frac{x(L-x)^2}{2}\right]dx=\frac{\kappa L^4}{8}. ∬ D x ρ ( x , y ) d A = κ ∫ 0 L [ x 2 ( L − x ) + 2 x ( L − x ) 2 ] d x = 8 κ L 4 .
交換坐標的對稱性給出相同的 y y y 矩。因此 x ‾ = y ‾ = 3 L / 8 \overline x=\overline y=3L/8 x = y = 3 L /8 。密度離原點越遠越大,使中心偏離均勻密度的位置 ( L / 3 , L / 3 ) (L/3,L/3) ( L /3 , L /3 ) 。對稱性給出坐標相等,但不給出它們的數值。
結果滿足 x ‾ , y ‾ ≥ 0 \overline x,\overline y\ge0 x , y ≥ 0 和 x ‾ + y ‾ = 3 L / 4 ≤ L \overline x+\overline y=3L/4\le L x + y = 3 L /4 ≤ L 。更一般地,三角形不等式在非負加權下仍成立:對 x ≥ 0 , y ≥ 0 , x + y ≤ L x\ge0,y\ge0,x+y\le L x ≥ 0 , y ≥ 0 , x + y ≤ L 乘以 ρ / M \rho/M ρ / M 積分即可證明相同的質心不等式。這是計算出的中心位於三角形內的結構檢查,獨立於算術計算。
2. 在利用徑向對稱性前推導球坐標縮放
使用坐標
T ( r , ϕ , θ ) = ( r sin ϕ cos θ , r sin ϕ sin θ , r cos ϕ ) , T(r,\phi,\theta)=(r\sin\phi\cos\theta,r\sin\phi\sin\theta,r\cos\phi), T ( r , ϕ , θ ) = ( r sin ϕ cos θ , r sin ϕ sin θ , r cos ϕ ) ,
其中 ϕ \phi ϕ 從正 z z z 軸量起,0 ≤ ϕ ≤ π 0\le\phi\le\pi 0 ≤ ϕ ≤ π ,θ \theta θ 是一整周的方位角。導數向量的長度為 1 , r , r sin ϕ 1,r,r\sin\phi 1 , r , r sin ϕ ,它們兩兩垂直,在內部具有相應的正定向。其標量三重積為 r 2 sin ϕ r^2\sin\phi r 2 sin ϕ ,得到
d V = r 2 sin ϕ d r d ϕ d θ . dV=r^2\sin\phi\,dr\,d\phi\,d\theta. d V = r 2 sin ϕ d r d ϕ d θ .
明確核算行列式:按 ( r , ϕ , θ ) (r,\phi,\theta) ( r , ϕ , θ ) 次序,導數列為
T r = ( sin ϕ cos θ , sin ϕ sin θ , cos ϕ ) , T_r=(\sin\phi\cos\theta,\sin\phi\sin\theta,\cos\phi), T r = ( sin ϕ cos θ , sin ϕ sin θ , cos ϕ ) ,
T ϕ = ( r cos ϕ cos θ , r cos ϕ sin θ , − r sin ϕ ) , T θ = ( − r sin ϕ sin θ , r sin ϕ cos θ , 0 ) . T_\phi=(r\cos\phi\cos\theta,r\cos\phi\sin\theta,-r\sin\phi),\qquad
T_\theta=(-r\sin\phi\sin\theta,r\sin\phi\cos\theta,0). T ϕ = ( r cos ϕ cos θ , r cos ϕ sin θ , − r sin ϕ ) , T θ = ( − r sin ϕ sin θ , r sin ϕ cos θ , 0 ) .
當 r > 0 r>0 r > 0 、0 < ϕ < π 0<\phi<\pi 0 < ϕ < π 時,其混合積為 r 2 sin ϕ > 0 r^2\sin\phi>0 r 2 sin ϕ > 0 。在開參數域 ( 0 , ∞ ) × ( 0 , π ) × ( − π , π ) (0,\infty)\times(0,\pi)\times(-\pi,\pi) ( 0 , ∞ ) × ( 0 , π ) × ( − π , π ) 上,映射對其切縫後的像是一對一的,光滑逆為 r = x 2 + y 2 + z 2 r=\sqrt{x^2+y^2+z^2} r = x 2 + y 2 + z 2 、ϕ = arccos ( z / r ) \phi=\arccos(z/r) ϕ = arccos ( z / r ) 及通常的切縫平面方位角。因此,已證明的有限維換元定理適用於該域內每個緊的截斷參數長方體。
還須獨立於待證明公式估計遺漏的目標區域。半徑 ε \varepsilon ε 的球包含在體積 8 ε 3 8\varepsilon^3 8 ε 3 的立方體中。在半徑 R R R 的球內,與軸的夾角至多為 δ < π / 2 \delta<\pi/2 δ < π /2 的兩個極帽包含在 ∣ x ∣ , ∣ y ∣ ≤ R sin δ |x|,|y|\le R\sin\delta ∣ x ∣ , ∣ y ∣ ≤ R sin δ 、∣ z ∣ ≤ R |z|\le R ∣ z ∣ ≤ R 的長方體中,體積為 8 R 3 sin 2 δ 8R^3\sin^2\delta 8 R 3 sin 2 δ 。總角寬 2 η 2\eta 2 η 、0 < η < π / 2 0<\eta<\pi/2 0 < η < π /2 的方位角接縫包含在 ∣ x ∣ , ∣ z ∣ ≤ R |x|,|z|\le R ∣ x ∣ , ∣ z ∣ ≤ R 、∣ y ∣ ≤ R sin η |y|\le R\sin\eta ∣ y ∣ ≤ R sin η 中,體積為 8 R 3 sin η 8R^3\sin\eta 8 R 3 sin η 。三種上界均隨截斷參數趨於零。乘上連續目標被積函數的上界,即控制其遺漏積分;下文的參數積分估計則控制另一側。這完成了非奇異區域逼近的論證,沒有用球座標體積公式來證明該公式自身。
原點、極點和角度接縫是坐標邊界的例外。對球上有界連續被積函數,在無奇點的截斷參數區域上計算並取極限;被省略的參數積分界趨於零。對於有界被積函數,零附近的徑向因子由常數乘 ∫ 0 ε r 2 d r = ε 3 / 3 \int_0^\varepsilon r^2dr=\varepsilon^3/3 ∫ 0 ε r 2 d r = ε 3 /3 控制,極冠因子為 ∫ 0 ε sin ϕ d ϕ = 1 − cos ε → 0 \int_0^\varepsilon\sin\phi\,d\phi=1-\cos\varepsilon\to0 ∫ 0 ε sin ϕ d ϕ = 1 − cos ε → 0 ,收縮的方位角接縫則具有趨零的角寬。該公式不允許忽略這些位置上的任意無界奇異性。
對半徑為 R > 0 R>0 R > 0 、密度為 ρ = k r 2 \rho=k r^2 ρ = k r 2 且 k > 0 k>0 k > 0 的球,
M = k ∫ 0 R r 4 d r ∫ 0 π sin ϕ d ϕ ∫ 0 2 π d θ = 4 π k R 5 5 . M=k\int_0^Rr^4dr\int_0^\pi\sin\phi\,d\phi\int_0^{2\pi}d\theta=\frac{4\pi kR^5}{5}. M = k ∫ 0 R r 4 d r ∫ 0 π sin ϕ d ϕ ∫ 0 2 π d θ = 5 4 π k R 5 .
反射對稱使每個一階矩為零,所以中心為原點。對於轉動慣量,到 z z z 軸的垂直距離平方為 r 2 sin 2 ϕ r^2\sin^2\phi r 2 sin 2 ϕ ,不是完整的 r 2 r^2 r 2 。因此
I z = k ∫ 0 R r 6 d r ∫ 0 π sin 3 ϕ d ϕ ∫ 0 2 π d θ = 8 π k R 7 21 , I_z=k\int_0^Rr^6dr\int_0^\pi\sin^3\phi\,d\phi\int_0^{2\pi}d\theta
=\frac{8\pi kR^7}{21}, I z = k ∫ 0 R r 6 d r ∫ 0 π sin 3 ϕ d ϕ ∫ 0 2 π d θ = 21 8 π k R 7 ,
其中使用了 ∫ 0 π sin 3 ϕ d ϕ = 4 / 3 \int_0^\pi\sin^3\phi\,d\phi=4/3 ∫ 0 π sin 3 ϕ d ϕ = 4/3 。比率 I z / ( M R 2 ) = 10 / 21 I_z/(MR^2)=10/21 I z / ( M R 2 ) = 10/21 無量綱。這裡 k k k 的單位為質量/長度的五次方,因為它乘以 r 2 r^2 r 2 後必須得到體密度;若把它當作普通的常數密度,單位將出錯。
3. 曲面面積包含斜率,不只是投影面積
對平面區域上的圖像 z = h ( x , y ) z=h(x,y) z = h ( x , y ) ,若 h ∈ C 1 h\in C^1 h ∈ C 1 ,參數化 r ( x , y ) = ( x , y , h ( x , y ) ) r(x,y)=(x,y,h(x,y)) r ( x , y ) = ( x , y , h ( x , y )) 的切向量叉積為 ( − h x , − h y , 1 ) (-h_x,-h_y,1) ( − h x , − h y , 1 ) 。其標量面積元為
d S = 1 + h x 2 + h y 2 d x d y . dS=\sqrt{1+h_x^2+h_y^2}\,dx\,dy. d S = 1 + h x 2 + h y 2 d x d y .
因此常數面密度 σ \sigma σ 對應的投影坐標質量密度為 σ 1 + h x 2 + h y 2 \sigma\sqrt{1+h_x^2+h_y^2} σ 1 + h x 2 + h y 2 ,一般不是常數 σ \sigma σ 。真實曲面上均勻和投影上均勻是不同模型。
對半徑為 R R R 的圓盤上方的拋物面圖像 h = x 2 + y 2 h=x^2+y^2 h = x 2 + y 2 ,面積為
A = 2 π ∫ 0 R r 1 + 4 r 2 d r = π 6 [ ( 1 + 4 R 2 ) 3 / 2 − 1 ] . A=2\pi\int_0^Rr\sqrt{1+4r^2}\,dr
=\frac\pi6[(1+4R^2)^{3/2}-1]. A = 2 π ∫ 0 R r 1 + 4 r 2 d r = 6 π [( 1 + 4 R 2 ) 3/2 − 1 ] .
相同圖像上常向量場 ( 0 , 0 , 1 ) (0,0,1) ( 0 , 0 , 1 ) 的向上通量則為 ∬ D 1 d A = π R 2 \iint_D1\,dA=\pi R^2 ∬ D 1 d A = π R 2 ,因為該場與定向面積向量的點積為一。曲面面積累積大小,通量累積帶符號法向分量。兩者的差異來自幾何,而不是積分中少了常數。
4. 通過真實極限處理奇異頂點
圓錐圖像 h ( x , y ) = a x 2 + y 2 h(x,y)=a\sqrt{x^2+y^2} h ( x , y ) = a x 2 + y 2 ,其中 a > 0 a>0 a > 0 ,在原點不可微。遠離原點時 h x 2 + h y 2 = a 2 h_x^2+h_y^2=a^2 h x 2 + h y 2 = a 2 。在去心圓環 ε ≤ r ≤ R \varepsilon\le r\le R ε ≤ r ≤ R 上,正則曲面面積為 π ( R 2 − ε 2 ) 1 + a 2 \pi(R^2-\varepsilon^2)\sqrt{1+a^2} π ( R 2 − ε 2 ) 1 + a 2 。令 ε ↓ 0 \varepsilon\downarrow0 ε ↓ 0 ,得到側面積 π R 2 1 + a 2 \pi R^2\sqrt{1+a^2} π R 2 1 + a 2 。
這個結論不能僅靠“一點的面積為零”。半徑 δ < ε \delta<\varepsilon δ < ε 之間任何更小圓環的貢獻都有上界 π ε 2 1 + a 2 \pi\varepsilon^2\sqrt{1+a^2} π ε 2 1 + a 2 。因此 Cauchy 尾部判據給出有限的反常曲面面積,且遺漏貢獻的極限由同一量控制並趨於零。另一種無界面積因子可能使該極限論證失效,正如無界體積被積函數在去心點處可能失效。
5. 對稱性、界和量綱是獨立檢查
若 E E E 及其密度在 x ↦ − x x\mapsto-x x ↦ − x 下不變,則被積函數 x ρ x\rho x ρ 變號,而區域和體積元不變,迫使積分為零。若密度破壞該對稱性,僅區域對稱還不夠。
對半徑為 R R R 的球內任意非負質量,0 ≤ x 2 + y 2 ≤ R 2 0\le x^2+y^2\le R^2 0 ≤ x 2 + y 2 ≤ R 2 蘊含 0 ≤ I z ≤ M R 2 0\le I_z\le MR^2 0 ≤ I z ≤ M R 2 。若算出的轉動慣量在此區間外,無論如何作符號化簡都與模型矛盾。質心是加權平均:整個密度乘一個正常數,會使質量和轉動慣量乘同一因子,但不改變質心。單位分析、正性、對稱性和縮放規律應分別檢查;符合其中一項不證明所有算術都正確。
分層應用與建模練習
加權中心。 重算三角形的質量和一階矩,再只用非負加權證明所得中心必滿足三角形不等式。
徑向轉動慣量。 對球密度 k r 2 kr^2 k r 2 ,推導 I z I_z I z 中的 r r r 因子及角度因子,解釋為甚麼把 r 2 r^2 r 2 當作到軸的距離平方是錯誤的。
坐標約定。 指出球坐標公式中的極角是哪一個,並說明不同角度約定如何改變所寫正弦因子和坐標範圍。
曲面與投影。 對圖像上的均勻面密度,指出投影上的等效密度。比較拋物面的標量面積與恆定竪直場的通量。
奇點檢查。 解釋為甚麼圓錐頂點不使算出的極限面積失效,並寫出僅憑“忽略該點”而缺失的估計。
積分界。 對支撐在半徑為 R R R 的球內的非負密度,證明 0 ≤ I z ≤ M R 2 0\le I_z\le MR^2 0 ≤ I z ≤ M R 2 。說明徑向例子滿足該界,並在保持普通密度不變時描述所有長度縮放對質量、質心和轉動慣量的影響。
完整解答
質量為 κ L 3 / 3 \kappa L^3/3 κ L 3 /3 ,每個一階矩為 κ L 4 / 8 \kappa L^4/8 κ L 4 /8 ,故中心為 ( 3 L / 8 , 3 L / 8 ) (3L/8,3L/8) ( 3 L /8 , 3 L /8 ) 。每個點滿足 x , y ≥ 0 x,y\ge0 x , y ≥ 0 和 x + y ≤ L x+y\le L x + y ≤ L ;將各不等式乘 ρ / M ≥ 0 \rho/M\ge0 ρ / M ≥ 0 後積分,其總積分為一。所得界為 x ‾ , y ‾ ≥ 0 \overline x,\overline y\ge0 x , y ≥ 0 和 x ‾ + y ‾ ≤ L \overline x+\overline y\le L x + y ≤ L ,計算得到的坐標對滿足這些界。這裡不要求密度恆定。
體積元貢獻 r 2 sin ϕ r^2\sin\phi r 2 sin ϕ ,密度貢獻 k r 2 kr^2 k r 2 ,到軸距離平方貢獻 r 2 sin 2 ϕ r^2\sin^2\phi r 2 sin 2 ϕ ,乘積為 k r 6 sin 3 ϕ kr^6\sin^3\phi k r 6 sin 3 ϕ 。積分分別給出 R 7 / 7 R^7/7 R 7 /7 、4 / 3 4/3 4/3 和 2 π 2\pi 2 π ,故得 8 π k R 7 / 21 8\pi kR^7/21 8 π k R 7 /21 。完整的 r 2 r^2 r 2 是到原點而非到 z z z 軸的距離平方;其中包括不貢獻於到該軸距離的竪直位移。
此處 ϕ ∈ [ 0 , π ] \phi\in[0,\pi] ϕ ∈ [ 0 , π ] 從正 z z z 軸向下量取,θ \theta θ 是完整方位角。若仰角從赤道平面量取,其範圍為 [ − π / 2 , π / 2 ] [-\pi/2,\pi/2] [ − π /2 , π /2 ] ,體積因子使用仰角的餘弦。若重新命名角度卻不相應改變其幾何定義與範圍,可能把三角因子放在錯誤變量上。
投影密度為 σ 1 + h x 2 + h y 2 \sigma\sqrt{1+h_x^2+h_y^2} σ 1 + h x 2 + h y 2 。拋物面面積為 π [ ( 1 + 4 R 2 ) 3 / 2 − 1 ] / 6 \pi[(1+4R^2)^{3/2}-1]/6 π [( 1 + 4 R 2 ) 3/2 − 1 ] /6 ,向上竪直通量為 π R 2 \pi R^2 π R 2 。當 R > 0 R>0 R > 0 時,除中心外面積因子大於一,故標量面積大於投影面積。通量使用定向叉積與場的點積,而非叉積的大小。
曲面在避開頂點的每個圓環上正則,其精確面積公式有有限極限。遺漏面積的貢獻由常數乘 ε 2 \varepsilon^2 ε 2 控制,並趨於零。若沒有控制奇點鄰域的估計,僅憑該點大小為零不能排除無界的累積貢獻。
將 0 ≤ x 2 + y 2 ≤ R 2 0\le x^2+y^2\le R^2 0 ≤ x 2 + y 2 ≤ R 2 乘非負密度並積分得到這些界。徑向例子的比率 10 / 21 10/21 10/21 嚴格介於零與一之間。若所有長度按 s > 0 s>0 s > 0 縮放,且保持所搬運的普通質量/體積密度分布不額外乘因子,則體積和質量按 s 3 s^3 s 3 縮放,質心坐標按 s s s 縮放,轉動慣量按 s 5 s^5 s 5 縮放。若密度規律本身隨縮放改變,必須另行計入其乘子;不檢查模型就不能既固定有量綱系數又假設數值密度分布不變。對規律 ρ = k r 2 \rho=kr^2 ρ = k r 2 ,固定 k k k 時半徑按 s s s 縮放也會使密度乘 s 2 s^2 s 2 ,故質量按 s 5 s^5 s 5 、轉動慣量按 s 7 s^7 s 7 縮放,與顯式公式一致。
總結與術語
質心 是按質量加權的位置,而不只是幾何中點。轉動慣量 以到軸的垂直距離平方加權質量。面密度 乘真實曲面面積,而通量 使用有定向的法向分量。奇異坐標極限 需要鄰域的定量估計。量綱、對稱性、正性和縮放規律是每個積分應用的重要檢查。