Particle energy spectrum inversion method and system based on double regular coefficients

Abstract

The invention discloses a particle energy spectrum inversion method and a particle energy spectrum inversion system based on double regular coefficients, so as to achieve the purpose of one-time inversion of electron and proton energy spectrums by using fewer detection devices. The detector is a single detector and detects electrons and protons at the same time, and the method comprises the following steps: generating a first inversion model based on the detector; wherein the first inversion model comprises a local relaxation coefficient and a biregular coefficient, and the biregular coefficient comprises a first regularization coefficient corresponding to electrons and a second regularization coefficient corresponding to protons; presetting a local relaxation coefficient set of local relaxation coefficients and a biregular coefficient set of biregular coefficients, performing iterative calculation on the first inversion model by adopting a Boosted-Gold algorithm based on the biregular coefficients on the basis of the local relaxation coefficient set and the biregular coefficient set, and recording an energy spectrum norm and a residual norm which meet a preset convergence condition; and generating an L curve based on the spectral norm and the residual norm, and taking an inflection point of the L curve as an optimal solution for inversion.

Description

Technical Field

[0001] The present application relates to the technical field of space high-energy particle detection, and in particular to a particle energy spectrum inversion method and system based on double canonical coefficients. Background Art

[0002] Charged particle radiation in the space environment is mainly composed of electrons and protons. Both types of particles use the mechanism of ionization excitation in the material after the particles lose energy when interacting with the sensor material after incidence. Without taking special measures, both electrons and protons will respond in the proton-electron load, that is, there is contamination between each other. The inversion of high-energy electron energy spectrum faces multiple technical difficulties. First, the interference between electron and proton signals is one of the main challenges. The proton signal may be confused with the electron signal in the detector, affecting the inversion accuracy. Secondly, the energy deposition of high-energy electrons of a single energy in the detector has the uncertainty of Gaussian distribution, which leads to the overlap of the deposition energy spectra of electrons of different energies, increasing the complexity of the inversion.

[0003] Due to the limitation of accelerator beam conditions, it is impossible to collect data covering all energy ranges, and GEANT4 Monte Carlo simulation is required. In the prior art, the basis of particle inversion based on Monte Carlo simulation is to use different detectors to collect electron and proton signals respectively and then perform inversion separately. Summary of the invention

[0004] The purpose of the present application is to provide a particle energy spectrum inversion method and system based on double canonical coefficients, so as to achieve the purpose of inverting electron and proton energy spectra at one time by using fewer detection devices.

[0005] In a first aspect, the present application provides a particle energy spectrum inversion method based on dual canonical coefficients, wherein the detector is a single detector, and the detector detects electrons and protons simultaneously, and the method comprises: Based on the detector, a first inversion model is generated; wherein the first inversion model includes a local relaxation coefficient and a dual regularization coefficient, and the dual regularization coefficient includes a first regularization coefficient corresponding to electrons and a second regularization coefficient corresponding to protons; Preset a local relaxation coefficient set of the local relaxation coefficient and a dual regularization coefficient set of the dual regularization coefficient, and based on the local relaxation coefficient set and the dual regularization coefficient set, use a Boosted-Gold algorithm based on the dual regularization coefficient to iteratively calculate the first inversion model, and record an energy spectrum norm and a residual norm that meet the preset convergence conditions; An L curve is generated based on the spectral norm and the residual norm, and an inflection point of the L curve is used as an inversion optimal solution.

[0006] In a possible embodiment, the iterative calculation of the first inversion model based on the local relaxation coefficient set and the dual regularization coefficient set using the Boosted-Gold algorithm with dual regularization coefficients includes: Traversing the set of biregular coefficients; When traversing the set of dual regular coefficients, each time the current dual regular coefficient is selected, the local relaxation coefficient is iteratively calculated; Iteratively updating the particle energy spectrum in the first inversion model according to a preset preliminary solution based on the current dual regularization coefficient and the current local relaxation coefficient calculated in each iteration; Repeat the previous step until the energy spectrum norm and the residual norm that meet the preset convergence condition are obtained when the current regularization coefficient is selected, and record the energy spectrum norm and the residual norm that meet the preset convergence condition when the current regularization coefficient is selected.

