CN103076622A - Generation method of random signals for spectrum stabilization - Google Patents

Generation method of random signals for spectrum stabilization Download PDF

Info

Publication number
CN103076622A
CN103076622A CN2012104265996A CN201210426599A CN103076622A CN 103076622 A CN103076622 A CN 103076622A CN 2012104265996 A CN2012104265996 A CN 2012104265996A CN 201210426599 A CN201210426599 A CN 201210426599A CN 103076622 A CN103076622 A CN 103076622A
Authority
CN
China
Prior art keywords
random
function
spectrum
velocity
space
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012104265996A
Other languages
Chinese (zh)
Other versions
CN103076622B (en
Inventor
黄洪全
丁卫撑
王超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN201210426599.6A priority Critical patent/CN103076622B/en
Publication of CN103076622A publication Critical patent/CN103076622A/en
Application granted granted Critical
Publication of CN103076622B publication Critical patent/CN103076622B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a generation method of random signals for spectrum stabilization, which comprises the following steps that firstly, one energy spectrum which is measured in radiometry is selected as a reference energy spectrum, and in the condition that a two-norm error is set, the variance of a generating function and the scale space of a basis function are adjusted to map function space of the energy spectrum; secondly, stable change and unstable random fluctuation of random velocity simulation test conditions are adopted; then a random velocity parameter adjustment method is adopted to adjust parameters of the function space; and finally, random sampling is performed in accordance with a time slice segmentation method, so as to generate random numbers. The generation method can provide various rich and realistic random signals for spectrum stabilization research and is an effective method.

Description

A kind of spectrum stabilization production method of random signal
Technical field
A kind of spectrum stabilization production method of random signal.
Background technology
In carrying out radioactivity survey, because detector and interlock circuit (such as pulse amplifier, pulse height analyzer, high-voltage power supply etc.) thereof are subjected to the impact of environment temperature, and the instability of instrument self (such as the fatigue effect of components and parts, catabiosis etc.), cause energy spectrometer to be compared with " standard spectrum " (to the same object) under the established condition in the spectrum of the instrument under certain condition, its peak position or spectrum shape change, and namely spectrum is floated.In energy spectrometer, generally all to take spectrum-stabilizing device or the method for certain form, adopt " reference source " spectrum stabilization method and adopt the computer software spectrum stabilization technology etc. of "dead" reference source such as energy spectrometer.For the validity of authentication processing method or instrument spectrum stabilization, usually to float the spectrum of complexity and analyze, if produce such radioactivity random signal by simulation, will greatly reduce workload and overcome environmental limit.The present invention produces the randomness measuring condition by effective model, and simulation produces and form the random signal that complicated spectrum is floated, for the research of spectrum-stabilizing device or method provides effective random signal.
Summary of the invention
The object of the invention is to disclose the production method that a kind of spectrum stabilization is used random signal.This invention can provide for the research of spectrum-stabilizing device or method effective random signal.
The present invention is achieved by the following technical solutions, and 1. concrete steps following～4..
1. a certain power spectrum that has recorded in the selective emission measurement is set up the Maps of Function Spaces of this power spectrum as the reference power spectrum: under two norm errors of setting, the adjustment of the metric space by generating function variance and basis function realizes the Maps of Function Spaces of power spectrum.
2. adopt steady change and the astable random fluctuation of random velocity simulation test condition:
At first set the smooth change speed of an initial velocity simulation test condition, then adopt the astable random fluctuation of the acceleration simulation test condition that produces at random; Random velocity is that initial velocity adds that the accumulated time of acceleration obtains; Acceleration adopts arma modeling to obtain: be zero white Gaussian noise input arma modeling system with average, the at random output valve of system is acceleration.
3. adopt random velocity parameter adjustment method that the parameter of function space is adjusted: with required random velocity regulating the speed as function space basis function average and variance.
4. produce random number by the random sampling of timeslice split plot design:
At first will be divided into the continuous time period test duration; Next asks for the variable quantity of per time period inner function spatial parameter; Then the represented power spectrum of function space new in per time period is carried out normalization as probability density function, and produce random number by sampling.
The invention has the beneficial effects as follows: under two norm errors of setting, a certain power spectrum of selecting is finished Maps of Function Spaces by the adjustment of the metric space of generating function variance and basis function, can realize by the precision of expectation the parametrization of this power spectrum; Adopting steady change and the astable random fluctuation of random velocity simulation test condition, is the randomness characteristics that combine test condition in the radioactivity survey; Adopting random velocity parameter adjustment method that the parameter of function space is adjusted, is that randomness and the parameterized power spectrum with test condition contacts; Produce random number by the random sampling of timeslice split plot design at last, realized the combination of time and space randomness.Therefore, the spectrum that the model that the present invention adopts and method can Simulation of Complex is floated process, for the research of spectrum-stabilizing device or method provides random signal rich and varied and true to nature.
Description of drawings
Fig. 1 is the process flow diagram of the inventive method.
Embodiment
For the purpose, technical scheme and the advantage that make invention is clearer, referring to the accompanying drawing embodiment that develops simultaneously, the present invention is described in further detail.
For the research of giving spectrum-stabilizing device or method provides random signal rich and varied and true to nature, the invention provides the production method that a kind of spectrum stabilization is used random signal.Fig. 1 has shown the flow process of the method for the invention.1. concrete steps following～4. go on foot described.
1. a certain power spectrum that has recorded in the selective emission measurement is as the reference power spectrum, set up the Maps of Function Spaces of this power spectrum: under two norm errors of setting, the adjustment of the metric space by generating function variance and basis function realizes the Maps of Function Spaces of power spectrum, comprises following steps A～D.
The power spectrum of for simplicity, establishing selection is y( n) ( n=1 N), wherein NBe Zong Dao location number, tale is N t
A sets two norm errors on demand ε 0, select initial generating function variance σ 2And the structure generating function, set the initial gauges space V j, wherein j is integer, makes up V jThe basis function of metric space:
Two norm errors ε 0Decide variance by required precision σ 2Decide metric space by the smooth degree of power spectrum V jUsually be chosen as j greater than 1;
If the initial generating function of selecting is Gauss type function, this generating function can be expressed as
Figure 791227DEST_PATH_IMAGE001
(1)；
V jThe basis function of metric space is by generating function
Figure 871310DEST_PATH_IMAGE002
The integer translation obtain
Figure 866948DEST_PATH_IMAGE003
(2)；
It is basis function
Figure 847411DEST_PATH_IMAGE004
For
Figure 168671DEST_PATH_IMAGE005
(3)；
In formula (1) and the formula (3)
Figure 470470DEST_PATH_IMAGE006
Need to prove that generating function may be selected to be other function of non-Gaussian.
B is in the initial gauges space V jInteriorly set up selected power spectrum y( n) Maps of Function Spaces:
At first, with basis function
Figure 473061DEST_PATH_IMAGE007
Carry out discretize
Figure 384996DEST_PATH_IMAGE008
nBe integer (4),
Try to achieve y( n) V j On Maps of Function Spaces be
Figure 611578DEST_PATH_IMAGE009
(5)，
In the formula (5) c Jk Provided by following formula
Figure 666253DEST_PATH_IMAGE010
(6)；
Then, trying to achieve tale by following formula is N tMaps of Function Spaces
Figure 269273DEST_PATH_IMAGE011
For:
Figure 958749DEST_PATH_IMAGE012
(7)。
C calculates two norm errors ε:
Figure 559495DEST_PATH_IMAGE013
(8)。
D works as The time, illustrating that mapping error meets the demands, 1. this step finishes; Otherwise, continue to process by following (a), (b) step:
(a) reduce the generating function variance σ 2, but
Figure 304914DEST_PATH_IMAGE015
, wherein
Figure 117405DEST_PATH_IMAGE016
Be the lower limit variance, return above steps A ~ D and proceed;
(b) the initial gauges space is set as V J-1, and reset initial generating function variance σ 2, return above steps A ~ D and proceed.
2. adopt steady change and the astable random fluctuation of random velocity simulation test condition, comprise following steps A～C.
A at first sets an initial velocity v 0 , in order to the smooth change speed of simulation test condition.
Then B adopts the astable random fluctuation of the acceleration simulation test condition that produces at random:
Acceleration adopts arma modeling to obtain, and being about to average is zero white Gaussian noise input arma modeling system, and the at random output valve of system is acceleration; The step that produces acceleration is undertaken by following (a)～(b):
(a) establishing acceleration is x( n), usefulness ARMA ( p, q) this stochastic process of model approximate description gets
Figure 685790DEST_PATH_IMAGE017
(9)
In the formula (9) e( n) be that average is zero white Gaussian noise, all the other parameters set up on their own as required;
(b) to the coefficient of formula (9) initial setting With
Figure 597562DEST_PATH_IMAGE019
Revise, with the ARMA that obtains having causality and minimum phase ( p, q) the model process:
Formula (9) is become following form
Figure 323948DEST_PATH_IMAGE020
(10)
If ARMA ( p, q) there is following distribution in the pole and zero of model
Figure 63233DEST_PATH_IMAGE021
(11)
Figure 314217DEST_PATH_IMAGE022
(12)
Figure 62731DEST_PATH_IMAGE023
(13)
(14)
Then the pole and zero inverting that unit circle is outer is in unit circle, obtain revised ARMA ( p, q) the model process is
Figure 760002DEST_PATH_IMAGE025
(15)
In the formula (15)
Figure 498282DEST_PATH_IMAGE026
,
Figure 784906DEST_PATH_IMAGE027
Use acceleration x( n) the astable random fluctuation of random value simulation test condition.
The C random velocity adds that by initial velocity the accumulated time of acceleration obtains:
If t random velocity constantly is v( t), then random velocity provides by following formula
Figure 423567DEST_PATH_IMAGE028
(16)
T is time constant in the formula, for simplicity, and order x( n)=0, n＜0.
3. adopt random velocity parameter adjustment method that the parameter of function space is adjusted: required random velocity as the regulating the speed of function space basis function average and variance, is provided by formula (17) and formula (18)
Figure 770235DEST_PATH_IMAGE029
(17)
(18)
In the formula
Figure 135062DEST_PATH_IMAGE031
With
Figure 192010DEST_PATH_IMAGE032
Representative function space basis function average is regulated the speed and is regulated the speed with variance respectively, With Be respectively the velocity factor of average and variance.
4. produce random number by the random sampling of timeslice split plot design, comprise following steps A～C.
A at first will be divided into the continuous time period test duration:
If the test duration is t, it is divided into continuous time period sequence by time constant T, i.e. 0 ~ T, T ~ 2T, 2T ~ 3T ...,
Figure 362463DEST_PATH_IMAGE035
,
Figure 398553DEST_PATH_IMAGE036
For simplicity, they are designated as respectively
Figure 87023DEST_PATH_IMAGE037
, ,
Figure 943476DEST_PATH_IMAGE039
...
Secondly B asks for the variable quantity of per time period inner function spatial parameter, comprises following steps (a) and (b):
(a) in each time period
Figure 912701DEST_PATH_IMAGE040
,
Figure 772072DEST_PATH_IMAGE038
,
Figure 692493DEST_PATH_IMAGE039
...,
Figure 600406DEST_PATH_IMAGE041
In ask for the function space mean parameter uVariable quantity be
Figure 673404DEST_PATH_IMAGE042
(b) in each time period
Figure 454409DEST_PATH_IMAGE037
,
Figure 612858DEST_PATH_IMAGE038
,
Figure 638976DEST_PATH_IMAGE039
...,
Figure 832060DEST_PATH_IMAGE041
In ask for the function space parameter variance σVariable quantity be
Figure 783967DEST_PATH_IMAGE043
Adopt the timeslice method to make things convenient for calculating, also guaranteed the real-time update of parameter simultaneously.
Then C carries out normalization as probability density function to the represented power spectrum of function space new in per time period, and produces random number by sampling, comprises following steps (a) and (b):
(a) at first ask for interior new power spectrum of per time period
Figure 101816DEST_PATH_IMAGE044
In time period, function space V j Basis function by following correction
Figure 679428DEST_PATH_IMAGE045
(21),
Then can try to achieve power spectrum y( n) at new function space V j On be mapped as
(22)，
In the formula (22) c Jk With 1. behind the EOS c Jk Value is equal,
Figure 348361DEST_PATH_IMAGE047
Be new power spectrum;
(b) within per time period to the new power spectrum generation random number of sampling
Will New power spectrum in time period
Figure 82279DEST_PATH_IMAGE047
Normalization is as probability density function, and the discrete distribution functional form that is expressed as:
Figure 922059DEST_PATH_IMAGE048
（23）
Wherein N n For
Figure 717233DEST_PATH_IMAGE047
The nThe counting of location, road, For Tale, and provided by following formula
（24）；
In the formula (24) NFor
Figure 727466DEST_PATH_IMAGE047
Zong Dao location number;
Produce random number by formula (25) sampling again X F , try to achieve obedience F( x) random number that distributes X F , this random number is obeyed exactly
Figure 507203DEST_PATH_IMAGE047
Spectral distribution
Figure 964730DEST_PATH_IMAGE044
Random number in time period, wherein εFor obeying [0-1] equally distributed random number;
Figure 592151DEST_PATH_IMAGE051
（25）
Except adopting above method sampling, also can add sampling by Gaussian function and produce random number, here no longer narration.
By 1. above～4. the step,
Figure 476930DEST_PATH_IMAGE037
,
Figure 58478DEST_PATH_IMAGE038
,
Figure 319695DEST_PATH_IMAGE039
... random number Deng generation in each time period is exactly the random signal of studying for spectrum stabilization.
Can find out with the production method of random signal from above-mentioned spectrum stabilization, the present invention combines the variation of radioactivity survey condition to the impact of power spectrum, the power spectrum of selecting is finished Maps of Function Spaces by the adjustment of the metric space of generating function variance and basis function, realize the parametrization of this power spectrum by the precision of expectation; Combine the randomness characteristics of test condition in the radioactivity survey, adopt steady change and the astable random fluctuation of random velocity simulation test condition; Randomness and the parameterized power spectrum of test condition are contacted, adopt random velocity parameter adjustment method that the parameter of function space is adjusted; Last binding time and space randomness produce the spectrum stabilization random signal by the random sampling of timeslice split plot design.Therefore, the spectrum that the model that the present invention adopts and method can Simulation of Complex is floated the random signal in the process, for the research of spectrum-stabilizing device or method provides random signal rich and varied and true to nature.
In the embodiment of the invention described above; the production method of spectrum stabilization in the radioactivity survey with random signal had been described in detail; but it should be noted that; the above only is one embodiment of the present of invention; the present invention can select the generating function except Gaussian function equally, and the spectrum stabilization research that also can be local spectral coverage provides random signal, and is within the spirit and principles in the present invention all; any modification of doing, be equal to replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (5)

1. a spectrum stabilization is characterized in that concrete steps are as follows with the production method of random signal:
1. a certain power spectrum that has recorded in the selective emission measurement is set up the Maps of Function Spaces of this power spectrum as the reference power spectrum;
2. adopt steady change and the astable random fluctuation of random velocity simulation test condition;
3. adopt random velocity parameter adjustment method that the parameter of function space is adjusted;
4. produce random number by the random sampling of timeslice split plot design.
2. spectrum stabilization according to claim 1 is characterized in that with the production method of random signal, described 1. in Maps of Function Spaces, refer to setting under two norm errors, the adjustment of the metric space by generating function variance and basis function realizes.
3. spectrum stabilization according to claim 2 is with the production method of random signal, it is characterized in that, described 2. middle steady change and the astable random fluctuation of adopting random velocity simulation test condition, refer at first set the smooth change speed of an initial velocity simulation test condition, then adopt the astable random fluctuation of the acceleration simulation test condition that produces at random; Random velocity is that initial velocity adds that the accumulated time of acceleration obtains; Acceleration adopts arma modeling to obtain: be zero white Gaussian noise input arma modeling system with average, the at random output valve of system is acceleration.
4. spectrum stabilization according to claim 3 is with the production method of random signal, it is characterized in that, described 3. middle employing random velocity parameter adjustment method is adjusted the parameter of function space, refers to required random velocity regulating the speed as function space basis function average and variance.
5. spectrum stabilization according to claim 4 is characterized in that with the production method of random signal, describedly produces random number by the random sampling of timeslice split plot design in 4., refers to: at first will be divided into the continuous time period test duration; Next asks for the variable quantity of per time period inner function spatial parameter; Then the represented power spectrum of function space new in per time period is carried out after the normalization as probability density function, and sampling produces random number.
CN201210426599.6A 2012-10-31 2012-10-31 A kind of production method of spectrum stabilization stochastic signal Expired - Fee Related CN103076622B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210426599.6A CN103076622B (en) 2012-10-31 2012-10-31 A kind of production method of spectrum stabilization stochastic signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210426599.6A CN103076622B (en) 2012-10-31 2012-10-31 A kind of production method of spectrum stabilization stochastic signal

Publications (2)

Publication Number Publication Date
CN103076622A true CN103076622A (en) 2013-05-01
CN103076622B CN103076622B (en) 2016-08-17

Family

ID=48153209

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210426599.6A Expired - Fee Related CN103076622B (en) 2012-10-31 2012-10-31 A kind of production method of spectrum stabilization stochastic signal

Country Status (1)

Country Link
CN (1) CN103076622B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105717140A (en) * 2016-01-29 2016-06-29 成都理工大学 Chang'e 2 gamma spectrum and nonlinear spectrum drift calibration method based on energy correspondence
CN105785424A (en) * 2016-02-25 2016-07-20 中国人民解放军63973部队 Cadmium zinc telluride detector gamma spectrum full-energy peak non-linear fitting algorithm
CN109086555A (en) * 2015-08-28 2018-12-25 易良碧 Using the simulation spectrum curve emulation mode of Monte Carlo method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1045186A (en) * 1989-02-22 1990-09-05 成都地质学院 Digital spectrum-stabilizing device with natural low-energy peak source for referencial use
EP0416970A2 (en) * 1989-08-30 1991-03-13 Schlumberger Limited Nuclear spectroscopy signal stabilization and calibration method and apparatus
US20090166542A1 (en) * 2007-12-29 2009-07-02 Ming Ruan Radiation detection system using a multichannel spectrometer and method thereof
CN102073060A (en) * 2009-11-24 2011-05-25 成都理工大学 Simulation method for random properties of nuclear signals
CN102298153A (en) * 2010-06-23 2011-12-28 成都理工大学 Method for decomposing multiple spectral peaks during radioactive measurement

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1045186A (en) * 1989-02-22 1990-09-05 成都地质学院 Digital spectrum-stabilizing device with natural low-energy peak source for referencial use
EP0416970A2 (en) * 1989-08-30 1991-03-13 Schlumberger Limited Nuclear spectroscopy signal stabilization and calibration method and apparatus
US20090166542A1 (en) * 2007-12-29 2009-07-02 Ming Ruan Radiation detection system using a multichannel spectrometer and method thereof
CN102073060A (en) * 2009-11-24 2011-05-25 成都理工大学 Simulation method for random properties of nuclear signals
CN102298153A (en) * 2010-06-23 2011-12-28 成都理工大学 Method for decomposing multiple spectral peaks during radioactive measurement

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张海勇等: "非平稳随机信号的参数模型分析方法", 《系统工程与电子技术》 *
罗智中: "随机信号AR频谱模型的本质研究", 《振动与冲击》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109086555A (en) * 2015-08-28 2018-12-25 易良碧 Using the simulation spectrum curve emulation mode of Monte Carlo method
CN105717140A (en) * 2016-01-29 2016-06-29 成都理工大学 Chang'e 2 gamma spectrum and nonlinear spectrum drift calibration method based on energy correspondence
CN105785424A (en) * 2016-02-25 2016-07-20 中国人民解放军63973部队 Cadmium zinc telluride detector gamma spectrum full-energy peak non-linear fitting algorithm
CN105785424B (en) * 2016-02-25 2019-02-12 中国人民解放军63973部队 A kind of tellurium-zinc-cadmium detector gamma spectrum full energy peak Nonlinear Quasi hop algorithm

Also Published As

Publication number Publication date
CN103076622B (en) 2016-08-17

Similar Documents

Publication Publication Date Title
Zhou et al. Estimation of the dynamic states of synchronous machines using an extended particle filter
CN108563867A (en) A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM
CN105352487A (en) Calibration method for precision of attitude measurement system
WO2019181032A1 (en) Wind flow sensing system and method for determining velocity fields of wind flow
CN109448124B (en) Water quality simulation method and device for river channel
CN106842306A (en) The staggered-mesh finite difference analogy method and device of a kind of global optimization
KR20190017454A (en) Device and method for generating location estimation model
Zhang et al. Measurement data fitting based on moving least squares method
CN109633570A (en) A kind of random Chaff Clouds modeling method of quick arbitrary shape
Shin et al. The local ensemble transform Kalman filter (LETKF) with a global NWP model on the cubed sphere
Nan et al. Groundwater parameter estimation using the ensemble Kalman filter with localization
EP4127456A1 (en) System and method for fast wind flow measurement by lidar in a complex terrain
CN103076622A (en) Generation method of random signals for spectrum stabilization
CN117252863A (en) Quick detection and analysis method for geographic information abnormal data
CN108966343A (en) A kind of self calibration localization method based on ultra wide band Location-Unknown anchor node
CN104105049A (en) Room impulse response function measuring method allowing using quantity of microphones to be reduced
Li et al. Uncertainty modelling and analysis of volume calculations based on a regular grid digital elevation model (DEM)
EP4127457A1 (en) System and method for wind flow turbulence measurement by lidar in a complex terrain
CN115979245B (en) Magnetic induction network positioning method based on self-calibration estimation
CN106815607A (en) A kind of isogram based on anti-distance weighting interpolation inverse function is as data extraction method
CN106528956A (en) Ray tracing model-based method for predicting field intensity by data interpolation method
KR20190042540A (en) Device and method for generating geomagnetic sensor based location estimation model using artificial neural networks
CN114019445B (en) Two-dimensional arrival angle measurement method based on position clustering dynamic sparse reconstruction
Minár et al. Towards exactness in geomorphometry
Duan et al. SimStr: A C++ program for offshore seismic towed-streamer positioning data simulation

Legal Events

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

Granted publication date: 20160817

Termination date: 20171031