专利名称：一种带通信号调制宽带转换器采样的重构方法
技术领域：
本发明涉及信号采样领域，特别涉及频域稀疏信号调制宽带转换器采样的重构方法。
背景技术：
Shannon/Nyquist采样理论指出，对信号的采样频率至少为信道带宽的两倍。近年来随着压缩感知(也被称为“压缩采样”)技术的发展，人们发现当信号是稀疏的或者是可压缩的，贝1J可以以低于Nyquist采样率的速率对信号进行采样，由此提出的以低于Nyquist速率的采样速率对模拟信号进行采样的方法称为子Nyquist采样方法。Mishali 等人在论文“From theory to practice:Sub-Nyquist sampling ofsparse wideband analog signals^Ci^jAL, IEEE Journa丄 of Selected Topics in bignal Processing, vol. 4，no. 2，pp. 375-391，April 2010)中提出的调制宽带转换器(MWC)采样就是ー种子Nyquist采样方法。MWC的结构如图2所示,具体来讲,频率范围为
的信号x(t)同时输入m个通道。在第i个通道，信号x(t)与周期为Tp的混频函数Pi (t)相乘，接着通过ー个截止频率为Iパ2TS)的低通滤波器，滤波后的信号以fs = 1/TS的速率进行采样。由于1/TS足够小，所以现有的商用ADC能用来完成采样过程。另外，该采样系统的ー个重要特点就是m个通道的采样速率之和仍远小于Nyquist采样率，即mfs〈〈fNYQ。现在来分析采样序列yi [n]与未知信号X (t)之间的关系。令fp = I/Tp, fs = I/Ts,Fp = [-fp/2, +fp/2],Fs = [_f■ノ2，+fs/2]。考虑第 i 个通道。由于 PiU)是周期イ■号，所以其Fourier展开式为
IMt
00/-It/；.(/) = ^ ' CjlC ^(I)
/=—00
I c T-プと"其中c =—J ' P1(I)C '"ぬ。对于MWC，有如下表达式y (f) = Az (f), f G Fs (2)其中y (f)是mX I维向量,其第i个元素为Yi[n]的离散时间Fourier变换，即yiif) = Yi(ej27tfTs')，I 彡 i 彡 m。未知矩阵 z (f) = [Z1 (f), z2(f),..., zL(f)]T 的长度为 L =
2L0+1,￡o = (/nyq +^)/(2/^) 一I’「o']表示大于或等于a且与a最接近的整数,其中Zi (f) =X(f+ (i-L0-l) fp), I ^ i ^ L, f G Fs (3)mXL矩阵A (被称为感知矩阵)第i行第I列元素an为aa =C,-,z0+i-/3 I 彡 i 彡 m, I 彡 I 彡 2L0+1 (4)由式3可知，只要求得z(f)，就能得到原信号。z(f)的求解依赖于式2的求解，即在已知y(f)和A的情况下求解z(f)。式2所示的模型称为重构模型。由于mXL矩阵A中，m<L，因此，该问题是ー个欠定问题。当已知z(f)为稀疏的情况下，可以采用压缩感知理论中的信号重构方法进行求解。Mishali等人在文献“From theory topractice:Sub-Nyquist sampling of sparse wideband analog signals，，中给出 J 一种具体的求解方法，表明了重构的可行性。然而，Mishali等人给出的重构方法针对的是低通信号，即频率范围为
的信号。而在实际中，待采样信号往往是带通信号，即频率范围为[fstart，fmd]的信号，其中fstart > O0如果将带通信号当成频率范围为
的低通信号来处理，则重构时感知矩阵A以及未知向量z(f)的维度将较大，重构复杂度较高。本发明提出ー种针对带通信号的MWC采样重构方法，相比于上述处理过程，感知矩阵和待求未知向量维度将降低，从而具有重构复杂度低的优势。
发明内容
本发明公开了ー种带通信号调制宽带转换器采样的重构方法，其包括以下步骤(a)构造感知矩阵k*，使得感知矩阵k*中的元素为
、卜,んー丨 _,，1</</; ,1</</./2
Ic、_L「卜 \<i<m,L!2 + \<l <L
_7]其中，厶=认マす1、-1も ニ 2f^— I，L = 2(L2-Li+1)，f加和し
分别为待采样频段的起始频率和终止频率，fs = 1/TS和fp = 1/TP ；(b)构造频域未知向量其中ド/十！-卿’2J,Fs
[X(f + (i-L + L2)fp), 772 + 1 </</../€/ ；其中X(f)为带通信号X (t)的Fourier变换，<(/)代表的是X (f)中以(i-l_L2)fp (I彡i彡L/2吋)或(i_L+L2)fp (L/2+1彡i彡L吋)为中心、以fs为宽度的频段移到基带对应的分量；(C)获得频域未知向量的值，其中根据对带通信号X (t)进行调制宽带转换器采样后得到的采样序列Ji [n]，可得XV2机)=2>,ネ(/-/ん)+ 2>"1(/ —//於)，f G Fs
I^-L2I=L1其中ム=1Lart丁 -UL2= 1Ljtjs -I；
2fp 2fp将该式写成矩阵形式，为y (f) = A*z* (f), f G Fs其中バわ是mX I维向量，其第i个元素为タ,(/) = )^バ嗔)，I ( i ( m，根据已知的y(f)和感知矩阵广获得未知的Z(f)；(d)根据映射关系得到原带通信号，其中根据X(f)与z*(f)的对应关系，得到X(f)，将其做反Fourier变换得到原带通信号x (t)。
其中，在步骤(C)中，采用压缩感知领域的正交匹配追踪法来获得频域未知向量z*(f)。本发明可取得以下有益效果本发明提出的针对带通信号的MWC采样的重构模型，相比于将带通信号当成低通信号来处理的MWC采样重构模型，未知量个数更少，感知矩阵维数更小，因此，所需的计算复杂度更低，更加适合于实时应用。
图I为本发明的重构方法流程图。图2为调制宽带转换器采样框图。图3为本发明重构中未知向量与原信号频谱的对应关系。
具体实施例方式下面结合附图和具体实施例对本发明的带通信号调制宽带转换器采样的重构方法作进ー步详细说明。图I是本发明的带通信号调制宽带转换器采样的重构方法的流程图，如图所示，本发明的重构方法包括以下步骤(I)构造感知矩阵
271 h正如在背景技术中介绍的，系数C;7 = ニp乃(/> Tp A，由此构造感知矩阵A%
if
使得感知矩阵A*中的元素为
*_|ca2+i-/5 \<i<m,\<l<LI2じ'其中，A= lfst;+L — I，L2 = 2マ + ハ-1，L = 2 (L2-L1+l)，fstart 和 fend
J pJ p
分别为待采样频段的起始频率和终止频率，fs和fp的定义与背景技术中的相同。(2)构造未知向量构造频域未知向量/(/)= [Zj(/),Z*(/),...,Z*(/)f ,其中如=ド(ハ(トレ他)，l^L,2J,Fs⑷
\x(f + (i-L + L7)fp), L/2 + l<i<L,feFs其中X(f)为带通信号x(t)的Fourier变换。由式(6)可知，:^(/)代表的是X(f)中以(I-I-L2) fp (I彡i彡L/2吋)或(i-L+L2) fp (L/2+1彡i彡L吋)为中心、以fs为宽度的频段移到基带对应的分量。(3)获得未知向量的值根据对带通信号x(t)进行调制宽带转换器(MWC)采样后得到的采样序列yi[n]，可得
Y1 {e)2wK) = CilX(/ -Ifp) + ^jCilX(/-IfpIf ^ Fs (7)
I=-L2I=L^其中LpL2与式(5)中定义相同。将式(7)写成矩阵形式，为y(f) = A*z* (f), f G Fs (8)其中，y(f)是mX I维向量，其第i个元素为ス.(/)=ぢレプ2桃)，I彡i彡m。在该式中，要解决的问题是根据已知的y(f)和A%获得未知的z*(f)。采用压缩感知领域的已有方法，如正交匹配追踪法，可以求解得到z*(f)。(4)根据映射关系得到原带通信号
求出z*(f)后，根据式6中X(f)与z*(f)的对应关系，可以得到X(f)，将其做反Fourier变换可得到x (t)，此即完成了原带通信号x (t)的重构，即得到原带通信号x(t)。下面以ー个具体例子来说明本发明的重构方法。设射频带通信号X (t)的频率范围为[500MHz，1000MHz]。对该带通信号x (t)进行如图2所示的调制宽带转换器(MWC)采样，其中采样通道数取m = 24，各个Pi (t)周期的倒数fp = 20MHz，每个通道采样频率fs = 20MHz，各个通道的低通滤波器的截止频率为IOMHz。根据对带通信号x(t)进行MWC采样后得到的采样序列yi[n]，I < i <m，本发明采用如下步骤进行重构(I)构造感知矩阵由在本例子中给出的上述參数，可得到^!= lj: +わ-1 = 26 ,
ムJ P
しニ 2^^+fs _1 = 51感知矩阵A*中元素为
Vp
* I ら:—z， \<i<m.}< I <26
11~ ， 1</<ot,27</<52
I-/ 其中じ"=—Ip XOc Tp 治，Tp = l/fp = 0. 5uSo TPJo(2)构造未知向量未知向量Z*(/)= [Z:(/),Z;(/),.. ,<(/)f的长度为L = 52，各个元素为 , [ズ(/. + 20(/-52)). l</<26,/'e[-20.20]MHzZi (/) = <'
1 [ Xi f + 20(/ -1)), 27 < i < L,./. e [-20，20]MHz其与原信号对应关系示意图如图3所示。由图可知，原信号Fourier变换X(f)与z*(f)之间存在一对ー的映射关系。(3)获得未知向量的值各个通道MWC采样后数字序列的离散Fourier变换为W2吼)=X CilXif - 20/) +— 20/), f G [-20，20] MHz
I=-L2I=L1将其写成矩阵形式
y(/)24xl = K4x52Z (f)52xl, f G [-20，20]MHz根据压缩感知领域的已有方法，如正交匹配追踪法，求得z*(f)。(4)根据映射关系得到原带通信号求出z*(f)后，根据X(f)与z*(f)的对应关系，可以得到X(f)，将其做反Fourier变换可得到X(t)，此即完成了原带通信号x(t)的重构，即得到原带通信号x(t)。作为比较，若将该带通信号作为低通信号(频率范围变成
MHz)处理，按相同的MWC进行采样，则根据公式2进行重构时，未知矩阵z (f) = [Zl(f),z2(f),...,zL(f)]
T的长度为L = 2LQ+1，4=「(/胃+./:)/—1 = 50所以L=101。感知矩阵A维度为
24X101。由此可知，本发明重构方法中感知矩阵A*维度要小于按低通信号处理的感知矩阵 A的维度，而且本发明未知向量z*(f)长度小于按低通信号处理的未知向量z (f)的长度，因此按本发明的重构方法进行重构计算的复杂度更低。
权利要求
1.一种带通信号调制宽带转换器采样的重构方法，其包括以下步骤 Ca)构造感知矩阵A%使得感知矩阵A*中的元素为
2.根据权利要求I所述的方法，其特征在于，在步骤(c)中，采用压缩感知领域的正交匹配追踪法来获得频域未知向量z*(f)。
全文摘要
本发明公开了一种带通信号调制宽带转换器采样的重构方法，其包括构造感知矩阵，构造频域未知向量，获得频域未知向量的值，和根据映射关系得到原带通信号。
文档编号H04L25/02GK102801665SQ20121029912
公开日2012年11月28日 申请日期2012年8月21日 优先权日2012年8月21日
发明者杨小牛, 郑仕链 申请人:中国电子科技集团公司第三十六研究所