季 朋，李 晖，陈 梅，戴震宇

(1.贵州大学计算机科学与技术学院,贵州 贵阳 550025;2.贵州省先进计算与医疗信息服务工程实验室，贵州 贵阳 550025)

0 引 言

快速傅里叶变换(FFT)[1]算法是数字信号处理中的核心算法，在频谱分析、数字通信、图像处理中有着重要的应用。快速傅里叶变换在天文学方面也有大量应用，一个典型的应用场景是对脉冲星信号进行相干消色散处理[2]。脉冲星辐射的电磁波在传播过程中，由于受到星际介质的干扰，不同频率的电磁波经过星际介质后，会引起观测到的脉冲星轮廓展宽，所以需要进行消色散处理。消色散处理通常使用相干消色散算法。相干消色散算法的基本原理[3]是将基带脉冲星数字信号进行快速傅里叶变换，然后将频域上各频点的信号乘以星际介质函数CHIRP即频率的响应函数，最终逆傅里叶变换到时域，从而实现将不同频率成分的脉冲星型号对齐到某一频点。傅里叶变换的快慢直接决定相关消色散的处理效率。

为了满足海量信号数据的存储与分析需求，本文采用Greenplum分布式数据库[4-5]存储脉冲星信号。在数据库中对数据计算分析通常使用SQL语句，但SQL语句不便于实现复杂的算法，例如快速傅里叶变换等科学分析算法，而用户自定义函数(UDF)可以满足这一需求。本文实现了UDF形式的快速傅里叶变换算法，为了提高算法效率，对算法的执行做了进一步的优化。

在Greenplum数据库集群中有多个segment实例，每个实例能够并行执行一部分任务。……