CN107169232B - Method and device for extracting circuit model parameters of acoustic wave filter - Google Patents
Method and device for extracting circuit model parameters of acoustic wave filter Download PDFInfo
- Publication number
- CN107169232B CN107169232B CN201710434138.6A CN201710434138A CN107169232B CN 107169232 B CN107169232 B CN 107169232B CN 201710434138 A CN201710434138 A CN 201710434138A CN 107169232 B CN107169232 B CN 107169232B
- Authority
- CN
- China
- Prior art keywords
- value
- harmonic oscillator
- circuit model
- acoustic wave
- objective function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Oscillators With Electromechanical Resonators (AREA)
Abstract
The invention relates to a method and a device for extracting circuit model parameters of an acoustic wave filter. The method comprises the following steps: performing multi-physical field simulation on a physical model of the acoustic wave filter, and extracting a first admittance value of each harmonic oscillator from a multi-physical field simulation result; respectively optimizing the initial values of the set circuit model parameters of each harmonic oscillator, and calculating a second admittance value of each harmonic oscillator according to the optimized parameter values of the circuit model parameters of each harmonic oscillator; respectively inputting the first admittance value and the second admittance value of each harmonic oscillator into the objective function to obtain the objective function value of each harmonic oscillator; if the objective function value of each harmonic oscillator meets a preset convergence condition, extracting the parameter value of the circuit model parameter of each harmonic oscillator corresponding to the objective function value, and taking the extracted parameter values of the circuit model parameters of all harmonic oscillators as the parameter values of the circuit model parameters of the acoustic wave filter.
Description
Technical Field
The present invention relates to the field of filter technologies, and in particular, to a method for extracting parameters of an acoustic wave filter circuit model, an apparatus for extracting parameters of an acoustic wave filter circuit model, a computer-readable storage medium, and a computer device.
Background
The acoustic wave filter (such as a bulk acoustic wave filter and a surface acoustic wave filter) based on a Micro-Electro-Mechanical System (MEMS) process is small in size, excellent in indexes such as insertion loss and power capacity and has wide prospects in miniaturized terminal equipment and System equipment. Especially, mobile communication is continuously updated, system devices and terminals are multi-frequency and multi-mode products, and the acoustic wave filter module with miniaturization and excellent performance is indispensable in order to eliminate interference between systems and ensure communication quality.
Simulation is an important link in the design of an acoustic wave filter, and is generally divided into circuit simulation and field simulation, wherein the former is to perform optimization analysis on an original equivalent circuit model of the acoustic wave filter to enable the characteristics to meet index requirements, so that a basis is provided for modeling of a physical model of the latter, namely, a real physical structure (physical model) of the acoustic wave filter is established according to the circuit model, and the characteristics of the acoustic wave filter are obtained through simulation calculation of multiple physical fields (an electric field, a magnetic field, an acoustic wave field and the like).
In the design process of the acoustic wave filter, it is important to extract corresponding circuit model parameters from the multi-physical field simulation result, because the corresponding circuit model parameters can be compared with corresponding element values in the original equivalent circuit model, differences between the circuit model parameters and the original equivalent circuit model can be found out, the corresponding structures of the acoustic wave filter can be further modified in sequence until the circuit model parameters and the element values are close to each other, and therefore the final structure size of the acoustic wave filter can be determined.
Simulation and test show that harmonic oscillators contained in the acoustic wave filter have secondary modes near a main mode, and have interference effect on the extraction of circuit model parameters, so that the difficulty of accurately proposing a parameter algorithm is increased, and an effective acoustic wave filter circuit model parameter extraction method is not provided in the conventional technology so far.
Disclosure of Invention
Based on the method, aiming at the problem that the traditional technology lacks a method for effectively extracting the parameters of the circuit model of the acoustic wave filter, the method and the device for extracting the parameters of the circuit model of the acoustic wave filter are provided, so that the parameters of the circuit model of the acoustic wave filter can be effectively extracted.
A method for extracting parameters of a circuit model of an acoustic wave filter comprises the following steps:
a method for extracting parameters of a circuit model of an acoustic wave filter is characterized by comprising the following steps:
performing multi-physical field simulation on a physical model of an acoustic wave filter, and extracting a first admittance value of each harmonic oscillator from a multi-physical field simulation result, wherein the acoustic wave filter comprises a plurality of harmonic oscillators;
respectively optimizing the initial values of the set circuit model parameters of each harmonic oscillator, and calculating a second admittance value of each harmonic oscillator according to the optimized parameter values of the circuit model parameters of each harmonic oscillator;
respectively inputting the first admittance value and the second admittance value of each harmonic oscillator into a preset objective function to obtain an objective function value of each harmonic oscillator, wherein the objective function value is a difference value between the second admittance value and the first admittance value;
and if the objective function value of each harmonic oscillator meets a preset convergence condition, extracting the parameter value of the circuit model parameter of each harmonic oscillator corresponding to the objective function value, and taking the extracted parameter values of the circuit model parameters of all harmonic oscillators as the parameter values of the circuit model parameters of the acoustic wave filter.
According to the method for extracting the circuit model parameters of the acoustic wave filter, the objective function value is calculated according to the multi-physical-field simulation result of the physical model of the acoustic wave filter and the theoretical calculation result of the admittance value, the circuit model parameters of the acoustic wave filter are effectively extracted according to the convergence condition that the objective function value meets the convergence condition, the frequency characteristics of each harmonic oscillator in the acoustic wave filter are accurately fitted through the circuit model parameters extracted by the method, and a basis is provided for the structure optimization of the acoustic wave filter.
In one embodiment, the method for extracting the parameters of the acoustic wave filter circuit model further comprises the following steps: and if the objective function value of the harmonic oscillator does not meet the preset convergence condition, re-optimizing the parameter value of the circuit model parameter of the harmonic oscillator, and re-obtaining the objective function value of the harmonic oscillator according to the optimized parameter value.
In one embodiment, the objective function is:
wherein, f (C)i0,Cim,Lim,Rim) Is the value of the objective function, Ci0Is the static capacitance of the ith harmonic oscillator, CimIs the ith harmonicEquivalent capacitance of the vibrator, LimIs the equivalent inductance of the ith harmonic oscillator, RimIs the resistance of the ith harmonic oscillator, k is the kth frequency point, n is the total number of frequency points,
is the phase value of the second admittance value of the ith harmonic oscillator,
is the phase value, omega, of the first admittance value of the ith harmonic oscillatorkIs the angular frequency of the k-th frequency point, Ai(ωk) Is the amplitude, A ', of the second admittance value of the ith harmonic oscillator'i(ωk) Is the amplitude of the first admittance value of the ith harmonic oscillator.
According to the set target function, an ideal result can be obtained, and the accuracy of circuit model parameter extraction is improved.
In one embodiment, before performing the multi-physical field simulation on the physical model of the acoustic wave filter, the method further comprises the following steps: obtaining an equivalent circuit model of the acoustic wave filter, and constructing a physical model of the acoustic wave filter according to the equivalent circuit model, wherein the model of each harmonic oscillator in the equivalent circuit model is a BVD circuit model.
In one embodiment, the circuit model parameters include static capacitance, equivalent inductance, and resistance; before optimizing the initial values of the circuit model parameters of each set harmonic oscillator, the method further comprises the following steps: and correspondingly setting the initial value of the static capacitor, the initial value of the equivalent inductor and the initial value of the resistor of each harmonic oscillator according to the value of the static capacitor, the value of the equivalent inductor and the value of the resistor in the BVD circuit model of each harmonic oscillator. According to the method, the initial value is set, an ideal result can be obtained, and the accuracy of circuit model parameter extraction is improved.
In one embodiment, extracting the first admittance value of each resonator from the result of the multi-physics simulation comprises the steps of: and extracting first admittance values of each harmonic oscillator at different frequency points in the pass band from the result of the multi-physical field simulation.
In one embodiment, the initial values and/or parameter values of the circuit model parameters are optimized by a numerical optimization algorithm. According to the method, the circuit model parameters are searched, and the convergence speed of searching calculation is improved.
An extraction device for circuit model parameters of an acoustic wave filter, comprising:
the first admittance value obtaining module is used for carrying out multi-physical field simulation on a physical model of the acoustic wave filter, and extracting a first admittance value of each harmonic oscillator from a multi-physical field simulation result, wherein the acoustic wave filter comprises a plurality of harmonic oscillators;
the second admittance value obtaining module is used for respectively optimizing the initial values of the set circuit model parameters of each harmonic oscillator and calculating the second admittance values of each harmonic oscillator according to the optimized parameter values of the circuit model parameters of each harmonic oscillator;
an objective function value obtaining module, configured to input the first admittance value and the second admittance value of each resonator into a preset objective function, respectively, to obtain an objective function value of each resonator, where the objective function value is a difference between the second admittance value and the first admittance value;
and the circuit model parameter extraction module is used for extracting the parameter values of the circuit model parameters of each harmonic oscillator corresponding to the objective function values when the objective function values of the harmonic oscillators meet the preset convergence condition, and taking the extracted parameter values of the circuit model parameters of all the harmonic oscillators as the parameter values of the circuit model parameters of the acoustic wave filter.
The extraction device for the acoustic wave filter circuit model parameters calculates the objective function value according to the multi-physical field simulation result of the acoustic wave filter physical model and the theoretical calculation result of the admittance value, effectively extracts the circuit model parameters of the acoustic wave filter according to the convergence condition that the objective function value meets the convergence condition, more accurately fits the frequency characteristic of each harmonic oscillator in the acoustic wave filter through the circuit model parameters extracted by the extraction device, and provides a basis for the structure optimization of the acoustic wave filter.
A computer-readable storage medium, on which a computer program is stored, which program, when being executed by a processor, is adapted to carry out the steps of the method of any of the preceding claims. The computer readable storage medium calculates an objective function value according to a multi-physical field simulation result of a physical model of the acoustic wave filter and a theoretical calculation result of an admittance value, effectively extracts a circuit model parameter of the acoustic wave filter according to the objective function value meeting a convergence condition, more accurately fits the frequency characteristic of each harmonic oscillator in the acoustic wave filter through the circuit model parameter extracted in the mode, and provides a basis for the structure optimization of the acoustic wave filter.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of any of the methods described above when executing the program. The computer equipment calculates the objective function value according to the multi-physical field simulation result of the physical model of the acoustic wave filter and the theoretical calculation result of the admittance value, effectively extracts the circuit model parameters of the acoustic wave filter according to the convergence condition that the objective function value meets the convergence condition, more accurately fits the frequency characteristic of each harmonic oscillator in the acoustic wave filter through the circuit model parameters extracted in the mode, and provides a basis for the structure optimization of the acoustic wave filter.
Drawings
FIG. 1 is a schematic flow chart of a method for extracting parameters of a circuit model of an acoustic wave filter according to an embodiment;
FIG. 2 is a schematic diagram of an equivalent circuit of a bulk acoustic wave harmonic oscillator BVD;
FIG. 3 is a schematic diagram of a four-unit ladder structure bulk acoustic wave filter;
FIG. 4 is a frequency response curve calculated by simulation of multiple physical fields of a four-unit ladder structure bulk acoustic wave filter and an extracted circuit model;
FIG. 5 is a schematic structural diagram of an apparatus for extracting parameters of a circuit model of an acoustic wave filter according to an embodiment;
fig. 6 is a schematic structural diagram of a computer device according to an embodiment.
Detailed Description
In order to further explain the technical means and effects of the present invention, the following description of the present invention with reference to the accompanying drawings and preferred embodiments will be made for clarity and completeness.
As shown in fig. 1, in one embodiment, there is provided a method for extracting parameters of a circuit model of an acoustic wave filter, including the steps of:
s110, performing multi-physical field simulation on a physical model of the acoustic wave filter, and extracting a first admittance value of each harmonic oscillator from a multi-physical field simulation result, wherein the acoustic wave filter comprises a plurality of harmonic oscillators;
s120, respectively optimizing the initial values of the set circuit model parameters of the harmonic oscillators, and calculating second admittance values of the harmonic oscillators according to the optimized parameter values of the circuit model parameters of the harmonic oscillators;
s130, respectively inputting the first admittance value and the second admittance value of each harmonic oscillator into a preset objective function to obtain an objective function value of each harmonic oscillator, wherein the objective function value is a difference value between the second admittance value and the first admittance value;
and S140, if the objective function value of each harmonic oscillator meets a preset convergence condition, extracting the parameter value of the circuit model parameter of each harmonic oscillator corresponding to the objective function value, and taking the extracted parameter values of the circuit model parameters of all harmonic oscillators as the parameter values of the circuit model parameters of the acoustic wave filter.
The Acoustic Wave filter includes a Bulk Acoustic Wave (BAW) filter, a Surface Acoustic Wave (SAW) filter, and the like. The physical model of the acoustic wave filter reflects the true physical structure of the acoustic wave filter. In one embodiment, before performing the multi-physical field simulation on the physical model of the acoustic wave filter, the method may further include the steps of: obtaining an equivalent circuit model of the acoustic wave filter, and constructing a physical model of the acoustic wave filter according to the equivalent circuit model, wherein the model of each harmonic oscillator in the equivalent circuit model is a BVD (Butterworth-Van Dyke, Butterworth-Vandyk) circuit model.
The equivalent circuit model is formed by abstracting the actual circuitApproximately reflecting the electrical characteristics of the actual circuit. Since the acoustic wave filter includes a plurality of resonators, the basic unit of the equivalent circuit model of the acoustic wave filter is the BVD circuit model of each resonator. Taking the BVD circuit model of the bulk acoustic wave harmonic oscillator as an example, as shown in FIG. 2, it includes a static capacitance C formed by electrodes0Equivalent capacitance C associated with mechanical vibrationsmEquivalent inductance L associated with mechanical vibrationsmAnd a resistance R associated with the lossesm. The topological structure of the four-unit bulk acoustic wave filter is shown in fig. 3, wherein a harmonic oscillator 1 and a harmonic oscillator 4 are connected in series, and a harmonic oscillator 2 and a harmonic oscillator 3 are connected in parallel.
A real physical structure (physical model) of the acoustic wave filter is established according to an equivalent circuit model of the acoustic wave filter, and then the characteristics of the acoustic wave filter are obtained through simulation calculation of multiple physical fields (an electric field, a magnetic field, an acoustic wave field and the like), wherein the establishment of the physical model and the simulation of the multiple physical fields can be realized by adopting the existing mode in the prior art. The admittance values of each harmonic oscillator at each frequency point can be extracted from the result of the multi-physical-field simulation, for example, the admittance corresponding to the ith harmonic oscillator is y'i(ω), where ω represents angular frequency, indicating admittance yi' is a function of frequency.
In the equivalent circuit model of BVD, each harmonic oscillator has two resonant frequencies: series resonant frequency fsAnd parallel resonant frequency fp. In the acoustic wave filter structure, the series resonance frequency of the harmonic oscillator (such as harmonic oscillators 1 and 4 in fig. 3) in the series path is positioned in the pass band of the acoustic wave filter, and the parallel resonance frequency is positioned outside the pass band; the two resonant frequencies of the parallel resonators (such as resonators 2 and 3 in fig. 3) are exactly opposite to each other, i.e. the series frequency is outside the pass band, and the parallel frequency is inside the pass band. Therefore, in one embodiment, extracting the first admittance value of each resonator from the result of the multi-physical field simulation comprises the steps of: and extracting first admittance values of each harmonic oscillator at different frequency points in the pass band from the result of the multi-physical field simulation. When multi-physical-field simulation is carried out, the admittance value is calculated only by taking points in a section of continuous spectrum containing the pass band, namely, the multi-physical-field simulation result is extractedDifferent frequency points omega of each harmonic oscillator in the pass band1,…ωk,…,ωnAnd admittance value y 'of'i(ωk)。
In one embodiment, the circuit model parameters include static capacitance Ci0And an equivalent capacitance CimEquivalent inductance LimAnd a resistance RimAnd the like. Presetting initial values of circuit model parameters of each harmonic oscillator, optimizing each initial value to obtain optimized parameter values, and calculating the optimized parameter values by adopting an admittance value calculation method to obtain theoretically calculated admittance values.
Taking the BVD circuit model as an example, theoretically, the admittance of the ith harmonic oscillator is:
wherein, yi(ω) is the admittance of the ith harmonic oscillator; omega is angular frequency;
is an imaginary unit; ci0The static capacitance is the static capacitance of the ith harmonic oscillator; l isimThe equivalent inductance of the ith harmonic oscillator; cimThe equivalent capacitance of the ith harmonic oscillator; rimIs the resistance of the ith harmonic oscillator.
Admittance yiThe amplitude and phase of (ω) are:
amplitude:
phase position:
extracting circuit model parameters Ci0、Cim、LimAnd RimThe mathematical problem of (2) is to solve a multivariate nonlinear function yi(ω) its value at each frequency point and according to the multiple physics simulation result yi′(ωk) BetweenThe difference value satisfies a preset convergence condition, which may be that the obtained objective function value is a minimum value, for example, by setting a minimum threshold, and if the objective function value is smaller than the minimum threshold, the minimum value of the objective function is found. Taking the objective function value as a minimum value as an example, the process of solving the minimum value of the objective function value is as follows: searching the final value of the circuit model parameter when the target value is minimum, namely, giving C according to the initial value of the set circuit model parameteri0、Cim、LimAnd RimFinding C from reasonable initial valuei0、Cim、LimAnd RimTo minimize the value of the objective function.
If the objective function value of a harmonic oscillator meets the convergence condition, the parameter value of the circuit model parameter corresponding to the theoretically calculated admittance value (second admittance value) of the harmonic oscillator can be extracted according to the objective function value meeting the convergence condition, the extracted parameter value of the circuit model parameter is the circuit model parameter value of the harmonic oscillator extracted through the multi-physical field simulation result, and each harmonic oscillator extracts the circuit model parameter according to the method, so that the circuit model parameters of all harmonic oscillators in the acoustic wave filter can be obtained.
In one embodiment, the method for extracting the parameters of the acoustic wave filter circuit model further comprises the following steps: and if the objective function value of the harmonic oscillator does not meet the preset convergence condition, re-optimizing the parameter value of the circuit model parameter of the harmonic oscillator, and re-obtaining the objective function value of the harmonic oscillator according to the optimized parameter value. If the objective function value of one harmonic oscillator does not meet the convergence condition, continuing to search the parameter values, then recalculating the second admittance value according to the searched parameter values, recalculating the objective function value according to the recalculated second admittance value and the first admittance value of the harmonic oscillator, and sequentially circulating until the objective function value of the harmonic oscillator meets the convergence condition, wherein the parameter value corresponding to the objective function value meeting the convergence condition is the circuit model parameter value of the harmonic oscillator.
Whether a desired result is achieved and the speed of convergence of the search calculation is increased is related to three factors,first, parameters (C) of the circuit model are giveni0、Cim、LimAnd Rim) An initial value of (d); second, optimizing the initial value and the parameter value; and thirdly, an objective function form. The invention provides corresponding preferable schemes aiming at three factors respectively, and the detailed description is provided below.
With respect to circuit model parameters (C)i0、Cim、LimAnd Rim) Initial value of (d): when designing the BVD circuit model of each resonator, each element is assigned with a corresponding element value, for example, a corresponding element value is set for a circuit model parameter such as a static capacitance, an equivalent inductance, and a resistance of each resonator, so in order to obtain an ideal result and increase the convergence speed of the search calculation, in an embodiment, before optimizing the initial values of the set circuit model parameters of each resonator, the method further includes the steps of: and correspondingly setting the initial value of the static capacitor, the initial value of the equivalent inductor and the initial value of the resistor of each harmonic oscillator according to the value of the static capacitor, the value of the equivalent inductor and the value of the resistor in the BVD circuit model of each harmonic oscillator. Namely, the corresponding element value in the BVD circuit model is taken as Ci0、Cim、LimAnd RimThe initial value of (c).
The optimization method for the initial values and the parameter values comprises the following steps: in order to obtain the desired result and increase the convergence speed of the search calculation, in one embodiment, the initial values and/or parameter values of the circuit model parameters may be optimized by an existing numerical optimization algorithm, for example, the numerical optimization algorithm is a conjugate gradient optimization algorithm.
For the objective function: due to admittance value yi(ωk) Is a complex number with an amplitude Ai(ωk) And phase
Two parameters, namely, the difference between the theoretically calculated admittance value and the multiphysics simulation admittance value cannot be simply used as a target function, and a target function which takes a value in a real number domain needs to be established as the basis for judging the optimization algorithm, and the real number is realThe domain values are real numbers for all values of the objective function.
In one embodiment, the objective function is:
wherein, f (C)i0,Cim,Lim,Rim) Is the value of the objective function, Ci0Is the static capacitance of the ith harmonic oscillator, CimIs the equivalent capacitance of the ith harmonic oscillator, LimIs the equivalent inductance of the ith harmonic oscillator, RimIs the resistance of the ith harmonic oscillator, k is the kth frequency point, n is the total number of frequency points,
is the phase value of the second admittance value of the ith harmonic oscillator,
is the phase value, omega, of the first admittance value of the ith harmonic oscillatorkIs the angular frequency of the k-th frequency point, Ai(ωk) Is the amplitude, A ', of the second admittance value of the ith harmonic oscillator'i(ωk) Tan () is a tangent trigonometric function and lg () is a logarithmic function with a base 10, which is the magnitude of the first admittance value of the ith harmonic oscillator.
Equation (4) represents that the objective function is obtained by taking the square after the tangent values of the admittance phase angles at the corresponding frequency points are subtracted, taking the square after the logarithm values with the amplitude of the base 10 are subtracted, and summing the two values about n frequency points after the two values are added. It should be noted that the present invention is not limited to the objective function represented by the formula (4), and equivalent modifications on the formula (4), such as modifying the logarithmic base number or increasing the coefficient, are within the protection scope of the present invention.
In addition, in the pair Ci0、Cim、LimAnd RimIn the value searching process, the range is limited, which is beneficial to shortening the calculation time and accelerating the convergence speed.
As shown in fig. 4, the frequency response of the acoustic wave filter obtained by direct multi-physics simulation of the four-element bulk acoustic wave filter shown in fig. 3 and the frequency response curve obtained by extracting the circuit parameters and using the corresponding circuit model are shown, and as can be seen from fig. 4, the two are highly matched.
Based on the same inventive concept, the invention also provides a device for extracting the circuit model parameters of the acoustic wave filter, and the following describes the specific implementation mode of the device in detail with reference to the attached drawings.
As shown in fig. 5, an apparatus for extracting parameters of a circuit model of an acoustic wave filter includes:
a first admittance value obtaining module 110, configured to perform multi-physical field simulation on a physical model of an acoustic wave filter, and extract a first admittance value of each resonator from a result of the multi-physical field simulation, where the acoustic wave filter includes a plurality of resonators;
a second admittance value obtaining module 120, configured to optimize the initial values of the circuit model parameters of the resonators, and calculate second admittance values of the resonators according to the optimized parameter values of the circuit model parameters of the resonators;
an objective function value obtaining module 130, configured to input the first admittance value and the second admittance value of each resonator into a preset objective function respectively, so as to obtain an objective function value of each resonator, where the objective function value is a difference value between the second admittance value and the first admittance value;
and the circuit model parameter extraction module 140 is configured to, when the objective function value of each resonator meets a preset convergence condition, extract a parameter value of a circuit model parameter of each resonator corresponding to the objective function value, and use the extracted parameter value of the circuit model parameter of all the resonators as the parameter value of the circuit model parameter of the acoustic wave filter.
In an embodiment, if the objective function value of the resonator does not satisfy the preset convergence condition, the second admittance value obtaining module 120 re-optimizes the parameter value of the circuit model parameter of the resonator, and the objective function value obtaining module 130 re-obtains the objective function value of the resonator according to the optimized parameter value.
In one embodiment, the objective function is:
wherein, f (C)i0,Cim,Lim,Rim) Is the value of the objective function, Ci0Is the static capacitance of the ith harmonic oscillator, CimIs the equivalent capacitance of the ith harmonic oscillator, LimIs the equivalent inductance of the ith harmonic oscillator, RimIs the resistance of the ith harmonic oscillator, k is the kth frequency point, n is the total number of frequency points,
is the phase value of the second admittance value of the ith harmonic oscillator,
is the phase value, omega, of the first admittance value of the ith harmonic oscillatorkIs the angular frequency of the k-th frequency point, Ai(ωk) Is the amplitude, A ', of the second admittance value of the ith harmonic oscillator'i(ωk) Is the amplitude of the first admittance value of the ith harmonic oscillator.
In an embodiment, before performing the multi-physical field simulation on the physical model of the acoustic wave filter, the first admittance value obtaining module 110 is further configured to obtain an equivalent circuit model of the acoustic wave filter, and construct the physical model of the acoustic wave filter according to the equivalent circuit model, where a model of each harmonic oscillator in the equivalent circuit model is a BVD circuit model.
In one embodiment, the circuit model parameters include static capacitance, equivalent inductance, and resistance; before the second admittance value obtaining module 120 optimizes the initial values of the set circuit model parameters of each resonator, the second admittance value obtaining module is further configured to set the initial values of the static capacitance, the equivalent inductance, and the resistance of each resonator according to the values of the static capacitance, the equivalent inductance, and the resistance in the BVD circuit model of each resonator.
In one embodiment, the first admittance value obtaining module 110 extracts the first admittance values of the respective harmonic oscillators at different frequency points within the pass band from the result of the multi-physical field simulation.
In one embodiment, the second admittance value obtaining module 120 optimizes the initial values and/or parameter values of the circuit model parameters by a numerical optimization algorithm.
In one embodiment, a computer-readable storage medium is provided, having stored thereon a computer program which, when executed by a processor, performs the steps of:
performing multi-physical field simulation on a physical model of an acoustic wave filter, and extracting a first admittance value of each harmonic oscillator from a multi-physical field simulation result, wherein the acoustic wave filter comprises a plurality of harmonic oscillators;
respectively optimizing the initial values of the set circuit model parameters of each harmonic oscillator, and calculating a second admittance value of each harmonic oscillator according to the optimized parameter values of the circuit model parameters of each harmonic oscillator;
respectively inputting the first admittance value and the second admittance value of each harmonic oscillator into a preset objective function to obtain an objective function value of each harmonic oscillator, wherein the objective function value is a difference value between the second admittance value and the first admittance value;
and if the objective function value of each harmonic oscillator meets a preset convergence condition, extracting the parameter value of the circuit model parameter of each harmonic oscillator corresponding to the objective function value, and taking the extracted parameter values of the circuit model parameters of all harmonic oscillators as the parameter values of the circuit model parameters of the acoustic wave filter.
In one embodiment, the computer program when executed by the processor further performs the steps of: and if the objective function value of the harmonic oscillator does not meet the preset convergence condition, re-optimizing the parameter value of the circuit model parameter of the harmonic oscillator, and re-obtaining the objective function value of the harmonic oscillator according to the optimized parameter value.
In one embodiment, the objective function is:
wherein, f (C)i0,Cim,Lim,Rim) Is the value of the objective function, Ci0Is the static capacitance of the ith harmonic oscillator, CimIs the equivalent capacitance of the ith harmonic oscillator, LimIs the equivalent inductance of the ith harmonic oscillator, RimIs the resistance of the ith harmonic oscillator, k is the kth frequency point, n is the total number of frequency points,
is the phase value of the second admittance value of the ith harmonic oscillator,
is the phase value, omega, of the first admittance value of the ith harmonic oscillatorkIs the angular frequency of the k-th frequency point, Ai(ωk) Is the amplitude, A ', of the second admittance value of the ith harmonic oscillator'i(ωk) Is the amplitude of the first admittance value of the ith harmonic oscillator.
In one embodiment, the computer program when executed by the processor further performs the steps of: before the multi-physical field simulation is carried out on the physical model of the acoustic wave filter, the method also comprises the following steps: obtaining an equivalent circuit model of the acoustic wave filter, and constructing a physical model of the acoustic wave filter according to the equivalent circuit model, wherein the model of each harmonic oscillator in the equivalent circuit model is a BVD circuit model.
In one embodiment, the circuit model parameters include static capacitance, equivalent inductance, and resistance; the computer program when executed by the processor further implements the steps of: before optimizing the initial values of the circuit model parameters of each set harmonic oscillator, the method further comprises the following steps: and correspondingly setting the initial value of the static capacitor, the initial value of the equivalent inductor and the initial value of the resistor of each harmonic oscillator according to the value of the static capacitor, the value of the equivalent inductor and the value of the resistor in the BVD circuit model of each harmonic oscillator.
In one embodiment, the computer program when executed by the processor further performs the steps of: the method for extracting the first admittance value of each harmonic oscillator from the result of the multi-physical-field simulation comprises the following steps: and extracting first admittance values of each harmonic oscillator at different frequency points in the pass band from the result of the multi-physical field simulation.
In one embodiment, the computer program when executed by the processor further performs the steps of: and optimizing the initial values and/or parameter values of the circuit model parameters through a numerical optimization algorithm.
Other technical features of the computer-readable storage medium are the same as those of the method for extracting the circuit model parameters of the acoustic wave filter, and are not described herein again.
As shown in fig. 6, in one embodiment, there is provided a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the following steps when executing the program:
performing multi-physical field simulation on a physical model of an acoustic wave filter, and extracting a first admittance value of each harmonic oscillator from a multi-physical field simulation result, wherein the acoustic wave filter comprises a plurality of harmonic oscillators;
respectively optimizing the initial values of the set circuit model parameters of each harmonic oscillator, and calculating a second admittance value of each harmonic oscillator according to the optimized parameter values of the circuit model parameters of each harmonic oscillator;
respectively inputting the first admittance value and the second admittance value of each harmonic oscillator into a preset objective function to obtain an objective function value of each harmonic oscillator, wherein the objective function value is a difference value between the second admittance value and the first admittance value;
and if the objective function value of each harmonic oscillator meets a preset convergence condition, extracting the parameter value of the circuit model parameter of each harmonic oscillator corresponding to the objective function value, and taking the extracted parameter values of the circuit model parameters of all harmonic oscillators as the parameter values of the circuit model parameters of the acoustic wave filter.
In one embodiment, the processor when executing the program further performs the steps of: and if the objective function value of the harmonic oscillator does not meet the preset convergence condition, re-optimizing the parameter value of the circuit model parameter of the harmonic oscillator, and re-obtaining the objective function value of the harmonic oscillator according to the optimized parameter value.
In one embodiment, the objective function is:
wherein, f (C)i0,Cim,Lim,Rim) Is the value of the objective function, Ci0Is the static capacitance of the ith harmonic oscillator, CimIs the equivalent capacitance of the ith harmonic oscillator, LimIs the equivalent inductance of the ith harmonic oscillator, RimIs the resistance of the ith harmonic oscillator, k is the kth frequency point, n is the total number of frequency points,
is the phase value of the second admittance value of the ith harmonic oscillator,
is the phase value, omega, of the first admittance value of the ith harmonic oscillatorkIs the angular frequency of the k-th frequency point, Ai(ωk) Is the amplitude, A ', of the second admittance value of the ith harmonic oscillator'i(ωk) Is the amplitude of the first admittance value of the ith harmonic oscillator.
In one embodiment, the processor when executing the program further performs the steps of: before the multi-physical field simulation is carried out on the physical model of the acoustic wave filter, the method also comprises the following steps: obtaining an equivalent circuit model of the acoustic wave filter, and constructing a physical model of the acoustic wave filter according to the equivalent circuit model, wherein the model of each harmonic oscillator in the equivalent circuit model is a BVD circuit model.
In one embodiment, the circuit model parameters include static capacitance, equivalent inductance, and resistance; the processor, when executing the program, further implements the steps of: before optimizing the initial values of the circuit model parameters of each set harmonic oscillator, the method further comprises the following steps: and correspondingly setting the initial value of the static capacitor, the initial value of the equivalent inductor and the initial value of the resistor of each harmonic oscillator according to the value of the static capacitor, the value of the equivalent inductor and the value of the resistor in the BVD circuit model of each harmonic oscillator.
In one embodiment, the processor when executing the program further performs the steps of: the method for extracting the first admittance value of each harmonic oscillator from the result of the multi-physical-field simulation comprises the following steps: and extracting first admittance values of each harmonic oscillator at different frequency points in the pass band from the result of the multi-physical field simulation.
In one embodiment, the processor when executing the program further performs the steps of: and optimizing the initial values and/or parameter values of the circuit model parameters through a numerical optimization algorithm.
Other technical features of the computer device are the same as those of the method for extracting the circuit model parameters of the acoustic wave filter, and are not described herein again.
The extraction of the parameters of the circuit model of the acoustic wave filter is of great importance in the design of the acoustic wave filter, and the parameters of the acoustic wave filter can be judged to be adjusted and the amplitude can be adjusted by obtaining the result, so that the design difficulty of the acoustic wave filter is reduced. According to the invention, the curve fitting and numerical optimization are carried out on the multi-physical-field simulation result to obtain the minimum value, the circuit model parameters such as the capacitance, the equivalent inductance and the resistance corresponding to each harmonic oscillator of the acoustic wave filter are efficiently and accurately extracted, the corresponding acoustic wave filter structure is analyzed and optimized, and the product simulation design time is shortened. In practical application, the frequency response of the circuit model of the acoustic wave filter obtained by the extracted circuit model parameters is highly consistent with the simulation result of the multiple physical fields, which shows that the method is accurate and effective. In addition, the invention is very suitable for programming by using the existing commercial mathematical software, and the program is simple and efficient.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
The technical features of the embodiments described above may be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the embodiments described above are not described, but should be considered as being within the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (9)
1. A method for extracting parameters of a circuit model of an acoustic wave filter is characterized by comprising the following steps:
performing multi-physical field simulation on a physical model of an acoustic wave filter, and extracting a first admittance value of each harmonic oscillator from a multi-physical field simulation result, wherein the acoustic wave filter comprises a plurality of harmonic oscillators;
respectively optimizing the initial values of the set circuit model parameters of each harmonic oscillator, and calculating a second admittance value of each harmonic oscillator according to the optimized parameter values of the circuit model parameters of each harmonic oscillator;
respectively inputting the first admittance value and the second admittance value of each harmonic oscillator into a preset objective function to obtain an objective function value of each harmonic oscillator, wherein the objective function value is a difference value between the second admittance value and the first admittance value;
if the objective function value of each harmonic oscillator meets a preset convergence condition, extracting the parameter value of the circuit model parameter of each harmonic oscillator corresponding to the objective function value, and taking the extracted parameter values of the circuit model parameters of all harmonic oscillators as the parameter values of the circuit model parameters of the acoustic wave filter;
the preset objective function is as follows:
wherein, f (C)i0,Cim,Lim,Rim) Is the value of the objective function, Ci0Is the static capacitance of the ith harmonic oscillator, CimIs the equivalent capacitance of the ith harmonic oscillator, LimIs the equivalent inductance of the ith harmonic oscillator, RimIs the resistance of the ith harmonic oscillator, k is the kth frequency point, n is the total number of frequency points,
is the phase value of the second admittance value of the ith harmonic oscillator,
is the phase value, omega, of the first admittance value of the ith harmonic oscillatorkIs the angular frequency of the k-th frequency point, Ai(ωk) Is the amplitude of the second admittance value, A, of the ith harmonic oscillatori′(ωk) The amplitude of the first admittance value of the ith harmonic oscillator;
the objective function can achieve equivalent deformation of the objective function by modifying the logarithmic base or increasing the coefficient.
2. The method for extracting parameters of an acoustic wave filter circuit model according to claim 1, further comprising the steps of:
and if the objective function value of the harmonic oscillator does not meet the preset convergence condition, re-optimizing the parameter value of the circuit model parameter of the harmonic oscillator, and re-obtaining the objective function value of the harmonic oscillator according to the optimized parameter value.
3. The method for extracting parameters of an acoustic wave filter circuit model according to any one of claims 1 to 2, wherein before performing the multi-physical field simulation on the physical model of the acoustic wave filter, the method further comprises the steps of:
obtaining an equivalent circuit model of the acoustic wave filter, and constructing a physical model of the acoustic wave filter according to the equivalent circuit model, wherein the model of each harmonic oscillator in the equivalent circuit model is a BVD circuit model.
4. The method for extracting parameters of an acoustic wave filter circuit model according to claim 3, wherein the parameters of the circuit model include static capacitance, equivalent inductance, and resistance;
before optimizing the initial values of the circuit model parameters of each set harmonic oscillator, the method further comprises the following steps:
and correspondingly setting the initial value of the static capacitor, the initial value of the equivalent inductor and the initial value of the resistor of each harmonic oscillator according to the value of the static capacitor, the value of the equivalent inductor and the value of the resistor in the BVD circuit model of each harmonic oscillator.
5. The method for extracting parameters of an acoustic wave filter circuit model according to any one of claims 1 to 2, wherein the step of extracting the first admittance values of each resonator from the result of the multi-physical field simulation comprises the steps of:
and extracting first admittance values of each harmonic oscillator at different frequency points in the pass band from the result of the multi-physical field simulation.
6. The method for extracting parameters of an acoustic wave filter circuit model according to any one of claims 1 to 2, wherein the initial values and/or parameter values of the circuit model parameters are optimized by a numerical optimization algorithm.
7. An apparatus for extracting parameters of a circuit model of an acoustic wave filter, comprising:
the first admittance value obtaining module is used for carrying out multi-physical field simulation on a physical model of the acoustic wave filter, and extracting a first admittance value of each harmonic oscillator from a multi-physical field simulation result, wherein the acoustic wave filter comprises a plurality of harmonic oscillators;
the second admittance value obtaining module is used for respectively optimizing the initial values of the set circuit model parameters of each harmonic oscillator and calculating the second admittance values of each harmonic oscillator according to the optimized parameter values of the circuit model parameters of each harmonic oscillator;
an objective function value obtaining module, configured to input the first admittance value and the second admittance value of each resonator into a preset objective function, respectively, to obtain an objective function value of each resonator, where the objective function value is a difference between the second admittance value and the first admittance value;
the circuit model parameter extraction module is used for extracting the parameter values of the circuit model parameters of the harmonic oscillators corresponding to the objective function values when the objective function values of the harmonic oscillators meet the preset convergence condition, and taking the extracted parameter values of the circuit model parameters of all the harmonic oscillators as the parameter values of the circuit model parameters of the acoustic wave filter;
the objective function is:
wherein, f (C)i0,Cim,Lim,Rim) Is the value of the objective function, Ci0Is the static capacitance of the ith harmonic oscillator, CimIs the equivalent capacitance of the ith harmonic oscillator, LimIs the equivalent inductance of the ith harmonic oscillator, RimIs the resistance of the ith harmonic oscillator, k is the kth frequency point, n is the total number of frequency points,
is the phase value of the second admittance value of the ith harmonic oscillator,
is the phase value, omega, of the first admittance value of the ith harmonic oscillatorkIs the angular frequency of the k-th frequency point, Ai(ωk) Is the amplitude of the second admittance value, A, of the ith harmonic oscillatori′(ωk) The amplitude of the first admittance value of the ith harmonic oscillator;
the objective function can achieve equivalent deformation of the objective function by modifying the logarithmic base or increasing the coefficient.
8. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 6.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method of any of claims 1 to 6 are implemented when the program is executed by the processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710434138.6A CN107169232B (en) | 2017-06-09 | 2017-06-09 | Method and device for extracting circuit model parameters of acoustic wave filter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710434138.6A CN107169232B (en) | 2017-06-09 | 2017-06-09 | Method and device for extracting circuit model parameters of acoustic wave filter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107169232A CN107169232A (en) | 2017-09-15 |
| CN107169232B true CN107169232B (en) | 2020-04-14 |
Family
ID=59825165
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710434138.6A Active CN107169232B (en) | 2017-06-09 | 2017-06-09 | Method and device for extracting circuit model parameters of acoustic wave filter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107169232B (en) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111010143A (en) * | 2019-11-20 | 2020-04-14 | 天津大学 | Bulk acoustic wave filter, method for manufacturing the same, and duplexer |
| CN111898333B (en) * | 2020-06-24 | 2021-12-10 | 北京智芯仿真科技有限公司 | Method and device for extracting frequency response frequency points and calculating response curve of integrated circuit |
| CN113609815B (en) * | 2021-09-16 | 2022-05-13 | 芯和半导体科技(上海)有限公司 | Circuit simulation optimization method and device, computer equipment and storage medium |
| CN114444429B (en) * | 2021-12-22 | 2023-09-01 | 北京无线电测量研究所 | Chip filter manufacturing method and system |
| CN114117690B (en) * | 2022-01-27 | 2022-06-14 | 深圳飞骧科技股份有限公司 | Method for simulating surface acoustic wave filter, related device and storage medium |
| CN115048885B (en) * | 2022-08-12 | 2022-11-15 | 阿里巴巴（中国）有限公司 | Circuit design parameter adjusting method and device, electronic equipment and storage medium |
| CN115495995B (en) * | 2022-10-25 | 2023-07-07 | 深圳飞骧科技股份有限公司 | Method, system, equipment and storage medium for fitting process parameters through simulation test |
| CN118052191A (en) * | 2024-04-16 | 2024-05-17 | 无锡频岢微电子有限公司 | Acoustic model parameter lookup table construction method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101237244A (en) * | 2008-03-11 | 2008-08-06 | 浙江大学 | Establishment method for antenna equivalent circuit model |
| CN103022624A (en) * | 2012-11-20 | 2013-04-03 | 深圳光启创新技术有限公司 | Harmonic oscillator as well as cavity filter and electromagnetic wave equipment thereof |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10007109B2 (en) * | 2013-06-20 | 2018-06-26 | The University Of North Carolina At Charlotte | Wavelength discriminating imaging systems and methods |
-
2017
- 2017-06-09 CN CN201710434138.6A patent/CN107169232B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101237244A (en) * | 2008-03-11 | 2008-08-06 | 浙江大学 | Establishment method for antenna equivalent circuit model |
| CN103022624A (en) * | 2012-11-20 | 2013-04-03 | 深圳光启创新技术有限公司 | Harmonic oscillator as well as cavity filter and electromagnetic wave equipment thereof |
Non-Patent Citations (3)
| Title |
|---|
| X波段交叉耦合型SIR带通滤波器的设计;郭婧; 汪常娥;《微波学报》;20140721;第161-164页 * |
| 基于参数提取和空间映射的滤波器优化设计;陈建忠;《华南理工大学学报(自然科学版)》;20111015;第39卷(第10期);第32-36、73页 * |
| 基于导纳参数的快速滤波器设计方法;李刚;《2015年全国微波毫米波会议论文集》;20150530;第741-743页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107169232A (en) | 2017-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107169232B (en) | Method and device for extracting circuit model parameters of acoustic wave filter | |
| US11036910B2 (en) | Element removal design in microwave filters | |
| JP6906606B6 (en) | Hierarchical cascade in 2D finite element simulation of sonic filter devices | |
| Yarman et al. | High precision LC ladder synthesis part II: Immittance synthesis with transmission zeros at DC and infinity | |
| KR101811528B1 (en) | Simulating effects of temperature on acoustic microwave filters | |
| CN107562990B (en) | BVD model-based BAW filter design method, device and equipment | |
| Uzunov et al. | Extraction of modified butterworth—Van Dyke model of FBAR based on FEM analysis | |
| CN103562736A (en) | Simple and minimally invasive methods and systems for sensing and computing load impedance | |
| RU133668U1 (en) | BAND ACTIVE FILTER | |
| Fan et al. | A sequentially coupled filter design approach using the reflected group delay method and the implicit space mapping technique | |
| Savin et al. | Mode matching technique for tunable shielded cylindrical metal-dielectric resonator | |
| Razinkin et al. | Elliptic filters with quarter-wave coupling | |
| WO2019201454A1 (en) | Method for determining element values of filter structure | |
| CN116579284B (en) | Method, system and related equipment for designing frequency band of surface acoustic wave filter | |
| Zahradnik et al. | Degree of equiripple narrow bandpass FIR filter | |
| Lehtovuori et al. | Improving handset antenna performance by wideband matching optimization | |
| Shrivastava et al. | Application of Newton iterative process in approximation of fractional differentiator and integrator | |
| Yildiz et al. | Design of multiband matching ladders without mutual coupling using parametric representation of Brune functions | |
| Biswal et al. | Design Analysis and Stability Improvement of Fractional-Order Allpass and Bandstop Filters | |
| Krishna | Realization of Fractional order Multivibrator using different approximation techniques | |
| Wen et al. | A Novel Exponential Transform Method for Evaluating the Small Signal Stability with Enhanced Efficiency and Accuracy | |
| KR20210085689A (en) | METHOD OF DETERMINING PARAMETERS OF mBVD MODEL | |
| Class et al. | Patent application title: NETWORK SYNTHESIS DESIGN OF MICROWAVE ACOUSTIC WAVE FILTERS | |
| Hadad | Finding ‘at-once’all dominantly excited modes in a layered medium by fourier-padé approximation of the green's function | |
| CN117879526A (en) | Model correction method of bulk acoustic wave resonator |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20200108 Address after: 510730 Guangdong city of Guangzhou province Guangzhou economic and Technological Development Zone Jinbi Road No. 6 Applicant after: Jingxin Communication Technology (Guangzhou) Co., Ltd. Address before: 510663 Shenzhou Road 10, Guangzhou Science City, Guangzhou economic and Technological Development Zone, Guangzhou, Guangdong Applicant before: Jingxin Communication System (China) Co., Ltd. Applicant before: Jingxin Communication Technology (Guangzhou) Co., Ltd. Applicant before: Jingxin Communication System (Guangzhou) Co., Ltd. Applicant before: TIANJIN COMBA TELECOM SYSTEMS CO., LTD. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |