买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

摘要：本发明公开了一种基于FFT与apFFT的精密测量雷达二次测频方法，包括以下步骤：将发射波与回波混频后处理得到差频信号序列SIFn；将序列SIFn做apFFT后得到幅度最大值|Yk|及次大值|Yk+r|；求出首次估测频率；计算频偏Δk，若Δk∈[‑0.5,‑0.25]∪[0.25,0.5]，将fIFFIRST作为最终估测频率fIFFINAL；若Δk∈‑0.25,0.25，则将差频信号序列SIFn后N个点加窗得到SWIFn；求SWIFn的能量谱Gx；找出Gx中的最大值Gk1并求出相邻能量谱值；计算差拍信号测试频率并将其作为最终估测频率。本发明结合FFT插值算法和apFFT算法，利用apFFT频谱泄漏较小的优点，在任意可测量的频率都有较好的测量精度。

主权项：1.一种基于FFT与apFFT的精密测量雷达二次测频方法，其特征在于：包括以下步骤：1将发射波与回波混频后经过低通滤波器，通过模数转换得到差频信号序列SIFn；2将序列SIFn做apFFT后得到幅度最大值|Yk|及次大值|Yk+r|，其中r＝±1，k为幅值最大值所对应的位置序列号；若离散谱峰值在主瓣中心右侧即离散谱次大值为Yk+1，r＝1；若离散谱峰值在主瓣中心左侧即离散谱次大值为Yk‑1，r＝‑1；3若r＝1，则将|Yk|，|Yk+r|带入下式：其中，N为FFT点数，fs为采样频率，k为谱线峰值最大值所对应的位置序列号；4若r＝‑1，则|Yk|，|Yk+r|带入下式：5利用步骤3或4求出的频率作为首次估测频率fIFFIRST；6计算频偏Δk，其中Δk为首次测量频率所对应的序列号与k之差，则若Δk∈[‑0.5,‑0.25]∪[0.25,0.5]，将fIFFIRST作为最终估测频率fIFFINAL；若Δk∈‑0.25,0.25，则转入步骤7；7将差频信号序列SIFn后N个点加窗得到SWIFn；8求SWIFn的能量谱Gx；9找出Gx中的最大值Gk1，其中k1为能量谱最大谱线所对应的位置序列号，并求出相邻能量谱值；10将k1的值与步骤9中的相邻能量谱值带入下式：其中，n∈[‑2,2],其中G1k1+i谱序列号为k1+i的能量谱值，并根据下式计算差拍信号测试频率fIFSECEND：最终将得到fIFSECEND作为最终估测fIFFINAL输出。

