CN103684348A - Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter - Google Patents

Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter Download PDF

Info

Publication number
CN103684348A
CN103684348A CN201310753123.8A CN201310753123A CN103684348A CN 103684348 A CN103684348 A CN 103684348A CN 201310753123 A CN201310753123 A CN 201310753123A CN 103684348 A CN103684348 A CN 103684348A
Authority
CN
China
Prior art keywords
read
filter
low pass
pass filter
order iir
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
CN201310753123.8A
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.)
National University of Defense Technology
Original Assignee
National University of Defense 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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201310753123.8A priority Critical patent/CN103684348A/en
Publication of CN103684348A publication Critical patent/CN103684348A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Complex Calculations (AREA)

Abstract

The invention relates to a multiplication removal rapid algorithm on the basis of a second-order IIR (Infinite Impulse Response) low pass filter, which comprises the following steps: (1) setting filter coefficients, wherein in a mode (1), b0 is equal to 1, b1 is equal to 21, b2 is equal to 1 and a1 is equal to 0, in a mode (2), b0 is equal to 1, b1 is equal to 21, b2 is equal to 1 and a1 is equal to minus one, N1 is a nonnegative integer and N 2 is an integer which is greater than 1; (2) setting an output updating equation of the second-order IIR low pass filter; (3) and in the filtering calculation process, firstly, carrying out shift operation, then using a result for the additive operation and as coefficients of the rest of feedforward quantities or feedback quantities are equal to 1, directly reading the coefficients to carry out additive operation. The multiplication removal rapid algorithm has the advantages that the principle is simple; the total calculated amount can be greatly reduced; rapid operation is implemented; and the like.

Description

Based on second order IIR low pass filter, remove to take advantage of fast algorithm
Technical field
The present invention is mainly concerned with digital filtering technique field, refers in particular to a kind ofly for second order IIR wave digital lowpass filter, to go to take advantage of fast algorithm and filter coefficient method for designing thereof.
Background technology
Digital filter is being brought into play very important effect in the various application of Digital Signal Processing.It is processed to reach the object of frequency filtering by sampled data signal is performed mathematical calculations, and has that precision is high, flexibility is large, reliability is high, be easy to large-scale integrated, can realize the advantages such as parallel processing.Digital filter is to extract that useful information is extremely important, method very flexibly, is the important content of modern signal processing.
Digital filter is divided into two large classes: infinite impulse response (IIR) and finite impulse response (FIR) filter.
From performance, iir filter can obtain higher frequency selectivity with lower exponent number, and memory cell used is few, and economy and efficiency are high; FIR filter can only just can reach high selectivity with higher exponent number, makes cost higher, and higher-order number also causes signal delay larger.Yet, FIR filter can obtain strict linear phase, and the phase-frequency response of iir filter is nonlinear, if will meet linear phase requirement, iir filter carries out phase place calibration with regard to adding all-pass network, and this will increase joint number and the complexity of filter greatly.
From structure, iir filter adopts recursive structure, and the position of limit on z-plane is indefinite, and when limit is not in unit circle, system is by unstable; In addition, the round off processing of meeting to sequence in calculating process, this finite word length effect may cause limit cycles oscillations in recursive structure.And FIR filter adopts onrecurrent structure, no matter existence and stability problem not in theory or in the limited precision operations of reality, arithmetic eror is also less.In addition, FIR filter can utilize fast fourier transform algorithm to realize quick computing, and under identical exponent number, FIR filter arithmetic speed can be more faster than iir filter.
As mentioned above, each tool advantages and disadvantages of dissimilar digital filter; So, in actual applications, need, according to application demand and filter feature, select suitable filter type.
Certain system for real-time signal processing needs a plurality of digital filters to realize the low-pass filtering function of different passbands, requires it to take few physical resource of trying one's best, shorten signal delay as far as possible, and reduces amount of calculation and realize quick computing.Application demand does not claim to linear phase, considers again and reduces physical resource occupancy and shorten signal delay, can realize above-mentioned filter function with low order iir filter; But, also need design for the fast algorithm of low order iir filter, meet the requirement of quick computing.
Second order iir filter system is expressed as follows with difference equation:
y(n)＝-a 1y(n-1)-a 2y(n-2)+b 0x(n)+b 1x(n-1)+b 2x(n-2)
Above formula is also the recursion equation that in filtering, filtering output variable is upgraded, and each filtering output of second order iir filter is upgraded computational process and mainly comprised 5 multiplyings and 4 sub-addition computings; Therefore, the efficient realization of multiply-add operation, is the key point that realizes the quick computing of iir filter.In addition, add operation relatively, multiplying need take more physical resource, will expend the more work period, so the amount of calculation of IIR filtering depends primarily on the number of times of multiplying in filtering algorithm.So in minimizing filtering calculating, the number of times of multiplying, even removes multiplying, can greatly reduce amount of calculation, realizes quick computing.This is also for iir filter Fast Algorithm Design provides a good thinking.
The filter coefficient that there is no at present open method design meets fast algorithm requirement, can, by the multiplying in removal or equivalent substitute filtering computational process, realize the quick computing of low order IIR low-pass filtering.
Summary of the invention
The technical problem to be solved in the present invention is just: the technical problem existing for prior art, the invention provides a kind of principle simple, can greatly reduce overall amount of calculation, realize quick computing based on second order IIR low pass filter, remove to take advantage of fast algorithm.
For solving the problems of the technologies described above, the present invention by the following technical solutions:
Based on second order IIR low pass filter remove to take advantage of a fast algorithm, step is:
(1) filter coefficient is set: during mode (), be: b 0=1, b 1=2 1, b 2=1, a 1=0, during mode (two), be: b 0=1, b 1=2 1, b 2=1, a 1=-1,
Figure BDA0000451172560000022
wherein, N 1for nonnegative integer, N 2for being greater than 1 integer;
(2) set second order IIR low pass filter output renewal equation:
Mode (one): y ( n ) = y ( n - 2 ) - 2 - N 1 y ( n - 2 ) + x ( n ) + 2 x ( n - 1 ) + x ( n - 2 )
Mode (two): y ( n ) = y ( n - 1 ) - 2 - N 2 y ( n - 2 ) + x ( n ) + 2 x ( n - 1 ) + x ( n - 2 )
Wherein, feedforward amount x (n-1) coefficient is 2, and feedback quantity y (n-2) coefficient is 2 -N(wherein N equals N 1or N 2), be all 2 integral number power; Based on binary system, x (n-1) is moved to left 1 and obtains 2x (n-1), by y (n-1) move to right N for position obtain 2 -Ny (n-2);
(3) in filtering computational process, read x (n-1), y (n-2) first does shift operation, x (n-1) is moved to left 1 and obtains 2x (n-1), y (n-1) N that moves to right is obtained to 2 for position -Ny (n-2), then result is used for to add operation; And the coefficient of all the other feedforward amounts or feedback quantity is all 1, therefore directly read them, do add operation.
As a further improvement on the present invention: when mode (one), idiographic flow is:
(1.1) after starting, parameter N assignment;
(1.2) variable k is set, initialize k=3;
Read: x (1), x (2), assignment: y (1)=x (1), y (2)=x (2);
(1.3) read y (k-2), assignment: y (k)=y (k-2);
Displacement: y (k-2) the N position that moves to right obtains 2 -Ny (k-2);
Cumulative: y (k)=y (k)-2 -Ny (k-2);
(1.4) read x (k), cumulative: y (k)=y (k)+x (k);
Read x (k-1), displacement: x (k-1) moves to left 1 and obtains 2x (k-1);
Cumulative: y (k)=y (k)+2x (k-1);
Read: x (k-2), cumulative: y (k)=y (k)+x (k-2);
(1.5) k=k+1, judges whether to finish, and returns to (1.3) circulate as do not finished.
As a further improvement on the present invention: when mode (two), idiographic flow is:
(2.1) after starting, parameter N assignment;
(2.2) variable k is set, initialize k=3;
Read: x (1), x (2), assignment: y (1)=x (1), y (2)=x (2);
(2.3) read y (k-1), assignment: y (k)=y (k-1);
Read y (k-2), displacement: y (k-2) the N position that moves to right obtains 2 -Ny (k-2);
Cumulative: y (k)=y (k)-2 -Ny (k-2);
(2.4) read x (k), cumulative: y (k)=y (k)+x (k);
Read x (k-1), displacement: x (k-1) moves to left 1 and obtains 2x (k-1);
Cumulative: y (k)=y (k)+2x (k-1);
Read: x (k-2), cumulative: y (k)=y (k)+x (k-2);
(2.5) k=k+1, judges whether to finish, and returns to (2.3) circulate as do not finished.
Compared with prior art, the invention has the advantages that: the present invention proposes a kind of configuration mode of filter coefficient, make corresponding second order IIR low pass filter possess special nature, can be used for realizing filter fast algorithm; According to this account form, each filtering output is upgraded computational process and is only included 2 shift operations and 4 sub-addition computings, and the multiplying in filtering calculating is realized by shift operation equivalence.Because the shared physical resource of shift operation and amount of calculation thereof are all far less than multiplying, shift operation is efficient more than multiplying, so can realize quick computing.
Accompanying drawing explanation
Fig. 1 is second order IIR low pass filter filtering performance and the table of comparisons of filter factor when mode (one) after application the inventive method.
Fig. 2 is second order IIR low pass filter filtering performance and the table of comparisons of filter factor when mode (two) after application the inventive method.
Fig. 3 be application during the inventive method in mode () second order IIR low pass filter remove to take advantage of the flow chart of fast algorithm.
Fig. 4 be application during the inventive method in mode (two) second order IIR low pass filter remove to take advantage of the flow chart of fast algorithm.
Embodiment
Below with reference to Figure of description and specific embodiment, the present invention is described in further details.
Because second order IIR low pass filter system transter is expressed as follows:
H ( z ) = b 0 + b 1 z - 1 + b 2 z - 2 1 + a 1 z - 1 + a 2 z - 2
Said system comprises 5 filter coefficient b 0, b 1, b 2, a 1, a 2, make [b 0b 1b 2]=[1 2 1], make system be fixed as z zero point 1,2=-1, can be by changing filter factor a 1, a 2, the pole location of configuration-system, carrys out the filtering performance of selective system.
Configuration mode based on above-mentioned filter coefficient, then according to filter factor a 1, a 2numerical characteristics, for second order IIR low pass filter, the present invention just configures its filter coefficient by mode as shown in following table l.
Table 1:
Figure BDA0000451172560000042
Wherein, N 1for nonnegative integer, N 2for being greater than the integer of l.By changing N 1or N 2value, configurable system pole location, and then choose corresponding filtering performance parameter.Performance parameter mainly refers to the cut-off frequency of second order IIR low pass filter, also comprises gain, adjusting time of filter etc.As scheme as shown in l and Fig. 2 the performance parameter of filter and N 1or N 2value corresponding.In addition the filter system that, in table, listed filter factor is corresponding is all progressive stable.The filter factor of configuration in a manner described, its feature is: the numerical value of filter coefficient can be expressed as 2 integral number power or 2 integral number powers and poor.This feature is to realize the key of removing multiplication algorithm.
Therefore,, in order to realize the quick computing of filter, core of the present invention is filtering output renewal equation.The following is second order IIR low pass filter output renewal equation proposed by the invention.
Mode (one): y ( n ) = y ( n - 2 ) - 2 - N 1 y ( n - 2 ) + x ( n ) + 2 x ( n - 1 ) + x ( n - 2 )
Mode (two): y ( n ) = y ( n - 1 ) - 2 - N 2 y ( n - 2 ) + x ( n ) + 2 x ( n - 1 ) + x ( n - 2 )
Wherein, feedforward amount x (n-1) coefficient is 2, and feedback quantity y (n-2) coefficient is 2 -N(wherein N equals N 1or N 2), be all 2 integral number power, so can replace with " shift operation " " multiplying " of coefficient; That is: based on binary system, x (n-1) the l position that moves to left is obtained to 2x (n-1), y (n-1) IV that moves to right is obtained to 2 for position -Ny (n-2).
So, in filtering computational process, read x (n-1), y (n-2) first does shift operation, then by result for add operation; And the coefficient of all the other feedforward amounts or feedback quantity is all 1, therefore directly read them, do add operation.According to this account form, each filtering output is upgraded computational process and is only included 2 shift operations and 4 sub-addition computings, and the multiplying in filtering calculating is realized by shift operation equivalence.Because the shared physical resource of shift operation and amount of calculation thereof are all far less than multiplying, shift operation is efficient more than multiplying, so can realize quick computing.
According to above-mentioned recursion equation, obtain the flow chart that the present invention removes to take advantage of fast algorithm, respectively as shown in Figure 3 and Figure 4.
Mode (one): after (1.1) start, parameter N assignment; (1.2) variable k is set, initialize k=3; Read: x (1), x (2), assignment: y (1)=x (1), y (2)=x (2); (1.3) read y (k-2), assignment: y (k)=y (k-2); Displacement: y (k-2) the N position that moves to right obtains 2 -Ny (k-2); Cumulative: y (k)=y (k)-2 -Ny (k-2); (1.4) read x (k), cumulative: y (k)=y (k)+x (k); Read x (k-1), displacement: x (k-1) moves to left 1 and obtains 2x (k-1); Cumulative: y (k)=y (k)+2x (k-1); Read: x (k-2), cumulative: y (k)=y (k)+x (k-2); (1.5) k=k+1, judges whether to finish, and returns to (1.3) circulate as do not finished.
Mode (two): after (2.1) start, parameter N assignment; (2.2) variable k is set, initialize k=3; Read: x (1), x (2), assignment: y (1)=x (1), y (2)=x (2); (2.3) read y (k-1), assignment: y (k)=y (k-1); Read y (k-2), displacement: y (k-2) the N position that moves to right obtains 2 -Ny (k-2); Cumulative: y (k)=y (k)-2 -Ny (k-2); (2.4) read x (k), cumulative: y (k)=y (k)+x (k); Read x (k-1), displacement: x (k-1) moves to left 1 and obtains 2x (k-1); Cumulative: y (k)=y (k)+2x (k-1); Read: x (k-2), cumulative: y (k)=y (k)+x (k-2); (2.5) k=k+1, judges whether to finish, and returns to (2.3) circulate as do not finished.
According to above-mentioned algorithm flow chart, the implementation platform that can provide based on application system, programming realizes second order IIR low pass filter and removes to take advantage of fast algorithm.
The present invention is directed to second order IIR low pass filter, according to the numerical characteristics of filter factor, select specific filter coefficient, particularity based on coefficient, utilize " shift operation " equivalence to replace " multiplying " in filtering computational process, will greatly reduce overall operand, improve operation efficiency; On the other hand, sum up the filtering performance of second order IIR low pass filter and the relation of its filter factor, show that form is for designing filter coefficient, make coefficient meet above-mentioned fast algorithm requirement.
Because second order iir filter ssystem transfer function is:
H ( z ) = Y ( z ) X ( z ) = b 0 + b 1 z - 1 + b 2 z - 2 1 + a 1 z - 1 + a 2 z - 2
Under unit step response, system output final value:
y ( ∞ ) = lim z → 1 ( z - 1 ) Y ( z ) = lim z → 1 [ ( z - 1 ) b 0 + b 1 z - 1 + b 2 z - 2 1 + a 1 z - 1 + a 2 z - 2 z z - 1 ] = b 0 + b 1 + b 2 1 + a 1 + a 2
Second order IIR low-pass filter coefficients substitution above formula by the invention described above, can obtain respective filter system gain as follows:
Mode (one): A = y ( ∞ ) = b 0 + b 1 + b 2 1 + a 1 + a 2 = 2 N 1 + 2
Mode (two): A = y ( ∞ ) = b 0 + b 1 + b 2 1 + a 1 + a 2 = 2 N 2 + 2
As can be seen here, above-mentioned second order IIR low pass filter, the feature of its system gain is: yield value can be expressed as 2 integral number power, this feature provides foundation for utilizing shift operation to realize system gain size adjustment.
In conjunction with as depicted in figs. 1 and 2, for going to take advantage of the table of comparisons of fast algorithm second order IIR low pass filter filtering performance and filter factor, the systematic sampling frequency f of classifying as s *during=100kHz, each organizes the cut-off frequency f of the corresponding second order IIR of filter coefficient low pass filter c *, the adjusting time
Figure BDA0000451172560000062
and gain A *.
In different sample frequencys, (establish it for f s) under, the cut-off frequency of establishing filter is f c, the adjusting time is t s, gain is for A, has following relation:
f c f s = f c * f s * , t s * f s = t s * * f s * , A = A *
Above-mentioned relation formula and the table of comparisons are the foundations that design second order IIR low pass filter removes to take advantage of fast algorithm filter coefficient.Based on above-mentioned relation, and by consulting the table of comparisons, can select suitable filter coefficient according to performance of filter parameter request.So, carry out filter coefficient design, be to determine performance of filter parameter according to filtering performance requirement, then according to the process of performance parameter selective filter coefficient.
Designing filter is applied to real system, need know in advance the sample frequency of application system, and according to using Location of requirement design object parameter, then design on this basis corresponding filter.Existing oneself clear and definite design object is second order IIR low pass filter, with the cut-off frequency as a token of parameter of filtering performance, on this basis designing filter coefficient.
The concrete application example of take is example, and oneself knows system proportion f s d, desired design cut-off frequency is f c dsecond order IIR low pass filter, be met the filter coefficient of taking advantage of fast algorithm requirement, its step is as follows:
Step1: calculation expectation equivalent limit frequency:
Figure BDA0000451172560000064
Step2: comparison f tables look-up c *with f c * d, select and f c * dthe f of approximate size c *as the actual equivalent limit frequency of filter;
Step3: obtain filter coefficient: in the table of comparisons, f c *the filter coefficient in being expert at be required coefficient;
Step4: the actual cut-off frequency of calculating filter:
Figure BDA0000451172560000065
Step5: calculate other filtering performance parameter: the adjusting time
Figure BDA0000451172560000066
gain A=A *
The second order IIR low pass filter being designed by said method, the cut-off frequency of its actual cut-off frequency and expectation is not quite identical, reason is: for making filter factor have described feature, to meet the requirement of taking advantage of fast algorithm, in filter coefficient design procedure Step2, equivalent limit frequency has adopted approximation.But the actual cut-off frequency of this filter, near desired value, can meet filtering performance requirement substantially, corresponding filter system is progressive stable; In addition, above-mentioned by the filter coefficient of method design, meet fast algorithm requirement.Therefore, the filter designing as stated above, being applicable to does not have be strict with but need filter to realize the application system of quick computing to design of filter index.
The digital filter of realizing based on the inventive method may be based on software implementation platform or hardware implementation platform.Software implementation platform refers to as realizing specific filter function develops respective filter software, necessary software support environment; Hardware implementation platform refers to as realizing the necessary hardware circuit of specific filter and chip (as FPGA, dsp processor etc.).Based on above-mentioned two class platforms, all can realize the second order IIR low pass filter that the present invention proposes; Its specific implementation, the software and hardware condition providing according to real application systems is determined.But, no matter based on which kind of implementation platform, select which kind of filter implementation, its key and common ground are the algorithm of filter.
Below be only the preferred embodiment of the present invention, protection scope of the present invention is also not only confined to above-described embodiment, and all technical schemes belonging under thinking of the present invention all belong to protection scope of the present invention.It should be pointed out that for those skilled in the art, some improvements and modifications without departing from the principles of the present invention, should be considered as protection scope of the present invention.

