CN105337588A - Combined integral filtering system - Google Patents

Combined integral filtering system Download PDF

Info

Publication number
CN105337588A
CN105337588A CN201510854733.6A CN201510854733A CN105337588A CN 105337588 A CN105337588 A CN 105337588A CN 201510854733 A CN201510854733 A CN 201510854733A CN 105337588 A CN105337588 A CN 105337588A
Authority
CN
China
Prior art keywords
combined integral
tau
link
signal
time delay
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
CN201510854733.6A
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.)
Donghua University
National Dong Hwa University
Original Assignee
Donghua University
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 Donghua University filed Critical Donghua University
Priority to CN201510854733.6A priority Critical patent/CN105337588A/en
Publication of CN105337588A publication Critical patent/CN105337588A/en
Pending legal-status Critical Current

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

The invention relates to a combined integral filtering system. The combined integral filtering system is characterized in that a first combined integral link and a time delay link are included, and signals from the first combined integral link and the time delay link are superposed to a second combined integral link; the overall transfer function of the combined integral filtering system is shown in the specification, the time delay constants tau1, tau2, tau3 and the constant K are adjusted according to requirements, and low-pass, high-pass and band-pass filtering is achieved. The filtered-out signals are stable, the signals with only one frequency are contained, noise is low, and the signals are smooth.

Description

A kind of combined integral filtering system
Technical field
The present invention relates to industrial control field, particularly relate to a kind of combined integral filtering system.
Background technology
Environment residing for Industry Control object is comparatively severe, and interference source is a lot, as strong electromagnetic, ambient temperature frequent variations, especially signal input channel, sampled value will be used for data processing and system fading margin, and the accuracy of sampled value is significant for the raising of robot control system(RCS) performance.Therefore in order to reduce the impact of interference on sampled value, the performance of guarantee system, first will carry out signal filtering to sampled value before data processing of being everlasting, to improve authenticity, the accuracy of signal, eliminate or reduce various interference and noise, improving the reliability of computer system further.
Filtering technique is often divided into two classes: analog filtering and digital filtering technique.Digital filtering technique overcomes the shortcoming of a lot of analog filtering, and mainly realize with programmed algorithm, reliability is high, good stability, can take different filtering methods according to signal type difference, very convenient.The goal in research of digital filtering algorithm can use up the interference removing other possible greatly after the signal digital process that can make to collect after digital filtering algorithm process, reverts to its actual value.
Enter into digital times now, people extremely pay close attention to Digital Signal Processing.Along with computer and microelectric technique develop, signal processing theory and algorithm obtain primary study, and be widely used in the field such as voice image processing, digital communication, analysis of spectrum, pattern recognition, automatically control, filter acts on significantly in digital signal processing.Digital filter be input as a Serial No., after mathematical operation rule, export as another Serial No..Change the relative scale of frequency content contained by input signal by operation relation or eliminate the signal of some frequency, by useful signal constituents extraction out, realizing the object weakening interference component.
Digital filter obtains general application, digital filter characteristic difference, hinders its application in other respects.As iir digital filter still has the excellent characteristic of analog filter, this sacrifices linear phase for cost.Namely utilize the methods such as Butterworth, Chebyshev and ellipse to complete digital filter design, make its width product characteristic close to ideal filter, the filter characteristic of acquisition is nonlinear often.Finite Impulse Response filter can use window function, frequency sampling, weighted least squares at any amplitude of design and frequency characteristic filter.
Classical digital filter is not enough:
(1) given technical indicator, select a kind of in multi-design methods of comforming, classical technology implies a lot of filter design method;
(2) iir filter transfer function limit is at unit circle Arbitrary distribution, and available lower-order number obtains multiple choices, and memory cell is few, high efficiency.Transfer function limit may be in outside unit circle, causes system unstable.Iir digital filter phase place is nonlinear, the non-linear decision selectivity quality of phase place;
(3) the transfer function limit of FIR filter is in initial point always, and selectivity is more, has come with high exponent number, and under design of filter index same case, the exponent number required by FIR filter than iir filter height 5-10 doubly, spends high, postpones a lot.
Summary of the invention
Technical problem to be solved by this invention is to provide a kind of combined integral filtering system, can improve filter effect.
The technical solution adopted for the present invention to solve the technical problems is: provide a kind of combined integral filtering system, comprise the first combined integral link and time delay process, described first combined integral link and time delay process Signal averaging to the second combined integral link out; The overall transfer function of described combined integral filtering system is wherein, for pure lag system, τ itime delay constant; Regulate time delay constant τ according to demand 1, τ 2, τ 3and constant K realizes low pass, high pass and bandpass filtering.
Described first combined integral link is low pass filter, and its transfer function is
The high pass filter that described first combined integral link and time delay process are formed, its transfer function is g h ( s ) = e - τ 2 s - 1 τ 1 s ( 1 - e - τ 1 s ) .
The transfer function of described second combined integral link is
Beneficial effect
Owing to have employed above-mentioned technical scheme, the present invention compared with prior art, has following advantage and good effect:
The present invention adopts combined integral and superposes with Time Delay, determines Selecting parameter principle, designs low pass, high pass, band pass filter respectively, realizes industrial collection signal filtering.Based on above reason, the special system performance analysis method for combined integral process of research and applying at other field of combined integral link have important theory and realistic meaning.
The present invention adopts the signal gone out based on the advanced filtering algorithm filtering of combined integral system very stable, and the signal only containing a kind of frequency, noise is few, signal smoothing.The Butterworth filter of contrast iir digital filter design, in filtered signal, noise is more, adopts combined integral system filter effect better.
Accompanying drawing explanation
Fig. 1 is the advanced filtering algorithm block diagram of combined integral system of the present invention.
Fig. 2 is the concrete block diagram of combined integral filtering system of the present invention.
Embodiment
Below in conjunction with specific embodiment, set forth the present invention further.Should be understood that these embodiments are only not used in for illustration of the present invention to limit the scope of the invention.In addition should be understood that those skilled in the art can make various changes or modifications the present invention, and these equivalent form of values fall within the application's appended claims limited range equally after the content of having read the present invention's instruction.
Embodiments of the present invention relate to a kind of combined integral filtering system, adopt combined integral and superpose with Time Delay, determine Selecting parameter principle, design low pass, high pass, band pass filter respectively, realize industrial collection signal filtering.
Principle of the present invention is:
(1) a kind of new filtering algorithm put forward on the basis of digital filtering algorithm, by constructing system model, regulates time delay constant τ according to demand 1, τ 2, τ 3and constant K realizes filtering, is easily understood, easy to operate.
(2) define: combined integral is to liking a class open loop stabilization, and its transfer function is made up of two or more integration Object with Time Delay and can expresses with following form:
Σ i = 1 n k i sG i ( s ) ( 1 - e - τ 1 i s ) e - τ 2 i s - - - ( 1 )
G is () is a Stable Polynomials not comprising integral element; formula is included in the systems stabilisation in combined integral system; In order to ensure the continuity of this process, time delay constant τ 1i, τ 2i, must equation be met: τ 2i2 (i-1)+ τ 1 (i-1), 1 < i < n, n is positive integer.Combined integral process comprises non-Left half-plane pole and zero simultaneously.There is anomalous integral minimum phase characteristic simultaneously.Integral process is classified: divide without time-lag process and time-lag process by time lag type; Single combined integral process and many combined integrals process is divided into by system complexity; Low order process and Advanced Process is divided into by systematic education.
Analysis on Mechanism modeling is adopted to obtain the transfer function model of the combined integral process in industrial process:
G p ( s ) = k &tau; s ( 1 - e - &tau; s ) - - - ( 2 )
G p ( s ) = k &tau; 1 s ( 1 - e - &tau; 1 s ) e - &tau; 2 s - - - ( 3 )
G p ( s ) = k &tau; s ( T s + 1 ) ( 1 - e - &tau; s ) - - - ( 4 )
G p ( s ) = k &tau; 1 s ( T s + 1 ) ( 1 - e - &tau; 1 s ) e - &tau; 2 s - - - ( 5 )
G p ( s ) = k &tau; 1 s ( T S + 1 ) ( 1 - e - &tau; 1 s ) e - &tau; 2 s + k &tau; 3 s ( T s + 1 ) ( 1 - e - &tau; 3 s ) e - &tau; 4 s - - - ( 6 )
Suppose Bounded Index signal f (t), its Laplace is transformed to F (s), is multiplied by combined integral link g (s):
g ( s ) = 1 &tau; s ( 1 - e - &tau; s ) - - - ( 7 )
Obtain transfer function:
G ( s ) = F ( s ) 1 &tau; s ( 1 - e - &tau; s ) = 1 &tau; ( F ( s ) s - F ( s ) e - &tau; s s ) - - - ( 8 )
Laplace inverse transformation is asked to obtain to formula (8):
G ( t ) = L - 1 ( G ( s ) ) = 1 &tau; &lsqb; &Integral; 0 t f ( u ) d u - &Integral; 0 t - &tau; f ( u ) d u &rsqb; = &Integral; t - &tau; t f ( u ) d u &tau; - - - ( 9 )
It is that the τ time period (t-τ, t) carries out integration that formula (9) to may be interpreted as signal in length, carries out τ average to gained integration divided by time span, and being equivalent to signal in length is the Mean Filtering Algorithm that the τ time period carries out.Therefore the combined integral link that formula (7) defines is equal to a mean filter, has phase same-action with mean filter.Mean filter is the special case of low pass filter, and releasing this combined integral link is thus a low pass filter.There is the ability of certain restraint speckle, because its essence is equivalent to mean filter, some can be suppressed thus completely with the cyclical signal of specific period.
As shown in Figure 1, this combined integral filtering system, comprises the first combined integral link and time delay process, described first combined integral link and time delay process Signal averaging to the second combined integral link out; The overall transfer function of described combined integral filtering system is:
G ( s ) = K ( e - &tau; 1 s - 1 - e - &tau; 2 s &tau; 1 s ) &CenterDot; 1 - e - &tau; 3 s &tau; 3 s - - - ( 10 )
Wherein, for pure lag system, τ itime delay constant; Regulate time delay constant τ according to demand 1, τ 2, τ 3and constant K realizes low pass, high pass and bandpass filtering.
Described first combined integral link is low pass filter, and its transfer function is
Derived by formula (8) and formula (9) and know, the mean filter in time range that this link completes [t-τ, t].This link realizes the effect of low pass signal filtering, time delay constant τ 1controlled, determine the item number participating in mean filter.τ 1larger, low-pass filter effect is better, and the high fdrequency component of elimination is more, and filtered signal is more level and smooth.
The high pass filter that described first combined integral link and time delay process are formed, its transfer function is g h ( s ) = e - &tau; 2 s - 1 &tau; 1 s ( 1 - e - &tau; 1 s ) .
As seen from Figure 2, former composite signal by after a time delay process with low-pass filtering after Signal averaging, because time delay can not change the characteristic of signal, just signal phase is had an impact.Known by low-pass filtering link, system by and signal passes through after link, only remaining low frequency signal in signal, then superpose with this link, namely deduct low frequency signal with original signal, so realize effect and be equal to High frequency filter effect, so link achieves the effect of high-pass filtering.
The transfer function of described second combined integral link is
G'(s) the equivalent g of effect is equal to 1the effect of (s) link, different unlike the item number participating in low-pass filtering, i.e. parameter τ 1, τ 2vary in size, and also different for the signal superposed.G 1s remaining in composite signal after () link is mixed signal medium and low frequency signal, with g hs () link superposes, remaining is high-frequency signal in signal, again superposes a low frequency filtering link, as long as Selecting parameter is reasonable, can realize the effect of bandpass filtering.
Added a Proportional coefficient K herein, because original signal is by having decay after two filtering links, amplitude is less, after increasing this link, by the amplification K of filtered signal doubly, can reduce the error of signal and filtering signal, with the contrast of original signal, can check than being easier to.
Adopt the signal gone out based on the advanced filtering algorithm filtering of combined integral system very stable, the signal only containing a kind of frequency, noise is few, signal smoothing.The Butterworth filter of contrast iir digital filter design, in filtered signal, noise is more, adopts combined integral system filter effect better.

