CN105324762B - 滤波器系数群计算设备以及滤波器系数群计算方法 - Google Patents
滤波器系数群计算设备以及滤波器系数群计算方法 Download PDFInfo
- Publication number
- CN105324762B CN105324762B CN201480035761.9A CN201480035761A CN105324762B CN 105324762 B CN105324762 B CN 105324762B CN 201480035761 A CN201480035761 A CN 201480035761A CN 105324762 B CN105324762 B CN 105324762B
- Authority
- CN
- China
- Prior art keywords
- fourier transform
- coefficient group
- filter coefficient
- frequency
- filter
- 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
Links
- 238000000205 computational method Methods 0.000 title claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 6
- 238000000034 method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000003801 milling Methods 0.000 claims description 2
- 230000006870 function Effects 0.000 description 25
- 238000001228 spectrum Methods 0.000 description 23
- 238000005070 sampling Methods 0.000 description 22
- 230000005540 biological transmission Effects 0.000 description 6
- 230000005236 sound signal Effects 0.000 description 6
- 238000000605 extraction Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03G—CONTROL OF AMPLIFICATION
- H03G5/00—Tone control or bandwidth control in amplifiers
- H03G5/02—Manually-operated control
- H03G5/025—Equalizers; Volume or gain control in limited frequency bands
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03G—CONTROL OF AMPLIFICATION
- H03G5/00—Tone control or bandwidth control in amplifiers
- H03G5/02—Manually-operated control
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0211—Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
- H03H17/0213—Frequency domain filters using Fourier transforms
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Circuit For Audible Band Transducer (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Tone Control, Compression And Expansion, Limiting Amplitude (AREA)
- Complex Calculations (AREA)
Abstract
滤波器系数群计算设备配置为具有:对通过输入装置输入的频率特性执行傅里叶逆变换的装置;对通过傅里叶逆变换获得的数字序列执行短时傅里叶变换的装置;将频率越高窗口长度越短的窗函数应用至通过短时傅里叶变换获得的频域信号的装置;对在应用窗函数之后的频域信号应用短时傅里叶逆变换的装置;对通过短时傅里叶变换获得的数字序列进行重叠相加的装置;以及将在重叠相加之后的数字序列确定为形成具有通过输入装置输入的频率特性的滤波器的滤波器系数群的装置。
Description
技术领域
本发明涉及一种滤波器系数群计算设备以及一种用于计算滤波器的系数的滤波器系数群计算方法。
背景技术
用于改变输入的音频信号的频率特性的其中一种装置为图形均衡器。一般地,在图形均衡器中,多个调节点需要逐个进行调节。因此,存在的问题在于当有很多调节点时,用户的操作负担很大。
出于该原因,日本专利临时公开第H11-112254A号(在本文中为“专利文件1”)和第2004-62503A号(在本文中为“专利文件2”)提出了能够克服以上问题的装置。在专利文件1中和专利文件2中公开的装置检测用户利用手写笔或手指在均衡器的屏幕上描绘的路径(在本文中为“触摸操作路径”),并且基于所检测的触摸操作路径设定用于每个调节点的增益(每个频带的中心频率)。用户能够通过跟踪均衡器的屏幕的一个操作来设定所有调节点的增益,而不需要分别地和独立地逐个操作多个调节点。然而,在图形均衡器中,由于每个频带的中心频率和Q值(调节频带宽度)是固定的,所以调节灵活性较低。因此,利用在专利文件1或专利文件2中公开的图形均衡器难以准确地设定用户想要的滤波器的频率特性(在本文中为“滤波器特性”)。
与图形均衡器相似的一个装置是参数均衡器。在参数均衡器中,除了用于每个频带的中心频率的增益之外,每个频带的中心频率和Q值可以调节。日本专利临时公开第H5-175773A号(在本文中为“专利文件3”)公开了参数均衡器的示例。在专利文件3中公开的参数均衡器提前准备了Q值与用于每个频带的中心频率的增益绝对值成比例地改变的多个滤波器系数群,由用户从提前准备的多个滤波器系数群中读出适于增益调节操作的滤波器系数群,并且使用读出的滤波器系数群将滤波器特性调节为总体上平缓的和连续的。可以考虑到,用户想要的滤波器特性可以通过应用在专利文件3中公开的技术来准确地设定。然而,在参数均衡器中,由于纹波保持处于滤波器特性,因此存在的问题在于当音乐等在回放时会检测到不必要的混响成分。
发明内容
本申请是考虑到以上情况而发明的,并且本发明的目的是提供一种滤波器系数群计算设备以及一种适用于计算滤波器系数群的滤波器系数群计算方法,所述滤波器系数群用于获得用户的输入被精确的反映并且混响成分被消除的滤波器特性。
本发明的实施方案的滤波器系数群计算设备设置有用户在其上输入频率特性的输入装置,并且配置为计算形成具有通过输入装置输入的频率特性的滤波器的滤波器系数群,所述滤波器系数群计算设备包括:傅里叶逆变换装置,其配置为对通过输入装置输入的频率特性执行傅里叶逆变换;短时傅里叶变换装置,其配置为对通过傅里叶逆变换获得的数字串执行短时傅里叶变换;开窗装置,其配置为利用窗口长度随着频率的增加而缩短的函数对通过短时傅里叶变换获得的频域信号执行开窗;短时傅里叶逆变换装置,其配置为对在开窗之后的频域信号执行短时傅里叶逆变换;重叠相加装置,其配置为对通过短时傅里叶逆变换获得的数字串执行重叠相加;以及滤波器系数群确定装置,其配置为将在重叠相加之后的数字串确定为形成具有通过所述输入装置输入的频率特性的滤波器的滤波器系数群。
在本实施方案中,通过利用窗口长度随着频率的增加而缩短的函数执行开窗并且之后执行短时傅里叶逆变换和重叠相加而获得了滤波器系数群,所述滤波器系数群用于获得用户的输入被精确的反映并且混响成分被消除的滤波器特性。
例如,所述函数定义窗口长度随着频率的增加而指数缩短的关系。
所述滤波器系数群计算设备可以包括减少装置,其配置为减少包括在通过所述滤波器系数群确定装置确定的滤波器系数群中的系数的数量。
所述滤波器系数群计算设备可以包括最小相位转换装置,其配置为对在通过所述傅里叶逆变换装置进行的傅里叶逆变换之后的数字串执行最小相位转换。在这种情况下,所述短时傅里叶变换装置对在最小相位转换之后的数字串执行短时傅里叶变换。
例如,所述输入装置为在其上能够进行触摸操作的触摸面板。
此外,本发明的实施方案的滤波器系数群计算方法是这样一种方法,其用于计算具有由用户输入的频率特性的滤波器的滤波器系数群,所述滤波器系数群计算方法包括:傅里叶逆变换步骤,对由用户输入的频率特性执行傅里叶逆变换;短时傅里叶变换步骤,对通过傅里叶逆变换获得的数字串执行短时傅里叶变换;开窗步骤,利用窗口长度随着频率的增加而缩短的函数对通过短时傅里叶变换获得的频域信号执行开窗;短时傅里叶逆变换步骤,对在开窗之后的频域信号执行短时傅里叶逆变换;重叠相加步骤,对通过短时傅里叶逆变换获得的数字串执行重叠相加;以及滤波器系数群确定步骤,将在重叠相加之后的数字串确定为形成具有由用户输入的频率特性的滤波器的滤波器系数群。
在本实施方案中,通过利用窗口长度随着频率的增加而缩短的函数执行开窗并且之后执行短时傅里叶逆变换和重叠相加而获得了滤波器系数群,所述滤波器系数群用于获得用户的输入被精确的反映并且混响成分被消除的滤波器特性。
附图说明
图1为示出本发明的实施方案的声音处理装置的配置的框图。
图2示出了在设置到本发明的实施方案的声音处理装置的输入接口单元上显示的屏幕的示例。
图3为示出通过设置到本发明的实施方案的声音处理装置的滤波器系数群计算单元所执行的滤波器系数群计算过程的流程图。
图4示出了根据从设置到本发明的实施方案的滤波器系数群计算单元的坐标输入单元输入的取样坐标点所计算的功率谱。
图5示出了在图3中的处理步骤S16中的最小相位转换之后的滤波器系数群(最小相位转换)。
图6示出了通过在图3中的处理步骤S17中的STFT获得的频谱图(短时傅立叶变换)。
图7示出了在图3中的处理步骤S18中使用的对图形所用的函数(开窗)。
图8示出了在图3中的处理步骤S18中的开窗之后的频谱图(开窗)。
图9示出了在图5中显示的滤波器系数群(虚线)以及在图3中的步骤S20中的重叠相加之后的滤波器系数群(重叠相加)(实线)。
图10示出了在图3中的处理步骤S20中的重叠相加之后的滤波器系数群(重叠相加)给到FIR滤波器单元时FIR滤波器单元的功率谱。
具体实施方式
在下文中,将参考附图来描述本发明的实施方案。应当注意的是,在下文中,声音处理装置表示为本发明的一个实施方案并进行解释。
[声音处理装置1的总体配置]
图1为本实施方案的声音处理装置1的配置的框图。如图1所示,声音处理装置1包括:输入接口单元12、滤波器系数群计算单元14和FIR(有限冲激响应)滤波器单元16。
在FIR滤波器中,以有限时间长度表示的冲激响应本身是滤波器的系数。所以,滤波器系数群的确定基本上与FIR滤波器的滤波器特性的设定是同义的。因此,滤波器系数群计算单元14根据由用户对输入接口单元12的操作来计算滤波器系数群。滤波器系数群计算单元14通过给出所计算的滤波器系数群来设定FIR滤波器单元16的滤波器特性。
对于FIR滤波器单元16,通过将以不可逆的或可逆的压缩格式的编码信号解码所产生的音频信号从声源输入。在本实施方案中,音频信号取样频率例如为44.1kHz。FIR滤波器单元16利用由用户通过输入接口单元12和滤波器系数群计算单元14设定的特性来滤波输入的音频信号,以修改音频信号的频率特性。经频率特性修改之后的音频信号经由在图中未示出的功率放大器通过扬声器输出。这使得用户能够听到反映对输入接口单元12的操作结果的音乐等(音乐等的音质得到修正或声像得到改进)。
[输入接口单元12的配置]
输入接口单元12为接收并处理由用户使用手写笔或手指进行的触摸操作的触摸面板。图2示出了在输入接口单元12上显示的屏幕的示例。如图2所示,在输入接口单元12的显示屏幕中显示了允许用户设定FIR滤波器单元16的滤波器特性的滤波器特性设定屏幕12a。滤波器特性设定屏幕12a为在其上画出根据用户的操作的功率谱的屏幕,其中纵轴为功率(单位:dB),以及横轴为频率(单位:Hz)。功率是频率与自身相乘的乘积。此外,人耳听觉特性对于频率是呈对数的。在横轴上的频率利用对数显示以适应于人耳听觉特性。
在用户的手指等触摸滤波器特性设定屏幕12a时,输入接口单元12在触摸部分处显示图标12b。当用户利用手指描绘滤波器特性设定屏幕12a时,输入接口单元12检测触摸操作路径,并且画出检测的触摸操作路径,同时在用户的手指触摸的部分处显示图标12b。如在图2中示例的,根据用户的操作的功率谱由此显示在滤波器特性设定屏幕12a上。如上所述,用户能够通过在滤波器特性设定屏幕12a上的触摸操作来输入想要的滤波器特性(在以上的示例中为功率谱)。本实施方案通过采用触摸面板装置作为输入接口单元12而能够实现直观的滤波器设计。
输入接口单元12不限于触摸面板装置,而可以是通过比如鼠标的指向装置来接收和处理操作的装置。在这种情况下,输入接口单元12画出在滤波器特性设定屏幕12a内作出的拖动操作路径作为对于FIR滤波器单元16设定的滤波器特性(功率谱)。
[滤波器系数群计算单元14的结构和滤波器系数群计算流程]
如图1所示,滤波器系数群计算单元14包括坐标计算单元14a、幅度谱计算单元14b、传递函数计算单元14c、傅里叶逆变换(IFFT:快速傅里叶逆变换)单元14d、实部提取单元14e、最小相位转换单元14f、短时傅里叶变换(STFT:短时傅里叶变换)单元14g、开窗单元14h、短时傅里叶逆变换(ISTFT:短时傅里叶逆变换)单元14i、重叠相加单元14j和抽头数量减少单元14k。滤波器系数群计算单元14计算滤波器系数群并将所计算的滤波器系数群给到FIR滤波器单元16。图3示出了通过滤波器系数群计算单元14执行的滤波器系数群计算过程的流程图。
[图3中的S11(坐标点取样)]
坐标输入单元14a对用户在滤波器特性设定屏幕12a上的触摸操作路径进行取样,并且将在滤波器特性设定屏幕12a内的取样的坐标点(在下文中称为“取样坐标点”)输出至幅度谱计算单元14b。在本实施方案中,取样坐标点的数量是8192个,其与FIR滤波器的坐标的数量(抽头数量)相同。应当注意,以上提到的取样坐标点的数量是一个示例。一般地,相对于用户的操作的FIR滤波器单元16的设定精确性随着取样坐标点的数量(抽头数量)的增加而提高。
[图3中的S12(幅度谱计算)]
幅度谱计算单元14b根据用户从坐标输入单元14a输入的取样坐标点的操作来计算功率谱,并根据所计算的功率谱来计算幅度谱|F(ω)|。图4示出了根据从坐标输入单元14a输入的取样坐标点计算的功率谱的示例。图4(a)示出了在全频率范围内的功率谱,而图4(b)示出了在全频率范围的一部分内并且扩大了的功率谱。如从这些图中可见的,在该阶段的功率谱不具有平滑和连续的特性但具有纹波在整个频率范围内保持的特性。纹波是由于在输入接口单元12上的触摸操作期间手抖动、屏幕的分辨率等所引起的。如上所述,滤波器特性设定屏幕12a的横轴(频率)利用对数示出以适应于人耳听觉特性。因此,如果在对应于滤波器特性设定屏幕12a的高频部分的频率范围内将频率轴考虑为线性标度,则微小的手抖动等引起大的纹波。
[图3中的S13(传递函数计算)]
传递函数计算单元14c根据由幅度谱计算单元14b计算的幅度谱|F(ω)|来计算传递函数F(ω)。即,传递函数计算单元14c将对应于根据用户的操作的功率谱设定为传递函数。
[图3中的S14(傅里叶逆变换)]
IFFT单元14d利用IFFT将由传递函数计算单元14c计算的传递函数F(ω)从频域取样序列变换为时域取样序列。
[图3中的S15(实部提取)]
实部提取单元14e从作为IFFT单元14b进行的IFFT结果而获得的取样序列中提取实部。此处提取的实部是给FIR滤波器单元16的对应于滤波器系数群的取样序取。
[图3中的S16(最小相位转换)]
最小相位转换单元14f对由实部提取单元14e提取的实部(滤波器系数群)执行最小相位转换,以有效地消除混响成分。图5示出了在最小相位转换后的滤波器系数群的示例(冲激响应)。在图5中,纵轴是功率(单位:dB),而横轴是取样。应当注意的是,在本实施方案中单位是取样的轴可以被替换为时间轴。
如图5所示,在最小相位转换后的滤波器系数群包括许多产生混响成分的若干系数。在本实施方案修改的示例中,线性相位转换可以代替最小相转换而执行。
[图3中的S17(短时傅里叶变换)]
STFT单元14g利用STFT将在通过最小相位转换单元14f的最小相位转换后的系数群变换为频域。下列为用于STFT的条件。
FFT长度(单位:取样):256
重叠长度(单位:取样):240
窗函数:汉明
图6示出通过STFT单元14g进行的STFT获得的频谱图的示例。在图6中,纵轴是频率(单位:Hz),横轴是取样,并且颜色的明暗指示功率(单位:dB)。如在图6中的示例所示,在该阶段的频谱图具有在从低范围到高范围上分布的混响成分。
[图3中的S18(开窗)]
开窗单元14h利用预义的函数对通过STFT单元14g进行的STFT获得的频率成分执行开窗。图7了示出在开窗过程中使用的对图形所用的函数。在图7中,纵轴是窗口长度(单位:取样),而横轴是频率(单位:Hz)。如在图7中的示例所示,该函数定义了在频率和窗口长度之间的关系。在下文中,该函数将被称为“频率窗口长度函数”。频率窗口长度函数定义了窗口长度随着频率的增加而对数缩短的关系。此外,由于如果窗口长度太短,则数据因过度的数据缩减而不能恢复,所以在特定频率以上的范围中的频率窗口长度函数的窗口长度被限定为特定值(14个取样是最小的窗口长度)。图8示出了在开窗之后的频谱图的示例。如在图8中所示,通过利用频率窗口长度函数执行开窗,混响成分被消除,同时特别地在低范围处保留了功率。
[图3中的S19(短时傅里叶逆变换)]
ISTFT单元14i利用ISTFT将在通过开窗单元14h进行的加窗之后的频率成分变换为时域取样序列。
[图3中的S20(重叠相加)]
重叠相加单元14j对作为通过ISTFT单元14i进行的IFFT的结果而获得的取样序列(滤波器系数群)执行重叠相加,以消除时域信号的不连续性。在本实施方案中,重叠长度是240个取样。图9示出了在图5中显示的滤波器系数群(虚线)以及在通过重叠相加单元14j进行的重叠相加之后的滤波器系数群的示例(实线)。此外,分别在图10(a)和图10(b)中示出了与在图4(a)和图4(b)中示出的那些相似的图。在图10中,在图4中示出的功率谱以虚线显示,而在重叠相加之后的滤波器系数群给到FIR滤波器单元16时,FIR滤波器单元16的功率谱以实线显示。
在本实施方案中,如从利用频率窗口长度函数的开窗之前的滤波器系数群(图9中的虚线)与在利用频率窗口长度函数的开窗之后的滤波器系数群(图9中的实线)之间的比较可以看出,通过利用窗口长度随着频率的增加而缩短的函数开窗,然后执行短时傅里叶逆变换和重叠相加,混响成分(对应于混响成分的滤波器系数)大幅度地衰减。此外,如在图10中的示例所示,在重叠相加之后的功率谱(实线)具有平滑且连续的特性,同时精确地反映用户的触摸操作路径,并且在输入接口单元12上的触摸操作期间的手抖动、屏幕的分辨率等引起的纹波在整个频率范围内都被消除。
如上所述,在本实施方案中,计算滤波器特性精确地反映用户想要的滤波器特性(通过在滤波器特性设定屏幕12a上的触摸操作输入的滤波器特性)的滤波器系数群。此外,所计算的滤波器系数群是消除了由在操作期间的手抖动、屏幕的分辨率等、以及不必要的混响成分所引起的纹波的群。通过将这样的滤波器系数群给到FIR滤波器单元16,可以获得用户的输入被精确地反映并且混响成分被消除的滤波器特性。
[图3中的S21(抽头数量减少)]
抽头数量减少单元14k减少包括在通过重叠相加单元14j进行的重叠相加之后的滤波器系统群中的系数的数量(抽头数量)。在本实施方案中,滤波器系数的绝对值大于0处的有限部分被提取,并且朝着正无穷大执行取整操作(部分处理),使得所提取的有限部分的长度变为2的幂。抽头数量因此减少。在本实施方案中,提取了大约2200个取样(见图9),并且取整为4096个取样。即,抽头数量从8192减少为4096。
应当注意,有限部分不限于以上提到的部分,而是可以通过使用小到足以被忽略的值作为阈值。此外,当提取有限部分时,不仅阈值还有滤波器系数的包络的形状可以包括在用于做出决定的信息中。此外,在本实施方案中,由于前提是通过执行FFT来执行在频域中的乘积运算,因此执行取整运算而使得所提取的有限部分的长度变为2的幂。通过在时域执行卷积运算,取整运算可以省略。
以上为本发明的示例性实施方案的描述。本发明的实施方案不限于以上所解释的实施方案,并且在本发明的技术概念范围内的各种变化都是可以的。例如,在说明书中说明的示例性实施方案的适当组合和/或从说明书中显而易见的示例性实施方案也被包括在本发明的实施方案中。
Claims (9)
1.一种滤波器系数群计算设备,其设置有用户在其上输入频率特性的输入装置,并且配置为计算形成具有通过输入装置输入的频率特性的滤波器的滤波器系数群,其特征在于,所述滤波器系数群计算设备包括:
傅里叶逆变换装置,其配置为对通过输入装置输入的频率特性执行傅里叶逆变换;
短时傅里叶变换装置,其配置为对通过傅里叶逆变换获得的数字串执行短时傅里叶变换;
开窗装置,其配置为利用窗口长度随着频率的增加而缩短的函数对通过短时傅里叶变换获得的频域信号执行开窗;
短时傅里叶逆变换装置,其配置为对在开窗之后的频域信号执行短时傅里叶逆变换;
重叠相加装置,其配置为对通过短时傅里叶逆变换获得的数字串执行重叠相加;以及
滤波器系数群确定装置,其配置为将在重叠相加之后的数字串确定为形成具有通过所述输入装置输入的频率特性的滤波器的滤波器系数群。
2.根据权利要求1所述的滤波器系数群计算设备,
其中,函数定义了窗口长度随着频率的增加而指数缩短的关系。
3.根据权利要求1所述的滤波器系数群计算设备,
进一步包括减少装置,其配置为减少包括在通过所述滤波器系数群确定装置确定的滤波器系数群中的系数的数量。
4.根据权利要求1至3中任一项所述的滤波器系数群计算设备,
进一步包括最小相位转换装置,其配置为对在通过傅里叶逆变换装置进行的傅里叶逆变换之后的数字串执行最小相位转换,以及
其中,短时傅里叶变换装置对在通过最小相位转换装置进行的最小相位转换之后的数字串执行短时傅里叶变换。
5.根据权利要求1所述的滤波器系数群计算设备,
其中,输入装置为在其上能够进行触摸操作的触摸面板。
6.一种滤波器系数群计算方法,其用于计算具有由用户输入的频率特性的滤波器的滤波器系数群,其特征在于,所述滤波器系数群计算方法包括:
傅里叶逆变换步骤,对由用户输入的频率特性执行傅里叶逆变换;
短时傅里叶变换步骤,对通过傅里叶逆变换获得的数字串执行短时傅里叶变换;
开窗步骤,利用窗口长度随着频率的增加而缩短的函数,对通过短时傅里叶变换获得的频域信号执行开窗;
短时傅里叶逆变换步骤,对在开窗之后的频域信号执行短时傅里叶逆变换;
重叠相加步骤,对通过短时傅里叶逆变换获得的数字串执行重叠相加;以及
滤波器系数群确定步骤,将在重叠相加之后的数字串确定为形成具有由用户输入的频率特性的滤波器的滤波器系数群。
7.根据权利要求6所述的滤波器系数群计算方法,
其中,函数定义了窗口长度随着频率的增加而指数缩短的关系。
8.根据权利要求6所述的滤波器系数群计算方法,
进一步包括减少步骤,其减少包括在滤波器系数群确定步骤中确定的滤波器系数群中的系数的数量。
9.根据权利要求6至8中的任一项所述的滤波器系数群计算方法,
进一步包括最小相位转换步骤,其对在傅里叶逆变换步骤中的傅里叶逆变换之后的数字串执行最小相位转换,以及
其中,在短时傅里叶变换步骤中,对在通过最小相位转换步骤进行的最小相位转换之后的数字串执行短时傅里叶变换。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013132178A JP6216550B2 (ja) | 2013-06-25 | 2013-06-25 | フィルタ係数群演算装置及びフィルタ係数群演算方法 |
JP2013-132178 | 2013-06-25 | ||
PCT/JP2014/065197 WO2014208319A1 (ja) | 2013-06-25 | 2014-06-09 | フィルタ係数群演算装置及びフィルタ係数群演算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105324762A CN105324762A (zh) | 2016-02-10 |
CN105324762B true CN105324762B (zh) | 2017-11-28 |
Family
ID=52141665
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201480035761.9A Active CN105324762B (zh) | 2013-06-25 | 2014-06-09 | 滤波器系数群计算设备以及滤波器系数群计算方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9559658B2 (zh) |
EP (1) | EP3015996B1 (zh) |
JP (1) | JP6216550B2 (zh) |
CN (1) | CN105324762B (zh) |
WO (1) | WO2014208319A1 (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016050780A1 (en) * | 2014-10-02 | 2016-04-07 | Sony Corporation | Method, apparatus and system |
CN105931649A (zh) * | 2016-03-31 | 2016-09-07 | 欧仕达听力科技(厦门)有限公司 | 基于频谱分析的超低时延音频处理方法与系统 |
JP6661866B2 (ja) * | 2017-10-30 | 2020-03-11 | アンリツ株式会社 | フィルタ係数算出装置を備えた信号発生装置及び信号発生方法 |
CN111835495B (zh) * | 2020-09-16 | 2020-12-08 | 南昌大学 | 参考信号的检测方法、系统、可读存储介质及电子设备 |
CN113611292B (zh) * | 2021-08-06 | 2023-11-10 | 思必驰科技股份有限公司 | 用于语音分离、识别的短时傅里叶变化的优化方法及系统 |
CN115982527B (zh) * | 2023-03-21 | 2023-07-07 | 西安电子科技大学 | 一种基于fpga的时频域变换算法实现方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4896285A (en) * | 1987-03-23 | 1990-01-23 | Matsushita Electric Industrial Co., Ltd. | Calculation of filter factors for digital filter |
US5051942A (en) * | 1987-12-18 | 1991-09-24 | Matsushita Electric Industrial Co., Ltd. | Method and apparatus for calculating filter factors |
JP2005223887A (ja) * | 2004-01-06 | 2005-08-18 | Pioneer Electronic Corp | 音響特性調整装置 |
CN102231280A (zh) * | 2011-05-06 | 2011-11-02 | 山东大学 | 卷积语音信号的频域盲分离排序算法 |
CN102760435A (zh) * | 2012-07-03 | 2012-10-31 | 合肥工业大学 | 一种语音信号频域盲解卷积方法 |
CN102866010A (zh) * | 2012-09-28 | 2013-01-09 | 苏州大学 | 一种信号的谱峭度滤波方法及相关装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05175773A (ja) | 1991-12-20 | 1993-07-13 | Alpine Electron Inc | オーディオイコライザ |
JPH11112254A (ja) | 1997-10-01 | 1999-04-23 | Sony Corp | 遠隔制御装置ならびに制御装置 |
JP2004062503A (ja) | 2002-07-29 | 2004-02-26 | Sony Corp | 電子機器、オーディオ機器および機器操作処理方法 |
US7853342B2 (en) * | 2005-10-11 | 2010-12-14 | Ejamming, Inc. | Method and apparatus for remote real time collaborative acoustic performance and recording thereof |
WO2007086993A1 (en) * | 2006-01-30 | 2007-08-02 | Exxonmobil Upstream Research Company | Method for spatial filtering of electromagnetic survey data |
US9031834B2 (en) * | 2009-09-04 | 2015-05-12 | Nuance Communications, Inc. | Speech enhancement techniques on the power spectrum |
FR2961938B1 (fr) * | 2010-06-25 | 2013-03-01 | Inst Nat Rech Inf Automat | Synthetiseur numerique audio ameliore |
-
2013
- 2013-06-25 JP JP2013132178A patent/JP6216550B2/ja active Active
-
2014
- 2014-06-09 WO PCT/JP2014/065197 patent/WO2014208319A1/ja active Application Filing
- 2014-06-09 EP EP14818513.5A patent/EP3015996B1/en active Active
- 2014-06-09 US US14/898,882 patent/US9559658B2/en active Active
- 2014-06-09 CN CN201480035761.9A patent/CN105324762B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4896285A (en) * | 1987-03-23 | 1990-01-23 | Matsushita Electric Industrial Co., Ltd. | Calculation of filter factors for digital filter |
US5051942A (en) * | 1987-12-18 | 1991-09-24 | Matsushita Electric Industrial Co., Ltd. | Method and apparatus for calculating filter factors |
JP2005223887A (ja) * | 2004-01-06 | 2005-08-18 | Pioneer Electronic Corp | 音響特性調整装置 |
CN102231280A (zh) * | 2011-05-06 | 2011-11-02 | 山东大学 | 卷积语音信号的频域盲分离排序算法 |
CN102760435A (zh) * | 2012-07-03 | 2012-10-31 | 合肥工业大学 | 一种语音信号频域盲解卷积方法 |
CN102866010A (zh) * | 2012-09-28 | 2013-01-09 | 苏州大学 | 一种信号的谱峭度滤波方法及相关装置 |
Also Published As
Publication number | Publication date |
---|---|
US20160126915A1 (en) | 2016-05-05 |
EP3015996A4 (en) | 2017-02-22 |
US9559658B2 (en) | 2017-01-31 |
JP2015007850A (ja) | 2015-01-15 |
EP3015996B1 (en) | 2018-05-30 |
WO2014208319A1 (ja) | 2014-12-31 |
EP3015996A1 (en) | 2016-05-04 |
CN105324762A (zh) | 2016-02-10 |
JP6216550B2 (ja) | 2017-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105324762B (zh) | 滤波器系数群计算设备以及滤波器系数群计算方法 | |
CN107731223B (zh) | 语音活性检测方法、相关装置和设备 | |
US10014002B2 (en) | Real-time audio source separation using deep neural networks | |
CN104423707B (zh) | 使用分段和组合的触觉转换系统 | |
CN106658284B (zh) | 频域中的虚拟低音的相加 | |
CN108735202A (zh) | 用于小占用资源关键词检索的卷积递归神经网络 | |
CN106463106A (zh) | 用于音频接收的风噪声降低 | |
WO2012050705A1 (en) | Automatic equalization using adaptive frequency-domain filtering and dynamic fast convolution | |
CN110070884B (zh) | 音频起始点检测方法和装置 | |
CN109819375A (zh) | 调节音量的方法与装置、存储介质、电子设备 | |
CN113299313B (zh) | 音频处理方法、装置及电子设备 | |
EP4254408A1 (en) | Speech processing method and apparatus, and apparatus for processing speech | |
CN106653049A (zh) | 时域中的虚拟低音的相加 | |
WO2020228226A1 (zh) | 一种纯音乐检测方法、装置及存储介质 | |
WO2022247494A1 (zh) | 音频信号补偿方法及装置、耳机、存储介质 | |
Shankar et al. | Efficient two-microphone speech enhancement using basic recurrent neural network cell for hearing and hearing aids | |
CN107481727A (zh) | 一种基于电音基调控制的音频信号处理方法及系统 | |
WO2022256577A1 (en) | A method of speech enhancement and a mobile computing device implementing the method | |
JP6556463B2 (ja) | フィルタ生成装置、フィルタ生成方法およびフィルタ生成プログラム | |
CN101221767B (zh) | 人声语音加强装置与应用于其上的方法 | |
CN116027911B (zh) | 一种基于音频信号的无接触手写输入识别方法 | |
CN110085214B (zh) | 音频起始点检测方法和装置 | |
CN107251574B (zh) | 相位控制信号生成设备、相位控制信号生成方法及计算机可读介质 | |
CN115273822A (zh) | 音频处理方法、装置、电子设备及介质 | |
US20230245668A1 (en) | Neural network-based audio packet loss restoration method and apparatus, and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |