CN106887238B - 一种基于改进独立向量分析算法的声信号盲分离方法 - Google Patents

一种基于改进独立向量分析算法的声信号盲分离方法 Download PDF

Info

Publication number
CN106887238B
CN106887238B CN201710117865.XA CN201710117865A CN106887238B CN 106887238 B CN106887238 B CN 106887238B CN 201710117865 A CN201710117865 A CN 201710117865A CN 106887238 B CN106887238 B CN 106887238B
Authority
CN
China
Prior art keywords
separation
signal
frequency domain
signals
analysis algorithm
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.)
Active
Application number
CN201710117865.XA
Other languages
English (en)
Other versions
CN106887238A (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.)
University of Chinese Academy of Sciences
Original Assignee
Shanghai Institute of Microsystem and Information Technology of CAS
University of Chinese Academy of Sciences
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 Shanghai Institute of Microsystem and Information Technology of CAS, University of Chinese Academy of Sciences filed Critical Shanghai Institute of Microsystem and Information Technology of CAS
Priority to CN201710117865.XA priority Critical patent/CN106887238B/zh
Publication of CN106887238A publication Critical patent/CN106887238A/zh
Application granted granted Critical
Publication of CN106887238B publication Critical patent/CN106887238B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/27Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the analysis technique
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0272Voice signal separating

