一种基于FPGA的实时产生混沌雷达发射信号系统及方法

一种基于fpga的实时产生混沌雷达发射信号系统及方法
技术领域
1.本发明涉及混沌雷达信号处理技术领域，特别涉及一种基于fpga的实时产生混沌雷达发射信号系统及方法。


背景技术：

2.噪声雷达采用随机或伪随机信号作为发射波形，目前已经应用于军事和民用等多个领域(参考文献[1]：专利cn104777461a，一种载频随机跳变的宽带混沌雷达信号产生方法及系统；参考文献[2]：krzysztof kulpa，signal processing in noise waveform radar，isbn：9781608076611)。
[0003]
在噪声雷达中，随机或伪随机发射信号的产生是一项关键技术。相比于普通的热噪声信号，混沌信号更加容易产生和控制。混沌信号是一种由确定性系统产生的伪随机信号。混沌信号对初值非常敏感，初值的细微变化会导致后续信号差异巨大；而且混沌信号具有非周期和不可预测的特点。采用混沌信号作为雷达的发射信号具有低截获概率特性、有效的频谱利用等优势。
[0004]
现场可编程门阵列(field programmable gate array，fpga)具有大容量、高速度、高集成度等一系列优点，在雷达控制、数字信号处理等方面的应用也越来越多。因此，基于fpga实时产生雷达发射所需的混沌信号将带来新的技术变革。


技术实现要素：

[0005]
本发明的目的在于克服现有技术缺陷，提出了一种基于fpga的实时产生混沌雷达发射信号系统及方法。
[0006]
为了实现上述目的，本发明提出了一种基于fpga的混沌雷达发射信号实时产生系统，所述系统包括：基带iq路波形产生模块、dac模块、iq调制器、发射上变频模块和频率源；其中，
[0007]
所述基带iq路波形产生模块，通过fpga实现，用于分别产生混沌雷达的基带i路波形数据和基带q路波形数据，并输入dac模块；
[0008]
所述dac模块，用于对基带i路波形数据和基带q路波形数据分别进行数模转换，得到基带i路模拟波形和基带q路模拟波形，并输入iq调制器；
[0009]
所述iq调制器，用于对基带i路模拟波形和基带q路模拟波形结合中频本振信号进行iq调制，产生中频混沌信号并输入发射上变频模块；
[0010]
所述发射上变频模块，用于对中频混沌信号和射频本振信号进行混频，产生所需要的射频混沌信号；
[0011]
所述频率源，用于提供fpga的工作时钟，用于提供中频本振信号至iq调制器，还用于提供射频本振信号至发射上变频模块。
[0012]
作为上述系统的一种改进，所述基带iq路波形产生模块包括：混沌序列产生器、带通滤波器、两个变频器和两个低通滤波；其中，两个变频器为并联，每个变频器和一个低通
滤波器串联；
[0013]
所述混沌序列产生器，用于产生混沌序列并输入带通滤波器；
[0014]
所述带通滤波器，用于对混沌序列进行带通滤波处理，得到带限混沌序列并同时输入两个变频器；
[0015]
所述两个变频器，用于对带限混沌序列进行数字下变频，得到两路正交的基带混沌信号，并分别输入各自串联的低通滤波器；
[0016]
所述两个低通滤波器，用于对两路正交的基带混沌信号分别进行低通滤波处理，得到混沌雷达的基带i路波形数据和基带q路波形数据。
[0017]
作为上述系统的一种改进，所述混沌序列为：
[0018]
x(n+1)＝k
·
x(n)
·
(1
‑
x(n))
[0019]
其中，x(n)表示logistic混沌信号映射的当前值，x(n+1)表示logistic混沌信号映射的下一个迭代值，k为logistic混沌信号映射的参数，n为logistic混沌信号映射序列的序号，n＝0,1,2,
…
,n
p
‑
1，n
p
为每个发射脉冲对应的logistic混沌序列长度：
[0020]
n
p
＝t
p
·
f
sr
[0021]
其中，t
p
为发射脉冲信号时长，f
sr
为dac模块的采样率。
[0022]
作为上述系统的一种改进，所述带限混沌序列x
b
(n)为：
[0023][0024]
其中，h
b
(k)为带通滤波器系数，k＝0,1,2,
…
,k
‑
1，k为带通滤波器阶数，h
b
(k)根据带通滤波器的通带范围[0.5
‑
b
r
/f
sr
,0.5+b
r
/f
sr
]计算得到，b
r
为雷达发射的混沌波形的带宽。
[0025]
作为上述系统的一种改进，所述两路正交的基带混沌信号为：
[0026][0027]
其中，f0为进行混频的中频本振信号，f0＝f
sr
/4。
[0028]
作为上述系统的一种改进，所述混沌雷达的基带i路波形数据和基带q路波形数据为：
[0029][0030]
其中，滤波器系数h
l
(m)，m＝0,1,2，
…
，m
‑
1，m为低通滤波器阶数，h
l
(m)根据低通滤波器通带带宽为[0,b
r
/f
sr
]计算得到。
[0031]
一种基于fpga的混沌雷达发射信号实时产生方法，基于上述系统实现，所述方法包括：
[0032]
所述基带iq路波形产生模块分别产生混沌雷达的基带i路波形数据和基带q路波形数据，并输入dac模块；
[0033]
所述dac模块对基带i路波形数据和基带q路波形数据分别进行数模转换，得到基
带i路模拟波形和基带q路模拟波形，并输入iq调制器；
[0034]
所述iq调制器对基带i路模拟波形和基带q路模拟波形结合中频本振信号进行iq调制，产生中频混沌信号并输入发射上变频模块；
[0035]
所述发射上变频模块对中频混沌信号和射频本振信号进行混频，产生所需要的射频混沌信号；
[0036]
所述基带iq路波形产生模块的工作时钟、所述中频本振信号和射频本振信号均由频率源提供。
[0037]
与现有技术相比，本发明的优势在于：
[0038]
1、基于fpga实时计算产生混沌雷达发射信号，每个脉冲的波形均为类似噪声的伪随机信号，且每个脉冲的发射波形均不相同，发射波形对敌方而言即为噪声信号，难以截获和破译，可有效提高雷达的低检测概率特性、低截获概率特性以及电子抗干扰能力；
[0039]
2、基于fpga实时计算产生混沌雷达发射信号，发射的波形为一种由确定性系统产生的伪随机信号，对于己方为确定性波形，在己方雷达信号处理机可以使用混沌模型重构发射的伪随机混沌信号，非常方便的对雷达回波信号进行脉冲压缩处理，而不需要而不需要像随机噪声模拟产生的方式那样将每个脉冲的发射信号经过延迟线再采集回来做定标，可以减少一路数据采集；
[0040]
3、该系统使用通用的fpga和dac，以及iq调制器，可移植性高，方便移植到其他雷达系统中。
附图说明
[0041]
图1是本发明实施例1的基于fpga的实时产生混沌雷达发射信号系统的组成框图；
[0042]
图2是本发明实施例1的基带iq路波形产生模块实时产生混沌雷达发射信号的流程图；
[0043]
图3是本发明实施例1的基带iq路波形产生模块生成logistic混沌序列的状态机框图；
[0044]
图4(a)是本发明实施例1的基带iq路波形产生模块生成的混沌雷达基带i路的波形图；
[0045]
图4(b)是本发明实施例1的基带iq路波形产生模块生成的混沌雷达基带q路的波形图；
[0046]
图5(a)是本发明实施例1的dac模块采样的混沌雷达基带i路波形图；
[0047]
图5(b)是本发明实施例1的dac模块采样的混沌雷达基带q路波形图；
[0048]
图6是本发明实施例1的信号处理机重构发射波形对dac模块采样的基带回波波形脉冲压缩的结果图。
具体实施方式
[0049]
本发明的目的在于基于fpga实时产生混沌雷达发射信号，每个发射脉冲的波形均为类似噪声的伪随机混沌信号，且每个脉冲的发射波形均不相同。发射波形对敌方而言即为噪声信号，难以截获和破译，可有效提高雷达的低截获概率特性以及电子抗干扰能力。同时，发射的波形为由确定性的混沌模型产生的伪随机信号，对于己方为确定性波形，在己方
信号处理机可以使用混沌模型重构发射的伪随机混沌信号，非常方便的对雷达回波信号进行脉冲压缩处理，而不需要像随机噪声模拟产生的方式那样将每个脉冲的发射信号经过延迟线再采集回来做定标(参考文献[2]：krzysztof kulpa，signal processing in noise waveform radar，isbn：9781608076611)，可以减少雷达系统的一路数据采集。
[0050]
下面结合附图和实施例对本发明的技术方案进行详细的说明。
[0051]
实施例1
[0052]
本发明的实施例1提出了一种基于fpga的混沌雷达发射信号实时产生系统，所述系统包括：基带iq路波形产生模块、dac模块、iq调制器、发射上变频模块和频率源；其中，
[0053]
基带iq路波形产生模块，通过fpga实现，用于分别产生混沌雷达的基带i路波形数据和基带q路波形数据，并输入dac模块；
[0054]
基带iq路波形产生模块包括：混沌序列产生器、带通滤波器、两个变频器和两个低通滤波；其中，两个变频器为并联，每个变频器和一个低通滤波器串联；
[0055]
所述混沌序列产生器，用于产生混沌序列并输入带通滤波器；
[0056]
所述带通滤波器，用于对混沌序列进行带通滤波处理，得到带限混沌序列并同时输入两个变频器；
[0057]
所述两个变频器，用于对带限混沌序列进行数字下变频，得到两路正交的基带混沌信号，并分别输入各自串联的低通滤波器；
[0058]
所述两个低通滤波器，用于对两路正交的基带混沌信号分别进行低通滤波处理，得到混沌雷达的基带i路波形数据和基带q路波形数据。
[0059]
dac模块，用于对基带i路波形数据和基带q路波形数据分别进行数模转换，得到基带i路模拟波形和基带q路模拟波形，并输入iq调制器；
[0060]
iq调制器，用于对基带i路模拟波形和基带q路模拟波形结合中频本振信号进行iq调制，产生中频混沌信号并输入发射上变频模块；
[0061]
发射上变频模块，用于对中频混沌信号和射频本振信号进行混频，产生所需要的射频混沌信号；
[0062]
频率源，用于提供fpga的工作时钟，用于提供中频本振信号至iq调制器，还用于提供射频本振信号至发射上变频模块。
[0063]
需要说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
[0064]
在本实施例中，雷达系统发射信号的时长t
p
＝10us，带宽b
r
＝100mhz，dac芯片的采样率为f
sr
＝2500msps。基带iq路波形产生模块在fpga内实时产生混沌雷达发射信号的方法流程图如附图2所示。
[0065]
可选的，上述技术方案中基于fpga实时产生混沌雷达发射信号的方法进一步包括：
[0066]
步骤101)在fpga内利用状态机产生混沌序列。
[0067]
在本实施例中，以logistic混沌信号为例，产生的混沌序列表达式为：
[0068]
x(n+1)＝k
·
x(n)
·
(1
‑
x(n))
ꢀꢀ
(1)
[0069]
其中，x(n)表示logistic映射的当前值，x(n+1)表示logistic映射的下一个迭代值，k为logistic映射的参数，n为logistic映射序号。
[0070]
通常，logistic映射的输入和输出均分布在(0,1)之间，在fpga内实现首先要把(0,1)上的小数写成二进制，表达式如下：
[0071][0072]
其中，c
i
为二进制的0或1。
[0073]
取前l位表示，则有：
[0074][0075]
其中，是一个由l位二进制表示的整数，和小数x(n)一一对应。把x(n+1)和x(n)写成上面形式，可以得到：
[0076]
x(n+1)＝k
·
x(n)
·
(2
l
‑1‑
x(n))
·2‑
(l
‑
1)
ꢀꢀ
(4)
[0077]
其中，x(0)＝round(2
l
‑1·
x(0))，round表示对变量取整。x(0)为logistic混沌序列的初始值，可以由外部通过串口或网口等接口注入到fpga中。在本实施例中，设定k＝0.875，可以看作：k＝1
‑
0.125。在fpga中可以通过移位来实现，0.125即为向右移3bit。
[0078]
上式(4)在fpga中以状态机进行实现，每产生一个logistic混沌序列数值需要4个状态，如附图3所示。
[0079]
根据发射脉冲信号时长t
p
和dac的采样率f
sr
进行计算，得到每个发射脉冲信号所需要的logistic混沌序列的长度为：
[0080]
n
p
＝t
p
·
f
sr
ꢀꢀ
(5)
[0081]
即上式n＝0,1,2,
…
,n
p
‑
1。在本实施例中，n
p
＝25000。
[0082]
步骤102)对产生的混沌序列x(n)进行带限滤波处理，进一步包括：
[0083]
步骤102
‑
1)设定带通滤波器通带范围为[0.5
‑
b
r
/f
sr
,0.5+b
r
/f
sr
]，滤波器阶数为k，计算得到滤波器系数h
b
(k)，k＝0,1,2,
…
,k
‑
1。b
r
为雷达发射的混沌波形的带宽。
[0084]
步骤102
‑
2)对伪随机混沌序列x(n)进行带通滤波处理，得到带限混沌序列x
b
(n)如下所示：
[0085][0086]
在本实施例中，带通滤波器的通带为[0.42,0.58]。设定带通滤波器阶数k＝31，计算得到带限滤波器系数h
b
(k)，然后在fpga中使用ip核实现式(6)对伪随机混沌序列x(n)的带通滤波处理。
[0087]
步骤103)对x
b
(n)进行数字下变频处理，得到i/q两路对应的序列，如下式所示。
[0088][0089]
使用工程应用中的参数提出一种简单设计方式。使用中频信号f0＝f
sr
/4，这样选择参数使得c o(sπ2f0n/f
sr
)信号在一个周期内的数值为{0,
‑
1,0,1}，
‑
sin(2πf0n/f
sr
)信号的数值为{1,0,0,
‑
1}。在fpga内部实现，数值{1,0,
‑
1}分别对应着不变、置零、取反操作。
[0090]
步骤104)对i(n)和q(n)进行低通滤波处理，进一步包括：
[0091]
步骤104
‑
1)设定低通滤波器通带带宽为[0,b
r
/f
sr
]，滤波器阶数为m，计算得到滤波器系数h
l
(m)，m＝0,1,2，
…
，m
‑
1。
[0092]
步骤104
‑
2)对i(n)和q(n)进行低通滤波处理，得到序列i
l
(n)和q
l
(n)如下所示：
[0093][0094]
在本实施例中，低通滤波器的通带为[0,0.08]，阶数m＝31，计算得到低通滤波器系数h
l
(m)。在fpga中使用ip核实现式(8)对i(n)和q(n)的低通滤波处理。经过上述步骤，得到的基带发射波形i路和q路及频谱如附图4(a)、(b)所示。
[0095]
步骤105)fpga按照雷达系统的工作时序将基带发射i路和q路波形同时输出给dac芯片。
[0096]
在实际雷达系统中，将dac输出的i路和q路波形使用采样率为250msps的adc进行采样，得到波形如附图5(a)、(b)所示。在信号处理机使用混沌模型重构发射信号，对adc采样的波形进行脉冲压缩处理，得到结果如附图6所示，可以看到脉冲压缩效果很好。
[0097]
需要说明的是，本实施例以logistic映射为例生成伪随机波形，对于不同算法的混沌波形，均可使用本系统，只是对于波形形成的算法有所区别。
[0098]
实施例2
[0099]
基于实施例1的系统，本申请的实施例2提出了一种基于fpga的混沌雷达发射信号实时产生方法，具体步骤如下：
[0100]
fpga产生混沌雷达基带i路和q路波形数据，同时输送给2路dac芯片；
[0101]
2路dac芯片将i和q两路波形数据分别进行数模转换，产生基带i路和q路模拟波形，输送给iq调制器；
[0102]
iq调制器将输入的i路和q路模拟波形结合中频本振信号lo2进行iq调制，产生中频混沌信号，输送给发射上变频模块；
[0103]
发射上变频模块将中频混沌信号和射频本振lo1进行混频，产生所需要的射频混沌信号。
[0104]
频率源提供fpga的工作时钟，iq调制器使用的中频本振lo2，发射上变频使用的射频本振lo1。
[0105]
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。