CN110212935B - 一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法 - Google Patents

一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法 Download PDF

Info

Publication number
CN110212935B
CN110212935B CN201910375497.8A CN201910375497A CN110212935B CN 110212935 B CN110212935 B CN 110212935B CN 201910375497 A CN201910375497 A CN 201910375497A CN 110212935 B CN110212935 B CN 110212935B
Authority
CN
China
Prior art keywords
noise
bandwidth
signal
frequency
discriminator
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
CN201910375497.8A
Other languages
English (en)
Other versions
CN110212935A (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.)
No63686 Troops Pla
Original Assignee
No63686 Troops Pla
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 No63686 Troops Pla filed Critical No63686 Troops Pla
Publication of CN110212935A publication Critical patent/CN110212935A/zh
Application granted granted Critical
Publication of CN110212935B publication Critical patent/CN110212935B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06Receivers
    • H04B1/16Circuits
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B14/00Transmission systems not characterised by the medium used for transmission
    • H04B14/02Transmission systems not characterised by the medium used for transmission characterised by the use of pulse modulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明涉及一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法,基于统计信号处理技术,在深入分析PCM/FM接收机鉴频器非相干解调原理的基础上,并区分两种不同的情况,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽等因素,结合工程实际进行合理的近似和简化,提出新的误码率计算公式。然后,通过载波定标的方式,确定指定误码率下的信号电平,对设备的误码率进行统计和分析。本发明能够方便快速准确地评估遥测和通信设备性能,便于系统设计和系统评估工作。本发明提出的方法可以推广到航空遥测、移动通信、无线网络等领域,具有较强的普适性。

Description

