基于混合压缩采样的非合作超宽带信号的重构方法

基于混合压缩采样的非合作超宽带信号的重构方法
【技术领域】
[0001] 本发明属于信号处理技术领域，特别涉及一种非合作超宽带信号的重构方法，可 用于超宽带信号的低速获取。
【背景技术】
[0002] 现代意义上的超宽带技术，出现于20世纪60年代。GeraldF.Ross博士开创性的 论证了超宽带技术在雷达和通讯方面的可行性。经过50多年的发展，超宽带技术广泛应用 于无线电通讯、雷达、精确定位、成像等领域，并且与常规无线电技术相比，超宽带技术具有 抗衰落、容量大、功率低等优点。
[0003] 在常规无线电技术的应用中，根据Nyquist采样定理，获得超宽带信号需要非常 高的采样频率。高频率的信号采样增加了采样系统设计的困难，而且离散化后产生的大量 数据也给存储与传输处理带来了巨大压力。因此需要储存量更大的硬件设备和超高计算能 力的计算机来处理这些数据，关键问题在于这些硬件资源是有限的。
[0004] 由于超宽带信号需要非常高的采样频率，这就需要超高速率、高分辨率的数模转 换器ADC对模拟信号进行采样。现有的一些文献中设计出了速度可达数百兆/每秒、分辨 率超过10位的A/D转换器（ADC)，但是当所需要的采样速度高达千兆/每秒时，由于A/D转 化过程中量化速度的限制，ADC芯片的高速度和高精度两个指标难以同时得到满足。当A/ D转换速率高达千兆/每秒时，现有A/D转换器的量化分辨率会降低几个量级，因此这些类 型的ADC无法满足高分辨率的需求。
[0005] 为了解决以上问题，近年来在信号处理领域兴起的压缩感知理论吸引了诸多学者 的关注，该理论指出：只要信号在某一个空间上具有稀疏性，可以利用观测矩阵以远低于奈 奎斯特的采样速率对其进行观测，并使用优化计算方法可高概率的实现信号的重构。这样 大大降低了传感器的采样成本，从而设计出具有高概率重构信号的低速采样结构，例如模 拟信息转换器AIC的设计。
[0006] 但是压缩感知理论有其可行性前提：已知信号在某个变换域中存在稀疏的表示。 在非合作性的环境下，若待处理信号缺乏足够先验信息，对信号进行稀疏表示是很困难的。 因此，要在非合作性的环境下用压缩感知的原理观测信号，就必须首先获取信号的先验信 息。

【发明内容】

