CN109959907B - 一种对高度数据进行处理的方法 - Google Patents
一种对高度数据进行处理的方法 Download PDFInfo
- Publication number
- CN109959907B CN109959907B CN201711418045.0A CN201711418045A CN109959907B CN 109959907 B CN109959907 B CN 109959907B CN 201711418045 A CN201711418045 A CN 201711418045A CN 109959907 B CN109959907 B CN 109959907B
- Authority
- CN
- China
- Prior art keywords
- covariance
- value
- measurement
- state
- next moment
- 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.)
- Active
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
- 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
一种对高度数据进行处理的方法,包括:a、已知初始时刻的状态值和此状态对应的协方差值;b、根据状态转移矩阵,得到下一时刻k的状态预测值;c、根据测量矩阵、下一时刻的状态预测值和下一时刻k的测量值,计算得出当前时刻的测量余量;d、根据状态转移模型的噪声协方差,计算下一时刻的预测值对应的协方差;e、根据测量模型中的噪声方差以及上一步得到的预测值对应的协方差,计算得出测量余量协方差,进而得出增益矩阵;f、根据上一步的增益矩阵,得出下一时刻k的估计值和对应的协方差;g、根据上一步得出的数据,重复步骤a‑f进行迭代计算。本申请使无人机飞控***在高度解算时进行数据融合算法时提供一个稳定性强、精度高一个数据源。
Description
技术领域
本发明属于无人机雷达高度表领域,具体说是一种对高度数据进行处理的方法。
背景技术
目前,我国的油气管道总里程已近15万公里,纵横万里的油气管道为国民经济和居民的日常生活提供着重要的能源。所以能源传输的安全性成了重中之重,为了保证能源传输的安全,必须对它们进行定期巡查。尤其是对于长距离输油输气管道,但对于这些测量任务常常是周期短、任务重、质量高,区域内地形、地貌复杂多变,有时人员和测量仪器无法到达。
传统的人工巡线方法不仅工作量大而且条件艰苦,特别是对山区、河流、沼泽以及一些无人区等地区的输油输气管道的巡检;或者是在特殊时期,例如,水灾、火灾、地震、滑坡泥石流、夜晚期间巡线检查,所花时间长、人力成本大、危险性大。
无人机具有成本低廉、方便运输等特点,这些特点使得无人机很适合用于对输油输气管线巡检。无人机在对输油输气管线巡检时无论是摄像测绘还是照相测绘都需要对飞行高度有一个精确的测量。
发明内容
针对现有技术存在的上述缺点和不足,本发明提供了一种对高度数据进行处理的方法,可以提供实时、稳定的高度值输出。
一方面,本发明提供一种对高度数据进行处理的方法,包括滤波处理的步骤,包括:
a、已知初始时刻的状态值X(k-1),和此状态对应的协方差值P(k-1);
b、根据状态转移矩阵A,得到下一时刻k的状态预测值X(k|k-1);
c、根据测量矩阵H、下一时刻k的状态预测值X(k|k-1)和下一时刻k的测量值Z(k),计算得出当前时刻的测量余量Y(k);
d、根据状态转移模型的噪声协方差Q(k),计算下一时刻k的预测值对应的协方差P(k|k-1);
e、根据测量模型中的噪声方差R(k),以及上一步得到的预测值对应的协方差P(k|k-1),计算得出测量余量协方差S(k),进而得出增益矩阵K(k);
f、根据上一步的增益矩阵K(k),得出下一时刻k的估计值X(k)和对应的协方差P(k);
g、根据上一步得出的数据,重复步骤a-f进行迭代计算。
具体的,状态转移矩阵A为:T为帧与帧之间的时间间隔;下一时刻的状态预测值为:
X(k|k-1)=AX(k-1|k-1)。
更具体的,测量矩阵H为H=[1 0],当前k时刻的测量余量为,
Y(k)=Z(k)-H(k)X(k|k-1)。
更具体的,噪声协方差Q(k)为:
amax为最大的加速度值;下一时刻k的预测值对应的协方差P(k|k-1)为:
P(k|k-1)=AP(k-1|k-1)AT+Q(k)。
更具体的,余量协方差S(k):S(k)=H(k)P(k|k-1)HT(k)+R(k),其中R为测量对应的噪声协方差;
进而得出增益矩阵K(k)为:
K(k)=P(k|k-1)HT(k)S-1(k)。
作为更具体的,下一时刻k的估计值X(k)为:
X(k|k)=X(k|k-1)+K(k)Y(k)
对应的协方差P(k)为:
P(k|k)=(I-K(k)H(k))P(k|k-1)。
作为更具体的,本申请还包括数据平滑处理的步骤,如下:
本发明由于采用以上技术方案,能够取得如下的技术效果:
1、在上下扫频的过程中对于运动中产生的多普勒频移可以很好的进行补偿,解算高度准确,且不受设备运动的影响。
2、本***明显提高数据稳定性及数据连续性,对于解算时产生的数据跳变可以很好的抑制、剔除。
3、本发明设备能通过多种方式交互解算结果,可通过多种通信接口将实时解算高度发送至无人机飞控***,供飞控***进行使用,由于本设备结算高度的精确度高、稳定性好,所以可以用于与其他高度解算设备进行数据融合提高解算高度的准确性。不仅如此,本设备还能通过无线传输模块将实时解算高度下发至地面站,供后期使用。
附图说明
本发明共有附图4幅:
图1为发射波形示意图;
图2为实施例中一种对高度数据进行处理的方法中硬件部分结构框图;
图3为实施例中一种对高度数据进行处理的方法流程图;
图4为实施例中滤波处理的流程图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步的具体说明。本发明所使用的毫米波雷达的工作频率在24GHz,60GHz或者77GHz,采用FMCW调频连续波体制。使用FMCW体制雷达主要是因为辐射功率小、测距精度高、设备相对简单、具有良好的电子对抗(ECM)性能等。工作原理是用回波信号和发射信号的一部分进行相干混频,得到差频信号,然后对差频信号进行处理和检测得到相应的距离信息。通常情况下线性调频使用锯齿波或者三角波。由于使用三角波调制可以抵消动目标在测量时的多普勒效应,所以发明使用三角波调制。本实施例设计的雷达高度表高度测量范围为20~200m,距离分辨率为0.75m,雷达中心频率为24.125GHz,发射波形使用对称三角波,周期为20ms,带宽为200M,发射波形如图1所示。
实施例1
本实施例提供一种对高度数据进行处理的方法,包括:
AD数据输入:通过Arm芯片的ADC模块采集一路I、Q数据;
去直流:对I、Q数据中的每一个数据分别减去各自数据总和的均值;
FFT计算:对去直流后的I、Q数据进行FFT变换,将时域数据转换成频域数据,并对计算后的复数取模;
非相干积累:采集3个周期的FFT计算数据,并对N个周期内每个点的数据求和后取平均值得到当前点的新数值;由于累积3个周期信号,在周期时间很短的情况下,有效目标幅度和位置应该是不会发生变化的,但是噪声在每个周期内都是会不断发生变化。又由于无人机在飞行过程中,有可能在某一时刻由于飞行器角度变化问题导致有效目标的幅度较小,此时在进行CFAR门限检测就有可能检测不到有效目标或者解算错误。采用三个周期的非相干累积,有效目标幅度经过累积、平均后,幅度达到相对稳定,同时对干扰信号进行累积、平均后,将会大幅降低干扰信号对于有效信号检测的影响,降低虚警率,此时再进行CFAR门限检测可以大幅提高有效目标被检测出的概率,对于增强数据稳定性和连续性起到非常重要的作用。
门限检测:对新数值进行CFAR门限检测,找出过门限点幅度最大的一个点;
解算距离:根据门限检测模块得出上扫频过门限最大点坐标为up_max_index,该点对应的频率为fd_up,下扫频过门限最大点坐标为down_max_index,该点对应的频率为fd_down,进而计算得出高度H,
其中C为光速,C=3*108;T为周期,T=2ms;B为调制带宽,B=300MHz;
数据处理模块,对解算出的高度H进行数据处理,以达到过滤、稳定数据的目的。
实施例2
作为对实施例1的补充,上述方法,还包括:无人机雷达高度表在石油管线巡线中的信号处理***,具体为:
射频前端:通过PLL芯片调制射频芯片将发射信号从发射天线发送出去,并且将接收天线接收到的回波信号和发射信号的一部分进行相干混频,得到中频信号通过GPIO接口传输至Arm芯片的ADC模块;
Arm芯片:通过SPI总线配置PLL芯片相关参数;ADC模块接收到的中频信号通过DMA通道送入CPU进行高度解算,CPU将解算出的距离值通过接口与无人机飞控***和地面站进行数据交互;
电源模块:为本***提供电源;
外部电路,包括晶振,为pll芯片提供基准时钟;RS232接口:Arm芯片输出TTL电平转换成232电平,可用MAXIM公司的max3232芯片进行电平转换;IIC接口,使用Arm芯片内部的IIC总线直接输出;CAN总线接口:可用MICROCHIP公司的MCP2515芯片进行电平转换;WiFi模块:采用LSD4WF-2MD05107。
CPU使用RS232、CAN总线及IIC总线与飞控***进行数据交互,CPU通过2.4GHzWiFi模块与地面站进行数据交互。
实施例3
作为对实施例1的补充,数据处理模块中对高度H进行数据处理,包括滤波处理模块,已知初始时刻的状态值X(k-1),和此状态对应的协方差值P(k-1);根据状态转移矩阵A为得到下一时刻k的状态预测值X(k|k-1):
X(k|k-1)=AX(k-1|k-1)。
根据测量矩阵H=[1 0]、下一时刻k的状态预测值X(k|k-1)和下一时刻k的测量值Z(k),该测量值Z(k)为解算距离后的高度H,得出当前时刻的测量余量Y(k);
Y(k)=Z(k)-H(k)X(k|k-1);
根据状态转移模型的噪声协方差amax为最大的加速度值,计算下一时刻k的预测值对应的协方差P(k|k-1),P(k|k-1)=AP(k-1|k-1)AT+Q(k);
根据测量模型中的噪声方差R(k),以及上一步得到的预测值对应的协方差P(k|k-1),计算得出测量余量协方差S(k),S(k)=H(k)P(k|k-1)HT(k)+R(k),本***测量模型中的噪声方差一般采用方差为0.75米,R=0.75;进而得出增益矩阵K(k),K(k)=P(k|k-1)HT(k)S-1(k);
根据增益矩阵K(k),得出下一时刻k的估计值X(k)和对应的协方差P(k);
下一时刻k的估计值X(k)为:
X(k|k)=X(k|k-1)+K(k)Y(k)
对应的协方差P(k)为:
P(k|k)=(I-K(k)H(k))P(k|k-1)。
数据处理模块中对高度H进行数据处理,包括数据平滑处理模块,累积N包解算高度,进行取平均处理,如下:
本发明的多种通信接口,不仅可以与地面站进行无线传输,而且可以与无人机飞控***实时进行数据交互,使无人机飞控***在高度解算时进行数据融合算法时提供一个稳定性强、精度高一个数据源。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (7)
1.一种对高度数据进行处理的方法,其特征在于,包括滤波处理的步骤,
a、已知初始时刻的状态值X(k-1),和此状态对应的协方差值P(k-1);
b、根据状态转移矩阵A,得到下一时刻k的状态预测值X(k|k-1);
c、根据测量矩阵H、下一时刻k的状态预测值X(k|k-1)和下一时刻k的测量值Z(k),计算得出当前时刻的测量余量Y(k);
d、根据状态转移模型的噪声协方差Q(k),计算下一时刻k的预测值对应的协方差P(k|k-1);
e、根据测量模型中的噪声方差R(k),以及上一步得到的预测值对应的协方差P(k|k-1),计算得出测量余量协方差S(k),进而得出增益矩阵K(k);
f、根据上一步的增益矩阵K(k),得出下一时刻k的估计值X(k)和对应的协方差P(k);
g、根据上一步得出的数据,重复步骤a-f进行迭代计算。
2.根据权利要求1所述一种对高度数据进行处理的方法,其特征在于,状态转移矩阵A为:T为帧与帧之间的时间间隔;下一时刻的状态预测值为:
X(k|k-1)=AX(k-1|k-1)。
3.根据权利要求1所述一种对高度数据进行处理的方法,其特征在于,测量矩阵H为H=[1 0],当前k时刻的测量余量为,
Y(k)=Z(k)-H(k)X(k|k-1)。
4.根据权利要求1所述一种对高度数据进行处理的方法,其特征在于,噪声协方差Q(k)为:
amax为最大的加速度值;下一时刻k的预测值对应的协方差P(k|k-1)为:
P(k|k-1)=AP(k-1|k-1)AT+Q(k)。
5.根据权利要求4所述一种对高度数据进行处理的方法,其特征在于,余量协方差S(k):S(k)=H(k)P(k|k-1)HT(k)+R(k),其中R为测量对应的噪声协方差;
进而得出增益矩阵K(k)为:
K(k)=P(k|k-1)HT(k)S-1(k)。
6.根据权利要求4所述一种对高度数据进行处理的方法,其特征在于,下一时刻k的估计值X(k)为:
X(k|k)=X(k|k-1)+K(k)Y(k)
对应的协方差P(k)为:
P(k|k)=(I-K(k)H(k))P(k|k-1)。
7.根据权利要求2所述一种对高度数据进行处理的方法,其特征在于,本方法还包括数据平滑处理的步骤,如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711418045.0A CN109959907B (zh) | 2017-12-25 | 2017-12-25 | 一种对高度数据进行处理的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711418045.0A CN109959907B (zh) | 2017-12-25 | 2017-12-25 | 一种对高度数据进行处理的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109959907A CN109959907A (zh) | 2019-07-02 |
CN109959907B true CN109959907B (zh) | 2023-10-20 |
Family
ID=67020648
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711418045.0A Active CN109959907B (zh) | 2017-12-25 | 2017-12-25 | 一种对高度数据进行处理的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109959907B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5617099A (en) * | 1996-01-22 | 1997-04-01 | Hughes Aircraft Company | Adaptive filtering of matched-filter data |
JP2009014596A (ja) * | 2007-07-06 | 2009-01-22 | Mitsubishi Electric Corp | 目標追尾装置及び目標追尾方法 |
CN102721951A (zh) * | 2012-05-04 | 2012-10-10 | 西安电子科技大学 | 一种高机动目标跟踪方法 |
CN104180801A (zh) * | 2014-08-20 | 2014-12-03 | 广州海格通信集团股份有限公司 | 基于ads-b***航迹点的预测方法和*** |
-
2017
- 2017-12-25 CN CN201711418045.0A patent/CN109959907B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5617099A (en) * | 1996-01-22 | 1997-04-01 | Hughes Aircraft Company | Adaptive filtering of matched-filter data |
JP2009014596A (ja) * | 2007-07-06 | 2009-01-22 | Mitsubishi Electric Corp | 目標追尾装置及び目標追尾方法 |
CN102721951A (zh) * | 2012-05-04 | 2012-10-10 | 西安电子科技大学 | 一种高机动目标跟踪方法 |
CN104180801A (zh) * | 2014-08-20 | 2014-12-03 | 广州海格通信集团股份有限公司 | 基于ads-b***航迹点的预测方法和*** |
Non-Patent Citations (1)
Title |
---|
自适应卡尔曼定位算法在分布式MIMO***应用;牟亚南;冀小平;凌宗平;邢伟;;电声技术(第Z2期);140-144 * |
Also Published As
Publication number | Publication date |
---|---|
CN109959907A (zh) | 2019-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11506512B2 (en) | Method and system using tightly coupled radar positioning to improve map performance | |
CN107607923B (zh) | 基于lfmcw雷达的振动监测***与信号处理方法 | |
CN106405556B (zh) | 车辆目标信息探测识别***及其信号处理方法 | |
CN110609287B (zh) | 一种双频雷达散射计及同时测量海面风场和流场的方法 | |
KR100795497B1 (ko) | 레이더를 이용한 파랑 계측 시스템 및 방법 | |
CN104391281A (zh) | 提高天波雷达海面船舶目标跟踪定位精度的方法 | |
WO2021077287A1 (zh) | 一种检测方法、检测装置以及存储介质 | |
CN110730913B (zh) | 退化可视环境的分布式多节点低频雷达***的方法和设备 | |
CN104252010B (zh) | 一种无线电探空仪及其气象数据测量方法 | |
US9310481B2 (en) | Wide band clear air scatter doppler radar | |
CN103914983A (zh) | 一种交通卡口管理雷达装置及实现方法 | |
CN113108767A (zh) | 一种无人机载雷达水文信息实时监测方法 | |
CN110988884B (zh) | 一种基于高频地波雷达的中纬度电离层探测方法 | |
CN104007286A (zh) | 一种用于河道表面流速的自动监测*** | |
CN106483495B (zh) | 一种室内运动标签的定位和测速方法 | |
CN107271725A (zh) | 激光测风雷达走航测量速度校准的方法 | |
CN101762812A (zh) | 一种机载单站无源定位***对目标辐射波长的测定方法 | |
CN113721253B (zh) | 基于fmcw激光雷达的运动物体速度检测方法 | |
CN107843895A (zh) | 一种双多普勒雷达三维风场反演方法 | |
Wang et al. | Highly accurate noncontact water level monitoring using continuous-wave Doppler radar | |
CN111766599A (zh) | 一种实时测量海浪高度的激光雷达装置及方法 | |
JP6567665B2 (ja) | 船舶の航路(パス)の全地点における各々のドリフト(浮動)ベクトルを推定する方法 | |
WO2019114299A1 (zh) | 一种地表冻融状态延迟多普勒图仿真模型的建立方法 | |
CN117907983B (zh) | 激光雷达测距测速的方法、激光雷达以及信号处理设备 | |
CN109959907B (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 |