[發(fā)明專利]一種改進(jìn)的ROMP水聲信道估計算法

在審

• 申請?zhí)枺?/td>: 202010576842.7
• 申請日: 2020-06-22
• 公開（公告）號: CN111800362A
• 公開（公告）日: 2020-10-20
• 發(fā)明（設(shè)計）人: 王好賢;李尊琦;周志權(quán)
• 申請（專利權(quán)）人: 哈爾濱工業(yè)大學(xué)（威海）
• 主分類號: H04L27/26
• 分類號: H04L27/26;H04B7/0413;H04B13/02
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種改進(jìn)的ROMP水聲信道估計算法，其特征在于如下步驟：

第一步：根據(jù)ROMP(正則化正交匹配追蹤)算法的正交化原則需要選取合適的正則化系數(shù)r，根據(jù)誤碼率性能選取合適的正則化系數(shù)r，并保證算法收斂，令表示u內(nèi)積向量。J₀表示篩選原子集，正交化原則體現(xiàn)為式(1)；

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

第二步：設(shè)導(dǎo)頻的接收端信號為觀測向量y，導(dǎo)頻信號矩陣為測量矩陣X，稀疏度為K，初始化迭代次數(shù)t＝1，殘差r₀＝y(tǒng)，Λ_t表示t次迭代的索引集合，X_t表示按索引Λ_t選出的X矩陣的列集合；

第三步：通過式(2)計算X與殘差向量的內(nèi)積的絕對值；

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

第四步：若0≤||Λ||₀≤1/2K則選擇u(i)中K/2個最大的值，將這些值對應(yīng)X的序列序號j構(gòu)成集合J中；若1/2K＜||Λ||₀≤3/4K則選擇u(i)中K/4個最大的值，將這些值對應(yīng)X的序列序號j構(gòu)成集合J中；不然，則選擇u(i)中1個最大的值，將這些值對應(yīng)X的序列序號j構(gòu)成集合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，則停止運(yùn)算，并根據(jù)X_t求出最終結(jié)果h，若不滿足上述條件，則跳轉(zhuǎn)至第三步繼續(xù)運(yùn)算。