CN105629219A - 均衡测距精度与运算量的方法 - Google Patents
均衡测距精度与运算量的方法 Download PDFInfo
- Publication number
- CN105629219A CN105629219A CN201511008974.5A CN201511008974A CN105629219A CN 105629219 A CN105629219 A CN 105629219A CN 201511008974 A CN201511008974 A CN 201511008974A CN 105629219 A CN105629219 A CN 105629219A
- Authority
- CN
- China
- Prior art keywords
- czt
- data
- distance accuracy
- target
- computing amount
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
均衡测距精度与运算量的方法,设定数据长度为Data_N,FFT变换点数为FFT_N,所需要分析的数据为Czt_sig,进行CZT所得到的点数为Czt_Pp_M,根据门限检测模块,得到目标数目N_mubiao及各个目标在频谱中的位置k1;如目标数目N_mubiao<=N1,则直接对该点与前后两个点进行CZT计算峰值位置;如目标数目N_mubiao>N1,则对判决出的前N1个目标采用CZT计算峰值位置,对其他数据采用二项式拟合计算峰值位置,根据峰值位置来获得差频频率。该方法既能保证在主要目标的距离分辨力和分辨率,又能在现有的硬件基础上,满足系统的实时性要求。
Description
技术领域
本发明属于测距精度技术领域,具体说是一种均衡测距精度与运算量的方法。
背景技术
为提高测距精度,则需要提高频谱峰值的估计精度,可采用的方法有对需要分析的频谱区间进行插值算法、CZT即线性调频Z变换或者小波分析。
二项式拟合属于插值算法的一种,如果二项式拟合估计频谱峰值,频谱估计的精度不高,采用CZT或者小波分析的方法,则运算量很大,尤其是当目标数目较多,如果对每个频谱峰值都进行CZT,由于运算量过大,甚至会造成系统死机。
发明内容
针对现有技术存在上述缺点或者不足,本发明提供了一种均衡测距精度与运算量的方法,既能保证在主要目标的距离分辨力和分辨率,又能在现有的硬件基础上,满足系统的实时性要求。
为实现上述目的,本发明的技术方案是,均衡测距精度与运算量的方法,设定数据长度为Data_N,FFT变换点数为FFT_N,所需要分析的数据为Czt_sig,进行CZT所得到的点数为Czt_Pp_M,根据门限检测模块,得到目标数目N_mubiao及各个目标在频谱中的位置k1;
如目标数目N_mubiao<=N1,则直接对该点与前后两个点进行CZT计算峰值位置;如目标数目N_mubiao>N1,则对判决出的前N1个目标采用CZT计算峰值位置,对其他数据采用二项式拟合计算峰值位置,根据峰值位置来获得差频频率。
进一步的,所述CZT计算峰值位置的步骤具体为:
S1.设置根据目标所在的位置k1,计算CZT变换的起始相位角;
S2.求出抽样点的矢量长度Czt_A和螺线的伸展率Czt_W:
S3.计算中间变量信号L点序列Czt_g(n)和L点序列Czt_h(n);
S4.进而求取Czt_g(n)和Czt_h(n)的圆周卷积Czt_q(n);
S5.最后求得所需的CZT后的频谱信息:
进一步的,步骤S1中起始相位角为:
Czt_qs_jd=2*π*(k1-1)/Data_N。
进一步的,步骤S2中抽样点的矢量长度Czt_A和螺线的伸展率Czt_W为:
Czt_A=Czt_A0*exp(Czt_qs_jd)
Czt_W=Czt_W0*exp(-j*(k2-k1)/(2*Czt_Pp_M)2*π/Data_N)。
进一步的,计算中间变量信号L点序列Czt_g(n)和L点序列Czt_h(n)为:
进一步的,步骤S4中先求得Czt_g(n)和Czt_h(n)的频谱Czt_G_P和Czt_H_P,将其频谱相乘然后取逆FFT变换IFFT得到圆周卷积Czt_q(n)。
进一步的,步骤S5中所需的CZT后的频谱信息为:
进一步的,所述的二项式拟合的具体为:采用y=-a*(x-b)2+c来实现,其中a、b、c为需要确定的系数,x为所需处理的数据在频谱中的位置,y为实际测得的数据,设需要处理的数据为(x1,y1),(x2,y2),(x3,y3);其中x3=x2+1=x1+2;则
y1=-a*(x1-b)2+c①
y2=-a*(x2-b)2+c=-a*(x1+1-b)2+c②
y3=-a*(x3-b)2+c=-a*(x1+2-b)2+c③
公式②减去①化简可得:
y2-y1=-a*(2*x1-2*b+1)④
公式③减去①化简可得:
y3-y1=-2*a*(2*x1-2*b+2)⑤
将④带入⑤,化简可得:
a=y2-1/2*(y1+y3)⑥
b=x1+1/2-(y1-y2)/(2*a)⑦
根据公式⑦和三个数据点的大小及位置,得出更精确的频谱峰值点位置b,该点对应的信号频率f为
f=b*fs/N_FFT。
作为更进一步的,CZT采用递推算法。
本发明由于采用以上技术方案,能够取得如下的技术效果:既能保证在主要目标的距离分辨力和分辨率,又能在现有的硬件基础上,满足系统的实时性要求。
本发明可通过设定所需要进行CZT分析的目标数目来均衡测距精度和运算量增加之间的矛盾。
附图说明
本发明共有附图3幅:
图1为本发明流程框图;
图2为CZT的原理框图;
图3为CZT的实现框图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步的具体说明。
实施例1
均衡测距精度与运算量的方法,设定数据长度Data_N为996,FFT变换点数FFT_N为1024,所需要分析的数据为Czt_sig,进行CZT所得到的点数为Czt_Pp_M,根据门限检测模块,得到目标数目N_mubiao及各个目标在频谱中的位置k1(目标频谱的峰值点位置);
如目标数目N_mubiao<=4,则直接对该点与前后两个点进行CZT计算峰值位置;如目标数目N_mubiao>4,则对判决出的前N1个目标采用CZT计算峰值位置,对其他数据采用二项式拟合计算峰值位置,根据峰值位置来获得差频频率。
所述CZT计算峰值位置的步骤具体为:
S1:根据目标所在的位置k1,计算出CZT变换的起始相位角
Czt_qs_jd=2*π*(k1-1)/Data_N
初始化抽样点的矢量长度Czt_A0=1,螺线的伸展率Czt_W0=1;在本设计中设定需扩展的序列长度Czt_L=1024.
S2:计算抽样点的矢量长度Czt_A和螺线的伸展率Czt_W:
Czt_A=Czt_A0*exp(Czt_qs_jd)
Czt_W=Czt_W0*exp(-j*(k2-k1)/(2*Czt_Pp_M)2*π/Data_N)
S3:计算中间变量信号L点序列Czt_g(n)和L点序列Czt_h(n):
(2-12)
为减少运算量,实现过程中采用如图3所示的递推算法来计算。
S4:求取Czt_g(n)和Czt_h(n)的圆周卷积Czt_q(n),可先求得求得Czt_g(n)和Czt_h(n)的频谱Czt_G_P和Czt_H_P,将其频谱相乘然后取IFFT(逆FFT变换)来得到。
S5:最后根据下式求得所需要的CZT后的频谱信息:
实施例2
本实施例以三角波的测距为例进行说明,数据长度Data_N,FFT变换点数FFT_N,采样率为fs,如果门限检测出目标的数目N_mubiao小于N1个,则只进行CZT变换来提高频谱计算精度及目标的距离分辨力。如果目标数目N_mubiao大于N1,则对其中的重要的N1个目标进行CZT变换提高频谱计算精度及目标的距离分辨力,而对其余目标则采用二次多项式拟合来计算频谱的峰值点来提高频谱计算精度,进而提高距离计算精度。
本设计中多项式采用y=-a*(x-b)2+c来实现,其中a、b、c为需要确定的系数,x为所需要处理的数据在频谱中的位置,y为实际测得的数据。假设需要处理的三个数据为(x1,y1),(x2,y2),(x3,y3);其中x3=x2+1=x1+2;则
y1=-a*(x1-b)2+c①
y2=-a*(x2-b)2+c=-a*(x1+1-b)2+c②
y3=-a*(x3-b)2+c=-a*(x1+2-b)2+c③
公式②减去①化简可得:
y2-y1=-a*(2*x1-2*b+1)④
公式③减去①化简可得:
y3-y1=-2*a*(2*x1-2*b+2)⑤
将④带入⑤,化简可得:
a=y2-1/2*(y1+y3)⑥
b=x1+1/2-(y1-y2)/(2*a)⑦
根据公式⑦和三个数据点的大小及位置,得出更精确的频谱峰值点位置b,该点对应的信号频率f为
f=b*fs/N_FFT。
为了减少运算量,CZT采用递推算法,CZT变换的原理框图如图2所示,实现框图如图3所示。
在申请中,根据系统要求来设定运算量,得到一个周期数据最多能对N1+1个目标进行CZT,由于数据拟合运算量远小于CZT,目标的距离最近在频谱中位置越靠近零频谱线,故当目标大于N1时,可采用对距离零频位置最小的N1个目标采用CZT,而对其余目标均采用二项式拟合。CZT可以提高系统的距离分辨力(分辨出两个不同距离目标的能力)及分辨率(分辨出单个目标的距离精度),而采用二次多项式拟合只能提高系统的距离分辨率。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (10)
1.均衡测距精度与运算量的方法,其特征在于,设定数据长度为Data_N,FFT变换点数为FFT_N,所需要分析的数据为Czt_sig,进行CZT所得到的点数为Czt_Pp_M,根据门限检测模块,得到目标数目N_mubiao及各个目标在频谱中的位置k1;
如目标数目N_mubiao<=N1,则直接对该点与前后两个点进行CZT计算峰值位置;如目标数目N_mubiao>N1,则对判决出的前N1个目标采用CZT计算峰值位置,对其他目标采用二项式拟合计算峰值位置,根据峰值位置来获得差频频率。
2.根据权利要求1所述的均衡测距精度与运算量的方法,其特征在于,所述CZT计算峰值位置的步骤具体为:
S1.根据目标所在的位置k1,计算CZT变换的起始相位角;
S2.求出抽样点的矢量长度Czt_A和螺线的伸展率Czt_W;
S3.计算中间变量信号L点序列Czt_g(n)和L点序列Czt_h(n);
S4.进而求取Czt_g(n)和Czt_h(n)的圆周卷积Czt_q(n);
S5.最后求得所需的CZT后的频谱信息。
3.根据权利要求2所述的均衡测距精度与运算量的方法,其特征在于,步骤S1中起始相位角为:
Czt_qs_jd=2*π*(k1-1)/Data_N。
4.根据权利要求2所述的均衡测距精度与运算量的方法,其特征在于,步骤S2中抽样点的矢量长度Czt_A和螺线的伸展率Czt_W为:
Czt_A=Czt_A0*exp(Czt_qs_jd)
Czt_W=Czt_W0*exp(-j*(k2-k1)/(2*Czt_Pp_M)2*π/Data_N)。
5.根据权利要求2所述的均衡测距精度与运算量的方法,其特征在于,计算中间变量信号L点序列Czt_g(n)和L点序列Czt_h(n)为:
6.根据权利要求2所述的均衡测距精度与运算量的方法,其特征在于,步骤S4中先求得Czt_g(n)和Czt_h(n)的频谱Czt_G_P和Czt_H_P,将其频谱相乘然后取逆FFT变换IFFT得到圆周卷积Czt_q(n)。
7.根据权利要求2所述的均衡测距精度与运算量的方法,其特征在于,步骤S5中所需的CZT后的频谱信息为:
8.根据权利要求1所述的均衡测距精度与运算量的方法,其特征在于,所述的二项式拟合的具体为:采用y=-a*(x-b)2+c来实现,其中a、b、c为需要确定的系数,x为所需处理的数据在频谱中的位置,y为实际测得的数据,设需要处理的数据为(x1,y1),(x2,y2),(x3,y3);其中x3=x2+1=x1+2;则
y1=-a*(x1-b)2+c①
y2=-a*(x2-b)2+c=-a*(x1+1-b)2+c②
y3=-a*(x3-b)2+c=-a*(x1+2-b)2+c③
公式②减去①化简可得:
y2-y1=-a*(2*x1-2*b+1)④
公式③减去①化简可得:
y3-y1=-2*a*(2*x1-2*b+2)⑤
将④带入⑤,化简可得:
a=y2-1/2*(y1+y3)⑥
b=x1+1/2-(y1-y2)/(2*a)⑦
根据公式⑦和三个数据点的大小及位置,得出更精确的频谱峰值点位置b,该点对应的信号频率f为
f=b*fs/N_FFT。
9.根据权利要求1-7任一项所述的均衡测距精度与运算量的方法,其特征在于,CZT采用递推算法。
10.根据权利要求1所述的均衡测距精度与运算量的方法,其特征在于,N1的取值为3-4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511008974.5A CN105629219B (zh) | 2015-12-29 | 2015-12-29 | 均衡测距精度与运算量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511008974.5A CN105629219B (zh) | 2015-12-29 | 2015-12-29 | 均衡测距精度与运算量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105629219A true CN105629219A (zh) | 2016-06-01 |
CN105629219B CN105629219B (zh) | 2018-03-06 |
Family
ID=56044335
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511008974.5A Active CN105629219B (zh) | 2015-12-29 | 2015-12-29 | 均衡测距精度与运算量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105629219B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107390210A (zh) * | 2017-06-22 | 2017-11-24 | 电子科技大学 | 一种物位测量中差拍信号的数字处理方法 |
CN112099038A (zh) * | 2020-09-17 | 2020-12-18 | 上海波汇科技有限公司 | 一种基于fmcw激光雷达的多物体识别方法及装置 |
CN114942053A (zh) * | 2022-04-28 | 2022-08-26 | 国家石油天然气管网集团有限公司华南分公司 | 一种基于相位估值的高精度雷达物位计测距方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5282154A (en) * | 1992-06-01 | 1994-01-25 | Thomson Consumer Electronics, Inc. | System for determining and repairing instability in an IIR filter suitable for use in deghosting apparatus |
CN1520082A (zh) * | 2003-01-23 | 2004-08-11 | 上海贝尔阿尔卡特股份有限公司 | 数字通信系统中载频同步方法和通信装置 |
US7436349B1 (en) * | 2006-06-26 | 2008-10-14 | Sandia Corporation | Controlling data collection to support SAR image rotation |
EP2650695A1 (en) * | 2012-08-02 | 2013-10-16 | Institute of Electronics, Chinese Academy of Sciences | Imaging method for synthetic aperture radar in high squint mode |
CN103823215A (zh) * | 2014-03-03 | 2014-05-28 | 中国科学院电子学研究所 | 线性调频连续波雷达测距方法 |
CN104819841A (zh) * | 2015-05-05 | 2015-08-05 | 西安交通大学 | 基于内置编码信息的单传感柔性角度域平均方法 |
-
2015
- 2015-12-29 CN CN201511008974.5A patent/CN105629219B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5282154A (en) * | 1992-06-01 | 1994-01-25 | Thomson Consumer Electronics, Inc. | System for determining and repairing instability in an IIR filter suitable for use in deghosting apparatus |
CN1520082A (zh) * | 2003-01-23 | 2004-08-11 | 上海贝尔阿尔卡特股份有限公司 | 数字通信系统中载频同步方法和通信装置 |
US7436349B1 (en) * | 2006-06-26 | 2008-10-14 | Sandia Corporation | Controlling data collection to support SAR image rotation |
EP2650695A1 (en) * | 2012-08-02 | 2013-10-16 | Institute of Electronics, Chinese Academy of Sciences | Imaging method for synthetic aperture radar in high squint mode |
CN103823215A (zh) * | 2014-03-03 | 2014-05-28 | 中国科学院电子学研究所 | 线性调频连续波雷达测距方法 |
CN104819841A (zh) * | 2015-05-05 | 2015-08-05 | 西安交通大学 | 基于内置编码信息的单传感柔性角度域平均方法 |
Non-Patent Citations (3)
Title |
---|
X XU等: ""Research on a kind of high precision and fast signal processing algorithm for FM/CW laser radar"", 《PROCEEDINGS OF SPIE-THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING》 * |
周旭辉等: ""基于正交多项式拟合的风廓线雷达峰谱识别"", 《现代雷达》 * |
钱克矛等: ""频谱校正的线性调频Z变换方法"", 《振动工程学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107390210A (zh) * | 2017-06-22 | 2017-11-24 | 电子科技大学 | 一种物位测量中差拍信号的数字处理方法 |
CN112099038A (zh) * | 2020-09-17 | 2020-12-18 | 上海波汇科技有限公司 | 一种基于fmcw激光雷达的多物体识别方法及装置 |
CN112099038B (zh) * | 2020-09-17 | 2023-11-10 | 上海波汇科技有限公司 | 一种基于fmcw激光雷达的多物体识别方法及装置 |
CN114942053A (zh) * | 2022-04-28 | 2022-08-26 | 国家石油天然气管网集团有限公司华南分公司 | 一种基于相位估值的高精度雷达物位计测距方法 |
CN114942053B (zh) * | 2022-04-28 | 2024-05-31 | 国家石油天然气管网集团有限公司华南分公司 | 一种基于相位估值的高精度雷达物位计测距方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105629219B (zh) | 2018-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103323822B (zh) | 一种估计通道误差的方法及装置 | |
CN107229048B (zh) | 一种高分宽幅sar动目标速度估计与成像方法 | |
CN113340191B (zh) | 时间序列干涉sar的形变量测量方法及sar系统 | |
US10009594B2 (en) | Stereo image processing device and stereo image processing method | |
CN103885061B (zh) | 基于改进相位梯度自聚焦的双基sar运动补偿方法 | |
CN105974414A (zh) | 基于二维自聚焦的高分辨聚束sar自聚焦成像方法 | |
CN103630878B (zh) | 一种基于迭代自适应方法的频带融合方法及装置 | |
CN103616687B (zh) | 分段线性估计的多项式拟合isar包络对齐方法 | |
CN104007318B (zh) | 获取信号时频函数的方法 | |
CN105629219A (zh) | 均衡测距精度与运算量的方法 | |
CN102621549A (zh) | 多基线/多频段干涉相位解缠频域快速算法 | |
CN108776342A (zh) | 一种高速平台sar慢速动目标检测与速度估计方法 | |
CN112098970A (zh) | 一种用于交通微波检测的速度解模糊算法及相关设备 | |
Moreau et al. | Estimation of power spectral density from laser Doppler data via linear interpolation and deconvolution | |
CN106842196B (zh) | 一种用于sar子孔径成像的方位模糊抑制方法 | |
CN103728619A (zh) | 基于变重频技术的机载大斜视条带sar成像方法 | |
CN102207546B (zh) | 一种子孔径高度计回波视配准装置及方法 | |
CN101806893A (zh) | 用于合成孔径雷达点目标成像质量评估的自适应二维插值方法 | |
CN103760534A (zh) | 一种星载sar数据的电离层色散效应校正方法 | |
CN110988873B (zh) | 基于能量中心提取的单通道sar舰船速度估计方法及系统 | |
KR101627419B1 (ko) | 이동 노드의 위치 추정 방법 및 그 장치 | |
JP2010127771A (ja) | 合成開口ソーナー、合成開口ソーナーの位相誤差補正方法及びプログラム | |
CN113933804A (zh) | 一种图像域方位向多通道sar误差校正方法 | |
CN106199607A (zh) | 一种麦克风阵列的声源方向定位方法及装置 | |
CN105572648A (zh) | 一种合成孔径雷达回波数据距离徙动校正方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |