天津大学 精密仪器与光电子工程学院，天津 300072
从辐射传输方程出发，应用球谐函数方法对辐射率进行多项式展开，并利用球谐函数的正交性和递推性，在二维笛卡儿坐标下，导出了三阶展开的球谐函数微分方程组（P3近似），改进了以往成像文献中忽略各项异性因子的P3近似，并用有限元方法对二维圆域均匀和非均匀两种情况做了数值模拟.与漫射近似模型相比较， P3近似能更准确地描述光源附近及吸收较强情况下边界的光辐射分布情况.
辐射传输方程 球谐函数 扩散方程 有限元方法 蒙特卡洛模拟 Radiative transfer Spherical harmonics Diffusion equation Finite element method Monte Carlo simulation
