CN103941280B - Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure - Google Patents

Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure Download PDF

Info

Publication number
CN103941280B
CN103941280B CN201310629796.2A CN201310629796A CN103941280B CN 103941280 B CN103941280 B CN 103941280B CN 201310629796 A CN201310629796 A CN 201310629796A CN 103941280 B CN103941280 B CN 103941280B
Authority
CN
China
Prior art keywords
formation system
gauss
digital
impulse response
gaussian
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.)
Expired - Fee Related
Application number
CN201310629796.2A
Other languages
Chinese (zh)
Other versions
CN103941280A (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 CN201310629796.2A priority Critical patent/CN103941280B/en
Publication of CN103941280A publication Critical patent/CN103941280A/en
Application granted granted Critical
Publication of CN103941280B publication Critical patent/CN103941280B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Radiation (AREA)

Abstract

The invention discloses a kind of digital core pulse Gauss manufacturing process based on Impulse invariance procedure, comprise: according to the circuit of simulation Gauss formation system, obtain the differential equation of circuit system input signal and output signal, the differential equation is derived by frequency-domain analysis, obtains the unit impulse response of simulation Gaussian Systems; The frequency response of analysis mode Gauss formation system, and frequency that simulation Gauss formation system is sampled is determined according to the spectral range of frequency response; Simulation Gauss formation system Impulse invariance procedure is converted to numeric field by analog domain, obtains the unit impulse response of digital Gaussian formation system, containing the forming parameter that can adjust waveform widths in described impulse response; Forming parameter in selected impulse response, and the unit impulse response of corresponding digital Gaussian formation system is determined according to selected forming parameter value; The unit impulse response of the core pulse signal and digital Gaussian formation system that are input to digital Gaussian formation system is carried out Convolution sums computing, realizes digital core pulse Gauss and be shaped.

Description

Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure
Technical field
The Gauss that the present invention relates to digital core pulse in radioactivity survey is shaped, and particularly relates to a kind of digital core pulse Gauss manufacturing process based on Impulse invariance procedure.
Background technology
In nuclear spectrum measurement system, reduce pulse pile-up for convenience of subsequent conditioning circuit process and analysis and improve signal to noise ratio (S/N ratio), the signal shaping that detector need be exported be suitable waveform.The signal to noise ratio (S/N ratio) of Gaussian waveform is high, top is more smooth, ballistic deficit is less, and therefore the normal signal shaping exported by detector is Gauss or accurate Gaussian waveform.Simulation Gauss is shaped and simulation Sallen-Key wave filter can be used to realize, but the formation system realized due to mimic channel, changing over shape parameter must adjust hardware; Digital forming does not then need adjustment System hardware, has higher dirigibility and stability, and digital Gaussian shaping Algorithm becomes the focus of Nuclear signal processing research in recent years.
Summary of the invention
For solving the problems of the technologies described above, the object of this invention is to provide a kind of digital core pulse Gauss manufacturing process based on Impulse invariance procedure.The method overcome the deficiency that core pulse simulation Gauss is shaped, solve the digital Gaussian shaping demand of core pulse, for the digitizing of nuclear instrument is laid a good foundation.
Object of the present invention is realized by following technical scheme:
Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure, comprising:
According to the circuit of simulation Gauss formation system, obtain the differential equation of circuit system input signal and output signal, the differential equation is derived by frequency-domain analysis, obtain the unit impulse response of simulation Gaussian Systems;
The frequency response of analysis mode Gauss formation system, and frequency that simulation Gauss formation system is sampled is determined according to the spectral range of frequency response;
Simulation Gauss formation system Impulse invariance procedure is converted to numeric field by analog domain, obtains the unit impulse response of digital Gaussian formation system, containing the forming parameter that can adjust waveform widths in described impulse response;
Forming parameter in selected impulse response, and the unit impulse response of corresponding digital Gaussian formation system is determined according to selected forming parameter value;
The unit impulse response of the core pulse signal and digital Gaussian formation system that are input to digital Gaussian formation system is carried out Convolution sums computing, realizes digital core pulse Gauss and be shaped.
Compared with prior art, one or more embodiment of the present invention can have the following advantages by tool:
Effectively overcome the deficiency of simulation Gauss formation system, be shaped digital core pulse signal the Gaussian waveform that is as the criterion, and can according to needing flexible parameter value to shaping waveform in using, to meet different measurement demands, the energy resolution of raising system, in digital core energy depressive spectroscopy, the real-time implementation of Gaussian particle filter is laid a good foundation.
Other features and advantages of the present invention will be set forth in the following description, and, partly become apparent from instructions, or understand by implementing the present invention.Object of the present invention and other advantages realize by structure specifically noted in instructions, claims and accompanying drawing and obtain.
Accompanying drawing explanation
Accompanying drawing is used to provide a further understanding of the present invention, and forms a part for instructions, with embodiments of the invention jointly for explaining the present invention, is not construed as limiting the invention.In the accompanying drawings:
Fig. 1 is the process flow diagram of digital Gaussian shaping implementation procedure;
Fig. 2 is the circuit theory diagrams of simulation Gauss formation system;
Fig. 3 is that Gauss corresponding to different parameters value is shaped waveform.
Embodiment
Easy understand, according to technical scheme of the present invention, do not changing under connotation of the present invention, one of ordinary skill in the art can propose multiple frame mode of the present invention and method for making.Therefore following embodiment and accompanying drawing are only illustrating of technical scheme of the present invention, and should not be considered as of the present invention all or be considered as restriction or the restriction of technical solution of the present invention.
Below in conjunction with embodiment and accompanying drawing, the present invention is described in further detail.
As shown in Figure 1, be the digital core pulse Gauss manufacturing process based on Impulse invariance procedure, the method comprises the following steps:
Step 10 is according to the circuit of simulation Gauss formation system, obtain circuit system input signal f (t) and the differential equation outputing signal y (t), the differential equation is derived by frequency-domain analysis, obtains unit impulse response h (t) of simulation Gaussian Systems;
Namely the circuit of above-mentioned simulation Gauss formation system is the circuit of simulation Sallen-Key wave filter, according to the circuit theory (as shown in Figure 2) of simulation Sallen-Key wave filter, the differential equation listed between input signal f (t) and output signal y (t) is:
R 2C 2y"(t)+RCy′(t)+y(t)=2f(t)(1)
The unit impulse response of simulation Gaussian Systems is derived with frequency domain analysis
f ( ω ) = ∫ - ∞ ∞ f ( t ) e - jωt dt - - - ( 2 )
Y ( ω ) = ∫ - ∞ ∞ y ( t ) e - jωt dt - - - ( 3 )
d n y ( t ) dt n ↔ L j ω n Y ( ω ) - - - ( 4 )
According to the frequency response function H (ω) that must simulate Gauss's formation system is:
H ( ω ) = Y ( ω ) F ( ω ) = 2 R 2 C 2 ( jω ) 2 + RCjω + 1 - - - ( 5 )
Carry out inverse Fourier transform to it to obtain:
h ( t ) = 2 e - ( 1 2 FC ) t sin ( B 2 RC ) t 3 RC u ( t ) - - - ( 6 )
The frequency response of step 20 analysis mode Gauss formation system, and frequency that simulation Gauss formation system is sampled is determined according to the spectral range of frequency response;
Guarantee that sample frequency is greater than 2 times of simulation Gauss formation system frequency response bandwidth, when changing to numeric field to avoid analog domain, spectral aliasing occurs.
Simulation Gauss formation system is converted to numeric field by analog domain by step 30, obtains unit impulse response h (n) of digital Gaussian formation system, containing the forming parameter that can adjust waveform widths in described impulse response;
Unit impulse response h (n) of described digital Gaussian formation system just in time equals the sampled value of impulse response h (t) of analog filter, i.e. h (n)=h (t) | t=nT, T is the sampling period.
h ( n ) = 2 ae - a 2 n sin 3 a 2 n 3 u ( n ) - - - ( 7 )
Wherein, T RC = a
Step 40 selectes forming parameter suitable in impulse response, namely determines the value of RC and T, and determines unit impulse response h (n) of corresponding digital Gaussian formation system according to selected forming parameter value;
The unit impulse response of core pulse signal and digital Gaussian formation system that step 50 is input to digital Gaussian formation system carries out Convolution sums computing, realizes digital core pulse Gauss and is shaped, that is:
y ( n ) = f ( n ) * h ( n ) = Σ k = - ∞ ∞ f ( m ) h ( n - m ) - - - ( 8 )
Fig. 3 is the waveform of parameter Gaussian particle filter when getting different value.As can be seen from the figure parameter value can affect width and the waveform shape of shaped pulse, can select suitable parameter value according to actual needs in practical application, to meet different measurement demands.
Although the embodiment disclosed by the present invention is as above, the embodiment that described content just adopts for the ease of understanding the present invention, and be not used to limit the present invention.Technician in any the technical field of the invention; under the prerequisite not departing from the spirit and scope disclosed by the present invention; any amendment and change can be done what implement in form and in details; but scope of patent protection of the present invention, the scope that still must define with appending claims is as the criterion.

Claims (4)

1., based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure, it is characterized in that, described method comprises:
According to the circuit of simulation Gauss formation system, obtain the differential equation of circuit system input signal and output signal, the differential equation is derived by frequency-domain analysis, obtain the unit impulse response of simulation Gaussian Systems;
The frequency response of analysis mode Gauss formation system, and frequency that simulation Gauss formation system is sampled is determined according to the spectral range of frequency response;
Simulation Gauss formation system Impulse invariance procedure is converted to numeric field by analog domain, obtains the unit impulse response of digital Gaussian formation system, containing the forming parameter that can adjust waveform widths in described impulse response;
Forming parameter in selected impulse response, and the unit impulse response of corresponding digital Gaussian formation system is determined according to selected forming parameter value;
The unit impulse response of the core pulse signal and digital Gaussian formation system that are input to digital Gaussian formation system is carried out Convolution sums computing, realizes digital core pulse Gauss and be shaped.
2. the digital core pulse Gauss manufacturing process based on Impulse invariance procedure according to claim 1, is characterized in that, described frequency of sampling to simulation Gauss formation system is greater than 2 times of simulation Gauss formation system frequency response bandwidth.
3. the digital core pulse Gauss manufacturing process based on Impulse invariance procedure according to claim 1, it is characterized in that, described forming parameter affects width and the waveform shape of shaped pulse, and parameter value is less, and waveform is wider.
4. the digital core pulse Gauss manufacturing process based on Impulse invariance procedure according to claim 1, it is characterized in that, described Impulse invariance procedure is: make the unit impulse response of digital Gaussian formation system equal to simulate the sampled value of Gauss's formation system unit impulse response, thus obtain the unit impulse response of digital Gaussian formation system.
CN201310629796.2A 2013-11-29 2013-11-29 Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure Expired - Fee Related CN103941280B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310629796.2A CN103941280B (en) 2013-11-29 2013-11-29 Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310629796.2A CN103941280B (en) 2013-11-29 2013-11-29 Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure

Publications (2)

Publication Number Publication Date
CN103941280A CN103941280A (en) 2014-07-23
CN103941280B true CN103941280B (en) 2016-01-20

Family

ID=51189026

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310629796.2A Expired - Fee Related CN103941280B (en) 2013-11-29 2013-11-29 Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure

Country Status (1)

Country Link
CN (1) CN103941280B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105301627B (en) * 2015-11-23 2018-10-09 重庆大学 A kind of energy spectrum analysis method, energy spectrum analysis system and gamma-ray detection system
CN112327347B (en) * 2020-10-29 2022-11-18 中广核久源(成都)科技有限公司 Digital nuclear pulse forming system with adjustable curvature
CN114252899B (en) * 2022-03-02 2022-05-20 四川新先达测控技术有限公司 Cascade impulse convolution forming method and device for kernel signal

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201251629Y (en) * 2008-07-14 2009-06-03 成都理工大学 Digital random pulse real-time discriminating and storing system
US8315836B2 (en) * 2009-12-14 2012-11-20 Mitsubishi Electric Corporation Radiation measuring device and diagnostic method thereof

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201251629Y (en) * 2008-07-14 2009-06-03 成都理工大学 Digital random pulse real-time discriminating and storing system
US8315836B2 (en) * 2009-12-14 2012-11-20 Mitsubishi Electric Corporation Radiation measuring device and diagnostic method thereof

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Study of time-domain digital pulse shaping algorithms for nuclear signals;ZHOU Jianbin et al.;《Nuclear Science and Techniques》;20121231;第23卷;150-155 *
基于高斯函数分布的数字化能谱模拟软件在X荧光分析中的应用;李飞等;《核技术》;20110930;第34卷(第9期);663-665 *
无限冲激响应数字滤波器的设计与 MATLAB 仿真;赵普渡;《工程控制计算机》;20111231;第24卷(第10期);97-98、100 *
模拟滤波器的数字化;赵伟刚,陈洪亮;《技术交流》;20011231;第17卷(第12期);46-47、36 *

Also Published As

Publication number Publication date
CN103941280A (en) 2014-07-23

Similar Documents

Publication Publication Date Title
CN103837884B (en) Based on the digital core pulse signal trapezoidal shaping algorithm of time-domain analysis
CN104181577B (en) Beam position and phase measurement system and method based on full digitalization technology
CN103257271A (en) Device and method for detecting micro grid harmonic wave and inter-harmonics based on STM32F107VCT6
CN103941280B (en) Based on the digital core pulse Gauss manufacturing process of Impulse invariance procedure
CN103888104B (en) Method and system for designing FIR digital filter
CN103675891B (en) Based on the digital core pulse Gauss manufacturing process of Bilinear transformation method
CN103777228B (en) Based on the digital core pulse signal Gauss manufacturing process of iir filter
Franck Efficient algorithms for arbitrary sample rate conversion with application to wave field synthesis
CN107276560A (en) A kind of FRI sparse samplings kernel function construction method and circuit
CN109946512A (en) A kind of dynamic power analysis method for improving frequency domain interpolation
CN103853930A (en) Numerical simulation method and device for earthquake vector wave field
CN104502698A (en) Method and system for measuring frequency of electric power signal
CN105373708B (en) A kind of Time-Frequency Analysis Method of the improvement generalized S-transform based on parameter optimization
CN106483563A (en) seismic energy compensation method based on complementary set empirical mode decomposition
Liu et al. A resampling method based on filter designed by window function considering frequency aliasing
CN104242877B (en) Core pulse signal digital Gaussian manufacturing process based on simulation CR RC circuits
CN105334381A (en) Method and device for measuring AC active power
Saini et al. Power spectral density analysis of speech signal using window techniques
Nilsen Recursive time-frequency reassignment
CN108983321B (en) Method for extracting periodic components of solar black number and geomagnetic Ap index based on synchronous compression wavelet transform
CN103777221B (en) Based on the digital core pulse signal Gauss manufacturing process of window function metht
CN108427032B (en) Frequency spectrum decomposition method and frequency-time inversion method
CN114252899B (en) Cascade impulse convolution forming method and device for kernel signal
CN108984472A (en) A kind of stable method for solving of Time domain electric field integral equation
CN103941279A (en) Nuclear pulse signal digital Gaussian forming method based on double parameters of pulse width and amplitude width

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160120

Termination date: 20161129