CN107942307B - 调频连续波着陆雷达波束的地面入射角估计方法 - Google Patents

调频连续波着陆雷达波束的地面入射角估计方法 Download PDF

Info

Publication number
CN107942307B
CN107942307B CN201710952931.5A CN201710952931A CN107942307B CN 107942307 B CN107942307 B CN 107942307B CN 201710952931 A CN201710952931 A CN 201710952931A CN 107942307 B CN107942307 B CN 107942307B
Authority
CN
China
Prior art keywords
frequency modulation
positive
echo
frequency
negative
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
Application number
CN201710952931.5A
Other languages
English (en)
Other versions
CN107942307A (zh
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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201710952931.5A priority Critical patent/CN107942307B/zh
Publication of CN107942307A publication Critical patent/CN107942307A/zh
Application granted granted Critical
Publication of CN107942307B publication Critical patent/CN107942307B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了一种调频连续波雷达波束地面入射角的估计方法，主要解决现有技术只能估计地面回波信号在天线上的入射角度，无法估计天线发射信号相对地面的入射角度的问题。其实现方案是：1、获得基于三角调制的调频连续波雷达分别在正、负调频周期的回波信号；2、分别对正调频和负调频回波信号进行解调频、FFT及非相参积累3、分别估计正、负调频阶段的噪声功率；4、分别估计正、负调频阶段回波包络的重心和回波包络的宽度5、计算回波重心对应的距离和回波宽度；6、利用回波宽度与回波重心对应距离的比值估计波束在地面的入射角。本发明能估计出雷达波束相对地面的入射角，可用于校正重心法估计出的波束重心，实现波束重心的无偏估计。

Description

调频连续波着陆雷达波束的地面入射角估计方法
技术领域
本发明属于数据处理技术领域，特别涉及一种调频连续波着陆雷达波束的地面入射角估计方法，可用于对雷达波束重心的校正。
背景技术
随着雷达技术的迅速发展和对调频连续波雷达的深入研究，调频连续波雷达已广泛应用于测距、测速、气象预测等领域。调频连续波雷达具有距离分辨率高，不存在距离测量盲区，系统容易实现等优点。线性频率调制的调频连续波信号作为一种比较成熟的低截获概率信号，目前已受到广泛的研究和应用。而在实际应用中，由于雷达天线波束在地面等角度间隔的投影面积不同以及天线波束在地面的入射角不同，调频连续波雷达接收到的回波信号的能量并不是关于天线波束中心轴线对称的，用重心法估计得到的频谱中心频率并非天线波束中心对应的频率。为实现无偏估计，需要对重心估计的结果进行校正，而重心校正因子由波束在地面的入射角唯一确定。因此，如何估计出天线波束在地面的入射角具有重要的现实意义。
现有的估计波束入射角的方法只能估计地面回波信号在天线上的入射角度，无法估计天线发射信号相对地面的入射角度。当调频连续波雷达运用于着陆器时，雷达接收到的面目标回波信号能量并不是关于天线波束中心轴线对称的，若只估计地面回波信号在天线上的入射角，则无法对着陆器雷达测速测距的结果进行校正，测速测距结果精度降低；同时为获得着陆器的姿态辅助信息，也需要对天线发射信号相对地面的入射角进行估计；着陆器在着陆的过程中需要观测地面的起伏情况，从而决定着陆姿态，在着陆器平台不动的情况下，只估计地面回波信号在天线上的入射角也无法获得地面的起伏情况。
发明内容
本发明的目的在于提出了一种调频连续波着陆雷达波束的地面入射角估计方法，以解决现有技术只能估计地面回波信号在天线上的入射角度，无法估计天线发射信号相对地面的入射角度的问题，从而实现对着陆器雷达的测速测距结果进行校正，获得着陆器姿态的辅助信息，及在着陆器平台不动的情况下获得地面的起伏情况。
为了实现上述目的，本发明的技术方案包括如下：
(1)获取基于三角调制的调频连续波雷达分别在正、负调频周期的回波信号
Figure GDA0002958386080000011
和
Figure GDA0002958386080000012
(2)设共有N个调频周期，每个调频周期有M个采样点，对正、负调频周期的回波信号分别进行解调频处理并排列，得到正、负调频周期的差频信号二维矩阵：s+[tr(m),ta(n)]、s-[tr(m),ta(n)]，其中tr(m)表示调频周期内第m个采样点的采样时间，ta(n)表示第n个调频周期的起始时间，m＝1、2、···、M,n＝1、2、···、N；
(3)分别对正、负调频周期的差频信号二维矩阵进行FFT处理和非相参积累，得到一维的正、负调频周期差频信号频谱s'+[fr(m)]、s'-[fr(m)]：
3a)对正、负调频周期的差频信号二维矩阵进行距离维和方位维的FFT处理，得到正、负调频周期差频信号频谱s+[fr(m),fa(n)]、s-[fr(m),fa(n)]，其中fr(m)表示距离维第m个采样点对应的频率，fa(n)表示方位维第n个调频周期对应的频率；
3b)分别对经过FFT处理后得到的正、负调频周期差频信号频谱进行非相参积累，积累的个数为调频周期数，得到非相参积累后的一维的正、负调频周期差频信号频谱s'+[fr(m)]、s'-[fr(m)]为：
Figure GDA0002958386080000021
Figure GDA0002958386080000022
其中，N为调频周期数；
(4)分别估计一维的正、负调频周期差频信号频谱s'+[fr(m)]、s'-[fr(m)]的噪声功率P1和P2：
4a)在一维的正调频周期差频信号频谱s'+[fr(m)]的峰值两端各取Np个采样点，进行一次重心估计，得出正频率回波包络的重心
Figure GDA0002958386080000023
和回波包络的宽度
Figure GDA0002958386080000024
4b)取
Figure GDA0002958386080000025
所在位置左侧的Np1个采样点估计正调频周期噪声功率P1；
4c)在一维的负调频周期差频信号频谱s'-[fr(m)]的峰值两端各取Np2个采样点，进行一次重心估计，得出负频率回波包络的重心
Figure GDA0002958386080000026
和回波包络的宽度
Figure GDA0002958386080000027
4d)取
Figure GDA0002958386080000028
所在位置右侧的Np3个采样点估计负调频周期噪声功率P2；
(5)利用正、负调频周期噪声功率P1和P2，分别计算正、负调频周期的回波包络重心
Figure GDA0002958386080000031
Figure GDA0002958386080000032
和回波包络宽度
Figure GDA0002958386080000033
5a)利用重心法，分别以
Figure GDA0002958386080000034
为中心、
Figure GDA0002958386080000035
为新窗口宽度，在窗口内重新计算正调频周期的回波包络重心
Figure GDA0002958386080000036
和回波包络宽度
Figure GDA0002958386080000037
5b)利用重心法，分别以
Figure GDA0002958386080000038
为中心、
Figure GDA0002958386080000039
为新窗口宽度，在窗口内重新计算负调频周期的回波包络重心
Figure GDA00029583860800000310
和回波包络宽度
Figure GDA00029583860800000311
(6)利用正、负调频周期的回波包络重心
Figure GDA00029583860800000312
和回波包络宽度
Figure GDA00029583860800000313
计算回波重心对应的距离
Figure GDA00029583860800000314
和回波包络宽度
Figure GDA00029583860800000315
(7)计算回波包络宽度
Figure GDA00029583860800000316
与回波重心对应的距离
Figure GDA00029583860800000317
的比值
Figure GDA00029583860800000318
Figure GDA00029583860800000319
利用比值
Figure GDA00029583860800000320
估计天线波束在地面的入射角
Figure GDA00029583860800000321
Figure GDA00029583860800000322
其中，α(θ)为在不同入射角条件下已知理论比值。
本发明将估计得到的频谱宽度和回波重心对应距离的比值与理论值进行比较，用估计值最接近的理论值对应的入射角作为雷达波束在地面的入射角估计值，解决了现有方法只能估计地面回波信号相对天线入射角的问题，从而可实现对着陆器雷达的测速测距结果校正，获得着陆器姿态的辅助信息，及在着陆器平台不动的情况下获得地面的起伏情况。
附图说明
下面结合附图和具体实施方式对本发明做进一步详细说明。
图1是本发明的实现流程图；
图2是不同入射角条件下回波包络宽度与回波重心对应距离的理论比值；
图3是用本发明计算雷达波束相对地面的入射角估计值与理论值的对比图。
具体实施方式
参照图1，本发明的实现步骤如下：
步骤1:获得基于三角调制的调频连续波雷达分别在正、负调频周期的回波信号
Figure GDA00029583860800000323
和
Figure GDA00029583860800000324
雷达发射一个三角调制的调频连续波信号，经过时延τ后接收到目标反射的回波信号，该回波信号包含正调频周期信号
Figure GDA0002958386080000041
和负调频周期信号
Figure GDA0002958386080000042
其分别表示如下：
Figure GDA0002958386080000043
Figure GDA0002958386080000044
其中，Et为雷达发射信号的幅度，f0为中心频率，μ为调制斜率，
Figure GDA0002958386080000045
为随机初相位，τ为回波延时，R(t)为目标到雷达的距离，G为天线功率增益，λ为波长，σ为雷达截面积，t为时间，j为虚数单位。
步骤2:分别对正、负调频周期回波信号进行解调频处理。
解调频处理就是将雷达发射信号与回波信号进行共轭相乘，得到差频信号的过程，其实现如下：
2a)对正调频周期的回波信号进行解调频处理：设正调频周期发射信号为
Figure GDA0002958386080000046
接收到的回波信号为
Figure GDA0002958386080000047
则解调频后得到的差频信号
Figure GDA0002958386080000048
表示为：
Figure GDA0002958386080000049
其中
Figure GDA00029583860800000410
表示对正调频周期发射信号
Figure GDA00029583860800000411
取共轭；
2b)对负调频周期的回波信号进行解调频处理：设负调频周期发射信号为
Figure GDA00029583860800000412
接收到的回波信号为
Figure GDA00029583860800000413
则解调频后得到的差频信号
Figure GDA00029583860800000414
表示为：
Figure GDA00029583860800000415
其中
Figure GDA00029583860800000416
表示对负调频周期发射信号
Figure GDA00029583860800000417
取共轭；
2c)设共有N个调频周期，每个调频周期有M个采样点，对正、负调频周期的回波信号分别进行解调频处理并排列，得到正、负调频周期的差频信号二维矩阵：s+[tr(m),ta(n)]、s-[tr(m),ta(n)]，其中tr(m)表示调频周期内第m个采样点的采样时间，ta(n)表示第n个调频周期的起始时间，m＝1、2、···、M,n＝1、2、···、N。
步骤3:分别对正、负调频周期差频信号进行FFT处理和非相参积累。
解调频后得到的正、负周期差频信号中包含目标的速度和距离等信息，为获取差频信号中的有用信息，需先得到差频信号的频谱，即需要分别对正、负调频周期差频信号进行距离维和方位维的FFT处理，再对正、负调频周期差频信号频谱进行非相参积累，其步骤如下：
3a)对正、负调频周期差频信号进行距离维的FFT处理，得到正调频周期的距离维差频频谱s+[fr(m),ta(n)]和负调频周期的距离维差频频谱s-[fr(m),ta(n)]：
Figure GDA0002958386080000051
Figure GDA0002958386080000052
其中，fr(m)表示距离维第m个采样点对应的频率；
3b)对得到的正、负调频周期的距离维差频频谱进行方位维的FFT处理，得到正调频周期差频信号频谱s+[fr(m),fa(n)]和负调频周期差频信号频谱s-[fr(m),fa(n)]：
Figure GDA0002958386080000053
Figure GDA0002958386080000054
其中，fa(n)表示方位维第n个调频周期对应的频率；
3c)对经过FFT处理后得到的正、负调频周期差频信号频谱分别进行非相参积累，积累的个数为调频周期数，非相参积累后得到一维的正调频周期差频信号频谱s'+[fr(m)]和一维的负调频周期差频信号频谱s'-[fr(m)]表示如下：
Figure GDA0002958386080000055
Figure GDA0002958386080000056
步骤4:分别估计正、负调频阶段的噪声功率P1和P2。
为了充分利用整个波束主瓣内的回波信号，采用重心法来对频谱包络重心进行估计。为提高估计精度，需要先分别估计出正、负调频阶段差频信号频谱s'+[fr(m)]、s'-[fr(m)]的噪声功率P1和P2，具体步骤如下：
4a)为了估计正调频阶段噪声功率，首先在峰值两端各取Np个采样点，进行一次重心估计：假设正调频周期差频频谱s'+[fr(m)]的峰值位置为
Figure GDA0002958386080000057
计算正频率回波包络的重心
Figure GDA0002958386080000058
和回波包络的宽度
Figure GDA0002958386080000059
计算公式如下：
Figure GDA0002958386080000061
Figure GDA0002958386080000062
4b)取
Figure GDA0002958386080000063
所在位置左侧的Np1个采样点估计正调频阶段噪声功率P1：
Figure GDA0002958386080000064
4c)为了估计负调频阶段噪声功率，首先在峰值两端各取Np2个采样点，进行一次重心估计：假设负调频周期差频频谱s'-[fr(m)]的峰值位置为
Figure GDA0002958386080000065
计算负频率回波包络的重心
Figure GDA0002958386080000066
和回波包络的宽度
Figure GDA0002958386080000067
计算公式如下：
Figure GDA0002958386080000068
Figure GDA0002958386080000069
4d)取
Figure GDA00029583860800000610
所在位置右侧的Np3个采样点估计负调频阶段噪声功率P2：
Figure GDA00029583860800000611
步骤4中的各采样点数数值相等，即Np＝Np1＝Np2＝Np3。
步骤5:分别估计正、负调频阶段回波包络的重心
Figure GDA00029583860800000612
和正、负回波包络的宽度
Figure GDA00029583860800000613
Figure GDA0002958386080000071
由于在不同入射角条件下，回波包络宽度和回波重心对应距离的比值随入射角单调变化，因此，利用这一特性来估计雷达波束的地面入射角，需要先计算出回波包络的重心和回波包络宽度，其具体实现步骤如下：
5a)估计正频率的回波包络重心
Figure GDA0002958386080000072
和回波包络宽度
Figure GDA0002958386080000073
根据估计噪声功率时计算出的正频率回波包络的重心
Figure GDA0002958386080000074
和回波包络的宽度
Figure GDA0002958386080000075
以
Figure GDA0002958386080000076
为中心、
Figure GDA0002958386080000077
为新窗口宽度，在窗口内重新计算正频率的回波包络重心
Figure GDA0002958386080000078
和回波包络宽度
Figure GDA0002958386080000079
计算公式如下：
Figure GDA00029583860800000710
Figure GDA00029583860800000711
5b)估计负频率的回波包络重心
Figure GDA00029583860800000712
和回波包络宽度
Figure GDA00029583860800000713
根据估计噪声功率时计算出的负频率回波包络的重心
Figure GDA00029583860800000714
和回波包络的宽度
Figure GDA00029583860800000715
以
Figure GDA00029583860800000716
为中心、
Figure GDA00029583860800000717
为新窗口宽度，在窗口内重新计算负频率的回波包络重心
Figure GDA00029583860800000718
和回波包络宽度
Figure GDA00029583860800000719
计算公式如下：
Figure GDA00029583860800000720
Figure GDA0002958386080000081
5c)根据5a)中计算得出的正频率的回波包络重心
Figure GDA0002958386080000082
和回波包络宽度
Figure GDA0002958386080000083
分别以
Figure GDA0002958386080000084
为中心、
Figure GDA0002958386080000085
为新窗口宽度重复步骤5a)，得到最终正调频周期的回波重心
Figure GDA0002958386080000086
和回波包络的宽度
Figure GDA0002958386080000087
5d)根据5b)中计算得出的负频率的回波包络重心
Figure GDA0002958386080000088
和回波包络宽度
Figure GDA0002958386080000089
分别以
Figure GDA00029583860800000810
为中心、
Figure GDA00029583860800000811
为新窗口宽度重复步骤5b)，得到最终负调频周期的回波重心
Figure GDA00029583860800000812
和回波包络的宽度
Figure GDA00029583860800000813
步骤6:计算回波重心对应的距离
Figure GDA00029583860800000814
和回波宽度
Figure GDA00029583860800000815
根据步骤5得到正、负调频周期的回波重心
Figure GDA00029583860800000816
和回波包络宽度
Figure GDA00029583860800000817
计算回波重心对应的距离和回波宽度：
Figure GDA00029583860800000818
Figure GDA00029583860800000819
式中c为光速，μ为调频率。
步骤7:利用回波宽度
Figure GDA00029583860800000820
与回波包络重心对应距离
Figure GDA00029583860800000821
的比值估计波束在地面的入射角
Figure GDA00029583860800000822
7a)计算在不同的入射角条件下回波宽度
Figure GDA00029583860800000823
与回波重心对应距离
Figure GDA00029583860800000824
的比值
Figure GDA00029583860800000825
Figure GDA00029583860800000826
7b)利用比值
Figure GDA00029583860800000827
估计出波束的地面入射角：
将不同入射角条件下回波宽度与回波包络重心对应距离的理论比值记为α(θ)，得出如图2所示的对应关系，从图2可见，理论比值α(θ)与波束的地面入射角成单调关系，故将估计比值
Figure GDA00029583860800000828
与理论比值α(θ)进行比较，把与估计比值最接近的理论值对应的入射角作为雷达波束在地面的入射角的估计值
Figure GDA00029583860800000829
Figure GDA0002958386080000091
此时，估计出了雷达天线波束在地面的入射角，进而可实现对重心法估计得到的测速测距结果进行校正，获得着陆器姿态的辅助信息，及在着陆器平台不动的情况下获得地面的起伏情况。
本发明的效果可通过以下仿真进一步验证。
1.实验场景：
设线性调频连续波发射信号中心频率f0＝35GHz，带宽B＝4MHz，脉冲重复周期1ms，采样频率fs＝20MHz，随机初相位
Figure GDA0002958386080000092
目标到雷达的距离范围为R∈[5000m,6000m]，天线波束在地面的入射角范围为θ∈[20°,40°]，仿真时间为[0,153.6s]，数据更新周期为0.128s,天线方向图使用sinc函数。
对发射信号的幅度做归一化处理后，生成线性调频连续波雷达回波信号。将正、负调频周期的回波信号进行解调频处理并排列为二维矩阵后，对正调频和负调频差频信号分别进行距离维和方位维的FFT及非相参积累处理。
在估计正、负调频阶段噪声功率时，采样点数Np＝Np1＝Np2＝Np3＝50。
2.实验内容与结果：
根据CE-3天线实测的雷达回波数据，计算出在不同入射角条件下回波宽度与回波包络重心对应距离的理论比值，得到单调曲线，如图2所示。
根据本发明的实现步骤，对三角调制的线性调频连续波雷达进行仿真，得出雷达波束相对地面的入射角估计值，结果如图3所示。
从图3可以看出，雷达波束相对地面的入射角估计值随着理论值上下波动，入射角估计值相对理论值的误差均值为-1.6097°，均方根误差为2.3959°，可见采用本发明可以较为准确的估计出雷达波束相对地面的入射角，从而用于校正重心法估计出的波束重心，实现波束重心的无偏估计。

Claims (10)

1.调频连续波着陆雷达波束的地面入射角估计方法，包括：
(1)获取基于三角调制的调频连续波雷达分别在正、负调频周期的回波信号
Figure FDA0002958386070000011
和
Figure FDA0002958386070000012
(2)设共有N个调频周期，每个调频周期有M个采样点，对正、负调频周期的回波信号分别进行解调频处理并排列，得到正、负调频周期的差频信号二维矩阵：s+[tr(m),ta(n)]、s-[tr(m),ta(n)]，其中tr(m)表示调频周期内第m个采样点的采样时间，ta(n)表示第n个调频周期的起始时间，m＝1、2、···、M,n＝1、2、···、N；
(3)分别对正、负调频周期的差频信号二维矩阵进行FFT处理和非相参积累，得到一维的正、负调频周期差频信号频谱s'+[fr(m)]、s'-[fr(m)]：
3a)对正、负调频周期的差频信号二维矩阵进行距离维和方位维的FFT处理，得到正、负调频周期差频信号频谱s+[fr(m),fa(n)]、s-[fr(m),fa(n)]，其中fr(m)表示距离维第m个采样点对应的频率，fa(n)表示方位维第n个调频周期对应的频率；
3b)分别对经过FFT处理后得到的正、负调频周期差频信号频谱进行非相参积累，积累的个数为调频周期数，得到非相参积累后的一维的正、负调频周期差频信号频谱s'+[fr(m)]、s'-[fr(m)]为：
Figure FDA0002958386070000013
Figure FDA0002958386070000014
其中，N为调频周期数；
(4)分别估计一维的正、负调频周期差频信号频谱s'+[fr(m)]、s'-[fr(m)]的噪声功率P1和P2：
4a)在一维的正调频周期差频信号频谱s'+[fr(m)]的峰值两端各取Np个采样点，进行一次重心估计，得出正频率回波包络的重心
Figure FDA0002958386070000015
和回波包络的宽度
Figure FDA0002958386070000016
4b)取
Figure FDA0002958386070000017
所在位置左侧的Np1个采样点估计正调频周期噪声功率P1；
4c)在一维的负调频周期差频信号频谱s'-[fr(m)]的峰值两端各取Np2个采样点，进行一次重心估计，得出负频率回波包络的重心
Figure FDA0002958386070000021
和回波包络的宽度
Figure FDA0002958386070000022
4d)取
Figure FDA0002958386070000023
所在位置右侧的Np3个采样点估计负调频周期噪声功率P2；
(5)利用正、负调频周期噪声功率P1和P2，分别计算正、负调频周期的回波包络重心
Figure FDA0002958386070000024
Figure FDA0002958386070000025
和回波包络宽度
Figure FDA0002958386070000026
5a)利用重心法，分别以
Figure FDA0002958386070000027
为中心、
Figure FDA0002958386070000028
为新窗口宽度，在窗口内重新计算正调频周期的回波包络重心
Figure FDA0002958386070000029
和回波包络宽度
Figure FDA00029583860700000210
5b)利用重心法，分别以
Figure FDA00029583860700000211
为中心、
Figure FDA00029583860700000212
为新窗口宽度，在窗口内重新计算负调频周期的回波包络重心
Figure FDA00029583860700000213
和回波包络宽度
Figure FDA00029583860700000214
(6)利用正、负调频周期的回波包络重心
Figure FDA00029583860700000215
和回波包络宽度
Figure FDA00029583860700000216
计算回波重心对应的距离
Figure FDA00029583860700000217
和回波包络宽度
Figure FDA00029583860700000218
(7)计算回波包络宽度
Figure FDA00029583860700000219
与回波重心对应的距离
Figure FDA00029583860700000220
的比值
Figure FDA00029583860700000221
Figure FDA00029583860700000222
利用比值
Figure FDA00029583860700000223
估计天线波束在地面的入射角
Figure FDA00029583860700000224
Figure FDA00029583860700000225
其中，α(θ)为在不同入射角条件下已知理论比值。
2.根据权利要求1所述的方法，其中步骤(1)中正、负调频周期的回波信号
Figure FDA00029583860700000226
和
Figure FDA00029583860700000227
分别表示如下：
Figure FDA00029583860700000228
Figure FDA00029583860700000229
其中，Et为雷达发射信号的幅度，f0为中心频率，μ为调制斜率，
Figure FDA00029583860700000230
为随机初相位，τ为回波延时，R(t)为目标到雷达的距离，G为天线功率增益，λ为波长，σ为雷达截面积，t为时间，j为虚数单位。
3.根据权利要求1所述的方法，其中步骤3a)中分别对正、负调频周期的差频信号二维矩阵进行距离维和方位维的FFT处理，得到正、负调频周期差频信号频谱s+[fr(m),fa(n)]、s-[fr(m),fa(n)]，按如下步骤进行：
(3a1)对正、负调频周期差频信号二维矩阵进行距离维的FFT处理，得到距离维差频频谱s+[fr(m),ta(n)]、s-[fr(m),ta(n)]：
Figure FDA0002958386070000031
Figure FDA0002958386070000032
(3a2)对得到的距离维差频频谱进行方位维的FFT处理，得到正、负调频周期差频信号频谱s+[fr(m),fa(n)]、s-[fr(m),fa(n)]：
Figure FDA0002958386070000033
Figure FDA0002958386070000034
4.根据权利要求1所述的方法，其中步骤4a)中计算正频率回波包络的重心
Figure FDA0002958386070000035
和回波包络的宽度
Figure FDA0002958386070000036
按如下公式进行：
Figure FDA0002958386070000037
Figure FDA0002958386070000038
其中，
Figure FDA0002958386070000039
为一维的正调频周期差频信号频谱s'+[fr(m)]的峰值位置，Np为在一维的正调频周期差频信号频谱峰值两端取采样点的个数。
5.根据权利要求1所述的方法，其中步骤4b)中估计正调频周期噪声功率P1，按如下公式进行：
Figure FDA0002958386070000041
其中，Np1为
Figure FDA0002958386070000042
所在位置左侧采样点的个数。
6.根据权利要求1所述的方法，其中步骤4c)中计算负频率回波包络的重心
Figure FDA0002958386070000043
和回波包络的宽度
Figure FDA0002958386070000044
按如下公式进行：
Figure FDA0002958386070000045
Figure FDA0002958386070000046
其中，
Figure FDA0002958386070000047
为一维的负调频周期差频信号频谱s'-[fr(m)]的峰值位置，Np2为在一维的负调频周期差频信号频谱峰值两端取采样点的个数。
7.根据权利要求1所述的方法，其中步骤4d)中估计负调频周期噪声功率P2，按如下公式进行：
Figure FDA0002958386070000048
其中，Np3为
Figure FDA0002958386070000049
所在位置右侧采样点的个数。
8.根据权利要求1所述的方法，其中步骤5a)中计算正调频周期的回波包络重心
Figure FDA00029583860700000410
和回波包络宽度
Figure FDA00029583860700000411
按如下公式进行：
Figure FDA0002958386070000051
Figure FDA0002958386070000052
9.根据权利要求1所述的方法，其中步骤5b)中计算负调频周期的回波包络重心
Figure FDA0002958386070000053
和回波包络宽度
Figure FDA0002958386070000054
按如下公式进行：
Figure FDA0002958386070000055
Figure FDA0002958386070000056
10.根据权利要求1所述的方法，其中步骤(6)中计算回波重心对应的距离
Figure FDA0002958386070000057
和回波包络宽度
Figure FDA0002958386070000058
其计算公式如下：
Figure FDA0002958386070000059
Figure FDA00029583860700000510
其中，c为光速，μ为调制斜率。
CN201710952931.5A 2017-10-13 2017-10-13 调频连续波着陆雷达波束的地面入射角估计方法 Active CN107942307B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710952931.5A CN107942307B (zh) 2017-10-13 2017-10-13 调频连续波着陆雷达波束的地面入射角估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710952931.5A CN107942307B (zh) 2017-10-13 2017-10-13 调频连续波着陆雷达波束的地面入射角估计方法

