CN113534145A - 一种基于线性调频连续波体制的高度表测高方法及系统 - Google Patents

一种基于线性调频连续波体制的高度表测高方法及系统 Download PDF

Info

Publication number
CN113534145A
CN113534145A CN202111083383.XA CN202111083383A CN113534145A CN 113534145 A CN113534145 A CN 113534145A CN 202111083383 A CN202111083383 A CN 202111083383A CN 113534145 A CN113534145 A CN 113534145A
Authority
CN
China
Prior art keywords
signal
frequency
difference frequency
chirp
modulation period
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.)
Pending
Application number
CN202111083383.XA
Other languages
English (en)
Inventor
李贤文
刘辉
宫旭泽
周震
谷明亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Andawell Aviation Equipment Co Ltd
Original Assignee
Beijing Andawell Aviation Equipment Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Andawell Aviation Equipment Co Ltd filed Critical Beijing Andawell Aviation Equipment Co Ltd
Priority to CN202111083383.XA priority Critical patent/CN113534145A/zh
Publication of CN113534145A publication Critical patent/CN113534145A/zh
Pending legal-status Critical Current

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
    • 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/882—Radar or analogous systems specially adapted for specific applications for altimeters
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
    • 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)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种基于线性调频连续波体制的高度表测高方法及系统。该方法包括获取由地面反射的回波信号；根据回波信号和当前发射的发射信号确定差频信号；对差频信号进行预处理；对预处理后的差频信号进行FFT运算，当FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内时，调整当前发射的发射信号的调制周期；反之，对FFT运算后的差频信号进行Chirp_Z变换，当Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率时，调整当前发射的发射信号的调制周期；反之，固定当前发射的发射信号的调制周期，并根据调制周期确定当前高度值。本发明能够在保证测高表测高精度的基础上，提高测高表的实时性。

Description

一种基于线性调频连续波体制的高度表测高方法及系统
技术领域
本发明涉及航空航天行业领域，特别是涉及一种基于线性调频连续波体制的高度表测高方法及系统。
背景技术
目前市场上基于线性调频连续波体制的高度表在军用和民用领域有广泛的应用，其测高原理基本相同，系统通常采用锯齿波调制，将发射信号与回波信号之间的时延转化为差频信号频率值，利用FFT等数字信号处理方法对差频信号进行求取，从而求出飞行高度。由于FFT的栅栏效应导致测量数据的误差，通过增加FFT的点数可以提高系统的测高精度，但势必会计算量从而降低系统的实时性能。
因此，为了解决上述问题，亟需一种新的高度表测高方法或系统。
发明内容
本发明的目的是提供一种基于线性调频连续波体制的高度表测高方法及系统，能够在保证测高表测高精度的基础上，提高测高表的实时性。
为实现上述目的，本发明提供了如下方案：
一种基于线性调频连续波体制的高度表测高方法，包括：
获取由地面反射的回波信号；所述回波信号为发射天线发射的发射信号到达地面，由地面反射回来的信号；所述发射信号为通过锯齿波调制并随时间作线性变换的调制信号；
根据回波信号和当前发射的发射信号确定差频信号；
对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；
若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；
若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；
若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；
若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值。
可选地，所述获取由地面反射的回波信号，具体包括：
利用公式
Figure DEST_PATH_IMAGE001
确定回波信号的瞬时频率；
其中，
Figure 399378DEST_PATH_IMAGE002
为回波信号的瞬时频率，
Figure DEST_PATH_IMAGE003
为载波信号频率，k为调频斜率，
Figure 16217DEST_PATH_IMAGE004
，
Figure DEST_PATH_IMAGE005
为锯齿波的调制周期，
Figure 665504DEST_PATH_IMAGE006
为调频带宽，
Figure DEST_PATH_IMAGE007
为获取回波信号的时间，
Figure 717642DEST_PATH_IMAGE008
为延迟时间，
Figure DEST_PATH_IMAGE009
。
可选地，所述根据回波信号和当前发射的发射信号确定差频信号，具体包括：
利用公式
Figure 292718DEST_PATH_IMAGE010
确定差频信号；
其中，
Figure DEST_PATH_IMAGE011
为差频信号。
可选地，所述若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值，具体包括：
利用公式
Figure 220354DEST_PATH_IMAGE012
确定当前高度值；
其中，
Figure DEST_PATH_IMAGE013
为Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率时的发射信号的调制周期，H为当前高度值，C为光速。
一种基于线性调频连续波体制的高度表测高系统，包括：
回波信号获取模块，用于获取由地面反射的回波信号；所述回波信号为发射天线发射的发射信号到达地面，由地面反射回来的信号；所述发射信号为通过锯齿波调制并随时间作线性变换的调制信号；
差频信号确定模块，用于根据回波信号和当前发射的发射信号确定差频信号；
预处理模块，用于对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
第一判断模块，用于对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；
调制周期第一调整模块，用于若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；
第二判断模块，用于若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；
调制周期第二调整模块，用于若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；
当前高度值确定模块，用于若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值。
一种基于线性调频连续波体制的高度表测高系统，应用于所述的一种基于线性调频连续波体制的高度表测高方法，包括：接收机、发射机、调制器、混频器、信号预处理装置、信号处理装置以及接口单元；
所述调制器与发射机连接；所述调制器用于生成发射信号，并将所述调制信号发送至发射机；所述发射信号通过锯齿波调制并随时间作线性变换的调制信号；
所述混频器分别与接收机和发射机连接；所述混频器用于根据接收机接收的回波信号和发射机发射的发射信号生成差频信号；
所述信号预处理装置与所述混频器连接；所述信号预处理装置用于对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
所述信号处理装置分别与信号预处理装置、调制器以及接口单元连接；所述信号处理装置用于对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值；并将当前高度值发送至接口单元。
可选地，所述信号预处理装置包括：滤波器、放大器以及模数转换器；
所述滤波器、放大器以及模数转换器依次连接。
可选地，所述发射机包括：发射天线和功分器；
所述功分器用于将发射信号分成两路，一路至发射天线发射出去，另一路作为本振信号进入混频器。
可选地，所述接收机包括：接收天线。
可选地，还包括：数模转换器；
所述数模转换器分别与调制器和信号处理装置连接。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明所提供的一种基于线性调频连续波体制的高度表测高方法及系统，对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内，当FFT运算后的差频信号的幅值对应的频率处于设定频率范围内时，对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；当Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率时，固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值。并且，当不满足判断条件时，则重新调整发射信号的调制周期，即通过调整调制周期，使防止FFT的栅栏效应导致测量数据的误差。能够在保证测高表测高精度的基础上，提高测高表的实时性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明所提供的一种基于线性调频连续波体制的高度表测高方法流程示意图；
图2为锯齿波调制下FMCW雷达时频曲线图；
图3为一种基于线性调频连续波体制的高度表测高系统结构示意图；
图4为一种基于线性调频连续波体制的高度表测高系统原理示意图；
图5为信号处理装置原理示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种基于线性调频连续波体制的高度表测高方法及系统，能够在保证测高表测高精度的基础上，提高测高表的实时性。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明所提供的一种基于线性调频连续波体制的高度表测高方法流程示意图，如图1所示，本发明所提供的一种基于线性调频连续波体制的高度表测高方法，包括：
S101，获取由地面反射的回波信号；所述回波信号为发射天线发射的发射信号到达地面，由地面反射回来的信号；所述发射信号为通过锯齿波调制并随时间作线性变换的调制信号；
S101具体包括：
利用公式
Figure 922599DEST_PATH_IMAGE001
确定回波信号的瞬时频率；
其中，
Figure 704611DEST_PATH_IMAGE002
为回波信号的瞬时频率，
Figure 202586DEST_PATH_IMAGE003
为载波信号频率，k为调频斜率，
Figure 820781DEST_PATH_IMAGE004
，
Figure 670925DEST_PATH_IMAGE005
为锯齿波的调制周期，
Figure 556710DEST_PATH_IMAGE006
为调频带宽，
Figure 443895DEST_PATH_IMAGE007
为获取回波信号的时间，
Figure 798653DEST_PATH_IMAGE008
为延迟时间，
Figure 295231DEST_PATH_IMAGE009
。
如图2所示，锯齿波的调制信号
Figure 536987DEST_PATH_IMAGE014
为：
Figure DEST_PATH_IMAGE015
；
将锯齿波调制电压信号加载到压控振荡器（VCO）上，控制其产生等幅连续调频信号，设等幅连续调频信号
Figure 31292DEST_PATH_IMAGE016
为：
Figure DEST_PATH_IMAGE017
；
其中，
Figure 279870DEST_PATH_IMAGE018
表示
Figure DEST_PATH_IMAGE019
的幅度，
Figure 399048DEST_PATH_IMAGE020
为载波信号频率，对
Figure DEST_PATH_IMAGE021
的相位进行微分，得到发射信号的瞬时频率函数
Figure 525005DEST_PATH_IMAGE022
为：
Figure DEST_PATH_IMAGE023
；
经过时延
Figure 81888DEST_PATH_IMAGE024
后接收到回波信号的频率函数
Figure DEST_PATH_IMAGE025
为：
Figure 191664DEST_PATH_IMAGE026
；
其中
Figure DEST_PATH_IMAGE027
，C为光速，H为相对高度；
S102，根据回波信号和当前发射的发射信号确定差频信号；
S102具体包括：
利用公式
Figure 108674DEST_PATH_IMAGE010
确定差频信号；
其中，
Figure 43132DEST_PATH_IMAGE011
为差频信号。
S103，对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
S104，对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；
S105，若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；
S105具体包括：
当FFT运算后的差频信号的幅值对应的频率大于设定频率范围时，增大调制信号的调制周期；
当FFT运算后的差频信号的幅值对应的频率小于设定频率范围时，减小调制信号的调制周期；
S106，若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；
S106具体包括：
生成
Figure 974178DEST_PATH_IMAGE028
序列；
M为Z平面上的频谱采样点数，
Figure DEST_PATH_IMAGE029
为任意复数，L为
Figure 565391DEST_PATH_IMAGE028
序列的点数，n为
Figure 849873DEST_PATH_IMAGE028
序列的变量：
Figure 543898DEST_PATH_IMAGE030
；
对序列
Figure 770480DEST_PATH_IMAGE028
进行快速傅里叶变换；
k为FFT的点数，
Figure DEST_PATH_IMAGE031
；
将差频信号
Figure 199056DEST_PATH_IMAGE032
与系数
Figure DEST_PATH_IMAGE033
相乘得到序列
Figure 83966DEST_PATH_IMAGE034
，并对
Figure 789754DEST_PATH_IMAGE034
序列末端补零；
Figure DEST_PATH_IMAGE035
；
其中，a为任意复数；
对序列
Figure 577450DEST_PATH_IMAGE034
进行快速傅里叶变换：
Figure 368689DEST_PATH_IMAGE036
；
将序列
Figure DEST_PATH_IMAGE037
与
Figure 905893DEST_PATH_IMAGE038
进行频域相乘：
Figure DEST_PATH_IMAGE039
；
对
Figure 246613DEST_PATH_IMAGE040
进行IFFT变换：
Figure DEST_PATH_IMAGE041
；
将
Figure 424785DEST_PATH_IMAGE042
和
Figure DEST_PATH_IMAGE043
进行相乘，并在频谱上对M点进行采样，提取差频信号的频率值：
Figure 421428DEST_PATH_IMAGE044
；
S107，若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；
S107具体包括：
当Chirp_Z变换后的差频信号的最大幅值对应的频率大于恒定差频频率时，增大调制信号的调制周期；
当Chirp_Z变换后的差频信号的最大幅值对应的频率小于恒定差频频率时，减小调制信号的调制周期；
S108，若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值。
S108具体包括：
利用公式
Figure 84360DEST_PATH_IMAGE012
确定当前高度值；
其中，
Figure DEST_PATH_IMAGE045
为Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率时的发射信号的调制周期，H为当前高度值，C为光速。
一种基于线性调频连续波体制的高度表测高系统，包括：
回波信号获取模块，用于获取由地面反射的回波信号；所述回波信号为发射天线发射的发射信号到达地面，由地面反射回来的信号；所述发射信号为通过锯齿波调制并随时间作线性变换的调制信号；
差频信号确定模块，用于根据回波信号和当前发射的发射信号确定差频信号；
预处理模块，用于对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
第一判断模块，用于对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；
调制周期第一调整模块，用于若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；
第二判断模块，用于若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；
调制周期第二调整模块，用于若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；
当前高度值确定模块，用于若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值。
图3为一种基于线性调频连续波体制的高度表测高系统结构示意图，如图3所示，本发明所提供的一种基于线性调频连续波体制的高度表测高系统，应用于所述的一种基于线性调频连续波体制的高度表测高方法，包括：接收机、发射机、调制器、混频器、信号预处理装置、信号处理装置以及接口单元；
所述调制器与发射机连接；所述调制器用于生成发射信号，并将所述调制信号发送至发射机；所述发射信号通过锯齿波调制并随时间作线性变换的调制信号；
所述混频器分别与接收机和发射机连接；所述混频器用于根据接收机接收的回波信号和发射机发射的发射信号生成差频信号；
所述信号预处理装置与所述混频器连接；所述信号预处理装置用于对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
所述信号处理装置分别与信号预处理装置、调制器以及接口单元连接；所述信号处理装置用于对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值；并将当前高度值发送至接口单元。
所述信号预处理装置包括：滤波器、放大器以及模数转换器；
所述滤波器、放大器以及模数转换器依次连接。
所述发射机包括：发射天线和功分器；
所述功分器用于将发射信号分成两路，一路至发射天线发射出去，另一路作为本振信号进入混频器。
所述接收机包括：接收天线。
本发明所提供的一种基于线性调频连续波体制的高度表测高系统，还包括：数模转换器；
所述数模转换器分别与调制器和信号处理装置连接。
如图4和图5所示，信号处理装置为FPGA，通过控制外部D/A转换器产生锯齿波调制信号，然后该信号对发射机进行调频，产生一个锯齿波线性调频的连续波，最后发射机对该信号进行功率放大等操作并有发射天线辐射出去，并产生相干本振信号送给接收机。
接收天线将回波信号接收送至混频器与当前的发射信号进行混频，经过滤波放大处理后进行A/D转换并将差频信号送至FPGA内部进行处理。
差频信号进入FPGA内部后，首先在DRAM中进行数据的存储，然后对其进行FFT运算，通过判断最大幅值所对应的频率，控制外部D/A转换调整锯齿波的调制周期。当最大幅值所对应的频率在指定频率范围内时，对差频信号进行Chirp-Z变换。
如图5所示，FPGA通过创建DDS的IP核产生线性调频信号，并放入内部存储器DRAM中，分别给系数1、系数2和h（n）使用。
在FPGA内部差频信号与系数1进行乘法运算产生
Figure 374527DEST_PATH_IMAGE046
信号，然后对其进行FFT运算，h（n）做FFT运算后与y（n）的FFT结果进行相乘，其结果通过IFFT运算产生g（n）信号。
将g（n）和系数2进行相乘，然后在频谱上对M点进行采样，提取差拍信号频率值，判断差拍信号最大幅度值对应的频点是否等于指定的频率值。
当最大幅度值所对应的频点大于指定频率时，则通过控制D/A转换器增大锯齿波调制信号的调制周期，当最大幅度值所对应的频点小于指定频率时，则减小锯齿波调制信号的调制周期，当最大幅度值所对应的频点等于指定频率时，则计算此时锯齿波调制信号的周期并根据公式求出此时的高度值。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种基于线性调频连续波体制的高度表测高方法，其特征在于，包括：
获取由地面反射的回波信号；所述回波信号为发射天线发射的发射信号到达地面，由地面反射回来的信号；所述发射信号为通过锯齿波调制并随时间作线性变换的调制信号；
根据回波信号和当前发射的发射信号确定差频信号；
对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；
若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；
若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；
若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；
若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值。
2.根据权利要求1所述的一种基于线性调频连续波体制的高度表测高方法，其特征在于，所述获取由地面反射的回波信号，具体包括：
利用公式
Figure 674348DEST_PATH_IMAGE001
确定回波信号的瞬时频率；
其中，
Figure 559959DEST_PATH_IMAGE002
为回波信号的瞬时频率，
Figure 522099DEST_PATH_IMAGE003
为载波信号频率，k为调频斜率，
Figure 901259DEST_PATH_IMAGE004
，
Figure 999665DEST_PATH_IMAGE005
为锯齿波的调制周期，
Figure 671824DEST_PATH_IMAGE006
为调频带宽，
Figure 121260DEST_PATH_IMAGE007
为获取回波信号的时间，
Figure 304110DEST_PATH_IMAGE008
为延迟时间，
Figure 257023DEST_PATH_IMAGE009
。
3.根据权利要求2所述的一种基于线性调频连续波体制的高度表测高方法，其特征在于，所述根据回波信号和当前发射的发射信号确定差频信号，具体包括：
利用公式
Figure 100083DEST_PATH_IMAGE010
确定差频信号；
其中，
Figure 36815DEST_PATH_IMAGE011
为差频信号。
4.根据权利要求3所述的一种基于线性调频连续波体制的高度表测高方法，其特征在于，所述若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值，具体包括：
利用公式
Figure 23356DEST_PATH_IMAGE012
确定当前高度值；
其中，
Figure 830775DEST_PATH_IMAGE013
为Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率时的发射信号的调制周期，H为当前高度值，C为光速。
5.一种基于线性调频连续波体制的高度表测高系统，其特征在于，包括：
回波信号获取模块，用于获取由地面反射的回波信号；所述回波信号为发射天线发射的发射信号到达地面，由地面反射回来的信号；所述发射信号为通过锯齿波调制并随时间作线性变换的调制信号；
差频信号确定模块，用于根据回波信号和当前发射的发射信号确定差频信号；
预处理模块，用于对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
第一判断模块，用于对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；
调制周期第一调整模块，用于若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；
第二判断模块，用于若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；
调制周期第二调整模块，用于若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；
当前高度值确定模块，用于若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值。
6.一种基于线性调频连续波体制的高度表测高系统，应用于权利要求1-4任意一项所述的一种基于线性调频连续波体制的高度表测高方法，其特征在于，包括：接收机、发射机、调制器、混频器、信号预处理装置、信号处理装置以及接口单元；
所述调制器与发射机连接；所述调制器用于生成发射信号，并将所述调制信号发送至发射机；所述发射信号通过锯齿波调制并随时间作线性变换的调制信号；
所述混频器分别与接收机和发射机连接；所述混频器用于根据接收机接收的回波信号和发射机发射的发射信号生成差频信号；
所述信号预处理装置与所述混频器连接；所述信号预处理装置用于对所述差频信号进行预处理；所述预处理包括：滤波、放大以及模数转换；
所述信号处理装置分别与信号预处理装置、调制器以及接口单元连接；所述信号处理装置用于对预处理后的差频信号进行FFT运算，判断FFT运算后的差频信号的幅值对应的频率是否处于设定频率范围内；若FFT运算后的差频信号的幅值对应的频率未处于设定频率范围内，则调整当前发射的发射信号的调制周期；若FFT运算后的差频信号的幅值对应的频率处于设定频率范围内，则对FFT运算后的差频信号进行Chirp_Z变换，并判断Chirp_Z变换后的差频信号的最大幅值对应的频率是否等于恒定差频频率；若Chirp_Z变换后的差频信号的最大幅值对应的频率不等于恒定差频频率，则调整当前发射的发射信号的调制周期；若Chirp_Z变换后的差频信号的最大幅值对应的频率等于恒定差频频率，则固定当前发射的发射信号的调制周期，并根据当前发射的发射信号的调制周期确定当前高度值；并将当前高度值发送至接口单元。
7.根据权利要求6所述的一种基于线性调频连续波体制的高度表测高系统，其特征在于，所述信号预处理装置包括：滤波器、放大器以及模数转换器；
所述滤波器、放大器以及模数转换器依次连接。
8.根据权利要求6所述的一种基于线性调频连续波体制的高度表测高系统，其特征在于，所述发射机包括：发射天线和功分器；
所述功分器用于将发射信号分成两路，一路至发射天线发射出去，另一路作为本振信号进入混频器。
9.根据权利要求6所述的一种基于线性调频连续波体制的高度表测高系统，其特征在于，所述接收机包括：接收天线。
10.根据权利要求6所述的一种基于线性调频连续波体制的高度表测高系统，其特征在于，还包括：数模转换器；
所述数模转换器分别与调制器和信号处理装置连接。
CN202111083383.XA 2021-09-16 2021-09-16 一种基于线性调频连续波体制的高度表测高方法及系统 Pending CN113534145A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111083383.XA CN113534145A (zh) 2021-09-16 2021-09-16 一种基于线性调频连续波体制的高度表测高方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111083383.XA CN113534145A (zh) 2021-09-16 2021-09-16 一种基于线性调频连续波体制的高度表测高方法及系统

