CN113068112B - 声场重现中仿真系数向量信息的获取算法及其应用 - Google Patents

声场重现中仿真系数向量信息的获取算法及其应用 Download PDF

Info

Publication number
CN113068112B
CN113068112B CN202110226269.1A CN202110226269A CN113068112B CN 113068112 B CN113068112 B CN 113068112B CN 202110226269 A CN202110226269 A CN 202110226269A CN 113068112 B CN113068112 B CN 113068112B
Authority
CN
China
Prior art keywords
filter
gain
algorithm
sound field
coefficient vector
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
CN202110226269.1A
Other languages
English (en)
Other versions
CN113068112A (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.)
Shenzhen Yueersheng Acoustics Co ltd
Original Assignee
Shenzhen Yueersheng Acoustics Co ltd
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 Shenzhen Yueersheng Acoustics Co ltd filed Critical Shenzhen Yueersheng Acoustics Co ltd
Priority to CN202110226269.1A priority Critical patent/CN113068112B/zh
Publication of CN113068112A publication Critical patent/CN113068112A/zh
Application granted granted Critical
Publication of CN113068112B publication Critical patent/CN113068112B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
    • H04S7/30—Control circuits for electronic adaptation of the sound field
    • H04S7/302—Electronic adaptation of stereophonic sound system to listener position or orientation
    • H04S7/303—Tracking of listener position or orientation
    • H04S7/304—For headphones
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
    • H04S7/30—Control circuits for electronic adaptation of the sound field
    • H04S7/302—Electronic adaptation of stereophonic sound system to listener position or orientation
    • H04S7/303—Tracking of listener position or orientation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S2400/00—Details of stereophonic systems covered by H04S but not provided for in its groups
    • H04S2400/01—Multi-channel, i.e. more than two input channels, sound reproduction with two speakers wherein the multi-channel information is substantially preserved
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S2400/00—Details of stereophonic systems covered by H04S but not provided for in its groups
    • H04S2400/11—Positioning of individual sound objects, e.g. moving airplane, within a sound field
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S2420/00—Techniques used stereophonic systems covered by H04S but not provided for in its groups
    • H04S2420/01—Enhancing the perception of the sound image or of the spatial distribution using head related transfer functions [HRTF's] or equivalents thereof, e.g. interaural time difference [ITD] or interaural level difference [ILD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明属于耳机技术领域，提供了一种声场重现中仿真系数向量信息的获取算法及其应用。所述声场重现中仿真系数向量信息的获取算法，采用大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法，应用的声场范围更大，可以达到‑105度～+105度，有更大的声场动态范围和高仿真效果。本发明的算法，既保证了真实性，也实现了平滑过渡；基于最小二乘‑多项式拟合的系数估计算法，可以实现更小粒度的平滑。

Description

声场重现中仿真系数向量信息的获取算法及其应用
技术领域
本发明属于耳机技术领域，尤其涉及一种声场重现中仿真系数向量信息的获取算法及其应用。
背景技术
HRTF(Head Related Transfer Function，头部相关传递函数)，是一种音效定位算法，目前已初步应用到头戴耳机领域。
围绕如何更准确的获得头部相关传递函数(HRTF)，采取的一些测量或仿真方法。目前SONY公司的专利公开了一种HRTF测量方法、HRTF测量装置和程序专利(申请号201780075539.5)，该专利可以实现HRTF测量装置显示表示用户应面对的目标方向的图像。在用户的正面方向与目标方向匹配的情况下，HRTF测量装置通过扬声器输出测量声音并且基于由佩戴在用户耳朵上的麦克风获取的测量声音的结果来测量HRTF。该专利技术的重点是基于参考图像的方法，可以获得具有个体差异的头部相关传递函数；本技术可以应用于例如测量头部相关传递函数等的装置，针对个体差异的HRTF测量技术，旨在针对每个人的个体声学结构差异引起的频响差异，进行个性化测量，使声场回放更逼真。但因为批量采集信息难以实现，所以很难大规模的应用在相关产品上。
关于基于HRTF技术的应用，目前应用集中在两个方面：一为录音重放；一为以头戴式耳机为载体的环绕音效重现技术。目前哈曼公司公开了一种具有头部跟踪的双耳头戴式耳机呈现(申请号201611243763.4)，该发明公开一种可增强由头戴式耳机和其他声音系统发出的声音的再现的声音增强系统(SES)。所述SES通过模拟期望的声音系统而不包括通常与声音系统模拟相关联的不想要的伪像来改进声音再现。所述SES借助于通过从直接和间接头部相关传递函数(HRTF)导出的一组一个或多个双耳呈现滤波器变换声音系统输出来促进此类改进。所述双耳呈现滤波器的参数基于佩戴所述头戴式耳机的用户的头部跟踪角度来更新，以便呈现稳定的立体声像。所述头部跟踪角度可根据从安装在头戴式耳机组件中的数字陀螺仪获得的传感器数据来确定。该技术方案在-45度～+45度范围内，采用限定坡型滤波器(shelf filter)和陷波滤波器(notch filter)两种拓扑结构类型和数量，根据陀螺仪采集的角度信息，固定中心频率不变，只更新增益和延时2个参数，以此提高滤波器切换的平滑性，避免听觉上的“咔哒”声，算法实现简单，容易实现滤波器的平滑切换，计算量小，在计算资源有限的情况下，工程上也容易实现。但是其缺点是这一切换算法和策略与实际声音的传播特性相差较大，听感上与实际的声场特性会有一定差别；主要原因是声源传播到人耳(或仿真人头，实际测量HRTF时多用仿真耳)的传播特性，其频响曲线是较复杂的，要用10几段级联的滤波器才能逼近。y为准确描述滤波器的波形特征，其相关参数，除前文提到的中心频率和增益外，通常还包括滤波器的拓扑类型、抬升增益、品质因数等多个参数。这些参数在声源发生角度或位置变化时，通常都会跟随变化，特别是大角度变化时，这些参数变化会很剧烈，呈现非线性变化特点，在实际中会出现若干“奇异点”，是不连续和非线性的，上述两种滤波器的拓扑结构在克服“奇异点”方面效果并不理想。例如滤波器拓扑类型，由一种类型衍化到另一种类型，这些非线性变化特点，是影响听感的重要因素。并且该专利记述的技术方案主要应用于-45度～+45度范围内，而在用户实际佩戴使用双耳头戴式耳机过程中，头部及其上半身配合的转动幅度经常在-45度～+45度范围以外，所以上述专利公开的音效切换方案在重现声场音效的真实感受上是有较大局限的。
目前声场重现中，多是用耳机模仿一定距离和角度的外部声源，如扬声器在声场变化下的环绕声效果。在获得环绕音效体验的同时，由于距离和角度变化带来频率响应的显著变化，同时也失去了原有耳机高保真音频的效果，不能不说是一种遗憾。本发明声场重现中仿真系数向量信息的获取算法，弥补了现有技术中在小幅度角度幅度内进行系数拟合的缺陷。
发明内容
本发明实施例提供一种声场重现中仿真系数向量信息的获取算法及其应用，旨在解决目前在高仿真声场重现过程中，在小幅度角度幅度内进行系数拟合的缺陷。
本发明实施例是这样实现的，本发明提供了一种声场重现中仿真系数向量信息的获取算法，采用大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法，具体包括以下步骤：
(1)大角度范围的离散工程测量，采用HRTF(头部相关传递函数)的测量方法，在-105～+105度范围内，从初始零位开始，以15度为测量步长，测得测量步长节点的滤波器参数再由数字信号处理软件(Digital Signal Processing，简写为DSP软件)拟合出参数型级联的多段滤波器，并将所述滤波器参数以测量的角度为单位保存，供出现从一种拓扑类型衍化到另一种拓扑类型的节点(即奇异点)时角度切换时调用，作为小角度范围插值平滑过渡的起算数据；
(2)小角度范围插值平滑过渡，采用插值算法，在每个15度测量步长的扇形区间内，再等间隔采样，采用基于最小二乘法-多项式拟合的逼近算法，估计出多项式系数；通过系数拟合估计滤波器参数，可以对角度变化实现连续的平滑过渡；当遇到拓扑结构的奇异点时，则重新起算多项式系数，然后通过系数拟合估计滤波器参数；结合对应扇形区域、滤波器的频段号、滤波器拓扑类型信息获得多项式拟合出的仿真系数向量组。
步骤(2)中所述插值算法，具体步骤如下：在-105～+105度范围内，对每个15度测量步长的扇形区间，从初始零位开始区分左右编成L_0～L_6，R_0～R_6索引号，以每个扇形区间的每级滤波器的拓扑类型-type为基准参数，对中心频率-f0、增益-gain、抬升增益-boost、品质因数-Q参数分别进行二阶多项式拟合输出相应的系数；当遇到拓扑结构的奇异点时，则重新起算多项式系数。
所述的二阶多项式拟合，具体如下：
y＝a0+a1·x+a2·x2
其中：
x：为测得的角度值；
y：为需要获得的滤波器参数，为f0、gain、boost、Q中的任意一个；
其中a0，a1，a2为多项式系数；
A＝[a0,a1,a2]T，表示需要估计出的多项式系数向量；其中T为滤波器拓扑类型。
考虑多参数下，步骤(2)中所述拟合出的仿真系数向量组包括以下7个维度：
coeff＝[bin_num,band_num,fil_T,A_f0,A_Q,A_boost,A_gain]
其中：
bin_num，表示扇形区域的编号，即；L_0～L_6，R_0～R_6，共14个元素；
band_num，表示多段滤波器的频段号；选用最多级数，作为参考编号；本发明中为band1～band12，即12个元素；
fil_T,表示滤波器拓扑类型；本发明中包括3个滤波器拓扑类型元素：
highshelf--高架滤波器，对应编码：1；
lowshelf--低架滤波器，对应编码：2；
peak--峰值滤波器，对应编码：3；
A_f0,A_Q,A_boost,A_gain，分别对应滤波器的中心频率-f0，品质因数-Q值，抬升增益-boost，和增益-gain对应的系数向量，各3个元素，即：
A＝[a0,a1,a2]T。
步骤(1)中所述滤波器参数包括基准参数拓扑类型-type，和中心频率-f0、增益-gain、抬升增益-boost、品质因数-Q参数。
步骤(1)中，所述从一种拓扑类型衍化到另一种拓扑类型的节点(即奇异点)，通过观察相邻扇形区域连接点多段滤波器的波形变化，如果确认某一段波形的拓扑结构已经从一种形式转换到另一种形式，则可确定在某个扇形区域内是否存在奇异点，再通过在该扇形区域内等间隔采样的办法，进一步缩小和估计出奇异点可能所在的范围区间。
步骤(2)中所述多项式拟合优选为二阶多项式拟合；在实际运用过程中，多项式拟合的阶数，根据需要平滑扇形区域的大小设定；当平滑区域较大，参数动态变化区间大，需要更多的阶数拟合，当平滑区域较小，参数动态变化范围较小，可以用较少的阶数拟合。本专利中根据实验结果，采用二阶多项式拟合。
其中，本发明通过采用相对较大离散角度间隔的工程测量方法与小角度范围插值平滑过渡相结合的算法；既保证了真实性，也实现了平滑过渡；可以实现更小粒度的平滑。
步骤(2)中所述插值算法可以对角度变化实现连续的平滑过渡；比申请号为201611243763.4专利(具有头部跟踪的双耳头戴式耳机呈现)中提到的固定值法，有更好的角度分辨率；实际应用中可根据需要选用平滑的粒度大小(即角度间隔)。
上述声场重现中仿真系数向量信息的获取算法在声场重现头戴设备中应用。
现有技术的算法方案中采用限定拓扑结构类型和数量，固定中心频率，只更新增益和延时2个参数的方法，算法实现简单，容易实现滤波器的平滑切换，计算量小，在计算资源有限的情况下，工程上也容易实现。其缺点是这一切换算法和策略与实际声音的传播特性相差较大，听感上与实际的声场特性会有一定差别；主要原因是声源传播到人耳的传播特性，其频响曲线是较复杂的，要用10几段级联的滤波器才能逼近。y为准确描述滤波器的波形特征，其相关参数，除前文提到的中心频率、增益外，通常还包括滤波器的拓扑类型、抬升增益、品质因数等多个参数。这些参数在声源发生角度或位置变化时，通常都会跟随变化，特别是大角度变化时，这些参数变化会很剧烈，呈现非线性变化特点，甚至出现奇异点。如滤波器拓扑类型，由一种类型衍化到另一种类型。这些非线性变化特点，是影响听感的重要因素。所以现有的获取声场重现的仿真系数向量信息的算法，在重现声场音效的真实感受上是有局限的。
本发明所达到的有益效果在于：
1.本发明的算法应用的声场范围更大，可以达到-105度～+105度，采用大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法，有更大的声场动态范围和高仿真效果。
2.本发明采用大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法。既保证了真实性，也实现了平滑过渡；基于最小二乘-多项式拟合的系数估计算法，可以实现更小粒度的平滑。比申请号为201611243763.4的专利(具有头部跟踪的双耳头戴式耳机呈现)中提到的以索引表的形式预置在内存中的固定值法，有更好的角度分辨率。
具体实施方式
下面将结合本发明实施例对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明，若本发明实施例中有涉及方向性指示(诸如上、下、左、右、前、后、顶、底……)，则该方向性指示仅用于解释在某一特定姿态下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
在本申请中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系，除非另有明确的限定。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本申请中的具体含义。
需要说明的是，当元件被称为“固定于”或“设置于”另一个元件，它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件，它可以是直接连接到另一个元件或者可能同时存在居中元件。
另外，若本发明实施例中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
目前，在声场重现技术领域，现有技术的切换算法和策略与实际声音的传播特性相差较大，听感上与实际的声场特性会有一定差别；造成这种差别的主要缺陷包括以下：仿真系数参数在声源发生角度或位置变化时，通常都会跟随变化，特别是大角度变化时，这些参数变化会很剧烈，呈现非线性变化特点，在实际中会出现若干奇异点从而导致“咔哒”声。为了解决上述技术问题，本发明提出了一种声场重现中仿真系数向量信息的获取算法。
实施例一
基于HRTF的测量方法，其滤波器相关参数随角度的变化，在较大范围时呈非线性特性，有较大的动态变换范围；在较小范围内，呈近似连续的线性变化特性。本发明基于上述认知，采用相对大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法。
本发明提供一种声场重现中仿真系数向量信息的获取算法及其应用；其中所述声场重现中仿真系数向量信息的获取算法，包括以下步骤：
一种声场重现中仿真系数向量信息的获取算法，采用大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法，具体包括以下步骤：
(1)大角度范围的离散工程测量，采用HRTF(头部相关传递函数)的测量方法，在-105～+105度范围内，从初始零位开始，以15度为测量步长，测得测量步长节点的滤波器参数再由数字信号处理软件(Digital Signal Processing，简写为DSP软件)拟合出参数型级联的多段滤波器，并将所述滤波器参数以测量的角度为单位保存，供出现奇异点(即从一种拓扑类型衍化到另一种拓扑类型的节点)角度切换时调用，作为小角度范围插值平滑过渡的起算数据；
(2)小角度范围插值平滑过渡，采用插值算法，在每个15度测量步长的扇形区间内，再等间隔采样，采用基于最小二乘法-多项式拟合的逼近算法，估计出多项式系数(a0，a1，a2)；通过系数拟合估计滤波器参数(y)，可以对角度变化实现连续的平滑过渡；当遇到拓扑结构的奇异点时，则重新起算多项式系数(a0，a1，a2)，然后通过系数拟合估计滤波器参数；结合对应扇形区域、滤波器的频段号、滤波器拓扑类型信息获得多项式拟合出的仿真系数向量组。
步骤(2)中所述插值算法，具体步骤如下：在-105～+105度范围内，对每个15度的扇形区间，从初始零位开始区分左右编成L_0～L_6，R_0～R_6索引号，以每个扇形区间的每级滤波器的拓扑类型-type为基准参数，对中心频率-f0、增益-gain、抬升增益-boost、品质因数-Q参数，分别进行二阶多项式拟合输出相应的系数(a0，a1，a2)；当遇到拓扑结构的奇异点时，则重新起算多项式系数(a0，a1，a2)。
所述的二阶多项式拟合，具体如下：
y＝a0+a1·x+a2·x2
其中：
x：为测得的角度值；
y：为需要获得的滤波器参数，为f0、gain、boost、Q中的任意一个。
其中a0，a1，a2为多项式系数；
A＝[a0,a1,a2]T，表示需要估计出的多项式系数向量；其中T为滤波器拓扑类型。
考虑多参数下，步骤(2)中所述拟合出的仿真系数向量组包括以下7个维度：
coeff＝[bin_num,band_num,fil_T,A_f0,A_Q,A_boost,A_gain]
其中：
bin_num，表示扇形区域的编号，即；L_0～L_6，R_0～R_6，共14个元素；
band_num，表示多段滤波器的频段号；选用最多级数，作为参考编号；本发明中为band1～band12，即12个元素；
fil_T,表示滤波器拓扑类型；本发明中包括3个滤波器拓扑类型元素：
highshelf--高架滤波器，对应编码：1；
lowshelf--低架滤波器，对应编码：2；
peak--峰值滤波器，对应编码：3；
A_f0,A_Q,A_boost,A_gain，分别对应滤波器的中心频率-f0，品质因数-Q值，抬升增益-boost，和增益-gain对应的系数向量，各3个元素，即：
A＝[a0,a1,a2]T。
步骤(1)中所述滤波器参数包括基准参数拓扑类型-type，和中心频率-f0、增益-gain、抬升增益-boost、品质因数-Q参数。
上述步骤(2)中所述拟合出的仿真系数向量组为“7维系数向量”。
步骤(1)中，所述奇异点，通过观察相邻扇形区域连接点多段滤波器的波形变化，如果确认某一段波形的拓扑结构已经从一种形式转换到另一种形式，则可确定在某个扇形区域内是否存在奇异点，再通过在该扇形区域内等间隔采样的办法，进一步缩小和估计出奇异点可能所在的范围区间。
步骤(2)中所述多项式拟合优选为二阶多项式拟合；在实际运用过程中，多项式拟合的阶数，根据需要平滑扇形区域的大小设定；当平滑区域较大，参数动态变化区间大，需要更多的阶数拟合，当平滑区域较小，参数动态变化范围较小，可以用较少的阶数拟合。本专利中根据实验结果，采用二阶多项式拟合。
其中，本发明通过采用相对较大离散角度间隔的工程测量方法与小角度范围插值平滑过渡相结合的算法；既保证了真实性，也实现了平滑过渡；可以实现更小粒度的平滑。
步骤(2)中所述插值算法可以对角度变化实现连续的平滑过渡；比申请号为201611243763.4专利(具有头部跟踪的双耳头戴式耳机呈现)中提到的固定值法，有更好的角度分辨率；实际应用中可根据需要选用平滑的粒度大小(角度间隔)。
上述声场重现中仿真系数向量信息的获取算法在声场重现头戴设备中应用。
现有技术算法方案中采用限定拓扑结构类型和数量，固定中心频率，只更新增益和延时2个参数的方法，算法实现简单，容易实现滤波器的平滑切换，计算量小，在计算资源有限的情况下，工程上也容易实现。其缺点是这一切换算法和策略与实际声音的传播特性相差较大，听感上与实际的声场特性会有一定差别；主要原因是声源传播到人耳(或仿真人头，实际测量HRTF时多用仿真耳)的传播特性，其频响曲线是较复杂的，要用10几段级联的滤波器才能逼近。y为准确描述滤波器的波形特征，其相关参数，除前文提到的中心频率、增益外，通常还包括滤波器的拓扑类型、抬升增益、品质因数等多个参数。这些参数在声源发生角度或位置变化时，通常都会跟随变化，特别是大角度变化时，这些参数变化会很剧烈，呈现非线性变化特点，甚至出现奇异点。如滤波器拓扑类型，由一种类型衍化到另一种类型。这些非线性变化特点，是影响听感的重要因素。所以现有的获取声场重现的仿真系数向量信息的算法，在重现声场音效的真实感受上是有局限的。
本发明所达到的有益效果在于：
1.本发明的算法应用的声场范围更大，可以达到-105度～+105度，采用大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法，有更大的声场动态范围和高仿真效果。
2.本发明采用大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法。既保证了真实性，也实现了平滑过渡；基于最小二乘-多项式拟合的系数估计算法，可以实现更小粒度的平滑。比申请号为201611243763.4的专利(具有头部跟踪的双耳头戴式耳机呈现)中提到的以索引表的形式预置在内存中的固定值法，有更好的角度分辨率。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种声场重现中仿真系数向量信息的获取算法，其特征在于：采用大角度范围的离散工程测量与小角度范围插值平滑过渡相结合的算法，具体包括以下步骤：
(1)大角度范围的离散工程测量，采用头部相关传递函数的测量方法，在-105～+105度范围内，从初始零位开始，以15度为测量步长，测得测量步长节点的滤波器参数再由数字信号处理软件拟合出参数型级联的多段滤波器，并将所述滤波器参数以测量的角度为单位保存，供出现从一种拓扑类型衍化到另一种拓扑类型的节点时角度切换时调用，作为小角度范围插值平滑过渡的起算数据；
(2)小角度范围插值平滑过渡，采用插值算法，在每个15度测量步长的扇形区间内，再等间隔采样，采用基于最小二乘法-多项式拟合的逼近算法，估计出多项式系数；通过系数拟合估计滤波器参数，可以对角度变化实现连续的平滑过渡；当遇到拓扑结构的奇异点时，则重新起算多项式系数，然后通过系数拟合估计滤波器参数；结合对应扇形区域、滤波器的频段号、滤波器拓扑类型信息获得多项式拟合出的仿真系数向量组；
步骤(2)中所述插值算法，具体步骤如下：在-105～+105度范围内，对每个15度测量步长的扇形区间，从初始零位开始区分左右编成L_0～L_6，R_0～R_6索引号，以每个扇形区间的每级滤波器的拓扑类型-type为基准参数，对中心频率-f0、增益-gain、抬升增益-boost、品质因数-Q参数分别进行二阶多项式拟合输出相应的系数；当遇到拓扑结构的奇异点时，则重新起算多项式系数；
所述的二阶多项式拟合，具体如下：
y＝a0+a1·x+a2·x2
其中：
x：为测得的角度值；
y：为需要获得的滤波器参数，为f0、gain、boost、Q中的任意一个；
其中a0，a1，a2为多项式系数；
A＝[a0,a1,a2]T，表示需要估计出的多项式系数向量；其中T为滤波器拓扑类型；
步骤(2)中所述拟合出的仿真系数向量组包括以下7个维度：
coeff＝[bin_num,band_num,fil_T,A_f0,A_Q,A_boost,A_gain]
其中：
bin_num，表示扇形区域的编号，即；L_0～L_6，R_0～R_6，共14个元素；
band_num，表示多段滤波器的频段号；选用最多级数，作为参考编号；为band1～band12，即12个元素；
fil_T,表示滤波器拓扑类型；包括3个滤波器拓扑类型元素：
highshelf--高架滤波器，对应编码：1；
lowshelf--低架滤波器，对应编码：2；
peak--峰值滤波器，对应编码：3；
A_f0,A_Q,A_boost,A_gain，分别对应滤波器的中心频率-f0，品质因数-Q值，抬升增益-boost，和增益-gain对应的系数向量，各3个元素，即：
A＝[a0,a1,a2]T；
步骤(1)中所述滤波器参数包括基准参数拓扑类型-type，和中心频率-f0、增益-gain、抬升增益-boost、品质因数-Q参数；
步骤(2)中所述拟合出的仿真系数向量组为7维系数向量。
2.根据权利要求1所述的声场重现中仿真系数向量信息的获取算法，其特征在于：步骤(1)中，所述从一种拓扑类型衍化到另一种拓扑类型的节点，通过观察相邻扇形区域连接点多段滤波器的波形变化，如果确认某一段波形的拓扑结构已经从一种形式转换到另一种形式，则可确定在某个扇形区域内是否存在奇异点，再通过在该扇形区域内等间隔采样的办法，进一步缩小和估计出奇异点可能所在的范围区间。
3.权利要求1～2任一项所述声场重现中仿真系数向量信息的获取算法在声场重现头戴设备中应用。
CN202110226269.1A 2021-03-01 2021-03-01 声场重现中仿真系数向量信息的获取算法及其应用 Active CN113068112B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110226269.1A CN113068112B (zh) 2021-03-01 2021-03-01 声场重现中仿真系数向量信息的获取算法及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110226269.1A CN113068112B (zh) 2021-03-01 2021-03-01 声场重现中仿真系数向量信息的获取算法及其应用

Publications (2)

Publication Number Publication Date
CN113068112A CN113068112A (zh) 2021-07-02
CN113068112B true CN113068112B (zh) 2022-10-14

Family

ID=76559455

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110226269.1A Active CN113068112B (zh) 2021-03-01 2021-03-01 声场重现中仿真系数向量信息的获取算法及其应用

Country Status (1)

Country Link
CN (1) CN113068112B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105872940A (zh) * 2016-06-08 2016-08-17 北京时代拓灵科技有限公司 一种虚拟现实声场生成方法及系统

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7085393B1 (en) * 1998-11-13 2006-08-01 Agere Systems Inc. Method and apparatus for regularizing measured HRTF for smooth 3D digital audio
US8229143B2 (en) * 2007-05-07 2012-07-24 Sunil Bharitkar Stereo expansion with binaural modeling
US9491560B2 (en) * 2010-07-20 2016-11-08 Analog Devices, Inc. System and method for improving headphone spatial impression
WO2015134658A1 (en) * 2014-03-06 2015-09-11 Dolby Laboratories Licensing Corporation Structural modeling of the head related impulse response
WO2016165776A1 (en) * 2015-04-17 2016-10-20 Huawei Technologies Co., Ltd. Apparatus and method for driving an array of loudspeakers with drive signals
US9860666B2 (en) * 2015-06-18 2018-01-02 Nokia Technologies Oy Binaural audio reproduction
US9918177B2 (en) * 2015-12-29 2018-03-13 Harman International Industries, Incorporated Binaural headphone rendering with head tracking
CN110036655B (zh) * 2016-12-12 2022-05-24 索尼公司 Hrtf测量方法、hrtf测量装置和存储介质

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105872940A (zh) * 2016-06-08 2016-08-17 北京时代拓灵科技有限公司 一种虚拟现实声场生成方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
虚拟制造系统中立体声音的合成方法研究;张卫波等;《现代制造工程》;20050618(第06期);第40-42页 *

Also Published As

Publication number Publication date
CN113068112A (zh) 2021-07-02

Similar Documents

Publication Publication Date Title
US7715568B2 (en) Binaural sound reproduction apparatus and method, and recording medium
CN107018460B (zh) 具有头部跟踪的双耳头戴式耳机呈现
KR100636213B1 (ko) 실시간 주파수 특성 보정 방법 및 그를 적용한 사운드재생 장치
JP6790654B2 (ja) フィルタ生成装置、フィルタ生成方法、及びプログラム
US20210250686A1 (en) Method and apparatus for processing an audio signal based on equalization filter
CN110024419A (zh) 用于不对称听觉传输音频再现的增益相位均衡(gpeq)滤波器和调谐方法
CN113099359B (zh) 一种基于hrtf技术的高仿真声场重现的方法及其应用
CN113068112B (zh) 声场重现中仿真系数向量信息的获取算法及其应用
CN106973355B (zh) 环绕立体声实现方法和装置
US20210297802A1 (en) Signal processing device, signal processing method, and program
CN115802274A (zh) 音频信号处理方法、电子设备和计算机可读存储介质
CN110313188A (zh) 头外定位处理装置、头外定位处理方法以及头外定位处理程序
JP2010011083A (ja) バイノーラル収音再生システム
JP2755081B2 (ja) 音像定位制御方法
US20230040821A1 (en) Processing device and processing method
US20210377684A1 (en) Processing device, processing method, reproducing method, and program
US20060245305A1 (en) System comprising sound reproduction means and ear microphones
CN113905323B (zh) 适用于服务型机器人播放音频时的感知声源高度修正方法
EP3855765B1 (en) Processing device, processing method, reproduction method, and program
DK180449B1 (en) A method and system for real-time implementation of head-related transfer functions
US20240080618A1 (en) Out-of-head localization processing device, out-of-head localization processing method, and computer-readable medium
JP2023047706A (ja) フィルタ生成装置、及びフィルタ生成方法
JP2023047707A (ja) フィルタ生成装置、及びフィルタ生成方法
Lazzarini et al. New csound opcodes for binaural processing
JP2023024040A (ja) 処理装置、及び処理方法

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
GR01 Patent grant
GR01 Patent grant