[0007] In a possible embodiment, generating the first inversion model includes: Obtaining an electron response function and a proton response function of the detector; determining a response matrix based on the electron response function and the proton response function; Determining the basic inversion model based on the response matrix and the detection value of the detector; Performing transformation based on the basic inversion model to generate the first inversion model; The preset preliminary solution of the first inversion model is determined.

[0008] In a possible embodiment, the response matrix is ​​as follows: , (1) in, To measure the total number of energy channels, is the total number of electron incident energy channels, is the total number of proton incident energy channels; For electrons in the channel The response function of For protons in the channel The response function of is the number of the measurement energy level, is the number of the incident energy channel, The value range is 1 to The natural number of The value range is 1 to The natural number of The basic inversion model is as follows: (2) in, is the response matrix, is the electron spectrum or proton spectrum, is the measurement error, is the measured value.

[0009] In a possible embodiment, the first inversion model is formulated as follows: (3) Among them, A T is the inverse matrix of A; μ is the local relaxation coefficient, which is used to adjust the update amplitude of the solution. Each element is ; is the first regularization coefficient, is the second regularization coefficient; is the particle energy spectrum, is the electron energy spectrum, Proton spectrum.

[0010] In a second aspect, the present application provides a particle energy spectrum inversion system based on dual canonical coefficients, wherein the detector is a single detector, and the detector detects electrons and protons simultaneously, including: A model unit, configured to generate a first inversion model based on the detector; wherein the first inversion model comprises a local relaxation coefficient and a dual regularization coefficient, and the dual regularization coefficient comprises a first regularization coefficient corresponding to electrons and a second regularization coefficient corresponding to protons; A processing unit, used for presetting a local relaxation coefficient set of the local relaxation coefficient and a dual regularization coefficient set of the dual regularization coefficient, and based on the local relaxation coefficient set and the dual regularization coefficient set, using a Boosted-Gold algorithm based on the dual regularization coefficient to iteratively calculate the first inversion model, and recording an energy spectrum norm and a residual norm that meet a preset convergence condition; A result unit is used to generate an L curve based on the spectral norm and the residual norm, and take the inflection point of the L curve as the best inversion solution.

[0011] In a possible embodiment, the processing unit performs iterative calculation on the first inversion model, including: Traversing the set of biregular coefficients; When traversing the set of dual regular coefficients, each time the current dual regular coefficient is selected, the local relaxation coefficient is iteratively calculated; Iteratively updating the particle energy spectrum in the first inversion model according to a preset preliminary solution based on the current dual regularization coefficient and the current local relaxation coefficient calculated in each iteration; Repeat the previous step until the energy spectrum norm and the residual norm that meet the preset convergence condition are obtained when the current regularization coefficient is selected, and record the energy spectrum norm and the residual norm that meet the preset convergence condition when the current regularization coefficient is selected.

[0012] In a possible embodiment, the model unit generates the first inversion model, including: Obtaining an electron response function and a proton response function of the detector; determining a response matrix based on the electron response function and the proton response function; Determining the basic inversion model based on the response matrix and the detection value of the detector; Performing transformation based on the basic inversion model to generate the first inversion model; The preset preliminary solution of the first inversion model is determined.

[0013] In a possible embodiment, the response matrix is ​​as follows: , (1) in, To measure the total number of energy channels, is the total number of electron incident energy channels, is the total number of proton incident energy channels; For electrons in the channel The response function of For protons in the channel The response function of is the number of the measurement energy level, is the number of the incident energy channel, The value range is 1 to The natural number of The value range is 1 to The natural number of The basic inversion model is as follows: (2) in, is the response matrix, is the electron spectrum or proton spectrum, is the measurement error, is the measured value.

[0014] In a possible embodiment, the first inversion model is formulated as follows: (3) Among them, A T is the inverse matrix of A; μ is the local relaxation coefficient, which is used to adjust the update amplitude of the solution. Each element is ; is the first regularization coefficient, is the second regularization coefficient; is the particle energy spectrum, is the electron energy spectrum, Proton spectrum.

[0015] The particle energy spectrum inversion method and system based on double regularization coefficients provided in the embodiments of the present application use a single detector in combination with the particle energy spectrum inversion method based on double regularization coefficients, thereby reducing the number of detectors and being able to invert electron and proton energy spectra at one time, thereby improving the inversion efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 A flow chart of a particle energy spectrum inversion method based on double canonical coefficients provided in an embodiment of the present application; Figure 2 A flowchart of iterative calculation based on dual regular coefficients provided in an embodiment of the present application; Figure 3 A flowchart for generating a first inversion model provided in an embodiment of the present application; Figure 4 A schematic diagram showing a mapping relationship between the incident energy spectrum and the measured energy spectrum of protons and electrons in a detector provided in an embodiment of the present application; Figure 5 A schematic diagram of an L curve provided in an embodiment of the present application; Figure 6 A more specific flow chart of a particle energy spectrum inversion method based on double canonical coefficients provided in an embodiment of the present application; Figure 7 A structural diagram of a particle energy spectrum inversion system based on double regularization coefficients provided in an embodiment of the present application. DETAILED DESCRIPTION

[0017] In order to make the purpose, technical solution and advantages of the embodiments of the present application clearer, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of this application.

[0018] like Figure 1 As shown, an embodiment of the present application provides a particle energy spectrum inversion method based on dual canonical coefficients, wherein the detector is a single detector, and the detector detects electrons and protons simultaneously, and the method includes: 101. Generate a first inversion model based on the detector; wherein the first inversion model includes a local relaxation coefficient and a dual regularization coefficient, and the dual regularization coefficient includes a first regularization coefficient corresponding to electrons and a second regularization coefficient corresponding to protons.

[0019] 102. A local relaxation coefficient set of a preset local relaxation coefficient and a dual regularization coefficient set of a dual regularization coefficient are preset. Based on the local relaxation coefficient set and the dual regularization coefficient set, a Boosted-Gold algorithm based on the dual regularization coefficient is used to iteratively calculate the first inversion model, and the energy spectrum norm and the residual norm that meet the preset convergence conditions are recorded.

[0020] 103. Generate an L curve based on the spectral norm and the residual norm, and take the inflection point of the L curve as the best inversion solution.

[0021] In a possible embodiment, Figure 2 As shown, in step 102, based on the local relaxation coefficient set and the dual regularization coefficient set, the Boosted-Gold algorithm with the dual regularization coefficient is used to iteratively calculate the first inversion model, including: 201. Traverse the set of double regular coefficients.

[0022] 202. When traversing the set of double regular coefficients, each time the current double regular coefficient is selected, the local relaxation coefficient is iteratively calculated.

[0023] 203. Based on the current bi-regular coefficient and the current local relaxation coefficient calculated in each iteration, the particle energy spectrum in the first inversion model is iteratively updated according to a preset preliminary solution.

[0024] 204. Repeat the previous step 203 until the energy spectrum norm and residual norm that meet the preset convergence condition are obtained when the current regularization coefficient is selected, and record the energy spectrum norm and residual norm that meet the preset convergence condition when the current regularization coefficient is selected.

[0025] In a possible embodiment, Figure 3 As shown, the first inversion model is generated, including: 301. Obtaining an electron response function and a proton response function of the detector; 302. Determine a response matrix based on the electron response function and the proton response function; 303. Determine a basic inversion model based on the response matrix and the detection value of the detector; 304. Perform transformation based on the basic inversion model to generate a first inversion model; 304. Determine a preset preliminary solution of the first inversion model.

[0026] In a possible embodiment, the response matrix is ​​as follows: , (1) in, To measure the total number of energy channels, is the total number of electron incident energy channels, is the total number of proton incident energy channels; For electrons in the channel The response function of For protons in the channel The response function of . is the number of the measurement energy level, is the number of the incident energy channel, The value range is 1 to The natural number of The value range is 1 to The natural number of The basic inversion model is as follows: (2) in, is the response matrix, is the electron spectrum or proton spectrum, is the measurement error, is the measured value.

[0027] In a possible embodiment, the first inversion model has the following formula: (3) Among them, A T is the inverse matrix of A; μ is the local relaxation coefficient, which is used to adjust the update amplitude of the solution. Each element is ; is the first regularization coefficient, is the second regularization coefficient; is the particle energy spectrum, is the electron energy spectrum, Proton spectrum.

[0028] It should be noted that is the total number of energy channels to be measured, and its value is , ,and Under certain circumstances, In step 203, the specific iteration is based on the initial solution (denoted as ) as the starting point, Iterate.

[0029] In order to better understand the particle energy spectrum inversion method based on double regular coefficients provided in this application, combined with Figures 4 to 6 The detailed instructions are as follows: like Figure 4 As shown, the proton spectrum It is incident on the high-energy electron detector, and the interfering proton response probability , which produces interfering proton counts ;Electron Spectroscopy It is incident on a high-energy electron detector, and according to the electron response probability (Probability of correct response + probability of incorrect response), which generates an electronic count . Instrument measurement results is the sum of interfering proton counts and electron counts, that is: = + ; Therefore, in high-energy electron detectors, there is a set of mapping relationships between the incident energy spectrum and the measured energy spectrum of protons and electrons.

[0030] In view of the response characteristics of protons and electrons in the instrument, a single instrument proton-electron matrix joint inversion method is used for proton-electron spectrum inversion, and a Boosted-Gold algorithm based on double regularization coefficients is proposed. In real situations, the equations can be simplified, as shown in Formula 2.

[0031] In formula 2 Used to adjust the update amplitude of the solution to avoid unstable convergence caused by excessive updates.

[0032] Based on considerations in practical applications, there may be the following two problems: First, the measured data is often affected by noise, which may cause the inversion result to overfit the noise; second, when using a single instrument to perform proton-electron spectrum analysis simultaneously, the equations for spectrum analysis are underdetermined, that is, the number of unknowns exceeds the number of equations. At this time, the solution to the equations is prone to numerical instability, resulting in unreliable results. Therefore, in order to obtain a stable solution, the problem is regularized, that is, an additional term is added to the problem to suppress overfitting and falling into the local optimal solution: . At the same time, considering that the solution of the system of equations It is composed of the energy spectrum of electrons and the energy spectrum of protons, so two regularization coefficients are introduced for the actual situation and The purpose is to more flexibly control the stability and smoothness of each part of the solution by acting on the electron energy spectrum and the proton energy spectrum respectively. The iterative formula obtained is the first inversion model, as shown in Formula 3: (3) Regularization coefficient corresponding to the expected result and The L-curve method is used to determine the regularization parameter. The L-curve method is a graphical method for selecting regularization parameters, which is widely used in dealing with underdetermined or ill-conditioned problems. The L-curve is an image drawn on a logarithmic scale. Figure 5 As shown in the figure, the horizontal axis is the norm of the solution and the vertical axis is the norm of the fitting error. At the "corner" of the L curve, the regularization parameter reaches a balance point. The regularization parameter corresponding to the corner point is usually considered to be the best because it achieves the best balance between the norm of the solution (representing the complexity or smoothness of the solution) and the fitting error.

[0033] Based on the above description, combined with the method steps, Figure 6 A flowchart showing a more specific inversion process includes: 601. Obtain an electron response function and a proton response function of a single detector.

[0034] 602. Determine a response matrix.

[0035] The method comprises determining a response matrix based on the electron response function and the proton response function. In a possible embodiment, the response matrix is ​​as shown in Formula 1.

[0036] 603. Determine the basic inversion model.

