CN111797541A - Method for quantifying coupling uncertainty and calculating global sensitivity of field line - Google Patents
Method for quantifying coupling uncertainty and calculating global sensitivity of field line Download PDFInfo
- Publication number
- CN111797541A CN111797541A CN202010696773.3A CN202010696773A CN111797541A CN 111797541 A CN111797541 A CN 111797541A CN 202010696773 A CN202010696773 A CN 202010696773A CN 111797541 A CN111797541 A CN 111797541A
- Authority
- CN
- China
- Prior art keywords
- polynomial
- hyperbolic
- chaotic polynomial
- truncation
- sensitivity
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Complex Calculations (AREA)
Abstract
The invention relates to a field line coupling uncertainty quantification and global sensitivity calculation method, which combines an adaptive hyperbolic truncation scheme with a minimum angle regression method, sparsely processes a generalized chaotic polynomial, establishes a proxy model to analyze uncertainty of radiation sensitivity of a multi-conductor transmission line, and analyzes global sensitivity of an input variable by combining a global sensitivity analysis method.
Description
Technical Field
The invention relates to the field of multi-conductor transmission line radiation sensitivity, in particular to a field line coupling uncertainty quantification and global sensitivity calculation method.
Background
With the continuous development of science and technology, various electrical systems and electronic devices are widely applied to real life, resulting in the increasingly complex electromagnetic environment in which the electrical devices are located. The transmission line is used as a carrier for transmitting electric energy and signals and is an indispensable part in the electric and electronic equipment; the problem of electromagnetic compatibility of transmission lines is also becoming increasingly serious due to the complexity of the electromagnetic environment. Radiation sensitivity is also gaining increasing attention from researchers as a class of important issues in electromagnetic compatibility.
Due to the possible changes of the electric and electronic equipment and the complexity of the electromagnetic environment, the problem of the radiation sensitivity of the transmission line has strong uncertainty, and the research on the problem of the uncertainty of the radiation sensitivity of the multi-conductor transmission line is more and more important. Different uncertain input variables obey corresponding random distribution, and the analysis of the influence degree of the different uncertain input variables on the whole system is also necessary.
The traditional Monte Carlo method has high calculation cost and low calculation efficiency, does not carry out global sensitivity analysis on uncertain input variables, and cannot accurately analyze the influence degree of each uncertain input variable on the whole system.
Disclosure of Invention
The invention aims to provide a method for accurately analyzing the influence degree of each uncertain input variable on the whole system.
In order to achieve the purpose, the invention provides the following scheme:
a field line coupling uncertainty quantification and global sensitivity calculation method, the method comprising:
determining an input variable according to an electromagnetic environment and an incident field where the multi-conductor transmission line is located;
determining a distribution type corresponding to the input variable according to the input variable;
determining an orthogonal base corresponding to the input variable according to the distribution type corresponding to the input variable;
establishing a generalized chaotic polynomial model of the radiation sensitivity of the multi-conductor transmission line by combining a transmission line theory according to the orthogonal substrate of the input variable to obtain an expansion of the generalized chaotic polynomial;
performing hyperbolic truncation processing on the expanded form of the generalized chaotic polynomial by using a self-adaptive hyperbolic truncation method to obtain the hyperbolic truncated generalized chaotic polynomial;
processing the hyperbolic truncated generalized chaotic polynomial by using a minimum angle regression method to obtain a self-adaptive sparse chaotic polynomial;
calculating induced current or induced voltage of the multi-conductor transmission line radiation sensitivity by using the self-adaptive sparse chaotic polynomial to obtain probability distribution of the induced current or the induced voltage at different frequency points;
converting the self-adaptive sparse chaotic polynomial into an incremental summation form, and analyzing the global sensitivity index of the input variable by using a global sensitivity analysis method; the global sensitivity indicators include: total sensitivity index and first order sensitivity index.
Optionally, the expansion of the generalized chaotic polynomial is:
wherein the content of the first and second substances,mixed orthogonal polynomials of order n, which are multidimensional standard random variablesInfinity is the highest power of the mixed orthogonal polynomial,and phiiAre respectively connected withAndin response to this, the mobile terminal is allowed to,for the generalized mixtureCoefficient of expansion of chaos polynomial, phiiAnd (xi) is the product of the one-dimensional orthogonal polynomial basis functions corresponding to the random variables.
Optionally, the hyperbolic truncation processing is performed on the expansion of the generalized chaotic polynomial by using a self-adaptive hyperbolic truncation method to obtain a hyperbolic truncated generalized chaotic polynomial, and the method specifically includes:
obtaining an expansion of the generalized chaotic polynomial;
setting the highest truncation order P of the self-adaptive hyperbolic truncation method, and then performing primary sparse processing on the expansion of the generalized chaotic polynomial to obtain the primarily sparse generalized chaotic polynomial; the ith polynomial maximum truncation order piSatisfies the following conditions:where k is the dimension of the random variable, lkAnd the order of the random variable in the k dimension is shown, q is a norm, and n is the maximum value of the dimension of the random variable.
Optionally, the method for determining the norm q includes:
setting the range of the preliminary norm q as [ a, b ] and the step length as s;
setting an error threshold;
the preliminary norm q is gradually increased by the step length s and is according to the formulaCalculating out a leave-one-out cross validation error;
stopping calculation when the cross validation error is smaller than the error threshold or exceeds the maximum value of the range of the preliminary norm q to obtain the norm q; wherein the content of the first and second substances,LOOto leave one out cross validation errors, M (x)(i)) For hyperbolic truncation model after setting norm q value at ith sample point x(i)Response value at Point, MPC\i(x(i)) Is a chaotic polynomial model after P-order truncation at the ith sample point x(i)The value of the response of (c) to (d),is the mean of the response values for the L sample points.
Optionally, the hyperbolic truncated generalized chaotic polynomial is processed by using a minimum angle regression method to obtain a self-adaptive sparse chaotic polynomial, which specifically includes:
acquiring the hyperbolic truncated generalized chaotic polynomial;
converting the hyperbolic truncated generalized chaotic polynomial into a Y ═ theta X model; wherein, theta is a set of coefficients of the hyperbolic truncated generalized chaotic polynomial, and X is a set of the hyperbolic truncated generalized chaotic polynomial;
processing the Y-theta X model by using a minimum angle regression method to obtain a coefficient of a polynomial;
and substituting the coefficients of the polynomial into a self-adaptive sparse chaotic polynomial to obtain the self-adaptive sparse chaotic polynomial.
Optionally, the calculating the induced current or the induced voltage of the multi-conductor transmission line radiation sensitivity by using the adaptive sparse chaotic polynomial to obtain the probability distribution of the induced current or the probability distribution of the induced voltage at different frequency points specifically includes:
acquiring the self-adaptive sparse chaotic polynomial;
converting the adaptive sparse chaotic polynomial intoWhereinFor the sparsely processed polynomial coefficient, phii'is a polynomial after sparse processing, and p' is the number of terms of the polynomial after sparse processing;
obtaining the mean value of the multi-conductor transmission line radiation sensitivity induction current or the mean value of the multi-conductor transmission line radiation sensitivity induction voltage according to the formula after the hyperbolic truncation and the generalized chaotic polynomial conversion:wherein the content of the first and second substances,a constant term of an expansion of the truncated generalized chaotic polynomial;
obtaining the variance of the multi-conductor transmission line radiation sensitivity induction current or the variance of the multi-conductor transmission line radiation sensitivity induction voltage according to the formula after the hyperbolic truncation and the generalized chaotic polynomial conversion:
optionally, the converting the adaptive sparse chaotic polynomial into an incremental summation form, and analyzing the global sensitivity index of the input variable by using a global sensitivity analysis method specifically includes:
acquiring the self-adaptive sparse chaotic polynomial;
converting the adaptive sparse chaotic polynomial into an expansion in an incremental summation form;
carrying out recursive calculation on the expansion of the incremental summation form by using an integral method, solving the corresponding coefficient of the decomposition term, and obtaining the expansion of the incremental summation form with known coefficient;
taking variance from both sides of the expansion of the known incremental summation form of the coefficient at the same time to obtain a variance decomposition formula;
obtaining a global sensitivity index according to the variance decomposition formula; the global sensitivity index includes a sensitivity index and a total sensitivity index.
Optionally, the expansion of the incremental summation form is:
optionally, the sensitivity index is:
The total sensitivity index is:
wherein S isiIs an index of first-order sensitivity.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects: the method combines the self-adaptive hyperbolic truncation method with the minimum angle regression method, sparsely processes the generalized chaotic polynomial, establishes a proxy model to analyze the uncertainty problem of the radiation sensitivity of the multi-conductor transmission line, analyzes the global sensitivity of the input variable by combining the Sobol method, and analyzes and calculates the uncertainty problem of the radiation sensitivity of the multi-conductor transmission line with lower calculation cost and higher efficiency on the premise of ensuring the calculation accuracy.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a schematic flow chart of a field line coupling uncertainty quantification and global sensitivity calculation method provided by the present invention;
FIG. 2 is a schematic diagram showing the positional relationship between a multi-conductor transmission line and an incident field;
FIG. 3 is a hyperbolic truncation diagram of two-dimensional input variables and different q-norms;
fig. 4 shows the truncation effect of hyperbolic truncation when P is 5 and n is 2;
FIG. 5 is a schematic diagram of a two-dimensional variable least angle regression method;
FIG. 6 is a 3+1 multi-conductor transmission line model;
FIG. 7 is a diagram illustrating comparison results between the calculation results of the mean and standard deviation of the conventional generalized chaotic polynomial with different truncation orders and the Monte Carlo method of 20000 orders;
FIG. 8 is a pair of different truncation orders at different frequency pointsLOOAn influence schematic;
FIG. 9 is a schematic view ofLOOAn error curve;
fig. 10 is a graph illustrating the comparison result of the standard deviation of the mean of the sparse chaotic polynomial with the monte carlo method when q is 0.8 and P is 15;
FIG. 11 is a comparison result of AS-PC calculating # 2 transmission line induced current probability distribution at different frequency points;
FIG. 12 is a graph comparing an upper limit of the induced current curve with a simulation curve of 20000 Monte Carlo;
FIG. 13 is a graph comparing total sensitivity indices;
FIG. 14 is a first order sensitivity index comparison plot;
FIG. 15 shows the degree of influence of various parameters in the [10MHz,1GHz ] frequency range on the radiation sensitivity of a multi-conductor transmission line.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a field line coupling uncertainty quantification and global sensitivity calculation method, which can be used for analyzing and calculating the uncertainty of the radiation sensitivity of a multi-conductor transmission line with lower calculation cost and higher efficiency on the premise of ensuring the calculation accuracy.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
As shown in fig. 1, a flow chart of a field line coupling uncertainty quantification and global sensitivity calculation method provided by the present invention is shown, and the method specifically includes:
step 101: the input variables are determined from the electromagnetic environment in which the multi-conductor transmission line is located and the incident field.
Step 102: and determining the distribution type corresponding to the input variable according to the input variable.
Step 103: and determining the orthogonal basis corresponding to the input variable according to the distribution type corresponding to the input variable.
Step 104: and establishing a generalized chaotic polynomial model of the radiation sensitivity of the multi-conductor transmission line by combining a transmission line theory according to the orthogonal substrate of the input variable to obtain an expansion of the generalized chaotic polynomial.
Step 105: and performing hyperbolic truncation processing on the expanded form of the generalized chaotic polynomial by using a self-adaptive hyperbolic truncation method to obtain the hyperbolic truncated generalized chaotic polynomial.
Step 106: and processing the hyperbolic truncated generalized chaotic polynomial by using a minimum angle regression method to obtain a self-adaptive sparse chaotic polynomial.
Step 107: and calculating the induced current or induced voltage of the multi-conductor transmission line radiation sensitivity by using the self-adaptive sparse chaotic polynomial to obtain the probability distribution of the induced current or the probability distribution of the induced voltage at different frequency points.
Step 108: converting the self-adaptive sparse chaotic polynomial into an incremental summation form, and analyzing the global sensitivity index of the input variable by using a global sensitivity analysis method; the global sensitivity indicators include: total sensitivity index and first order sensitivity index.
The scheme integrally combines an adaptive hyperbolic truncation method and a minimum angle regression method, sparsely processes a generalized chaotic polynomial, establishes a proxy model to analyze the uncertainty problem of the radiation sensitivity of the multi-conductor transmission line, and combines a Sobol method to analyze the global sensitivity of an input variable.
The positional relationship between the multi-conductor transmission line and the incident field is shown in fig. 2, and since the uncertainty in the multi-conductor radiation sensitivity uncertainty problem is concentrated in the incident field, the present invention determines the input variables as the elevation angle θ, the azimuth angle ψ, the polarization angle η, and the level amplitude E of the incident plane wave in the implementation process0(ii) a L is the length of the transmission line, R is the termination impedance of the transmission line, and the number of transmission lines is S.
And (3) taking the original model as Y (xi), and carrying out expansion on the original model by using a generalized chaotic polynomial, wherein the expansion result is as follows:
wherein the content of the first and second substances,mixed orthogonal polynomials of order n, which are multidimensional standard random variablesInfinity is the highest power of the mixed orthogonal polynomial,and phiiAre respectively connected withAndin response to this, the mobile terminal is allowed to,coefficient of expansion being said generalized chaotic polynomial, ΦiAnd (xi) is the product of the one-dimensional orthogonal polynomial basis functions corresponding to the random variables.
For the generalized chaotic polynomial, the following describes the processing procedure of the adaptive hyperbolic truncation method and the minimum angle regression method.
The method for performing hyperbolic truncation processing on the expanded form of the generalized chaotic polynomial by using the self-adaptive hyperbolic truncation method specifically comprises the following steps:
obtaining an expansion of the generalized chaotic polynomial;
setting the highest truncation order P of the self-adaptive hyperbolic truncation method, and then performing primary sparse processing on the expansion of the generalized chaotic polynomial to obtain the primarily sparse generalized chaotic polynomial; the ith polynomial maximum truncation order piSatisfies the following conditions:where k is the dimension of the random variable, lkAnd the order of the random variable in the k dimension is shown, q is a norm, and n is the maximum value of the dimension of the random variable.
According to the highest truncation order piSatisfying the formula, it can be seen that when q is 1, pmaxThe truncation scheme is the conventional truncation scheme, i.e. PMeanwhile, the hyperbolic truncation scheme is to perform sparse processing on the high-order effect of the model on the basis of the traditional truncation scheme, and when q is used<1, the reserved polynomials are all positioned under a hyperbolic curve or curved surface. As shown in fig. 3, when n is 2, the norm q is 1, 0.75, and 0.5, and the hyperbolic truncation diagram corresponds to a schematic diagram, it can be intuitively seen that, for input variables with different dimensions, the hyperbolic truncation scheme can effectively reduce the influence of the model high order effect and keep the low order effect below the curve, and as the q value decreases, the penalty of the hyperbolic truncation scheme on the model high order effect is more obvious.
Fig. 4 is a truncation effect of a conventional truncation order P of 5 and an input variable dimension n of 2 hyperbolic truncation, and in order to show the influence of a related sparse processing method on a model, Q is a sparse coefficient: q is N '/N, where N' is a polynomial expansion term number after the sparse processing, N is a polynomial expansion term number without the sparse processing, and the smaller the sparse coefficient Q is, the lower the required computation cost is, for a model in which an input variable N is 4 and a hyperbolic truncation norm Q is 0.8, a polynomial term pair of a traditional truncation scheme polynomial term number and a hyperbolic truncation scheme polynomial term number with different truncation orders is as shown in table 1:
TABLE 1
As can be seen from table (1), the hyperbolic truncation scheme can effectively reduce the number of terms of the polynomial, and as the truncation order increases, the sparse coefficient also decreases, so that the hyperbolic truncation scheme can effectively perform preliminary sparse processing on the model.
The norm q is an important parameter introduced by the invention, and a specific determination method comprises the following steps:
setting the range of the preliminary norm q as [ a, b ] and the step length as s;
setting an error threshold;
the preliminary norm q is gradually increased by the step length s and is according to the formulaCalculating out a leave-one-out cross validation error;
stopping calculation when the cross validation error is smaller than the error threshold or exceeds the maximum value of the range of the preliminary norm q to obtain the norm q; wherein the content of the first and second substances,LOOto leave one out cross validation errors, M (x)(i)) For hyperbolic truncation model after setting norm q value at ith sample point x(i)Response value at Point, MPC\i(x(i)) Is a chaotic polynomial model after P-order truncation at the ith sample point x(i)The value of the response of (c) to (d),is the mean of the response values for the L sample points.
When the range of the hyperbolic truncation norm q is selected, as can be seen from fig. 3, along with the decrease of q, while punishing a high-order effect, the influence of hyperbolic truncation on a low-order effect is gradually increased, when the value of q is too small, the low-order effect of the model is punished by a hyperbolic curve, which may greatly affect the precision of the model, and if a small value is included in the selection range of the norm q, this may cause unnecessary iterative computation of the program, which may cause a waste of the computation cost, so generally when the range of q is selected, a and b are often set to 0.5 and 1, respectively, so that the norm q of the model may be computed more quickly, and the computation cost may be saved.
The chaos polynomial adopting the hyperbolic truncation scheme effectively reduces the influence of high-order effects among variables in the model, and as can be seen from table 1, compared with the generalized chaos polynomial, the polynomial term number after the hyperbolic truncation processing is reduced by more than half, but as the truncation order is increased, the polynomial term number after the hyperbolic truncation processing still has a certain sparse processing space due to the fact that the basic term number of the polynomial is very much. The minimum angle regression method is an important sparse algorithm and can effectively carry out sparse processing on the model, so that on the basis of the self-adaptive hyperbolic truncation scheme, in order to further sparsify the number of terms of the chaotic polynomial, the minimum angle regression method is adopted to carry out sparse processing on the hyperbolic truncated polynomial, the method integrates the advantages of the forward selection algorithm and the forward gradient algorithm, greedy iteration is not needed to be carried out like the forward selection algorithm, and fitting calculation can be completed only by fewer iteration times than the forward gradient algorithm.
Processing the hyperbolic truncated generalized chaotic polynomial by using a minimum angle regression method to obtain a self-adaptive sparse chaotic polynomial, which specifically comprises the following steps:
acquiring the hyperbolic truncated generalized chaotic polynomial;
converting the hyperbolic truncated generalized chaotic polynomial into a Y ═ theta X model; wherein, theta is a set of coefficients of the hyperbolic truncated generalized chaotic polynomial, and X is a set of the hyperbolic truncated generalized chaotic polynomial;
selecting independent variable X most relevant to dependent variable Y by cosine similarity methodj(ii) a Independent variable XjBelonging to a set of the hyperbolic truncated generalized chaotic polynomials;
calculating Y and XjResidual γ of from XjIs advanced in the direction of thetajAt this time, the residual γ is Y- θj*XjProceed until another variable X appearstThe degree of correlation with the residual y being equal to XjEqual to the residual gamma, i.e. the residual gamma is located at XjAnd XtThen continues forward along the bisector until the next argument appears that is the same as the residual correlation. In short, the variable most correlated to the current residual γ is continuously sought.
Fig. 5 is a schematic diagram of the principle of the minimum angle regression method using two-dimensional input variables as an example, which specifically includes:
1. finding the variable X most related to the dependent variable Y through cosine similarity1。
2. Let variable X1Move in the current direction until another variable X appears2So that Y-X1*θ1Can divide X equally1And X2I.e. X2Correlation with residual gamma and X1The residual gamma is the same as the correlation degree of the residual gamma, and the residual gamma is positioned at X1And X2On the bisector of (2), in X, as illustrated in FIG. 41Is moved by theta along the current direction1At time of multiplication, the residual γ at this time is Y- θ1*X1And gamma bisects X1And X2Angle of (e), at this time theta1I.e. variable X1The coefficient of (a).
3. Similarly to the second step, the process continues to proceed along the direction of the angular bisector obtained in the step 2 until a new independent variable X appears3The correlation with the residual γ is the same.
4. Updating the coefficients and moving eligible arguments from the candidate set to the active set.
5. Repeating the steps until all the variables are iterated completely, namely all important polynomials are selected, and the obtained theta is the coefficient of the polynomials. And substituting the coefficients of the polynomial into a self-adaptive sparse chaotic polynomial to obtain the self-adaptive sparse chaotic polynomial.
After an AS-PC agent model of the multi-conductor transmission line radiation sensitivity is established, the mean value and the variance of induced current or induced voltage are calculated, and the self-adaptive sparse chaotic polynomial is converted intoWhereinFor the sparsely processed polynomial coefficient, phii'is a polynomial after sparse processing, and p' is the number of terms of the polynomial after sparse processing;
obtaining the mean value of the multi-conductor transmission line radiation sensitivity induction current or the mean value of the multi-conductor transmission line radiation sensitivity induction voltage according to the formula after the hyperbolic truncation and the generalized chaotic polynomial conversion:wherein the content of the first and second substances,a constant term of an expansion of the truncated generalized chaotic polynomial;
obtaining the variance of the multi-conductor transmission line radiation sensitivity induction current or the variance of the multi-conductor transmission line radiation sensitivity induction voltage according to the formula after the hyperbolic truncation and the generalized chaotic polynomial conversion:
according to the self-adaptive sparse chaotic polynomial, the global sensitivity index of each input variable is calculated by combining a global sensitivity analysis method, and the influence degree of each input variable on the system is quantized, wherein the specific process comprises the following steps:
acquiring the self-adaptive sparse chaotic polynomial;
converting the adaptive sparse chaotic polynomial into an expansion in an incremental summation form;
carrying out recursive calculation on the expansion of the incremental summation form by using an integral method, solving the corresponding coefficient of the decomposition term, and obtaining the expansion of the incremental summation form with known coefficient;
taking variance from both sides of the expansion of the known incremental summation form of the coefficient at the same time to obtain a variance decomposition formula;
obtaining a global sensitivity index according to the variance decomposition formula; the global sensitivity index includes a sensitivity index and a total sensitivity index.
The expansion of the incremental summation form is:
wherein S isiIs an index of first-order sensitivity.
In order to more intuitively show the technical effect that the present invention can achieve, the following takes the number S of the transmission lines as 3 as an example, and further details the specific scheme of the present invention.
A 3+1 multi-conductor transmission line model with the ground as a reference conductor is established as shown in fig. 6, the length L of the transmission line is 1m, the radius r of the transmission line is 0.4mm, the distance between the transmission lines is 1cm, the height h from the ground is 2cm, and the loads of the source end and the load end of the transmission line are both 50 Ω.
Setting random input variables as an elevation angle theta, an azimuth angle psi, a polarization angle eta and a level amplitude E of incident plane waves, enabling the elevation angle theta to obey uniform distribution on an interval [0,0.5 pi ], enabling the azimuth angle psi to obey uniform distribution on the interval [ -pi, pi ], enabling the polarization angle eta to obey uniform distribution on the interval [0,2 pi ], enabling the level amplitude E to obey normal distribution with the mean value of 1V/m and the standard deviation of 0.2V/m, and selecting corresponding orthogonal bases according to different distributions obeyed by the random input variables to construct a generalized chaotic polynomial model.
Under the conventional truncation scheme, the current I is induced at the far end of the second transmission line in FIG. 6R2For example, let the truncation order P be 5, 10, 15, 20, respectively, the mean and standard deviation of the generalized chaotic polynomial model of the multi-conductor transmission line radiation sensitivity are calculated, and the frequency range of the incident field is set to [10MHz,1GHz []The calculated mean and standard deviation were compared with the calculation results of the 20000 MC methods, and the comparison results are shown in fig. 7.
From the comparison results of fig. 7, the induced current I is radiated sensitively to the multi-conductor transmission lineR2The mean and standard deviation of the low-order truncation is not good under the traditional truncation scheme, but the mean and standard deviation calculation result is closer to the calculation result of the MC as the truncation order increases. Then, for different truncation orders, three frequency points of 50MHz, 80MHz and 100MHz are selectedLOOFor comparison, as shown in fig. 8. As can be seen in FIG. 7, the low order truncations are, for example, 5 th and 10 th orderLOORelatively large, with calculation results not ideal, and with 15 th order truncationLOOWith truncation of order 20LOOLower and substantially consistent, it is proved that the accuracy of the calculation is higher as the truncation order P increases, which indicates that the 15 th order truncation can reach higher accuracy already in the case of the conventional truncation.
The expansion of the generalized chaotic polynomial in the 15-order conventional truncation scheme has 3876 terms, so that the generalized chaotic polynomial still has a large sparse processing space, and in order to further reduce the calculation cost and improve the calculation efficiency, the adaptive hyperbolic truncation scheme introduced in the foregoing is adopted to process the polynomial. At 50MHzInduced current I ofR2For example, let the selection interval of hyperbolic truncated norm q be [0.5, 1%]The step length s is 0.05,LOOthe threshold was 0.05, which was calculated to give the value shown in FIG. 9LOOThe error curve of (2). Table 2 provides sparse coefficients corresponding to different norms q when P ═ 15:
TABLE 2
q | 0.5 | 0.55 | 0.6 | 0.65 | 0.7 | 0.75 | 0.8 | 0.85 | 0.9 | 0.95 | 1 |
N’ | 263 | 358 | 480 | 662 | 873 | 1133 | 1476 | 1893 | 2404 | 3095 | 3876 |
Q | 6.79% | 9.24% | 12.4% | 17.1% | 22.5% | 29.2% | 38.1% | 48.8% | 62.0% | 79.9% | 1 |
As can be seen from fig. 9, when q is 0.8,LOOthe threshold is reached and it can be seen that when q is reached>At the time of 0.8, the temperature of the alloy is higher,LOOis small, and combines the sparse coefficient Q in the table with respect to Q>0.8 in the case ofLOOAlmost the same time, the sparse coefficient with q equal to 0.8 is the smallest, and the calculation cost is the lowest, so that the method with q equal to 0.8 meets the expectation of the invention, and can complete the preliminary sparse processing on the polynomial. The invention proceeds to further analytical studies with q-norm set to 0.8.
After the selection of the hyperbolic truncated norm q is completed, according to the introduction in the foregoing, the present invention combines the minimum angle regression method introduced in the foregoing to perform deeper sparse processing on the model. It has been demonstrated in the foregoing that under the conventional truncation scheme, better accuracy can be achieved when P-15, in order to verify that dual truncation is usedWhen the curved truncation is combined with the minimum angle regression method, the ideal precision can be achieved by taking the induced currents at 50MHz, 80MHz and 100MHz as examples, and the induced currents with different truncation orders are usedLOOFor comparison, as shown in fig. 10. It can be seen that when q is 0.8, P is 15 and P is 20 compared to the calculation of the low order truncationLOOThe variation is small and basically consistent, and the model adopted by the invention can achieve ideal precision when q is 0.8 and P is 15. In order to further verify the accuracy of the combination of the hyperbolic truncation and the minimum angle regression method, the invention calculates the probability distribution of the radiation sensitivity induced current of the multi-conductor transmission line at the frequency points of 50Mhz, 80Mhz and 100Mhz, and calculates the probability distribution curve of the radiation sensitivity induced current of the multi-conductor transmission line by combining the hyperbolic truncation and the minimum angle regression method, as shown in fig. 11, from the result of the probability distribution curve in the graph, it can be seen that the method of combining the hyperbolic truncation and the minimum angle regression method can effectively calculate the induced current I at different frequency points on the premise of ensuring the calculation accuracyR2The calculation time of the method combining the curved truncation and the minimum angle regression method, the calculation time of the generalized chaotic polynomial and the sparse coefficient Q are further analyzed, as shown in table 3:
TABLE 3
Calculating time/s | Sparse coefficient Q | |
AS-PC | 6.39 | 3.25% |
Generalized chaotic polynomial | 13.69 | 100% |
20000 times of MC | 2001.54 | nan |
The CPU main frequency of the computer used by the invention is 2.3GHz, the operating memory is 8GB, and table 3 shows that compared with the computation time of 20000 Monte Carlo methods, the computation time of the generalized chaotic polynomial is greatly shortened, but the hyperbolic truncation and the minimum angle regression method are combined to further effectively compress the computation time, so that the computation cost is saved while the computation precision is ensured, and the computation efficiency is improved. The calculation results are combined, and the hyperbolic truncation and minimum angle regression method adopted by the method can effectively make the generalized chaotic polynomial sparse and accurately calculate the probability distribution of the multi-conductor transmission line radiation sensitivity induced current at different frequency points. In practical applications, the induced current of the radiation sensitivity of the multi-conductor transmission line is limited in the interval [ mu-3 sigma, mu +3 sigma]Can be used as an induced current IR3Approximate upper and lower limits, with high quantile confidence (e.g., 0.95, 0.99, 0.995, etc.) may be used to demonstrate that μ -3 σ may be effective as the sense current IR2But this is not within the scope of the present invention's study, which was validated in conjunction with the calculation of 20000 monte carlo method, as shown in fig. 12: the upper limit mu +3 sigma obtained by hyperbolic truncation and minimum angle regression method calculation well covers most Monte Carlo simulation curves, and proves that the mu +3 sigma can be effectively used as the upper limit of the variation range of the induced current.
Through the calculation and analysis, the hyperbolic truncation and minimum angle regression method adopted by the invention can quickly and accurately calculate related statistical characteristic parameters (such as mean standard deviation probability distribution) and the like in the uncertainty problem of the radiation sensitivity of the multi-conductor transmission line on the premise of ensuring the calculation precision. Next, in order to further analyze the influence degree of different input variables on the model in the multi-conductor transmission line radiosensitivity system, the present invention will be analyzed and studied in combination with the global sensitivity calculation method described in the foregoing.
The response current of the right end of the No. 2 transmission line at 50MHz is taken as an index, the total sensitivity index and the first-order sensitivity index comparison graph of each input variable are obtained by calculation, as shown in fig. 13 and fig. 14, and it can be seen from the two graphs that the total sensitivity index and the first-order sensitivity index calculated based on the method of the invention are basically consistent with the result calculated by 20000 times of MC method, the influence degree of each different input variable on the model is also consistent, the total sensitivity index and the first-order sensitivity index of the polarization angle eta are both kept at a higher level, which are important factors influencing the radiation sensitivity of the multi-conductor transmission line at the frequency point, and the sensitivity index of the elevation angle theta is lower, and the whole influence degree on the model is not large.
By combining the analysis, the method adopted by the invention has higher calculation speed and higher efficiency, and proves that the method adopted by the invention is effective in calculating the total sensitivity index and the first-order sensitivity index of the radiation sensitivity of the multi-conductor transmission line.
In order to more intuitively represent the influence degrees of different random variables at different frequency points on the whole model, the total sensitivity index of the influence degrees of each parameter on the multi-conductor transmission radiation sensitivity at the frequency band [10MHz,1GHz ] is calculated, as shown in fig. 15, it can be seen that the influence degree of the level amplitude E is kept at a lower level on the frequency band [10MHz,1GHz ], and the influence on the whole model is not large, although the influence at the higher frequency band, such as [900MHz,1GHz ] is increased, the influence degree of the polarization angle η at the same-frequency band is far greater than the level amplitude E, when the frequency is higher than 200MHz, the influence degree of the elevation angle θ on the model is also increased, the azimuth angle ψ becomes an important factor influencing the model at about 200MHz, but the influence degree at the high-frequency band is greatly reduced.
By combining the above analysis, in practical engineering application, when the position and frequency range of the radiation source in the surrounding environment are known, the polarization angle of the radiation source received by the system should be considered when designing an electrical system in the environment, so as to reasonably and effectively adjust the position distribution of the transmission line and avoid the occurrence of unnecessary electromagnetic compatibility problems.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (9)
1. A method for field line coupling uncertainty quantification and global sensitivity calculation, the method comprising:
determining an input variable according to an electromagnetic environment and an incident field where the multi-conductor transmission line is located;
determining a distribution type corresponding to the input variable according to the input variable;
determining an orthogonal base corresponding to the input variable according to the distribution type corresponding to the input variable;
establishing a generalized chaotic polynomial model of the radiation sensitivity of the multi-conductor transmission line by combining a transmission line theory according to the orthogonal substrate of the input variable to obtain an expansion of the generalized chaotic polynomial;
performing hyperbolic truncation processing on the expanded form of the generalized chaotic polynomial by using a self-adaptive hyperbolic truncation method to obtain the hyperbolic truncated generalized chaotic polynomial;
processing the hyperbolic truncated generalized chaotic polynomial by using a minimum angle regression method to obtain a self-adaptive sparse chaotic polynomial;
calculating induced current or induced voltage of the multi-conductor transmission line radiation sensitivity by using the self-adaptive sparse chaotic polynomial to obtain probability distribution of the induced current or the induced voltage at different frequency points;
converting the self-adaptive sparse chaotic polynomial into an incremental summation form, and analyzing the global sensitivity index of the input variable by using a global sensitivity analysis method; the global sensitivity indicators include: total sensitivity index and first order sensitivity index.
2. The field line coupling uncertainty quantization and global sensitivity calculation method of claim 1, wherein the expansion of the generalized chaotic polynomial is:
wherein the content of the first and second substances,mixed orthogonal polynomials of order n, which are multidimensional standard random variablesInfinity is the highest power of the mixed orthogonal polynomial,and phiiAre respectively connected withAndin response to this, the mobile terminal is allowed to,coefficient of expansion being said generalized chaotic polynomial, ΦiAnd (xi) is the product of the one-dimensional orthogonal polynomial basis functions corresponding to the random variables.
3. The field line coupling uncertainty quantification and global sensitivity calculation method according to claim 1, wherein hyperbolic truncation processing is performed on the expansion of the generalized chaotic polynomial by using a self-adaptive hyperbolic truncation method to obtain a hyperbolic truncated generalized chaotic polynomial, and specifically the method comprises:
obtaining an expansion of the generalized chaotic polynomial;
setting the highest truncation order P of the self-adaptive hyperbolic truncation method, and then performing primary sparse processing on the expansion of the generalized chaotic polynomial to obtain the primarily sparse generalized chaotic polynomial; the ith polynomial maximum truncation order piSatisfies the following conditions:q is more than 0 and less than or equal to 1, wherein k is the dimension of a random variable, lkAnd the order of the random variable in the k dimension is shown, q is a norm, and n is the maximum value of the dimension of the random variable.
4. The method for field line coupling uncertainty quantification and global sensitivity calculation of claim 3, wherein the method for determining the norm q comprises:
setting the range of the preliminary norm q as [ a, b ] and the step length as s;
setting an error threshold;
the preliminary norm q is gradually increased by the step length s and is according to the formulaCalculating out a leave-one-out cross validation error;
when the cross-validation error is less than the error threshold orStopping calculation when the maximum value of the range of the preliminary norm q is exceeded, and obtaining the norm q; wherein the content of the first and second substances,LOOto leave one out cross validation errors, M (x)(i)) For hyperbolic truncation model after setting norm q value at ith sample point x(i)Response value at Point, MPC\i(x(i)) Is a chaotic polynomial model after P-order truncation at the ith sample point x(i)The value of the response of (c) to (d),is the mean of the response values for the L sample points.
5. The field line coupling uncertainty quantification and global sensitivity calculation method according to claim 1, wherein the hyperbolic truncated generalized chaotic polynomial is processed by using a minimum angle regression method to obtain an adaptive sparse chaotic polynomial, specifically comprising:
acquiring the hyperbolic truncated generalized chaotic polynomial;
converting the hyperbolic truncated generalized chaotic polynomial into a Y ═ theta X model; wherein, theta is a set of coefficients of the hyperbolic truncated generalized chaotic polynomial, and X is a set of the hyperbolic truncated generalized chaotic polynomial;
processing the Y-theta X model by using a minimum angle regression method to obtain a coefficient of a polynomial;
and substituting the coefficients of the polynomial into a self-adaptive sparse chaotic polynomial to obtain the self-adaptive sparse chaotic polynomial.
6. The field line coupling uncertainty quantification and global sensitivity calculation method according to claim 1, wherein the calculating of the induced current or induced voltage of the multi-conductor transmission line radiation sensitivity using the adaptive sparse chaotic polynomial to obtain the probability distribution of the induced current or the probability distribution of the induced voltage at different frequency points specifically comprises:
acquiring the self-adaptive sparse chaotic polynomial;
will be self-supportingAdaptive sparse chaotic polynomial conversionWhereinIs a sparse processed polynomial coefficient of phi'iIs a polynomial after sparse processing, and p' is the number of terms of the polynomial after sparse processing;
obtaining the mean value of the multi-conductor transmission line radiation sensitivity induction current or the mean value of the multi-conductor transmission line radiation sensitivity induction voltage according to the formula after the hyperbolic truncation and the generalized chaotic polynomial conversion:wherein the content of the first and second substances,a constant term of an expansion of the truncated generalized chaotic polynomial;
obtaining the variance of the multi-conductor transmission line radiation sensitivity induction current or the variance of the multi-conductor transmission line radiation sensitivity induction voltage according to the formula after the hyperbolic truncation and the generalized chaotic polynomial conversion:
7. the field line coupling uncertainty quantification and global sensitivity calculation method according to claim 1, wherein the transforming the adaptive sparse chaotic polynomial into an incremental summation form, and analyzing the global sensitivity index of the input variable by using a global sensitivity analysis method specifically comprises:
acquiring the self-adaptive sparse chaotic polynomial;
converting the adaptive sparse chaotic polynomial into an expansion in an incremental summation form;
carrying out recursive calculation on the expansion of the incremental summation form by using an integral method, solving the corresponding coefficient of the decomposition term, and obtaining the expansion of the incremental summation form with known coefficient;
taking variance from both sides of the expansion of the known incremental summation form of the coefficient at the same time to obtain a variance decomposition formula;
obtaining a global sensitivity index according to the variance decomposition formula; the global sensitivity index includes a sensitivity index and a total sensitivity index.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010696773.3A CN111797541B (en) | 2020-07-20 | 2020-07-20 | Method for quantifying coupling uncertainty and calculating global sensitivity of field line |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010696773.3A CN111797541B (en) | 2020-07-20 | 2020-07-20 | Method for quantifying coupling uncertainty and calculating global sensitivity of field line |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111797541A true CN111797541A (en) | 2020-10-20 |
CN111797541B CN111797541B (en) | 2022-08-26 |
Family
ID=72807902
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010696773.3A Active CN111797541B (en) | 2020-07-20 | 2020-07-20 | Method for quantifying coupling uncertainty and calculating global sensitivity of field line |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111797541B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113312841A (en) * | 2021-05-28 | 2021-08-27 | 合肥工业大学 | Variable norm equivalent source near-field acoustic holography algorithm with acoustic source sparsity adaptivity |
CN114019250A (en) * | 2021-10-31 | 2022-02-08 | 吉林大学 | Random non-uniform transmission line electromagnetic radiation analysis method |
CN115048809A (en) * | 2022-07-11 | 2022-09-13 | 中国核动力研究设计院 | Uncertainty quantification method for multi-physics multi-scale coupling system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106155980A (en) * | 2015-04-01 | 2016-11-23 | 国家电网公司 | Overhead transmission line field wire coupling modeling method based on inhomogeneous excitation electric field |
CN108416148A (en) * | 2018-03-13 | 2018-08-17 | 西北核技术研究所 | A kind of high-altitude electromagnetic pulse field wire coupling uncertainty acquisition methods based on polynomial chaos expression |
CN110031406A (en) * | 2019-03-21 | 2019-07-19 | 北京航空航天大学 | A kind of atmosphere-ocean coupled system optics shallow water remote sensing signal Sensitivity Analysis |
CN110188963A (en) * | 2019-06-05 | 2019-08-30 | 中国矿业大学(北京) | Return the building energy consumption prediction technique with gradient lift scheme based on minimum angular convolution |
CN111141974A (en) * | 2019-12-26 | 2020-05-12 | 吉林大学 | Multi-conductor transmission line radiation sensitivity analysis method |
-
2020
- 2020-07-20 CN CN202010696773.3A patent/CN111797541B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106155980A (en) * | 2015-04-01 | 2016-11-23 | 国家电网公司 | Overhead transmission line field wire coupling modeling method based on inhomogeneous excitation electric field |
CN108416148A (en) * | 2018-03-13 | 2018-08-17 | 西北核技术研究所 | A kind of high-altitude electromagnetic pulse field wire coupling uncertainty acquisition methods based on polynomial chaos expression |
CN110031406A (en) * | 2019-03-21 | 2019-07-19 | 北京航空航天大学 | A kind of atmosphere-ocean coupled system optics shallow water remote sensing signal Sensitivity Analysis |
CN110188963A (en) * | 2019-06-05 | 2019-08-30 | 中国矿业大学(北京) | Return the building energy consumption prediction technique with gradient lift scheme based on minimum angular convolution |
CN111141974A (en) * | 2019-12-26 | 2020-05-12 | 吉林大学 | Multi-conductor transmission line radiation sensitivity analysis method |
Non-Patent Citations (4)
Title |
---|
MOURAD LARBI .ETC: "Variability Impact of Many Design Parameters : The Case of a Realistic Electronic Link", 《IEEE TRANSACTIONS ON ELECTROMAGNETIC COMPATIBILITY》 * |
于全毅等: "基于广义混沌多项式法的多道题传输线辐射敏感度分析方法", 《电工技术学报》 * |
王天皓等: "基于混沌多项式展开式的线束串扰统计模型", 《吉林大学学报(工学版)》 * |
胡军等: "基于多项式混沌的全局敏感度分析", 《计算物理》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113312841A (en) * | 2021-05-28 | 2021-08-27 | 合肥工业大学 | Variable norm equivalent source near-field acoustic holography algorithm with acoustic source sparsity adaptivity |
CN113312841B (en) * | 2021-05-28 | 2022-09-13 | 合肥工业大学 | Variable norm equivalent source near-field acoustic holography algorithm with acoustic source sparsity adaptivity |
CN114019250A (en) * | 2021-10-31 | 2022-02-08 | 吉林大学 | Random non-uniform transmission line electromagnetic radiation analysis method |
CN114019250B (en) * | 2021-10-31 | 2023-09-29 | 吉林大学 | Random non-uniform transmission line electromagnetic radiation analysis method |
CN115048809A (en) * | 2022-07-11 | 2022-09-13 | 中国核动力研究设计院 | Uncertainty quantification method for multi-physics multi-scale coupling system |
CN115048809B (en) * | 2022-07-11 | 2023-08-29 | 中国核动力研究设计院 | Uncertainty quantification method for multi-physical multi-scale coupling system |
Also Published As
Publication number | Publication date |
---|---|
CN111797541B (en) | 2022-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111797541B (en) | Method for quantifying coupling uncertainty and calculating global sensitivity of field line | |
CN111209347B (en) | Method and device for clustering mixed attribute data | |
Xu et al. | Asymptotic properties of Pearson׳ s rank-variate correlation coefficient in bivariate normal model | |
Akrou et al. | Enhanced modal tracking for characteristic modes | |
CN111141974A (en) | Multi-conductor transmission line radiation sensitivity analysis method | |
CN115480206A (en) | Off-grid DOA estimation method | |
CN114520755A (en) | Improved characteristic value-to-angle loading information source number estimation method | |
Forrester et al. | A Fuchsian matrix differential equation for Selberg correlation integrals | |
CN115906741A (en) | Radio frequency circuit optimization design method based on high-performance calculation | |
CN115577779A (en) | Bare state information determination method and device for multi-body system in superconducting quantum chip layout | |
CN115343673A (en) | Broadband signal DOA estimation method based on Dirichlet process prior | |
CN115577778A (en) | Method and device for determining equivalent coupling strength between quantum devices in superconducting quantum chip layout | |
Yucel | Uncertainty quantification for electromagnetic analysis via efficient collocation methods | |
Seo | An IE-FFT algorithm to analyze PEC objects for MFIE formulation | |
Ewe et al. | Solving mixed dielectric/conducting scattering problem using adaptive integral method | |
Wang et al. | Radiated Susceptibility Analysis of Multiconductor Transmission Lines Based on Polynomial Chaos. | |
CN114417585A (en) | Antenna broadband method, system, equipment and medium | |
Kirichenko | Radiation from a plasma layer with a strong longitudinal irregularity | |
CN115659905B (en) | Method and device for determining coupling strength between quantum devices in superconducting quantum chip layout | |
Meng et al. | A wideband 2-d fast multipole algorithm with a novel diagonalization form | |
Moroz et al. | The CORDIC method of calculating the exponential function | |
Guo et al. | Synthesis of Linear Antenna Array for Wireless Power Transmission | |
Chen et al. | An Improved Algorithm for High-Precision Frequency Estimation | |
Chen et al. | A novel PO solver for uncertainty EM computation of electrically large targets | |
Saçın et al. | Neural network modeling of antennas on package for 5G applications |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |