买专利,只认龙图腾
首页 专利交易 科技果 科技人才 科技服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

【发明授权】声源探测装置及方法、记录有声源探测程序的记录介质_松下电器(美国)知识产权公司_201810158701.6 

申请/专利权人:松下电器(美国)知识产权公司

申请日:2018-02-26

公开(公告)日:2020-11-24

公开(公告)号:CN108540898B

主分类号:H04R1/40(20060101)

分类号:H04R1/40(20060101)

优先权:["20171109 JP 2017-216735","20170303 US 62/466498"]

专利状态码:有效-授权

法律状态:2020.11.24#授权;2020.02.04#实质审查的生效;2018.09.14#公开

摘要:本公开涉及声源探测装置及方法、记录有声源探测程序的记录介质,对探测对象的声源的方向进行探测的声源探测装置1具备:相关矩阵算出部,其算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号;学习部,其以使得对预先存储在存储部中的多个第2相关矩阵分别乘以权重后的线性和等于第1相关矩阵的方式,通过学习来算出权重,所述多个第2相关矩阵是根据麦克风阵列的阵列排列而算出的按方向的相关矩阵;以及空间谱算出部,其使用通过学习部算出的权重,算出表示按方向的声压强度的、观测信号的空间谱。

主权项:1.一种声源探测装置,对探测对象的声源的方向进行探测,所述声源探测装置具备:存储部;和处理电路,其在运行中执行包括以下的工作:算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号,以使得对预先存储在所述存储部中的多个第2相关矩阵分别乘以权重后的线性和等于所述第1相关矩阵的方式,通过学习来算出所述权重,所述多个第2相关矩阵是根据所述麦克风阵列的阵列排列而算出的按方向的相关矩阵,使用所算出的所述权重,算出所述观测信号的空间谱,所述空间谱表示按方向的声压强度,所述工作还包括:选择第1元素和第2元素,所述第1元素是构成所述第1相关矩阵的元素中的一个元素,所述第2元素是构成所述多个第2相关矩阵的各第2相关矩阵的元素中的、位于与所述第1元素对应的位置的元素,逐次地切换要选择的所述第1元素以及所述第2元素,在所述权重的算出中,通过逐次地反复进行如下处理来算出所述权重:以使得对所述第2元素乘以第1权重后的第1元素线性和等于所述第1元素的方式,将所述第1权重更新为在所述权重的算出中算出的第2权重;以使得对下一个选择出的所述第2元素乘以更新后的所述第2权重后的第2元素线性和等于下一个选择出的所述第1元素的方式,将所述第2权重更新为在所述权重的算出中算出的第3权重。

