立体角.docVIP

立体角 一个锥面所围成的空间部分称为“立体角”。 立体角是以锥的顶点为心，半径为1的球面被锥面所截得的面积来度量的，度量单位称为“立体弧度”。和平面角的定义类似。在平面上我们定义一段弧微分S与其矢量半径r的比值为其对应的圆心角记作dθ ds/r；所以整个圆周对应的圆心角就是2π；与此类似，定义立体角为曲面上面积微元ds与其矢量半径的二次方的比值为此面微元对应的立体角记作dΩ ds/r^2；由此可得，闭合球面的立体角都是4π。 球面度 （steradian，符号∶sr）是 立体角 的 国际单位 。它可算是三维的 弧度 。其英文字是希腊语「立体」（stereos）和弧度（radian）的混合。 例如一个球的球面度是4π，一个椭球的球面度是4π。）这正像平面角是圆的弧长与半径的比。 立体角的国际制单位是球面度（steradian）。 更严密的，立体角是面S对点P的面积分： [编辑] 常见几何体的立体角 [编辑] 圆锥，球冠 Section of cone 1 and spherical cap 2 inside a sphere. In this figure θ? ?a/2 and r? ?1. 顶角为2θ的圆锥的立体角为一个单位球的球冠。 （上面结果由下式得到，参见 surface element in spherical polars） 应该注意阿基米德在2200年前不用微积分证明了球冠的表面积与半径为球冠边沿到球冠最低点的距离的圆的面积相等。球冠边沿到球冠最低点的距离为 显然，在单位圆中球冠立体角为 相关的维基共享资源： 立体角 当 θ? ?π/2，球冠变为有着立体角 2π的半球。 当 θ? ?π，立体角涵盖整个球体，球冠变为有着立体角 4π的球，我们将4π称为全方位立体角。 [编辑] 任意四面体的立体角 对于任意一个四面体OABC,其中O,A,B,C分别为四面体的四个顶点。下面给出一个公式，计算从O点观察三角形ABC的立体角Ω的方便简单的公式。令α ∠BOC，β ∠AOC，γ ∠AOB 均为各自平面内两条直线的夹角，可以采用平面三角形的余弦公式计算求得 ，S （α+β+γ）/3有 参见L Huilier tan Ω/4 sqrt[tan s/2 tan s/2-α/2 tan s/2-β/2 tan s/2-γ/2 ] 附上相应的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*180./3.14159 end