具体实施方式

为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。

需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

本发明提供的超分辨多径条件下时延估计方法的核心是在源信号和多径数量都未知情况下通过重建潜在信号模型对多径数量及时延值的联合估计的方法，使得重建的潜在信号与源信号相匹配，匹配的相似程度越大说明重建信号越准确，相应的多径时延值就越精确。如何准确的重建潜在信号求解多径的时延值，要完成这一核心目标，需要解决问题如下：

(1)信号模型是多径时延准确估计的前提和基础，复杂的无线环境下，信号系统模型精确建立的难度较高；

(2)无线通信中，数字调制信号波形对多径时延的超分辨显示会产生较大程度的影响；

(3)无线传播多径数量未知的情况下，获取准确的多径时延估计结果的难度系数较高。

本发明通过以下技术方案解决上述技术问题：

本发明提供的超分辨多径数量及时延联合估计方法是一种源信号与多径数量未知的时延估计方法，通过消除源信号波形和干扰的影响，提出了一种新的互谱函数，利用互谱函数的超分辨显示提取时延参数信息，根据时延参数信息重建潜在信号与互谱函数相匹配。

如图1所示，本发明提供了一种超分辨多径数量及时延联合估计方法，包括如下步骤：

S1、建立通用的多径接收信号离散时间系统模型；

S2、利用信号的频域变换，构造互谱函数以消除未知源信号波形的影响；

S3、对构造的互谱函数波形的超分辨显示进行谱峰搜索；

S4、根据互谱函数局部峰值确定多径相对时延及对应幅度值，重建潜在信号与互谱函数进行匹配对应；

S5、构造代价函数，通过最大似然估计对多径数量及相应的时延值进行联合估计，得出时间延迟值。

具体实施时，作为本发明优选的实施方式，信号在传播的过程中由于道路中的建筑物或移动物体遮挡等原因会造成信号有多条路径传播，如图2所示，其中发送端到接收端都有一条主径，也就是直达路径，但由于反射的影响，信号在发送端发出传播的过程中分别产生了2条和3条多径数，信号会由于多径的影响造成接收信号产生偏差，为了能准确地提取出主径的源信号，在接收端由2个接收站对发送的数字调制信号进行时间同步采集，多个接收站则类似，以此类推。则所述步骤S1的具体实现过程如下：

S11、将接收到的离散时域信号进行建模，如下：

其中，n＝1,2,...,K，n为离散采样点，K表示总采样点数，是第一路源信号到接收站之间的信道衰减系数，是第二路源信号到接收站之间的信道衰减系数，J₁为第一路接收信号的多径数量，J₂为第二路接收信号的多径数量，是第一路接收信号的每一条多径数相对应的时延值，是第二路接收信号的每一条多径数相对应的时延值，ω₁(t)是第一路接收的加性零均值噪声，ω₂(t)是第二路接收的加性零均值噪声；

S12、对离散时域信号作频域变换，表达式如下：

其中，k＝1,2,...,K，R₁(k)是r₁(n)的离散傅立叶变换，R₂(k)是r₂(n)的离散傅立叶变换，S(k)是s(n)的离散傅立叶变换，W₁(k)是ω₁(n)的离散傅立叶变换，W₂(k)是ω₂(n)的离散傅立叶变换。

具体实施时，作为本发明优选的实施方式，由步骤S12中的表达式得到的离散傅里叶变换可知，接收信号的多径时延信息存在于频率差中，由于源信号是未知的，受波形的影响，在频域范围内无法计算出时延，要消除信号波形产生的影响，需要将信号本身的波形幅值化为1以降低偏差。为了能求解出时延信息，所述步骤S2的具体实现过程如下：

S21、设置一组时间延迟点D^l，l＝1,2,...,N，其中N是等间隔采样的总点数，对步骤S12中的表达式作归一化处理，得到互谱函数如下：

其中，H表示共轭；

S22、对傅里叶的采样平均值进行离散来抑制噪声的部分，令将互谱函数重写为：

S23、基于步骤S22的互谱函数，设置的时间延迟点D^l最接近真实的时间延迟点时会有较大的谱峰，即当时间延迟点D^l＝τ_i-λ_j时，φ^l的幅值最大，即：

S24、若设置的时间延时点D^l不接近真实的时延值，则表示频率差依然存在，复正弦波形在几个完整周期的平均值为零。

具体实施时，作为本发明优选的实施方式，接收信号中由于存在随机噪声的干扰，所以会导致互谱函数中出现多个局部的最大值φ^l，因此很难确定哪些局部最大值是所需要的多径时延值，因为不知道接收信号中到底是由多少条多径进行传播的，因此提出了一种谱峰搜索的方法，本实施例中，人为设置M个峰值，但此处没有区分真实路径和干扰信息，所述步骤S3的具体实现过程如下：

S31、令m＝1,2,...,M，搜索|φ^l|的相应M峰的信息，利用函数得到相应的时延信息和幅度信息

S32、定义间隔宽度Δl，并将|φ^l|的值在[l^*-Δl,l^*+Δl]范围内设置为零；

S33、若m＝M，则表示多径数与所设值相匹配，否则m＝m+1，返回执行步骤S31。

具体实施时，作为本发明优选的实施方式，所述步骤S4的具体实现过程如下：

在步骤S3中可知，互谱函数由复正弦信号组成，通过加权和组合，由于每个接收站的路径数都是未知的，所以需要恢复一组具有普适性的信号来处理所有的可能性，在每个假设中，m时从1,2,...,M间延迟数得到的参数，如下：

由于已知幅度和时延信息，则将重构的信号表示为：

具体实施时，作为本发明优选的实施方式，本发明提供的超分辨多径时延估计方法，有机地创新性地提出最大似然与多径条件相结合，使目标函数与所构造的信号之间有尽可能大的相似性，所述步骤S5的具体实现过程如下：

S51、定义代价函数，如下：

S52、由上述代价函数得出，代价函数输出越小，代表所重构的信号与互谱函数有极大的相似性，即重构的信号接近于含有真实时延的谱，所以通过使代价函数输出最小化来搜索每个数的最优值，表达式如下：

S53、遍历每个代价函数的可能取值，得到最小的值，即得到时延的估计，如下：

实施例

为了验证本发明方法的可行性，本发明测试了2个接收站接收的多径数字调制信号，取源信号为BPSK信号，BPSK信号在持续的码元周期内保持不变的状态，因此对于超分辨的情况，也就是时延长度小于码元持续长度时，由于信号在码元周期内状态保持不变，而时延长度较小时，对持续长度的码元影响较小，不易得出时延的准确估计结果。

源信号取BPSK信号也具有一定的代表性，其采样频率为f_s＝2×10⁷Hz，采样间隔Ts＝1/fs，载波频率为f_c＝0.2f_s，定义码速率为R＝1×10⁶bps，码元内采样点数为f_s/R＝20，本发明分别对两个接收站接收信号的数量进行调整，考虑了四种不同的多径数量情况，分别是：

当第一路接收信号的多径数量为1条，第二路接收信号的多径数量也为1条；

当第一路接收信号的多径数量为1条，第二路接收信号的多径数量为2条；

当第一路接收信号的多径数量为1条，第二路接收信号的多径数量为3条；

当第一路接收信号的多径数量为2条，第二路接收信号的多径数量也为2条。使得多径数目最终为1、2、3、4这4种情况；

第一种情况：两路信号的时间延迟差和衰减系数分别是τ₁₁＝0，α₁₁＝1和τ₂₁＝-148，α₂₁＝1；

第二种情况：是在第一种情况下在第二路增加了τ₂₂＝-153，α₂₂＝0.8的第2路多径信号；

第三种情况：是在第二种情况下在第二路继续增加了τ₂₃＝-157，α₂₃＝0.7的第3路多径信号；

第四种情况：是在第二种情况下在第一路继续增加了τ₁₂＝-9，α₁₂＝0.9的多径信号。

在多径数量为2及以上情况下，其相对时延值Δτ₂₁＝|τ₂₂-τ₂₁|＝5＜20；Δτ₃₁＝|τ₂₃-τ₂₂|＝4＜20；Δτ₃₂＝|τ₂₃-τ₂₁|＝9＜20；Δτ₄₁＝|τ₁₂-τ₁₁|＝9＜20，时延值小于码元持续时间的为超分辨情况。

如图3所示，为这4种情况下的互谱函数的超分辨显示，其时延真值已在图中标出，由图可知，其峰值与时延真值基本吻合，误差较小。可见本发明方法是可行的。

为验证本发明的广泛性，考虑到低信噪比情况下时延估计的准确性判别，在信噪比为5dB时，分别做了以上的多径数量的相同操作，结果如图4所示，可见在低信噪比情况下本发明方法同样适用。在互谱函数的超分辨显示中可以清晰地看到，除了有想要得到的多径谱峰外，还会在其他地方产生局部的小峰值，由于多径数是未知的，仅仅根据峰值情况来判断多径数量是行不通的，本发明借助构造的潜在信号，令M＝6，构造了6组潜在信号，分别对应多径数量从1到6的6种不同的情况，由构造的信号分别从多径数量为1开始与互谱函数匹配对应。本发明利用两信号相减做比较，两信号的匹配度高，也就是相似程度越高，那么相应的得到相减后的结果就越小，最终遍历完全部的构造信号，最小取值的结果就为实际的多径数量，如图5所示，四条曲线分别对应多径数为1，2，3，4的情况，可以很清晰的看出在对应的多径数量时有最小的取值，其结果与预先设定的多径数量吻合。相应的，就可以得到多径所对应的时延值。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。