全文数据:声源探测装置及方法、记录有声源探测程序的记录介质技术领域[0001]本公开涉及声源探测装置、声源探测方法以及记录有声源探测程序的记录介质。背景技术[0002]例如在日本特开2014-56181号公报中,提出了能够根据由多个麦克风单元获得的多个声响信号高精度地推定声源的方向的声源方向推定装置。在该日本特开2014-56181号公报中,通过使用基于多个声响信号的噪音信号的相关矩阵进行噪音应对措施,从而根据多个声响信号高精度地推定声源的方向。发明内容[0003]发明要解决的技术问题[0004]然而,在专利文献1中,基于作为观测信号的由多个麦克风单元获得的多个声响信号来算出噪音信号的相关矩阵。因此,在噪声源和探测对象的声源同时存在、或噪声的等级比探测对象的声源的等级高的情况下,难以准确地求出仅是噪声分量的相关矩阵。也就是说,存在如下问题:在基于由多个麦克风单元获得的多个声响信号的信号相位差进行声源探测的方式中,在存在比探测对象的声源高的声压级的噪声的情况下,会因该噪声的影响而导致不能对探测对象的声源进行检测即探测。[0005]本公开是鉴于上述的情况而完成的,其目的在于提供一种能够更切实地探测处于探测对象范围内的探测对象的声源的方向的声源探测装置。[0006]用于解决问题的技术方案[0007]本公开的一技术方案的声源探测装置,对探测对象的声源的方向进行探测,所述声源探测装置具备:存储部;和处理电路,其在运行中执行包括以下的工作:算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号,以使得对预先存储在所述存储部中的多个第2相关矩阵分别乘以权重后的线性和等于所述第1相关矩阵的方式,通过学习来算出所述权重,所述多个第2相关矩阵是根据所述麦克风阵列的阵列排列而算出的按方向的相关矩阵,使用所算出的所述权重,算出表示按方向的声压强度的、所述观测信号的空间谱。[0008]此外,这些中的一部分的具体技术方案既可以通过系统、方法、集成电路、计算机程序或计算机可读取的CD-ROM等记录介质来实现,也可以通过系统、方法、集成电路、计算机程序和记录介质的任意组合来实现。[0009]发明的效果[0010]根据本公开,能够实现能更切实地探测处于探测对象范围内的探测对象的声源的方向的声源探测装置等。附图说明[0011]图1是表示实施方式1的声源探测系统的构成的一例的图。[0012]图2是表示实施方式1的麦克风阵列与存在探测对象的声源的声源方向之间的位置关系的说明图。[0013]图3是在图2所示的位置关系下麦克风阵列观测的观测信号的空间谱图。[0014]图4是表示图1所示的声源探测装置的详细构成的一例的图。[0015]图5是实施方式1的选择部的选择方法的说明图。[0016]图6是表示实施方式1的非线性函数部的构成的一例的图。[0017]图7是表示实施方式1的声源探测装置的声源探测处理的流程图。[0018]图8是表示图7所示的声源探测处理的详细情况的流程图。[0019]图9是比较例的空间谱图。[0020]图10是实施方式1的空间谱图。[0021]图11是表示实施方式2的声源探测系统的构成的一例的图。具体实施方式[0022]本公开的一技术方案的声源探测装置,对探测对象的声源的方向进行探测,所述声源探测装置具备:存储部;和处理电路,其在运行中执行包括以下的工作:算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号,以使得对预先存储在所述存储部中的多个第2相关矩阵分别乘以权重后的线性和等于所述第1相关矩阵的方式,通过学习来算出所述权重,所述多个第2相关矩阵是根据所述麦克风阵列的阵列排列而算出的按方向的相关矩阵,使用所算出的所述权重,算出表示按方向的声压强度的、所述观测信号的空间谱。[0023]根据该构成,能够更切实地探测处于探测对象范围内的探测对象的声源的方向。进一步,因为使用通过学习算出的权重来算出观测信号的空间谱,所以能够实现抗噪声性以及对声音变化的跟随性优良的源探测装置。[0024]在此,例如也可以是,所述工作还包括:选择第1元素和第2元素,所述第1元素是构成所述第1相关矩阵的元素中的一个元素,所述第2元素是构成所述多个第2相关矩阵的各第2相关矩阵的元素中的、位于与所述第1元素对应的位置的元素,逐次地切换要选择的所述第1元素以及所述第2元素,在所述权重的算出中,通过所述逐次地反复进行如下处理来算出所述权重:以使得对所述第2元素乘以第1权重后的第1元素线性和等于所述第1元素的方式,将所述第1权重更新为在所述权重的算出中算出的第2权重;以使得对下一个选择出的所述第2元素乘以更新后的所述第2权重后的第2元素线性和等于下一个选择出的所述第1元素的方式,将所述第2权重更新为在所述权重的算出中算出的第3权重。[0025]由此,能够通过学习来算出使第1相关矩阵与多个第2相关矩阵的对应的各矩阵元素同时相等的权重,因此,基于由包括三个以上的麦克风单元的麦克风阵列收集到的声响信号,能够更切实地探测处于探测对象范围内的探测对象的声源的方向。[0026]另外,例如也可以是,在所述选择中,仅从构成所述第1相关矩阵以及所述第2相关矩阵的除对角成分之外的元素中的、由所述对角成分划分开的两组的多个元素的一方的组的多个元素中,选择所述第1元素以及所述第2元素。[0027]由此,能够削减运算量,因此能够更高速地探测处于探测对象范围内的探测对象的声源的方向。[0028]另外,例如也可以是,在所述权重的算出中,通过使用最小均方算法即LMS算法或独立成分分析即ICA,根据作为所述线性和与所述第1相关矩阵之差的误差和所述第2相关矩阵,算出所述权重。[0029]由此,能够在互相消除方向间的影响的同时算出按方向的强度,因此能够实现抗噪声性优良的声源探测装置。[0030]另外,例如也可以是,在所述权重的算出中,保持权重算出对所述多个第2相关矩阵分别乘以所保持的权重后的线性和,算出作为所述线性和与所述第1相关矩阵之差的误差,根据所述误差与所述第2相关矩阵之积来算出权重更新量,通过对所述保持的权重加上所述权重更新量来更新所述权重。[0031]在此,例如也可以是,在所述权重的更新中,通过使用LMS算法或ICA,根据所述误差以及所述第2相关矩阵来算出所述权重更新量。[0032]另外,例如也可以是,在所述权重的算出中,进一步,使用预定的非线性函数,对所述误差附加非线性性质,在所述更新量的算出中,根据附加了非线性性质的所述误差以及所述第2相关矩阵来算出所述权重更新量。[0033]由此,能够对所算出的误差赋予非线性性质并抑制方向间相互影响,因此能够实现抗噪声性优良的声源探测装置。[0034]另外,本公开的一技术方案的声源探测方法,对探测对象的声源的方向进行探测,所述声源探测方法包括:算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号,以使得对预先存储在存储部中的多个第2相关矩阵分别乘以权重后的线性和等于所述第1相关矩阵的方式,通过学习来算出所述权重,所述多个第2相关矩阵是根据所述麦克风阵列的阵列排列而算出的按方向的相关矩阵,使用所算出的所述权重,算出表示按方向的声压强度的、所述观测信号的空间谱。[0035]另外,本公开的一技术方案的记录介质,是存储有用于使计算机执行对探测对象的声源的方向进行探测的声源探测方法的程序的计算机可读取的非瞬时性记录介质,所述程序在由计算机执行时,使所述计算机执行包括以下的方法:算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号,以使得对预先存储在存储部中的多个第2相关矩阵分别乘以权重后的线性和等于所述第1相关矩阵的方式,通过学习来算出所述权重,所述多个第2相关矩阵是根据所述麦克风阵列的阵列排列而算出的按方向的相关矩阵,使用所算出的所述权重,算出表示按方向的声压强度的、所述观测信号的空间谱。[0036]此外,这些中的一部分的具体技术方案既可以通过系统、方法、集成电路、计算机程序或计算机可读取的CD-ROM等记录介质来实现,也可以通过系统、方法、集成电路、计算机程序和记录介质的任意组合来实现。[0037]以下,参照附图对本公开的一技术方案的声源探测装置进行具体说明。此外,以下说明的实施方式都表示本公开的一个具体例子。以下的实施方式中示出的数值、形状、材料、构成要素、构成要素的配置位置等是一例,并非限定本公开的意思。另外,关于以下的实施方式中的构成要素中的未记载在表示最上位概念的独立权利要求中的构成要素,作为任意的构成要素来说明。另外,在所有的实施方式中,也可以组合各实施方式的内容。[0038]实施方式I[0039]图1是表示实施方式1的声源探测系统1000的构成的一例的图。[0040]声源探测系统1000被用于对探测对象的声源的方向进行探测。在本实施方式中,声源探测系统1000如图1所示具备声源探测装置1、麦克风阵列200和频率分析部300。[0041][麦克风阵列200][0042]麦克风阵列200包括互相分离地配置的两个以上的麦克风单元,观测即收集从所有的方向到来的声波并输出将该声波变换成电信号后的声响信号。在本实施方式中,以下,设为麦克风阵列200由三个麦克风单元即麦克风单元201、202、203构成来进行说明。麦克风单元201、麦克风单元202以及麦克风单元203例如是对声压的灵敏度高的无指向性的麦克风元件,分离地配置换言之是配置于不同的位置)。在此,麦克风单元201输出将收集到的声波变换成电信号后的时域信号即声响信号mlη。同样地,麦克风单元202输出将收集到的声波变换成电信号后的时域的信号即声响信号m2η,麦克风单元203输出将收集到的声波变换成电信号后的时域的信号即声响信号m3η。[0043]图2是表示实施方式1中的麦克风阵列200与存在探测对象的声源S的声源方向之间的位置关系的说明图。图3是在图2所示的位置关系下麦克风阵列200观测的观测信号的空间谱图。图2中示出了由麦克风单元201、麦克风单元202以及麦克风单元203在θ=〇度的轴上排列成一列的阵列排列构成的麦克风阵列200的结构。另外,图2中不出了相对于麦克风阵列200在θ=θ8的方向上存在探测对象的声源S且不存在成为妨碍声的声源的情况。在该情况下,作为声源探测装置1的探测结果的空间谱成为如图3所示那样。即,在作为探测结果的图3所示的空间谱中,表示最高强度的角度成为0s。[0044][频率分析部300][0045]频率分析部300将在两个以上的麦克风单元各自中观测到的声响信号变换成频域的信号,作为频谱信号进行输出。更具体而言,频率分析部300对从麦克风阵列200输入的声响信号进行频率分析,输出作为频域信号的频谱信号。此外,在频率分析中,使用快速傅里叶变换(FastFourierTransform:FFT或离散傅里叶变换(DiscreteFourierTransfom=DFT等将时间信号变换成各频率分量的振幅信息和相位信息的方法即可。[0046]在本实施方式中,频率分析部300由进行快速傅里叶变换的FFT301、FFT302以及FFT303构成。FFT301将从麦克风单元201输出的声响信号mlη作为输入,使用快速傅里叶变换进行从时域向频域的变换而输出频谱信号SmlωFT302将从麦克风单元202输出的声响信号m2η作为输入,使用快速傅里叶变换进行从时域向频域的变换而输出频谱信号Sm2ωFT303将从麦克风单元203输出的声响信号m3η作为输入,使用快速傅里叶变换进行从时域向频域的变换而输出频谱信号Sm3ω。[0047][声源探测装置1][0048]图4是表示图1所示的声源探测装置1的详细构成的一例的图。[0049]声源探测装置1对探测对象的声源的方向进行探测。在本实施方式中,声源探测装置1如图1以及图4所示,具备相关矩阵算出部10、存储部20、选择部30、学习部40、空间谱算出部100以及输出部110。此外,如果构成麦克风阵列200的麦克风单元的个数为2个,则声源探测装置1也可以不具备选择部30。另外,也可以设为声源探测装置1具备麦克风阵列200以及频率分析部300。以下,对各构成要素进行说明。[0050]〈相关矩阵算出部10[0051]相关矩阵算出部10算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由麦克风阵列200收集到的声响信号。[0052]在本实施方式中,相关矩阵算出部10根据频率分析部300所输出的频谱,算出作为第1相关矩阵的观测相关矩阵Rxω。更具体而言,相关矩阵算出部10使用下述的(式1以及(式2,将来自FFT301的频谱信号Smlω、来自FFT302的频谱信号Sm2co以及来自FFT303的频谱信号Sm3ω作为输入,算出观测相关矩阵Rxω。[0053]在此,构成观测相关矩阵Rxω的各元素Xljω保存有相位差信息,该相位差信息是针对从存在于实际环境中的多个声源向各麦克风单元到来的多个声波的相位差信息。例如,(式1所示的元素X12ω表示针对向麦克风单元201以及麦克风单元202到来的声波的相位差信息。另外,例如式1所示的元素X13ω表示针对向麦克风单元201以及麦克风单元203到来的声波的相位差信息。(式2所示的(·Γ表示复共辄。[0056]此外,在本实施方式中,在作为麦克风单元201〜203示出的各麦克风单元的声压灵敏度特性大致相等且均匀的情况下,能够通过式3来表示观测相关矩阵Rxω的各元素Xljω。(式3所示的观测相关矩阵RxCO的各元素XljCO相当于对式2中的分母的标准化项进行了省略的各元素。[0057]Xijω=Smiω*Smjω…试3[0058]〈存储部20[0059]存储部20预先存储根据麦克风阵列200的阵列排列而算出的按方向的相关矩阵即多个第2相关矩阵。[0060]在本实施方式中,存储部20由存储器等构成,预先存储有作为第2相关矩阵的、各探测方向9的参照相关矩阵Rrθ,ω。在图4所示的例子中,存储部20预先存储有例如方向数N=ISO个的0兰Θ兰180的范围内的参照相关矩阵Rr仇,ω〜RrΘΝ,ω。[0061]参照相关矩阵Rrθ,ω表示与各方向Θ的声波对应的麦克风单元间的相位差,因此,如果声源的方向Θ和麦克风阵列200的麦克风单元排列即阵列排列已决定,则能够在理论上算出该参照相关矩阵。以下,以图2所示的麦克风阵列200的阵列排列为例,对算出参照相关矩阵Rrθ,ω的方法进行说明。[0062]在图2中,如上所述,示出了构成麦克风阵列200的麦克风单元201〜203呈直线状排列的阵列排列的例子。另外,图2中也示出了在方向0S上存在声源S这一位置关系。[0063]对于从声源S向麦克风单元201〜203的声波的到来时刻,若以中央的麦克风单元202为基准,则对于麦克风单元201而言提前时间τ,对于麦克风单元203而言延迟时间τ。时间τ能够使用以下的式4来算出。在式4中,L表示麦克风单元间距离,c表示音速。[0064]t=L·cos9sc···式4[0065]并且,对于表示与来自方向Θ的声源的声波对应的麦克风单元201〜203的相位差关系的方向矢量,若以位于中央的麦克风单元202的位置为基准,则可使用式5来表现。[0066][0067]因此,若换成声源位于方向Θ上时来说,则相对于方向Θ的参照相关矩阵Rrθ,ω能够根据式2、(式3以及式5的关系而定义成式6所示那样。在式6中,(·)η表示复共辄转置。[0068][0069]如此,能够算出相对于方向Q1-Qn例如N=180的参照相关矩阵Rr^1,ω〜RrΘν,ω〇[0070]〈选择部30[0071]选择部30选择第1元素以及第2元素,并逐次地切换要选择的第1元素以及第2元素,所述第1元素是构成第1相关矩阵的元素中的一个元素,所述第2元素是构成多个第2相关矩阵的各第2相关矩阵的元素中的位于与第1元素对应的位置的元素。在此,选择部30仅从构成第1相关矩阵以及第2相关矩阵的除对角成分之外的元素中的、由对角成分划分开的两组的多个元素的一方的组的多个元素中,选择第1元素以及第2元素即可。[0072]在本实施方式中,选择部30将来自相关矩阵算出部10的观测相关矩阵RxCO和来自存储部20的参照相关矩阵Rr0,ω作为输入,选择并输出观测相关矩阵RxCO以及多个参照相关矩阵Rr0,ω中的对应的相关矩阵的元素。选择部30例如如图4所示,具备矩阵元素选择部31和矩阵元素选择部32-1〜矩阵元素选择部32-Ν。此外,图4中示出了设置有被输入对于方向以的参照相关矩阵RrQ1,ω的矩阵元素选择部32-1和被输入对于方向ΘΝ的参照相关矩阵RrΘΝ,ω的矩阵元素选择部32-Ν这两个的情况下的例子,但并不限于此。在方向数N=180的情况下,设置有被输入对于方向Q1-Qn的参照相关矩阵Rr^1,ω〜RrΘΝ,ω的N个矩阵元素选择部32-1〜矩阵元素选择部32-Ν。[0073]以下,使用图5对选择部30的选择方法的一例进行具体说明。[0074]图5是实施方式1中的选择部30的选择方法的说明图。[0075]如图5所示,矩阵元素选择部31选择构成从相关矩阵算出部10输入的观测相关矩阵Rxω的元素也称为矩阵元素)中的一个元素,作为相位差信号Xω进行输出。矩阵元素选择部32-mm为1以上且N以下的自然数选择构成从存储部20输入的参照相关矩阵Rr0m,ω的元素中的、与矩阵元素选择部31所选择出的元素相同行且相同列的元素,作为相位差信号r9m,ω进行输出。[0076]此外,通常,相关矩阵的对角成分成为1而在信号处理上没有意义。另外,在相关矩阵中对行编号和列编号进行了调换的XU和Xjl,在相位旋转相反的关系下作为信息来说是相同的。考虑到这些情况,选择部30从构成参照相关矩阵Rrθ,ω以及观测相关矩阵Rxω的相关矩阵的除对角成分之外的元素中的、由对角成分划分开的两组的多个元素的一方的组的多个元素中选择元素并进行输出即可。也就是说,选择部30选择参照相关矩阵Rrθ,ω以及观测相关矩阵RxCO的除相关矩阵的对角成分之外的上三角矩阵或下三角矩阵的元素并进行输出即可。由此,声源探测装置1能够削减运算量。[0077]进一步,选择部30从运算量的削减等的观点出发,也可以间隔地选择上三角矩阵或下三角矩阵的元素。[0078]〈学习部40[0079]学习部40以使得对预先存储在存储部20中的多个第2相关矩阵分别乘以权重后的线性和(日文“線形和”)等于第1相关矩阵的方式,通过学习来算出该权重。在此,学习部40通过使用LMS算法(最小均方差)或ICAIndependentComponentAnalysis,独立成分分析),根据作为该线性和与第1相关矩阵之差的误差和第2相关矩阵来算出权重。更具体而言,学习部40以使得对由选择部30选择出的第2元素乘以第1权重后的第1元素线性和等于由选择部30选择出的第1元素的方式,将第1权重更新为通过学习而算出的第2权重。接着,学习部40以使得对下一个由选择部30选择出的第2元素乘以更新后的第2权重后的第2元素线性和等于下一个由选择部30选择出的第1元素的方式,将第2权重更新为通过学习而算出的第3权重。学习部40通过逐次反复进行这些更新,由此通过学习来算出该权重。[0080]在本实施方式中,学习部40如图1以及图4所示,具备保持部50、线性和算出部60、误差算出部70、非线性函数部80以及权重更新部90。此外,非线性函数部80并不是必须的结构,学习部40也可以不具备非线性函数部80。[0081]《保持部50》[0082]保持部50保持由权重更新部90更新的权重。保持部50保持有要对各参照相关矩阵Rrθ,ω乘以的权重。换言之,该权重对于构成参照相关矩阵Rr^1,ω〜RrΘΝ,ω各自的相关矩阵的各元素而言是共同的。[0083]另外,权重是以Θ以及ω为变量的函数,但通过ω作为常数来对待,其能够作为一次项的系数来对待。以下,将权重称为权重系数aθ,ω来进行说明。[0084]在本实施方式中,权重系数aθ,ω是对各方向Θ的参照相关矩阵Rrθ,ω乘以的系数。在图4中,作为一例,示出了与例如180个的05Θ5180的范围内的参照相关矩阵RrΘ,ω对应的方向Θ^ΘνN=180的权重系数a⑼,ω〜aΘΝ,ω。[0085]保持部50保持由权重更新部90更新的权重系数aθ,ω。也就是说,权重系数aθ,ω是基于由权重更新部90算出的权重更新量来对值进行更新的学习系数。另外,保持部50将保持的权重系数aθ,ω向空间谱算出部100输出。[0086]《线性和算出部60》[0087]线性和算出部60算出对多个第2相关矩阵分别乘以保持部50保持的权重后的线性和。[0088]在本实施方式中,线性和算出部60如图4所示具备信号乘法运算部61-1〜信号乘法运算部61-Ν和信号加法运算部62。[0089]信号乘法运算部61-1对由矩阵元素选择部32-1选择出的参照相关矩阵ω的元素He1,ω乘以方向0:的权重系数ω并输出到信号加法运算部62。同样地,信号乘法运算部61-Ν对由矩阵元素选择部32-Ν选择出的参照相关矩阵RrΘΝ,ω的元素rΘΝ,ω乘以方向ΘΝ的权重系数aΘΝ,ω并输出到信号加法运算部62。如此,信号乘法运算部61-1〜信号乘法运算部61-N分别将按各方向Q1WQn对参照相关矩阵Rrθ,ω乘以权重系数aθ,ω而得到的信号输出到信号加法运算部62。[0090]信号加法运算部62将对从信号乘法运算部61-1〜信号乘法运算部61-Ν输出的信号进行了加法运算而得到的推定相位差信号xrco输出到误差算出部70。更具体而言,信号加法运算部62使用式7,算出从信号乘法运算部61-1〜信号乘法运算部61-Ν输出的信号的线性和来作为推定相位差信号xrco。[0091][0092]《误差算出部70》[0093]误差算出部70算出作为由线性和算出部60算出的线性和与第1相关矩阵之差的误差。在本实施方式中,误差算出部70如图4所示,具备信号减法运算部71。[0094]信号减法运算部71通过从来自矩阵元素选择部31的相位差信号Xω减去来自信号加法运算部62的推定相位差信号χΓω,算出误差信号eco。更具体而言,信号减法运算部71使用式8来算出误差信号eω。[0095]eω=χω-xrω…(式8[0096]《非线性函数部80》[0097]非线性函数部80使用预定的非线性函数,对误差附加非线性性质。更具体而言,非线性函数部80将从信号减法运算部71输入的误差信号eω变换成通过具备非线性输入输出特性的函数即非线性函数附加了非线性性质的信号。非线性函数例如是双曲线函数,但并不限于此。只要是具有能够对信号振幅给予限制的非线性输入输出特性的非线性函数,则可以是任何函数。这是因为:即使在因外部干扰而导致相位差失常、误差信号eω暂时变大,也能够抑制对在后述的权重更新部90中学习的权重更新量产生的影响。[0098]图6是表示实施方式1的非线性函数部80的构成的一例的图。非线性函数部80如图6所示具备实部提取部801、虚部提取部802、非线性性质追加部803、非线性性质追加部804、虚数单位乘法运算部805以及信号加法运算部806。[0099]实部提取部801提取所输入的误差信号eω的实数部,输出到非线性性质追加部803。虚部提取部802提取所输入的误差信号eω的虚数部,输出到非线性性质追加部804。[0100]非线性性质追加部803对从实部提取部801输入的误差信号eω的实数部的信号振幅通过非线性函数附加非线性性质,并输出到信号加法运算部806。非线性性质追加部804对从虚部提取部802输入的误差信号eω的虚数部的信号振幅通过非线性函数附加非线性性质,并输出到虚数单位乘法运算部805。[0101]虚数单位乘法运算部805为了使从非线性性质追加部804输入的信号恢复为虚数而乘以虚数单位j,并输出到信号加法运算部806。信号加法运算部806将对作为实部信号的从非线性性质追加部803输入的信号和作为虚部信号的从虚数单位乘法运算部805输入的信号进行了加法运算后的附加了非线性性质的复信号fe«输出到权重更新部90。[0102]试9中示出附加了非线性性质的复信号feco的一例。(式9是对非线性函数使用双曲线函数tanh·的情况下的例子。real·表示实数部,imag·表示虚数部,j表示虚数单位。[0103]feω=tanhrealeω+j·tanhimageω···ΐζ9[0104]《权重更新部90》[0105]权重更新部90通过使用LMSLeastMeanSquare,最小均方)算法或ICAIndependentComponentAnalysis,根据误差以及第2相关矩阵来算出权重更新量,对保持部50保持的权重加上该权重更新量而作为保持部50保持的权重。另外,在声源探测装置1具备非线性函数部80的情况下,权重更新部90根据通过非线性函数部80附加了非线性性质的误差以及第2相关矩阵来算出权重更新量,对保持部50保持的权重加上该权重更新量而作为保持部50保持的权重。[0106]在本实施方式中,权重更新部90被输入从非线性函数部80输入的复信号feω和从选择部30输入的N个相位差信号rQ1,ω〜rΘΝ,ω。并且,权重更新部90算出与要对N个相位差信号rθι,ω〜rΘΝ,ω乘上的权重系数aQ1,ω〜aΘΝ,ω对应的权重更新量Aa9i,ω〜Aa%,ω。[0107]例如,在声源探测装置1不具备非线性函数部80的情况下,权重更新部90使用(式10,算出权重更新量AaQ^co〜Aa0N,co。另一方面,在声源探测装置1具备非线性函数部80的情况下,权重更新部90使用式11,算出权重更新量ΔaQ1,ω〜ΔaΘΝ,ω。[0108]Δa9k,ω=realβ·eω·r9k,ω…(式10[0109]Δa9k,ω=realβ·feω·r9k,ω…(式11[0110]此外,(式1〇以及式11示出了使用LMS算法来算出权重更新量的情况。β是控制更新速度的参数。另外,对于相关矩阵,在元素Η^ω*ηιω存在相位反转的关系。因此,在式10以及式11中,因复共辄的关系而虚部被消除,因此设置了real·的部分。[0111]并且,权重更新部90如下述的式12所示,使用所算出的权重更新量,更新保持部50所保持的权重系数a9k,ω。[0112]a9k,ω=a9k,ω+Δa9k,ω…(式12[0113]《空间谱算出部100》[0114]空间谱算出部100使用由学习部40算出的权重,算出表示按方向的声压强度的、观测信号的空间谱。[0115]在本实施方式中,空间谱算出部100将保持部50保持的、由权重更新部90通过学习而更新的权重系数aQ1,ω〜aθΝ,ω作为输入,算出空间谱ρΘ,输出到输出部11〇。[0116]更具体而言,空间谱算出部100如下述的(式13所示,通过将保持部50保持的权重系数aθ,ω针对频率ω来计算和或者平均值,能够得到空间谱ρΘ。原理将在后面叙述,但这是因为:权重系数aθ,ω作为表示各方向Θ及频率的声波的强度的函数来处理。[0117][0118][声源探测装置1的工作][0119]对如上述那样构成的声源探测装置1所进行的声源探测处理进行说明。[0120]图7是表示实施方式1的声源探测装置1的声源探测处理的流程图。[0121]首先,声源探测装置1进行观测信号的相关矩阵算出处理SlO。更具体而言,声源探测装置1算出观测信号的相关矩阵即观测相关矩阵Rx«,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列200收集到的声响信号。[0122]接着,声源探测装置1进行对参照相关矩阵分别乘以的权重的学习处理(S20。更具体而言,声源探测装置1以使得对预先存储在存储部20中的多个参照相关矩阵Rrθ,ω分别乘以权重系数aθ,ω后的线性和等于观测相关矩阵RxCO的方式,通过学习来算出权重,所述多个参照相关矩阵Rrθ,ω是根据麦克风阵列的阵列排列而算出的按方向的相关矩阵。[0123]接着,声源探测装置1进行观测信号的空间谱算出处理S30。更具体而言,声源探测装置1使用在步骤S20中算出的权重,算出表示按方向的声压强度的、观测信号的空间谱P[0124]图8是表示图7所示的声源探测处理的详细情况的流程图。对与图7同样的要素标注相同的标号。[0125]更详细而言,首先,在步骤SlO中,麦克风阵列200取得时刻t的声响信号SlOl。接着,频率分析部300进行通过步骤SlOl取得的声响信号的频率分析S102,变换成作为频域信号的频谱信号。然后,声源探测装置1根据在步骤S102中变换得到的频谱信号,算出时刻t的观测信号的相关矩阵即观测相关矩阵Rx«S103。[0126]接着,在步骤S20中,首先,对声源探测装置1设定预定次数Nt来将作为进行权重的学习处理的次数S201。接着,声源探测装置1选择观测相关矩阵RxCO以及参照相关矩阵Rr0,ω的对应的矩阵元素,输出相位差信号χω以及相位差信号Γθ,ωS202。接着,声源探测装置1根据相位差信号χω、相位差信号rθ,ω以及权重系数a0,ω,算出误差信号eωS203。接着,声源探测装置1算出对误差信号eω附加了非线性性质的复信号feωS204〇接着,声源探测装置1根据通过步骤S204算出的复信号feω和通过步骤S203算出的相位差信号rθ,ω,算出权重系数aθ,ω的权重更新量Δaθ,ω,更新权重系数a0,ωS205。然后,声源探测装置1判断通过步骤S202选择出的观测相关矩阵Rxω以及参照相关矩阵Rrθ,ω的矩阵元素是否轮完一次S206。在轮完一次的情况下S206:是),判断进行权重系数aθ,ω的学习处理的次数是否达到了预定次数NtS207。在达到了预定次数Nt的情况下(S207:是),声源探测装置1进入接下来的步骤S30。此外,在步骤S206或步骤S207中未轮完一次的情况S206:否或未达到预定次数Nt的情况下(S207:否),返回到步骤S202,反复进行处理。[0127]接着,在步骤S30中,声源探测装置1根据通过步骤S20中的学习而更新的权重系数a0,ω,算出观测信号的空间谱ρθS301。[0128]接着,声源探测装置1在步骤S40中,将时刻t更新为例如时刻t+Δt,在步骤S50中,判断是否结束声源探测处理。此外,在不结束声源探测处理的情况下(S50:否),返回到步骤S10,算出时刻t+At的观测信号的相关矩阵即观测相关矩阵RxCO。[0129]如此,声源探测装置1以使得对多个参照相关矩阵Rrθ,ω分别乘以权重系数aθ,ω后的线性和等于观测相关矩阵RxCO的方式,反复进行对应的各矩阵元素的学习,直到学习到关于全部矩阵元素的权重系数aθ,ω。进一步,声源探测装置1也可以按预定次数Nt来反复进行学习。[0130]例如,若是3行3列的参照相关矩阵Rrθ,ω以及观测相关矩阵RxCO、且预定次数Nt为3次,则会针对上三角矩阵或下三角矩阵的3个元素进行3次学习处理,因此总共会进行9次学习处理。[0131]如此,能够学习到使对多个参照相关矩阵Rr0,co分别乘以权重系数θ,ω后的线性和与观测相关矩阵Rxco更加相等的权重系数aθ,ω。[0132][工作的原理][0133]接着,对能够通过学习来算出使对多个参照相关矩阵Rrθ,ω分别乘以权重系数aθ,ω后的线性和等于观测相关矩阵Rxω的权重系数aθ,ω的原理进行说明。另外,对能够使用所获得的权重系数aθ,ω来算出空间谱ρΘ的原理也进行说明。[0134]对于基于来自麦克风阵列200的信号而观测的观测相关矩阵RxCO、即作为相关矩阵算出部10的输出的观测相关矩阵RxW,如下述的式14所示,可知:能够用与存在于方向Θ上的空间的声源对应的相关矩阵Rsθ,ω和强度uθ,ω的线性和来近似。Rsθ,ω是因声波的到来方向产生的麦克风单元间的相位差信息,表示方向信息。强度uθ,ω表示声波的强度。并且,通过求出与各方向Θ的声波对应的强度uθ,ω,能够导出空间谱ρΘ。[0135][0136]在式14中,观测相关矩阵RxCO是能够观测的相关矩阵,是已知数。另一方面,强度uθ,ω以及相关矩阵Rsθ,ω是未知数。在此,相关矩阵Rsθ,ω是按方向Θ的相关矩阵,其矩阵元素是声波到来方向为方向Θ时的麦克风单元间的相位差。着眼于该情况,相关矩阵Rsθ,ω使用作为已知信息的麦克风阵列的麦克风单元排列、方向Θ以及音速c,能够替换为理论值。此外,上述的(式4、(式5以及式6是将相关矩阵Rsθ,ω替换成使用已知信息预先计算出的、作为理论值的参照相关矩阵Rrθ,ω而得到的。[0137]另外,通过使在声源探测装置1中作为空间谱求出的未知数等于权重系数aΘ,ω、即使权重系数aθ,ω等于试14的强度uθ,ω,由此试14能够替换为试15。[0138][0139]因此,对于算出(式15,因为观测相关矩阵RxCO是观测值、且参照相关矩阵Rrθ,ω是已知的理论值,所以成为求出权重系数aθ,ω的问题。此外,这样的问题也被称为半盲识别(Semi-blindIdentification,日文K同定”)的问题。[0140]在此,与通常的声响信号的识别不同之处在于:观测相关矩阵RxCO和参照相关矩阵Rr0,ω是矩阵、且权重系数aθ,ω是一次项的系数;与观测信号和参照信号相当的信号是表现相位差的转子、且总是振幅1的复数。[0141]根据观测相关矩阵Rxω和参照相关矩阵Rrθ,ω是矩阵、且权重系数aθ,ω是一次项的系数这一点,可知:对于权重系数aθ,ω,会成为求出对观测相关矩阵Rxco和参照相关矩阵Rrθ,ω的对应的各矩阵元素而言共同成立的值。也就是说,会成为求出满足用矩阵的元素改写了(式15而得到的式16的权重系数aθ,ω。在式16中,Xljω是观测相关矩阵Rxω的矩阵元素。rijθ,ω表示参照相关矩阵Rrθ,ω的矩阵元素。[0142][0143]在本实施方式中,将式16改写为式17,使用使作为推定误差的误差信号θω最小化的LMS或ICAIndependentComponentAnalysis等学习方式,由此求出a9,co。此夕卜,学习方式并不限定于这些方式。[0144][0145]更具体而言,在式17中,为了算出对Xljω以及rijθ,ω的矩阵元素共同成立的权重系数aθ,ω,通过选择部30依次选择矩阵元素来进行权重系数的学习。并且,信号乘法运算部61-1、……、61-Ν对应于式17的右边第2项的乘法运算,信号加法运算部62对应于式17的右边的Σ,信号减法运算部71对应于式17的右边的减法运算。[0146]另外,根据与观测信号和参照信号相当的信号是表现相位差的转子、且常是振幅1的复数这一点,能对误差信号eω附加非线性性质而添加独立分量分析ICA的效果以抑制方向间相互影响。[0147]在本实施方式中,在如图6所示那样分解成实部和虚部之后,适用上述的(式9那样的非线性函数。如此,能够将作为声波到来方向的方向Θ的差异作为独立的分量来学习,因此,能够获得不容易受到不同方向的影响的收敛动作。[0148]基于以上那样的考虑,使用式11以及式12来进行权重系数的更新。并且,通过使用学习后的权重系数a0,ω并使用式13,能够算出作为声源探测装置1的输出的空间谱P⑻。[0149][效果][0150]如上所述,根据本实施方式的声源探测装置1,能够基于由作为麦克风阵列200的多个元件的麦克风单元观测到的声响信号的观测相关矩阵Rxω来求出空间谱ρΘ。更具体而言,预先准备能够根据麦克风阵列200的阵列排列而计算为理论值的按方向的参照相关矩阵Rrθ,ω,以使得对表示各个方向的参照相关矩阵Rrθ,ω分别乘以权重系数aΘ,ω后的线性和等于观测相关矩阵Rxω的方式,通过学习来算出权重系数aθ,ω。并且,使用所获得的权重系数aθ,ω来算出空间谱ρθ。由此,不进行根据运算量大的相关矩阵和方向矢量导出空间谱的计算,而能够逐次地推定与成为妨碍声的声源以及探测对象的声源的方向对应的强度来作为权重系数aθ,ω,因此,能够基于按频率分析帧单位的毫秒〜秒级这样的间隔在麦克风单元中观测到的声响信号的观测相关矩阵Rx«来求出空间谱ρΘ。也就是说,根据本实施方式,可知能够实现在对声音变化的跟随性上优良的声源探测装置1。[0151]另外,根据本实施方式的声源探测装置1,能够在使方向间的影响互相消除的同时算出按方向的强度。例如,将θΐ〜的角度范围设为应检测的探测范围的方向,将θω+1〜0咖角度范围设为存在妨碍声的非探测范围的方向。并且,将式15变形成式18那样,以使得应检测的探测范围来到左边,存在妨碍声的非探测范围来到右边。[0152][0153]于是,可知:(式18的左边相当于与作为声源探测结果获得的空间谱对应的相关矩阵。可知:(式18的右边的第1项相当于被观测的全部方向的声波混合存在的观测相关矩阵,试18的右边的第2项相当于表示妨碍声分量的相关矩阵。另外,可知:通过式18的右边的减法运算,妨碍声分量的相关矩阵被从观测相关矩阵Rxω中减掉,可获得消除效果。对于这种情况,可知:由于成为对各方向Θ的各分量互相消除干涉的工作,因此抗噪声性能提高。另外,可知:由于针对与全部方向对应的权重系数aθ,ω同时进行求解,因此在对声音变化的跟随性上也优良。[0154]因此,本实施方式的声源探测装置1通过根据探测范围的权重系数aθ,ω来算出空间谱PΘ,能够实现在抗噪声性能和对声音变化的跟随性上优良的声源探测。[0155]如上所述,根据本实施方式的声源探测装置1,能够更切实地探测处于探测对象范围内的探测对象的声源的方向。进而,本实施方式的声源探测装置1通过使用权重系数aΘ,ω算出空间谱pΘ,能够发挥抗噪声性以及对声音变化优良的跟随性。[0156]在此,使用图9以及图10,对本实施方式的声源探测装置1的效果进行说明。[0157]图9是比较例的空间谱图。在图9中,作为比较例示出了在存在探测对象的声源S和在声源S的附近成为声源S的妨碍声的声源Nl以及声源Ν2的情况下,使用专利文献1的技术算出了空间谱时的图。[0158]在图9所示的空间谱中,作为妨碍声的声源Nl的强度,不仅呈现在存在声源Nl的方向上,还呈现为其随着从声源Nl的方向(角度)离开而衰减。作为妨碍声的声源Ν2的强度也呈现与声源Nl同样的特性。因此,如图9所示,在声源Nl和声源Ν2的声压级比探测对象的声源S的声压级高的情况下,会成为声源S的强度的峰值被作为妨碍声的2个声源Nl和声源Ν2的强度的峰值埋没的状态。因此,在比较例中,无法检测探测对象的声源S的存在即声源S的强度的峰值,所以无法探测声源S的方向。[0159]另一方面,图10是实施方式1的空间谱图。在图10中,也示出了在存在探测对象的声源S和在声源S的附近成为声源S的妨碍声的声源Nl以及声源Ν2的情况下,本实施方式的声源探测装置1算出了空间谱时的图。声源探测装置1使用权重系数aθ,ω算出空间谱ρΘ,因此能够对各方向9的各分量互相消除干涉。因此,如图10所示,与声源Nl和声源Ν2的声压级比探测对象的声源S的声压级高还是低无关地,会成为声源S的强度的峰值、作为妨碍声的2个声源Nl和声源Ν2的强度的峰值独立出现的状态。也就是说,能够同时独立地探测作为妨碍声的声源Nl、声源Ν2以及探测对象的声源S的强度的峰值。[0160]因此,根据本实施方式的声源探测装置1,可知能够更切实地探测处于探测对象范围内的探测对象的声源的方向。[0161]此外,对于相关矩阵算出部10算出的观测相关矩阵RxCO以及存储部20所存储的各探测方向9的参照相关矩阵Rrθ,ω,也可以将在运算中使用的相关矩阵的上三角矩阵的元素或任意选择出的元素用矢量的形式来实现。该情况下,选择部30依次选择矢量的元素并输出即可。[0162]另外,在本实施方式中,在各方向Θ的参照相关矩阵Rrθ,ω和权重系数aθ,ω,作为方向数N为180个进行了说明,但并不限于此。根据声源探测装置1的用途、麦克风阵列的规模或运算规模,方向数N既可以多也可以少,并没有特别限制。另外,要设定的角度间隔既可以均匀,也可以不均。[0163]另外,在本实施方式中,在各频率ω的观测相关矩阵RxCO、参照相关矩阵RrΘ,ω以及权重系数aθ,ω,没有对频率ω的范围进行特别限制,但也可以根据探测对象的声源所包含的频率成分来对频率ω的范围设置限制。[0164]实施方式2[0165]在实施方式1中,对使用学习后的权重系数aθ,ω来算出空间谱ρΘ的情况进行了说明,但并不限于此。也可以使用学习后的权重系数aθ,ω来算出从所指定的方向到来的声响信号波形。以下,将该情况作为实施方式2来进行说明。[0166]图11是表示实施方式2的声源探测系统1000Α的构成的一例的图。声源探测系统1000Α相当于利用了声源探测装置的麦克风装置。对与图1以及图4同样的要素标注相同的标号,省略详细的说明。[0167]图11所示的声源探测系统1000Α相当于实施方式1的声源探测系统1000,声响信号谱算出部100Α、输出部IlOA以及IFFT120的构成不同。[0168][声响信号谱算出部100Α][0169]声响信号谱算出部100Α将保持部50所保持的权重系数aθ,ω、与来自麦克风单元201的声响信号mlη对应的频谱信号SmlCo、以及为了进行信号取得而指定的方向即方向θ〇作为输入,算出输出的声响信号谱Υω。[0170]更具体而言,声响信号谱算出部100Α使用式19,算出声响信号谱Υω。[0171]Yco=a9〇,coSmlω…(式19[0172]此外,从声源探测的角度分辨率的观点出发,根据麦克风阵列200的尺寸或麦克风单元数,也可以将与进行指定的方向相邻的角度的权重系数如式20所示那样进行加法运算来使用。[0173][0174]对于式19以及式20中的权重系数aθ,ω,如上述工作原理中所述那样,表现与各方向Θ的声波对应的强度,因此表现Θ方向的谱相对于全部方向的谱的强度比率。因此,通过乘以全部方向的频谱Smlω,能够算出与从指定的方向θ〇到来的声波对应的声响信号谱Υω。[0175][IFFT120][0176]IFFTInverseFastFourierTransform,快速傅里叶逆变换)120算出对由声响信号谱算出部IOOA算出的声响信号谱Yω进行了快速傅里叶逆变换而得到的声响信号波形yn,并输出到输出部110Α。[0177][效果][0178]如上所述,根据本实施方式的声源探测系统1000A,使用在抗噪声性优良的声源探测装置中通过学习算出的权重系数aθ,ω,能够仅输出所指定的特定方向的声响信号波形yη。由此,能够实现仅提取特定方向的声音的麦克风装置的功能。[0179]以上,基于实施方式以及变形例对本公开的一个或多个技术方案的声源探测装置等进行了说明,但本公开并不限定于这些实施方式等。在不脱离本公开主旨的范围内,在本实施方式中实施本领域技术人员能想到的各种变形而得到的方式、组合不同实施方式中的构成要素而构成的方式,也可以包含在本公开的一个或多个技术方案的范围内。例如,以下这样的情况也包含在本公开中。[0180]1上述的声源探测装置等具体也可以是由微处理器、R0M、RAM、硬盘单元、显示器单元、键盘、鼠标等构成的计算机系统。在所述RAM或硬盘单元中存储有计算机程序。通过所述微处理器按照所述计算机程序进行工作,各构成要素实现其功能。在此,计算机程序是为了实现预定功能而组合多个表示针对计算机的指令的命令码而构成的。[0181]2构成上述的声源探测装置等的构成要素的一部分或全部,也可以包括一个系统LSILargeScaleIntegration:大规模集成电路)。系统LSI是将多个构成部集成到一个芯片上而制造出的超多功能LSI,具体是包括微处理器、R0M、RAM等而构成的计算机系统。在所述RAM中存储有计算机程序。通过所述微处理器按照所述计算机程序进行工作,系统LSI实现其功能。[0182]3构成上述的声源探测装置等的构成要素的一部分或全部也可以包括能够装卸于各装置的IC卡或单体模块。所述IC卡或所述模块是包括微处理器、R0M、RAM等的计算机系统。所述IC卡或所述模块也可以包括上述的超多功能LSI。通过微处理器按照计算机程序进行工作,所述IC卡或所述模块实现其功能。该IC卡或该模块可以具有防篡改性能。[0183]产业上的可利用性[0184]本公开能够利用于使用了多个麦克风单元的声源探测装置,特别是,能够利用于如下的声源探测装置,该声源探测装置能够更切实地探测位于距声源探测装置较远位置的遥控直升机或无人机等到达麦克风单元的声音比周围的声音小的声源的方向。

权利要求:1.一种声源探测装置,对探测对象的声源的方向进行探测,所述声源探测装置具备:存储部;和处理电路,其在运行中执行包括以下的工作:算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号,以使得对预先存储在所述存储部中的多个第2相关矩阵分别乘以权重后的线性和等于所述第1相关矩阵的方式,通过学习来算出所述权重,所述多个第2相关矩阵是根据所述麦克风阵列的阵列排列而算出的按方向的相关矩阵,使用所算出的所述权重,算出所述观测信号的空间谱,所述空间谱表示按方向的声压强度。2.根据权利要求1所述的声源探测装置,所述工作还包括:选择第1元素和第2元素,所述第1元素是构成所述第1相关矩阵的元素中的一个元素,所述第2元素是构成所述多个第2相关矩阵的各第2相关矩阵的元素中的、位于与所述第1元素对应的位置的元素,逐次地切换要选择的所述第1元素以及所述第2元素,在所述权重的算出中,通过所述逐次地反复进行如下处理来算出所述权重:以使得对所述第2元素乘以第1权重后的第1元素线性和等于所述第1元素的方式,将所述第1权重更新为在所述权重的算出中算出的第2权重;以使得对下一个选择出的所述第2元素乘以更新后的所述第2权重后的第2元素线性和等于下一个选择出的所述第1元素的方式,将所述第2权重更新为在所述权重的算出中算出的第3权重。3.根据权利要求2所述的声源探测装置,在所述选择中,仅从构成所述第1相关矩阵以及所述第2相关矩阵的除对角成分之外的元素中的、由所述对角成分划分开的两组的多个元素的一方的组的多个元素中,选择所述第1元素以及所述第2元素。4.根据权利要求1〜3中任一项所述的声源探测装置,在所述权重的算出中,通过使用最小均方算法或独立成分分析,根据作为所述线性和与所述第1相关矩阵之差的误差和所述第2相关矩阵,算出所述权重。5.根据权利要求1〜3中任一项所述的声源探测装置,在所述权重的算出中,保持权重,算出对所述多个第2相关矩阵分别乘以所保持的权重后的线性和,算出作为所述线性和与所述第1相关矩阵之差的误差,根据所述误差与所述第2相关矩阵之积来算出权重更新量,通过对所述所保持的权重加上所述权重更新量来更新所述权重。6.根据权利要求5所述的声源探测装置,在所述权重的更新中,通过使用最小均方算法或独立成分分析,根据所述误差以及所述第2相关矩阵来算出所述权重更新量。7.根据权利要求5或6所述的声源探测装置,在所述权重的算出中,进一步,使用预定的非线性函数,对所述误差附加非线性性质,在所述更新量的算出中,根据附加了非线性性质的所述误差以及所述第2相关矩阵来算出所述权重更新量。8.—种声源探测方法,对探测对象的声源的方向进行探测,所述声源探测方法包括:算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号,以使得对预先存储在存储部中的多个第2相关矩阵分别乘以权重后的线性和等于所述第1相关矩阵的方式,通过学习来算出所述权重,所述多个第2相关矩阵是根据所述麦克风阵列的阵列排列而算出的按方向的相关矩阵,使用所算出的所述权重,算出所述观测信号的空间谱,所述空间谱表示按方向的声压强度。9.一种记录介质,是存储有用于使计算机执行对探测对象的声源的方向进行探测的声源探测方法的程序的计算机可读取的非瞬时性记录介质,所述程序在由计算机执行时,使所述计算机执行包括以下的方法:算出观测信号的相关矩阵即第1相关矩阵,所述观测信号是由包括互相分离地配置的两个以上的麦克风单元的麦克风阵列收集到的声响信号,以使得对预先存储在存储部中的多个第2相关矩阵分别乘以权重后的线性和等于所述第1相关矩阵的方式,通过学习来算出所述权重,所述多个第2相关矩阵是根据所述麦克风阵列的阵列排列而算出的按方向的相关矩阵,使用所算出的所述权重,算出所述观测信号的空间谱,所述空间谱表示按方向的声压强度。

百度查询: 松下电器(美国)知识产权公司 声源探测装置及方法、记录有声源探测程序的记录介质

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。