基于cuda并行框架的简化球谐波模型快速求解方法-生物医学工程专业论文.docxVIP

万方数据 万方数据 西安电子科技大学 学位论文独创性声明 秉承学校严谨的学风和优良的科学道德，本人声明所呈交的论文是我个人 在导师指导下进行的研究工作及取得的研究成果。尽我所知，除了文中特别加以 标注和致谢中所罗列的内容以外，论文中不包含其他人已经发表或撰写过的研究 成果；也不包含为获得西安电子科技大学或其它教育机构的学位或证书而使用过 的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中做了明确的 说明并表示了谢意。 申请学位论文与资料若有不实之处，本人承担一切法律责任。 本人签名： 日期 西安电子科技大学 关于论文使用授权的说明 本人完全了解西安电子科技大学有关保留和使用学位论文的规定，即：研 究生在校攻读学位期间论文工作的知识产权单位属西安电子科技大学。学校有权 保留送交论文的复印件，允许查阅和借阅论文；学校可以公布论文的全部或部分 内容，可以允许采用影印、缩印或其它复制手段保存论文。同时本人保证，毕业 后结合学位论文研究课题再撰写的文章一律署名单位为西安电子科技大学。 （保密的论文在解密后遵守此规定） 本学位论文属于保密，在 年解密后适用本授权书。 本人签名： 日期 导师签名： 日期 摘 要 摘 要 光学分子成像作为一种迅速发展的医学影像技术，由于具有无创性、高灵敏 度、高时间分辨率、成像价格低等优点而得到了广泛重视和应用，近年来已经成 为了影像领域的研究热点。在光学分子成像领域中，光子在生物组织中传输过程 的精确建模与快速求解具有重要的研究意义。本文围绕生物组织中光传输问题展 开研究，通过详尽地分析各种光传输前向模型，特别是简化球谐波模型（simplified spherical harmonics，SPN）的特点，针对该模型的数值求解方法进行了深入的研究， 并探讨了有限元数值求解方法移植到图形处理器（graphics processing unit，GPU） 并行执行的可行性。 对于多维复杂几何结构，随着近似阶数 N 的增大，SPN 模型将会面临计算量较 大的问题。同时考虑有限元算法作为一种离散化近似求解的数值方法，具有天然 的可并行性。因此，本文提出了一种基于 GPU 高性能计算的并行有限元算法，用 于加速求解基于 SPN 方程的光传输过程。本方法采用由 NVIDIA 公司推出的通用 并行计算框架——统一计算设备架构（compute unified device architecture，CUDA） 作为并行计算解决方案，使其求解结果能够在确保较高精度的同时，又缩短了大 量的计算时间。 在充分考虑 CUDA 平台 GPU 的硬件特点后，本文又对有限元算法程序中计算 量较大的线性方程组求解部分进行了相应的算法优化设计，包括精简优化 GPU 核 心函数代码和对程序结构、存储器、数据通信进行合理优化，充分利用 GPU 资源， 进一步提高并行加速性能。通过与蒙特卡罗仿真结果对比，验证了基于 CUDA 并 行框架的 SPN 快速求解算法的准确性；并通过与并行优化前的有限元算法进行对 比，获取优化后并行算法的加速比，结果表明该算法获得了较好的加速比和并行 加速性能。 关键词：简化球谐波模型 并行有限元算法 图形处理器 统一计算设备架构 ABSTRACT ABSTRACT Optical molecular imaging is a biomedical imaging technology rapidly developing that has received comprehensive study and application due to its advantages such as noninvasiveness, high sensitivity, high time resolution and low cost. The most important research for optical molecular imaging is the accurate modeling of optical transmission and fast solution in biological tissue. The main work of this dissertation is focused on the photon transport problem based on the analysis on the different light transport models. We studied the photon transport issue based on simplified spherical harmonics equations (SPN) in optical