[发明专利]一种四核并行的大点数脉压数据处理方法

摘要

本发明公开了一种四核并行的大点数脉压算法，利用DIT‑FFT原理，将N点的大点数FFT以及逆FFT各通过两次蝶形运算分解为四个N/4点的组合，分配四个核并行计算，显著缩短了大点数脉压算法的运算时间，提高了系统的实时性，有较好的工程应用价值。

主权项

1.一种四核并行的大点数脉压算法，其特征在于，包括以下步骤：S1、读取输入的脉冲数据，并分配执行并行脉压的四个dsp核，记为核0，核1，核2，核3；S2、核0取点数编号除4余1的数据，生成序列3；核1取点数编号除4余3的数据，生成序列4；核2取点数编号除4余2的数据，生成序列5；核3取点数编号除4余0的数据，生成序列6；S3、四核分别对序列3、序列4、序列5、序列6做原长度四分之一的FFT；S4、四核同步，完成后进入S5；S5、调用SIMD指令，核0将序列4的FFT结果通过_complex_mpysp函数复乘旋转因子，并通过_daddsp函数与序列3的FFT结果相加，得到序列1的前半部分；核1将序列4的FFT结果复乘取反的旋转因子，并加上序列3的FFT结果，得到序列1的后半部分；核2将序列6的FFT结果复乘旋转因子，并与序列5的FFT结果相加，得到序列2的前半部分；核3将序列6的FFT结果复乘取反的旋转因子，并加上序列5的FFT结果，得到序列2的后半部分；S6、四核同步，完成后进入S7；S7、调用SIMD指令，核0将序列2的前半部分复乘旋转因子，并与序列1的前半部分相加，得到原输入数据FFT结果的第一个1/4部分；核1将序列2的后半部分复乘旋转因子，并与序列1的后半部分相加，得到原输入数据FFT结果的第二个1/4部分；核2将序列2的前半部分复乘取反的旋转因子，并与序列1的前半部分相加，得到原输入数据FFT结果的第三个1/4部分；核3将序列2的后半部分复乘旋转因子，并与序列1的后半部分相加，得到原输入数据FFT结果的第四个1/4部分；S8、四核同步，完成后进入S9；S9、调用SIMD指令，将输入数据FFT结果按长度等分为四个部分，核0、核1、核2、核3分别将这四个部分复乘匹配滤波器；S10、四核同步，完成后进入S11；S11、将S9中的复乘结果作为逆FFT的输入数据，核0取点数编号除4余1的数据，生成序列3；核1取点数编号除4余3的数据，生成序列4；核2取点数编号除4余2的数据，生成序列5；核3取点数编号除4余0的数据，生成序列6；S12、四核分别对序列3、序列4、序列5、序列6做原长度四分之一的逆FFT；S13、四核同步，完成后进入S14；S14、调用SIMD指令，核0将序列4的逆FFT结果通过_complex_mpysp函数复乘旋转因子，并通过_daddsp函数与序列3的逆FFT结果相加，得到序列1的前半部分；核1将序列4的逆FFT结果复乘取反的旋转因子，并加上序列3的逆FFT结果，得到序列1的后半部分；核2将序列6的逆FFT结果复乘旋转因子，并与序列5的逆FFT结果相加，得到序列2的前半部分；核3将序列6的逆FFT结果复乘取反的旋转因子，并加上序列5的逆FFT结果，得到序列2的后半部分，所有结果再乘以1/2；S15、四核同步，完成后进入S16；S16、调用SIMD指令，核0将序列2的前半部分复乘旋转因子，并与序列1的前半部分相加，得到原输入数据逆FFT结果的第一个1/4部分；核1将序列2的后半部分复乘旋转因子，并与序列1的后半部分相加，得到原输入数据逆FFT结果的第二个1/4部分；核2将序列2的前半部分复乘取反的旋转因子，并与序列1的前半部分相加，得到原输入数据逆FFT结果的第三个1/4部分；核3将序列2的后半部分复乘旋转因子，并与序列1的后半部分相加，得到原输入数据逆FFT结果的第四个1/4部分，所有结果再乘以1/2；S17、四核同步，完成后得到当前脉冲的脉压结果，返回S1。