CN101561293A - Weak characteristic information recovery system based on stochastic resonance - Google Patents
Weak characteristic information recovery system based on stochastic resonance Download PDFInfo
- Publication number
- CN101561293A CN101561293A CNA2009100687840A CN200910068784A CN101561293A CN 101561293 A CN101561293 A CN 101561293A CN A2009100687840 A CNA2009100687840 A CN A2009100687840A CN 200910068784 A CN200910068784 A CN 200910068784A CN 101561293 A CN101561293 A CN 101561293A
- Authority
- CN
- China
- Prior art keywords
- mimic channel
- signal
- bistable
- multiplier
- recovery
- 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
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The invention discloses a weak characteristic information recovery system based on stochastic resonance. The system consists of a sensor, a bistable analog circuit, a recovery analog circuit and a filter module, wherein practical sensed signals of the sensor as primary data are sent to the bistable analog circuit and the like in the system. The recovery analog circuit consists of three constructed analog multiplier chips, a direct current power supply and a feedback amplifier, wherein an output end of a third multiplier is connected to the feedback amplifier to carry out feedback output for the signals. The system adopts the principle that: the signals practically sensed by the sensor are full submerged by noise and cannot be distinguished, and the signal to noise ratio of the primary weak signals input by the system is greatly enhanced through the processing such as the stochastic resonance, waveform recovery, filtration and the like. The system solves the problem of fine detection of weak signals in strong noise, breaks through the limitation that the prior low-pass filter distorts the information detection for the weak signals in the strong noise, and can realize the functions of multi-period, multi-frequency and multi-component complex weak characteristic information identification and signal reduction, and the like.
Description
Technical field
The invention belongs to detection technique and signal Processing, be specifically related to a kind of method of utilizing the accidental resonance technology under strong noise background, the multicycle random waveform to be recovered.
Technical background
Traditional low-pass filter can make signal distortion largely (shown in Fig. 1～5) occur after to the signal filtering under the certain noise intensity, and Fig. 1 is the primary standard sine wave, and its amplitude is 0.2V, and frequency is 100Hz; Fig. 2 is that effective value is the white noise of 0.35V; Fig. 3 is the mixed waveform signal figure behind Fig. 2 very noisy stacking diagram 1 sine wave; Fig. 4 is this mixed signal signal graph after through the 500Hz low-pass filter; Fig. 5 is this mixed signal signal graph after through the 1500Hz low-pass filter.Can find out that by Fig. 1～5 bigger distortion can occur after the signal filtering under the certain noise intensity, this will bring great difficulty to the extraction than weak signal, so on engineering, need to study a kind of method and apparatus that can carry out the signal recovery at present.
The accidental resonance notion is to propose when people such as Benzi studied ancient meteorological glacier problem in 1981, after that, accidental resonance amplifies, detects the peculiar advantages in aspect such as identification, transmission reduction in enhancing of signal Processing with it, more and more be applied to multidisciplinary among.In recent years, the accidental resonance technology also causes in engineering fields such as machinery, electronics, information widely to be paid close attention to, and begin to be widely used in engineering aspects such as fault diagnosis, image processing, target following, in the signal Processing field, accidental resonance is mainly used to study characteristics such as the detection, amplification, transmission, recovery of faint information.For accidental resonance, it is single frequency sinusoidal signal s (t)=Asin (2 π f that its signal resonance process can be regarded as in driving force
0T) and under the effect of white noise n (t), particle equation U (x)=-1/2ax
2+ 1/4bx
4The bistable potential well in shake motion; And the process that its signal recovers can be regarded as through the particle behind the bistable equation U (x) by h (x)=U ' (x)=-ax+bx
3The result who moves in the system of decision.Two procedures systems are called " bistable system U (x) " and " recovery system h (x) ".See that on the whole resonance and rejuvenation come down to the motion result that particle passes through " bistable system U (x) " and " recovery system h (x) " in turn.
Therefore, thereby resonance and recovery principle how accurately to use accidental resonance recover to detect the original signal that is mingled in the certain noise without distortion, to obtain real feeble signal, and this signal to be carried out that subsequent analysis uses be the difficult point place of breaking through this technology, and the method also has extremely strong practical value to the input analysis under strong noise background simultaneously.
Summary of the invention
At the problems referred to above, the purpose of this invention is to provide a kind of Weak characteristic information recovery system, for meticulous information analysis and signal Processing provide a kind of new information detection technology based on accidental resonance.
Below in conjunction with accompanying drawing technical scheme of the present invention is illustrated.Weak characteristic information recovery system based on accidental resonance comprises sensor, bistable mimic channel, recovers mimic channel and filtration module etc.Detected noisy original weak signal inputs to the bistable mimic channel through sensor, realizes that the accidental resonance of signal extracts, and recovery mimic channel and follow-up filtration module thereof are finished recovery and the low-pass filtering to signal.Original weak signal in the sensor measured signal is flooded fully by noise and can not be differentiated (as Figure 11 or Figure 16 or shown in Figure 21), and this measured signal is sent into the bistable mimic channel as original sampling data.The effect of bistable mimic channel is that signals and associated noises is carried out accidental resonance.The signal of bistable mimic channel output carries out the signal recovery by recovering mimic channel.Recover mimic channel and mainly form (Fig. 8) by three multiplier chips, amplifier and direct supplys.The effect that recovers mimic channel is that the signal that will produce Stochastic Resonance Phenomenon carries out the denoising recovery, and its circuit structure is: two input ends of first multiplier link to each other with the output terminal of bistable mimic channel respectively; An input end of second multiplier links to each other with the output terminal of first multiplier, and another input end also links to each other with the output terminal of bistable mimic channel; An input end of the 3rd multiplier still links to each other with the output terminal of bistable mimic channel, and another links to each other with direct supply.The output terminal of second multiplier and the 3rd multiplier is received feedback amplifier simultaneously, and the signal of the output terminal of second multiplier and the 3rd multiplier being exported by feedback amplifier feeds back output.This kind connected mode is finished as previously mentioned " recovery system " middle h's (x)-ax+bx
3Connect, and the signal that produces Stochastic Resonance Phenomenon is recovered, the output signal input filtration module that recovers mimic channel subsequently carries out low-pass filtering, finally draws the original weak signal (as Figure 13 or Figure 18 or shown in Figure 23) that is mingled in the noise.Figure 13, Figure 18 are the same substantially with Fig. 9, Figure 14 and Figure 19 waveform respectively with Figure 23 as can be seen simultaneously, illustrate that effect of the present invention is very good.
The signal of sensor actual measurement can be noisy vibration signal, force signal etc., after the bistable mimic channel carries out accidental resonance, signal enters bistable state, signal enters the recovery mimic channel afterwards, after 3 multiplier acting in conjunction outputs through recovering mimic channel, obtain pure signal through the output of filtration module low-pass filtering again.
Description of drawings
Fig. 1 is the primary standard sinusoidal pattern.
Fig. 2 is that effective value is the white noise of 0.35V.
Fig. 3 is the mixed waveform signal figure behind Fig. 2 very noisy stacking diagram 1 sine wave.
Fig. 4 is the signal graph behind the mixed signal process 500Hz low-pass filter.
Fig. 5 is the signal graph behind the mixed signal process 1500Hz low-pass filter.
Fig. 6 is the catenation principle figure of system of the present invention.
Fig. 7 is a bistable mimic channel schematic diagram.
Fig. 8 is for recovering the mimic channel schematic diagram.
Fig. 9 is an initial sinusoids signal waveform curve map.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 200mV.Sinusoidal wave amplitude is 0.2V, and the cycle is 0.01s.
Figure 10 is the oscillogram of original very noisy.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.Noise intensity effective value 0.35V.
Figure 11 is the oscillogram behind Figure 10 very noisy stacking diagram 9 sine waves.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.
Figure 12 is the sinusoidal signal squiggle figure through system recovery output.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.
Figure 13 is the sine waveform curve map of Figure 12 through exporting after the filtering.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 200mV.
Figure 14 is original FM signal squiggle figure.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 200mV.
Figure 15 is the oscillogram of original very noisy.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.Noise intensity effective value 0.35V.
The oscillogram of Figure 16 after for Figure 15 very noisy stack Figure 14 frequency-modulated wave.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.
Figure 17 is the fm waveform curve map through system recovery output.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.
Figure 18 is the fm waveform curve map of Figure 17 through exporting after the filtering.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 200mV.
Figure 19 is original amplitude-modulated signal squiggle figure.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 200mV.
Figure 20 is the oscillogram of original very noisy.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.Noise intensity effective value 0.35V.
The oscillogram of Figure 21 after for Figure 20 very noisy stack Figure 19 modulated wave.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.
Figure 22 is the modulated wave squiggle figure through system recovery output.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 500mV.
Figure 23 is the modulated wave squiggle figure of Figure 22 through exporting after the filtering.Among the figure: horizontal ordinate is a time shaft, and its scale is 10ms; Ordinate is a magnitude of voltage, and its scale is 200mV.
Embodiment
Below in conjunction with accompanying drawing and by embodiment technical scheme of the present invention is described further.
The connectivity scenario of device of the present invention is: be connected in series (as shown in Figure 6) by described order by sensor 1, bistable mimic channel 2, recovery mimic channel 3, filtration module 4.The mimic channel that bistable mimic channel 2 is made up of operational amplifier, multiplier etc. is realized (as shown in Figure 7), and wherein s (t) is the input of bistable system, and x (t) is the output of bistable system.Recovering mimic channel 3 can realize (as shown in Figure 8) by the mimic channel that analog multiplier, feedback amplifier etc. is formed, and its circuit is connected to: two input ends of first multiplier 5-1 link to each other with the output terminal of bistable mimic channel 2 respectively; The input end of second multiplier 5-2 links to each other with the output terminal of first multiplier 5-1, and another input end also links to each other with the output terminal of bistable mimic channel 2; The input end of the 3rd multiplier 5-3 links to each other with the output terminal of bistable mimic channel 2, another input end links to each other with direct supply 6, the output terminal of second multiplier 5-2 and the 3rd multiplier 5-3 is received feedback amplifier 7 simultaneously, and by feedback amplifier 7 it is fed back output.Wherein h (t) is that the input of recovery mimic channel also is the output of bistable mimic channel, and y (t) is the output that recovers mimic channel.
Its parameter of bistable mimic channel 2 (Fig. 7) is: R=15K Ω, and R1=R2=150K Ω, C=180pf, the multiplier in the bistable mimic channel is identical with 3 multiplier architectures in recovering mimic channel, and the product of its coefficient is taken as 0.001.The accidental resonance state of signal is regulated generation by voltage divider K1, K2, and its intrinsic standoff ratio is respectively 0.4 and 0.2.Recover mimic channel 3 as shown in Figure 8, first is adjustable respectively to the 3rd multiplier 5-1～5-3 parameter, and its coefficient is respectively 0.05,0.1,0.02 in order; Direct supply 6 voltages are 1.5V, and used resistance is 15K Ω in the recovery mimic channel 3.Original signal, noise signal, sensor 1, the output signal of recovering all steps such as mimic channel 3 and filtration module 4 are carried out signals collecting by the NI data acquisition system (DAS) and are uploaded to and carry out waveform on the PC and show.
Embodiment 1: the signal that sensor 1 is surveyed is made up of a sine wave signal (as shown in Figure 9) and very noisy (as shown in figure 10), and sine wave signal is flooded fully by noise and can not differentiate (as shown in figure 11).The amplitude of sine wave signal is 0.2V, and the cycle is 0.01s, and the intensity of noise is 0.35V, carries out accidental resonance by 2 pairs of these signals and associated noises of bistable mimic channel, and its output signal is imported recovery mimic channel 3 subsequently and carried out the waveform recovery.Output signal after the recovery (Figure 12) input filtration module carries out 200Hz (butterworth) low-pass filtering (for embodying the integrality of waveform to greatest extent, the passband of low-pass filter is made as 5KHz) and draws net result, as shown in figure 13.
Embodiment 2: the signal that sensor 1 is surveyed is made up of a FM signal (as shown in figure 14) and very noisy (as shown in figure 15), and the frequency-modulated wave signal is flooded fully by noise and can not differentiate (as shown in figure 16).In the frequency-modulated wave signal, carrier amplitude is 0.2V, and carrier cycle is 0.01s; The frequency of modulation signal is 0.1s, and its degree of modulation is that 5 (the FM signal expression formula is U
s=0.2cos (200 π+cos20 π t) t).The intensity of noise is 0.35V.Carry out accidental resonance by 2 pairs of these signals and associated noises of bistable mimic channel, its output signal inputs to recovery mimic channel 3 subsequently and carries out the waveform recovery.Output signal after the recovery (Figure 17) input filtration module carries out 200Hz (butterworth) low-pass filtering (for embodying the integrality of waveform to greatest extent, the passband of low-pass filter is made as 5KHz) and draws net result, as shown in figure 18.
Embodiment 3: the signal that sensor 1 is surveyed is made up of an amplitude-modulated signal (as shown in figure 19) and very noisy (as shown in figure 20), and amplitude-modulated signal is flooded fully by noise and can not differentiate (as shown in figure 21).In amplitude-modulated signal, carrier amplitude is 0.2V, and carrier cycle is 0.01s; The frequency of modulation signal is 0.1s, and its degree of modulation is 1, and (the amplitude-modulated signal expression formula is U
s=(0.2+0.2cos20 π t) cos200 π t).The intensity of noise is 0.35V.Carry out accidental resonance by 2 pairs of these signals and associated noises of bistable mimic channel, its output signal is imported recovery mimic channel 3 subsequently and is carried out the waveform recovery.Output signal after the recovery (Figure 22) input filtration module 4 carries out 200Hz (butterworth) low-pass filtering (for embodying the integrality of waveform to greatest extent, the passband of low-pass filter is made as 5KHz) and draws net result, as shown in figure 23.
The effect of filtration module 4 is mainly used in to carrying out low-pass filtering by the signal that recovers mimic channel 3 outputs, so that the signal to noise ratio (S/N ratio) of output signal is further enhanced.
Characteristics of the present invention and useful effect are: this system has solved the meticulous test problems of weak signal in the very noisy, Break through in the past low pass filter and can realize the multicycle, many to the limitation of the information function of weak signal in the very noisy The functions such as frequency, multicomponent complicated Weak characteristic information identification and signal reduction.
Claims (3)
1. based on the Weak characteristic information recovery system of accidental resonance, has sensor (1), bistable mimic channel (2), recover mimic channel (3) and filtration module (4), it is characterized in that by sensor (1), bistable mimic channel (2), recovering mimic channel (3) and filtration module (4) is connected in series by described order, signal by sensor (1) actual measurement is sent into bistable mimic channel (2) successively as raw data, recover mimic channel (3) and filtration module (4), draw the original weak signal with practical value to be extracted by filtration module (4), recover mimic channel (3) and be made up of 3 analog multiplier chips and direct supply (6) and feedback amplifier (7), two input ends of first multiplier (5-1) chip link to each other with the output terminal of bistable mimic channel (2); An input end of second multiplier (5-2) chip links to each other with the output terminal of first multiplier 5-1, and another input end is connected to the output terminal of bistable mimic channel (2); An input end of the 3rd multiplier (5-3) chip links to each other with the output terminal of bistable mimic channel (2), another input end links to each other with direct supply (6), second multiplier (5-2) received feedback amplifier (7) simultaneously with the output terminal of the 3rd multiplier (5-3), feeds back output by feedback amplifier (7).
2. according to the described Weak characteristic information recovery system based on accidental resonance of claim 1, the coefficient value scope that it is characterized in that described recovery mimic channel (3) is to described first to the 3rd analog multiplier chip (5-1～5-3) be 0～1.
3. according to the described Weak characteristic information recovery system of claim 1, it is characterized in that described filtration module (4) adopts low-pass filter based on accidental resonance.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100687840A CN101561293B (en) | 2009-05-11 | 2009-05-11 | Weak characteristic information recovery system based on stochastic resonance |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2009100687840A CN101561293B (en) | 2009-05-11 | 2009-05-11 | Weak characteristic information recovery system based on stochastic resonance |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101561293A true CN101561293A (en) | 2009-10-21 |
| CN101561293B CN101561293B (en) | 2011-09-07 |
Family
ID=41220191
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2009100687840A Expired - Fee Related CN101561293B (en) | 2009-05-11 | 2009-05-11 | Weak characteristic information recovery system based on stochastic resonance |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101561293B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101860347A (en) * | 2010-06-02 | 2010-10-13 | 天津大学 | Stochastic resonance signal recovery method based on signal classification |
| CN105492929A (en) * | 2014-02-27 | 2016-04-13 | 富士电机株式会社 | Radiation detection device, radiation dosimetry processing method, and radiation dosimetry processing program |
| CN107246912A (en) * | 2017-06-22 | 2017-10-13 | 西北工业大学 | A kind of marine riser vortex-induced vibration monitoring method based on accidental resonance |
| CN108306838A (en) * | 2017-01-11 | 2018-07-20 | 联发科技（新加坡）私人有限公司 | Demodulator and the method for modulating input signal for solving modulation |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6008642A (en) * | 1997-08-25 | 1999-12-28 | The United States Of America As Represented By The Secretary Of The Navy | Stochastic resonance detector for weak signals |
| US6285249B1 (en) * | 2000-01-21 | 2001-09-04 | The United States Of America As Represented By The Secretary Of The Navy | Controlled stochastic resonance circuit |
| CN1595179A (en) * | 2004-07-09 | 2005-03-16 | 天津大学 | Cascade double stabilization system for information detection |
-
2009
- 2009-05-11 CN CN2009100687840A patent/CN101561293B/en not_active Expired - Fee Related
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101860347A (en) * | 2010-06-02 | 2010-10-13 | 天津大学 | Stochastic resonance signal recovery method based on signal classification |
| CN101860347B (en) * | 2010-06-02 | 2012-11-14 | 天津大学 | Stochastic resonance signal recovery method based on signal classification |
| CN105492929A (en) * | 2014-02-27 | 2016-04-13 | 富士电机株式会社 | Radiation detection device, radiation dosimetry processing method, and radiation dosimetry processing program |
| CN105492929B (en) * | 2014-02-27 | 2017-06-13 | 富士电机株式会社 | Radiation detecting apparatus and radiation dose measurement processing method |
| CN108306838A (en) * | 2017-01-11 | 2018-07-20 | 联发科技（新加坡）私人有限公司 | Demodulator and the method for modulating input signal for solving modulation |
| CN108306838B (en) * | 2017-01-11 | 2020-12-11 | 联发科技（新加坡）私人有限公司 | Demodulator and method for demodulating an amplitude-modulated input signal |
| CN107246912A (en) * | 2017-06-22 | 2017-10-13 | 西北工业大学 | A kind of marine riser vortex-induced vibration monitoring method based on accidental resonance |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101561293B (en) | 2011-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101561293B (en) | Weak characteristic information recovery system based on stochastic resonance | |
| CN102967759A (en) | Weak signal detection circuit based on chaotic oscillator | |
| CN104545887A (en) | Method and device for identifying artifact electrocardiograph waveforms | |
| CN102410781A (en) | Detection device and detection method for frequency of artillery launchings | |
| CN106323448A (en) | High-frequency vibration distributed optical fiber sensing system based on random sampling | |
| CN107153139A (en) | A kind of optical fibre gyro detection circuit board noise characteristic method of testing | |
| CN109297705A (en) | Epicyclic gearbox vibration signal method for diagnosing faults based on MED and fuzzy entropy | |
| CN102914432A (en) | Vibration and acoustic transmission signal regulation instrument for detecting mechanical faults | |
| CN203054093U (en) | Apparatus for extracting DC (direct current) signal from mixed signal by using synchronous sampling | |
| CN107193374A (en) | A kind of detection means and detection method of the intentional gesture motion of active | |
| CN102670182B (en) | Pulse wave acquisition quality real-time analysis device | |
| CN104035012A (en) | Partial-discharge ultrahigh-frequency-signal detection conditioning circuit | |
| CN106770522A (en) | A kind of electronic tongue system and detection method for detecting fruit juice purity | |
| CN103713159A (en) | Feedback method of closed-loop micromechanics accelerometer | |
| CN103472485B (en) | Array frequency domain electromagnetic detection system quality of data monitoring device and monitoring method | |
| CN106097723A (en) | Vehicle based on bias magnetic field is delayed unloading the method for work of reality line detector | |
| CN201426110Y (en) | Stochastically resonated faint feature information retrieval device | |
| CN108896454A (en) | A kind of multichannel wear particle detection method and device based on time-division multiplex technology | |
| CN101726542B (en) | Method and system for monitoring and identifying ultrasonic bonding quality in package and interconnection of chips on line | |
| CN104714127B (en) | The anti-traction current interference test device of track circuit based on virtual instrument | |
| CN102818577A (en) | Automatic test device for frequency characteristic of laser gyro electromechanical assembly | |
| CN104132690A (en) | Multi-sensor test system | |
| CN106813784A (en) | A kind of real-time microwave pulse chirp detection means and its detection method | |
| CN202994625U (en) | Weak optical signal detection circuit of liquid particle counter | |
| CN108983065A (en) | A kind of voltage amplifier circuit, detection circuit and its circuit detecting method |
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 | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110907 Termination date: 20120511 |