CN109637550B - Method and system for controlling elevation angle of sound source - Google Patents
Method and system for controlling elevation angle of sound source Download PDFInfo
- Publication number
- CN109637550B CN109637550B CN201811613487.5A CN201811613487A CN109637550B CN 109637550 B CN109637550 B CN 109637550B CN 201811613487 A CN201811613487 A CN 201811613487A CN 109637550 B CN109637550 B CN 109637550B
- Authority
- CN
- China
- Prior art keywords
- frequency
- altitude
- sound source
- calculating
- hrtf
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000008447 perception Effects 0.000 claims abstract description 32
- 230000003595 spectral effect Effects 0.000 claims description 18
- 230000006870 function Effects 0.000 claims description 8
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000001228 spectrum Methods 0.000 abstract description 10
- 230000000694 effects Effects 0.000 abstract description 4
- 210000003128 head Anatomy 0.000 description 6
- 230000010354 integration Effects 0.000 description 3
- 230000005236 sound signal Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000003454 tympanic membrane Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
- G10L19/02—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders
-
- 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/307—Frequency adjustment, e.g. tone control
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Acoustics & Sound (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
The application provides a sound source altitude angle control method and system, wherein the method comprises the following steps: calculating a correlation coefficient of the amplitude value of each frequency point in the HRTF database along with the change of the altitude angle; obtaining at least one center frequency by calculating a correlation coefficient; obtaining at least one filter bank by calculating at least one central frequency, then convolving the at least one filter bank with HRTFs of different altitude angles, and modeling to obtain a sub-band energy model; and inputting the sound source with the given altitude angle into the sub-band energy model, calculating the sub-band energy ratio of the target altitude angle signal to the given altitude angle signal, and then adjusting the sub-band energy of the sound source with the given altitude angle to enable the sound source perception angle to be the target altitude angle. According to the invention, the division of the sound source frequency band is more in accordance with the characteristic that the HRTF frequency spectrum changes along with the altitude angle, so that the frequency band with high altitude angle correlation has higher frequency resolution, and a better sound source altitude angle control effect is realized.
Description
Technical Field
The invention relates to the technical field of three-dimensional audio, in particular to a method and a system for controlling a height angle of a sound source.
Background
The elevation angle control method for replaying a sound source in the three-dimensional audio technology mainly modifies the frequency spectrum of the sound source, so that the human auditory system can sense the change of the elevation angle of the sound source.
Since a Head Related Transfer Function (HRTF) represents a process of a sound wave from a sound source to an eardrum of a human ear through physiological structures such as a Head, an auricle, a trunk and the like, a conventional altitude control method generally modifies a signal spectrum by convolving the sound signal with HRTFs of different altitude angles, thereby controlling the altitude of the sound source. However, the structure of the HRTF is complex, and it is difficult to meet the real-time requirement of a three-dimensional audio playback system, so a more simplified method is to modify the band energy of a sound source to achieve the control of the elevation angle. The principle is that according to the characteristic that the human auditory system is insensitive to the frequency spectrum details of the sound signal, an auditory filter bank (such as a Mel filter) is utilized to filter the sound signal to obtain the total energy of each frequency band, and the control of the sound source height angle is realized by adjusting the energy of the corresponding frequency band. However, the method based on the auditory filter is not perfect for dividing the frequency band, for example, the mel filter is designed according to the pitch characteristics, the frequency resolution thereof is reduced along with the increase of the frequency, and the study on the HRTF shows that the characteristic of the signal spectrum changing along with the height angle is different from the frequency and is not consistent with the characteristic of the auditory filter such as the mel filter.
Disclosure of Invention
In order to solve the above problems, the present invention provides a method and a system for controlling elevation angle of a sound source.
In order to achieve the above purpose, the embodiments of the present application adopt the following technical solutions:
in a first aspect, the present application provides a method for controlling an elevation angle of a sound source, including: calculating a correlation coefficient of the amplitude value of each frequency point in the head related transfer function HRTF database along with the change of the altitude angle; obtaining at least one center frequency by calculating the correlation coefficient; obtaining at least one filter bank by calculating the at least one central frequency, then convolving the at least one filter bank with HRTFs of different altitude angles, and modeling to obtain a sub-band energy model; and inputting the sound source with a given altitude angle into the sub-band energy model, calculating the sub-band energy ratio of the target altitude angle signal to the given altitude angle signal, and then adjusting the sub-band energy of the sound source with the given altitude angle to enable the sound source perception angle to be the target altitude angle.
In another possible implementation, the calculating a correlation coefficient of the amplitude value of each frequency point in the HRTF database varying with the altitude angle includes: calculating a correlation coefficient of the amplitude value of each frequency point along with the change of the altitude angle by utilizing Fisher distribution, wherein the calculation formula is as follows:
wherein,
representing the spectral amplitude value, μ, of the HRTF at elevation angle i of the jth subject in the HRTF databaseiRepresents the mean of all subject HRTF spectral amplitude values at altitude i, μ represents the mean of all altitude and all subject HRTF spectral amplitude values, N represents the total number of samples for altitude, and M represents the total number of samples for subject.
In another possible implementation, the obtaining at least one center frequency by calculating the correlation coefficient specifically includes: converting the correlation coefficient from a frequency domain to an auditory perception domain to obtain the correlation coefficient of the perception domain; fitting the correlation coefficients of the perception domains into a continuous curve, then integrating to obtain an integral value, and equally dividing the continuous curve into at least two regions; and acquiring boundary points of each region and the region boundary, and converting the at least one boundary point into a frequency domain to obtain at least one center frequency.
In another possible implementation, the method for converting the correlation coefficient from the frequency domain to the auditory perception domain to obtain the correlation coefficient of the perception domain is as follows:
wherein f represents the frequency value corresponding to the correlation coefficient, ERBNNumber represents the number of the auditory perception domain to which this frequency corresponds.
In another possible implementation, the method for converting the at least one intersection point onto the frequency domain is: and carrying out ERB inverse operation on the at least one boundary point to obtain at least one corresponding center frequency.
In another possible implementation, the at least one center frequency is calculated to obtain at least one filter bank, where the filter bank is generated by:
wherein f iscThe center frequency of each sub-filter is represented,
indicating the initial phase and b the sub-filter bandwidth.
In another possible implementation, the convolving the at least one filter bank with HRTFs of different elevation angles and modeling to obtain a subband energy model specifically includes:
the modeling adopts a polynomial regression mode to carry out modeling, and the specific method comprises the following steps:
An＝{a0,n,a1,n,…,aK,n,}
wherein A isn＝{a0,n,a1,n,…,aK,nExpressed as coefficients of respective polynomials for the nth frequency band, K is expressed as an order of polynomial modeling, AnEach term coefficient represented can be obtained by the minimum mean square error of the modeling result and the band energy of the actual HRTF.
In a second aspect, the present application provides a sound source elevation angle control system, comprising: the generating unit is used for calculating a correlation coefficient of the amplitude value of each frequency point in the head related transfer function HRTF database along with the change of the altitude angle; and obtaining at least one center frequency by calculating the correlation coefficient; obtaining at least one filter bank by calculating the at least one central frequency, then convolving the at least one filter bank with HRTFs of different altitude angles, and modeling to obtain a sub-band energy model; and the processing unit is used for inputting the sound source with a given altitude angle into the sub-band energy model, calculating the sub-band energy ratio of the target altitude angle signal to the given altitude angle signal, and then adjusting the sub-band energy of the sound source with the given altitude angle to enable the sound source perception angle to be the target altitude angle.
In another possible implementation, the generating unit is specifically configured to: calculating a correlation coefficient of the amplitude value of each frequency point along with the change of the altitude angle by utilizing Fisher distribution, wherein the calculation formula is as follows:
wherein,
representing the spectral amplitude value, μ, of the HRTF at elevation angle i of the jth subject in the HRTF databaseiRepresents the mean of all subject HRTF spectral amplitude values at altitude i, μ represents the mean of all altitude and all subject HRTF spectral amplitude values, N represents the total number of samples for altitude, and M represents the total number of samples for subject.
In another possible implementation, the modeling manner is specifically: the modeling adopts a polynomial regression mode to carry out modeling, and the specific method comprises the following steps:
An＝{a0,n,a1,n,…,aK,n,}
wherein A isn＝{a0,n,a1,n,…,aK,nExpressed as coefficients of respective polynomials for the nth frequency band, K is expressed as an order of polynomial modeling, AnEach term coefficient represented can be obtained by the minimum mean square error of the modeling result and the band energy of the actual HRTF.
The invention provides a sound source altitude angle control method, which enables a frequency band with high altitude angle correlation to have higher frequency resolution and realizes better sound source altitude angle control effect by enabling the division of a sound source frequency band to better conform to the characteristic that an HRTF frequency spectrum changes along with the altitude angle.
Drawings
The drawings that accompany the detailed description can be briefly described as follows.
Fig. 1 is a flowchart of a sound source elevation angle control method according to an embodiment of the present disclosure;
fig. 2 is a flowchart of a method for generating a center frequency of a subband filter according to an embodiment of the present disclosure;
fig. 3 is a block diagram of a sound source elevation angle control system according to an embodiment of the present disclosure.
Detailed Description
Reference will now be made in detail to embodiments of the present application, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the drawings are exemplary only for the purpose of explaining the present application and are not to be construed as limiting the present application.
Fig. 1 is a flowchart of a sound source elevation angle control method according to an embodiment of the present disclosure. The sound source altitude control method shown in fig. 1 specifically includes the following steps:
step S102, calculating a correlation coefficient of the amplitude value of each frequency point in the head related transfer function HRTF database along with the change of the altitude angle.
Specifically, calculating a correlation coefficient of the amplitude value of each frequency point changing along with the altitude angle in the HRTF database, and calculating the correlation coefficient of the amplitude value of each frequency point changing along with the altitude angle by using Fisher distribution, wherein the calculation formula is as follows:
wherein,
representing the spectral amplitude value, μ, of the HRTF at elevation angle i of the jth subject in the HRTF databaseiRepresents the mean of all subject HRTF spectral amplitude values at altitude i, μ represents the mean of all altitude and all subject HRTF spectral amplitude values, N represents the total number of samples for altitude, and M represents the total number of samples for subject.
Step S104, at least one center frequency is obtained by calculating the correlation coefficient.
Fig. 2 is a flowchart of a method for generating a center frequency of a subband filter according to an embodiment of the present disclosure. As shown in fig. 2, in an embodiment, the method for generating the center frequency of the subband filter includes the following specific steps:
step S202, converting the relation number from the frequency domain to the auditory perception domain to obtain the correlation coefficient of the perception domain.
Specifically, the correlation coefficient of the amplitude value of each frequency point generated in step S102, which varies with the altitude angle, is converted from the frequency domain to the auditory perception domain, so as to obtain the correlation coefficient of the perception domain. The method for converting the frequency domain correlation coefficient into the perceptual domain correlation coefficient is as follows:
wherein f represents the frequency value corresponding to the correlation coefficient, ERBNNumber represents the number of the auditory perception domain to which this frequency corresponds.
Step S204, the correlation coefficients of the perception domain are fitted into a continuous curve, integration is carried out to obtain an integral value, and the continuous curve is equally divided into at least two areas.
Specifically, the correlation coefficients of the sensing domain are fitted into a continuous curve, the integral value of the continuous curve is obtained by means of numerical integration, the continuous curve fitted into the sensing domain is divided into a plurality of regions, and the regions are equally divided.
Preferably, the present invention divides the curve into 25 regions, although equally dividing into other numbers of regions is also possible.
Step S206, obtaining the boundary point of each region and the boundary of the regions, and converting at least one boundary point to a frequency domain to obtain at least one center frequency.
Specifically, 24 boundary points among the 25 equally divided regions are obtained, and then the 24 boundary points are subjected to ERB inverse operation to obtain corresponding 24 frequency values, namely, the center frequency.
Step S106, at least one filter bank is obtained by calculating at least one central frequency, and then at least one filter bank is convoluted with the HRTFs with different altitude angles and modeled to obtain a sub-band energy model.
Preferably, the filter selected by the present application is a Gamma tone filter.
Specifically, the 24 frequency values obtained in step S206 are processed by a Gamma tone filter, so as to obtain a filter bank of 24 channels. The method for generating the filter bank comprises the following steps:
wherein f iscThe center frequency of each sub-filter is represented,
indicating the initial phase and b the sub-filter bandwidth.
Convolving the filter banks of the 24 channels with the HRTFs of different elevation angles to obtain the energy of 24 sub-bands corresponding to the different elevation angles; and finally, modeling the energy of 24 sub-bands in a polynomial regression mode to obtain a sub-band energy model. The modeling expression mode of the polynomial is as follows:
An＝{a0,n,a1,n,…,aK,n,}
wherein A isn＝{a0,n,a1,n,…,aK,nExpressed as coefficients of respective polynomials for the nth frequency band, K is expressed as an order of polynomial modeling, AnEach term coefficient represented can be obtained by the minimum mean square error of the modeling result and the band energy of the actual HRTF.
Step S108, inputting the sound source with the given altitude angle into the sub-band energy model, calculating the sub-band energy ratio of the target altitude angle signal and the given altitude angle signal, and then enabling the sound source perception angle to be the target altitude angle by adjusting the sub-band energy of the sound source with the given altitude angle.
Specifically, for a sound source with a given known elevation angle, the sound source is input into a sub-band energy model, the ratio of the sub-band energy of a target elevation angle signal and the sub-band energy of the known elevation angle signal is calculated, and then the frequency spectrum energy of the known sound source is adjusted to enable the perception angle of the sound source to be the target elevation angle.
According to the invention, the division of the sound source frequency band is more in accordance with the characteristic that the HRTF frequency spectrum changes along with the altitude angle, so that the frequency band with high altitude angle correlation has higher frequency resolution, and a better sound source altitude angle control effect is realized.
Fig. 3 is a block diagram of a sound source elevation angle control system according to an embodiment of the present disclosure. As shown in fig. 3, the sound source altitude control system includes: a generating unit 301 and a processing unit 302.
The generating unit is used for calculating a correlation coefficient of the amplitude value of each frequency point in the head related transfer function HRTF database along with the change of the altitude angle; and obtaining at least one center frequency by calculating the correlation coefficient; and obtaining at least one filter bank by calculating at least one central frequency, then convolving the at least one filter bank with HRTFs of different altitude angles, and modeling to obtain a sub-band energy model.
In one implementation, calculating a correlation coefficient of the amplitude value of each frequency point along with the change of the altitude angle in the HRTF database, and calculating the correlation coefficient of the amplitude value of each frequency point along with the change of the altitude angle by using Fisher distribution, wherein the calculation formula is as follows:
wherein,
representing the spectral amplitude value, μ, of the HRTF at elevation angle i of the jth subject in the HRTF databaseiRepresents the mean of all subject HRTF spectral amplitude values at altitude i, μ represents the mean of all altitude and all subject HRTF spectral amplitude values, N represents the total number of samples for altitude, and M represents the total number of samples for subject.
In one embodiment, firstly, the correlation coefficient of the amplitude value of each frequency point changing with the altitude angle is converted from the frequency domain to the auditory perception domain, and the correlation coefficient of the perception domain is obtained. The method for converting the frequency domain correlation coefficient into the perceptual domain correlation coefficient is as follows:
wherein f represents the frequency value corresponding to the correlation coefficient, ERBNNumber represents the number of the auditory perception domain to which this frequency corresponds.
Then, the correlation coefficients of the perception domain are fitted into a continuous curve, the integral value of the perception domain is obtained in a numerical integration mode, the perception domain is fitted into the continuous curve to be divided into 25 regions, 24 boundary points among the 25 equally divided regions are obtained, and then the 24 boundary points are subjected to ERB inverse operation to obtain corresponding 24 frequency values, namely the central frequency.
And then, processing the obtained 24 frequency values through a Gamma tone filter to obtain a filter bank with 24 channels. The method for generating the filter bank comprises the following steps:
wherein f iscThe center frequency of each sub-filter is represented,
indicating the initial phase and b the sub-filter bandwidth.
And finally, convolving the filter banks of the 24 channels with the HRTFs of different altitude angles to obtain the energy of 24 sub-bands corresponding to the different altitude angles, and modeling the energy of the 24 sub-bands in a polynomial regression mode to obtain a sub-band energy model. The modeling expression mode of the polynomial is as follows:
An＝{a0,n,a1,n,…,aK,n,}
wherein A isn＝{a0,n,a1,n,…,aK,nExpressed as coefficients of respective polynomials for the nth frequency band, K is expressed as an order of polynomial modeling, AnEach term coefficient represented can be obtained by the minimum mean square error of the modeling result and the band energy of the actual HRTF.
And the processing unit is used for inputting the sound source with the given altitude angle into the sub-band energy model, calculating the sub-band energy ratio of the target altitude angle signal to the given altitude angle signal, and then enabling the sound source perception angle to be the target altitude angle by adjusting the sub-band energy of the sound source with the given altitude angle.
The method comprises the steps of inputting a sound source with a given known elevation angle into a sub-band energy model, calculating the ratio of sub-band energy of a target elevation angle signal to sub-band energy of the known elevation angle signal, and adjusting the frequency spectrum energy of the known sound source to enable the sensing angle to be the target elevation angle.
According to the invention, the division of the sound source frequency band is more in accordance with the characteristic that the HRTF frequency spectrum changes along with the altitude angle, so that the frequency band with high altitude angle correlation has higher frequency resolution, and a better sound source altitude angle control effect is realized.
Finally, the description is as follows: the above embodiments are only used to illustrate the technical solutions of the present application, but not to limit the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present application.
Claims (10)
1. A sound source elevation angle control method, comprising:
calculating a correlation coefficient of the amplitude value of each frequency point in the head related transfer function HRTF database along with the change of the altitude angle;
obtaining at least one center frequency by calculating the correlation coefficient;
obtaining at least one filter bank by calculating the at least one central frequency, then convolving the at least one filter bank with HRTFs of different altitude angles, and modeling to obtain a sub-band energy model;
and inputting the sound source with a given altitude angle into the sub-band energy model, calculating the sub-band energy ratio of the target altitude angle signal to the given altitude angle signal, and then adjusting the sub-band energy of the sound source with the given altitude angle to enable the sound source perception angle to be the target altitude angle.
2. The method as claimed in claim 1, wherein the calculating of the correlation coefficient of the amplitude value of each frequency point in the HRTF database with the change of the altitude angle comprises: calculating a correlation coefficient of the amplitude value of each frequency point along with the change of the altitude angle by utilizing Fisher distribution, wherein the calculation formula is as follows:
wherein,
representing the spectral amplitude value, μ, of the HRTF at elevation angle i of the jth subject in the HRTF databaseiRepresents the mean of all subject HRTF spectral amplitude values at altitude i, μ represents the mean of all altitude and all subject HRTF spectral amplitude values, N represents the total number of samples for altitude, and M represents the total number of samples for subject.
3. The method according to claim 2, wherein the obtaining at least one center frequency by calculating the correlation coefficient specifically comprises:
converting the correlation coefficient from a frequency domain to an auditory perception domain to obtain the correlation coefficient of the perception domain;
fitting the correlation coefficients of the perception domains into a continuous curve, then integrating to obtain an integral value, and equally dividing the continuous curve into at least two regions;
and acquiring boundary points of each region and the region boundary, and converting the at least one boundary point into a frequency domain to obtain at least one center frequency.
4. The method according to claim 3, wherein the method for converting the correlation coefficient from the frequency domain to the auditory perception domain to obtain the correlation coefficient of the perception domain comprises:
wherein f represents the frequency value corresponding to the correlation coefficient, ERBNNumber represents the number of the auditory perception domain to which this frequency corresponds.
5. The method of claim 3, wherein the method for transforming the at least one intersection point into the frequency domain comprises: and performing ERB inverse operation on the at least one junction point through an equivalent rectangular bandwidth to obtain at least one corresponding center frequency.
6. The method according to claim 1, wherein said calculating the at least one center frequency results in at least one filter bank, and wherein the filter bank is generated by:
wherein f iscThe center frequency of each sub-filter is represented,
indicating the initial phase and b the sub-filter bandwidth.
7. The method of claim 1, wherein convolving the at least one filter bank with HRTFs of different elevation angles and modeling the at least one filter bank to obtain a subband energy model specifically comprises:
the modeling adopts a polynomial regression mode to carry out modeling, and the specific method comprises the following steps:
An＝{a0,n,a1,n,…,aK,n,}
where θ represents the elevation angle of the HRTF, An＝{a0,n,a1,n,…,aK,nExpressed as coefficients of respective polynomials for the nth frequency band, K is expressed as an order of polynomial modeling, AnEach term coefficient represented can be obtained by the minimum mean square error of the modeling result and the band energy of the actual HRTF.
8. A sound source elevation angle control system comprising a generating unit and a processing unit, characterized by comprising:
the generating unit is used for calculating a correlation coefficient of the amplitude value of each frequency point in the head related transfer function HRTF database along with the change of the altitude angle; and
obtaining at least one center frequency by calculating the correlation coefficient; and
obtaining at least one filter bank by calculating the at least one central frequency, then convolving the at least one filter bank with HRTFs of different altitude angles, and modeling to obtain a sub-band energy model;
and the processing unit is used for inputting the sound source with a given altitude angle into the sub-band energy model, calculating the sub-band energy ratio of the target altitude angle signal to the given altitude angle signal, and then adjusting the sub-band energy of the sound source with the given altitude angle to enable the sound source perception angle to be the target altitude angle.
9. The system according to claim 8, wherein the generating unit is specifically configured to: calculating a correlation coefficient of the amplitude value of each frequency point along with the change of the altitude angle by utilizing Fisher distribution, wherein the calculation formula is as follows:
wherein,
representing the spectral amplitude value, μ, of the HRTF at elevation angle i of the jth subject in the HRTF databaseiRepresents the mean of all subject HRTF spectral amplitude values at altitude i, μ represents the mean of all altitude and all subject HRTF spectral amplitude values, N represents the total number of samples for altitude, and M represents the total number of samples for subject.
10. The system according to claim 8, wherein the modeling approach is specifically:
the modeling adopts a polynomial regression mode to carry out modeling, and the specific method comprises the following steps:
An＝{a0,n,a1,n,…,aK,n,}
where θ represents the elevation angle of the HRTF, An＝{a0,n,a1,n,…,aK,nExpressed as coefficients of respective polynomials for the nth frequency band, K is expressed as an order of polynomial modeling, AnEach term coefficient represented can be obtained by the minimum mean square error of the modeling result and the band energy of the actual HRTF.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811613487.5A CN109637550B (en) | 2018-12-27 | 2018-12-27 | Method and system for controlling elevation angle of sound source |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811613487.5A CN109637550B (en) | 2018-12-27 | 2018-12-27 | Method and system for controlling elevation angle of sound source |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109637550A CN109637550A (en) | 2019-04-16 |
| CN109637550B true CN109637550B (en) | 2020-11-24 |
Family
ID=66078516
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811613487.5A Active CN109637550B (en) | 2018-12-27 | 2018-12-27 | Method and system for controlling elevation angle of sound source |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109637550B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2620796A (en) * | 2022-07-22 | 2024-01-24 | Sony Interactive Entertainment Europe Ltd | Methods and systems for simulating perception of a sound source |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5802180A (en) * | 1994-10-27 | 1998-09-01 | Aureal Semiconductor Inc. | Method and apparatus for efficient presentation of high-quality three-dimensional audio including ambient effects |
| EP1058481A2 (en) * | 1999-05-29 | 2000-12-06 | Central Research Laboratories Limited | A method of modifying one or more original head related transfer functions |
| JP2005005949A (en) * | 2003-06-11 | 2005-01-06 | Matsushita Electric Ind Co Ltd | Transfer function interpolation method |
| CN101263739A (en) * | 2005-09-13 | 2008-09-10 | Srs实验室有限公司 | Systems and methods for audio processing |
| CN101483797A (en) * | 2008-01-07 | 2009-07-15 | 昊迪移通(北京)技术有限公司 | Head-related transfer function generation method and apparatus for earphone acoustic system |
| CN106162499A (en) * | 2016-07-04 | 2016-11-23 | 大连理工大学 | The personalized method of a kind of related transfer function and system |
| CN107205207A (en) * | 2017-05-17 | 2017-09-26 | 华南理工大学 | A kind of approximate acquisition methods of virtual sound image based on middle vertical plane characteristic |
| CN107480100A (en) * | 2017-07-04 | 2017-12-15 | 中国科学院自动化研究所 | Head-position difficult labor modeling based on deep-neural-network intermediate layer feature |
| WO2018026963A1 (en) * | 2016-08-03 | 2018-02-08 | Hear360 Llc | Head-trackable spatial audio for headphones and system and method for head-trackable spatial audio for headphones |
-
2018
- 2018-12-27 CN CN201811613487.5A patent/CN109637550B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5802180A (en) * | 1994-10-27 | 1998-09-01 | Aureal Semiconductor Inc. | Method and apparatus for efficient presentation of high-quality three-dimensional audio including ambient effects |
| EP1058481A2 (en) * | 1999-05-29 | 2000-12-06 | Central Research Laboratories Limited | A method of modifying one or more original head related transfer functions |
| JP2005005949A (en) * | 2003-06-11 | 2005-01-06 | Matsushita Electric Ind Co Ltd | Transfer function interpolation method |
| CN101263739A (en) * | 2005-09-13 | 2008-09-10 | Srs实验室有限公司 | Systems and methods for audio processing |
| CN101483797A (en) * | 2008-01-07 | 2009-07-15 | 昊迪移通(北京)技术有限公司 | Head-related transfer function generation method and apparatus for earphone acoustic system |
| CN106162499A (en) * | 2016-07-04 | 2016-11-23 | 大连理工大学 | The personalized method of a kind of related transfer function and system |
| WO2018026963A1 (en) * | 2016-08-03 | 2018-02-08 | Hear360 Llc | Head-trackable spatial audio for headphones and system and method for head-trackable spatial audio for headphones |
| CN107205207A (en) * | 2017-05-17 | 2017-09-26 | 华南理工大学 | A kind of approximate acquisition methods of virtual sound image based on middle vertical plane characteristic |
| CN107480100A (en) * | 2017-07-04 | 2017-12-15 | 中国科学院自动化研究所 | Head-position difficult labor modeling based on deep-neural-network intermediate layer feature |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109637550A (en) | 2019-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6934979B2 (en) | Methods and devices for rendering audio field representations for audio playback | |
| KR102149214B1 (en) | Audio signal processing method and apparatus for binaural rendering using phase response characteristics | |
| JP4718559B2 (en) | Method and apparatus for individualizing HRTFs by modeling | |
| Smith et al. | Bark and ERB bilinear transforms | |
| JP6454027B2 (en) | Audio signal processing apparatus and method for binaural rendering | |
| EP3069530B1 (en) | Method and device for compressing and decompressing sound field data of an area | |
| US20130136282A1 (en) | System and Method for Spectral Personalization of Sound | |
| CN101695148B (en) | Multi-channel wide dynamic range compressing system for digital hearing aid | |
| Zhong et al. | Head-related transfer functions and virtual auditory display | |
| CN113889125A (en) | Audio generation method and device, computer equipment and storage medium | |
| Sakamoto et al. | Sound-space recording and binaural presentation system based on a 252-channel microphone array | |
| Choueiri | Optimal crosstalk cancellation for binaural audio with two loudspeakers | |
| CA2983359C (en) | An audio signal processing apparatus and method | |
| CN103731780A (en) | Digital hearing-aid signal processing method | |
| US20170272881A1 (en) | Audio signal processing apparatus and method for modifying a stereo image of a stereo signal | |
| US9084050B2 (en) | Systems and methods for remapping an audio range to a human perceivable range | |
| CN109637550B (en) | Method and system for controlling elevation angle of sound source | |
| CN104967952B (en) | A kind of personalized method based on HRTF structural models and subjective feedback | |
| Ramos et al. | A parallel approach to HRTF approximation and interpolation based on a parametric filter model | |
| Li et al. | Modeling perceived externalization of a static, lateral sound image | |
| CN113766396A (en) | Loudspeaker control | |
| Marelli et al. | Efficient approximation of head-related transfer functions in subbands for accurate sound localization | |
| CN106569780A (en) | Real-time audio processing method and system for multi-channel digital audio signal | |
| CN109688531B (en) | Method for acquiring high-sound-quality audio conversion information, electronic device and recording medium | |
| JPWO2020066542A1 (en) | Acoustic object extraction device and acoustic object extraction method |
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 |