CN110726981A - 一种针对固定目标的雷达信号相位差检测方法 - Google Patents
一种针对固定目标的雷达信号相位差检测方法 Download PDFInfo
- Publication number
- CN110726981A CN110726981A CN201911309725.8A CN201911309725A CN110726981A CN 110726981 A CN110726981 A CN 110726981A CN 201911309725 A CN201911309725 A CN 201911309725A CN 110726981 A CN110726981 A CN 110726981A
- Authority
- CN
- China
- Prior art keywords
- radar
- log2
- signal
- phase difference
- matrix
- 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
Images
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种针对固定目标的雷达信号相位差检测方法,通过采用不同宽度的高斯窗对雷达信号进行加权处理,从不同时间‑频率尺度减少干扰信号对有用频率成分的频谱干扰,通过离散频谱峰值相位差的加权平均,获得雷达信号相位差的高精度计算结果,克服了传统雷达信号相位差检测方法难以同时满足抗噪性和检测精度要求的不足。
Description
技术领域
本发明属于信号处理技术领域,尤其涉及一种针对固定目标的雷达信号相位差检测方法。
背景技术
随着微波作为载波,根据发射和发射信号的相位差对固定目标的位置进行测量,具有不受天气情况影响,测点布置灵活等优点。
在实际应用中,由于被测固定目标的环境往往比较复杂,可能存在大量的悬浮的颗粒物质、测量通道的不平整等因素,导致固定目标的位置测量结果准确性较低。目前的针对固定目标的雷达信号相位差检测方法在实际运用中难以同时满足抗干扰和检测精度的需求,因此需要探索新的雷达信号相位差高精度检测方法,对固定目标的位置测量技术进行深入、全面的研究,以满足实际应用要求。
本发明提出的方法,克服了传统雷达信号相位差检测方法难以同时满足抗噪性和检测精度要求的不足。采用不同宽度的高斯窗对雷达信号进行加权处理,减少干扰信号对有用频率成分的频谱干扰,通过离散频谱峰值相位差的加权平均,为雷达信号相位差高精度检测提供了一条有效的途径。
发明内容
本发明提供了一种针对固定目标的雷达信号相位差检测方法,克服了传统雷达信号相位差检测方法难以同时满足抗噪性和检测精度要求的不足。采用不同宽度的高斯窗对雷达信号进行加权处理,从不同时间-频率尺度减少干扰信号对有用频率成分的频谱干扰,通过离散频谱峰值相位差的加权平均,获得雷达信号相位差的高精度计算结果。
本发明的技术方案是构造一种针对固定目标的雷达信号相位差检测方法,所述步骤如下:
步骤1:发射天线所发射的雷达驱动信号x和接收天线所接收的雷达反馈信号y送入模数转换器,分别得到离散化采样后长度为N的雷达驱动信号x1和雷达反馈信号y1,其中采样频率为fs,N为2的整数次幂,根据采样定理频率fs应大于雷达驱动信号最高频率的2倍;
步骤2:利用高斯窗函数构建加窗系数矩阵W,所述加窗系数矩阵W的列数为N,行数为log2(N),其中加窗系数矩阵W的第i行为长度为N、系数为σ(i),第i行高斯窗系数σ(i)=i/log2(N),加窗系数矩阵W表达式为
步骤3:运用W对长度为N的离散采样序列x1和y1进行加窗处理,获得加窗后的信号矩阵x2和y2;
步骤4:对信号矩阵x2和y2的每一行进行离散傅里叶变换,分别获得加窗离散频谱矩阵X2和Y2,分别计算加窗离散频谱矩阵X2和Y2的模值矩阵,记为MX2和MY2;
步骤5:搜索MX2和MY2中每行的最大峰值谱线,设MX2和MY2的第i行中的峰值谱线分别为第k列和第j列,记为MX2(i, k)和MY2(i, j);
步骤6:分别计算加窗离散频谱矩阵X2中第i行、第k列的相位角AX(i, k),加窗离散频谱矩阵Y2中第i行、第j列的相位角AY(i, j);
步骤7:计算加窗离散频谱矩阵X2和Y2中第i行峰值谱线对应的相位角差值AXY(i)=AX(i, k)-AY(i, j);
步骤8:计算加窗系数矩阵W第i行的权系数q(i),其中当i<0.5 log2(N)时,q(i)=i,当i=0.5 log2(N)时,q(i)=i,当i>0.5 log2(N)时,q(i)=log2(N)-i;
优选地,所述步骤2中的加窗系数矩阵中的高斯窗Gi(t)的表达式为:
优选地,所述步骤3中的加窗处理方法为:将长度为N的离散采样序列x1和y1分别扩充成列数为N,行数为log2(N)的矩阵,其中矩阵中的每一行均分别为x1和y1,扩充后的矩阵表达式分别为:
然后将xJ与yJ分别与加窗系数矩阵W进行点乘,获得加窗后的信号矩阵x2和y2。
优选地,所述雷达驱动信号x和雷达反馈信号y在送入模数转换器之前经过放大过滤处理。
本发明有益效果:
本克服了传统雷达信号相位差检测方法难以同时满足抗噪性和检测精度要求的不足。采用不同宽度的高斯窗对雷达信号进行加权处理,从不同时间-频率尺度减少干扰信号对有用频率成分的频谱干扰,通过离散频谱峰值相位差的加权平均,获得雷达信号相位差的高精度计算结果。
附图说明
图1为本发明中雷达信号相位差检测的原理框图。
具体实施方式
下面结合实施例对本发明作进一步详细的描述,但本发明的实施方式不限于此。本发明的实施方式并不受所述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、 替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
本发明提出了一种针对固定目标的雷达信号相位差检测方法。本领域工程技术人员可根据本发明所公开的检测方法编写程序,将所编写的程序下载到测深装置的微控制器中,实现高精度雷达信号相位差检测,从而实现固定目标位置信息的测量,检测方法的原理框图如图1所示。
实施例1:
一种针对固定目标的雷达信号相位差检测方法,所述步骤如下:
步骤1:发射天线所发射的雷达驱动信号x和接收天线所接收的雷达反馈信号y送入模数转换器,分别得到离散化采样后长度为N的雷达驱动信号x1和雷达反馈信号y1,其中采样频率为fs,N为2的整数次幂,根据采样定理频率fs应大于雷达驱动信号最高频率的2倍;
步骤2:利用高斯窗函数构建加窗系数矩阵W,所述加窗系数矩阵W的列数为N,行数为log2(N),其中加窗系数矩阵W的第i行为长度为N、系数为σ(i),第i行高斯窗系数σ(i)=i/log2(N),加窗系数矩阵W表达式为
步骤3:运用W对长度为N的离散采样序列x1和y1进行加窗处理,获得加窗后的信号矩阵x2和y2;
步骤4:对信号矩阵x2和y2的每一行进行离散傅里叶变换,分别获得加窗离散频谱矩阵X2和Y2,分别计算加窗离散频谱矩阵X2和Y2的模值矩阵,记为MX2和MY2;
步骤5:搜索MX2和MY2中每行的最大峰值谱线,设MX2和MY2的第i行中的峰值谱线分别为第k列和第j列,记为MX2(i, k)和MY2(i, j);
步骤6:分别计算加窗离散频谱矩阵X2中第i行、第k列的相位角AX(i, k)=angle(X2(i,k)),加窗离散频谱矩阵Y2中第i行、第j列的相位角AY(i, j) =angle(Y2(i, k));
步骤7:计算加窗离散频谱矩阵X2和Y2中第i行峰值谱线对应的相位角差值AXY(i)=AX(i, k)-AY(i, j);
步骤8:计算加窗系数矩阵W第i行的权系数q(i),其中当i<0.5 log2(N)时,q(i)=i,当i=0.5 log2(N)时,q(i)=i,当i>0.5 log2(N)时,q(i)=log2(N)-i;
步骤9:计算相位差加权平均结果qXY=[q(i)*AXY(i)]/(q(1)+ q(2)+…+q(log2(N))),作为雷达信号相位差的高精度计算结果。
本实施例步骤4中所述对信号矩阵x2和y2的每一行进行离散傅里叶变换和步骤5中对所述搜索MX2和MY2中每行的最大峰值谱线,关于怎样进行对对信号矩阵x2和y2的每一行进行离散傅里叶变换和怎样搜索MX2和MY2中每行的最大峰值谱线,在此不做解释,因为已经有文献公开了此技术,参见:庞浩、 李东霞、 俎云霄、王赞.应用 FFT 进行电力系统谐波分析的改进算法.中国电机工程学报,其中采样频率fs应大于雷达驱动信号最高频率的2倍,是根据采样定量所得,在此不做解释。
实施例2:
在实施例1的基础上,所述步骤2中的加窗系数矩阵中的高斯窗Gi(t)的表达式为:
实施例3:
在实施例1的基础上,所述步骤3中的加窗处理方法为,将长度为N的离散采样序列x1和y1分别扩充成列数为N,行数为log2(N)的矩阵,其中矩阵中的每一行均分别为x1和y1,扩充后的矩阵表达式分别为:
然后将xJ与yJ分别与加窗系数矩阵W进行点乘,获得加窗后的信号矩阵x2和y2。
实施例4:
在实施例1的基础上,所述雷达驱动信号x和雷达反馈信号y在送入模数转换器之前经过放大过滤处理。
表1 测量结果对比
表1提供的数据是在传统的微波窄束测深装置KJC-901的微控制器中未使用本发明所提供的方法和使用本发明所提供的方法的测量结果。
由表1可看出,本发明提供的针对固定目标的雷达信号相位差检测方法,大大的提高了测量精度。
Claims (4)
1.一种针对固定目标的雷达信号相位差检测方法,其特征在于,所述步骤如下:
步骤1:发射天线所发射的雷达驱动信号x和接收天线所接收的雷达反馈信号y送入模数转换器,分别得到离散化采样后长度为N的雷达驱动信号x1和雷达反馈信号y1,其中采样频率为fs,N为2的整数次幂,采样频率fs应大于雷达驱动信号最高频率的2倍;
步骤2:利用高斯窗函数构建加窗系数矩阵W,所述加窗系数矩阵W的列数为N,行数为log2(N),其中加窗系数矩阵W的第i行为长度为N、系数为σ(i),第i行高斯窗系数σ(i)=i/log2(N),加窗系数矩阵W表达式为
步骤3:运用W对长度为N的离散采样序列x1和y1进行加窗处理,获得加窗后的信号矩阵x2和y2;
步骤4:对信号矩阵x2和y2的每一行进行离散傅里叶变换,分别获得加窗离散频谱矩阵X2和Y2,分别计算加窗离散频谱矩阵X2和Y2的模值矩阵,记为MX2和MY2;
步骤5:搜索MX2和MY2中每行的最大峰值谱线,设MX2和MY2的第i行中的峰值谱线分别为第k列和第j列,记为MX2(i, k)和MY2(i, j);
步骤6:分别计算加窗离散频谱矩阵X2中第i行、第k列的相位角AX(i, k)=angle(X2(i,k)),加窗离散频谱矩阵Y2中第i行、第j列的相位角AY(i, j) =angle(Y2(i, k));
步骤7:计算加窗离散频谱矩阵X2和Y2中第i行峰值谱线对应的相位角差值AXY(i)=AX(i, k)-AY(i, j);
步骤8:计算加窗系数矩阵W第i行的权系数q(i),其中当i<0.5 log2(N)时,q(i)=i,当i=0.5 log2(N)时,q(i)=i,当i>0.5 log2(N)时,q(i)=log2(N)-i;
步骤9:计算相位差加权平均结果qXY=[q(i)*AXY(i)]/(q(1)+ q(2)+…+q(log2(N))),作为雷达信号相位差的高精度计算结果。
4.根据权利要求1所述一种针对固定目标的雷达信号相位差检测方法,其特征在于,所述雷达驱动信号x和雷达反馈信号y在送入模数转换器之前经过放大过滤处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911309725.8A CN110726981B (zh) | 2019-12-18 | 2019-12-18 | 一种针对固定目标的雷达信号相位差检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911309725.8A CN110726981B (zh) | 2019-12-18 | 2019-12-18 | 一种针对固定目标的雷达信号相位差检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110726981A true CN110726981A (zh) | 2020-01-24 |
CN110726981B CN110726981B (zh) | 2020-07-10 |
Family
ID=69226019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911309725.8A Active CN110726981B (zh) | 2019-12-18 | 2019-12-18 | 一种针对固定目标的雷达信号相位差检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110726981B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112649678A (zh) * | 2020-12-24 | 2021-04-13 | 广州山锋测控技术有限公司 | 天馈线测量方法、装置、天馈线测量器件和测试仪 |
CN114217288A (zh) * | 2022-02-22 | 2022-03-22 | 湖南纳雷科技有限公司 | 雷达的回波信号的码片间的高度相参的同步方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050203730A1 (en) * | 2004-03-11 | 2005-09-15 | Yoshirou Aoki | Weight function generating method, reference signal generating method, transmission signal generating apparatus, signal processing apparatus and antenna |
CN109283506A (zh) * | 2018-10-23 | 2019-01-29 | 中北大学 | 基于稀疏傅立叶变换的弱目标信号检测装置及方法 |
CN109343020A (zh) * | 2018-11-16 | 2019-02-15 | 南京理工大学 | 一种基于改进窗函数的s变换时频分析方法 |
CN109886280A (zh) * | 2019-02-21 | 2019-06-14 | 西安微电子技术研究所 | 一种基于核相关滤波的异源图像目标匹配方法 |
-
2019
- 2019-12-18 CN CN201911309725.8A patent/CN110726981B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050203730A1 (en) * | 2004-03-11 | 2005-09-15 | Yoshirou Aoki | Weight function generating method, reference signal generating method, transmission signal generating apparatus, signal processing apparatus and antenna |
CN109283506A (zh) * | 2018-10-23 | 2019-01-29 | 中北大学 | 基于稀疏傅立叶变换的弱目标信号检测装置及方法 |
CN109343020A (zh) * | 2018-11-16 | 2019-02-15 | 南京理工大学 | 一种基于改进窗函数的s变换时频分析方法 |
CN109886280A (zh) * | 2019-02-21 | 2019-06-14 | 西安微电子技术研究所 | 一种基于核相关滤波的异源图像目标匹配方法 |
Non-Patent Citations (1)
Title |
---|
黄斌 等: "广义S变换与短时傅里叶变换在地震时频分析中的对比研究", 《中国煤炭地质》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112649678A (zh) * | 2020-12-24 | 2021-04-13 | 广州山锋测控技术有限公司 | 天馈线测量方法、装置、天馈线测量器件和测试仪 |
CN112649678B (zh) * | 2020-12-24 | 2024-05-28 | 广州山锋测控技术有限公司 | 天馈线测量方法、装置、天馈线测量器件和测试仪 |
CN114217288A (zh) * | 2022-02-22 | 2022-03-22 | 湖南纳雷科技有限公司 | 雷达的回波信号的码片间的高度相参的同步方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110726981B (zh) | 2020-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842128B (zh) | 运动目标的声学跟踪方法及装置 | |
CN110275158B (zh) | 基于贝叶斯压缩感知的宽带雷达回波信号参数估计方法 | |
CN103630888B (zh) | 基于对称三角lfmcw雷达的高精度实时微波测速测距装置 | |
CN110726981B (zh) | 一种针对固定目标的雷达信号相位差检测方法 | |
CN103616679B (zh) | 基于差波束调制和波形分析的pd雷达测距测角方法 | |
CN107561508B (zh) | 一种用于匀加速运动目标的相参积累检测方法 | |
CN106896358A (zh) | 一种基于位置旋转变换的高速目标相参积累检测方法 | |
CN105334266A (zh) | 一种岩石声发射源定位方法 | |
CN103513244B (zh) | 一种基于动态规划的多帧相参积累目标检测前跟踪方法 | |
CN109343020B (zh) | 一种基于改进窗函数的s变换时频分析方法 | |
CN104297740B (zh) | 基于相位分析的雷达目标多普勒谱估计方法 | |
CN110208785B (zh) | 基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法 | |
CN106405501B (zh) | 一种基于相位差回归的单声源定位方法 | |
CN109188362B (zh) | 一种麦克风阵列声源定位信号处理方法 | |
CN103983903A (zh) | 利用示波器检测射频信号的变电站全站局部放电定位方法 | |
CN107356908B (zh) | 一种捷变频信号相参积累方法 | |
CN106950528A (zh) | 一种基于线性调频信号的波达方向估计方法 | |
CN102833016B (zh) | 一种在多径环境下针对同源信号时差的测定方法 | |
CN105223481A (zh) | 基于差值能量函数的局部放电特高频信号起始时刻确定方法 | |
CN111352083B (zh) | 一种高频地波雷达多接收通道增益自动校准方法及装置 | |
CN104731762A (zh) | 基于循环移位的立方相位信号参数估计方法 | |
CN108020818A (zh) | 一种在噪声下基于滑动dft的正弦脉冲信号测距方法 | |
CN116381629A (zh) | 基于实时动态卷积的雷达大型目标模拟系统以及方法 | |
CN111856400B (zh) | 一种水下目标声源定位方法及系统 | |
RU2380720C2 (ru) | Способ определения азимутальных и угломестных пеленгов источников радиоизлучения с повышенным быстродействием |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |