CN109889231B

CN109889231B - 基于随机解调和有限新息率的脉冲串信号欠采样方法

Info

Publication number: CN109889231B
Application number: CN201910105414.3A
Authority: CN
Inventors: 付宁; 黄国兴; 乔立岩
Original assignee: Harbin Institute of Technology
Current assignee: Harbin Institute of Technology
Priority date: 2019-02-01
Filing date: 2019-02-01
Publication date: 2021-03-02
Anticipated expiration: 2039-02-01
Also published as: CN109889231A

Abstract

本发明提供基于随机解调和有限新息率的脉冲串信号欠采样方法，信号处理技术领域。本发明首先生成被测脉冲串信号和伪随机序列；然后对被测脉冲串信号和伪随机序列进行随机解调：采用乘法器将被测脉冲串信号和伪随机序列进行混频，得到混频后的信号；采用低通滤波器对混频后的信号进行滤波，得到滤波后的信号；采用数据采集卡对滤波后的信号均匀采样得到样本；最后进行参数估计：利用获得的样本估计被测脉冲串信号中的幅值参数和时延参数。本发明解决了现有脉冲串信号欠采样技术需要根据具体的输入信号单独设计采样结构，导致系统的通用性较差的问题。本发明可用于脉冲串信号欠采样。

Description

基于随机解调和有限新息率的脉冲串信号欠采样方法

技术领域

本发明涉及脉冲串信号欠采样方法，属于信号处理技术领域。

背景技术

脉冲串信号广泛存在于雷达、通信和生物医学等领域。对于该类信号，现有的采样设备大多是基于经典的奈奎斯特采样定理，也就是要求采样率大于或等于信号带宽的两倍。然而，随着现代技术的发展，脉冲串信号中脉冲的时宽越来越窄，对应的信号带宽越来越宽，导致奈奎斯特采样率也越来越高，这无疑给基于奈奎斯特采样定理的采样系统设计带来了严峻的挑战。较高的采样频率要求ADC(Analog-to-Digital Converter模数转换器)器件具有较宽的模拟带宽及采样率，同时，高采样率将产生大量的采样数据，这给系统的存储及传输资源也增加了额外的负担。因而，奈奎斯特采样定理逐渐成为脉冲串信号采样系统设计的瓶颈，制约着信号处理技术的发展。

早在2002年，瑞士联邦理工学院的Vetterli研究团队提出了一种利用信号参数化特性的欠采样理论——有限新息率(Finite Rate of Innovation，简称FRI)采样理论。该理论不同于奈奎斯特采样定理，它指出：对于某些可以由有限数目的自由参量完全表示的参数稀疏信号，如脉冲串信号，只要选用合适的采样核对其滤波后，并以高于或等于信号新息率的速率进行均匀采样，就可以利用采集到的样本完全重构出未知参数。FRI采样理论的采样速率仅由信号的新息率，即单位时间内的自由参数数目决定，一般情况下，信号的新息率要远小于信号的带宽，这样就能大大降低系统采样速率。

考虑如下形式的脉冲串信号：

其中，h(t)是波形已知的脉冲，又可称之为基函数；a_l(l＝1,2,...,L)为幅值参数，t_l(l＝1,2,...,L)为时延参数，t_l∈[0,T)。显然，信号x(t)中的未知参数为2L个幅值和时延参数

对公式(1)所示的脉冲串信号进行FRI采样，一般的手段是通过获取其傅里叶系数来实现的。求该信号x(t)的连续时间傅里叶变换可得：

其中，H(ω)是基函数h(t)的连续时间傅里叶变换。将频谱进行离散化处理，即令

公式(2)可以表示为：

显然，未知参数

可以通过获取K≥2L个非零傅里叶系数

得到。

获取傅里叶系数最简单的方法是采用Sinc采样核(即理想低通滤波器)，然而，该方法只能够获取一段连续的傅里叶系数子集，重构效果较差；Tur等人引入一种单通道预采样滤波过程来获取任意离散的傅里叶系数，但是这种滤波器的频响特性在实际应用中是难以实现的；Gedalyahu等人利用多通道混频积分结构直接采集输入信号的傅里叶系数，这种结构中每个通道只能采集一个傅里叶系数，硬件系统规模庞大，而且各通道之间的同步是一个很大的问题；为了获取离散分布的傅里叶系数，以色列理工学院的Eldar等人设计了一种多通道欠采样雷达接收机。然而，由于不同脉冲串信号频谱的差异性，以上的FRI采样方法需要根据具体的输入信号单独设计采样结构，效率较低，且系统的通用性较差。综上所述，研究适用于各种类型的脉冲串信号的有限新息率采样方法，是非常有意义的。

发明内容

本发明为解决现有脉冲串信号欠采样技术需要根据具体的输入信号单独设计采样结构，导致系统的通用性较差的问题，提供了基于随机解调和有限新息率的脉冲串信号欠采样方法。

本发明所述基于随机解调和有限新息率的脉冲串信号欠采样方法，通过以下技术方案实现：

步骤一、生成被测脉冲串信号x(t)和伪随机序列p(t)；

步骤二、对被测脉冲串信号x(t)和伪随机序列p(t)进行随机解调：首先，采用乘法器将被测脉冲串信号x(t)和伪随机序列p(t)进行混频，得到混频后的信号y(t)；然后，采用低通滤波器对混频后的信号y(t)进行滤波，得到滤波后的信号g(t)；最后，采用数据采集卡对滤波后的信号g(t)均匀采样得到样本g(n)；

步骤三、参数估计：利用步骤二中获得的样本g(n)估计被测脉冲串信号x(t)中的幅值参数和时延参数。

本发明最为突出的特点和显著的有益效果是：

本发明所涉及的基于随机解调和有限新息率的脉冲串信号欠采样方法，利用随机解调中的扩频技术来获取任意类型脉冲串信号的频谱信息。在模拟预处理阶段，采用伪随机序列对被测脉冲串信号进行混频，从而将任意频带的频谱信息扩展到基带，采用低通滤波器即可获取一组傅里叶系数。解决了由于不同脉冲串信号频谱的差异性导致现有FRI采样方法通用性差的问题。且本发明方法精度较高；仿真实验结果表明，在无噪声环境下，参数估计精确接近100％；在不同信噪比的高斯白噪声环境下，参数估计精度平均高于LPF-FRI方法和滤波器组-FRI方法约5％。

附图说明

图1为脉冲串信号欠采样系统结构框图；

图2为实施例中不同方法在噪声情况下重构效果对比曲线图。

具体实施方式

具体实施方式一：结合图1对本实施方式进行说明，本实施方式给出的基于随机解调和有限新息率的脉冲串信号欠采样方法，具体包括以下步骤：

步骤一、生成被测脉冲串信号x(t)和伪随机序列p(t)；

步骤三、参数估计：利用步骤二中获得的时域样本g(n)估计被测脉冲串信号x(t)中的未知幅值参数和时延参数。

具体实施方式二：本实施方式与具体实施方式一不同的是，步骤一中所述的伪随机序列p(t)中随机±1的跳变频率f_p需要满足f_p≥max{f_min,f_cut-f_max}，其中f_min和f_max分别表示被测脉冲串信号x(t)的最小频率和最大频率，f_cut为低通滤波器的截止频率。

伪随机序列p(t)按照以下公式生成：

其中，α_n∈{+1,-1}(n∈Z)为随机±1序列，Q为时域区间t∈[0,T)内伪随机序列p(t)的分隔段数。

其他步骤及参数与具体实施方式一相同。

具体实施方式三：本实施方式与具体实施方式一不同的是，步骤一中所述脉冲串信号x(t)具体为：

其中，h(t)是波形已知的脉冲，也称之为基函数；t∈[0,T)为观测时间，T为x(t)的持续时间长度；a_l为幅值参数，t_l为时延参数，t_l∈[0,T)；l＝1,2,...,L；L为脉冲的数目，则幅值参数与时延参数的数目均为L。

其他步骤及参数与具体实施方式一或二相同。

具体实施方式四：本实施方式与具体实施方式二不同的是，步骤二中所述混频后的信号y(t)的具体时域形式为：

y(t)＝x(t)·p(t) (5)

其中，t∈[0,T)为观测时间，T为x(t)的持续时间长度；

其他步骤及参数与具体实施方式二相同。

具体实施方式五：本实施方式与具体实施方式四不同的是，步骤二中所述采用低通滤波器对混频后的信号y(t)进行滤波，得到滤波后的信号g(t)的具体过程包括：

求得y(t)的傅里叶系数：

其中，

f_max是被测脉冲串信号x(t)的最大频率，

表示对数值向下取整；X[m]为被测脉冲串信号x(t)的傅里叶系数；P[k-m]为伪随机序列p(t)的傅里叶系数；j表示虚数单位；m＝-M,…,M；k∈Z，Z表示整数；

所述滤波后的信号g(t)对应的傅里叶系数为：

其中，

f_cut为低通滤波器的截止频率。

其他步骤及参数与具体实施方式四相同。

具体实施方式六：本实施方式与具体实施方式二不同的是，步骤二中所述均匀采样得到样本g(n)的具体操作为：

以f_s≥2f_cut为采样率对滤波后的信号g(t)均匀采样，得到样本

其中，

其他步骤及参数与具体实施方式二至五相同。

具体实施方式七：本实施方式与具体实施方式五或六不同的是，步骤三中所述参数估计具体包括以下步骤：

步骤三一、构建观测向量：对样本g[n]求离散时间傅里叶变换，得到傅里叶系数G[k]，然后构建观测向量：

g＝[G[-K],G[1-K],…,G[K]]^T (8)

步骤三二、观测时域区间离散化：将观测时间[0,T)平均划分为N个均匀的网格，

即量化间隔为△＝T/N；那么，任意时间变量能够近似为t≈n△，n＝0,1,…,N-1；任意被测脉冲串信号的时延参数能够近似为t_l≈n_l△，n_l为幅值向量s中的非零项，n_l∈Z；

步骤三三、构建观测矩阵：按照以下公式，构建观测矩阵Φ：

Φ＝PHΨ (12)

其中，

f_max是被测脉冲串信号x(t)的最大频率，

表示对数值向下取整；X[m]为被测脉冲串信号x(t)的傅里叶系数；

f_cut为低通滤波器的截止频率；P[m]为伪随机序列p(t)的傅里叶系数；H[m]为基函数h(t)的傅里叶系数，P、H、Ψ为中间矩阵；

步骤三四、求幅值参数向量的解：幅值参数向量s的求解能够转换为求解一个最小L0范数(L0范数指向量中非零元素的个数)下的优化问题：

利用正交匹配追踪算法对公式(13)进行求解，得到幅值参数向量的估计值

并从

中找出非零项

l＝1,…,L；L为脉冲的数目；

步骤三五、估计时延和幅值参数：被测脉冲串信号中的幅值参数a_l和时延参数t_l能够估计为：

其中，

为幅值参数的估计值，

为时延参数的估计值。

其他步骤及参数与具体实施方式五或六相同。

具体实施方式八：本实施方式与具体实施方式七不同的是，步骤三四中所述幅值参数向量s具体为：

s＝[s[0],s[1],…,s[N-1]]^T (15)

其中，

其他步骤及参数与具体实施方式一至七相同。

实施例

采用以下实施例验证本发明的有益效果：

本实施例按照如图1所示脉冲串信号欠采样系统进行，被测脉冲串信号的参数设置如下：基函数采用高斯脉冲的频移形式，

其中α＝3e-4，载频f_c＝2kHz；脉冲个数L＝4；信号持续时间长度T＝0.1秒，时延参数在区间[0,0.1)秒内随机选取；幅值参数在[0.5,1.5]内随机选取。伪随机序列p(t)的在时域区间[0,0.1)秒内的分段个数为Q＝1000。在采样阶段，低通滤波器的截止频率设置为f_cut＝500Hz，数据采集卡的采样率设置为f_s＝2f_cut＝1kHz；在参数估计阶段，模拟时间轴的最小量化单位为△＝0.1毫秒，因此划分的网格数为N＝T/△＝1000。

为了从数值上评估该方法的性能，采用均方误差作为评价指标，为了方便比较对均方误差取对数：

其中，MSE表示均方误差；L是脉冲的数量，t_l是真实的时延参数，

是估计的时延参数。由于幅值参数的误差与时延参数的误差成正比，因此只需要用时延参数的均方误差来衡量该方法的性能。

(一)在无噪声环境下，验证本发明方法的有效性，实验结果如表1所示。从表中可以看出，对于L＝4个脉冲的脉冲串信号，本发明方法能够实现时延参数的精确估计，准确率接近100％。

表1 无噪声情况下时延参数的估计值(秒)

方法	t<sub>1</sub>	t<sub>2</sub>	t<sub>3</sub>	t<sub>4</sub>
					真实值	0.02	0.05	0.07	0.09
本发明方法	0.02	0.05	0.07	0.09

(二)在不同信噪比的高斯白噪声环境下(信噪比SNR由-20dB增加到80dB)，采用本发明方法和最近提出的低通滤波器-有限新息率LPF-FRI方法、滤波器组-FRI方法(基于4路多通道系统)进行比较，重复实验100次后取平均结果，仿真实验结果如图2所示。从图2中可以看出，本发明所述基于随机解调和有限新息率的脉冲串信号欠采样方法在噪声情况下有很好的表现，参数估计精度要优于LPF-FRI方法和滤波器组-FRI方法，平均约高出5％。

本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。