CN113922850B - 一种LoRa上行***中分布式MU-MIMO信道估计方法及传输方法 - Google Patents
一种LoRa上行***中分布式MU-MIMO信道估计方法及传输方法 Download PDFInfo
- Publication number
- CN113922850B CN113922850B CN202111042442.9A CN202111042442A CN113922850B CN 113922850 B CN113922850 B CN 113922850B CN 202111042442 A CN202111042442 A CN 202111042442A CN 113922850 B CN113922850 B CN 113922850B
- Authority
- CN
- China
- Prior art keywords
- nodes
- node
- transmission
- signals
- sampling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 46
- QVFWZNCVPCJQOP-UHFFFAOYSA-N chloralodol Chemical compound CC(O)(C)CC(C)OC(O)C(Cl)(Cl)Cl QVFWZNCVPCJQOP-UHFFFAOYSA-N 0.000 title claims abstract description 42
- 239000011159 matrix material Substances 0.000 claims abstract description 66
- 238000005070 sampling Methods 0.000 claims abstract description 65
- 238000001514 detection method Methods 0.000 claims description 11
- 108010003272 Hyaluronate lyase Proteins 0.000 claims description 4
- 238000005562 fading Methods 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0452—Multi-user MIMO systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Radio Transmission System (AREA)
Abstract
本发明提供了一种LoRa上行***中分布式MU‑MIMO信道估计方法及传输方法,其中:对于K个节点传输的信号y同时到达接收端时,接收端首先在接收到的信号中检测前导码,在检测到前导码后,继续检测不同节点的符号边缘是否对齐。当符号边缘对齐时,计算预存的哈达玛矩阵阶数,使信号y中chiirp符号数量,以采样率对信号y进行采样,遍历哈达玛序列与各信号相乘求解h ji。当符号边缘不对齐时,以采样率开始采样,建立方程求解h ji。通过特殊采样率降低了传输方程中未知数的个数,并基于哈达玛序列提出一种求解信道信息的方法,在获取信道信息后采用信道均衡算法分离混合信号,从而解决了在同频异步传输的LoRa网络中实现MU‑MIMO传输的问题。
Description
技术领域
本发明属于无线通信技术领域,涉及LoRa上行***,具体涉及一种LoRa上行***中分布式MU-MIMO信道估计方法。
背景技术
为了提升LoRa***整体吞吐量、充分利用频带资源,但由于商用网关的限制,一个解调器只能解调一路数据,在同一频率上只能使用一个SF进行传输,这极大浪费了频带资源,限制了传输吞吐量和传输效率。现有方法通过联合时域特征来实现并发解码,将***吞吐量提升为LoRaWAN的3倍;也有方法利用累积频谱系数区分没有时间偏移的冲突信号,但需要冲突信号具有不同的功率。虽然高速率并非LoRa的设计初衷,但较低的吞吐量将会限制***传输效率的提升。通过采用MU-MIMO传输可以提升***吞吐量,提升传输效率。但LoRa节点的分布式布置和异步传输的特点是LoRa实现MU-MIMO的一大挑战。
由于LoRa采用CSS调制,采用不同SF调制的信号彼此正交,因此不同节点可以在同一频率上使用不同SF进行传输。但是在有数据率约束的场景下,如果只有较小的SF才能满足速率需求,而且如果使用同一频率的节点都需要选择同一SF才能满足速率需求时,在标准LoRaWAN下无法正常进行传输,此时LoRa***的吞吐量将受到限制。
为了将MU-MIMO技术应用在LoRa中,需要在同一频率上设置了相同SF的LoRa节点可靠地进行传输。LoRa由于其分布式的特性,各个节点的传输不具有同步性,因此无法直接采用MU-MIMO进行传输,而现有方法中提出的盲分布式MU-MIMO传输方案可以实现在不具同步性的分布式节点上启用MU-MIMO传输,该方案虽然可以直接在LoRa***上使用,但在一些特定情况下存在限制,例如在具有K个单天线发送端和一个具有R根天线接收端的MU-MIMO***中,信号同时到达网关时,即使前导码数量足够多,但是节点传输的前导码对应相等,建立的传输方程组对应的向量组线性相关,必然存在多余的方程,因此无法对信道信息进行求解。此外,还要求接收端天线至少比发送端的天线总数多1根,这对接收端天线数量也进行了限制。
发明内容
针对现有技术存在的不足,本发明的目的在于,提供一种LoRa上行***中分布式MU-MIMO信道估计方法,以解决现有技术中尚没有在同频异步传输的LoRa网络中实现MU-MIMO传输的技术问题。
为了解决上述技术问题,本发明采用如下技术方案予以实现:
一种LoRa上行***中分布式MU-MIMO信道估计方法,包括以下步骤:
步骤1,对于K个节点传输的信号y同时到达接收端时,由于K个节点都是新加入***的节点,因此有个信道/>未知;传输1个符号会产生/>个未知数和R个方程,因此采用的编码矩阵至少为K+1阶;
其中,R个方程中的第一个方程中,第一个接收天线接收到的数据为:
式12;
R个方程中的其它方程中,接收天线接收到的数据的表达式在式12的基础上依次类推;
j表示第j个接收天线;
i表示第i个节点;
表示第i个节点到第j个接收天线的信道;
s表示chirp;
表示第/>个采样点;
步骤2,接收端首先在接收到的信号中检测前导码,如果检测到前导码,则进行步骤3,如果未检测到前导码会持续检测,直到检测到前导码;
步骤3,在检测到前导码后,继续检测不同节点的符号边缘是否对齐;
当符号边缘对齐时,执行步骤4,当符号边缘不对齐时,执行步骤5;
步骤4,当符号边缘对齐时,执行步骤4.1至步骤4.4;
步骤4.1,首先计算预存的哈达玛矩阵阶数,然后接收等于哈达玛矩阵阶数/>的采样点;
步骤4.2,接收端预存阶数为K+1的哈达玛矩阵,以采样率对信号y进行采样;
其中,B表示信号带宽;SF表示扩频因子;
步骤4.3,每接收K+1个采样点就遍历哈达玛矩阵中的行来与相乘,遍历到了第一个节点所使用的序列,由于哈达玛矩阵与编码矩阵不同行相乘求和为0,该序列与、/>、…、/>相乘求和均为0,因此,只剩下一项/>,除以即可求出信道信息;
重复该过程,当完整遍历完一遍哈达玛矩阵后即可求得全部信道信息;
步骤4.4,对于K个节点中有个节点新加入***,/>个节点稳定传输的情况,/>个稳定传输的节点传输过程中的信道/>已知,并且数据部分以采样率/>采样后转换为已知符号,再利用哈达玛矩阵通过控制节点传输功率编码为已知序列,又由于/>个新加入的节点传输的符号已知,因此只需估计/>个信道/>,由于节点总数为K,因此依然需要编码矩阵至少为K+1阶;
其中,R个方程中的第一个方程中,第一个接收天线接收到的数据为:
式13;
R个方程中的其它方程中,接收天线接收到的数据的表达式在式13的基础上依次类推;
对于个稳定传输的节点,其信道已知但序列并未与新加入节点的序列对齐,由于哈达玛矩阵与编码矩阵不同行错位相乘求和结果为0,因此依然可以区分不同节点传输的信号;对于/>个新加入的节点,利用节点对应的序列与/>相乘求和,最后只剩下一项,而该项即为该节点对应的信道信息与/>的乘积;重复该过程能够求出所有未知信道;
接收天线数大于等于发射天线数,即需满足;
步骤5,当符号边缘不对齐时,执行步骤5.1至步骤5.2;
步骤5.1,以采样率开始采样;
步骤5.2,在1个节点已稳定传输,另一个节点新加入***,接收端使用2根接收天线的情况下,建立如式14所示的传输方程,通过该传输方程可以求解出所有未知量;
式14;
式中:
和/>分别对应第一个接收天线两个时刻的采样点,这两个采样点分别对应两路信号各自的采样点;
和/>则为另一个接收天线上的两个采样点;
和/>由于是已加入***的节点的信道,因此为已知量;
、/>、/>、/>均为未知量,共计4个;
表示表示节点传输的数据;
步骤6,在获取到全部参数后,信道矩阵H可以直接获得。
本发明还保护一种LoRa上行***中多节点MU-MIMO传输方法,包括以下步骤:
步骤一,利用多个节点给网关发送信号,不同节点的信号经过瑞利衰落信道进行传输,所有信号与噪声组合之后到达接收天线;
步骤二,对步骤一得到的信号,首先需要检测前导码;
步骤三,接收端进行分布式MU-MIMO信道估计,获取信道信息;
所述的分布式MU-MIMO信道估计的方法采用如上所述的LoRa上行***中分布式MU-MIMO信道估计方法;
步骤四,将步骤三得到的信道信息,再将信道矩阵H代入信道均衡算法中即可解出MU-MIMO***中所有节点传输的信号,最后对各路信号分别使用LoRa解码算法,求解出各个节点传输的数据,从而完成了一个完整的MU-MIMO传输过程。
本发明与现有技术相比,具有如下技术效果:
(Ⅰ)本发明的算法实现了在LoRa***中启用MU-MIMO传输的可能性,为LoRa启用MU-MIMO提出了一种可用的解决方案。
(Ⅱ)本发明的算法通过特殊采样率降低了传输方程中未知数的个数,并基于哈达玛序列提出一种求解信道信息的方法,在获取信道信息后采用信道均衡算法分离混合信号,从而解决了在同频异步传输的LoRa网络中实现MU-MIMO传输的问题。
(Ⅲ)在LoRa***中采用MU-MIMO时,信号之间的干扰远远低于未采用MU-MIMO的***,这极大降低了节点间由于传输功率导致的传输错误,提高了***的鲁棒性。
附图说明
图1是LoRa上行***中分布式MU-MIMO信道估计的流程示意图。
图2是LoRa上行***中多节点MU-MIMO传输方法的流程示意图。
图3是信号同时到达接收端及采样结果图。
图4是新节点异步到达接收端及采样结果图。
图5是相同SF并发传输时信号检测算法比较图。
图6是不同SF并发传输时误比特率随SF变化图。
图7是相同SF并发传输时误比特率随SF变化图。
图8是不同SF并发传输时单个节点吞吐量随SNR变化图。
图9是相同SF并发传输时单个节点吞吐量随SNR变化图。
以下结合实施例对本发明的具体内容作进一步详细解释说明。
具体实施方式
需要说明的是,本发明中:
LoRa,全称为Long Range Radio,即远距离无线电。
CSS,全称为Chirp Spread Spectrum,即线性调频扩频调制。
MU-MIMO,全称为Multi-User Multiple-Input Multiple-Output,即多用户-多输入多输出。
SVD,全称为Singular Value Decomposition,即奇异值分解。
SNR,全称为Signal to Noise Ratio,即信噪比。
MMSE,全称为Minimum Mean-Squared Error,即最小均方误差。
MRC,全称为Maximum-Ratio Combining,即最大比合并。
ZF,全称为Zero-Forcing,即迫零算法。
为了解决背景技术中的这些问题,本发明提出一种新颖的MU-MIMO信道估计方法,该方法通过一种特殊的采样率降低了MU-MIMO传输过程中产生的未知数个数,利用正交编码的正交特性提出一种求解信道信息的方法,从而为接收端使用信道均衡算法分离各路信号时提供信道信息,从而实现 MU-MIMO 传输。
基于信道均衡算法的MIMO需要在接收端获取信道信息,而在实际***中,遵循ALOHA协议的LoRa***并不能保证节点有序传输,在异步传输的***中,由于不同节点传输的前导码部分和数据部分很可能叠加在一起,因此很难获取信道信息。这也是异步传输的LoRa***中实现MIMO的关键问题。
为了解决上述问题,提出一种信道估计方法,用于在异步传输状态下获取LoRa节点的信道,从而为后续利用信道均衡算法分离混合信号提供信道信息。
为了实现分布式MU-MIMO传输,现有方法中通过增加数据样本联合多个传输方程求解未知信息,其问题在于当其它节点加入MU-MIMO***时,已经处于稳定状态的节点传输的数据未知,这会增加所需求解的未知数个数。当大量数据包快速到达时,所需求解的未知数数量过多而方程有限,可能无法完成解码任务。
为了降低异步传输时的未知数的个数,通过采用不同于正常LoRa解码时所采用的采样间隔,而使用采样间隔为/>对传输的信号进行采样,以一个具有任意循环移位C的up-chirp为例,其表达式为:
式1;
式中,表示时间,/>为频率,表示在一个chirp周期内,一个up-chirp的频率随时间的增加而线性增加。up-chirp 的频率可以表示为:
式2;
式2中,C为循环移位值,可表示为:
式3;
为一个信息向量,表示一个chirp传输的SF个比特,也就是所需要传输数据的二进制形式,式3中/>为/>中的一个元素。
式2中,为一个chirp的周期,可表示为:
式4;
一个标准的up-chirp,其频率由线性增长至,其中/>为中心频率,B为信号带宽。
对式1以采样间隔为进行采样得:
式5;
式中,SF表示扩频因子,B表示信号带宽,表示采样点,/>表示频率,/>表示一个chirp的周期,C表示循环移位值。
根据式5可以发现,通过采样间隔采样消除了循环移位C,而LoRa前导码作为一类特殊的符号,仅将循环移位C置为0,对于循环移位C不为0的载荷信号部分,使用上述采样间隔采样后与前导码采样后得到的是相同的结果,因此这种方法可以视为延长了已知序列的长度。
对于LoRa中的前导码部分和数据部分而言,通过使用采样间隔采样消除了具有不同循环移位的信号之间的差异,从而将前导码部分和数据部分归并为了一类相同的数据,这有效减少了MU-MIMO传输方程中未知数的个数,减少了参数估计时所需样本数。
在将前导码和数据归为一类数据之后,对于信号同时到达网关或一路信号正在传输数据,另一路信号新加入并传输前导码的情况都可以采用上述方式归类为一种MU-MIMO传输形式。以K发R收的MU-MIMO为例,K个节点放置在任意位置,且在任意时刻进行传输,网关上的R根接收天线接收到信号后以采样间隔进行采样。
经过以为采样间隔采样的R根接收天线上接收到的信号在不考虑噪声的情况下可以简化表示为:
式6;
式中,表示第/>个采样点,/>表示节点i到接收天线j的信道,共有K个节点和R个接收天线,/>为采样后的chirp,其值为:
式7;
式中,表示第/>个采样点,式7在式5的基础上增加了一项功率控制系数/>,表示一个功率随时间变化的序列,该序列采用哈达玛(Hadamard)矩阵中的行序列,即哈达玛序列,以保证不同节点传输信号的正交性。一个/>阶哈达玛矩阵如下所示:
式8;
该矩阵可通过2阶哈达玛矩阵逐次求得,2阶哈达玛矩阵为:
式9;
由式8和式9可获得一个4阶哈达玛矩阵如式10所示,在该矩阵中,不同行序列具有正交性,即不同行序列对应相乘求和结果为0,同一序列对应相乘求和结果为矩阵阶数,矩阵中不同行序列错位相乘求和的结果也为0。
式10;
在实际情况下,由于在节点处信号幅值只能通过节点的发射功率来控制,因此所使用的编码矩阵中的元素均应为正值。为了防止载荷部分功率过高导致前导码检测的错误,需要控制编码矩阵中的元素取值范围在(0,1]区间内。所以采用的编码矩阵如式11所示:
式11;
该矩阵通过将哈达玛矩阵中的-1替换为0.5获得。使用编码矩阵编码的数据流对应使用与编码矩阵同阶的哈达玛矩阵解码。但需注意的是,编码矩阵中的第一行,即全为1的行不可用于编码,因为编码矩阵中的元素全为正值。解码时如果依然使用全为1的行与接收信号相乘,并不能消去其它干扰项,从而造成解码错误,因此,一个Q阶编码矩阵中只有Q-1行序列可用于编码。此外,编码矩阵与哈达玛矩阵对应的行错位奇数位时,相乘求和结果的符号会变化,需要乘以-1处理。
遵从上述技术方案,以下给出本发明的具体实施例,需要说明的是本发明并不局限于以下具体实施例,凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。
实施例1:
本实施例给出一种LoRa上行***中分布式MU-MIMO信道估计方法,如图1所示,包括以下步骤:
步骤1,对于K个节点传输的信号y同时到达接收端时,由于K个节点都是新加入***的节点,因此有个信道/>未知;传输1个符号会产生/>个未知数和R个方程,因此采用的编码矩阵至少为K+1阶;
其中,R个方程中的第一个方程中,第一个接收天线接收到的数据为:
式12;
R个方程中的其它方程中,接收天线接收到的数据的表达式在式12的基础上依次类推;
j表示第j个接收天线;
i表示第i个节点;
表示第i个节点到第j个接收天线的信道;
s表示chirp;
表示第/>个采样点;
本实施例中,两个节点传输的信号如图3所示,图中信号的前8个符号为前导码部分,后4个符号为载荷部分,图中红色圆圈标出的是采样后的序列。
步骤2,接收端首先在接收到的信号中检测前导码,如果检测到前导码,则进行步骤3,如果未检测到前导码会持续检测,直到检测到前导码;
步骤3,在检测到前导码后,继续检测不同节点的符号边缘是否对齐;
当符号边缘对齐时,执行步骤4,当符号边缘不对齐时,执行步骤5;
步骤4,当符号边缘对齐时,执行步骤4.1至步骤4.4;
步骤4.1,首先计算预存的哈达玛矩阵阶数,然后接收等于哈达玛矩阵阶数/>的采样点;
步骤4.2,接收端预存阶数为K+1的哈达玛矩阵,以采样率对信号y进行采样;
其中,B表示信号带宽;SF表示扩频因子;
步骤4.3,每接收K+1个采样点就遍历哈达玛矩阵中的行来与相乘,遍历到了第一个节点所使用的序列,由于哈达玛矩阵与编码矩阵不同行相乘求和为0,该序列与、/>、…、/>相乘求和均为0,因此,只剩下一项/>,除以即可求出信道信息;
重复该过程,当完整遍历完一遍哈达玛矩阵后即可求得全部信道信息;
由于前导码和数据采样后变为了一类数据并使用了已知序列对信号功率进行编码,每加入一个节点就会产生R个信道,而对于每根接收天线来说,接收信号是独立的;在哈达玛矩阵的阶数足够大的情况下,就能够解出传输到每根天线上节点的信道信息,因此解决了所有信号同时到达网关时信道信息的无法求解的问题。
步骤4.4,对于K个节点中有个节点新加入***,/>个节点稳定传输的情况,/>个稳定传输的节点传输过程中的信道/>已知,并且数据部分以采样率/>采样后转换为已知符号,再利用哈达玛矩阵通过控制节点传输功率编码为已知序列,又由于/>个新加入的节点传输的符号已知,因此只需估计/>个信道/>,由于节点总数为K,因此依然需要编码矩阵至少为K+1阶;
其中,R个方程中的第一个方程中,第一个接收天线接收到的数据为:
式13;
R个方程中的其它方程中,接收天线接收到的数据的表达式在式13的基础上依次类推;
对于个稳定传输的节点,其信道已知但序列并未与新加入节点的序列对齐,由于哈达玛矩阵与编码矩阵不同行错位相乘求和结果为0,因此依然可以区分不同节点传输的信号;对于/>个新加入的节点,利用节点对应的序列与/>相乘求和,最后只剩下一项,而该项即为该节点对应的信道信息与/>的乘积;重复该过程能够求出所有未知信道;
接收天线数大于等于发射天线数,即需满足;
本实施例中,两个节点传输的信号如图4所示,第二个信号相比第一个信号增加了一个符号周期的时间延迟,图中红色圆圈标出的是采样后的序列。
在上述这种情况下即使可以求解出信道,也无法在信道均衡算法中求解传输的数据。因为受MIMO传输本身的限制,发射天线数量决定了采用信道均衡算法求解数据时未知数的个数,接收天线数量决定了每次传输所能建立的方程个数,因此,接收天线数不可少于发射天线数,即需满足/>;
步骤5,当符号边缘不对齐时,执行步骤5.1至步骤5.2;
对于各路传输信号异步到达且符号边缘未对齐的情况,如果没有信道信息就无法使用信道均衡算法进行信号检测,需要对该情况下的信道进行测量。符号边缘未对齐时,由于采样后序列也未对齐,因此接收信号并非是采样后序列的线性组合,而是由已知采样点与未知采样点线性组合而成;由于已知采样点减少了传输方程中未知数的个数,且在新节点加入***时已加入***并稳定传输的节点信道已知,因此也能够使用增加采样点的方式求解未知信息。
步骤5.1,以采样率开始采样;
步骤5.2,在1个节点已稳定传输,另一个节点新加入***,接收端使用2根接收天线的情况下,建立如式14所示的传输方程,通过该传输方程可以求解出所有未知量;
式14;
式中:
和/>分别对应第一个接收天线两个时刻的采样点,这两个采样点分别对应两路信号各自的采样点;
和/>则为另一个接收天线上的两个采样点;
和/>由于是已加入***的节点的信道,因此为已知量;
、/>、/>、/>均为未知量,共计4个;
表示表示节点传输的数据;
步骤6,在获取到全部参数后,信道矩阵H可以直接获得。
实施例2:
本实施例给出一种LoRa上行***中多节点MU-MIMO传输方法,如图2所示,包括以下步骤:
步骤一,利用多个节点给网关发送信号,不同节点的信号经过瑞利衰落信道进行传输,所有信号与噪声组合之后到达接收天线;
步骤二,对步骤一得到的信号,首先需要检测前导码;
前导码是接收端检测信号的依据,在检测到前导码后就可以进行信道估计,因此检测前导码是LoRa是否能启用MU-MIMO的前提条件;
步骤三,接收端进行分布式MU-MIMO信道估计,获取信道信息;
所述的分布式MU-MIMO信道估计的方法采用实施例1中的LoRa上行***中分布式MU-MIMO信道估计方法;
步骤四,将步骤三得到的信道信息,再将信道矩阵H代入信道均衡算法中即可解出MU-MIMO***中所有节点传输的信号,最后对各路信号分别使用LoRa解码算法,求解出各个节点传输的数据,从而完成了一个完整的MU-MIMO传输过程。
性能测试:
第一,信号检测算法比较:
在测试MU-MIMO的各种性能前,首先对各种信号检测算法进行对比,包括三种信道均衡算法,还加入一种采用预编码的方法SVD作为对比。
***中采用2个单天线节点进行发送,接收端使用4根接收天线进行接收,SF均设置为7,并在SNR分别为-5 dB、-10 dB、-15 dB、-20 dB、-25 dB、-30 dB 的情况下进行传输。统计各次传输的误比特率,测试结果如图5所示。图中MMSE具有最优的性能,尽管MRC在低SNR时与MMSE相似,但当SNR超过-15 dB 时,MRC的误比特率相对较高;ZF的性能与MMSE非常接近,但在低SNR时误比特率略高于MMSE。因此本发明将采用MMSE作为主要分析的MU-MIMO信号检测算法。
第二,SF和SNR对误比特率和吞吐量的影响:
为了研究采用MU-MIMO的LoRa***中影响***吞吐量的主要参数,讨论具有2个单天线发送端和一个4根天线接收端的MU-MIMO***在不同SNR时,误比特率随SF的变化情况,变化曲线如图6和图7所示。由图6可知,当采用不同SF并发传输时,MU-MIMO中单个节点的误比特率随SF的增加逐渐降低。这与节点单独传输时的情况非常相似,也验证了对于不同SF并发传输的情况,MU-MIMO可以有效分离各路信号。对于相同SF并发传输的情况,如图7所示,MU-MIMO中单个节点的误比特率在高于-15 dB 时几乎不随SF的变化而变化,而只与网关接收到节点传输信号的SNR相关,在低于-15 dB 时误比特率随SF的增加逐渐降低。
在获取了各SF的误比特率之后,可以计算得出该MU-MIMO***中单个节点的吞吐量变化情况,吞吐量可表示为:
式15;
式中,和/>分别表示以 SF 进行传输时信噪比为SNR的吞吐量和误比特率,/>表示chirp的周期。
由于该MU-MIMO***可以将混合信号分解为独立的数据流,当计算***吞吐量时只需将所有节点的吞吐量相加即可,因此,这里只给出单个节点的吞吐量。图8和图9中给出了在2×4的MU-MIMO***中单个节点吞吐量的变化情况。
由图8和图9可知MU-MIMO中单个节点的吞吐量随SNR的增大逐渐上升,此外,在使用SF7、SF8和SF9进行传输时,吞吐量受SNR的变化影响较大;而在使用SF10、SF11和SF12进行传输时,吞吐量受SNR的变化影响较小。从该图中还可以看出采用各SF传输时,SF越大吞吐量达到饱和时对应的SNR越低。
Claims (2)
1.一种LoRa上行***中分布式MU-MIMO信道估计方法,其特征在于,包括以下步骤:
步骤1,对于K个节点传输的信号y同时到达接收端时,由于K个节点都是新加入***的节点,因此有个信道/>未知;传输1个符号会产生/>个未知数和R个方程,因此采用的编码矩阵至少为K+1阶;
其中,R个方程中的第一个方程中,第一个接收天线接收到的数据为:
式12;
R个方程中的其它方程中,接收天线接收到的数据的表达式在式12的基础上依次类推;
j表示第j个接收天线;
i表示第i个节点;
表示第i个节点到第j个接收天线的信道;
s表示chirp;
表示第/>个采样点;
步骤2,接收端首先在接收到的信号中检测前导码,如果检测到前导码,则进行步骤3,如果未检测到前导码会持续检测,直到检测到前导码;
步骤3,在检测到前导码后,继续检测不同节点的符号边缘是否对齐;
当符号边缘对齐时,执行步骤4,当符号边缘不对齐时,执行步骤5;
步骤4,当符号边缘对齐时,执行步骤4.1至步骤4.4;
步骤4.1,首先计算预存的哈达玛矩阵阶数,然后接收等于哈达玛矩阵阶数/>的采样点;
步骤4.2,接收端预存阶数为K+1的哈达玛矩阵,以采样率对信号y进行采样;
其中,B表示信号带宽;SF表示扩频因子;
步骤4.3,每接收K+1个采样点就遍历哈达玛矩阵中的行来与相乘,遍历到了第一个节点所使用的序列,由于哈达玛矩阵与编码矩阵不同行相乘求和为0,该序列与/>、/>、…、/>相乘求和均为0,因此,只剩下一项/>,除以/>即可求出信道信息;
重复该过程,当完整遍历完一遍哈达玛矩阵后即可求得全部信道信息;
步骤4.4,对于K个节点中有个节点新加入***,/>个节点稳定传输的情况,个稳定传输的节点传输过程中的信道/>已知,并且数据部分以采样率/>采样后转换为已知符号,再利用哈达玛矩阵通过控制节点传输功率编码为已知序列,又由于/>个新加入的节点传输的符号已知,因此只需估计/>个信道/>,由于节点总数为K,因此依然需要编码矩阵至少为K+1阶;
其中,R个方程中的第一个方程中,第一个接收天线接收到的数据为:
式13;
R个方程中的其它方程中,接收天线接收到的数据的表达式在式13的基础上依次类推;
对于个稳定传输的节点,其信道已知但序列并未与新加入节点的序列对齐,由于哈达玛矩阵与编码矩阵不同行错位相乘求和结果为0,因此依然可以区分不同节点传输的信号;对于/>个新加入的节点,利用节点对应的序列与/>相乘求和,最后只剩下一项,而该项即为该节点对应的信道信息与/>的乘积;重复该过程能够求出所有未知信道;
接收天线数大于等于发射天线数,即需满足;
步骤5,当符号边缘不对齐时,执行步骤5.1至步骤5.2;
步骤5.1,以采样率开始采样;
步骤5.2,在1个节点已稳定传输,另一个节点新加入***,接收端使用2根接收天线的情况下,建立如式14所示的传输方程,通过该传输方程可以求解出所有未知量;
式14;
式中:
和/>分别对应第一个接收天线两个时刻的采样点,这两个采样点分别对应两路信号各自的采样点;
和/>则为另一个接收天线上的两个采样点;
和/>由于是已加入***的节点的信道,因此为已知量;
、/>、/>、/>均为未知量,共计4个;
表示表示节点传输的数据;
步骤6,在获取到全部参数后,信道矩阵H可以直接获得。
2.一种LoRa上行***中多节点MU-MIMO传输方法,其特征在于,包括以下步骤:
步骤一,利用多个节点给网关发送信号,不同节点的信号经过瑞利衰落信道进行传输,所有信号与噪声组合之后到达接收天线;
步骤二,对步骤一得到的信号,首先需要检测前导码;
步骤三,接收端进行分布式MU-MIMO信道估计,获取信道信息;
所述的分布式MU-MIMO信道估计的方法采用如权利要求1所述的LoRa上行***中分布式MU-MIMO信道估计方法;
步骤四,将步骤三得到的信道信息,再将信道矩阵H代入信道均衡算法中即可解出MU-MIMO***中所有节点传输的信号,最后对各路信号分别使用LoRa解码算法,求解出各个节点传输的数据,从而完成了一个完整的MU-MIMO传输过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111042442.9A CN113922850B (zh) | 2021-09-07 | 2021-09-07 | 一种LoRa上行***中分布式MU-MIMO信道估计方法及传输方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111042442.9A CN113922850B (zh) | 2021-09-07 | 2021-09-07 | 一种LoRa上行***中分布式MU-MIMO信道估计方法及传输方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113922850A CN113922850A (zh) | 2022-01-11 |
CN113922850B true CN113922850B (zh) | 2023-10-31 |
Family
ID=79234039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111042442.9A Active CN113922850B (zh) | 2021-09-07 | 2021-09-07 | 一种LoRa上行***中分布式MU-MIMO信道估计方法及传输方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113922850B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114553268B (zh) | 2022-02-23 | 2022-12-27 | 广东工业大学 | 一种基于MIMO的非相干LoRa通信*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112671680A (zh) * | 2021-01-20 | 2021-04-16 | 哈尔滨工业大学 | 一种改进误码性能的基于多倍带宽采样的LoRa解调方法 |
CN112804182A (zh) * | 2021-04-12 | 2021-05-14 | 中国人民解放军国防科技大学 | 一种多载波扩频方法及其应用的通信方法和相关装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10812130B2 (en) * | 2016-10-18 | 2020-10-20 | University Of Washington | Backscatter systems, devices, and techniques utilizing CSS modulation and/or higher order harmonic cancellation |
-
2021
- 2021-09-07 CN CN202111042442.9A patent/CN113922850B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112671680A (zh) * | 2021-01-20 | 2021-04-16 | 哈尔滨工业大学 | 一种改进误码性能的基于多倍带宽采样的LoRa解调方法 |
CN112804182A (zh) * | 2021-04-12 | 2021-05-14 | 中国人民解放军国防科技大学 | 一种多载波扩频方法及其应用的通信方法和相关装置 |
Non-Patent Citations (4)
Title |
---|
Aiju Thomas ; N V Eldhose.Performance Evaluation of Chirp Spread Spectrum as used in LoRa Physical Layer.2019 IEEE International Conference on System, Computation, Automation and Networking (ICSCAN).2019,全文. * |
Tallal Elshabrawy ; Joerg Robert.Analysis of BER and Coverage Performance of LoRa Modulation under Same Spreading Factor Interference.2018 IEEE 29th Annual International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC).2018,全文. * |
TD-SCDMA集群***中的联合功率控制;周围;舒鹏;;计算机应用(第02期);全文 * |
基于物联网的电缆局部放电分布式监测***设计与实现;谭笑;陈杰;刘洋;胡丽斌;贾勇勇;高旭泽;;电测与仪表;第57卷(第15期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113922850A (zh) | 2022-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9838230B2 (en) | Pre-coding in a faster-than-Nyquist transmission system | |
TWI443989B (zh) | 用於無線網路中通訊的方法,設備與系統及其製造物 | |
KR20050009233A (ko) | 채널 변동들을 보상하는 다중 안테나 모바일 통신시스템에서 신호들을 전송하는 방법 및 장치 | |
US8831128B2 (en) | MIMO communication system signal detection method | |
US7729458B2 (en) | Signal decoding apparatus, signal decoding method, program, and information record medium | |
KR101080624B1 (ko) | 다중-입력 다중 출력(mimo) 통신 링크의 성능 평가 | |
CN113922850B (zh) | 一种LoRa上行***中分布式MU-MIMO信道估计方法及传输方法 | |
US8553803B2 (en) | Estimation of the signal to interference and noise ratio (SINR) in multiple input multiple output (MIMO) systems | |
Uthansakul et al. | Multipath signal effect on the capacity of MIMO, MIMO-OFDM and spread MIMO-OFDM | |
US20090262869A1 (en) | Decision-Feedback Detection for Block Differential Space-Time Modulation | |
Haider et al. | Mimo network and the alamouti, stbc (space time block coding) | |
CN104618297A (zh) | Ofdm-dcsk通信*** | |
George et al. | Performance analysis of transmission modes in 4G LTE/LTE-A downlink systems | |
Singh et al. | BER performance analysis of OFDM-MIMO system using GNU Radio | |
Waliullah et al. | Study the BER performance comparison of MIMO systems using BPSK modulation with ZF and MMSE equalization | |
Anil et al. | Implementation of MIMO OFDM NOMA system using iterative algorithm | |
Raghavendra et al. | Co-channel interference mitigation for 3G LTE MIMO-OFDM systems | |
Miatton | Interference alignment at intermediate SNR with perfect or noisy CSI | |
Le et al. | Efficient algorithm for blind detection of orthogonal space-time block codes | |
Hassan | A Simple Time–Space Diversity Scheme Using Extended Shaping Pulse | |
Singh et al. | BER-SNR Evaluation using Phase Shift Keying with OSTBC in MIMO Environment | |
Avendi | Differential modulation and non-coherent detection in wireless relay networks | |
Nguyen et al. | Recursive receivers for diversity channels with correlated flat fading | |
Dawi et al. | Performance of DSTM MIMO Systems with 2, 4 and 8 Transmit Antennas Using Extensions of the Weyl Group | |
Hu | Joint communication and sensing in massive MIMO architectures with low-precision converters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |