CN106483550B - A kind of simulation spectrum curve emulation mode - Google Patents

A kind of simulation spectrum curve emulation mode Download PDF

Info

Publication number
CN106483550B
CN106483550B CN201510539124.1A CN201510539124A CN106483550B CN 106483550 B CN106483550 B CN 106483550B CN 201510539124 A CN201510539124 A CN 201510539124A CN 106483550 B CN106483550 B CN 106483550B
Authority
CN
China
Prior art keywords
spectrum curve
nuclear
practical
simulation
curve
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.)
Active
Application number
CN201510539124.1A
Other languages
Chinese (zh)
Other versions
CN106483550A (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.)
Array Microelectronics Ltd
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 CN201510539124.1A priority Critical patent/CN106483550B/en
Priority to CN201811086619.3A priority patent/CN109239763B/en
Priority to CN201811085398.8A priority patent/CN109271707B/en
Priority to CN201811085406.9A priority patent/CN109086555B/en
Publication of CN106483550A publication Critical patent/CN106483550A/en
Application granted granted Critical
Publication of CN106483550B publication Critical patent/CN106483550B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Complex Calculations (AREA)
  • Measurement Of Radiation (AREA)

Abstract

The present invention relates to a kind of simulation spectrum curve emulation modes, include the following steps:Step S1 obtains practical nuclear spectrum curve graph;Step S2 handles to obtain simulation spectrum curve practical nuclear spectrum curve graph;The simulation spectrum curve emulation mode of the present invention is by carrying out Curves Recognition to practical nuclear spectrum curve graph and spectrum curve quantizing, to obtain the numerical value (i.e. the counting rate of the energy level of nuclear spectrum and each energy level) of spectrum curve each point, again by this group of numerical value of the random direct sampling of Monte Carlo method to obtain the random number about each nuclear level, to simulate the randomness of nuclear decay process, statistical disposition finally is carried out to the random number again and obtains simulation spectrum curve, by inverting comparative simulation spectrum curve and practical spectrum curve to determine the reliability and accuracy of imitative nuclear signal generator.

Description