Claims (3)

  1. Based on second order IIR low pass filter remove to take advantage of a fast algorithm, it is characterized in that,
    (1) filter coefficient is set; During mode (), be: b 0=1, b 1=2 1, b 2=1, a 1=0, during mode (two), be: b 0=1, b 1=2 1, b 2=1, a 1=-1,
    Figure FDA0000451172550000012
    wherein, N 1for nonnegative integer, N 2for being greater than 1 integer;
    (2) set second order IIR low pass filter output renewal equation:
    Mode (one): y ( n ) = y ( n - 2 ) - 2 - N 1 y ( n - 2 ) + x ( n ) + 2 x ( n - 1 ) + x ( n - 2 )
    Mode (two): y ( n ) = y ( n - 1 ) - 2 - N 2 y ( n - 2 ) + x ( n ) + 2 x ( n - 1 ) + x ( n - 2 )
    Wherein, feedforward amount x (n-1) coefficient is 2, and feedback quantity y (n-2) coefficient is 2 -N, wherein N equals N 1or N 2, be all 2 integral number power; Based on binary system, x (n-1) is moved to left 1 and obtains 2x (n-1), by y (n-1) move to right N for position obtain 2 -Ny (n-2);
    (3) in filtering computational process, read x (n-1), y (n-2) first does shift operation, x (n-1) is moved to left 1 and obtains 2x (n-1), y (n-1) N that moves to right is obtained to 2 for position -Ny (n-2), then result is used for to add operation; And the coefficient of all the other feedforward amounts or feedback quantity is all 1, therefore directly read them, do add operation.
  2. According to claim 1 based on second order IIR low pass filter remove to take advantage of fast algorithm, it is characterized in that, when mode (one), idiographic flow is:
    (1.1) after starting, parameter N assignment;
    (1.2) variable k is set, initialize k=3;
    Read: x (1), x (2), assignment: y (1)=x (1), y (2)=x (2);
    (1.3) read y (k-2), assignment: y (k)=y (k-2);
    Displacement: y (k-2) the N position that moves to right obtains 2 -Ny (k-2);
    Cumulative: y (k)=y (k)-2 -Ny (k-2);
    (1.4) read x (k), cumulative: y (k)=y (k)+x (k);
    Read x (k-1), displacement: x (k-1) moves to left 1 and obtains 2x (k-1);
    Cumulative: y (k)=y (k)+2x (k-1);
    Read: x (k-2), cumulative: y (k)=y (k)+x (k-2);
    (1.5) k=k+1, judges whether to finish, and returns to (1.3) circulate as do not finished.
  3. According to claim 1 based on second order IIR low pass filter remove to take advantage of fast algorithm, it is characterized in that, when mode (two), idiographic flow is:
    (2.1) after starting, parameter N assignment;
    (2.2) variable k is set, initialize k=3;
    Read: x (1), x (2), assignment: y (1)=x (1), y (2)=x (2);
    (2.3) read y (k-1), assignment: y (k)=y (k-1);
    Read y (k-2), displacement: y (k-2) the N position that moves to right obtains 2 -Ny (k-2);
    Cumulative: y (k)=y (k)-2 -Ny (k-2);
    (2.4) read x (k), cumulative: y (k)=y (k)+x (k);
    Read x (k-1), displacement: x (k-1) moves to left 1 and obtains 2x (k-1);
    Cumulative: y (k)=y (k)+2x (k-1);
    Read: x (k-2), cumulative: y (k)=y (k)+x (k-2);
    (2.5) k=k+1, judges whether to finish, and returns to (2.3) circulate as do not finished.
