主权项 |
适用于正交频分复用接收机抑制测距仪脉冲干扰的方法,其特征在于:步骤一:过采样;以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>Σ</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>[</mo><mi>k</mi><mo>]</mo><mo>=</mo><munderover><mo>Σ</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>[</mo><mi>k</mi><mo>]</mo><mo>+</mo><mi>n</mi><mo>[</mo><mi>k</mi><mo>]</mo><mo>+</mo><msub><mi>s</mi><mi>I</mi></msub><mo>[</mo><mi>k</mi><mo>]</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>[</mo><mi>k</mi><mo>]</mo><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>x</mi><mo>[</mo><mi>k</mi><mo>]</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>x</mi><mo>[</mo><mi>k</mi><mo>]</mo><mo><</mo><mi>μ</mi><mo>*</mo><mi>x</mi><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>x</mi><mo>[</mo><mi>k</mi><mo>]</mo><mo>></mo><mi>μ</mi><mo>*</mo><mi>x</mi><mo>[</mo><mi>k</mi><mo>]</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>Σ</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mi>Φ</mi><mo>`</mo><mo>[</mo><mi>k</mi><mo>]</mo><msub><msup><mi>x</mi><mo>′</mo></msup><mi>I</mi></msub><mo>[</mo><mi>k</mi><mo>]</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>。 |