一种信号重建方法、装置及系统与流程

本发明涉及信号处理技术领域，特别涉及一种信号重建方法、装置及系统。

背景技术：

脉冲uwb(ultrawideband，超宽带)通信系统通过发送和接收亚纳秒级的脉冲信号来传输信息。由于脉冲uwb信号的带宽非常大，按照传统的香农采样定理，对其进行数字化采样，需要设计采样率高达ghz的adc(analog-to-digitalconverter，模数转换器)。例如，当脉冲宽度为0.5纳秒时，脉冲uwb信号的带宽约为2ghz，这时，基带信号处理算法要求adc采样率大于4ghz。就目前的半导体集成电路技术工艺水平而言，设计这样高采样率，而又低成本、低功耗的adc，是非常困难的。因而，在奈奎斯特采样条件难以满足的条件下，研究对脉冲uwb信号进行欠奈奎斯特采样，然后通过后续的数字信号处理算法估计信号参数，并基于估计的信号参数重建信号，具有重要的实用价值。

技术实现要素：

本发明实施例提供了一种信号重建方法、装置及系统，能够基于欠奈奎斯特采样重建信号。所述技术方案如下：

第一方面，提供了一种信号重建方法，所述方法包括：

在对接收到的超宽带信号进行带通滤波后，对带通滤波后的超宽带信号进行欠采样，得到离散时间信号；

基于所述离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息，所述超宽带信号基于所述冲激串信号产生；

基于估算出的冲激串信号的幅度和时移信息，计算所述冲激串信号；

基于计算得到的所述冲激串信号，对所述接收到的超宽带信号进行重建。

可选地，所述带通滤波采用的带宽大于或者等于所述超宽带信号的新息率，欠采样频率大于或等于所述带通滤波采用的带宽的两倍。

可选地，所述基于所述离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息，包括：

通过傅里叶变换将所述离散时间信号变换到频域；

对频域信号进行迫零均衡处理；

基于迫零均衡处理后的频域信号，采用所述旋转不变性技术算法估算所述冲激串信号的幅度和时移信息。

可选地，所述基于迫零均衡处理后的频域信号，采用所述旋转不变性技术算法估算所述冲激串信号的幅度和时移信息，包括：

基于所述迫零均衡处理后的频域信号，获取设定数量的、且连续的频率采样点；

基于所述设定数量的、且连续的频率采样点，确定所述超宽带信号的入射角信息；

基于所述超宽带信号的入射角信息，确定所述冲激串信号的时移信息和幅度。

第二方面，提供了一种信号重建装置，所述装置包括：

采样模块，用于在对接收到的超宽带信号进行带通滤波后，对带通滤波后的超宽带信号进行欠采样，得到离散时间信号；

估算模块，用于基于所述离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息，所述超宽带信号基于所述冲激串信号产生，并基于估算出的冲激串信号的幅度和时移信息，计算所述冲激串信号；

重建模块，用于基于计算得到的所述冲激串信号，对所述接收到的超宽带信号进行重建。

可选地，带通滤波采用的带宽大于或者等于所述超宽带信号的新息率，欠采样频率大于或等于带通滤波采用的带宽的两倍。

可选地，所述估算模块用于，

通过傅里叶变换将所述离散时间信号变换到频域；

对频域信号进行迫零均衡处理；

基于迫零均衡处理后的频域信号，采用所述旋转不变性技术算法估算所述冲激串信号的幅度和时移信息。

可选地，所述估算模块用于，

基于所述迫零均衡处理后的频域信号，获取设定数量的、且连续的频率采样点；

基于设所述定数量的、且连续的频率采样点，确定所述超宽带信号的入射角信息；

基于所述超宽带信号的入射角信息，确定所述冲激串信号的时移信息和幅度。

第三方面，提供了一种信号重建装置，所述装置包括处理器和存储器，所述存储器中存储有至少一条指令，所述指令由所述处理器加载并执行以实现如第一方面所述的信号重建方法所执行的操作。

第四方面，提供了一种信号重建系统，所述系统包括：

带通滤波器，用于对接收到的超宽带信号进行带通滤波；

信号重建装置，用于对带通滤波后的超宽带信号进行欠采样，得到离散时间信号；基于所述离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息，所述超宽带信号基于所述冲激串信号产生；基于估算出的冲激串信号的幅度和时移信息，计算所述冲激串信号；基于计算得到的所述冲激串信号，对所述接收到的超宽带信号进行重建。

本发明实施例提供的技术方案带来的有益效果是：通过在对接收到的超宽带信号进行带通滤波后，对带通滤波后的超宽带信号进行欠采样，得到离散时间信号；基于所述离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息；基于估算出的所述冲激串信号的幅度和时移信息，计算所述冲激串信号；由于所述超宽带信号基于所述冲激串信号产生，因此基于计算得到的所述冲激串信号，能够对所述接收到的超宽带信号进行重建；旋转不变性技术算法在重建信号时抗噪声性能较优，所生成的解能准确估计出信号的时移信息和幅度，具有较好的工程实践性和实用价值。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例提供的一种信号重建系统的结构示意图；

图2是本发明实施例提供的一种信号重建装置的结构示意图；

图3是本发明实施例提供的信号重建系统的工作流程示意图；

图4和图5是本发明实施例提供的一种信号重建方法的流程图；

图6是本发明实施例提供的带通滤波后的uwb信号对应的频域信号的示意图；

图7是本发明实施例提供的带通采样频率的范围的示意图；

图8是本发明实施例提供的相邻采样频率窗的示意图；

图9是本发明实施例提供的一种信号重建装置的结构示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。

为便于对本发明实施例的理解，首先对涉及的名词进行解释。

旋转不变性技术算法：指esprit(estimatingsignalparameterviarotationalinvariancetechniques，基于旋转不变技术的信号参数估计)算法，其原理主要是基于信号的旋转不变性估算信号参数。

脉冲信号的时移信息：指脉冲在时间轴上的位置。

脉冲信号的幅度：指脉冲的振幅。

超宽带信号的入射角信息：指超宽带信号入射至天线阵列的方位方向。

新息率：指单位时间内信号的自由参量的个数。其中，信号可以由有限数目的自由参量表示。

图1示出了本发明实施例提供的一种信号重建系统。参见图1，系统包括：带通滤波器11和信号重新装置10。

带通滤波器11，用于对接收到的超宽带信号进行带通滤波。

信号重新装置10，用于，对带通滤波后的超宽带信号进行欠采样，得到离散时间信号；基于离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息；基于估算出的冲激串信号的幅度和时移信息，计算冲激串信号；基于计算得到的冲激串信号，对接收到的超宽带信号进行重建。

图2示出了本发明实施例提供的一种信号重建装置，参见图2，该信号重新装置10包括采样模块12、估算模块13和重建模块14。

采样模块12，用于对带通滤波后的超宽带信号进行欠采样，得到离散时间信号。

估算模块13，用于基于离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息，超宽带信号基于冲激串信号产生，并基于估算出的冲激串信号的幅度和时移信息，计算冲激串信号。

重建模块14，用于基于估算出的冲激串信号的幅度和时移信息，对接收到的超宽带信号进行重建。

需要说明的是，信号重新装置10可以由集成在一起的控制芯片实现。

参见图3，uwb信号s(t)由脉冲uwb通信系统90发送。脉冲uwb通信系统90包括冲激串产生器90a和高斯单脉冲产生器90b。冲激串产生器90a产生冲激串信号x(t)。高斯单脉冲产生器90b基于冲激串产生器90a产生的冲激串信号x(t)，产生uwb信号s(t)。其中，uwb信号s(t)由高斯脉冲的一阶导数生成。

其中，带通滤波器11(bpf)的中心频率为fc，脉冲uwb信号s(t)的频谱在频率fc处具有最大信噪比，其噪声主要是高斯白噪声e(t)。

在uwb信号的一个数据帧内，共有p个高斯脉冲，帧长度t＝p×tf，tf为高斯脉冲周期，包含2p个自由参数幅度{di}和时移{ti}，i＝1,2,…,p，则uwb信号的信息率为ρ＝2p/t。新息率由单位时间内的脉冲个数来获取。在本实施例中，带通滤波器11的带宽b大于或者等于uwb信号s(t)的新息率ρ，即b≥ρ。采样模块12欠采样频率fs大于或等于带通滤波器11的带宽b的两倍。参见图3，uwb信号s(t)经过带通滤波器11带通滤波输出信号h(t)，并以欠采样频率fs采样后，得到离散时间信号y[n]。离散时间信号y[n]经傅里叶变换(dft)变换至频域后，输出频域信号y[k]；频域信号y[k]通过迫零均衡处理，输出冲激串信号x[k]。

