[发明专利]基于DSP芯片的可变长度FFT计算方法

权利要求书

1.一种基于DSP芯片的可变长度的FFT计算方法，其特征在于步骤为：

(1)判断存储在片外存储器中待运算数据的长度N，若N<N₁，转入执行步骤(2)；若N₁<N<TH₁，转入执行步骤(3)；若TH₁<N<TH₂，转入执行步骤(4)；其中N＝2^k，N₁为DSP芯片内FFT执行装置能够直接支持的最大点数，TH₁为片上存储器能够存储的最大点数，TH₂＝N₁*N₁；

(2)直接以DMA方式从片外存储器中读取初始数据，所述FFT执行装置执行N点一维FFT运算后将运算结果以DMA方式写回片外存储器中，退出计算；

(3)将所述待运算数据视为N₁*N₂的二维数据矩阵并按行顺序读取到片上存储器，在二维数据矩阵中行数为N₁、列数为N₂，其中N＝N₁*N₂；所述FFT执行装置对片上存储器中数据矩阵执行二维FFT运算后得到最终结果矩阵，按列顺序读取所述最终结果矩阵并按行顺序写回片外存储器，退出计算；

(4)将所述待运算数据视为N₁*N₂的二维矩阵并按列分为多个列数据块，每个列数据块包含TH₁个数据点，依次读取每个列数据块到片上存储器中并由所述FFT执行装置进行列方向的FFT计算，将计算结果顺序写回片外存储器；将所述计算结果按行分为多个行数据块，每个行数据块包含TH₁个数据点，依次读取每个行数据块到片上存储器中，由所述FFT执行装置进行行方向的FFT计算，按列顺序读取计算结果并按行顺序写回片外存储器，退出计算。

2.根据权利要求1所述的基于DSP芯片的可变长度的FFT计算方法，其特征在于，所述步骤(3)的具体步骤为：

(3.1)将片外存储器中N点所述待运算数据视为N₁*N₂的二维数据矩阵，按行顺序读取所述二维数据矩阵并以相同顺序存储到片上存储器中，得到初始数据矩阵；

(3.2)所述FFT执行装置读取片上存储器中的所述初始数据矩阵并执行列方向上的N₂次N₁点FFT运算，进行旋转因子补偿后得到列方向运算结果矩阵，并写入片上存储器中对应位置；

(3.3)所述FFT执行装置从片上存储器中顺序读取所述列方向运算结果矩阵并执行行方向上的N₁次N₂点FFT运算，得到最终运算结果矩阵并顺序写入到片上存储器中对应位置；

(3.4)按列顺序从片上存储器中读取所述最终结果矩阵并以行顺序写回片外存储器中，得到最终结果矩阵，退出计算。

3.根据权利要求2所述的基于DSP芯片的可变长度的FFT计算方法，其特征在于，所述步骤(3.2)的具体步骤为：

(3.2.1)以间隔寻址的方式从片上存储器中读取所述初始数据矩阵的一列数据至所述FFT执行装置中；

(3.2.2)由所述FFT执行装置进行N₁点FFT运算，得到当前列数据的蝶形运算结果；对所述蝶形运算结果进行一级补偿旋转因子计算，得到当前列数据的计算结果；以间隔寻址的方式读取下一列数据至所述FFT执行装置中，返回执行步骤(3.2.2)，总计执行N₂次后得到列方向运算结果矩阵；

(3.2.3)将所述列方向运算结果矩阵的每列数据以间隔寻址的方式写入片上存储器的对应位置中，转入执行步骤(3.3)。