[0007] 本发明的目的在于提出一种基于混合压缩采样的非合作超宽带信号的重构方法， 以解决现有技术在非合作性环境下的，因为缺乏对信号的足够先验，而不能直接应用压缩 感知理论通过观测数据有效重构原信号的问题。
[0008] 本发明的技术关键是：对输入的模拟超宽带信号分两个方面进行观测，一是采用 模拟信息转换器AIC对信号的随机投影观测，二是在变换域基函数上的投影观测，以获得 原始超宽带信号的部分先验信息，再利用压缩感知的方法对信号进行精确重构。具体实现 步骤包括如下：
[0009] (1)将超宽带信号s(t)进行随机投影观测，得到随机观测系数向量y和随机观测 矩阵A;将超宽带信号s(t)在频率变换域的n个正交基函数Odt)上分别进行投影观测， 得到n个正交投影观测系数ai，1 = 1，2,. . .，n;
[0010] (2)将投影观测系数&1与对应的正交基函数进行线性组合，得到粗略重构信号s'⑴：
[0011]
[0012] (3)对粗略重构信号s'（t)进行K-L分解，得到信号自相关矩阵的特征向量矩阵 V，并把该特征向量矩阵V作为信号的稀疏矩阵，利用随机观测系数向量y和随机观测矩阵 A，根据压缩感知理论，求解优化问题，得到优化系数向量i:
[0013]
[0014] 其中，x为待求的系数向量，A为正则化参数，| | | |2表示向量的2范数，| | | |i表示向量的1范数；
[0015] (4)根据优化系数向量i和特征向量矩阵V，得到最终精确重构信号为：
[0016]互= 7-4。
[0017] 本发明与现有技术相比有如下优点：
[0018] 本发明由于将基于压缩感知的随机投影观测和频率域的正交投影观测两种不同 采样结构相结合，进行优势互补，解决了在非合作性环境下不能利用压缩感知理论对信号 进行压缩处理的问题，不仅很大程度的降低采样率，而且提高了对信号重构的精确度。
[0019] 仿真结果表明：与传统的Nyquist采样方法相比本发明的采样方法可以大大降低 信号的采样速度，并保证信号的高精度恢复。
【附图说明】
[0020] 图1为本发明的实现原理图；
[0021] 图2为本发明中随机投影观测的示意图；
[0022] 图3为本发明中正交投影观测的示意图；
[0023] 图4为仿真实验的对比结果；
【具体实施方式】
[0024] 以下结合附图对本发明的实施力和效果做进一步详细描述。
[0025]参照图1，本发明的具体实现步骤如下：
[0026] 步骤1:对超宽带信号s(t)进行混合观测
[0027] (1)将超宽带信号s(t)进行随机投影观测，得到随机观测系数向量y和随机观测 矩阵A:
[0028] 参照图2,本步骤的具体实现如下：
[0029](la)将超宽带信号s(t)与伪随机序列函数p(t)进行调制后进入积分器，积分后 的信号使用低速的数模转换器ADC进行均匀的采样，随机投影采样数据向量.小〃]表示为：
[0030]
[0031] 其中，At为随机投影观测结构中数模转换器ADC的采样时间间隔，m为随机投影 观测值的下标，M为总观测点数。
[0032] (lb)对随机投影采样数据向量通过数模转换器ADC量化，得到最终输出的随 机观测系数向量y;
[0033] (lc)根据随机投影观测结构中的伪随机序列函数p(t)的幅度只有{+1，_1}两个 值的特性，超宽带信号与伪随机序列函数P(t)进行随机调制，等效为对超宽带的离散信号 乘以对角矩阵D:
[0034]
[0035] 其中￡i= {+1，_1}表示随机序列p(t)中第i个调制符号的幅度，N为随机序列 调制符号的总数。假定，在随机投影观测中数模转换器ADC的采样率为超宽带信号Nyquist 采样率的1/3,则积分器对应的矩阵形式为：
[0036]
[0037] 因此，对超宽带信号进行随机观测的随机观测矩阵A为：
[0038] A=HD
[0039] (2)将超宽带信号s(t)在频率变换域的n个正交基函数Odt)上分别进行投影 观测，得到n个正交投影观测系数ai，1 = 1，2,. . .，n。
[0040] 参照附图3,本步骤的具体实现如下：
[0041] (2a)将超宽带信号s(t)与频率域的基函数(t)进行调制后通过积分器进行积 分，再使用数模转换器ADC采样得到正交投影采样数据尽：
[0042]
[0043] 其中T为正交投影结构中数模转换器ADC的采样间隔，n为所选取的频率域基函 数的个数；
[0044] (2b)将正交投影采样数据4经过数模转换器ADC进行量化，得到n个正交投影观 测系数4。
[0045] 步骤2 :信号的粗略重构
[0046] 将步骤（2b)中的投影观测系数与对应的正交基函数〇dt)进行线性组合，得 到粗略重构信号S'（t):
[0047]
[0048] 步骤3 :基于压缩感知理论的优化系数向量求解
[0049] 对步骤二中粗略重构信号s'（t)进行K-L分解，得到信号自相关矩阵的特征向量 矩阵V，并把该特征向量矩阵V作为信号的稀疏矩阵，利用步骤（lb)中得到的随机观测系数 向量y和步骤（1C)中得到的随机观测矩阵A，根据压缩感知理论，求解优化问题，得到优化 系数向量文：
[0050]
[0051] 其中，x为待求的系数向量，人为正则化参数，| | | |2表示向量的2范数，| | |Ii
[0052] 表示向量的1范数。
[0053] 步骤4:信号的精确重构
[0054] 根据步骤三中得到的优化系数向量i和特征向量矩阵V，得到最终精确重构信号 为：
[0055] 5=7-4。
[0056] 本发明的效果可以通过以下仿真实验说明：
[0057] 1?仿真参数
[0058] 超宽带信号选择高斯调制的正弦脉冲信号，中心频率：;^= 2GHz;高斯脉冲在低 于最大值40dB处截断；截断的信号时长T约为2ns;脉冲在-10dB处的带宽B约为1. 6GHz; 采样频率fs= 50GHz;数模转换器ADC的采样量化位数为8位。
[0059] 在频率变换域正交投影观测，在输入信号的带宽内均匀设置投影基函数，投影基 函数的数量满足〃确保在粗略重构信号时获得原超宽带信号的大部分信息。
[0060] 2.仿真内容与结果
[0061] 用本发明的方法和传统Nyquist采样方法分别的对上述的信号进行采样，用本发 明的方法首先对信号进行粗略重构，再对信号进行精确的重构，精确重构的信号与Nyquist 采样的信号对比结果如图4,其中
[0062] 图4(a)是仿真使用的高斯脉冲信号及频谱，其在-10dB的带宽大于1GHz，是一个 标准的超宽带信号；
[0063] 图4(b)是用本发明方法粗略重构后的信号与Nyquist采样方法的信号的对比图。 从图4 (b)可见本发明粗略重构的信号只有波形与Nyquist采样方法的信号保持一致，说明 它含有原超宽带信号的大部分信息，但粗略重构的信号信噪比SRN只有40dB，只能粗略的 重构信号；
[0064] 图4(c)是用本发明精确重构后的信号与Nyquist采样方法的信号的对比图。从 图4(c)可见本发明精确重构的信号与Nyquist采样方法的信号完全一致，精确重构的信号 的信噪比SRN达到64dB，表明本发明能精确的重构信号。
[0065] 在上述仿真试验中，用传统Nyquist的采样点数为Nny(iusit= 95。用本发明方法在 频率域正交投影观测中采样点的个数为q=「32^ = 4，随机投影观测中的采样点数为n2 = 10,总的采样点数仅为1^=1^+%= 14。对比这些参数表明：本发明方法对信号有压缩 功成，其f曰5压缩比为H=Nny(iuist/Nmixed= 6. 8。
[0066] 通过本仿真实验的结果表明本发明的方法可以大大降低信号的采样速度，并保证 信号的高精度恢复。
【主权项】
1. 一种基于混合压缩采样的非合作超宽带信号的重构方法，包括以下步骤： (1) 将超宽带信号S(t)进行随机投影观测，得到随机观测系数向量y和随机观测矩阵 A ;将超宽带信号s (t)在频率变换域的η个正交基函数O1U)上分别进行投影观测，得到 η个正交投影观测系数七，1 = 1，2,. . .，η ; (2) 将投影观测系数a,与对应的正交基函数进行线性组合，得到粗略重构信号s'（t):(3) 对粗略重构信号s'（t)进行K-L分解，得到信号自相关矩阵的特征向量矩阵V，并 把该特征向量矩阵V作为信号的稀疏矩阵，利用随机观测系数向量y和随机观测矩阵A，根 据压缩感知理论，求解优化问题，得到优化系数向量h其中，X为待求的系数向量，λ为正则化参数，Μ M2表示向量的2范数，Μ M1表示 向量的1范数； (4) 根据优化系数向量f和特征向量矩阵V，得到最终精确重构信号为： S=V xO2. 根据权利要求1所述的基于混合压缩采样的非合作超宽带信号的重构方法，其中所 述步骤（1)中的对超宽带信号s(t)进行随机投影观测，得到随机观测系数向量 y和随机观 测矩阵A，按照如下步骤进行： (la) 将超宽带信号s (t)与伪随机序列函数p (t)进行调制后进入积分器，积分后的信 号使用低速的数模转换器ADC进行均匀的采样，随机投影采样数据向量表示为：其中，Δ t为随机投影观测结构中数模转换器ADC的采样时间间隔，m为随机投影观测 值的下标，M为总观测点数。 (lb) 对随机投影采样数据向量通过数模转换器ADC量化，得到最终输出的随机观 测系数向量y; (lc) 根据随机投影观测结构中的伪随机序列函数p(t)的幅度只有{+1，-1}两个值的 特性，超宽带信号与伪随机序列函数P (t)进行随机调制，等效为对超宽带的离散信号乘以 对角矩阵D :其中ε i= {+1，-1}表示随机序列p (t)中第i个调制符号的幅度，N为随机序列调制 符号的总数。假定，在随机投影观测中数模转换器ADC的采样率为超宽带信号Nyquist采 样率的1/3,则积分器对应的矩阵形式为：因此，对超宽带信号进行随机观测的随机观测矩阵A为： A = HD03.根据权利要求1所述的基于混合压缩采样的非合作超宽带信号的重构方法，其中所 述步骤（1)中的对超宽带信号s (t)进行正交投影观测，得到η个正交投影观测系数ai，l = 1，2,. . .，n，按照如下步骤进行： (Id)将超宽带信号s (t)与频率域的基函数O1U)进行调制后通过积分器进行积分， 再使用数模转换器ADC采样得到正交投影采样数据矣： U1 = ￡ Λ'(/)Φ;(/)?// /=1,2,...,/7 其中T为正交投影结构中数模转换器ADC的采样间隔，η为所选取的频率域基函数的 个数； (If)将正交投影采样数据A经过数模转换器ADC进行量化，得到η个正交投影观测系 数七。
【专利摘要】本发明公开了一种基于混合压缩采样的非合作超宽带信号的重构方法，主要解决现有技术无法利用有效先验对非合作性信号进行重构的问题，其实现步骤为：(1)对输入的模拟信号进行随机投影观测和在变换域基函数上的投影观测；(2)利用投影观测系数与基函数的线性组合实现对信号的粗略重构；(3)利用主成分分析法对粗略重构信号进行K-L分解，得到信号自相关矩阵的特征向量；(4)以特征向量矩阵作为信号稀疏的域，在压缩感知理论框架下实现信号的精细重构。本发明能够在非合作环境下，实现对超宽带信号的低速采样和精确重构，可用于超宽带信号的低速获取。
【IPC分类】H03M7/30
【公开号】CN104883192
【申请号】CN201510340783
【发明人】林杰, 徐同飞, 石光明, 贺玉高, 赵光辉, 王晓甜
【申请人】西安电子科技大学
【公开日】2015年9月2日
【申请日】2015年6月18日