CN102591204A - Realization method of Gauss system - Google Patents

Realization method of Gauss system Download PDF

Info

Publication number
CN102591204A
CN102591204A CN2012100450610A CN201210045061A CN102591204A CN 102591204 A CN102591204 A CN 102591204A CN 2012100450610 A CN2012100450610 A CN 2012100450610A CN 201210045061 A CN201210045061 A CN 201210045061A CN 102591204 A CN102591204 A CN 102591204A
Authority
CN
China
Prior art keywords
omega
gauss
infin
order
cascade
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.)
Pending
Application number
CN2012100450610A
Other languages
Chinese (zh)
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.)
HARBIN HIGHTECH PRECISION GAGE CO Ltd
Original Assignee
HARBIN HIGHTECH PRECISION GAGE CO Ltd
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 HARBIN HIGHTECH PRECISION GAGE CO Ltd filed Critical HARBIN HIGHTECH PRECISION GAGE CO Ltd
Priority to CN2012100450610A priority Critical patent/CN102591204A/en
Publication of CN102591204A publication Critical patent/CN102591204A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Complex Calculations (AREA)

Abstract

The invention discloses a realization method of a Gauss system, belonging to the measurement technology. Finite terms of the expansion of a Gauss prototype function are intercepted as an approximant of a Gauss function; a finite number of first-order systems with the same characteristics are cascaded; the system characteristics approximate the Gauss characteristics; when the cascade series of the first-order systems is gradually increased, the error of first-order system characteristics is gradually decreased compared with that of Gauss system characteristics; along with the increase of the cascade series, the approximation precision to the Gauss system is improved; the method is also suitable for a second-order system; with the instruction of the method, a Gauss system can be realized through the cascade approximation of the first-order systems and second-order systems with the same characteristics; and compared with the traditional system, the Gibbs phenomenon generated by the cut-off frequency is eliminated.

Description