示例性地，估算模块13用于，通过傅里叶变换将离散时间信号变换到频域；对频域信号进行迫零均衡处理；基于迫零均衡处理后的频域信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息。

示例性地，估算模块13用于，基于迫零均衡处理后的频域信号，获取设定数量的、且连续的频率采样点；基于设定数量的、且连续的频率采样点，确定超宽带信号的入射角信息；基于超宽带信号的入射角信息，确定冲激串信号的时移信息和幅度。

本发明实施例通过在对接收到的超宽带信号进行带通滤波后，对带通滤波后的超宽带信号进行欠采样，得到离散时间信号；基于离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息；基于估算出的冲激串信号的幅度和时移信息，计算冲激串信号；由于超宽带信号基于冲激串信号产生，因此基于计算得到的冲激串信号，能够对接收到的超宽带信号进行重建；旋转不变性技术算法在重建信号时抗噪声性能较优，所生成的解能准确估计出信号的时移信息和幅度，具有较好的工程实践性和实用价值。

对上述实施例提供的信号重建装置及系统未进行详尽描述的地方请参照下文方法实施例。

图4示出了本发明实施例提供的一种信号重建方法，参见图4，该方法包括如下步骤。

步骤301、在对接收到的超宽带信号进行带通滤波后，对带通滤波后的超宽带信号进行欠采样，得到离散时间信号。

步骤302、基于离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息。

其中，超宽带信号基于冲激串信号产生。

步骤303、基于估算出的冲激串信号的幅度和时移信息，计算冲激串信号。

步骤304、基于计算得到的冲激串信号，对接收到的超宽带信号进行重建。

本发明实施例通过在对接收到的超宽带信号进行带通滤波后，对带通滤波后的超宽带信号进行欠采样，得到离散时间信号；基于离散时间信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息；基于估算出的冲激串信号的幅度和时移信息，计算冲激串信号；由于超宽带信号基于冲激串信号产生，因此基于计算得到的冲激串信号，能够对接收到的超宽带信号进行重建；旋转不变性技术算法在重建信号时抗噪声性能较优，所生成的解能准确估计出信号的时移信息和幅度，具有较好的工程实践性和实用价值。

图5示出了本发明实施例提供的一种信号重建方法，参见图5，该方法包括如下步骤。

步骤401、对接收到的超宽带信号进行带通滤波。

示例性地，可以通过图1或图3示出的带通滤波器对接收到的超宽带信号进行带通滤波。

步骤402、对带通滤波后的超宽带信号进行欠采样，得到离散时间信号。

下面介绍一下欠采样采样频率的确定方式。假设通过傅里叶变换获得带通滤波后的uwb信号h(t)对应的频域信号x(f)。参见图6中(a)所示，x(f)满足如下条件。

其中fl，fh分别为x(f)的最低频率和最高频率。

按照传统的采样定理进行采样，采样后的频谱是原信号频谱的周期搬移，其周期与采样频率fs相同，即周期为fs。图6中(b)绘出了以fs为采样频率对uwb信号进行采样后，采样的带通uwb信号负频率向右搬移后的频谱。其中，图6中(b)示出了相邻两个周期向右搬移的频谱。为了避免频谱周期性搬移后发生混叠，必须保证：

fs≥2(fh-fl)(2)

(k-1)fs-fl≤fl(3)

以及，

kfs-fh≤fh(4)

其中k为周期序号，k＝1,2,...。由于fl＝fh-b，结合公式(2)-(4)得到：

并且，

结合公式(5)-(6)得到带通采样频率fs的范围为：

其中，可以将公式(7)示出的频率范围用图7表示，其中非阴影区域表示满足限制条件，而阴影区域表示不满足限制条件，图中实线是最小采样频率，可以由公式(8)给出：

公式(8)中m是不大于fh/b的最大整数。可以发现，当fh/b为整数时，最小采样频率为2b。即，欠采样频率fs大于或等于带通滤波器11的带宽b的两倍。当带通滤波器的带宽b等于uwb信号s(t)的新息率ρ时，欠采样频率fs为新息率的两倍。

下面介绍一下旋转不变性技术算法的原理。

参见图8，考虑一个单一的复指数其复幅度是d，时间是t，f0是频率分辨率。这个信号有以下特征：

从式子(9)可以看出，下一个频率采样值是当前频率采样值的相移的结果。这个相移可以表示为：在单位圆上的一个的转动。

基于此，将由冲激串频域信号x[k]组成的矢量信号x[k]分解为由时移频率响应成分s[k]和噪声成分w[k]组成的频率窗口矢量模型，如式子(10)所示。

式子(10)中，矩阵s的p个列：s＝[s(t1)s(t2)…s(tp)]是长为m的频率窗口时间复指数矢量。矢量d由复指数的幅度di组成。矩阵是表示s[k]的单独复指数分量的相邻频率采样值之间相位移动的对角矩阵，如式子(11)所示。

式子(11)中，ti为时移(时间估计)。由于复指数的时间完全描述了这个旋转矩阵，时间估计ti可以通过而得到。需要解决的是，在仅已知冲激串频域信号x[k]的情况下对矩阵进行求解。矩阵的解法原理如下。

旋转不变性技术算法的实质在于其在交错的子空间之间的转动特性，这些子空间被调用以得到时移ti的估计。在离散频域信号的情况下，这个特性依赖于在两个频率上交错的相同区间上的信号观测值。

首先，对未知的时移频率响应矩阵s进行分析，考虑在长为m的频率窗矢量中，两个交叠的、长为m-1的子窗口。参见图8，这两个交叠的、长为m-1的子窗口分别为sm-1[k-1]、sm-1[k]。考虑时移频率响应成分s[k]是由复指数的和组成的信号(参见式子(10))，将时移频率响应成分s[k]表示为式子(12)。

这里，sm-1[k]是s[k]的长为m-1的子窗口，即：

除了它的频率窗口时间矢量的长度是m-1外，矩阵sm-1和矩阵s以相同的方式构造，记为sm-1(t)：sm-1＝[sm-1(t1)sm-1(t2)…sm-1(tp)]。

s[k]是在频率kf0由复指数的和构成的标量信号。定义如下矩阵：

这里，s1对应非交错的窗口，s2对应交错的窗口，即：

这两个窗口频率矢量的矩阵有如下关系：两个矩阵各自生成一个虽然不同、但相互联系的m-1维子空间。因此，可以基于在两个频率上交错的相同区间上的信号观测值s1和s2，求解出矩阵

步骤403、通过傅里叶变换将离散时间信号变换到频域。

假设离散时间信号y[n]为，

y[n]＝<s(t),h(t-nts)>,n＝0,1,...,n-1(15)

式中：h(t)是带通滤波器的冲激响应函数；ts是采样时间间隔；n≥2p。p为脉冲的个数。该离散时间信号y[n]的傅里叶变换为：

设uwb信号s(t)的傅里叶变换为s(f)，将离散时间信号y[n]写成逆傅里叶变换形式：

交换式子(17)积分顺序后：

将公式(18)中[fc-b/2,fc+b/2]段积分写为分段求和后：

式中，为小于等于x的最大整数；带通滤波器h(t)的中心频率fc＝lf0，f0＝fs/n为离散傅里叶变换对应的频率间隔；记由n≥2p可知λ≥p；k＝l-λ,...,l+λ。

将uwb脉冲串信号在频率kf0上的傅里叶变换s(kf0)记为s[k]：

步骤404、对频域信号进行迫零均衡处理。

其中，迫零均衡处理是为了得到冲激串信号x(t)的离散傅里叶变换x[k]。在本实施例中，参见式子(21)，迫零均衡处理方式包括：uwb脉冲信号的频域信号除以高斯脉冲的频域响应，即uwb脉冲信号的傅里叶变换s[k]除以高斯脉冲的频率响应g[k]，得到冲激串信号x(t)的离散傅里叶变换x[k]。

式中g[k]为高斯二阶脉冲信号g(t)(高斯脉冲的二阶导数生成)的离散傅里叶变换。需要说明的是，高斯脉冲的频率响应g[k]是已知信号，不需要求解。

