Method, device, electronic device and storage medium for removing filter loading phase

Abstract

The present application provides a method, device, electronic device and storage medium for removing the loading phase of a filter, and relates to the technical field of filters, wherein the method comprises obtaining initial scattering parameters of the filter; obtaining the loading phase corresponding to the initial reflection parameter for each initial reflection parameter; constructing a correction matrix based on the loading phase; correcting the initial scattering parameters by means of the correction matrix to obtain the corresponding de-loading phase scattering parameters; wherein the following steps are performed to obtain the loading phase corresponding to the initial reflection parameter: performing vector fitting on the initial reflection parameters to obtain a first zero-pole set and a second zero-pole set; determining multiple target zero-pole sets caused by the loading phase based on the first zero-pole set and the second zero-pole set; calculating a phase factor based on the multiple target zero-pole sets; and calculating the corresponding loading phase based on the phase factor.

Description

Technical Field

[0001] The present application relates to the field of filter technology, and in particular to a method, device, electronic device and storage medium for removing a filter loading phase. Background Art

[0002] With the continuous advancement of microwave technology, microwave filters, as an important RF device, have been widely used in the fields of communications, radar, etc. In the design process of microwave filters, accurately extracting the modeling parameters of the filter is a very critical step. Among them, S parameters (scattering parameters) are one of the basic parameters to describe the characteristics of microwave networks, and their accurate measurement is crucial for subsequent filter modeling.

[0003] In the prior art, a variety of methods have been proposed to eliminate the influence of the loading phase on the S parameters. One of the methods is to determine the amount of loading phase removal by curve polynomial fitting of the phase shift and time delay of the actual loaded physical structure. However, due to the small actual loading phase value, the extraction accuracy of this method is not ideal. Another method is to use an optimization algorithm to solve the loading phase. Although this type of method can improve the accuracy of the solution through different optimization strategies, due to the optimization iteration, the algorithm efficiency is relatively low and may fall into the problem of local optimal solution. Another loading phase removal method based on vector fitting lacks adaptive boundary division in its method. Under large detuning conditions, the effectiveness of this method will be affected, and its efficiency will also be affected. Summary of the invention

[0004] The purpose of the embodiments of the present application is to provide a method, device, electronic device and storage medium for removing the filter loading phase, so as to improve the effectiveness of the filter loading phase removal, especially in the face of a large detuning state, so as to determine a more accurate S parameter.

[0005] In a first aspect, the present invention provides a method for removing a filter loading phase, the method comprising:

[0006] Acquire initial scattering parameters of the filter, the initial scattering parameters including initial reflection parameters and initial transmission parameters; for each initial reflection parameter, acquire a loading phase corresponding to the initial reflection parameter; construct a correction matrix based on the loading phase; correct the initial scattering parameters by the correction matrix to obtain the corresponding scattering parameters after removing the loading phase;

[0007] For each initial reflection parameter, the following steps are performed to obtain a loading phase corresponding to the initial reflection parameter:

[0008] The frequency is normalized and the initial reflection parameters are vector-fitted to obtain a first number of zero poles as a first zero-pole set; the frequency is normalized and the initial reflection parameters are vector-fitted to obtain a second number of zero poles as a second zero-pole set; based on the first zero-pole set and the second zero-pole set, multiple target zero poles caused by the loading phase are determined; based on the multiple target zero poles, a phase factor is calculated; based on the phase factor, a corresponding loading phase is calculated.

[0009] In an optional implementation, the step of determining a plurality of target zero-poles caused by the loading phase based on the first zero-pole set and the second zero-pole set specifically includes:

[0010] Determine the intersection between the first zero-pole set and the second zero-pole set; determine the non-repeating zero-pole points between the first zero-pole set and the intersection as target zero-pole points.

[0011] In an optional implementation, vector fitting is performed in the following manner to obtain a preset number of zeros and poles:

[0012] For the residue matrix in the following formula To solve:

[0013] ;

[0014] ;

[0015] in, For the The initial reflection parameters at the frequency point, , is the number of filter ports, is the total number of frequency points, is the first coefficient matrix, is the second coefficient matrix, Based on The diagonal matrix formed;

[0016] Determine the residual matrix corresponding to the residual Is it less than the preset threshold?

[0017] If not, recalculate the poles to update the first coefficient matrix and the second coefficient matrix, and return to the step of solving the residue matrix;

[0018] If so, a preset number of poles and residues are output to obtain a preset number of zero poles.

[0019] In an optional implementation, the first coefficient matrix and the second coefficient matrix are constructed in the following manner:

[0020] The first coefficient matrix in the second matrix is ​​obtained by Line Column element values:

[0021] Calculate the The complex frequency of the frequency point is The difference between the extreme points;

[0022] Calculate the ratio between 1 and the difference as the first Line Column element values:

[0023] Among them, a preset number is randomly generated within a preset range The initial point of ;

[0024] The complex frequency of each frequency point can be obtained by the following method :

[0025] ;

[0026] ;

[0027] in, For the The complex frequency of the frequency point, is the center frequency, is the starting frequency, is the cut-off frequency, is a unit imaginary number, is bandwidth;

[0028] And in the second coefficient matrix To the right of the column, add the preset vertical matrix to obtain the first coefficient matrix. The preset vertical matrix includes element values, each element value is 1.

[0029] In an optional implementation, the phase factor is calculated in the following manner:

[0030] ;

[0031] in, For the The zero point of the target zero pole, For the The poles of the target zero poles.

[0032] In an optional embodiment, the loading phase is calculated by :

[0033] ;

[0034] in, Used for argument calculations.

[0035] In an optional implementation, for each frequency point, the initial scattering parameter , the deloading phase scattering parameters are generated by :

[0036] ;

[0037] .

[0038] In a second aspect, the present invention provides a device for removing a filter loading phase, the device comprising:

[0039] An acquisition module, used to obtain initial scattering parameters of the filter;

[0040] The processing module is used to obtain, for each initial scattering parameter, a loading phase corresponding to the initial scattering parameter, wherein, for each initial scattering parameter, the following steps are performed to obtain the loading phase corresponding to the initial scattering parameter:

[0041] Performing vector fitting on the scattering parameters to obtain a first number of zero poles as a first zero pole set; performing vector fitting on the scattering parameters to obtain a second number of zero poles as a second zero pole set; determining multiple target zero poles caused by the loading phase based on the first zero pole set and the second zero pole set; calculating a phase factor based on the multiple target zero poles; and calculating the loading phase based on the phase factor;

[0042] A correction module, for constructing a correction matrix based on a loading phase;

[0043] The deloading phase module is used to correct the initial scattering parameters through a correction matrix to obtain corresponding deloading phase scattering parameters.

[0044] In a third aspect, the present invention provides an electronic device comprising: a processor, a memory and a bus, the memory storing machine-readable instructions executable by the processor, and when the electronic device is running, the processor and the memory communicate through the bus, and the processor executes the machine-readable instructions to perform the steps of the method for removing the filter loading phase in any of the aforementioned implementations.

[0045] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the method for removing the filter loading phase in any of the aforementioned embodiments are executed.

[0046] The present application provides a filter loading phase removal method, device, electronic device and storage medium, wherein the method includes obtaining the initial scattering parameters of the filter, the initial scattering parameters including the initial reflection parameters and the initial transmission parameters; for each initial reflection parameter, obtaining the loading phase corresponding to the initial reflection parameter; constructing a correction matrix based on the loading phase; correcting the initial scattering parameters by the correction matrix to obtain the corresponding scattering parameters of the de-loading phase; wherein, for each initial reflection parameter, the following steps are performed to obtain the loading phase corresponding to the initial reflection parameter: vector fitting the initial reflection parameters to obtain a first number of zero poles as a first zero pole set; vector fitting the initial reflection parameters to obtain a second number of zero poles as a second zero pole set; based on the first zero pole set and the second zero pole set, determining multiple target zero poles caused by the loading phase; calculating the phase factor based on the multiple target zero poles; and calculating the corresponding loading phase based on the phase factor. By analyzing the influence of the loading phase from a mathematical perspective and then correcting the S parameter, the calculation efficiency is high, the time consumption is less, and the loading phase removal accuracy is high. It can be applied to the loading phase removal under various working conditions, especially for large detuning conditions. It can also show good removal ability. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments of the present application will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present application and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0048] Figure 1 A flowchart of a method for removing a filter loading phase provided in an embodiment of the present application;

[0049] Figure 2 A flowchart of steps for determining a loading phase provided in an embodiment of the present application;

[0050] Figure 3 A flowchart of a vector fitting step provided in an embodiment of the present application;

[0051] Figure 4 A schematic diagram of the structure of a filter loading phase removal device provided in an embodiment of the present application;

[0052] Figure 5 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application;

[0053] Figure 6 A schematic diagram of transmission parameters of a three-cavity waveguide filter provided in an embodiment of the present application;

[0054] Figure 7 A zero-pole distribution diagram provided in an embodiment of the present application;

[0055] Figure 8 Another zero-pole distribution diagram provided in an embodiment of the present application;

[0056] Fig. 9 The first embodiment provided in this application Schematic diagram of the fitting;

[0057] Fig.10 The first embodiment provided in this application Schematic diagram of the fitting;

[0058] Fig.11 The second embodiment provided in this application Schematic diagram of the fitting;

[0059] Fig.12 The second embodiment provided in this application Schematic diagram of the fitting;

[0060] Fig.13 The first embodiment provided in this application Schematic diagram of the real part of ;

[0061] Fig.14 The first embodiment provided in this application Schematic diagram of the imaginary part of ;

[0062] Fig.15 The second embodiment provided in this application Schematic diagram of the real part of ;

[0063] Fig.16 The second embodiment provided in this application Schematic diagram of the imaginary part of ;

[0064] Fig.17 The first embodiment provided in this application Schematic diagram of the fitting error;

[0065] Fig.18 The first embodiment provided in this application Schematic diagram of the fitting error;

[0066] Fig.19 The second embodiment provided in this application Schematic diagram of the fitting error;

[0067] Fig. 20 The second embodiment provided in this application Schematic diagram of the fitting error. DETAILED DESCRIPTION

[0068] First, the application scenario of the present application is described. The technical solution of the present application is applicable to the loading phase removal of S parameters in filter modeling simulation, and can be used to form a coupling matrix to complete the design and optimization of the filter.

[0069] In the design process of microwave filters, accurately extracting the modeling parameters of the filter is a very critical step. Among them, S parameters (scattering parameters) are one of the basic parameters that describe the characteristics of microwave networks, and their accurate measurement is crucial for subsequent filter modeling.

[0070] In the process of filter design and debugging, since the port of the filter is connected to a transmission line, an additional phase shift is introduced, which causes the port plane to be not a parameter plane. Therefore, it is very important to eliminate the loading phase in the S parameter, which directly affects the extraction accuracy of the coupling matrix.

[0071] Based on this, the present application provides a method, device, electronic device and storage medium for removing filter loading phase.

[0072] The technical solutions in the embodiments of the present application will be described below in conjunction with the drawings in the embodiments of the present application.

[0073] Figure 1 Flow chart of a method for removing filter loading phase provided in an embodiment of the present application. Figure 1 As shown, the present application provides a method for removing a filter loading phase, comprising:

[0074] S1. Obtaining initial scattering parameters of the filter, where the initial scattering parameters include initial reflection parameters and initial transmission parameters.

[0075] Here, the target filter can be modeled and the initial scattering parameters can be obtained through simulation. Ports can include input ports and output ports for connecting with other devices such as signal sources, loads, filters, etc.

[0076] The reference plane is used to measure multiple parameters of the filter, such as S parameters (scattering parameters). The S parameters here include the dispersion effects of input and output high-order modes and the phase shift of the port itself.

[0077] For example, taking a two-port three-cavity waveguide filter as an example, including an input port and an output port, the measured S parameters may include reflection parameters corresponding to the sampling frequency point. and , transmission parameters and .

