Abstract:
The present invention discloses a method for parameter identification and parameter optimization of microspeakers. Measurement procedures for identifying electromechanical constants of microspeaker and a GUI are developed to facilitate estimation of electroacoustic parameters of the microspeaker under test. In light of the thus identified microspeaker parameters, a parameter optimization procedure is carried out to obtain the design that attains the best acoustic performance with minimum harmonic distortion.

Description:
BACKGROUND OF THE INVENTION 
       [0001]    1. Field of the Invention 
         [0002]    The present invention concerns a method for parameter identification and parameter optimization of a speaker, in particular a method for parameter identification and parameter optimization of a microspeaker. 
         [0003]    2. Description of the Related Art 
         [0004]    Micro speakers have been extensively used in electronic products recently as important components in mobile phones, digital cameras, personal digital assistants and MPEG3 display devices. In order to achieve the best performance and the minimum harmonic distortion of microspeakers, it is necessary to estimate relevant electroacoustic parameters thereof. 
         [0005]    Speaker parameters refer to the physical properties which affect the performance of a speaker mechanically and acoustically such as the resonance frequency, frequency response, mechanical system quality factor and electrical system quality factor. However, conventional parameter identification tools for the electroacoustic system of speakers, such as the procedure proposed by R. H. Small in “Closed-Box Loudspeaker Systems Part 1: Analysis”,  Journal of the Audio Engineering Society,  1972, can only apply to large-size speakers. When the conventional parameter identification tools are applied to a microspeaker, the results will be incorrect because the volume of the microspeaker is too small to perform a precise measurement. Currently, the technology of microspeakers focuses on the fields of the design, assemblage, and impedance measurement, and appropriate system and method for the evaluation and analysis of microspeakers have not appeared yet. 
         [0006]    Accordingly, the present invention proposes a method for parameter identification and parameter optimization of electroacoustic systems of microspeakers. 
       SUMMARY OF THE INVENTION 
       [0007]    To attain the advantages of the present method and overcome the disadvantages of the conventional method in accordance with the purpose of the invention as embodied and broadly described herein, the present invention provides a means for parameter identification of microspeakers. Measurement procedures are required to identify the electromechanical constants of the microspeaker under test. 
         [0008]    Another objective of the present invention is to provide a method for parameter optimization of microspeakers with the aid of an optimization algorithm. An optimal parameter design of a micro speaker under limiting conditions can be obtained to achieve the best acoustic performance and minimum harmonic distortion. 
         [0009]    In the present invention, an external circuit serves as the front-end to measure the impedance frequency response of the microspeaker. The front-end comprising a passive circuit and a signal analyzer is capable of measuring the impedance frequency response of the microspeaker as a dedicated impedance analyzer. 
         [0010]    To achieve the abovementioned objective, the present invention proposes a method for parameter identification of a microspeaker, wherein firstly, the impedance frequency response of a microspeaker is measured; next, the microspeaker is placed inside a test box to measure its impedance frequency response; next, a first simulation circuit is used to simulate the peak value of the impedance frequency response curve, and a second simulation circuit is used to simulate the peak value of the inside-test box impedance frequency response curve; then, the transfer functions of the first simulation circuit and the second simulation circuit are calculated to obtain the parameters of the microspeaker. 
         [0011]    In the measurement of impedance frequency response, a voltage is input into a passive circuit, which comprises the microspeaker and a load with known impedance, and then, the voltage and the obtained voltage drop over the load are input to a signal analyzer to calculate the impedance frequency response of the microspeaker. 
         [0012]    In the microspeaker parameter optimization of the present invention, parameter identification is performed for at least one micro speaker firstly; next, a target parameter and at least one limit parameter that is used as a limiting condition, are selected from parameters; then, the target parameter is optimized under the limiting condition with an optimization algorithm. 
         [0013]    These and other objectives of the present invention will become obvious to those of ordinary skill in the art after reading the following detailed description of preferred embodiments. 
         [0014]    It is understood that both the foregoing general description and the following detailed description are exemplary, and are intended to provide further explanation of the invention as claimed. 
     
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0015]    The accompanying drawings are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification. The drawings illustrate embodiments of the invention and, together with the description, serve to explain the principles of the invention. In the drawings: 
           [0016]      FIG. 1  is a diagram showing an impedance frequency response measurement device according to one embodiment of the present invention. 
           [0017]      FIG. 2  is a diagram showing a test box according to the present invention. 
           [0018]      FIG. 3  is a flowchart of the method for parameter identification and parameter optimization according to the present invention. 
           [0019]      FIG. 4  is a diagram showing the results of measuring the impedance frequency response of a microspeaker respectively disposed inside and outside a test box. 
           [0020]      FIG. 5  is a diagram showing the axial sound pressure frequency response functions before and after the parameter optimization of the axial sound pressure sensitivity of a microspeaker. 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
       [0021]    Reference will now be made in detail to the preferred embodiments of the present invention, examples of which are illustrated in the accompanying drawings. Wherever possible, the same reference numbers are used in the drawings and the description to refer to the same or like parts. 
         [0022]    The present invention utilizes a front-end device to measure the impedance frequency response of a microspeaker and utilizes a test-box method to obtain the impedance curve of the microspeaker. The electromechanical parameters of the microspeaker are calculated according to the impedance curve. After the electromechanical parameters have been identified, the performances of the microspeaker are evaluated, including: sound-pressure sensitivity, efficiency, total harmonic distortion, and inter-modulation distortion. Then, the analysis and design for optimizing the electromechanical parameters of the microspeaker are undertaken to obtain the best output performance of the microspeaker. 
         [0023]    Refer to  FIG. 1  a diagram showing an impedance frequency response measurement device according to one embodiment of the present invention. In this embodiment, a signal generator  10 , a microspeaker  12  and a load  14  with a known impedance form a measurement circuit, and a resistor is used as the load herein. The signal generator  10  outputs an alternating voltage e s , one branch of the positive pole of the signal generator  10  is connected to a signal analyzer  16 , such as a first channel ch 1  of the spectrum analyzer; via the resistor  14 , the other branch of the positive pole of the signal generator  10  is connected to the microspeaker  12 ; and the negative pole of the signal generator  10  is also connected to microspeaker  12 . When a current flows through the resistor  14 , there is a voltage drop e occurring, and thus, the voltage over the microspeaker  12  is e s −e. Therefore, once the voltage drop over the resistor  14  is obtained and input into the signal analyzer  16  via a second channel ch 2  thereof, the signal analyzer  16  can calculate the impedance frequency response Z of the microspeaker  12  according to the equation 
         [0000]    
       
         
           
             
               Z 
               = 
               
                 
                   R 
                    
                   
                     
                       
                         e 
                         s 
                       
                       - 
                       e 
                     
                     e 
                   
                 
                 = 
                 
                   R 
                    
                   
                     ( 
                     
                       
                         1 
                         
                           H 
                            
                           
                             ( 
                             f 
                             ) 
                           
                         
                       
                       - 
                       1 
                     
                     ) 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    wherein H(f) is the impedance frequency response of the resistor  14 , and R is the impedance of the resistor  14 . 
         [0024]    After the impedance frequency response of the microspeaker has been obtained, the parameters of the microspeaker can be measured. Limited by the size of the microspeaker, the parameter identification is undertaken with a test-box method in the present invention, as shown in  FIG. 2 . The test box must be an airtight chamber, and no air leakage is permitted. Refer to  FIG. 3  a flowchart of the parameter identification and parameter optimization process according to the present invention. In Step  10 , the microspeaker is respectively disposed inside and outside the test box, and the measurement device shown in  FIG. 1  is used to measure the inside-test box impedance frequency response and the outside-test box impedance frequency response respectively. The results are shown in  FIG. 4 , and the red curve represents the outside-test box impedance frequency response curve of the microspeaker, and the blue curve represents the inside-test box impedance frequency response curve of the microspeaker. 
         [0025]    Next, the process proceeds to Step  12 . A first simulation circuit, which comprises a resistor, an inductor and a capacitor, is used to simulate the peak value of the outside-test box impedance frequency response curve of the microspeaker. A second simulation circuit, which also comprises a resistor, an inductor and a capacitor, is used to simulate the peak value of the inside-test box impedance frequency response curve of the microspeaker. The objective of the abovementioned simulation is to utilize a curve-fitting method to identify the mechanical system quality factor Q MS  and the closed-box system electrical quality factor Q EC . The simulation steps comprise selecting appropriate resistance R, inductance M, and capacitance C so that the peak value of the frequency response curve of the first simulation circuit comprising said resistor, said inductor and said capacitor is the same as the peak value of the outside-test box impedance frequency response curve of the microspeaker; comparing the coefficients of the second order transfer function 
         [0000]    
       
         
           
             1 
             
               
                 j 
                  
                 
                     
                 
                  
                 ω 
                  
                 
                     
                 
                  
                 M 
               
               + 
               R 
               + 
               
                 1 
                 
                   j 
                    
                   
                       
                   
                    
                   ω 
                    
                   
                       
                   
                    
                   C 
                 
               
             
           
         
       
     
         [0000]    of the inductor, resistor and capacitor with 
         [0000]    
       
         
           
             
               1 
               
                 
                   s 
                   2 
                 
                 + 
                 
                   2 
                    
                   
                       
                   
                    
                   ζ 
                    
                   
                       
                   
                    
                   
                     ω 
                     s 
                   
                    
                   s 
                 
                 + 
                 
                   ω 
                   s 
                   2 
                 
               
             
             , 
           
         
       
     
         [0000]    the resonance frequency ω S  and the mechanical system quality factor Q MS  of the microspeaker are then identified by utilizing Equation (1) to (3), as shown in Step  14 , wherein Equations (1) to (3) are respectively expressed by: 
         [0000]    
       
         
           
             
               
                 
                   
                     ω 
                     s 
                   
                   = 
                   
                     2 
                      
                     
                         
                     
                      
                     π 
                      
                     
                         
                     
                      
                     
                       f 
                       s 
                     
                   
                 
               
               
                 
                   ( 
                   1 
                   ) 
                 
               
             
             
               
                 
                   
                     Q 
                     MS 
                   
                   = 
                   
                     1 
                     
                       2 
                        
                       
                           
                       
                        
                       ζ 
                     
                   
                 
               
               
                 
                   ( 
                   2 
                   ) 
                 
               
             
             
               
                 
                   
                     Q 
                     ES 
                   
                   = 
                   
                     
                       Q 
                       MS 
                     
                      
                     
                       ( 
                       
                         
                           R 
                           E 
                         
                         
                           R 
                           ES 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   ( 
                   3 
                   ) 
                 
               
             
           
         
       
     
         [0026]    Similarly, the inside-test box resonance frequency f c  of the microspeaker and the closed-box system electrical quality factor Q EC  are obtained via comparing the coefficients of the second order transfer function of the second simulation circuit. After the closed-box system electrical quality factor Q EC  and the mechanical system quality factor Q MS  have been identified, the equivalent volume of the microspeaker can be calculated via the equation 
         [0000]    
       
         
           
             
               
                 V 
                 AS 
               
               = 
               
                 
                   V 
                   T 
                 
                  
                 
                   ( 
                   
                     
                       
                         
                           f 
                           c 
                         
                         
                           f 
                           s 
                         
                       
                        
                       
                         
                           Q 
                           EC 
                         
                         
                           Q 
                           ES 
                         
                       
                     
                     - 
                     1 
                   
                   ) 
                 
               
             
             , 
           
         
       
     
         [0000]    wherein V T  is the volume of the test box. The mechanical mass of the vibrating diaphragm M MD , the mass of the mechanical system of the vibrating diaphragm and air load M MS  and the mechanical compliance of the vibrating diaphragm suspension C MS  can be identified with Equations (4) to (6), which are expressed by: 
         [0000]    
       
         
           
             
               
                 
                   
                     C 
                     MS 
                   
                   = 
                   
                     
                       V 
                       AS 
                     
                     
                       
                         ρ 
                         0 
                       
                        
                       
                         c 
                         2 
                       
                        
                       
                         S 
                         D 
                         2 
                       
                     
                   
                 
               
               
                 
                   ( 
                   4 
                   ) 
                 
               
             
             
               
                 
                   
                     M 
                     MS 
                   
                   = 
                   
                     1 
                     
                       
                         ω 
                         S 
                         2 
                       
                        
                       
                         C 
                         MS 
                       
                     
                   
                 
               
               
                 
                   ( 
                   5 
                   ) 
                 
               
             
             
               
                 
                   
                     M 
                     MD 
                   
                   = 
                   
                     
                       M 
                       MS 
                     
                     - 
                     
                       2 
                        
                       
                         M 
                         1 
                       
                     
                   
                 
               
               
                 
                   ( 
                   6 
                   ) 
                 
               
             
           
         
       
     
         [0000]    wherein ρ 0  is the air density; S D  is the effective area of the vibrating diaphragm; M 1  is the low-frequency air load impedance. The mechanical resistance of the vibrating diaphragm R MS  and the motor constant B 1  can be obtained with Equations (7) and (8), which are respectively expressed by: 
         [0000]    
       
         
           
             
               
                 
                   
                     R 
                     MS 
                   
                   = 
                   
                     
                       
                         ω 
                         S 
                       
                        
                       
                         M 
                         MS 
                       
                     
                     
                       Q 
                       MS 
                     
                   
                 
               
               
                 
                   ( 
                   7 
                   ) 
                 
               
             
             
               
                 
                   Bl 
                   = 
                   
                     
                       
                         
                           ω 
                           S 
                         
                          
                         
                           R 
                           E 
                         
                          
                         
                           M 
                           MS 
                         
                       
                       
                         Q 
                         ES 
                       
                     
                   
                 
               
               
                 
                   ( 
                   8 
                   ) 
                 
               
             
           
         
       
     
         [0027]    The other important parameters, such as the acoustic compliance of vibrating diaphragm suspension C AS , the acoustic mass of the vibrating diaphragm and air load M AS , the acoustic resistance of suspension loss R AS , the capacitance driving the total displacement mass C MES , the inductance driving the mechanical compliance L CES , the acoustic resistance of suspension loss and electrical loss R AT , the total mechanical resistance of suspension loss and electrical loss R MT , and the mechanical mass of the vibrating diaphragm M MD , are respectively identified by: 
         [0000]    
       
         
           
             
               
                 
                   
                     C 
                     AS 
                   
                   = 
                   
                     
                       S 
                       D 
                       2 
                     
                      
                     
                       C 
                       MS 
                     
                   
                 
               
               
                 
                   ( 
                   9 
                   ) 
                 
               
             
             
               
                 
                   
                     M 
                     AS 
                   
                   = 
                   
                     
                       M 
                       MS 
                     
                     
                       S 
                       D 
                       2 
                     
                   
                 
               
               
                 
                   ( 
                   10 
                   ) 
                 
               
             
             
               
                 
                   
                     R 
                     AS 
                   
                   = 
                   
                     
                       R 
                       MS 
                     
                     
                       S 
                       D 
                       2 
                     
                   
                 
               
               
                 
                   ( 
                   11 
                   ) 
                 
               
             
             
               
                 
                   
                     C 
                     MES 
                   
                   = 
                   
                     
                       M 
                       MS 
                     
                     
                       Bl 
                       2 
                     
                   
                 
               
               
                 
                   ( 
                   12 
                   ) 
                 
               
             
             
               
                 
                   
                     L 
                     CES 
                   
                   = 
                   
                     
                       Bl 
                       2 
                     
                      
                     
                       C 
                       MS 
                     
                   
                 
               
               
                 
                   ( 
                   13 
                   ) 
                 
               
             
             
               
                 
                   
                     R 
                     AT 
                   
                   = 
                   
                     
                       R 
                       AE 
                     
                     + 
                     
                       R 
                       AS 
                     
                   
                 
               
               
                 
                   ( 
                   14 
                   ) 
                 
               
             
             
               
                 
                   
                     R 
                     MT 
                   
                   = 
                   
                     
                       R 
                       AT 
                     
                     + 
                     
                       S 
                       D 
                       2 
                     
                   
                 
               
               
                 
                   ( 
                   15 
                   ) 
                 
               
             
             
               
                 
                   
                     M 
                     MD 
                   
                   = 
                   
                     
                       M 
                       MS 
                     
                     - 
                     
                       2 
                        
                       
                         S 
                         D 
                         2 
                       
                        
                       
                         M 
                         A 
                       
                     
                   
                 
               
               
                 
                   ( 
                   16 
                   ) 
                 
               
             
           
         
       
     
         [0000]    wherein R AE  is the acoustic resistance of electrical loss, and M A  is the acoustic mass. The equivalent coil resistance and the equivalent coil inductance of the speaker can be identified with the following equations: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             Z 
                             E 
                           
                            
                           
                             ( 
                             
                               j 
                                
                               
                                   
                               
                                
                               ω 
                             
                             ) 
                           
                         
                         ≈ 
                         
                           
                             
                               ( 
                               
                                 j 
                                  
                                 
                                     
                                 
                                  
                                 ω 
                               
                               ) 
                             
                             n 
                           
                            
                           
                             L 
                             E 
                           
                         
                       
                        
                       
                         
 
                       
                       ⇒ 
                       
                         R 
                         E 
                         ′ 
                       
                     
                     = 
                     
                       
                         [ 
                         
                           
                             L 
                             e 
                           
                           
                             cos 
                              
                             
                               ( 
                               
                                 n 
                                  
                                 
                                     
                                 
                                  
                                 
                                   π 
                                   / 
                                   2 
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                        
                       
                         ω 
                         n 
                       
                     
                   
                   , 
                   
                     
                       L 
                       E 
                     
                     = 
                     
                       
                         [ 
                         
                           
                             L 
                             e 
                           
                           
                             sin 
                              
                             
                               ( 
                               
                                 n 
                                  
                                 
                                     
                                 
                                  
                                 
                                   π 
                                   / 
                                   2 
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                        
                       
                         ω 
                         
                           n 
                           - 
                           1 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   17 
                   ) 
                 
               
             
           
         
       
     
         [0000]    The values of n and L e  can be worked out with the measured value Z VC  and the following equation: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       Z 
                       E 
                     
                     = 
                     
                       
                         Z 
                         VC 
                       
                       - 
                       
                         R 
                         E 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       n 
                       = 
                       
                         
                           
                             1 
                             90 
                           
                            
                           
                             
                               tan 
                               
                                 - 
                                 1 
                               
                             
                              
                             
                               [ 
                               
                                 
                                   Im 
                                    
                                   
                                     ( 
                                     
                                       Z 
                                       E 
                                     
                                     ) 
                                   
                                 
                                 
                                   Re 
                                    
                                   
                                     ( 
                                     
                                       Z 
                                       E 
                                     
                                     ) 
                                   
                                 
                               
                               ] 
                             
                           
                         
                         = 
                         
                           
                             
                               ln 
                                
                               
                                  
                                 
                                   Z 
                                   2 
                                 
                                  
                               
                             
                             - 
                             
                               ln 
                                
                               
                                  
                                 
                                   Z 
                                   1 
                                 
                                  
                               
                             
                           
                           
                             
                               ln 
                                
                               
                                   
                               
                                
                               
                                 ω 
                                 2 
                               
                             
                             - 
                             
                               ln 
                                
                               
                                   
                               
                                
                               
                                 ω 
                                 1 
                               
                             
                           
                         
                       
                     
                     , 
                     
                       
                         L 
                         e 
                       
                       = 
                       
                         
                            
                           
                             Z 
                             E 
                           
                            
                         
                         
                           ω 
                           n 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   18 
                   ) 
                 
               
             
           
         
       
     
         [0028]    The calculation of the abovementioned parameters can be implemented with software having calculation function, such as Matlab GUI. After the outside-test box impedance response frequency of the microspeaker, the inside-test box impedance response frequency of the microspeaker and the size of the test box have been input, Matlab can automatically calculate the values of the abovementioned parameters. Therefore, the parameter identification method of the present invention can be presented in the form of a computer program. 
         [0029]    Further, the present invention proposes an optimization method for the parameters of microspeakers. Since microspeakers are limited in volume and thickness, and the elements of a microspeaker are separately fabricated before assembled, it is hard to ensure that the elements are perfect matching, and the acoustic volume and quality of the microspeaker is hard to achieve the best performance. Thus, an optimization method is needed to fully achieve the designed performance of microspeakers. In the optimization method of the present invention, a target parameter and a limit parameter (used as a limiting condition) are selected from parameters; under the limiting condition, an optimization algorithm is used to perform optimization and find the maximum or minimum of the target parameter, as shown in Step  16 . For example, when the target parameter is the axial sound pressure sensitivity p sens   1V , it is the value of the sound pressure sensitivity at the axial distance of 1 meter and under an input voltage e g =1 V rms . The limiting condition may be the displacement of the vibrating diaphragm, the density of magnetic flux, the acoustic compliance, the resonance frequency, etc. The aim of the optimization is to obtain the maximum sound pressure sensitivity. 
         [0030]    Refer to  FIG. 5  for the comparison between the frequency response functions of axial sound pressure before and after the parameter optimization, wherein the red curve is the measurement result of the frequency response functions of axial sound pressure after the parameter optimization, and the blue curve is the simulation result of the frequency response functions of axial sound pressure without the parameter optimization, and the black curve is the measurement result of the frequency response functions of axial sound pressure without the parameter optimization. In the parameter optimization, a Sequential Quadratic Programming method is used to optimize the sensitivity of axial sound pressure; the impedance of the microspeaker, the velocity of the vibrating diaphragm and the frequency response function of axial sound pressure are respectively obtained with the following equations: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       Z 
                       VC 
                     
                      
                     
                       ( 
                       s 
                       ) 
                     
                   
                   = 
                   
                     
                       R 
                       E 
                     
                     + 
                     
                       
                         L 
                         E 
                       
                        
                       s 
                        
                       
                           
                       
                        
                       
                         •R 
                         E 
                         ′ 
                       
                     
                     + 
                     
                       
                         R 
                         ES 
                       
                        
                       
                         
                           
                             ( 
                             
                               1 
                               / 
                               
                                 Q 
                                 MS 
                               
                             
                             ) 
                           
                            
                           
                             ( 
                             
                               s 
                               / 
                               
                                 ω 
                                 s 
                               
                             
                             ) 
                           
                         
                         
                           
                             
                               ( 
                               
                                 s 
                                 / 
                                 
                                   ω 
                                   s 
                                 
                               
                               ) 
                             
                             2 
                           
                           + 
                           
                             
                               ( 
                               
                                 1 
                                 / 
                                 
                                   Q 
                                   MS 
                                 
                               
                               ) 
                             
                              
                             
                               ( 
                               
                                 s 
                                 / 
                                 
                                   ω 
                                   s 
                                 
                               
                               ) 
                             
                           
                           + 
                           1 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   19 
                   ) 
                 
               
             
             
               
                 
                   
                     U 
                     D 
                   
                   = 
                   
                     
                       
                         
                           S 
                           D 
                         
                          
                         
                           e 
                           g 
                         
                       
                       Bl 
                     
                      
                     
                       
                         R 
                         AE 
                       
                       
                         R 
                         AT 
                       
                     
                      
                     
                       
                         
                           ( 
                           
                             1 
                             / 
                             
                               Q 
                               TS 
                             
                           
                           ) 
                         
                          
                         
                           ( 
                           
                             s 
                             / 
                             
                               ω 
                               s 
                             
                           
                           ) 
                         
                       
                       
                         
                           
                             ( 
                             
                               s 
                               / 
                               
                                 ω 
                                 s 
                               
                             
                             ) 
                           
                           2 
                         
                         + 
                         
                           
                             ( 
                             
                               1 
                               / 
                               
                                 Q 
                                 TS 
                               
                             
                             ) 
                           
                            
                           
                             ( 
                             
                               s 
                               / 
                               
                                 ω 
                                 s 
                               
                             
                             ) 
                           
                         
                         + 
                         1 
                       
                     
                   
                 
               
               
                 
                   ( 
                   20 
                   ) 
                 
               
             
             
               
                 
                   
                     p 
                      
                     
                       ( 
                       
                         r 
                         = 
                         
                           1 
                            
                           m 
                         
                       
                       ) 
                     
                   
                   = 
                   
                     
                       
                         ρ 
                         0 
                       
                       
                         2 
                          
                         
                             
                         
                          
                         π 
                       
                     
                      
                     
                       
                         Ble 
                         g 
                       
                       
                         
                           S 
                           D 
                         
                          
                         
                           R 
                           E 
                         
                          
                         
                           M 
                           AS 
                         
                       
                     
                      
                     
                       
                         
                           ( 
                           
                             s 
                             / 
                             
                               ω 
                               s 
                             
                           
                           ) 
                         
                         2 
                       
                       
                         
                           
                             ( 
                             
                               s 
                               / 
                               
                                 ω 
                                 s 
                               
                             
                             ) 
                           
                           2 
                         
                         + 
                         
                           
                             ( 
                             
                               1 
                               / 
                               
                                 Q 
                                 TS 
                               
                             
                             ) 
                           
                            
                           
                             ( 
                             
                               s 
                               / 
                               
                                 ω 
                                 s 
                               
                             
                             ) 
                           
                         
                         + 
                         1 
                       
                     
                   
                 
               
               
                 
                   ( 
                   21 
                   ) 
                 
               
             
           
         
       
     
         [0000]    wherein R AT =R AE +R AS . It can be observed from  FIG. 5  that the frequency response curve after the parameter optimization is smoother, i.e. the acoustic balance and the sensitivity of axial sound pressure of the microspeaker are better. 
         [0031]    In summary, the present invention provides a method of utilizing an external electronic circuit to measure the impedance frequency response of a microspeaker. The simple external electronic circuit serves as the front-end and replaces the conventional impedance analyzer. Further, the present invention proposes a method for parameter identification of a microspeaker, wherein the parameters of a microspeaker are identified via measurement procedures for identifying electromechanical constants. After the parameters of the microspeaker have been calculated, the optimal parameter design can be obtained so that the microspeaker can achieve the best acoustic performance with minimum harmonic distortion. 
         [0032]    Those embodiments described above are to clarify the present invention to enable the person skilled in the art to understand, make and use the present invention. However, it is not intended to limit the scope of the present invention. Any equivalent modification and variation according to the spirit of the present invention is to be also included within the claims of the present invention stated below.