CN103208996A - 准循环码的频域编码方法 - Google Patents
准循环码的频域编码方法 Download PDFInfo
- Publication number
- CN103208996A CN103208996A CN2013101340097A CN201310134009A CN103208996A CN 103208996 A CN103208996 A CN 103208996A CN 2013101340097 A CN2013101340097 A CN 2013101340097A CN 201310134009 A CN201310134009 A CN 201310134009A CN 103208996 A CN103208996 A CN 103208996A
- Authority
- CN
- China
- Prior art keywords
- matrix
- row
- quasi
- coding
- expression
- 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.)
- Granted
Links
Images
Abstract
Description
技术领域
本发明涉及一种准循环码的频域编码方法,属于信息论及编码技术领域。
背景技术
1、现有技术的技术方案。
假设准循环码的生成矩阵G=[I|P],其中,G是ek×en的准循环矩阵,e表示准循环阵中每个小循环方阵的行数或者列数,每个准循环矩阵由k×n个小循环方阵构成,k表示准循环矩阵的行数,n表示准循环矩阵的列数,G由I和P两部分组成。I是ek×ek的单位阵,P是ke×(n-k)e的准循环阵。要编码的消息为消息序列m=[m0,m1,…,mi,…,mk],其中mi是一个长为e的向量。编码后的码字为C=[C0,C1,…,Ci,…,Cn],其中Ci是一个长为e的向量。对于二进制准循环码,生成矩阵G,消息序列m以及码字C中的元素均为二元域GF(2)中的元素。对于非二进制准循环码,生成矩阵G,消息序列m以及码字C中的元素均为GF(2r)中的元素,其中GF(2r)表示二元域GF(2)的r次方扩域。其传统的编码方式均为:
C=m×G
2、现有技术的缺点
现有准循环码编码技术编码的计算复杂度很高,对于二进制准循环码以及非二进制准循环码而言,其编码复杂度均为O(k(n-k)e2),即其编码的复杂度在k(n-k)e2的数量级上。
发明内容
本发明的目的是为了解决在常规准循环码编码过程较复杂的缺陷,提出了一种能大幅度降低编码运算复杂度的准循环码的频域编码方法。
一种准循环码的频域编码方法,包括以下几个步骤:
步骤一、判断准循环码的类别;
步骤二、针对非二进制准循环码进行编码;
步骤三、针对二进制准循环码进行编码。
本发明的优点为,频域编码相比于常规编码方式复杂度大大降低,具体表现如下:
(1)对于非二进制准循环码,常规的编码方式将长度为ke的消息编码成长度为ne的码字运算复杂度为O(k(n-k)e2),而对于频域编码其运算复杂度为O(ek(n-k))。
(2)对于二进制准循环码,常规的编码方式将长度为ke的消息编码成长度为ne的码字运算复杂度为O(k(n-k)e2),而对于频域编码其运算复杂度为O(ek(n-k)log2e)。
附图说明
图1是本发明的方法流程图;
图2是本发明非二进制准循环码编码方法流程图;
图3是本发明二进制准循环码编码方法流程图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种准循环码的频域编码方法,流程如图1所示,包括以下几个步骤:
步骤一、判断准循环码的类别。
准循环码的编码过程是将一段消息序列与准循环码生成矩阵相乘最终变为一段信息序列,设定消息序列为m,其中,m=[m0,m1,...,mi,...,me]。
设n、e、k为正整数,ke为编码前消息序列m的长度,ne为编码后信息序列的长度,其中ke<ne。
(1)当(ne,ke)线性分组码Cqc中元素在GF(2)中,符合以下条件:
1)Cqc中的码字由n组e个比特组成。
2)对Cqc中的码字的n组比特分别进行循环移位得到的码字仍然是Cqc中的码字。
则此线性分组码Cqc为二进制准循环码,转入步骤三。
(2)当(ne,ke)线性分组码Cqc中的元素在GF(2r)中,且满足上述1)与2)的条件,则此线性分组码Cqc为非二进制准循环码,转入步骤二。
步骤二、针对非二进制准循环码进行编码。
具体包括:
对非二进制准循环码,假设准循环码的生成矩阵G(ek×en)是由k×n组的非二进制准循环矩阵组成,其中e+1为2的r次幂,即e=2r-1。其中,e表示准循环阵中每个小循环方阵的行数或者列数,每个准循环矩阵由k×n个小循环方阵构成,k表示准循环矩阵的行数,n表示准循环矩阵的列数。
假设G=[Wi,j],其中Wi,j=circ(ai,j,0,ai,j,1,…,ai,j,e-1)是一个e×e的准循环矩阵,且第一行为ai,j,0,ai,j,1,…,ai,j,e-1,其中0≤i<k,0≤j<n,Wi,j表示G矩阵中位于第i行j列的小矩阵,circ(ai,j,0,ai,j,1,…,ai,j,e-1)表示Wi,j中第一行的元素分别为ai,j,0,ai,j,1,…,ai,j,e-1。
其中:V=[α-ij],V-1=[αij],α为GF(2r)中的本原元,[αij]表示第i行j列的元素为αij的矩阵,diag表示对角矩阵。
将中第0行至第ek-1行中的元素分别置于第πrow(0)行至第πrow(ek-1)行中,然后再将第0列至第en-1列中的元素分别置于第πcol(0)列至第πcol(en-1)列中,得到矩阵其中Dt=[di,j,t],di,j,t为矩阵Dt中的第i行j列的元素,对于0≤i<k,0≤j<n,0≤t<e,有
需要编码的消息序列m=[m0,m1,…,mi,…,me],其中mi是一个长为k的子序列,即mi=[mi,t],mi,t表示m中第i个子序列mi的第t个元素,且0≤t<k。
其中,编码后的码字为C=[C0,C1,…,Ci,…,Ce],其中Ci是一个长为n的子序列,即Ci=[ci,t],ci,t表示C中第i个子序列Ci的第t个元素,且0≤t<n。
编码过程结束。
步骤三、针对二进制准循环码进行编码。
具体包括:
对于二进制准循环码,假设其生成矩阵G(ek×en)是由k×n组的二进制准循环矩阵组成,其中e=2r-1。其中,e为准循环码中循环块的大小,ek为编码前的消息长度, en为编码后的信息长度。
假设G=[Wi,j],其中Wi,j=circ(ai,j,0,ai,j,1,…,ai,j,e-1)是一个e×e的准循环矩阵,其中0≤i<k,0≤j<n。在GF(2r)中有λ组共轭组,分别表示为Ψ0,Ψ1,…,Ψλ-1。Ψi定义为其中ηi为第i个共轭组中元素的个数,ti表示第i个共轭组中的第一个元素。为GF(2r)的子域GF(2ηi)中的一组基。共轭组中的所有元素构成了GF(2r)中的所有元素。
转换方式与非二进制准循环码编码中G的转换方式相同,与步骤二1)相同。
需要编码的消息序列m=[m0,m1,…,me],其中mi是一个长为k的向量。
其中,编码后的码字为C=[C0,C1,…,Ci,…,Ce],其中Ci是一个长为n的子序列,即Ci=[ci,t],ci,t表示C中第i个子序列Ci的第t个元素,且0≤t<n。
至此,编码过程结束。
实施实例
对于在无线通信领域利用到的一组64进制的(4095,2160)准循环码,其中的元素为64进制,码长为4095位,信息位长为2160位,准循环码中循环块大小e为63。利用针对非二进制准循环码编码方法,可将一个长为2160位的64进制消息编码为长为4095位的64进制信息。相比于传统的编码方法,该编码方法的计算复杂度仅为1.59%。
对于在航天领域利用到的一组2进制的(4095,2160)准循环码,其中的元素为2进制,码长为4095位,信息位长为2160位,准循环码中循环块大小e为63。利用针对二进制准循环码编码方法,可将一个长为2160位的二进制消息编码为长为4095位的二进制信息。相比于传统的编码方法,该编码方法的计算复杂度仅为传统编码方法的9.52%。
对于在存储和航空领域利用到的一组2进制的(8176,7154)准循环码,其中的元素为2进制,码长为8176位,信息位长为7154位,准循环码中循环块大小e为511。利用针对二进制准循环码编码方法,可将一个长为7154位的二进制消息编码为长为8176位的二进制信息。相比于传统的编码方法,该编码方法的计算复杂度仅为传统编码方法的1.77%。
Claims (3)
1.一种准循环码的频域编码方法,其特征在于,包括以下几个步骤:
步骤一、判断准循环码的类别;
设消息序列为m,其中,m=[m0,m1,…,mi,…,me],其中mi是一个长为k的向量;
设n、e、k为正整数,ke为编码前消息序列m的长度,ne为编码后信息序列的长度,其中ke<ne;
(1)当(ne,ke)线性分组码Cqc中元素在GF(2)中,其中,GF(2)表示二元域,若符合以下条件:
1)Cqc中的码字由n组e个比特组成;
2)对Cqc中的码字的n组比特分别进行循环移位得到的码字仍然是Cqc中的码字;
则此线性分组码Cqc为二进制准循环码,转入步骤三;
(2)当(ne,ke)线性分组码Cqc中的元素在GF(2r)中,且满足上述1)与2)的条件,则此线性分组码Cqc为非二进制准循环码,转入步骤二,其中GF(2r)表示二元域GF(2)的r次方扩域;
步骤二、针对非二进制准循环码进行编码;
步骤三、针对二进制准循环码进行编码;
2.根据权利要求1所述的一种准循环码的频域编码方法,其特征在于,所述的步骤二具体包括以下几个步骤:
具体包括:
对非二进制准循环码,假设准循环码的生成矩阵G(ek×en)是由k×n组的非二进制准循环矩阵组成,其中e+1为2的r次幂,即e=2r-1,e表示准循环阵中每个小循环方阵的行数或者列数,每个准循环矩阵由k×n个小循环方阵构成,k表示准循环矩阵的行数,n表示准循环矩阵的列数;
假设G=[Wi,j],其中Wi,j=circ(ai,j,0,ai,j,1,…,ai,j,e-1)是一个e×e的准循环矩阵,且第一行为ai,j,0,ai,j,1,…,ai,j,e-1,其中0≤i<k,0≤j<n,Wi,j表示G矩阵中位于第i行j列的小矩阵,circ(ai,j,0,ai,j,1,…,ai,j,e-1)表示Wi,j中第一行的元素分别为ai,j,0,ai,j,1,…,ai,j,e-1;
将G转化为矩阵具体为:
其中:V=[α-ij],V-1=[αij],α为GF(2r)中的本原元,[αij]表示第i行j列的元素为αij的矩阵,diag表示对角矩阵;
设πrow和πcol分别为长为ek,en的序列,πrow中的第p个元素其中0≤p<ek,(p)e表示p/e的余数,表示p/e的商;0≤q<en,其中0≤q<en,(q)e表示q/e的余数,表示q/e的商;
将中第0行至第ek-1行中的元素分别置于第πrow(0)行至第πrow(ek-1)行中,然后再将第0列至第en-1列中的元素分别置于第πcol(0)列至第πcol(en-1)列中,得到矩阵其中Dt=[di,j,t],di,j,t为矩阵Dt中的第i行j列的元素,对于0≤i<k,0≤j<n,0≤t<e,有
需要编码的消息序列m=[m0,m1,…,mi,…,me],其中mi是一个长为k的子序列,即mi=[mi,t],mi,t表示m中第i个子序列mi的第t个元素,且0≤t<k;
其中,编码后的码字为C=[C0,C1,…,Ci,…,Ce],其中Ci是一个长为n的子序列,即Ci=[ci,t],ci,t表示C中第i个子序列Ci的第t个元素,且0≤t<n,α为GF(2r)中的本原元;
编码过程结束。
3.根据权利要求1所述的一种准循环码的频域编码方法,其特征在于,所述的步骤三具体包括以下几个步骤:
具体包括:
对于二进制准循环码,假设其生成矩阵G(ek×en)是由k×n组的二进制准循环矩阵组成,其中e=2r-1;其中,e为准循环码中循环块的大小,ek为编码前的消息长度,en为编码后的信息长度;
假设G=[Wi,j],其中Wi,j=circ(ai,j,0,ai,j,1,…,ai,j,e-1)是一个e×e的准循环矩阵,其中0≤i<k,0≤j<n;在GF(2r)中有λ组共轭组,分别表示为Ψ0,Ψ1,…,Ψλ-1;Ψi定义为其中ηi为第i个共轭组中元素的个数,ti表示第i个共轭组中的第一个元素;为GF(2r)的子域GF(2ηi)中的一组基;共轭组中的所有元素构成了GF(2r)中的所有元素;
其中:V=[α-ij],V-1=[αij],α为GF(2r)中的本原元,[αij]表示第i行j列的元素为αij的矩阵,diag表示对角矩阵;
设πrow和πcol分别为长为ek,en的序列,πrow中的第p个元素其中0≤p<ek,(p)e表示p/e的余数,表示p/e的商;0≤q<en,其中0≤q<en,(q)e表示q/e的余数,表示q/e的商;
将中第0行至第ek-1行中的元素分别置于第πrow(0)行至第πrow(ek-1)行中,然后再将第0列至第en-1列中的元素分别置于第πcol(0)列至第πcol(en-1)列中,得到矩阵其中Dt=[di,j,t],di,j,t为矩阵Dt中的第i行j列的元素,对于0≤i<k,0≤j<n,0≤t<e,有
需要编码的消息序列m=[m0,m1,…,me],其中mi是一个长为k的向量;
此时中的元素并不满足共轭关系,共轭关系指的是对于序列中的第t个子序列和第(2t)e个子序列,其中0≤t<e,都存在即子序列中每个元素进行平方后都等于中的对应元素,其中(2t)e表示(2t)/e的余数;
其中,编码后的码字为C=[C0,C1,…,Ci,…,Ce],其中Ci是一个长为n的子序列,即Ci=[ci,t],ci,t表示C中第i个子序列Ci的第t个元素,且0≤t<n;
编码过程结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310134009.7A CN103208996B (zh) | 2013-04-17 | 2013-04-17 | 准循环码的频域编码方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310134009.7A CN103208996B (zh) | 2013-04-17 | 2013-04-17 | 准循环码的频域编码方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103208996A true CN103208996A (zh) | 2013-07-17 |
CN103208996B CN103208996B (zh) | 2016-06-22 |
Family
ID=48756099
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310134009.7A Expired - Fee Related CN103208996B (zh) | 2013-04-17 | 2013-04-17 | 准循环码的频域编码方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103208996B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017041232A1 (zh) * | 2015-09-08 | 2017-03-16 | 广东超算数据安全技术有限公司 | 一种二进制循环码的编解码框架 |
CN107026654A (zh) * | 2016-02-02 | 2017-08-08 | 中国科学院声学研究所 | 一种准循环多进制低密度奇偶校验码的快速频域编码方法 |
CN111447158A (zh) * | 2020-03-30 | 2020-07-24 | 电子科技大学 | 一种基于有限域傅里叶变换的信道编码特征提取方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1964200A (zh) * | 2006-11-21 | 2007-05-16 | 华为技术有限公司 | 一种线性码扩展编码及解码方法以及数据传输方法及装置 |
US20100058150A1 (en) * | 2008-08-26 | 2010-03-04 | Hiroyuki Yamagishi | Coding Apparatus and Method |
CN101699808A (zh) * | 2009-11-12 | 2010-04-28 | 上海交通大学 | 差分编码空时频调制方法 |
CN102457286A (zh) * | 2010-10-21 | 2012-05-16 | 航天信息股份有限公司 | 准循环ldpc码编码方法、装置及校验矩阵生成方法 |
-
2013
- 2013-04-17 CN CN201310134009.7A patent/CN103208996B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1964200A (zh) * | 2006-11-21 | 2007-05-16 | 华为技术有限公司 | 一种线性码扩展编码及解码方法以及数据传输方法及装置 |
US20100058150A1 (en) * | 2008-08-26 | 2010-03-04 | Hiroyuki Yamagishi | Coding Apparatus and Method |
CN101699808A (zh) * | 2009-11-12 | 2010-04-28 | 上海交通大学 | 差分编码空时频调制方法 |
CN102457286A (zh) * | 2010-10-21 | 2012-05-16 | 航天信息股份有限公司 | 准循环ldpc码编码方法、装置及校验矩阵生成方法 |
Non-Patent Citations (1)
Title |
---|
王祖林: "循环冗余校验码的查表生成算法及其实现", 《北京航空航天大学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017041232A1 (zh) * | 2015-09-08 | 2017-03-16 | 广东超算数据安全技术有限公司 | 一种二进制循环码的编解码框架 |
CN107026654A (zh) * | 2016-02-02 | 2017-08-08 | 中国科学院声学研究所 | 一种准循环多进制低密度奇偶校验码的快速频域编码方法 |
CN107026654B (zh) * | 2016-02-02 | 2019-06-18 | 中国科学院声学研究所 | 一种准循环多进制低密度奇偶校验码的快速频域编码方法 |
CN111447158A (zh) * | 2020-03-30 | 2020-07-24 | 电子科技大学 | 一种基于有限域傅里叶变换的信道编码特征提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103208996B (zh) | 2016-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107832837B (zh) | 一种基于压缩感知原理的卷积神经网络压缩方法及解压缩方法 | |
CN103957015B (zh) | 用于ldpc码解码的非均匀量化编码方法及其在解码器的应用 | |
CN100414841C (zh) | 一种ldpc编码的实现方法 | |
CN101902228B (zh) | 快速循环冗余校验编码方法及装置 | |
DE602008003456D1 (de) | Erzeugung von paritätsprüfmatrizen | |
CN101834613A (zh) | 一种ldpc码的编码方法及编码器 | |
CN103746774A (zh) | 一种高效数据读取的容错编码方法 | |
CN101159435B (zh) | 基于移位矩阵分级扩展的低密度校验码校验矩阵构造方法 | |
PH12019502212A1 (en) | Method for performing encoding on basis of parity check matrix of low density parity check (ldpc) code in wireless communication system and terminal using same | |
CN101699770A (zh) | 一种高速ldpc码编码器及其编码方法 | |
CN103208996A (zh) | 准循环码的频域编码方法 | |
CN103220083A (zh) | 一种对上行控制信息进行编码的方法和装置 | |
CN109644006A (zh) | 编码数据和解码数据的装置及方法 | |
CN107026654B (zh) | 一种准循环多进制低密度奇偶校验码的快速频域编码方法 | |
CN103236860A (zh) | 用于生成ldpc码校验矩阵的方法、及该ldpc码编码方法 | |
CN101465655B (zh) | 极短码长低密度奇偶校验码的编码方法 | |
CN114285417A (zh) | 一种全局耦合低密度奇偶校验码的编码方法 | |
CN101577554B (zh) | 多码长多码率的低密度奇偶校验码的编码方法 | |
CN101882972B (zh) | 一种Raptor码的解码方法 | |
CN103368585B (zh) | 一种ldpc码校验矩阵的构造方法 | |
CN101789795B (zh) | 基于多码率原模图ldpc码的编码方法及编码器 | |
CN101471671A (zh) | Ldpc编码中校验比特的计算方法和装置 | |
CN103401566A (zh) | 参数化的bch纠错码的并行编码方法及装置 | |
CN109639290A (zh) | 一种半随机分组叠加编码及译码方法 | |
CN101834612A (zh) | 一种ldpc码的编码方法及编码器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160622 Termination date: 20170417 |