CN106059980B - 一种基于快速跳频的多载波扩频方法 - Google Patents
一种基于快速跳频的多载波扩频方法 Download PDFInfo
- Publication number
- CN106059980B CN106059980B CN201610365626.1A CN201610365626A CN106059980B CN 106059980 B CN106059980 B CN 106059980B CN 201610365626 A CN201610365626 A CN 201610365626A CN 106059980 B CN106059980 B CN 106059980B
- Authority
- CN
- China
- Prior art keywords
- spread spectrum
- signal
- ffh
- frequency
- channel
- 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
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000013507 mapping Methods 0.000 claims description 21
- 230000007480 spreading Effects 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 5
- 239000000969 carrier Substances 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 3
- ONUFESLQCSAYKA-UHFFFAOYSA-N iprodione Chemical compound O=C1N(C(=O)NC(C)C)CC(=O)N1C1=CC(Cl)=CC(Cl)=C1 ONUFESLQCSAYKA-UHFFFAOYSA-N 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 27
- 238000005562 fading Methods 0.000 abstract description 13
- 238000004891 communication Methods 0.000 abstract description 7
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000004088 simulation Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2626—Arrangements specific to the transmitter only
- H04L27/2627—Modulators
- H04L27/2628—Inverse Fourier transform modulators, e.g. inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details 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/69—Spread spectrum techniques
- H04B1/713—Spread spectrum techniques using frequency hopping
- H04B1/7143—Arrangements for generation of hop patterns
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2649—Demodulators
- H04L27/265—Fourier transform demodulators, e.g. fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Discrete Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明属于通信技术领域,尤其涉及一种基于多载波直接序列扩频码分多址(MC‑DS‑CDMA)以及FFH‑OFDM的多载波扩频技术。一种基于快速跳频的多载波扩频方法,该方法兼顾了MC‑DS‑CDMA和FFH‑OFDM共同的特点,避免了传统的MC‑DS‑CDMA无频率分集增益的弱点,在抗衰落、抗干扰、码分多址、通讯保密性方面都具有一定的优势。
Description
技术领域
本发明属于通信技术领域,尤其涉及一种基于多载波直接序列扩频码分多址(MC-DS-CDMA)以及FFH-OFDM的多载波扩频技术。
背景技术
随着OFDM和扩频通信的发展,由此出现了一种全新通信技术——多载波扩频技术。该技术将OFDM技术和扩频技术结合在了一起,很好的兼顾了两种技术在抗多径、抗干扰、码分多址、保密性方面的优势,成为目前通信技术研究的热点。根据其结构和扩频形式的不同可以大致划分成四类型:MC-CDMA、MC-DS-CDMA、MT-CDMA、二维扩频技术。其中MC-DS-CDMA是一种时域扩频技术,该技术首先对输入信号进行串并变换,然后在时间域上实现扩频,把扩频输出结果分别调制到不同的正交载波上。从其扩频方式可以看出,由于MC-DS-CDMA是一种在时域上进行的多载波扩频技术,因此在不进行前向纠错的情况下,无频域分集增益。
FFH-OFDM技术作为一种全新的OFDM衍生技术,从形式上来看,FHH将传统OFDM单一子载波上传输的数据扩散到所有子载波上进行传输,克服了传统OFDM无频率分集的弱点。从本质上看,FFH是一种预编码,通过这种预编码实现了对DFT变换矩阵因子的重新排序。并行信号中的任意一路信号遍历了所有的子载波,实现了在OFDM符号周期内信号的快速跳变。由于FFH-OFDM具有以上的特点,使其可以对来自各子信道的衰退和噪声干扰进行均匀的混合,在频率选择性衰退和高斯白噪的环境中性能优于传统的OFDM。但是普通的FFH-OFDM不包含扩频,因此与传统的OFDM相比在码分多址、保密性方面、抗干扰方面并不具有优势。
发明内容
本发明在MC-DS-CDMA的基础上结合FFH-OFDM提出了一种基于快速跳频的多载波扩频方法,该方法兼顾了MC-DS-CDMA和FFH-OFDM共同的特点,避免了传统的MC-DS-CDMA无频率分集增益的弱点,在抗衰落、抗干扰、码分多址、通讯保密性方面都具有一定的优势。
一种基于快速跳频的多载波扩频方法,具体为:
S1、时域直接扩频,具体为:已调信号首先经过串并变换后变为多路并行的信号,对所述多路并行信号用相同的扩频序列进行直接扩频得到时域扩频后的信号其中,dk(t)为时域信号,N为扩频序列长度,an为扩频序列,M为并行子载波个数,gk(t)为时域直接扩频后的输出信号,所述gk(t)的集合G(t)=(g1(t),g2(t),g3(t),...,gM(t))T,为脉冲信号,其持续周期为Tc,n=1,...,N,k=1,…,M;
S2、进行FFH映射,具体为:
S21、S1所述G(t)进入FFH映射模块后,FFH映射模块根据预先设定好的跳频矩阵对信号进行快跳频得到
将S1所述gk(t)代入所述zv(t)中得到其中,uq,k为FFH跳频矩阵UM中第q行第k列的元素,所述v表示采样时间点,u表示输入数据序号,k=1,…,M,v=1,…,M,所述fM为不同时间点的起始瞬时频率,FM=(f1,f2,…,fM)T为(1,2,…,M)这M个整数的乱序排列,*H表示共轭转置,DM为M行M列的离散傅里叶变换矩阵,
为DM中第q行第k列的元素;
S22、利用S221所述FM产生出不相同的FFH矩阵和IFFH矩阵,完成不一样的跳频映射;
S3、进行多载波正交调制,将信号zv(t)通过离散傅里叶反变换(IDFT)调制到相互正交子载波上得到f(t),具体为:
设任意两个正交子载波之间的频率间隔为△f=1/Tc,有M个正交子载波,该系统的频带宽度为M△f,则IDFT输出信号f(t)可以表示为
S4、进行多载波正交解调,具体为:
在接收端,设信号已经经过了精确的信号同步和信道均衡,同时接收信号不受射频中频模块的影响,可得到经过离散傅里叶变换后的信号为
其中,n(t)为均值为0,方差为δ2的高斯白噪声,该噪声是由于在无线信道中传输引起的,Hp为第p个子载波上的信道频率响应;
S5、信道均衡,即信道经过了准确的估计,得到信道频率响应的Hp,经过迫零信道均衡以后可以得到
S6、进行IFFH映射,具体为:将S5所述送入IFFH映射模块进行解调,IFFH映射模块输出为其中,是IFFH解调矩阵为第k行第q列的元素,将S4所述代入S5所述送入IFFH映射模块进行解调得到:
S7、时域直接扩频的解扩,具体为:将S6所述进行直接扩频的解扩,利用与发送端事先约定的扩频码进行相关解扩,得最终的接收信号
本发明的有益效果是:
具有二维扩频特性:MC-DS-CDMA通过直接扩频实现了在时间域上的扩频,是一种时域多载波扩频技术。FFH-OFDM通过FFH预编码实现了在所有子载波上的快速跳频,是一种频域上的扩频。将两种技术结合起来,实现了时频域上的双扩频。
抗多径能力强:改进的多载波扩频技术利用FFH实现了频域上的扩频,在一个符号周期内信号遍历性的在所有正交子载波上进行跳频,这种符号周期内的快速跳频可以认为是一种频率分集,FFH将所有子信道的衰落和噪声进行均匀的混合,在抗多径频率选择性衰退方面具有优良的性能。
抗干扰能力强:基于快速跳频的多载波扩频技术采用直接扩频的方式实现了时域上的扩频,因此具有直接扩频技术抗干扰能力强的特点。同时,在技术中信号遍历了所有正交子载波。当其中某一个载波受到窄带噪声干扰,可以丢弃受到干扰载波上携带的信息,利用剩余载波上的分集信息进行信号提取,可以大大提高还技术抗干扰的能力。
具有码分多址能力:基于快速跳频的多载波扩频技术是在MC-DS-CDMA基础上发展起来的,因此具有码分多址的能力,可以实现多用户共享同一频带资源。
保密性强:基于快速跳频的多载波扩频技术中包含直接扩频模块,在未知扩频码的情况下很难恢复出传输信号,因此信号在传输过程中不易被截获、窃听或侦查,具有很好的保密性和安全性。
附图说明
图1为改进的多载波扩频技术原理图。
图2为平坦衰落的高斯信道下误码率随信噪比变换曲线。
图3为五条可分辨路径的多径高斯信道下误码率曲线。
图4为九条可分辨路径的多径高斯信道下误码率曲线。
具体实施方式
下面结合附图,对本发明作进一步地详细描述。
本发明在原有的MC-DS-CDMA技术基础上,结合FFH-OFDM技术提出了一种基于快速跳频的多载波扩频技术,该方案的具体步骤如下所示:
发射机部分:
步骤一、对原始信号进行调制,调制方式可以为PSK或者QAM等。
步骤二、串变并,将已调信号进行串并转换,变为M路并行信号。
步骤三、时域直接扩频,采用相同的扩频序列对M路并行信号采用相同的扩频码分别进行直接扩频,如所示。
步骤四、FFH映射,FFH映射模块根据预先设定好的跳频矩阵对信号进行快跳频,映射方式如所示。
步骤五、插入导频,在这里采用的是块状导频。
步骤六、多载波正交调制,利用离散傅里叶反变换(IDFT)将经过FFH跳频以后的信号分别调制到相互正交子载波上,调制方式如所示。
步骤七、并变串,并加入循环前缀(CP),通过射频模块经天线辐射到无线信道中。
接收机部分:
步骤一、去掉CP,串变并,送入IDFT模块。
步骤二、多载波正交解调,对接收信号进行离散傅里叶变换(DFT),得到M路并行信号,如所示。
步骤三、信道均衡,利用导频信号对信道进行估计得到信道的频率响应为Hv,采用迫零法对DFT变换输出信号进行信道均衡,如
所示。
步骤四、IFFH映射,IFFH映射模块根据预先设定好的解调矩阵对信号进行快跳频的解调,映射方式如所示。
步骤五、时域直接扩频的解扩,利用与发送端事先约定的扩频码进行相关解扩。解扩方式如所示。
步骤六、解调,按照PSK或者QAM的解调规律进行接收信号的解调。
仿真一、为了分析改进的效果,在这里利用MATLAB仿真平台进行了验证。选用了多径高斯信道模型,其功率谱延迟包络服从负指数分布。仿真过程中,采用FFH-OFDM和MC-DS-CDMA作为对照组。系统的参数设定为:符号间隔为0.05us,子载波个数为256,CP长度为16,CP时长为0.8us,OFDM符号长度为13.6us,信号的有效带宽为20MHZ,调制方式选用QPSK,改进的多载波扩频技术使用WASH码作为扩频码,扩频码长度为16。
图2为平坦衰落的高斯信道下误码率随SNR的变换曲线,此时信道的最大延迟时间小于符号周期,接收信号中不会出现ISI。从图中可以看出,改进技术与MC-DS-CDMA的性能相似,前者略优于后者,当SNR大于8dB以后误码率都接近于零。而FFH-OFDM的性能最差,当SNR大于12dB以后误码率才接近于零。因此,在平坦衰落的环境中,改进技术的性能优于MC-DS-CDMA以及FFH-OFDM。
图3为多径高斯白噪声信道下误码率随SNR的变换曲线,其中可分辨路径数为5条,包括1条期望路径和4条反射路径。相邻两条路径的时延间隔为0.15us,最大时延为0.6us。模拟了频率选择性衰落情况下系统的情况,此时多径时延没有超过CP长度(CP时间长度为0.8us)。从图中可以看出,频率选择性衰落使得三种技术的误码率都有所上升。由于FFH对各信道的衰落具有混合的作用,改进技术的优势显现出来,当SNR高于12dB时误码率趋于零,明显优于MC-DS-CDMA。同时还可以发现当SNR大于16dB时,FFH-OFDM的误码率也低于MC-DS-CDMA。说明FFH跳频映射在频率选择性衰落的信道下性能更好。
图4为多径高斯白噪声信道下误码率随SNR的变换曲线,其中可分辨路径数为9条,包括1条期望路径和8条反射路径。相邻两条路径的时延间隔为0.2us,最大时延为1.6us。此时多径时延超过了CP长度,CP无法起到保护作用,会造成ICI的产生。当SNR较低时,误码率随着噪声功率的下降而下降。但当SNR高于4dB时,由于ICI的影响,误码率不随信噪比的升高而发生变化,在10-2附近保持不变。就总体而言,改进技术的误码率依然低于FFH-OFDM和MC-DS-CDMA。
Claims (1)
1.一种基于快速跳频的多载波扩频方法,其特征在于,包括如下步骤:
S1、时域直接扩频,具体为:已调信号首先经过串并变换后变为多路并行的信号,对所述多路并行信号用相同的扩频序列进行直接扩频得到时域扩频后的信号其中,dk(t)为时域信号,N为扩频序列长度,an为扩频序列,M为并行子载波个数,gk(t)为时域直接扩频后的输出信号,所述gk(t)的集合G(t)=(g1(t),g2(t),g3(t),...,gM(t))T,为脉冲信号,其持续周期为Tc,n=1,...,N,k=1,...,M;
S2、进行FFH映射,具体为:
S21、S1所述G(t)进入FFH映射模块后,FFH映射模块根据预先设定好的跳频矩阵对信号进行快跳频得到将S1所述gk(t)代入所述zv(t)中得到
其中,uq,k为FFH跳频矩阵UM中第q行第k列的元素,所述v表示采样时间点,u表示输入数据序号,k=1,...,M,v=1,...,M,所述fM为不同时间点的起始瞬时频率,FM=(f1,f2,...,fM)T为(1,2,...,M)这M个整数的乱序排列,*H表示共轭转置,DM为M行M列的离散傅里叶变换矩阵, 为DM中q行第k列的元素;
S22、利用S21 所述FM产生出不相同的FFH矩阵和IFFH矩阵,完成不一样的跳频映射;
S3、进行多载波正交调制,将信号zv(t)通过离散傅里叶反变换(IDFT)调制到相互正交子载波上得到f(t),具体为:
设任意两个正交子载波之间的频率间隔为△f=1/Tc,有M个正交子载波,该系统的频带宽度为M△f,则IDFT输出信号f(t)可以表示为
S4、进行多载波正交解调,具体为:
在接收端,设信号已经经过了精确的信号同步和信道均衡,同时接收信号不受射频中频模块的影响,可得到经过离散傅里叶变换后的信号为
其中,n(t)为均值为0,方差为δ2的高斯白噪声,该噪声是由于在无线信道中传输引起的,Hp为第p个子载波上的信道频率响应;
S5、信道均衡,即信道经过了准确的估计,得到信道频率响应的Hp,经过迫零信道均衡以后可以得到
S6、进行IFFH映射,具体为:将S5所述送入IFFH映射模块进行解调,IFFH映射模块输出为其中,是IFFH解调矩阵为第k行第q列的元素,将S4所述代入S5所述送入IFFH映射模块进行解调得到:
S7、时域直接扩频的解扩,具体为:将S6所述进行直接扩频的解扩,利用与发送端事先约定的扩频码进行相关解扩,得最终的接收信号
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610365626.1A CN106059980B (zh) | 2016-05-27 | 2016-05-27 | 一种基于快速跳频的多载波扩频方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610365626.1A CN106059980B (zh) | 2016-05-27 | 2016-05-27 | 一种基于快速跳频的多载波扩频方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106059980A CN106059980A (zh) | 2016-10-26 |
CN106059980B true CN106059980B (zh) | 2019-03-29 |
Family
ID=57175979
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610365626.1A Expired - Fee Related CN106059980B (zh) | 2016-05-27 | 2016-05-27 | 一种基于快速跳频的多载波扩频方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106059980B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106998219B (zh) * | 2017-05-10 | 2019-04-23 | 厦门大学 | 基于N选m组合选择器调制的多载波水声跳频通信系统 |
CN112152663B (zh) * | 2020-09-25 | 2021-08-31 | 中国电子科技集团公司第五十四研究所 | 一种基于时分多址的抗多径抗干扰信号接收方法 |
CN112152664B (zh) * | 2020-09-25 | 2021-07-27 | 中国电子科技集团公司第五十四研究所 | 一种基于时分多址的抗多径抗干扰信号发送方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1758576A (zh) * | 2005-09-29 | 2006-04-12 | 北京交通大学 | 一种空时分组码mc-ds-cdma系统上行链路发射与接收方法 |
CN101534279A (zh) * | 2009-04-10 | 2009-09-16 | 北京大学 | 直接序列扩频级联频域正交码分复用收发装置、方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040208624A1 (en) * | 2002-04-05 | 2004-10-21 | Universite Laval | Fast frequency hopping spread spectrum for code division multiple access communications networks (FFH-CDMA) |
WO2006118433A1 (en) * | 2005-05-04 | 2006-11-09 | Samsung Electronics Co., Ltd. | Adaptive data multiplexing method in ofdma system and transmission/reception apparatus thereof |
US9294240B2 (en) * | 2011-02-25 | 2016-03-22 | Qualcomm Incorporated | Multi-carrier operations with fast frequency hopping |
-
2016
- 2016-05-27 CN CN201610365626.1A patent/CN106059980B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1758576A (zh) * | 2005-09-29 | 2006-04-12 | 北京交通大学 | 一种空时分组码mc-ds-cdma系统上行链路发射与接收方法 |
CN101534279A (zh) * | 2009-04-10 | 2009-09-16 | 北京大学 | 直接序列扩频级联频域正交码分复用收发装置、方法及系统 |
Non-Patent Citations (1)
Title |
---|
Performance Evaluation of Fast Frequency Hopping OFDM over UWB Channels;Poramate Tarasak等;《IEEE》;20091231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106059980A (zh) | 2016-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105530217B (zh) | 基于加权分数傅里叶变换的gfdm系统的信号发射和接收方法 | |
Choi et al. | On channel estimation and detection for multicarrier signals in fast and selective Rayleigh fading channels | |
WO2003081823A1 (en) | Cdma system with frequency domain equalization | |
CN101534278B (zh) | 时频扩展的正交频分复用收发装置、方法及系统 | |
CN101394385B (zh) | 基于时域处理联合信道估计的提高正交频分复用系统性能的方法 | |
CN103428154A (zh) | 基于Vector OFDM的双选择性信道的变换域复用方法 | |
CN106059980B (zh) | 一种基于快速跳频的多载波扩频方法 | |
CN102255707A (zh) | 针对sc-fdma mu- mimo系统的基于互信息的空间频率调度方法 | |
CN103647740B (zh) | 基于拉马努金和的正交非均匀多载波间距的多载波调制解调方法 | |
CN112822133B (zh) | 一种多用户正交频分复用差分混沌键控通信方法 | |
Sun et al. | MMSE detection for space-time coded MC-CDMA | |
CN103780529B (zh) | 通信系统及其信号发送方法与装置、信号接收方法与装置 | |
Sohl et al. | Comparison of localized, interleaved and block-interleaved FDMA in terms of pilot multiplexing and channel estimation | |
Haab et al. | Filter bank multi-carrier spread spectrum with biorthogonal signaling for high speed data transmission through HF skywave channels | |
CN102487368A (zh) | Per-tone均衡器的设计方法及实现装置 | |
CN101789810A (zh) | 最小均方误差两级扩频ofdm自适应码分配方法及系统 | |
Kushwah et al. | Performance Estimation of 2* 2 MIMO-MC-CDMA Using Convolution Code in Different Modulation Technique | |
Banelli et al. | Spectral efficiency of multicarrier schemes for 5G | |
Yu et al. | Performance of an MC-CDMA system based on the complex wavelet packet and pre-equalisation technique | |
Ghanim et al. | Achievement of novel multicarrier cdma with various modulation techniques | |
Domínguez-Bolaño et al. | Experimental evaluation of 5G modulation schemes in quasi-static scenarios | |
Ghanem et al. | Image Transmission on MC-CDMA System over Rayleigh and AWGN Channel–comparison study | |
CN101753501A (zh) | 一种td-scdma下行链路结合正交频分复用的方法 | |
Smirani et al. | An Uplink LTE-A Channel Estimation Method Based On Connexionist System | |
Chang et al. | A novel architecture of single-carrier block transmission DS-CDMA |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20190329 |
|
CF01 | Termination of patent right due to non-payment of annual fee |