一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法
技术领域
本发明涉及一种准确评估PCM/FM接收机鉴频器非相干解调性能的新方法,属于信息与通信工程领域。
背景技术
PCM/FM是一种常用的靶场遥测调制体制,自上世纪70年代至今,已有近50年的历史。其调制体制简单而性能优越,地位至今无法替代。
在一般工程应用中,PCM/FM大都采用基于鉴频器的非相干解调方法。目前,在遥测设备的研制和测试过程中,PCM/FM的理论误码率一般按照FSK非相干解调的理论误码率进行计算,即:
Figure GDA0002804992350000011
然而,事实上,这个公式并不符合PCM/FM体制的情况。首先,这个公式只适用于FSK匹配滤波器方式的解调过程。而PCM/FM体制基本均采用基于鉴频器的解调方案。二者有一定区别。上述公式的根本问题是未能体现PCM/FM中调制频偏的影响。在FM体制中,调制频偏越大,解调输出的信噪比也就越高,而上述公式完全忽略了这一因素。
其次,采用鉴频器方案的PCM/FM解调的理论误码特性,在很大程度上依赖解调器的具体设计,如前置滤波器、中频带宽、检后低通滤波器等等。而上述公式未能体现这些影响。实际测试结果也表明:上述公式和实际的理论误码率存在着显著的差异。
本方法基于统计信号处理技术,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽等因素,结合工程实际进行合理的近似和简化,能够准确地评估通信系统的接收解调性能。
发明内容
本发明所要解决的技术问题是针对上述现有技术提供一种准确评估PCM/FM接收机鉴频器非相干解调性能的新方法,能够方面准确地评估遥测设备的误码率性能,从而便于遥测和通信系统的系统设计和系统评估。
本发明解决上述问题所采用的技术方案为:一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法,该方法基于统计信号处理和随机过程分析,建立基于鉴频器的PCM/FM信号接收解调数学模型,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽因素,开发出定量的评估和定标方法,具体包括以下步骤:
步骤一、建立基于鉴频器的PCM/FM信号接收解调数学模型
记中频带通滤波器的带宽为BIF,鉴频器的归一化增益系数为1/2π,低通滤波器的带宽为BL
接收信号的时域表达式为:
xFM(t)=Ac·cos[2πfct+φs(t)]+n(t)
其中:
Ac为接收信号的幅度
fc为接收信号的频率;
φs(t)为数据调制造成的相位动态变化;
m(t)为基带的数据信号,不妨设其幅度为1和-1;
Kf为最大调制频偏;
中频带通滤波器输出的噪声表示为:
n(t)=nI(t)cos2πfct-nQ(t)sin2πfct
=r(t)cos[2πfct+ψ(t)]
其中:
r(t)为噪声的幅度;ψ(t)为噪声的相位,单边带噪声功率谱密度为N0/2,功率为Ni(t)=N0BIF,nI(t)、nQ(t)分别为噪声中的同相分量和正交分量;
步骤二、鉴频器输出噪声功率谱密度的求解
信号和噪声叠加后,总的相位表示为:
Figure GDA0002804992350000021
鉴频器的输出表示为:
xd(t)=(1/2π)dθ(t)/dt
=(1/2π)dφs(t)/dt+nd(t)
=Kfm(t)+nd(t)
其中,nd(t)是鉴频器输出的噪声,表示为:
Figure GDA0002804992350000031
当信噪比较大时,认为
Figure GDA0002804992350000032
在[-π,π]间均匀分布,则上式进一步简化为:
Figure GDA0002804992350000033
由上式可以看出,鉴频器输出的噪声过程是输入噪声过程对时间的微分,则nd(t)的噪声功率谱密度可表示为:
Nd0=(1/Ac)2N0f2
步骤三、根据滤波器带宽选取误码率计算公式
噪声功率的计算分为两种情况:
(1)中频带宽小于或等于两倍低通滤波器带宽
鉴频器输出噪声功率通过对步骤二中获得的噪声功率谱密度积分取得:
Figure GDA0002804992350000034
则:
Figure GDA0002804992350000035
每比特码元的能量Eb由载波功率来计算,即
Figure GDA0002804992350000036
令Rb=1/Tb,h=2Kf/Rb,其中,h为调制指数,Tb为传输数据的码元宽度,Rb为传输数据的码元速率;
则误码率表示为:
Figure GDA0002804992350000037
其中
Figure GDA0002804992350000038
(2)中频带宽大于或等于两倍低通滤波器带宽
鉴频器输出噪声功率同样通过对步骤二中获得的噪声功率谱密度的积分来计算,其积分区间由低通滤波器带宽BL决定:
Figure GDA0002804992350000041
当BL大于或等于1/Tb,满足无码间干扰的要求,则BL=1/Tb,上式重写为:
Nd=(2/3)(1/Ac)2N0(1/Tb)3
Figure GDA0002804992350000042
则误码率可表示为:
Figure GDA0002804992350000043
令L=4/3h2,误码率公式可进一步表示为:
Figure GDA0002804992350000044
步骤四、测试信号定标和测试过程
根据系统测试指定的误码率,由步骤三中的公式计算出所需要的Eb/N0,再根据下列公式计算出载波的信号功率与噪声功率谱密度之比:
S/N0=Eb/N0+10logRb
此处,Rb为传输数据的码元速率,
按上式计算得到的S/N0,加上系统的解调性能指标要求Ds,将相应的信号输入至通信系统或遥测系统中。
与现有技术相比,本发明的优点在于:
1)本发明能够方便快速准确地评估遥测和通信设备性能,便于系统设计和系统评估工作。
2)本发明提出的方法可以推广到航空遥测、移动通信、无线网络等领域,具有较强的普适性。
附图说明
图1为PCM/FM体制鉴频解调原理示意图。
图2为三种调制体制理论误码率比较示意图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本实施例中的一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法,该方法基于统计信号处理和随机过程分析,建立基于鉴频器的PCM/FM信号接收解调数学模型,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽等因素,开发出定量的评估和定标方法,提高了通信系统性能评估的准确性。
所述方法具体包括以下步骤:
步骤一、建立基于鉴频器的PCM/FM信号接收解调数学模型
典型的PCM/FM体制接收机的接收解调原理如图1所示。
信号和噪声叠加后,通过一个中频带通滤波器,滤除带外干扰。经过限幅器,克服信道传输过程中的幅度抖动。经过鉴频器,对信号的相位进行微分(实际是检测频率的变化),再通过低通滤波器,获得原始的调制信号。经判决后,得到原始1和0的信息。
记中频带通滤波器的带宽为BIF,鉴频器的归一化增益系数为1/2π,低通滤波器的带宽为BL
接收信号的时域表达式为:
xFM(t)=Ac·cos[2πfct+φs(t)]+n(t)
其中:
Ac为接收信号的幅度
fc为接收信号的频率;
φs(t)为数据调制造成的相位动态变化。
中频带通滤波器可以视为以fc为中心的窄带滤波器,从而中频带通滤波器输出的噪声可表示为:
n(t)=nI(t)cos2πfct-nQ(t)sin2πfct
=r(t)cos[2πfct+ψ(t)]
其中:
r(t)为噪声的幅度;ψ(t)为噪声的相位,其单边带噪声功率谱密度为N0/2,功率为Ni(t)=N0BIF。nI(t)、nQ(t)分别为噪声中的同相分量和正交分量。
步骤二、鉴频器输出噪声功率谱密度的求解
信号和噪声叠加后,总的相位可表示为:
Figure GDA0002804992350000061
这个式子的表示比较复杂,难于运算和推导。鉴于在实际应用中,一般仅考虑信噪比较大的情况。在这种情况下,上式可近似简化为:
Figure GDA0002804992350000062
鉴频器的作用是对输入信号的相位进行微分,因此鉴频器的输出可表示为:
xd(t)=(1/2π)dθ(t)/dt
=(1/2π)dφs(t)/dt+nd(t)
=Kfm(t)+nd(t)
其中,nd(t)是鉴频器输出的噪声,m(t)为基带的数据信号,设其幅度为1和-1;
Kf为最大调制频偏;
Figure GDA0002804992350000063
上式仍难于求解。为了进一步处理,需进行一定程度的近似。当信噪比较大时,可以认为
Figure GDA0002804992350000064
在[-π,π]间均匀分布,则上式可进一步简化为:
Figure GDA0002804992350000065
由上式可以看出,鉴频器输出的噪声过程是输入噪声过程对时间的微分。微分运算相当于将输入噪声过程通过一个频率响应为j2πf的滤波器。如果输入是正态过程,输出仍为正态过程。则nd(t)的噪声功率谱密度可表示为:
Nd0=(1/Ac)2N0f2
此处值得注意的是,由于输入噪声经过正负频率轴上的两个中频滤波器(在物理上实为一个),此处输入的噪声功率谱密度应该按N0而不是N0/2计算。
由上式可以看出,鉴频器输出的噪声不再是白噪声,其功率谱密度呈抛物线型。这就造成了PCM/FM非相干检测的理论误码公式要比相干检测方式更为复杂。
步骤三、根据滤波器带宽选取误码率计算公式
鉴频器输出的信号可表示为:
Figure GDA0002804992350000071
当传输数据为1时,鉴频器输出的信号符合正态分布,其均值为Kf,方差
Figure GDA0002804992350000077
(Nd为噪声功率),其表示式为:
Figure GDA0002804992350000072
由于受噪声的干扰,当x≤0时,判决器会错误地将1判决为0。错误概率为
Figure GDA0002804992350000073
其中,Q(·)函数的定义为:
Figure GDA0002804992350000074
当数据中1和0的概率相等时,解调过程中的误码率为:
Pe=Q(Kf/σ)
然而,噪声功率的计算却要区分两种情况。
噪声功率的计算区分两种情况。
(1)中频带宽小于或等于两倍低通滤波器带宽
中频滤波器等效到基带之后,对噪声的作用等效为一个区间为[-BIF/2,BIF/2]的低通滤波器。如果这个等效低通滤波器的带宽小于鉴频器后固有的硬件的低通滤波器(区间为[-BL,BL]),则鉴频器的输出噪声功率实际上取决于中频滤波器带宽BIF
鉴频器输出噪声功率可通过对步骤二中获得的噪声功率谱密度的积分来计算:
Figure GDA0002804992350000075
则:
Figure GDA0002804992350000076
由于鉴频器的输出只能反映输入信号相位的变化,无法反映功率的变化,因此,每比特码元的能量Eb不能由Kf计算,只能由载波功率来计算,即
Figure GDA0002804992350000081
令Rb=1/Tb,h=2Kf/Rb,其中,h一般被称为调制指数。
则误码率可进一步表示为:
Figure GDA0002804992350000082
Figure GDA0002804992350000083
则上式表示为:
Figure GDA0002804992350000084
考虑到,FSK相干检测方式的误码率为
Figure GDA0002804992350000085
此处,称L为PCM/FM非相干解调相对于FSK相干解调的恶化因子。为达到某个指定的误码率指标,记PCM/FM非相干解调所需的Eb/N0为Eb/N0|uncoherentCPFSK,记FSK相干解调所需的Eb/N0为Eb/N0|coherentFSK
由于函数
Figure GDA0002804992350000086
较易计算和查表,则在其基础上,将某指定误码率所需的Eb/N0加上10log10L即可得到PCM/FM非相干检测所需的Eb/N0
它们的关系是:
Eb/N0|uncoherentCPFSK=Eb/N0|coherentFSK+10log10L
(2)中频带宽大于或等于两倍低通滤波器带宽
在这种情况下,鉴频器的输出噪声功率实际上取决于低通滤波器带宽BL
鉴频器输出噪声功率同样可通过对步骤二中获得的噪声功率谱密度的积分来计算,但其积分区间由低通滤波器带宽BL决定:
Figure GDA0002804992350000087
假设传输数据的码元宽度为Tb。一般来说,只要BL大于或等于1/Tb,即可满足无码间干扰的要求。此处假定BL=1/Tb,则上式可重写为:
Nd=(2/3)(1/Ac)2N0(1/Tb)3
Figure GDA0002804992350000091
则误码率可表示为:
Figure GDA0002804992350000092
令L=4/3h2,误码率公式可进一步表示为:
Figure GDA0002804992350000093
步骤四、测试信号定标和测试过程
通信系统或遥测系统的解调性能指标一般定义为:为达到某个指定的误码率,所需要的实际信号功率和理论信号功率的差值,一般用分贝(dB)表示。这个差值越小,说明系统的解调性能越优秀。此处用Ds表示这个指标。
根据系统测试指定的误码率,由步骤三中的公式计算出所需要的Eb/N0,再根据下列公式计算出载波的信号功率与噪声功率谱密度之比:
S/N0=Eb/N0+10logRb
此处,Rb为传输数据的码元速率。
按上式计算得到的S/N0,加上系统的解调性能指标要求Ds,将相应的信号输入至通信系统或遥测系统中。如果系统的误码率低于指定的误码率,则表明该系统满足设计要求,如图2所示。
除上述实施例外,本发明还包括有其他实施方式,凡采用等同变换或者等效替换方式形成的技术方案,均应落入本发明权利要求的保护范围之内。