全文数据：基于FFT与apFFT的精密测量雷达二次测频方法技术领域[0001]本发明涉及一种精密测量雷达测频方法，尤其涉及一种基于FFT与apFFT的精密测量雷达二次测频方法。背景技术[0002]目前，精密测量雷达测频测距方法已经很成熟了，测量方法越来越多，传统方法主要包括RIFE法、能量重心法和FFT细化法，RIFE法是一种插值算法，FFT即快速傅氏变换算法。传统方法计算量普遍偏小，但是都存在一个问题:仅仅在某些频偏所估测的频率与FFT谱线峰值对应频率之差下测频误差较小而在剩余区间测频误差较大。目前较为先进的测量雷达测频方法多基于apFFT，apFFT是全相位FFT的简称，是对传统FFT进行改进而得到的，它具有FFT不具备的“相位不变性”等优良性能;但是，目前采用的apFFT算法也存在相位模糊的问题，从而导致测量精度不高、尤其在某些频率测量精度较低的问题。发明内容[0003]本发明的目的就在于为了解决上述问题而提供一种在任意频率都具有较高测量精度的基于FFT与apFFT的精密测量雷达二次测频方法。[0004]本发明通过以下技术方案来实现上述目的：[0005]—种基于FFT与apFFT的精密测量雷达二次测频方法，包括以下步骤：[0006]1将发射波与回波混频后经过低通滤波器，通过模数转换得到差频信号序列SIFη；[0007]⑵将序列SIF⑹做apFFT后得到幅度最大值|Υ⑹I及次大值IYk+rI，其中r=±l，k为幅值最大值所对应的位置序列号;若离散谱峰值在主瓣中心右侧即离散谱次大值为Yk+1，r=l;若离散谱峰值在主瓣中心左侧即离散谱次大值为Yk-l，r=-l;[0008]3若r=I，则将IY⑹I，IYk+rI带入下式：[0009][0010]其中，N为FFT点数，fs为采样频率，k为谱线峰值最大值所对应的位置序列号；[0011]⑷若r=-l，则IY⑹I，IYk+rI带入下式：[0012][0013]5利用步骤⑶或⑷求出的频率作为首次估测频率fIFFIRST;[0014]⑹计算频偏Λk，其中Δk为首次测量频率所对应的序列号与k之差，则若八1^£[-0·5，-0·25]U[0·25，0·5]，将fiFFiRST作为最终估测频率fiFFiNAL;若Δke-〇·25，0·25，则转入步骤⑺；[0015]⑺将差频信号序列SIFη后N个点加窗得到SWIFη;[0016]⑻求SWIFη的能量谱GX;[0017]⑼找出GX中的最大值GIu，其中Iu为能量谱最大谱线所对应的位置序列号，并求出相邻能量谱值；[0018]10将1^的值与步骤⑼中的相邻能量谱值带入下式：[0019][0020]其中，ne[-2,2]，其中G1Mi谱序列号为h+i的能量谱值，[0021]并根据下式计算差拍信号测试频率flFSECEND:[0022][0023]最终将得到fIFSECEND作为最终估测fIFF皿L输出。[0024]作为优选，所述步骤⑷、⑸中，N=512,fs=512kHz。[0025]本发明的有益效果在于：[0026]本发明结合FFT插值算法和apFFT算法，利用apFFT频谱泄漏较小的优点，对于频偏较大的信号序列，插值方向判定准确，有较好的测量精度，而在频偏绝对值较小时则采用传统能量重心算法，测频误差小;利用本方法在任意可测量的频率都有较好的测量精度。附图说明[0027]图1-1是信噪比SNR=5dB时利用本发明所述基于FFT与apFFT的精密测量雷达二次测频方法获得的均值误差绝对值与频率的对照关系示意图；[0028]图1-2是信噪比SNR=5dB时利用本发明所述基于FFT与apFFT的精密测量雷达二次测频方法获得的均方根误差与频率的对照关系示意图；[0029]图2-1是信噪比SNR=5dB时利用传统apFFT插值测频方法获得的均值误差绝对值与频率的对照关系示意图；[0030]图2-2是信噪比SNR=5dB时利用传统apFFT插值测频方法获得的均方根误差与频率的对照关系示意图；[0031]图3-1是信噪比SNR=5dB时利用传统能量重心测频方法获得的均值误差绝对值与频率的对照关系示意图；[0032]图3-2是信噪比SNR=5dB时利用传统能量重心测频方法获得的均方根误差与频率的对照关系示意图；[0033]图4-1是信噪比SNR=3dB时利用本发明所述基于FFT与apFFT的精密测量雷达二次测频方法获得的均值误差绝对值与频率的对照关系示意图；[0034]图4-2是信噪比SNR=3dB时利用本发明所述基于FFT与apFFT的精密测量雷达二次测频方法获得的均方根误差与频率的对照关系示意图。具体实施方式[0035]下面结合附图和实施例对本发明作进一步说明：[0036]本发明所述基于FFT与apFFT的精密测量雷达二次测频方法，包括以下步骤：[0037]1将发射波与回波混频后经过低通滤波器，通过模数转换得到差频信号序列SIFη；[0038]2将序列SIF⑹做apFFT后得到幅度最大值|Υ⑹I及次大值IYk+rI，其中r=±l，k为幅值最大值所对应的位置序列号;若离散谱峰值在主瓣中心右侧即离散谱次大值为Yk+1，r=l;若离散谱峰值在主瓣中心左侧即离散谱次大值为Yk-l，r=-l;[0039]3若r=l，则将IY⑹I，IYk+rI带入下式：[0040][0041]其中，N为FFT点数，fs为采样频率，k为谱线峰值最大值所对应的位置序列号；[0042]⑷若r=-1，则IY⑹I，IYk+rI带入下式：[0043][0044]5利用步骤⑶或⑷求出的频率作为首次估测频率fIFFIRST;[0045]⑹计算频偏Δk，其中Δk为首次测量频率所对应的序列号与k之差，则若八1^£[-0·5，-0·25]U[0·25，0·5]，将fiFFiRST作为最终估测频率fiFFiNAL;若Δke-〇·25，0·25，则转入步骤⑺；[0046]⑺将差频信号序列SIFη后N个点加窗得到SWIFη;[0047]⑻求SWIFη的能量谱GX;[0048]⑼找出GX中的最大值GIu，其中Iu为能量谱最大谱线所对应的位置序列号，并求出相邻能量谱值；[0049]10将1^的值与步骤⑼中的相邻能量谱值带入下式：[0050][0051]其中，ne[-2,2]，其中Gdlu+i谱序列号为h+i的能量谱值，[0052]并根据下式计算差拍信号测试频率flFSECEND:[0053][0054]最终将得到flFSECEND作为最终估测fIFFINAL输出。[0055]实施例：[0056]图1-1至图4-2的仿真条件同为FFT点数为N=512,采样频率设置为fs=512kHz，差拍信号频率设置为fi「e[30000hz，31000hz]，两两相差10Hz，做1000次蒙特卡洛。[0057]由图1-1、图1-2可知，利用本发明所述基于FFT与apFFT的精密测量雷达二次测频方法获得的均值误差绝对值最高不超过80且多数低于40，均方根误差最大不超过120且多数低于50;对应地，由图2-1、图2-2可知，利用传统apFFT插值测频方法获得的均值误差绝对值最高超过180，均方根误差最大超过220;由图2-1、图2-2可知，利用传统能量重心测频方法获得的均值误差绝对值最高超过80但多数高于50,均方根误差最大不到100但多数高于60;可以看出本发明所述基于FFT与apFFT的精密测量雷达二次测频方法的测量性能卓越。[0058]由图4-1、图4-2可知，本发明所述基于FFT与apFFT的精密测量雷达二次测频方法应用在信噪比SNR=3的环境中时，其测量特性相较于SNR=5时会增大，但相比于传统测频方法仍然具有较高的测量精度。综上可以证明，本发明基于FFT与apFFT的精密测量雷达二次测频方法有着良好的测频性能。[0059]上述实施例只是本发明的较佳实施例，并不是对本发明技术方案的限制，只要是不经过创造性劳动即可在上述实施例的基础上实现的技术方案，均应视为落入本发明专利的权利保护范围内。