Publications (2)

Publication Number Publication Date
CN107942307A CN107942307A (zh) 2018-04-20
CN107942307B true CN107942307B (zh) 2021-04-20

Family

ID=61935282

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710952931.5A Active CN107942307B (zh) 2017-10-13 2017-10-13 调频连续波着陆雷达波束的地面入射角估计方法

Country Status (1)

Country Link
CN (1) CN107942307B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109839623B (zh) * 2019-02-14 2020-09-11 北京遥感设备研究所 一种地外天体着陆测量雷达面目标回波信号测距处理方法
CN110988858B (zh) * 2019-11-11 2021-12-07 西安空间无线电技术研究所 一种双波束微波着陆雷达高精度测距方法及系统
CN113608205B (zh) * 2021-06-25 2024-04-09 西安空间无线电技术研究所 一种超高精度地外天体微波着陆雷达测距测速方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103630888A (zh) * 2013-02-27 2014-03-12 中国科学院电子学研究所 基于对称三角lfmcw雷达的高精度实时微波测速测距装置
CN104076353A (zh) * 2014-06-30 2014-10-01 西安空间无线电技术研究所 一种面目标回波波束中心速度测量方法
CN105629224A (zh) * 2016-02-29 2016-06-01 珠海上富电技股份有限公司 调频连续波雷达高精度测距方法
WO2016138238A1 (en) * 2015-02-27 2016-09-01 Texas Instruments Incorporated Gesture recognition using frequency modulated continuous wave radar
CN106646447A (zh) * 2017-01-18 2017-05-10 武汉雷博合创电子技术有限公司 基于线性调频连续波的雷达目标长时间积累检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103630888A (zh) * 2013-02-27 2014-03-12 中国科学院电子学研究所 基于对称三角lfmcw雷达的高精度实时微波测速测距装置
CN104076353A (zh) * 2014-06-30 2014-10-01 西安空间无线电技术研究所 一种面目标回波波束中心速度测量方法
WO2016138238A1 (en) * 2015-02-27 2016-09-01 Texas Instruments Incorporated Gesture recognition using frequency modulated continuous wave radar
CN105629224A (zh) * 2016-02-29 2016-06-01 珠海上富电技股份有限公司 调频连续波雷达高精度测距方法
CN106646447A (zh) * 2017-01-18 2017-05-10 武汉雷博合创电子技术有限公司 基于线性调频连续波的雷达目标长时间积累检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
一种应用于地形匹配雷达高度表中的回波跟踪与测高算法;马洪等;《宇航学报》;20051130;第26卷(第6期);第762-767页 *
月球着陆系统微波测距雷达波束中心估计算法研究;易予生等;《中国宇航学会深空探测技术专业委员会第七届学术年会论文集》;20101231;第327-330页 *
月面探测器回波模拟方法研究;朱保安;《中国优秀硕士学位论文全文数据库 信息科技辑》;20131215(第S2期);正文全文 *