Claims (1)

1.一种准确评估PCM/FM接收机鉴频器非相干解调性能的方法,其特征在于:该方法基于统计信号处理和随机过程分析,建立基于鉴频器的PCM/FM信号接收解调数学模型,综合考虑信道滤波器带宽、噪声功率谱密度、解调滤波器带宽因素,开发出定量的评估和定标方法,具体包括以下步骤:
步骤一、建立基于鉴频器的PCM/FM信号接收解调数学模型
记中频带通滤波器的带宽为BIF,鉴频器的归一化增益系数为1/2π,低通滤波器的带宽为BL
接收信号的时域表达式为:
xFM(t)=Ac·cos[2πfct+φs(t)]+n(t)
其中:
Ac为接收信号的幅度
fc为接收信号的频率;
φs(t)为数据调制造成的相位动态变化;
n(t)中频带通滤波器输出的噪声,表示为:
n(t)=nI(t)cos2πfct-nQ(t)sin2πfct
=r(t)cos[2πfct+ψ(t)]
其中:
r(t)为噪声的幅度;ψ(t)为噪声的相位,单边带噪声功率谱密度为N0/2,功率为Ni(t)=N0BIF,nI(t)、nQ(t)分别为噪声中的同相分量和正交分量;
步骤二、鉴频器输出噪声功率谱密度的求解
信号和噪声叠加后,总的相位表示为:
Figure FDA0002804992340000011
鉴频器的输出表示为:
xd(t)=(1/2π)dθ(t)/dt
=(1/2π)dφs(t)/dt+nd(t)
=Kfm(t)+nd(t)
其中,m(t)为基带的数据信号,设其幅度为1和-1;Kf为最大调制频偏,nd(t)是鉴频器输出的噪声,表示为:
Figure FDA0002804992340000021
当信噪比较大时,认为
Figure FDA0002804992340000022
在[-π,π]间均匀分布,则上式进一步简化为:
Figure FDA0002804992340000023
由上式可以看出,鉴频器输出的噪声过程是输入噪声过程对时间的微分,则nd(t)的噪声功率谱密度可表示为:
Nd0=(1/Ac)2N0f2,其中f为接收机接收到的频率范围;
步骤三、根据滤波器带宽选取误码率计算公式
噪声功率的计算分为两种情况:
(1)中频带宽小于或等于两倍低通滤波器带宽
鉴频器输出噪声功率通过对步骤二中获得的噪声功率谱密度积分取得:
Figure FDA0002804992340000024
则:
Figure FDA0002804992340000025
每比特码元的能量Eb由载波功率来计算,即
Figure FDA0002804992340000026
令Rb=1/Tb,h=2Kf/Rb,其中,h为调制指数,Tb为传输数据的码元宽度,Rb为传输数据的码元速率;
则误码率表示为:
Figure FDA0002804992340000027
其中
Figure FDA0002804992340000031
(2)中频带宽大于或等于两倍低通滤波器带宽
鉴频器输出噪声功率同样通过对步骤二中获得的噪声功率谱密度的积分来计算,其积分区间由低通滤波器带宽BL决定:
Figure FDA0002804992340000032
当BL大于或等于1/Tb,满足无码间干扰的要求,则BL=1/Tb,上式重写为:
Nd=(2/3)(1/Ac)2N0(1/Tb)3
Figure FDA0002804992340000033
则误码率可表示为:
Figure FDA0002804992340000034
令L=4/3h2,误码率公式可进一步表示为:
Figure FDA0002804992340000035
步骤四、测试信号定标和测试过程
根据系统测试指定的误码率,由步骤三中的公式计算出所需要的Eb/N0,再根据下列公式计算出载波的信号功率与噪声功率谱密度之比:
S/N0=Eb/N0+10logRb
此处,Rb为传输数据的码元速率,
按上式计算得到的S/N0,加上系统的解调性能指标要求Ds,将相应的信号输入至通信系统或遥测系统中。
CN201910375497.8A 2018-11-06 2019-05-07 一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法 Expired - Fee Related CN110212935B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2018113145168 2018-11-06
CN201811314516 2018-11-06

Publications (2)

Publication Number Publication Date
CN110212935A CN110212935A (zh) 2019-09-06
CN110212935B true CN110212935B (zh) 2021-05-11

Family

ID=67785582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910375497.8A Expired - Fee Related CN110212935B (zh) 2018-11-06 2019-05-07 一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法

Country Status (1)

Country Link
CN (1) CN110212935B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5711009A (en) * 1996-04-01 1998-01-20 The United States Of America As Represented By The Secretary Of The Navy Method and apparatus for minimizing the effects of frequency selective fading on a PCM/FM data signal
US7512395B2 (en) * 2006-01-31 2009-03-31 International Business Machines Corporation Receiver and integrated AM-FM/IQ demodulators for gigabit-rate data detection
CN101640654B (zh) * 2009-07-27 2011-11-09 北京航空航天大学 一种用于卫星通信系统的超低码速率psk解调器
CN101741786B (zh) * 2009-12-18 2012-12-26 中国人民解放军理工大学 数字通信系统超宽带接收机及其信号处理方法
CN103501203B (zh) * 2013-09-26 2015-03-18 北京空间飞行器总体设计部 一种适用于遥感卫星激光星地通信链路测试系统
CN104486288B (zh) * 2014-12-01 2017-09-15 北京理工大学 一种适用于pcm/fm遥测接收机的载波频偏抑制方法

