EP1665884A2 - Audio reproduction system - Google Patents
Audio reproduction systemInfo
- Publication number
- EP1665884A2 EP1665884A2 EP04780339A EP04780339A EP1665884A2 EP 1665884 A2 EP1665884 A2 EP 1665884A2 EP 04780339 A EP04780339 A EP 04780339A EP 04780339 A EP04780339 A EP 04780339A EP 1665884 A2 EP1665884 A2 EP 1665884A2
- Authority
- EP
- European Patent Office
- Prior art keywords
- signal
- transducer
- providing
- actuator
- compensator
- 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.)
- Withdrawn
Links
- 230000026683 transduction Effects 0.000 claims abstract description 70
- 238000010361 transduction Methods 0.000 claims abstract description 70
- 230000001419 dependent effect Effects 0.000 claims abstract description 65
- 230000033001 locomotion Effects 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims description 395
- 230000008569 process Effects 0.000 claims description 252
- 230000006870 function Effects 0.000 claims description 189
- 238000005259 measurement Methods 0.000 claims description 117
- 230000009021 linear effect Effects 0.000 claims description 72
- 230000005291 magnetic effect Effects 0.000 claims description 63
- 230000005236 sound signal Effects 0.000 claims description 54
- 230000004044 response Effects 0.000 claims description 45
- 230000036962 time dependent Effects 0.000 claims description 32
- 238000001514 detection method Methods 0.000 claims description 26
- 230000008878 coupling Effects 0.000 claims description 20
- 238000010168 coupling process Methods 0.000 claims description 20
- 238000005859 coupling reaction Methods 0.000 claims description 20
- 230000003287 optical effect Effects 0.000 claims description 12
- 239000002131 composite material Substances 0.000 claims description 7
- 238000012512 characterization method Methods 0.000 claims description 4
- 238000005183 dynamical system Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 description 109
- 230000000694 effects Effects 0.000 description 75
- 238000010586 diagram Methods 0.000 description 50
- 230000003750 conditioning effect Effects 0.000 description 36
- 230000001939 inductive effect Effects 0.000 description 33
- 239000000523 sample Substances 0.000 description 33
- HBBOZFUQJDYASD-LPHOMBEVSA-N alpha-L-Fucp-(1->3)-[beta-D-Galp-(1->4)]-beta-D-GlcpNAc Chemical compound O[C@H]1[C@H](O)[C@H](O)[C@H](C)O[C@H]1O[C@H]1[C@H](O[C@H]2[C@@H]([C@@H](O)[C@@H](O)[C@@H](CO)O2)O)[C@@H](CO)O[C@@H](O)[C@@H]1NC(C)=O HBBOZFUQJDYASD-LPHOMBEVSA-N 0.000 description 28
- 238000013459 approach Methods 0.000 description 26
- 239000000725 suspension Substances 0.000 description 24
- 230000003071 parasitic effect Effects 0.000 description 22
- 238000004088 simulation Methods 0.000 description 21
- 230000009467 reduction Effects 0.000 description 20
- 238000001228 spectrum Methods 0.000 description 19
- 230000008859 change Effects 0.000 description 17
- 238000005070 sampling Methods 0.000 description 16
- 230000002829 reductive effect Effects 0.000 description 14
- 230000003446 memory effect Effects 0.000 description 12
- 238000012360 testing method Methods 0.000 description 12
- 239000003570 air Substances 0.000 description 11
- 238000004422 calculation algorithm Methods 0.000 description 11
- 238000009826 distribution Methods 0.000 description 11
- 239000013598 vector Substances 0.000 description 11
- 238000013461 design Methods 0.000 description 10
- 230000004069 differentiation Effects 0.000 description 10
- 230000003068 static effect Effects 0.000 description 10
- 239000004065 semiconductor Substances 0.000 description 9
- 238000006243 chemical reaction Methods 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 8
- 230000003247 decreasing effect Effects 0.000 description 7
- 230000001934 delay Effects 0.000 description 7
- 238000009795 derivation Methods 0.000 description 7
- 230000001965 increasing effect Effects 0.000 description 7
- 230000009022 nonlinear effect Effects 0.000 description 7
- 230000036961 partial effect Effects 0.000 description 7
- 238000012546 transfer Methods 0.000 description 7
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 6
- 230000001133 acceleration Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 6
- 238000001914 filtration Methods 0.000 description 6
- 238000007493 shaping process Methods 0.000 description 6
- 238000004804 winding Methods 0.000 description 6
- 230000003321 amplification Effects 0.000 description 5
- 239000003990 capacitor Substances 0.000 description 5
- 238000003199 nucleic acid amplification method Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 230000000717 retained effect Effects 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 5
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 4
- 229910052782 aluminium Inorganic materials 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000006467 substitution reaction Methods 0.000 description 4
- 230000009897 systematic effect Effects 0.000 description 4
- 241000239290 Araneae Species 0.000 description 3
- 101100381996 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) BRO1 gene Proteins 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 3
- 238000009413 insulation Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 229910052742 iron Inorganic materials 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 229910001209 Low-carbon steel Inorganic materials 0.000 description 2
- 241000721701 Lynx Species 0.000 description 2
- 101100286571 Mus musculus Pyhin1 gene Proteins 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 230000003190 augmentative effect Effects 0.000 description 2
- 230000001364 causal effect Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 239000004020 conductor Substances 0.000 description 2
- 238000012885 constant function Methods 0.000 description 2
- 229910052802 copper Inorganic materials 0.000 description 2
- 239000010949 copper Substances 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000007429 general method Methods 0.000 description 2
- 230000001976 improved effect Effects 0.000 description 2
- 230000005415 magnetization Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000003121 nonmonotonic effect Effects 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- NLZUEZXRPGMBCV-UHFFFAOYSA-N Butylhydroxytoluene Chemical compound CC1=CC(C(C)(C)C)=C(O)C(C(C)(C)C)=C1 NLZUEZXRPGMBCV-UHFFFAOYSA-N 0.000 description 1
- 238000006842 Henry reaction Methods 0.000 description 1
- 230000005483 Hooke's law Effects 0.000 description 1
- 238000004566 IR spectroscopy Methods 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 239000012080 ambient air Substances 0.000 description 1
- 230000004323 axial length Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000005520 electrodynamics Effects 0.000 description 1
- 230000005288 electromagnetic effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003090 exacerbative effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000002847 impedance measurement Methods 0.000 description 1
- 238000007373 indentation Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000000696 magnetic material Substances 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000000873 masking effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910001172 neodymium magnet Inorganic materials 0.000 description 1
- 239000003973 paint Substances 0.000 description 1
- 230000005624 perturbation theories Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 229910052702 rhenium Inorganic materials 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 229910052720 vanadium Inorganic materials 0.000 description 1
- 238000013022 venting Methods 0.000 description 1
- 230000003245 working effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R29/00—Monitoring arrangements; Testing arrangements
- H04R29/001—Monitoring arrangements; Testing arrangements for loudspeakers
- H04R29/003—Monitoring arrangements; Testing arrangements for loudspeakers of the moving-coil type
Definitions
- FIG. 1 is a block diagram illustrating a typical audio reproduction system 100. As is seen in step 101 , an electrical audio signal, which may be digital or analog, is provided to a signal analysis shaping system 102.
- 5 signal analysis shaping system 102 is based on a speaker enclosure and a preference model. Thereafter, a modified version of the analog signal 103 is provided to a power switch or switches 104 that activate a transducer 105 contained in the speaker enclosure 106.
- a conventional speaker assembly there are generally a plurality of transducers which are typically voice coil transducers. Transducers are also commonly o referred to as drivers. However, many types of devices can be utilized as transducers in a speaker system.
- a conventional signal processing system also provides for standard audio amplification.
- Signal analysis shaping system 102 can be described functionally as illustrated in FIG. 2, which is a flow chart thereof for standard audio amplification.
- the input signal, 5 which may be in either analog or digital format, is provided to the signal processing system via step 201.
- the signal is adjusted to correct for speaker enclosure effects, via step 202. This may comprise correctional adjustments for frequency response due to resonances, anti-resonances and phase errors created in multi-transducer systems within speaker enclosures. o Conventional approaches may also include correctional adjustments of frequency response due to resonances, anti-resonances and phase errors arising from room and environmental distortions, which is accomplished in step 203. For example, adjustments may involve de-peaking of resonances to try to flatten the frequency response.
- the input signal is also adjusted for user preferences, in terms of frequency amplitude adjustment, which is accomplished in step 204.
- step 205 may be performed, in which the input signal may be adjusted for each transducer of the speaker system, for example, sending only the high frequency signal to the tweeter, and the low frequencies to the woofer or subwoofers. Following the completion of all correctional adjustments, the signal is sent to an output amplifier in step 206.
- a problem with the foregoing system is that there are frequency dependent errors as well as phase dependent errors which are not corrected, as well as errors due to the non-linear distortion of the transducer which reduce the effectiveness of the other corrections.
- FIG. 3 is an illustration of a typical voice coil transducer 300.
- the frame 301 holds the cone, or diaphragm 302.
- the diaphragm 302 is acted upon by voice coil 303 which acts as a motor, causing the diaphragm 302 to vibrate and create pressure waves in the ambient air.
- Voice coil 303 is comprised of a coil of wire wound around a tube or former. Voice coil 303 receives an electrical current, which is acted upon by the static magnetic field developed by the permanent magnet 304 and iron assembly 305 in the annular gap 306 in which voice coil 303 rides.
- the additional magnetic field from voice coil 303 which is induced by the external current driven through voice coil 303, interacts with the static magnetic field due to the permanent magnet 304 and iron assembly 305 within the annular gap 306, causing the voice coil 303 to move forward (toward the listener, to the right in FIG.
- a dome 309 acts as a dust cap and as a diffuser for high frequency sound.
- Nonlinear effects are an intrinsic part of the design of voice coil transducers. Nonlinearities in the motor factor in a voice coil transducer result from the fact that the coil and the region of uniform static magnetic field are limited in size, coupled with the fact that the coil moves relative to the static field. The actual size of the static magnetic field region, and its size relative to the voice coil, represent engineering and economic compromises. For a voice coil in a transducer, a stronger field results in a larger motor factor, and hence a larger motive force per given coil current magnitude. As the field falls off away from the annular gap 306, the motive force is reduced.
- the motive force per unit coil current is defined as the motor factor, and depends on the geometry of the coil and on the shape and position of the coil with respect to the static magnetic field configuration, the latter being generated by the permanent magnet or magnets and guided by the magnetic pole structures.
- This motor factor is usually denoted as the Bl factor, and is a function of x , the outward displacement of the coil/diaphragm assembly away from its equilibrium position (which the transducer relaxes to after the driving audio signal ceases).
- x is positive when the coil/diaphragm assembly is displaced from equilibrium in the direction of the listener, i.e. towards the front of the speaker.
- FIG. 4 represents data for actual large signal (LS) parameters of a transducer from a small desktop stereo system, model name: Spin70, manufactured by Labtec.
- the large signal parameters shown in FIG. 4 were obtained using a commercially available laser metrology system (Klippel GMBH).
- the magnitude of Bl is shown by curve 401 as a function of the displacement x of the coil/diaphragm assembly from the no-signal equilibrium position, which is indicated in FIG. 4 by a zero on the horizontal axis; at that position, no elastic restoring force is applied to the coil/diaphragm assembly.
- the unit for Bl is Newton I Ampere (or N I A). The highly non-constant nature of the Bl factors of commercial voice coil transducers is recognized in the current art.
- the cone suspension is axially symmetric and typically includes two parts: a corrugated suspension near the coil, typically referred to as the spider 307, and the surround 308 connecting the large end of cone 302 to the frame 301 of the speaker. These two suspensions together act as an effective spring, which provides a restoring force to the coil/diaphragm assembly and determines the equilibrium position of the assembly to which it relaxes when not being driven.
- This effective spring restoring force is again a highly non-constant function of coil/cone axial position x ; that is to say, the effective spring stiffness varies significantly as a function of x .
- curve 402 shows a plot of K , the spring stiffness, as a function of x for the speaker transducer mentioned above.
- Spring stiffness K is expressed in units of N I mm (i.e. Newton per millimeter).
- the mechanical equation of motion for the transducer can be approximated as a second order ODE (ordinary differential equation) in the position x of the coil/diaphragm assembly, treated as if it were a rigid piston.
- R ms represents the effective drag coefficient experienced by the assembly, mainly due to air back pressure and suspension friction
- K ⁇ x) is the position dependent effective spring stiffness due to the elastic suspension
- Bl(x) is the position dependent motor factor
- i(t) is the time dependent voice-coil current, which responds to the input audio signal and constitutes the control variable.
- x is used as the term for acceleration and x is used as the term for velocity.
- the second order differential equation (1 ) would be straightforward to solve, but for the nonlinearities in the elastic restoring force and in the motor force terms; these nonlinearities stem from the x dependence of K(x) and R (x) , and they preclude a closed-form analytical solution in the general case.
- the Ohmic resistance of the coil is R e .
- the coil's effective inductance, L e (x) is a function of x because it depends upon the instantaneous position of the coil relative to the magnetic pole structure and its airgap.
- curve 403 shows a typical plot of the position dependence of coil inductance L e (x) at low audio frequencies.
- the units of L e are mH (milli-Henries), and the values of L e shown in curve 403 have been multiplied by a factor 10 to render the graph more readable.
- Prior art includes a number of approaches for controlling the nonlinearities in audio transducers. These approaches include classic control methods based on negative feedback of a motional signal, as well as more recent methods based on system modeling and state estimation. It may seem apparent that a negative feedback system would be advantageous for reducing the nonlinear response of a voice coil transducer, and descriptions of several examples of such feedback systems do exist. Nevertheless, none of these prior techniques appear to have made any significant impact on commercial audio practice.
- Such feedback systems include ones based upon signals from microphones (U.S. Patent No. 6,122,385, U.S. Patent Application 2003/0072462A1 ), extra coils in the speakers (U.S. Patent Nos. 6,104,817, 4,335,274, 4,243,839, 3,530,244 and U.S. Patent Application 2003/0072462A1 ), piezoelectric accelerometers (U.S. Patent Application 2002/015906 A1, U.S. Patent Nos. 6,104,817, 5,588,065, 4,573,189) or back EMF (BEMF) (U.S. Patent Nos. 5,542,001 , 5,408,533).
- microphones U.S. Patent No. 6,122,385, U.S. Patent Application 2003/0072462A1
- extra coils in the speakers U.S. Patent Nos. 6,104,817, 4,335,274, 4,243,839, 3,530,244 and U.S. Patent Application 2003/0072462A1
- each step in the audio reproduction process is treated independently - by concentration on either amplifier design (drive), transducer design, or enclosure design - because there is little point in having a full-system control loop with such a large non-linear element, the transducer, running open-loop within the system. Accordingly, there are several factors described above that significantly affect the ability to provide accurate sound from a conventional audio reproduction system. Some of the issues can be addressed by improving the circuitry through digital means; but even with the digital circuitry to handle the signal shaping, the transducer itself has significant nonlinearities that can never be addressed adequately by shaping the input signal to the transducer. Therefore, what is needed is a system that controls the transducer in such a manner that optimum linear sound is provided. Such a system should also be easy to implement, cost effective, and easily adaptable to existing systems.
- the present invention provides a control system for a transducer to provide linear sound, and the present invention also provides an integrated audio reproduction system.
- a process for controlling an audio reproduction system that includes a sound transducer.
- the process includes preparing a model of a portion of the audio reproduction system and providing a control circuit having first and second inputs.
- the control circuit is configured as a function of the model.
- An audio signal is provided to the first input of the control circuit, and a position indication signal of the sound transducer is provided to a second input of the control circuit.
- the control circuit generates an output signal that is a function of both the position indication signal and the audio signal.
- the position indication signal is generated using an electrical characteristic of the system.
- the electrical characteristic is the impedance of a coil that is part of the sound transducer.
- the position indication signal is provided using an optical means.
- the optical technique involves directing an infrared light source to a portion of the sound transducer using an infrared light emitting diode.
- the electrical characteristic which is measured by the system to give an indication of a position involves measuring the capacitance of a coil of the sound transducer with respect to a structure of the sound transducer.
- preparing a model of a portion of the audio reproduction system comprises preparing a model of the sound transduction portion of the audio reproduction system.
- the model is prepared by determining at least one operational parameter of the speaker transducer as a function of a position of a coil and diaphragm of the speaker transducer with respect to another portion of the speaker transducer.
- the operational parameter can be an impedance of the coil, a motor factor of the coil and diaphragm, or the operational parameter may be a stiffness of a spring coupled to the diaphragm.
- the control circuit generates an output signal which compensates the system with respect to a spring stiffness of a spring support of a voice coil transducer.
- the control circuit can generate an output signal which compensates the system with respect to a motor factor of the voice coil transducer.
- a model may be made of the signal conditioning portion of the audio reproduction system.
- the control circuit generates an output signal for compensation of the system with respect to a back electromotive force of a driver of the sound transducer.
- the control circuit may be utilized to generate an output signal to compensate the system with respect to an impedance of the driver of the sound transducer.
- a process is provided for controlling an audio reproduction system including a sound transducer, which process comprises: providing a model of a portion of the audio reproduction system, providing a control circuit configured according to the model, and providing to the control circuit a signal indicative of a position of the sound transducer.
- the portion of the audio reproduction system which is modeled may be the signal conditioning portion of the system.
- the portion of the audio reproduction system modeled may be the sound transducer.
- the control circuit may be utilized to condition an audio signal as a function of a back electromotive force of a driver of the sound transducer. Additionally, the control circuit may be used to condition the audio signal as a function of an impedance of a driver of the sound transducer.
- the model of the audio reproduction system is a model of the sound transducer
- the control circuit may be utilized to condition an audio signal as a function of a spring stiffness of a spring support of the audio transducer.
- the control circuit may be utilized to condition the audio signal as a function of a motor factor of the coil and diaphragm assembly.
- FIG. 1 is a block diagram illustrating a typical audio reproduction system
- FIG. 2 is a flow chart depicting the functionality of a signal analysis shaping system
- FIG. 3 is an illustration of a typical voice coil transducer
- FIG. 4 graphically illustrates curves of Large Signal (LS) data for the actual parameters of the transducer from a Spin70 desktop stereo system manufactured by Labtec
- FIG. 5 illustrates relationships between the main areas of the present invention, grouped under three different headings: control systems, instrumentation, and audio reproduction
- FIG. 6 is a block diagram of an audio reproduction system in accordance with the three processes identified in the context of the present invention
- FIG. 1 is a block diagram illustrating a typical audio reproduction system
- FIG. 2 is a flow chart depicting the functionality of a signal analysis shaping system
- FIG. 3 is an illustration of a typical voice coil transducer
- FIG. 4 graphically illustrates curves of Large Signal (LS) data for the actual parameters of the transducer from a Spin70 desktop stereo system manufactured by Labtec
- FIG. 7 is a flow chart illustrating the process of feedback linearization in accordance with the present invention
- FIG. 8 is a block diagram of the main portion of a sound reproduction system, including a control system for controlling the operation of the sound reproduction system in accordance with the present invention
- FIG. 9 is a block diagram of the feedback linearization process using the control law of equation (34), which only linearizes the transduction component of the signal conditioning process, and without an electronically restored linear restoring force
- FIG. 10 is a block diagram of the feedback linearization process using the control law given by equation (40), which provides transduction correction along with a linear spring constant (suspension stiffness) that is electronically added;
- FIG. 40 linear spring constant
- FIG. 11 is a block diagram of the feedback linearization process for the control law correcting for spring, motor factor and BEMF nonlinearities, including an electronically restored linear spring and an electronically restored contribution to the linear drag force term;
- FIG. 12 is a block diagram of the feedback linearization process for the control law implementing all four corrections: spring, motor factor, BEMF and inductive, and also implementing two numerical Low Pass Filters: one between the position-indicator variable measurement and the sensor inversion, and another after the computation of the fully corrected coil voltage and before it is fed as input to the coil;
- FIG. 13 illustrates a process of applying a state variable feedback law based on a plurality of measurements of one or a plurality of state variables;
- FIG. 14 illustrates Power Spectrum Distribution simulation curves showing the effect of the transduction corrections (spring stiffness and motor factor correction) upon harmonic distortion for a single 100 Hz tone input, both with and without BEMF and nonlinear inductance in the physical model of the Labtec Spin 70 transducer;
- FIG. 15 illustrates Power Spectrum Distribution simulation curves for a single 100 Hz tone input, showing the reduction in distortion as a function of the delay in the correction loop;
- FIG. 16 illustrates simulated waveforms of the coil/diaphragm axial position versus time in the presence of a single-tone excitation, both with and without electronically restored effective spring stiffness, showing that without such restoration the cone may drift from its equilibrium position and reach its limit of excursion;
- FIG. 17 is a graph of suspension restoring force due to an electronically implemented linear spring without including the effect of the transducer motor-factor, Bl(x) ;
- FIG. 18 is a graph of the simulated phase lag between coil voltage and coil current as a function of audio frequency at low frequencies, which is almost entirely due to BEMF;
- FIG. 19 illustrates the simulated power spectrum distribution curves for the two- tone (60 Hz and 3 kHz ) intermodulation and harmonic distortion test for the 3" Audax speaker transducer, showing the forest of intermodulation peaks near the 3 kHz main peak.
- FIG. 20 is a block diagram of a control loop, including a digital controller, an amplifier, and a transducer with position sensor;
- FIG. 21 is a flow diagram of an offline calibration process for determining S as a function of position for an audio transducer using a ramped DC-voltage drive;
- FIG. 22 illustrates voltage plotted versus time for two full sweeps of the S calibration ramped DC voltage drive, including thirty-two steps of equal duration per sweep from highest to lowest or lowest to highest voltage value;
- FIG. 20 is a block diagram of a control loop, including a digital controller, an amplifier, and a transducer with position sensor;
- FIG. 21 is a flow diagram of an offline calibration process for determining S as a function of position for an audio transducer using a ramped DC-voltage drive;
- FIG. 22 illustrates voltage plotted versus time for two full sweeps of the S calibration ramped DC voltage drive, including thirty-two steps of equal duration per sweep from highest to lowest or lowest to highest voltage value;
- FIG. 23 is a general block diagram depicting an audio transducer with a controller
- FIG. 24 illustrates a plot of suspension stiffness K in Newton/mm together with a plot of Bl in Newton/amp, both of which are plotted against L e for the same Labtec Spin 70 transducer data
- FIG. 25 illustrates the S parameter plotted as a function of L e for the same Labtec Spin 70 transducer data
- FIG. 26 shows a curve that illustrates the variation of L e with position at 43k ⁇ z for a Labtec Spin70 transducer
- FIG. 27 and FIG. 28 illustrate, respectively, magnitude and phase parts of Bode plots of V rali0 for progressively larger values of L e ;
- FIG. 30 illustrate, respectively, magnitude and phase parts of Bode plots of V rati0 for progressively larger values of R.
- FIG. 31 is a block diagram for a circuit that, together with parameter estimation, measures transducer coil inductance via a supersonic probe tone and reference RL circuit
- FIG. 32 shows a curve illustrating the functional relation C (tic (x) for the mechanically moved, non-driven set of measurements of a speaker transducer
- FIG. 33 shows a curve that illustrates the variation of C pamsitic with V coil for driven measurement
- C Mc is measured in arbitrary units obtained using the method described in Detail 12
- FIG. 34 illustrates in cross-section a cell-phone speaker transducer
- FIG. 35 shows a cross-section of a portion of a speaker transducer and illustrates geometrical details of voice coil undergoing canting and its associated magnetic assembly
- FIG. 36 illustrates an audio transducer undergoing canting
- FIG. 37 is a cross-sectional view of a speaker transducer that includes an IR- LED diode and an associated PIN diode, mounted on the back side of an audio transducer of the type shown in FIG. 3, as part of an optical position detection system
- FIG. 38 is a block diagram showing in more detail an embodiment of the generalized control system shown in FIG. 8
- FIG. 39 is a block diagram of an embodiment of an audio reproduction system in accordance with one aspect of the present invention
- FIG. 40 illustrates a process flow used to linearize the transconductance component of the signal conditioning process and the transduction process of an audio transducer
- FIG. 41 illustrates the structure, in one embodiment, of the Software Control Program that is used both for obtaining data during calibration and for operating in normal mode
- FIG. 42 shows an overall flow diagram of a calibration of S and x versus f(x)
- FIG. 43 shows the details of HW and ISR operations for the S calibration in step 11504 of FIG. 42
- FIG. 44 shows a flow chart detailing the steps of mainline S calibration loop 11505
- FIG. 45 illustrates an overall flow diagram of normal mode of operation (NM, module 111104 of FIG. 41);
- FIG. 46 illustrates the operations of process 11203 of FIG.
- FIG. 47 shows a flow diagram of the ISR 11303 of FIG. 46;
- FIG. 48 shows the operation of the Wait Loop and Command Parser 11204;
- FIG. 50 shows a flow chart illustrating the details of operations performed by the DSP software in program 111208 in order to reduce the order of the approximate polynomial interpolating functions for S , x , Bl , and L e as functions of x ir for the specified rms and maximum error values, while maintaining 'Best Fit';
- FIG. 51 shows the details of the operations within step 111305 of FIG. 50;
- FIG. 52 shows a block diagram of a potential divider circuit;
- FIG. 53 shows a block diagram of the Z e (x) detection system using the probe tone 12101;
- FIG. 54 shows a block diagram of a control circuit for transducer linearization, which includes the Z e (x) detection circuit 12200;
- FIG. 55 shows a circuit diagram of the summing circuit 12202;
- FIG. 56 shows a circuit diagram of the potential divider 12203 and the high pass filter 12204;
- FIG. 57 shows a circuit diagram of the full wave bridge detector circuit 12205;
- FIG. 58 shows a circuit diagram of the low pass filter 12206;
- FIG. 59 shows the details of the circuit of the audio amplifier 12303;
- FIG. 60 shows a partial schematic and a partial block diagram of the capacitance detector and speaker arrangement, together with the DSP used for correction;
- FIG. 61 shows the input from speaker 13100 and details of the oscillator circuit 13208;
- FIG. 62 shows the detailed circuitry of the frequency to voltage converter 13210;
- FIG. 63 shows an overall block diagram of the IR-LED method for detecting a position-indicator state variable
- FIG. 64 shows a schematic diagram of IR-LED detection circuit 14400
- FIG. 65 shows a portion near 3 kHz of the FFT power spectrum distribution of the SPL (sound pressure level) wave-pattern picked up by a microphone in the acoustic near-field, with both corrected and uncorrected spectra depicted
- FIG. 66 shows a low-frequency portion of the same power spectrum distribution shown in FIG. 65, displaying multiple harmonics of the 60 Hz tone, with spectra depicted both with and without correction.
- An enabling invention in the area of control engineering 501 was the linearization method for dynamical equations 504 used in modeling physical systems to be controlled, such as actuators and transducers. This method relies on finding the control equation for the non-linear part of the dynamical equation and substituting this into the full equation.
- the application of this method to a second order differential equation 505 shows that a non-linear second order ordinary differential equation can be linearized by solving the control equation for the non-linear first order differential equation, provided the second order and first order differential terms are linear.
- This is a general method for linearizing such differential equations, and covers the application to the control of all actuators and transducer systems that can be modeled in full, or in part, by such an equation.
- the application of the linearizing method 505 to an equation with nonlinearities dependent on one state variable 506 shows that only one state variable is required for linearization.
- the application of 506 relies on positional sensing. That is to say, neither the velocity, nor the acceleration, nor the instantaneous driving force state variables are required in order to linearize the process.
- Position dependent sensing and feedback linearization can be used with many classes of non-linear motors and actuators. In the present work it was discovered that there are multiple processes in a sound reproduction system, that each process can influence the performance of other processes, that each process has non-linearities that must be considered in the design of a control loop, and that each control loop must have a sufficient number of state measurements which must be measured with sufficient discrimination against noise and with sufficient speed to control the process.
- Control of multiple processes with multiple control loops can be effected if the criteria for sufficiency is met for each control loop. It has been discovered that for the correction of non-linear transduction a necessary condition for control is a positional state measurement, in distinction to the motional measurements of prior art.
- the positional state measurement must be of sufficiently low noise and latency and of sufficiently high bandwidth to effect the control while not adding unacceptable noise to, nor engendering instability in, the sound output.
- Multiple positional measurements can be used to estimate the positional state for the purpose of transducer linearization.
- a control system approach that is based on measurement of the state of the processes in the time domain is utilized.
- the sufficiency of state measurements is based on modeling and measurement of the processes.
- Modeling of the processes in the frequency domain can also give parameters that can be reduced to the time domain.
- time domain methods can be used to measure the state of the system at each instant in time, even as the system becomes very non-linear. No assumptions need be made about the relationships of the transfer function, the input and the output.
- the signals that are used to measure state variables can come from a plurality of sensors throughout the system. Multiple state measurements are used to estimate the state of the overall system, not just the state of the output. Then, for example, amongst other properties, the instantaneous forward transduction can be estimated from a model and a measurement of the state. Thus the measurement of signals from different parts of the system is used for modeling the system response.
- the method and system comprise providing a model of at least a portion of the audio transducer system and utilizing a control engineering technique in the time domain to control an output of the audio transducer system based upon the model.
- a method to determine, in real time, the nonlinear parameters of the transducer from measurement of internal state parameters of the transducer is provided.
- the electrical properties of the voice coil can be used as a measure of positional state and a predictor of the major non-linearities of the transducer.
- "Real time" in this context means with sufficiently low latency to effect control.
- the present invention relates generally to an audio reproduction system.
- an audio transducer can be defined as: signal conditioning/ transduction/sound conditioning.
- FIG. 6 illustrates a block diagram of audio reproduction system 1100 in accordance with these processes.
- a signal conditioning process 1102 takes an audio signal 1101 (digital or analog) and performs signal conversion, amplification, filtering and frequency partitioning to provide a drive signal 1103.
- the drive signal 1103 is provided to a transduction process 1104.
- the transduction process 1104 typically utilizes a plurality of transducers, and results in diaphragm motion 1105, which drives an air load.
- Nonlinear effects resulting from sound conditioning are much smaller in normal operating conditions, and are thus neglected in the physical model described in this section, and in the control model based upon it and described in Detail 2. But these nonlinear acoustical effects, along with other higher-order effects described and then neglected in this section, can in principle also be linearized, via separate control loops according to the 'modular' approach to linearization disclosed as part of this invention. All of the effects mentioned above vary with time and circumstances. They are nonlinear and thus distort the sound wave shape, in both amplitude and phase, relative to the input audio information.
- the common method is to convert the audio signal to a voltage level, and then use this voltage to drive the impedance of the voice coil, providing current through the coil. This current then results in coil/diaphragm motion (electromechanical transduction).
- the signal conditioning may utilize a linear amplifier, in which one voltage signal is converted to another with greater driving power.
- Other options include converting the audio signal into a pulse width modulated (PWM) drive signal; thus a drive voltage is produced only during the pulse time period, thereby modulating the average current flow.
- PWM pulse width modulated
- the BEMF itself is not only dependent upon coil position, but also modulated by coil current, which introduces yet another type of nonlinearity.
- Other nonlinear response effects arise when a plurality of transducers are employed to cover a wide frequency range and the drive signal is partitioned by filters into low, medium, and high frequency ranges.
- the sound conditioning process includes the radiation of sound waves (pressure waves) from the diaphragm; reflections of the support and enclosure system (speaker enclosure) which generate multiple interfering pressure waves; and the effects of room acoustics, including noise, furniture, audience and other sound sources.
- the pressure waves present in the enclosure influence the motion of the diaphragm and the attached voice coil, thereby influencing also the signal conditioning by back-reacting upon the coil circuit.
- This back-reaction arises because the coil motion feeds into the BEMF, as well as into the coil impedance (through the latter's dependence upon coil position).
- the three processes can be described by a mathematical model, comprising a system of coupled equations specifying the rate of change (evolution) of each of a complete set of state variables, such as coil current and coil position, at any given time, in terms of the state vector at the same and all previous times.
- Such equations are termed "integro-differential equations", and are nonlinear in the case at hand.
- the model equations are usually approximated as having no "memory", in the sense that the rates of change of state variables are taken to be wholly determined by
- a nonlinear process can be very complex, and the number of terms kept in the evolution equations, as well as the decision whether or not to include memory effects, and if so which ones, can vary depending on the degree of approximation required in the control methodology.
- simplifying the approximations to the most basic mechanisms of the three processes yields several coupled "ordinary" nonlinear differential equations.
- using approximations is a compromise, and that beyond a certain point, enlarging or truncating the list of modeled effects does not alter the fundamentals of the invention.
- the most basic functionality of the signal conditioning process 1102 is transconductance, that is to say: the conversion of a voltage signal 1101 containing the audio information (audio program) into a current 1103 in the voice coil.
- the transduction process 1104 the basic functionality is the conversion of coil current to diaphragm motion (or motions) 1105; this conversion includes both electrodynamic and elasto-acoustic aspects.
- the basic functionality of the sound conditioning process 1106 is the conversion of diaphragm motion into acoustic radiation and subsequently perceived sound 1107. This can be thought of as the acoustic side of the "elastoacoustic transduction".
- i d ( t ) g,(t - ⁇ ⁇ , t - ⁇ 2 , x( ⁇ l ), x( ⁇ 2 )) i( ⁇ x )i( ⁇ 2 ) (8) is an EMF voltage term described in more detail below.
- the transduction process is governed by the mechanical equation of motion for the coil/diaphragm assembly treated as a rigid piston; including friction, acoustic loss and magnetic (Lorentz) force terms.
- Equation (14) is an oversimplification.
- a transducer voice coil is characterized by a frequency-dependent complex effective impedance, which we denote Z e ( ⁇ ,x) to indicate that it also depends upon coil position; it also implicitly depends upon other, more slowly varying parameters, such as temperature.
- the effective coil impedance Z e ( ⁇ ,x) characterizes one aspect of the relation between voltage signal V coil (t) applied to the voice-coil circuit on the one hand, and the coil current i(t) caused by this voltage, on the other.
- This voltage-current relation, or functional is nonlinear, and furthermore involves electrodynamical memory effects (distributed delays) as described above.
- this relation can be expanded in a functional series of the type known in the literature as a Volterra series.
- the multivariate coefficient-functions of this Volterra series depend on coil position and motion within the magnetic-circuit airgap. Current-nonlinear effects, i.e. deviations from linearity of the voltage-current functional, were found to be measureable.
- equation (17) In deriving equation (17) an approximation was made, namely, only linear terms in velocity x(t) were retained. This is a reasonable approximation for the physical regimes in which most speakers operate.
- equation (17) was obtained from equations (11)-(13) by dropping all EMF terms that are quadratic in the state-vector components (/(/) , x(t)) .
- the second (velocity dependent) term on the right hand side of equation (17) is the BEMF due to coil motion; the other two terms comprise the EMF due to the overall effective coil impedance.
- the subtracted impedance Z sub e ( ⁇ ,x) has both resistive and reactive components; the former is attributable to eddy-current dissipation inside the magnetic poles (and also in the coil former, in case that is made of aluminum).
- the reactive component of Z sub e ( ⁇ ,x) is known in prior art as L e (x) , with the frequency dependence often left implicit, as it was in equations (14)-(15) above.
- the subtracted effective coil impedance Z e sub ( ⁇ , x) is determined by the geometries of coil solenoid, metallic former (if any) and pole structure, as well as by the material composition within the magnetic structure (which includes the poles as well as one or more permanent magnets).
- nonlinearities arise in all of the processes involved in converting audio information into a sound wave.
- a control system such as the one described in the present invention, corrects for these distortions by applying a linearizing filter that predistorts the voltage V coil (t) applied across the coil so that it is no longer linear with the audio program signal V mdi0 (t) .
- V coil (t) applied across the coil so that it is no longer linear with the audio program signal V mdi0 (t) .
- the control paradigm used in accordance with the present invention seeks to simplify the control system by decomposing the overall control problem into reasonably independent modular parts, each of which controls a single process or sub-process. Any set of sub- processes which has already been controlled (i.e.
- FIG. 7 is a flow chart that illustrates the process of linearization in accordance with the present invention.
- a model of a portion of the audio reproduction system is provided in step 1301.
- a control engineering technique is utilized in the time domain to control an output of the audio transducer system based upon the model, via step 1302.
- the present invention controls an audio reproduction system including all three processes shown in FIG. 6. But it is not necessary that the method and system be applied to each process, but rather, that the method be available for control as the need arises.
- the model provided in step 1301 covers those processes that are appropriate to any particular implementation of the audio reproduction system.
- FIG. 8 is a block diagram of the main portion of a sound reproduction system and a control system for controlling the operation of the sound reproduction system in accordance with the present invention.
- An audio signal 1401 is input to a controller 1402, which contains algorithms based on a control model, which in turn is based on a physical model (such as the one described by equations (6)-(16) of this section) of the processes within the audio transducer system.
- the modeled processes may include the signal conditioning process 1102, the voice coil transduction process 1104, and the sound conditioning process 1106, as discussed above.
- the state variables 1403 from the sound reproduction processes are input to the controller 1402 from a measurement system 1404.
- the measurement system 1404 consists of a sensor conditioner 1405 and one or more sensors, 1406a, 1406b, and 1406c, which take measurements of variables from the sound reproduction system.
- the sensor conditioner 1405 amplifies and converts the signals from the sensors 1406a, 1406b, and 1406c to the state variables 1403, which are provided to the controller 1402.
- Sensor 1406a may, for example, measure a variable such as current from the drive amplifier 1407.
- Sensor 1406b may, for example, measure an internal circuit parameter, such as parasitic capacitance, of the transducer 1408.
- sensor 1406b could electronically measure the impedance of one of the voice coils of transducer 1408, or it could optically measure an indicator of voice coil position.
- Sensor 1406c may, for example, measure a variable from the acoustic environment, such as sound pressure by using a microphone.
- the controller 1402 modifies the audio input 1401, converts it back to an analog voltage, and thus outputs a compensated analog audio signal on line 1409 to the amplifier 1407.
- the amplifier 1407 outputs a drive signal on line 1410 to the transducer 1408.
- the audio transducer state variables that are measured and fed back to controller 1402 are generalized coordinates of the transducer dynamical system. These generalized coordinates usually vary nonlinearly with the position of the voice coil/diaphragm assembly with respect to the transducer frame, and thus, with suitable calibrations, serve to provide controller 1402 with estimates of recent values of that position. Controller 1402 then uses these real-time position estimates to suitably modify the input audio voltage signal before applying it across the voice coil. Multiple position- indicating signals can be fed to the controller, as depicted in FIG. 8; they are derived from one or more position-indicating generalized coordinates.
- the advantage of measuring and feeding back values for multiple generalized coordinates is that these coordinates may be chosen in such a way that the configuration space of their joint values is approximately a one dimensional differentiable manifold, where the coil/diaphragm position is a continuous and differentiable function on this manifold.
- each of the selected generalized coordinates is also a continuous and differentiable function of coil/diaphragm position
- the mapping between a tuple of simultaneously measured generalized coordinates and the corresponding position is both invertible and differentiable, allowing the use of the tuple to compute the audio signal modification within the controller DSP.
- One embodiment of this computation based on a single generalized coordinate that is derived from infrared optical measurements, is described in detail in Detail 10. It will be readily apparent to those skilled in the art that additional and different sensors may be utilized, and different signal conditioners may be used to recover state variables and internal parameters from the sensor signals and provide control signals to the system.
- Additional sensors may include, for example: accelerometers, additional transducer coils, or new coil-circuit elements.
- Such sensors can provide analog measurements of various voltages appearing in the transconductance equation (14), or of other voltages that allow the estimation of various terms and state variables in either equation (14) or the mechanical (transduction process) equation (15).
- State variables and parameters must be identified for each of the sound reproduction processes, and a sufficient set of them must be measured to effect control. It has been discovered that measurements not usually regarded as state variables can be used effectively in controlling the audio reproduction processes.
- Some measurable variables can be measured by reference to other variables through known functional dependencies; for instance, temperature can be inferred from coil resistance and a lookup table.
- Internal parameters and other variables not listed in the above list include, for example: V(t) voice coil voltage, i(t) voice coil current, R e voice coil resistance, L e voice coil inductance, Z e complex voice coil impedance, C paras i t i c voice coil/magnet parasitic capacitance, BEMF back-EMF, ⁇ complex phase angle of voice coil impedance, T e voice coil temperature.
- There are other internal parameters such as Bl and K , respectively the motor factor and suspension stiffness. These parameters may be difficult to measure directly, although they can be extracted from measurements of other variables via parameter estimation methods.
- the voice-coil voltage V(t) and voice coil current i(t) are considered internal variables, rather than stimuli, because the full audio transduction process according to the present invention includes creating V(t) and i(t) as internal variables.
- DETAILED DESCRIPTION 2 CONTROL MODEL
- the present invention is described in the context of controlling part of or all of an audio reproduction system using a control model.
- the control model is based upon the physical models for one or more of the three processes in the audio reproduction system; these processes, and physical models for their main components, were described above (Detail 1 ).
- the control model is based on the physical models expressed by the electromechanical evolution equations (14) and (15), but with terms non-linear in velocity and/or current neglected.
- the electrical circuit equation (18) describes the transconductance component of the signal conditioning process; whereas the mechanical equation of motion (19) describes the transduction process.
- a modular control model was developed in the context of the present invention, including separate corrections of nonlinearities in the transduction and signal- conditioning processes based on the measurement of a minimum of one position- indicator state variable during operation.
- an implementation of this control model removes a significant and adjustable portion of the audio distortions caused by the nonlinearities in equations (18) and (19). Furthermore, the control model removes nonlinearities in a modular way. Specifically, as described in the remainder of this section, this control model linearizes either the BEMF voltage term in the transconductance equation (18), or it linearizes the effective voice-coil inductance term in equation (18), or it linearizes the suspension stiffness and/or motor drive factor in the mechanical transduction equation (19); or it linearizes any combination of these.
- the particular combination of modular control laws implemented in the controller is determined by user preferences. And all modular control laws are based upon a single state measurement of position, or of a position-indicating variable.
- the linearizations are performed in a controller, such as that described in connection with FIG. 8.
- the control model treats the motor factor R/(x) , the effective coil inductance
- L e (x) and the suspension stiffness K(x) as functions of x(t) , the current axial position of the coil/diaphragm assembly. These three functions cause most of the nonlinearities, and thus distortions, of audio transducers, as explained above.
- the motor factor R/(x) determines the motive force term in equation (19) as well as the BEMF term in equation (18);
- L e (x) determines the inductive EMF term in equation (18); while
- K(x) determines the elasto-acoustic restoring force in equation (19).
- these three functions are derived from calibration measurements on the system, which yield the functional dependence of Bl , L e and K upon x ; these functions can, for instance, be obtained from commercially available transducer test equipment such as a Klippel GMBH laser metrology system.
- the functional dependences R (x) and L e (x) are entirely obtained from such a laser metrology system, while K(x) is obtained by combining knowledge of R (x) and L e (x) with ramped DC-drive calibration runs, as fully described in Details 5 and 10 below.
- the three functions RJ(x) , L e (x) and K(x) must be combined with approximants to a function mapping the measured position-indicator state variable onto the actual position , as described in Details 4,5, and 10 below, in order to provide the controller DSP with an estimate for the values of Bl , L e and K(x) at the present moment t .
- the controller estimates the BEMF term by multiplying the estimated present value of Bl(x(t)) by an estimate for the present velocity x(t) ; the latter may be obtained either from a numerical differentiation of the recent history of discrete position measurements, or from an independent velocity measurement. In one embodiment of the present invention, velocity is estimated via numerical differentiation of estimated position, as described in Detail 10 below.
- the inductive control law is not neglected, and describe another type of modular control law in the context of the present invention, namely a control law that corrects for the inductive EMF term in equation (18).
- the inductive control law partially linearizes the transductance sub-process. Specifically, the inductive control law addresses the nonlinearity, and thus distortion, caused by the position dependence of the effective coil inductance L e (x) .
- the BEMF term is temporarily ignored in the transconductance equation (18); later in this section, all four of the modular control laws described in the context of this invention (BEMF, inductive, spring and motor factor) will be combined.
- equation (23) can be viewed as a linear first-order ordinary differential equation for the unknown function i(t) .
- the nonlinear effects in the transconductance equation (18) can be partially eliminated in a modular manner by the control laws given by equations (20) and (22), leaving approximately linear effects for the back-EMF and inductive EMF, respectively.
- the BEMF and inductive EMF corrections have little overlap in frequency; that is to say, the BEMF has significantly lower frequency content than the inductive EMF. Therefore, the order of application of the two separate modular control laws thus far described in this section, equation (20) for BEMF and equation (22) for the inductive term, should not greatly matter in terms of amount of distortion reduction, in case the user elects to implement both of these control laws.
- the correction of the nonlinear electromechanical effects in the mechanical (transduction) equation of motion (19) is based upon a derivation similar to, but different from, the standard control theory derivation of a control equation presented in the Background section above as prior art.
- One practical problem with the mechanical equation (19) as a starting point for a control model is that the inertia term involves the coil/diaphragm acceleration x . This term increases rapidly with frequency, eventually becoming too large to be considered in a compensation system. However, because the acoustical radiation efficiency of the cone also increases with frequency, the inertia noncompensation is balanced by the radiating efficiency, within limits.
- the controller linearizing the transduction process should cause the transducer output x(t) to be proportional to the audio input. Equating x(t) with V mdi0 (t) in equation (26) and solving for u(t) , and assuming that the function ⁇ (x) defined in (27) is nonsingular, we obtain: where w(t) is the generator or reference (in our case the audio program input V mdio (t) to the uncorrected transducer), and R e u(t) is the actual voltage input to the voice coil in the controlled (corrected) transducer if the signal conditioning process is ignored. Substituting and rearranging terms in equations (27), (28) and (30), provides:
- Equation (33) is a linear differential equation with constant coefficients. Note o that from the above a general method of linearizing this form of nonlinear dynamical equation is presented, and any further linear terms can be added to the equation without changing the validity of the linearization approach.
- Equation (34) provides a correction for the open loop non-linear transfer function of the speaker transducer, provided that the dependencies of S(x) and R(x) on x are 0 known and that real-time measurements or estimates of x are made available to the controller during transducer operation. The validity of equation (34) as a control law can be simulated when applied to a full physical model of an actual transducer.
- control models of equations (20) and (22) above suitably combined, eliminate or reduce only those nonlinearities arising from the transconductance component of the signal conditioning process, but do not correct either of the other two processes (transduction or sound conditioning).
- all of the above control laws can, and have been, applied together, or in various partial combinations, in the context of the present invention.
- the transduction control law of equation (34) can be subdivided into "spring correction" and "motor factor” modular units; e.g. if only the first term on the right hand side of equation (34) is used, this represents a control law which only linearizes the elastic restoring force.
- the number of modular control laws described by the above equations can actually be counted as four: BEMF, inductive, spring, and motor factor. If a choice is made to simultaneously implement all of these modular corrections: the BEMF correction (equation (20)), the inductive correction (equation (22)), and the transduction corrections (equation (34)), this can for example be done as follows.
- the last term of equation (20) is added to the voltage given by the right-hand side of equation (34); then the new overall voltage, u x (t) , still in the digital domain, is numerically differentiated (as described in Detail 10 below), and this numerical derivative is combined with u x (t) itself in accordance with equation (22).
- u(t) S(x)- ⁇ - + wB(x) (40) R/ (x)
- FIG. 9, FIG. 10, FIG. 11 and FIG. 12 are process block diagrams depicting the workings of various possible combinations of control laws as applied to the overall three- state system, or to parts thereof, in the context of the present invention. What follows is a detailed description of these diagrams.
- FIG. 9 shows the feedback linearization process 20400 with the control law of equation (34), which only linearizes the transduction component of the signal conditioning process, without an electronically restored linear restoring force.
- the audio signal, 20401 is input to a Linear Compensation Process module 20402 (henceforth abbreviated as LCP).
- LCP 20402 multiplies w by the compensation function B(z) , where z 20411 is the estimated present value of the position variable.
- the present value of position variable z 20411 is obtained from the transduction module 20408 of the three-state overall transducer system, via a two step process: first the position indicator state variable f(x) 20413 is measured by the positional sensor module
- a sensor inversion module 20414 which estimates actual position x via an interpolation method as described in Details 5 and 10.
- Actual position x 20409 and actual velocity x 20410 are fed from the output of transduction module 20408 back into the input of the transconductance module 20406, via the physical system itself (not as measured data).
- the estimated x value, z 20411 is fed into to the LCP 20402 and also to an S -lookup module 20415.
- module 20415 S(z) * S(x) 20416, as well as the LCP output B(z)w 20403, are both fed as inputs to a summing junction 20404, the output 20405 of which is the corrected audio signal (V coil of equation (34)).
- This corrected audio signal 20405 is provided as input to the transconductance module 20406 of the three-state transducer system.
- the current output I coil 20407 of the transconductance module 20406 is provided as input to the transduction module 20408.
- FIG. 10 shows the feedback linearization process 20500 for the control law given by equation (40); again only transduction corrections are made, but now a linear spring constant (suspension stiffness) is electronically added, as explained above and in Detail
- the LCP 20502 multiplies w by the compensation function B(z) , where z 20514 is the estimated current value of the position variable.
- Value z 20514 is obtained from the transduction module 20508 of the three-state overall transducer system, via a two step process as in FIG. 9: the positional sensor module 20511 outputs the measured position indicator state variable f(x) 20512, and measured state variable f(x) 20512 is fed as input to a sensor inversion module 20513, which estimates actual position x via the interpolation method.
- Actual position x 20510 and velocity x 20509 are fed back from the output of the transduction module 20508 to the input of the transconductance module 20506 via the physical system itself.
- the estimated x value, z 20514 is this time fed into three modules: to the LCP 20502, to an S -lookup module 20516, and to a new 'Electronically Restored Linear Spring'(henceforth ERLS) module 20517.
- FIG. 11 shows the feedback linearization process 20600 for the control law given by equation (39) alone, without the inductive correction (36); i.e. for a control law correcting for spring, motor factor and BEMF nonlinearities, including an electronically restored linear spring and electronically restored contribution to the linear drag force term, as explained above.
- the LCP 20602 multiplies w by the compensation function B(z) , where z 20622 is the estimated present value of the position variable.
- Value z 20622 is obtained from the transduction module 20610 of the three-state overall transducer system, via a two step process as in the previous figures: the positional sensor module 20613 outputs the measured position indicator state variable f(x) 20614, which is then fed as input to a sensor inversion module 20615. Sensor inversion module 20615 estimates actual position x via the interpolation method. And as in previous figures, the actual position x 20612 and velocity x 20611 are fed back by the actual physical system from the output of the transduction module 20610 to the input of the transconductance module 20608.
- the estimated x value, z 20622 is now fed into four modules: to the LCP 20602; to the S -lookup module 20618; to an ERLS module 20620; and finally, to a BEMF-computation module 20616, which applies a numerical differentiation operation D to z 20622.
- the output 20619 of the module 20618, as well as output 20621 of module 20620 and output 20603 of the LCP 20602, are summed in the summing junction 20604.
- FIG. 12 shows the feedback linearization process 20900 for the control law given by equations (36) and (39), i.e.
- the audio signal, V mdi0 w 20901, is input to an LCP module 20902.
- the LCP module
- LPF2 would typically roll off at 1 - 2 kHz .
- the output 20925 of LPF2 20924 is fed to the sensor inversion module 20914.
- Sensor inversion module 20914 again estimates actual position x via the interpolation method, in the digital domain; while the actual position x 20911 and velocity i 20912 are fed via the physical transducer plant, back from the transduction module 20910 to the transconductance module 20908.
- the estimated x value, now called z f 20921 is fed into the following three modules: to the LCP 20902, to the ERLS module 20920, and to the BEMF-computation module 20915.
- the S - lookup module 20917 receives its input this time from the filtered, but not inverted, positional indicator variable measurement result 20925.
- the output 20905 of summing junction 20904 is passed to an inductive-correction module 20927, which again applies a numerical differentiation operation D , this time to the numerical output voltage 20926 of the summing junction 20904.
- the output 20906 of the inductive- correction module 20927 is provided along with numerical output voltage 20926 multiplied by R e to a second summing junction 20928, whose output 20907 is fed to the low pass filter LPF1 20922.
- the low pass filter LPF1 20922 implements a (partial) correction for the voice coil inductance at equilibrium.
- the output 20923 of LPF1 20922 is finally fed as the corrected analog voltage V coil to the transconductance module 20908 of the three-state transducer system.
- the physical transducer plant provides the analog output current I coil 20909, output by the transconductance module 20908, as input to the transduction module 20910.
- the present invention requires at least one state variable to be measured in operation for any given run. In the control diagrams depicted in FIG. 9, FIG. 10, FIG. 11 and FIG.
- the process of applying a state variable feedback law based on a plurality of measurements of one or many state variables is depicted in FIG. 13.
- the process 21000 begins with one or several measurements of a state variable or variables from a plurality of sensors, 21001 through 21002.
- a transducer's coil/diaphragm displacement may be measured both via the parasitic capacitance method (Details 7 and 12 below) and the IR method (Details 8 and 13 below).
- the respective state variable measurement signals, 21003 through 21004 are passed from the sensors to the state estimation module 21005, which synthesizes the desired partial or full state variable estimate, 21006, which in general is a vector state variable.
- This state variable estimate 21006 is in turn used in the application of the control law 21007 in place of the actual state variable.
- none of the sensors, 21001 through 21002 can measure its intended state variable exactly.
- the measurement is always corrupted to some extent by factors including nonlinearities in the measurement, measurement noise, quantization noise, systematic errors, etc.
- the task of the state estimation module 21005 is to mitigate these corrupting effects. This task may include all or some of the following ingredients: inverting the nonlinearities of the sensors to provide a more linear response to the measurements 21001 through 21002; adaptation to minimize the sensitivity of the state variable estimate 21006 to parametric uncertainties in the measurement, such as uncertainty in gain; filtering the measurement signals 21003 through 21004 to minimize the effects of noise; or fusing multiple measurements of a state variable into one state variable estimate 21006.
- many engineering objectives are taken into consideration in the design of the state estimation module 21005.
- the tradeoffs include such desirable properties as simplicity of design, overall reduction in the effects of noise in the system, minimization of the order of the state variable estimator, and cost of implementation.
- one possible method by which to invert the nonlinearities in any of the measurements 21001 to 21002 is via a lookup table based upon offline calibration runs; another possible method, also based upon offline calibration, is via polynomial expansion. The latter is the method used in one embodiment of the present invention, as described in Detail 10 below.
- Noise reduction may be accomplished by filtering, for example by using finite impulse response (FIR) or infinite impulse response (IIR) digital filters, or else analog filters.
- FIR finite impulse response
- IIR infinite impulse response
- the structure of an IIR noise reduction and data fusion filter, and its coefficient values may be determined by trial and error or by analysis.
- a positional estimation filter could be designed via Kalman filtering techniques, in which a stochastic model of the input signal and state measurement noise is combined with a model of the transconductance and transduction dynamics (such as equations (18)-(19) above) to resolve the order and coefficient values of the estimation filter.
- Kalman filtering techniques in which a stochastic model of the input signal and state measurement noise is combined with a model of the transconductance and transduction dynamics (such as equations (18)-(19) above) to resolve the order and coefficient values of the estimation filter.
- DETAILED DESCRIPTION 3 JUSTIFICATION OF APPROXIMATIONS
- these mathematical models are used in the context of the present invention: in actual physical implementation, and in simulation.
- the chosen collection of one or more of the four basic control laws is implemented within DSP hardware and software, which control the transducer in order to linearize sound.
- nonlinear effective spring "constant” nonlinear motor factor
- nonlinear effective voice coil inductance nonlinear motor BEMF
- computer simulations based upon the transducer-plus-controller model were used in the present work to study the effect of all of these nonlinearities, thereby elucidating the merits of implementing partial correction for a subset of the nonlinearities. For instance, it was found via simulation that transconductance nonlinearities (BEMF and inductive) are responsible for significant audio distortions at various important frequency ranges, which led to the inclusion of corrections for these effects in the control law (equations (20) and (22) above).
- FIG. 14 shows curves 4100 of simulated Power Spectral Density (PSD) which illustrate the effect of the transduction corrections alone (spring stiffness and motor factor correction, equation (34)) both with and without BEMF and nonlinear inductance in the system.
- PSD Power Spectral Density
- the vertical axis is a measure of PSD in relative dB units.
- the curves of FIG. 14 were generated by simulating the performance of a particular transducer (that of the Labtec Spin 70 speaker) using a single 100 Hz tone; each curve clearly shows that the highest power is in the fundamental 100 Hz tone, but that significant power is also present in the various harmonics of this tone. Overall, the curves of FIG. 14 shows that even at frequencies where BEMF is significant, introduction of corrections for spring and force constant greatly improve the system performance.
- Curve 4103 depicts the simulated PSD with no BEMF voltage term modelled, with linear (i.e. position independent) inductive EMF voltage term modelled, and with no correction incorporated in the modelling; the harmonics, and power present at non-harmonic frequencies, are an artifact of the finite time windowing used to perform the FFT (Fast Fourier Transform) in the simulation.
- Curve 4101 shows the PSD when the position-dependent (nonlinear) BEMF and position-dependent inductive EMF voltage terms both modelled, but still with no correction; the harmonics, as well as the general diffuse high-frequency content of the power spectrum, are seen to be enhanced by nonlinearity-caused distortions.
- Curve 4102 again depicting the PSD with nonlinear
- curve 4104 depicts the PSD with no BEMF and with linear inductive EMF, as in curve 4103, but with the difference that the transduction correction is applied. It is inevitable that there will be some delay between measuring and reading the sensor output, and sending out the command to compensate for the position-dependent nonlinear spring stiffness and motor factor (and for any other nonlinearities for which terms are included in the controller). Using model-based simulation, it was possible to determine that the existence of this delay, while somewhat degrading the performance of the control algorithm, did not cause a significant problem, nor did it render the algorithm ineffective.
- the vertical axis is a measure of relative PSD magnitudes in dB.
- the curves of FIG. 15 depict the simulated PSD of the transducer-cone velocity, again for a 100 Hz audio input tone. In obtaining these simulation results, it was important to keep the amount of the nonlinearities the same for all the cases that were considered. This was achieved by suitably scaling the driving force as the time delay was varied. It is clear, from the curves of FIG. 15, that longer delays in the correction loop will increase distortion. However, for a 100 Hz tone, even at 200 ⁇ sec delay, the distortion is seen to be less than that of the uncorrected system.
- Curve 4201 depicts the PSD with no correction
- curve 4202 depicts the PSD with transduction correction but for the ideal case of no delay
- curves 4203 and 4204 show the PSD curves with correction modelled and with simulated delays in the amounts of 100 ⁇ sec and 200 /sec , respectively. While a complete nonlinear spring cancellation will reduce the distortion in a speaker's acoustic output, it will also remove the restoring force that was provided by the mechanical spring in the uncorrected speaker transducer, as discussed in Detail 2 above. In order to keep the speaker cone centered near its equilibrium position and place the mechanical resonance of the speaker at the desirable frequency, linear stiffness can be added electronically, as seen in Detail 2 above.
- FIG. 16 displays a plot
- FIG. 17 shows the spring force due to an electronically implemented linear spring without including the effect of the transducer motor-factor, Bl(x) .
- FIG. 18 depicts the simulated phase lag between coil voltage and coil current at low audio frequencies, which is almost entirely due to BEMF. At high frequencies this phase lag would be mainly due to the inductive term in the electrical circuit equation (18).
- FIG. 19 is a simulated version of spectral plot results 4600 of the two-tone intermodulation and harmonic distortion test for which actual, physical implementation results are reported in Detail 14 below. The two input tones are at 60 Hz and 3/ ⁇ Hz , and the portion of the simulated power spectrum distributions (PSDs) shown in the curves of FIG. 19 are in the vicinity of 3 kHz .
- PSDs power spectrum distributions
- the curves (4601 through 4603) clearly show the forest of intermodulation peaks, spaced uniformly 60 Hz apart and with decreasing power level away from the 3kHz main peak.
- the simulation shows the intermodulation peaks to be significantly suppressed when all four linearizing-filter corrections are applied (i.e. with the combined correction law given by equations (36)-(39)).
- Curve 4601 shows the simulated uncorrected PSD
- curve 4602 shows the dramatic intermodulation reduction when the corrections ' are applied, with 10// sec simulated delay.
- curve 4603 shows the simulated PSD with corrections and with the longer simulated delay of 50 //sec . It is seen that while the larger delay increases distortions, even the corrected spectrum with the higher simulated delay value is still less distorted than the uncorrected spectrum with no delay at all. It will be clear to those skilled in the art that simulation of any particular implementation of the linearization and control methods described in this disclosure provides valuable information for practically implementing such systems for any particular application; and, furthermore, that the simulations developed here can be greatly expanded to cover many such systems and applications.
- DETAILED DESCRIPTION 4 STATE MEASUREMENT THEORY
- the present invention is described in the context of controlling an audio reproduction system, in part, by a model requiring real time measurement of at least one position-dependent state variable of the speaker transducer.
- one such state variable is the axial position x of the coil/diaphragm assembly. Real-time values of the state variable x are needed during transducer operation in order to effect the linearization of the transconductance and transduction processes, as set out in Detail 2.
- a position-indicator state variable i.e. a variable which varies monotonically (but, in general, nonlinearly) with x within the range of possible diaphragm excursions.
- This disclosure discusses in detail three main choices of f(x) measurement techniques: an optical method using IR detection; a method using the effective impedance, or inductance, of the voice coil; and a method that uses the parasitic capacitance between the voice coil and the magnet assembly of the transducer.
- the above-mentioned three methods are referred to as the IR method, the Z e (or L e ) method and the C method, respectively.
- Other choices of position-indicator state variables could be made, depending on the application.
- the IR method is fully described in Details 8 and 13.
- the Z e method is fully described in Details 6 and 11.
- the C method is fully described in Details 7 and 12.
- the position information derived by Z e and C methods is generated using internal electronic parameters of the transducer.
- the IR method is based on an external measurement of position.
- the respective variables must be monotonic, but not necessarily linear, with position.
- position indicators there are other possible position indicators according to the present invention, which are measurable from internal electronic circuit parameters of the transducer that are not constant during transducer operation, but instead vary monotonically with x .
- K(x) , RJ(x) , and L e (x) are commonly presented as the parameters most responsible for the nonlinearities in the operation of such a transducer.
- variables S and B can be redefined as functions of x ir ,L e , Z e or C parasitic , depending on the positional-detection method selected.
- i ⁇ t) S ⁇ r (x ir ) + wB lr (x ⁇ r ) (41 )
- i ⁇ t) S L ⁇ L e )+ wB L ⁇ L e ) (41a)
- i(t) S z (Z e )+ wB z ⁇ Z e ) (42)
- i(t) S c (C parasitlc ) + wB c (C paras ⁇ (42a)
- the transduction control law, equation (34), has been used to illustrate the use of nonlinear position indicators for linearization corrections.
- the same indicators can be used for some of the other corrections that can be added in a modular fashion to any particular implementation.
- These combinations of the modular control laws, described in the context of the present invention, are given by the control equations (20), (22), and (36)-(39) in Detail 2 above.
- the motor factor R (x) can be stored in the controller as a function of the nonlinear state variable f(x) , while the instantaneous velocity x can be obtained not by measuring a motional state variable, but rather via numerical differentiation of the position, which in turn is obtained from f(x) via the stored inverse functional relation f ⁇ l .
- All controller-stored functions, whether having the form of polynomials, look-up tables or splines, or some combination of the these, will be computed, based upon calibration or characterization of the transducer, 'offline'; i.e. before actual transducer operation.
- L e (x) can be characterized as a function of the position-indicator variable f(x) , while the time derivative of the voltage can again be computed numerically.
- Information from other external measurement apparata not utilized in the context of this invention, such as accelerometers, microphones, voltages from additional coils and/or additional transducers, can also be used to provide additional state variables, and thus can be used to add precision to, or reduce the noise, for positional or motional estimates.
- FIG. 20 shows a block diagram of a control loop 6100.
- the control loop 6100 includes a digital controller 6101 , an amplifier 6102, and a transducer 6103 with position sensor 6104 (illustrated graphically) that outputs a measurement of a signal which is indicative of a state variable that is a monotonic, and generally nonlinear, function of position, f(x) 6105.
- This nonlinear state variable could be an internal circuit parameter or a signal from an external position-sensing device.
- the nonlinear state variable serves as a measure of position in the control system according to the present invention.
- Values for S can be measured directly from the control loop 6100.
- the spring force term S can be output independently simply by outputting a DC value - because for a DC signal, the only force in the correction equation is the static (spring- force to motor-factor ratio) term S(x) , and the numerical value of S can thus be measured.
- FIG 21 is a flow diagram of a process for determining S as a function of position of a transducer.
- FIG 22 shows the voltage waveform 6206, the current from which is utilized to move the cone of transducer 6103 and thus to determine and plot S as a function of .
- control parameter S used in the control loop 6100 is the transducer-coil current in voltage units - which is taken to be V . This procedure is approximately correct (in the case of a voltage controlled amplifier assumed here) to the extent that the non-Ohmic EMF terms in the coil circuit, including the effective coil inductance and BEMF voltage terms, are neglected. This is a justifiable approximation for sufficiently slow ramping, i.e. long ramp-times and settling times. The ramp is made slow relative to audio signal timescales, because it is undesirable to put out audio information in the ramp.
- the current into the coil is proportional to voltage by Ohm's law, to a good approximation.
- waveform 6206 shown in FIG. 22 included thirty-two steps of equal duration per each sweep from highest to lowest or lowest to highest voltage value.
- the output voltage was zero.
- the voltage increment or decrement was 1 / 16 th of the zero- to-peak amplitude of the waveform 6206, which was 0.25 volt. This value was before amplification.
- the amplitude of the ramp-sweep voltage signal fed to the voice coil of transducer 6103 was about 20 times higher. This amplitude is determined, for each speaker transducer, by the need to cover the full excursion of the coil/diaphragm motion that is encountered in normal operation.
- each thirty-two step sweep was completed over a one-second time interval, and two such full sweeps are shown in FIG.
- FIG. 22 only shows half the number of DC voltage steps per sweep as were actually used for the case of the 3" Audax speaker transducer.
- a table of the V(n) outputs, and the corresponding measured values of the nonlinear position-indicator state variable f( ⁇ personally) is created.
- This table is then polynomial-fitted to yield an approximate polynomial interpolating formula for the function S o f ⁇ l , or (more generally) a new look-up table for interpolation of this function; in general both approaches could be used, for example via a polynomial spline (piecewise polynomial) and interpolation.
- the values of V(n) in the table can be either actual voltage values, or values in the numerical format used by controller 6101.
- the output values of V(n) could be fixed format digital words that are output to a digital-to-analog converter (DAC).
- the off- resonance choice of tone frequency provides a relatively simple relation between the measured SPL and the motor factor Bl , which in turn is inversely related to R .
- the deduced values of R can then be tabulated against corresponding measurements of f(x) , for a stairway-ramped voltage signal 6206, in a manner similar to that used in the
- the low-amplitude tone is applied after that DC level has been held a sufficient time to allow electromechanical relaxation of the transducer to a steady state current and mechanical equilibrium.
- the frequency of the tone is fixed for each stairway-ramped voltage sweep, but can be varied from sweep to sweep.
- the foregoing approach is complicated by two factors. Firstly, the speaker's acoustic transfer function (diaphragm motion to SPL) is not a priori known for realistic speaker enclosures; and secondly, the suspension stiffness still affects the conversion of SPL values to R values, through the xicide -dependent elastic resonance frequency, for tone frequencies low enough so that coil-inductance effects do not spoil the simple Ohmic conversion of voltage to coil current.
- this initial stage need only be performed once per given speaker, since drifts in the motor-factor function Bl (x) are almost entirely multiplicative, stemming from temperature dependence of the airgap magnetic field, and thus hardly affect the ratio B(x) .
- a stairway-ramped voltage sweep of the type described above is performed, in which the position-indicator nonlinear state variable f(x) and the actual position x are simultaneously measured. The latter is measured via a Klippel-type laser, which returns a voltage known to vary linearly with actual position to a high accuracy.
- FIG. 23 is a general block diagram of a system 6300 depicting an audio transducer 6304 with the digital controller 6301.
- Digital controller 6301 received two inputs: the audio voltage signal w 6302 (also referred to as V mdh ; see Detail 2), and the most recent measurement of the position-indicator nonlinear state variable f(x) 6303. This nonlinear state variable is measured in the transducer 6304.
- Digital controller 6301 combines the audio input with the measured value of f(x) to compute the corrected
- control law may be that given by equation
- the voltage V coil is output in analog form 6305 by digital controller 6301, and provided to the amplifier 6306.
- the output voltage from amplifier 6306 is provided to transducer 6304.
- Equation (46) the three quantities S , B and S' are all expressed as interpolated polynomials in the measured position- indicator nonlinear state variable f(x) , as described above.
- a suitable choice of the residual linear spring coefficient k in equation (46) is also important in order to tune the resonant properties of the transducer appropriately for the given program material: a low effective spring stiffness will yield a low resonant frequency, and vice versa.
- parameterized linearization-filter functions characterizing the given transducer which are measured and estimated using in-operation measurements of at least one nonlinear position-indicator state variable, augmented by preliminary (characterization) calibration runs in which this nonlinear state variable is measured simultaneously with a more linear position-indicating variable (such as the Klippel-GMBH laser metrology system).
- the nonlinear position-indicator variable measured in operation can be a voltage output from an optical device, as is the case in one embodiment of the present invention and as is described in Details 8 and 13 below; or it could be an output from the internal electronic parameter measurements, as described in Details 6,7,11 and 12.
- These measurements could be augmented by an external measurement of sound pressure level during characterization runs, as described above.
- the S and R parameters which are needed by the controller to implement the transduction-process portion of the linearizing control law, can be matched to the program material by adjusting the parameter q governing the electronic spring force compensation, as described in equations (39), (40) and (46).
- curve 5101 is a plot of K in Newton/mm and curve 5102 is a plot of Bl in Newton/amp, both of which are plotted against L e for the same data from FIG. 4.
- This new mapping provides a basis of a correction scheme. Because the inductance of the voice coil is a function of its position, by measuring the inductance the position of the voice coil is determined. Thus L e provides an inductive position detector. From the definition of S and R in Detail 2, it can be seen that S is a function of x (determined by the functions K(x) and R/(x) ) and can thus be expressed and plotted as a function of
- FIG. 25 displays S plotted as a function of L e for the same Labtec Spin 70 transducer data as in FIG. 4. Similarly R can be plotted versus L e .
- the use of the voice coil inductance, L e as a position estimator can be generalized as a method by considering that we are in fact using the effective complex voice-coil impedance Z e ( ⁇ ,x) , defined in Detail 1 above, to provide the estimate f(x) .
- the effective complex voice-coil impedance Z e ( ⁇ ,x) is measured electronically at some suitably chosen supersonic probe-tone frequency.
- the reactive component of Z e ( ⁇ ,x) that is L e
- L e is also a state variable that depends monotonically upon .
- the variation of L e with position at 43kHz is shown in FIG. 26 for a Labtec Spin70 transducer.
- the impedance Z e ( ⁇ ,x) depends not only on coil position x and probe tone frequency ⁇ , but also on the temperature distribution in various components of the transducer; the most significant such dependence is upon the average instantaneous voice-coil temperature, T coil .
- a Z e method which involves electronically measuring Z e ( ⁇ ,x) , for a range of values of coil/diaphragm position , using a suitably chosen supersonic probe-tone frequency ⁇ , and encoding the resulting function Z e (x) via a polynomial fit to the measured data.
- the polynomial fit can be used during speaker operation to dynamically calculate the current value of x(t) from the electronically measured values of Z e ; the calculated x value is input into a correction (any of the linearizing-filter control laws described in Detail 2 above).
- the fitted function is used to generate and store a Look-Up Table (LUT).
- LUT Look-Up Table
- This implementation utilizes a potential divider circuit to measure the overall (complex) effective coil impedance, Z e ( ⁇ ,x) , at the particular probe tone frequency of 43kHz , with no attempt at either theoretical modeling of the trivariate complex function Z e ( ⁇ ,x, T coil ) , or at separating the real (resistive) component of Z e from its imaginary (reactive or inductive) component.
- L e (x) curve shows a typical prior-art L e (x) curve, coil inductance versus coil position, obtained by polynomial fitting of data at audio frequencies; the impedance measurements upon which FIG. 4 was based ignore the resistive component of Z e .
- the inductance changes monotonically with position, and measurement of this inductance thus yields a suitable substitute for the coil position itself in the control model of the present invention.
- FIG. 31 A method for measuring the coil inductance is illustrated by the block diagram in FIG. 31.
- a supersonic probe tone (“carrier signal") is applied via input line 7401 to the voice coil of transducer 7402.
- a reference R L circuit 7403 is placed in series with the voice coil.
- the supersonic signal is then injected into the voice coil of the transducer 7402 in addition to the audio signal, and the voltage across the voice coil of the transducer 7402 and the reference R L circuit 7403 is measured.
- Reference R L circuit 7403 may be implemented using a resistor and a coil in series. Alternatively, a coil or a resistor may be used to implement circuit 7403.
- the measured voltage signals are sent via summing junction 7404 and summing junction 7405 through filter 7406 and filter 7407, respectively, and the ratio of the output of the filters is then determined in either the analog or digital domain.
- the filter 7406 and filter 7407 are band pass filters implemented about the frequency of the carrier signal. Envelope detection via envelope detector 7408 and envelope detector 7409 is used to extract the signal due to changes in L e . The ratio of the voltages coming out of the envelope detector 7408 and detector
- R e ' is the resistive component of coil impedance at the probe tone frequency, including both the Ohmic coil resistance R e and the lossy effective coil impedance component due to eddy currents.
- R ref and L ref are the respective series resistance and inductance of the reference R L circuit 7403; and s is the Laplace variable. Because the ratio of the two voltages is taken, the signals that are close in frequency to that of the carrier, and thus cannot be rejected by the band-pass filter 7406 and filter 7407, will not introduce significant error in L e determination. As long as L ref and R re/ are chosen
- FIG. 27 shows the Bode plot of the transfer function V ratt0 given in equation (47), while FIG. 28 shows the corresponding phase Bode plots.
- the ordinate in FIG. 27 is the magnitude of V raii0 , in dB units, while the ordinate in FIG. 28 is the phase of V rati0 , in degrees; in both plots, the abscissa represents angular frequency in units of radians per second.
- the family of Bode plots is for progressively larger values of L e , with the highest L e value resulting in the curve 7201 and curve 7204, while the lowest value results in the curves 7202 and
- V rati0 to changes in L e is clearly a function of the probe tone frequency. The higher this frequency, the more sensitive V rati0 will be to L e variations. To reduce the effect of the common mode in-band noise, which is present in the voltage across the voice coil (i.e. (L e s + R e ') - i ) and in the voltage across the reference rcuit , p i
- FIG. 29 shows a series of Bode plots for the magnitude 7300 of V rati0
- FIG. 30 shows the corresponding plots for the phase of V rati0 V rati0 7303.
- Each plot-pair is for one of a decreasing sequence of R e values, and thus corresponds to a sequence of decreasing voice-coil temperatures, for example, the magnitude plot for the highest R e value 7301; the magnitude plot for lowest R e 7302; the phase plot for highest R 7304; and finally, the phase plot for lowest R e 7305. Because FIG. 27, FIG. 28, FIG. 29 and FIG.
- V rati0 can be utilized to accurately determine L e in the presence of thermal changes to R e .
- the use of a carrier signal at 150 kHz will significantly reduce the thermal effects upon L e measurement.
- DETAILED DESCRIPTION 7 C THEORY - Parasitic Capacitance and CANT Dynamics An important aspect of the present invention is described in the context of a digital control system that linearizes audio reproduction using a position-indicator state variable, f(x) , which is monotonic in position.
- the parasitic capacitance C parasitic between the voice coil and the body of a transducer can be used to give such a position-indicator state variable. This method applies to many other classes of nonlinear actuators and motors.
- the parasitic capacitance C parasitic between the voice coil of a transducer and the body of the transducer is largely determined by the relative positions of the voice coil and the magnetic pole pieces and central core.
- this capacitance with position is relatively straightforward and robust (reproducible).
- the voice coil 303 fits about a central core 310 that is part of the iron assembly 305.
- the variation in the parasitic capacitance depends largely on the overlap of the voice coil 303 with the central core 310, and, to some extent, with the outer pole piece 311 as well. More precisely, the parasitic capacitance is between the voice coil-copper wire and the entire magnetic circuit, each regarded as a single, equipotential, electrical conductor. C parasitit .
- C parasitic is an easily measurable internal circuit parameter of the transducer which is, at the same time, a state variable which depends monotonically upon axial coil position x .
- C parasitic capacitance As the coil moves deeper into the magnetic airgap, the capacitive contact areas between the metallic surfaces of coil and poles on the one hand, and between former and poles on the other hand, increases; and thus so does the value of the parasitic capacitance.
- Detailed measurements of C parasitic have been made as a function of x for the transducer of the Labtec Spin70 speaker, the large signal parameters of which are given by the curves depicted in FIG. 4. This transducer is of the type shown in cross- section in FIG. 3.
- the C parasitic measurements made were of two types: driven and non- driven. In the non-driven class of measurements, the voice coil was not driven, i.e.
- FIG. 32 shows the functional relation C ' para i t i c ix) for the mechanically moved, non-driven set of measurements.
- C parasitic is measured in arbitrary units obtained using the method described in Detail 12. While it is not possible to directly compare FIG. 32 and FIG. 33, it is known that V cail is monotonic with . It will be appreciated that the qualitative behaviors of the two curves agree for x values corresponding to a coil displaced outward from its equilibrium position. For the lower portion of the x range, however, the voltage driven variation in C parasitic (x) function is no longer monotonic. As illustrated in FIG.
- C parasith ( ⁇ ) curve displayed in FIG. 32 This lower portion of the position range corresponds to a coil/diaphragm assembly at mechanical equilibrium or displaced inward from equilibrium.
- the non-monotonicity in C p ⁇ msitic (x) displayed in FIG. 33 is understood as resulting from canting of the coil/diaphragm assembly as it moves into the airgap; the canting, in turn, results from magnetic torques on the incomplete wire-turns that terminate the coil solenoid on its outward end.
- transducer models used in this disclosure typically assume perfect azimuthal symmetry (i.e. invariance under rotations about the axis of symmetry) of both5 the transducer's geometry and its dynamics; this assumption is also made in most prior art models.
- azimuthal symmetry i.e. invariance under rotations about the axis of symmetry
- cant tilt
- FIG 35 shows the voice-coil and magnetic assembly for a canting transducer 300 in more detail.
- FIG. 35 illustrates the tilted voice coil 303, showing its dimensions h and r and variable tilt-angle ⁇ (mechanical connection of the coil to former and diaphragm assembly not shown), the core pole 310, and outer pole 311, both made of low-carbon steel in the case of the Labtec and similar speakers (typically 1008 or 1010 steel), and a permanent magnet 304 (sometimes one of several permanent magnets in the magnetic assembly).
- ⁇ denotes torque
- i(t) is the coil current, time independent in the DC case
- R/(x) is the transducer motor factor
- N the total number of windings in the voice coil.
- This magnetic torque is opposed by an elastic torque, caused by the elastic restoring forces acting to counter the canting.
- l ⁇ 2 p el ⁇ stic ( ⁇ ) K(x) the relevant torsional spring constant - i.e.
- h is the coil's height (defined above equation (48))
- K(x) is the coil/diaphragm suspension stiffness recognized in prior art
- p el ⁇ stic (x) is a dimensionless elastic ratio modulus , characteristic of the coil/diaphragm assembly.
- the p el ⁇ stic ratio modulus is expected to be significantly larger than unity, as speaker diaphragms are designed to resist canting while allowing axial motion.
- the elastic restoring torque is simply: el a stic * ⁇ Pel st i ⁇ x) K(x) ⁇ (t) (51 ) where ⁇ (t) represents the canting (or tilt) angle, in radian units, as a function of time.
- ⁇ (t) represents the canting (or tilt) angle, in radian units, as a function of time.
- equations (52)-(53) now yield the predicted fractional increase in parasitic capacitance due to canting: _ ° ⁇ parasitic ⁇ ,-, , , y 1 Bl(x) C paras i t i c ""** ⁇ °" ⁇ h * S Q R e N Pdastic (x) K(x)
- equation (54) only holds when the voltage V coil is of the sign corresponding to an inward magnetic Lorentz force acting on the coil; when V coil has the opposite sign, the fractional winding is too far from the airgap's magnetic field to result in significant canting, and ⁇ C parasitic becomes approximately zero.
- speckle-removing filters create signal delay.
- the bandwidth of the Klippel GMBH laser-based metrology system is on the order 1 kHz , which is too low for controlling a mid-range audio transducer.
- FIG. 37 illustrates the detection system 14200.
- An IR-LED 14201 , and a PIN diode detector 14202 are secured to a transducer frame 14203.
- a region 14204, consisting of reflective material or coating, such as white paint, is sprayed or placed on the back side of the transducer cone 14205.
- the IR-LED 14201 illuminates reflecting region 14204 with infrared light 14206.
- the electrical resistance of PIN diode detector 14202 changes with the position- dependent intensity variations of infrared light scattered from reflecting region 14204 on the back side of cone 14205. Due to the use of an area illuminator with finite emittance, a relatively widely illuminated region, and a finite-area detector with finite acceptance angle, the position information derived via this IR-LED method is quite linear with x over most of the cone's excursion.
- the IR-LED derived positional measure f(x) can be calibrated by comparing LED measurements against the laser output from a metrology instrument such as Klippel GMBH.
- FIG. 38 shows.a block diagram of a more specific embodiment of the generalized control system shown in FIG. 8.
- a DSP based controller 10101 consists of a DSP processor and software system 10102 and an interface system 10103 consisting of analog input/output and user interface software. Audio input is provided to DSP based controller 10101 through a signal-matching network 10104 that filters the audio input and provides the correct level of input to the interface system 10103. The audio input is acted on by the control routines in the DSP based controller 10101 and is output to a second signal-matching network 10105.
- the signal from the signal-matching network 10105 is provided to a power amplifier 10106.
- the output of power amplifier 10106 drives a speaker transducer 10107.
- a position sensor 10108, or sensors, is used to provide a position indication signal, indicating the position of the coil/diaphragm assembly of the speaker transducer 10107 to sensor signal conditioner 10109.
- Such position sensors could be, for example, the Z e detector of Detail 6, or IR detector described in Details 8 and 13, or C detector described in Details 7 and 12.
- Sensor signal conditioning system 10109 is used to amplify and filter the positional signal and match it to the level required for the interface system 10103.
- FIG. 39 is a block diagram of a particular embodiment of an audio reproduction system 15100 that includes a DSP based controller 10101.
- a Personal Computer (PC) 15101 which could be an eMachines T1742, is used as a control and user input environment for the DSP based controller 10101.
- the DSP based controller 10101 is implemented using a M67 DSP board 15102 and a A4D4 I/O board 15103 both manufactured by Alternative Integration Inc. (Simi Valley, CA).
- the M67 DSP board 15102 is a motherboard for the A4D4 I/O board 15103.
- the M67 DSP board 15102 contains a 106MHz TMS320C6701 floating point DSP manufactured by Texas Instruments and has been modified to add an inverter (74LS14) between JP14 pin 34 to
- the A4D4 I/O board 15103 consists of four 16 bit analog-to-digital converters (ADCs) and four 16 bit digital-to-analog converters (DACs) with interface circuitry to the M67 DSP board 15102.
- a Lynx L22 card 15104 manufactured by Lynx Studio Technology, Inc (Newport Beach, CA) installed on the PC 15101 provides an audio signal 15105 that is input to the A4D4 I/O board 15103.
- the Lynx L22 card 15104 receives input via Cool Edit Pro software 15106 (version 2) installed on PC 15101.
- the Cool Edit Pro software 15106 generates a '.wav' type digital sound file from a music source, which could be a CD player 15107 also installed on the PC 15101.
- the corrected analog audio signal 15108 is output from the A4D4 I/O board 15103, and provided as an input to a 20:1 attenuator
- Output from the attenuator 15109 is provided as input to a Marchand PM224 amplifier 15110 with internal jumpers set to give a DC coupled amplifier.
- the Marchand PM224 amplifier 15110 is manufactured by Marchand Electronics Inc (Webster NY).
- the Marchand PM224 amplifier 15110 is used to drive a 3" transducer 15111 manufactured by Audax (Westlake Village, CA).
- the embodiment of audio reproduction system shown in FIG. 39 uses the IR method of position sensing.
- An IR detector 15112 the operation of which is described in Details 8 and 13, is used both to measure the position of the coil/diaphragm assembly of the 3" transducer 15111 and to match the signal to the input stage of the A4D4 I/O board 15103.
- the output 15113 of the IR detector 15112 is an input to the A4D4 I/O board 15103.
- FIG. 40 shows the process flow used to linearize the transconductance component of the signal conditioning process and the transduction process of a given audio transducer, based upon the control model given by equations (36)-(39) in Detail 2 above.
- FIG. 40 applies also for the case that only a subset of these corrections are applied.
- the first step 111001 entails measuring large signal (LS) transducer parameters. This step yields coefficients of polynomial interpolations for the functions R (x) and L e (x) .
- the measurements are performed using a Klippel GMBH laser metrology system, with procedure as detailed in Klippel System Manual dated May 2, 2002.
- a software control program is invoked.
- a third step 111002 a software control program is invoked.
- the invoked software control program is run in 'Calibrate' mode in order to calibrate the functional relation between coil/diaphragm position x and the position- indicator nonlinear state variable, f(x) , which in one embodiment of the present invention is the voltage output of the IR circuitry: x ir - f(x) .
- the software control program collects corresponding values of x as measured to an approximation by the Klippel laser, and f(x) , in relation to the corresponding values of voltage outputs as described in Detail 5 so that the dependence of f(x) with x and the dependence of S with f(x) can be determined.
- step 111003 An example of the software control program used in step 111003 is provided by FIG. 41, FIG.42, FIG.43, and FIG. 44.
- the data obtained from steps 111001 and 111003 are used to find Best Fit coefficients for lowest order polynomials of S , x , Bl and L e as functions of x ir , as indicated by step 111004.
- 'Best Fit' is defined as that curve which is of the lowest order and which does not exceed specified rms and maximum errors, subject to substantial weighting in the mid section of the range of the (x) variable. More details and specifics on 'Best Fit' are provided later in this section.
- the user then inserts the polynomial coefficients obtained from step 111004 into the Software Control Program - step 111005.
- FIG. 41 shows the structure of one embodiment of the Software Control Program that is used both for obtaining data during calibration 111003, and for operating in normal mode 111007 in which linearized sound is produced.
- the initialization process 111101 places the system in a known state.
- the software control system can then be selected to operate in calibration mode 111103, which consists of an S and an x calibration process, or to operate in the normal mode 111104.
- calibration mode 111103 which consists of an S and an x calibration process, or to operate in the normal mode 111104.
- the first time around the user needs to select the calibration mode 111103, as indicated in 111003.
- the system can be selected for normal operating mode 111104, in which the software controls the sound reproduction process through an Interrupt Service Routine (ISR) 111106.
- ISR Interrupt Service Routine
- the system stops the program 111107.
- FIG. 45, FIG. 46, FIG. 47 and FIG. 48 cover the normal operation in detail, while FIG. 42, FIG. 43 and FIG. 44 cover the calibration mode in detail; all these figures are described later in this section.
- FIG. 49, FIG. 50 and FIG. 51 show the process of obtaining Best Fit Coefficients for S , x , Bl , and L e .
- operation 111204 the 'polyfit' function supplied with Matlab is utilized in order to fit two polynomials: S and x , each a different polynomial function of the corresponding position-indicator variable, f(x) .
- Bl and L e are provided from step 111001 as functions of the corresponding laser measurement x rather than as functions of x ir
- operation 111205 entails composing the functional relationships Bl and L e with the function f ⁇ l to yield the functions RJ ° f ⁇ l and L e ⁇ - f ⁇ l respectively, in accordance with the notation introduced in Detail 5 above.
- each of the monomial terms in the high-order polynomial approximant for that parameter is checked to see whether its maximal absolute value can exceed the tolerance divided by a significance factor such as ten. Those monomial terms which can exceed this bound in absolute value, are retained; while those that cannot exceed it, are discarded. This procedure results in a significant reduction in the order of the polynomial approximations to Bi f 1 and L e o f ⁇ l , especially for the former function.
- step 111202 an attempt is made, using the 'Best Fit' approach, to reduce the orders of all 4 polynomials: S , x , Bl , and L e .
- the approach is to specify a given amount of root mean square (rms ) error, and a corresponding maximum amount of error 111207, and then to run the 'Best Fit' polynomial order reduction program 111208 so as to fit polynomials of the lowest order possible without exceeding the specified errors.
- the polynomial coefficients are initialized to those obtained from the operation 111206.
- step 111210 coefficients are chosen from one of the following sets. Six rms error values were run: 0.1%, 0.3%, 0.5%, 1.0%, 3.0% and 5.0%. And for each of these rms error values, the maximum desired value was set at 5 times the rms value. The results for the case with rms error 1.0% was chosen as a compromise between low error magnitude and low online computation requirements: smaller values of errors yield higher orders of coefficients, which require a higher amount of online computation. 5 FIG.
- the DSP software in program 111208 provides details of the operations performed by the DSP software in program 111208 in order to reduce the order of the approximate polynomial interpolating functions for S , x , Bl , and L e as functions of x ir for the specified rms and maximum error values, while maintaining 'Best Fit'.
- the user specifies the range of f(x) , the midsection of f(x) , and a weight for this o midsection.
- the following set was chosen, in units of volts for the IR circuit output voltage: a range of [-0.8 to 0.8]; midsection [-0.3 to 0.3]; and a weight of 10 for the midsection, with the rest of the range being assigned the weight of 1.
- the high weight value (10) chosen for the midsection was motivated by the need to accommodate three requirements: (a) to emphasize a better fit in this5 predominantly linear section; (b) to account for the fact that the outer section is much larger; and (c) to account for the fact that there are more points in the outer sections than indicated by mere proportion, due to more predominant nonlinearity of S in the outer section (since coil DC voltage, rather than position x , was ramped in equal step- sizes, as shown e.g. in FIG. 22).
- This weighting results in a better fit in the linear region 0 compared with the fit obtained by a non-weighted approach.
- someone skilled in the art will recognize that other choices for range, midsection and weights are possible within the framework of this invention.
- Step 111302 is a programming maintenance function (file name specifications).
- step 111303 the operations for polynomial order reduction are repeated for S(x ir ) ,5 x(x ir ) , Bl(x ir ) and L e (x ir ) , with a reduced set of coefficients determined one curve at a time.
- S the first curve for polynomial reduction, although the process could have equally well began with x , Bl , or L e with identical overall results.
- the coefficients for the next curve are supplied 111304.0
- the operations within step 111305 are detailed in FIG. 51.
- the above Y orig (p) values are then used in step 111403 to compute new coefficients, and in module 111404 to compute errors.
- Y orig values Y orig . , Y o g2 , ...
- Y orig33 are calculated for 33 points ⁇ l , ..., ⁇ 33 distributed uniformly over the above range. It will be readily recognized that that the number of points used can be changed within the framework of this invention.
- the 'Best Fit' coefficients are computed as described here and based on a weighted least-squares curve fitting approach used in signal processing [P.M. Embree and Damon Danieli, C++ Algorithms for Digital Signal Processing, Second Ed., 1999, Prentice Hall].
- a Jk W j * p , where j is the data point index, k is the power index, and w, is the weight for the point p j .
- j ranges from 1 to N , where N is the number of data points chosen over the range above, while k ranges from 0 through M , with M being the reduced order for which best fit coefficients are being derived.
- the data point index starts with 1
- the power or order index starts with 0.
- the error between Y orig and Y new is computed, squared, and weighted by corresponding weights.
- the total is divided by a weighted divisor, i.e a number obtained by taking the total points in the mid section, multiplying it by 10, and adding to it the number of data points outside of the mid section. Taking the square root of the divided result yields the rms value.
- the maximum magnitude of error between points of Y orig and Y new is also determined in step 111404.
- step 111405 For the error test in step 111405, if either the rms error or the maximum magnitude error exceeds corresponding specified value, the control goes to the 'Yes' branch; else it goes to the 'No' branch, to reduce the order further (step 111402) by repeating the above process.
- step 111406 checks whether the polynomial order has been reduced; only if the answer is 'Yes' on this latter test, does the program declare 'Pass' and output the lowest order b vector that had both the rms and maximum magnitude not exceeding corresponding specified values. Otherwise, it declares 'Fail' and outputs the original coefficients.
- the program passes control to the calling program 111306 which tests if any more curves need to be processed for reduction of order while obtaining 'Best Fit'.
- the steps of FIG. 50 and FIG. 51 have been implemented in a program, written in Matlab, which uses the function developed by Tymphany Corporation for module
- FIG. 45, FIG. 46, FIG. 47 and FIG. 48 cover the normal mode of operation 111104.
- FIG. 45 shows an overall flow diagram of normal mode of operation 111104.
- an initialization process 11201 receives the user inputs such as the sampling frequency and the initial audio volume level.
- Step 11201 initializes the Digital-to-Analog converter (DAC), enables Analog-to- Digital converter (ADC) and DAC triggers, and initializes and sets up the ISR 11203.
- Step 11202 enables the ISR, sets the sampling rate of the real time clock, and enables it. The enabling of the sampling clock spawns the process: execute normal mode HW & ISR operations 11203. The software then enters a wait loop and command parser
- FIG. 46 shows the operations of process 11203 that are spawned as a result of enabling the sampling clock and ISR in 11202. These elements are spawned in parallel with the mainline operation. Note that the three processes: Sampling Clock 11301 , ADC Convert 11302, and the ISR 11303 are activated essentially in parallel. However, ADC convert 11302 starts on the rising edge of sampling clock 11301 , while ISR 11303 starts on the falling edge of the sampling clock 11301. Moreover, when the falling edge of sampling clock 11301 occurs, the ISR 11303 uses the most recently converted sample from ADC convert 11302.
- the Sampling Clock 11301 is typically set at 48 kHz, although any frequency above the Nyquist frequency for audio (typically above 40kHz) can be chosen.
- the sampling clock 11301 runs as an autonomous hardware loop, operating until powered down, or disabled by the software control program.
- the ADC Convert module 11302 samples and converts an analog stream representing the sensor measurement of the position-indicator state variable and the audio source.
- the ISR 11303 operates on the converted data provided by ADC convert 11302.
- FIG. 47 shows a flow diagram of the ISR 11303.
- the software control passes from the wait loop and command parser 11204 to step 11401.
- Step 11401 limits the value of the word to be sent to the
- the DAC can be an onboard DAC as it is with the innovative Integration A4D4, or a serial-port-based off-board DAC.
- the analog signal that is created is the corrected audio signal V coil , and is fed to a power amplifier 10106.
- the ISR module 11303 uses IR sensor data f(x) from module 11403 and audio data from module
- Module 11405 A digital filter 11404 is used to minimize sensor noise in the measurement of f(x) .
- Module 11406 computes S , B , and L e corrections from the filtered value of f(x) 11404, as described below.
- the input f(x) read from ADC in module 11403 is scaled to volts by dividing the value of f(x) by 3,276.7.
- the divisor 3,276.7 was chosen because of the DAC resolution.
- the onboard DACs of the innovative Integration M67 are
- 32767 counts / 10 volts. If an off board 1V DAC is used, the divisor would be 32,767
- S S 0 + S X f ⁇ x) + S 2 (f(x)) 2 + - + S 5 (f(x)) 5 - k f- 1 (f(x))
- Equation (61) is a digital filter designed to estimate the velocity of the transducer needed for the BEMF correction.
- Equation (62) calculates the required BEMF correction.
- the BEMF correction comprises two components: the removal of the nonlinear BEMF and the replacement with a linear BEMF.
- the equations incorporate a multiplier for each term to allow for fine adjustment of the correction.
- Equation (63) and (64) implements the above components of the audio correction. It will be appreciated that there are many different ways of discretizing the numerical differentiation operation of the control diagrams FIG. 11 and FIG. 12, and that the implementation of these numerical differentiations used in one embodiment of the invention, and shown in equations (61) and (64), represent but one possible choice.
- Digital filters may be added to equation (64) for smoothing, equalizing and noise reduction.
- module 11410 executes a return from ISR, which passes the software control to the wait loop 11204; and the process then repeats, unless stopped by a 'Stop' command to the wait loop 11204 which resides in the normal mode 111104.
- FIG. 42, FIG. 43 and FIG. 44 show flow diagrams of S and x versus f(x) calibration 111103.
- FIG. 42, FIG. 43 and FIG. 44 illustrate S and x calibration 111103.
- FIG. 42 shows the overall flow diagram of S and x versus f(x) calibration.
- An array is initialized with S values that will be used as S drive for calibration.
- the magnitude of the S drive should be large enough to drive the transducer close to its maximal and minimal x excursions.
- the operations in FIG. 42 are similar to those in FIG. 45.
- the diagram shows the mainline S calibration loop 11505.
- FIG. 43 shows the details of HW and ISR operations for S calibration 11504. It depicts Sampling Clock 11601 and ADC Convert 11602, which are similar to corresponding modules in FIG. 46; the same description applies, and is thus not repeated.
- Modules 11604 and 11605 limit and convert the digital values to analog waveform.
- Module 11606 tests whether the data is to be collected.
- the mainline S calibration loop 11505 sets and clears the flag 'Collect_data'. If this flag is set, the data collection is done by the module 11607, and a sample count is tallied.
- module 11608 reads the S value from the array, to be used in the variable 'dacvalue'. If the flag is not set, these two modules are bypassed.
- Module 11609 executes the return from ISR.
- FIG. 44 shows the details of mainline S calibration loop 11505.
- Module 11701 checks whether any value of S is left with which to operate the loop. If there is one, it executes the path comprising modules 11702 through 11707 to send out the S value via the ISR 11603, and to collect the corresponding value of f(x) and x as follows. Module 11702 executes a wait of 100 milliseconds to allow the transients in the transducer to attenuate.
- Module 11703 sets the 'Collect_data' flag which signals the ISR 11603 to collect data.
- Module 11704 allows 1 millisecond to collect samples, which at 48 kSPS collects 48 samples. These samples suffice to give a good reading of f(x) , the IR data, and x, the laser data.
- Module 11706 performs averaging, and the module 11707 stores S , f(x) and x for offline curve fitting. As long as there is an S value to be covered, the process continues. To ensure reliable calibration, the values in the arrays are such that each point of S is covered at least twice, each at very different instances of time.
- the calibration of S is started at 0, and increased in steps until an upper limit is reached, and then decreased in steps until a lower (negative) limit is reached. Again it is increased until the top limit is reached. From the top limit, it is decreased in steps until the negative limit is reached. From the negative value, S is increased in steps until it returns to 0. Thus it forms a W pattern.
- the mainline loop for S commences a termination procedure, as shown in the module 11708.
- the sampling clock is disabled, which stops the operations of ADC convert 11602 and the ISR 11603.
- FIG. 48 illustrates the details of the Wait Loop and Command Parser 11204, shown in FIG. 45, which is abbreviated below as WLCP.
- the system enters into the 11801 step of WLCP from Enable ISR Setup and Enable Sampling Clock 11202; in step
- step 11801 it is determined whether Normal Mode operation should stop. If 'Yes', system enters into step 11803, in which Interrupt is disabled and the HW is put into a known state; then system is passed out of WLCP and into User Mode Select 111102. But if the answer to the 'Stop?' query (step 11801) is 'No', the DSP passes to step 'Command?' 11802, in which the WLCP checks to see whether User has entered a keyboard command since the last check (checks are spaced several microseconds apart during the Wait Loop).
- User keyboard response 'b' causes the DSP to enter the mode 'Adjust Linear BEMF' 11805, from which it is again returned to 'Stop?' query 11801.
- the following are the remainder of the allowed keyboard responses, and their effects.
- Response '+' puts the DSP into mode 'Increase Volume' 11806, from which it returns to 'Stop?' query 11801 ; similarly, response '-' puts DSP into mode 'Decrease Volume' 11809, and thence to 'Stop?' query 11801.
- Response 'u' puts DSP into 'Uncorrected Audio Mode' 11807, and thence to 'Stop?' query 11801.
- Response T puts DSP into mode 'Adjust dL/dx Correction' 11808, and thence to 'Stop?' query 11801.
- Response 'o' puts DSP into mode 'Adjust Offset' 11810, and thence to 'Stop?' query 11801.
- Response 'j' puts DSP into mode 'Adjust dL/dx Offset' 11811 , and thence to 'Stop?' query 11801.
- Response 'm' puts DSP into mode 'Mute On' 11812, and thence to 'Stop?' query 11801.
- Response 'k' puts DSP into mode 'Adjust Linear Spring' 11813, and thence to 'Stop?' query 11801.
- Response T puts DSP into mode 'Turn IR
- DETAILED DESCRIPTION 11 Z. METHODS AND CIRCUITS
- the present invention is described, in one aspect, in the context of controlling an audio reproduction system, in part by a system, consisting of methods and electronic circuits, which provide at least one position-indicator transducer state variable derived from effective circuit parameters of the transducer during operation.
- the position-indicator state variable f(x) utilized in this embodiment of the invention is an output voltage derived from the functional dependence of the effective complex coil impedance Z e ( ⁇ ,x) upon coil/diaphragm position , at some fixed supersonic probe frequency ⁇ .
- FIG. 52 shows a block diagram of a potential divider circuit 12100.
- An exciting signal, a probe tone 12101 at a fixed frequency and fixed amplitude, is connected across a potential divider consisting of the transducer voice coil 12102 and a reference impedance Z re/ 12103.
- the magnitude of the output voltage 12104 across the reference impedance 12103 is a fraction of the magnitude of probe tone voltage 12101, depending on the relative impedances of the transducer voice coil 12102 and the reference impedance 12103.
- the input signal to the voice coil will include audio information (program material) together with the probe tone. It is therefore necessary to separate the probe tone and program material in frequency, so that the probe tone measurement is not interfered with by the audio drive signal.
- the Nyquist criterion suggests that the probe tone 12101 should have a frequency of at least twice the audio frequency bandwidth, to avoid aliasing with the program material.
- a probe tone having a frequency of 43kHz has been found to be particularly desirable. However, many other frequency values could be used.
- a desirable implementation utilizes a potential divider measurement system that is filtered to separate out the contributions of the audio program material and of the ultrasonic probe tone frequency.
- the filtered probe tone 12101 is then envelope-detected and reduced to an audio frequency signal, which varies as Z e (x) changes due to the voice-coil motion created by the transducer in response to the audio input signal.
- FIG. 53 shows a block diagram of the Z e (x) detection system 12200.
- the probe tone 12101 is added to the audio drive signal 12201 in a summing circuit 12202.
- the summed signal excites a potential divider 12203, which includes the transducer voice coil 12102.
- FIG. 54 shows a block diagram of a control circuit for transducer linearization, which includes the Z e (x) detection circuit 12200 (FIG. 53).
- An incoming audio signal 12301 is converted into digital form and input to a DSP, for example, using the mixed signal device 12302 which may be, for example, implemented by an Analog
- the Z e (x) signal 12207 is also provided as an input to and converted by the mixed signal device 12302.
- the DSP core runs the linearization algorithm, with the Z e (x) signal 12207 as the positional signal.
- the corrected audio signal 12305 is an input signal to amplifier 12303, which produces the audio drive signal
- the sine wave generator 12304 preferably has a low impedance output, for example below 1.0 Ohm.
- FIG. 55 shows a circuit diagram of the summing circuit 12202.
- the audio drive signal 12201 is provided as an input to filter 12401 which isolates the probe tone 12101 from the low impedance of the audio amplifier output.
- the filter 12401 is composed of resistive, capacitive, and inductive elements, as indicated in FIG. 55.
- the probe tone 12101 is provided to a capacitor 124C4, which in turn is connected to the summing point 12402. Capacitor 124C4 decouples the audio drive signal at the summing point 12402 from the low impedance output of the sine wave generator 12304.
- the signal at the summing point 12402 is provided at output terminal 12403 which is connected to an input of the potential divider circuit 12203.
- FIG. 56 shows the circuits of the potential divider 12203 and the high pass filter 12204.
- the summed output 12403 excites the potential divider 12203, which includes the voice coil 12501 of the transducer being used in the audio system, and a reference inductor 12502.
- the proportional excitation across the reference inductor 12502 is input to the capacitor 125C1 of the high pass filter 12204.
- the high pass filter 12204 may be, for example, a standard 2nd order Butterworth filter, designed to discriminate against the audio signal and pass the 43kHz probe tone.
- Operational amplifier 12504 may be, for example, a National Semiconductor part LM741.
- the filter has as its output the filtered 43kHz signal 12503.
- FIG. 57 shows the circuit of the full wave bridge detector circuit 12205. This is a standard circuit that rectifies the filtered 43kHz signal 12503 and outputs a full wave rectified signal 12601.
- Operational amplifiers 1260A1 and 1260A2 may be implemented by National Semiconductor part LM741 devices.
- One skilled in the art will recognize that many different circuit arrangements could be used for the full wave bridge detection circuit 12205 and that the standard circuit shown here is only one example.
- FIG. 58 shows the circuit of the low pass filter 12206.
- the first part of the low pass filter incorporating the operational amplifier 1270A1
- the second part of the filter is an inverting amplifying stage that includes an operational amplifier 1270A2 and a variable resistance 127VR1 that produces a DC offset in the output signal. This offset is set to reduce the DC offset in the magnitude of the probe tone that has been detected.
- the gain of the inverting amplifying stage is set to enhance the signal significance when it is converted to digital form.
- One skilled in the art will recognize that many different circuit arrangements could be used for the filter, gain and offset circuit, and that the rather straightforward circuit shown in FIG. 58 can be modified without changing the essence of the design.
- Operational amplifiers 1270A1 and 1270A2 may be National Semiconductor part number LM741.
- FIG. 59 shows the circuit of the audio amplifier 12303 of FIG. 54 in more detail.
- the corrected audio signal 12305 received from the ADI-21992 EZ-KIT 12302 is a positive unipolar signal and must be offset to a signal oscillating about zero for output as audio.
- the requisite offset is achieved by utilizing an inverting operational amplifier 1280A1, which may be, for example, a National Semiconductor part LM741, in a unity gain stage, with offset provided by a variable resistor 128VR1 connected to a positive voltage.
- a power operational amplifier 1280A2 for example a National Semiconductor part LM575, is used to amplify the corrected audio signal 12305 and drive the speaker
- FIG. 54 is sensitive to changes in output impedance of the audio amplifier 12303.
- some types of amplifiers exhibit large crossover distortion effects, which in effect are a change in output impedance.
- This change in output impedance can cause noise in the Z e (x) measurement.
- transducers driven with large currents there can be considerable heating effects in the coil.
- someone skilled in the art would recognize that a more complex circuit is required to separate out these two effects for the full range of transducers, but that this would not materially change the invention detailed here.
- DETAILED DESCRIPTION 12 C METHODS AND CIRCUITS
- the present invention is described, in one aspect, in the context of controlling an audio reproduction system, in part by a system, consisting of methods and electronic circuits, which provide at least one position-indicator transducer state variable derived from effective circuit parameters of the transducer during operation.
- the position-indicator state variable, f(x) utilized in this embodiment of the invention is an output voltage derived from the internal parasitic capacitance C parasitic between the transducer voice-coil and the transducer magnetic pole structure.
- the method utilizes the functional dependence C parasitic (x) of this capacitance upon the axial position of the transducer's coil/diaphragm assembly as a positional sensor.
- FIG. 34 shows a schematic cross section of a typical cell phone speaker or receiver 13100; actual three-dimensional speaker geometry is a figure of revolution about the central horizontal axis of symmetry (not shown).
- Speaker 13100 consists of a transducer and integral acoustic venting.
- a voice coil 13101 is mounted on the diaphragm 13102.
- Coil 13101 is positioned in the gap between a neodymium magnet 13103 and a magnetic base plate 13104.
- a plastic surround 13105 supports the diaphragm 13102 and a faceplate 13106.
- the surround and faceplate have acoustic vents 13107 which tune the frequency response of the speaker 13100.
- the depth, indicated in FIG. 34 by D1 is typically 2 mm.
- the main difference between this type of transducer assembly and the transducer shown in FIG. 3 is the single surrounding support of the relatively flat diaphragm 13102. This means that the system is resistant to the tilt ("canting") that can complicate capacitance position-sensing methods in other transducers as described in Detail 7.
- the preferred method of detecting the variation of capacitance with coil/diaphragm axial position, C pamsiiic (x) is to place the capacitance within an oscillator circuit.
- FIG. 60 shows a schematic of the capacitance detector and speaker arrangement, together with the DSP used for correction.
- An analog audio signal, provided over input line 13201 is digitized by DSP based mixed-signal controller 13202.
- Mixed-signal controller 13202 is embodied by a AD21992 chip which includes an ADC (analog to digital converter).
- the output of the DSP based controller 13202 is connected to a standard DAC (digital to analog converter) 13203.
- the output of the DAC 13203 is amplified by a DC-connected audio amplifier 13204.
- the output of amplifier 13204 has a drive connection 13205a to one terminal of the voice coil 13101 of the speaker 13100.
- the magnetic base plate 13104 of the speaker 13100 has a connection 13207 to one input of an oscillator circuit 13208 (detailed in FIG. 61). Another input to the oscillator circuit 13208 is connected to the drive connections 13205a and 13205b of the coil 13101 through blocking capacitors 13209a and 13209b, respectively.
- the output of oscillator circuit 13208 is connected to a frequency to voltage converter 13210, which converts the variable frequency received from the oscillator circuit 13208, and also amplifies and level-shifts the varying voltage output.
- the output 13404 from the frequency to voltage converter 13210 which is a measure of C peritic (x) (abbreviated as C p (x) in the Figure), and hence the position-indicator state variable f(x) , is input into the mixed signal DSP controller 13202.
- C p (x) abbreviated as C p (x) in the Figure
- f(x) the position-indicator state variable
- the purpose of the DSP functionality within the controller 13202 is to furnish the DAC 13203 with a digital signal such that the output of DAC 13203, after amplification by amplifier 13204, will feed the speaker-transducer voice coil with a voltage signal including both the audio program and a pre-distortion calculated to cancel out a significant portion of the nonlinearities introduced by the transducer in the course of its normal uncorrected operation.
- FIG. 61 shows the input from speaker 13100 and the detail of the oscillator circuit 13208.
- the audio amplifier drive signal connections 13205a and 13205b are decoupled using 60pF capacitors 13209a and 13209b connected to the ground of the oscillator circuit 13208.
- the parasitic capacitance between the voice coil 13101 and the base plate 13104 is part of the R C oscillator created by the circuit, with the resistance values shown and an LF411 operational amplifier 13303 (available, for example, from National Semiconductor).
- the parasitic capacitance between the voice coil 13101 (FIG.
- the base plate 13104 and the base plate 13104 is part of the R C oscillator created by the circuit, with the resistance values shown and an LF411 operational amplifier 13303 (available, for example, from National Semiconductor).
- the electrical connection to magnetic base plate is indicated by reference character 13207.
- the values of the variable parasitic capacitance C p (x) typically ranges between 2pF and lO Efor the above-mentioned type of speaker, and thus the oscillator circuit must be physically close to the speaker to avoid the effects of environmental sources of further stray capacitance. Such further stray capacitance would reduce the sensitivity of the system.
- the oscillator output signal (at terminal 13304) is a square wave of varying frequency between I MHz and 2 MHz .
- FIG. 62 shows the detailed circuitry of the frequency to voltage converter 13210.
- Frequency to voltage converter 13210 consists of two parts: a frequency to pulse converter circuit 13401, and a low pass filter, amplifier and level shifter circuit 13402.
- the frequency to pulse converter 13401 consists of a mono-stable multi-vibrator circuit
- the mono-stable multi-vibrator circuit 13407 takes the square wave output signal 13304 received from the oscillator circuit 13208, which has a constant rms value, and converts it to a pulse train that is provided on line 13403.
- the pulse train 13403 has an rms value varying with frequency, which is a function of the transducer coil/core capacitance C , which in turn varies with coil/diaphragm position .
- the low pass filter, amplifier and level shifter circuit 13402 converts the pulse train on line 13403 to a varying analog voltage output provided on line 13404.
- the low pass filter, amplifier and level shifter circuit 13402 includes an operational amplifier 1340A1 , which receives the output signal on line 13403 and, using a gain of 10 as determined by resistor values, low-pass-filters and offsets the signal 13403; and operational amplifier 1340A2, which has a gain of unity and implements a second-order Butterworth filter.
- operational amplifiers may be embodied, for example, as National Semiconductor part number LM741 , or equivalent.
- Resistor 134VR1 is adjusted such that the coil/diaphragm equilibrium position produces a zero output voltage.
- Operational amplifier 1340A2 receives, at its input terminal 13406, the output signal provided at output terminal 13405 of operational amplifier 1340A1 , and then converts that signal to a voltage which is provided on line 13404 to mixed signal DSP 13202.
- the capacitance dependent voltage output 13404 is also a position sensitive signal (since C p depends on x ).
- the functional dependence C p (x) is monotonic, and C p can thus be used as a position-indicator nonlinear state variable in lieu of the position variable x itself in a feedback linearization control law.
- IR METHODS AND CIRCUITS The present invention is described, in one aspect, in the context of controlling an audio reproduction system, in part by a system, consisting of methods and electronic circuits, which provide at least one position-indicator transducer state variable.
- the position-indicator state variable, f(x) utilized in this embodiment of the invention is an output voltage from an optical IR-LED system, as discussed in Detail 8. This embodiment is called the IR method.
- FIG. 63 shows an overall block diagram of a system 14100 for implementing the IR-LED method for detecting a position-indicator state variable.
- IR light 14206 is emitted by an IR-LED 14201.
- the IR light 14106 is scattered off a reflecting region 14204 on the back side of the transducer cone.
- the scattered IR light 14104 is detected by a PIN diode detector 14202.
- a detection circuit 14106 supplies current to the IR-LED 14201 and detects the photo-current flowing in the PIN diode 14202.
- the electronic circuit 14106 converts the photo-current flowing in the PIN diode 14202 to a positional signal, the present value of the position-indicator transducer state variable f(x) 14107.
- FIG. 64 shows an embodiment of the circuit schematic of IR-LED detection circuit 14106 of FIG. 63.
- the IR-LED 14201 and PIN diode 14202 are both connected into the circuit with a short (less than 1 meter) shielded cable (not shown) that extends from the circuit board which includes the remaining electronics to the frame 14203 of the transducer on which the IR-LED 14201 and PIN diode 14202 are supported.
- the IR- LED 14201 may be implemented by a SLI-0308CP purchased from Jameco Electronics in Belmont, California and PIN diode 14202 may be implemented by a IRD500 purchased from Jameco Electronics in Belmont, California.
- the detector configuration used in the IR-LED detection circuit 14106 is operated in the "reversed biased" mode of operation. In this mode the PIN diode 14202 is biased by an external direct voltage.
- this voltage is 6V , though it may be as high as 40 to 60V .
- the PIN diode 14202 operates as a leaky diode, with the leakage current depending upon the intensity of the light striking the device's active area.
- a silicon PIN diode of the type described above will typically leak nearly 1 mA of current per 2mW of light striking it, which constitutes a high quantum efficiency.
- Low cost IR LEDs of which the one mentioned above is an example, will produce sufficient power for this application. It should be noted that a PIN photodiode has both the speed and the sensitivity required for the position detection described herein, and is available at a low cost.
- the IR-LED detection circuit 14106 is configured as a transimpedance amplifier.
- Resistor 144R5 which converts the PIN diode 14202 current into a voltage is connected from the output to the input of an inverting operational amplifier 1440P1.
- the amplifier 1440P1 thus acts as a buffer, and produces an output voltage proportional to the PIN diode current.
- the zero balance meaning that the cone of the transducer is at the rest position, is set by a variable resistor 144VR2.
- the transimpedance amplifier 1440P1 is followed by another high gain amplifier 1440P2.
- a variable resistor 144VR3 is used to set the gain of the amplifier in order to match the input range of the A/D converter which receives the voltage f(x) , which in one embodiment was + 1.00 volt.
- the IR-LED 14201 and PIN diode 14202 are epoxied side-by-side onto the transducer frame 14203, with both diodes pointing at a reflecting region 14204 on the transducer cone 14205.
- Reflecting region 14204 should subtend a sufficient angle such that, as the transducer cone moves, the PIN diode 14202 detector admittance cone is always pointed within the region.
- the diodes are preferably inclined towards each other and pointed towards the axis of the transducer at approximately a right angle to the direction of motion, or towards the curve of the cone.
- the PIN diode output is not completely linear with cone position and therefore requires calibration by comparison with a metrology system.
- the position-indicator variable, f(x) and the degree of its non-linearity, can be varied by changing the positions and orientations of the two diodes relative to each other and to the transducer cone.
- the circuit 14400 is prone to saturation and interference from ambient light. Hence prior to operation the diodes must be shielded from external light, either by masking or by the speaker cabinet. All adjustable resistors in the circuit are put at the center of their resistive ranges.
- the circuit board is connected to the diodes with a shielded cable, and powered.
- the IR LED current resistor 144VR1 is adjusted until the output is approximately at ground potential.
- the transducer voice coil (not shown) is connected to a low power, low frequency AC source (for example, 20 -60 kHz ), and the power to the voice coil is adjusted to give maximal Peak-to-Peak motion, while avoiding excursions large enough to cause the cone to hit its encasement.
- a low power, low frequency AC source for example, 20 -60 kHz
- DETAILED DESCRIPTION 14 IR RESULTS
- FIG. 65 shows a portion near 3 kHz of the FFT power spectrum distribution of the SPL (sound pressure level) wave-pattern picked up by a microphone in the acoustic near-field; both corrected spectra which is indicated by reference character 1521, and uncorrected spectra which is indicated by reference character 1522 are shown, and it is clearly seen that the powers in the 60 Hz -spaced lattice of intermodulation frequency peaks, are significantly reduced when the correction is applied.
- FIG. 66 shows the low- frequency portion of the same power spectrum distribution, showing multiple harmonics of the 60Hz tone; again, spectra are depicted both with and without correction, and again, significant reduction in the magnitude of the harmonic distortion peaks can be seen.
Landscapes
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Control Of Linear Motors (AREA)
Abstract
Description
Claims
Applications Claiming Priority (11)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/638,201 US20060104451A1 (en) | 2003-08-07 | 2003-08-07 | Audio reproduction system |
US10/638,224 US20050031134A1 (en) | 2003-08-07 | 2003-08-07 | Position detection of an actuator using infrared light |
US10/638,232 US20050031139A1 (en) | 2003-08-07 | 2003-08-07 | Position detection of an actuator using impedance |
US10/638,243 US20050031140A1 (en) | 2003-08-07 | 2003-08-07 | Position detection of an actuator using a capacitance measurement |
US10/638,198 US20050031137A1 (en) | 2003-08-07 | 2003-08-07 | Calibration of an actuator |
US10/638,200 US20050031138A1 (en) | 2003-08-07 | 2003-08-07 | Method of measuring a cant of an actuator |
US10/637,960 US20050031131A1 (en) | 2003-08-07 | 2003-08-07 | Method of modifying dynamics of a system |
US10/638,197 US20050031133A1 (en) | 2003-08-07 | 2003-08-07 | Process for position indication |
US10/638,170 US20050031132A1 (en) | 2003-08-07 | 2003-08-07 | Control system |
US10/638,237 US20050031117A1 (en) | 2003-08-07 | 2003-08-07 | Audio reproduction system for telephony device |
PCT/US2004/025485 WO2005015955A2 (en) | 2003-08-07 | 2004-08-06 | Audio reproduction system |
Publications (1)
Publication Number | Publication Date |
---|---|
EP1665884A2 true EP1665884A2 (en) | 2006-06-07 |
Family
ID=34140020
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP04780339A Withdrawn EP1665884A2 (en) | 2003-08-07 | 2004-08-06 | Audio reproduction system |
Country Status (2)
Country | Link |
---|---|
EP (1) | EP1665884A2 (en) |
WO (1) | WO2005015955A2 (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101621728B (en) * | 2009-06-25 | 2013-03-06 | 北京卓锐微技术有限公司 | Method and device for calibrating sensitivity of microphone |
DE102018002290A1 (en) | 2017-03-27 | 2018-09-27 | Sound Solutions International Co., Ltd. | A system and method for applying a sound signal to a multi-coil electrodynamic acoustic transducer |
DE102018002289A1 (en) | 2017-03-27 | 2018-09-27 | Sound Solutions International Co., Ltd. | A method for avoiding a deviation of a diaphragm of an electrodynamic acoustic transducer |
CN110460945A (en) * | 2018-05-08 | 2019-11-15 | 西安中兴新软件有限责任公司 | A kind of audio channel detecting circuit and detection method |
CN111565353B (en) * | 2020-03-10 | 2021-05-28 | 南京大学 | A Loudspeaker Nonlinear Parameter Identification Method with Adaptive Multi-step Size |
CN114735020A (en) * | 2022-03-29 | 2022-07-12 | 重庆邮电大学 | Steering control method and device for audio environment interaction |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6104817A (en) * | 1996-12-12 | 2000-08-15 | Ding; Chih-Shun | Speaker and amplifier system |
-
2004
- 2004-08-06 WO PCT/US2004/025485 patent/WO2005015955A2/en active Application Filing
- 2004-08-06 EP EP04780339A patent/EP1665884A2/en not_active Withdrawn
Non-Patent Citations (1)
Title |
---|
See references of WO2005015955A3 * |
Also Published As
Publication number | Publication date |
---|---|
WO2005015955A2 (en) | 2005-02-17 |
WO2005015955A3 (en) | 2007-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20050031139A1 (en) | Position detection of an actuator using impedance | |
US20050031134A1 (en) | Position detection of an actuator using infrared light | |
US20060104451A1 (en) | Audio reproduction system | |
US20050031132A1 (en) | Control system | |
US20050031137A1 (en) | Calibration of an actuator | |
US20050031133A1 (en) | Process for position indication | |
US20050031140A1 (en) | Position detection of an actuator using a capacitance measurement | |
US20050031131A1 (en) | Method of modifying dynamics of a system | |
US20050031138A1 (en) | Method of measuring a cant of an actuator | |
US20050031117A1 (en) | Audio reproduction system for telephony device | |
KR101647315B1 (en) | Arrangement and method for converting an input signal into an output signal and for generating a predefined transfer behavior between said input signal and said output signal | |
EP3047657B1 (en) | Systems and methods for protecting a speaker from overexcursion | |
US10667040B1 (en) | System and method for compensating for non-linear behavior for an acoustic transducer based on magnetic flux | |
CN102843633A (en) | Control of loudspeaker output | |
CN111885475B (en) | System and method for compensating for nonlinear behavior of an acoustic transducer | |
Klippel | Measurement of large-signal parameters of electrodynamic transducer | |
US3988541A (en) | Method and apparatus for frequency compensation of electro-mechanical transducer | |
US11012786B2 (en) | Armature-based acoustic receiver having improved output and method | |
Sturtzer et al. | Comparison between voltage and current driving methods of a micro-speaker | |
Klippel | Diagnosis and remedy of nonlinearities in electrodynamical transducers | |
US20030118193A1 (en) | Method and system for digitally controlling a speaker | |
Klippel | Mechanical overload protection of loudspeaker systems | |
Schurer | Linearization of electroacoustic transducers | |
EP1665884A2 (en) | Audio reproduction system | |
EP3476135A1 (en) | Method for simulating total harmonic distortion of a loudspeaker |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
17P | Request for examination filed |
Effective date: 20060208 |
|
AK | Designated contracting states |
Kind code of ref document: A2 Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PL PT RO SE SI SK TR |
|
AX | Request for extension of the european patent |
Extension state: AL HR LT LV MK |
|
DAX | Request for extension of the european patent (deleted) | ||
PUAK | Availability of information related to the publication of the international search report |
Free format text: ORIGINAL CODE: 0009015 |
|
RIC1 | Information provided on ipc code assigned before grant |
Ipc: H04R 29/00 20060101ALI20070223BHEP Ipc: H04R 3/00 20060101AFI20070223BHEP |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION HAS BEEN WITHDRAWN |
|
18W | Application withdrawn |
Effective date: 20080625 |