The implementation method of a kind of Gauss system
Technical field
The invention belongs to the measuring system technical field, relate generally to the implementation method of a kind of Gauss system.
Background technology
In order to realize the undistorted transmission of signal, require system to have infinitely-great bandwidth and linear phase characteristic in theory, promptly desirable instrument should satisfy two conditions: 1) amplitude versus frequency characte is the constant with frequency-independent, 2) phase-frequency characteristic and frequency are linear.In fact, theoretic desirable instrument is impossible realize physically.Generally speaking; The bandwidth of actual signal all is limited; Even the spectrum distribution of signal can be extended to infinitely great frequency place, but owing to the contained energy of each frequency component along with increasing of frequency reduces, can ignore by the distortion that higher frequency components causes.System is transmitted undistortedly for assurance, if system can satisfy 1) amplitude versus frequency characte is a constant in certain frequency range, 2) phase-frequency characteristic and frequency are linear in this frequency range, and then this system also can be used as the instrument with undistorted characteristic.
For Design of Measurement System, people have carried out sufficient research to traditional desirable lowpass system model with rectangle amplitude-frequency response etc., and its time domain specification, frequency domain characteristic etc. had deep understanding.These ideal system models are non-causal systems; Physically can not realize; But available method of approaching is similar to realization; The system that is realized has produced Gibbs' effect because the decay at cut-off frequency place sharply changes, and promptly has overshoot and oscillatory occurences at the cutoff frequency place, and this is disadvantageous for the stable of measuring system.Desirable Gauss's system model has level and smooth step response characteristic, and the time frequency range long-pending minimum, therefore, can be used as the target of optimization system design.
Summary of the invention
The objective of the invention is in order to solve above-mentioned existing problem; On the basis of further investigation measuring system theory; With system with Gaussian response characteristic target as the optimization system design; Propose the method that a kind of Gauss system realizes, the system that is realized does not exist because the decay at cut-off frequency place sharply changes and has produced Gibbs' effect.
The objective of the invention is to realize like this:
The implementation method of a kind of Gauss system, exponential function e xExpansion
e x = Σ n = 0 ∞ x n n ! = 1 + x + x 2 2 ! + . . . + x n n ! + . . . , x ∈ ( - ∞ , + ∞ ) ,
Have in view of the above
e x 2 = Σ n = 0 ∞ ( x 2 ) n n ! = 1 + x 2 + ( x 2 ) 2 2 ! + . . . + ( x 2 ) n n ! + . . . , x ∈ ( - ∞ , + ∞ ) ,
The different power items of the above-mentioned function of intercepting
Figure BSA00000675174400013
expansion are as Gaussian function approximating function formula; There is following relational expression to set up for object function
(1) for first-order system
lim n → ∞ [ 1 1 + ( β n Ω / Ω c ) 2 ] n = e - ( αΩ / Ω c ) 2 ,
Employing has cascade such as sensor, hardware filtering device, software filter and amplifying circuit, the detecting circuit of first-order system characteristic and approaches the Gauss system that realizes;
(2) for second-order system
lim n → ∞ [ 1 1 + ( β n Ω / Ω c ) 2 + ( β n Ω / Ω c ) 4 2 ] n = e - ( αΩ / Ω c ) 2 ,
Employing has cascade such as sensor, hardware filtering device, software filter and amplifying circuit, the detecting circuit of second-order system characteristic and approaches the Gauss system that realizes.
The invention has the advantages that: the implementation method of Ben Gaosi system is compared with traditional system, has eliminated because the Gibbs phenomenon that cut-off frequency produces.
Description of drawings
Fig. 1 is the amplitude versus frequency characte of n level first-order system cascade system and Gauss system
Fig. 2 is the amplitude versus frequency characte deviation curve of n level first-order system cascade system and Gauss system
Fig. 3 is the system responses characteristic of first-order system prototype
Fig. 4 is 8 grades of first-order system cascade system response characteristics
Fig. 5 is the amplitude versus frequency characte of n level second-order system cascade system and Gauss system
Fig. 6 is the amplitude versus frequency characte deviation curve of n level second-order system cascade system and Gauss system
Fig. 7 is the system responses characteristic of second-order system prototype
Fig. 8 is 4 grades of second-order system cascade system response characteristics
Fig. 9 is the sensor signal processing system schematic diagram
Figure 10 is a joint RC active filter circuit
Figure 11 is-the 2RC filtering circuit
Figure 12 realizes the schematic diagram of Gauss system for cascade
Embodiment:
For geometric measurement; When the inductance type micro-displacement sensor that uses is first-order system or second-order system; For such sensing system, when forming measuring system according to the basic physical model of geometric sense multimetering instrument, the system how to realize ideal with Gaussian response characteristic; Seem particularly important, carry out detailed analysis and study in the face of this problem down.
According to e xExpansion
e x = Σ n = 0 ∞ x n n ! = 1 + x + x 2 2 ! + . . . + x n n ! + . . . , x ∈ ( - ∞ , + ∞ ) ,
Have
e x 2 = Σ n = 0 ∞ ( x 2 ) n n ! = 1 + x 2 + ( x 2 ) 2 2 ! + . . . + ( x 2 ) n n ! + . . . , x ∈ ( - ∞ , + ∞ ) ,
The different power items of the above-mentioned function expansion of intercepting are as Gaussian function approximating function formula; According to central-limit theorem; For object function:
Figure BSA00000675174400034
should have following relational expression to set up
For first-order system
lim n → ∞ [ 1 1 + ( β n Ω / Ω c ) 2 ] n = e - ( αΩ / Ω c ) 2 - - - ( 1 )
The cascades such as sensor, amplifying circuit, detecting circuit, hardware filtering device and software filter that employing has the first-order system characteristic approach the Gauss system that realizes;
For second-order system
lim n → ∞ [ 1 1 + ( β n Ω / Ω c ) 2 + ( β n Ω / Ω c ) 4 2 ] n = e - ( αΩ / Ω c ) 2 - - - ( 2 )
The cascades such as sensor, amplifying circuit, detecting circuit, hardware filtering device and software filter that employing has the second-order system characteristic approach the Gauss system that realizes.
Instance 1 first-order system cascade approaches the Gauss system that realizes
1.1 the first-order system cascade approaches the central limit theorem of Gauss system
By aforementioned analysis, suppose that used sensor is a first-order system, and hope that the characteristic of other part all is a single order, when a lot of identical first-order system cascades, square can be described as of the amplitude versus frequency characte function of its system
| H ( Ω / Ω c ) | n = ( 1 1 + ( β n Ω / Ω c ) 2 ) n
So, can use the first-order system cascade to realize the Gauss system, will solve the represented limiting equation of formula (1)
lim n → ∞ [ 1 1 + ( β n Ω / Ω c ) 2 ] n = e - ( αΩ / Ω c ) 2 ,
Whether set up, this is the key that realizes the Gauss system.Proof procedure is following.
Proof: order
G ( Ω / Ω c ) = ( 1 1 + ( β n Ω / Ω c ) 2 ) n
At cutoff frequency Ω=Ω cThe damping capacity at place is G (Ω/Ω cDuring)=p, have
( 1 1 + ( β n ) 2 ) n = p
1 + ( β n ) 2 = p - 1 n
β n = p - 1 n - 1 - - - ( 3 )
If
q=(β nΩ/Ω c) 2
Then
G ( Ω / Ω c ) = ( 1 1 + q ) n
If
β n 2 = p - 1 n - 1 = 1 t
Then
q = 1 t · ( Ω / Ω c ) 2
n = ln 1 p ln ( 1 + 1 t ) - - - ( 4 )
When n →+∞, promptly t →+during ∞,
lim t → + ∞ G ( Ω / Ω c ) = lim t → ∞ ( 1 1 + 1 t ( Ω / Ω c ) 2 ) ln 1 p ln ( 1 + 1 t )
= lim t → ∞ ( 1 + 1 t ( Ω / Ω c ) 2 ) - t t · ln 1 p ln ( 1 + 1 t )
= lim t → ∞ ( 1 + 1 t ( Ω / Ω c ) 2 ) - t · ln 1 p t · ln ( 1 + 1 t )
= lim t → ∞ ( 1 + 1 t ( Ω / Ω c ) 2 ) - t ln 1 p ln ( 1 + 1 t ) t
= lim t → ∞ ( 1 + 1 t ( Ω / Ω c ) 2 ) 1 1 t ( Ω / Ω c ) 2 · 1 t ( Ω / Ω c ) 2 · ( - t ) · ln 1 p ln ( 1 + 1 t ) t
= lim t → ∞ ( 1 + 1 t ( Ω / Ω c ) 2 ) 1 1 t ( Ω / Ω c ) 2 · ( - 1 ) · ( Ω / Ω c ) 2 · ln 1 p ln ( 1 + 1 t ) t - - - ( 5 )
Turn to according to the limit
Figure BSA00000675174400055
formula (5):
lim t → + ∞ G ( Ω / Ω c ) = e - ln 1 p · ( Ω / Ω c ) 2 = e - ( αΩ / Ω c ) 2
In the formula--- α = Ln 1 p
The first-order system cascade approaches the central limit theorem of Gauss system and must demonstrate,prove.Formula (1) shows the cascade of infinite a plurality of first-order systems, and its system performance is just ad infinitum approached gaussian characteristics.Also promptly, limited cascade with first-order system of identical characteristics, its system performance is that of gaussian characteristics is approximate.
According to With
Figure BSA00000675174400059
When making p=0.707, get α=0.5887, the cascade progression n=1 for first-order system, has β at 2,4,8,16 o'clock 1=0.6436, β 2=0.435, β 4=0.3008, β 8=0.2104, β 16=0.148, make the cascaded functions of first-order system
Figure BSA000006751744000510
Gauss's system function Amplitude-frequency response, as shown in Figure 1.
The amplitude error of first-order system cascade systems at different levels and Gauss system amplitude-frequency response is as shown in Figure 2, and its amplitude peak deviation is listed in table 1.
The amplitude peak deviation that the cascade of table 1n level first-order system approaches the Gauss system
Figure BSA000006751744000512
Can know that by Fig. 2 and table 1 when the cascade progression of first-order system increased gradually, its system performance and Gauss's system performance ratio error were mutually reducing gradually, this means the increase along with cascade progression, with the approximation accuracy that improves the Gauss system.
1.2 single order simulation system prototype and system responses characteristic thereof
When getting 8 grades of first-order system cascades,
G 8 ( Ω ) = 1 1 + ( β 8 Ω / Ω c ) 2
By the amplitude square principle of design:
G 8(Ω)=|H a8(jΩ)| 2=H a8(jΩ)H a8(-jΩ)
With s=j Ω substitution following formula, then have:
G 8 ( s ) = H a 8 ( s ) H a 8 ( - s ) = 1 1 - ( β 8 2 / Ω c 2 ) s 2
Its limit is:
Figure BSA00000675174400063
Figure BSA00000675174400064
Choose the limit s of the left half-plane on s plane 2As the limit of simulation system function H (s), establishing gain constant is K 0, obtain
H a 8 ( s ) = K 0 s - s 2 - - - ( 6 )
By H A8(s) | S=0=H (j Ω) | Ω=0 solves
Figure BSA00000675174400066
H A8(s) be design Gauss and approach the H of system β 8(s) single order simulation system prototype.
Single order simulation system prototype H A8(s) unit impulse response
h 8 ( t ) = L - 1 [ K 0 s - s 2 ]
= K 0 e s 2 t - - - ( 7 )
Unit-step response
y ( t ) = L - 1 [ 1 s K 0 s - s 2 ] (8)
= - s 2 K 0 ( 1 - e s 2 t )
Make the unit impact response and the unit-step response curve of single order simulation system according to formula (7) and formula (8), as shown in Figure 3.
1.3 single order simulation system prototype cascade and system responses characteristic thereof
8 grades of single orders are simulated when system-level, the unit impulse response of system
h 8 ( t ) = L - 1 [ ( K 0 s - s 2 ) 8 ]
= K 0 8 7 ! t 7 e s 2 s 2 t - - - ( 9 )
Unit-step response
y ( t ) = ∫ 0 t h 8 ( x ) dx
= K 0 8 7 ! ∫ 0 t x 7 e s 2 x dx
Then
y ( t ) = K 0 8 7 ! ( s 2 7 t 7 e s 2 t - 7 s 2 6 t 6 e s 2 t + 42 s 2 5 t 5 e s 2 t - 210 s 2 4 t 4 e s 2 t + (10)
840 s 2 3 t 3 e s 2 t - 2520 s 2 2 t 2 e s 2 t + 5040 s 2 te s 2 t - 5040 e s 2 t + 5040 ) / s 2 8
According to formula (9) and formula (10), make 8 grades of first-order system cascaded system unit impact response curves and unit-step response curve respectively, as shown in Figure 4.
Can know that by Fig. 4 the unit impulse response of the cascade system of 8 grades of first-order systems and unit-step response family curve very approach the resonse characteristic of Gauss system.
Instance 2 second-order system cascades approach the Gauss system that realizes
2.1 the second-order system cascade approaches the central limit theorem of Gauss system
When supposing that used sensor is second-order system, if will use a plurality of second-order systems to realize the Gauss system, then demand is demonstrate,proved the limit
lim n → ∞ [ 1 1 + ( βΩ / Ω c ) 2 + ( βΩ / Ω c ) 4 2 ] n = e - ( αΩ / Ω c ) 2
Set up.Issued a certificate below.
Proof: order
G ( Ω / Ω c ) = ( 1 1 + ( β n Ω / Ω c ) 2 + 1 2 ( β n Ω / Ω c ) 4 ) n
Make cutoff frequency Ω=Ω cDamping capacity G (Ω/the Ω at place c)=p when (0<p<1), then has
( 1 + ( β n ) 2 + 1 2 ( β n ) 4 ) - n = p
1 + ( β n ) 2 + 1 2 ( β n ) 4 = p - 1 n
β n = 2 × p - 1 n - 1 - 1 - - - ( 11 )
If
q=(β nΩ/Ω c) 2
Then
G ( Ω / Ω c ) = ( 1 1 + q + q 2 / 2 ) n
If β n 2 = 2 × p - 1 n - 1 - 1 = 2 × p - 1 n - 1 - 1 = 1 t
Then
q = 1 t · ( Ω / Ω c ) 2
n = ln 1 p ln [ ( 1 + 1 t ) 2 + 1 2 ] - - - ( 12 )
When n →+∞, promptly t →+during ∞,
lim t → + ∞ G ( Ω Ω c ) = lim t → ∞ ( 1 1 + 1 t ( Ω Ω c ) 2 + 1 2 t 2 ( Ω Ω c ) 4 ) ln 1 p ln [ ( 1 t + 1 ) 2 + 1 2 ]
= lim t → ∞ ( 1 + 1 t ( Ω Ω c ) 2 + 1 2 t 2 ( Ω Ω c ) 4 ) - t t · ln 1 p ln [ ( 1 t + 1 ) 2 + 1 2 ]
= lim t → ∞ ( 1 + 1 t ( Ω Ω c ) 2 + 1 2 t 2 ( Ω Ω c ) 4 ) - t · ln 1 p t · ln [ ( 1 t + 1 ) 2 + 1 2 ]
= lim t → ∞ ( 1 + 1 t ( Ω Ω c ) 2 + 1 2 t 2 ( Ω Ω c ) 4 ) - t · ln 1 p ln [ ( 1 t + 1 ) 2 + 1 2 ] t - - - ( 13 )
Can derive obtains by
Figure BSA00000675174400091
:
lim t → + ∞ ln [ ( 1 + 1 t ) 2 + 1 2 ] t = lim t → ∞ ln ( 1 + 1 + 2 t 2 t 2 ) 2 t 2 1 + 2 t · 1 + 2 t 2 t
= lim t → + ∞ ln ( e 1 + 2 t 2 t )
= lim t → + ∞ 1 + 2 t 2 t
= lim t → + ∞ 1 + 1 2 t
= 1
Therefore, formula (13) turns to
lim t → + ∞ G ( Ω / Ω c ) = lim t → ∞ ( 1 + 1 t ( Ω / Ω c ) 2 + 1 2 t 2 ( Ω / Ω c ) 4 ) - t · ln 1 p
= lim t → ∞ ( [ 1 + 1 t ( Ω / Ω c ) 2 ] · [ 1 + 1 2 t 2 ( Ω / Ω c ) 4 1 + 1 t ( Ω / Ω c ) 2 ] ) - t · ln 1 p
= lim t → ∞ ( 1 + 1 t ( Ω / Ω c ) 2 ) - t · ln 1 p · lim t → ∞ ( 1 + 1 2 t 2 ( Ω / Ω c ) 4 1 + 1 t ( Ω / Ω c ) 2 ) - t · ln 1 p - - - ( 14 )
For the subitem in the formula (14)
lim t → ∞ ( 1 + 1 2 t 2 ( Ω Ω c ) 4 1 + 1 t ( Ω Ω c ) 2 ) - t · ln 1 p
Have
1 + 1 2 t 2 ( Ω Ω c ) 4 > 1 + 1 2 t 2 ( Ω Ω c ) 4 1 + 1 t ( Ω Ω c ) 2 > 1 , t > 0
And
lim t → ∞ ( 1 + 1 2 t 2 ( Ω Ω c ) 4 ) - t · ln 2 2 = lim t → ∞ ( 1 + 1 2 t 2 ( Ω Ω c ) 4 ) 2 t 2 · ( Ω Ω c ) - 4 · 1 2 t 2 · ( Ω Ω c ) 4 · ( - t · ln 2 2 )
= lim t → ∞ e - ln 2 4 t · ( Ω Ω c ) 4
= 1
So Lim t → ∞ ( 1 + 1 2 t 2 ( Ω Ω c ) 4 1 + 1 t ( Ω Ω c ) 2 ) - t · Ln 2 2 = 1
Therefore, formula (14) turns to
lim t → + ∞ G ( Ω Ω c ) = lim t → ∞ ( 1 + 1 t ( Ω Ω c ) 2 ) - t · ln 1 p
= lim t → ∞ ( 1 + 1 t ( Ω Ω c ) 2 ) - t · ( Ω Ω c ) - 2 · ( Ω Ω c ) 2 · ln 1 p
= e - ln 1 p ( Ω Ω c ) 2
= e - ( α Ω Ω c ) 2
In the formula--- α = Ln 1 p
To sum up Lim n → ∞ ( 1 1 + ( β n Ω / Ω c ) 2 + 1 2 ( β n Ω / Ω c ) 4 ) n = e - ( α Ω Ω c ) 2
The second-order system cascade approaches the central limit theorem of Gauss system and must demonstrate,prove.Formula (2) shows the cascade of infinite a plurality of second-order systems, and its system performance is just ad infinitum approached gaussian characteristics.Also promptly, limited cascade with second-order system of identical characteristics, its system performance is gaussian characteristics-individual being similar to.
According to
Figure BSA000006751744001011
With
Figure BSA000006751744001012
When making p=0.707, α=0.5887 is arranged then, work as n=1,2,4,8,16 o'clock, β 1=0.5935, β 2=0.4172, β 4=0.2945, β 8=0.2082, β 16=0.1472, it is more as shown in Figure 5 to map.
The amplitude error of second-order system cascade systems at different levels and Gauss system amplitude-frequency response is as shown in Figure 6, and its amplitude peak deviation is listed in table 2.
The amplitude peak deviation that the cascade of table 2n level second-order system approaches the Gauss system
Can know that by Fig. 6 when the cascade progression of second-order system increased gradually, its characteristic and Gauss's system performance ratio error were mutually reducing gradually, this means the increase along with cascade progression, with the approximation accuracy that improves the Gauss system; And the cascade of second-order system is faster than the cascade velocity of approch of first-order system.
2.2 second order simulation system prototype and system responses characteristic thereof
When getting 4 grades of second-order system cascades,
G 4 ( Ω ) = 1 1 + ( β 4 Ω / Ω c ) 2 + 1 2 ( β 4 Ω / Ω c ) 4
By the amplitude square principle of design [89,90]:
G 4(Ω)=|H a4(jΩ)| 2
=H a4(jΩ)H a4(-jΩ)
With s=j Ω substitution following formula, then have:
G 4 ( s ) = H a 4 ( s ) H a 4 ( - s )
= 1 1 - ( β 4 2 / Ω c 2 ) s 2 + ( β 4 4 / 2 Ω c 4 ) s 4
Its limit is: s 1 = Ω c β 4 ( 1.0987 + j 0.4551 ) , s 2 = Ω c β 4 ( - 1.0987 + j 0.4551 ) , s 3 = Ω c β 4 ( - 1.0987 - j 0.4551 ) , s 4 = Ω c β 4 ( 1.0987 - j 0.4551 ) . Choose the limit s of the left half-plane on S plane 2, s 3As the limit of simulation system function H (s), establishing gain constant is K 0, obtain
H a 4 ( s ) = K 0 ( s - s 2 ) ( s - s 3 ) - - - ( 15 )
By H A4(s) | S=0=H (j Ω) | Ω=0Condition obtain
Figure BSA000006751744001110
H A4Be design Gauss and approach the H of system β 4Second order simulation system prototype.
The unit impulse response of second order simulation system prototype
X(s)=1
Y ( s ) X ( s ) = K 0 ( s - s 2 ) ( s - s 3 )
Y ( s ) = K 0 s 2 - s 3 [ 1 s - s 2 - 1 s - s 3 ]
h ( t ) = K 0 s 2 - s 3 ( e s 2 t - e s 3 t ) - - - ( 16 )
Unit-step response
X ( s ) = 1 s
Y ( s ) X ( s ) = K 0 ( s - s 2 ) ( s - s 3 )
Y ( s ) = 1 s · K 0 ( s - s 2 ) ( s - s 3 )
= K 0 s 2 s 3 1 s + K 0 s 2 ( s 2 - s 3 ) 1 s - s 2 - K 0 s 3 ( s 2 - s 3 ) 1 s - s 3
Then
y ( t ) = K 0 s 2 s 3 + K 0 s 2 ( s 2 - s 3 ) e s 2 t - K 0 s 3 ( s 2 - s 3 ) e s 3 t - - - ( 17 )
According to formula (16) and formula (17), it is as shown in Figure 7 to map.
2.3 second order simulation system prototype cascade and system responses characteristic thereof
The cascade of 4 grades of second order simulation system prototypes is
Figure BSA00000675174400129
its system unit impulse response
Y ( s ) X ( s ) = [ K 0 ( s - s 2 ) ( s - s 3 ) ] 4
Y ( s ) = K 0 4 ( s - s 2 ) 4 ( s - s 3 ) 4
Then
h ( t ) = L - 1 [ K 0 4 ( s - s 2 ) 4 ( s - s 3 ) 4 ]
= K 0 4 t 3 e s 2 t 3 ! * t 3 e s 3 t 3 !
= K 0 4 81 t 3 e s 2 t ( 6 - 6 e s 3 t + e s 3 t s 3 3 t 3 - 3 e s 3 t s 3 2 t 2 - 6 e s 3 t s 3 t ) / s 3 4 - - - ( 18 )
Unit-step response
y ( t ) = ∫ 0 t h ( x ) dx - - - ( 19 )
According to formula (18) and formula (19), it is as shown in Figure 8 to map.
Can know that by Fig. 8 the cascaded system unit impulse response of 4 grades of second-order systems and unit-step response family curve very approach the resonse characteristic of Gauss system.
The physics realization of 3 Gauss systems
For the electronic system of geometric sense multimetering system, sensor signal need be passed through processing such as amplification, detection, filtering, sampling usually, and what influence the system responses characteristic mainly is the realization of wave filter.Typical sensor signal processing system schematic diagram is as shown in Figure 9.
For first-order system, approach the Gauss system that realizes with cascades such as sensor, hardware filtering device, software filter and amplifying circuit with first-order system characteristic, detecting circuits; For second-order system, approach the Gauss system that realizes with cascades such as sensor, hardware filtering device, software filter and amplifying circuit with second-order system characteristic, detecting circuits.
Realize that with first-order system the Gauss system is an example, on the basis that the principle that the cascade of Gauss system is realized is studied, consider to influence the key component of system performance with the realization of RC wave filter.The RC wave filter is a low-pass filter commonly used in the instrument, a joint RC active filter [92]Shown in figure 10.
For a joint RC filtering circuit, its ssystem transfer function
H f 1 = V out V in = jRCω 1 + jRCω
Wherein,
Figure BSA00000675174400134
R and C represent resistance and electric capacity respectively.
The amplitude-frequency function square do
| H f 1 | 2 = 1 1 + ( ω c / ω ) 2
2RC filter circuit construction form is shown in figure 11.
For the 2RC filtering circuit, its ssystem transfer function
Hf 2 = V out V in = Hf 1 2 = ( jRCω 1 + jRCω ) 2
Wherein,
Figure BSA00000675174400137
Figure BSA00000675174400138
R and C represent resistance and electric capacity respectively.
The amplitude-frequency function square do
| Hf 2 | 2 = 1 1 + ( ω c / ω ) 2
For n joint 2RC filtering circuit, have
Hf n = V out V in = Hf 1 n = ( jRCω 1 + JRCω ) n
The amplitude-frequency function square do
| H f n | 2 = [ 1 1 + ( ω c / ω ) 2 ] n
According to formula (1), the n level cascade of then available 2RC wave filter approaches the realization Gaussian filter.Gaussian filter can approach realization through the cascade of multistage Butterworth filter, so just changes into our the very design of familiar Butterworth filter to the design of the Gaussian filter of complicacy [90]Then, the cascade of sensing system, amplification, detection, filtering system according to formula (1), can realize having the system of gaussian characteristics, and it is simpler to use the system architecture that this method obtains, and is prone to realize, and is shown in figure 12.
According to above-mentioned analysis, in the practical application, can obtain satisfied effect with 4 joint 2RC wave filter cascades.The application cascade approaches the sensor and the electronic system that realize the Gauss system and in many cover multimetering instruments, is used widely, and is the basis of supporting geometric sense multimetering technical application.