Also Published As

Publication number Publication date
CN107942307A (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN1078949C (zh) 减轻脉冲多普勒雷达中距离一多普勒模糊的设备和方法
CN102998672B (zh) 基于相干化处理的步进频率isar成像方法
CN108387900B (zh) 一种直升机载旋转式合成孔径雷达的振动误差补偿方法
CN108802716B (zh) 基于重心校正的调频连续波着陆雷达测距方法
CN100541227C (zh) 子孔径雷达高度计
CN103454637B (zh) 基于调频步进频的太赫兹逆合成孔径雷达成像方法
CN107942307B (zh) 调频连续波着陆雷达波束的地面入射角估计方法
CN108535719B (zh) 基于多普勒频谱重心校正的调频连续波着陆雷达测速方法
CN106291524A (zh) 基于人体模型的lfmcw雷达探测运动人体目标的方法
CN106842166A (zh) 一种适用于线性调频连续波雷达系统的解速度模糊方法
CN102707269B (zh) 一种机载雷达距离走动校正方法
CN108490443B (zh) 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法
CN108132461B (zh) 抑制调频连续波着陆雷达直流泄露的方法
CN107966686B (zh) 一种基于线性调频的时间依赖性频控阵目标探测方法
CN109507669A (zh) 基于相参累积的地面微弱运动目标参数估计方法
CN108061882A (zh) 基于修正牛顿迭代的isar横向定标与越多普勒走动校正方法
CN102914775A (zh) 一种改进的双基地合成孔径雷达nlcs成像算法
CN110879391B (zh) 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法
CN114527430A (zh) 一种频率分组编码的捷变频抗干扰信号相参积累方法
CN103105610A (zh) 基于非均匀采样的dpc-mab sar成像方法
CN109884337B (zh) 一种利用高频地波雷达探测海面风向的方法
CN103064084A (zh) 基于距离频域的解模糊方法
CN109738893B (zh) 一种双基合成孔径雷达回波数据生成方法和装置
CN109856636B (zh) 曲线合成孔径雷达自适应三维成像方法
CN107271996A (zh) 一种机载cssar地面运动目标成像方法

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