[0078] S2. For each initial reflection parameter, obtain a loading phase corresponding to the initial reflection parameter.

[0079] Due to the relatively insignificant change in the coupling parameters of the filter cavity with frequency under large detuning conditions, the zeros and poles that characterize the ideal S parameters are not all near the origin. This results in the dilemma of being unable to correctly classify the zeros and poles. Therefore, an improved zero-pole classification step is provided here to calculate the loading phase.

[0080] Specifically, in step S2, only the reflection parameters of the input port are and reflection parameters of the output port , calculate the loading phase.

[0081] like Figure 2 As shown, for each initial reflection parameter, the following steps are performed to obtain the loading phase corresponding to the initial reflection parameter:

[0082] S20. Perform vector fitting on the initial reflection parameters to obtain a first number of zeros and poles as a first zero-pole set.

[0083] S22. Perform vector fitting on the initial reflection parameters to obtain a second number of zero-pole points as a second zero-pole set.

[0084] In step S20 and step S22, a VF algorithm may be used to perform vector fitting on the initial reflection parameters of each frequency point to obtain a preset number of zero poles. The preset number here includes a first number and a second number. The second number is greater than the first number.

[0085] S24. Determine a plurality of target zero-pole points caused by the loading phase based on the first zero-pole set and the second zero-pole set.

[0086] In step S24, based on the first zero-pole set and the second zero-pole set, a step of determining a plurality of target zero-poles caused by the loading phase specifically includes:

[0087] Determine the intersection between the first zero-pole set and the second zero-pole set. Determine the non-repeated zero-pole points between the first zero-pole set and the intersection as target zero-pole points.

[0088] Here, the zero poles in the intersection between the first zero pole set and the second zero pole set are the zero poles after the loading phase is removed. The multiple target zero poles that are not repeated between the first zero pole set and the intersection are the zero poles with the loading phase.

[0089] In a specific embodiment, the first number is , the second quantity is ,in, is the filter order.

[0090] S26. Calculate a phase factor based on multiple target zero poles.

[0091] Here, the phase factor can be calculated as follows:

[0092] ;

[0093] in, For the The zero point of the target zero pole, For the The poles of the target zero poles, is the number of target zeros and poles, is the complex frequency.

[0094] S28. Calculate the corresponding loading phase based on the phase factor.

[0095] Here, the loading phase can be calculated as follows :

[0096] ;

[0097] in, Used for argument calculations.

[0098] That is to say, taking the two-port three-cavity waveguide filter as an example, the reflection parameters can be obtained respectively Corresponding and reflection parameters of the output port Corresponding .

[0099] S3. Construct a correction matrix based on the loading phase.

[0100] Correction Matrix It can be expressed as:

[0101] .

[0102] Taking a two-port three-cavity waveguide filter as an example, the correction matrix .

[0103] S4. Correct the initial scattering parameters by using a correction matrix to obtain corresponding optimized scattering parameters.

[0104] In step S4, the initial scattering parameters corresponding to each frequency point are , the deloading phase scattering parameters are generated by :

[0105] ;

[0106] .

[0107] in, For the The initial scattering parameters corresponding to the frequency point are: For the The optimized scattering parameters corresponding to each frequency point.

[0108] The scattering parameters of the unloading phase after removing the loading phase are used to form a coupling matrix, thereby completing the design and optimization of the filter.

[0109] The present application provides a method for removing the loading phase of a filter. By analyzing the influence of the loading phase from a mathematical perspective and then correcting the S parameters, compared with the curve polynomial fitting and optimization algorithm solution methods in the prior art, the method has high computational efficiency, less time consumption, and high loading phase removal accuracy. It can be applied to loading phase removal under various working conditions, and can also exhibit good removal capabilities under large detuning conditions.

[0110] like Figure 3 As shown, in a specific embodiment of the present application, vector fitting can be performed in the following manner to obtain a preset number of zeros and poles:

[0111] S2000, generate initial poles.

[0112] When determining the initial poles, the imaginary part of the poles needs to be evenly distributed in [-1,1.1], and the real part is 1% of the imaginary part. For example, the initial poles can be , .

[0113] S2002, assemble the matrix.

[0114] ;

[0115] ;

[0116] in, For the The initial reflection parameters at the frequency point, , is the number of filter ports, is the total number of frequency points, is the first coefficient matrix, is the second coefficient matrix, Based on The resulting diagonal matrix. , the dimension is .

[0117] In step S2002, the first coefficient matrix and the second coefficient matrix may be constructed in the following manner:

[0118] Specifically, the first coefficient matrix in the second coefficient matrix is ​​obtained by the following method: Line Column element values:

[0119] Calculate the The complex frequency of the frequency point is Calculate the ratio between 1 and the difference as the value of the first Line The element value of the column. Among them, a preset number of random The initial point of .

[0120] That is, the second coefficient matrix The row can be represented as .

[0121] It should be noted that physical frequency normalization can be performed on each frequency point before vector fitting. Specifically, the complex frequency of each frequency point can be obtained by the following method: :

[0122] ;

[0123] ;

[0124] in, For the The complex frequency of the frequency point, is the center frequency, is the starting frequency, is the cut-off frequency, is a unit imaginary number, For bandwidth.

[0125] And in the second coefficient matrix To the right of the column, add the preset vertical matrix to obtain the first coefficient matrix. The preset vertical matrix includes element values, each element value is 1.

[0126] That is, the first coefficient matrix can be expressed as:

[0127] .

[0128] S2004, based on the least squares method of the residue matrix To solve.

[0129] Assumptions , based on the least squares method, solve the following formula to obtain the residue matrix:

[0130] .

[0131] S2006. Determine the residual corresponding to the residue matrix Is it less than the preset threshold?

[0132] In a specific embodiment, the preset threshold may be 1e-4. for The corresponding L2 norm.

[0133] S2008: If not, recalculate the poles to update the first coefficient matrix and the second coefficient matrix.

[0134] In step S2008, Solve the eigenvalue to obtain a new pole, and then return to execute step S2002 again. , the dimension is .

[0135] S2010: If yes, output a preset number of poles and residues to obtain a preset number of zero poles.

[0136] Furthermore, it is also possible to determine whether a preset number of iterations has been reached, and if so, the result can be output.

[0137] The present application provides a method for removing the loaded phase of a filter. Based on the classification of zeros and poles based on mathematical principles, the method can effectively distinguish the zeros and poles of the loaded phase from the zeros and poles of the real S parameters, thereby effectively removing the loaded phase of the S parameters. The method has strong robustness and high accuracy even under large detuning conditions.

[0138] In a specific embodiment of the present application, the three-cavity waveguide filter is simulated to obtain the response of the three-cavity waveguide filter in a certain state. Figure 6 shown.

[0139] The VF algorithm is used to perform vector fitting on the initial reflection parameters of each frequency point to obtain the first number of zero poles. Figure 7 The distribution of zero poles in .

[0140] The VF algorithm is used to perform vector fitting on the initial reflection parameters of each frequency point to obtain the second number of zero poles, which can be seen in Figure 8 The distribution of zero poles in .

[0141] It can be seen that the zero poles generated by the loading phase are close to the zero poles generated by the inherent S parameters (S parameters after removing the loading phase), and cannot be directly distinguished using a circle with a fixed radius. If the optimization method is used to solve the loading phase, it will be difficult to converge to the global optimal solution, and the convergence will be slow. If the existing loading phase removal method is used to remove the loading phase, the final fitting S parameter result is as follows Fig. 9 and Fig.10 As shown, the curves do not overlap highly. (The coupling matrix corresponding to the S parameters after removing the loading phase is extracted, and the S parameters represented by the coupling matrix are the final fitted S parameters. Under the premise that the coupling matrix is ​​correctly extracted, the fitted S parameters are equal to the S parameters after removing the loading phase.)