Claims (1)

1. the implementation method of a Gauss system is characterized in that: exponential function e xExpansion
e x = Σ n = 0 ∞ x n n ! = 1 + x + x 2 2 ! + . . . + x n n ! + . . . , x ∈ ( - ∞ , + ∞ ) ,
Have in view of the above
e x 2 = Σ n = 0 ∞ ( x 2 ) n n ! = 1 + x 2 + ( x 2 ) 2 2 ! + . . . + ( x 2 ) n n ! + . . . , x ∈ ( - ∞ , + ∞ ) ,
The different power items of the above-mentioned function of intercepting
Figure FSA00000675174300013
expansion are as Gaussian function
Figure FSA00000675174300014
approximating function formula; There is following relational expression to set up for object function
(1) for first-order system
lim n → ∞ [ 1 1 + ( β n Ω / Ω c ) 2 ] n = e - ( αΩ / Ω c ) 2 ,
The cascades such as sensor, amplifying circuit, detecting circuit, hardware filtering device and software filter that employing has the first-order system characteristic approach the Gauss system that realizes;
(2) for second-order system
lim n → ∞ [ 1 1 + ( β n Ω / Ω c ) 2 + ( β n Ω / Ω c ) 4 2 ] n = e - ( αΩ / Ω c ) 2 ,
The cascades such as sensor, amplifying circuit, detecting circuit, hardware filtering device and software filter that employing has the second-order system characteristic approach the Gauss system that realizes.
CN2012100450610A 2012-02-27 2012-02-27 Realization method of Gauss system Pending CN102591204A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012100450610A CN102591204A (en) 2012-02-27 2012-02-27 Realization method of Gauss system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012100450610A CN102591204A (en) 2012-02-27 2012-02-27 Realization method of Gauss system

