The present invention refers to an audio signal processing method for enhancing the quality and/or other characteristic of an audio signal. This method corrects a non-linearity of electro-acoustic transducers in an audio chain by taking into account also a non-linear psychoacoustical characteristics of the human ear by adding non-linearities in the audio chain in a controlled manner. Furthermore, the present invention relates to a device/apparatus for the implementation of said method and audio chain configured to correct the non-linearity of electroacoustic transducers, taking into account also the non-linear psychoacoustical characteristics of the human ear. The audio chain contains at least one apparatus for the implementation of the audio signal processing method.
Nowadays, the audio chain before the electroacoustic converter displays impeccable features. It is not known why some audio chain components with greater distortions produce better sound than components with lower distortion. Some amplifiers have incorporated vacuum tubes in order to sound better, whereas other employ a small feedback loop to intensify non-linearities of the components. The audio chain distortions before the electroacoustic transducer do not mean that it will sound better or worse. Two different electroacoustic transducers that sound good on their audio chains, will sound worse when they swap places. One of the reasons thereof is that the audio chain before the electroacoustic transducer has non-linearities that reduce its non-linearities which makes it sound better than on the other audio chain.
Technical problem that gets solved with the present invention is a method and an apparatus for audio signal processing in audio chain, that correct non-linearity of electroacoustic transducers in audio chain, taking into account also the non-linear psychoacoustic feature of the human ear.
Non-linearities of the electroacoustic transducers have been known for some time now. Non-linear distortions characterize the entire electroacoustic reproduction chain, from the sound recording process on the sound recording medium all the way to the reproduction of sound from the sound recording medium, the amplifier and the loudspeaker itself. There are many publications documenting these non-linearities. The application of non-linearities in musical instruments in order to change the sound has also been known for a while. People do not perceive some non-linearities in the sound, whereas others are perceived, even though they have the same acoustic energy as described in the article Amplifier Musicality—A Study of Amplifier Harmonic Distortion Spectrum Analysis by Jean Hiraga. The document U.S. Pat. No. 5,133,015 discloses the process and the apparatus for audio signal processing, more precisely, the technique that permits various audio signal distortion grades comprising the audio signal distortion to a certain grade. The document US2011255701 discloses the electronic circuit and the audio enhancement method, particularly the electronic circuit that can introduce a predictive and controllable harmonic distortion that increases with an increased signal amplitude. The document US2015249889 discloses the system and the method for digital audio signal processing by extending the loudspeaker frequency response and reducing or eliminating non-linear loudspeaker distortion. An audio signal can be extended by applying a digital linear filter, based on a modified loudspeaker frequency response. A non-linear distortion of a loudspeaker can be cancelled or reduced by a digital non-linear filter based on a reverse parametric model of the loudspeaker.
Most of the known conventional approaches related to audio signal processing with the view to enhancing the quality and/or other characteristics of audio signal do not take into consideration also the non-linear psychoacoustical characteristic of the human ear.
The present invention relates to an audio signal processing method and apparatus in an audio chain that correct a non-linearity of electroacoustic transducers in the audio chain, taking into consideration also a non-linear psychoacoustical characteristics of the human ear by adding non-linearities in the audio chain in a controlled manner, in order to obtain a better acoustic image and more details when reproducing the sound by using approximation of the quadratic and a fifth degree polynomial function in some range.
According to the present invention, a method comprises approximating of the non-linear psychoacoustical characteristic of the human ear by a fifth-degree polynomial and adding of at least one non-linear element in front of at least one electroacoustic transducer in the audio chain, whereby the non-linear element has a function of adding a non-linearity in the audio chain that corrects the non-linearity of at least one electroacoustic transducer and/or the non-linearity of the approximated psychoacoustical characteristic of the human ear for a pressure change by the human ear up to pΔ.
The audio chain for implementing of said audio signal processing method, according to the present invention, is configured to correct the non-linearity of electroacoustic transducers in the audio chain, taking into account also the non-linear psychoacoustical characteristic of the human ear. Said audio chain contains at least one apparatus for implementing of the audio signal processing method. The aforementioned apparatus has the function of adding the non-linearity to the audio chain that corrects the non-linearity of at least one electroacoustic transducer and/or the non-linearity of the approximate psychoacoustical characteristic of the human ear for the pressure change by the human ear up to pΔ.
The method of the present invention, the apparatus and the audio chain reduce limitations of the electroacoustic transducers as well as of the human ear by adding non-linearities that, ultimately reduce non-linearities of an entire audio chain with the human ear, i.e. adding non-linearities to the audio chain so that an audio chain characteristic reduces the non-linearity of the human ear polynomial approximation to the change of pressures pΔ=±1 Pa.
In the following, the invention shall be described in detail with reference to the drawings, wherein:
with asymptotes;
with asymptotes;
A method of the present invention takes into consideration one non-linearity of an electroacoustic transducer and non-linearity of the human ear.
According to the present invention, an audio signal processing method in an audio chain, which corrects the non-linearity of the electroacoustic transducers in the audio chain, taking into account also the non-linear psychoacoustical characteristic of the human ear, comprises approximating the psychoacoustical characteristics of the human ear by a fifth degree polynomial function, and adding of at least one non-linear element 4 in front of at least one electroacoustic transducer in the audio chain, said non-linear element 4 has a function to add a non-linearity in the audio chain that corrects the non-linearity of at least one electroacoustic transducer and/or the non-linearity of the approximated psychoacoustical characteristic of the human ear for a pressure change by the human ear up to pΔ. According to the present method, the non-linear element 4 reduces the non-linearity of the electroacoustic transducer by applying a quadratic non-linearity which is an inverse function of ax+bx2 where x is a relative membrane excursion or a relative force on a membrane of the electroacoustic transducer, a and b are positive constants.
According to the one embodiment of the invention, the non-linear element 4 reduces the non-linearity of the psychoacoustical characteristic of the human ear x−a x2−b x3−c x4−d x5 by applying the function which reduces at least two times the non-linearities introduced by the members x2, x3 and x4, wherein the constants
stay within the tolerances±30% for each constant and x is a relative pressure by the human ear.
According to the other embodiment of the invention, the non-linear element 4 reduces the non-linearity of the psychoacoustical characteristic of the human ear by applying the hyperbolic function
where x is the relative pressure by the human ear.
According to another embodiment of the invention, the non-linear element 4 reduces the non-linearity of the psychoacoustical characteristic of the human ear by applying the function x1.5, where x is the relative pressure by the human ear.
The present method will be further described in more detail and in accordance with the embodiment of the audio chain according to the present invention.
The non-linearity within the electroacoustic transducer is defined by an adiabatic process defined as:
0pVn=const [1]
Said non-linearity within the electroacoustic transducer affects the quality of sound. In the case of the electroacoustic transducer that produces sound by moving the membrane, the air by the membrane changes the pressure by adiabatic process. The volume of air being compressed is unknown. However, changes in air pressure can be measured. A larger volume of air being compressed requires a greater membrane excursion for the same pressure and vice versa. As the air pressure changes by adiabatic process, the same membrane excursion in the direction that increases the pressure, will create greater pressure change that the excursion in the opposite direction. We will consider two ideal cases. In both cases the mass of the membrane is negligibly small, and it is rigid. In the first case, the membrane excursion is linear and the volume of the compressed air changes linearly with the membrane excursion. We will use the adiabatic process of air. The initial air pressure is atmospheric pressure. Adiabatic equation for air is:
pV1.4=const [2]
As the membrane moves, the volume changes, which changes the air pressure adiabatically:
Air pressure by the membrane is:
where V0 is the initial volume we compress, and VΔ the volume change that occurs by moving the membrane. VΔ has the negative sign because the volume decreases as the membrane moves forward. The initial conditions will be: const=p0, V0=1, and the volume change VΔ=d, where p0 is atmospheric pressure and dr relative membrane excursion. Consequently, we can write:
If we expand the function into Taylor series according to the relative excursion d, the first five members are:
p=p0(1+1.4x+1.68x2+1.904x3+2.0944x4+ . . . ), [6]
p=p0+pΔ, [7]
where pΔis the pressure change:
pΔ=p0(1.4x+1.68x2+1.904x3+2.0944x4+ . . . ). [8]
For the pressure change of pΔ=1 Pa the relative membrane excursion is:
If we put it in Taylor series, the members after the quadratic member are negligible:
p0(1.904x3+2.0944x4+ . . . )≈0. [10]
The greatest non-linearity at normal loudness is the quadratic function of the pressure change
pΔ≈p0(1.4x+1.68x2). [11]
In the second case we have the force on the electroacoustic transducer membrane and the air volume that changes linearly with the membrane excursion. For easier calculation, we will use an isothermal process defined for ideal gas as:
pV=const. [12]
The force on the membrane is the sum of the forces on both sides of the membrane. Since we listen to the sound only from one side of the membrane, we will monitor the pressure on that side. The force for the membrane surface is:
F=A0(p1−p2), [13]
where p1 is the pressure on the side of the membrane facing us, p2 is the pressure on the opposite side of the membrane and A0 is the surface for the membrane that is constant. Pressure p1, p2 is:
where V0 is the initial volume we compress, and VΔ the volume change that occurs by moving the membrane. The initial conditions will be const=p0, V0=1 and VΔ=d, where p0 is atmospheric pressure, and dr the relative membrane excursion in the direction of listening. We get the equations for p1, p2:
The force on the membrane is:
If we assume that the relative force is Fr=F/(A0 p0), then it is:
and the relative membrane excursion is:
The pressure on the listening side is then p1=p0/(1−d) which results in:
Developed into Taylor series per the relative force Fr, we get the pressure on the side of the membrane facing us:
wherein the pressure p1 on the side of the membrane facing us is disclosed as:
p1=p0+pΔ [21]
and the pressure change pΔ on the side of listening is:
For the pressure change of pΔ=1 Pa the relative membrane excursion is:
For such a small relative force we can ignore the impact of the bigger members of Taylor series:
The greatest non-linearity at normal loudness is the quadratic function of pressure change
In both cases, we can approximate the air pressure change on the membrane by quadratic function ax+bx2 where x is the relative membrane excursion in the first case or relative pressure on membrane in the second case. If we consider a normal loudness with the pressure change±1 Pa by the human ear, the pressure on the membrane is greater, because the pressure decreases with the distance. The smaller the surface of the electroacoustic transducer membrane, other parameters being identical, the greater the pressure on it by the same loudness at the same distance. Assuming that, at 2 meters from the electroacoustic transducer, the pressure difference is ±1 Pa and the electroacoustic transducer has a surface 1.272π cm2 and an ideal dispersion in all directions without sound reflection, then the acoustic power at the membrane is equal to the power at the spherical surface at some distance from the membrane. By a sphere at a 2 meters distance, this is 4·22π m2, which makes 160000π cm2. A sound power is:
P=I·A=const [26]
where P is a power, I is an intensity and A is a surface area. If intensity I is proportional to the square of the pressure change I∝p12 then we can write p12A1=p22A2 meaning that the pressure on the membrane in the direction of listening is
As the pressure on the membrane increases, the electroacoustic transducer works in a non-linear area, influencing the quality of sound we hear. For the calculated loudness pΔ=p0(ax+bx2)
And the ratio of the quadratic component bx2 to the linear component ax is
In the first case is a=1.4, b=1.68 and pΔ=314.96 Pa, the quadratic component is 0.27% of the linear component, which is not to be ignored. In the second case is a=½, b=¼, pΔ=314.96 Pa and the quadratic component is 0.31% of the linear component, which is also not to be ignored. To reduce the quadratic non-linearity of the electroacoustic transducer in the chain before it, we incorporate the non-linear element that corrects the non-linearity of the audio chain behind it:
y=a(x+bx2) [30]
where a and b are positive constants. The easiest way to correct the non-linearity of the electroacoustic transducer is by using the non-linear element that approximates the inverse function x+bx2 which makes:
Developed into Taylor series, we get x−bx2+2bx2+2b2x3−5b3x5+ . . .
We will take the first two members of Taylor series:
y−1≈x−bx2, [32]
And we will ignore the remaining members, because their impact is negligible when x is very small. To obtain the characteristics of the non-linear element and the audio chain after it, in a(x+bx2) we replace x with x−bx2 and get a(x−2b2x3+b3x4), where |−2b2x3+b3x4|<<|bx2| is when x is very small. That way we reduced distortions by low values x, which is the case by listening of the audio chain at normal loudness, where the pressure change by the human ear is up to pΔ=±1 Pa. If the electroacoustic transducer has a smaller membrane surface, a greater pressure will be on the membrane for the same loudness at the same distance. This will increase the adiabatic distortion of the electroacoustic transducer. It is sufficient to adjust the non-linear element to reduce at least three times the quadratic non-linearity of the electroacoustic transducer to feel a significant enhancement of sound.
SET (Single Ended Triode) tube amplifiers are known to have a non-linearity greater than 1% at rated power and are not audible to the human ear. Jean Hiraga wrote an article that received a lot of attention and criticism called Amplifier Musicality—A Study of Amplifier Harmonic Distortion Spectrum Analysis where he describes the harmonic structure of the non-linearity of various amplifiers and subjectively evaluates their sound. In addition to not hearing the non-linearity of SET tube amplifiers, their non-linearity overrides details of sound that we no longer hear. If we assume that the human ear has a similar non-linearity and we do not hear it, then we would not hear it even if the non-linearity were in a part of the audio chain. It is known that the frequency sine wave f1 and the same one with added frequencies f2, f3, f4, f5, f6 that are 2, 3, 4, 5, 6 time greater than f1 where the amplitudes are: f1 at 0 db, f2 at −40 db, f3 at −50 db, f4 at −60 db, f5 at −70 db and f6 at −80 db will sound the same to the human ear (
x−ax2−bx3−cx4−dx5 [33]
Where a, b, c and d are real positive numbers and x is the relative pressure by the human ear. To determine the values of a, b, c, and d, we add non-linearities to the audio signal until we have reached the distortion of the harmonic structure of the human ear we hear. To determine the coefficient a, we use the non-linearity x+ax2 which, with an approximation of the characteristic of the human ear, gives:
x−(2a2+b)x3−(a3+3ab+c)x4−(3a2b+4ac+d)x5− . . . [34]
where we removed the member x2 and disturbed the harmonic structure of the human ear. To determine the coefficient b, we use the non-linearity x+bx3 which, with an approximation of the characteristic of the human ear, gives:
x−ax2−(2ab+c)x4−(3b2+d)x5− . . . [35]
where we removed the member x3 and disturbed the harmonic structure of the human ear. To determine the coefficient c, we use the non-linearity x+cx4 which, with an approximation of the characteristic of the human ear, gives:
x−ax2−bx3−(2ac+d)x5− . . . [36]
where we removed the member x4 and disturbed the harmonic structure of the human ear. To determine the coefficient d, we use the non-linearity x+dx5 which, with an approximation of the characteristic of the human ear, gives:
x−ax2−bx3−cx4− . . . [37]
where we removed the member x5 and disturbed the harmonic structure of the human ear. The members
within the tolerances±30% for each member were obtained through hearing tests. Approximated function of the psychoacoustic feature of the human ear is:
By applying the Lagrange-Bürmann formula, we get the following inverse function of the approximation of the human ear:
Since the coefficient of the x5 member of the approximated function of the psychoacoustical characteristic of the human ear is very small, it can be ignored, as well as the bigger members. In order to hear enough details, it is necessary to reduce at least two times the non-linearities introduced by the x2, x3 and x4 members of the approximated psychoacoustic characteristics of the human ear. The inverse function of the approximation of the psychoacoustic feature of the human ear can be derived using the hyperbolic curves
where a=0.00372, b=0.06061, c=0.002484 and d=0.01313 (
which, when developed into Taylor series, makes the first five members:
In order to see how the non-linearity of the human ear decreases, in the approximate psychoacoustic feature of the human ear
we replace x with
and obtain the first five members:
we reduced at least two times the non-linearities introduced by the x2, x3 and x4 members of the approximated psychoacoustic characteristics of the human ear.
According to the present invention, an apparatus for the implementation of the method comprises at least one non-linear element 4 in the audio chain that has the function of adding the non-linearity to the audio chain that corrects the non-linearity of at least one electroacoustic transducer and/or the non-linearity of the approximate psychoacoustical characteristic of the human ear for the pressure change by the human ear to pΔ.
A method for audio signal processing in the audio chain carried out by using the apparatus 19 illustrated in
and the amplifier/attenuator 5 after the non-linear element 4 for the positive value b. The hyperbolic non-linear element 4 is derived from the signal inverter 11, a source 12 of the value of the constant 1, a signal adder 13, a signal scaler 14 and the signal multiplier 10. At the signal adder 13 output is 1−x where the signal further enters the signal scaler 14 that splits the signal x÷(1−x) which the signal multiplier 10 multiplies by x and
is obtained. The total transfer function of the circuit on the
By adjusting the values a and b we can obtain any function
where c and d are arbitrary positive values.
and amplifier/attenuator 5 after the non-linear element 4 having positive value b. The hyperbolic non-linear element 4 is derived from the source 12 of the value of the constant 1, the signal adder 13, the signal scaler 14, the signal multiplier 10 and the signal inverter 11. At the signal adder 13 output is
where the signal further enters the signal scaler 14 that splits the signal x÷(1+x) which the signal multiplier 10 multiplies by x and
is obtained. The total transfer function of the circuit on the
By adjusting the values a and b we can obtain any function
where c and d are arbitrary positive values.
The audio signal processing method in audio chain as illustrated in
The other embodiment of the apparatus 19 and of the method within audio chain is illustrated in
The audio signal processing method in audio chain shown on the
According to the method of the present invention, the apparatus 19 reduces by two times the non-linearity of the approximated psychoacoustical characteristic of the human ear and/or by 3 times the quadratic non-linearity of the electroacoustic transducer, and the pressure change by the human ear up to pΔ=±1 Pa.
Furthermore, according to the method of the present invention, the audio signal can be processed either in an analogue format or in a digital format.
The present invention relates also to a computer program adapted to run on a processor and to perform the method steps according to the present invention when carried out on a computer device.
The outputs of three branches of non-linear parts of the signal 7 enter the adder 8 made of a resistor network that converts the non-linear output voltage signals 7, as well as the audio signal after the input stage 23, into a sum of currents that make up the output audio signal 9, where the output inverting stage 26 converts them into the output voltage 9a.
The inverse psychoacoustic feature of the human ear can be approximated also by other functions and derivations of the non-linear element 4 can be performed by applying non-linearities of electronic elements such as diodes, transistors and vacuum tubes.
by non-linearity x1.5, which corresponds to the current/voltage characteristic of the vacuum diode I=k·U1.5. The approximation on the
In order to see how the non-linearity of the human ear decreases, in the approximate psychoacoustic feature of the human ear
we replace x with
and obtain the first five members:
we reduced at least two times the non-linearities introduced by the x2, x3 and x4 members of the approximated psychoacoustic characteristics of the human ear.
The implementation of the non-linear element 4 by applying vacuum diodes is illustrated in
Number | Date | Country | Kind |
---|---|---|---|
P20190292A | Feb 2019 | HR | national |
P20191903A | Oct 2019 | HR | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/HR2019/000027 | 10/22/2019 | WO |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2020/165616 | 8/20/2020 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
20060133620 | Lashkari | Jun 2006 | A1 |
20110257982 | Smithers | Oct 2011 | A1 |
20150146890 | Yang | May 2015 | A1 |
20170061982 | Pakarinen | Mar 2017 | A1 |
Number | Date | Country |
---|---|---|
0841836 | May 1998 | EP |
2018075442 | Apr 2018 | WO |
Entry |
---|
International Preliminary Report on Patentability dated Oct. 3, 2021 issued in International Patent Application PCT/HR2019/000027. |
International Search Report and Written Opinion of the Searching Authority dated Feb. 17, 2020 issued in International Patent Application PCT/HR2019/000027. |
Klippel, Wolfgang. “The Mirror Filter—A New Basis for Reducing Nonlinear Distortion and Equalizing Response in Woofer Systems”, J. Audio Eng. Soc., vol. 40, No. 9, pp. 675-691, (Sep. 30, 1992). |
Klippel, Wolfgang. “Active Compensation of Transducer Nonlinearities”, AES 23rd International Conference, Copenhagen, Denmark, May 23-25, 2003, pp. 1-17. |
Wee-Tong, Lim et al, “Synthesis of Polynomial-Based Nonlineal Device and Harmonic Shifting Technique for Virtual Bass System”, Circuits and Systems, 2009. ISCAS 2009. IEEE International Symposium on, IEEE, Piscataway, NJ, USA, May 24, 2009, pp. 1871-1874. |
Number | Date | Country | |
---|---|---|---|
20220141610 A1 | May 2022 | US |