权利要求：1.一种基于FFT与apFFT的精密测量雷达二次测频方法，其特征在于:包括以下步骤：1将发射波与回波混频后经过低通滤波器，通过模数转换得到差频信号序列SIFη;⑵将序列SIFη做apFFT后得到幅度最大值|YkI及次大值|Yk+rI，其中r=±l，k为幅值最大值所对应的位置序列号;若离散谱峰值在主瓣中心右侧即离散谱次大值为Yk+I，r=l;若离散谱峰值在主瓣中心左侧即离散谱次大值为Yk-l，r=-l;⑶若r=l，则将|YkI，|Yk+rI带入下式：其中，N为FFT点数，fs为采样频率，k为谱线峰值最大值所对应的位置序列号；⑷若r=-l，则|YkI，|Yk+rI带入下式：⑶利用步骤⑶或⑷求出的频率作为首次估测频率fIFFIRST;6计算频偏△k，其中△k为首次测量频率所对应的序列号与k之差，则若△keίο·5，-O·25]U[0·25，O·5]，将fiFFiRST作为最终估测频率fiFFiNAL;若Δke-〇·25，O·25，则转入步骤⑵；⑵将差频信号序列SIFη后N个点加窗得到SWIFη;⑶求SWIFη的能量谱GX;⑶找出GX中的最大值GIu，其中Iu为能量谱最大谱线所对应的位置序列号，并求出相邻能量谱值；10将匕的值与步骤⑶中的相邻能量谱值带入下式：其中，ne[-2,2]，其中GiOu+i谱序列号为la+i的能量谱值，并根据下式计算差拍信号测试频率flFSECEND:最终将得到flFSEGEND作为最终估测fIFFINAL输出。2.根据权利要求1所述的基于FFT与apFFT的精密测量雷达二次测频方法，其特征在于：所述步骤⑷、（5中，N=512，fs=512kHz。

百度查询： 成都雷尼尔科技有限公司 基于FFT与apFFT的精密测量雷达二次测频方法