Also Published As

Publication number Publication date
CN110212935A (zh) 2019-09-06

Similar Documents

Publication Publication Date Title
Wong et al. Specific emitter identification using convolutional neural network-based IQ imbalance estimators
Sills Maximum-likelihood modulation classification for PSK/QAM
US4381546A (en) System for the quantitative measurement of impairments in the communication channel of a quadrature amplitude modulation data communication system
Berni et al. On the utility of chirp modulation for digital signaling
US20090066578A1 (en) System and method for cumulant-based geolocation of cooperative and non-cooperative rf transmitters
Eduardo et al. Experimental evaluation of performance for spectrum sensing: Matched filter vs energy detector
CN104869096B (zh) 基于Bootstrap的BPSK信号盲处理结果可信性检验方法
CN111800359A (zh) 一种通信信号调制方式的识别方法、装置、设备及介质
CN112737992A (zh) 一种水声信号调制方式自适应类内识别方法
Wong et al. Emitter identification using CNN IQ imbalance estimators
CN108683622A (zh) 一种基于信号平方谱平坦系数的调制识别方法
CN110212935B (zh) 一种准确评估pcm/fm接收机鉴频器非相干解调性能的方法
US9106485B1 (en) System and method for FSK demodulation
KR102097367B1 (ko) 직교성을 갖는 동기 워드를 이용한 독립적인 패킷 검출 방법 및 그 수신 장치
US20100265139A1 (en) System and method for cumulant-based geolocation of cooperative and non-cooperative RF transmitters
Maleki et al. Channel-aware m-ary distributed detection: Optimal and suboptimal fusion rules
US7860191B2 (en) Optimal two-layer coherent demodulation for D-PSK (Differential Phase Shift Keying)
Kikuzuki et al. Automatic standard classification method for the 2.4 GHz ISM band
Chen et al. Universal classifier synchronizer demodulator
CN113765838A (zh) 一种dpsk信号的解调方法、装置、设备及存储介质
Song et al. Quantifying the accuracy of high SNR BER approximation of MPSK and MDPSK in fading channels
CN113875202B (zh) 射频信号的灵敏且稳健的帧同步
US20100211990A1 (en) Packet Detection Method for Wireless Communication Device and Related Device
CN115695124B (zh) 一种uqpsk相干解调方法和系统
Mouton Automatic Modulation Classification through clustering in the I/Q plane

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210511

CF01 Termination of patent right due to non-payment of annual fee