Claims (4)

1. a combined integral filtering system, is characterized in that, comprises the first combined integral link and time delay process, described first combined integral link and time delay process Signal averaging to the second combined integral link out; The overall transfer function of described combined integral filtering system is G ( s ) = K ( e - &tau; 1 s - 1 - e - &tau; 2 s &tau; 1 s ) &CenterDot; 1 - e - &tau; 3 s &tau; 3 s , Wherein, for pure lag system, τ itime delay constant; Regulate time delay constant τ according to demand 1, τ 2, τ 3and constant K realizes low pass, high pass and bandpass filtering.
2. combined integral filtering system according to claim 1, is characterized in that, described first combined integral link is low pass filter, and its transfer function is g 1 ( s ) = 1 &tau; 1 s ( 1 - e - &tau; 1 s ) .
3. combined integral filtering system according to claim 1, is characterized in that, the high pass filter that described first combined integral link and time delay process are formed, and its transfer function is
4. combined integral filtering system according to claim 1, is characterized in that, the transfer function of described second combined integral link is g &prime; ( s ) = K &tau; 3 s ( 1 - e - &tau; 3 s ) .
CN201510854733.6A 2015-11-30 2015-11-30 Combined integral filtering system Pending CN105337588A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510854733.6A CN105337588A (en) 2015-11-30 2015-11-30 Combined integral filtering system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510854733.6A CN105337588A (en) 2015-11-30 2015-11-30 Combined integral filtering system

