CN106483550A - A kind of simulation spectrum curve emulation mode - Google Patents
A kind of simulation spectrum curve emulation mode Download PDFInfo
- Publication number
- CN106483550A CN106483550A CN201510539124.1A CN201510539124A CN106483550A CN 106483550 A CN106483550 A CN 106483550A CN 201510539124 A CN201510539124 A CN 201510539124A CN 106483550 A CN106483550 A CN 106483550A
- Authority
- CN
- China
- Prior art keywords
- spectrum curve
- nuclear
- simulation
- actual
- random number
- 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
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 173
- 238000004088 simulation Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 claims abstract description 83
- 238000005070 sampling Methods 0.000 claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 35
- 230000003595 spectral effect Effects 0.000 claims description 12
- 230000008034 disappearance Effects 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000008439 repair process Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000000052 comparative effect Effects 0.000 abstract description 3
- 238000009826 distribution Methods 0.000 description 13
- 238000013461 design Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000005654 stationary process Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Radiation (AREA)
- Complex Calculations (AREA)
Abstract
The present invention relates to a kind of simulation spectrum curve emulation mode, comprise the steps：Step S1, obtains actual nuclear spectrum curve chart；Step S2, is processed to actual nuclear spectrum curve chart to obtain simulation spectrum curve；The simulation spectrum curve emulation mode of the present invention by carrying out Curves Recognition and spectrum curve quantizing to actual nuclear spectrum curve chart, thus obtaining 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 the random number with regard to each nuclear level is obtained by this group numerical value of the random direct sampling of DSMC, thus simulating the randomness of nuclear decay process, finally more described random number is carried out with statistical disposition obtain simulating spectrum curve, determine reliability and the accuracy of imitative nuclear signal generator by inverting comparative simulation spectrum curve and actual spectrum curve.
Description
Technical field
The present invention relates to a kind of nuclear energy field, more particularly, to a kind of simulation spectrum curve emulation mode.
Background technology
Nuclear decay process is random generation in time, and decay process release ray (energy) is also random
, the analysis but time interval that it is occurred and energy value take statistics, can learn that nuclear decay process has following
Characteristic：Approximate obedience exponential in the time interval that nuclear decay occurs；Nuclear decay process externally discharges
The approximate Gaussian distributed of energy (i.e. power spectrum).
Above characteristic is had based on nuclear decay process, the imitative nuclear signal generator of traditional approach is to obey difference
The random number of distribution simulating nuclear signal feature, that is, to obey exponential random number simulation nuclear signal in the time
Interval statistical characteristic；Statistical property in amplitude for the nuclear signal is simulated with random numbers of Gaussian distribution.But only to take
It is inapt for simulating statistical property on amplitude for the nuclear signal from the random number of Gauss distribution, and simulation obtains
High phase distribution curve and actual spectrum curve there is larger error it is impossible to reflection nuclear signal spy exactly
Property；Meanwhile, the power spectrum statistical property of each nucleic is not quite similar, therefore the amplitude to variety classes nucleic
Simulated behavior then need produce different parameters random numbers of Gaussian distribution matching, this is unpractiaca,
It is difficult in the operating process of border.
In view of a variety of drawbacks that traditional imitative nuclear signal generator exists, set forth herein one kind completely newly side
Method is in order to solve problem above.
Content of the invention
It is an object of the invention to provide a kind of simulation spectrum curve emulation mode, to realize nuclear energy spectral line is carried out
Emulation.
In order to solve above-mentioned technical problem, the invention provides a kind of simulation spectrum curve emulation mode, it is special
Levy and be, comprise the steps：Step S1, obtains actual nuclear spectrum curve chart；And step S2, to reality
Border nuclear spectrum curve chart is processed to obtain simulation spectrum curve.
Further, described simulation spectrum curve emulation mode also includes：Step S3, compares simulation by inverting
Spectrum curve and actual spectrum curve, to obtain the error between simulation spectrum curve and actual spectrum curve.
Further, in step S2, actual nuclear spectrum curve chart is processed to obtain the side of simulation spectrum curve
Method includes：Step S21, carries out Curves Recognition and spectrum curve quantizes to actual nuclear spectrum curve chart, with
Obtain the numerical value of spectrum curve each point；Step S22, by this group number of the random direct sampling of DSMC
It is worth to the random number with regard to each nuclear level, to simulate the randomness of nuclear decay process；And step S23, right
Described random number carries out statistical disposition and obtains described simulation spectrum curve.
Further, in described step S21, actual nuclear spectrum curve chart is carried out with Curves Recognition and by spectrum curve
The method of the numerical value to obtain spectrum curve each point of quantizing includes：Actual nuclear energy is respectively set a song to music the image of line chart
Again actual spectrum curve figure is shown after filtered denoising, and the actual spectrum curve according to display
Each key point that figure copies this nuclear energy spectral curve obtains spectrum curve data, to set up spectrum curve data base；
Or by actual nuclear energy respectively set a song to music line chart image is filtered, noise reduction pretreatment, Curves Recognition, curvilinear characteristic carry
Take, and carry out interpolation processing to improve and to repair each point data of disappearance spectrum curve, can be set a song to music with setting up
Line data base.
Further, in step S22 by this group numerical value of the random direct sampling of DSMC obtain with regard to
The random number of each nuclear level is included with the method simulating the randomness of nuclear decay process：Nuclear signal time statistics is special
Property simulation；And nuclear signal amplitude statistics simulated behavior.
Further, the method for described nuclear signal time statistical property simulation includes：By obeying exponential
Random number realizes the simulation of nuclear signal time statistical property, wherein
The random number of exponential by (0,1] equally distributed random number by inverse function method conversion obtain, and
(0,1] uniform random number is suitable to try to achieve by linear congruential method.
Further, the method for described nuclear signal amplitude statistics simulated behavior includes：
By to actual nuclear spectrum Curves Recognition and digitized obtains each energy level amplitude and counting rate, then by covering
Special Caro method direct sampling simultaneously exports described random number；Wherein
Identification digitized process to actual nuclear energy spectral curve include：
Step S221, is filtered to actual spectrum curve figure, noise reduction；
Step S222, by maximum kind between split-run obtain threshold value, and spectrum curve figure is carried out at binaryzation
Reason, then be coordinate by the numerical value that pixel scanning method extracts each point on spectrum curve；
Step S223, is repaired to spectrum curve and is quantized.
Further, described by DSMC direct sampling and the method that exports described random number, that is,
By each point value on DSMC direct sampling spectrum curve and curve, a series of to obtain
Random number, thus to simulate the randomness of nuclear decay process.
Further, the method in described step S221, actual spectrum curve figure being filtered, that is, to actual energy
Line chart of setting a song to music carries out Wiener filtering process, to filter the Gaussian noise of spectrum curve in figure.
Further, in described step S223, spectrum curve is repaired and the method that quantized is included：Pass through
Cubic spline interpolation filling up during spectrum curve feature extraction the data point of disappearance, and by sitting
Target ratio expands stretches to obtain the numerical value of each point on spectrum curve figure.
The invention has the beneficial effects as follows, the simulation spectrum curve emulation mode of the present invention is by actual nuclear spectrum
Curve chart carries out Curves Recognition and spectrum curve quantizes, thus obtaining the numerical value of spectrum curve each point (i.e.
The energy level of nuclear spectrum and the counting rate of each energy level), then pass through this group number of the random direct sampling of DSMC
Value is to obtain the random number with regard to each nuclear level, thus simulating the randomness of nuclear decay process, finally again to institute
State random number and carry out statistical disposition and obtain simulating spectrum curve, by inverting comparative simulation spectrum curve with actual
Spectrum curve is to determine reliability and the accuracy of imitative nuclear signal generator.
Brief description
The present invention is further described with reference to the accompanying drawings and examples.
Fig. 1 is the theory 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 described step S2 in actual nuclear spectrum curve chart is processed to obtain simulation energy
The method flow diagram of spectral curve；
Fig. 4 is generation n=10000 involved in the present invention (0,1) uniform random number scattergram；
Fig. 5 is the exponential random-number distribution figure of the present invention；
Fig. 6 is that 1000 group squares that above exponential random number span is evenly dividing of the present invention are gone forward side by side
The cartogram of row statistics；
Fig. 7 is the spectrum curve characteristic pattern extracting of the present invention；
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 this spectrum curve that finally gives of the present invention simulates design sketch；
Figure 11 shows the design sketch of the simulation random generating process of nuclear signal；
Figure 12 shows using DSMC direct sampling final effect figure.
Specific embodiment
In conjunction with the accompanying drawings, the present invention is further detailed explanation.These accompanying drawings are the schematic diagram of simplification,
The basic structure of the present invention is only described in a schematic way, therefore it only shows the composition relevant with the present invention.
As shown in figure 1, the simulation spectrum curve emulation mode of the present invention is by entering to actual nuclear spectrum curve chart
Spectrum curve simultaneously quantizes by row Curves Recognition, thus obtaining numerical value (the i.e. nuclear spectrum of spectrum curve each point
Energy level and the counting rate of each energy level), then obtained by this group numerical value of the random direct sampling of DSMC
With regard to the random number of each nuclear level, thus simulating the randomness of nuclear decay process, finally again to described random number
Carry out statistical disposition to obtain simulating spectrum curve, by inverting comparative simulation spectrum curve and actual spectrum curve
To determine reliability and the accuracy of imitative nuclear signal generator.
The specific embodiment 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, comprise the steps：
Step S1, obtains actual nuclear spectrum curve chart；
Step S2, is processed to actual nuclear spectrum curve chart to obtain simulation spectrum curve.
Optionally, described simulation spectrum curve emulation mode also includes：
Step S3, compares simulation spectrum curve and actual spectrum curve by inverting, can be set a song to music with obtaining simulation
Error between line and actual spectrum curve.
Further, as shown in figure 3, being processed to actual nuclear spectrum curve chart in described step S2 to obtain
The method of simulation spectrum curve includes：
Step S21, carries out Curves Recognition and spectrum curve quantizes, to obtain to actual nuclear spectrum curve chart
The numerical value of spectrum curve each point；Step S22, is obtained by this group numerical value of the random direct sampling of DSMC
To the random number with regard to each nuclear level, to simulate the randomness of nuclear decay process；And step S23, to described
Random number carries out statistical disposition and obtains described simulation spectrum curve.
Specifically, in described step S21, actual nuclear spectrum curve chart is carried out with Curves Recognition and by spectrum curve
The method of the numerical value to obtain spectrum curve each point of quantizing includes：
Actual nuclear energy is respectively set a song to music line chart the filtered denoising of image after actual spectrum curve figure is carried out
Display, and obtain and can set a song to music according to each key point that the actual spectrum curve figure of display copies this nuclear energy spectral curve
Line number evidence, to set up spectrum curve data base；Or by actual nuclear energy respectively set a song to music line chart image filtered, fall
Make an uproar pretreatment, Curves Recognition, curvilinear characteristic extract, and carry out interpolation processing with improve and repair disappearance energy
Each point data of spectral curve, to set up spectrum curve data base.
Wherein, closed by this group numerical value of the random direct sampling of DSMC in described step S22
Random number in each nuclear level is included with the method simulating the randomness of nuclear decay process：The nuclear signal time counts
Simulated behavior and nuclear signal amplitude statistics simulated behavior.
The method of described nuclear signal time statistical property simulation includes：Random number by obedience exponential is real
Existing nuclear signal time statistical property simulation, the wherein random number of exponential by (0,1] equally distributed random
Number is converted by inverse function method and obtains, and (0,1] uniform random number is suitable to try to achieve by linear congruential method.
Specifically, tried to achieve 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 choosing1Referred to as seed, has one to be fixed the generation quality of random number
Ring, its value is respectively 1～216Choose between=65535.For the ease of using on computers, generally 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 random number, that is, exponential random number can be realized by inverse function method, specifically
Process is as follows：
If the distribution function of stochastic variable X obeys exponential：
F (x)=1-e-ax, x >=0 (3)
Wherein, a is a time constant, and e is the nature truth of a matter.
Permissible by above formula, F (x) ∈ [0,1), and monotone decreasing in domain of definition, therefore function F (x) is in 0～+∞
Must there is inverse function, seek its inverse function：
Due to 0 ＜ 1-F (x)≤1, therefore above formula can be reduced to
Can be learnt by meeting by formula (5) (0,1] equally distributed Random sampling obtain obeying exponential with
Machine number x.
TakeThe exponential being produced by inverse function method by the equal distribution random numbers of above unit is random
Number scattergram such as Fig. 5.Above exponential random number span is evenly dividing with 1000 group squares and carries out
Statistics, final cartogram is as shown in Figure 6.
The method of described nuclear signal amplitude statistics simulated behavior includes：By to actual nuclear spectrum Curves Recognition simultaneously
Digitized obtains each energy level amplitude and counting rate, then by DSMC direct sampling and export described with
Machine number；Wherein identification digitized process to actual nuclear energy spectral curve include：
Step S221, is filtered to actual spectrum curve figure, noise reduction；Step S222, by between maximum kind
Split-run obtains threshold value, and spectrum curve figure is carried out binary conversion treatment, then is extracted by pixel scanning method
On spectrum curve, the numerical value of each point is coordinate；Step S223, is repaired to spectrum curve and is quantized.
Specifically, described by DSMC direct sampling and the method that exports described random number, that is, lead to
Cross each point value on DSMC direct sampling spectrum curve and curve, a series of random to obtain
Random number, thus to simulate the randomness of nuclear decay process.
The method in described step S221, actual spectrum curve figure being filtered, that is, to actual spectrum curve figure
Carry out Wiener filtering process, to filter the Gaussian noise of spectrum curve in figure, thus reducing noise as far as possible
The interference bringing.
The specific implementation process of the method for described nuclear signal amplitude statistics simulated behavior is as follows：
In described step S221, actual spectrum curve figure is filtered, the specific implementation step of noise reduction as follows：
By Wiener filtering, actual spectrum curve figure is filtered, noise reduction process, i.e. described Wiener filter
For a kind of linear filter, and also it is a kind of based on minimum mean square error criterion, the optimum to stationary process
Estimator.
Hypothesis Wiener filter input signal is s (t), superimposed noise n (t).Output signal x (t) passes through wave filter g (t)
Obtained by following convolution algorithm：
X (t)=g (t) * (s (t)+n (t)) (6)
For signal x (t) estimating it is desirable to be 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 wave filter output.
X (t) is write as convolution integral, that is,
Square error 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 x (t) and s (t) from phase
Close function.The final purpose of Wiener filtering is exactly to seek optimum g (t) so that E (e2) minimum.
In described step S222 by maximum kind between split-run obtain threshold value, and spectrum curve figure is carried out two-value
Change is processed, then is coordinate by the numerical value that pixel scanning method extracts each point on spectrum curve；
The specific algorithm process of maximum variance between clusters is as follows：
If the gray value of piece image is 1～m, wherein the pixel number for i for the gray value is ni, N represents figure
As pixel sum, then the probability that gray value occurs for i is：
Making gray value be more than threshold value k is c1Group, i.e. c1={ 1～k }, that gray value is more than threshold value k is then c2Group,
c2={ k+1～m }, then C1And C2The probability occurring is respectively：
It is calculated C1And C2Gray average be：
Wherein,So can obtain：
μr=ω1·μ1+ω2·μ2(16)
The variances sigma between two groups thus can be calculated2For：σ2(k)=ω1(μ1-μr)2+ω2(μ2-μr)2(17)
Formula (16) is substituted into formula (17) can obtain：σ2(k)=ω2ω2(μ2-μ1)2
So optimal threshold T*=Arg max { σ2(k) }, 0≤k ＜ m-1 (18)
Try to achieve segmentation threshold T*=0.6353.
That in described step S223, spectrum curve is repaired and quantized comprises the following steps that：
It is the number extracting each point on nuclear energy spectral curve after the filtered denoising of actual nuclear spectrum curve chart, binaryzation
Value is coordinate, need to extract spectrum curve feature, and curve is quantized.Detailed process is as follows：
First, Straight Line Identification, that is, pass through to scan the row and column of nuclear energy spectral curve binary map, identify nuclear spectrum
The straight line of in figure；
Secondly, fixed point, is judged horizontal stroke, the vertical coordinate of spectrum curve place coordinate system by the straight line identifying, and
Positioning initial point, typically from top to bottom, scans from left to right, and the straight line identifying is just horizontal, vertical seat
Mark；
3rd, spectrum curve feature extraction.For reduce image in frame and coordinate pair curve impact, need by
Frame filters.Again the point (two that pixel is 0 line by line or is scanned by column by pixel scan method after filtering frame
In value image, black is 0, and white is 1).
Finally, curve quantizes.After extracting curve, arrived by calculating the spectrum curve available point scanning
The row of scan origin and stringer distance determine that this pixel, in the position of in figure, is sat finally by being multiplied by expand
Target scale factor obtains the coordinate figure of this 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, can be seen that from Fig. 7 and Fig. 8, the spectrum curve figure of the simulation obtaining is than proper energy spectral curve
Lack in the data of some points.For truly reflecting actual spectrum curve characteristic as far as possible, need to scarce
The data lost carries out filling up reparation.
Specifically, described being filled up by cubic spline interpolation is lacked during spectrum curve feature extraction
The data point lost, and the numerical value obtaining each point on spectrum curve figure is stretched by the ratio expansion of coordinate, to realize
The data of disappearance is filled up effectively and repairs.
Described cubic spline interpolation is filling up during spectrum curve feature extraction the data point of disappearance
Specific algorithm is as follows：
Piecewise function S (x) on interval of definition [a, b], if meet：
1. S (x) is in each subinterval [xi, xi+1] on be a cubic polynomial function；
2. S (x) has continuous second dervative on entirely interval [a, b].
Then S (x) is called with regard to a=x on interval [a, b]0＜ x1＜ ... ＜ xnOne cubic spline function of=b.
Thus cubic spline interpolation problem is：N+1 node x of given function g (x)0, x1..., xnObtain function
y0, y1..., yn, seek cubic spline function S (x) so as to meet：
S(xj)=yj, j=0,1 ..., n (19)
Wherein, function S (x) is referred to as the cubic spline functions of g (x).
If S (x) is three sample spline interpolation functions of f (x), then must be fulfilled for following condition：
1. interpolation condition, that is,
S(xj)=yj, j=0,1 ..., n-1
2. the condition of continuity, that is,
3. the first derivative condition of continuity, that is,
4. the second dervative condition of continuity, that is,
Design sketch after cubic spline interpolation, as shown in figure 9, can be seen by its partial enlarged drawing,
Data point after cubic spline interpolation is more smooth, more approaches actual value.
Actual nuclear spectrum curve simulation effect, that is, to actual nuclear spectrum curve chart after above-mentioned image procossing,
Obtain this spectrum curve simulation design sketch eventually as shown in Figure 10.
Specifically, by each point value on DSMC direct sampling spectrum curve and curve, to obtain
Obtain a series of random randoms number, thus 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 actual spectrum curve using DSMC direct sampling final effect figure
Figure obtains energy level and this array of counting rate after quantizing after, random sampling procedure counting obtains.This figure
It is suitable to prove by the reasonability of Monte-Carlo step and accuracy).
Obtain simulating the numerical value of each point on spectrum curve and curve by numbers above image processing process
(abscissa is road location Channel, and vertical coordinate is counting rate Count), more directly taken out with DSMC
This group data of sample just can obtain a series of random energy levels, and (energy level such as but not limited to passes through multichannel analyzer
Quantization obtains, and described road location refers to what the energy of nuclear decay process release obtained after multichannel analyzer quantization)
Random number, thus to simulate the randomness of nuclear decay process.Finally more described random number is counted, can
To obtain simulating spectrum curve figure, so on the one hand can be with the reliability of checking system and accuracy, the opposing party
Face can also inverting in multichannel analyzer, in order to demarcate the accuracy of multichannel analyzer.
Using DSMC, probability P (A)=p (unknown) that certain event A is occurred, it is simulated meter
Calculate, circular：
(1) carry out n times and repeat independent sampling test, calculating event A frequency is nA.
Introduce stochastic variable Xi, represent event A frequency in i ＆ lt test, order
Then have
(2) event A that calculates repeats the occurrence frequency f in independent sampling test in n timesN, it is
(3) when N is fully big, with probability fNEstimated value as probability P (A)=pFor
(4) require estimated valueFor the unbiased esti-mator of probability P (A)=p, that is,
And direct sampling, that is, to nuclear signal, the characteristic on time and amplitude is with two groups of obedience different distributions
Random number simulating, and random number is discrete, discontinuous.For the sampling of discrete random sequence,
Direct sampling method is ideal.
The concrete sampling process of discrete distribution direct sampling method is as follows：
If the span of discrete random variable X is Xi(i=0,1,2,3 ...), its probability distribution is
P (X=Xi)=Pi(i=0,1,2,3 ...).Wherein Pi>=0,
(1) produce equally distributed random number r on (0,1) interval；
(2) try to achieve positive integer n=0,1,2... so that r meets
(3) sample value extracting discrete random variable X is X=Xn.And work as 0<r≤P0When X=X0；
(4) repeat step (1), (2), (3) are until extracting n sample value.
Due to producing (0,1) if equally distributed random number r is in intervalProbability be
I.e. eventThe probability occurring is equal to event X=XnThe probability occurring.
Again because random number r obeys being uniformly distributed on (0,1), its probability density function is
Its distribution function is as follows：
Therefore it is X=x that the random number r producing draws sample valuenProbability be
It follows that (X=X is drawn into by direct sampling methodn) probability be equivalent to random number XnIn random number sequence
Row X1, X2... XnThe frequency occurring.
Can be by following proof for direct sampling method reliability：
If X is discontinuous variable, its probability distribution is Pi=P { X=Xi, wherein i=1,2 ....X respectively with
PiObtain Xi, thenEvent | X-E (X) | >=ε represents that stochastic variable X obtains all satisfactions
Inequality | Xi- E (X) | the probable value X of >=εi, then
∵
∴
∴
Due to event X=XiThe probability that (i=0,1,2 ... N) occurs is pi(0<pi<1), then X ≠ XiGeneral
Rate is then 1-pi, and each X=XiThe probability occurring is constant, and each sampling resultses are with other each extractions
Result is unrelated.Therefore X=XiIndividual event is a bernoulli experiment, then sampling n time, then weigh primary for n
Nu Li tests.If making event A (X=Xi) number of times that occurs is nA, i.e. nA～B (n, p).Due to X1, X2...,
XnIt is n stochastic variable that is separate and obeying the 0-1 distribution that parameter is p, and
HaveD (Xi)=p (1-p), i=1,2 ... n.Giving any ε ＞ 0 then has
Be can be derived from by (4.31) formula
And
Therefore can be derived from
Abbreviation obtains
I.e. when the frequency n extracting is bigger, after sampling, the number of times of event A appearance and the frequency ratio of population of samples are got over
The probability occurring close to event A.
Extracting its error of random number by direct sampling is：
OrderTherefore
I.e.It is the unbiased esti-mator of p,
I.e. the number of times of sampling n is bigger, estimated valueCloser to theoretical value p.
With the above-mentioned desirable embodiment according to the present invention for enlightenment, by above-mentioned description, related work
Personnel can carry out various change and modification completely in the range of without departing from this invention technological thought.
The technical scope of this invention is not limited to content in description it is necessary to according to right
To determine its technical scope.
Claims (10)
1. a kind of simulation spectrum curve emulation mode is it is characterised in that comprise the steps：
Step S1, obtains actual nuclear spectrum curve chart；
Step S2, is processed to actual nuclear spectrum curve chart to obtain simulation spectrum curve.
2. simulation spectrum curve emulation mode according to claim 1 is it is characterised in that described simulation
Spectrum curve emulation mode also includes：
Step S3, compares simulation spectrum curve and actual spectrum curve by inverting, can be set a song to music with obtaining simulation
Error between line and actual spectrum curve.
3. according to claim 1 and 2 simulation spectrum curve emulation mode it is characterised in that
In step S2, actual nuclear spectrum curve chart is processed and included with the method obtaining simulation spectrum curve：
Step S21, carries out Curves Recognition and spectrum curve quantizes, to obtain to actual nuclear spectrum curve chart
The numerical value of spectrum curve each point；
Step S22, is obtained with regard to each nuclear level by this group numerical value of the random direct sampling of DSMC
Random number, to simulate the randomness of nuclear decay process；
Step S23, carries out statistical disposition and obtains described simulation spectrum curve to described random number.
4. according to claim 3 simulation spectrum curve emulation mode it is characterised in that
In described step S21, actual nuclear spectrum curve chart is carried out Curves Recognition and by spectrum curve quantize with
The method obtaining the numerical value of spectrum curve each point includes：
Actual nuclear energy is respectively set a song to music line chart the filtered denoising of image after actual spectrum curve figure is carried out
Display, and obtain and can set a song to music according to each key point that the actual spectrum curve figure of display copies this nuclear energy spectral curve
Line number evidence, to set up spectrum curve data base；Or
By actual nuclear energy respectively set a song to music line chart image is filtered, noise reduction pretreatment, Curves Recognition, curvilinear characteristic
Extract, and carry out interpolation processing to improve and to repair each point data of disappearance spectrum curve, to set up power spectrum
Diagram database.
5. according to claim 4 simulation spectrum curve emulation mode it is characterised in that
Obtained with regard to each nuclear level by this group numerical value of the random direct sampling of DSMC in step S22
Random number included with the method simulating the randomness of nuclear decay process：
Nuclear signal time statistical property is simulated；And nuclear signal amplitude statistics simulated behavior.
6. according to claim 5 simulation spectrum curve emulation mode it is characterised in that
The method of described nuclear signal time statistical property simulation includes：Random number by obedience exponential is real
Existing nuclear signal time statistical property simulation, wherein
The random number of exponential by (0,1] equally distributed random number by inverse function method conversion obtain, and
(0,1] uniform random number is suitable to try to achieve by linear congruential method.
7. simulation spectrum curve emulation mode according to claim 6 is it is characterised in that described core is believed
The method of number amplitude statistics simulated behavior includes：
By to actual nuclear spectrum Curves Recognition and digitized obtains each energy level amplitude and counting rate, then by covering
Special Caro method direct sampling simultaneously exports described random number；Wherein
Identification digitized process to actual nuclear energy spectral curve include：
Step S221, is filtered to actual spectrum curve figure, noise reduction；
Step S222, by maximum kind between split-run obtain threshold value, and spectrum curve figure is carried out at binaryzation
Reason, then be coordinate by the numerical value that pixel scanning method extracts each point on spectrum curve；
Step S223, is repaired to spectrum curve and is quantized.
8. simulation spectrum curve emulation mode according to claim 7 is it is characterised in that described pass through
DSMC direct sampling the method exporting described random number, that is,
By each point value on DSMC direct sampling spectrum curve and curve, a series of to obtain
Random number, thus to simulate the randomness of nuclear decay process.
9. simulation spectrum curve emulation mode according to claim 8 is it is characterised in that described step
The method in S221, actual spectrum curve figure being filtered, carries out Wiener filtering to actual spectrum curve figure
Process, to filter the Gaussian noise of spectrum curve in figure.
10. according to claim 9 simulation spectrum curve emulation mode it is characterised in that
In described step S223, spectrum curve is repaired and is quantized including：By cubic spline interpolation
To fill up during spectrum curve feature extraction the data point of disappearance, and expanded by the ratio of coordinate stretch with
Obtain the numerical value of each point on spectrum curve figure.
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| 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 |
| CN201510539124.1A CN106483550B (en) | 2015-08-28 | 2015-08-28 | A kind of simulation spectrum curve emulation mode |
| 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 |
|---|---|---|---|
| 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 |
| CN201811086619.3A Division CN109239763B (en) | 2015-08-28 | 2015-08-28 | Simulation method for simulating nuclear decay process by using simulated energy spectrum curve |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106483550A true CN106483550A (en) | 2017-03-08 |
| CN106483550B CN106483550B (en) | 2018-09-18 |
Family
ID=58234736
Family Applications (4)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811085398.8A Active CN109271707B (en) | 2015-08-28 | 2015-08-28 | Simulation energy spectrum curve simulation method for simulating nuclear energy spectrum line |
| CN201811086619.3A Active CN109239763B (en) | 2015-08-28 | 2015-08-28 | Simulation method for simulating nuclear decay process by using simulated energy spectrum curve |
| CN201811085406.9A Active CN109086555B (en) | 2015-08-28 | 2015-08-28 | Simulation method for simulating energy spectrum curve by adopting Monte Carlo method |
| 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 |
|---|---|---|---|
| CN201811085398.8A Active CN109271707B (en) | 2015-08-28 | 2015-08-28 | Simulation energy spectrum curve simulation method for simulating nuclear energy spectrum line |
| CN201811086619.3A Active CN109239763B (en) | 2015-08-28 | 2015-08-28 | Simulation method for simulating nuclear decay process by using simulated energy spectrum curve |
| CN201811085406.9A Active CN109086555B (en) | 2015-08-28 | 2015-08-28 | Simulation method for simulating energy spectrum curve by adopting Monte Carlo method |
Country Status (1)
| Country | Link |
|---|---|
| CN (4) | CN109271707B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108416355A (en) * | 2018-03-09 | 2018-08-17 | 浙江大学 | A kind of acquisition method of the industry spot creation data based on machine vision |
| CN111276195A (en) * | 2020-01-20 | 2020-06-12 | 深圳大榆树科技有限公司 | Method for calculating number of compounds in gel energy spectrum by using maximum inter-class variance method |
| CN112462676A (en) * | 2021-01-27 | 2021-03-09 | 泛华检测技术有限公司 | Device capable of simulating overlapped nuclear pulse signal generation and control method thereof |
| CN112462675A (en) * | 2021-01-27 | 2021-03-09 | 泛华检测技术有限公司 | Cascaded dual-index nuclear pulse signal generating device and control method thereof |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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)
| 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)
| 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 |
| CN102177444B (en) * | 2008-10-10 | 2014-07-09 | 皇家飞利浦电子股份有限公司 | 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 |
-
2015
- 2015-08-28 CN CN201811085398.8A patent/CN109271707B/en active Active
- 2015-08-28 CN CN201811086619.3A patent/CN109239763B/en active Active
- 2015-08-28 CN CN201811085406.9A patent/CN109086555B/en active Active
- 2015-08-28 CN CN201510539124.1A patent/CN106483550B/en active Active
Patent Citations (6)
| 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)
| Title |
|---|
| 万东阳 等: "《数字核信号图形化识别和处理研究》", 《第十六届全国核电子学与核探测技术学术年会》 * |
| 谭承君 等: "《基于随机抽样的核脉冲信号发生器的研究》", 《原子能科学技术》 * |
| 霍建文 等: "《任意分布的高速仿核信号发生器》", 《核电子学与探测技术》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108416355A (en) * | 2018-03-09 | 2018-08-17 | 浙江大学 | A kind of acquisition method of the industry spot creation data based on machine vision |
| CN111276195A (en) * | 2020-01-20 | 2020-06-12 | 深圳大榆树科技有限公司 | Method for calculating number of compounds in gel energy spectrum by using maximum inter-class variance method |
| CN111276195B (en) * | 2020-01-20 | 2023-09-19 | 深圳大榆树科技有限公司 | Method for calculating number of gel energy spectrogram compounds by using maximum inter-class variance method |
| CN112462676A (en) * | 2021-01-27 | 2021-03-09 | 泛华检测技术有限公司 | Device capable of simulating overlapped nuclear pulse signal generation and control method thereof |
| CN112462675A (en) * | 2021-01-27 | 2021-03-09 | 泛华检测技术有限公司 | Cascaded dual-index nuclear pulse signal generating device and control method thereof |
| CN112462675B (en) * | 2021-01-27 | 2021-05-07 | 泛华检测技术有限公司 | Cascaded dual-index nuclear pulse signal generating device and control method thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109239763A (en) | 2019-01-18 |
| CN109086555B (en) | 2022-04-19 |
| CN109271707B (en) | 2022-04-19 |
| CN106483550B (en) | 2018-09-18 |
| CN109239763B (en) | 2022-04-19 |
| CN109086555A (en) | 2018-12-25 |
| CN109271707A (en) | 2019-01-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106483550A (en) | A kind of simulation spectrum curve emulation mode | |
| CN106483842A (en) | A kind of high-precision intelligent is imitated nuclear signal and system and its method of work is occurred | |
| CN106483551B (en) | A kind of imitative nuclear signal generator and its working method | |
| CN104134204B (en) | Image definition evaluation method and image definition evaluation device based on sparse representation | |
| Venema et al. | A stochastic iterative amplitude adjusted Fourier transform algorithm with improved accuracy | |
| CN109784410B (en) | Characteristic extraction and classification method for ship radiation noise signals | |
| CN110059845B (en) | Metering device clock error trend prediction method based on time sequence evolution gene model | |
| CN108181613B (en) | Non-uniform quantitative sorting method for sequence difference values of PRI (pulse repetition index) jittering signals | |
| CN108681689B (en) | Frame rate enhanced gait recognition method and device based on generation of confrontation network | |
| CN108776582A (en) | A kind of true random number sequence production method based on quantum tunneling effect | |
| Szarka III et al. | On the equivalence of the Bernoulli and geometric CUSUM charts | |
| CN112989361A (en) | Model security detection method based on generation countermeasure network | |
| CN109597123B (en) | Effective signal detection method and system | |
| CN116432703B (en) | Pulse height estimation method, system and terminal based on composite neural network model | |
| CN106483867A (en) | There is system processor, system and method for work in imitative nuclear signal | |
| CN112183260A (en) | One-way valve fault diagnosis method based on total variation noise reduction and RQA | |
| CN114553606B (en) | Industrial control network intrusion detection method and system | |
| CN107886113B (en) | Electromagnetic spectrum noise extraction and filtering method based on chi-square test | |
| CN114972330A (en) | Workpiece surface roughness detection optimization method based on improved histogram homogenization algorithm | |
| CN111008356A (en) | WTSVD algorithm-based background-subtracted gamma energy spectrum set analysis method | |
| Rao et al. | An Automatic CADI’s Ionogram Scaling Software Tool for Large Ionograms Data Analytics | |
| CN111931412A (en) | Underwater target noise LOFAR spectrogram simulation method based on generative countermeasure network | |
| CN112508862A (en) | Method for enhancing magneto-optical image of crack by improving GAN | |
| CN110988802A (en) | Radar radiation source identification system based on signal scale decomposition | |
| CN114638851B (en) | Image segmentation method, system and storage medium based on generation countermeasure network |
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 |