CN107612658B - 一种基于b类构造格型码的高效编码调制与译码方法 - Google Patents
一种基于b类构造格型码的高效编码调制与译码方法 Download PDFInfo
- Publication number
- CN107612658B CN107612658B CN201710979281.3A CN201710979281A CN107612658B CN 107612658 B CN107612658 B CN 107612658B CN 201710979281 A CN201710979281 A CN 201710979281A CN 107612658 B CN107612658 B CN 107612658B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- vector
- bit
- decoding
- complex signal
- 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
Images
Landscapes
- Error Detection And Correction (AREA)
Abstract
本发明公开了一种基于B类构造格型码的高效编码调制与译码方法,涉及格型码的编译码技术领域;其包括如下步骤:1)在编码调制端,构建基于B类构造格型码的消息空间,在消息空间将给定的n+k‑1维q元原始消息划分后进行编码调制得到n维复信号x,将n维复信号x传输至译码端得到待处理的n维复信号y;2)在译码端,将待处理的n维复信号y进行逐位量化实现分块后并将其进行分级译码得到k维q元向量和n‑1维q元向量;3)将k维q元向量和n‑1维q元向量拼接得到n+k‑1维q元消息完成译码;本发明解决了现有B类构造格型码无法高效实施提高编码效率的问题,达到了高效地实现基于B类构造格型消息空间的信息传输的效果。
Description
技术领域
本发明涉及格型码的编译码技术领域,尤其是一种基于B类构造格型码的高效编码调制与译码方法。
背景技术
格型码是一类可无限逼近高斯信道理论容量log2(1+SNR)的编码结构,有着良好的代数结构;格型码可应用在物理层网络编码中,尤其是双向中继系统和多用户中继传输系统等,所传输的消息都能调制映射到格型空间中的一个元素,用户将调制后的n维复数信号同时传送至中继节点,中继节点通过所接收的由多个用户传输信号自然线性叠加而成的混合信号,直接计算出一组由多用户原始消息线性组合而成的新消息,并将新消息基于格型码重新调制再广播出去;除了基线格型码结构,目前经典的格型码构造方法有A类构造(Construction A)、B类构造(Construction B)、D类构造(Construction D)等等。A类构造格型的结构简单,跟纠错码有着一一映射的关系,因此对要传输的消息进行调制比较容易,对格型空间中的任一点找到其对应的陪集代表也易于实现;基于B类构造方法的格型结构比较复杂,不易于实现,对其的译码研究比较少;目前格型码主要有三种:整数Z上的格型划分、高斯整数Z[i]上的格型划分、爱森斯坦整数Z[ω]上的格型划分,其中高斯整数爱森斯坦整数由于Z[ω]有更好的代数结构,相比Z[i]有0.167dB的形状增益(shaping gain)。现有对A类构造格型研究较多,而B类构造格型除自身具有调制编码结合和良好的代数结构的特点外,与A类构造格型相比,在编码增益相同的情况下可以提供更高的编码速率,所以需要一种编码调制与译码方法利用B类构造格型实现高效调制、编码和译码。
发明内容
本发明的目的在于:本发明的目的在于提供了一种基于B类构造格型码的高效编码调制与译码方法,解决了现有B类构造方法由于格型结构复杂无法高效实施提高传输编码效率的问题。
本发明采用的技术方案如下:
一种基于B类构造格型码的高效编码调制与译码方法,包括如下步骤:
步骤1:在编码调制端,构建基于B类构造格型码的消息空间,在消息空间将给定的n+k-1维q元原始消息划分后进行编码调制得到n维复信号x,将n维复信号x传输至译码端得到待处理的n维复信号y;
步骤2:在译码端,将待处理的n维复信号y进行逐位量化实现分块后并进行分级译码得到k维q元向量和n-1维q元向量;
步骤3:将k维q元向量和n-1维q元向量拼接得到n+k-1维q元消息完成译码。
优选地,所述步骤1包括如下步骤:
步骤1.1:根据给定模值为q的爱森斯坦素数β得到伽罗华域GF(q)=Z[ω]/βZ[ω],其中GF(q)中元素均等价于Z[ω]/βZ[ω]中陪集;由伽罗华域GF(q)下(n,k)纠错码C基于B类构造生成的格型码Λ/Λ′为包含qn+k-1个消息的消息空间,每个消息均取值为伽罗华域GF(q)=Z[ω]/βZ[ω];根据伽罗华域GF(q)下(n,k)纠错码C得到其k×n维q元生成矩阵G;
步骤1.2:给定n+k-1维q元原始消息m=(m1,m2,...,mn+k-1)并从中选择前k个q元符号组成向量(m1,m2,...,mk),将其乘以k×n维生成矩阵G得到n维q元向量c1,将该向量中n个GF(q)元素逐位映射至其所对应的陪集代表,形成n维复数向量x1=(x11,x12,...,x1n),其中每个元素属于Z[ω];
步骤1.3:根据B类构造格型的定义得到(n-1)×n维生成矩阵M=[In–1 1](其中In–1表示(n-1)×(n-1)的单位阵、1表示取值均为1的n-1维列向量);
步骤1.4:从n+k-1维q元原始消息m=(m1,m2,...,mn+k-1)中选择后n-1个q元符号组成向量(mk+1,mk+2,...,mk+n-1),将其乘以(n-1)×n维生成矩阵M得到n维q元向量c2,将该向量中n个GF(q)元素逐位映射至其所对应的陪集代表后乘以爱森斯坦素数β,形成n维复数向量x2=(x21,x22,...,x2n),其中每个元素属于Z[ω];
步骤1.5:将步骤1.2中n维复数向量x1=(x11,x12,...,x1n)与步骤1.4中n维复数向量x2=(x21,x22,...,x2n)相加得到(x11+x21,x12+x22,...,x1n+x2n)后将其中每个元素进行模值缩减得到n维复数向量(x1,x2,...,xn)即n维复信号x,该向量为n+k-1维q元原始消息(m1,m2,...,mn+k-1)于消息空间中所对应的陪集代表,将n维复信号x传输至译码端得到待处理的n维复信号y。
优选地,所述步骤2包括如下步骤:
步骤2.1:在译码端,对待处理的n维复信号y=(y1,y2,...,yn)进行逐位量化,并将量化所得爱森斯坦整数映射至其于伽罗华域GF(q)=Z[ω]/βZ[ω]中所属的陪集,得到n维q元向量(a1,a2,...,an);
步骤2.2:基于(n,k)纠错码C的译码方法,将步骤2.1中所得n维q元向量(a1,a2,...,an)映射到(n,k)纠错码C中的一个码字得到n维q元向量u;将n维q元向量u中n个GF(q)元素逐位映射至其所对应的陪集代表,形成n维复数向量w;
步骤2.3:将待处理的n维复信号y减去步骤2.2所得n维复数向量w再逐位除以爱森斯坦素数β,得到新的n维待处理复信号y′=(y-w)/β;
步骤2.4:将n维待处理复信号y′=(y1′,y2′,...,yn′)进行逐位量化得到n维复向量w′=(w1′,w2′,...,wn′),并将量化所得爱森斯坦整数wj′映射至其于伽罗华域GF(q)=Z[ω]/βZ[ω]中所属的陪集得到n维q元向量v=(v1,v2,...,vn);
步骤2.5:对n维q元向量v=(v1,v2,...,vn),计算其中各元素之和若结果不等于GF(q)中0元素,则逐位计算y′与w′间元素的欧氏距离,找出欧氏距离最大所对应的元素位置i,即并将vi的取值重新定义为∑1≤j≤n,j≠ivj,完成更新n维q元向量v;
优选地,所述步骤3包括如下步骤:
综上所述,由于采用了上述技术方案,本发明的有益效果是:
1.本发明基于B类构造格型码构建消息空间,在消息空间将给定的n+k-1维q元原始消息划分后,进行编码调制得到n维复信号x进行传输可以高效的将原始n+k-1维原始消息编码调制到B类构造格型消息空间中;通过将接收复信号分块后分级译码可以高效地将复向量译码还原至原始n+k-1维消息,从而高效的实现基于B类构造格型消息空间的信息传输;
2.本发明针对(n,k)纠错码编码后较复杂的B类构造格型进行分块,使复杂的格型变得有规律可循,有利于加快译码速度;
3.本发明通过陪集代表能直接对复杂的B类构造格型构建出消息空间,且能直接将调制和编码结合,简化流程,进一步实现高效调制和编码;
4.本发明通过分级译码来减少时间复杂度,先对格型进行量化分块,加快量化和映射,通过第一次译码将其确定为具体哪一块,再在该块里进行第二次译码,降低时间复杂度,快速地对接收消息进行译码,有利于实现高效译码。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1是本发明译码流程图;
图2是本发明编码调制流程图。
具体实施方式
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
下面结合图1-2对本发明作详细说明。
实施例1
一种基于B类构造格型码的高效编码调制与译码方法,包括如下步骤:
步骤1:在编码调制端,构建基于B类构造格型码的消息空间,在消息空间将给定的n+k-1维q元原始消息m=(m1,m2,...,mn+k-1)划分,分别乘以根据(n,k)纠错码C得到k×n维生成矩阵G和根据B类构造格型的定义得到的(n-1)×n维生成矩阵M后再乘爱森斯坦素数β得到对应的向量,将向量映射至消息空间中得到对应的qn+k-1个陪集代表,并将qn+k-1陪集代表作为n维复信号x传输至译码端得到待处理的n维复信号y;
步骤1包括如下步骤:
步骤1.1:给定模值为q的爱森斯坦素数β,根据爱森斯坦素数β得到伽罗华域GF(q)=Z[ω]/βZ[ω],其中GF(q)中元素均等价于Z[ω]/βZ[ω]中陪集;给定任意基于伽罗华域GF(q)=Z[ω]/βZ[ω]的(n,k)纠错码C;构建消息空间Λ/Λ′(其中Λ为由(n,k)纠错码C基于B类构造生成的n维Z[ω]-格型,Λ′=(β2Z[ω])n),该消息空间Λ/Λ′共包含qn+k-1个消息;给定任意n+k-1维q元原始消息m=(m1,m2,...,mn+k-1),其中每个符号mj均取值为伽罗华域GF(q)=Z[ω]/βZ[ω];
步骤1.2:根据(n,k)纠错码C得到其k×n维q元生成矩阵G;
步骤1.3:从n+k-1维q元原始消息m=(m1,m2,...,mn+k-1)中选择前k个q元符号组成向量(m1,m2,...,mk),将其乘以k×n维生成矩阵G得到n维q元向量c1,将该向量中n个GF(q)元素逐位映射至其所对应的陪集代表,形成n维复数向量x1=(x11,x12,...,x1n),其中每个元素属于Z[ω];
步骤1.4:根据B类构造格型的定义得到一个(n-1)×n维生成矩阵M=[In–1 1](其中In–1表示(n-1)×(n-1)的单位阵、1表示取值均为1的n-1维列向量);
步骤1.5:从n+k-1维q元原始消息m=(m1,m2,...,mn+k-1)中选择后n-1个q元符号组成向量(mk+1,mk+2,...,mk+n-1),将其乘以(n-1)×n维生成矩阵M得到n维q元向量c2,将该向量中n个GF(q)元素逐位映射至其所对应的陪集代表后乘以爱森斯坦素数β,形成n维复数向量x2=(x21,x22,...,x2n),其中每个元素属于Z[ω];
步骤1.6:将步骤1.3中n维复数向量x1=(x11,x12,...,x1n)与步骤1.5中n维复数向量x2=(x21,x22,...,x2n)相加得到(x11+x21,x12+x22,...,x1n+x2n),对每个元素x1j+x2j按如下方式进行模值缩减:首先将x1j+x2j表示成k1+k2ω的形式,其中k1、k2均为整数。若k1>q/2,令k1′=k1-q,若k1<-q/2,令k1′=k1+q;若k2>q/2,则令k2′=k2-q;若k1<-q/2,则令k2′=k2+q;若k1 2+k2 2-k1k2<k1′2+k2′2-k1′k2′,则设xj=k1+k2ω,否则设xj=k1′+k2′ω,从而生成(x1,x2,...,xn)即n维复信号x,该向量为n+k-1维q元原始消息(m1,m2,...,mn+k-1)于消息空间Λ/Λ′中所对应的陪集代表,将qn+k-1陪集代表作为n维复信号x传输至译码端得到待处理的n维复信号y。
步骤2:在译码端,对待处理的n维复信号y进行第一次译码得到n维q元向量u,并将n维q元向量u中n个GF(q)元素逐位映射至其所对应的陪集代表得到n维复数向量w;
步骤2包括如下步骤:
步骤2.1:在译码端,对待处理的n维复信号y=(y1,y2,...,yn)进行逐位量化,即对每一元素yj,计算与其欧氏距离最近的爱森斯坦整数,以及将该爱森斯坦整数映射至其于伽罗华域GF(q)=Z[ω]/βZ[ω]中所属的陪集,记为aj,从而形成n维q元向量(a1,a2,...,an);
步骤2.2:将步骤2.1中所得n维q元向量(a1,a2,...,an)映射到(n,k)纠错码C中的一个码字得到n维q元向量u;
步骤2.3:将步骤2.2中n维q元向量u中n个GF(q)元素逐位映射至其所对应的陪集代表,形成n维复数向量w。
步骤3:在译码端,将n维复信号y减去步骤2中n维复数向量w后逐位除爱森斯坦素数β得到n维待处理复信号y′,对y′进行第二次译码得到n维q元向量v;
步骤3包括如下步骤:
步骤3.1:在译码端,将n维复信号y减去步骤2.3所得n维复数向量w再逐位除以爱森斯坦素数β,得到新的n维待处理复信号y′=(y-w)/β;
步骤3.2:将n维待处理复信号y′=(y1′,y2′,...,yn′)进行逐位量化,即对每一元素yj′,计算与其欧氏距离最近的爱森斯坦整数wj′以及该爱森斯坦整数wj′在GF(q)=Z[ω]/βZ[ω]中所属的陪集,记为aj′,从而形成n维复向量w′=(w1′,w2′,...,wn′)以及n维q元向量v=(v1,v2,...,vn);
步骤3.3:对n维q元向量v=(v1,v2,...,vn),计算其中各元素之和若结果不等于GF(q)中0元素,则逐位计算y′与w′间元素的欧氏距离,找出欧氏距离最大所对应的元素位置i,即并将vi的取值重新定义为∑1≤j≤n,j≠ivj,完成更新n维q元向量v。
步骤4:在译码端,将步骤2中n维q元向量u乘以k×n维生成矩阵G对应的n×k维q元矩阵H得到k维q元向量将步骤3所得n维q元向量v删除最后一位元素,得到n-1维q元向量将k维q元向量和n-1维q元向量拼接得到n+k-1维q元消息完成译码;
步骤4包括如下步骤:
步骤4.1:基于1.2中(n,k)纠错码C的q元k×n维生成矩阵G,计算n×k维q元矩阵H(其中满足G乘以H等于k×k维的单位阵);
本发明基于B类构造格型码构建消息空间,在消息空间将给定的n+k-1维q元原始消息划分后,进行编码调制得到n维复信号x进行传输可以高效的将原始n+k-1维原始消息编码调制到B类构造格型消息空间中;通过将接收复信号分块后分级译码可以高效地将复向量译码还原至原始n+k-1维消息,从而高效的实现基于B类构造格型消息空间的信息传输。
Claims (1)
1.一种基于B类构造格型码的高效编码调制与译码方法,其特征在于:包括如下步骤:
步骤1:在编码调制端,构建基于B类构造格型码的消息空间,在消息空间将给定的n+k-1维q元原始消息划分后进行编码调制得到n维复信号x,将n维复信号x传输至译码端得到待处理的n维复信号y;
步骤2:在译码端,将待处理的n维复信号y进行逐位量化实现分块后并进行分级译码得到k维q元向量和n-1维q元向量;
步骤3:将k维q元向量和n-1维q元向量拼接得到n+k-1维q元消息完成译码;
所述步骤1包括如下步骤:
步骤1.1:根据给定模值为q的爱森斯坦素数β得到伽罗华域,其中中元素均等价于中陪集;由伽罗华域下(n,k)纠错码C基于B类构造生成的格型码为包含个消息的消息空间,每个消息均取值为伽罗华域;根据伽罗华域下(n,k)纠错码C得到其k×n维q元生成矩阵G;
步骤1.2:给定n+k-1维q元原始消息并从中选择前k个q元符号组成向量,将其乘以k×n维生成矩阵G得到n维q元向量,将该向量中n个元素逐位映射至其所对应的陪集代表,形成n维复数向量,其中每个元素属于;
步骤1.4:从n+k-1维q元原始消息中选择后n-1个q元符号组成向量,将其乘以(n-1)×n维生成矩阵M得到n维q元向量,将该向量中n个元素逐位映射至其所对应的陪集代表后乘以爱森斯坦素数β,形成n维复数向量 ,其中每个元素属于;
步骤1.5:将步骤1.2中n维复数向量与步骤1.4中n维复数向量 相加得到后将其中每个元素进行模值缩减得到n维复数向量即n维复信号x,该向量为n+k-1维q元原始消息于消息空间中所对应的陪集代表,将n维复信号x传输至译码端得到待处理的n维复信号y;
所述步骤2包括如下步骤:
步骤2.2:基于(n,k)纠错码C的译码方法,将步骤2.1中所得n维q元向量映射到(n,k)纠错码C中的一个码字得到n维q元向量u;将n维q元向量u中n个元素逐位映射至其所对应的陪集代表,形成n维复数向量w;
所述步骤3包括如下步骤:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710979281.3A CN107612658B (zh) | 2017-10-19 | 2017-10-19 | 一种基于b类构造格型码的高效编码调制与译码方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710979281.3A CN107612658B (zh) | 2017-10-19 | 2017-10-19 | 一种基于b类构造格型码的高效编码调制与译码方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107612658A CN107612658A (zh) | 2018-01-19 |
CN107612658B true CN107612658B (zh) | 2020-07-17 |
Family
ID=61078713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710979281.3A Active CN107612658B (zh) | 2017-10-19 | 2017-10-19 | 一种基于b类构造格型码的高效编码调制与译码方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107612658B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101082592B1 (ko) * | 2010-07-20 | 2011-11-10 | 한국해양연구원 | 이진 부호를 사용하는 물리계층 네트워크 코딩 방법 |
CN102918590A (zh) * | 2010-03-31 | 2013-02-06 | 韩国电子通信研究院 | 编码方法和装置、以及解码方法和装置 |
CN102957435A (zh) * | 2011-08-17 | 2013-03-06 | 国家广播电影电视总局广播科学研究院 | Ldpc码校验矩阵的构造方法、装置和编码方法及系统 |
CN105846961A (zh) * | 2016-05-06 | 2016-08-10 | 电子科技大学 | 一种dmr协议网格码快速译码方法及译码装置 |
CN107070586A (zh) * | 2016-12-12 | 2017-08-18 | 北京邮电大学 | 一种基于低密度格码的中继协作编码方法及系统 |
-
2017
- 2017-10-19 CN CN201710979281.3A patent/CN107612658B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102918590A (zh) * | 2010-03-31 | 2013-02-06 | 韩国电子通信研究院 | 编码方法和装置、以及解码方法和装置 |
KR101082592B1 (ko) * | 2010-07-20 | 2011-11-10 | 한국해양연구원 | 이진 부호를 사용하는 물리계층 네트워크 코딩 방법 |
CN102957435A (zh) * | 2011-08-17 | 2013-03-06 | 国家广播电影电视总局广播科学研究院 | Ldpc码校验矩阵的构造方法、装置和编码方法及系统 |
CN105846961A (zh) * | 2016-05-06 | 2016-08-10 | 电子科技大学 | 一种dmr协议网格码快速译码方法及译码装置 |
CN107070586A (zh) * | 2016-12-12 | 2017-08-18 | 北京邮电大学 | 一种基于低密度格码的中继协作编码方法及系统 |
Non-Patent Citations (2)
Title |
---|
Lattice Network Codes Based on Eisenstein Integers;Qifu Sun, Jinhong Yuan,Tao Huang,Kenneth W.Shum;《IEEE TRANSACTIONS ON COMMUNICATIONS 》;20130707;全文 * |
基于有限域的有结构低密度格码的构造;李家云;《CNKI》;20150430;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107612658A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100808664B1 (ko) | 패리티 검사행렬 저장 방법 및 이를 이용한 블록 저밀도패리티 검사 부호화 방법 및 장치 | |
CN101459430B (zh) | 低密度生成矩阵码的编码方法及装置 | |
Böcherer et al. | Matching dyadic distributions to channels | |
US7743309B2 (en) | Distributed source coding with context | |
RU2450442C2 (ru) | Способ и устройство для кодирования и декодирования канала в системе связи с использованием кодов с низкой плотностью проверок на четность | |
KR101496182B1 (ko) | 최소거리가 확장된 극 부호 생성 방법 및 장치 | |
CN105656604A (zh) | 一种比特交织极化编码调制方法及装置 | |
US10826648B2 (en) | Methods of converting or reconverting a data signal and method and system for data transmission and/or data reception | |
CN108964669B (zh) | 基于度分解和交替乘子法的ldpc码二次规划译码方法 | |
US11038525B2 (en) | Methods of converting or reconverting a data signal and method and system for data transmission and/or data reception | |
CN107124251B (zh) | 一种基于任意内核的极化码编码方法 | |
CN111050291B (zh) | 一种基于系统在线喷泉码的高编码符号度值广播方法 | |
US10523480B1 (en) | K-bit enumerative sphere shaping of multidimensional constellations | |
Shi et al. | Zigzag decodable online fountain codes with high intermediate symbol recovery rates | |
Nguyen et al. | Integrating sparsity into Fulcrum codes: Investigating throughput, complexity and overhead | |
CN111464300A (zh) | 一种适用于连续变量量子密钥分发的高速后处理方法 | |
US20100031116A1 (en) | Method for encoding low density parity check codes using result of checking previously specified parity bits | |
CN103368586B (zh) | 面向深空探测多媒体业务的独立窗不等保护喷泉编码方法 | |
CN107612658B (zh) | 一种基于b类构造格型码的高效编码调制与译码方法 | |
CN110601699B (zh) | 码率动态可变的多元ldpc码实现方法 | |
JP6744982B2 (ja) | OvXDMシステムに適用される一種類のファストデコード方法、装置及びOvXDMシステム | |
CN101854179A (zh) | 一种应用于ldpc译码的5比特量化方法 | |
EP3767457A1 (en) | Data communication | |
CN107592958B (zh) | 不同码长下使用一个多比率指数表和一个扩展因子的多比率的ldpc编码和解码 | |
MX2014012121A (es) | Codificador de revision de paridad de baja densidad que tiene una longitud de 64800 y un indice de codigo de 7/15 y metodo de codificacion de revision de paridad de baja densidad que utiliza el mismo. |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211108 Address after: 410000 Room 401, building 5, phase II, Depu Wuhe Enterprise Park, No. 1 and 3, luolongtang Road, Changsha Economic and Technological Development Zone, Hunan Province Patentee after: Hunan Wanban Shangpin Information Technology Co.,Ltd. Address before: No. 30 Xueyuan Road, Haidian District, Beijing 100083 Patentee before: University OF SCIENCE AND TECHNOLOGY BEIJING |