一种基于压缩感知的LoRa信道估计方法技术

本发明专利技术涉及一种基于压缩感知的LoRa信道估计方法，属于无线传感器网络领域。该方法包括：首先，将LoRa线性调频扩频技术和传输特性与信号因多径衰落、多普勒频移衰减的影响相结合，构建符合LoRa的信道模型；然后，在信道估计时，提出一种基于循环替换导频搜索的算法来设计导频图案，使得LoRa信道模型中的测量矩阵的互相关性最小；最后，发射机通过发送经过导频设计的子载波，当接收机收到子载波后，采用基于自适应步长的稀疏度自适应匹配追踪算法，通过导频子载波重构出信道的状态信息，恢复出LoRa信道的时频域特征值。本发明专利技术能够减小重构过程的复杂度，提高信道估计的整体性能。

A LoRa Channel Estimation Method Based on Compressed Sensing

The invention relates to a LoRa channel estimation method based on compressed sensing, which belongs to the field of wireless sensor networks. The method includes: Firstly, LoRa LFM spread spectrum technology and transmission characteristics are combined with the influence of multipath fading and Doppler frequency shift attenuation to construct LoRa-compliant channel model; secondly, in channel estimation, an algorithm based on cyclic replacement pilot search is proposed to design pilot pattern, which minimizes the cross-correlation of measurement matrix in LoRa channel model; After that, the transmitter transmits the pilot-designed subcarriers. When the receiver receives the subcarriers, the sparseness adaptive matching pursuit algorithm based on the adaptive step size is adopted to reconstruct the channel state information through the pilot subcarriers and recover the time-frequency domain eigenvalues of LoRa channel. The invention can reduce the complexity of the reconstruction process and improve the overall performance of channel estimation.

【技术实现步骤摘要】
一种基于压缩感知的LoRa信道估计方法
本专利技术属于无线传感器网络
，涉及一种基于压缩感知的LoRa信道估计方法。
技术介绍
在无线传感器网络中，由于空间环境的复杂性，发射端的信号在传输的过程中会受到各种障碍物的遮挡、反射、折射等现象，导致接收端接收到的信号是多路信号的叠加，同时信号在传输的过程中能量也在不断的衰减。这就造成了信道的多径衰落，使得整个通信系统的性能受到很大的影响，为了实现对无线信道传输中多径衰落的补偿，就需要对无线信道的信道特征进行准确的估计。随着通信与网络技术的发展，物联网得到了广泛的关注，低功耗广域网(LPWAN)应运而生。其中最为代表的LoRa技术更是得到了广泛的运用，因此对LoRa通信系统的信道估计，准确的获取其信道的特征有着一定的现实意义。在目前信道估计中，首先通过在时域、频域或者空间发送训练序列，再根据接收端接收到的信号通过线性重构的方法来获取信道的特征参数。但是这些数据都是会占据一定的信道资源，降低了频带的利用率。随着压缩感知理论的提出，在无线信道的估计中利用尽量少的训练序列来高精度的估计信道的冲激响应，从而提高系统的频谱利用率。
技术实现思路
有鉴于此，本专利技术的目的在于提供一种基于压缩感知的LoRa信道估计方法，减小基于导频训练序列的LoRa信道估计重构过程的计算复杂度，从而提高信道估计的性能。为达到上述目的，本专利技术提供如下技术方案：一种基于压缩感知的LoRa信道估计方法，具体包括以下步骤：S1：建立LoRa信道模型：根据LoRa特性和调制方式和信号在传输过程中的衰减特性，建立LoRa信道模型；S2：设计导频方案：根据压缩感知原理和建立的LoRa信道模型，在满足模型测量矩阵互相关最小化准则的基础上，得出循环替换导频搜索算法的方案；S3：信道重构：基于自适应步长的稀疏度自适应匹配追踪的重构算法，在接收机接收到经过压缩采样的导频信号后，利用该算法进行信道重构，以当前重构状态(信号能量)来自适应的调整步长，使得重建误差快速收敛，减小重构过程的计算复杂度。进一步，所述步骤S1具体包括以下步骤：S11：LoRa信号在经过线性调制扩频和匹配滤波器后，在某一时刻t，输出的时域信号表示为：其中，A'表示发射信号的幅值，f0为载波中心频率，k表示频率斜率，S＝B2/k表示信号的时间带宽积，B表示信号带宽；S12：在LoRa信号传输时需要配置载波频率f，扩频因子SF，信道带宽Bw和编码率CR,根据不同的参数配置，得到扩频调制后的码片速率r’表示为：S13：根据步骤S12的码片速率，假设一个LoRa数据包的前导码长度为npreamble，有效数据长度为PL，则数据包的传播时间Tpacket表示为：其中，Tpreamble表示前导码传输时间、Tpayload表示前导码传输时间；H＝0表示使用报文头，H＝1表示不使用报文头；DE表示当sx1278芯片中寄存器LowDataRateOptimize位设置为1时，DE＝1，否则DE＝0；S14：假设发送机到接收机的传输距离为D，则LoRa信号在发送过程中每一单位路径d消耗的时间为考虑多径效应与多普勒效应的影响，得到单位长度下的传播时间表示为：t'＝td+tj+αltp其中，tj表示多径时延，tp表示多普勒效应下的时延，αl表示多普勒调频率；则LoRa信道的冲激响应表示为：其中δ(·)表示阶跃函数，hl(t)表示路径的幅值；S15：假设通信中有L条路径，多普勒调频率不变，上行信道的集合Hup＝[h1,h2,…,hK1]，K1表示上行信道数，得到上行传输中的接收信号表示为：假设通信中有L条路径，多普勒调频率不变，下行信道的集合Hdown＝[h1，h1，…，hK2]，K2表示下行信道数，得到下行传输中的接收信号表示为：其中，n(t)表示高斯白噪声；考虑所有子载波，定义接收信号向量Y，发送信号向量为X，噪声向量为n，则系统的输入输出关系为：Y＝XH+n。进一步，根据压缩感知原理在基于导频序列的信道估计中，在获得导频信号的输入输出关系后，需要得到求解信道状态参数的函数，其步骤如下：S161：假设信号的载波数为N，需要的导频个数为P，则假设导频集合为Λ＝{k1，k2，…，kp}，它是所有子载波集合Ω＝{1，2，...，N}的子集，导频发送的信号表示为：[X(k1)X(k2)…X(kP)]T，导频接收信号表示为：[Y(k1)Y(k2)…Y(kP)]T，噪声信号表示为：[n(1)n(2)…n(P)]T，则发送和接收信号之间用矩阵表示为：其中，FP×L为离散傅里叶矩阵，其表示为：其中，ω＝e-j2π/N；则YP＝XPFP×Lh+nP＝Ah+nP，A＝XPFP×L表示测量矩阵；S162：由最小均方差估计，当均方差最小时，得到信道状态信息矩阵h的表达式为：h＝(AHA)-1AHy。进一步，所述步骤S2具体包括：测量矩阵A要保证稀疏的信道状态矩阵h能从观测值y中恢复出来，需要测量矩阵满足互相关最小化准则，得出近视最优的导频方案算法，其步骤如下：S21：根据不同的导频图案，将测量矩阵的互相关性表达式为：则最后输出的导频模式pout表示为：pout＝argming(p)S22：初始化：循环次数T，导频数目P，随机生成M种导频图案则对应的测量矩阵互相关性为S23：并行替换，在第s(s≤kp)次并行替换时，分别将每种导频图案的第s-1次替换后的图案中的第s个元素进行替换，第一次替代是针对初始备选节点；即对于导频序列为的备选图案固定导频位置将所有子载波的集合Ω与的差集中的元素置于ps-1，然后根据互相关性表达式计算导频图案的互相关性，共计算N-kp次；从计算得到的组合中，再选择M个互相关性最小的导频图案作为第s次替换的备选图案，供第s+1次替换使用，记做S24：循环选优，经过kp次并行替换后，得到M个备选导频方案m＝1，2，…M，如果此时循环次数小于T，则继续步骤S22，否则结束循环，得到一个相关性最小的导频方案作为最终结果pout。进一步，所述步骤S3具体包括：当接收信号是在近似最优导频方案下的信号，通过导频信号重构信道状态时，采用自适应步长的稀疏度自适应匹配追踪的重构算法进行信道重构，其具体包括以下步骤：输入：设置初始的接收导频信号，其中信号的子载波为YP，测量矩阵A，残差阈值ε，附加阈值Γ，初始步长sI，存放测量矩阵A与残差r内积值的集合S31：初始化h＝[0，0，K，0]T，hold＝[0，0，K，0]T，rtemp＝[0，0，K，0]T，索引集候选集残差r0＝Yp，最终支持集的大小LF＝s＝sI，最终支持集迭代索引i＝1；S32：在信号参数经过第i-1次重构之后，比较信道参数的实际残差与设定残差阈值是否接近，即||ri-1||≥ε，如果满足，则直接结束，得到所要估计的信道参数h，不满足则继续执行以下步骤；S33：计算测量矩阵各列与残差直接的值，存在SP中，SP＝||AHri-1||；S34：在步骤S33中计算得到的SP中，选择出前LF个最大的值，并且将这些值对应A的列序号存入到索引结合Bi中；S35：将步骤S34中得到的索引集合和最终支持集合合并到候选支持集中，得到Ci＝Bi∪Fi-1：S36：计算出的值，并在其中选择出前LF个最大元素的索引，得到最终支持集F本文档来自技高网...

【技术保护点】
1.一种基于压缩感知的LoRa信道估计方法，其特征在于，该方法具体包括以下步骤：S1：建立LoRa信道模型：根据LoRa特性和调制方式和信号在传输过程中的衰减特性，建立LoRa信道模型；S2：设计导频方案：根据压缩感知原理和建立的LoRa信道模型，在满足模型测量矩阵互相关最小化准则的基础上，得出循环替换导频搜索算法的方案；S3：信道重构：基于自适应步长的稀疏度自适应匹配追踪的重构算法，在接收机接收到经过压缩采样的导频信号后，利用该算法进行信道重构，以当前重构状态来自适应的调整步长，使得重建误差快速收敛，减小重构过程的计算复杂度。

【技术特征摘要】
1.一种基于压缩感知的LoRa信道估计方法，其特征在于，该方法具体包括以下步骤：S1：建立LoRa信道模型：根据LoRa特性和调制方式和信号在传输过程中的衰减特性，建立LoRa信道模型；S2：设计导频方案：根据压缩感知原理和建立的LoRa信道模型，在满足模型测量矩阵互相关最小化准则的基础上，得出循环替换导频搜索算法的方案；S3：信道重构：基于自适应步长的稀疏度自适应匹配追踪的重构算法，在接收机接收到经过压缩采样的导频信号后，利用该算法进行信道重构，以当前重构状态来自适应的调整步长，使得重建误差快速收敛，减小重构过程的计算复杂度。2.根据权利要求1所述的一种基于压缩感知的LoRa信道估计方法，其特征在于，所述步骤S1具体包括以下步骤：S11：LoRa信号在经过线性调制扩频和匹配滤波器后，在某一时刻t，输出的时域信号表示为：其中，A'表示发射信号的幅值，f0为载波中心频率，k表示频率斜率，S＝B2/k表示信号的时间带宽积，B表示信号带宽；S12：在LoRa信号传输时需要配置载波频率f，扩频因子SF，信道带宽Bw和编码率CR,根据不同的参数配置，得到扩频调制后的码片速率r’表示为：S13：根据步骤S12的码片速率，假设一个LoRa数据包的前导码长度为npreamble，有效数据长度为PL，则数据包的传播时间Tpacket表示为：其中，Tpreamble表示前导码传输时间，Tpayload表示数据负载传输时间；H＝0表示使用报文头，H＝1表示不使用报文头；DE表示当芯片中寄存器LowDataRateOptimize位设置为1时，DE＝1，否则DE＝0；S14：假设发送机到接收机的传输距离为D，则LoRa信号在发送过程中每一单位路径d消耗的时间为考虑多径效应与多普勒效应的影响，得到单位长度下的传播时间表示为：t'＝td+tj+αltp其中，tj表示多径时延，tp表示多普勒效应下的时延，αl表示多普勒调频率；则LoRa信道的冲激响应表示为：其中δ(·)表示阶跃函数，hl(t)表示路径的幅值；S15：假设通信中有L条路径，多普勒调频率不变，上行信道的集合Hup＝[h1,h2,…,hK1]，K1表示上行信道数，得到上行传输中的接收信号表示为：假设通信中有L条路径，多普勒调频率不变，下行信道的集合Hdown＝[h1,h1,…,hK2]，K2表示下行信道数，得到下行传输中的接收信号表示为：其中，n(t)表示高斯白噪声；考虑所有子载波，定义接收信号向量Y，发送信号向量为X，噪声向量为n，则系统的输入输出关系为：Y＝XH+n。3.根据权利要求2所述的一种基于压缩感知的LoRa信道估计方法，其特征在于，根据压缩感知原理在基于导频序列的信道估计中，在获得导频信号的输入输出关系后，需要得到求解信道状态参数的函数，其步骤如下：S161：假设信号的载波数为N，需要的导频个数为P，则假设导频集合为Λ＝{k1,k2,…,kP}，它是所有子载波集合Ω＝{1,2,...,N}的子集，导频发送的信号表示为：[X(k1)X(k2)…X(kP)]T，导频接收信号表示为：[Y(k1)Y(k2)…Y(kP)]T，噪声信号表示为：[n(1)n(2)…n(P)]T，则发送和接收信号之间用矩阵表示为：其中，FP×L为离散傅里叶矩阵，其表示为：其中，ω＝e-j2π/N；则YP＝XPFP×Lh+nP＝Ah+nP，A＝XPFP×L表示测量矩阵；S162：由最小均方差估计，当均方差最小时，得到信道状态信息矩阵h的表达式为：h＝(AHA)-1AHy。4.根据权利要求3所述的一种基于压缩感知的LoRa信道估计方法，其特征在于，所述步骤S2具体包括：测量矩阵A要保证稀疏的信道状态矩阵h能从观测值y中恢复出来，需要测量矩阵满足...

【专利技术属性】
技术研发人员：谢昊飞，丁凡，王平，黄代维，
申请(专利权)人：重庆邮电大学，
类型：发明
国别省市：重庆,50