CN101707579B - 连续相位调制频域均衡调制解调方法 - Google Patents

连续相位调制频域均衡调制解调方法 Download PDF

Info

Publication number
CN101707579B
CN101707579B CN2009101545692A CN200910154569A CN101707579B CN 101707579 B CN101707579 B CN 101707579B CN 2009101545692 A CN2009101545692 A CN 2009101545692A CN 200910154569 A CN200910154569 A CN 200910154569A CN 101707579 B CN101707579 B CN 101707579B
Authority
CN
China
Prior art keywords
continuous phase
symbolic blocks
modulation
symbol
frequency domain
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.)
Expired - Fee Related
Application number
CN2009101545692A
Other languages
English (en)
Other versions
CN101707579A (zh
Inventor
赵民建
钟杰
付涛
吕春苗
彭岑
郑雅敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2009101545692A priority Critical patent/CN101707579B/zh
Publication of CN101707579A publication Critical patent/CN101707579A/zh
Application granted granted Critical
Publication of CN101707579B publication Critical patent/CN101707579B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明公开了一种连续相位调制频域均衡调制解调方法,该方法可用于频率选择性衰落信道中。连续相位调制信号的产生包括映射、调制、***循环前缀;连续相位调制信号的解调包括接收经历了频率选择性衰落信道下的信号,移除循环前缀、正交分解、信道估计、频域信道均衡、变换回时域后维特比译码。本发明通过正交分解的方法,消除CPM自身的相关性,后将其变换到频域做信道估计与均衡,大大的降低了复杂度,尤其是在频率选择性衰落信道下。本发明通过频域均衡的方法可以有效的抵消信道的影响,消除信道影响带来的码间干扰,因此可以工作在频率选择性衰落信道下。

Description

连续相位调制频域均衡调制解调方法
技术领域
本发明涉及信号调制解调技术,特别地,涉及一种连续相位调制频域均衡调制解调方法。
背景技术
CPM(连续相位调制)信号具有很好的恒包络特性,这使得它经过高效的非线性功放放大而不受影响,功率效率高。同时,CPM调制信号具有连续相位的特性,没有相位突变,从而没有高频分量。因此,CPM信号频谱紧凑,其频带利用率较高。在同样的码元速率下,CPM信号的带宽比普通FSK和PSK调制都要小一些。CPM调制信号由于其良好的恒包络特性、高功率效率、高频谱效率,得到了广泛的应用。
随着通信技术的发展,以及通信业务的需求变得更加多样化,对通信的带宽要求越来越高。在无线信道传输中,宽带信号比窄带信号经历的衰落要大很多,为频率选择性信道,因此,抗衰落技术的选取十分重要。在衰落信道条件下,CPM最佳的时域均衡技术是采用维特比算法的最大似然序列估计(MLSE)方法,MLSE在信道弥散长度变长的情况下有非常高的复杂度,即使是使用降低复杂度的延迟判决反馈序列估计方法(DDFSE),仍然具有很高的复杂度,这样,寻求CPM的频域均衡显得很有意义。
发明内容
本发明的目的是克服现有技术的不足,提供一种连续相位调制频域均衡调制解调方法,该方法可用于频率选择性衰落信道中。
包括如下步骤:
1)从数据源获得二进制序列,将二进制序列映射成调制符号序列,将调制符号或训练字符号组成符号块,按照符号块长度逐块送入连续相位调制器;
2)每个符号块在完成连续相位调制之后,***循环前缀波形,得到连续相位调制信号;
3)将连续相位调制信号在通信通道上发送,连续相位调制信号经历频率选择性衰落;
4)接收衰落后的连续相位调制信号,移除循环前缀,得到经历衰落的训练字符号块和数据符号块,用存储于本地的正交基对每个训练字符号块或数据符号块做正交分解;
5)正交分解之后,把训练字符号块变换到频域做信道估计,数据符号块则根据信道估计结果,在频域上做均衡,去除信道干扰的影响后再变换回时域;
6)均衡后变换回时域的数据符号块送进维特比译码器解调得到源序列。
所述的从数据源获得二进制序列,将二进制序列映射成调制符号序列,将调制符号或训练字符号组成符号块,按照符号块长度逐块送入连续相位调制器的步骤为:
1)将从数据源获得的二进制序列an映射为数据符号序列In,具体是将1映射为+1,0映射为-1;
2)每252个数据符号In和4个归零符号组成一个数据符号块,数据符号块的长度为256;
3)每隔20个数据符号块***一个训练字符号块,以训练字符号块为起点组成帧结构,训练字符号块仍然是由长度为252的训练字符号和4个归零符号组成;
4)每个符号块依次进入连续相位调制器,在每个符号块的起点相位为0,同时每个符号块的第231(256-25)个符号处相位也为0,连续相位调制信号波形表达式为其中
Figure G2009101545692D00022
( t ) = ∫ - ∞ t ( τ ) dτ 为成形函数,In为数据符号块或训练序列符号块,
Figure G2009101545692D00024
h=0.5。
所述的每个符号块在完成连续相位调制之后,***循环前缀波形,得到连续相位调制信号的步骤为:
1)循环前缀长度的大小要大于信道最大的弥散长度,设定频率选择性衰落信道的弥散长度最大为25个符号,故循环前缀的长度设定为25个符号;
2)在第229(256-25-2)个符号处将相位归零,这将用掉2个相位归零符号;
3)在每个符号块的末尾,需要加上另外2个归零符号,使相位从任意状态归0;
4)对于每个符号块,将最后25个符号的连续相位调制信号波形复制到符号块最前面形成循环前缀。
所述的将连续相位调制信号在通信通道上发送,连续相位调制信号经历频率选择性衰落的步骤为:
1)连续相位调制信号
Figure G2009101545692D00031
经历频率选择性衰落信号后表示为: r ( t ) = Σ l = 0 M D - 1 h l s ( t - lT ) + z ( t ) , 其中z(t)是零均值复高斯噪声,方差为N0;MD是信道的最大弥散长度,hl为信道第L径的增益。
所述的接收衰落后的连续相位调制信号,移除循环前缀,得到经历衰落的训练字符号块和数据符号块,用存储于本地的正交基对每个训练字符号块或数据符号块做正交分解的步骤为:
1)接收经历了频率选择性衰落信道的连续相位调制信号r(t),对r(t)先移除循环前缀;
2)加上倾斜相位
Figure G2009101545692D00033
使相位网格状态时不变,得到r′(t);
3)存储于本地的正交基为f0(t)=1, f 1 ( t ) = 1 1 - 4 π 2 ( e jπt T - j 2 π ) , 0≤t<T,其中T为符号周期;将加入了倾斜相位的接收信号r′(t)在符号周期内与两个正交基作相关,得到两组向量r0(n)和r1(n)。
所述的正交分解之后,把训练字符号块变换到频域做信道估计,数据符号块则根据信道估计结果,在频域上做均衡,去除信道干扰的影响后再变换回时域的步骤为:
1)正交分解将训练字符号块和数据符号块同时分解为两路向量;在两路向量中分别提取出训练字符号块的向量段r0 T(n)、r1 T(n);
2)将r0 T(n)、r1 T(n)和存储于本地的训练字符号块向量s0 T(n)、s1 T(n)一起变换到频域得到R0 T(k)、R1 T(k)和S0 T(k)、S1 T(k);
3)根据最小均方误差准则,得到正交分解后信道频域的估计值为:
H ~ 0 = H 0 * | H 0 | 2 + N 0 , H ~ 1 = H 1 * | H 1 | 2 + N 0 ( H 0 = R 0 T / S 0 T , H 1 = R 1 T / S 1 T );
4)将数据符号块正交分解后得到的向量r0(n)和r1(n)变换到频域上得到R0,R1,根据信道估计结果,在频域去除信道的影响,得到两路均衡后的输出
Figure G2009101545692D00039
Figure G2009101545692D000310
S 0 ~ = H 0 ~ * R 0 , S 1 ~ = H 1 ~ * R 1 ;
5)将
Figure G2009101545692D000313
Figure G2009101545692D000314
变换到时域,送进解调器解调。
本发明通过正交分解的方法,消除CPM自身的相关性,后将其变换到频域做信道估计与均衡,大大的降低了复杂度,尤其是在频率选择性衰落信道下。本发明通过频域均衡的方法可以有效的抵消信道的影响,消除信道影响带来的码间干扰,因此可以工作在频率选择性衰落信道下。
附图说明
图1是依据本发明的实施例的一种示范性通信***的方框图;
图2是依据本发明的实施例中的符号块结构图;
图3是依据本发明的实施例中的循环前缀添加示意图;
图4是依据本发明的实施例中的CPM信号正交分解方框图;
图5是依据本发明的实施例中的CPM频域均衡方框图。
具体实施方式
本发明实现了CPM信号利用频域均衡技术,在频率选择性衰落信道下的调制解调。
图1示出了本发明的实施方案的一种示范性通信***。发送部分由数据源101、映射器102、成帧103、CPM调制器104、加入循环前缀105模块相连,完成CPM信号的产生过程;接收部分接收经历频率选择性衰落信道的CPM信号,并进行解调处理,由循环前缀移除107、正交分解108、信道估计110、频域均衡109和解调111模块组成,最终解调出二进制流,恢复数据源。
连续相位调制频域均衡调制解调方法包括如下步骤:
1)从数据源获得二进制序列,将二进制序列映射成调制符号序列,将调制符号或训练字符号组成符号块,按照符号块长度逐块送入连续相位调制器;
2)每个符号块在完成连续相位调制之后,***循环前缀波形,得到连续相位调制信号;
3)将连续相位调制信号在通信通道上发送,连续相位调制信号经历频率选择性衰落;
4)接收衰落后的连续相位调制信号,移除循环前缀,得到经历衰落的训练字符号块和数据符号块,用存储于本地的正交基对每个训练字符号块或数据符号块做正交分解;
5)正交分解之后,把训练字符号块变换到频域做信道估计,数据符号块则根据信道估计结果,在频域上做均衡,去除信道干扰的影响后再变换回时域;
6)均衡后变换回时域的数据符号块送进维特比译码器解调得到源序列。
所述的从数据源获得二进制序列,将二进制序列映射成调制符号序列,将调制符号或训练字符号组成符号块,按照符号块长度逐块送入连续相位调制器的步骤为:
1)将从数据源获得的二进制序列an映射为数据符号序列In,具体是将1映射为+1,0映射为-1;
2)每252个数据符号In和4个归零符号组成一个数据符号块,数据符号块的长度为256,如图2所示,可以看出,在一个符号块中,数据符号的个数为252个,分为两段,第一段229个(201),第二段23个(203),中间和最后分别有一组归零符号(202、204);
3)每隔20个数据符号块***一个训练字符号块,以训练字符号块为起点组成帧结构,训练字符号块仍然是由长度为252的训练字符号和4个归零符号组成;
4)每个符号块依次进入连续相位调制器,在每个符号块的起点相位为0,同时每个符号块的第231(256-25)个符号处相位也为0,连续相位调制信号波形表达式为其中
Figure G2009101545692D00052
q ( t ) = ∫ - ∞ t g ( τ ) dτ 为成形函数,In为数据符号块或训练序列符号块,h=0.5。
图3示意了在本发明实施例中循环前缀的添加方式,由于CPM信号是在相位上进行调制的,故需要在CPM调制之后添加循环前缀,具体是将图3中304、305模块复制形成301得到循环前缀。而为了保证CPM信号相位的连续性,需要加入归零符号使相位状态归零,见图3中303、305。
所述的每个符号块在完成连续相位调制之后,***循环前缀波形,得到连续相位调制信号的步骤为:
1)循环前缀长度的大小要大于信道最大的弥散长度,设定频率选择性衰落信道的弥散长度最大为25个符号,故循环前缀的长度设定为25个符号;
2)在第229(256-25-2)个符号处将相位归零,这将用掉2个相位归零符号;
3)在每个符号块的末尾,需要加上另外2个归零符号,使相位从任意状态归0;
4)对于每个符号块,将最后25个符号的连续相位调制信号波形复制到符号块最前面形成循环前缀,参见图3中301。
所述的将连续相位调制信号在通信通道上发送,连续相位调制信号经历频率选择性衰落的步骤为:
1)连续相位调制信号经历频率选择性衰落信号后表示为: r ( t ) = Σ l = 0 M D - 1 h l s ( t - lT ) + z ( t ) , 其中z(t)是零均值复高斯噪声,方差为N0;MD是信道的最大弥散长度,hl为信道第L径的增益。
图4示出了本发明实施例中的CPM信号正交分解方框图,可以采用Gram-Schmidt正交分解法得到正交基,在我们实施例中正交基的数目为2个,因此,采用两路相关的方式进行正交分解,得到两路向量。将接收到的CPM信号分解之后,消除CPM信号自身的相关性,可以更好的进行信道均衡。
所述的接收衰落后的连续相位调制信号,移除循环前缀,得到经历衰落的训练字符号块和数据符号块,用存储于本地的正交基对每个训练字符号块或数据符号块做正交分解的步骤为:
1)接收经历了频率选择性衰落信道的连续相位调制信号r(t),对r(t)先移除循环前缀;
2)加上倾斜相位使相位网格状态时不变,得到r′(t);
3)存储于本地的正交基为f0(t)=1, f 1 ( t ) = 1 1 - 4 π 2 ( e jπt T - j 2 π ) , 0≤t<T,其中T为符号周期;将加入了倾斜相位的接收信号r′(t)在符号周期内与两个正交基作相关,得到两组向量r0(n)和r1(n)。
图5示出了CPM信号频域均衡的过程。CPM信号正交分解之后,得到两路向量r0(n)和r1(n),将这两路向量变换到频域(501、504),根据信道估计的结果和MMSE准则,在频域上均衡消除信道的影响(502、505),将均衡后的频域数据变换到时域进行解调(503、506),完成频域均衡。
所述的正交分解之后,把训练字符号块变换到频域做信道估计,数据符号块则根据信道估计结果,在频域上做均衡,去除信道干扰的影响后再变换回时域的步骤为:
1)正交分解将训练字符号块和数据符号块同时分解为两路向量;在两路向量中分别提取出训练字符号块的向量段r0 T(n)、r1 T(n);
2)将r0 T(n)、r1 T(n)和存储于本地的训练字符号块向量s0 T(n)、s1 T(n)一起变换到频域得到R0 T(k)、R1 T(k)和S0 T(k)、S1 T(k);
3)根据最小均方误差准则,得到正交分解后信道频域的估计值为:
H ~ 0 = H 0 * | H 0 | 2 + N 0 , H ~ 1 = H 1 * | H 1 | 2 + N 0 ( H 0 = R 0 T / S 0 T , H 1 = R 1 T / S 1 T );
4)将数据符号块正交分解后得到的向量r0(n)和r1(n)变换到频域上得到R0,R1,根据信道估计结果,在频域去除信道的影响,得到两路均衡后的输出
Figure G2009101545692D00076
S 0 ~ = H 0 ~ * R 0 , S 1 ~ = H 1 ~ * R 1 ;
5)将
Figure G2009101545692D000710
变换到时域,送进解调器解调。
本发明是一个CPM调制使用频域均衡方法对抗频率选择性衰落信道影响的***。这里所述的方法与线路,对于CPM调制解调单体部件可以是完全传统的,我们要求将它与频域信道均衡的结合作为发明进行保护,频域信道均衡包括CPM信号正交分解、频域上的信道估计与频域均衡。以上所述仅为特定应用场合的具体实施方式,但本发明的真实精神和范围不局限于此,任何熟悉本领域的技术人员可以修改单体部件的具体方法,实现不同应用场合的CPM频域均衡调制解调***。本发明仅由后附权利要求书及其等效技术方案来限定,我们要求将这些作为本发明来保护。

Claims (4)

1.一种连续相位调制频域均衡调制解调方法,其特征在于包括如下步骤:
1)从数据源获得二进制序列,将二进制序列映射成调制符号序列,将调制符号和训练字符号组成符号块,按照符号块长度逐块送入连续相位调制器;
2)每个符号块在完成连续相位调制之后,***循环前缀波形,得到连续相位调制信号;
3)将连续相位调制信号在通信通道上发送,连续相位调制信号经历频率选择性衰落;
4)接收衰落后的连续相位调制信号,移除循环前缀,得到经历衰落的训练字符号块和数据符号块,用存储于本地的正交基对每个训练字符号块和数据符号块做正交分解;
5)正交分解之后,把训练字符号块变换到频域做信道估计,数据符号块则根据信道估计结果,在频域上做均衡,去除信道干扰的影响后再变换回时域;
6)均衡后变换回时域的数据符号块送进维特比译码器解调得到源序列。
2.根据权利要求1所述的一种连续相位调制频域均衡调制解调方法,其特征在于,所述的从数据源获得二进制序列,将二进制序列映射成调制符号序列,将调制符号和训练字符号组成符号块,按照符号块长度逐块送入连续相位调制器步骤为:
1)将从数据源获得的二进制序列an映射为数据符号序列In,具体是将1映射为+1,0映射为-1;
2)每252个数据符号In和4个归零符号组成一个数据符号块,数据符号块的长度为256;
3)每隔20个数据符号块***一个训练字符号块,以训练字符号块为起点组成帧结构,训练字符号块仍然是由长度为252的训练字符号和4个归零符号组成;
4)每个符号块依次进入连续相位调制器,在每个符号块的起点相位为0,同时每个符号块的第231个符号处相位也为0,连续相位调制信号波形表达式为
Figure FSB00000639035300011
其中
Figure FSB00000639035300012
为成形函数,In为数据符号块和训练序列符号块,
Figure FSB00000639035300014
h=0.5,T为符号周期。
3.根据权利要求1所述的一种连续相位调制频域均衡调制解调方法,其特征在于,所述的每个符号块在完成连续相位调制之后,***循环前缀波形,得到连续相位调制信号步骤为:
1)循环前缀长度的大小要大于信道最大的弥散长度,设定频率选择性衰落信道的弥散长度最大为25个符号,故循环前缀的长度设定为25个符号;
2)在第229个符号处将相位归零,这将用掉2个相位归零符号;
3)在每个符号块的末尾,需要加上另外2个归零符号,使相位从任意状态归0;
4)对于每个符号块,将最后25个符号的连续相位调制信号波形复制到符号块最前面形成循环前缀。
4.根据权利要求1所述的一种连续相位调制频域均衡调制解调方法,其特征在于,所述的将连续相位调制信号在通信通道上发送,连续相位调制信号经历频率选择性衰落步骤为:连续相位调制信号
Figure FSB00000639035300021
经历频率选择性衰落信号后表示为:
Figure FSB00000639035300022
其中z(t)是零均值复高斯噪声,方差为N0;MD是信道的最大弥散长度,hl为信道第L径的增益。
CN2009101545692A 2009-11-12 2009-11-12 连续相位调制频域均衡调制解调方法 Expired - Fee Related CN101707579B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101545692A CN101707579B (zh) 2009-11-12 2009-11-12 连续相位调制频域均衡调制解调方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101545692A CN101707579B (zh) 2009-11-12 2009-11-12 连续相位调制频域均衡调制解调方法

Publications (2)

Publication Number Publication Date
CN101707579A CN101707579A (zh) 2010-05-12
CN101707579B true CN101707579B (zh) 2012-05-30

Family

ID=42377772

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101545692A Expired - Fee Related CN101707579B (zh) 2009-11-12 2009-11-12 连续相位调制频域均衡调制解调方法

Country Status (1)

Country Link
CN (1) CN101707579B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107026808B (zh) * 2017-03-23 2019-10-11 西安电子科技大学 基于混沌交织的多进制连续相位调制频域均衡方法
CN108156106B (zh) * 2017-12-19 2020-02-18 浙江大学 适用于可见光通信的连续相位调制信号的传输与重建方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1113058A (zh) * 1993-12-21 1995-12-06 日本电气株式会社 用于移动卫星通信终端的发射机装置
CN101262240A (zh) * 2008-04-25 2008-09-10 浙江大学 一种易于硬件实现的全数字频率变换的方法及其装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1113058A (zh) * 1993-12-21 1995-12-06 日本电气株式会社 用于移动卫星通信终端的发射机装置
CN101262240A (zh) * 2008-04-25 2008-09-10 浙江大学 一种易于硬件实现的全数字频率变换的方法及其装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈 晨等.单载波频域均衡技术的定时同步研究.<<浙 江 大 学 学 报(工学版)>>.2007, *

Also Published As

Publication number Publication date
CN101707579A (zh) 2010-05-12

Similar Documents

Publication Publication Date Title
CN106911374B (zh) 一种低复杂度软输出空间调制检测方法
CN106941465B (zh) 超奈奎斯特率块的传输方法、发射机、接收机及***
CN101662434B (zh) 利用导频格式设计的WiMAX信道估计方法
CN105814855A (zh) 超Nyquist发送***中的预编码
WO2006125343A1 (fr) Emetteur, recepteur et procedes de multiplexage orthogonal a repartition frequentielle et temporelle
CN110011947A (zh) 一种超奈奎斯特传输***中基于分解矩阵的干扰消除调制方法
CN113381951A (zh) 时变频选衰落信道下mftn联合信道估计与均衡方法
KR20070089558A (ko) 직교 주파수 분할 다중화 기반의 이동 통신 시스템에서다이버시티 이득을 향상시키기 위한 선 부호화 방법 및이를 이용한 송신 장치 및 방법
CN110290083B (zh) 基于低秩矩阵估计的多载波差分混沌***解调方法
CN107332606A (zh) 基于双采样的leo***差分空时正交频分复用编码方法
US8189709B2 (en) Frequency domain equalization method for continuous phase modulated signals
CN101707579B (zh) 连续相位调制频域均衡调制解调方法
CN101325438B (zh) 基于gmc***的4天线分集发送装置、发送方法及接收装置、接收方法
CN101895492A (zh) 一种单载波频域均衡技术的过采样接收方法
CN111314027A (zh) 一种双向中继信道物理层网络编码的符号检测方法
CN103259758A (zh) 一种基于压缩扩展的ofdm通信***及方法
Mohr et al. Forney observation models for faster-than-Nyquist signaling on nonlinear satellite links
CN109639301A (zh) 一种基于置信度估计的ftn均衡方法
CN104135455A (zh) 一种通信***迭代接收方法
CN105187354A (zh) 一种基于pts技术抑制ofdm通信信号峰平比的方法
CN102404274A (zh) 降低ofdm信号峰平比的双曲正切压扩变换方法
CN101945062B (zh) 一种单载波无线通信***中的接收处理方法
CN116668247B (zh) 一种超奈奎斯特***Cholesky预编码方法
Valluri et al. Receiver design for UW-GFDM systems
CN104980321A (zh) 一种基于判决反馈的ofdm信道质量估计算法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Zhao Minjian

Inventor after: Zhong Jie

Inventor after: Fu Tao

Inventor after: Lv Chunmiao

Inventor after: Peng Cen

Inventor after: Zheng Yamin

Inventor before: Zhao Minjian

Inventor before: Fu Tao

Inventor before: Lv Chunmiao

Inventor before: Peng Cen

Inventor before: Zheng Yamin

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: ZHAO MINJIAN FU TAO LV CHUNMIAO PENG CEN ZHENG YAMIN TO: ZHAO MINJIAN ZHONG JIE FU TAO LV CHUNMIAO PENG CEN ZHENG YAMIN

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120530

Termination date: 20141112

EXPY Termination of patent right or utility model