CN108663707A - A kind of repeatedly two-way S-K smoothing processings system and method - Google Patents

A kind of repeatedly two-way S-K smoothing processings system and method Download PDF

Info

Publication number
CN108663707A
CN108663707A CN201810285498.9A CN201810285498A CN108663707A CN 108663707 A CN108663707 A CN 108663707A CN 201810285498 A CN201810285498 A CN 201810285498A CN 108663707 A CN108663707 A CN 108663707A
Authority
CN
China
Prior art keywords
buffer
filters
reversed
filtering
repeatedly
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
CN201810285498.9A
Other languages
Chinese (zh)
Other versions
CN108663707B (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 CN201810285498.9A priority Critical patent/CN108663707B/en
Publication of CN108663707A publication Critical patent/CN108663707A/en
Application granted granted Critical
Publication of CN108663707B publication Critical patent/CN108663707B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36Measuring spectral distribution of X-rays or of nuclear radiation spectrometry

Abstract

The present invention relates to a kind of repeatedly two-way S-K smoothing processings system and methods, first buffer, forward direction number SK filters, the second buffer and reversed number SK filters are sequentially connected, the reversed number SK filters are connected with the first buffer, and first buffer, forward direction number SK filters, the second buffer and reversed number SK filters are respectively connected with control register and parameter register.It is provided by the invention that repeatedly two-way S-K smoothing processings system and method can be identical as the effect of multiple Gaussian smoothing, compared with Gaussian smoothing, the execution efficiency higher of the present invention, algorithm calculating process is simple, it is easy to use, it is easy to software and hardware realization, present invention can apply to nuclear instrument and actual measurement nuclear energy spectrum analysis, to solve the problem of that characteristics such as Electronics noice and statistic fluctuation in existing nuclear spectrum measurement influence the qualitative and quantitative analysis of low content nucleic.

Description

A kind of repeatedly two-way S-K smoothing processings system and method
Technical field
The present invention relates to nuclear spectrum analysis fields, and in particular to a kind of repeatedly two-way S-K smoothing processings system and method.
Background technology
With the progress and development of science, the application field of Nuclear Instruments And Industrial Gauges is also more and more extensive, almost covers agricultural, work The fields such as industry, health care, geology and mine, environmental protection, space flight, teaching, scientific research, customs.In nuclear decay and nuclear spectrum measurement Statistic fluctuation characteristic, can make count it is smaller when, real Gaussian Profile desired value can be submerged in noise, though When counting larger, this statistic fluctuation characteristic also can make measured spectrum data burr phenomena occur, influence later data analysis;For The presence of reliable recognition low content nucleic characteristic peak, while carrying out qualitative analysis to peak energy amount and position needs research new Energy spectrum analysis method.
S-K filters propose that low pass S-K filters are in nuclear electronics in nineteen fifty-five by R.P.Sallen and E.L.Key It is widely used in the filtering forming of core pulse signal, the output of approximate Gaussian waveform can be obtained under less series, When by nuclear spectrum data back by digital S-K filters, spectral peak can occur to drift about to the left.In live data analysis, for The good energy spectrum analysis tool of scale, spectral peak drift are unfavorable for data analysis.
Invention content
Technical problem to be solved by the invention is to provide a kind of repeatedly two-way S-K smoothing processings system and methods, solve The characteristics such as Electronics noice and statistic fluctuation in existing nuclear spectrum measurement, influence asking for the qualitative and quantitative analysis of low content nucleic Topic.
The technical solution that the present invention solves above-mentioned technical problem is as follows:A kind of repeatedly two-way S-K smoothing processings system, packet It includes:
First buffer is kept in filtering for the temporary input traffic when filtering beginning through reversed digital Data flow after SK filters, and output stream after filtering;
Second buffer, for the temporary data flow after positive number SK filters;
Positive number SK filters, for forward direction SK filtering;
Reversed number SK filters, filter for reversed SK;
Control register, for controlling filter times and data flow direction;
Parameter register, for storing filtering parameter;
First buffer, forward direction number SK filters, the second buffer and reversed number SK filters are sequentially connected, The reversed number SK filters are connected with the first buffer, first buffer, forward direction number SK filters, the second caching Device and reversed number SK filters are respectively connected with control register and parameter register.
Further, the output signal of second buffer inputs reversed number SK filters after direction is corrected.
A kind of repeatedly two-way S-K smoothing processing methods, include the following steps:
S1, by data stream to be filtered to the first buffer;
S2, by the data flow of the first buffer by being sent to the second buffer after positive number SK filter filterings Caching;
S3, the data flow cached in the second buffer is corrected into line direction, then passes through reversed number SK filter filterings After be sent to the first buffer;
S4, filter times are compared with the filtering parameter stored in parameter register, if filter times are less than filtering and join Number, then enter step S2, if filter times are more than filtering parameter, enters step S5;
S5, the data flow cached in the first buffer is exported.
The beneficial effects of the invention are as follows:It is provided by the invention repeatedly two-way S-K smoothing processings system and method can with it is multiple The effect of Gaussian smoothing is identical, and compared with Gaussian smoothing, execution efficiency higher of the invention, algorithm calculating process is simple, It is easy to use, it is easy to software and hardware realization, it is existing to solve present invention can apply to nuclear instrument and actual measurement nuclear energy spectrum analysis The characteristics such as Electronics noice and statistic fluctuation in nuclear spectrum measurement, the problem of influencing the qualitative and quantitative analysis of low content nucleic.
Description of the drawings
Fig. 1 is principle of the invention figure;
Fig. 2 is flow chart of the present invention;
Fig. 3 is that cadmium rice energy translation measures actual measurement spectral line region of interest enlarged drawing in the embodiment of the present invention;
Fig. 4 is that the multiple Gaussian smoothing of region of interest spectral line and multiple two-way S-K smoothing processings are flat in the embodiment of the present invention Slip type analysis figure;
Fig. 5 is that the cadmium rice region of interest spectral line after SNIP backgrounds is deducted in the embodiment of the present invention.
Specific implementation mode
The principle and features of the present invention will be described below with reference to the accompanying drawings, and the given examples are served only to explain the present invention, and It is non-to be used to limit the scope of the present invention.
As shown in Figure 1, a kind of smooth (referred to as MBSK) processing systems of repeatedly two-way S-K, including:
First buffer is kept in filtering for the temporary input traffic when filtering beginning through reversed digital Data flow after SK filters, and output stream after filtering;
Second buffer, for the temporary data flow after positive number SK filters;
Positive number SK filters, for forward direction SK filtering;
Reversed number SK filters, filter for reversed SK;
Control register, for controlling filter times and data flow direction;
Parameter register, for storing filtering parameter;
First buffer, forward direction number SK filters, the second buffer and reversed number SK filters are sequentially connected, The reversed number SK filters are connected with the first buffer, first buffer, forward direction number SK filters, the second caching Device and reversed number SK filters are respectively connected with control register and parameter register.
The output signal of second buffer inputs reversed number SK after the correction of direction first-in last-out by pressing data and filters Wave device.
As shown in Fig. 2, a kind of repeatedly two-way S-K smoothing processing methods, include the following steps:
S1, by data stream to be filtered to the first buffer;
S2, by the data flow of the first buffer by being sent to the second buffer after positive number SK filter filterings Caching;
S3, the data flow cached in the second buffer is corrected into line direction, then passes through reversed number SK filter filterings After be sent to the first buffer;
S4, filter times are compared with the filtering parameter stored in parameter register, if filter times are less than filtering and join Number, then enter step S2, if filter times are more than filtering parameter, enters step S5;
S5, the data flow cached in the first buffer is exported.
In embodiments of the present invention, using multiple two-way S-K filtering techniques, basic principle is:Make forward filtering first, Nuclear energy modal data is first passed through into digital S-K filters by sequence from the front to the back;Then inverse filtering is remake, i.e., by previous filter Gamma-spectrometric data stream after wave is reversed, is allowed to from the front to the back through digital S-K filters.The primary two-way filtered nuclear energy of SK of operation Modal data, spectrum peak position is constant, facilitates and carries out data analysis with the energy calibration information that proper energy is composed.
Cadmium rice is measured using 2048 X-ray Fluorescence Analyzers, region of interest is intercepted after measurement (roads 1605-1665), spectral line is as shown in Figure 3.
In embodiments of the present invention, after suitable filtering parameter is arranged, multiple two-way S-K smoothing processings can be with multiple Gauss The effect of smoothing processing is identical, and compared with Gaussian smoothing, execution efficiency higher of the invention, algorithm calculating process is simple, uses It is convenient, it is easy to software and hardware realization.
Multiple Gaussian smoothing parameter setting is as follows:Window width is 7, standard deviation 1.25, filter times 55;It is repeatedly double The filtering parameter smooth to S-K be:K=1.02, a=1, filter times 10.Using MATLAB2016b, record respectively above-mentioned The Gaussian smoothing of parameter and multiple run time smooth two-way S-K, can obtain:Gaussian smoothing run time 0.0122s, it is repeatedly double To S-K smoothness run time 0.0022s, it is known that the method for the invention efficiency of code execution higher.Experiment is it is found that two kinds smooth Data height overlap, the multiple two-way S-K of region of interest is smoothly done into linear fit with data after Gaussian smoothing, such as Fig. 4 institutes Show, goodness of fit R can be obtained2It is 0.99999, as shown in figure 5, cadmium rice original measurement spectral line and the multiple two-way S-K of warp is smooth Spectral line afterwards uses SNIP method background corrections, then intercepts the spectral line of region of interest respectively.
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all the present invention spirit and Within principle, any modification, equivalent replacement, improvement and so on should all be included in the protection scope of the present invention.

Claims (3)

1. a kind of repeatedly two-way S-K smoothing processings system, which is characterized in that including:
First buffer is kept in filtering and is filtered through reversed number SK for the temporary input traffic when filtering beginning Data flow after wave device, and output stream after filtering;
Second buffer, for the temporary data flow after positive number SK filters;
Positive number SK filters, for forward direction SK filtering;
Reversed number SK filters, filter for reversed SK;
Control register, for controlling filter times and data flow direction;
Parameter register, for storing filtering parameter;
First buffer, forward direction number SK filters, the second buffer and reversed number SK filters are sequentially connected, described Reversed number SK filters are connected with the first buffer, first buffer, forward direction number SK filters, the second buffer and Reversed number SK filters are respectively connected with control register and parameter register.
2. repeatedly two-way S-K smooth processing units according to claim 1, which is characterized in that second buffer Output signal inputs reversed number SK filters after direction is corrected.
3. a kind of repeatedly two-way S-K smoothing processing methods, which is characterized in that include the following steps:
S1, by data stream to be filtered to the first buffer;
S2, the data flow of the first buffer is delayed by being sent to the second buffer after positive number SK filter filterings It deposits;
S3, the data flow cached in the second buffer is corrected into line direction, then by being passed after reversed number SK filter filterings It send to the first buffer;
S4, filter times are compared with the filtering parameter stored in parameter register, if filter times are less than filtering parameter, S2 is entered step, if filter times are more than filtering parameter, enters step S5;
S5, the data flow cached in the first buffer is exported.
CN201810285498.9A 2018-04-02 2018-04-02 Multi-time bidirectional S-K smoothing processing system and method Active CN108663707B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810285498.9A CN108663707B (en) 2018-04-02 2018-04-02 Multi-time bidirectional S-K smoothing processing system and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810285498.9A CN108663707B (en) 2018-04-02 2018-04-02 Multi-time bidirectional S-K smoothing processing system and method

Publications (2)

Publication Number Publication Date
CN108663707A true CN108663707A (en) 2018-10-16
CN108663707B CN108663707B (en) 2020-12-08

Family

ID=63782744

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810285498.9A Active CN108663707B (en) 2018-04-02 2018-04-02 Multi-time bidirectional S-K smoothing processing system and method

Country Status (1)

Country Link
CN (1) CN108663707B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110135402A (en) * 2019-06-04 2019-08-16 东华理工大学 A kind of nuclear spectrum data smoothing method based on multi-parameters model
CN113189634A (en) * 2021-03-02 2021-07-30 四川新先达测控技术有限公司 Gaussian-like forming method
CN114511477A (en) * 2022-01-20 2022-05-17 上海盈方微电子有限公司 NEON optimization method suitable for Gaussian smoothing

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102302363A (en) * 2011-06-29 2012-01-04 深圳市理邦精密仪器股份有限公司 Device and method for automatically identifying fetal heart rate baseline
US20130285746A1 (en) * 2011-09-23 2013-10-31 Tensorcom, Inc. Differential Source Follower having 6dB Gain with Applications to WiGig Baseband Filters
CN103941279A (en) * 2014-05-08 2014-07-23 东华理工大学 Nuclear pulse signal digital Gaussian forming method based on double parameters of pulse width and amplitude width
CN105806338A (en) * 2016-03-17 2016-07-27 孙红星 GNSS/INS integrated positioning and directioning algorithm based on three-way Kalman filtering smoother
CN105989410A (en) * 2015-03-05 2016-10-05 成都理工大学 Overlap kernel pulse separation method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102302363A (en) * 2011-06-29 2012-01-04 深圳市理邦精密仪器股份有限公司 Device and method for automatically identifying fetal heart rate baseline
US20130285746A1 (en) * 2011-09-23 2013-10-31 Tensorcom, Inc. Differential Source Follower having 6dB Gain with Applications to WiGig Baseband Filters
US8803596B2 (en) * 2011-09-23 2014-08-12 Tensorcom, Inc. Differential source follower having 6dB gain with applications to WiGig baseband filters
CN103941279A (en) * 2014-05-08 2014-07-23 东华理工大学 Nuclear pulse signal digital Gaussian forming method based on double parameters of pulse width and amplitude width
CN105989410A (en) * 2015-03-05 2016-10-05 成都理工大学 Overlap kernel pulse separation method
CN105806338A (en) * 2016-03-17 2016-07-27 孙红星 GNSS/INS integrated positioning and directioning algorithm based on three-way Kalman filtering smoother

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
C.ARELLANO ET AL.: "Modeling of Opto-Electronics in Complex Photonic Integrated Circuits", 《PROC SPIE》 *
ZHOU JIANBIN ET AL.: "Improvement of digital S-K filter and its application in nuclear signal processing", 《NUCLEAR SCIENCE AND TECHNIQUES》 *
纪跃波 等: "零相位数字滤波器", 《重庆大学学报(自然科学版)》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110135402A (en) * 2019-06-04 2019-08-16 东华理工大学 A kind of nuclear spectrum data smoothing method based on multi-parameters model
CN113189634A (en) * 2021-03-02 2021-07-30 四川新先达测控技术有限公司 Gaussian-like forming method
CN114511477A (en) * 2022-01-20 2022-05-17 上海盈方微电子有限公司 NEON optimization method suitable for Gaussian smoothing

Also Published As

Publication number Publication date
CN108663707B (en) 2020-12-08

Similar Documents

Publication Publication Date Title
CN108663707A (en) A kind of repeatedly two-way S-K smoothing processings system and method
DE112008003554B4 (en) Multi-channel momentum analyzer for use with multi-channel spectrometers
CN106568805A (en) Highly-integrated Langmuir probe diagnosis system and method
CN108983279B (en) Low-background spectrum resolving method based on sodium iodide detector
US5132540A (en) Dead time correction and regulation method and apparatus for nuclear spectroscopy
CN108181648B (en) The quality control method and device of continuous record aliased seismic acquisition data
WO2021174901A1 (en) Molecular omics data structure implementation method based on data independent acquisition mass spectrum
CN108009122B (en) Improved HHT method
CN107219547B (en) Method for deducting gamma energy spectrum full spectrum background by utilizing successive approximation SNIP (selective non-ip) calculation
CN107797154A (en) A kind of electric imaging logging image scale method and device
CN104156565A (en) System state analysis method and analysis device based on offline CAN bus data
CN109634944B (en) Network loss data cleaning method based on multi-dimensional space-time analysis
CN105629290A (en) Method of forming Mexico hat wavelet by digital nuclear pulse signal
Petriska et al. CDBTools-evaluate positron annihilation coincidence Doppler broadening spectrum
CN106656200A (en) Program counter compression method and hardware circuit thereof
CN107220215B (en) Method for calculating initial parameters of particle spectrum fitting function based on moment estimation method
CN111413725B (en) System and method for realizing gamma-gamma digital coincidence measurement by using virtual instrument technology
CN101464420A (en) Portable high-speed multi-channel energy spectrometer
CN109753983A (en) Image classification method, device and computer readable storage medium
CN116319383A (en) Traffic matrix prediction method, device, computer equipment and storage medium
Landsberger et al. NADA92: An automated, user-friendly program for neutron activation data analysis
CN115310472A (en) Nuclear pulse peak sequence-based one-dimensional convolution neural network nuclide identification method
CN109117873A (en) A kind of user behavior analysis method based on Bayesian Classification Arithmetic
CN112035159B (en) Configuration method, device, equipment and storage medium of audit model
EP3676640B1 (en) Methods and systems for calibration of particle detectors

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