A kind of simulation spectrum curve emulation mode
Technical field
The present invention relates to a kind of nuclear energy field more particularly to a kind of simulation spectrum curve emulation modes.
Background technology
Nuclear decay process occurs at random in time, and it is also random, but right that decay process, which discharges ray (energy), Its time interval occurred and energy value take statistics analysis, it can be seen that nuclear decay process has following characteristic:Nuclear decay is occurring Time interval on approximate obey exponential distribution;The approximate Gauss that obeys of energy (i.e. power spectrum) that nuclear decay process externally discharges divides Cloth.
Based on nuclear decay process, there are the above characteristics, and the imitative nuclear signal generator of traditional approach is to obey different distributions Random number simulates nuclear signal feature, i.e., to obey exponential distribution random number simulation nuclear signal in time interval statistical property;With Random numbers of Gaussian distribution simulates statistical property of the nuclear signal in amplitude.But core is only simulated with the random number of Gaussian distributed Statistical property of the signal in amplitude be it is inapt, the high phase distribution curve simulated and practical spectrum curve there is compared with Big error cannot accurately reflect nuclear signal characteristic;At the same time, the power spectrum statistical property of each nucleic is not quite similar, therefore Then need the random numbers of Gaussian distribution for generating different parameters matching the amplitude characteristic simulation of variety classes nucleic, this is not It is actual, it is difficult to realize in the operating process of border.
In view of various drawbacks existing for traditional imitative nuclear signal generator, set forth herein a kind of completely new approach to Solve problem above.
Invention content
The object of the present invention is to provide a kind of simulation spectrum curve emulation modes, are emulated to nuclear energy spectral line with realizing.
In order to solve the above technical problem, the present invention provides a kind of simulation spectrum curve emulation modes, which is characterized in that Include the following steps:Step S1 obtains practical nuclear spectrum curve graph;And step S2, at practical nuclear spectrum curve graph Reason is to obtain simulation spectrum curve.
Further, the simulation spectrum curve emulation mode further includes:Step S3 compares simulation spectrum curve by inverting With practical spectrum curve, the error between spectrum curve and practical spectrum curve is simulated with acquisition.
Further, practical nuclear spectrum curve graph is handled to obtain the method packet of simulation spectrum curve in step S2 It includes:Step S21 carries out Curves Recognition to practical nuclear spectrum curve graph and spectrum curve quantizes, each to obtain spectrum curve The numerical value of point;Step S22 is obtained by Monte Carlo method this group of numerical value of random direct sampling about the random of each nuclear level Number, to simulate the randomness of nuclear decay process;And step S23, statistical disposition is carried out to the random number and obtains the simulation Spectrum curve.
Further, Curves Recognition is carried out to practical nuclear spectrum curve graph in the step S21 and spectrum curve quantizes Include in the method for obtaining the numerical value of spectrum curve each point:Practical nuclear energy is respectively set a song to music after the filtered denoising of image of line chart Practical spectrum curve figure is shown again, and copies each key of the nuclear energy spectral curve according to the practical spectrum curve figure of display Point obtains spectrum curve data, to establish spectrum curve database;Or by practical nuclear energy respectively set a song to music line chart image it is filtered, drop Pretreatment, Curves Recognition, curvilinear characteristic extraction, and progress interpolation processing make an uproar to improve and repair each point of missing spectrum curve Data, to establish spectrum curve database.
Further, it is obtained about each nuclear energy by Monte Carlo method this group of numerical value of random direct sampling in step S22 Grade random number include in the method for simulating the randomness of nuclear decay process:Nuclear signal time statistical property is simulated;And core letter Number amplitude statistics simulated behavior.
Further, the method for the nuclear signal time statistical property simulation includes:By the random number for obeying exponential distribution Realize the simulation of nuclear signal time statistical property, wherein
The random number of exponential distribution by (0,1] equally distributed random number converts to obtain by inverse function method, and (0,1] Even distribution random numbers are suitable for acquiring by linear congruential method.
Further, the method for the nuclear signal amplitude statistics simulated behavior includes:
By to practical nuclear spectrum Curves Recognition and digitizing and obtaining each energy level amplitude and counting rate, then pass through Monte Carlo Method direct sampling simultaneously exports the random number;Wherein
Identification and digitized process to practical nuclear energy spectral curve include:
Step S221 is filtered practical spectrum curve figure, noise reduction;
Step S222 finds out threshold value by split plot design between maximum kind, and by spectrum curve figure progress binary conversion treatment, then by Pixel scanning method extracts the numerical value i.e. coordinate of each point on spectrum curve;
Step S223, repairs spectrum curve and is quantized.
Further, described by Monte Carlo method direct sampling and the method that exports the random number, i.e.,
It is a series of random to obtain by each point value on Monte Carlo method direct sampling spectrum curve and curve Number, to simulate the randomness of nuclear decay process.
Further, the method practical spectrum curve figure being filtered in the step S221, i.e., to practical spectrum curve Figure carries out Wiener filtering processing, to filter out the Gaussian noise in spectrum curve figure.
Further, the method repaired and quantized to spectrum curve in the step S223 includes:Pass through sample three times Interpolation method is stretched by the expansion of the ratio of coordinate to obtain to fill up the data point lacked during spectrum curve feature extraction Obtain the numerical value of each point on spectrum curve figure.
The invention has the advantages that the simulation spectrum curve emulation mode of the present invention passes through to practical nuclear spectrum curve graph Carry out Curves Recognition simultaneously spectrum curve is quantized, to obtain spectrum curve each point numerical value (i.e. the energy level of nuclear spectrum and respectively The counting rate of energy level), then by this group of numerical value of the random direct sampling of Monte Carlo method with obtain about each nuclear level with Machine number finally carries out statistical disposition to the random number again and obtains simulating and can set a song to music to simulate the randomness of nuclear decay process Line, by inverting comparative simulation spectrum curve and practical spectrum curve to determine the reliability and accurately of imitative nuclear signal generator Property.
Description of the drawings
Present invention will be further explained below with reference to the attached drawings and examples.
Fig. 1 is the functional block diagram of the simulation spectrum curve emulation mode of the present invention;
Fig. 2 is the simulation spectrum curve emulation mode flow chart of the present invention;
Fig. 3 be the present invention the step S2 in practical nuclear spectrum curve graph is handled to obtain simulation spectrum curve Method flow diagram;
Fig. 4 is generation n=10000 according to the present invention (0,1) uniform random number distribution maps;
Fig. 5 is the exponential distribution random-number distribution figure of the present invention;
Fig. 6 is being evenly dividing 1000 group squares to the above exponential distribution random number value range and counting for the present invention Statistical chart;
Fig. 7 is the spectrum curve characteristic pattern of the present invention extracted;
Fig. 8 is the spectrum curve Preliminary Simulation design sketch of the present invention;
Fig. 9 is the design sketch after the cubic spline interpolation of the present invention;
Figure 10 is that the spectrum curve that finally obtains of the present invention simulates design sketch;
Figure 11 shows the design sketch of the simulation random generating process of nuclear signal;
Figure 12 is shown using Monte Carlo method direct sampling final effect figure.
Specific implementation mode
In conjunction with the accompanying drawings, the present invention is further explained in detail.These attached drawings are simplified schematic diagram, only with Illustration illustrates the basic structure of the present invention, therefore it only shows the composition relevant to the invention.
As shown in Figure 1, the simulation spectrum curve emulation mode of the present invention to practical nuclear spectrum curve graph by carrying out curve It identifies and spectrum curve quantizes, to obtain numerical value (the i.e. meter of the energy level of nuclear spectrum and each energy level of spectrum curve each point Digit rate), then by Monte Carlo method this group of numerical value of random direct sampling to obtain the random number about each nuclear level, to The randomness of nuclear decay process is simulated, finally carrying out statistical disposition to the random number again obtains simulation spectrum curve, by anti- Comparative simulation spectrum curve is drilled with practical spectrum curve to determine the reliability and accuracy of imitative nuclear signal generator.
The specific implementation mode of the present invention is as illustrated in the examples below.
As shown in Fig. 2, a kind of simulation spectrum curve emulation mode of the present invention, includes the following steps:
Step S1 obtains practical nuclear spectrum curve graph;
Step S2 handles to obtain simulation spectrum curve practical nuclear spectrum curve graph.
Optionally, the simulation spectrum curve emulation mode further includes:
Step S3 compares simulation spectrum curve and practical spectrum curve by inverting, to obtain simulation spectrum curve and reality Error between the spectrum curve of border.
Further, as shown in figure 3, being handled practical nuclear spectrum curve graph to obtain simulation power spectrum in the step S2 The method of curve includes:
Step S21 carries out Curves Recognition to practical nuclear spectrum curve graph and spectrum curve quantizes, to obtain to set a song to music The numerical value of line each point;Step S22 is obtained by Monte Carlo method this group of numerical value of random direct sampling about each nuclear level Random number, to simulate the randomness of nuclear decay process;And step S23, random number progress statistical disposition is obtained described Simulate spectrum curve.
Specifically, carrying out Curves Recognition to practical nuclear spectrum curve graph in the step S21 and spectrum curve quantizing Include in the method for obtaining the numerical value of spectrum curve each point:
Practical nuclear energy is respectively set a song to music and is again shown practical spectrum curve figure after the filtered denoising of image of line chart, And each key point that the nuclear energy spectral curve is copied according to the practical spectrum curve figure of display obtains spectrum curve data, to establish energy Spectral curve database;Or by practical nuclear energy respectively set a song to music line chart image is filtered, noise reduction pretreatment, Curves Recognition, curvilinear characteristic Extraction, and interpolation processing is carried out to improve and repair each point data of missing spectrum curve, to establish spectrum curve database.
Wherein, it is obtained about each core by Monte Carlo method this group of numerical value of random direct sampling in the step S22 The random number of energy level includes in the method for simulating the randomness of nuclear decay process:Nuclear signal time statistical property is simulated and nuclear signal Amplitude statistics simulated behavior.
The method of nuclear signal time statistical property simulation includes:Random number by obeying exponential distribution realizes core letter The simulation of number time statistical property, the wherein random number of exponential distribution by (0,1] equally distributed random number become by inverse function method Get in return, and (0,1] uniform random number be suitable for acquired by linear congruential method.
Specifically, acquired by linear congruential method (0,1] method of uniform random number is as follows:
The recurrence formula of linear congruential method is as follows:
xi+1≡λxi+c(mod M) (1)
Wherein λ, c are constant.The initial x chosen1Referred to as seed has certain influence, value to the generation quality of random number Respectively 1~216It is chosen between=65535.For the ease of using on computers, usually take
M=2S, wherein S is binary maximum possible number of significant digit in computer.
Fig. 4 be take 10000 (0,1] random-number distribution situation
The production method of exponential distribution random number, i.e. exponential distribution random number can be realized that detailed process is such as by inverse function method Under:
If the distribution function of stochastic variable X obeys exponential distribution:
F (x)=1-e-ax, x >=0 (3)
Wherein, a is a time constant, and e is the nature truth of a matter.
Can be with by above formula, F (x) ∈ [0,1), and the monotone decreasing in domain, therefore function F (x) must have in 0~+∞ Inverse function seeks its inverse function:
Due to 0 < 1-F (x)≤1, above formula can be reduced to
By formula (5) it can be seen that by meet (0,1] equally distributed Random sampling obtains obeying exponential distribution random number x.
It takesThe exponential distribution random-number distribution generated by inverse function method by the equal distribution random numbers of the above unit Figure such as Fig. 5.1000 group squares are evenly dividing to the above exponential distribution random number value range and are counted, final statistical chart is such as Shown in Fig. 6.
The method of the nuclear signal amplitude statistics simulated behavior includes:By to practical nuclear spectrum Curves Recognition and digitizing Each energy level amplitude and counting rate are obtained, then by Monte Carlo method direct sampling and exports the random number;Wherein to reality The identification of nuclear energy spectral curve and digitized process includes:
Step S221 is filtered practical spectrum curve figure, noise reduction;Step S222 is asked by split plot design between maximum kind Go out threshold value, and spectrum curve figure is subjected to binary conversion treatment, then extracts the number of each point on spectrum curve by pixel scanning method Value is coordinate;Step S223, repairs spectrum curve and is quantized.
Specifically, described by Monte Carlo method direct sampling and the method that exports the random number, that is, pass through cover it is special A series of each point value on Caro method direct sampling spectrum curve and curve, to obtain random random numbers, to mould The randomness of nucleoid decay process.
The method being filtered to practical spectrum curve figure in the step S221 is tieed up practical spectrum curve figure It receives and is filtered, to filter out the Gaussian noise in spectrum curve figure, to reduce the interference of grass as far as possible.
The specific implementation process of the method for the nuclear signal amplitude statistics simulated behavior is as follows:
Practical spectrum curve figure is filtered in the step S221, the specific implementation step of noise reduction it is as follows:
Practical spectrum curve figure is filtered by Wiener filtering, noise reduction process, i.e., the described Wiener filter is one kind Linear filter, and a kind of still optimal estimation device based on minimum mean square error criterion, to stationary process.
Assuming that Wiener filter input signal is s (t), superimposed noise n (t).Output signal x (t) passes through filter g (t) It is obtained by following convolution algorithm:
X (t)=g (t) * (s (t)+n (t)) (6)
For the signal x (t) estimated, it is expected that being equal to s (t).
Its error is:E (t)=s (t+d)-x (t) (7)
Variance is:e2(t)=s2(t+d)-2s(t+d)x(t)+x2(t) (8)
Wherein s (t+d) is desired filter output.
Write x (t) as convolution integrals, i.e.,
Square error, which can be calculated, is:
Wherein RsIt is the auto-correlation function of s (t), RxIt is the auto-correlation function of x (t), RxsIt is the auto-correlation of x (t) and s (t) Function.The final purpose of Wiener filtering is exactly to seek optimal g (t) so that E (e2) minimum.
Threshold value is found out by split plot design between maximum kind in the step S222, and spectrum curve figure is carried out at binaryzation It manages, then extracts the numerical value i.e. coordinate of each point on spectrum curve by pixel scanning method;
The specific algorithm process of maximum variance between clusters is as follows:
If the gray value of piece image is 1~m, the pixel number that wherein gray value is i is ni, N expression image pixels Point sum, then the probability that gray value is i appearance is:
It is c to enable gray value be more than threshold value k1Group, i.e. c1={ 1~k }, gray value are then c more than threshold value k's2Group, c2={ k+ 1~m }, then C1And C2The probability of appearance is respectively:
C is calculated1And C2Gray average be:
μr1·μ12·μ2 (16)
Thus the variances sigma between two groups can be calculated2For:σ2(k)=ω11r)222r)2 (17)
Formula (16), which is substituted into formula (17), to be obtained:σ2(k)=ω2ω221)2
So optimal threshold T*=Arg max { σ2(k) }, (18) 0≤k < m-1
Acquire segmentation threshold T*=0.6353.
Spectrum curve is repaired and is quantized in the step S223 and is as follows:
It is to extract the numerical value of each point on nuclear energy spectral curve to sit after the filtered denoising of practical nuclear spectrum curve graph, binaryzation Mark, need to extract spectrum curve feature, and curve is quantized.Detailed process is as follows:
First, Straight Line Identification identifies straight in nuclear energy spectrogram that is, by scanning the row and column of nuclear energy spectral curve binary map Line;
Secondly, fixed point, cross, the ordinate of coordinate system where judging spectrum curve by the straight line identified, and origin is positioned, Generally from top to bottom, it scans from left to right, the straight line identified is just horizontal, ordinate;
Third, spectrum curve feature extraction.To reduce the influence of frame and coordinate pair curve in image, frame need to be filtered It removes.Filter out after frame again by pixel scan method line by line or scan by column point that pixel is 0 (black is 0 in bianry image, 1) white is.
Finally, curve quantizes.After extracting curve, by calculating the spectrum curve available point scanned to scan origin Row and stringer distance determine position of the pixel in figure, be somebody's turn to do finally by the scale factor for expanding coordinate is multiplied by The coordinate value of pixel.
Final extraction spectrum curve characteristic effect is as shown in Figure 7.
And spectrum curve Preliminary Simulation effect is as shown in Figure 8.
Further, it can be seen from Fig. 7 and Fig. 8, the spectrum curve figure of obtained simulation is compared with proper energy spectral curve in certain points Data lack.Truly to reflect practical spectrum curve characteristic as far as possible, need that the data of missing fill up repairing It is multiple.
Specifically, described fill up the number lacked during spectrum curve feature extraction by cubic spline interpolation Strong point, and stretched by the expansion of the ratio of coordinate to obtain the numerical value of each point on spectrum curve figure, it is effective to the data of missing to realize It fills up and repairs in ground.
The cubic spline interpolation fills up the specific of the data point lacked during spectrum curve feature extraction Algorithm is as follows:
Piecewise function S (x) on interval of definition [a, b], if meeting:
1. S (x) is in each subinterval [xi, xi+1] on be a cubic polynomial function;
2. S (x) has continuous second dervative on entire section [a, b].
Then S (x) is referred to as on section [a, b] about a=x0< x1< ... < xnA cubic spline function of=b.To three Secondary spline interpolation problem is:The n+1 node x of given function g (x)0, x1..., xnObtain function y0, y1..., yn, ask one three times Spline function S (x) makes its satisfaction:
S(xj)=yj, j=0,1 ..., n (19)
Wherein, function S (x) is known as the cubic spline functions of g (x).
If S (x) is the sample spline interpolation function three times of f (x), the following conditions are must satisfy:
1. interpolation condition, i.e.,
S(xj)=yj, j=0,1 ..., n-1
2. the condition of continuity, i.e.,
3. the first derivative condition of continuity, i.e.,
4. the second dervative condition of continuity, i.e.,
By the design sketch after cubic spline interpolation, as shown in figure 9, can see by its partial enlarged view, cubic spline Data point after interpolation is more smooth, more approaches actual value.
Practical nuclear spectrum curve simulation effect finally obtains practical nuclear spectrum curve graph after above-mentioned image procossing It is as shown in Figure 10 that the spectrum curve simulates design sketch.
Specifically, by each point value on Monte Carlo method direct sampling spectrum curve and curve, to obtain a system Random random number is arranged, to simulate the randomness of nuclear decay process.
Figure 11 shows the design sketch of the simulation random generating process of nuclear signal;
Figure 12 shows that (this figure is by practical spectrum curve figure numerical value using Monte Carlo method direct sampling final effect figure After obtaining energy level and counting rate this array after change, what random sampling procedure and counting obtained.This figure is suitable for proving by Meng Teka The reasonability and accuracy of sieve sampling).
Numerical value (the abscissa of each point on simulation spectrum curve and curve has been obtained by numbers above image processing process For road location Channel, ordinate is counting rate Count), then can be obtained with this group of data of Monte Carlo method direct sampling (energy level is such as, but not limited to quantified by multichannel analyzer a series of random energy levels, and the road location refers to nuclear decay What the energy of process release obtained after multichannel analyzer quantifies) random number, to simulate the randomness of nuclear decay process.Most The random number is counted again afterwards, simulation spectrum curve figure can be obtained, on the one hand can verify the reliable of system in this way Property and accuracy, on the other hand can also inverting in multichannel analyzer, to demarcate the accuracy of multichannel analyzer.
Simulation calculating is carried out to probability P (A)=p (unknown) that certain event A occurs using Monte Carlo method, it is specific to count Calculation method:
(1) n times are carried out and repeat independent sampling experiment, calculating event A frequencies are nA
Introduce stochastic variable Xi, indicate event A frequencies in ith experiment, enable
(2) occurrence frequency fs of the event A in n times repeat independent sampling experiment is calculatedN, it is
(3) when N is fully big, with probability fNEstimated value as probability P (A)=pFor
(4) estimated value is requiredFor the unbiased esti-mator of probability P (A)=p, i.e.,
And direct sampling, the i.e. characteristic to nuclear signal on time and amplitude are with the random of two groups of obedience different distributions Number is simulated, and random number be it is discrete, it is discontinuous.Sampling for discrete random sequence, the non-convention of direct sampling method Think.
The specific sampling process of discrete distribution direct sampling method is as follows:
If the value range of discrete random variable X is Xi(i=0,1,2,3 ...), probability distribution is
P (X=Xi)=Pi(i=0,1,2,3 ...).Wherein Pi>=0,
(1) equally distributed random number r on (0,1) section is generated;
(2) positive integer n=0,1,2... are acquired so that r meets
(3) sample value for extracting discrete random variable X is X=Xn.And when 0<r≤P0When X=X0
(4) step (1), (2), (3) are repeated until extracting n sample value.
Due to generation (0,1) if equally distributed random number r is in sectionProbability be
That is eventThe probability of appearance is equal to event X=XnThe probability of generation.
Again because random number r obeys being uniformly distributed on (0,1), probability density function is
Its distribution function is as follows:
Therefore it is X=x that the random number r generated, which draws sample value,nProbability be
It follows that being drawn into (X=X by direct sampling methodn) probability be equivalent to random number XnIn random number sequence X1, X2... XnThe frequency of appearance.
It can be by following proof for direct sampling method reliability:
If X is discrete random variable, probability distribution Pi=P { X=Xi, wherein i=1,2 ....X is respectively with PiIt takes Obtain Xi, thenEvent | X-E (X) | >=ε indicates that stochastic variable X acquirements are all and meets inequality | Xi- E (X) | >=ε's can It can value Xi, then
Due to event X=XiThe probability that (i=0,1,2 ... N) occurs is pi(0<pi<1), then X ≠ XiProbability then be 1- pi, and each X=XiThe probability of generation is constant, and sampling results are unrelated with other each times extraction results every time.Therefore X=Xi Individual event is a bernoulli experiment, then n times of sampling, then be n again Bernoulli trials.If enabling event A (X=Xi) occur Number is nA, i.e. nA~B (n, p).Due to X1, X2..., XnIt is n independently of each other and obeys the random of the 0-1 distributions that parameter is p Variable, and
HaveD (Xi)=p (1-p), i=1,2 ... n.Giving arbitrary ε > 0 then has
It can be derived from by (4.31) formula
Therefore it can be derived from
Abbreviation obtains
I.e. when the frequency n of extraction is bigger, the frequency ratio of event A occurs after sampling number and population of samples is closer to thing The probability that part A occurs.
By direct sampling extraction random number, its error is:
It enablesTherefore
I.e.It is the unbiased esti-mator of p,
That is the number of sampling n is bigger, estimated valueCloser to theoretical value p.
It is enlightenment with above-mentioned desirable embodiment according to the present invention, through the above description, relevant staff is complete Various changes and amendments can be carried out without departing from the scope of the technological thought of the present invention' entirely.The technology of this invention Property range is not limited to the contents of the specification, it is necessary to determine its technical scope according to right.

Claims (1)

1. a kind of simulation spectrum curve emulation mode, which is characterized in that include the following steps:
Step S1 obtains practical nuclear spectrum curve graph;
Step S2 handles to obtain simulation spectrum curve practical nuclear spectrum curve graph;
The simulation spectrum curve emulation mode further includes:
Step S3 compares simulation spectrum curve and practical spectrum curve by inverting, to obtain simulation spectrum curve and practical energy Error between spectral curve;
Wherein, practical nuclear spectrum curve graph is handled in step S2 and includes in the method for obtaining simulation spectrum curve:
Step S21 carries out Curves Recognition to practical nuclear spectrum curve graph and spectrum curve quantizes, each to obtain spectrum curve The numerical value of point;
Step S22 obtains the random number about each nuclear level by Monte Carlo method this group of numerical value of random direct sampling, with Simulate the randomness of nuclear decay process;
Step S23 carries out statistical disposition to the random number and obtains the simulation spectrum curve;
Wherein, Curves Recognition is carried out to practical nuclear spectrum curve graph in the step S21 and spectrum curve quantizes to obtain The method of the numerical value of spectrum curve each point, including:
Practical nuclear energy is respectively set a song to music and is again shown practical spectrum curve figure after the filtered denoising of image of line chart, and root Each key point that the nuclear energy spectral curve is copied according to the practical spectrum curve figure of display obtains spectrum curve data, can be set a song to music with establishing Line database;
Wherein, in step S22 by Monte Carlo method this group of numerical value of random direct sampling obtain about each nuclear level with Method of the machine number to simulate the randomness of nuclear decay process, including:
Nuclear signal time statistical property is simulated;And nuclear signal amplitude statistics simulated behavior.
CN201510539124.1A 2015-08-28 2015-08-28 A kind of simulation spectrum curve emulation mode Active CN106483550B (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201510539124.1A CN106483550B (en) 2015-08-28 2015-08-28 A kind of simulation spectrum curve emulation mode
CN201811086619.3A CN109239763B (en) 2015-08-28 2015-08-28 Simulation method for simulating nuclear decay process by using simulated energy spectrum curve
CN201811085398.8A CN109271707B (en) 2015-08-28 2015-08-28 Simulation energy spectrum curve simulation method for simulating nuclear energy spectrum line
CN201811085406.9A CN109086555B (en) 2015-08-28 2015-08-28 Simulation method for simulating energy spectrum curve by adopting Monte Carlo method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510539124.1A CN106483550B (en) 2015-08-28 2015-08-28 A kind of simulation spectrum curve emulation mode

Related Child Applications (3)

Application Number Title Priority Date Filing Date
CN201811086619.3A Division CN109239763B (en) 2015-08-28 2015-08-28 Simulation method for simulating nuclear decay process by using simulated energy spectrum curve
CN201811085406.9A Division CN109086555B (en) 2015-08-28 2015-08-28 Simulation method for simulating energy spectrum curve by adopting Monte Carlo method
CN201811085398.8A Division CN109271707B (en) 2015-08-28 2015-08-28 Simulation energy spectrum curve simulation method for simulating nuclear energy spectrum line

Publications (2)

Publication Number Publication Date
CN106483550A CN106483550A (en) 2017-03-08
CN106483550B true CN106483550B (en) 2018-09-18

Family

ID=58234736

Family Applications (4)

Application Number Title Priority Date Filing Date
CN201811085406.9A Active CN109086555B (en) 2015-08-28 2015-08-28 Simulation method for simulating energy spectrum curve by adopting Monte Carlo method
CN201811086619.3A Active CN109239763B (en) 2015-08-28 2015-08-28 Simulation method for simulating nuclear decay process by using simulated energy spectrum curve
CN201811085398.8A Active CN109271707B (en) 2015-08-28 2015-08-28 Simulation energy spectrum curve simulation method for simulating nuclear energy spectrum line
CN201510539124.1A Active CN106483550B (en) 2015-08-28 2015-08-28 A kind of simulation spectrum curve emulation mode

Family Applications Before (3)

Application Number Title Priority Date Filing Date
CN201811085406.9A Active CN109086555B (en) 2015-08-28 2015-08-28 Simulation method for simulating energy spectrum curve by adopting Monte Carlo method
CN201811086619.3A Active CN109239763B (en) 2015-08-28 2015-08-28 Simulation method for simulating nuclear decay process by using simulated energy spectrum curve
CN201811085398.8A Active CN109271707B (en) 2015-08-28 2015-08-28 Simulation energy spectrum curve simulation method for simulating nuclear energy spectrum line

Country Status (1)

Country Link
CN (4) CN109086555B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416355B (en) * 2018-03-09 2021-07-30 浙江大学 Industrial field production data acquisition method based on machine vision
CN111276195B (en) * 2020-01-20 2023-09-19 深圳大榆树科技有限公司 Method for calculating number of gel energy spectrogram compounds by using maximum inter-class variance method
CN111273337B (en) * 2020-02-27 2022-11-29 成都工贸职业技术学院 Nuclear energy spectrum processing method based on Monte Carlo pulse interpolation method
CN111553111B (en) * 2020-04-30 2023-03-28 成都航空职业技术学院 Digital imitation nuclear signal generator based on MCNP
CN112462676B (en) * 2021-01-27 2021-05-07 泛华检测技术有限公司 Device capable of simulating overlapped nuclear pulse signal generation and control method thereof
CN112462675B (en) * 2021-01-27 2021-05-07 泛华检测技术有限公司 Cascaded dual-index nuclear pulse signal generating device and control method thereof
CN114422041B (en) * 2021-12-23 2023-03-24 中国原子能科学研究院 Nuclear signal simulation method, device, terminal and storage medium
CN114241846A (en) * 2021-12-24 2022-03-25 西安恒律模训科技发展有限公司 Simulation gamma radioactive nuclide identification training method and system thereof

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201233446Y (en) * 2008-07-14 2009-05-06 成都理工大学 Arbitrary nuclear power spectrum generator
CN102073060A (en) * 2009-11-24 2011-05-25 成都理工大学 Simulation method for random properties of nuclear signals
CN102353972A (en) * 2011-07-01 2012-02-15 成都理工大学 Multimode digital multichannel spectrometer
CN102902958A (en) * 2012-09-19 2013-01-30 四川大学 Digital nuclear signal processing and multi-parameter nuclear information acquisition method based on image recognition
CN104316954A (en) * 2014-09-28 2015-01-28 中国石油大学(华东) Nuclear physics experiment simulation system and method for carrying out energy spectrum measurement experiment and intensity measurement experiment by using the nuclear physics experiment simulation system

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4633088A (en) * 1985-04-08 1986-12-30 Packard Instrument Co., Inc. Reverse sum quench measurement using a liquid scintillation counter
AUPO427796A0 (en) * 1996-12-20 1997-01-23 University Of Queensland, The Computer simulation of magnetic resonance spectra employing homotopy
JP4309733B2 (en) * 2003-09-29 2009-08-05 株式会社東芝 Combustion calculation method and combustion calculation program
US7411188B2 (en) * 2005-07-11 2008-08-12 Revera Incorporated Method and system for non-destructive distribution profiling of an element in a film
US8698090B2 (en) * 2008-10-10 2014-04-15 Koninklijke Philips N.V. Practical SPECT calibration method for quantification of nuclides with high-energy contributions
CN201622351U (en) * 2009-11-18 2010-11-03 成都理工大学 Nuclear signal random characteristic simulator
WO2011123837A2 (en) * 2010-04-01 2011-10-06 University Of Georgia Research Foundation, Inc. Method and system using computer simulation for the quantitative analysis of glycan biosynthesis
CN102298652B (en) * 2010-06-23 2013-02-27 成都理工大学 Method for simulating energy spectrum drift during radioactive measurement
JP2012037305A (en) * 2010-08-05 2012-02-23 Fujita Corp Sequential nonlinear earthquake response analysis method for foundation and storage medium with analysis program stored thereon
DE102011055075B4 (en) * 2010-11-05 2013-04-18 Stefan Brühl A method and apparatus for the preliminary planning of a medical ion beam therapy and apparatus for performing a medical ion beam therapy
CN102928866B (en) * 2011-08-09 2015-05-20 中国辐射防护研究院 Method for measuring spectrum and accumulated dose of neutrons by utilizing passive detector
CN103091701B (en) * 2011-10-28 2015-09-30 中国原子能科学研究院 Multipurpose flight time equipment for measuring quality of cold neutron beam
CN102916683B (en) * 2012-10-18 2016-09-14 成都理工大学 A kind of Parameter-adjustablenuclear nuclear pulse simulation method
CN103076622B (en) * 2012-10-31 2016-08-17 成都理工大学 A kind of production method of spectrum stabilization stochastic signal
WO2014090297A1 (en) * 2012-12-12 2014-06-19 Qatar Foundation System and method for the simulation of metabolic profiles
US10593436B2 (en) * 2013-11-21 2020-03-17 Terrapower, Llc Method and system for generating a nuclear reactor core loading distribution
CN103853929B (en) * 2014-03-17 2016-06-15 东华理工大学 A kind of based on the low resolution gamma spectrum inversion analysis system and method covering card response matrix

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201233446Y (en) * 2008-07-14 2009-05-06 成都理工大学 Arbitrary nuclear power spectrum generator
CN102073060A (en) * 2009-11-24 2011-05-25 成都理工大学 Simulation method for random properties of nuclear signals
CN102353972A (en) * 2011-07-01 2012-02-15 成都理工大学 Multimode digital multichannel spectrometer
CN102353972B (en) * 2011-07-01 2013-04-10 成都理工大学 Multimode digital multichannel spectrometer
CN102902958A (en) * 2012-09-19 2013-01-30 四川大学 Digital nuclear signal processing and multi-parameter nuclear information acquisition method based on image recognition
CN104316954A (en) * 2014-09-28 2015-01-28 中国石油大学(华东) Nuclear physics experiment simulation system and method for carrying out energy spectrum measurement experiment and intensity measurement experiment by using the nuclear physics experiment simulation system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《任意分布的高速仿核信号发生器》;霍建文 等;《核电子学与探测技术》;20140331;第34卷(第3期);全文 *
《基于随机抽样的核脉冲信号发生器的研究》;谭承君 等;《原子能科学技术》;20141031;第48卷;全文 *
《数字核信号图形化识别和处理研究》;万东阳 等;《第十六届全国核电子学与核探测技术学术年会》;20121231;全文 *

Also Published As

Publication number Publication date
CN109086555A (en) 2018-12-25
CN109239763A (en) 2019-01-18
CN109239763B (en) 2022-04-19
CN109271707B (en) 2022-04-19
CN109271707A (en) 2019-01-25
CN106483550A (en) 2017-03-08
CN109086555B (en) 2022-04-19

Similar Documents

Publication Publication Date Title
CN106483550B (en) A kind of simulation spectrum curve emulation mode
CN109884891B (en) High-precision intelligent nuclear simulation signal generation system and working method thereof
CN106443379B (en) The recognition methods of partial discharge of transformer fault type and device
CN104134204B (en) Image definition evaluation method and image definition evaluation device based on sparse representation
CN106483551B (en) A kind of imitative nuclear signal generator and its working method
CN110879980B (en) Nuclear magnetic resonance spectrum denoising method based on neural network algorithm
Venema et al. A stochastic iterative amplitude adjusted Fourier transform algorithm with improved accuracy
CN109270345A (en) A kind of detection method of mains by harmonics signal
CN107729289B (en) Polynomial phase signal self-adaptive time-frequency transformation method based on genetic optimization
CN102184451B (en) Multi-wavelet fusion feature based neuron action feature extraction method
US10396911B1 (en) Noise analysis to reveal jitter and crosstalk&#39;s effect on signal integrity
CN113486750A (en) Oil-immersed transformer partial discharge signal denoising method based on improved VMD algorithm and wavelet packet denoising algorithm
CN112858782B (en) Harmonic responsibility quantification method and system under influence of multi-user interaction of power system
CN106483867B (en) The working method of imitative nuclear signal generating system processor
CN116150594B (en) Method for identifying switch element characteristics in spectrum test data
CN114118177B (en) Bit plane noise reduction method and system based on singular spectrum analysis and storage medium
CN106501293A (en) A kind of method that XRF power spectrum background is deducted based on Wavelet Transform
CN114021424B (en) PCA-CNN-LVQ-based voltage sag source identification method
CN115936196A (en) Monthly rainfall model prediction method based on time sequence convolution network
CN115374605A (en) Hydroelectric generating set runout data noise elimination method based on SVMD and permutation entropy
JP7151561B2 (en) Pulse train signal period estimation device, pulse train signal period estimation method, and pulse train signal period estimation program
CN107622035B (en) Polynomial phase signal self-adaptive time-frequency transformation method based on simulated annealing
CN112054803A (en) Communication signal sorting method based on compressed sensing
CN111008356A (en) WTSVD algorithm-based background-subtracted gamma energy spectrum set analysis method
CN112788057A (en) LDoS attack detection method based on FSWT time-frequency distribution

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20180821

Address after: 610000 1 East three road, Erxian bridge, Chengdu, Sichuan

Applicant after: Chengdu University of Technology

Address before: 213000 Hohai University, Xinbei District, Changzhou, Jiangsu

Applicant before: Yi Liangbi

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20240218

Address after: 201100 Room 201-202, Building E, No. 1618, Yishan Road, Minhang District, Shanghai

Patentee after: ARRAY MICROELECTRONICS LTD.

Country or region after: China

Address before: 610000 1 East three road, Erxian bridge, Chengdu, Sichuan

Patentee before: Chengdu University of Technology

Country or region before: China

TR01 Transfer of patent right