[發(fā)明專利]基于四項(xiàng)加權(quán)分?jǐn)?shù)傅里葉變換的窄帶干擾抑制的方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201210120739.7
• 申請(qǐng)日: 2012-04-23
• 公開（公告）號(hào): CN102624660A
• 公開（公告）日: 2012-08-01
• 發(fā)明（設(shè)計(jì)）人: 沙學(xué)軍;李勇;梅林;吳宣利;房宵杰;邱昕;吳瑋
• 申請(qǐng)（專利權(quán)）人: 哈爾濱工業(yè)大學(xué)
• 主分類號(hào): H04L25/03
• 分類號(hào): H04L25/03;H04L27/26
• 代理公司: 哈爾濱市松花江專利商標(biāo)事務(wù)所 23109
• 代理人: 楊立超

1.一種基于四項(xiàng)加權(quán)分?jǐn)?shù)傅里葉變換的窄帶干擾抑制的方法，其特征在于：所述方法的具體過(guò)程為：

步驟一、α階4-WFRFT變換：對(duì)長(zhǎng)度為N的數(shù)字信號(hào)X₀＝{x₁，x₂，...，x_N}進(jìn)行α階4-WFRFT變換，變換到分?jǐn)?shù)域信號(hào)S＝{s₁，s₂，...s_N}，其中s_i表示第i個(gè)分?jǐn)?shù)域采樣點(diǎn)信號(hào)；

步驟二、確定強(qiáng)干擾點(diǎn)的位置：

設(shè)定分?jǐn)?shù)域受到強(qiáng)干擾的點(diǎn)的位置為第i處，設(shè)i處信號(hào)為s_i+n，其中n是強(qiáng)干擾信號(hào)，

si=A0xi+A1yi+A2xN-i+1+A3yN-i+1=[A0,A1,A2,A3]xiyixN-i+1yN-i+1;]]>

A₀～A₃是加權(quán)系數(shù)；

步驟三、系數(shù)替換：

步驟三(一)、利用矩陣廣義逆分離出分?jǐn)?shù)域信號(hào)點(diǎn)s_N-i+1的變換前信號(hào)點(diǎn)[.]^H表示矩陣的轉(zhuǎn)置，分離過(guò)程如下：

記A＝[A₀，A₁，A₂，A₃]，其中A₀～A₃是加權(quán)系數(shù)，由于A是行滿秩矩陣，根據(jù)矩陣的廣義逆可得，A^-1＝A⁺＝A^H(AA^H)^-1，其中A^H表示A的共軛轉(zhuǎn)置，

于是s^N-i+1=[xN-i+1,yN-i+1,xi,yi]H=A-1sN-i+1;]]>

步驟三(二)、根據(jù)s_N-i+1和s_i的表達(dá)式，將第N-i+1點(diǎn)的變換前的信號(hào)點(diǎn)乘上相應(yīng)的加權(quán)系數(shù)得到s_i，將已有的加權(quán)系數(shù)A＝[A₀，A₁，A₂，A₃]經(jīng)過(guò)一定的變換可得到相應(yīng)的加權(quán)系數(shù)，定義如下的置換矩陣：

H=0010000110000100]]>

利用置換矩陣H和矩陣A的乘積得到相應(yīng)的加權(quán)系數(shù)A′＝(HA^H)^H＝AH；

步驟三(三)、將步驟三(二)得到的加權(quán)系數(shù)A′和第i個(gè)點(diǎn)的被加權(quán)函數(shù)的乘積，作為新的分?jǐn)?shù)域的點(diǎn)去代替受強(qiáng)干擾的點(diǎn)；

上述變換系數(shù)替換的整個(gè)過(guò)程由下述公式表達(dá)：

s_i＝AH[A^H(AA^H)^-1s_N-i+1]；

步驟四、-α階4-WFRFT變換：將消除干擾后的變換域信號(hào)，經(jīng)過(guò)4-WFRFT逆變換，得到處理后的時(shí)域輸出信號(hào)。