CN111337890A - 一种lfmcw雷达目标回波信号模拟方法 - Google Patents

一种lfmcw雷达目标回波信号模拟方法 Download PDF

Info

Publication number
CN111337890A
CN111337890A CN202010100663.6A CN202010100663A CN111337890A CN 111337890 A CN111337890 A CN 111337890A CN 202010100663 A CN202010100663 A CN 202010100663A CN 111337890 A CN111337890 A CN 111337890A
Authority
CN
China
Prior art keywords
signal
frequency
control word
doppler
dds
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
Application number
CN202010100663.6A
Other languages
English (en)
Other versions
CN111337890B (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202010100663.6A priority Critical patent/CN111337890B/zh
Publication of CN111337890A publication Critical patent/CN111337890A/zh
Application granted granted Critical
Publication of CN111337890B publication Critical patent/CN111337890B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4052Means for monitoring or calibrating by simulation of echoes
    • G01S7/4056Means for monitoring or calibrating by simulation of echoes specially adapted to FMCW
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

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

本发明公开了一种LFMCW雷达目标回波信号模拟方法,首先,通过Matlab产生所需的对称三角LFMCW雷达对应的基带信号、多普勒DDS频率控制字、延时控制字、上变频DDS频率控制字、杂波与噪声信号,并存入FPGA芯片的ROM中;其次,对基带信号进行多普勒调制;然后,对多普勒调制后的信号进行延时处理;在系统中加入杂波和噪声;最后,对信号进行上变频,输出模拟回波信号。本发明采用软硬结合的方法,使得控制更加灵活,系统更加集成化,能更方便地产生符合要求的回波信号。

Description

一种LFMCW雷达目标回波信号模拟方法
技术领域
本发明属于雷达检测与信号处理技术领域,具体涉及一种LFMCW雷达目标回波信号模拟方法。
背景技术
在现代雷达系统的设计、研发、调试及使用过程中,雷达系统的测试是一个十分重要的环节,但是如果都采用真实目标来为待测雷达创造检验环境,不但会耗费大量的人力、物力与财力,而且不一定都能实现。另外在真实的外场测试中,由于各种情况都可能存在干扰,将会使得测试的可重复性差且控制较为复杂。雷达回波模拟技术的实现解决了外场验证中存在的大量问题,缩短了开发时间、节约了雷达系统的开发经费,并且提高了可测试目标样本的数量,降低了风险系数,从而更有助于完善雷达系统。对于雷达回波信号模拟系统的设计,一般采用全软件实现或全硬件实现。全硬件实现的硬件复杂度高,对硬件处理数据的性能要求高,且灵活性低,难以进行多种波形的模拟;全软件实现方法实时性不强,主要运用在雷达系统的分析及论证中。本发明采用软硬件结合的方法,结合计算机软件和现场可编程门阵列(Field Programmable Gate Array,FPGA)能够方便地对模拟的回波参数进行控制,缩短硬件的数据处理时间,加快系统处理速度,提高系统的灵活性与实时性,增大系统容量,提高目标回波模拟的有效性。
发明内容
发明目的:提供了一种LFMCW雷达目标回波信号模拟方法,能模拟出包含所需目标信息的回波信号,缩短硬件的数据处理时间,提高系统的灵活性和实时性,提高系统集成化,增大系统容量。
技术方案:本发明所述的一种LFMCW雷达目标回波信号模拟方法,包括以下步骤:
(1)通过Matlab产生所需的对称三角LFMCW雷达对应的基带信号、多普勒DDS频率控制字、延时控制字、上变频DDS频率控制字、杂波与噪声信号,并存入FPGA芯片的ROM中;
(2)对基带信号进行多普勒调制;
(3)对多普勒调制后的信号进行延时处理;
(4)在系统中加入杂波和噪声;
(5)对信号进行上变频,输出模拟回波信号。
进一步地,所述步骤(2)包括以下步骤:
(21)根据式子
Figure BDA0002386417690000021
计算出多普勒频移,其中fd为所求多普勒频率,v为模拟目标的速度,c为无线电波在均匀介质中传播的速度,可认为为光速,f0为LFMCW雷达发射信号载波的频率,对目标速度的模拟就是将速度对应的多普勒频移加载到基带信号的载频上;
(22)通过式
Figure BDA0002386417690000022
计算多普勒DDS的频率控制字,其中fout为DDS输出的正余弦波信号频率,M为频率控制字,fclk为时钟频率,N为相位累加器的位数,根据fout=fd即可求得频率控制字M1,将频率控制字M1存入FPGA中;
(23)使用相位累加器和正余弦查询表构造DDS模块,控制频率控制字M1作为相位累加器的输入端,控制DDS输出信号的频率,相位累加器的输出作为正余弦查找表的输入,其中余弦查找表中存入的是一个完整周期的余弦信号,正弦查找表中存入的为一个完整周期正弦函数取负数的结果,余弦查询表和正弦查询表的输出分别为I路与Q路的多普勒载波信号,设为
Figure BDA0002386417690000023
Figure BDA0002386417690000024
各模块时钟均为系统时钟;
(24)将基带信号SI和SQ分别与多普勒载波信号
Figure BDA0002386417690000025
Figure BDA0002386417690000026
通过有符号数乘法器进行交互相乘,得到
Figure BDA0002386417690000027
Figure BDA0002386417690000028
(25)通过有符号数减法器与加法器,将S11(t)减去S22(t)得
Figure BDA0002386417690000029
将S12(t)加上S21(t)得
Figure BDA00023864176900000210
Figure BDA00023864176900000211
Figure BDA00023864176900000212
即为多普勒调制后的基带信号,他们的载频为fd
进一步地,所述步骤(3)包括以下步骤:
(31)通过公式
Figure BDA0002386417690000031
计算出所需距离对应的延时,这里R为模拟目标相对于雷达天线的距离,根据延时和系统时钟频率可计算出该延时对应的系统时钟脉冲的个数D,D即为延时控制字,将延时控制字存入FPGA中;
(32)创建一个D位计数器,延时控制字决定计数器的最大位数,计数器的输出使能端初值设为0,当计数器计数小于D时,输出使能端始终为0,当计数器计数到D时,在下一个时钟周期来临时输出使能端变为1;
(33)将多普勒已调信号
Figure BDA0002386417690000033
Figure BDA0002386417690000034
分别通过两个同步FIFO存储器,同步FIFO存储器有四个输入端,分别为信号输入端、写使能端、读使能端和时钟端口,其中写使能端和读使能端分别控制FIFO存储器的写入和读出,均为高电平有效,FIFO存储器的写使能端与系统的复位模块输出端相连,复位模块输出1时系统工作,即系统开始工作后FIFO模块即可同步写入,读使能端连接到D位计数器的输出端,如此,通过控制延时控制字D便可以控制FIFO存储器延迟D个时钟脉冲读出。
进一步地,所述步骤(4)实现过程如下:
采用Matlab生成满足需要的杂波和噪声序列,采用Q8方法为杂波和噪声分别定标,将定标后的杂波和噪声信号分别存入FPGA的ROM中,在系统时钟的控制下对ROM中的杂波和噪声信号进行读取,并将读取出的数据通过有符号数加法器与延时后的信号进行相加。
进一步地,所述步骤(5)包括以下步骤:
(51)采用CIC插值滤波器对信号进行P倍插值,来提高信号的采样频率,插值后信号的带宽由B变为B/P,B/P范围以外的P-1个频谱为镜像频谱;
(52)将CIC插值滤波器的输出信号通过一个截止频率为B/P的低通滤波器,滤除插值产生的镜像频谱;
(53)通过式
Figure BDA0002386417690000032
计算上变频DDS的频率控制字,根据fout=fI,fI为上变频载波频率,即可求得频率控制字M2,将频率控制字M2存入FPGA中;
(54)使用相位累加器和正余弦查询表构造DDS模块,控制频率控制字M2作为相位累加器的输入端,控制DDS输出信号的频率,相位累加器的输出作为正余弦查找表的输入,其中余弦查找表和正弦查找表中存入的分别是一个完整周期的余弦信号和正弦信号,余弦查询表和正弦查询表的输出分别为I路与Q路的上变频载波信号,设为Iduc和Qduc,各模块时钟均为系统时钟;
(55)通过有符号数乘法器,将经过各项调制后的I路信号与I路载波信号Iduc相乘,得到SIf,经过各项调制后的Q路信号与Q路载波信号Qduc相乘,得到SQf,通过有符号数减法器用SIf减去SQf得SU,SU即为上变频的输出信号。
有益效果:与现有技术相比,本发明的有益效果:1、采用软硬件结合的方法对LFMCW雷达回波信号进行模拟,能够减少硬件的计算量,使控制更加灵活,更方便的产生所需的回波信号;2、采用FPGA对LFMCW雷达回波信号硬件处理,使得系统的体积更小,运算速度更快,系统容量更大,可运用此方法实现多种情况下的目标回波信号模拟。
附图说明
图1为本发明的流程图;
图2为基带信号存储模块原理图;
图3为DDS模块原理图;
图4为多普勒调制模块原理图;
图5为延时模块原理图;
图6为背景信号模块原理图;
图7为插值模块原理图;
图8为上变频原理图;
图9为目标一差拍信号正调制段对应频谱图;
图10为目标一差拍信号负调制段对应频谱图;
图11为目标二差拍信号正调制段对应频谱图;
图12为目标二差拍信号负调制段对应频谱图。
具体实施方式
下面结合附图对本发明作进一步详细描述。
基于FPGA芯片的线性调频连续波(LFMCW)雷达回波信号模拟系统的实现方法,主要包括软件实现和硬件实现两大部门。软件实现主要包括通过Matlab产生所需的对称三角LFMCW雷达对应的基带信号、多普勒DDS频率控制字、延时控制字、上变频DDS频率控制字、杂波与噪声信号;硬件实现主要包括基带信号存储模块、多普勒调制模块、延时模块、背景信号模块、上变频模块、复位模块及时钟模块。具体包括以下步骤:
1、通过Matlab产生所需的对称三角LFMCW雷达对应的基带信号、多普勒DDS频率控制字、延时控制字、上变频DDS频率控制字、杂波与噪声信号,并存入FPGA芯片的ROM中。
采用Matlab仿真出I路与Q路两路基带信号,由于LFMCW雷达信号的正调制段与负调制段的数字处理过程相同,可将正调制段与负调制段的基带信号合并为一个信号,并采用Q8法对合并后的基带信号定标,将定标后的I路与Q路基带信号通过mif文件分别存入到FPGA的两个ROM中,创建系统时钟与地址模块对ROM进行重复读取,设I路基带信号为SI,Q路基带信号为SQ
2、对基带信号进行多普勒调制。
(1)根据式子
Figure BDA0002386417690000051
计算出多普勒频移,其中fd为所求多普勒频率,v为模拟目标的速度,c为无线电波在均匀介质中传播的速度,可认为为光速,f0为LFMCW雷达发射信号载波的频率,对目标速度的模拟就是将速度对应的多普勒频移加载到基带信号的载频上。
(2)通过式
Figure BDA0002386417690000052
计算多普勒DDS的频率控制字,其中fout为DDS输出的正余弦波信号频率,M为频率控制字,fclk为时钟频率,N为相位累加器的位数,常取N=24,根据fout=fd即可求得频率控制字M1,将频率控制字M1存入FPGA中。
(3)使用相位累加器和正余弦查询表构造DDS模块,控制频率控制字M1作为相位累加器的输入端,控制DDS输出信号的频率,相位累加器的输出作为正余弦查找表的输入,其中余弦查找表中存入的是一个完整周期的余弦信号,正弦查找表中存入的为一个完整周期正弦函数取负数的结果,余弦查询表和正弦查询表的输出分别为I路与Q路的多普勒载波信号,设为
Figure BDA0002386417690000053
Figure BDA0002386417690000054
各模块时钟均为系统时钟。
(4)将基带信号SI和SQ分别与多普勒载波信号
Figure BDA0002386417690000061
Figure BDA0002386417690000062
通过有符号数乘法器进行交互相乘,得到
Figure BDA0002386417690000063
Figure BDA0002386417690000064
(5)通过有符号数减法器与加法器,将S11(t)减去S22(t)得
Figure BDA0002386417690000065
将S12(t)加上S21(t)得
Figure BDA0002386417690000066
Figure BDA0002386417690000067
Figure BDA0002386417690000068
即为多普勒调制后的基带信号,他们的载频为fd
3、对多普勒调制后的信号进行延时处理。
(1)通过公式
Figure BDA0002386417690000069
计算出所需距离对应的延时,这里R为模拟目标相对于雷达天线的距离,根据延时和系统时钟频率可计算出该延时对应的系统时钟脉冲的个数D,D即为延时控制字,将延时控制字存入FPGA中。
(2)创建一个D位计数器,延时控制字决定计数器的最大位数,计数器的输出使能端初值设为0,当计数器计数小于D时,输出使能端始终为0,当计数器计数到D时,在下一个时钟周期来临时输出使能端变为1。
(3)将多普勒已调信号
Figure BDA00023864176900000610
Figure BDA00023864176900000611
分别通过两个同步FIFO存储器,同步FIFO存储器有四个输入端,分别为信号输入端、写使能端、读使能端和时钟端口,其中写使能端和读使能端分别控制FIFO存储器的写入和读出,均为高电平有效,FIFO存储器的写使能端与系统的复位模块输出端相连,复位模块输出1时系统工作,即系统开始工作后FIFO模块即可同步写入,读使能端连接到D位计数器的输出端,如此,通过控制延时控制字D便可以控制FIFO存储器延迟D个时钟脉冲读出。
4、在系统中加入杂波和噪声。
采用Matlab生成满足需要的杂波和噪声序列,采用Q8方法为杂波和噪声分别定标,将定标后的杂波和噪声信号分别存入FPGA的ROM中,在系统时钟的控制下对ROM中的杂波和噪声信号进行读取,并将读取出的数据通过有符号数加法器与延时后的信号进行相加。
5、对信号进行上变频,输出模拟回波信号。
(1)采用CIC插值滤波器对信号进行P倍插值,来提高信号的采样频率,插值后信号的带宽由B变为B/P,B/P范围以外的P-1个频谱为镜像频谱。
(2)将CIC插值滤波器的输出信号通过一个截止频率为B/P的低通滤波器,滤除插值产生的镜像频谱。
(3)通过式
Figure BDA0002386417690000071
计算上变频DDS的频率控制字,根据fout=fI即可求得频率控制字M2,将频率控制字M2存入FPGA中。
(4)使用相位累加器和正余弦查询表构造DDS模块,控制频率控制字M2作为相位累加器的输入端,控制DDS输出信号的频率,相位累加器的输出作为正余弦查找表的输入,其中余弦查找表和正弦查找表中存入的分别是一个完整周期的余弦信号和正弦信号,余弦查询表和正弦查询表的输出分别为I路与Q路的上变频载波信号,设为Iduc和Qduc,各模块时钟均为系统时钟。
(5)通过有符号数乘法器,将经过各项调制后的I路信号与I路载波信号Iduc相乘,得到SIf,经过各项调制后的Q路信号与Q路载波信号Qduc相乘,得到SQf,通过有符号数减法器用SIf减去SQf得SU,SU即为上变频的输出信号。
本发明在基于FPGA实现的基础上,采用软件与硬件相结合的方法,由软件产生部分回波模拟所需要的信号及参数并预先存入到FPGA中,再通过FPGA对数据进行处理,实现LFMCW雷达回波信号的模拟。
假设LFMCW雷达回波信号模拟的参数如表1所示。
表1目标参数
Figure BDA0002386417690000072
Figure BDA0002386417690000081
参考表1的参数,对目标一和目标二进行LFMCW目标回波信号模拟,并验证模拟回波的正确性,目标回波模拟的总体流程图如图1所示。首先在Matlab中生成基带信号并进行定标,设基带信号为S(t),当数据S(t)≥0时,S(t)量化后的数据S(k)为:
S(k)=fix(S(t)·2W-1)
当数据S(t)<0时,S(t)量化后的数据S(k)为:
S(k)=fix(S(t)·2W-1+2W)
这里W为量化位宽,本文中W=8,fix函数为向零方向取整。
将基带信号存入FPGA内部存储器的方法为:建立一个数据位宽为8bit、数据长度(设为L)等于基带信号定标点数的mif文件,将定标后的信号存入mif文件中并保存该mif文件。调用ROM存储器的IP核,将该ROM的存储文件设置为保存的mif文件,并保存。建立一个地址循环模块address.v,该模块的作用是循环输出0~L-1,将地址读取模块的输出作为ROM核的地址端的输入。地址读取模块与ROM模块的触发时钟均为同一个时钟,该时钟的频率等于Matlab中基带信号的采样频率。当时钟触发后,ROM模块即可完成每一个时钟上升沿读取一个数据,如此循环。基带信号存储模块的原理图如图2所示。
对基带信号进行多普勒调制,其中DDS模块的原理图如图3所示。其中相位累加器可在每个时钟周期来临时将频率控制字所决定的相位增量M累加一次,当计数大于2N时则自动溢出,只保留后面的N位数字至累加器中。正弦查询表ROM中存放了一个保存一个完整正/余弦周期的mif文件,相位累加器的输出作为ROM的地址端,通过控制频率控制字的大小即可调节读取ROM中一个完整周期的正/余弦波的时间,即控制了输出信号的频率。
多普勒调制模块的原理图如图4所示。以目标远离雷达的方向为正方向,则
Figure BDA0002386417690000091
其中f0为载频,则由DDS模块产生多普勒载波信号
Figure BDA0002386417690000092
Figure BDA0002386417690000093
可以表示为
Figure BDA0002386417690000094
Figure BDA0002386417690000095
其中
Figure BDA0002386417690000096
可以看出,DDS余弦查找表中存入的是一个完整周期的余弦信号,而DDS正弦查找表中存入的为一个完整周期正弦函数取负数的结果。DDS模块产生的多普勒信号需要将其与I路、Q路的基带信号一一对应相乘。由于基带信号的数值有正负之分,所以该乘法模块对应的为有符号乘法器,可以通过调用乘法器的IP核,并将输入输出信号均设置为signed形式,即可完成有符号数信号的相乘。
再对信号进行延时处理,延时模块的原理图如图5所示。写使能端控制信号的存储,当该端口为高电平时,FIFO存储器开始按照数据的先后顺序存储数据;读使能端控制信号的输出,当该端口为高电平时FIFO存储器由先入先出的方式开始读取数据。通过控制FIFO存储器的读使能端与写使能端即可控制信号在FIFO存储器中的缓存时间。定时计数器的最大计数值可以称为延时控制字。设置延时控制字时还需要考虑系统各个硬件模块工作产生的延时大小。最终的延时控制字WordR应为延时τ所对应的时钟周期数减去硬件模块工作产生的延时周期数。
背景信号模块主要包括距离电压系数调制、杂波和噪声调制。背景信号模块参考基带信号存储模块,以杂波为例,背景信号模块的原理图如图6所示。
数字上变频模块的主要是由三部分组成:插值、滤波、正交上变频。在FPGA硬件实现中,CIC插值滤波器可以通过调用IP核生成。配置主要有:
Filter Type:设置CIC滤波器的模式为抽取Dcimator或插值Interpolator;
Number ofStages:设置级联的CIC滤波器的级数,根据插值因子和系统资源来决定。本文插值因子为5,选择CIC滤波器的级数为4;
Differential delay:设置差分延时,该值与输出数据位宽又关,本文输出数据位宽采取阶段方式,查表可知差分延时设置为1;
Rate change factor:设置抽取/插值因子,本文中设置插值因子为5。
经过插值后的信号还需进行低通滤波,可以通过调用IP核生成滤波器。调用滤波器IP核需要设置的主要参数有:滤波器的类型、输入信号位宽、输出信号位宽、滤波器、阶数、窗类型、采样速率及截止频率。采用的滤波器类型为低通滤波器,输入信号位宽为16bit,输出信号位宽截取信号的高16位,也为16bit,窗类型为Blackman窗,阶数设置为127阶,采样速率为360MHz,截止频率为10MHz。插值滤波及上变频的原理图如图7和图8所示。
为验证模拟出的回波信号是否包含正确的距离及速度信息,采用Signal Tap抓取FPGA输出的信号,并联合Matlab读取。I路延时后的信号与模拟出的回波信号相比形式相同,只缺少距离电压系数、杂波与上变频载频部分,但都包含了目标的距离与速度信息,所以可以利用I路延时后的信号与I路基带信号来求出该LFMCW雷达的差拍信号,由此计算出目标的距离与速度信息。
当目标参数为距离R=150m,速度v=90m/s时,由Signal Tap抓取出I路延时后的信号SI_delay,将SI_delay作为回波信号,I路基带信号SI作为发射信号。将回波信号不变,发射信号取共轭,相乘即可求得差拍信号,即S_chapai=SI·SI_delay *,其中(·)*表示取共轭。
通过对差拍信号进行频谱分析可以得出正调制段对应的差拍信号频谱峰值fb +,负调制段对应的差拍信号频谱峰值fb -
Figure BDA0002386417690000101
Figure BDA0002386417690000102
通过fb +和fb -即可以求出该回波信号对应的目标的距离信息与速度信息,即
Figure BDA0002386417690000103
Figure BDA0002386417690000104
求得的目标一对应的差拍信号频谱如图9和图10所示。
由图9和图10可以看出目标一正调制段对应的差拍信号频率为fb +=0.5156MHz,负调制段对应的差拍信号频率为fb -=0.4688MHz。经过Rife算法修正后,正调制段对应的差拍信号频率为fb +=0.52172MHz,负调制段对应的差拍信号频率为fb -=0.47853MHz,带入fb +和fb -的计算公式中求得距离R=150.0371m,速度V=92.5342m。
同理当目标参数为距离R=600m,速度V=200m/s时,采用同样的方法求出差拍信号,并通过正调制段与负调制段对应的差拍信号频率值求出估计的目标参数,结果如图11和图12所示。由图11和图12可以看出目标二正调制段对应的差拍信号频率为fb +=2.063MHz,负调制段对应的差拍信号频率为fb-=1.969MHz。经过Rife算法修正后,正调制段对应的差拍信号频率为fb +=2.0455MHz,负调制段对应的差拍信号频率为fb-=1.9523MHz,带入fb +和fb-的计算公式中求得距离R=599.671m,速度V=199.6344m/s。
由目标一和目标二的验证结果可知,估算出的目标参数与目标的真实参数十分接近,FPGA芯片模拟出的LFMCW雷达目标回波信号满足设计要求。
综上可知,该LFMCW雷达回波模拟系统模拟可对目标的回波信号进行准确模拟,模拟出的回波信号包含正确的距离、速度信息,满足设计要求。
通过以上多组验证可以得出,可对目标的回波信号进行准确模拟,并且可以灵活地控制或改变基带信号、频率控制字、延时控制字等,提高了系统的灵活性和实时性,系统的硬件处理部分仅通过FPGA即可完成,提高了硬件的处理速度,减小了系统的体积,集成度更高。

Claims (5)

1.一种LFMCW雷达目标回波信号模拟方法,其特征在于,包括以下步骤:
(1)通过Matlab产生所需的对称三角LFMCW雷达对应的基带信号、多普勒DDS频率控制字、延时控制字、上变频DDS频率控制字、杂波与噪声信号,并存入FPGA芯片的ROM中;
(2)对基带信号进行多普勒调制;
(3)对多普勒调制后的信号进行延时处理;
(4)在系统中加入杂波和噪声;
(5)对信号进行上变频,输出模拟回波信号。
2.根据权利要求1所述的一种LFMCW雷达目标回波信号模拟方法,其特征在于,所述步骤(2)包括以下步骤:
(21)根据式子
Figure FDA0002386417680000011
计算出多普勒频移,其中fd为所求多普勒频率,v为模拟目标的速度,c为无线电波在均匀介质中传播的速度,可认为为光速,f0为LFMCW雷达发射信号载波的频率,对目标速度的模拟就是将速度对应的多普勒频移加载到基带信号的载频上;
(22)通过式
Figure FDA0002386417680000012
计算多普勒DDS的频率控制字,其中fout为DDS输出的正余弦波信号频率,M为频率控制字,fclk为时钟频率,N为相位累加器的位数,根据fout=fd即可求得频率控制字M1,将频率控制字M1存入FPGA中;
(23)使用相位累加器和正余弦查询表构造DDS模块,控制频率控制字M1作为相位累加器的输入端,控制DDS输出信号的频率,相位累加器的输出作为正余弦查找表的输入,其中余弦查找表中存入的是一个完整周期的余弦信号,正弦查找表中存入的为一个完整周期正弦函数取负数的结果,余弦查询表和正弦查询表的输出分别为I路与Q路的多普勒载波信号,设为
Figure FDA0002386417680000013
Figure FDA0002386417680000014
各模块时钟均为系统时钟;
(24)将基带信号SI和SQ分别与多普勒载波信号
Figure FDA0002386417680000015
Figure FDA0002386417680000016
通过有符号数乘法器进行交互相乘,得到
Figure FDA0002386417680000017
Figure FDA0002386417680000018
(25)通过有符号数减法器与加法器,将S11(t)减去S22(t)得
Figure FDA0002386417680000021
将S12(t)加上S21(t)得
Figure FDA0002386417680000022
Figure FDA0002386417680000023
即为多普勒调制后的基带信号,他们的载频为fd
3.根据权利要求1所述的一种LFMCW雷达目标回波信号模拟方法,其特征在于,所述步骤(3)包括以下步骤:
(31)通过公式
Figure FDA0002386417680000024
计算出所需距离对应的延时,这里R为模拟目标相对于雷达天线的距离,根据延时和系统时钟频率可计算出该延时对应的系统时钟脉冲的个数D,D即为延时控制字,将延时控制字存入FPGA中;
(32)创建一个D位计数器,延时控制字决定计数器的最大位数,计数器的输出使能端初值设为0,当计数器计数小于D时,输出使能端始终为0,当计数器计数到D时,在下一个时钟周期来临时输出使能端变为1;
(33)将多普勒已调信号
Figure FDA0002386417680000025
Figure FDA0002386417680000026
分别通过两个同步FIFO存储器,同步FIFO存储器有四个输入端,分别为信号输入端、写使能端、读使能端和时钟端口,其中写使能端和读使能端分别控制FIFO存储器的写入和读出,均为高电平有效,FIFO存储器的写使能端与系统的复位模块输出端相连,复位模块输出1时系统工作,即系统开始工作后FIFO模块即可同步写入,读使能端连接到D位计数器的输出端,如此,通过控制延时控制字D便可以控制FIFO存储器延迟D个时钟脉冲读出。
4.根据权利要求1所述的一种LFMCW雷达目标回波信号模拟方法,其特征在于,所述步骤(4)实现过程如下:
采用Matlab生成满足需要的杂波和噪声序列,采用Q8方法为杂波和噪声分别定标,将定标后的杂波和噪声信号分别存入FPGA的ROM中,在系统时钟的控制下对ROM中的杂波和噪声信号进行读取,并将读取出的数据通过有符号数加法器与延时后的信号进行相加。
5.根据权利要求1所述的一种LFMCW雷达目标回波信号模拟方法,其特征在于,所述步骤(5)包括以下步骤:
(51)采用CIC插值滤波器对信号进行P倍插值,来提高信号的采样频率,插值后信号的带宽由B变为B/P,B/P范围以外的P-1个频谱为镜像频谱;
(52)将CIC插值滤波器的输出信号通过一个截止频率为B/P的低通滤波器,滤除插值产生的镜像频谱;
(53)通过式
Figure FDA0002386417680000031
计算上变频DDS的频率控制字,根据fout=fI,fI为上变频载波频率,即可求得频率控制字M2,将频率控制字M2存入FPGA中;
(54)使用相位累加器和正余弦查询表构造DDS模块,控制频率控制字M2作为相位累加器的输入端,控制DDS输出信号的频率,相位累加器的输出作为正余弦查找表的输入,其中余弦查找表和正弦查找表中存入的分别是一个完整周期的余弦信号和正弦信号,余弦查询表和正弦查询表的输出分别为I路与Q路的上变频载波信号,设为Iduc和Qduc,各模块时钟均为系统时钟;
(55)通过有符号数乘法器,将经过各项调制后的I路信号与I路载波信号Iduc相乘,得到SIf,经过各项调制后的Q路信号与Q路载波信号Qduc相乘,得到SQf,通过有符号数减法器用SIf减去SQf得SU,SU即为上变频的输出信号。
CN202010100663.6A 2020-02-18 2020-02-18 一种lfmcw雷达目标回波信号模拟方法 Active CN111337890B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010100663.6A CN111337890B (zh) 2020-02-18 2020-02-18 一种lfmcw雷达目标回波信号模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010100663.6A CN111337890B (zh) 2020-02-18 2020-02-18 一种lfmcw雷达目标回波信号模拟方法

Publications (2)

Publication Number Publication Date
CN111337890A true CN111337890A (zh) 2020-06-26
CN111337890B CN111337890B (zh) 2023-05-09

Family

ID=71181688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010100663.6A Active CN111337890B (zh) 2020-02-18 2020-02-18 一种lfmcw雷达目标回波信号模拟方法

Country Status (1)

Country Link
CN (1) CN111337890B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112187278A (zh) * 2020-10-15 2021-01-05 西北工业大学 一种声纳信号模拟器延时控制的声纳信号模拟方法及装置
CN112904285A (zh) * 2021-02-03 2021-06-04 北京航空航天大学杭州创新研究院 一种基于fpga芯片的信号采集方法及装置
CN115166640A (zh) * 2022-06-08 2022-10-11 中国计量大学 一种时域Frank码信号产生方法
CN115425967A (zh) * 2022-08-31 2022-12-02 北京北方华创微电子装备有限公司 相位同步装置和方法、射频电源、半导体工艺设备
CN116884455A (zh) * 2023-07-12 2023-10-13 武汉贞坤电子有限公司 用于多目标模拟的双口ram读出方法及电路

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101082667A (zh) * 2006-06-01 2007-12-05 北京航空航天大学 一种毫米波捷变频雷达目标模拟器
CN102023293A (zh) * 2010-09-14 2011-04-20 中国兵器工业第二○六研究所 多波束比幅测角雷达目标中频回波仿真系统及其控制方法
RU2014109008A (ru) * 2014-03-07 2015-09-20 Открытое Акционерное Общество "Уральское проектно-конструкторское бюро "Деталь" Имитатор радиолокационной цели при зондировании преимущественно длительными сигналами
CN104965197A (zh) * 2015-06-30 2015-10-07 南京理工大学 基于fpga的雷达回波信号模拟器
RU2625567C1 (ru) * 2016-07-11 2017-07-17 Акционерное общество "Уральское проектно-конструкторское бюро "Деталь" Устройство для имитации ложной радиолокационной цели при зондировании сигналами с линейной частотной модуляцией
CN108089164A (zh) * 2017-10-30 2018-05-29 北京航天福道高技术股份有限公司 一种宽带目标模拟器信号处理单元
CN108828540A (zh) * 2018-07-25 2018-11-16 华清瑞达(天津)科技有限公司 雷达回波模拟装置
CN109765535A (zh) * 2019-03-05 2019-05-17 上海志良电子科技有限公司 超高速目标雷达回波的模拟方法及模拟器
CN109884603A (zh) * 2019-03-14 2019-06-14 中国人民解放军海军航空大学青岛校区 基于pxi机箱的雷达目标模拟系统及方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101082667A (zh) * 2006-06-01 2007-12-05 北京航空航天大学 一种毫米波捷变频雷达目标模拟器
CN102023293A (zh) * 2010-09-14 2011-04-20 中国兵器工业第二○六研究所 多波束比幅测角雷达目标中频回波仿真系统及其控制方法
RU2014109008A (ru) * 2014-03-07 2015-09-20 Открытое Акционерное Общество "Уральское проектно-конструкторское бюро "Деталь" Имитатор радиолокационной цели при зондировании преимущественно длительными сигналами
CN104965197A (zh) * 2015-06-30 2015-10-07 南京理工大学 基于fpga的雷达回波信号模拟器
RU2625567C1 (ru) * 2016-07-11 2017-07-17 Акционерное общество "Уральское проектно-конструкторское бюро "Деталь" Устройство для имитации ложной радиолокационной цели при зондировании сигналами с линейной частотной модуляцией
CN108089164A (zh) * 2017-10-30 2018-05-29 北京航天福道高技术股份有限公司 一种宽带目标模拟器信号处理单元
CN108828540A (zh) * 2018-07-25 2018-11-16 华清瑞达(天津)科技有限公司 雷达回波模拟装置
CN109765535A (zh) * 2019-03-05 2019-05-17 上海志良电子科技有限公司 超高速目标雷达回波的模拟方法及模拟器
CN109884603A (zh) * 2019-03-14 2019-06-14 中国人民解放军海军航空大学青岛校区 基于pxi机箱的雷达目标模拟系统及方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
任丽香 等: "一种通用雷达信号模拟器的设计与实现" *
卢丁丁 等: "一种雷达中频信号回波模拟系统的设计" *
崔佩佩 等: "基于LabVIEW FPGA模块的雷达模拟回波多普勒频移调制" *
王建超 等: "一种宽带LFMCW雷达信号模拟技术研究" *
靳澎 等: "基于PCI总线的雷达目标回波模拟" *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112187278A (zh) * 2020-10-15 2021-01-05 西北工业大学 一种声纳信号模拟器延时控制的声纳信号模拟方法及装置
CN112904285A (zh) * 2021-02-03 2021-06-04 北京航空航天大学杭州创新研究院 一种基于fpga芯片的信号采集方法及装置
CN112904285B (zh) * 2021-02-03 2023-12-15 北京航空航天大学杭州创新研究院 一种基于fpga芯片的信号采集方法及装置
CN115166640A (zh) * 2022-06-08 2022-10-11 中国计量大学 一种时域Frank码信号产生方法
CN115166640B (zh) * 2022-06-08 2024-05-07 中国计量大学 一种时域Frank码信号产生方法
CN115425967A (zh) * 2022-08-31 2022-12-02 北京北方华创微电子装备有限公司 相位同步装置和方法、射频电源、半导体工艺设备
CN115425967B (zh) * 2022-08-31 2023-11-14 北京北方华创微电子装备有限公司 相位同步装置和方法、射频电源、半导体工艺设备
CN116884455A (zh) * 2023-07-12 2023-10-13 武汉贞坤电子有限公司 用于多目标模拟的双口ram读出方法及电路

Also Published As

Publication number Publication date
CN111337890B (zh) 2023-05-09

Similar Documents

Publication Publication Date Title
CN111337890B (zh) 一种lfmcw雷达目标回波信号模拟方法
US7719446B2 (en) Method and apparatus for computing interpolation factors in sample rate conversion systems
US8022863B1 (en) Spatially variant apodization (SVA) filter
US9124244B2 (en) Method and device for generating a filter coefficient in real time
CN103019117B (zh) 基于PXIe接口的数字化仪
CN109765535B (zh) 超高速目标雷达回波的模拟方法及模拟器
CN106502309B (zh) 基于da归零保持函数的时域交织任意波形合成装置及方法
US4011438A (en) Simplified digital moving target indicator filter
CN108427102A (zh) 一种雷达基带杂波生成装置及方法
CN201226025Y (zh) 一种脉冲多普勒雷达信号处理器
WO2007084687A1 (en) Asynchronous sample rate conversion using a digital simulation of an analog filter
CN104267385A (zh) 具有预失真功能的线性调频信号发生器
JP2762080B2 (ja) 平方根の概算装置
CN107994954A (zh) 一种模拟多普勒变化的正交调制信号生成装置
CN109976660B (zh) 基于线性插值的任意信号采样率重构方法和数据采样系统
CN103259602B (zh) 信号产生方法及信号产生系统
CN114185014B (zh) 一种应用于雷达信号处理的并行卷积方法及装置
CN116466308A (zh) 回波动态调整方法、装置、设备、介质及雷达回波模拟器
CN113687320B (zh) 一种雷达回波的模拟方法、装置、终端设备和存储介质
CN115685108A (zh) 脉冲伪码体制引信体目标模拟系统及其方法
KR102064232B1 (ko) 데이터 레이트 변환 장치 및 방법
Bokov et al. Generation of Radar Ground Clutter Echoes with Jakes' Doppler Spectrum on FPGA
Shehata et al. Design and implementation of lfmcw radar signal processor for slowly moving target detection using fpga
CN110061719A (zh) 一种基于zynq-7010的自适应滤波器
Mehra et al. Reconfigurable Area and Speed Efficient Interpolator Using DALUT Algorithm

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