Publications (1)

Publication Number Publication Date
CN113534145A true CN113534145A (zh) 2021-10-22

Family

ID=78123214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111083383.XA Pending CN113534145A (zh) 2021-09-16 2021-09-16 一种基于线性调频连续波体制的高度表测高方法及系统

Country Status (1)

Country Link
CN (1) CN113534145A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114070352A (zh) * 2021-11-02 2022-02-18 中国电子科技集团公司第二十九研究所 一种用于无线电高度表的调频周期切换方法
CN116087942A (zh) * 2023-04-06 2023-05-09 烟台初心航空科技有限公司 航空高度表调制信号的产生方法

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CH244084A (de) * 1943-10-27 1946-08-31 Patelhold Patentverwertung Verfahren und Einrichtung zur Beobachtung eines Gegenstandes von einer Stelle aus, an der Schwingungen gesendet und empfangen werden.
JPH11183601A (ja) * 1997-12-19 1999-07-09 Fujitsu Ten Ltd レーダ装置の信号処理方法
JP2012068122A (ja) * 2010-09-24 2012-04-05 Casio Comput Co Ltd 高度計
CN102707275A (zh) * 2012-05-25 2012-10-03 北京航空航天大学 一种线性调频连续波雷达高度表的数字化处理方法
US20160223665A1 (en) * 2015-02-04 2016-08-04 Honeywell International Inc. Systems and methods for measuring velocity with a radar altimeter
CN106840310A (zh) * 2017-02-17 2017-06-13 水利部南京水利水文自动化研究所 一种连续调频波制雷达水位计测量方法
EP3208633A1 (en) * 2016-02-19 2017-08-23 Honeywell International Inc. Method and system for fmcw radar altimeter system height measurement resolution improvement
US20180074187A1 (en) * 2015-02-11 2018-03-15 Honeywell International Inc. Velocity and attitude estimation using an interferometric radar altimeter
CN107907878A (zh) * 2017-11-08 2018-04-13 零八电子集团有限公司 高精度获取fmcw雷达距离测量值的方法
CN109061623A (zh) * 2018-06-25 2018-12-21 南京信息工程大学 一种应用于无人机的平面集成式微波测高雷达及测量方法
CN109581369A (zh) * 2019-01-25 2019-04-05 西安电子科技大学 非均匀多通道恒差拍频率体制的雷达高度表
CN111025278A (zh) * 2019-12-31 2020-04-17 北京中科飞鸿科技股份有限公司 一种无线电高度表
CN212903232U (zh) * 2020-10-09 2021-04-06 北京安达维尔航空设备有限公司 一种机载无线电罗盘及高度表模块

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CH244084A (de) * 1943-10-27 1946-08-31 Patelhold Patentverwertung Verfahren und Einrichtung zur Beobachtung eines Gegenstandes von einer Stelle aus, an der Schwingungen gesendet und empfangen werden.
JPH11183601A (ja) * 1997-12-19 1999-07-09 Fujitsu Ten Ltd レーダ装置の信号処理方法
JP2012068122A (ja) * 2010-09-24 2012-04-05 Casio Comput Co Ltd 高度計
CN102707275A (zh) * 2012-05-25 2012-10-03 北京航空航天大学 一种线性调频连续波雷达高度表的数字化处理方法
US20160223665A1 (en) * 2015-02-04 2016-08-04 Honeywell International Inc. Systems and methods for measuring velocity with a radar altimeter
US20180074187A1 (en) * 2015-02-11 2018-03-15 Honeywell International Inc. Velocity and attitude estimation using an interferometric radar altimeter
EP3208633A1 (en) * 2016-02-19 2017-08-23 Honeywell International Inc. Method and system for fmcw radar altimeter system height measurement resolution improvement
CN106840310A (zh) * 2017-02-17 2017-06-13 水利部南京水利水文自动化研究所 一种连续调频波制雷达水位计测量方法
CN107907878A (zh) * 2017-11-08 2018-04-13 零八电子集团有限公司 高精度获取fmcw雷达距离测量值的方法
CN109061623A (zh) * 2018-06-25 2018-12-21 南京信息工程大学 一种应用于无人机的平面集成式微波测高雷达及测量方法
CN109581369A (zh) * 2019-01-25 2019-04-05 西安电子科技大学 非均匀多通道恒差拍频率体制的雷达高度表
CN111025278A (zh) * 2019-12-31 2020-04-17 北京中科飞鸿科技股份有限公司 一种无线电高度表
CN212903232U (zh) * 2020-10-09 2021-04-06 北京安达维尔航空设备有限公司 一种机载无线电罗盘及高度表模块

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
L. KUANG 等: ""A new algorithm for long range FMCW radar altimeter"", 《2013 IEEE CHINA SUMMIT AND INTERNATIONAL CONFERENCE ON SIGNAL AND INFORMATION PROCESSING》 *
廖志佳: ""FMCW雷达高度表信号处理方法研究与实现"", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *
邢慕增: ""直升机无线电测高受外吊挂物影响研究"", 《中国科技信息》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114070352A (zh) * 2021-11-02 2022-02-18 中国电子科技集团公司第二十九研究所 一种用于无线电高度表的调频周期切换方法
CN114070352B (zh) * 2021-11-02 2023-03-14 中国电子科技集团公司第二十九研究所 一种用于无线电高度表的调频周期切换方法
CN116087942A (zh) * 2023-04-06 2023-05-09 烟台初心航空科技有限公司 航空高度表调制信号的产生方法

Similar Documents

Publication Publication Date Title
US7583222B2 (en) Method for using pulse compression in weather radar
EP3144701B1 (en) Method and device for generating non-linear frequency modulation signal
CN113534145A (zh) 一种基于线性调频连续波体制的高度表测高方法及系统
CN103048650B (zh) 基于步进频雷达的回波模拟方法及系统
CN104076352A (zh) 低截获测速方法及雷达装置
KR101289315B1 (ko) 넓은 고도 범위와 낮은 측정 오차를 갖는 fmcw 전파고도계 및 그의 고도 측정방법
CN105301592B (zh) 采用汽车防撞雷达系统的多目标识别算法
Sundaresan et al. Real time implementation of FMCW radar for target detection using GNU radio and USRP
CN106840310B (zh) 一种连续调频波制雷达水位计测量方法
KR101239166B1 (ko) Ｆｍｃｗ 근접 센서
Griffiths et al. Digital generation of high time-bandwidth product linear FM waveforms for radar altimeters
CN109001729B (zh) 太赫兹成像中调频连续波线性度实时校准方法及其系统
CN113093157B (zh) 基于微波光子稳相传输链路的分布式接收阵列通道误差标定方法、系统
CN101266293A (zh) 激光合成孔径雷达成像距离向相位补偿方法
CN109257051A (zh) 一种数字信道化射频存储系统
CN104991231A (zh) 一种中频对消技术
CN113055108A (zh) 一种混频器群时延测量方法与装置
US5187484A (en) Linearizing a swept-frequency radar
CN113433540A (zh) 一种用于fmcw调制的相位差检测设备
JPH04357485A (ja) パルス・ドップラーレーダ装置
CN209946381U (zh) 一种隧道安全距测量系统
US10267666B2 (en) Fill level measurement device
CN106842159A (zh) 一种雷达目标探测中信息量的计算方法
JPH0318784A (ja) Ｆｍ―ｃｗ測距方法
CN114217148B (zh) 基于prf谱线的脉冲调制波相位噪声定标装置及方法

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20211022

RJ01 Rejection of invention patent application after publication