Publications (1)

Publication Number Publication Date
CN102591204A true CN102591204A (en) 2012-07-18

Family

ID=46480040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012100450610A Pending CN102591204A (en) 2012-02-27 2012-02-27 Realization method of Gauss system

Country Status (1)

Country Link
CN (1) CN102591204A (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08137673A (en) * 1994-11-15 1996-05-31 Oki Electric Ind Co Ltd Method for generating random number
CN101782940A (en) * 2008-10-17 2010-07-21 霍夫曼-拉罗奇有限公司 Determination of melting temperatures by equation-less methods
CN202043084U (en) * 2010-12-30 2011-11-16 哈尔滨理工大学 Gauss filter implementing circuit

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08137673A (en) * 1994-11-15 1996-05-31 Oki Electric Ind Co Ltd Method for generating random number
CN101782940A (en) * 2008-10-17 2010-07-21 霍夫曼-拉罗奇有限公司 Determination of melting temperatures by equation-less methods
CN202043084U (en) * 2010-12-30 2011-11-16 哈尔滨理工大学 Gauss filter implementing circuit

Similar Documents

Publication Publication Date Title
CN101476906B (en) Anti-strong fixed interference digital signal processing system of vortex street flowmeter
Upadhyay Class of Recursive Wideband Digital Differentiators and Integrators.
CN102495280B (en) Anti-noise wide-range frequency measurement method and phase locking frequency meter
Vashi et al. Laplace and Sumudu transforms and their application
CN102591204A (en) Realization method of Gauss system
CN103954351B (en) The measuring method of the micro-angular oscillation of a kind of space vehicle
CN103532530A (en) Pulse peak detection device
CN101982782A (en) Chaotic frequency measuring system and method in similar micro-photon source navigation system
CN201773147U (en) Circuit for calculating virtual value of commercial electricity voltage
CN104504272A (en) Digital integrator implementing method based on ideal reconstruction filter impulse response fitting
Kim The solution of Laguerre’s equation by using G-transform
Ahmed Design Analysis of High Pass FIR Filters Using Hanning, Bartlett and Kaiser Windows
CN102710237B (en) First-order digital low-pass filtering method, filter and electronic equipment
US7046053B2 (en) Total harmonic distortion standard
Zhang et al. Signal measurement of projectile penetration overload based on charge sensor
KR100898567B1 (en) Method for calculating integrated value of pulse input signal and Integrator for the method
CN104914136A (en) Differential signal interweaving control soil water content sensor
Cui et al. Flow velocity measurement by cross-correlation with tailored modulation
Banuelos-Saucedo et al. Fast response lock-in amplifier
Adams et al. Approximate realization of fractional-order 2-D IIR frequency-planar filters
Hasler et al. Reconfigurable Analog PDE computation for Baseband and RF Computation
CN204679442U (en) The soil moisture content sensor that differential signal intertexture controls
Vandenbussche et al. Multiplicative finite impulse response filters: implementations and applications using field programmable gate arrays
CN105680826A (en) Method for estimating parameters of analog filter
Xu et al. Distributed arithmetic FIR filter for electrical resistance tomography system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
ASS Succession or assignment of patent right

Owner name: JIANGSU HAITAI PRECISE MEASUREMENT AND CONTROL TEC

Free format text: FORMER OWNER: HARBIN HIGHTECH PRECISION GAGE CO., LTD.

Effective date: 20150203

C41 Transfer of patent application or patent right or utility model
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Yuan Yibao

Inventor after: Xu Lianhu

Inventor after: Pu Weiying

Inventor before: Yuan Yibao

Inventor before: Xu Lianhu

Inventor before: Pu Weiying

Inventor before: Wang Lei

COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 150060 HARBIN, HEILONGJIANG PROVINCE TO: 226010 NANTONG, JIANGSU PROVINCE

Free format text: CORRECT: INVENTOR; FROM: YUAN YIBAO XU LIANHU PIAO WEIYING WANG LEI TO: YUAN YIBAO XU LIANHU PIAO WEIYING

TA01 Transfer of patent application right

Effective date of registration: 20150203

Address after: 226010 workshop No. 11, bureau of Optics and electricity, Nantong Development Zone, Jiangsu, China

Applicant after: Jiangsu Haitai Precision Control Technology Co.

Address before: 150060 Heilongjiang city of Harbin province Harbin Development Zone haping Road District No. 18 North Road in Yantai

Applicant before: Harbin Hightech Precision Gage Co., Ltd.

C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120718