首页> 中国专利> CMMB中LDPC码的译码算法及部分并行译码器

CMMB中LDPC码的译码算法及部分并行译码器

摘要

本发明公开了一种CMMB中LDPC码的译码算法,包括以下步骤:(a)对LDPC码的现有校验矩阵进行置换,得到由18个256×9216的子矩阵组成的新校验矩阵,且每个子矩阵的每一列仅包含一个“1”;(b)将新校验矩阵应用于部分并行译码器中,控制单元控制整个部分并行译码器进行系统初始化;(c)信息更新单元接收变量节点存储单元的信息,与校验节点存储单元的信息共同作用,完成所有变量节点和校验节点的信息更新,直至CMMB中LDPC码的译码完成。本发明减少了译码过程中的信息存储量,有效地提高了信息收敛速度,降低了迭代次数,提高了译码效率;同时,还降低了硬件的实现难度。

著录项

  • 公开/公告号CN102291153A

    专利类型发明专利

  • 公开/公告日2011-12-21

    原文格式PDF

  • 申请/专利权人 电子科技大学;

    申请/专利号CN201110157038.6

  • 发明设计人 武畅;李玉柏;谭太秋;

    申请日2011-06-13

  • 分类号H03M13/11(20060101);

  • 代理机构

  • 代理人

  • 地址 610000 四川省成都市高新区(西区)西源大道2006号

  • 入库时间 2023-12-18 04:08:41

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2015-08-05

    未缴年费专利权终止 IPC(主分类):H03M13/11 授权公告日:20140402 终止日期:20140613 申请日:20110613

    专利权的终止

  • 2014-04-02

    授权

    授权

  • 2012-02-08

    实质审查的生效 IPC(主分类):H03M13/11 申请日:20110613

    实质审查的生效

  • 2011-12-21

    公开

    公开

说明书

技术领域

本发明涉及一种译码器,具体地说,是涉及一种CMMB中LDPC码的部分并 行译码算法及译码器。

背景技术

低密度奇偶校验(LDPC)码是一种具有逼近shannon限的优秀纠错码,具 有极强的纠错和检错能力。近年来,低密度奇偶校验码的优异性能及其良好的 应用前景已引起了研究人员的高度重视,现已成功应用于多个行业标准,如第 二代欧洲数字电视广播标准(DVB-S2)、地面数字电视广播标准(CDTTB)、我国 2006年提出的中国移动多媒体广播标准(CMMB)。

CMMB系统的信道编码采用LDPC码作为内码,此LDPC码的现有校验矩阵是 经过精心设计构造出来的一种高度结构化的稀疏矩阵,具有巧妙的准循环特点:

(1)1/2码率的校验矩阵H为4608×9216,可以划分为256个18×9216 的行子矩阵,其中下一个子矩阵可通过上一个子矩阵向右循环移36位而得到; 也可以划分为256个4608×36的列子矩阵,其中后一个子矩阵可通过前一个子 矩阵向下移18位而得到;这种校验矩阵的行重为6,列重为3。

(2)3/4码率的校验矩阵H为2304×9216,可以划分为256个9×9216的 行子矩阵,其中下一个子矩阵可通过上一个子矩阵向右循环移36位而得到;也 可以划分为256个2304×36的列子矩阵,其中后一个子矩阵可通过前一个子矩 阵向下移9位而得到;这种校验矩阵的行重为12,列重为3。

在校验矩阵中只有极少数矩阵元素为“1”,绝大多数矩阵元素为零,上述 所谓校验矩阵的行重即是指矩阵的每行中“1”的个数,而所谓的列重即是指矩 阵的每列中“1”的个数。

CMMB中LDPC码的标准校验矩阵很大,为了便于分析和阐述原理,此处构造 一个和标准校验矩阵具有类似结构的较小的矩阵Ha以便分析。

上述矩阵Ha是用与标准校验矩阵相似的构造方法构造的规则为9×18的矩 阵,其行重为6,列重为3,准循环特性为:每隔3行,“1”的位置向右循环移 6列,每隔6列,“1”的位置向下循环移3行。因此,在上述矩阵Ha中,可以 将矩阵的第4行和第7行看成为第1行循环移位生成的两行,同理,第5行和 第8行看成为第2行循环移位生成的两行,第6行和第9行看成为第3行循环 移位生成的两行。

每一个二进制LDPC码均可以用上述类似的一个M×N稀疏校验矩阵进行描 述,其中行数M表示校验节点数目,列数N表示变量节点数目(也可称为信息 节点)。LDPC(Low Density Parity Check Code)译码主要是通过对数似然信 息在校验节点和变量节点之间传递更新而迭代译码。在目前的研究中,LDPC译 码器的译码结构可以分为全并行结构、全串行结构和部分并行结构。其中,全 并行结构中所有的变量节点和校验节点同时更新,译码速度快,但是并行处理 单元太多,各处理单元与存储单元之间的连线复杂度随着码长增加而急剧增加, 给硬件设计造成了极大的困难;全串行结构每次只更新一行中的校验节点和变 量节点,一次迭代需要m次校验节点更新和变量节点更新,译码延迟大,不利 于在实时通信系统中应用;而部分并行结构是复杂度和译码速度的折中,可根 据需要调整并行量。

按照消息传递方式的不同,目前LDPC码的译码算法可分为两种,一是TPMP (two-phase message passing)算法;另一种是将校验节点更新和变量节点更 新融合在一起的TDMP(turbo decoding message passing)算法。下面分别对 这两种现有技术做简要介绍。

现有技术一

该方案采用校验节点更新和变量节点更新分开处理的TPMP算法,其大致过 程为:

步骤1.初始化,θn→m=L(cn),Λm→n=0。

步骤2.校验节点更新,也叫行更新。对每个m,n∈N(m),计算

Λmn=2tanh-1{ΠnN(m)\ntanh[θnm/2]}---(1)

步骤3.变量节点更新,也叫列更新。对每个n,m∈M(n),计算

θnm=L(cn)+ΣmM(n)\mΛmn---(2)

λn=L(cn)+ΣmM(n)Λmn---(3)

步骤4.itea=itea+1,如果itea<Max_it,转到步骤2,否则转到步骤5。

步骤5.根据λn的值判定第n比特的接收值,n=1,...,N。

现有技术一由于将变量节点更新和校验节点更新独立开来,因此,在处理 校验节点更新的时候,必须把所有的行都处理完,而在处理变量节点更新的时 候,必须把所有的列都处理完,不能在处理一种节点更新的过程中对另一种节 点进行更新。这种处理方法的缺点在于消息更新速度太慢,译码速度得不到提 高。而且,在处理消息更新的同时,还需要存储大量的中间变量。以CMMB标准 为例,在处理校验节点更新的时候,需要存储4608×6个更新变量;在处理变 量节点更新的时候,需要存储9216×3个更新变量,这将需要占用大量的硬件 资源,因此,对硬件设备的要求极高。

现有技术二

针对现有技术一的不足,近年来提出了一种将校验节点更新和变量节点更 新融合在一起的消息更新机制,即TDMP算法。该算法的实现过程将现有技术一 的步骤2和步骤3融合在了一起,该步骤具体如下:

ρn=λnm→n    (4)

Λmn=ΠnM(n)\nsign(ρn)max{minnM(n)\n(|ρn|-β),0}---(5)

λn=ρnm→n    (6)

其中,sign()是符号函数,min()和max()分别是求最小值和最大值的函数。 β是补偿最小和算法中的补偿因子,取值在0到1之间。

现有技术二在计算完一行的信息更新后,立即更新对应列的后验信息,更 新后的列又用于下一行的计算。这样在一次迭代中,消息不断的在行和列之间 传递,大大加快了迭代收敛的速度。和现有技术一相比,它能降低大约50%的迭 代次数。同时在中间计算的时候,不需要存储对应节点的更新信息,可减少9216 ×3个存储单元。校验节点更新信息的计算也可通过采用补偿最小和算法进行简 化,存储最小值、次小值以及每个变量更新后的符号。

该方案虽然有许多优点,但它的并行实现目前只能适用于特定的校验矩阵。 这样的校验矩阵有一个特点,即每次并行处理的子矩阵中每一列的列重必须小 于等于1,否则将会导致变量节点更新信息还未来得及使用便被覆盖,破坏消息 传递结构而导致译码错误。很显然,从前述构造的矩阵Ha可以看出,CMMB中LDPC 码的标准校验矩阵中列重并非小于等于1,因此,使用现有技术二来进行CMMB 中LDPC码的译码并不适合。

发明内容

本发明的目的在于提供一种CMMB中LDPC码的译码算法,解决现有CMMB中 LDPC码的译码过程中更新信息的存储量大、译码速度慢的问题,在实现对CMMB 标准中LDPC码的正常译码的同时,减少更新信息的存储量和迭代次数,提高译 码速度。

为了实现上述目的,本发明采用的技术方案如下:

CMMB中LDPC码的译码算法,包括以下步骤:

(a)对LDPC码的现有校验矩阵进行置换,得到由18个256×9216的子矩 阵组成的新校验矩阵,且每个子矩阵的每一列仅包含一个“1”;

(b)将新校验矩阵应用于部分并行译码器中,控制单元控制整个部分并行 译码器进行系统初始化;

(c)信息更新单元接收变量节点存储单元的信息,与校验节点存储单元的 信息共同作用,完成所有变量节点和校验节点的信息更新,直至LDPC码的译码 完成。

具体地说,所述步骤(a)中对LDPC码的现有校验矩阵进行置换的置换方 法为:(a1)将LDPC码的现有校验矩阵中第1行和由第1行循环移位生成的所 有行构成第1个子矩阵,该子矩阵占据新校验矩阵的第1行至第256行;(a2) 将LDPC码的现有校验矩阵中第2行和由第2行循环移位生成的所有行构成第2 个子矩阵,该子矩阵占据新校验矩阵的第257行至第512行;(a3)按照上述方 法,类推至LDPC码的现有校验矩阵的第18行,得到18个256×9216的子矩阵, 从而构成LDPC码的新校验矩阵。

所述步骤(c)的具体步骤包括:(c1)预处理单元对输入的校验节点信息 和变量节点信息进行处理;(c2)最小值次小值计算单元按照补偿最小和(OMS) 算法对预处理后的信息进行处理,并更新当前行的校验节点信息,然后由变量 节点更新单元完成信息更新,即更新与当前处理行的校验节点相连的变量节点 信息;(c3)判断是否完成所有变量节点的信息更新,如果是,则执行(c5), 反之,则执行(c4);(c4)变量节点存储单元将变量节点的更新信息传输至与 之相连的下一行校验节点,并执行(c1)~(c3);(c5)根据λn的值判定第n比 特的接收值,n=1,...,N。

更具体地说,所述步骤(c1)的具体方法为:首先从变量节点存储单元中 读取和当前处理行校验节点相连的变量节点信息,然后再从当前处理行的校验 节点存储单元中读取当前行的当前校验节点信息,两者相减,得到预处理单元 输出信息,将该信息输入最小值次小值计算单元,并存入对应的移位存储器中。

所述步骤(c2)中补偿最小和算法的具体方法为:首先对预处理单元的处 理结果求取最小值min1、次小值min2,以及求取处理结果的所有值的符号的乘 积signs;然后针对每一个校验节点的更新信息,取出它的符号signx与signs做异 或运算得到更新后的符号signf;最后将预处理单元的处理结果的绝对值与min1 进行比较,若相等,则它更新后的值为(signf)min2,反之,则它更新后的值为 (signf)min1。

所述步骤(c2)中变量节点更新单元完成变量节点的信息更新的具体方法 为:将经过最小值次小值计算单元处理后的校验节点的更新信息与移位寄存器 中存储的值相加得到与校验节点相连的变量节点的更新信息。

以上述译码算法为基础,本发明还提供了一种CMMB中LDPC码的部分并行 译码器,包括变量节点存储单元,与该变量节点存储单元进行信息交换的至少 一个信息更新单元,以及对变量节点存储单元和所有信息更新单元进行控制的 控制单元。

进一步地说,所述信息更新单元包括用于存储校验节点信息的校验节点存 储单元、用于对待更新信息进行预处理的预处理单元、用于对预处理单元的处 理结果进行补偿最小和算法的最小值次小值计算单元,以及用于完成变量节点 信息更新的变量节点更新单元和移位寄存器。

再进一步地说,所述校验节点存储单元和变量节点存储单元均由存储器构 成;控制单元由状态器构成,预处理单元由处理器构成,而最小值次小值计算 单元则由计算器构成。

本发明的设计原理在于:通过改变LDPC码的校验矩阵,来提高校验节点与 变量节点之间的信息收敛速度,减少校验节点与变量节点之间的迭代次数,进 而提高校验节点和变量节点之间的信息传递效率,最终达到提高译码速度的目 的。LDPC码的校验矩阵的改变,在于利用LDPC码的现有校验矩阵的准循环特性 进行巧妙置换,使构成的新校验矩阵中,每个子矩阵的列重为1,从而解决CMMB 中LDPC码的现有校验矩阵中每个子矩阵的列重超过1导致在译码过程中可能出 现更新信息尚未使用便被覆盖的问题。

与现有技术相比,本发明具有以下有益效果:

(1)本发明设计巧妙,原理简单,实现方便,有效地解决了现有CMMB中 LDPC码的译码算法信息存储量大、速度慢等问题;

(2)本发明中部分并行译码器的硬件结构简单,变量节点存储单元只需要 存储校验矩阵的列数个变量节点即可,因此,存储单元大大减少,节省了存储 空间,减小了硬件资源消耗;

(3)本发明中将校验节点存储单元集成于信息更新单元之中,大大减小了 硬件的布线难度,使校验节点信息的读写更加方便,硬件更加容易实现;

(4)本发明采用巧妙的置换方法,实现了新校验矩阵中每个子矩阵的列重 为1,从而有效地避免了译码过程中同一节点在先更新的信息在尚未使用之前便 被在后更新的信息覆盖的问题,保证了译码过程的顺利进行,提高了译码的准 确率;

(5)本发明中采用校验节点、变量节点、校验节点、变量节点交叉的顺序 进行信息更新,且在同一次迭代中,变量节点还会用更新后的信息参与到其他 校验节点更新中,这种信息传递机制更充分地利用了已知信息,大大提高了信 息的收敛速度,降低了迭代次数,提高译码吞吐率,最终实现了提高译码速度 的效果。

附图说明

图1为本发明中部分并行译码器的系统框图。

图2为图1中信息更新单元的系统框图。

图3为本发明中部分并行译码器的译码算法的流程示意图。

图4为本发明中部分并行译码器的译码算法中步骤(c)的流程示意图。

图5为现有技术中矩阵Ha的tanner图。

图6为本发明-实施例中矩阵Ha′的tanner图。

具体实施方式

下面结合附图和实施例对本发明作进一步说明,本发明的实施方式包括但 不限于下列实施例。

首先对本发明中各个字母及函数代表的含义进行解释,具体如下:

L(cn)表示随机变量cn的对数似然比,L(cn)=log{P(cn=0|rn)/P(cn=1|rn)};

集合N(m)={n;Hmn≠0}表示与校验节点m相连的所有变量节点;

集合M(n)={m;Hmn≠0}表示与变量节点n相连的所有校验节点;

N(m)\n表示N(m)中除了n以外的其他点的集合;

M(n)\m表示M(n)中除了m以外的其他点的集合;

Λm→n表示校验节点m传递给变量节点n的信息;

θn→m表示变量节点n传递给校验节点m的信息;

λn表示第n个变量节点的信息;

Max_it是设定的最大迭代次数;

itea表示迭代次数;

Ln表示接收到的第n个比特信道信息的对数似然比;

Λij表示第i行校验节点传递给第j个变量节点信息的对数似然比,i=1,...,M, M表示行数,j=1,...,ci,ci表示第i行的行重;

Ik表示第k行中所有1所在的位置;

集合N(m)={n;Hmn≠0}表示与校验节点m相连的所有变量节点;

N(m)\j表示N(m)中除了j以外的所有点的集合。

实施例

如图1和图2所示,采用32路并行译码,每个信息更新单元处理2行,因 此CMMB中LDPC码的部分并行译码器,由控制单元、变量节点存储单元和16个 信息更新单元。其中,信息更新单元是本发明中的核心器件,主要用于实现校 验节点和变量节点之间的信息更新,其信息更新的速度和准确率直接关系到整 个部分并行译码器的译码速度和准确率;变量节点存储单元主要用于存储变量 节点的更新信息;而控制单元则是整个部分并行译码器的控制中心,它向部分 并行译码器内部的各个功能器件发出控制信号,控制各功能器件执行各自的工 作。所有的信息更新单元分别采用双向总线与变量节点存储单元连接,以实现 两者之间的信息交换;而控制单元分别通过单向总线与变量节点存储单元和所 有的信息更新单元连接,以便向变量节点存储单元和所有的信息更新单元发出 控制信号,实现对各个功能器件的控制。

上述16个信息更新单元结构相同,均由校验节点存储单元、预处理单元、 移位存储器、最小值次小值计算单元和变量节点更新单元组成。其中,校验节 点存储单元用于存储校验节点信息;预处理单元用于对待更新的节点信息进行 预处理,并将处理结果输入最小值次小值计算单元和移位存储器中,以便进一 步处理;最小值次小值计算单元完成校验节点的信息更新;而变量节点更新单 元则是用于实现变量节点的信息更新。在上述结构中,预处理单元分别通过移 位存储器和最小值次小值计算单元与变量节点更新单元连接,而校验节点存储 单元则分别与预处理单元和最小值次小值计算单元连接。

下面以1/2码率的校验矩阵为例,对本发明中部分并行译码器的译码算法 进行详细说明。

假设S为并行处理的行数,D为子矩阵的个数,满足校验节点数目M=D×S, 则本实施例中部分并行译码器的总体译码过程可以按照以下步骤进行描述:

步骤1.对现有校验矩阵进行置换,得到本发明中所需的新校验矩阵;

步骤2.部分并行译码器初始化,λn=Ln,n=1,...,N;Λij=0,i=1,...,M;j=1,...,ci

步骤3.对所有的d,d从1到D,计算当前处理行的校验节点的更新信息, 每次并行计算S行,然后完成与当前处理行的校验节点相连的变量节点的信息 更新;

ρjs=λ(I(d-1)S+s,j)-Λ(d-1)S+s,j---(7)

Λ(d-1)S+s,j=ΠjN(m)\jsign(ρjs)max{minjN(m)\j(|ρjs|-β),0}---(8)

λ(I(d-1)S+s,j)=ρjs+Λ(d-1)S+s,j---(9)

步骤4.itea=itea+1,如果itea<Max_it,转到步骤2,否则转到步骤4。

步骤5.根据λn的值判定第n比特的接收值,n=1,...,N。

为了更加容易理解上述步骤1~步骤5,使本发明的实现更加容易,下面分 别对上述步骤1~步骤5进行详细说明。本发明进行CMMB中LDPC码的译码流程 示意图如图3和图4所示。

步骤1.对LDPC码的现有校验矩阵进行置换,得到由18个256×9216的子 矩阵组成的新校验矩阵,且每个子矩阵的列重均为1。

根据现有技术可知,LDPC码的现有校验矩阵具有巧妙的准循环特点,利用 该准循环特点对现有校验矩阵进行置换便可以得到本发明中所需的新校验矩 阵,具体的置换方法为:首先,将LDPC码的现有校验矩阵中第1行和由第1行 循环移位生成的所有行构成第1个子矩阵,该子矩阵占据新校验矩阵的第1行 至第256行;然后再将LDPC码的现有校验矩阵中第2行和由第2行循环移位生 成的所有行构成第2个子矩阵,该子矩阵占据新校验矩阵的第257行至第512 行;按照上述方法,依次类推至LDPC码的现有校验矩阵的第18行,即最后一 行,得到18个256×9216的子矩阵,这18个子矩阵按照顺序排列便形成了本 发明中LDPC码的新校验矩阵。

在形成LDPC码的新校验矩阵之后,必须验证其对译码结果的影响。由于LDPC 码的译码主要是通过信息在校验节点m和变量节点n之间传递而迭代译码,而 信息的传递则按照校验矩阵所描述的tanner图进行,因此,验证新校验矩阵对 译码结果是否存在影响的方法便是将LDPC码的新校验矩阵的tanner图和现有 校验矩阵的tanner图进行比较,如果tanner图中信息的传递结构相同则表明 新校验矩阵对译码结果无影响,反之,则说明新校验矩阵会对译码结果产生影 响,不能使用。为了方便对现有校验矩阵和新校验矩阵的tanner图进行比较, 本实施例中以现有技术中构造的矩阵Ha为原形,采用上述置换方法得到一个新 矩阵Ha′。

明显可以看出,矩阵Ha′由三个子矩阵H1、H2、H3构成,且每个子矩阵的 列重均为1。矩阵Ha和矩阵Ha′的tanner图分别为图5和图6。从图5可以看出, 在校验节点更新时,变量节点的信息通过与校验节点的连线传递到相应的校验 节点;变量节点更新时,校验节点的信息通过与变量节点的连线传递到相应的 变量节点。而从图5和图6的对比来看,两者只是校验节点的位置进行了交换, 变量节点和校验节点之间的连接关系并没有变化,由此便说明信息的传递结构 没有变化,所以从Ha到Ha′的变换对译码结果没有影响,进而推断采用本发明所 述的置换方法对LDPC码的现有校验矩阵进行置换,所得到的新校验矩阵对译码 结果也没有影响。

步骤2.将置换所得的新校验矩阵应用于部分并行译码器中,由控制单元发 出控制信号,控制整个部分并行译码器进行系统初始化。

控制单元由控制器构成。在系统初始化过程中,控制单元向系统输入先验 信息,并生成各存储单元的地址信息,此处的存储单元是指变量节点存储单元 和所有的校验节点存储单元。另外,在译码过程中,控制单元还进行迭代是否 终止的判断,即进行图4所示所有的变量节点的信息是否均已经更新的判断: 所有的变量节点的信息均已更新表明迭代终止,反之,则迭代继续。

步骤3.本步骤实质上可以细分为以下三个小步骤,分别为:

(I)预处理单元对输入的变量节点信息和校验节点信息处理;

预处理单元的具体工作过程为:预处理单元首先从变量节点存储单元中读 取和当前处理行校验节点相连的变量节点信息,然后再从校验节点存储单元中 读取当前行的校验节点信息,两者相减,即完成(7)式所描述的运算,所得结 果同时传输至最小值次小值计算单元和移位存储器中。

在本实施例中,LDPC码的新校验矩阵由18个256×9216的子矩阵排列构成, 因此,必然设置有18个校验节点存储单元,用于分别存储每个子矩阵的校验节 点信息。而每个子矩阵的行数为256行,因此在计算子矩阵的校验节点信息时, 每次并行计算的行数也为256行。

(II)最小值次小值计算单元按照补偿最小和算法对预处理单元的输出信 息进行处理;

最小值次小值计算单元完成(8)式所描述的补偿最小和算法,具体方法为: 由于预处理单元每次处理的结果均包含多个带符号的数值,因此,首先对预处 理单元的处理结果求取最小值min1、次小值min2,以及求取处理结果的所有值 的符号的乘积signs;然后针对每一个校验节点的更新信息,取出它的符号signx与 signs做异或运算得到更新后的符号signf;最后将预处理单元的处理结果的绝对 值与min1进行比较,若相等,则它更新后的值为(signf)min2,反之,则它更新后 的值为(signf)min1。本部分所述的符号是指数值的正号与负号。

(III)由变量节点更新单元对与当前处理行的校验节点相连的变量节点进 行信息更新;

具体地说,变量节点更新单元的功能在于完成(9)式的运算,其方法为: 由经过最小值次小值计算单元处理后的校验节点的更新信息与移位寄存器中存 储的值相加得到与该校验节点相连的变量节点的更新信息。

变量节点的更新信息存储于变量节点存储单元之中。

步骤4.判断是否完成所有的变量节点的信息更新;

在本步骤中,采用的具体判断方法是设定最大迭代次数Max_it,然后在每 完成一个变量节点的信息更新之后,相应的迭代次数加1,再判断该迭代次数是 否大于设定的最大迭代次数Max_it。如果大于设定的最大迭代次数Max_it,则 证明完成了所有的变量节点的信息更新,反之,则证明尚未完成所有的变量节 点的信息更新。

如果尚未完成所有的变量节点的信息更新,则由变量节点存储单元将当前 处理行的变量节点的更新信息传输至与之相连的下一行的校验节点,然后再重 复步骤3,直到完成所有的变量节点的信息更新。

步骤5.根据λn的值判定第n比特的接收值,n=1,...,N;

在所有变量节点的信息更新完成之后,只需要得出每个变量节点的接收值, 即可完成CMMB中LDPC码的译码,此步骤与现有技术一中的步骤5完全相同, 在此不再赘述。

按照上述实施例,便可很好地实现本发明。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号