[发明专利]正交追踪算法在GPU上的并行实现方法

权利要求书

1.一种正交追踪算法在GPU上的并行实现方法，其特征在于，包括以下步骤：

S1、在GPU上生成观测矩阵；

S2、使用正交追踪算法反复迭代估计原始信号，利用上述观测矩阵计算出所述原始信号对应的观测数据，并与真实观测数据进行比较，判断是否终止所述迭代操作；

S3、计算所述观测矩阵中与残差相关性最大的列，将其补充到基矩阵中，所述基矩阵为所述观测矩阵的一部分；

S4、利用最小二乘法在所述基矩阵中估计所述原始信号的非零元素，更新原始信号，继续步骤S2。

2.如权利要求1所述的方法，其特征在于，在步骤S1中，所述观测矩阵为对DCT离散余弦变换矩阵进行随机抽行得到的矩阵，其元素根据下面的公式计算：

Φ(m,n)=1N,h(m)=0,0≤n≤N-12Ncosπ(2n+1)h(m)2N,1≤h(m)≤N-1,0≤n≤N-1]]>

其中，Φ为观测矩阵，m为行数，m∈(0,1,2,…,N-1)，n为列数，N为正交追踪算法待恢复信号的长度，h为计算机生成的伪随机数序列，该序列中有M个元素，M为压缩感知中的观测数，且M<N。

3.如权利要求2所述的方法，其特征在于，GPU将所述元素的生成任务分配到64个线程中,其中第i个线程负责生成(Φ(i,0),Φ(i,1),Φ(i,2),...,Φ(i,N-1))，多个线程在多个处理器上并行完成观测矩阵Φ的生成。