DFT计算程序