CN101413824B - 一种基于随机传声器阵列的运动物体声场测量方法 - Google Patents
一种基于随机传声器阵列的运动物体声场测量方法 Download PDFInfo
- Publication number
- CN101413824B CN101413824B CN2008102388872A CN200810238887A CN101413824B CN 101413824 B CN101413824 B CN 101413824B CN 2008102388872 A CN2008102388872 A CN 2008102388872A CN 200810238887 A CN200810238887 A CN 200810238887A CN 101413824 B CN101413824 B CN 101413824B
- Authority
- CN
- China
- Prior art keywords
- mentioned
- array
- microphone
- random
- microphone array
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000013598 vector Substances 0.000 claims abstract description 61
- 238000005259 measurement Methods 0.000 claims abstract description 22
- 238000012216 screening Methods 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000010586 diagram Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000001629 suppression Effects 0.000 claims description 3
- 238000003491 array Methods 0.000 abstract description 7
- 238000013461 design Methods 0.000 abstract description 4
- 238000001914 filtration Methods 0.000 abstract description 2
- 238000012795 verification Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000001093 holography Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种基于随机传声器阵列的运动物体声场测量方法,属于噪声分析和控制技术领域。通过基于阵列单元的随机传声器阵列生成、根据基线向量统计筛选、仿真计算评价三个步骤设计二维随机传声器阵列布局,按照阵列布局搭建随机传声器阵列,对运动物体声场进行测量。本发明方法减少了传声器布局的任意性,降低布局设计中的搜索空间;应用基线统计指标预先滤除不满足要求的阵列布局,减少了性能验证的计算量,提高了搜索效率;通过仿真计算对不同频率下的随机传声器阵列性能进行验证,保证所得到的阵列满足测量使用要求。本发明方法使运动物体声场测量系统能够在较短时间内确定随机传声器阵列布局,实现对运动声源更准确的定位和分析。
Description
技术领域
本发明涉及一种基于随机传声器阵列的运动物体声场测量方法,属于噪声分析和控制技术领域。
背景技术
汽车、火车等复杂运动机械的噪声对环境影响很大。准确的识别和定位主要的发声部件所在,能够为噪声治理工作提供依据,对噪声测量以及治理工作都具有重要意义。由于运动的原因,对这类(多声源)运动型噪声,要比较精确地获取其声场,将各个声源分辨出来,并准确定位是比较困难的。
目前,世界上针对交通工具类的运动声源,对其声场进行测量分析研究的主要有声全息方法和阵列方法。全息法测量运动声场受多普勒效应影响很大,其工程实际应用受到了局限。与全息技术相比,阵列技术在运动声源识别研究中的应用更为广泛。其中随机阵列的布局方式通过传声器的不规则布置可以有效的克服规则阵列的旁瓣效应,提高声源识别的精确度;目前对于二维随机阵列布局的优化方法一般效率较低,需要进行大量计算才能找到性能较好的布局。
发明内容
本发明的目的是提出一种基于随机传声器阵列的运动物体声场测量方法,通过基于阵列单元的随机传声器阵列生成、根据基线向量统计筛选、仿真计算评价三个步骤进行二维随机传声器阵列布局的设计,按照计算得出的阵列布局搭建随机传声器阵列,对运动物体声场进行测量,从而提高声场测量的精确度。
本发明提出的基于随机传声器阵列的运动物体声场测量方法,包括以下步骤:
(1)设定测量用传声器之间的最短间隔距离为d,
其中c为声速,fmax为设定的最大测量分析频率;
(2)设上述测量用传声器所在的阵列平面A的长度为L1、高度为H1,待测运动物体所在平面R的长度为L2、高度为H2,阵列平面A与平面R之间的垂直距离为D,阵列平面A的长度L1和高度H1是上述最短间隔距离d的整数倍;
(3)将上述阵列平面A划分为n个大小形状相同的阵列单元,每个阵列单元的长度和宽度分别为k1×d、k2×d,其中k1、k2为正整数,d为上述最短间隔距离,将每个阵列单元划分为k1×k2个方格,每个方格的中心点为测量用传声器的预设置点,共有k1×k2预设置点;
(4)在上述k1×k2预设置点中,随机选定一个预设置点作为传声器放置位置,得到一个由n个传声器构成的随机传声器阵列,重复该过程N1次,生成N1个随机传声器阵列;
(5)计算上述N1个随机传声器阵列中每个随机传声器阵列的所有基线向量,将同一随机传声器阵列的所有基线向量放到同一点距坐标系下,得到基线向量的点距坐标图,设定点距坐标图中坐标轴的单位为最短间隔距离d,形成U个由直线x=±J和y=±J所构成的正方形坐标环,其中J=1,2,…U-1,U,U=max(片1/d,L1/d),计算基线向量终点在第1坐标环上的所有基线向量数M1,基线向量终点在第1坐标环上且互相不重合的所有基线向量总数K1,基线向量终点在第2坐标环上且互相不重合的所有基线向量总数K2,以及基线向量终点在第U坐标环上且互相不重合的所有基线向量总数KU;
(6)根据上述基线向量总数,对上述N1个随机传声器阵列进行筛选,其过程为:
(a)若上述阵列单元数n≤24,且上述基线向量数M1≤4,则进行步骤(b),若上述阵列单元数n>24,且上述基线向量总数K1≤4,则进行步骤(b);
(b)若上述阵列单元数n≥14,且上述基线向量数K2≥14,则进行步骤(c),若上述阵列单元数10≤n≤13,且上述基线向量数K2≥12,则进行步骤(c),若上述阵列单元数8≤n≤9,且上述基线向量数K2≥10,则进行步骤(c);
(c)若上述基线向量数KU≥4,则选择该随机传声器阵列;
若γf≥n且上述阵列单元数8≤n<16,或者γf≥16且n≥16,则该随机传声器阵列满足频率f下的筛选条件,进行下一频率的筛选。
在下一频率重复上述仿真步骤,若在对每一频率都满足上述筛选条件,则得到测量用的随机传声器阵列,否则对上述选择的N2个随机传声器阵列中的下一个随机传声器阵列重复上述筛选过程,直至得到一个满足上述筛选条件的随机传声器阵列。
(9)根据上述筛选得到的随机传声器阵列组成测量系统,采集运动物体发出的声压信号,计算得到待测运动物体所在平面R上任意一点s(ε,η)处在测量时间t1至t2内的声场特征函数:
其中,pi(t)为t时刻第i个传声器接收到的信号声压,ri(t,ε,η)为t时刻,任意点s(ε,η)与第i个传声器之间的距离。
上述方法中,阵列单元为2×2、2×3、2×4、3×3、3×4或2×6中的任何一种。
本发明提出的基于随机传声器阵列的运动物体声场测量方法,通过基于阵列单元的随机传声器阵列生成、根据基线向量统计筛选、仿真计算评价三个步骤进行二维随机传声器阵列布局的设计,按照计算得出的阵列布局搭建随机传声器阵列,对运动物体声场进行测量。通过划分阵列单元的方法进行随机传声器阵列生成,减少布局的任意性,降低布局设计中的搜索空间;应用基线统计指标预先滤除不满足要求的阵列布局,从而减少了性能验证的计算量,提高了搜索效率;通过仿真计算对不同频率下的随机传声器阵列性能进行验证,保证所得到的阵列满足测量使用要求。通过应用本发明方法,使运动物体声场测量系统能够在较短时间内确定随机传声器阵列布局,实现对运动声源更准确的定位和分析。
附图说明
图1是应用本发明方法进行运动物体声场测量的现场示意图。
图2是被测运动物体所在平面与阵列平面的空间关系示意图。
图3是一个由16个2×4式阵列单元组成的随机传声器阵列示例图。
图4是一个2×4式阵列单元结构示意图。
图5是一个随机传声器阵列的基线向量点距坐标图。
图6是仿真计算得出的声场特征函数分布示意图。
图7是不同形式的阵列单元结构示意图。
图1中,1是被测运动物体,2是激光接收器,3是传声器,4是传声器阵列支架,5是PC机,6是激光发射器,7是前置放大器接口箱,8是数字采集箱。
图7中,(a)是2×2式阵列单元,(b)是2×3式阵列单元,(c)是2×4式阵列单元,(d)是3×3式阵列单元,(e)是3×4式阵列单元,(f)是2×6式阵列单元。
具体实施方式
以下结合附图,详细介绍本发明提出的基于随机传声器阵列的运动物体声场测量方法,应用本发明方法进行运动物体声场测量的现场如附图1所示。
其中c为声速,fmax为设定的最大测量分析频率;
(2)如附图2所示,设上述测量用传声器所在的阵列平面A的长度为L1、高度为H1,待测运动物体所在平面R的长度为L2、高度为H2,阵列平面A与平面R之间的垂直距离为D,阵列平面A的长度L1和高度H1是上述最短间隔距离d的整数倍;两面垂直距离D取值范围为2~6米。
(3)将上述阵列平面A划分为n个大小形状相同的阵列单元,每个阵列单元的长度和宽度分别为k1×d、k2×d,其中k1、k2为正整数,d为上述最短间隔距离,将每个阵列单元划分为k1×k2个方格,每个方格的中心点为测量用传声器的预设置点,共有k1×k2预设置点;阵列单元可以按不同方向放置,划分得到的n个阵列单元应覆盖整个阵列平面的面积,或者除边角少量位置外的全部面积。如图3所示为一个16阵列单元的划分方案实例,图4所示为一个有2×4个传声器预设置点的阵列单元结构。
(4)在上述k1×k2预设置点中,随机选定一个预设置点作为传声器放置位置,得到一个由n个传声器构成的随机传声器阵列,如图3为一个16阵列单元的阵列布局示例,其中“”表示随机选定的传声器位置;重复该过程N1次,生成N1个随机传声器阵列;
(5)计算上述N1个随机传声器阵列中每个随机传声器阵列的所有基线向量,如图5所示,将同一随机传声器阵列的所有基线向量放到同一点距坐标系下,得到基线向量的点距坐标图,其中基线向量是从一个传声器位置指向另外一个传声器位置的有向线段,对于传声器数目为n的传声器阵列,其基线向量数目共有n(n-1)个。设定点距坐标图中坐标轴的单位为最短间隔距离d,形成U个由直线x=±J和y=±J所构成的正方形坐标环,其中J=1,2,…U-1,U,U=max(H1/d,L1/d);各坐标环均用数字表示,由x=±J及y=±J组成的坐标环称为第J环。计算基线向量终点在第1坐标环上的所有基线向量数M1,基线向量终点在第1坐标环上且互相不重合的所有基线向量总数K1,基线向量终点在第2坐标环上且互相不重合的所有基线向量总数K2,以及基线向量终点在第U坐标环上且互相不重合的所有基线向量总数KU;
(6)根据上述基线向量总数,对上述N1个随机传声器阵列进行筛选,其过程为:
(a)若上述阵列单元数n≤24,且上述基线向量数M1≤4,则进行步骤(b),若上述阵列单元数n>24,且上述基线向量总数K1≤4,则进行步骤(b);
(b)若上述阵列单元数n≥14,且上述基线向量数K2≥14,则进行步骤(c),若上述阵列单元数10≤n≤13,且上述基线向量数K2≥12,则进行步骤(c),若上述阵列单元数8≤n≤9,且上述基线向量数K2≥10,则进行步骤(c);
(c)若上述基线向量数KU≥4,则选择该随机传声器阵列;
(8)设上述选择的随机传声器阵列数有N2个,分别对N2个随机传声器阵列进行仿真计算,仿真计算模型如图2所示,仿真中的参数均取实际测量中的数值,仿真被测声源Q为一个点声源,位于待测运动物体所在平面R的中心,设定其声源特征函数其中为仿真实测声源强度,f为仿真被测声源的频率。
其中为在t时刻Q点距第i个传声器间的距离,θi(t)为在t时刻Q点到第i个传声器连线与Q点运动方向间的夹角,M为马赫数。
若γf≥n且上述阵列单元数8≤n<16,或者γf≥16且n≥16,则该随机传声器阵列满足频率f下的筛选条件,进行下一频率的筛选。
在下一频率重复上述仿真步骤,若在对每一频率都满足上述筛选条件,则得到测量用的随机传声器阵列,否则对上述选择的N2个随机传声器阵列中的下一个随机传声器阵列重复上述筛选过程,直至得到一个满足上述筛选条件的随机传声器阵列。
其中上述筛选步骤中仿真频率的选择方法为:将考察频率范围分为高频段、低频段和中频段,在每个频段中,各取1~2个特征频率。
(9)根据上述筛选得到的随机传声器阵列组成如图1所示的测量系统,采集运动物体发出的声压信号,计算得到待测运动物体所在平面R上任意一点s(ε,η)处在测量时间t1至t2内的声场特征函数:
其中,pi(t)为t时刻第i个传声器接收到的信号声压,ri(t,ε,η)为t时刻,任意点s(ε,η)与第i个传声器之间的距离。
如图7中的(a)、(b)、(c)、(d)、(e)、(f)、(g)所示,上述步骤(3)中所述的阵列单元结构可采用2×2、2×3、2×4、3×3、3×4或2×6中的任何一种。
Claims (2)
1.一种基于随机传声器阵列的运动物体声场测量方法,其特征在于该方法包括以下步骤:
(1)设定测量用传声器之间的最短间隔距离为d,
其中c为声速,fmax为设定的最大测量分析频率;
(2)设上述测量用传声器所在的阵列平面A的长度为L1、高度为H1,待测运动物体所在平面R的长度为L2、高度为H2,阵列平面A与平面R之间的垂直距离为D,阵列平面A的长度L1和高度H1是上述最短间隔距离d的整数倍;
(3)将上述阵列平面A划分为n个大小形状相同的阵列单元,每个阵列单元的长度和宽度分别为k1×d、k2×d,其中k1、k2为正整数,d为上述最短间隔距离,将每个阵列单元划分为k1×k2个方格,每个方格的中心点为测量用传声器的预设置点,共有k1×k2预设置点;
(4)在上述k1×k2预设置点中,随机选定一个预设置点作为传声器放置位置,得到一个由n个传声器构成的随机传声器阵列,重复该过程N1次,生成N1个随机传声器阵列;
(5)计算上述N1个随机传声器阵列中每个随机传声器阵列的所有基线向量,将同一随机传声器阵列的所有基线向量放到同一点距坐标系下,得到基线向量的点距坐标图,所述的基线向量是从一个传声器位置指向另外一个传声器位置的有向线段,设定点距坐标图中坐标轴的单位为最短间隔距离d,形成U个由直线x=±J和y=±J所构成的正方形坐标环,其中J=1,2,…U-1,U,U=max(H1/d,L1/d),计算基线向量终点在第1坐标环上的所有基线向量数M1,基线向量终点在第1坐标环上且互相不重合的所有基线向量总数K1,基线向量终点在第2坐标环上且互相不重合的所有基线向量总数K2,以及基线向量终点在第U坐标环上且互相不重合的所有基线向量总数KU;
(6)根据上述基线向量总数,对上述N1个随机传声器阵列进行筛选,其过程为:
(a)若上述阵列单元数n≤24,且上述基线向量数M1≤4,则进行步骤(b),若上述阵列单元数n>24,且上述基线向量总数K1≤4,则进行步骤(b);
(b)若上述阵列单元数n≥14,且上述基线向量数K2≥14,则进行步骤(c),若上述阵列单元数10≤n≤13,且上述基线向量数K2≥12,则进行步骤(c),若上述阵列单元数8≤n≤9,且上述基线向量数K2≥10,则进行步骤(c);
(c)若上述基线向量数KU≥4,则选择该随机传声器阵列;
若γf≥n且上述阵列单元数8≤n<16,或者γf≥16且n≥16,则该随机传声器阵列满足频率f下的筛选条件,进行下一频率的筛选;
在下一频率重复上述仿真步骤,若在对每一频率都满足上述筛选条件,则得到测量用的随机传声器阵列,否则对上述选择的N2个随机传声器阵列中的下一个随机传声器阵列重复上述筛选过程,直至得到一个满足上述筛选条件的随机传声器阵列;
(8)根据上述筛选得到的随机传声器阵列组成测量系统,采集运动物体发出的声压信号,计算得到待测运动物体所在平面R上任意一点s(ε,η)处在测量时间t1至t2内的声场特征函数:
其中,pi(t)为t时刻第i个传声器接收到的信号声压,ri(t,ε,η)为t时刻,任意点s(ε,η)与第i个传声器之间的距离。
2.如权利要求1所述的测量方法,其特征在于其中所述的阵列单元为2×2、2×3、2×4、3×3、3×4或2×6中的任何一种。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102388872A CN101413824B (zh) | 2008-12-04 | 2008-12-04 | 一种基于随机传声器阵列的运动物体声场测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102388872A CN101413824B (zh) | 2008-12-04 | 2008-12-04 | 一种基于随机传声器阵列的运动物体声场测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101413824A CN101413824A (zh) | 2009-04-22 |
CN101413824B true CN101413824B (zh) | 2010-04-07 |
Family
ID=40594465
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102388872A Expired - Fee Related CN101413824B (zh) | 2008-12-04 | 2008-12-04 | 一种基于随机传声器阵列的运动物体声场测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101413824B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101662714B (zh) * | 2009-07-28 | 2012-08-15 | 南京大学 | 基于时间反转的复杂声场定位拾声传声器阵列设计方法 |
CN103308154B (zh) * | 2013-06-15 | 2016-05-04 | 西安费斯达自动化工程有限公司 | 一种多孔金属复合结构阵列驻极体噪声检测器 |
CN104678359B (zh) * | 2015-02-28 | 2017-01-04 | 清华大学 | 一种声场识别的多孔声全息方法 |
CN104964739B (zh) * | 2015-07-08 | 2017-12-19 | 安徽合力股份有限公司 | 一种便携式近场噪声源识别的声阵列支架 |
CN106500828A (zh) * | 2016-09-29 | 2017-03-15 | 中国人民解放军军械工程学院 | 一种手持式传声器阵列及其测试自行装备驾驶舱内噪声分布的方法 |
CN106525224B (zh) * | 2016-10-26 | 2019-12-17 | 青岛理工大学 | 一种测量多源运动声场的非规则声阵列搭建方法 |
WO2018119910A1 (zh) * | 2016-12-29 | 2018-07-05 | 清华大学 | 一种运动声源的定位方法和定位装置 |
CN107167227B (zh) * | 2017-05-31 | 2019-07-19 | 合肥工业大学 | 一种多普勒效应恢复方法 |
CN107218996B (zh) * | 2017-05-31 | 2019-07-19 | 合肥工业大学 | 一种多普勒效应消除方法 |
US12010483B2 (en) * | 2021-08-06 | 2024-06-11 | Qsc, Llc | Acoustic microphone arrays |
CN114199368B (zh) * | 2021-11-30 | 2024-04-26 | 北京工商大学 | 一种全频带pp声强自动测量装置和测量方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1293350A (zh) * | 2000-12-15 | 2001-05-02 | 清华大学 | 高速运动物体表面声场分析方法 |
CN101290347A (zh) * | 2008-06-13 | 2008-10-22 | 清华大学 | 用规则声阵列和单摄像机获取静止声源声场图像的方法 |
CN101295017A (zh) * | 2008-06-13 | 2008-10-29 | 清华大学 | 基于规则声阵列和双目视觉的获取运动声场视频的方法 |
-
2008
- 2008-12-04 CN CN2008102388872A patent/CN101413824B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1293350A (zh) * | 2000-12-15 | 2001-05-02 | 清华大学 | 高速运动物体表面声场分析方法 |
CN101290347A (zh) * | 2008-06-13 | 2008-10-22 | 清华大学 | 用规则声阵列和单摄像机获取静止声源声场图像的方法 |
CN101295017A (zh) * | 2008-06-13 | 2008-10-29 | 清华大学 | 基于规则声阵列和双目视觉的获取运动声场视频的方法 |
Non-Patent Citations (1)
Title |
---|
罗禹贡 等.应用动态声全息方法识别轿车的行使噪声源.汽车工程25 6.2003,25(6),595-598. * |
Also Published As
Publication number | Publication date |
---|---|
CN101413824A (zh) | 2009-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101413824B (zh) | 一种基于随机传声器阵列的运动物体声场测量方法 | |
CN106339556B (zh) | 一种轨道交通桥梁和钢轨全空间噪声预测方法 | |
CN100442030C (zh) | 一种声场分离方法 | |
CN104237849B (zh) | 一种双五元十字阵被动声定位融合方法 | |
CN101566496B (zh) | 采用双面振速测量和等效源法分离声场的方法 | |
CN101403634B (zh) | 压力速度法声场分离方法 | |
CN105205461A (zh) | 一种用于模态参数识别的信号降噪方法 | |
CN101414000B (zh) | 基于随机传声器阵列和双目视觉的获取运动声场视频的方法 | |
CN102955152B (zh) | 一种海浪的sar 信号模拟方法 | |
CN103235286B (zh) | 一种对电噪声源的高精度定位方法 | |
CN110399680B (zh) | 一种浅海弹性结构辐射声场计算方法 | |
CN104596636A (zh) | 声场分离方法 | |
CN112613092B (zh) | 一种路基压实度空间分布的预测方法和预测装置 | |
CN104678359B (zh) | 一种声场识别的多孔声全息方法 | |
CN104407328A (zh) | 基于空间脉冲响应匹配的封闭空间声源定位方法及系统 | |
CN104656123A (zh) | 一种测定区域岩体等效波速的方法 | |
CN106646645A (zh) | 一种新的重力正演加速方法 | |
CN105260568A (zh) | 基于离散型卡尔曼滤波的超高层建筑风荷载反分析方法 | |
CN102901950A (zh) | 平面阵列识别声源三维坐标的方法 | |
CN103616071A (zh) | Patch近场声全息-声品质客观参量三维分布可视化方法 | |
CN105466661A (zh) | 基于改进卡尔曼滤波的超高层建筑风荷载反分析方法 | |
CN104316958A (zh) | 一种识别不同尺度地层断裂的相干处理方法 | |
CN104568024A (zh) | 振动式流量计特征信号提取方法 | |
CN101922964A (zh) | 一种可辨别噪声来源方向的环境噪声监测装置和方法 | |
CN107884743A (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: 20100407 Termination date: 20151204 |
|
EXPY | Termination of patent right or utility model |