Landscapes

  • Engineering & Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Quality & Reliability (AREA)
  • Complex Calculations (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提供一种基于改进独立向量分析算法的声信号盲分离方法,包括:步骤1,通过M个传感器采集N个源信号的卷积混合信号,以得到M个观测信号;步骤2,通过短时傅里叶变换将所述观测信号从时域转换到频域,并建立基于分离矩阵的频域分离信号模型;步骤3,逐频点对频域的所述观测信号和分离矩阵进行预处理;步骤4,采用基于自适应更新的自然梯度独立向量分析算法逐频点进行盲分离处理,以得到频域的分离信号;步骤5,对频域的分离信号进行逆短时傅里叶变换,以得到时域的分离信号。本发明考虑声信号的卷积混合,自适应更新,能够更精确地恢复源信号,更快达到收敛;算法对更新步长的初始值不敏感,鲁棒性好,稳定性高。

Description

一种基于改进独立向量分析算法的声信号盲分离方法
技术领域
本发明涉及盲源分离信号处理领域,特别是涉及一种基于改进独立向量分析算法的声信号盲分离方法。
背景技术
随着数字信号处理理论和技术的发展以及相关学科的不断深入,大量有效的盲分离算法不断被提出,使盲分离问题逐渐成为当今信息处理领域中最热门的研究课题之一。在无线通信、图像处理、地震信号处理、阵列信号处理和生物医学信号处理等领域得到了广泛的应用。盲源分离用于解决瞬时线性混合问题目前已经有很成熟的算法,而且也得到了很好的应用,例如在图像增强和瞬时混合信号分离。
在实际系统中,传感器接收到的信号往往是源信号经过不同时延的线性组合,即传感器得到的观测信号是源信号的卷积和,称为线性卷积混合模型。这种混合模型更接近实际。所以语音信号的混合是典型的卷积混合信号,但是现有算法解决卷积混合信号分离问题并不是很有效,因为混合语音分离问题有其自身的难点,例如很多卷积盲分离算法是针对平稳信号的,而语音信号是典型的非平稳信号,因而关于实际情况中声信号的盲分离是一个亟待解决的难题。
发明内容
针对上述现有技术的不足,本发明提供一种基于改进独立向量分析算法的声信号盲分离方法。
为实现上述目的,本发明采用如下技术方案:
一种基于改进独立向量分析算法的声信号盲分离方法,包括以下步骤:
步骤1,通过M个传感器采集N个源信号的卷积混合信号,以得到M个观测信号;
步骤2,通过短时傅里叶变换将所述观测信号从时域转换到频域,并建立基于分离矩阵的频域分离信号模型;
步骤3,逐频点对频域的所述观测信号和分离矩阵进行预处理;
步骤4,采用基于自适应自然梯度的独立向量分析算法逐频点进行盲分离处理,以得到频域的分离信号;
步骤5,对频域的分离信号进行逆短时傅里叶变换,以得到时域的分离信号。
进一步地,在所述步骤1中,将N个源信号表示为s(t)=[s1(t),s2(t),...,sN(t)]T,M个观测信号表示为x(t)=[x1(t),x2(t),...,xM(t)]T,则在每一时刻t有如下关系式:
Figure BDA0001236134000000021
在式(1)中,aij(l)表示第i个源信号到第j个传感器的冲激响应,L表示冲激响应的长度,且L<∞。
进一步地,在所述步骤2中,通过短时傅里叶变换将式(1)变换到频域,得到Xf(n)=AfSf(n),并建立频域分离信号模型Yf(n)=Wf·Xf(n);
其中,f表示频点,且f=1,2,…,F,F表示频点数目,
Figure BDA0001236134000000022
Figure BDA0001236134000000023
Af和Wf分别表示第f个频点对应的混合矩阵和分离矩阵。
进一步地,在所述步骤3中,所述预处理包括使观测信号Xf(n)零均值化和白化,并使分离矩阵Wf正交。
进一步地,在所述步骤4中,所述自适应自然梯度的独立向量分析算法通过以下步骤实现:
步骤41,构建目标函数为:
Figure BDA0001236134000000024
在式(2)中,q(yi)为分离信号yi的概率密度函数,p(y1…yN)为联合概率密度函数,I为单位矩阵;
步骤42,采用自然梯度更新算法得到分离矩阵Wf,即,令第k+1次迭代公式表示为:
Wf(k+1)=Wf(k)+μ(k)[I-f(yf)(yf)T]Wf(k) (3),
在式(3)中,
Figure BDA0001236134000000031
μ(k)表示更新步长,且μ(k)的表达式如下:
μ(k)=α|1-eβ·RE| (4),
在式(4)中,α和β是经验常数,RE为目标函数的相对误差;
步骤43,根据Yf(n)=Wf·Xf(n)逐频点得到频域的分离信号Yf(n)。
优选地,在所述步骤42中,α取0.7,β取0.3。
进一步地,在所述步骤5中,对Yf(n)进行逆短时傅里叶变换,以得到时域的分离信号。
由于采用了上述的技术方案,本发明与现有技术相比,具有以下的优点和积极效果:
本发明利用改进的独立向量分析算法进行声信号的盲分离,相比传统的盲源分离算法,本发明考虑声信号的卷积混合,采用步长自适应的频域盲分离算法,能够更精确地恢复源信号,更快达到收敛条件;而且本发明对更新步长的初始值不敏感,自适应更新,具有很好的鲁棒性,稳定性。本发明在无线通信,雷达信号处理以及声信号处理等领域具有很好的应用前景。
附图说明
图1为本发明一种基于改进独立向量分析算法的声信号盲分离方法的流程图;
图2为本发明一种基于改进独立向量分析算法的声信号盲分离方法的处理框架图;
图3为本发明中采用的自适应更新自然梯度独立向量分析算法的流程图;
图4为本发明一个实例中源信号的波形图;
图5为本发明一个实例中混合信号的波形图;
图6为本发明一个实例中得到的分离信号的波形图。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
在现有技术中,对于声音信号的卷积混合基本处理思想为:将观测信号通过STFT变换到频域,然后使用复数域上的瞬时盲源分离算法对每个频段依次进行分离,最后再通过短时傅里叶逆变换ISTFT将分离所得的信号变换回时域。
基于上述内容,本发明提供一种基于改进的独立向量分析(IVA)算法的声信号盲分离方法,如图1所示,本发明的方法包括以下步骤:步骤1,通过M个传感器采集N个独立的源信号的卷积混合信号,得到M个观测信号;步骤2,通过短时傅里叶变换(STFT)将观测信号从时域转换到频域,并建立基于分离矩阵的频域分离信号模型;步骤3,逐频点对频域的观测信号和分离矩阵进行预处理;步骤4,采用基于自然梯度的独立向量分析算法逐频点进行盲分离处理,以得到频域的分离信号;步骤5,对频域的分离信号进行逆短时傅里叶变换,以得到时域的分离信号。
下面结合图2和图3,对本发明各步骤进行详细描述:
在步骤1中,N个源信号表示为s(t)=[s1(t),s2(t),...,sN(t)]T,M个观测信号表示为x(t)=[x1(t),x2(t),...,xN(t)]T,其中,si(t)表示t时刻的第i(i=1,2,…,N)个源信号,xj(t)表示t时刻的第j(j=1,2,…,M)个观测信号。周知的是,观测信号是由N个源信号卷积混合而成,则在每一时刻t都有如下关系式:
Figure BDA0001236134000000041
式(1)写成向量形式为X(t)=AS(t),aij(t)为混合矩阵A的第i行、第j列的元素。
在式(1)中,aij(l)表示第i个源信号到第j个传感器的冲激响应,L表示冲激响应的长度,理论上L<∞。在一个实施例中,一源信号的波形见图4,一观测信号的波形图见图5。
在步骤2中,本发明采用STFT将观测信号从时域变换到频域,转换为复数瞬时混合问题,建立频域分离信号模型。具体来说,在频域将时域时间t用对应的n表示,n对应于每个短时傅里叶变换,得到频域观测信号Xf(n)=AfSf(n),则可建立对应的频域分离信号模型Yf(n)=Wf·Xf(n)。
其中,f表示频点,f=1,2,…,F,F为频点数目,
Figure BDA0001236134000000051
Figure BDA0001236134000000052
Af和Wf分别表示第f个频点对应的混合矩阵和分离矩阵。
在步骤3中,预处理包括使观测数据Xf(n)零均值化和白化,并使分离矩阵Wf正交,从而保持Yf(n)均值化和白化。在本发明中,白化处理采用PCA降维方法,目的是去除信号之间的相关性,即保证分离矩阵Wf的各行之间保持正交,以减少复杂度和更好的收敛性。
在步骤4中,本发明采用多元变量目标函数,选择两个函数之间的KL散度(相对熵)来度量信号的独立性。在第f个频点内,构建目标函数为:
Figure BDA0001236134000000053
在式(2)中,q(yi)为每个分离信号yi的概率密度函数,p(y1…yN)为联合概率密度。
本发明采用了自然梯度更新算法来更新分离矩阵Wf,令第k+1次迭代公式表示为:
Wf(k+1)=Wf(k)+μ(k)[I-f(yf)(yf)T]Wf(k) (3),
其中,I为单位矩阵;认为分离信号服从多元变量超高斯分布,
Figure BDA0001236134000000054
特别地,本发明采用了如下步长自适应更新公式:
μ(k)=α|1-eβ·RE| (4),
其中α和β是经验常数,一般α取0.7,β取0.3,收敛条件是目标函数的相对误差小于10-6。RE为目标函数的相对误差,第k次迭代时对应的RE即为目标函数在第k次迭代时的瞬时估计值与在第k-1次迭代时的瞬时估计值之差。利用本发明的自适应算法迭代约25次就能达到收敛,比未改进的IVA算法快了很多倍。
然后,根据Yf(n)=Wf·Xf(n)逐频点对频域观测信号进行源信号恢复,得到频域的分离信号Yf(n)。
在步骤5中,对频域的分离信号Yf(n)进行逆短时傅里叶变换,即可得到时域的分离信号。输出的分离信号语音可懂度高,其波形图见图5。
不难发现,本发明采用改进独立向量分析算法,充分利用了源信号在频段间的相关信息,同时采用自适应步长更新的自然梯度算法。相比现有的盲源分离算法,能够有效分离实际卷积混合的声音信号,并且算法结构简单,收敛更快,鲁棒性好,稳定性高。

Claims (6)

1.一种基于改进独立向量分析算法的声信号盲分离方法,其特征在于,包括以下步骤:
步骤1,通过M个传感器采集N个源信号的卷积混合信号,
以得到M个观测信号;
步骤2,通过短时傅里叶变换将所述观测信号从时域转换到频域,并建立基于分离矩阵的频域分离信号模型;
步骤3,逐频点对频域的所述观测信号和分离矩阵进行预处理;
步骤4,采用基于自适应更新的自然梯度独立向量分析算法逐频点进行盲分离处理,以得到频域的分离信号;
步骤5,对频域的分离信号进行逆短时傅里叶变换,以得到时域的分离信号;
其中,步骤4中所述自适应更新的自然梯度的独立向量分析算法通过以下步骤实现:
步骤41,构建目标函数为:
Figure FDA0002361989380000011
在式(2)中,q(yi)为分离信号yi的概率密度函数,p(y1...yN)为联合概率密度函数,I为单位矩阵,f表示频点,yf表示频域中第f个频点的分离信号;
步骤42,采用自然梯度更新算法得到分离矩阵Wf,即,令第k+1次迭代公式表示为:
Wf(k+1)=Wf(k)+μ(k)[I-f(yf)(yf)T]Wf(k) (3),
在式(3)中,
Figure FDA0002361989380000012
μ(k)表示更新步长,且μ(k)的表达式如下:
μ(k)=α|1-eβ·RE| (4),
在式(4)中,α和β是经验常数,RE为目标函数的相对误差;
步骤43,根据Yf(n)=Wf·Xf(n)逐频点得到频域的分离信号Yf(n)。
2.根据权利要求1所述的基于改进独立向量分析算法的声信号盲分离方法,其特征在于,在所述步骤1中,将N个源信号表示为s(t)=[s1(t),s2(t),...,sN(t)]T,M个观测信号表示为x(t)=[x1(t),x2(t),...,xM(t)]T,则在每一时刻t有如下关系式:
Figure FDA0002361989380000021
在式(1)中,aij(l)表示第i个源信号到第j个传感器的冲激响应,L表示冲激响应的长度,且L<∞。
3.根据权利要求2所述的基于改进独立向量分析算法的声信号盲分离方法,其特征在于,在所述步骤2中,通过短时傅里叶变换将式(1)变换到频域,得到Xf(n)=AfSf(n),并建立频域分离信号模型Yf(n)=Wf·Xf(n);
其中,f表示频点,且f=1,2,...,F,F表示频点数目,
Figure FDA0002361989380000022
Figure FDA0002361989380000023
Af和Wf分别表示第f个频点对应的混合矩阵和分离矩阵。
4.根据权利要求3所述的基于改进独立向量分析算法的声信号盲分离方法,其特征在于,在所述步骤3中,所述预处理包括使观测信号Xf(n)零均值化和白化,并使分离矩阵Wf正交。
5.根据权利要求1所述的基于改进独立向量分析算法的声信号盲分离方法,其特征在于,在所述步骤42中,α取0.7,β取0.3。
6.根据权利要求1所述的基于改进独立向量分析算法的声信号盲分离方法,其特征在于,在所述步骤5中,对Yf(n)进行逆短时傅里叶变换,以得到时域的分离信号。
CN201710117865.XA 2017-03-01 2017-03-01 一种基于改进独立向量分析算法的声信号盲分离方法 Active CN106887238B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710117865.XA CN106887238B (zh) 2017-03-01 2017-03-01 一种基于改进独立向量分析算法的声信号盲分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710117865.XA CN106887238B (zh) 2017-03-01 2017-03-01 一种基于改进独立向量分析算法的声信号盲分离方法

Publications (2)

Publication Number Publication Date
CN106887238A CN106887238A (zh) 2017-06-23
CN106887238B true CN106887238B (zh) 2020-05-15

Family

ID=59179008

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710117865.XA Active CN106887238B (zh) 2017-03-01 2017-03-01 一种基于改进独立向量分析算法的声信号盲分离方法

Country Status (1)

Country Link
CN (1) CN106887238B (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019016494A1 (en) * 2017-07-19 2019-01-24 Cedar Audio Ltd ACOUSTIC SOURCE SEPARATION SYSTEMS
CN107563300A (zh) * 2017-08-08 2018-01-09 浙江上风高科专风实业有限公司 基于预白化方法的降噪预处理技术
CN108198570B (zh) * 2018-02-02 2020-10-23 北京云知声信息技术有限公司 审讯时语音分离的方法及装置
CN108597531B (zh) * 2018-03-28 2021-05-28 南京大学 一种通过多声源活动检测来改进双通道盲信号分离的方法
CN109086686B (zh) * 2018-07-12 2022-09-30 西安电子科技大学 基于自适应动量因子的时变信道下的盲源分离方法
CN110890098B (zh) * 2018-09-07 2022-05-10 南京地平线机器人技术有限公司 盲信号分离方法、装置和电子设备
CN109243483B (zh) * 2018-10-17 2022-03-08 西安交通大学 一种含噪频域卷积盲源分离方法
CN109671447A (zh) * 2018-11-28 2019-04-23 广东工业大学 一种双通道欠定卷积混叠信号盲分离方法
CN111986695B (zh) * 2019-05-24 2023-07-25 中国科学院声学研究所 一种无重叠子带划分快速独立向量分析语音盲分离方法及系统
CN110473565A (zh) * 2019-07-04 2019-11-19 中国人民解放军63892部队 一种无需源数估计的独立向量分析信号分离方法
CN113223553B (zh) * 2020-02-05 2023-01-17 北京小米移动软件有限公司 分离语音信号的方法、装置及介质
CN111312276B (zh) * 2020-02-14 2023-01-17 北京声智科技有限公司 一种音频信号处理的方法、装置、设备和介质
CN111667846B (zh) * 2020-05-14 2023-06-16 北京声智科技有限公司 一种盲源分离方法及装置
CN111916075A (zh) * 2020-07-03 2020-11-10 北京声智科技有限公司 一种音频信号的处理方法、装置、设备及介质
CN112349292B (zh) * 2020-11-02 2024-04-19 深圳地平线机器人科技有限公司 信号分离方法和装置、计算机可读存储介质、电子设备
CN113470675A (zh) * 2021-06-30 2021-10-01 北京小米移动软件有限公司 音频信号处理方法及装置
CN116866123B (zh) * 2023-07-13 2024-04-30 中国人民解放军战略支援部队航天工程大学 一种无正交受限的卷积盲分离方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102231280A (zh) * 2011-05-06 2011-11-02 山东大学 卷积语音信号的频域盲分离排序算法
CN102592607A (zh) * 2012-03-30 2012-07-18 北京交通大学 一种使用盲语音分离的语音转换系统和方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8874439B2 (en) * 2006-03-01 2014-10-28 The Regents Of The University Of California Systems and methods for blind source signal separation
WO2007103037A2 (en) * 2006-03-01 2007-09-13 Softmax, Inc. System and method for generating a separated signal
CN103763230B (zh) * 2014-01-17 2017-03-01 南京信息工程大学 一种改进的自适应盲源分离方法
CN104866866A (zh) * 2015-05-08 2015-08-26 太原理工大学 改进的自然梯度变步长盲源分离算法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102231280A (zh) * 2011-05-06 2011-11-02 山东大学 卷积语音信号的频域盲分离排序算法
CN102592607A (zh) * 2012-03-30 2012-07-18 北京交通大学 一种使用盲语音分离的语音转换系统和方法

Also Published As

Publication number Publication date
CN106887238A (zh) 2017-06-23

Similar Documents

Publication Publication Date Title
CN106887238B (zh) 一种基于改进独立向量分析算法的声信号盲分离方法
US9668066B1 (en) Blind source separation systems
US8874439B2 (en) Systems and methods for blind source signal separation
CN110634502B (zh) 基于深度神经网络的单通道语音分离算法
Shin et al. Normalised subband adaptive filter with variable step size
Nesta et al. Convolutive underdetermined source separation through weighted interleaved ICA and spatio-temporal source correlation
Lee et al. Single-channel speech enhancement method using reconstructive NMF with spectrotemporal speech presence probabilities
CN109166590B (zh) 一种基于空域相关性的二维时频掩模估计建模方法
JP6099032B2 (ja) 信号処理装置、信号処理方法及びコンピュータプログラム
Tengtrairat et al. Single-channel separation using underdetermined blind autoregressive model and least absolute deviation
Kemiha et al. Complex blind source separation
JP6448567B2 (ja) 音響信号解析装置、音響信号解析方法、及びプログラム
Das et al. ICA methods for blind source separation of instantaneous mixtures: A case study
KR101802444B1 (ko) 독립 벡터 분석 및 반향 필터 재추정을 이용한 베이시안 특징 향상에 의한 강인한 음성 인식 장치 및 방법
CN113095113B (zh) 一种用于水下目标识别的小波线谱特征提取方法及系统
CN108280416A (zh) 一种小波跨尺度相关滤波的宽带水声信号处理方法
CN109658944A (zh) 直升机声信号增强方法及装置
Jyoshna et al. An Efficient Reference Free Adaptive Learning Process for Speech Enhancement Applications.
Peng et al. Reducing permutation error in subband-based convolutive blind separation
Malek Blind compensation of memoryless nonlinear distortions in sparse signals
JP2014048398A (ja) 音響信号解析装置、方法、及びプログラム
JP6644356B2 (ja) 音源分離システム、方法及びプログラム
CN110956978B (zh) 一种基于欠定卷积混叠模型的稀疏盲分离方法
US11297418B2 (en) Acoustic signal separation apparatus, learning apparatus, method, and program thereof
Adiloğlu et al. A general variational Bayesian framework for robust feature extraction in multisource recordings

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
TA01 Transfer of patent application right

Effective date of registration: 20171026

Address after: 200050 Changning Road, Shanghai, No. 865, No.

Applicant after: Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences

Applicant after: University of Chinese Academy of Sciences

Address before: 200050 Changning Road, Shanghai, No. 865, No.

Applicant before: Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant