CN102970647B - Simulating calculation method for nonlinear characteristics in loudspeaker vibration - Google Patents
Simulating calculation method for nonlinear characteristics in loudspeaker vibration Download PDFInfo
- Publication number
- CN102970647B CN102970647B CN201210462713.0A CN201210462713A CN102970647B CN 102970647 B CN102970647 B CN 102970647B CN 201210462713 A CN201210462713 A CN 201210462713A CN 102970647 B CN102970647 B CN 102970647B
- Authority
- CN
- China
- Prior art keywords
- voice coil
- function
- frequency
- loudspeaker voice
- nonlinear
- 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
Landscapes
- Audible-Bandwidth Dynamoelectric Transducers Other Than Pickups (AREA)
Abstract
The invention proposes a simulating calculation method for nonlinear characteristics in loudspeaker vibration. According to a relation that physical quantities such as a stiffness factor of a loudspeaker vibration system, an average magnetic induction intensity at a voice coil and an electrical inductance of the voice coil are changed with the position changes of the voice coil, each frequency component of voice coil vibration under the excitation of any single-frequency or multiple-frequency voltages is calculated, thereby obtaining nonlinear characteristics such as direct-current bias and harmonic distortion and intermodulation distortion in voice coil vibration. By the method disclosed by the invention, various linear characteristics in voice coil vibration and various linear and nonlinear characteristics of voice coil currents also can be obtained.
Description
Technical field
The invention belongs to loud speaker field, relate to the emulated computation method of nonlinear characteristic in a kind of loudspeaker vibration.Use the present invention also can obtain the multiple linear characteristic in voice coil vibrations simultaneously, the multiple line style of voice coil loudspeaker voice coil electric current and nonlinear characteristic.
Background technology
Loud speaker, when relatively high power works, because the physical quantity such as average magnetic induction intensity (B) and voice coil loudspeaker voice coil inductance (L) at vibrational system stiffness factor (K), voice coil loudspeaker voice coil place (x) changes with voice coil loudspeaker voice coil displacement, causes the Non-Linear Vibration of vibrational system.The Non-Linear Vibration of vibrational system causes the nonlinear distortion of loud speaker, shows sound impure, unpleasant, therefore studies and suppresses the nonlinear distortion in loudspeaker vibration to be main points in medium-to-high grade loudspeaker design.Loudspeaker vibration parts area is larger, motion state changes to some extent with position, for the ease of quantitative research, can the Vibration Condition at primary study voice coil loudspeaker voice coil place, obtain the various nonlinear characteristics in voice coil vibrations, such as: displacement bias (the balance point place of average displacement not necessarily when static of voice coil loudspeaker voice coil steady-state vibration, the distance departed from is displacement bias), each harmonic distortion is (when single frequency excitation, voice coil vibrations not only comprises driving frequency, also 2 overtones bands of driving frequency are comprised, 3 overtones bands etc., this phenomenon is called harmonic distortion), each intermodulation distortion is (when double frequency excitation, voice coil vibrations not only comprises driving frequency, also comprise driving frequency with frequently and difference frequency etc., this phenomenon is called intermodulation distortion) etc.Although the nonlinear characteristic of voice coil loudspeaker voice coil place vibration is not equal to the nonlinear characteristic of loud speaker entirety, have directly close contact, voice coil vibrations non-linear is the main source of Nonlinear Distortion of Loudspeakers.Understood the nonlinear vibration characteristics of loudspeaker voice coil by simulation calculation, effectively can instruct design and the improvement of loud speaker.
Existing correlation technique and method:
(1) method of sample preparation measurement.According to designing and producing a sample, measure the characteristic such as harmonic distortion, intermodulation distortion of sample, the distortion situation understanding design in this way is generally adopted by current electroacoustic enterprise.The shortcoming of this method has: 1. owing to needing sample preparation to each design improvement, therefore the cycle grows up to this height; 2. every detection is more consuming time, the detection of especially multi-frequency intermodulation distortion, and various combination of frequency total amount is difficult to complete detection very greatly; 3. detect and can find distortion situation, but can not illustrate it is the nonlinear distortion which partial design causes;
(2) paper " Tutorial:Loudspeaker Nonlinearities-Causes, Parameters, Symptoms " of Wolfgang Klippel is for the comprehensive elaboration that caused the principle of nonlinear distortion and feature to do of each parts.Based on this cover theory, Klippel measuring system can provide analysis and proposition qualitatively to the generation of the nonlinear distortion of loud speaker.But this cover system just rule of thumb provides the reason of some possibilities " producing 10% nonlinear distortion ", and cannot analyze concrete nonlinear distortion situation quantitatively;
(3) comprise circuit, magnetic circuit and vibrational system Holistic modeling to loud speaker, utilize nonlinear finite element method to calculate transient motion characteristic, recurrence calculation can obtain vibrational system nonlinear motion curve.The shortcoming of this method has: 1. model is large, and relate to multiple physical field coupling and calculate and dynamic calculation, need very large memory, computing time, overlength, needed special computing cluster just may realize; 2. the algorithm of multistep transient superposition easily produces error accumulation; 3. pair user's physics, mathematics standard require high, are unfavorable for that electroacoustic engineering teacher grasps.
Summary of the invention
The object of the invention is for the problems referred to above, propose physical quantity and the voice coil loudspeaker voice coil displacement functional relations (x) such as a kind of average magnetic induction intensity (B) according to known stiffness factor (K), voice coil loudspeaker voice coil place and voice coil loudspeaker voice coil inductance (L), calculate the method for the nonlinear characteristic under any single-frequency or multifrequency voltage drive in loudspeaker voice coil vibration.Utilize this method can calculate multiple nonlinear characteristic in the vibration of under the coupling of various predominant non-linearity factor loudspeaker voice coil and in voice coil loudspeaker voice coil electric current fast.
The object of the present invention is achieved like this:
The emulated computation method of nonlinear characteristic in loudspeaker vibration of the present invention, according to each component physical character of loud speaker, calculates under single-frequency or multifrequency voltage drive, the nonlinear characteristic in loudspeaker voice coil vibration.The method comprises the following steps:
(1) to the voltage drive function of any period, its frequency content and the range value of each frequency is determined;
(2) according to above-mentioned each frequency values, calculate each frequency multiplication within the scope of concern and frequently and the frequency values of difference frequency, remove and wherein repeat or too close part, obtain calculated rate set, within the scope of concern, generally comprise 20Hz to 5kHz, frequency sets is designated as:
;
(3) loudspeaker voice coil displacement function and voice coil loudspeaker voice coil current function are expressed as the form of the weighted sum of one group of basic function:
Coefficient
for unknown quantity, basic function set is:
; Wherein
it is the calculated rate set obtained in step (2);
(4) utilize the result of step (3), by time differentiate, obtain the time-derivative of voice coil loudspeaker voice coil speed, voice coil loudspeaker voice coil acceleration, voice coil loudspeaker voice coil electric current:
They are all be the function of independent variable with the two groups of unknown quantitys mentioned in step (3) and time;
(5) utilize the result of step (3), it is the function of independent variable that vibrational system stiffness factor (K), voice coil loudspeaker voice coil place average magnetic induction intensity (B) and voice coil loudspeaker voice coil inductance (L) physical quantity are write as with the two groups of unknown quantitys mentioned in step (3) and time:
；
(6) function in step (3), step (4), step (5) is substituted into nonlinear loudspeaker circuit equation and Nonlinear Equations of Motion, forming two with the two groups of unknown quantitys mentioned in step (3) and time is the differential equation of independent variable;
(7) two equations in step (6) to be multiplied with one group of trial function and to time variable integration, not to be comprised the time, about the equation group of the two groups of unknown quantitys mentioned in step (3); Wherein time integral interval should be taken as the common multiple in the cycle of each single-frequency obtained in step (2), or gets the value that a low-limit frequency obtained in step (2) do not comprise the cycle of 0;
(8) by the equation group obtained in Solving Nonlinear Systems of Equations Algorithm for Solving step (7);
(9) according to the voice coil loudspeaker voice coil displacement function of hypothesis in the two groups of coefficient values solved and step (3) and the form of voice coil loudspeaker voice coil current function, the last solution of voice coil loudspeaker voice coil displacement function and voice coil loudspeaker voice coil current function is obtained; The multiple linear and nonlinear characteristic in speaker loop current and voice coil vibrations is also calculated according to these two groups of coefficient values;
(10) choose the voltage excitation signals of multiple different frequency, repeat step (1) to the process of step (9), to obtain in voice coil loudspeaker voice coil electric current and voice coil vibrations various linear processes characteristic with the curve of frequency change.
Described loud speaker is coil-moving speaker.
Basic function set described in step (3) is respectively the SIN function of each frequency values obtained in step (2) by a constant function and a series of frequency, cosine function forms.
Nonlinear loudspeaker circuit equation described in step (6) includes but not limited to:
(1) coil is modeled as inductance series resistance, the basic circuit equation of formation:
Wherein t represents the time, and u (t) is driving voltage function, i (t) voice coil loudspeaker voice coil electric current, x (t) voice coil loudspeaker voice coil displacement function, constant Re voice coil loudspeaker voice coil D.C. resistance, constant l voice coil loudspeaker voice coil length;
(2) based on basic circuit equation, magnetic circuit current vortex phenomenon is modeled as inductance in parallel resistance, the circuit equation comprising eddy current of formation:
;
Wherein L2 and R2 is analog electrical vortex phenomenon device parameters, is respectively voice coil loudspeaker voice coil inductance and resistance, and i2 (t) is the voice coil loudspeaker voice coil electric current by resistance R2;
Separate this equation to need to add one group of unknown quantity, simultaneously many prescription journeys, copy the mode of process voice coil loudspeaker voice coil electric current to process;
Based on basic circuit equation, consider that voice coil loudspeaker voice coil electric current affects to some extent on magnetic flux density, the circuit equation of formation:
;
Nonlinear Equations of Motion described in step (6) includes but not limited to:
Be spend the spring of coefficient with voice coil loudspeaker voice coil change in displacement without a break by vibrating system modeling, form the equation of motion:
;
Consider that voice coil loudspeaker voice coil electric current affects to some extent on magnetic flux density, the equation of motion of formation:
；
Wherein Rms is the system mechanics damping in vibration, and m is vibrational system effective mass;
One group of trial function described in step (7) is elected as and is included but not limited to: A) identical with basic function set described in step (3); B) one group of unit pulse function, also referred to as: Dirac delta function, Dirac Delta function; C) each trial function is 1 in a segmentation, is 0 outside this segmentation;
Solving Nonlinear Systems of Equations algorithm described in step (8) includes but not limited to: Picard iterative method Picard iteration, confidence domain method Trust Region Method, the derivation algorithm of these algorithms;
Linear characteristic described in step (9) includes but not limited to: impedance curve comprises amplitude and phase place, circuit transfer function, voice coil loudspeaker voice coil Displacement Transfer Function; Nonlinear characteristic described in step (9) at least comprises that the Nonlinear Steady response of voice coil vibrations, vibration are biased, vibration peak, each order harmonics distortion, total harmonic distortion, the intermodulation distortion of each rank, overall intermodulation distortion, amplitude modulation distortion, frequency modulation distortion, and the Nonlinear Steady response of voice coil loudspeaker voice coil electric current, current offset, current peak, each order harmonics distortion, total harmonic distortion, the intermodulation distortion of each rank, overall intermodulation distortion, amplitude modulation distortion, frequency modulation distortion.
Owing to present invention employs above-mentioned technical scheme, therefore have the following advantages:
(1) contain the principal element causing nonlinear loudspeaker to move in calculating, the nonlinear situation that loud speaker shows under concrete excitation can be calculated very accurately.
(2) amount of calculation is very little, takies resource short for few computing time.
(3) not only can calculate the loudspeaker movements state under single frequency excitation, the motion state of loud speaker under double frequency or more complicated pumping signal can also be calculated.
(4) to various nonlinear physics amount, only need the functional relation determined, do not need expression formula that is clear and definite or special shape.
(5) obtain the solution of the equation of motion and circuit equation, therefore this method also can be used for the linear processes characteristic emulating speaker circuit part simultaneously, such as, emulate impedance curve (amplitude and phase place), current harmonics distortion, electric current intermodulation distortion etc. simultaneously.
(6) based on weighted residual method and confidential interval solution, obtain semi-analytical solution, computational accuracy is high.As long as select enough harmonic numbers just can obtain result accurately.
(7) this method directly obtains the steady state solution of motion and circuit, is not subject to the impact of transient state factor, does not also need time recursion to so there is no the problem of error accumulation.
Embodiment
Below the embodiment of this patent is elaborated.
(1) functional relation first needing the physical quantitys such as known vibrational system stiffness factor (K), voice coil loudspeaker voice coil place average magnetic induction intensity (B) and voice coil loudspeaker voice coil inductance (L) (x) to change with voice coil position, these functions are write respectively and are done: K (x), B(x) and L(x).Also need some the constant value physical quantitys knowing this loud speaker in addition, comprising: voice coil loudspeaker voice coil D.C. resistance
, voice coil loudspeaker voice coil length (l), this amount also often calculates together with B, is denoted as Bl), vibrational system mechanical damping
with vibrational system effective mass (m).These functional relations and constant can being calculated from design by methods such as finite elements, also can obtain these functional relations and constant by making sample measurement.Emulation or detect the form that the physical quantity that obtains is all generally discrete point with the variation relation of position, can fit to polynomial form and use, the method matching of spline curve and so on also can be adopted to use, and the form of this method to these functions does not limit.
(2) according to demand, determine a kind of voltage excitation signals that will calculate, this signal can be simple signal, and also can be the multiple-frequency signal of multiple simple signal combination, the frequency wherein comprised be designated as:
.
(3) displacement function of voice coil loudspeaker voice coil comprises each rank frequency multiplication and frequently and difference frequency component of each frequency of pumping signal.That is, to arbitrarily
and arbitrary integer
, the single-frequency components that in voice coil loudspeaker voice coil displacement function and voice coil loudspeaker voice coil current function may be containing frequency.Only need in Practical Calculation to get with all less, and the component (such as, getting 20Hz to 10kHz) within the scope of concern calculates.Remove the frequency too close with some repeated in said frequencies, the frequency set finally determined arranges from small to large and is designated as.
In order to follow-up use, also should determine that a time integral interval is 0 to T.T should be taken as the common multiple in the cycle of above-mentioned single-frequency, or gets a value much larger than the cycle of the low-limit frequency (not comprising 0) in said frequencies.
(4), after determining the frequency that may contain in voice coil loudspeaker voice coil displacement function and voice coil loudspeaker voice coil current function, voice coil loudspeaker voice coil displacement function can be write and do:
voice coil loudspeaker voice coil current function can be write as:
, wherein
undetermined coefficient, 4n+2 altogether.If row vector
Be denoted as w(t), order
, voice coil loudspeaker voice coil displacement function and current function can be write as simpler and clearer vector form:
(5) can obtain time differentiate above equation both sides:
Wherein due to w(t) be that known function (works as frequency set
after determining), therefore its time-derivative and second derivative are all the known functions that can write out.
(6) according to the result of step (4), the known conditions prepared in integrating step (1), K, B and L can be write as the function of X and time t:
。Generally speaking these functions are all nonlinear functions, and are often difficult to clearly write out expression formula, but the value of any given a group, just can draw these physical quantitys value at any time.
(7) for the basic nonlinear circuit equation of loud speaker and basic Nonlinear Equations of Motion:
It is wherein driving voltage function.Bring step (4) into above-mentioned two formulas to the equation in step (6) to obtain:
(7) two equations are all set up random time above, and therefore getting arbitrarily a function of time (being called " trial function "), is all 0 with the inner product (inner product of two functions of time is: two functions are multiplied, then quadratures on 0 to T to the time) of equation levoform.Such as, can amount of orientation w(t) in 2n+1 the function of time as trial function, ask inner product with above two equation levoforms respectively, do not comprised time variable, about common 4n+2 the equation of X and I.With matrix form write out into:
Because K (X, t), B (X, t) and L (X, t) not necessarily have simple expression formula, therefore these equations cannot direct solution.But, as long as there is the value of one group of X and I, just can know the functional value of each time point in integration item, use numerical integration method, such as Gauss – Legendre integral formula, just can calculate the approximation of each integration item, can show that above-mentioned equation group brings 4n+2 the margin value in left side after specific X and I into further.
(8) use Solving Nonlinear Systems of Equations algorithm, such as Picard iterative method (Picard iteration), confidence domain method (Trust Region Method) etc., can solve the equation group obtained in step (7).Wherein attempt each time needing surplus on the left of the middle equation group of calculation procedure (7).Specific algorithm is with reference to correlative theses.
(9), after solving coefficient X and I, the steady state solution of voice coil loudspeaker voice coil displacement and voice coil loudspeaker voice coil electric current is just obtained according to the equation in step (4).According to these coefficients and pumping signal information, multiple linear and nonlinear parameter can also be calculated.Such as, if pumping signal U (t) is simple signal, this Frequency and Amplitude is U, can obtain from X and I vector:
1, the direct current biasing of this frequency place vibration is a
0.
2, the fundamental component of this frequency place vibration is:
, the transfer function amplitude being activated to vibration is:
.
3, the n order harmonics distortion ratio of this frequency place vibration is:
, wherein
it is gross energy.
4, the total harmonic distortion ratio of this frequency place vibration is:
.
5, this frequency place loudspeaker impedance amplitude is:
, phase place is:
.
Each order harmonics distortion in electric current and total harmonic distortion etc. also can be obtained by the formula in similar 3 and 4.
If pumping signal is two-frequency signal, can also obtain n rank intermodulation distortion ratio is:
The frequency that wherein in denominator, coefficient is corresponding is one of the fundamental frequency in pumping signal (organize dual-frequency data if calculated, this frequency should be " fixing " frequency wherein more), is called f
1, another fundamental frequency in pumping signal is called f
2, then two frequencies that in molecule, coefficient is corresponding respectively:
with
.Similarly can also obtain overall intermodulation distortion etc.
In conjunction with the steady state solution of voice coil loudspeaker voice coil displacement, voice coil loudspeaker voice coil electric current and coefficient, displacement/current peak, amplitude modulation distortion, frequency modulation distortion etc. can also be obtained, various nonlinear loudspeaker characteristic quantity.
(10) generally can to the operation of a series of frequency repeated application this method step (2) to step (9), calculate a series of nonlinear characteristic value, obtain the relation of these characteristic values with frequency change, such as, use this method can obtain the vibration transfer function curve, impedance curve, direct current biasing curve, each order harmonics distortion curve, total harmonic distortion curve etc. of loud speaker to a series of simple signal.To a series of two-frequency signal (one of them frequency is constant, is called " fixing " frequency), each rank intermodulation distortion curve, overall intermodulation distortion curve etc. can be calculated.Utilize these characteristic curves, loudspeaker design and improvement can be instructed.
Above embodiment is only in order to illustrate the present invention and and unrestricted technical scheme described in the invention.Although this specification with reference to each embodiment above-mentioned to invention has been detailed description, but will be understood by those skilled in the art that, when not departing from the spirit and scope limited of following claims, many amendments, change or equivalence can be made and replace; And all do not depart from technical scheme and the improvement thereof of the spirit and scope of the present invention, it all should be covered by the scope of the present invention.
Claims (8)
1. the emulated computation method of nonlinear characteristic in loudspeaker vibration, according to each component physical character of loud speaker, calculates under single-frequency or multifrequency voltage drive, and the nonlinear characteristic in loudspeaker voice coil vibration, is characterized in that the method comprises the following steps:
(1) to the voltage drive function of any period, its frequency content and the range value of each frequency is determined;
(2) according to above-mentioned each frequency values, calculate each frequency multiplication within the scope of concern and frequently and the frequency values of difference frequency, remove and wherein repeat or too close part, obtain calculated rate set, generally comprise 20Hz to 5kHz within the scope of concern, frequency sets is designated as:
;
(3) loudspeaker voice coil displacement function and voice coil loudspeaker voice coil current function are expressed as the form of the weighted sum of one group of basic function:
Coefficient
for unknown quantity, basic function set is:
; Wherein
it is the calculated rate set obtained in step (2);
(4) utilize the result of step (3), by time differentiate, obtain the time-derivative of voice coil loudspeaker voice coil speed, voice coil loudspeaker voice coil acceleration, voice coil loudspeaker voice coil electric current:
They are all be the function of independent variable with the two groups of unknown quantitys mentioned in step (3) and time;
(5) utilize the result of step (3), it is the function of independent variable that vibrational system stiffness factor K, voice coil loudspeaker voice coil place average magnetic induction intensity B and voice coil loudspeaker voice coil inductance L physical quantity are write as with the two groups of unknown quantitys mentioned in step (3) and time:
(6) function in step (3), step (4), step (5) is substituted into the basic nonlinear circuit equation of loud speaker and basic Nonlinear Equations of Motion, forming two with the two groups of unknown quantitys mentioned in step (3) and time is the differential equation of independent variable;
(7) two differential equations in step (6) to be multiplied with one group of trial function and to time variable integration, not to be comprised the time, about the equation group of the two groups of unknown quantitys mentioned in step (3); Wherein time integral interval should be taken as the common multiple in the cycle of each single-frequency obtained in step (2), or gets the value that a low-limit frequency obtained in step (2) do not comprise the cycle of 0;
(8) by the equation group obtained in Solving Nonlinear Systems of Equations Algorithm for Solving step (7);
(9) according to the voice coil loudspeaker voice coil displacement function of hypothesis in the two groups of coefficient values solved and step (3) and the form of voice coil loudspeaker voice coil current function, the last solution of voice coil loudspeaker voice coil displacement function and voice coil loudspeaker voice coil current function is obtained; The multiple linear and nonlinear characteristic in speaker loop current and voice coil vibrations is also calculated according to these two groups of coefficient values;
(10) choose the voltage excitation signals of multiple different frequency, repeat step (1) to the process of step (9), to obtain in voice coil loudspeaker voice coil electric current and voice coil vibrations various linear processes characteristic with the curve of frequency change.
2. the emulated computation method of nonlinear characteristic in a kind of loudspeaker vibration according to claim 1, is characterized in that: described loud speaker is coil-moving speaker.
3. the emulated computation method of nonlinear characteristic in a kind of loudspeaker vibration according to claim 1, is characterized in that: the basic function set described in step (3) is respectively the SIN function of each frequency values obtained in step (2) by a constant function and a series of frequency, cosine function forms.
4. the emulated computation method of nonlinear characteristic in a kind of loudspeaker vibration according to claim 1, is characterized in that: the nonlinear loudspeaker circuit equation described in step (6) includes but not limited to:
(1) coil is modeled as inductance series resistance, the basic circuit equation of formation:
Wherein t represents the time, and u (t) is driving voltage function, and i (t) is voice coil loudspeaker voice coil electric current, and x (t) is voice coil loudspeaker voice coil displacement function, and constant Re is voice coil loudspeaker voice coil D.C. resistance, constant
for voice coil loudspeaker voice coil length;
(2) based on basic circuit equation, magnetic circuit current vortex phenomenon is modeled as inductance in parallel resistance, the circuit equation comprising eddy current of formation:
;
Wherein L2 and R2 is analog electrical vortex phenomenon device parameters, is respectively voice coil loudspeaker voice coil inductance and resistance, and i2 (t) is the voice coil loudspeaker voice coil electric current by resistance R2;
Separate this equation to need to add one group of unknown quantity, simultaneously many prescription journeys, copy the mode of process voice coil loudspeaker voice coil electric current to process;
(3) based on basic circuit equation, consider that voice coil loudspeaker voice coil electric current affects to some extent on magnetic flux density, the circuit equation of formation:
。
5. the emulated computation method of nonlinear characteristic in a kind of loudspeaker vibration according to claim 1, is characterized in that: the Nonlinear Equations of Motion described in step (6) includes but not limited to:
(1) be spend the spring of coefficient with voice coil loudspeaker voice coil change in displacement without a break by vibrating system modeling, form the equation of motion:
；
(2) consider that voice coil loudspeaker voice coil electric current affects to some extent on magnetic flux density, the equation of motion of formation:
；
Wherein Rms is the system mechanics damping in vibration, and m is vibrational system effective mass.
6. the emulated computation method of nonlinear characteristic in a kind of loudspeaker vibration according to claim 1, is characterized in that: one group of trial function described in step (7) is elected as and included but not limited to: A) identical with basic function set described in step (3); B) one group of unit pulse function, also referred to as: Dirac delta function, Dirac Delta function; C) each trial function is 1 in a segmentation, is 0 outside this segmentation.
7. the emulated computation method of nonlinear characteristic in a kind of loudspeaker vibration according to claim 1, it is characterized in that: the Solving Nonlinear Systems of Equations algorithm described in step (8) includes but not limited to: Picard iterative method Picard iteration, confidence domain method Trust Region Method, the derivation algorithm of these algorithms.
8. the emulated computation method of nonlinear characteristic in a kind of loudspeaker vibration according to claim 1, is characterized in that: the linear characteristic described in step (9) includes but not limited to: impedance curve comprises amplitude and phase place, circuit transfer function, voice coil loudspeaker voice coil Displacement Transfer Function; Nonlinear characteristic described in step (9) at least comprises that the Nonlinear Steady response of voice coil vibrations, vibration are biased, vibration peak, each order harmonics distortion, total harmonic distortion, the intermodulation distortion of each rank, overall intermodulation distortion, amplitude modulation distortion, frequency modulation distortion, and the Nonlinear Steady response of voice coil loudspeaker voice coil electric current, current offset, current peak, each order harmonics distortion, total harmonic distortion, the intermodulation distortion of each rank, overall intermodulation distortion, amplitude modulation distortion, frequency modulation distortion.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210462713.0A CN102970647B (en) | 2012-11-16 | 2012-11-16 | Simulating calculation method for nonlinear characteristics in loudspeaker vibration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210462713.0A CN102970647B (en) | 2012-11-16 | 2012-11-16 | Simulating calculation method for nonlinear characteristics in loudspeaker vibration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102970647A CN102970647A (en) | 2013-03-13 |
| CN102970647B true CN102970647B (en) | 2015-04-01 |
Family
ID=47800453
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210462713.0A Active CN102970647B (en) | 2012-11-16 | 2012-11-16 | Simulating calculation method for nonlinear characteristics in loudspeaker vibration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102970647B (en) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105246018B (en) * | 2015-10-29 | 2019-02-26 | 浙江中科电声研发中心 | The approximate calculation method of limited big baffle axial symmetry harmonic distortion of loudspeaker |
| CN105611473A (en) * | 2016-03-28 | 2016-05-25 | 歌尔声学股份有限公司 | Diaphragm-voice coil assembly and system and method for adjusting balance of vibrating voice coil |
| CN106068007B (en) * | 2016-06-07 | 2019-09-13 | 瑞声科技(新加坡)有限公司 | Nonlinear loudspeaker system identifying method |
| EP3476135A1 (en) * | 2016-06-24 | 2019-05-01 | Harman Becker Automotive Systems GmbH | Method for simulating total harmonic distortion of a loudspeaker |
| CN106060730B (en) * | 2016-06-29 | 2019-05-03 | 歌尔股份有限公司 | The method for reducing super-thin speaker distortion |
| CN108540242B (en) * | 2018-01-05 | 2021-06-25 | 中国传媒大学广州研究院 | Method and device for measuring harmonic distortion index of broadcast transmitter |
| CN109145514A (en) * | 2018-09-30 | 2019-01-04 | 浙江中科电声研发中心 | A kind of Numerical Simulation Analysis method of loudspeaker distortions |
| CN110225433B (en) * | 2019-05-16 | 2021-04-13 | 音王电声股份有限公司 | Nonlinear measurement and tone quality tuning method of loudspeaker system |
| CN111551848B (en) * | 2019-12-30 | 2022-09-02 | 瑞声科技(新加坡)有限公司 | Motor experience distortion index testing method, electronic equipment and storage medium |
| CN113286233B (en) * | 2021-05-20 | 2022-12-27 | 维沃移动通信有限公司 | Loudspeaker simulation method, loudspeaker simulation device, loudspeaker and electronic equipment |
| CN117412226B (en) * | 2023-02-21 | 2024-08-23 | 东莞市惟动智能科技有限公司 | Moving coil vibrator with parallel magnetic coils and push-pull nonlinear cancellation function and application thereof |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102118678A (en) * | 2011-04-02 | 2011-07-06 | 嘉兴中科声学科技有限公司 | Method and system for measuring speaker parameters by using current sensor |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20050023841A (en) * | 2003-09-03 | 2005-03-10 | 삼성전자주식회사 | Device and method of reducing nonlinear distortion |
-
2012
- 2012-11-16 CN CN201210462713.0A patent/CN102970647B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102118678A (en) * | 2011-04-02 | 2011-07-06 | 嘉兴中科声学科技有限公司 | Method and system for measuring speaker parameters by using current sensor |
Non-Patent Citations (3)
| Title |
|---|
| M. Faifer,."A Novel, Cost-Effective Method for Loudspeakers Parameters Measurement".《IEEE AUTOTESTCON 2008》.2008,第1-7页. * |
| Walter Frank,."Loudspeaker Nonlinearities- Analysis and Compensation".《Signals, Systems and Computers, 1992. 1992 Conference Record of The Twenty-Sixth Asilomar Conference on》.1992,参见第756-760页. * |
| 陆晓,."动圈式扬声器数值分析方法".《声学技术》.2012,第 31 卷(第 4 期),参见第407-412页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102970647A (en) | 2013-03-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102970647B (en) | Simulating calculation method for nonlinear characteristics in loudspeaker vibration | |
| Tronchin | The emulation of nonlinear time-invariant audio systems with memory by means of Volterra series | |
| Yeh | Digital implementation of musical distortion circuits by analysis and simulation | |
| CN102843633A (en) | Control of loudspeaker output | |
| CN113286233B (en) | Loudspeaker simulation method, loudspeaker simulation device, loudspeaker and electronic equipment | |
| CN110427650A (en) | The Numerical Simulation Analysis method of moving iron type loud-speaker fundamental characteristics | |
| CN110907827B (en) | Motor transient distortion measuring method and system | |
| Novak et al. | Nonparametric identification of nonlinear systems in series | |
| King et al. | Fractional derivative loudspeaker models for nonlinear suspensions and voice coils | |
| Saraswat et al. | Modeling and study of nonlinear effects in electrodynamic shakers | |
| Parker et al. | Physical modeling using recurrent neural networks with fast convolutional layers | |
| Massi et al. | Deep learning-based wave digital modeling of rate-dependent hysteretic nonlinearities for virtual analog applications | |
| Nielsen et al. | Estimation of optimal values for lumped elements in a finite element—Lumped parameter model of a loudspeaker | |
| TW201414991A (en) | Method for measuring electroacoustic parameters of transducer | |
| CN112487679A (en) | Method, system, terminal and medium for simulating acoustic characteristics of sound box by numerical values | |
| Zhang et al. | Parameter estimation of linear fractional-order system from laplace domain data | |
| Brunet et al. | State-space modeling and identification of loudspeaker with nonlinear distortion | |
| Middleton et al. | The application of Fourier neural operator networks for solving the 2D linear acoustic wave equation | |
| Novak | Measurement of loudspeaker parameters: A pedagogical approach | |
| Peng et al. | Finite linear representation of nonlinear structural dynamics using phase space embedding coordinate | |
| Holmes | Guitar effects: pedal emulation and identification | |
| Vanhatalo et al. | Neural Network-Based Virtual Analog Modeling | |
| Geddes et al. | Efficient loudspeaker linear and nonlinear parameter estimation | |
| Nemeth et al. | Measuring and simulating magnetic characteristics using epstein frame | |
| Munroe | Real Time Loudspeaker Control |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |