CN106595572A - 一种飞行器低空飞行高度测量方法及装置 - Google Patents
一种飞行器低空飞行高度测量方法及装置 Download PDFInfo
- Publication number
- CN106595572A CN106595572A CN201610916471.6A CN201610916471A CN106595572A CN 106595572 A CN106595572 A CN 106595572A CN 201610916471 A CN201610916471 A CN 201610916471A CN 106595572 A CN106595572 A CN 106595572A
- Authority
- CN
- China
- Prior art keywords
- signal
- sampled point
- peak value
- echo
- vertical coordinate
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000012545 processing Methods 0.000 claims abstract description 158
- 238000005070 sampling Methods 0.000 claims abstract description 118
- 238000005259 measurement Methods 0.000 claims abstract description 33
- 238000001914 filtration Methods 0.000 claims description 43
- 230000005540 biological transmission Effects 0.000 claims description 28
- 238000000691 measurement method Methods 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 14
- 241000196324 Embryophyta Species 0.000 description 29
- 238000010586 diagram Methods 0.000 description 10
- 239000000654 additive Substances 0.000 description 5
- 230000000996 additive effect Effects 0.000 description 5
- 230000003116 impacting effect Effects 0.000 description 4
- 238000003672 processing method Methods 0.000 description 4
- 230000004888 barrier function Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000004080 punching Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- 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
- G01C5/005—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels altimeters for aircraft
-
- 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)
Abstract
本发明实施例提供了一种飞行器低空飞行高度测量方法及装置，所述方法包括：发射信号；接收所述发射信号的回波信号，对所述回波信号进行信号处理，得到回波处理信号；根据所述回波处理信号获取第一峰值采样点和第二峰值采样点，并根据所述第一峰值采样点的横坐标和所述第二峰值采样点的横坐标计算出第一高度距离和第二高度距离。所述装置用于执行上述方法。本发明实施例提供的飞行器低空飞行高度测量方法及装置，根据第一峰值采样点的横坐标计算出飞行器距离植物顶端的高度，并根据第二峰值采样点的横坐标计算出飞行器距离地面的高度作为参考，提高了飞行器低空飞行时飞行高度测量的准确性。
Description
技术领域
本发明实施例涉及飞行高度测量技术领域，具体涉及一种飞行器低空飞行高度测量方法及装置。
背景技术
当飞行器在空中飞行时，通常需要测量飞行器距离地面的距离，尤其当飞行器贴近地面飞行时，若地面存在草本植物，为了安全飞行，需要较为准确的测量飞行器与植物的距离如：使用无人机在农田林地中工作时，就需要较为准确的测量无人机与植物的距离。
现有技术中，飞行器飞行高度测量一般采用气压测高方法，其原理为：测量飞行器当前位置气压与起飞地面位置的气压之差，根据计算得出的气压差计算出飞行器当前位置的高度。这种测量方法实现简单、成本较低，但在低空测量时误差较大，不能准确的测量出飞行器与植物之间的距离。
因此，如何提出一种方法，能够提高飞行器低空飞行时飞行高度测量的准确性，成为亟待解决的问题。
发明内容
针对现有技术中的缺陷，本发明实施例提供一种飞行器低空飞行高度测量方法。
一方面本发明实施例提供一种飞行器低空飞行高度测量方法，包括：
发射信号；
接收所述发射信号的回波信号，对所述回波信号进行信号处理，得到回波处理信号；
根据所述回波处理信号获取第一峰值采样点和第二峰值采样点，并根据所述第一峰值采样点的横坐标和所述第二峰值采样点的横坐标计算出第一高度距离和第二高度距离。
另一方面本发明实施例提供一种飞行器低空飞行高度测量装置，包括：
信号发射单元，用于发射信号；
信号处理单元，用于接收所述发射信号的回波信号，对所述回波信号进行信号处理，得到回波处理信号；
飞行高度测量单元，用于根据所述回波处理信号获取第一峰值采样点和第二峰值采样点，并根据所述第一峰值采样点的横坐标和所述第二峰值采样点的横坐标计算出第一高度距离和第二高度距离。
本发明实施例提供的飞行器低空飞行高度测量方法及装置，通过使用冲击雷达发射信号，并由接收机接收发射信号的回波信号，对回波信号进行处理后，找出回波处理信号中的第一峰值采样点和第二峰值采样点。根据第一峰值采样点的横坐标计算出飞行器距离植物顶端的高度，并根据第二峰值采样点的横坐标计算出飞行器距离地面的高度作为参考，提高了飞行器低空飞行时飞行高度测量的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中飞行器低空飞行高度测量方法的流程示意图；
图2为本发明实施例中冲击雷达发射的波形示意图；
图3为本发明实施例中的原始回波信号波形图；
图4为本发明实施例中均值去噪信号处理后的波形图；
图5为本发明实施例中中值滤波信号处理后的波形图；
图6为本发明实施例中匹配滤波信号处理后的波形图；
图7为本发明实施例中方差滤波信号处理后的波形图；
图8为本发明实施例中飞行器低空飞行高度测量方法的整体流程示意图；
图9为本发明实施例中第一峰值采样点和第二峰值采样点在原始回波信号波形图中对应位置的示意图；
图10为本发明实施例中第一峰值采样点和第二峰值采样点在方差滤波信号处理后的波形图中对应位置的示意图；
图11为本发明实施例中飞行器低空飞行高度测量装置结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例中飞行器低空飞行高度测量方法的流程示意图，如图1所示，本发明实施例提供的一种飞行器低空飞行高度测量方法包括以下步骤：
S1、发射信号。
具体地，本发明实施例中的飞行器低空飞行高度测量方法采用冲击雷达向地面发射信号，在飞行器上安装冲击雷达，图2为本发明实施例中冲击雷达发射的波形示意图。冲击雷达是一种超宽带雷达，其脉宽极窄，一般都在纳秒或皮秒量级，并且具有高距离分辨率和穿透能力强的优点。其中图2中的波形示意图是冲击雷达的一阶微分高斯脉冲发射信号的波形示意图。
S2、接收所述发射信号的回波信号，对所述回波信号进行信号处理，得到回波处理信号。
具体地，当冲击雷达向地面发射信号后，信号在介质中传播，遇到障碍物后反射，再经过介质传播后，由冲击雷达中的接收机接收到回波信号。由于接收到的回波信号会有杂波，对计算飞行器的飞行高度有影响，因此需要对接收机接收到的回波信号进行信号处理，得到回波处理信号，为了便于观察可以根据回波处理信号绘制回波处理信号波形图。具体可以将接收到的回波数据传输到计算机，由计算机软件进行信号处理并绘制回波处理信号波形图，当然也可以采用其他方法，本发明实施例不作具体限定。
S3、根据所述回波处理信号获取第一峰值采样点和第二峰值采样点，并根据所述第一峰值采样点的横坐标和所述第二峰值采样点的横坐标计算出第一高度距离和第二高度距离。
具体地，将回波信号处理好后，根据回波处理信号获取第一峰值采样点和第二峰值采样点。由于地面上的植物和地面对冲击雷达发射的信号会有集中的反射，反应在回波处理信号波形图上就会出现明显的波峰，其中植物顶端的回波是接收机最先收到的回波信号，因此第一峰值采样点对应植物顶端的回波，相应的第二峰值采样点对应地面的回波。再根据回波处理信号中的第一峰值采样点在所述回波处理信号波形中的横坐标计算出第一高度距离，根据第二峰值采样点在所述回波处理信号波形中的横坐标计算出第二高度距离，也就是飞行器距离植物顶端的高度和飞行器距离地面的高度。需要说明的是，其中第一高度距离是指飞行器距离最近目标物的距离，不仅仅指距离植物顶端的距离。
具体根据回波处理信号波形图中的第一峰值采样点的横坐标和第二峰值采样点的横坐标就计算第一高度距离和第二高度距离的方法如下：
假设冲击雷达的发射信号为s(t)，冲击雷达接收机接收到的回波信号如公式(1)所示：
x(t)＝s(t)*h(t)*δ(t-τ)+n(t) (1)
式中：h(t)为目标的冲击响应，*表示卷积运算，n(t)为加性噪声和杂波，τ为接收回波信号时刻与发射信号时刻之间的延时。
接收到回波信号后要对回波信号作信号处理，尽可能的去除加性噪声和杂波n(t)，保留发射信号由目标反射后的延迟回波部分，回波处理信号如公式(2)所示：
x’(t)＝x(t)-n(t)＝s(t)*h(t)*δ(t-τ) (2)
根据公式(2)分别计算出植物顶端的信号延时和地面的信号延时后，就可以根据公式(3)计算出第一高度距离和第二高度距离。由于在实际应用时获得的回波信号不是连续的，是由许多采样点组成的，因此具体计算延时τ可以通过回波处理信号波形图中采样点的横坐标乘以采样时间间隔得到，采样时间间隔是已知的，只要能够获取到回波处理信号波形图中采样点的横坐标，就可以计算出延时τ，进一步可以计算出第一高度距离和第二高度距离。公式(3)如下所示：
式中：v为发射信号在介质中的传播速度。将植物顶端的信号延时和地面的信号延时分别带入公式(3)就可以计算出飞行器距离植物顶端的高度和飞行器距离地面的高度，即第一高度距离和第二高度距离。
本发明实施例提供的飞行器低空飞行高度测量方法，通过使用冲击雷达发射信号，并由接收机接收发射信号的回波信号，对回波信号进行处理后，找出回波处理信号中的第一峰值采样点和第二峰值采样点。根据第一峰值采样点的横坐标计算出飞行器距离植物顶端的高度，并根据第二峰值采样点的横坐标计算出飞行器距离地面的高度作为参考，提高了飞行器低空飞行时飞行高度测量的准确性。
在上述实施例的基础上，所述根据所述回波处理信号波形图获取第一峰值采样点包括：预先设置门限值，从所述回波处理信号波形图中的第一个采样点开始依次将所有采样点的纵坐标的绝对值与所述门限值进行比对，获取第一个纵坐标的绝对值大于所述门限值的采样点作为所述第一峰值采样点。
具体地，根据实际情况预先设置门限值，从经过信号处理后的回波处理信号波形图中的第一个采样点开始，判断采样点的纵坐标的绝对值是否超过门限值，若超过则将所述采样点作为第一峰值采样点，并记录该采样点的横坐标，即记录该采样点是第几个采样点；若不超过则继续判断下一采样点直至寻找到纵坐标的绝对值超过门限值的采样点。
本发明实施例提供的飞行器低空飞行高度测量方法，通过合理设置门限值来获取第一峰值采样点，从而更加准确的计算出第一高度距离，即飞行器距离最近目标的距离。
在上述实施例的基础上，所述预先设置门限值包括：获取所述回波处理信号波形图中的最大幅值，将所述最大幅值乘以比例系数所得的值作为所述门限值，其中所述比例系数大于0小于1。
具体地，首先获取回波处理信号波形图中的最大幅值x_max，将x_max乘以一个比例系数i，得到的值作为门限值，即门限值为r＝i×x_max，其中比例系数0<i<1。比例系数的具体设置可以根据实际情况而定，本发明实施例中i取0.3，具体还可以根据环境和对精度的要求作出调整，本发明实施例不作具体限定。
本发明实施例提供的飞行器低空飞行高度测量方法，通过获取回波处理信号的最大幅值，再乘以比例系数来设置门限值，使得门限值可以根据不同的回波信号设置，提高了第一峰值采样点获取的准确性，从而更加准确的计算出第一高度距离。
在上述实施例的基础上，所述对所述回波信号进行信号处理，包括：对所述回波信号进行均值去噪信号处理或中值滤波信号处理或匹配滤波信号处理。
具体地，本发明实施例在采用预先设置门限的方法获取第一峰值采样点在前，对回波信号进行信号处理可以是均值去噪信号处理或中值滤波信号处理或匹配滤波信号处理。需要说明的是，具体对回波信号进行信号处理的方法可以是其中一种或多种结合，还可以是其他的信号处理方法，本发明实施例不做具体限定。
下面具体介绍几种信号处理方法，本发明实施例是将几种信号处理结合到一起依次对接收到的回波信号作均值去噪信号处理、中值滤波信号处理以及匹配滤波信号处理，得到回波处理信号。
由于在实际的应用中无法获得连续的回波信号x[t]，而是由不连续的采样点组成的x[nT]，T为采样时间间隔，一般记作x[n]，以下都用x[n]表示回波信号，图3为本发明实施例中的原始回波信号波形图。
(1)均值去噪信号处理：对接收到的回波信号x[n]进行均值去噪，这种去噪方法适用在飞行器静止或飞行速度较低且环境变化不大的情况。此时冲击雷达可获得连续连续N组相似的回波信号，他们的差别主要体现在信号中加性随机噪声。通过多组数据取平均的方法，可将信噪比改善约10logNdB，经过均值去噪信号处理后得到的信号为x1[n]，具体如公式(4)所示：
式中：x0i[n]为未处理的N组相似信号中第i组信号。
图4为本发明实施例中均值去噪信号处理后的波形图。
(2)中值滤波信号处理：由于回波信号中可能叠加有杂波，影响计算结果，而相较于频率很高的目标反射冲击雷达信号，杂波信号的频率较低，可以通过中值滤波的方法将变化较快的高频部分滤除，保留下低频杂波部分，再将其减去，就实现了去除杂波的目的。对经过均值去噪信号处理后得到的信号为x1[n]进行中值滤波信号处理，杂波信号表示为x2[n]，具体如公式(5)所示：
经过中值滤波信号处理后得到的信号为x3[n]，具体如公式(6)所示：
图5为本发明实施例中中值滤波信号处理后的波形图。
(3)匹配滤波信号处理：将中值滤波信号处理后得到的信号为x3[n]通过与冲击雷达的发射信号做互相关得到匹配滤波信号，经过匹配滤波信号处理后得到的信号为x4[n]，具体如公式(7)所示：
式中：s[n]为发射的一阶微分高斯脉冲的离散采样，即s[n]＝s[nT]，T为采样时间间隔。
图6为本发明实施例中匹配滤波信号处理后的波形图。
本发明实施例提供的飞行器低空飞行高度测量方法，通过对接收到的回波信号进行信号处理，得到更加准确的回波处理信号，从而提高了飞行器飞行高度测量的准确性，保证了飞行器飞行平稳和安全。
在上述实施例的基础上，所述根据所述回波处理信号获取第一峰值采样点和第二峰值采样点包括：
获取所述回波处理信号中的多个初始峰值采样点和多个波谷采样点，所述初始峰值采样点的纵坐标大于其前后采样点的纵坐标，所述波谷采样点是相邻所述初始峰值采样点之间的纵坐标最小的采样点，其中将第一个所述初始峰值采样点之前的第一个所述波谷采样点的纵坐标设置为零；
将纵坐标大于其之后某一个采样点纵坐标的第一预设倍数且大于其之前相邻的波谷采样点纵坐标的第二预设倍数的初始峰值采样点作为峰值采样点；
获取所有峰值采样点中纵坐标最大的采样点作为第一峰值采样点，纵坐标次最大的采样点作为第二峰值采样点。
具体地，从回波处理信号中的第一个采样点开始依次将所有采样点的纵坐标进行比对，将每一个采样点的纵坐标与其前后采样点的纵坐标作比对，找出纵坐标大于前后采样点的纵坐标的采样点，即为初始峰值采样点，每两个初始峰值采样点之间纵坐标最小采样点作为波谷采样点，并且第一个初始分值采样点前的波谷采样点的纵坐标设置为0。在依次将采样点的纵坐标进行比对时，若判断某个初始峰值采样点的纵坐标大于其之后的某个采样点的纵坐标的第一预设倍数，并且该初始峰值采样点的纵坐标大于与其相邻且在其前面的一个波谷采样点纵坐标的第二预设倍数时，则将该初始峰值采样点作为峰值采样点。依次对所有的采样点进行纵坐标的比较判断，找出满足条件的所有的峰值采样点，对找出的峰值采样点的纵坐标进行比对，将纵坐标最大的峰值采样点作为第一峰值采样点，将纵坐标第二大的峰值采样点作为第二峰值采样点。其中第一预设倍数和第二预设倍数的设置可以根据实际使用情况进行设置，可以相等也可以不相等，本发明实施例不作具体限定。
本发明实施例提供的飞行器低空飞行高度测量方法，通过对回波处理信号波形图中的采样点的纵坐标进行比对，并设置第一预设倍数和第二预设倍数，使得找到的第一峰值采样点和第二峰值采样点之后的采样点的能量都回落到比较小的值，从而获得更加准确的第一峰值采样点和第二峰值采样点，进一步使得第一高度距离和第二高度距离更加准确，提高了飞行器飞行高度计算的准确性。
在上述实施例的基础上，所述对所述回波信号进行信号处理，包括：对所述回波信号进行方差滤波信号处理。
具体地，本发明实施例在采用设置第一预设倍数和第二预设倍数找到满足初始峰值采样点之后的采样点的能量都回落到比较小的值的采样点来获取第一峰值采样点和第二峰值采样点的方法之前，为了更好的获得回波信号的能量值，还对回波信号进行了方差滤波信号处理。本发明实施例在进行方差处理之前还依次进行了均值去噪信号处理、中值滤波信号处理和匹配滤波信号处理，具体的处理方法同上述实施例一致，此处不再赘述。方差滤波处理的具体方法如下：
对经过匹配滤波信号处理后得到的信号x4[n]进行方差滤波信号处理，经过方差滤波信号处理后得到的信号为y[n]，即x4[n]及其左右各W个采样点数据求方差保存在y[n]中，x4[n]中发射信号较大的回波能量对应y[n]中的较大数值，具体如公式(8)所示：
式中:W为滤波的宽度，为x4[n-W]到x4[n+W]的均值。
图7为本发明实施例中方差滤波信号处理后的波形图。
在上述实施例的基础上，若判断获知所述第二峰值采样点的纵坐标和所述第一峰值采样点的纵坐标的比值小于预设比值，则令所述第一高度距离与所述第二高度距离相等。
具体地，获取到第一峰值采样点的纵坐标和第二峰值采样点的纵坐标，计算第二峰值采样点的纵坐标与第一峰值采样点的纵坐标的比值，判断所述比值是否小于预设比值，若小于，则认为地面上没有植物或植物高度可以忽略，令第一高度距离与第二高度距离相等。可以理解的是，由于本发明实施例中第一峰值采样点对应的是地面上的植物或障碍物对冲击雷达发射信号进行反射得到的回波信号中的峰值，第二峰值采样点对应的是地面对冲击雷达发射信号进行反射得到的回波信号中的峰值，第二峰值采样点要比第一峰值采样点对应的回波信号消耗的能量多，因此第二峰值采样点的纵坐标一般比第一峰值采样点的纵坐标小，所以预设比值一般大于0小于1。
图8为本发明实施例中飞行器低空飞行高度测量方法的整体流程示意图，下面结合图8具体介绍本发明实施例中的飞行器低空飞行高度测量方法，以便更好的理解本发明的技术方案，如图8所示，本发明实施例中飞行器低空飞行高度测量方法的具体流程如下：
R1、接收回波信号。飞行器上的冲击雷达向地面发射信号，经地面上的植物或障碍物反射后，冲击雷达的接收机会接收到相应的回波信号x[n]。
R2、均值去噪信号处理。对接收到的回波信号进行均值去噪信号处理，得到信号x1[n]。
R3、中值滤波信号处理。对均值去噪信号处理后得到的信号x1[n]进行中值滤波信号处理得到杂波信号x2[n]，再去除杂波，得到信号x3[n]＝x1[n]-x2[n]。
R4、匹配滤波信号处理。对中值滤波信号处理得到的信号x3[n]进行匹配滤波信号处理得到信号x4[n]，再进行步骤R5和R11。
R5、获取最大幅值。匹配滤波信号处理得到信号x4[n]后，获取信号x4[n]中的最大幅值，也就是采样点纵坐标的最大绝对值x_max。
R6、设置门限值。获取到最大幅值x_max后，设置门限值r＝i×x_max，其中i为比例系数，并且0<i<1。
R7、k＝1。设置循环次数k＝1，即从第一个采样点开始搜索。
R8、判断│x4[k]│>r。判断采样点的纵坐标的绝对值是否大于门限值，若不大于门限值执行步骤R9，否则执行步骤R10。
R9、k＝k+1。若采样点的纵坐标的绝对值不超过门限值，则将循环次数k累加1继续对下一个采样点进行判断，直至寻找到超过门限的采样点。
R10、计算第一高度距离。找到第一个纵坐标的绝对值大于门限值的采样点，并记录该采样点的横坐标kp，即记录该采样点是第几个采样点，根据公式(9)计算得到第一高度距离，即距离飞行器最近的目标距离飞行器的距离。
式中：hp为第一高度距离，kp为纵坐标绝对值大于门限值的采样点的横坐标，T为回波信号中采样点的采样时间间隔，v为冲击雷达发射的电磁波在介质中的传播速度，一般可以直接按光速计算。
R11、方差滤波信号处理。对步骤R4匹配滤波信号处理得到的信号x4[n]进行方差滤波信号处理得到信号y[n]。
R12、初始化。方差滤波信号处理后得到信号y[n]，开始初始化，令循环次数k＝1，当前散射点回波最大值y_max＝0，当前散射点回波的最大值对应的位置y_max_index＝0，最大波峰大小h[1]＝0，最大波峰位置h_index[1]＝0，第二大波峰大小h[2]＝0，第二大波峰位置h_index[2]＝0，两波峰间的最小值，即波谷采样点的纵坐标y_min＝0。
R13、判断y[k]>y_max且y[k]>4×y_min。判断当前采样点的纵坐标是否大于y_max并且大于y_min的4倍，若判断结果为否执行步骤R14，否则执行步骤R16。y_min对应的是波谷采样点的纵坐标，这里第二预设倍数为4，在获取第一峰值采样点时y_min＝0，此时的第二预设倍数没有作用，在获取第二峰值采样点时，需要先获取波谷采样点，并判断第二峰值采样点的纵坐标是否大于波谷采样点的纵坐标的第二预设倍数。
R14、判断y[k]<y_min。若判断结果为是，则执行步骤R15，否则执行步骤R17。
R15、y_min＝y[k]。将y[k]的值赋给y_min。
R16、y_max＝y[k]，y_max_index＝k。若判断获知y[k]>y_max且y[k]>4×y_min，将y[k]的值赋给y_max，k的值值赋给y_max_index。
R17、k＝k+1。循环次数k累加1，继续对下一个采样点进行比对判断。
R18、判断k>M。M为回波信号中数据的总采样点数，判断当前循环次数是否大于总采样点数M，即判断是否达到最后一个采样点，是则跳出循环执行R22，否则执行R19。
R19、判断y[k]<0.25×y_max。判断当前采样点是否小于当前y_max的0.25倍，即此时的散射点回波最大值y_max是否大于当前采样点纵坐标y[k]的4倍，即第一预设倍数为4，是则认为当前的回波能量已经回落到了相比前一波峰值较低的大小，认为y_max是第一峰值采样点或第二峰值采样点，执行R20，否则执行R13，进行下一次循环。可以理解的是，本发明实施例中第一预设倍数和第二预设倍数相等，都是4，实际应用时可以根据需要取其他值，本发明实施例不作具体限定。
R20、将h[1]、h[2]和y_max中最大的值存入h[1]，第二大的值存入h[2]，对应的坐标分别存入h_index[1]和h_index[2]。将最近获得的极值y_max与之前最大的两个波峰值h[1]、h[2]进行比较，保留两个最大的两个值分别存入h[1]、h[2]，并且h[1]>h[2]，并将对应的横坐标存入h_index[1]和h_index[2]，这样就获取到第一峰值采样点h[1]，以及对应的横坐标h_index[1]，并在之后的循环中获取到第二峰值采样点h[2]，以及对应的横坐标h_index[2]。
R21、y_max＝0，y_max_index＝0，y_min＝y[k]。继续寻找下一强散射点回波能量的峰值，即第二峰值采样点，此时将y_max和y_max_index重新置为零，将y[k]的值赋给y_min，寻找两个峰值采样点之间的波谷采样点，即y_min＝y[k]，返回步骤R13，重新循环，以获取第二峰值采样点。
R22、判断h[2]/h[1]<q。判断两个波峰纵坐标的比值即h[2]/h[1]是否小于预设比值q，其中0<q<1，若判断结果为是，执行步骤23，否则否执行步骤R24。
R23、忽略第一高度距离。若判断两个波峰纵坐标的比值小于预设比值，则认为回波信号中仅有唯一强散射点，可以认为没有植物或植物高度可忽略，即第一高度距离与第二高度距离相等，忽略通过门限值计算出的第一高度距离，以步骤R24计算出的第二高度距离为飞行器飞行高的的参考。
R24、计算第二高度距离。若判断获知h[2]/h[1]<q不成立，因为地面距飞行器比植物距飞行器更远，对应回波信号中横坐标更大，认为h_index[1]、h_index[2]中较小的为植物散射回波能量最大时刻的横坐标，即第一峰值采样点的横坐标，较大的为地面散射回波能量最大时刻的横坐标，即第二峰值采样点的横坐标。利用h_index[1]、h_index[2]中较大的值kg根据公式(10)计算第二高度距离，也就是地面距飞行器的距离。
式中：hg为第二高度距离，kg为h_index[1]、h_index[2]中较大的值，T为回波信号中采样点的采样时间间隔，v为冲击雷达发射的电磁波在介质中的传播速度，一般可以直接按光速计算。
需要说明的是，本发明实施例中是通过门限值来计算第一高度距离，实际应用时还可以在步骤R20获取到第一峰值采样点时，根据其横坐标采用计算第二峰值采样点的方法计算第一高度距离。也就是说本发明实施例中实际上根据两种方法获取到了两个第一峰值采样点，具体计算第一高度距离时可以采用其中任何一个第一峰值采样点的纵坐标来计算。
图9为本发明实施例中第一峰值采样点和第二峰值采样点在原始回波信号波形图中对应位置的示意图，本发明实施例中的比例系数i＝0.3，最大幅值x_max＝1.238，门限值r＝0.3*1.238＝0.3714，对应的循环次数为采样点数M＝512，采样间隔T＝26ps，信号在介质中的传播速度按光速算v＝300000000m/s，预设比值q＝0.15，采用第一预设倍数和第二预设倍数计算出的h[1]＝0.3881，h[2]＝0.1577，因此h[2]/h[1]>q，认为植物回波存在，即第一高度距离存在。如图9所示，本发明实施例中最终计算出的第一峰值采样点的横坐标为kp＝118，第二峰值采样点的横坐标为kg＝203，计算出的第一高度距离为0.4602米，第二高度距离为0.7917米，需要说明的是，图9中的第一峰值采样点是通过门限值计算得到的，第二峰值采样点是通过对经过方差滤波后的波形图中的采样点的纵坐标进行比较计算后得到的。图10为本发明实施例中第一峰值采样点和第二峰值采样点在方差滤波信号处理后的波形图中对应位置的示意图，通过对采样点的坐标进行依次循环比较得到的第一峰值采样点和第二峰值采样点在方差滤波后的波形中对应的位置如图10所示。需要说明的是，从图9和图10可以看出本发明实施例提供的两种获取第一峰值采样点的方法获取到的第一峰值采样点不完全一样，但是最终计算出的第一高度距离的差别不大，具体应用时可以采用其中任一方法获取到第一峰值采样点计算第一高度距离，本发明实施例不作具体限定。
本发明实施例提供的飞行器低空飞行高度测量方法，通过使用冲击雷达发射信号，并由接收机接收发射信号的回波信号，对回波信号进行处理后，找出回波处理信号波形图中的第一峰值采样点和第二峰值采样点，即对应地面上的障碍物或植物对冲击雷达发射信号的回波能量峰值和地面对冲击雷达发射信号的回波能量峰值。根据第一峰值采样点的横坐标计算出飞行器距离植物顶端的高度，并根据第二峰值采样点的横坐标计算出飞行器距离地面的高度作为参考，提高了飞行器低空飞行时飞行高度测量的准确性。
图11为本发明实施例中飞行器低空飞行高度测量装置结构示意图，如图11所示，在上述实施例的基础上，本发明提供的飞行高度测量装置包括：信号发射单元01，信号处理单元02，飞行高度测量单元03，其中：
信号发射单元01用于发射信号；信号处理单元02用于接收所述发射信号的回波信号，对所述回波信号进行信号处理，得到回波处理信号；飞行高度测量单元03用于根据所述回波处理信号获取第一峰值采样点和第二峰值采样点，并根据所述第一峰值采样点的横坐标和所述第二峰值采样点的横坐标计算出第一高度距离和第二高度距离。
具体地，本发明实施例提供的飞行器低空飞行高度测量装置由信号发射单元01也就是冲击雷达发射信号，由信号处理单元02接收信号发射单元01的发射信号对应的回波信号，并对接收到的回波信号进行信号处理得到回波处理信号，为了便于观察可以根据回波处理信号绘制回波处理信号波形图，具体可以将数据传输到计算机，由计算机进行信号处理以及绘制回波处理信号波形图，当然还可以使用其他方法，本发明实施例不作具体限定。飞行高度测量单元03根据信号处理单元02绘制的回波处理信号获取第一峰值采样点和第二峰值采样点，并根据第一峰值采样点的横坐标和第二峰值采样点的横坐标计算出第一高度距离和第二高度距离，也就是飞行器距离植物顶端的高度和飞行器距离地面的高度。
本发明实施例提供的飞行器低空飞行高度测量装置，通过使用信号发射单元发射信号，由信号处理单元接收发射信号的回波信号，并对回波信号进行处理，由飞行高度测量单元找出回波处理信号中的第一峰值采样点和第二峰值采样点，并根据第一峰值采样点的横坐标计算出飞行器距离植物顶端的高度，根据第二峰值采样点的横坐标计算出飞行器距离地面的高度作为参考，提高了飞行器低空飞行时飞行高度测量的准确性。
在上述实施例的基础上，所述飞行高度测量单元具体用于：预先设置门限值，从所述回波处理信号波形图中的第一个采样点开始依次将所有采样点的纵坐标的绝对值与所述门限值进行比对，获取第一个纵坐标的绝对值大于所述门限值的采样点作为所述第一峰值采样点。
具体地，飞行高度测量单元03根据实际情况预先设置门限值，从信号处理单元02绘制的回波处理信号波形图中的第一个采样点开始，判断采样点的纵坐标的绝对值是否超过门限值，若超过则将所述采样点作为第一峰值采样点，并记录该采样点的横坐标；若不超过则继续判断下一采样点直至寻找到纵坐标的绝对值超过门限值的采样点。需要说明的是，本发明实施例为了更好的获取到第一峰值采样点，在对回波信号进行信号处理时，没有进行方差滤波信号处理。
在上述实施例的基础上，所述飞行高度测量单元还用于：获取所述回波处理信号波形图中的最大幅值，将所述最大幅值乘以比例系数所得的值作为所述门限值，其中所述比例系数大于0小于1。
具体地，飞行高度测量单元03获取回波处理信号波形图中的最大幅值，并将获取到的最大幅值乘以一个比例系数i，得到的值作为门限值，其中比例系数0<i<1。其中比例系数的具体取值同上述实施例一致，此处不再赘述。
在上述实施例的基础上，所述信号处理单元具体用于：对所述回波信号进行均值去噪信号处理或中值滤波信号处理或匹配滤波信号处理。
具体地，信号处理单元02对接收到的回波信号进行信号处理的方法包括均值去噪信号处理或中值滤波信号处理或匹配滤波信号处理。具体进行均值去噪信号处理、中值滤波信号处理以及匹配滤波信号处理的方法同上述实施例一致，此处不再赘述。需要说明的是，具体对回波信号进行信号处理的方法可以是其中一种或多种结合，还可以是其他的信号处理方法，本发明实施例不做具体限定。
在上述实施例的基础上，所述飞行高度测量单元具体用于：获取所述回波处理信号中的多个初始峰值采样点和多个波谷采样点，所述初始峰值采样点的纵坐标大于其前后采样点的纵坐标，所述波谷采样点是相邻所述初始峰值采样点之间的纵坐标最小的采样点，其中将第一个所述初始峰值采样点之前的第一个所述波谷采样点的纵坐标设置为零；
将纵坐标大于其之后某一个采样点纵坐标的第一预设倍数且大于其之前相邻的波谷采样点纵坐标的第二预设倍数的初始峰值采样点作为峰值采样点；
获取所有峰值采样点中纵坐标最大的采样点作为第一峰值采样点，纵坐标次最大的采样点作为第二峰值采样点。
具体地，飞行高度测量单元03从回波处理信号中的第一个采样点开始依次将所有采样点的纵坐标进行比对，将每一个采样点的纵坐标与其前后采样点的纵坐标作比对，找出纵坐标大于前后采样点的纵坐标的采样点，即为初始峰值采样点，每两个初始峰值采样点之间纵坐标最小采样点作为波谷采样点，并且第一个初始分值采样点前的波谷采样点的纵坐标设置为0。在依次将采样点的纵坐标进行比对时，若判断某个初始峰值采样点的纵坐标大于其之后的某个采样点的纵坐标的第一预设倍数，并且该初始峰值采样点的纵坐标大于与其相邻且在其前面的一个波谷采样点纵坐标的第二预设倍数时，则将该初始峰值采样点作为峰值采样点。依次对所有的采样点进行纵坐标的比较判断，找出满足条件的所有的峰值采样点，对找出的峰值采样点的纵坐标进行比对，将纵坐标最大的峰值采样点作为第一峰值采样点，将纵坐标第二大的峰值采样点作为第二峰值采样点。其中第一预设倍数和第二预设倍数的设置与上述实施例一致，此处不再赘述。
在上述实施例的基础上，所述信号处理单元具体用于：对所述回波信号进行方差滤波信号处理。
具体地，在获取第二峰值采样点之前，信号处理单元02还对回波信号进行了方差滤波信号处理，具体进行方差滤波信号处理的方法同上述实施例一致，此处不再赘述。
在上述实施例的基础上，所述飞行高度测量单元还用于：若判断获知所述第二峰值采样点的纵坐标和所述第一峰值采样点的纵坐标的比值小于预设比值，则令所述第一高度距离与所述第二高度距离相等。
具体地，飞行高度测量单元03获取到第一峰值采样点的纵坐标和第二峰值采样点的纵坐标，计算第二峰值采样点的纵坐标与第一峰值采样点的纵坐标的比值，判断所述比值是否小于预设比值，若小于，则认为地面上没有植物或植物高度可以忽略，令第一高度距离与第二高度距离相等。其中预设比值的设置和上述实施例一致，此处不再赘述。
本发明实施例提供的装置用于执行上述方法，其具体的实施方式与方法的实施方式一致，此处不再赘述。
本发明实施例提供的飞行器低空飞行高度测量装置，通过使用信号发射单元发射信号，由信号处理单元接收发射信号的回波信号，并对回波信号进行处理，由飞行高度测量单元找出回波处理信号波形图中的第一峰值采样点和第二峰值采样点，并根据第一峰值采样点的横坐标计算出飞行器距离植物顶端的高度，根据第二峰值采样点的横坐标计算出飞行器距离地面的高度作为参考，提高了飞行器低空飞行时飞行高度测量的准确性。最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (14)
1.一种飞行器低空飞行高度测量方法，其特征在于，包括：
发射信号；
接收所述发射信号的回波信号，对所述回波信号进行信号处理，得到回波处理信号；
根据所述回波处理信号获取第一峰值采样点和第二峰值采样点，并根据所述第一峰值采样点的横坐标和所述第二峰值采样点的横坐标计算出第一高度距离和第二高度距离。
2.根据权利要求1所述的方法，其特征在于，所述根据所述回波处理信号获取第一峰值采样点包括：预先设置门限值，从所述回波处理信号中的第一个采样点开始依次将所有采样点的纵坐标的绝对值与所述门限值进行比对，获取第一个纵坐标的绝对值大于所述门限值的采样点作为所述第一峰值采样点。
3.根据权利要求2所述的方法，其特征在于，所述预先设置门限值包括：获取所述回波处理信号中的最大幅值，将所述最大幅值乘以比例系数所得的值作为所述门限值，其中所述比例系数大于0小于1。
4.根据权利要求2所述的方法，其特征在于，所述对所述回波信号进行信号处理，包括：对所述回波信号进行均值去噪信号处理或中值滤波信号处理或匹配滤波信号处理。
5.根据权利要求1所述的方法，其特征在于，所述根据所述回波处理信号获取第一峰值采样点和第二峰值采样点包括：获取所述回波处理信号中的多个初始峰值采样点和多个波谷采样点，所述初始峰值采样点的纵坐标大于其前后采样点的纵坐标，所述波谷采样点是相邻所述初始峰值采样点之间的纵坐标最小的采样点，其中将第一个所述初始峰值采样点之前的第一个所述波谷采样点的纵坐标设置为零；
将纵坐标大于其之后某一个采样点纵坐标的第一预设倍数且大于其之前相邻的波谷采样点纵坐标的第二预设倍数的初始峰值采样点作为峰值采样点；
获取所有峰值采样点中纵坐标最大的采样点作为第一峰值采样点，纵坐标次最大的采样点作为第二峰值采样点。
6.根据权利要求5所述的方法，其特征在于，所述对所述回波信号进行信号处理，包括：对所述回波信号进行方差滤波信号处理。
7.根据权利要求1所述的方法，其特征在于，若判断获知所述第二峰值采样点的纵坐标和所述第一峰值采样点的纵坐标的比值小于预设比值，则令所述第一高度距离与所述第二高度距离相等。
8.一种飞行器低空飞行高度测量装置，其特征在于，包括：
信号发射单元，用于发射信号；
信号处理单元，用于接收所述发射信号的回波信号，对所述回波信号进行信号处理，得到回波处理信号；
飞行高度测量单元，用于根据所述回波处理信号获取第一峰值采样点和第二峰值采样点，并根据所述第一峰值采样点的横坐标和所述第二峰值采样点的横坐标计算出第一高度距离和第二高度距离。
9.根据权利要求8所述的装置，其特征在于，所述飞行高度测量单元具体用于：预先设置门限值，从所述回波处理信号波形图中的第一个采样点开始依次将所有采样点的纵坐标的绝对值与所述门限值进行比对，获取第一个纵坐标的绝对值大于所述门限值的采样点作为所述第一峰值采样点。
10.根据权利要求9所述的装置，其特征在于，所述飞行高度测量单元还用于：获取所述回波处理信号中的最大幅值，将所述最大幅值乘以比例系数所得的值作为所述门限值，其中所述比例系数大于0小于1。
11.根据权利要求9所述的装置，其特征在于，所述信号处理单元具体用于：对所述回波信号进行均值去噪信号处理或中值滤波信号处理或匹配滤波信号处理。
12.根据权利要求8所述的装置，其特征在于，所述飞行高度测量单元具体用于：
获取所述回波处理信号中的多个初始峰值采样点和多个波谷采样点，所述初始峰值采样点的纵坐标大于其前后采样点的纵坐标，所述波谷采样点是相邻所述初始峰值采样点之间的纵坐标最小的采样点，其中将第一个所述初始峰值采样点之前的第一个所述波谷采样点的纵坐标设置为零；
将纵坐标大于其之后某一个采样点纵坐标的第一预设倍数且大于其之前相邻的波谷采样点纵坐标的第二预设倍数的初始峰值采样点作为峰值采样点；
获取所有峰值采样点中纵坐标最大的采样点作为第一峰值采样点，纵坐标次最大的采样点作为第二峰值采样点。
13.根据权利要求12所述的装置，其特征在于，所述信号处理单元具体用于：对所述回波信号进行方差滤波信号处理。
14.根据权利要求8所述的装置，其特征在于，所述飞行高度测量单元还用于：若判断获知所述第二峰值采样点的纵坐标和所述第一峰值采样点的纵坐标的比值小于预设比值，则令所述第一高度距离与所述第二高度距离相等。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610916471.6A CN106595572B (zh) | 2016-10-20 | 2016-10-20 | 一种飞行器低空飞行高度测量方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610916471.6A CN106595572B (zh) | 2016-10-20 | 2016-10-20 | 一种飞行器低空飞行高度测量方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106595572A true CN106595572A (zh) | 2017-04-26 |
| CN106595572B CN106595572B (zh) | 2020-07-03 |
Family
ID=58556227
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610916471.6A Expired - Fee Related CN106595572B (zh) | 2016-10-20 | 2016-10-20 | 一种飞行器低空飞行高度测量方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106595572B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108496130A (zh) * | 2017-05-31 | 2018-09-04 | 深圳市大疆创新科技有限公司 | 飞行控制方法、设备、控制终端及其控制方法、无人机 |
| CN109375199A (zh) * | 2017-08-09 | 2019-02-22 | 宏碁股份有限公司 | 距离检测装置及其距离检测方法 |
| CN109633525A (zh) * | 2018-11-29 | 2019-04-16 | 上海无线电设备研究所 | 一种快速的均匀圆阵空间谱测向实现方法 |
| CN109673280A (zh) * | 2019-01-10 | 2019-04-26 | 珠海格力电器股份有限公司 | 基于微波雷达的栽培方法、装置、设备及处理器 |
| CN111480092A (zh) * | 2017-12-07 | 2020-07-31 | 威力登激光雷达有限公司 | 用于高效多返回光检测器的系统和方法 |
| CN113116153A (zh) * | 2019-12-31 | 2021-07-16 | 佛山市顺德区美的饮水机制造有限公司 | 水杯高度的测量方法、出水控制方法和出水控制装置 |
Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN87104480A (zh) * | 1986-04-18 | 1988-02-24 | 森德斯特兰德数据控制公司 | 无源无线电高度计 |
| CN1735896A (zh) * | 2002-11-05 | 2006-02-15 | 半导体技术及器械公司 | 凸块高度测量的系统和方法 |
| CN1932549A (zh) * | 2006-09-29 | 2007-03-21 | 重庆邮电大学 | 一种低空飞行目标无源定位方法及系统 |
| CN201266242Y (zh) * | 2007-12-27 | 2009-07-01 | 中国人民解放军空军工程大学 | 低空无人驾驶飞行器的高度测量半实物仿真装置 |
| CN201803718U (zh) * | 2010-09-19 | 2011-04-20 | 航宇救生装备有限公司 | 一种中低空高度测量装置 |
| CN102387109A (zh) * | 2010-09-06 | 2012-03-21 | 电信科学技术研究院 | 一种削除信号峰值功率的方法及装置 |
| CN102597801A (zh) * | 2009-07-02 | 2012-07-18 | 鹦鹉股份有限公司 | 能区别从另一无人驾驶飞机发出的伪回声的无人驾驶飞机的超声遥测方法 |
| CN102650689A (zh) * | 2012-05-17 | 2012-08-29 | 中国路桥工程有限责任公司 | 一种步进频率脉冲雷达位移测量方法 |
| JP2013170860A (ja) * | 2012-02-20 | 2013-09-02 | Pasuko:Kk | 計測点抽出プログラム、計測点抽出方法及び計測点抽出装置 |
| CN203385413U (zh) * | 2013-06-18 | 2014-01-08 | 桂林理工大学 | 一种适用于低空无人机的高度测量系统 |
| CN203803115U (zh) * | 2013-02-08 | 2014-09-03 | 上海科斗电子科技有限公司 | 浮力飞行器高度控制装置 |
| CN104502812A (zh) * | 2014-11-26 | 2015-04-08 | 国家电网公司 | 局部放电采集方法和装置 |
| CN104838229A (zh) * | 2012-10-23 | 2015-08-12 | 通用电气公司 | 超声波测量设备和方法 |
| CN105068063A (zh) * | 2015-07-16 | 2015-11-18 | 同济大学 | 一种基于回波全波形分解的地物特征获取方法 |
| CN105319546A (zh) * | 2015-04-02 | 2016-02-10 | 沈阳天择智能交通工程有限公司 | 多目标雷达及其所采用的数字信号处理方法 |
| CN205049136U (zh) * | 2015-10-21 | 2016-02-24 | 零度智控（北京）智能科技有限公司 | 无人机飞行高度测量装置 |
| CN205121343U (zh) * | 2015-11-08 | 2016-03-30 | 杨珊珊 | 基于无人飞行器的植被数据采集装置及无人飞行器 |
| CN105891828A (zh) * | 2016-04-22 | 2016-08-24 | 西安电子科技大学 | 一种机载cssar雷达动目标的检测方法 |
-
2016
- 2016-10-20 CN CN201610916471.6A patent/CN106595572B/zh not_active Expired - Fee Related
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN87104480A (zh) * | 1986-04-18 | 1988-02-24 | 森德斯特兰德数据控制公司 | 无源无线电高度计 |
| CN1735896A (zh) * | 2002-11-05 | 2006-02-15 | 半导体技术及器械公司 | 凸块高度测量的系统和方法 |
| CN1932549A (zh) * | 2006-09-29 | 2007-03-21 | 重庆邮电大学 | 一种低空飞行目标无源定位方法及系统 |
| CN201266242Y (zh) * | 2007-12-27 | 2009-07-01 | 中国人民解放军空军工程大学 | 低空无人驾驶飞行器的高度测量半实物仿真装置 |
| CN102597801A (zh) * | 2009-07-02 | 2012-07-18 | 鹦鹉股份有限公司 | 能区别从另一无人驾驶飞机发出的伪回声的无人驾驶飞机的超声遥测方法 |
| CN102387109A (zh) * | 2010-09-06 | 2012-03-21 | 电信科学技术研究院 | 一种削除信号峰值功率的方法及装置 |
| CN201803718U (zh) * | 2010-09-19 | 2011-04-20 | 航宇救生装备有限公司 | 一种中低空高度测量装置 |
| JP2013170860A (ja) * | 2012-02-20 | 2013-09-02 | Pasuko:Kk | 計測点抽出プログラム、計測点抽出方法及び計測点抽出装置 |
| CN102650689A (zh) * | 2012-05-17 | 2012-08-29 | 中国路桥工程有限责任公司 | 一种步进频率脉冲雷达位移测量方法 |
| CN104838229A (zh) * | 2012-10-23 | 2015-08-12 | 通用电气公司 | 超声波测量设备和方法 |
| CN203803115U (zh) * | 2013-02-08 | 2014-09-03 | 上海科斗电子科技有限公司 | 浮力飞行器高度控制装置 |
| CN203385413U (zh) * | 2013-06-18 | 2014-01-08 | 桂林理工大学 | 一种适用于低空无人机的高度测量系统 |
| CN104502812A (zh) * | 2014-11-26 | 2015-04-08 | 国家电网公司 | 局部放电采集方法和装置 |
| CN105319546A (zh) * | 2015-04-02 | 2016-02-10 | 沈阳天择智能交通工程有限公司 | 多目标雷达及其所采用的数字信号处理方法 |
| CN105068063A (zh) * | 2015-07-16 | 2015-11-18 | 同济大学 | 一种基于回波全波形分解的地物特征获取方法 |
| CN205049136U (zh) * | 2015-10-21 | 2016-02-24 | 零度智控（北京）智能科技有限公司 | 无人机飞行高度测量装置 |
| CN205121343U (zh) * | 2015-11-08 | 2016-03-30 | 杨珊珊 | 基于无人飞行器的植被数据采集装置及无人飞行器 |
| CN105891828A (zh) * | 2016-04-22 | 2016-08-24 | 西安电子科技大学 | 一种机载cssar雷达动目标的检测方法 |
Non-Patent Citations (4)
| Title |
|---|
| 严利华: "《机载雷达原理与系统》", 31 May 2010, 航空工业出版社 * |
| 陈行禄，秦永元: "《信号分析与处理》", 31 July 1993 * |
| 马士忠: "《测距机》", 31 August 1994, 国防工业出版社 * |
| 龙建忠，王勇: "《电路系统分析与设计》", 31 October 2002, 四川大学出版社 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108496130A (zh) * | 2017-05-31 | 2018-09-04 | 深圳市大疆创新科技有限公司 | 飞行控制方法、设备、控制终端及其控制方法、无人机 |
| CN108496130B (zh) * | 2017-05-31 | 2021-06-11 | 深圳市大疆创新科技有限公司 | 飞行控制方法、设备、控制终端及其控制方法、无人机 |
| CN109375199A (zh) * | 2017-08-09 | 2019-02-22 | 宏碁股份有限公司 | 距离检测装置及其距离检测方法 |
| CN109375199B (zh) * | 2017-08-09 | 2022-12-06 | 宏碁股份有限公司 | 距离检测装置及其距离检测方法 |
| CN111480092A (zh) * | 2017-12-07 | 2020-07-31 | 威力登激光雷达有限公司 | 用于高效多返回光检测器的系统和方法 |
| CN109633525A (zh) * | 2018-11-29 | 2019-04-16 | 上海无线电设备研究所 | 一种快速的均匀圆阵空间谱测向实现方法 |
| CN109673280A (zh) * | 2019-01-10 | 2019-04-26 | 珠海格力电器股份有限公司 | 基于微波雷达的栽培方法、装置、设备及处理器 |
| CN113116153A (zh) * | 2019-12-31 | 2021-07-16 | 佛山市顺德区美的饮水机制造有限公司 | 水杯高度的测量方法、出水控制方法和出水控制装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106595572B (zh) | 2020-07-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106595572A (zh) | 一种飞行器低空飞行高度测量方法及装置 | |
| CN103339525B (zh) | 用于监测地势变化的方法和设备 | |
| CN103048658A (zh) | 基于径向加速度的RA-Singer-EKF机动目标跟踪算法 | |
| CN104035095B (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
| CN108307767B (zh) | 适用于全自动除草机的障碍物检测避障系统及方法 | |
| CN106872969B (zh) | 基于mtd脉冲积累及滑动处理的雷达目标角度估计方法 | |
| CN106500671A (zh) | 一种基于lm算法分解激光雷达波形确定海水深度的方法 | |
| CN104215938B (zh) | 一种对地形反弹干扰的极化检测识别方法和装置 | |
| CN104714225B (zh) | 一种基于广义似然比的动态规划检测前跟踪方法 | |
| CN102288949B (zh) | 基于最优处理器的机动目标检测方法 | |
| CN106125088B (zh) | 基于激光雷达测深系统确定海水深度的方法 | |
| CN108959705B (zh) | 一种预测亚热带森林生物量的方法 | |
| CN103901430B (zh) | 双基前视sar俯冲段雷达发射机轨道设计方法 | |
| CN106443792B (zh) | 一种水下多缆声学网络定位方法 | |
| CN110687521B (zh) | 一种车载激光雷达标定的方法 | |
| CN110133648B (zh) | 一种选取逆合成孔径雷达船只成像时窗的方法 | |
| CN104483668A (zh) | 一种高精度雷达信号检测和跟踪系统及其方法 | |
| CN105738887A (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
| CN102621536B (zh) | 基于relax的空中多机动目标检测与参数估计方法 | |
| CN105445699A (zh) | 一种非视距误差消除的测距方法及系统 | |
| CN109932720A (zh) | 基于雷达的运动目标识别方法及装置 | |
| CN107092015B (zh) | 一种激光雷达回波信号散斑噪声的滤除方法 | |
| CN110488273A (zh) | 一种基于雷达的车辆跟踪检测方法及装置 | |
| CN109143262A (zh) | 无人驾驶汽车自动控制装置及其控制方法 | |
| CN104075710B (zh) | 一种基于航迹预测的机动扩展目标轴向姿态实时估计方法 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200703 |