CN106842150A - 一种简化的雷达地形模拟方法 - Google Patents
一种简化的雷达地形模拟方法 Download PDFInfo
- Publication number
- CN106842150A CN106842150A CN201510880668.4A CN201510880668A CN106842150A CN 106842150 A CN106842150 A CN 106842150A CN 201510880668 A CN201510880668 A CN 201510880668A CN 106842150 A CN106842150 A CN 106842150A
- Authority
- CN
- China
- Prior art keywords
- echo
- loudspeaker
- angle
- radar
- signal
- 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
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/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
Abstract
本发明提供了一种简化的雷达地形模拟方法,其特征在于,包括以下步骤:1)在雷达的远场位置设置喇叭1和喇叭2,喇叭1和喇叭2分别位于模拟地形的最高角度和最低角度2)模拟地形中位于最高角度和最低角度之间的所有角度回波D1,…,Di,…,Dn,通过喇叭1和喇叭2的回波合成模拟地形回波,则第i角度的回波Di为:Di=S1 i+S2 i其中,S1 i为喇叭1在第i角度回波,S2 i为喇叭2第i角度回波,设第1角度回波信号为雷达发射脉冲延迟后回波信号D1=Lp(m),m为回波的延迟距离门,m=(T*(2*f0)),T表示回波相对于发射信号的延时时间,单位为秒,f0为信号的采样率,则第i角度的回波距离门为m-i+1,第i角度的回波信号为:Di=Lp(m-i+1)=S1 i+S2 i=aiLp(m-i+1)+biLp(m-i+1)。
Description
技术领域
本发明属于雷达技术信号处理方向,涉及一种简化的空馈方式进行雷达的地形模拟的一种方法。
背景技术
本发明结合了空馈模拟不同角度目标的原理,来实现对原有地形模拟方式的改进。空馈模拟不同角度目标方法以及过零点检测方法基本原理介绍如下:
根据单脉冲测角原理,在模拟不同角度的目标时,可以通过两个喇叭进行不同加权进行模拟某个虚拟角度。暗室空馈方式虚拟某角度目标模拟方法见图 1。而对于地形模拟,一般的方法是在主波束内的地形分成N个区域,每个区域用等效散射中心表示,则各散射中心可以用不同的喇叭来模拟回波信号(图 2)。当等效散射中心角度偏差由负到正变化时,天线的相位差也由负到正变化。因此,不同的角度的回波在被雷达天线接收后,相位差有一个过零点位置,该位置对应的距离延时就为模拟的视轴方向的距离。由此便完成了地形过零点模拟。
现有的空馈模拟过零点存在的主要缺点是通过较多喇叭模拟不同角度的回波,则需要根据不同雷达、及雷达与喇叭距离来相应增减喇叭个数,并且每个喇叭的延迟和相位必须严格保持相应关系,实现成本较高。
发明内容
提出一种简单的空馈方式进行地形过零点模拟的方法,实现一种实用的空馈方式雷达对地形目标模拟。
技术方案
本发明结合暗室空馈方式模拟不同方向目标方法与空馈方式模拟过零点方法,从而通过两个喇叭来虚拟所有雷达主波束内的回波,达到简化布阵的目的,从而实现了简化的地形模拟的方法。
本发明一种简化的雷达地形模拟方法,其特征在于,包括以下步骤:
1)在雷达的远场位置设置喇叭1和喇叭2,喇叭1和喇叭2分别位于模拟地形的最高角度和最低角度,
2)模拟地形中位于最高角度和最低角度之间的所有角度回波D1,…,Di,…,Dn,通过喇叭1和喇叭2的回波合成模拟地形回波,则第i角度的回波Di为:
Di=S1 i+S2 i (1)
其中S1 i为喇叭1在第i角度回波,S2 i为喇叭2第i角度回波,
设第1角度回波信号为雷达发射脉冲延迟后回波信号D1=Lp(m),m为回波的延迟距离门,m=(T*(2*f0)),T表示回波相对于发射信号的延时时间,单位为秒,f0为信号的采样率,则第i角度的回波距离门为m-i+1,第i角度的回波信号为:
Di=Lp(m-i+1)=S1 i+S2 i=aiLp(m-i+1)+biLp(m-i+1) (2)
ai为喇叭1对相对于第i角度回波强度S1 i的加权系数,ai范围在0到1,bi为喇叭2对相对于第i角度回波强度S2 i的加权系数,bi范围在0到1;
3)将所有虚拟回波对应喇叭1的信号合并,得到喇叭1的总的回馈信号为:
令K1 i=ai,为信号幅度的加权,则可化简为如下公式:
其中:conv(x,y)表示x与y的卷积,同理得喇叭2的总的回馈信号S2:
S2=conv(K2,Lp) (5)
则通过S1、S2即能模拟多个方向回波,进而模拟地形过零点仿真。
在实际条件下,馈源喇叭俯仰角应满足主波束内,即雷达与馈源间距应满足远场条件L≤d2/λ;
其中:λ——雷达波长(单位m);
——馈源相对于雷达相位中心的俯仰角间距,单位为弧度;
d——雷达天线直径(单位m);
L——馈源与雷达距离,(单位m)。
本文提出的简化方式的地形模拟,推广到雷达试验暗室,能有效的降低试验室为进行地形模拟所需的布阵中馈源的个数,同时也满足充分的验证雷达天线在不同地形条件下的地形检测能力的要求。
附图说明
图 1暗室空馈方式虚拟某角度目标回波;
图 2本发明方法模拟地形过零点;
图 3实例示意图;
图 4线性调频信号示意图;
图 5加权系数;
图 6上方与下方喇叭输出信号。
具体实施方式
参考图 3到图 6,一种简化的雷达地形模拟方法,其特征在于,包括以下步骤:
1)在雷达的远场位置设置喇叭1和喇叭2,喇叭1和喇叭2分别位于模拟地形的最高角度和最低角度。
2)模拟地形中位于最高角度1°和最低角度-1°之间的间隔0.02°所有101个角度回波D1,…,Di,…,D101,通过喇叭1和喇叭2的回波合成模拟地形回波,则第i角度的回波Di为:
Di=S1 i+S2 i (1)
S1 i为喇叭1在第i角度回波,S2 i为喇叭2第i角度回波;
设第1角度回波信号为雷达发射脉冲延迟后回波信号D1=Lp(m),此时,Lp为带宽为25MHz,脉宽为4.5us的线性调频信号。
其中,m为回波的延迟距离门,回波D1相对于发射信号的延时时间为11us,f0为100MHz,m=(11*10-6*(2*100*106)),m=2200,即从2200点后开始接收到D1回波;
令K1,K2为信号幅度的加权,加权参数为K1=[1 0.99 0.98 0.97 … 0],K2=[0 0.01 0.02 0.03 … 1],通过如下公式(5)和(6)进行S1、S2计算:
S1=conv(K1,Lp) (5)
S2=conv(K2,Lp) (6)
则通过S1、S2即能模拟多个方向回波,进而模拟地形过零点仿真。
已知雷达波长λ为0.0086m,雷达直径d为0.4m。则馈源与雷达间距为19m,角度为0.02弧度;
如 图 3,两个喇叭分别放在雷达主波束前方俯仰0.01和-0.01弧度的位置,设采样率为100MHz,发射信号为宽度为4.5us,带宽为25MHz的线性调频信号(图 4),模拟个斜面,在主波束内,斜面与雷达径向最近距离为1500m,长度为150m,斜面角度为雷达主波束宽度,斜面在主波束内共有101个散射中心。各个点的值如下表:
表 1地形各散射中心参数
地形回波序号i | 延时 | ||
1 | 10us+1000ns | 1 | 0 |
2 | 10us+990ns | 0.98 | 0.01 |
… | … | … | … |
101 | 10us | 0 | 1 |
因此,在加权系数为表中的K1、K2条件下,可以通过公式(5)和(6)计算出两个喇叭的回波信号(图 6)虚拟这101个方向的角度回波。从而得到上下喇叭的输出信号S1和S2(图 5),完成地形模拟。
Claims (2)
1.一种简化的雷达地形模拟方法,其特征在于,包括以下步骤:
1)在雷达的远场位置设置喇叭1和喇叭2,喇叭1和喇叭2分别位于模拟地形的最高角度和最低角度,
2)模拟地形中位于最高角度和最低角度之间的所有角度回波D1,…,Di,…,Dn,通过喇叭1和喇叭2的回波合成模拟地形回波,则第i角度的回波Di为:
Di=S1 i+S2 i (1)
其中,S1 i为喇叭1在第i角度回波,S2 i为喇叭2第i角度回波,
设第1角度回波信号为雷达发射脉冲延迟后回波信号D1=Lp(m),m为回波的延迟距离门,m=(T*(2*f0)),T表示回波相对于发射信号的延时时间,单位为秒,f0为信号的采样率,则第i角度的回波距离门为m-i+1,第i角度的回波信号为:
Di=Lp(m-i+1)=S1 i+S2 i=aiLp(m-i+1)+biLp(m-i+1) (2)
ai为喇叭1对相对于第i角度回波强度S1 i的加权系数,ai范围在0到1,bi为喇叭2对相对于第i角度回波强度S2 i的加权系数,bi范围在0到1;
将所有虚拟回波对应喇叭1的信号合并,得到喇叭1的总的回馈信号为:
令K1 i=ai,为信号幅度的加权,则可化简为如下公式:
其中:conv(x,y)表示x与y的卷积,同理得喇叭2的总的回馈信号S2:
S2=conv(K2,Lp) (5)
则通过S1、S2即能模拟多个方向回波,进而模拟地形过零点仿真。
2.根据权利要求1所述的简化的雷达地形模拟方法,其特征在于,包括以下步骤:
馈源喇叭俯仰角应满足主波束内,即雷达与馈源间距应满足远场条件L≤d2/λ;
其中:λ——雷达波长(单位m);
——馈源相对于雷达相位中心的俯仰角间距,单位为弧度;
d——雷达天线直径(单位m);
L——馈源与雷达距离,(单位m)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510880668.4A CN106842150B (zh) | 2015-12-03 | 2015-12-03 | 一种简化的雷达地形模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510880668.4A CN106842150B (zh) | 2015-12-03 | 2015-12-03 | 一种简化的雷达地形模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106842150A true CN106842150A (zh) | 2017-06-13 |
CN106842150B CN106842150B (zh) | 2020-04-28 |
Family
ID=59149889
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510880668.4A Active CN106842150B (zh) | 2015-12-03 | 2015-12-03 | 一种简化的雷达地形模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106842150B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050052714A1 (en) * | 2003-07-24 | 2005-03-10 | Zebra Imaging, Inc. | Enhanced environment visualization using holographic stereograms |
CN101482616A (zh) * | 2008-08-13 | 2009-07-15 | 中国科学院电子学研究所 | 一种地形测量方法 |
CN102707284A (zh) * | 2012-06-25 | 2012-10-03 | 西安煤航信息产业有限公司 | 一种基于机载InSAR的地面控制测量布点方法 |
CN104949629A (zh) * | 2015-04-20 | 2015-09-30 | 同济大学 | 一种基于激光测距的危岩变形多点信息提取及报警系统 |
-
2015
- 2015-12-03 CN CN201510880668.4A patent/CN106842150B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050052714A1 (en) * | 2003-07-24 | 2005-03-10 | Zebra Imaging, Inc. | Enhanced environment visualization using holographic stereograms |
CN101482616A (zh) * | 2008-08-13 | 2009-07-15 | 中国科学院电子学研究所 | 一种地形测量方法 |
CN102707284A (zh) * | 2012-06-25 | 2012-10-03 | 西安煤航信息产业有限公司 | 一种基于机载InSAR的地面控制测量布点方法 |
CN104949629A (zh) * | 2015-04-20 | 2015-09-30 | 同济大学 | 一种基于激光测距的危岩变形多点信息提取及报警系统 |
Non-Patent Citations (1)
Title |
---|
张朋 等: "合成孔径雷达成像三维地形目标模拟方法", 《系统仿真学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN106842150B (zh) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102346249B (zh) | 合成孔径雷达宽测绘带对地观测步进扫描模式的实现方法 | |
CN105259552A (zh) | 一种基于非线性调频信号的合成孔径雷达成像方法和装置 | |
CN105445701B (zh) | Ddma‑mimo雷达目标的单脉冲角度估计方法 | |
CN104035095A (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
CN107728109A (zh) | 一种非合作目标辐射噪声测量定位技术 | |
CN104614713A (zh) | 一种适合于艇载雷达系统的雷达回波信号模拟器 | |
CN106501802A (zh) | 高分辨率多维协同昆虫迁飞雷达测量仪 | |
CN107976660A (zh) | 弹载多通道雷达超低空目标分析与多径回波建模方法 | |
CN105738891B (zh) | 一种机载数字阵列雷达对弱机动目标角度进行跟踪的方法 | |
CN102914775B (zh) | 一种改进的双基地合成孔径雷达nlcs成像算法 | |
CN104391279A (zh) | 基于电离层传播特性的相径扰动抑制方法 | |
CN102435992A (zh) | 基于广义相干系数的合成聚焦成像方法 | |
CN107607936A (zh) | 一种高频天地波雷达海洋表面流反演方法 | |
CN105738887A (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
CN109520383B (zh) | 基于matlab的体目标回波模拟方法 | |
CN104076337B (zh) | 基于阵元幅相误差校正的机载雷达杂波抑制方法 | |
CN104635214A (zh) | 机载前视频率分集阵列雷达距离模糊杂波抑制方法 | |
CN103267964A (zh) | 一种基于低秩矩阵恢复的弹载导引头雷达σδ-stap方法 | |
CN104200110A (zh) | 基于去耦合的二维波束扫描测角方法 | |
CN102788970B (zh) | 机载数字阵列雷达近程杂波高效自适应抑制方法 | |
CN106291547B (zh) | 基于天线方向图辅助的多普勒模糊分量自适应抑制方法 | |
CN112098999A (zh) | 一种高动态雷达导引头掠海目标电磁信号建模方法 | |
CN106054187A (zh) | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 | |
CN109884621B (zh) | 雷达高度表回波相参积累方法 | |
CN104020465B (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 |