发明名称 适用于正交频分复用接收机抑制测距仪脉冲干扰的方法
摘要 本发明公开了一种适用于正交频分复用接收机抑制测距仪脉冲干扰的方法,属于脉冲干扰抑制技术领域。所述方法通过过采样,对过采样数据进行预处理,采用基于压缩感知的干扰重构方法,得到64路干扰信号,最后将干扰信号从接收信号中消除,得到的期望信号作为接收机干扰消除模块的输出。本发明利用正交匹配跟踪算法,运算速度快,程序上易于实现,能够满足实时处理干扰信号的要求;本发明通过迭代寻找最优化测量矩阵,重构干扰信号概率高,能有效消除干扰;本发明采取的压缩感知算法立足于信号的稀疏度自适应变化,贴合实际测距仪DME信号发生特点,重构成功率高。
申请公布号 CN103560991B 申请公布日期 2016.08.17
申请号 CN201310491566.4 申请日期 2013.10.18
申请人 北京航空航天大学 发明人 张学军;赵亮;白云
分类号 H04L27/26(2006.01)I;H04L25/03(2006.01)I 主分类号 H04L27/26(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 姜荣丽
主权项 适用于正交频分复用接收机抑制测距仪脉冲干扰的方法,其特征在于:步骤一:过采样;以4倍于L‑DACS1系统采样频率的频率进行采样,经过过采样,接收机收到的包含期望信号的接收信号X[k]时域表示如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><msub><mi>x</mi><mi>d</mi></msub><mo>+</mo><mi>n</mi><mo>+</mo><msub><mi>S</mi><mi>I</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>VN</mi><mrow><mi>f</mi><mi>f</mi><mi>t</mi></mrow></msub><mo>-</mo><mn>1</mn></mrow></munderover><mi>x</mi><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>VN</mi><mrow><mi>f</mi><mi>f</mi><mi>t</mi></mrow></msub><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>x</mi><mi>d</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>+</mo><mi>n</mi><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>+</mo><msub><mi>s</mi><mi>I</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000927623890000011.GIF" wi="1318" he="188" /></maths>其中,V为过采样因子,V=4;N<sub>fft</sub>为L‑DACS1系统中正交频分复用系统的FFT长度,N<sub>fft</sub>=64,所述的期望信号即OFDM系统信号x<sub>d</sub>、可加高斯白噪声信号n和测距仪DME干扰信号s<sub>I</sub>;步骤二:对64路经过过采样的接收信号进行预处理,包括以下步骤:A1、对接收信号X的某采样值x[k]进行阈值判断,阈值选择为μ*x[k],μ∈(0,0.6),将大于此阈值的接收信号幅值设置为0,得到削峰后的接收信号x1[k]:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>x</mi><mn>1</mn><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>x</mi><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>x</mi><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>&lt;</mo><mi>&mu;</mi><mo>*</mo><mi>x</mi><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>x</mi><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><mo>&gt;</mo><mi>&mu;</mi><mo>*</mo><mi>x</mi><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000927623890000012.GIF" wi="734" he="123" /></maths>A2、对各列正交的正态高斯随机矩阵进行优化,优化方式如下:第j次迭代时,j∈(1,2,…,200),首先利用下面的公式按列将测量矩阵的每一列单位化:<img file="FDA0000927623890000013.GIF" wi="412" he="87" />其中<img file="FDA0000927623890000014.GIF" wi="64" he="72" />表示第j次迭代中,测量矩阵Φ的第n列;每列单位化结束后,计算得到优化测量矩阵Φ<sup>j+1</sup>:Φ<sup>j+1</sup>=Φ<sup>j</sup>‑0.01Φ<sup>j</sup>(Φ<sup>jT</sup>Φ<sup>j</sup>‑1)其中,Φ<sup>j</sup>表示第j次迭代得到测量矩阵;A3、分块:使用分块向量G∈R<sup>N</sup>处理待处理信号x<sub>I</sub>及优化后的测量矩阵,得到具有块结构的信号x′<sub>I</sub>及最优测量矩阵Φ`,所述分块向量G如下:G=[1,…,1,2,..,2,..,m,…,m]<sup>T</sup>得到:<img file="FDA0000927623890000015.GIF" wi="691" he="210" />x<sup>T</sup>[1]……x<sup>T</sup>[m]分别表示待处理信号的第1、……m子块,m为子块的个数;<img file="FDA0000927623890000016.GIF" wi="670" he="170" />其中,Φ[1]……Φ[m]表示测量矩阵的第1、………m子块,m为子块的个数;A4、使用分块后的最优测量矩阵Φ`对具有块结构的信号x′<sub>I</sub>[k]进行压缩采样,得到:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>y</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mi>&Phi;</mi><mo>`</mo><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo><msub><msup><mi>x</mi><mo>&prime;</mo></msup><mi>I</mi></msub><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000927623890000021.GIF" wi="437" he="175" /></maths>其中,L为分块后的测量矩阵、分块后的待处理信号矩阵中块结构数量的和;步骤三:将步骤二经预处理后的64路信号采用基于压缩感知的干扰重构方法,得到64路干扰信号;步骤四、将步骤三经过压缩感知方法得到的64路干扰信号从接收信号中消除,得到的期望信号x<sub>d</sub>作为接收机干扰消除模块的输出,输入到FFT模块中以进一步处理:x<sub>d</sub>=X‑S<sub>I</sub>。
地址 100191 北京市海淀区学院路37号