CN105245477B - 一种低复杂度差分空间调制检测算法 - Google Patents
一种低复杂度差分空间调制检测算法 Download PDFInfo
- Publication number
- CN105245477B CN105245477B CN201510551207.2A CN201510551207A CN105245477B CN 105245477 B CN105245477 B CN 105245477B CN 201510551207 A CN201510551207 A CN 201510551207A CN 105245477 B CN105245477 B CN 105245477B
- Authority
- CN
- China
- Prior art keywords
- matrix
- moment
- detection
- difference
- transmitting
- 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
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 70
- 238000001514 detection method Methods 0.000 title claims abstract description 48
- 239000011159 matrix material Substances 0.000 claims abstract description 73
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 18
- 238000013507 mapping Methods 0.000 claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 230000004913 activation Effects 0.000 claims description 10
- 238000003032 molecular docking Methods 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 3
- 230000005540 biological transmission Effects 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Landscapes
- Radio Transmission System (AREA)
- Error Detection And Correction (AREA)
Abstract
本发明提供了一种LC‑DSM算法(低复杂度的差分检测算法)，首先，在差分空间调制***发射端应用比特填充方案与传输比特计算公式构造符号矩阵，经过差分变换与发送在接收端得到接收信号，再以差分最大似然检测算法的算式结构为基础，对接收矩阵采取分列检测的搜索方式，逐列进行信号判决，得到发射符号与发射天线的序号，最终根据比特与发射符号、比特与发射天线序号的对应关系逆映射为信息比特。该方法是将空间调制下低复杂度的检测算法与差分空间调制***下最大似然检测算法形成的算式结构相结合而得到的。研究表明通过使用LC‑DSM算法，在保证算法性能没有损失的前提下，相对于差分最大似然检测算法复杂度得到了大大的降低。
Description
技术领域
本发明涉及通信技术领域，尤其涉及无线通线***发射端的发射方法与信号检测方法，具体地说是一种用于空间调制***的低复杂度检测算法。
背景技术
多输入多输出***(MIMO：multiple input multiple output)技术通过使用多根发射和接收天线，能够提高无线通信***的容量与频谱利用率。但由于其***结构的限制，传输信道间的干扰与多天线发射严格同步性成为制约其发展的两个主要原因。在这样的情况下，空间调制(SM：spatial Modulation)作为一种新型的MIMO发射方式被提出，***在每一个发送时隙，只有一根发射天线被激活用来发送数据，这种特性使得空间调制***可以有效的避免信道间干扰和多天线同步的问题。但由于在高速移动的场景中很难获得准确的信道状态信息，性能损失变得严重。
基于上述问题，差分信号作为解决方案被引入到空间调制***当中。有学者提出了基于分组的差分空时相移键控(DSTSK：differential modulation for space-timeshift keying)，但DSTSK需要穷尽搜索线性分散矩阵并且只传输了单个符号。为了解决这两个问题，基于符号的差分空间调制(DSM：differential spatial modulation)被提出，通过被激活的发射天线来构造传输的符号矩阵，不仅降低了穷尽搜索分散矩阵的限制，同时由于多符号的传输增加了频谱利用率。但随着传输符号数的增多，DSM***的检测算法的复杂度也变得很高.
现有针对差分空间调制的的研究工作主要集中在探寻具有高性能低复杂度的检测方法上。差分最大似然检测算法(DSM-ML)作为性能最优的检测算法被提出，但由于其采用遍历搜索的方式使得其复杂度非常高。基于此，本发明结合差分最大似然检测算法的算式结构提出了一种低复杂度的差分检测算法。
发明内容
本发明所要解决的技术问题是提供一种新的算法，称作LC-DSM算法(低复杂度的差分检测算法)，该方法是将空间调制下低复杂度的检测算法与差分空间调制***下最大似然检测算法形成的算式结构相结合而得到的。研究表明通过使用LC-DSM算法，在保证算法性能没有损失的前提下，相对于差分最大似然检测算法复杂度得到了大大的降低。
为此，本发明提供以下技术方案：
一种低复杂度差分空间调制检测算法，首先，在差分空间调制***发射端应用比特填充方案(Yuli Yang,Aissa S.Bit-Padding Information Guided Channel Hopping,IEEE Communications Letters,2011,15,(2):163-165)与传输比特计算公式构造符号矩阵，经过差分变换与发送在接收端得到接收信号(接收信号为Nr×Nt的矩阵)，再以差分最大似然检测算法的算式结构为基础，对接收矩阵采取分列检测的搜索方式，逐列进行信号判决，得到发射符号与发射天线的序号，最终根据比特与发射符号、比特与发射天线序号的对应关系逆映射为信息比特。在采用上述技术方案的基础上，本发明还可采用以下进一步的技术方案：
所述差分空间调制***的发射端设为Nt×Nt维符号矩阵，Nt×Nt维符号矩阵可能个数的计算公式为：
A＝Nt！×(M)Nt………(1)
发射过程中发射一个符号矩阵所需要的信息比特个数为：
其中Nt代表发射天线的个数，M代表调制阶数，表示向下取整；
将信息比特每个为一组映射成一个Nt行Nt列的符号矩阵X(X中任一行任一列有且只有一个非零元素，满足不同时隙激活不同天线发射单一符号的要求)。
将得到的符号矩阵X进行差分变换St＝St-1Xt，其中St为t时刻经过差分变化之后得到的差分接收矩阵，Xt为t时刻发送的符号矩阵；将Xt加载在发射天线上进行发射，在接收端得到接收信号：
Yt＝HtSt+Vt………(3)；
式中Yt、Ht、Vt分别为t时刻的接收矩阵，信道增益矩阵与高斯白噪声。
差分最大似然检测算法公式为：
式中Θ为Nt！MNt种符号矩阵的集合；以差分最大似然检测算法的检测结构为基础，将在t-1时刻得到的接收信号视为t时刻检测发射信号所需要的信道增益矩阵，(4)式可化为：
式中分别表示第t个时刻被激活的发射天线与加载在天线上的调制符号，ht|l表示第t个时刻的第l根发射天线对应的信道增益矩阵，Q为调制符号的集合，s为符号矩阵中的候选调制符号；
将接收矩阵Yt按列取出，然后逐列进行信号检测，直至判决出所有的天线序号与调制符号，最终逆映射得到发射端的信息比特。
由于采用本发明的技术方案，本发明的有益效果为：
(1)在发射端舍弃传统的映射方法，转而采用比特填充方案，解决了由发射天线个数导致的接收端判决得到的天线映射方案无法逆映射的问题，同时增加了每个符号矩阵携带的信息比特数，提高了频谱利用率。
(2)本发明算法在接收端以差分最大似然检测算法的算式结构为基础，将接收信号矩阵按列拆分，再进行信号检测，将参与搜索的候选符号从符号矩阵转变为可在空间调制***下进行检测的单一符号向量，从而在检测方式上将基于矩阵的差分最大似然检测转化为基于向量的相干空间调制检测算法。
(3)本发明算法的复杂度与调至阶数无关，降低了调制阶数的影响，同时也降低了算法复杂度。
(4)在检测过程中信道矩阵会根据检测到的天线序号进行缩小，避免了多次利用无关的信道增益向量与候选符号相乘，进一步降低了复杂度。
附图说明
图1是差分空间调制***框图；
图2是本发明LC-DSM算法检测示意图；
图3是本发明LC-DSM算法与DSM-ML算法的性能比较图；
具体实施方式
本发明采用的一种低复杂度差分空间调制检测算法为：首先，在差分空间调制***发射端应用比特填充方案与传输比特计算公式构造符号矩阵，将信号经过差分变化后进行发送，在接收端得到接收信号(接收信号为Nr×Nt的矩阵)，再以差分最大似然检测算法的算式结构为基础，对接收矩阵采取分列检测的搜索方式，逐列进行信号判决，得到调制符号与发射天线的序号，最终逆映射为信息比特。
进一步的，所述差分空间调制***的发射端设为Nt×Nt维符号矩阵，Nt×Nt维符号矩阵可能个数的计算公式为：
A＝Nt！×(M)Nt………(1)
发射过程中发射一个符号矩阵所需要的信息比特个数为：
其中Nt代表发射天线的个数，M代表调制阶数，表示向下取整；
将信息比特每个为一组映射成一个Nt行Nt列的符号矩阵X(X中任一行任一列有且只有一个非零元素，满足不同时隙激活不同天线发射单一符号的要求)。
将得到的符号矩阵X进行差分变换St＝St-1Xt，其中St为t时刻经过差分变化之后得到的差分矩阵，Xt为t时刻发送的符号矩阵；将Xt加载在发射天线上进行发射，在接收端得到接收信号：
Yt＝HtSt+Vt………(3)；
式中Yt、Ht、Vt分别为t时刻的接收矩阵，信道增益矩阵与高斯白噪声。
差分最大似然检测算法公式为：
式中Θ为Nt！MNt种符号矩阵的集合；以差分最大似然检测算法的检测结构为基础，将在t-1时刻得到的接收信号视为t时刻检测发射信号所需要的信道增益矩阵，(4)式可化为：
式中分别表示第t个时刻被激活的发射天线与加载在天线上的调制符号，ht|l表示第t个时刻的第l根发射天线对应的信道增益矩阵，Q为调制符号的集合，s为符号矩阵中的候选调制符号；
将接收矩阵Yt按列取出，然后逐列进行信号检测，直至判决出所有的天线序号与调制符号，最终逆映射得到发射端的信息比特。
下面主要结合附图对本发明的具体实施进行详细描述。
图1表示的是差分空间调制***的原理框图。
我们考虑DSM***具有Nt根发射天线和Nr根接收天线。发射端的信息比特被分成两部分，一部映射成为发射天线激活方案，另一部分经过调制变成符号；然后将得到的符号加载在被激活的发射天线上，我们得到发射符号矩阵X；将得到的符号矩阵X进行差分变换之后得到发射矩阵St；将St加载在发射天线上进行发送，在接收端得到接收矩阵Yt；Yt经过LC-DSM算法运算后得到发射天线激活方案与调制符号，再经过逆映射得到最终的比特信息。图中Map表示信息比特的映射，delay表示进行差分变换，LC-DSM表示本专利提出的检测算法，Demap to bits表示将检测得到的符号与天线序号映射成信息比特。
LC-DSM算法在检测时采取分列检测的搜索方式，即对接收到的信号矩阵Yt进行逐列检测.假设Nt＝3,Nr＝2，调制方式为QPSK的差分空间调制***，算法检测示意如图2所示。
LC-DSM算法首先将接收矩阵Yt的第一列Yt|1单独取出，由于参与搜索的候选符号不再是符号矩阵而是单一的符号向量，从而在检测方式上将基于矩阵的差分最大似然检测转化为基于向量的相干空间调制检测算法.根据接收向量信号Yt|1检测出对应的天线序号与符号后，将Ht缩小为2×2的矩阵.取出Yt|2进行下一次信号检测得重复上一步直到发射天线与发射符号全部检测出来。
根据实验分析，当LC-DSM算法采取逐列搜索的方式时，假设Nt＝4，个比特可以映射16种天线激活方案，但根据矩阵维度实际有Nt！＝24种天线选择方案.如表1所示，
表1
记录了第16到24种天线选择方案。在这种情况下采用LC-DSM算法进行分列检测出的天线激活方案可能超出这16种情况，以至于天线序号无法逆映射为信息比特，如表1中第17到24中方案所示。
为了在发射端解决这个问题，本发明引入了比特填充方案。如表2所示，
表2
表2列出了发射天线数为2和3的情况。其中Nt！表示发射天线的激活方案。由上述分析可得，当Nt不为2的整数次幂时实际采用了种天线激活方案，在接收端进行分列检测时可能检测到的天线激活方案属于中的一种，这样会导致无法根据天线激活顺序进行逆映射。
在比特填充方案中首先根据(6)式计算出k值
2k＜Nt！＜2k+1 (6)
图3表示DSM-ML算法(差分最大似然检测算法)与LC-DSM算法在传输效率为2.5bit/s/Hz，Nt＝2，Nr＝2,3,4情况下的性能比较结果。
表3
从上面的图形可见，LC-DSM在低信噪比区域相比于DSM算法的性能略差，但在较高信噪比区域LC-DSM算法与DSM-ML算法的性能一致。
在分析算法复杂度时我们采用计算算法的实数次乘法。
LC-DSM算法的复杂度可由下面计算公式得到：
DSM-ML算法的复杂度可由下面计算公式得到：
表3给出了不同情况下DSM-ML算法与LC-DSM算法的复杂度比较结果。从表中可以看到，在Nt＝2，M＝4，Nr分别等于2、3、4的情况下，LC-DSM算法较DSM-ML算法复杂度降低了80％左右；在Nr＝2，M＝4，Nt分别等于2、3、4的情况下，随着Nt的增大，LC-DSM算法较DSM-ML算法的复杂度降低的更多；在Nt＝2，Nr＝2，M分别等于2、4、8的情况下由于LC-DSM算法的复杂度与调制阶数M无关，此时LC-DSM算法的低复杂度特性表现的更为明显。因此，根据以上分析得出LC-DSM算法在基本保持了原有DSM-ML性能，且在高信噪比性能优于DSM-ML算法；另外LC-DSM算法相比于DSM-ML算法大大降低了复杂度。
Claims (3)
1.一种低复杂度差分空间调制检测方法，其特征在于，首先，在差分空间调制***发射端应用比特填充方案与传输比特计算公式构造符号矩阵，将信息比特经过差分变化后进行发送，在接收端得到接收信号，差分空间调制***接收信号为Nr×Nt维的矩阵，其中Nr为接收天线的个数，Nt为发射天线的个数，再以差分最大似然检测算法的算式结构为基础，对接收矩阵采取分列检测的搜索方式，逐列进行信号判决，得到调制符号与发射天线的序号，最终逆映射为信息比特；
所述差分空间调制***发射信号为Nt×Nt的维符号矩阵，Nt×Nt维的符号矩阵的个数A的计算公式为：
A＝Nt！×(M)Nt………(1)；
发射过程中发射一个符号矩阵所需要的信息比特个数为：
其中Nt代表发射天线的个数，M代表调制阶数，表示向下取整；
将信息比特每个为一组映射成一个Nt行Nt列的符号矩阵X，X中任一行任一列有且只有一个非零元素，满足不同时隙激活不同天线发射单一符号的要求。
2.如权利要求1所述的一种低复杂度差分空间调制检测方法，其特征在于，
将得到的符号矩阵进行差分变换St＝St-1Xt，其中St为t时刻经过差分变化之后得到的差分矩阵，St-1为t-1时刻经过差分变化之后得到的差分矩阵，Xt为t时刻发送的符号矩阵；将Xt加载在发射天线上进行发射，在接收端得到接收信号：
Yt＝HtSt+Vt………(3)；
式中Yt、Ht、Vt分别为t时刻的接收矩阵，信道增益矩阵与高斯白噪声。
3.如权利要求1所述的一种低复杂度差分空间调制检测方法，其特征在于，差分最大似然检测算法公式为：
式中Θ为Nt！MNt种符号矩阵的集合，Yt-1为t-1时刻的接收矩阵；以差分最大似然检测算法的检测结构为基础，将在t-1时刻得到的接收信号视为t时刻检测发射信号所需要的信道增益矩阵，公式(4)化为：
式中分别表示第t个时刻被激活的发射天线与加载在天线上的调制符号，i代表第几根天线的索引值，ht|l表示第t个时刻的第l根发射天线对应的信道增益矩阵,Yt|i是指第t时刻接收信号矩阵的第i列，Q为调制符号的集合，S为符号矩阵中的候选调制符号，ht|lS代表ht|l矩阵与S符号相乘结果；Xt为t时刻发送的符号矩阵，式中Yt、Ht、Vt分别为t时刻的接收矩阵，信道增益矩阵与高斯白噪声；
将接收矩阵Yt按列取出，然后逐列进行信号检测，直至判决出所有的调制符号与发射天线的序号，最终根据信息比特与调制符号、信息比特与发射天线的序号对应关系逆映射为信息比特。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510551207.2A CN105245477B (zh) | 2015-09-01 | 2015-09-01 | 一种低复杂度差分空间调制检测算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510551207.2A CN105245477B (zh) | 2015-09-01 | 2015-09-01 | 一种低复杂度差分空间调制检测算法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105245477A CN105245477A (zh) | 2016-01-13 |
| CN105245477B true CN105245477B (zh) | 2018-09-07 |
Family
ID=55042981
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510551207.2A Active CN105245477B (zh) | 2015-09-01 | 2015-09-01 | 一种低复杂度差分空间调制检测算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105245477B (zh) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106027207B (zh) * | 2016-05-27 | 2018-12-21 | 华北电力大学（保定） | 一种差分空间调制方法 |
| CN107959519B (zh) * | 2016-10-17 | 2022-06-24 | 北京三星通信技术研究有限公司 | 一种差分空间调制传输方法、发射机和接收机 |
| CN106850475B (zh) * | 2017-02-13 | 2020-05-12 | 重庆邮电大学 | 一种应用于空间调制的低复杂度排序a*检测算法 |
| CN107483376A (zh) * | 2017-09-07 | 2017-12-15 | 西安电子科技大学 | 一种用于mimo‑ofdm***的信号检测方法 |
| CN108289013B (zh) * | 2018-01-05 | 2021-01-08 | 中国计量大学 | 一种基于补码技术的差分空间调制协作***抗干扰方法 |
| CN108449120B (zh) * | 2018-03-14 | 2020-05-19 | 中南大学 | 一种基于差分度量低复杂度qam-mimo检测方法 |
| CN108616475B (zh) * | 2018-04-17 | 2020-12-18 | 电子科技大学 | 一种用于ofdm***的子载波索引差分调制方法 |
| CN108989262B (zh) * | 2018-08-08 | 2021-04-06 | 中国计量大学 | 一种基于apsk调制的低复杂度非相干空间调制检测方法 |
| CN109412670B (zh) * | 2018-11-28 | 2022-01-18 | 中国计量大学 | Gsm-mbm***中基于松弛迭代的低复杂度检测方法 |
| CN109818663B (zh) * | 2018-12-21 | 2023-10-03 | 中国计量大学 | 一种低复杂度差分正交空间调制检测方法 |
| CN109547077B (zh) * | 2019-01-22 | 2020-10-13 | 重庆京东方智慧电子***有限公司 | 信号发送方法、信号接收方法、通信设备及存储介质 |
| CN110855328B (zh) * | 2019-10-25 | 2021-01-19 | 西安交通大学 | 一种基于天线分组的差分空间调制方法、设备及存储介质 |
| CN110855333B (zh) * | 2019-11-13 | 2021-05-11 | 中国计量大学 | 一种基于广义空间调制的差分传输方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104333434A (zh) * | 2014-08-31 | 2015-02-04 | 电子科技大学 | 一种低复杂度的空间调制检测方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100870557B1 (ko) * | 2001-05-25 | 2008-11-27 | 리전츠 어브 더 유니버시티 오브 미네소타 | 무선통신망에서의 공간-시간 부호화 전송 장치 및 방법 |
-
2015
- 2015-09-01 CN CN201510551207.2A patent/CN105245477B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104333434A (zh) * | 2014-08-31 | 2015-02-04 | 电子科技大学 | 一种低复杂度的空间调制检测方法 |
Non-Patent Citations (2)
| Title |
|---|
| 《Reduced complexity sphere decoder for spatial modulation detection receivers》;YOUNIS A;《IEEE》;20101230;全文 * |
| 《低复杂度最大似然MIMO信号检测算法》;赵新雪等;《计算机工程与设计》;20140116;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105245477A (zh) | 2016-01-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105245477B (zh) | 一种低复杂度差分空间调制检测算法 | |
| CN110098855B (zh) | 一种可见光mimo通信预编码及解码方法 | |
| CN105429709B (zh) | 一种mu‑miso可见光通信***迫零预编码矩阵确定方法 | |
| CN105827293B (zh) | 一种多用户广义空间调制***线性迫零预编码方法 | |
| CN107634787A (zh) | 一种大规模mimo毫米波信道估计的方法 | |
| CN108989262B (zh) | 一种基于apsk调制的低复杂度非相干空间调制检测方法 | |
| CN108540185A (zh) | 一种结合空时分组码的差分空间调制方法 | |
| CN105246158A (zh) | 基于高信噪比的能效最大化多天线中继***功率分配方法 | |
| CN106301496B (zh) | 基于天线选择和预编码的空间调制*** | |
| CN107743043A (zh) | 多用户毫米波***中基于带外空间信息的用户分组方法 | |
| CN105933042A (zh) | 一种lte***中基于分簇的自适应有限反馈新方法 | |
| CN109995404A (zh) | 空频域调制的一种差分调制和解调方法 | |
| CN106130613B (zh) | 一种基于酉空时码获得灵活发射分集的空间调制方法 | |
| CN104092516B (zh) | 一种适用于mu‑mimo***的非正交预编码码本设计方法 | |
| KR101055573B1 (ko) | 다중 사용자, 다중 안테나 무선 송출 시스템에서의 프리 코딩 장치 | |
| CN103326825B (zh) | 一种准正交空时分组码低复杂度译码方法 | |
| US11108458B2 (en) | Method and apparatus for combining plurality of radio signals | |
| Chi et al. | Sphere decoding for spatial permutation modulation MIMO systems | |
| CN107948114B (zh) | 一种mimo盲源信号分离***及信号分离方法 | |
| CN110190880A (zh) | 一种广义空间调制递增发射天线组合选择方法 | |
| CN102340377A (zh) | 一种非线性预编码处理方法、装置和*** | |
| CN113708813B (zh) | 一种基于波束赋形的多用户空间调制方法 | |
| CN109818663A (zh) | 一种低复杂度差分正交空间调制检测方法 | |
| CN105119869B (zh) | 基于空时矩阵星座图的空移键控通信方法 | |
| CN111901022B (zh) | 一种预编码辅助的信号发送和接收方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20231106 Address after: Room 407-10, floor 4, building 2, Haichuang science and technology center, Cangqian street, Yuhang District, Hangzhou City, Zhejiang Province, 311100 Patentee after: Zhejiang Zhiduo Network Technology Co.,Ltd. Address before: 315470 126 Guangming Road, Simen Town, Yuyao City, Ningbo, Zhejiang Patentee before: China Jiliang University |