球面度 立体角 ,常用字母Ω 表示,是一个物体对特定点的三维空间的角度,是平面角 在三维空间中的类比。它描述的是站在某一点的观察者测量到的物体大小的尺度。例如,对于一个特定的观察点,一个在该观察点附近的小物体有可能和一个远处的大物体有着相同的立体角。
锥体 的立体角大小定义为,以锥体的顶点为球心作球面 ,该锥体在球表面截取的面积与球半径平方之比,单位为球面度 (sr)。[ 1] 一個球體即是4π,半球 即是2π。
立體角 d Ω = sin θ d θ d φ {\displaystyle d\Omega =\sin \theta \,d\theta \,d\varphi } 以觀測點為球心,構造一個單位球面;任意物體投影到該單位球面上的投影面積,即為該物體相對於該觀測點的立體角。
因此,立體角是單位球面上的一塊面積,這和“平面角是單位圆上的一段弧长”類似。
在球坐标系 中,任意球面的極小面積為:
d A = ( r sin θ d φ ) ( r d θ ) = r 2 ( sin θ d θ d φ ) {\displaystyle dA=(r\sin \theta \,d\varphi )(rd\theta )=r^{2}(\sin \theta \,d\theta \,d\varphi )} 因此,極小立體角(單位球面上的極小面積)為:
d Ω = d A r 2 = sin θ d θ d φ {\displaystyle d\Omega ={\frac {dA}{r^{2}}}=\sin \theta \,d\theta \,d\varphi } 所以,立體角是投影面积与球半径平方值的比,这和“平面角是圆的弧长与半径的比”类似。 對極小立體角做曲面积分 即可得立體角:
Ω = ∬ S d Ω = ∬ S sin θ d θ d φ {\displaystyle \Omega =\iint _{S}d\Omega =\iint _{S}\sin \theta \,d\theta \,d\varphi } 任意定向曲面 S → {\displaystyle {\vec {S}}} 相對於某一個點 P {\displaystyle P} 的立體角,即為該曲面投影到以 P {\displaystyle P} 為球心的單位球面上的面積。 令 r → {\displaystyle {\vec {r}}} 為該單位球面上以 P {\displaystyle P} 為原點的極小面積的位置向量 ,可以得到以下公式:
Ω = ∬ S d A r 2 = ∬ S r → ⋅ d S → | r → | r 2 = ∬ S r → ⋅ d S → r 3 {\displaystyle \Omega =\iint _{S}{\frac {dA}{r^{2}}}=\iint _{S}{\frac {{\vec {r}}\cdot {\textrm {d}}{\vec {S}}}{\left|{\vec {r}}\right|\,r^{2}}}=\iint _{S}{\frac {{\vec {r}}\cdot {\textrm {d}}{\vec {S}}}{r^{3}}}} 立体角的国际制单位是球面度 (steradian,sr)。立体角有一个非国际制单位平方度 ,1 sr = (180/π)2 square degree 。
一个完整的球面对于球内任意一点的立体角为4π sr (对于球外任意一点的立体角为0 sr ):
∯ {\displaystyle \oiint } S {\displaystyle {\scriptstyle S}} sin θ d θ d φ {\displaystyle \ \sin \theta \,d\theta \,d\varphi } = ∫ 0 π sin θ d θ ∫ 0 2 π d φ = [ − cos θ ] 0 π ( 2 π ) = 4 π {\displaystyle \ =\int _{0}^{\pi }\sin \theta d\theta \int _{0}^{2\pi }d\varphi =[-\cos \theta ]_{0}^{\pi }(2\pi )=4\pi } 這個定理對所有封閉曲面皆成立,它也是高斯定律 的主要依據[ 2] 。
立體角在物理上有相當多的應用:
球體中圓錐(1)與球冠 (2)的切面。圖中 θ = A /2,r = 1。 顶角为2 θ {\displaystyle \theta } 的圆锥的立体角为一个单位球的球冠。
Ω = 2 π ( 1 − cos θ ) . {\displaystyle \Omega =2\pi \left(1-\cos {\theta }\right).\,\!} (上面结果由下式得到)
∫ 0 2 π ∫ 0 θ sin θ ′ d θ ′ d ϕ = 2 π ∫ 0 θ sin θ ′ d θ ′ = 2 π [ − cos θ ′ ] 0 θ = 2 π ( 1 − cos θ ) . {\displaystyle \int _{0}^{2\pi }\int _{0}^{\theta }\sin \theta '\ d\theta '\ d\phi =2\pi \int _{0}^{\theta }\sin \theta '\ d\theta '=2\pi \left[-\cos \theta '\right]_{0}^{\theta }\ =2\pi \left(1-\cos \theta \right).} 应该注意阿基米德 在2200年前不用微积分证明了球冠的表面积与半径为球冠边沿到球冠最低点的距离的圆的面积相等。球冠边沿到球冠最低点的距离为
2 r sin ( θ 2 ) . {\displaystyle 2r\sin \left({\frac {\theta }{2}}\right).\,} 显然,在单位圆中球冠立体角为
Ω = 4 π sin 2 ( θ 2 ) = 2 π ( 1 − cos θ ) . {\displaystyle \Omega =4\pi \sin ^{2}\left({\frac {\theta }{2}}\right)=2\pi \left(1-\cos {\theta }\right).\,\!} 当 θ = π/2,球冠变为有着立体角 2π的半球。
当 θ = π,立體角涵蓋整個球體,球冠变为有着立体角 4π的球,我們將4π稱為全方位立體角。
对于任意一个四面体 OABC,其中O,A,B,C分别为四面体的四个顶点。下面给出一个公式,计算从O点观察三角形ABC的立体角Ω的方便简单的公式。令α=∠BOC,β=∠AOC,γ=∠AOB(均为各自平面内两条直线的夹角,可以采用平面三角形 的余弦公式 计算求得), s = 1 2 ( α + β + γ ) . {\displaystyle s={\frac {1}{2}}(\alpha +\beta +\gamma ).} 有(参见L' Huilier' )
tan ( Ω 4 ) = tan ( s 2 ) tan ( s − α 2 ) tan ( s − β 2 ) tan ( s − γ 2 ) {\displaystyle \tan \left({\frac {\Omega }{4}}\right)={\sqrt {\tan \left({\frac {s}{2}}\right)\tan \left({\frac {s-\alpha }{2}}\right)\tan \left({\frac {s-\beta }{2}}\right)\tan \left({\frac {s-\gamma }{2}}\right)}}} 附上相应的Fortran 程序:
program solid_angle REAL X ( 4 ), Y ( 4 ), Z ( 4 ) !存放四个顶点的三维坐标,依次为A,B,C,O(程序中以P点代替,为观测点) X ( 1 ) = 1 X ( 2 ) = 0 X ( 3 ) = 0 X ( 4 ) = 0 Y ( 1 ) = 0 Y ( 2 ) = 1 Y ( 3 ) = 0 Y ( 4 ) = 0 Z ( 1 ) = 0 Z ( 2 ) = 0 Z ( 3 ) = 1 Z ( 4 ) = 0 !以单位1球体,在第一象限的八分之一球所含的四面体为例。程序可适用于任意四面体。 !分别计算6条棱的长度的平方 PA = ( X ( 4 ) - X ( 1 )) ** 2 + ( Y ( 4 ) - Y ( 1 )) ** 2 + ( Z ( 4 ) - Z ( 1 )) ** 2 PB = ( X ( 4 ) - X ( 2 )) ** 2 + ( Y ( 4 ) - Y ( 2 )) ** 2 + ( Z ( 4 ) - Z ( 2 )) ** 2 PC = ( X ( 4 ) - X ( 3 )) ** 2 + ( Y ( 4 ) - Y ( 3 )) ** 2 + ( Z ( 4 ) - Z ( 3 )) ** 2 AB = ( X ( 1 ) - X ( 2 )) ** 2 + ( Y ( 1 ) - Y ( 2 )) ** 2 + ( Z ( 1 ) - Z ( 2 )) ** 2 BC = ( X ( 3 ) - X ( 2 )) ** 2 + ( Y ( 3 ) - Y ( 2 )) ** 2 + ( Z ( 3 ) - Z ( 2 )) ** 2 AC = ( X ( 1 ) - X ( 3 )) ** 2 + ( Y ( 1 ) - Y ( 3 )) ** 2 + ( Z ( 1 ) - Z ( 3 )) ** 2 !应用余弦公式计算三个夹角 COSALPHA = ( PB + PC - BC ) / 2.0 / SQRT ( PB * PC ) COSBETA = ( PA + PC - AC ) / 2.0 / SQRT ( PA * PC ) COSGAMA = ( PB + PA - AB ) / 2.0 / SQRT ( PB * PA ) ALPHA = acos (( PB + PC - BC ) / 2.0 / SQRT ( PB * PC )) BETA = acos (( PA + PC - AC ) / 2.0 / SQRT ( PA * PC )) GAMA = acos (( PB + PA - AB ) / 2.0 / SQRT ( PB * PA )) s = ( alpha + beta + gama ) / 2. omiga = atan ( sqrt ( tan ( s / 2. ) * tan ( s / 2. - alpha / 2. ) * tan ( s / 2. - beta / 2. ) * tan ( s / 2. - gama / 2. ))) * 4. write ( * , * ) omiga * 18 0. / 3.14159 end
线性(平动)的量 角度(转动)的量 量纲 — L L2 量纲 — — — T 时间 : t s 位移积分 : A m s T 时间 : t s — 距离 : d , 位矢 : r , s , x , 位移 m 面积 : A m2 — 角度 : θ , 角移 : θ rad 立體角 : Ω rad2 , sr T−1 頻率 : f s−1 , Hz 速率 : v , 速度 : v m s−1 面積速率 : ν m2 s−1 T−1 頻率 : f s−1 , Hz 角速率 : ω , 角速度 : ω rad s−1 T−2 加速度 : a m s−2 T−2 角加速度 : α rad s−2 T−3 加加速度 : j m s−3 T−3 角加加速度 : ζ rad s−3 M 质量 : m kg ML2 轉動慣量 : I kg m2 MT−1 动量 : p , 冲量 : J kg m s−1 , N s 作用量 : 𝒮 , actergy : ℵ kg m2 s−1 , J s ML2 T−1 角动量 : L , 角衝量 : ι kg m2 s−1 作用量 : 𝒮 , actergy : ℵ kg m2 s−1 , J s MT−2 力 : F , 重量 : F g kg m s−2 , N 能量 : E , 功 : W kg m2 s−2 , J ML2 T−2 力矩 : τ , moment : M kg m2 s−2 , N m 能量 : E , 功 : W kg m2 s−2 , J MT−3 加力 : Y kg m s−3 , N s−1 功率 : P kg m2 s−3 , W ML2 T−3 rotatum : P kg m2 s−3 , N m s−1 功率 : P kg m2 s−3 , W