CN102023010B - 基于mems的小波域多传感器信息融合系统及融合方法 - Google Patents
基于mems的小波域多传感器信息融合系统及融合方法 Download PDFInfo
- Publication number
- CN102023010B CN102023010B CN2010105228251A CN201010522825A CN102023010B CN 102023010 B CN102023010 B CN 102023010B CN 2010105228251 A CN2010105228251 A CN 2010105228251A CN 201010522825 A CN201010522825 A CN 201010522825A CN 102023010 B CN102023010 B CN 102023010B
- Authority
- CN
- China
- Prior art keywords
- sigma
- wavelet
- signal
- module
- scale
- 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
Images
Landscapes
- Gyroscopes (AREA)
Abstract
本发明公开的一种基于MEMS的小波域多传感器信息融合系统及融合方法,包括依次连接的微机电陀螺阵列模块、小波域多尺度分析模块、多尺度多传感器信息融合模块、小波多尺度重构模块及组合微机电惯性测量模块,微机电陀螺阵列模块还通过数据统计性分析模块与多尺度多传感器信息融合模块相连接。通过对各个微机电陀螺的信息进行多尺度小波分解和统计分析,然后在多尺度上进行融合,最后将融合后各尺度上的微机电陀螺信息重构到原始尺度上。使微机电测量设备的精度有效提高,各微机电陀螺仪的协调性得以应用,同时也增强了该系统的可靠性。融合系统性能稳定、工作可靠、体积小、性价比高,可以为各种设备提供精确的角度测量。
Description
技术领域
本发明属于信号处理技术领域,具体涉及一种基于MEMS的小波域多传感器信息融合系统,本发明还涉及基于该系统的融合方法。
背景技术
目前国内外都在研究微机电陀螺,以及微机电惯性测量设备,微机电惯性测量设备主要是由多个微机电测量器件(陀螺仪、加速度计)通过信息融合技术组合而成的。
(1)微机电陀螺仪
微机电陀螺仪(Micro-Electro Mechanical System Gyroscope,MEMSGyro),用于敏感模拟坐标系相对理想坐标系的偏角或角速度,是各类惯性系统中的核心部件。由于陀螺仪的广阔应用前景,国内外都对MEMS陀螺进行了大量的研究工作,在多轴集成MEMS陀螺方面,国内外也很受重视,但迄今尚未有成熟的样机和产品。
(3)信息融合
信息融合是关于协同利用多传感器信息,涉及多级别、多方面、多层次信息检测、相关、估计和综合的一种多级别多自动智能信息处理的过程。它利用计算机技术对按时序获得的多个传感器的信息或已处理的数据在一定准则下加以分析、处理和综合,可以得到单一传感器所不能得到的信息,目的是获取更精确的信息。多传感器系统是信息融合技术的硬件基础,多传感器所给出的信息是融合系统的加工对象,对多源化信息的协调优化是融合的核心,通常所说的融合指的是多传感器信息融合。
(2)小波域多尺度分析
对具体的融合系统而言,它所接收的数据可以是单一层次上的信息,也可以是几个层次上的信息。融合的基本策略是同一层次上的信息进行融合,从而获得更高层次的融合信息,然后再汇入相应的信息融合层次。总的来说,信息融合本质上是一种由低层至顶层对多源信息进行融合,逐层抽象的信息处理过程。小波分析具有多分辨率特性,可以给同一层次上的信息进行多尺度分解,得到多层次上的信息。
现在大多MEMS陀螺的产品都没有加入信息融合技术,而是对陀螺的输出数据进行一些滤波处理,这种方法对精度的提高并不大,而且由于陀螺的随机误差不具备简单滤波所需要的条件,所以滤波效果不佳。本发明的途径在于采用多个单轴的MEMS陀螺集成高精度的微机电测量系统,该系统的输出可以看作是虚拟的MEMS陀螺仪,该系统的精度较原陀螺仪来说可以提高一个量级。
发明内容
本发明的目的是提供一种基于MEMS的小波域多传感器信息融合系统,解决了现有滤波系统精度低,滤波效果不佳的问题,使得系统融合的精度更高,稳定性更好。
本发明的另一目的是提供一种基于MEMS的小波域多传感器信息融合方法。
本发明所采用的技术方案是,一种基于MEMS的小波域多传感器信息融合系统,包括依次连接的微机电陀螺阵列模块、小波域多尺度分析模块、多尺度多传感器信息融合模块、小波多尺度重构模块及组合微机电惯性测量模块,微机电陀螺阵列模块还通过数据统计性分析模块与多尺度多传感器信息融合模块相连接。
本发明所采用的另一技术方案是,一种基于MEMS的小波域多传感器信息融合方法,采用基于MEMS的小波域多传感器信息融合系统,其结构为:包括依次连接的微机电陀螺阵列模块、小波域多尺度分析模块、多尺度多传感器信息融合模块、小波多尺度重构模块及组合微机电惯性测量模块,微机电陀螺阵列模块还通过数据统计性分析模块与多尺度多传感器信息融合模块相连接,
具体按照以下步骤实施:
步骤1:微机电陀螺阵列模块产生各个微机电陀螺仪的信号,传递给小波域多尺度分析模块及数据统计性分析模块,一方面,小波域多尺度分析模块对各个微机电陀螺仪的信号进行小波多尺度分析,得到第i+1-N尺度上的细节信号和第N层尺度上的近似信号,传递给多尺度多传感器信息融合模块;另一方面,数据统计性分析模块对收集到的各个微机电陀螺仪的信号进行统计分析,计算得到均值和标准差信号,传递给多尺度多传感器信息融合模块;
步骤2:多尺度多传感器信息融合模块将收集到的由小波域多尺度分析模块输入的第i+1-N尺度上的细节信号和第N层尺度上的近似信号,以及数据统计性分析模块输入的均值和标准差信号,进行多尺度多传感器信息融合,得到融合后的精确多尺度上的微机电陀螺仪信号,传递给小波多尺度重构模块;
步骤3:小波多尺度重构模块对步骤2得到的融合后的精确多尺度上的微机电陀螺仪信号进行小波多尺度重构,得到原始尺度上的微机电陀螺信号,传递给组合微机电惯性测量设备模块;
步骤4:组合微机电惯性测量设备模块将步骤3得到的原始尺度上的微机电陀螺信号整理成需要的形式输出,完成基于MEMS的小波域多传感器信息融合。
本发明的有益效果是,
(1)小波域多传感器信息融合方法,使微机电测量设备的精度有效提高,各微机电陀螺仪的协调性得以应用,同时也增强了该系统的可靠性。
(2)基于微机电陀螺的小波域多传感器信息融合系统性能稳定、工作可靠、体积小、性价比高,可以为各种设备提供精确的角度测量。
附图说明
图1是本发明基于MEMS的小波域多传感器信息融合系统的结构示意图;
图2是本发明融合方法中小波多尺度分析的工作流程图。
图中,1.微机电陀螺阵列模块,2.小波域多尺度分析模块,3.数据统计性分析模块,4.多尺度多传感器信息融合模块,5.小波多尺度重构模块,6.组合微机电惯性测量模块。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明基于MEMS的小波域多传感器信息融合系统的结构,如图1所示,包括依次连接的微机电陀螺阵列模块1、小波域多尺度分析模块2、多尺度多传感器信息融合模块4、小波多尺度重构模块5及组合微机电惯性测量模块6,微机电陀螺阵列模块1还通过数据统计性分析模块3与多尺度多传感器信息融合模块4相连接。
其中的,微机电陀螺阵列模块1,用于测量各个微机电陀螺的信息。小波域多尺度分析模块2,用于将各个微机电陀螺的信号进行多尺度小波分解。数据统计性分析模块3,用于将得到的各个微机电陀螺信号进行统计分析,得到均值和方差等统计特性。多尺度多传感器信息融合模块4,用于将多个微机电陀螺仪的信号在多尺度上进行融合。小波多尺度重构模块5,用于将融合后各尺度上的微机电陀螺信息重构到原始尺度上。组合微机电惯性测量模块6,用于按要求测量输出组合微机电陀螺的信号。
本发明基于MEMS的小波域多传感器信息融合方法,具体按照以下步骤实施:
步骤1:微机电陀螺阵列模块1产生各个微机电陀螺仪的信号,传递给小波域多尺度分析模块2及数据统计性分析模块3,
一方面,小波域多尺度分析模块2对各个微机电陀螺仪的信号进行小波多尺度分析,分析的过程如图2所示,具体按照以下步骤实施:
在尺度i+1上,对于输入陀螺仪的信号序列x(i+1,k)∈Rn×1,(k∈Z),其离散小波变换的分析形势和综合形式分别为:
其中,h(2k-l)为小波变换中的尺度系数,g(2k-l)为小波变换中的小波系数,选择合适的小波基函数,就可以得到相应的尺度系数和小波系数,由公式(1)可以得到尺度i上陀螺仪信息中的近似信号x(i,k)和细节信号d(i,k)。继续这个过程对尺度i上陀螺仪的近似信号x(i,k)进行离散小波分解,可以得到尺度i-1上陀螺仪信息中的近似信号x(i-1,k)和细节信号d(i-1,k),重复到最优的分解尺度N上,可以得到小波离散分解的近似信号x(N,k)和细节信号d(N,k)。
然后将多尺度分析的结果:第(i+1-N)尺度上的细节信号和第N层尺度上的近似信号送到多尺度多传感器信息融合模块4等待处理。
另一方面,数据统计性分析模块3将收集到的各个微机电陀螺仪的信号进行统计分析,计算出均值、标准差信号到多尺度多传感器信息融合模块4等待处理。
步骤2:多尺度多传感器信息融合模块4将收集到的由小波域多尺度分析模块2输入的第(i+1-N)尺度上的细节信号和第N层尺度上的近似信号,以及数据统计性分析模块3输入的均值和标准差信号,按照以下算法进行多尺度多传感器信息融合,得到融合后的精确多尺度上的微机电陀螺仪信号,传递给小波多尺度重构模块5:
假设得到不同尺度上的多传感器数据,这类数据可以视为被测量信息及干扰信息迭加的结果。传感器系统中,每个传感器的观测都是一致的。假设,传感器的测量模型:
Xi=X+Vi,i=1,2,…,m (2)
其中,X表示要估计的传感器的真值,Vi为均值为0,标准差为有限值的测量噪声,其与传感器的测量误差模型有关。近似认为传感器测量数据服从高斯分布。m个传感器的标准差分别为各传感器的加权因子分别为归一化权值W1,W2,…,Wi,该权值一般取为各个传感器的标准差满足:
融合后的X值和加权因子满足:
则由X1,X2,…,Xi彼此独立,时域融合后序列的标准差:
从上式可以看出,总的均方差是关于各加权因子的多元二次函数,因此必然存在最小值。该最小值的求取是加权因子W1,W2,…,Wi满足约束条件的多元函数极值求取。
根据多元函数求值理论,可求出总均方误差最小时所对应的加权因子为:
时域加权融合的最小均方误差为:
时域融合后的方差小于各传感器的方差。时域融合后的估计是各传感器测量值的线性函数。根据单一尺度上的最优加权信息融合公式,结合小波域随机信号的重构公式为:
XJ是经过J层分解后重构回去的融合数据,是第J层的尺度函数,ψj,k是从1到J层的小波函数。hn(J,k)是第J层上k个近似系数,gn(j,k)是1到J层的k个细节系数。经过J层分解后重构回去的融合结果:
式中,各传感器相应近似信号和细节信号的加权因子分别为Whi和Wgi,同理有权值的归一化条件。根据小波分解系数之间的相关性结论式等。hn,gn是独立不相关的,多尺度数据融合重构序列的方差:
根据相关函数和均方差之间的关系可知,相关函数值等于均方差加均值。这里均值取为相应变量的无偏估计。同样关于各加权因子的多元二次函数,结合加权因子的约束条件为:
根据拉格朗日方程组,求取多元函数的极值可得总均方误差最小时所对应的加权因子为:
此时所对应的最小均方误差为:
步骤3:小波多尺度重构模块5对步骤2得到的融合后的精确多尺度上的微机电陀螺仪信号进行小波多尺度重构,按照如图2小波域多尺度分析的反向过程,具体过程如下:
式中,x(i,l)和d(i,l)分别是尺度i上的近似信号和细节信号,h(2l-k)和g(2l-k)分别是尺度i上的尺度系数和小波系数,x(i+1,k)是经过小波重构得到的尺度i+1上的近似信号,即尺度i+1上陀螺仪的近似信息。对各个尺度上的信号进行多尺度的重构,得到原始尺度上的微机电陀螺信号,将其传送给组合微机电惯性测量设备模块6。
步骤4:组合微机电惯性测量设备模块6将步骤3得到的原始尺度上的微机电陀螺信号整理成需要的形式输出,以为后续设备提供高精度,高稳定性的组合微机电陀螺测量的信息,完成基于MEMS的小波域多传感器信息融合。
Claims (1)
1.一种基于MEMS的小波域多传感器信息融合方法,其特征在于,采用基于MEMS的小波域多传感器信息融合系统,其结构为:包括依次连接的微机电陀螺阵列模块(1)、小波域多尺度分析模块(2)、多尺度多传感器信息融合模块(4)、小波多尺度重构模块(5)及组合微机电惯性测量模块(6),所述的微机电陀螺阵列模块(1)还通过数据统计性分析模块(3)与多尺度多传感器信息融合模块(4)相连接,
具体按照以下步骤实施:
步骤1:微机电陀螺阵列模块(1)产生各个微机电陀螺仪的信号,传递给小波域多尺度分析模块(2)及数据统计性分析模块(3),一方面,小波域多尺度分析模块(2)对各个微机电陀螺仪的信号进行小波多尺度分析,得到第i+1-N尺度上的细节信号和第N层尺度上的近似信号,传递给多尺度多传感器信息融合模块(4);另一方面,数据统计性分析模块(3)对收集到的各个微机电陀螺仪的信号进行统计分析,计算得到均值和标准差信号,传递给多尺度多传感器信息融合模块(4);
所述的小波域多尺度分析模块(2)对各个微机电陀螺仪的信号进行小波多尺度分析,具体按照以下步骤实施:
在尺度i+1上,对于输入陀螺仪的信号序列x(i+1,k)∈Rn×1,(k∈Z),其离散小波变换的分析形势和综合形式分别为:
其中,h(2k-l)为小波变换中的尺度系数,g(2k-l)为小波变换中的小波系数,选择合适的小波基函数,得到相应的尺度系数和小波系数,由上式得到尺度i上陀螺仪信息中的近似信号x(i,k)和细节信号d(i,k),继续这个过程对尺度i上陀螺仪的近似信号x(i,k)进行离散小波分解,得到尺度i-1上陀螺仪信息中的近似信号x(i-1,k)和细节信号d(i-1,k),重复到最优的分解尺度N上,得到小波离散分解的近似信号x(N,k)和细节信号d(N,k);然后将多尺度分析的结果:第(i+1-N)尺度上的细节信号和第N层尺度上的近似信号送到多尺度多传感器信息融合模块4等待处理;
步骤2:多尺度多传感器信息融合模块(4)将收集到的由小波域多尺度分析模块(2)输入的第i+1-N尺度上的细节信号和第N层尺度上的近似信号,以及数据统计性分析模块(3)输入的均值和标准差信号,进行多尺度多传感器信息融合,得到融合后的精确多尺度上的微机电陀螺仪信号,传递给小波多尺度重构模块(5);
所述的多尺度多传感器信息融合,具体按照以下步骤实施:
假设,传感器的测量模型为:
Xi=X+Vi,i=1,2,…,m,
其中,X表示要估计的传感器的真值,Vi为均值为0,标准差为有限值的测量噪声,近似认为传感器测量数据服从高斯分布,m个传感器的标准差分别为各传感器的加权因子分别为归一化权值W1,W2,…,Wi,各个传感器的标准差满足:
融合后的X值和加权因子满足:
则由X1,X2,…,Xi彼此独立,时域融合后序列的标准差:
根据多元函数求值理论,求出总均方误差最小时所对应的加权因子为:
时域加权融合的最小均方误差为:
根据单一尺度上的最优加权信息融合公式,小波域随机信号的重构公式为:
XJ是经过J层分解后重构回去的融合数据,是第J层的尺度函数,ψj,k是从1到J层的小波函数,hn(J,k)是第J层上k个近似系数,gn(j,k)是1到J层的k个细节系数,经过J层分解后重构回去的融合结果:
式中,各传感器相应近似信号和细节信号的加权因子分别为Whi和Wgi,同理有权值的归一化条件,hn,gn是独立不相关的,多尺度数据融合重构序列的方差:
关于各加权因子的多元二次函数,结合加权因子的约束条件为:
根据拉格朗日方程组,求取多元函数的极值可得总均方误差最小时所对应的加权因子为:
此时所对应的最小均方误差为:
步骤3:小波多尺度重构模块(5)对步骤2得到的融合后的精确多尺度上的微机电陀螺仪信号进行小波多尺度重构,得到原始尺度上的微机电陀螺信号,传递给组合微机电惯性测量设备模块(6);
所述的小波多尺度重构,具体按照以下步骤实施:
式中,x(i,l)和d(i,l)分别是尺度i上的近似信号和细节信号,h(2l-k)和g(2l-k)分别是尺度i上的尺度系数和小波系数,x(i+1,k)是经过小波重构得到的尺度i+1上的近似信号,即尺度i+1上陀螺仪的近似信息;
步骤4:组合微机电惯性测量设备模块(6)将步骤3得到的原始尺度上的微机电陀螺信号整理成需要的形式输出,完成基于MEMS的小波域多传感器信息融合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105228251A CN102023010B (zh) | 2010-10-26 | 2010-10-26 | 基于mems的小波域多传感器信息融合系统及融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105228251A CN102023010B (zh) | 2010-10-26 | 2010-10-26 | 基于mems的小波域多传感器信息融合系统及融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102023010A CN102023010A (zh) | 2011-04-20 |
CN102023010B true CN102023010B (zh) | 2012-06-13 |
Family
ID=43864591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010105228251A Expired - Fee Related CN102023010B (zh) | 2010-10-26 | 2010-10-26 | 基于mems的小波域多传感器信息融合系统及融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102023010B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103162710B (zh) * | 2011-12-15 | 2015-11-11 | 洛阳理工学院 | 基于小波熵的mems陀螺故障检测系统及其检测方法 |
CN103363966B (zh) * | 2012-03-26 | 2016-01-20 | 北京星网宇达科技股份有限公司 | 一种低成本组合型陀螺仪 |
CN102654403A (zh) * | 2012-05-07 | 2012-09-05 | 扬州蓝剑电子系统工程有限公司 | 一种基于小波熵的微机电陀螺信息融合系统及方法 |
CN102706347B (zh) * | 2012-05-17 | 2015-03-25 | 南京航空航天大学 | 一种惯性传感器网络节点装置及其信息融合方法 |
CN108132053B (zh) * | 2017-11-24 | 2020-01-07 | 北京工商大学 | 一种行人轨迹构建方法、系统及惯性测量装置 |
CN108108766B (zh) * | 2017-12-28 | 2021-10-29 | 东南大学 | 基于多传感器数据融合的驾驶行为识别方法及系统 |
CN108716913A (zh) * | 2018-07-03 | 2018-10-30 | 深圳市中科金朗产业研究院有限公司 | 一种角速度测量装置及运动控制装置 |
CN110362779A (zh) * | 2019-06-11 | 2019-10-22 | 南京江岛环境科技研究院有限公司 | 一种多尺度环境数据融合方法 |
CN113434713B (zh) * | 2021-06-18 | 2024-03-12 | 广东翁源滃江源国家湿地公园管理处 | 一种基于一张图的林业大数据建设方法及系统 |
CN113807318B (zh) * | 2021-10-11 | 2023-10-31 | 南京信息工程大学 | 一种基于双流卷积神经网络和双向gru的动作识别方法 |
CN114018279B (zh) * | 2021-10-25 | 2023-08-22 | 北京航空航天大学 | 一种阵列式传感器多采样率数据融合姿态校正方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201626A (zh) * | 2007-12-10 | 2008-06-18 | 华中科技大学 | 机器人自主定位系统 |
-
2010
- 2010-10-26 CN CN2010105228251A patent/CN102023010B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201626A (zh) * | 2007-12-10 | 2008-06-18 | 华中科技大学 | 机器人自主定位系统 |
Non-Patent Citations (3)
Title |
---|
任亚飞,柯熙政,李树州.基于小波熵的组合定位系统数据融合.《仪器仪表学报》.2006,第27卷(第6期),1323-1325. * |
任亚飞,柯熙政.基于阿伦方差的微机电陀螺误差建模及其粒子滤波.《中国计量学院学报》.2009,第20卷(第2期),102-106. * |
柯熙政,任亚飞.多尺度多传感器融合算法在微机电陀螺数据处理中的应用.《兵工学报》.2009,第30卷(第7期),994-998. * |
Also Published As
Publication number | Publication date |
---|---|
CN102023010A (zh) | 2011-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102023010B (zh) | 基于mems的小波域多传感器信息融合系统及融合方法 | |
Park et al. | Displacement estimation using multimetric data fusion | |
Weng et al. | Output-only modal identification of a cable-stayed bridge using wireless monitoring systems | |
CN102279084B (zh) | 基于微惯性测量组合的输电导线舞动定位方法 | |
CN104407328B (zh) | 基于空间脉冲响应匹配的封闭空间声源定位方法及系统 | |
CN104123463A (zh) | 一种随机动态载荷的时域识别方法 | |
CN104133950B (zh) | 一种悬臂梁运行模态分析实验方法及装置 | |
CN101957443B (zh) | 声源定位方法 | |
CN101105126B (zh) | 基于微石英角速率传感器的随钻方位测量误差补偿方法 | |
CN102306206B (zh) | 自适应一致性数据融合方法 | |
CN104112072A (zh) | 基于小波阈值去噪的主成分分析的工作模态参数识别方法 | |
CN104165742B (zh) | 一种基于互谱函数的运行模态分析实验方法及装置 | |
CN101587007A (zh) | 识别柔性桥梁结构动力参数的惟输出小波基分析方法 | |
Kaloop et al. | Analysis of the dynamic behavior of structures using the high-rate GNSS-PPP method combined with a wavelet-neural model: Numerical simulation and experimental tests | |
CN103714045A (zh) | 面向异步多速率不均匀采样观测数据的信息融合估计方法 | |
CN104931040B (zh) | 基于机器学习的北斗ⅱ代导航系统电力铁塔变形监测设备安装和调试方法 | |
CN104921736A (zh) | 一种包含参数估计功能滤波模块的连续血糖监测设备 | |
CN104915534A (zh) | 基于序列学习的电力铁塔变形分析与决策方法 | |
CN106548031A (zh) | 一种结构模态参数识别方法 | |
CN101901209A (zh) | 基于改进emd和arma模型的结构响应分析方法 | |
CN106202977A (zh) | 一种基于盲源分离算法的低频振荡模式分析方法 | |
CN110398775B (zh) | 隧道突涌水灾害微震事件信号波动初至拾取方法及系统 | |
CN105066967A (zh) | 基于mems运动传感器的测量波浪方法 | |
Zou et al. | Virtual sensing of subsoil strain response in monopile-based offshore wind turbines via Gaussian process latent force models | |
CN103018729A (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: 20120613 Termination date: 20211026 |
|
CF01 | Termination of patent right due to non-payment of annual fee |