[0037] The basic inversion model is determined based on the response matrix and the detection value of the detector. In a possible embodiment, the basic inversion model is as shown in Formula 2.

[0038] 604. Generate a first inversion model and determine a preliminary solution.

[0039] The method comprises transforming the basic inversion model, generating a first inversion model and determining a preliminary solution, wherein the first inversion model includes a local relaxation coefficient. , regularization coefficient and .

[0040] 605. Traverse the set of double regular coefficients.

[0041] 606. Each time the current double regular coefficient is selected, the local relaxation coefficient is iteratively calculated. .

[0042] 607. Iteratively update the particle energy spectrum in the first inversion model.

[0043] Including based on the current bi-regular coefficient and and the current local relaxation coefficient calculated at each iteration , the particle energy spectrum in the first inversion model is iteratively updated according to the preset preliminary solution.

[0044] 608. Check whether the convergence meets the preset conditions.

[0045] If the convergence meets the preset conditions, the energy spectrum norm and residual norm that meet the preset convergence conditions when the current regularization coefficient is recorded; if the convergence does not meet the preset conditions, the iterative calculation of the local relaxation coefficient in step 606 is performed. , until the current regularization coefficient is selected, the energy spectrum norm and residual norm that meet the preset convergence conditions are obtained, and the energy spectrum norm and residual norm that meet the preset convergence conditions when the current regularization coefficient is selected are recorded.

[0046] 609. Check whether the traversal of the double regular coefficient set is completed.

[0047] If yes, proceed to step 610 ; if no, proceed to step 606 .

[0048] 610. Generate an L curve and analyze the L curve to determine the best inversion solution. norm and residual norm to generate the L curve.

[0049] 611. Output energy spectrum.

[0050] like Figure 7 As shown, the present application provides a particle energy spectrum inversion system based on dual canonical coefficients, wherein the detector is a single detector, and the detector detects electrons and protons simultaneously, including: A model unit 701 is used to generate a first inversion model based on the detector; wherein the first inversion model includes a local relaxation coefficient and a double regularization coefficient, and the double regularization coefficient includes a first regularization coefficient corresponding to electrons and a second regularization coefficient corresponding to protons; The processing unit 702 is used to preset a local relaxation coefficient set of local relaxation coefficients and a dual regularization coefficient set of dual regularization coefficients, and based on the local relaxation coefficient set and the dual regularization coefficient set, use the Boosted-Gold algorithm based on the dual regularization coefficients to iteratively calculate the first inversion model, and record the energy spectrum norm and residual norm that meet the preset convergence conditions; The result unit 703 is used to generate an L curve based on the spectral norm and the residual norm, and take the inflection point of the L curve as the best inversion solution.

[0051] In a possible embodiment, the processing unit 702 performs iterative calculation on the first inversion model, including: Traverse the set of biregular coefficients; When traversing the set of double regular coefficients, each time the current double regular coefficient is selected, the local relaxation coefficient is iteratively calculated; Iteratively updating the particle energy spectrum in the first inversion model according to a preset preliminary solution based on the current bi-regular coefficient and the current local relaxation coefficient calculated at each iteration; Repeat the previous step until the energy spectrum norm and residual norm that meet the preset convergence condition are obtained when the current regularization coefficient is selected, and record the energy spectrum norm and residual norm that meet the preset convergence condition when the current regularization coefficient is selected.

[0052] In a possible embodiment, the model unit 701 generates the first inversion model, including: Obtain the electron response function and proton response function of the detector; Determine a response matrix based on the electron response function and the proton response function; Determine a basic inversion model based on the response matrix and the detection value of the detector; Performing transformation based on the basic inversion model to generate a first inversion model; A preset preliminary solution of the first inversion model is determined.

[0053] In a possible embodiment, the response matrix is ​​as follows: , (1) in, To measure the total number of energy channels, is the total number of electron incident energy channels, is the total number of proton incident energy channels; For electrons in the channel The response function of For protons in the channel The response function of . is the number of the measurement energy level, is the number of the incident energy channel, The value range is 1 to The natural number of The value range is 1 to The natural number of The basic inversion model is as follows: (2) in, is the response matrix, is the electron spectrum or proton spectrum, is the measurement error, is the measured value.

