CN102270992A - Turbo译码器中前向递推概率获取方法 - Google Patents
Turbo译码器中前向递推概率获取方法 Download PDFInfo
- Publication number
- CN102270992A CN102270992A CN2011100721694A CN201110072169A CN102270992A CN 102270992 A CN102270992 A CN 102270992A CN 2011100721694 A CN2011100721694 A CN 2011100721694A CN 201110072169 A CN201110072169 A CN 201110072169A CN 102270992 A CN102270992 A CN 102270992A
- Authority
- CN
- China
- Prior art keywords
- forward recursion
- recursion probability
- probability
- state
- prime
- 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
Landscapes
- Error Detection And Correction (AREA)
Abstract
本发明公开了一种Turbo译码器中获取前向递推概率的方法,目的是减少前向递推概率存储容量,并减少时间开销。技术方案是在获取每个时刻的前向递推概率时,将此时刻所有状态的前向递推概率都减去0状态的前向递推概率,保证所有时刻非0状态的前向递推概率都是偶数,同时也保证所有时刻所有状态的前向递推概率不会溢出,避免了传统方法中通过减去前向递推概率中的最大值来防溢出时所带来的存储面积开销和逻辑延迟。采用本发明可使得前向递推概率存储器容量减小为(G-1)×(B-1)×W,且使得获取一个时刻的前向递推概率的时间开销只有τsub,极大地提高了译码速率,同时,通过实验模拟仿真,本发明也可保证与传统方法一样的误码率性能。
Description
技术领域:本发明涉及一种应用于turbo译码器中的前向递推概率获取方法,属于纠错码领域。
背景技术:turbo码是一种译码性能优越的前向纠错码,自1993年由C.Berrou,A.Glavieux和P.Thitimajshiwa三人在瑞士日内瓦举办的通信国际年会(ICC)上提出以来,立刻备受关注。目前,Turbo码已经成功地应用到磁介质光介质数据存储、多媒体和有线、无线、光纤、星载通信等多个领域。
Turbo译码过程如图1所示,首先对一帧长度为N的数据正向地计算前向递推概率,然后反向地计算后向递推概率,最后再反向地计算后验概率。后验概率根据前向、后向递推概率来计算,且后向递推概率和后验概率的计算顺序都是反向的,它们可以同时计算,因此Turbo译码过程必须先正向获取前向递推概率并保存,然后供反向计算后验概率时使用。
图2是3GPP LTE通信协议采用的Turbo译码网格图和编码网格图,译码是编码的逆过程,一个Turbo译码器必须有与之相应的Turbo编码器,译码网格图和编码网格图必须相同。图2中箭头方向是网格图的正向,与箭头相反的方向是网格图的反向,前向递推概率按照正向获取,后向递推概率按照反向获取。网格图中,圆圈表示每个时刻的状态,状态的个数由通信协议采用的Turbo编码器决定,状态的个数G等于2D,D是与Turbo译码器相应的编码器中寄存器单元的个数,那么每个状态可以用D位的二进制整数表示,如3GPP LTE协议编码器中延迟单元个数是3,那么状态的个数就是8个,即000、001、……、111,也可以表示为s0、s1、......、s7,箭头表示从k时刻到k+1时刻的状态跳转,称之为路径,每条跳转路径上的Xk/P1kP2k表示的是k时刻输入为Xk,输出为P1kP2k,Xk,P1k,P2k∈{0,1}。例如,网格图中k时刻的状态是s2,假设输入Xk为0,那么k时刻的状态就按照Xk/P1kP2k等于0/10路径跳转到k+1时刻的状态s5。图2中的虚线表示输入Xk为0,实线表示输入Xk为1。
k+1时刻状态s的前向递推概率与k时刻规约到状态s的路径有关,如图2中k+1时刻s0状态的前向递推概率就与k时刻的s0和s4状态有关,因为k时刻的s0和s4状态分别有一条路径规约到k+1时刻的s0状态。设规约到k+1时刻状态s在k时刻的状态分别为s1’和s2’,获取k+1时刻状态s的前向递推概率As,k+1的简化公式为:
公式一是对Turbo译码算法简化后的公式,称之为Max*-Log-MAP算法。其中的As,k+1是k+1时刻状态s的前向递推概率,s∈{s0、s1、......、sG-1},和分别是k时刻状态s1’和s2’的前向递推概率,fLUT是为了补偿公式简化时丢失的精度,称之为查找表LUT(Look-Up Table)函数,如图7中的曲线示意的就是fLUT函数。γs′→s,k是图2中每条路径的分支量度,它的计算公式为:
公式二中s′∈{s′1、s′1},系统信息校验信息和是从无线信道接收到的信息,m(Xk)是图2输入比特Xk经过反相调制(0调制到1;1调制到-1)后的信息,m(P1k)和m(P2k)是输出比特P1k和P2k经过调制后的信息,即m(Xk),m(P1k),m(P2k)∈{1,-1},Lc是信道置信度,先验信息L(Xk)是Turbo译码时的反馈信息。
对一帧长度为N的数据,根据图2网格图的正向重复计算N-1次,即完成前向递推概率的获取过程。前向递推概率存储器的宽度B有限,采用二进制补码时,前向递推概率存储器保存的前向递推概率值的范围是:-2B-1~2B-1-1。从公式一看出获取前向递推概率是不断叠加的过程,不可避免地会导致计算的前向递推概率值超出-2B-1~2B-1-1的范围,使前向递推概率存储器保存的前向递推概率值发生错误,降低正确性。
图3是传统获取前向递推概率的流程图,分为下列几个步骤。
第一步:初始化时刻计数值k=0。将k=0时刻s0状态的前向递推概率设置为0,将k=0时刻s1、s2、......、sG-1状态的前向递推概率 都设置成-2B-1,B是前向递推概率值的宽度,G等于2D,D为Turbo编码器中寄存器单元的个数。
查找逻辑中2选1查找逻辑的级数为假设一级2选1查找逻辑的延迟是τmax,完成第二步总的延迟是如图4所示,假设G等于8,那么查找8个前向递推概率中的最大值,需要3级2选1的查找逻辑,完成第二步总的延迟是3×τmax。
第四步:前向递推概率防溢出处理。将第二步得到的G个前向递推概率都减去第三步找出的最大前向递推概率Max Ak+1,得到k+1时刻G个新的前向递推概率设减法操作的延迟是τsub,传统方法获取一个时刻的前向递推概率的时间开销是3×τmax+τsub。
第六步:k=k+1。
第七步:判断k是否等于N-1,N为一帧数据长度如果不等于,转第二步;如果等于,结束。
传统获取前向递推概率的过程有两个缺点:一是前向递推概率存储器必须容量大,假设每个前向递推概率值的宽度是B,那么前向递推概率存储器总容量为G×B×N。尤其是随着通信协议的不断发展,N值越来越大,前向递推概率所需要的存储空间也越来越大。第二个缺点是时间开销大,主要是第二步查找G个前向递推概率最大值需要级逻辑来实现,导致总的时间开销是
针对传统计算方法的第一个缺点,相关学者提出滑窗算法以减少前向递推概率存储器容量,如图5所示。滑窗算法将一帧长度为N的待译码数据分成S个段,称之为“窗”,每个窗的长度为W,即N=SW。滑窗算法的流程如图6所示,分为下列步骤。
第一步:初始化窗计算值j=1。
第二步:初始化时刻计数值k=0;将k=0时刻s0状态的前向递推概率设置为0,将k=0时刻s1、s2、......、sG-1状态的前向递推概率 都设置成-2B-1,B是前向递推概率值的宽度,G等于2D,D为Turbo编码器中寄存器单元的个数。
第六步:保存第五步得到的G个防溢出处理后的前向递推概率到前向递推概率存储器。
第七步:k=k+1。
第八步:判断k是否等于W-1,如果不等于,转第三步;如果等于,跳到第九步。
第九步:j=j+1。
第十步:判断j是否等于S,如果不等于,转第二步;如果等于,结束。
采用滑窗算法后,前向递推概率存储总的容量减小为G×B×W,但是前向递推概率计算的时间开销仍然是
发明内容:
针对滑窗算法获取前向递推概率的两个缺点:一是需要的前向递推概率存储容量大,总容量仍为G×B×W,二是前向递推概率获取的时间开销大,总的时间开销是提出一种新的前向递推概率获取方法,既减小前向递推概率存储容量,又减少时间开销。
本发明的技术方案包括以下步骤:
第一步:初始化窗计算值j=1。
第二步:初始化时刻计数值k=0。将k=0时刻s0状态的前向递推概率设置为0,将k=0时刻s1、s2、......、sG-1状态的前向递推概率 都设置成-2B-1,B是前向递推概率值的宽度,G等于2D,D为Turbo译码器相应的Turbo编码器中寄存器单元的个数。
第三步:从窗j的k时刻开始,根据公式一获取k+1时刻G个状态的前向递推概率公式一中的fLUT因子量化成偶数QLUT。
对Turbo译码算法进行简化时,相关学者提出Max-Log-MAP算法,如公式一,s∈{s0、s1、......、sG-1},将获取前向递推概率简化为即fLUT因子等于0,这种方法的误码率性能仍然很高,在现代通信领域应用广泛。而传统的Max*-Log-MAP算法对fLUT因子进行奇偶量化,对误码率性能的提高不是很明显,但是增加了硬件复杂度。本发明采取只对fLUT因子进行偶数量化的方法,误码率性能介于Max-Log-MAP算法和传统的Max*-Log-MAP算法,既保持了高的误码率性能,又不会大幅增加硬件复杂度,对误码率性能需求高的应用场合仍然适用。
从前向递推概率值的宽度B考虑,k+1时刻前向递推概率的范围是-2B-1~2B-1-1,背景技术所述的传统的减去最大值的方法中,最大值也是G个前向递推概率中随机的一个值,两个B位的操作数进行减法操作,需要加入1位的符号保护位,那么减法的结果位数需要B+1位。本发明采取减去s0状态的前向递推概率,s0状态的前向递推概率也是G个前向递推概率中的一个值,只要减法过程中也加入1位保护位,那么本方法就符合Turbo译码算法基于随机概率的本质,通过模拟仿真结果也能证明这一点。
第五步:保存第四步得到的G-1个防溢出处理后的非零状态的前向递推概率到前向递推概率存储器,即保存k+1时刻的非最低位,它们的最低位都是0,不必保存,只需保存B-1位,而等于0,不用保存。在后续计算后验概率的过程中需要用到k+1时刻G个前向递推概率时,将未保存的0值补充完整,补0的过程不会增加硬件开销和计算延迟,也可保证计算的正确性。
第六步:k=k+1。
第七步:判断k是否等于W-1,如果不等于,转第三步;如果等于,跳到第八步。
第八步:j=j+1。
第九步:判断j是否等于S,如果不等于,转第二步;如果等于,结束。
采用本发明可以达到以下技术效果:
1、采用本发明后使得所有时刻s0状态的前向递推概率都等于0,且s1、s2、......、sG-1状态的前向递推概率都是偶数,前向递推概率存储器容量减小为(G-1)×(B-1)×W。
2、比传统计算方法中省略了查找G个前向递推概率最大值的逻辑,既节省了查找电路的面积,又节约了计算的时间。假设一级2选1查找逻辑的延迟是τmax,减法操作的延迟是τsub,那么,采用本发明计算一个时刻的前向递推概率的时间开销只有τsub,比传统方法节省了极大地提高了译码速率。
3、由于Turbo译码的过程中前向递推概率的获取和后向递推概率一样,本发明也可用于后向递推概率的获取。
附图说明
图1为传统的Turbo译码过程;
图2为3GPP LTE协议的Turbo译码(或编码)网格图;
图3为传统获取前向递推概率的流程图;
图4为传统前向递推概率计算逻辑结构;
图5为滑窗算法的Turbo译码过程示意图;
图6为采用滑窗算法获取前向递推概率的流程图;
图7为本发明采用的查找表LUT构建图;
图8为本发明总体流程图;
图9为采用本发明进行译码时的译码错误帧数模拟结果;
图10为采用本发明进行译码时的译码误码率模拟结果。
具体实施方式
如图8所示,本发明总体流程如下:
第一步:初始化窗计算值j=1。
第二步:初始化时刻计数值k=0。将k=0时刻s0状态的前向递推概率设置为0,将k=0时刻s1、s2、......、sG-1状态的前向递推概率 都设置成-2B-1,B是前向递推概率值的宽度,G等于2D,D为Turbo编码器中寄存器单元的个数。
图7中的曲线是fLUT的函数曲线,函数幅值为0.7-0。根据公式一可知:横座标|x-y|表示纵座标表示的函数值。通常定点硬件实现时选取2-q来量化fLUT,其中q是量化的小数位数。q越大,精度越高。图7中带圆圈的折线表示QLUT,量化方式是q=3,即采用[.FFF]的格式,选取的量化点为2-1、2-2、2-3。本发明选择其中的两个点2-1和2-2,即选择0.25和0.5这两个点,从图7可以看出:当|x-y|小于或等于0.125时,函数fLUT被量化成0.25加0.5,等于0.75;当|x-y|大于0.75且小于等于2时,函数fLUT被量化成0.25。QLUT利用下列分段函数表示:
|x-y|<=0.125 →QLUT=2-2+2-1=0.75;
0.125<|x-y|<=0.75 →QLUT=2-1=0.5;
0.75<|x-y|<=2 →QLUT=2-1=0.25;
2<|x-y| →QLUT=0;
采用q=3位小数位,在定点硬件中采用3位二进制表示QLUT={.110,.100,.010,.000}。那么,任意两个值QLUT1和QLUT2的查找表输出之差QLUT1-QLUT2也为偶数。对Turbo译码算法进行简化时,相关学者提出Max-Log-MAP算法,将前向递推概率的计算简化为s∈{s0、s1、......、sG-1},即fLUT因子等于0。而传统的Max*-Log-MAP算法对fLUT因子进行奇偶量化,对误码率性能的提高不是很明显。本发明只对fLUT因子进行偶数量化,因此,本发明的误码率性能介于Max-Log-MAP算法和传统的Max*-Log-MAP算法,对误码率性能需求高的应用场合仍然适用。
第四步:将第三步得到的G个前向递推概率(包括s0状态的前向递推概率)都减去状态s0的前向递推概率即得到G个新的前向递推概率本发明不采用传统的查找最大前向递推概率的方法,使得本发明比传统方法节省了的时间开销,极大地提高了译码速率。
由公式一得到防溢出具体计算过程如下式:
上式中首先是两个找最大值的操作,因此,由公式一可知max()的最终结果仍然还是或是只是选择的状态跳转路径不同。例如公式三中的最大值max()函数的结果不是就是即要么选择的是从状态s1’跳转到状态s的路径,要么选择的是从状态s2’跳转到状态s的路径。
在第二步中,将0时刻s0状态的前向递推概率设置为0,而0时刻s1、s2、......、sG-1状态的前向递推概率都设置成-2B-1,这些值都是偶数。k=1时刻的前向递推概率由k=0时刻的初始值(即0或-2B-1)和分支量度决定。以后每个时刻的前向递推概率都由前一时刻的前向递推概率和分支量度决定,即前向递推概率实际上是初始状态和各个时刻分支量度的叠加。将公式二代入到公式三中的
对于根据公式二都有:例如图2表示的3GPP LTE协议规定的网格图中,从k时刻跳转到k+1时刻的分支路径1/11与分支路径1/00的分支量度之差为:计算得到:任意选取两条分支路径的分支量度之差都等于2乘以一个数,即前向递推概率计算公式中的
上式说明每个时刻0状态的前向递推概率都是0,其他状态的前向递推概率都是偶数。那么,0状态的前向递推概率和其他状态前向递推概率的最低位都不用保存到前向递推概率存储器中。即传统的前向递推概率存储器由深度为N(采用滑窗后,深度为W),宽度为B的G个存储体构成,而采用本发明后,前向递推概率存储器只需要深度为W,宽度为B-1的G-1个存储体,即前向递推概率存储器面积为(G-1)×(B-1)×W。
第五步:保存第四步得到的G-1个防溢出处理后的非零状态的前向递推概率到前向递推概率存储器,即保存k+1时刻的非最低位,它们的最低位都是0,不必保存,即只需保存B-1位,而等于0,不用保存。在后续过程中需要用到k+1时刻G个前向递推概率时,将未保存的0值补充完整,补0的过程不会增加硬件开销和计算延迟,也可保证计算的正确性。
第六步:k=k+1。
第七步:判断k是否等于W-1,如果不等于,转第三步;如果等于,跳到第八步。
第八步:j=j+1。
第九步:判断j是否等于S,如果不等于,转第二步;如果等于,结束。
图9和图10是在Turbo译码中对比采用本发明和采用传统方法后译码性能的模拟结果。图9的纵座标是在高斯白噪声信道中不同信噪比下错误的帧数,图10的纵座标是在不同信噪比下的误码率,图9图10的横座标都是信噪比。模拟环境是3GPP LTE编码,帧长256,帧数10000,迭代5次。从图9和图10看出,采用本发明后译码性能曲线和传统方法的基本重叠,这说明本发明在减小前向递推概率存储容量,提高译码速率的同时可以保证译码的正确性。
Claims (1)
1.一种Turbo译码器中前向递推概率获取方法,其特征在于包括以下步骤:
第一步:初始化窗计算值j=1;
第二步:初始化时刻计数值k=0,将k=0时刻s0状态的前向递推概率设置为0,将k=0时刻s1、s2、......、sG-1状态的前向递推概率 都设置成-2B-1,B是前向递推概率值的宽度,G等于2D,D为Turbo译码器相应的Turbo编码器中寄存器单元的个数。
公式一中的As,k+1是k+1时刻状态s的前向递推概率,s∈{s0、s1、......、sG-1},和分别是k时刻状态s1’和s2’的前向递推概率,s1’和s2’是规约到k+1时刻状态s在k时刻的状态,fLUT是为了补偿公式简化时丢失的精度,称为查找表LUT函数,γs′→s,k是Turbo译码网格图中每条路径的分支量度,计算公式为:
公式二中s′∈{s′1、s′2},系统信息校验信息和是从无线信道接收到的信息,m(Xk)是Turbo译码网格图输入比特Xk经过反相调制后的信息,m(P1k)和m(P2k)是输出比特P1k和P2k经过调制后的信息,即m(Xk),m(P1k),m(P2k)∈{1,-1},Lc是信道置信度,先验信息L(Xk)是Turbo译码时的反馈信息;
第六步:k=k+1;
第七步:判断k是否等于W-1,如果不等于,转第三步;如果等于,跳到第八步;
第八步:j=j+1;
第九步:判断j是否等于S,如果不等于,转第二步;如果等于,结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110072169 CN102270992B (zh) | 2011-03-24 | 2011-03-24 | Turbo译码器中前向递推概率获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110072169 CN102270992B (zh) | 2011-03-24 | 2011-03-24 | Turbo译码器中前向递推概率获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102270992A true CN102270992A (zh) | 2011-12-07 |
CN102270992B CN102270992B (zh) | 2013-01-30 |
Family
ID=45053169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110072169 Expired - Fee Related CN102270992B (zh) | 2011-03-24 | 2011-03-24 | Turbo译码器中前向递推概率获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102270992B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106856406A (zh) * | 2016-11-18 | 2017-06-16 | 北京航空航天大学 | 一种译码方法中校验节点的更新方法及译码器 |
CN114567411A (zh) * | 2022-01-19 | 2022-05-31 | 北京理工大学 | 译码方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1246991A (zh) * | 1997-11-10 | 2000-03-08 | Ntt移动通信网株式会社 | 交织方法、交织装置以及存储交织模式产生程序的媒体 |
US6202178B1 (en) * | 1995-08-21 | 2001-03-13 | Alcatel N.V. | Complementary methods for interleaving and deinterleaving data frames and corresponding forward error correcting devices in a transmitter and receiver |
-
2011
- 2011-03-24 CN CN 201110072169 patent/CN102270992B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6202178B1 (en) * | 1995-08-21 | 2001-03-13 | Alcatel N.V. | Complementary methods for interleaving and deinterleaving data frames and corresponding forward error correcting devices in a transmitter and receiver |
CN1246991A (zh) * | 1997-11-10 | 2000-03-08 | Ntt移动通信网株式会社 | 交织方法、交织装置以及存储交织模式产生程序的媒体 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106856406A (zh) * | 2016-11-18 | 2017-06-16 | 北京航空航天大学 | 一种译码方法中校验节点的更新方法及译码器 |
CN106856406B (zh) * | 2016-11-18 | 2020-05-15 | 北京航空航天大学 | 一种译码方法中校验节点的更新方法及译码器 |
CN114567411A (zh) * | 2022-01-19 | 2022-05-31 | 北京理工大学 | 译码方法、装置、电子设备及存储介质 |
CN114567411B (zh) * | 2022-01-19 | 2024-03-15 | 北京理工大学 | 译码方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN102270992B (zh) | 2013-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104079382B (zh) | 一种基于概率计算的极化码译码器和极化码译码方法 | |
CN108847848B (zh) | 一种基于信息后处理的极化码的bp译码算法 | |
CN109586730B (zh) | 一种基于智能后处理的极化码bp译码算法 | |
CN103269229B (zh) | 一种ldpc-rs二维乘积码的混合迭代译码方法 | |
KR101583139B1 (ko) | 높은 처리량과 낮은 복잡성을 갖는 연속 제거 극 부호 복호 장치 및 그 방법 | |
Wu et al. | A maximum cosinoidal cost function method for parameter estimation of RSC turbo codes | |
CN105634508A (zh) | 一种低复杂度近性能限的Turbo译码器的实现方法 | |
Cyriac et al. | Polar code encoder and decoder implementation | |
CN107659318B (zh) | 一种自适应的极化码译码方法 | |
CN112332864A (zh) | 一种自适应有序移动剪枝列表的极化码译码方法及系统 | |
US10735141B2 (en) | System and a method for error correction coding using a deep neural network | |
CN106209312B (zh) | 一种利用软判决的循环码参数盲识别方法 | |
CN102270992B (zh) | Turbo译码器中前向递推概率获取方法 | |
Huang et al. | Latency reduced method for modified successive cancellation decoding of polar codes | |
CN103220007B (zh) | 一种自适应调整子码不可靠位数的tpc迭代译码算法 | |
Ahmed et al. | Efficient list‐sphere detection scheme for joint iterative multiple‐input multiple‐output detection | |
CN102832954B (zh) | 一种基于软信息平均最小值的Turbo码译码迭代停止方法 | |
RU2699833C1 (ru) | Способ ускоренного декодирования линейного кода | |
CN106059597A (zh) | 一种基于概率Turbo译码器的有符号概率计算单元 | |
CN107733570B (zh) | 基于代数交织器的星座映射方法和映射方式的搜索方法 | |
CN103138769B (zh) | 一种具有不等错误保护的编码方法 | |
CN110190925B (zh) | 一种数据处理方法及装置 | |
CN107682122B (zh) | 一种无线光通信多级编码调制系统的迭代解调译码方法 | |
Bobby et al. | Comparison of fast radix 2 ACS with adaptive fast radix 2 ACS in Viterbi Decoder | |
CN109347600A (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 | ||
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: 20130130 Termination date: 20150324 |
|
EXPY | Termination of patent right or utility model |