CN201310753123.8A 2013-12-31 2013-12-31 Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter Pending CN103684348A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310753123.8A CN103684348A (en) 2013-12-31 2013-12-31 Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310753123.8A CN103684348A (en) 2013-12-31 2013-12-31 Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter

Publications (1)

Publication Number Publication Date
CN103684348A true CN103684348A (en) 2014-03-26

Family

ID=50320841

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310753123.8A Pending CN103684348A (en) 2013-12-31 2013-12-31 Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter

Country Status (1)

Country Link
CN (1) CN103684348A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107766845A (en) * 2017-11-20 2018-03-06 苏州蓝珀医疗科技股份有限公司 A kind of breathing and BCG method for extracting signal based on light shock sensor

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101867355A (en) * 2009-04-17 2010-10-20 哈曼国际工业有限公司 Use the active noise control system of infinite impulse response filter
EP2315353A1 (en) * 2008-06-10 2011-04-27 Japan Science and Technology Agency Filter
CN103378821A (en) * 2012-04-12 2013-10-30 西门子公司 Filter system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2315353A1 (en) * 2008-06-10 2011-04-27 Japan Science and Technology Agency Filter
CN101867355A (en) * 2009-04-17 2010-10-20 哈曼国际工业有限公司 Use the active noise control system of infinite impulse response filter
CN103378821A (en) * 2012-04-12 2013-10-30 西门子公司 Filter system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈玙等: "一种IIR数字滤波器无乘法结构的实现算法", 《声学技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107766845A (en) * 2017-11-20 2018-03-06 苏州蓝珀医疗科技股份有限公司 A kind of breathing and BCG method for extracting signal based on light shock sensor

Similar Documents

Publication Publication Date Title
CN107636965B (en) Sparse cascaded integration comb filter
Shahein et al. A novel hybrid monotonic local search algorithm for FIR filter coefficients optimization
Lesnikov et al. A new paradigm in design of IIR digital filters
Soltan et al. Butterworth passive filter in the fractional-order
Choo et al. Complexity reduction of digital filters using shift inclusive differential coefficients
CN210142997U (en) Cascaded integrator-comb decimation filter
CN102891662B (en) A kind of general rate down-conversion, up conversion device and method
CN103684348A (en) Multiplication removal rapid algorithm on basis of second-order IIR (Infinite Impulse Response) low pass filter
Singh et al. Digital FIR Filter Designs
Abramovitch The multinotch, part II: Extra precision via Δ coefficients
CN108257628B (en) Lattice scale fractional order memristor for arbitrary order high-pass filtering
Owen et al. 384 TMAC/s FIR filtering on an Artix-7 FPGA using Prism signal processing
Sergiyenko et al. VHDL generation of optimized IIR filters
Poučki et al. Elliptic IIR filter sharpening implemented on FPGA
Karimov et al. Compact Fixed-Point Filter Implementation
Damian et al. A low area FIR filter for FPGA implementation
Mekhnache et al. Using new hybrid integrators to design IIR fractional-order integrators via impulse invariance method
Dolecek et al. Low power non-recursive comb-based decimation filter design
CN102868385B (en) Digital filter capable of configuring coefficient and implementation method
Chan et al. A new method for least-squares and minimax group-delay error design of allpass variable fractional-delay digital filters
CN116579284B (en) Method, system and related equipment for designing frequency band of surface acoustic wave filter
CN102710237A (en) First-order digital low-pass filtering method, filter and electronic equipment
CN103095251A (en) Improving method for high speed of digital filter
Mirković et al. Design of IIR digital filters with critical monotonic passband amplitude characteristic-A case study
Dam et al. Allpass VFD filter design

Legal Events

Date Code Title Description
PB01 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

Application publication date: 20140326

RJ01 Rejection of invention patent application after publication