步骤405、基于迫零均衡处理后的频域信号，采用旋转不变性技术算法估算冲激串信号的幅度和时移信息。

示例性地，本步骤405可以包括如下步骤4051-步骤4054。

步骤4051、基于迫零均衡处理后的频域信号，获取设定数量的、且连续的频率采样点。

其中，设定数量的、且连续的频率采样点可以包括，长为m的频率窗口矢量信号x[k]的n个数据记录。

步骤4052、基于设定数量的、且连续的频率采样点，确定超宽带信号的入射角信息。

示例性地，本步骤4052可以包括如下步骤。

第一步、基于长为m的频率窗口矢量信号x[k]的n个数据记录，构造一个汉克尔(hankel)矩阵h。

数据矩阵h的元素包括长为m的频率窗口矢量信号x[k]的n个数据记录。

第二步、对矩阵h进行奇异值分解，获取由矩阵h的p个主右奇异值矢量构成的矩阵v。

h＝uσv^h(23)

其中，u是m×p维矩阵，由h的p个主左奇异值矢量构成，∑是奇异值矩阵，也是uwb信号的入射角矩阵，v^h是l×p维矩阵，由h的p个主右奇异值矢量构成。

第二步、对矩阵v进行信号子空间和噪声子空间的分解。

具体地，参见式子(24)，分解方式包括：对矩阵v进行分块，前s列为信号子空间，后n列为噪声子空间。

v＝[vs|vn](24)

这里，vs是对应p个最大奇异值的右奇异矢量的矩阵。由于信号部分由用频率窗时间矢量s(t)作为模型的复指数的和构成(具体请参见步骤402的推导过程)，所以，矩阵s(式子(10))和vs生成相同的子空间。因此，存在一个可逆变换t将vs映射为s，即：s＝vst。

第三步、将信号子空间分成vs两个较小的m-1维的子空间v1和v2。

这里，与s1和s2类似，v1和v2分别对应非交错的和交错的子空间。由于s1和s2对应相同的子空间，映射关系对这些子空间同样成立：s1＝v1t，s2＝v2t。t为前述可逆变换。

第四步、求解v1和v2之间的关系矩阵ψ。

矩阵s的非交错的和交错的分量s1和s2通过子空间旋转相关联。由于矩阵v1和v2同样分别生成这些相关的子空间，一个相似的联系v1和v2的(旋转)关系肯定存在：v2＝v1ψ。本实施例中，求解矩阵ψ有如下两种方法：ls(最小二乘)和tls(总的最小二乘)。

ls解法中，通过使用ls准则解v2＝v1ψ得到子空间的旋转的估计：ψls＝(v1^hv1)^-1v1^hv2，这个ls解通过在一个最小二乘意义上从公式v2+e2＝v1ψ最小化误差得到，这里，e2是由v2和对应于s2的真实子空间之间的误差构成的矩阵。这个ls公式假定误差仅存在于对s2的估计中，而在v1和欲估计的相应于s1的真实子空间之间没有误差。

tls解法中，由于v1也是一个估计的子空间，一个更合适的公式是v2+e2＝(v1+e1)ψ。这里，e1是表示v1和对应于s1的真实子空间之间的误差的矩阵。这个问题的解称为tls解，是通过最小化两个误差矩阵的frobenius范数：||e1e2||f得到的。

将矩阵v1+e1和v2+e2并排放置后，可以得到方程

将估计的信号子空间矩阵v1和v2并排放置并进行奇异值分解：

其中，和是为m×p阶矩阵，其余为p阶方阵。由最小frobenius范数定理，可知：

其中σi表示矩阵[v1v2]的奇异值，及对角矩阵中的元素。由于两个误差矩阵的最小frobenius范数具有唯一解。交错子空间并排放置矩阵[v1v2]所张成的与脉冲相移相关的信号子空间矩阵为：

此信号子空间矩阵的零空间是的象空间。因此，从公式(25)可知，对某个p阶方阵s，有：

从方程和可看出因此必有：

第五步、求解ψ的特征值，并基于ψ的特征值确定入射角信息。

对矩阵s2：s2＝v2t＝v1ψt，那么由于ψ与有相同的特征值，ψ和的特征值就是矩阵的对角元素{ui}i＝1,...,p。

