[发明专利]一种改进的ROMP水声信道估计算法

权利要求书

1.一种改进的ROMP水声信道估计算法，其特征在于如下步骤：

第一步：根据ROMP(正则化正交匹配追踪)算法的正交化原则需要选取合适的正则化系数r，根据误码率性能选取合适的正则化系数r，并保证算法收敛，令表示u内积向量。J₀表示筛选原子集，正交化原则体现为式(1)；

|u_i|≤r|u_j|，i,j∈J₀ (1)

第二步：设导频的接收端信号为观测向量y，导频信号矩阵为测量矩阵X，稀疏度为K，初始化迭代次数t＝1，残差r₀＝y，Λ_t表示t次迭代的索引集合，X_t表示按索引Λ_t选出的X矩阵的列集合；

第三步：通过式(2)计算X与残差向量的内积的绝对值；

u＝abs[X^Tr_t-1] (2)

第四步：若0≤||Λ||₀≤1/2K则选择u(i)中K/2个最大的值，将这些值对应X的序列序号j构成集合J中；若1/2K＜||Λ||₀≤3/4K则选择u(i)中K/4个最大的值，将这些值对应X的序列序号j构成集合J中；不然，则选择u(i)中1个最大的值，将这些值对应X的序列序号j构成集合J中；

第五步：根据式(1)对J做一步滤除操作，并通过式(3)选出满足下式的J₀；

第六步：通过式(4)更新索引集合，其中a_j表示矩阵Α的第j列；

Λ_t＝Λ_t-1∪J₀，X_t＝X_t-1∪{a_j} j∈J₀ (4)

第七步：通过式(5)求最小二乘解；

第八步：通过式(6)更新残差；

第九步：令t＝t+1，若满足t＞K，或||Λ||₀≥K，或残差r_t＝0，则停止运算，并根据X_t求出最终结果h，若不满足上述条件，则跳转至第三步继续运算。