[0142] After classification by the loading phase removal method provided in the embodiment of the present application, the zero poles generated by the loading phase can be effectively identified, and the obtained loading phase is as follows: Fig.11 and Fig.12 shown.

[0143] It can be seen that Fig.11 The two curves completely overlap. Fig.12 Although the two curves do not completely overlap, they completely overlap above -35db, and differences are allowed in places with smaller db values.

[0144] The two can be presented in the form of plural contrast. Fig.13 and Fig.14 After removing the loading phase using existing technology The fitting of the real and imaginary parts. Fig.15 and 16 After the loading phase removal method provided in the embodiment of the present application is adopted The fitting of the real and imaginary parts.

[0145] Furthermore, the error value after removing the loading phase using the existing technology can be determined. , For example, The fitting error can be expressed as:

[0146] .

[0147] like Fig.17 As shown, it can be seen that the maximum fitting error reaches 0.9, indicating that there is a large phase error.

[0148] The fitting error can be expressed as:

[0149] 。

[0150] like Fig.18 As shown, it can be seen that the maximum fitting error reaches 1.6, indicating that there is a large phase error.

[0151] After adopting the loading phase removal method provided in the embodiment of the present application, The fitting error is Fig.19 As shown, the errors are This indicates that there is almost no phase error. The fitting error is Fig. 20 As shown, the errors are Whether from the theory or the numerical experimental results, it can be seen that the loading phase removal method provided in the embodiment of the present application can accurately extract the loading phase of the filter and has strong robustness.

[0152] like Figure 4 As shown, based on the same inventive concept, an embodiment of the present application further provides a filter loading phase removal device 400, the device comprising:

[0153] An acquisition module 4100 is used to obtain initial scattering parameters corresponding to a target reference plane of the filter, where the target reference plane is determined according to the number of ports of the filter;

[0154] The processing module 4200 is used to obtain, for each initial scattering parameter, a loading phase corresponding to the initial scattering parameter, wherein, for each initial scattering parameter, the following steps are performed to obtain the loading phase corresponding to the initial scattering parameter:

[0155] Performing vector fitting on the normalized scattering parameters to obtain a first number of zero poles as a first zero pole set; performing vector fitting on the normalized scattering parameters to obtain a second number of zero poles as a second zero pole set; determining multiple target zero poles caused by the loading phase based on the first zero pole set and the second zero pole set; calculating a phase factor based on the multiple target zero poles; and calculating the loading phase based on the phase factor;

[0156] A correction module 4300, for constructing a correction matrix based on the loading phase;

[0157] The deloading phase module 4400 is used to correct the initial scattering parameters through the correction matrix to obtain the corresponding optimized scattering parameters.

[0158] In a preferred embodiment, the processing module 4200 is specifically used for:

[0159] Determine the intersection between the first zero-pole set and the second zero-pole set; determine the non-repeating zero-pole points between the first zero-pole set and the intersection as target zero-pole points.

[0160] In a preferred embodiment, the processing module 4200 performs vector fitting in the following manner to obtain a preset number of zeros and poles:

[0161] For the residue matrix in the following formula To solve:

[0162] ;

[0163] ;

[0164] in, For the The initial reflection parameters at the frequency point, , is the number of filter ports, is the total number of frequency points, is the first coefficient matrix, is the second coefficient matrix, Based on The diagonal matrix formed;

[0165] Determine the residual matrix corresponding to the residual Is it less than the preset threshold?

[0166] If not, recalculate the poles to update the first coefficient matrix and the second coefficient matrix, and return to the step of solving the residue matrix;

[0167] If so, a preset number of poles and residues are output to obtain a preset number of zero poles.

[0168] In a preferred embodiment, the processing module 4200 constructs the first coefficient matrix and the second coefficient matrix in the following manner:

[0169] The first coefficient matrix in the second matrix is ​​obtained by Line Column element values:

[0170] Calculate the The complex frequency of the frequency point is The difference between the extreme points;

[0171] Calculate the ratio between 1 and the difference as the first Line Column element values:

[0172] Among them, a preset number is randomly generated within a preset range The initial point of ;

[0173] The complex frequency of each frequency point can be obtained by the following method :

[0174] ;

[0175] ;

[0176] in, For the The complex frequency of the frequency point, is the center frequency, is the starting frequency, is the cut-off frequency, is a unit imaginary number, is bandwidth;

[0177] And in the second coefficient matrix To the right of the column, add the preset vertical matrix to obtain the first coefficient matrix. The preset vertical matrix includes element values, each element value is 1.

[0178] In a preferred embodiment, the processing module 4200 calculates the phase factor in the following manner:

[0179] ;

[0180] in, For the The zero point of the target zero pole, For the The poles of the target zero poles.

[0181] In a preferred embodiment, the processing module 4200 calculates the loading phase by the following method: :

[0182] ;

[0183] in, Used for argument calculations.

[0184] In a preferred embodiment, the initial scattering parameters corresponding to each frequency point are The correction module 4300 generates the scattering parameters of the unloading phase by the following method: :

[0185] ;

[0186] .

[0187] See also Figure 5 , Figure 5 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. Figure 5As shown in , the electronic device 500 includes a processor 510 , a memory 520 and a bus 530 .

[0188] The memory 520 stores machine-readable instructions executable by the processor 510. When the electronic device 500 is running, the processor 510 communicates with the memory 520 through the bus 530. When the machine-readable instructions are executed by the processor 510, the steps of a method for removing a filter loading phase in the above-mentioned method embodiment can be executed. The specific implementation method can be found in the method embodiment, which will not be repeated here.

[0189] An embodiment of the present application also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of a method for removing a filter loading phase in the above-mentioned method embodiment can be executed. The specific implementation method can be found in the method embodiment, which will not be repeated here.

[0190] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0191] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some communication interfaces, and the indirect coupling or communication connection of the devices or units can be electrical, mechanical or other forms.

[0192] In addition, the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0193] Furthermore, the functional modules in the various embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.

[0194] It should be noted that if the function is implemented in the form of a software function module and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM), random access memory (RAM), disk or optical disk, and other media that can store program codes.

[0195] In this document, relational terms such as first and second, etc. are used merely to distinguish one entity or operation from another entity or operation, but do not necessarily require or imply any such actual relationship or order between these entities or operations.

[0196] The above description is only an embodiment of the present application and is not intended to limit the protection scope of the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for removing filter loading phase, characterized in that: The method comprises: Acquire initial scattering parameters of the filter, wherein the initial scattering parameters include initial reflection parameters and initial transmission parameters; For each initial reflection parameter, obtaining a loading phase corresponding to the initial reflection parameter; Constructing a correction matrix based on the loading phase; The initial scattering parameters are corrected by using a correction matrix to obtain corresponding unloading phase scattering parameters; For each initial reflection parameter, the following steps are performed to obtain a loading phase corresponding to the initial reflection parameter: Performing vector fitting on the initial reflection parameters to obtain a first number of zero poles as a first zero pole set; Performing vector fitting on the initial reflection parameters to obtain a second number of zero poles as a second zero pole set; Based on the first zero-pole set and the second zero-pole set, determining a plurality of target zero-pole points caused by the loading phase; Calculate the phase factor based on multiple target zero poles; Calculating a corresponding loading phase based on the phase factor; The step of determining a plurality of target zero-poles caused by the loading phase based on the first zero-pole set and the second zero-pole set specifically includes: Determine an intersection between the first pole-zero set and the second pole-zero set; Determine the non-repeated zero-pole points between the first zero-pole set and the intersection as target zero-pole points; Perform vector fitting to obtain a preset number of poles and zeros by: For the residue matrix in the following formula To solve: ； ； in, For the The initial reflection parameters at the frequency point, , is the number of filter ports, is the total number of frequency points, is the first coefficient matrix, is the second coefficient matrix, Based on The diagonal matrix formed; Determine the residual corresponding to the residue matrix Is it less than the preset threshold? If not, recalculate the poles to update the first coefficient matrix and the second coefficient matrix, and return to the step of solving the residue matrix; If yes, output a preset number of poles and residues to obtain a preset number of zero poles; The first coefficient matrix and the second coefficient matrix are constructed in the following way: The first coefficient matrix in the second matrix is ​​obtained by Line Column element values: Calculate the The complex frequency of the frequency point is The difference between the extreme points; Calculate the ratio between 1 and the difference as the first Line Column element values: Among them, a preset number is randomly generated within a preset range The initial point of ; The complex frequency of each frequency point can be obtained by the following method : ； ； in, For the The complex frequency of the frequency point, is the center frequency, is the starting frequency, is the cut-off frequency, is a unit imaginary number, is bandwidth; And in the second coefficient matrix To the right of the column, add a preset longitudinal matrix to obtain a first coefficient matrix, wherein the preset longitudinal matrix includes element values, each element value is 1.

2. The method according to claim 1, characterized in that The phase factor is calculated as follows: ； in, For the The zero point of the target zero pole, For the The poles of the target zero poles.

3. The method according to claim 2, characterized in that The loading phase is calculated by : ； in, Used for argument calculations.

4. The method according to claim 3, characterized in that Initial scattering parameters corresponding to each frequency point , the deloading phase scattering parameters are generated by : ； 。 5. A filter loading phase removal device, characterized in that: The device comprises: An acquisition module, used to obtain initial scattering parameters of the filter; The processing module is used to obtain, for each initial scattering parameter, a loading phase corresponding to the initial scattering parameter, wherein, for each initial scattering parameter, the following steps are performed to obtain the loading phase corresponding to the initial scattering parameter: Performing vector fitting on scattering parameters to obtain a first number of zero poles as a first zero pole set; performing vector fitting on scattering parameters to obtain a second number of zero poles as a second zero pole set; determining multiple target zero poles caused by loading phase based on the first zero pole set and the second zero pole set; calculating a phase factor based on the multiple target zero poles; calculating the loading phase based on the phase factor; the step of determining multiple target zero poles caused by loading phase based on the first zero pole set and the second zero pole set specifically includes: Determine the intersection between the first zero-pole set and the second zero-pole set; determine the non-overlapping zero-pole points between the first zero-pole set and the intersection as target zero-pole points; The processing module performs vector fitting to obtain a preset number of zeros and poles in the following way: For the residue matrix in the following formula To solve: ； ； in, For the The initial reflection parameters at the frequency point, , is the number of filter ports, is the total number of frequency points, is the first coefficient matrix, is the second coefficient matrix, Based on The diagonal matrix formed; Determine the residual corresponding to the residue matrix Is it less than the preset threshold? If not, recalculate the poles to update the first coefficient matrix and the second coefficient matrix, and return to the step of solving the residue matrix; If yes, output a preset number of poles and residues to obtain a preset number of zero poles; The processing module constructs the first coefficient matrix and the second coefficient matrix in the following manner: The first coefficient matrix in the second matrix is ​​obtained by Line Column element values: Calculate the The complex frequency of the frequency point is The difference between the extreme points; Calculate the ratio between 1 and the difference as the first Line Column element values: Among them, a preset number is randomly generated within a preset range The initial point of ; The complex frequency of each frequency point can be obtained by the following method : ； ； in, For the The complex frequency of the frequency point, is the center frequency, is the starting frequency, is the cut-off frequency, is a unit imaginary number, is bandwidth; And in the second coefficient matrix To the right of the column, add a preset longitudinal matrix to obtain a first coefficient matrix, wherein the preset longitudinal matrix includes element values, each element value is 1; A correction module, for constructing a correction matrix based on a loading phase; The deloading phase module is used to correct the initial scattering parameters through a correction matrix to obtain corresponding deloading phase scattering parameters.

6. An electronic device, characterized in that: include: A processor, a memory and a bus, wherein the memory stores machine-readable instructions executable by the processor, and when the electronic device is running, the processor communicates with the memory via the bus, and the processor executes the machine-readable instructions to perform the steps of the method for removing the filter loading phase as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for removing the filter loading phase according to any one of claims 1 to 4 are executed.

Images