在解出ψ的特征值后，就能得到入射角信息，入射角信息为∠ui。

步骤4053、基于超宽带信号的入射角信息，确定冲激串信号的时移信息和幅度。

在解出入射角信息后，就能得到时间估计是：同样，再由即可计算出冲激串信号的幅度di。

步骤4054、基于估算出的冲激串信号的幅度和时移信息，计算冲激串信号。

具体地，可以按照式子(10)对冲激串信号x(t)的离散傅里叶变换x[k]进行计算。

通过步骤403-步骤405实现了，基于离散时间信号，采用旋转不变性技术算法估算超宽带信号的幅度和时移信息。

步骤406、基于计算得到的冲激串信号，对接收到的超宽带信号进行重建。

在本实施例中，计算得到的冲激串信号具体是指冲激串信号x(t)的离散傅里叶变换x[k]。可以按照式子(21)对接收到的超宽带信号进行重建，即将计算得到的冲激串信号与高斯脉冲的频率响应g[k]作卷积运算，得到超宽带信号的傅里叶变换s[k]，完成重建。

需要说明的是：上述实施例提供的信号重建装置及系统在重建信号时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的信号重建装置及系统与信号重建方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。

参见图9，其示出了本发明实施例提供的一种信号重建装置的结构示意图。该装置可能是服务器或者终端等计算系统1800，具体来讲：

计算系统1800包括中央处理单元(cpu)1801、包括随机存取存储器(ram)1802和只读存储器(rom)1803的系统存储器1804，以及连接系统存储器1804和中央处理单元1801的系统总线1805。计算系统1800还包括帮助计算机内的各个器件之间传输信息的基本输入/输出系统(i/o系统)1806，和用于存储操作系统1813、应用程序1814和其他程序模块1815的大容量存储设备1807。

基本输入/输出系统1806包括有用于显示信息的显示器1808和用于用户输入信息的诸如鼠标、键盘之类的输入设备1809。其中显示器1808和输入设备1809都通过连接到系统总线1805的输入输出控制器1810连接到中央处理单元1801。基本输入/输出系统1806还可以包括输入输出控制器1810以用于接收和处理来自键盘、鼠标、或电子触控笔等多个其他设备的输入。类似地，输入输出控制器1810还提供输出到显示屏、打印机或其他类型的输出设备。

大容量存储设备1807通过连接到系统总线1805的大容量存储控制器(未示出)连接到中央处理单元1801。大容量存储设备1807及其相关联的计算机可读介质为计算系统1800提供非易失性存储。也就是说，大容量存储设备1807可以包括诸如硬盘或者cd-rom驱动器之类的计算机可读介质(未示出)。

不失一般性，计算机可读介质可以包括计算机存储介质和通信介质。计算机存储介质包括以用于存储诸如计算机可读指令、数据结构、程序模块或其他数据等信息的任何方法或技术实现的易失性和非易失性、可移动和不可移动介质。计算机存储介质包括ram、rom、eprom、eeprom、闪存或其他固态存储其技术，cd-rom、dvd或其他光学存储、磁带盒、磁带、磁盘存储或其他磁性存储设备。当然，本领域技术人员可知计算机存储介质不局限于上述几种。上述的系统存储器1804和大容量存储设备1807可以统称为存储器。

根据本公开的各种实施例，计算系统1800还可以通过诸如因特网等网络连接到网络上的远程计算机运行。也即计算系统1800可以通过连接在系统总线1805上的网络接口单元1811连接到网络1812，或者说，也可以使用网络接口单元1811来连接到其他类型的网络或远程计算机系统(未示出)。

上述存储器还包括一个或者一个以上的程序，一个或者一个以上程序存储于存储器中，被配置由cpu执行。所述一个或者一个以上程序包含用于进行图4和图5其中任一所提供的信号重建方法的指令。

本发明实施例还提供了一种非临时性计算机可读存储介质，当所述存储介质中的指令由计算系统的处理器执行时，使得计算系统能够执行图4和图5其中任一提供的信号重建方法。

一种包含指令的计算机程序产品，当其在计算机上运行时，使得计算机执行进行图4和图5其中任一提供的信号重建方法的指令。

本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。