CN107479035B - 一种基于fpga和dac的宽带线性调频信号实时产生方法 - Google Patents

一种基于fpga和dac的宽带线性调频信号实时产生方法 Download PDF

Info

Publication number
CN107479035B
CN107479035B CN201710572132.5A CN201710572132A CN107479035B CN 107479035 B CN107479035 B CN 107479035B CN 201710572132 A CN201710572132 A CN 201710572132A CN 107479035 B CN107479035 B CN 107479035B
Authority
CN
China
Prior art keywords
fpga
frequency modulation
linear frequency
dac
real
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
CN201710572132.5A
Other languages
English (en)
Other versions
CN107479035A (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.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
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 Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201710572132.5A priority Critical patent/CN107479035B/zh
Publication of CN107479035A publication Critical patent/CN107479035A/zh
Application granted granted Critical
Publication of CN107479035B publication Critical patent/CN107479035B/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/28Details of pulse systems
    • G01S7/282Transmitters

Abstract

本发明提出了一种基于FPGA和DAC的宽带线性调频信号实时产生方法,首先确定需要产生的宽带线性调频信号参数,然后确定DAC芯片的处理时钟速率,根据DAC芯片的处理时钟速率以及FPGA芯片的处理能力确定FPGA实时并行处理线性调频信号的路数以及时钟速率,并确定FPGA实时计算的角度量化位数。然后通过FPGA加法器计算每路线性调频信号的实时角度,求取每路线性调频信号实时角度的余弦值,最后通过FPGA控制DAC芯片实时输出宽带线性调频信号。应用本发明方法,信号的带宽和信号频率上限不受器件采样频率的限制,同时突破了传统方法FPGA内部存储器资源受限的瓶颈,不需要外部存储器,降低了FPGA程序设计的时序风险,提高了设计的可靠性。

Description

一种基于FPGA和DAC的宽带线性调频信号实时产生方法
技术领域
本发明涉及一种基于FPGA和DAC的宽带线性调频信号实时产生方法,属于雷达宽带线性调频信号产生领域。
背景技术
宽带线性调频信号广泛应用于遥感卫星信号处理类单机中,在对地探测雷达、星载SAR、散射计等遥感载荷中,宽带线性调频信号的性能是一项非常重要的指标。宽带线性调频信号比窄带线性调频信号有显著的优势。宽带线性调频信号增加了系统的信息容量,提供更多的关于目标频域信息,为目标识别提供更多的依据,降低多径效应影响,同时,宽带线性调频信号能增强雷达的抗干扰能力,有效地从杂波中检测出目标。雷达探测距离远近与信号的功率相关,探测距离越远需要的信号功率越大。受发射机和馈电系统功率的限制,宽带线性调频信号是增大雷达信号的功率一个很好的选择。速度分辨率取决于对目标的总的观测时间,角度分辨率取决于天线孔径和对目标的观测时间。宽带线性调频信号也是提高速度分辨率的一个有效方法。
采用DDS技术产生宽带线性调频信号己成为高分辨雷达系统中较为普遍的方法。DDS信号产生技术具有众多优点:极高的频率稳定度、超高速的频率转换时间(可达ns)、相位噪声低、信号形式灵活多样、集成度高。同时,它也存在着信号的带宽和信号频率上限受到器件采样频率的限制(常用DDS芯片AD9858,最高输出的频率为400MHz)。
随着高速DAC芯片的广泛引用,采用FPGA控制高速DAC芯片产生宽带线性调频信号也得到应用,常规方法是从FPGA内部或外部的存储器中读取宽带线性调频信号的波形文件实时输出,控制DAC产生宽带线性调频信号。该方法存在内部存储器资源受限、外部存储器硬件设计复杂、FPGA时序风险较大等问题。
发明内容
本发明的技术解决问题:克服现有技术的不足,提供一种基于FPGA和DAC的宽带线性调频信号实时产生方法,一方面信号的带宽和信号频率上限不受器件采样频率的限制,另一方面突破了传统方法FPGA内部存储器资源受限的瓶颈,不需要外部存储器,降低了FPGA程序设计的时序风险,提高了设计的可靠性。
本发明的技术解决方案是:一种基于FPGA和DAC的宽带线性调频信号实时产生方法,步骤如下:
(1)确定需要的宽带线性调频信号参数,所述宽带线性调频信号参数包括中心频率、时宽、带宽、占空比;
(2)结合宽带线性调频信号参数确定DAC芯片的处理时钟速率;
(3)根据DAC芯片的处理时钟速率以及FPGA芯片的处理能力确定FPGA实时并行处理宽带线性调频信号的路数以及时钟速率;根据宽带线性调频信号的精度要求确定FPGA实时计算的角度量化位数;
(4)通过FPGA加法器计算每路宽带线性调频信号的实时角度;
(5)求取每路宽带线性调频信号实时角度的余弦值;
(6)确定DAC芯片的MUX配比值,通过FPGA控制DAC芯片实时输出宽带线性调频信号。
所述步骤(2)中,DAC芯片的处理时钟速率的确定方法如下:
DAC芯片的处理时钟速率fDAC满足
Figure BDA0001349954240000021
n取能够满足fDAC≥2(fH-fL)的最大整数值,f0表示宽带线性调频信号的中心频率,fH表示宽带线性调频信号的最高频率,fL表示宽带线性调频信号的最低频率,fDMAX表示DAC芯片标称的最高工作频率。
所述步骤(3)中,FPGA实时并行处理宽带线性调频信号的路数以及时钟速率的确定方法如下:
FPGA实时并行处理宽带线性调频信号的时钟速率
Figure BDA0001349954240000031
FPGA与DAC的接口数据速率
Figure BDA0001349954240000032
其中MUX为DAC芯片的MUX值,MUX的选取需同时满足
Figure BDA0001349954240000033
以及
Figure BDA0001349954240000034
其中fmax和fmax_int分别表示FPGA芯片的最高工作速率和最高接口速率,k1表示FPGA芯片最高工作速率的降额因子,k2表示FPGA芯片最高接口速率的降额因子;
FPGA实时并行处理宽带线性调频信号的路数N=2×MUX。
所述步骤(4)中,计算FPGA处理的每路宽带线性调频信号实时角度的方法如下:
(4.1)利用如下标准线性调频信号表达式计算FPGA处理的第i路宽带线性调频信号实时角度的初值Θ(i-1):
Θ(i-1)=θ((i-1)T)×2p/2π1≤i≤N
θ(t)=2πfat+παt2,t=[0,τ)
其中fa表示宽带线性调频信号的起始频率,
Figure BDA0001349954240000035
表示宽带线性调频信号的调频率,B表示宽带线性调频信号带宽,τ表示宽带线性调频信号时宽,p表示FPGA实时计算的角度量化位数;T为DAC时钟周期;
(4.2)第m时刻,FPGA处理的第i路宽带线性调频信号实时角度Θ(Nm+i-1)为:
Θ(Nm+i-1)=NfaT2p+(N+2i-2)NαT22p-1+N2(m-1)αT22p+Θ(N(m-1)+i-1)
N为FPGA实时并行处理宽带线性调频信号的路数。
所述步骤(5)中求取FPGA处理的每路宽带线性调频信号实时角度余弦值的方法如下:
(5.1)在FPGA内部建立N个独立的“Sine/Cosine Look-Up Table”核,N为FPGA实时并行处理线性调频信号的路数,N个独立的“Sine/Cosine Look-Up Table”核与FPGA实时并行处理的N路线性调频信号一一对应;
(5.2)确定“Sine/Cosine Look-Up Table”核要求输入角度的量化位宽q,其中q≤16;
(5.3)为每个“Sine/Cosine Look-Up Table”核输入对应路的宽带线性调频信号实时角度的高q位,每个“Sine/Cosine Look-Up Table”核即输出对应路线性调频信号实时角度的余弦值。
所述步骤(5.3)中每个“Sine/Cosine Look-Up Table”核输出余弦值的位宽与DAC芯片的数字量化位数一致。
所述DAC芯片的MUX配比值为N:2,N为FPGA实时并行处理线性调频信号的路数。
所述DAC芯片后端设置有带通滤波器,用于保留目标奈奎斯特域的波形,滤除非目标奈奎斯特域的波形。
本发明与现有技术相比的有益效果是:
(1)相对于传统的使用DDS芯片产生宽带线性调频信号,本发明提出的方法可以控制实时角度计算的量化位数,选取量化位宽大于32bit,即可保证实时角度的量化误差小于传统的DDS芯片,提高了宽带线性调频信号的调频精度;另一方面DAC芯片的控制采用带通采样定理,降低了DAC芯片的处理时钟速率,同时在多个奈奎斯特域均可产生相应的宽带线性调频信号,突破了DDS芯片产生线性调频信号频率限制的瓶颈。以AD9858为例,最高输出的频率为400MHz,本发明提出的方法能够输出的最高频率由DAC芯片的输出带宽决定,最高输出频率可达6GHz。
(2)相对于传统的从FPGA内部或外部存储器中读取波形数据,控制DAC芯片实时产生线性调频信号的方法,本发明提出的方法仅需要确定数字量化后的角度初值以及递增常量,就可以实时产生宽带线性调频信号,提高了线性调频信号产生的自适应性;同时本发明仅通过FPGA芯片内部的加法器资源就可以实时地计算出宽带线性调频信号的角度信息,不使用FPGA内部的存储器资源,突破了传统方法FPGA内部存储器资源受限的瓶颈,不使用外部存储器资源,简化了硬件设计。
(3)在DAC控制方面,本发明提出的方法采用多路数据并行处理的方式,可以降低FPGA的处理速率,从而降低了FPGA程序设计的时序风险,提高了设计的可靠性。
附图说明
图1是本发明方法实时处理流程图;
图2是本发明实施例中的硬件实现框图;
图3是本发明实施例中的FPGA软件接口示意图;
图4是本发明实施例中宽带线性调频信号的频谱特性;
图5是本发明实施例中宽带线性调频信号的频率误差曲线。
具体实施方式
本发明提出了一种基于FPGA和DAC的宽带线性调频信号实时产生方法,利用线性调频信号频率变化的线性关系,仅通过FPGA芯片内部的加法器资源就可以实时地计算出线性调频信号的角度信息,大幅降低FPGA资源使用率的同时能够根据任务要求实时切换线性调频信号的控制参数;在DAC控制方面,采用多路数据并行处理的方式,可以降低FPGA的处理时钟速率,提高了设计的可靠性。
如图1所示,本发明方法步骤如下:
(1)确定需要产生的宽带线性调频信号参数,宽带线性调频信号参数包括中心频率、时宽、带宽、占空比。
(2)结合宽带线性调频信号参数和数字信号的采样定理确定DAC芯片的处理时钟速率。
DAC芯片的处理时钟速率的确定方法如下:
DAC芯片能够无失真地合成相应的宽带线性调频信号应该满足带通采样定理,即DAC芯片的处理时钟速率应满足
Figure BDA0001349954240000051
n取能够满足fDAC≥2(fH-fL)的最大整数值。fDAC为DAC芯片的处理时钟速率,标识f0表示宽带线性调频信号的中心频率,fH表示宽带线性调频信号的最高频率,fL表示宽带线性调频信号的最低频率,fDMAX表示DAC芯片标称的最高工作频率。目前主流的高速DAC芯片能够保证在第一、第二和第三奈奎斯特区域内合成高质量宽带信号,所以n的取值建议选取范围为0、1、2。
宽带线性调频信号的最高频率=宽带线性调频信号中心频率+1/2*带宽。
宽带线性调频信号的最低频率=宽带线性调频信号中心频率-1/2*带宽。
(3)根据DAC芯片的处理时钟速率以及FPGA芯片的处理能力确定FPGA实时并行处理线性调频信号的路数以及时钟速率;根据宽带线性调频信号的精度要求确定FPGA实时计算的角度量化位数。
FPGA实时并行处理宽带线性调频信号的路数以及时钟速率的确定方法如下:
FPGA实时并行处理宽带线性调频信号的速率
Figure BDA0001349954240000061
FPGA与DAC的接口数据速率
Figure BDA0001349954240000062
其中MUX为DAC芯片的MUX值,常用的高速DAC芯片的MUX通常为2:1或4:1。
DAC芯片的MUX值是可以选择的,MUX值的选取保证f和fint满足FPGA芯片对应速率的降额要求,即应该同时满足
Figure BDA0001349954240000063
Figure BDA0001349954240000064
其中fmax、fmax_int分别表示FPGA芯片的最高工作速率和最高接口速率,k1、k2分别表示两种速率的降额因子,降额因子建议小于0.8。
FPGA实时并行处理线性调频信号的路数N=2×MUX。
常用DDS芯片AD9858实时角度的量化位数为32位,采用FPGA实时计算角度信息,建议实时角度的量化位数p>32位,能够确保该方法的实时角度计算性能优于DDS芯片。
(4)通过FPGA加法器计算每路线性调频信号的实时角度。
具体方法如下:
标准的线性调频信号表达式为:
θ(t)=2πfat+παt2,t=[0,τ)
其中fa表示宽带线性调频信号的起始频率,
Figure BDA0001349954240000071
表示线性调频信号的调频率,B表示信号带宽,τ表示时宽。对角度值采用p位进行量化即:Θ=(θ/2π)×2p,可将θ量化至0~2p之间。FPGA实时并行处理N路数据,第m时刻(即第m+1个DAC时钟周期),各路数据的实时角度量化值为:
第1路:Θ(Nm)=fa×(NmT)×2p+α×(NmT)2×2p-1
第2路:Θ(Nm+1)=fa×(Nm+1)×T×2p+α×((Nm+1)T)2×2p-1
第i路:Θ(Nm+i-1)=fa×(Nm+i-1)×T×2p+α×((Nm+i-1)T)2×2p-1
FPGA实时并行处理N路数据为例,其他多路实时并行处理的推导用相同方法即可得出。以DAC时钟周期T的N倍为间隔进行离散化采样,分别计算各路角度的差值:
第1路:
ΔΘ1(m)=Θ(Nm)-Θ(Nm-N)=NfaT2p+N2αT22p-1+N2(m-1)αT22p
ΔΘ1(m+1)=Θ(Nm+N)-Θ(Nm)=NfaT2p+N2αT22p-1+N2mαT22p
两个角度变化量的差值:
ΔΘ1(m+1)-ΔΘ1(m)=N2αT22p
第2路:
ΔΘ2(m)=Θ(Nm+1)-Θ(Nm-N+1)=NfaT2p+(N+2)NαT22p-1+N2(m-1)αT22p
ΔΘ2(m+1)=Θ(Nm+N+1)-Θ(Nm+1)=NfaT2p+(N+2)NαT22p-1+N2mαT22p
相邻两个角度变化量的差值:
ΔΘ2(m+1)-ΔΘ2(m)=N2αT22p
第i路:
ΔΘi(m)=Θ(Nm+i-1)-Θ(Nm-N+i-1)=NfaT2p+(N+2i-2)NαT22p-1+N2(m-1)αT22p
ΔΘi(m+1)=Θ(Nm+N+i-1)-Θ(Nm+i-1)=NfaT2p+(N+2i-2)NαT22p-1+N2mαT22p
相邻两个角度变化量的差值:
ΔΘi(m+1)-ΔΘi(m)=N2αT22p
可以看出,相邻两个角度变化量的差值均N2αT22p,所以通过两级加法器就可以完成线性调频信号实时角度值的求取。
在实际的FPGA程序设计开始时,需要根据线性调频信号的参数以及DAC芯片的时钟周期计算出N路角度的初值:Θ(i-1)(i=1,2...N),以及上述公式中的其他常量N2αT22p、NfaT2p、(N+2i-2)NαT22p-1(i=1,2...N)等,然后求取实时角度数据,具体方法如下。
第i(i=1,2...N)路的实时数据为:
第1个DAC时钟周期:Θ(i-1):初值,根据标准的线性调频信号表达式获得,即Θ(i-1)=θ((i-1)T)×2p/2π1≤i≤N
第2个DAC时钟周期:Θ(N+i-1)=NfaT2p+(N+2i-2)NαT22p-1+Θ(i-1)
第3个DAC时钟周期:
Θ(2N+i-1)=NfaT2p+(N+2i-2)NαT22p-1+N2αT22p+Θ(N+i-1)
第m+1个DAC时钟周期:
Θ(Nm+i-1)=NfaT2p+(N+2i-2)NαT22p-1+N2(m-1)αT22p+Θ(N(m-1)+i-1)
(5)求取每路线性调频信号实时角度的余弦值。
方法如下:
(5.1)在FPGA内部建立N个独立的“Sine/Cosine Look-Up Table”核,N为FPGA实时并行处理线性调频信号的路数,N个独立的“Sine/Cosine Look-Up Table”核与FPGA实时并行处理的N路线性调频信号一一对应;
(5.2)确定“Sine/Cosine Look-Up Table”核要求输入角度的量化位宽q,其中q≤16;
(5.3)为每个“Sine/Cosine Look-Up Table”核输入对应路的线性调频信号实时角度的高q位,每个“Sine/Cosine Look-Up Table”核即输出对应路线性调频信号实时角度的余弦值,实时角度的余弦值y=cos(2πfat+παt2)。
每个“Sine/Cosine Look-Up Table”核输出的实时数据将作为DAC芯片的输入数字信号,故该核的输出位宽应该与DAC芯片的数字量化位数保持一致。
(6)确定DAC芯片的MUX配比值,通过FPGA控制DAC芯片实时输出宽带线性调频信号。其中DAC芯片的MUX配比值为N:2,N为FPGA实时并行处理线性调频信号的路数。
DAC芯片后端设置有带通滤波器,用于保留目标奈奎斯特域的波形,滤除非目标奈奎斯特域的波形。
实施例:
设需要产生的宽带线性调频信号中心频率f0=600MHz,带宽B=200MHz,时宽τ=50us。
宽带线性调频信号产生步骤如下:
(a1)确定宽带线性调频信号中心频率、时宽、带宽、占空比等参数。设计实现宽带线性调频信号的中心频率f0=600MHz,带宽B=200MHz,时宽τ=50us,占空比50%,信号最低频率fL=500MHz,最高频率fH=700MHz。
(a2)结合宽带线性调频信号的参数和数字信号的采样定理确定DAC芯片的处理时钟速率。
根据数字带通采样定理DAC工作频率
Figure BDA0001349954240000101
n取1,f0=600MHz,计算得出fDAC为800MHz,且满足fDAC≥2(fH-fL),最终使用DAC芯片输出的第2奈奎斯特域的波形。
(a3)根据DAC芯片的处理时钟速率以及FPGA芯片的处理能力确定FPGA实时处理数据的路数以及时钟速率;根据宽带线性调频信号的精度要求确定FPGA实时计算的角度量化位数。FPGA选用了Xilinx公司的XC4VSX55,DAC芯片选用了e2v公司的EV10DS130AG。
DAC芯片选取的工作速率为800MHz,通过4:1MUX配比可以将FPGA的接口数据速率fint降低为200MHz。FPGA接口采用ODDR设计,FPGA内部实时并行处理宽带线性调频信号的速率f可降低为100MHz。FPGA接口速率和处理速率均满足FPGA芯片的速率降额要求。实时角度计算的量化位数选取48bit,优于DDS芯片的32bit。FPGA实时并行处理线性调频信号的路数N=8。
FPGA输出的数字线性调频信号通过DAC芯片后会在三个奈奎斯特域产生带宽为200M的波形,因此DAC芯片后端需要增加带通滤波器,滤除第一和第三奈奎斯特域的波形,仅保留第二奈奎斯特域,中心频率为600MHz的波形。硬件设计框图如图2所示。
(a4)确定FPGA实时处理的初始角度及常量,通过FPGA的加法器求取实时角度数据。标准的线性调频信号表达式为:
θ(t)=2πfat+παt2,t=[0,τ)
其中fa表示宽带线性调频信号的起始频率,
Figure BDA0001349954240000102
表示线性调频信号的调频率,对角度值采用p位进行量化即:Θ=(θ/2π)×2p,可将θ量化至0~2p之间
FPGA实时并行处理8路数据,以DAC时钟周期的T的8倍为间隔进行离散化采样,在实际的FPGA程序设计开始时,需要根据线性调频信号的参数以及DAC芯片的时钟周期计算出8路角度的初值:Θ(i-1)(i=1,2...8),以及上述公式中的其他常量64αT22p、8faT2p、8(6+2i)αT22p-1(i=1,2...8)等,然后求取实时角度数据,具体方法如下。
第i(i=1,2...8)路的实时数据为:
第1个DAC时钟周期:Θ(i-1):初值
第2个DAC时钟周期:Θ(8+i-1)=8faT2p+8(6+2i)αT22p-1+Θ(i-1)
第3个DAC时钟周期:Θ(16+i-1)=8faT2p+8(6+2i)αT22p-1+64αT22p+Θ(7+i)
第m+1个DAC时钟周期:
Θ(8m+i-1)=8faT2p+8(6+2i)αT22p-1+64(m-1)αT22p+Θ(8m+i-9)
(a5)使用FPGA内部的“Sine/Cosine Look-Up Table”核,选取合适的量化位宽以及输出数据位宽,求取实时角度的余弦值。实施例采用8路数据并行处理,所以共需要建立8个独立的“Sine/Cosine Look-Up Table”,采用16位角度量化位宽,即“Sine/Cosine Look-Up Table”核的输入数据即为前述实时角度数据Θ的高16位。输出数据位宽为10bit,与DAC芯片EV10DS130AG的输入位宽一致。
(a6)结合DAC芯片的器件手册,选取合理的MUX配比值,控制DAC芯片实时输出宽带线性调频信号。实施例中DAC芯片采用4:1MUX配比,8路角度数据分别通过相对应的“Sine/Cosine Look-Up Table”核,计算出实时地8路余弦数字量,然后第1路与第5路、第2路与第6路、第3路与第7路、第4路与第8路通过四个ODDR核,将8路100MHz的数据合并为4路200M的数据,之后再经过单端转差分处理后最终输出给DAC芯片。FPGA软件接口示意图如图3所示。
实施例中DAC芯片EV10DS130AG有效位数8位,根据当前设计采用MATLAB进行了仿真分析,设计实例结果与仿真分析基本一致。图4给出了基于FPGA和DAC的宽带线性调频信号的频谱特性,图4中(a)为理论波形的频谱,图4中(b)为实施例的仿真波形频谱,可以看出实施例仿真波形的带宽与理想波形的带宽完全一致,实施例仿真波形的带外抑制方面不及理想波形,但能够达到50dBc,满足使用要求;图5给出了宽带线性调频信号的频率误差曲线,经计算,宽带线性调频信号的线性度<0.3%,宽带线性调频信号的频率一致性较好。实际应用中,增加实时角度计算的量化位数、选取有效位数更高的DAC芯片都能够提高宽带线性调频信号的带外抑制,减小频率误差。
应用本发明方法,宽带线性调频信号的时宽、带宽、调频斜率、调制模式均可实时切换和控制。相对于传统的使用FPGA芯片控制DDS芯片(常用芯片:AD公司的AD9858)产生线性调频信号,该方法一方面提高了线性调频信号的调频精度,另一方面突破了DDS芯片产生线性调频信号频率限制的瓶颈;相对于传统的从FPGA内部或外部存储器中读取波形数据,控制DAC芯片实时产生线性调频信号的方法,该方法一方面提高了线性调频信号产生的自适应性,另一方面突破了传统方法FPGA内部存储器资源受限的瓶颈。使用FPGA芯片实时地计算线性调频信号的角度信息,能够大幅度的降低FPGA的资源使用率,同时能够根据任务要求实时切换线性调频信号的控制参数。在DAC控制方面,采用多路数据并行处理的方式,大幅降低了FPGA的处理速率,从而降低了FPGA程序设计的时序风险,提高了设计的可靠性。本发明可以用于多种参数化的宽带线性调频信号的实时产生,具有广阔的发展应用前景。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。

Claims (8)

1.一种基于FPGA和DAC的宽带线性调频信号实时产生方法,其特征在于步骤如下:
(1)确定需要的宽带线性调频信号参数,所述宽带线性调频信号参数包括中心频率、时宽、带宽、占空比;
(2)结合宽带线性调频信号参数确定DAC芯片的处理时钟速率;
(3)根据DAC芯片的处理时钟速率以及FPGA芯片的处理能力确定FPGA实时并行处理宽带线性调频信号的路数以及时钟速率;根据宽带线性调频信号的精度要求确定FPGA实时计算的角度量化位数;
(4)通过FPGA加法器计算每路宽带线性调频信号的实时角度;
(5)求取每路宽带线性调频信号实时角度的余弦值;
(6)确定DAC芯片的MUX配比值,通过FPGA控制DAC芯片实时输出宽带线性调频信号。
2.根据权利要求1所述的一种基于FPGA和DAC的宽带线性调频信号实时产生方法,其特征在于:所述步骤(2)中,DAC芯片的处理时钟速率的确定方法如下:
DAC芯片的处理时钟速率fDAC满足
Figure FDA0002160775660000011
n取能够满足fDAC≥2(fH-fL)的最大整数值,f0表示宽带线性调频信号的中心频率,fH表示宽带线性调频信号的最高频率,fL表示宽带线性调频信号的最低频率,fDMAX表示DAC芯片标称的最高工作频率。
3.根据权利要求2所述的一种基于FPGA和DAC的宽带线性调频信号实时产生方法,其特征在于:所述步骤(3)中,FPGA实时并行处理宽带线性调频信号的路数以及时钟速率的确定方法如下:
FPGA实时并行处理宽带线性调频信号的时钟速率
Figure FDA0002160775660000012
FPGA与DAC的接口数据速率
Figure FDA0002160775660000021
其中MUX为DAC芯片的MUX值,MUX的选取需同时满足
Figure FDA0002160775660000022
以及
Figure FDA0002160775660000023
其中fmax和fmax_int分别表示FPGA芯片的最高工作速率和最高接口速率,k1表示FPGA芯片最高工作速率的降额因子,k2表示FPGA芯片最高接口速率的降额因子;
FPGA实时并行处理宽带线性调频信号的路数N=2×MUX。
4.根据权利要求1所述的一种基于FPGA和DAC的宽带线性调频信号实时产生方法,其特征在于:所述步骤(4)中,计算FPGA处理的每路宽带线性调频信号实时角度的方法如下:
(4.1)利用如下标准线性调频信号表达式计算FPGA处理的第i路宽带线性调频信号实时角度的初值Θ(i-1):
Θ(i-1)=θ((i-1)T)×2p/2π 1≤i≤N
θ(t)=2πfat+παt2,t∈[0,τ)
其中fa表示宽带线性调频信号的起始频率,
Figure FDA0002160775660000024
表示宽带线性调频信号的调频率,B表示宽带线性调频信号带宽,τ表示宽带线性调频信号时宽,p表示FPGA实时计算的角度量化位数;T为DAC时钟周期;
(4.2)第m时刻,FPGA处理的第i路宽带线性调频信号实时角度Θ(Nm+i-1)为:
Θ(Nm+i-1)=NfaT2p+(N+2i-2)NαT22p-1+N2(m-1)αT22p+Θ(N(m-1)+i-1)
N为FPGA实时并行处理宽带线性调频信号的路数。
5.根据权利要求1所述的一种基于FPGA和DAC的宽带线性调频信号实时产生方法,其特征在于:所述步骤(5)中求取FPGA处理的每路宽带线性调频信号实时角度余弦值的方法如下:
(5.1)在FPGA内部建立N个独立的“Sine/Cosine Look-Up Table”核,N为FPGA实时并行处理线性调频信号的路数,N个独立的“Sine/Cosine Look-Up Table”核与FPGA实时并行处理的N路线性调频信号一一对应;
(5.2)确定“Sine/Cosine Look-Up Table”核要求输入角度的量化位宽q,其中q≤16;
(5.3)为每个“Sine/Cosine Look-Up Table”核输入对应路的宽带线性调频信号实时角度的高q位,每个“Sine/Cosine Look-Up Table”核即输出对应路线性调频信号实时角度的余弦值。
6.根据权利要求5所述的一种基于FPGA和DAC的宽带线性调频信号实时产生方法,其特征在于:所述步骤(5.3)中每个“Sine/Cosine Look-Up Table”核输出余弦值的位宽与DAC芯片的数字量化位数一致。
7.根据权利要求1所述的一种基于FPGA和DAC的宽带线性调频信号实时产生方法,其特征在于:所述DAC芯片的MUX配比值为N:2,N为FPGA实时并行处理线性调频信号的路数。
8.根据权利要求1所述的一种基于FPGA和DAC的宽带线性调频信号实时产生方法,其特征在于:所述DAC芯片后端设置有带通滤波器,用于保留目标奈奎斯特域的波形,滤除非目标奈奎斯特域的波形。
CN201710572132.5A 2017-07-13 2017-07-13 一种基于fpga和dac的宽带线性调频信号实时产生方法 Active CN107479035B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710572132.5A CN107479035B (zh) 2017-07-13 2017-07-13 一种基于fpga和dac的宽带线性调频信号实时产生方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710572132.5A CN107479035B (zh) 2017-07-13 2017-07-13 一种基于fpga和dac的宽带线性调频信号实时产生方法

Publications (2)

Publication Number Publication Date
CN107479035A CN107479035A (zh) 2017-12-15
CN107479035B true CN107479035B (zh) 2020-04-10

Family

ID=60595634

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710572132.5A Active CN107479035B (zh) 2017-07-13 2017-07-13 一种基于fpga和dac的宽带线性调频信号实时产生方法

Country Status (1)

Country Link
CN (1) CN107479035B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110488228B (zh) * 2019-07-11 2021-12-07 中国科学院电子学研究所 线性调频信号生成方法、装置及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101576619A (zh) * 2009-05-07 2009-11-11 北京航空航天大学 基于fpga的uwb雷达信号模拟器及uwb雷达信号产生方法
CN105572641A (zh) * 2015-12-10 2016-05-11 武汉滨湖电子有限责任公司 一种在线重构式宽/窄带通用雷达源及其产生方法
CN106501779A (zh) * 2016-09-27 2017-03-15 哈尔滨工程大学 一种基于sca的软件化lpi雷达信号产生方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101576619A (zh) * 2009-05-07 2009-11-11 北京航空航天大学 基于fpga的uwb雷达信号模拟器及uwb雷达信号产生方法
CN105572641A (zh) * 2015-12-10 2016-05-11 武汉滨湖电子有限责任公司 一种在线重构式宽/窄带通用雷达源及其产生方法
CN106501779A (zh) * 2016-09-27 2017-03-15 哈尔滨工程大学 一种基于sca的软件化lpi雷达信号产生方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"宽带线性调频信号产生方法及实现研究";张俊;《中国优秀硕士学位论文全文数据库》;20150730;第3-4章 *

Also Published As

Publication number Publication date
CN107479035A (zh) 2017-12-15

Similar Documents

Publication Publication Date Title
CN108414966B (zh) 一种基于时间调制的宽带线性调频信号测向系统及方法
CN105024745B (zh) 多路宽带接收信号射频采样下的时延调整方法
CN104201991B (zh) 实现捷变频米波雷达的数字下变频系统
CN108614271A (zh) 一种带反馈校正的多通道超声波任意波形信号发生器
CN105404495B (zh) 用于调制宽带转换器的高速伪随机序列发生器及发生方法
CN110488228A (zh) 线性调频信号生成方法、装置及存储介质
CN104272595A (zh) 用于接收器的逐次逼近模-数转换架构布置
CN103929387B (zh) 基于fpga的大动态范围数字信道化接收机及工作方法
CN108627808A (zh) 雷达干扰器超宽带数字信号处理方法
CN108667483A (zh) 一种用于宽带信号的收发装置
CN106849950A (zh) 基于多速率并行采样的射频信号模数转换系统及方法
CN107479035B (zh) 一种基于fpga和dac的宽带线性调频信号实时产生方法
CN109639277A (zh) 一种基于iserdes和并行fir滤波的高速信号预处理方法
CN108107409A (zh) 一种mimo雷达导引头的中频信号源及其信号产生方法
CN204613394U (zh) 基于dds的四通道宽带雷达信号源
CN102063075A (zh) 中频采集卡板载dsp实时数字信号处理系统
CN113671456A (zh) 一种雷达信号处理仿真平台及仿真方法
CN111367196B (zh) W波段宽带可变分数延时方法及系统
CN202634451U (zh) 宽频无线电频谱监测站
Wang et al. Design and FPGA implementation of digital pulse compression for HF chirp radar based on modified orthogonal transformation
CN104242932B (zh) 一种具有非线性补偿功能的宽带信号发生器
CN101231711B (zh) 声表面波型小波神经网络器件
CN216013642U (zh) 一种雷达信号处理仿真平台
CN103439690B (zh) 一种超宽带雷达频率综合器装置
CN111552204A (zh) 一种基于fpga的旋转变压器模拟器

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