发明名称 一种OFDM雷达距离像处理的新方法
摘要 本发明提出一种新方法用来计算OFDM雷达系统的距离像。该方法特点在于直接作用于调制符号,这样能克服传统基带信号相关处理方法的缺点。同时为了验证本方法的可行性用MATLAB对一个特定的OFDM雷达系统进行建模与仿真。用传统相关性算法和本发明提出的新方法计算出OFDM雷达的距离像,并分析这两种方法获得距离像的性能。本方法针对调制符号进行处理,能避免传统基带信号相关性处理方法带来的一系列不利影响,例如OFDM相关性不理想带来的不可消除旁瓣和相关运算计算量大的缺点。
申请公布号 CN104062641B 申请公布日期 2016.09.07
申请号 CN201410324725.6 申请日期 2014.07.09
申请人 重庆大学 发明人 廖勇;陈欢;许锦;何娟
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 代理人
主权项 一种OFDM雷达距离像处理的方法,其信号处理过程为:首先对发射端信号选择调制方法,再将调制的信号经过IDFT/IFFT变换,在IDFT/IFFT模块输出提取出发射端的调制符号序列;该调制方法为:设雷达发射的信号是OFDM多载频相参脉冲串:<maths num="0001"><math><![CDATA[<mrow><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><msub><mi>&pi;f</mi><mn>0</mn></msub><mi>t</mi><mo>)</mo></mrow><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>mT</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000932896350000011.GIF" wi="662" he="140" /></maths>其中f<sub>0</sub>是发射中心载频,M是脉冲串个数,T<sub>r</sub>是脉冲串重复周期,x(t)为OFDM多载频信号的复包率,其表示为:<maths num="0002"><math><![CDATA[<mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>w</mi><mi>n</mi></msub><mi>r</mi><mi>e</mi><mi>c</mi><mi>t</mi><mrow><mo>(</mo><mfrac><mrow><mi>t</mi><mo>-</mo><mi>T</mi><mi>p</mi></mrow><mrow><mi>T</mi><mi>p</mi></mrow></mfrac><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mi>&Delta;</mi><mi>f</mi><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000932896350000012.GIF" wi="765" he="134" /></maths>其中,N是子载频个数,w<sub>n</sub>是n个子载频的频率加权系数,T<sub>p</sub>是脉冲宽度,rect()是矩形窗函数,Δf是相邻子载频的频率间隔,其中Δf应该满足Δf=1/T<sub>p</sub>以确保各个子载波的正交性,这里我们主要关注OFDM调制后的信号x(t),<maths num="0003"><math><![CDATA[<mrow><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>I</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><mi>n</mi><mi>&Delta;</mi><mi>f</mi><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000932896350000013.GIF" wi="589" he="134" /></maths>其中I(n)是包含了调制信息的序列;对应的频域形式X(f)用下面的式子表示:<maths num="0004"><math><![CDATA[<mrow><mi>X</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>I</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msqrt><mi>T</mi></msqrt><mfrac><mrow><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>&pi;</mi><mo>(</mo><mrow><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>n</mi></msub></mrow><mo>)</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><mo>(</mo><mi>&pi;</mi><mo>(</mo><mrow><mi>f</mi><mo>-</mo><msub><mi>f</mi><mi>n</mi></msub></mrow><mo>)</mo><mi>T</mi><mo>)</mo></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000932896350000014.GIF" wi="773" he="135" /></maths>接收天线将接收到的回波通过A/D变换,去掉循环前缀CP,进行串并转换,再进行FFT变换并对子载波进行解映射,再对信号进行解调,最后经并串转换得到接收信号Rx;则OFDM雷达回波r(t)用下面的式子表示:<maths num="0005"><math><![CDATA[<mrow><mi>r</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>&rho;</mi><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mn>2</mn><msub><mi>&pi;f</mi><mn>0</mn></msub><mo>(</mo><mrow><mi>t</mi><mo>-</mo><mi>&tau;</mi></mrow><mo>)</mo><mo>)</mo></mrow><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>MT</mi><mi>r</mi></msub><mo>-</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000932896350000015.GIF" wi="925" he="126" /></maths>其中ρ表示目标散射中心向后散射的强度,D是目标散射中心和雷达之间的径向距离,τ是延时,则τ=2D/c,其中c表示光速,得到计算雷达距离像的公式如下:<img file="FDA0000932896350000016.GIF" wi="470" he="103" />在接收端对接收信号进行采样,将采样后的信号经过DFT/FFT变换,在DFT/FFT模块输出提取出接收到的调制符号序列;将接收的调制符号序列和发射的调制符号序列进行相除获得采样信道脉冲响应,通过点除法可以得到采样后的信道传递函数:<maths num="0006"><math><![CDATA[<mrow><msub><mi>I</mi><mrow><mi>d</mi><mi>i</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>I</mi><mi>r</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mrow><mi>I</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000932896350000021.GIF" wi="325" he="141" /></maths>采样的信道脉冲响应与雷达距离像是相对应的。采样信道脉冲响应可以通过对I<sub>div</sub>(n)进行逆傅里叶变换得到:<maths num="0007"><math><![CDATA[<mrow><mi>h</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mi>D</mi><mi>F</mi><mi>T</mi><mrow><mo>(</mo><mo>{</mo><msub><mi>I</mi><mrow><mi>d</mi><mi>i</mi><mi>v</mi></mrow></msub><mo>(</mo><mi>n</mi><mo>)</mo><mo>}</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>I</mi><mrow><mi>d</mi><mi>i</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mi>n</mi><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>,</mo></mrow>]]></math><img file="FDA0000932896350000022.GIF" wi="1373" he="143" /></maths>最后根据采样信道脉冲响应和OFDM雷达距离像两者之间的关系,获得OFDM雷达的距离像,<maths num="0008"><math><![CDATA[<mrow><mi>h</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>I</mi><mrow><mi>d</mi><mi>i</mi><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>T</mi></mfrac><mi>n</mi><mi>t</mi><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000932896350000023.GIF" wi="662" he="143" /></maths>
地址 400044 重庆市沙坪坝区沙正街174号重庆大学A区通信与测控中心主教2609室