CN108132260A - A kind of power spectrum analytic method based on linear superposition model - Google Patents

A kind of power spectrum analytic method based on linear superposition model Download PDF

Info

Publication number
CN108132260A
CN108132260A CN201711351491.4A CN201711351491A CN108132260A CN 108132260 A CN108132260 A CN 108132260A CN 201711351491 A CN201711351491 A CN 201711351491A CN 108132260 A CN108132260 A CN 108132260A
Authority
CN
China
Prior art keywords
power spectrum
sup
lin
particle
model
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
Application number
CN201711351491.4A
Other languages
Chinese (zh)
Other versions
CN108132260B (en
Inventor
黄洪全
马兴科
彭亚娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN201711351491.4A priority Critical patent/CN108132260B/en
Publication of CN108132260A publication Critical patent/CN108132260A/en
Application granted granted Critical
Publication of CN108132260B publication Critical patent/CN108132260B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Molecular Biology (AREA)
  • Pathology (AREA)
  • Data Mining & Analysis (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Biochemistry (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Feedback Control In General (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of power spectrum analytic methods based on linear superposition model.For multielement and there are problems that the complicated power spectrum parsing for the destabilizing factor that power spectrum is caused to distort, the characteristic of Momentum profiles forming process, it is proposed a kind of linear superposition model (Linear Superposition), natural selection is introduced into algorithm, parameter optimization is carried out using the improved Particle Swarm Algorithm of Stochastic inertia weight.Each particle corresponds to a LIN SUP model, and calculate the fitness value of each LIN SUP model, the processes such as the assessment of initialization, particle quality by population and the update of particle " flight " speed and the iteration of position, search the LIN SUP models with " adaptive optimal control angle value ", i.e. global optimum LIN SUP models.The result shows that this method has faster convergence rate and higher search precision, the parsing of various power spectrums can be widely used for.

Description

A kind of power spectrum analytic method based on linear superposition model
Technical field
The present invention relates to a kind of power spectrum analytic methods based on linear superposition model.
Background technology
In radioactivity survey is carried out, due to radioactive source there are multielement, detector energy resolution ratio is relatively low and power spectrum not Stablize and (such as drift about) reason, power spectrum is caused to become complicated：Shape such as spectrum becomes strange, adjacent spectral peaks are overlapped, this gives The parsing (qualitative and quantitative analysis of such as ingredient) of power spectrum is made troubles.Once there is researcher to be calculated using gauss hybrid models iteration The method that method, genetic algorithm, curve-fitting method and wavelet transformation and neural network are combined, is fitted overlapping spectra. But at present for multielement and in the presence of the complex situations of destabilizing factor that power spectrum is caused to distort, the report of also rare power spectrum parsing Road.
The present invention is for above-mentioned multielement and there are problems that the complicated power spectrum parsing for the destabilizing factor that power spectrum is caused to distort, The characteristic of Momentum profiles forming process proposes a kind of linear superposition model (Linear Superposition), by natural selection It is introduced into algorithm, parameter optimization is carried out using the improved Particle Swarm Algorithm of Stochastic inertia weight.The result shows that this method has Faster convergence rate and higher search precision, can be widely used for the parsing of various power spectrums.
Invention content
It is an object of the invention to disclose a kind of power spectrum analytic method based on linear superposition model.For multielement and deposit In the parsing problem of the complicated power spectrum of destabilizing factor that power spectrum is caused to distort, there is faster convergence rate and higher search Precision compensates for the deficiency of current power spectrum analytic method.Be by it is following 1.~3. step realize.
1. the power spectrum parsed is normalized, the power spectrum that area is 1 is obtained.
2. the power spectrum after normalization is characterized with a kind of linear superposition model (Linear Superposition), herein Should " linear superposition model " be named as LIN-SUP models.
3. using the collective search technology of particle cluster algorithm, each particle corresponds to a LIN-SUP model, and calculates every The fitness value of one LIN-SUP model；Initialization, the assessment of particle quality and particle " flight " of the algorithm Jing Guo population The processes such as speed and the update of the iteration of position, search the LIN-SUP models with " adaptive optimal control angle value ", i.e. global optimum LIN-SUP models.
By above 1.~3. step acquire global optimum position, the LIN-SUP model parameters corresponding to the position are exactly institute Parse the solution of power spectrum.
The beneficial effects of the invention are as follows：
For multielement and there are problems that the complicated power spectrum parsing for the destabilizing factor that power spectrum is caused to distort, Momentum profiles shape Into the characteristic of process, a kind of linear superposition model (abbreviation LIN-SUP models) is proposed.It is searched using the population of Stochastic inertia weight Rope technology, natural selection is introduced into algorithm, each particle corresponds to a LIN-SUP model, and calculates each LIN- The fitness value of SUP models, the assessment of initialization, particle quality by population and particle " flight " speed and position The processes such as iteration update search the LIN-SUP models with " adaptive optimal control angle value ", i.e. global optimum LIN-SUP models, real Now to the parsing of complicated power spectrum.The result shows that this method has faster convergence rate and higher search precision, can use extensively In the parsing of various power spectrums.
Description of the drawings
Fig. 1 is the flow chart of the method for the present invention.
Specific embodiment
It elaborates below in conjunction with the accompanying drawings to the embodiment of the present invention：The present embodiment using technical solution of the present invention before It puts and is implemented, give detailed embodiment and process, but protection scope of the present invention is not limited to following embodiments.
Present embodiment assumes that the multielement complexity power spectrum obtained in radioactivity survey is F0(x), wherein x is road location, is adopted Power spectrum is parsed with this method by following specific steps 1.~3..
The step 1. power spectrum F to being parsed0(x) it is normalized, obtains the power spectrum F (x) that area is 1.
2. step characterizes the power spectrum F (x) after normalization with linear superposition model (Linear Superposition), Herein should " linear superposition model " be named as LIN-SUP models：
In formula (1)：M represents element species；Fi(x) (i=1 ... M) represents the corresponding power spectrum of the i-th dvielement；NiRepresent the The number of the sub- power spectrum of i dvielements；fij(x) j-th of sub- power spectrum of the i-th dvielement, weight a are representedij, and meet：
θ represents the set of LIN-SUP model parameters：
3. step uses the collective search technology of particle cluster algorithm, each particle corresponds to a LIN-SUP model, and counts Calculate the fitness value of each LIN-SUP model；Initialization of the algorithm Jing Guo population, the assessment of particle quality and particle The processes such as " flight " speed and the update of the iteration of position, search the LIN-SUP models with " adaptive optimal control angle value ", i.e., global Optimal L IN-SUP models.Specifically realized by following steps A~D.
A generates primary group：Model parameter θ is used for represent particleDimension space position；If particle number is Pnum；Being created on dimension space position has equally distributed initial population；Set each particleThe initial speed of dimension V is spent, it is correspondingThe speed of dimension space position.
B Fitness analysis：Calculate the fitness value y of each particlek(θk), as the following formula
Wherein yk(θk) represent the fitness value of k-th of particle, beginx and endx represent respectively power spectrum start channel location and Terminate road location；yk(θk) value is smaller, position is more excellent.
The speed of C more new particles and position, as the following formula
Wherein, Vij(t+1),Xij(t+1) speed and position of i-th of particle in t+1 iteration in jth dimension are represented respectively It puts；pijAnd pgjThe individual optimal value and global optimum of i-th of particle at the end of the t times iteration are represented respectively；c1And c2Respectively For Studying factors；r1And r2Uniform random number respectively in the range of [0,1]；β claims constraint factor, for the weight regulated the speed； ω is inertia weight, is obtained as the following formula；
Wherein N (0,1) is represented to obey the random number of standardized normal distribution, is uniformly distributed between rand (0,1) expressions (0,1) Random number, μmaxWith μminThe upper and lower bound of the parameter μ of normal distribution is represented respectively, and δ represents the side of random weighted mean Difference.
D sorts to population according to fitness value, worst half particle is replaced with half particle best in group, together When retain the history optimal value that original each individual is remembered.
When algorithm reaches stop condition, then stop search and export result；Otherwise B is returned to continue search for.
By above 1.~3. step acquire global optimum position, the LIN-SUP model parameters corresponding to the positionIt is exactly solved, completes the parsing of power spectrum.
For multielement and there are problems that the complicated power spectrum parsing for the destabilizing factor that power spectrum is caused to distort, Momentum profiles shape Into the characteristic of process, propose a kind of linear superposition model (Linear Superposition), natural selection be introduced into algorithm, Parameter optimization is carried out using the improved Particle Swarm Algorithm of Stochastic inertia weight.Each particle corresponds to a LIN-SUP model, And the fitness value of each LIN-SUP model is calculated, the assessment of initialization, particle quality by population and particle The processes such as " flight " speed and the update of the iteration of position, search the LIN-SUP models with " adaptive optimal control angle value ", i.e., global Optimal L IN-SUP models.The result shows that this method has faster convergence rate and higher search precision, can be widely used for The parsing of various power spectrums.
In embodiments of the invention described above, the parsing of multielement complexity power spectrum is described in detail, but need to illustrate , the foregoing is merely one embodiment of the present of invention, the present invention can equally carry out the power spectrum of other various rays Parsing, all within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in this hair Within bright protection domain.

Claims (1)

1. a kind of power spectrum analytic method based on linear superposition model, which is characterized in that step is as follows：
1. the power spectrum parsed is normalized, the power spectrum F (x) that area is 1 is obtained；
2. the power spectrum F (x) after normalization is characterized with a kind of linear superposition model (hereinafter referred to as LIN-SUP)：
In formula (1)：M represents element species；Fi(x) (i=1 ... M) represents the corresponding power spectrum of the i-th dvielement；NiRepresent the i-th class The number of first sub-prime power spectrum；fij(x) j-th of sub- power spectrum of the i-th dvielement, weight a are representedij, and meet：
θ represents the set of LIN-SUP model parameters：
3. using the collective search technology of particle cluster algorithm, global optimum's LIN-SUP models are found, method is：By LIN-SUP moulds Position of the parameter θ of type as particle in space by particle " flight " speed and the iteration renewal process of " position ", will be searched Parameter of the location parameter with " global optimum's fitness value " that rope arrives as global optimum's LIN-SUP models；Fitness value Calculating as the following formula：
Wherein yk(θk) representing the fitness value of k-th of particle, beginx and endx represent start channel location and the end of power spectrum respectively Road location；yk(θk) value is smaller, position is more excellent.
CN201711351491.4A 2017-12-15 2017-12-15 Energy spectrum analysis method based on linear superposition model Expired - Fee Related CN108132260B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711351491.4A CN108132260B (en) 2017-12-15 2017-12-15 Energy spectrum analysis method based on linear superposition model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711351491.4A CN108132260B (en) 2017-12-15 2017-12-15 Energy spectrum analysis method based on linear superposition model

Publications (2)

Publication Number Publication Date
CN108132260A true CN108132260A (en) 2018-06-08
CN108132260B CN108132260B (en) 2020-04-28

Family

ID=62390371

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711351491.4A Expired - Fee Related CN108132260B (en) 2017-12-15 2017-12-15 Energy spectrum analysis method based on linear superposition model

Country Status (1)

Country Link
CN (1) CN108132260B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109799541A (en) * 2019-01-25 2019-05-24 中国自然资源航空物探遥感中心 A kind of measurement spectrum drift of gamma spectrum and integral nonlinearity bearing calibration

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59228134A (en) * 1983-06-10 1984-12-21 Shimadzu Corp Device for decomposing and processing spectrum
CN101788507A (en) * 2010-02-03 2010-07-28 北京矿冶研究总院 Spectral analysis method
CN102298153A (en) * 2010-06-23 2011-12-28 成都理工大学 Method for decomposing multiple spectral peaks during radioactive measurement
CN102608649A (en) * 2012-03-02 2012-07-25 成都理工大学 Statistics distributed gamma or X ray energy spectrum unscrambling method
CN103401625A (en) * 2013-08-23 2013-11-20 西安电子科技大学 Particle swarm optimization algorithm based cooperative spectrum sensing optimization method
CN104239709A (en) * 2014-09-05 2014-12-24 西安奥华电子仪器有限责任公司 Method for determining yield of stratum element logging instrument by virtue of spectrum unfolding
CN104422704A (en) * 2013-08-21 2015-03-18 同方威视技术股份有限公司 Method of decomposing energy spectrum information of X-ray energy spectrum CT and corresponding reconstruction method
CN105203565A (en) * 2014-06-11 2015-12-30 成都理工大学 Energy spectrum overlapping peak analysis method
CN105989410A (en) * 2015-03-05 2016-10-05 成都理工大学 Overlap kernel pulse separation method
CN106405624A (en) * 2016-08-30 2017-02-15 天津大学 Medical CT oriented X-ray energy spectrum analysis method through reconstruction
CN107229787A (en) * 2017-05-24 2017-10-03 南京航空航天大学 A kind of gamma-ray spectrum analysis method based on approximation coefficient and deep learning

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59228134A (en) * 1983-06-10 1984-12-21 Shimadzu Corp Device for decomposing and processing spectrum
CN101788507A (en) * 2010-02-03 2010-07-28 北京矿冶研究总院 Spectral analysis method
CN102298153A (en) * 2010-06-23 2011-12-28 成都理工大学 Method for decomposing multiple spectral peaks during radioactive measurement
CN102608649A (en) * 2012-03-02 2012-07-25 成都理工大学 Statistics distributed gamma or X ray energy spectrum unscrambling method
CN104422704A (en) * 2013-08-21 2015-03-18 同方威视技术股份有限公司 Method of decomposing energy spectrum information of X-ray energy spectrum CT and corresponding reconstruction method
CN103401625A (en) * 2013-08-23 2013-11-20 西安电子科技大学 Particle swarm optimization algorithm based cooperative spectrum sensing optimization method
CN105203565A (en) * 2014-06-11 2015-12-30 成都理工大学 Energy spectrum overlapping peak analysis method
CN104239709A (en) * 2014-09-05 2014-12-24 西安奥华电子仪器有限责任公司 Method for determining yield of stratum element logging instrument by virtue of spectrum unfolding
CN105989410A (en) * 2015-03-05 2016-10-05 成都理工大学 Overlap kernel pulse separation method
CN106405624A (en) * 2016-08-30 2017-02-15 天津大学 Medical CT oriented X-ray energy spectrum analysis method through reconstruction
CN107229787A (en) * 2017-05-24 2017-10-03 南京航空航天大学 A kind of gamma-ray spectrum analysis method based on approximation coefficient and deep learning

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIN YU,ET AL: "A method for separation of overlapping absorption lines in intracavitygas detection", 《SENSORS AND ACTUATORS B》 *
吴和喜 等: "基于最小二乘法的航空能谱解析", 《核技术》 *
赵立林 等: "神经网络算法解析便携式高纯锗γ 谱", 《东华理工大学学报( 自然科学版)》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109799541A (en) * 2019-01-25 2019-05-24 中国自然资源航空物探遥感中心 A kind of measurement spectrum drift of gamma spectrum and integral nonlinearity bearing calibration

Also Published As

Publication number Publication date
CN108132260B (en) 2020-04-28

Similar Documents

Publication Publication Date Title
CN107871155B (en) Spectral overlapping peak decomposition method based on particle swarm optimization
CN113885555B (en) Multi-machine task allocation method and system for transmission line dense channel inspection
CN107392397A (en) A kind of short-term wind speed forecasting method, apparatus and system
CN111292525A (en) Traffic flow prediction method based on neural network
CN104765690A (en) Embedded software test data generating method based on fuzzy-genetic algorithm
CN102682219A (en) Method for forecasting short-term load of support vector machine
CN110427690A (en) A kind of method and device generating ATO rate curve based on global particle swarm algorithm
CN108320504B (en) Dynamic OD matrix estimation method based on monitoring data
CN108132260A (en) A kind of power spectrum analytic method based on linear superposition model
Diefenbacher et al. Refining Fast Calorimeter Simulations with a Schr\"{o} dinger Bridge
Li et al. Many-objective rapid optimization of reactor shielding design based on NSGA-III
CN105578472A (en) Wireless sensor network performance online optimization and planning method based on immune theory
Fülöp A method for approximating pairwise comparison matrices by consistent matrices
Tang et al. Intelligent vehicle lateral tracking control based on multiple model prediction
CN113656707A (en) Financing product recommendation method, system, storage medium and equipment
CN117575870A (en) Urban carbon balance optimization method based on SD-FLUS model
CN107276093A (en) Power system probability load flow calculation method based on scene reduction
Napolitani et al. Boltzmann-Langevin one-body dynamics for fermionic systems
CN113536685B (en) Neural network extrapolation-based wind speed probability model modeling method
CN105243446A (en) Electricity consumption combined forecasting method based on particle swarm optimization
CN113689694B (en) Traffic flow prediction method, device, equipment and readable storage medium
Huang et al. A Matrix Translation Model for Evacuation Path Optimization
Du et al. AGRU and convex optimization based energy management for plug-in hybrid electric bus considering battery aging
CN105469644A (en) Flight conflict resolution method and flight conflict resolution device
CN112634620A (en) Road network dynamic traffic distribution method based on Encoder-Decoder deep neural network

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200428

Termination date: 20201215

CF01 Termination of patent right due to non-payment of annual fee