[0054] In a possible embodiment, the first inversion model has the following formula: (3) Among them, A T is the inverse matrix of A; μ is the local relaxation coefficient, which is used to adjust the update amplitude of the solution. Each element is ; is the first regularization coefficient, is the second regularization coefficient; is the particle energy spectrum, is the electron energy spectrum, Proton spectrum.

[0055] The particle energy spectrum inversion method and system based on double regularization coefficients provided in the embodiments of the present application use a single detector in combination with the particle energy spectrum inversion method based on double regularization coefficients, thereby reducing the number of detectors and being able to invert electron and proton energy spectra at one time, thereby improving the inversion efficiency.

[0056] It should be understood that the terms used in this application, such as "unit", are a method for distinguishing different components at different levels or different functional divisions. However, if other terms can achieve the same purpose, they may also be used in this application to replace the above terms.

[0057] The terms such as "first" or "second" described in this application are only used to distinguish the relationship between the components, and do not limit them to be necessarily different. If other terms can achieve the same purpose, these other terms may be used in this application to replace the above terms.

[0058] Although the present application has been illustrated and described with specific embodiments, it should be appreciated that the above embodiments are only used to illustrate the technical solutions of the present application rather than to limit them. A person skilled in the art should understand that the technical solutions described in the above embodiments may be modified or some or all of the technical features thereof may be replaced by equivalents without departing from the spirit and scope of the present application. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present application. Therefore, this means that all such replacements and modifications within the scope of the present application are included in the appended claims.

Claims

1. A particle energy spectrum inversion method based on double canonical coefficients, wherein the detector is a single detector, and the detector detects electrons and protons simultaneously, characterized in that: include: Based on the detector, a first inversion model is generated; wherein the first inversion model includes a local relaxation coefficient and a dual regularization coefficient, and the dual regularization coefficient includes a first regularization coefficient corresponding to electrons and a second regularization coefficient corresponding to protons; Preset a local relaxation coefficient set of the local relaxation coefficient and a dual regularization coefficient set of the dual regularization coefficient, and based on the local relaxation coefficient set and the dual regularization coefficient set, use a Boosted-Gold algorithm based on the dual regularization coefficient to iteratively calculate the first inversion model, and record an energy spectrum norm and a residual norm that meet the preset convergence conditions; An L curve is generated based on the spectral norm and the residual norm, and an inflection point of the L curve is used as an inversion optimal solution.

2. The particle energy spectrum inversion method based on double canonical coefficients according to claim 1, characterized in that: The iterative calculation of the first inversion model based on the local relaxation coefficient set and the dual regularization coefficient set using the Boosted-Gold algorithm with dual regularization coefficients includes: Traversing the set of biregular coefficients; When traversing the set of dual regular coefficients, each time the current dual regular coefficient is selected, the local relaxation coefficient is iteratively calculated; Iteratively updating the particle energy spectrum in the first inversion model according to a preset preliminary solution based on the current dual regularization coefficient and the current local relaxation coefficient calculated in each iteration; Repeat the previous step until the energy spectrum norm and the residual norm that meet the preset convergence condition are obtained when the current regularization coefficient is selected, and record the energy spectrum norm and the residual norm that meet the preset convergence condition when the current regularization coefficient is selected.

3. The particle energy spectrum inversion method based on double canonical coefficients according to claim 2, characterized in that: The generating the first inversion model comprises: Obtaining an electron response function and a proton response function of the detector; determining a response matrix based on the electron response function and the proton response function; Determining the basic inversion model based on the response matrix and the detection value of the detector; Performing transformation based on the basic inversion model to generate the first inversion model; The preset preliminary solution of the first inversion model is determined.

4. The particle energy spectrum inversion method based on double canonical coefficients according to claim 3, characterized in that: The response matrix is ​​as follows: ， （1） in, To measure the total number of energy channels, is the total number of electron incident energy channels, is the total number of proton incident energy channels; For electrons in the channel The response function of For protons in the channel The response function of is the number of the measurement energy level, is the number of the incident energy channel, The value range is 1 to The natural number of The value range is 1 to The natural number of The basic inversion model is as follows: （2） in, is the response matrix, is the electron spectrum or proton spectrum, is the measurement error, is the measured value.

5. The particle energy spectrum inversion method based on double canonical coefficients according to claim 4, characterized in that: The first inversion model is formulated as follows: （3） Among them, A T is the inverse matrix of A; μ is the local relaxation coefficient, which is used to adjust the update amplitude of the solution. Each element is ; is the first regularization coefficient, is the second regularization coefficient; is the particle energy spectrum, is the electron energy spectrum, Proton spectrum.

6. A particle energy spectrum inversion system based on dual canonical coefficients, wherein the detector is a single detector, and the detector detects electrons and protons simultaneously, characterized in that: include: A model unit, configured to generate a first inversion model based on the detector; wherein the first inversion model comprises a local relaxation coefficient and a dual regularization coefficient, and the dual regularization coefficient comprises a first regularization coefficient corresponding to electrons and a second regularization coefficient corresponding to protons; A processing unit, used for presetting a local relaxation coefficient set of the local relaxation coefficient and a dual regularization coefficient set of the dual regularization coefficient, and based on the local relaxation coefficient set and the dual regularization coefficient set, using a Boosted-Gold algorithm based on the dual regularization coefficient to iteratively calculate the first inversion model, and recording an energy spectrum norm and a residual norm that meet a preset convergence condition; A result unit is used to generate an L curve based on the spectral norm and the residual norm, and take the inflection point of the L curve as the best inversion solution.

7. The particle energy spectrum inversion system based on double canonical coefficients according to claim 6, characterized in that: The processing unit performs iterative calculation on the first inversion model, including: Traversing the set of biregular coefficients; When traversing the set of dual regular coefficients, each time the current dual regular coefficient is selected, the local relaxation coefficient is iteratively calculated; Iteratively updating the particle energy spectrum in the first inversion model according to a preset preliminary solution based on the current dual regularization coefficient and the current local relaxation coefficient calculated in each iteration; Repeat the previous step until the energy spectrum norm and the residual norm that meet the preset convergence condition are obtained when the current regularization coefficient is selected, and record the energy spectrum norm and the residual norm that meet the preset convergence condition when the current regularization coefficient is selected.

8. The particle energy spectrum inversion system based on double canonical coefficients according to claim 7, characterized in that: The model unit generates the first inversion model, including: Obtaining an electron response function and a proton response function of the detector; determining a response matrix based on the electron response function and the proton response function; Determining the basic inversion model based on the response matrix and the detection value of the detector; Performing transformation based on the basic inversion model to generate the first inversion model; The preset preliminary solution of the first inversion model is determined.

9. The particle energy spectrum inversion system based on double canonical coefficients according to claim 8, characterized in that: The response matrix is ​​as follows: ， （1） in, To measure the total number of energy channels, is the total number of electron incident energy channels, is the total number of proton incident energy channels; For electrons in the channel The response function of For protons in the channel The response function of is the number of the measurement energy level, is the number of the incident energy channel, The value range is 1 to The natural number of The value range is 1 to The natural number of The basic inversion model is as follows: （2） in, is the response matrix, is the electron spectrum or proton spectrum, is the measurement error, is the measured value.

10. The particle energy spectrum inversion system based on double canonical coefficients according to claim 9, characterized in that: The first inversion model is formulated as follows: （3） Among them, A T is the inverse matrix of A; μ is the local relaxation coefficient, which is used to adjust the update amplitude of the solution. Each element is ; is the first regularization coefficient, is the second regularization coefficient; is the particle energy spectrum, is the electron energy spectrum, Proton spectrum.

Images