Publications (1)

Publication Number Publication Date
CN105337588A true CN105337588A (en) 2016-02-17

Family

ID=55287932

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510854733.6A Pending CN105337588A (en) 2015-11-30 2015-11-30 Combined integral filtering system

Country Status (1)

Country Link
CN (1) CN105337588A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105529022A (en) * 2015-12-18 2016-04-27 东华大学 Self-adaptive combined integration filtering apparatus and noise control system and method
CN107272404A (en) * 2017-07-13 2017-10-20 东华大学 A kind of combined integral controller
CN110221544A (en) * 2019-06-14 2019-09-10 四川烟叶复烤有限责任公司 A kind of tobacco redrying entrance moisture data filtering method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1035215A (en) * 1987-12-12 1989-08-30 德国Itt工业有限公司 Digital decimation filter
WO1993005573A1 (en) * 1991-09-09 1993-03-18 Quinton Instrument Company Ecg baseline wander filter and method
CN102780469A (en) * 2012-08-16 2012-11-14 钜泉光电科技(上海)股份有限公司 Cascade integrator comb filter and implementation method thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1035215A (en) * 1987-12-12 1989-08-30 德国Itt工业有限公司 Digital decimation filter
WO1993005573A1 (en) * 1991-09-09 1993-03-18 Quinton Instrument Company Ecg baseline wander filter and method
CN102780469A (en) * 2012-08-16 2012-11-14 钜泉光电科技(上海)股份有限公司 Cascade integrator comb filter and implementation method thereof

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李娜 等: "基于组合积分系统的滤波", 《电子设计工程》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105529022A (en) * 2015-12-18 2016-04-27 东华大学 Self-adaptive combined integration filtering apparatus and noise control system and method
CN107272404A (en) * 2017-07-13 2017-10-20 东华大学 A kind of combined integral controller
CN110221544A (en) * 2019-06-14 2019-09-10 四川烟叶复烤有限责任公司 A kind of tobacco redrying entrance moisture data filtering method

Similar Documents

Publication Publication Date Title
CN105222809B (en) The method of the boat magnetic disturbance penalty coefficient estimation of a kind of ground magnetic gradient robust
CN103674001B (en) A kind of optical fibre gyro denoising method based on strengthening self-adaptation time-frequency method
CN105913835B (en) Adaptive filtering method based on mathematical morphology
CN104579243B (en) A kind of Fast design method of DFT modulated filters group
CN105337588A (en) Combined integral filtering system
CN101425152B (en) Design method for FIR filter based on learning rate changing neural net
CN107276560B (en) A kind of FRI sparse sampling kernel function construction method and circuit
CN104253650B (en) The estimation unit of intrachannel nonlinear damage and method
CN103308829B (en) A kind of GIS single Partial discharge signal extracts and trigger instants method of adjustment
CN116088320B (en) Sparse double-space non-convex penalty self-adaptive Chirp modal cross aliasing decomposition method
CN106844291A (en) For controlling the method and system with collaborative simulation physical system
CN107122724A (en) A kind of method of the online denoising of sensing data based on adaptive-filtering
Badajos et al. Reduction of audio noise with lowpass Chebyshev Type II filter simulated using GNU octave
CN103152298B (en) Blind signal reconstruction method based on distribution-type compressed sensing system
CN104571401A (en) Implementing device of high-speed guiding filter on FPGA (field programmable gate array) platform
CN103580647B (en) A kind of filter construction
CN104539264A (en) Filtering method and filtering circuit applied to EPS system
CN102624660B (en) Four-weighted fractional Fourier transform-based narrowband interference suppression method
CN104504272B (en) The digital integrator implementation method being fitted based on perfect reconstruction filter impulse responses
CN103777221B (en) Based on the digital core pulse signal Gauss manufacturing process of window function metht
CN102227096A (en) Identification method for variable step-size least mean P-norm system in non-Gaussian environment
CN115240418A (en) Short-term traffic flow prediction method based on causal gating-low-pass graph convolution network
CN104618033B (en) A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method
CN104655136A (en) Multi-concave-point FIR filtering method applicable to laser strap-down inertial navigation system
CN109639258B (en) Design method of two-dimensional FIR notch filter based on Hopfield neural network

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160217