CN101340182B - Low-complexity implementing method and apparatus for FIR digital filter group - Google Patents

Low-complexity implementing method and apparatus for FIR digital filter group Download PDF

Info

Publication number
CN101340182B
CN101340182B CN2008101190648A CN200810119064A CN101340182B CN 101340182 B CN101340182 B CN 101340182B CN 2008101190648 A CN2008101190648 A CN 2008101190648A CN 200810119064 A CN200810119064 A CN 200810119064A CN 101340182 B CN101340182 B CN 101340182B
Authority
CN
China
Prior art keywords
filter
kernel
impulse response
finite impulse
group
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.)
Expired - Fee Related
Application number
CN2008101190648A
Other languages
Chinese (zh)
Other versions
CN101340182A (en
Inventor
彭克武
刘在爽
杨知行
宋健
符剑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN2008101190648A priority Critical patent/CN101340182B/en
Publication of CN101340182A publication Critical patent/CN101340182A/en
Application granted granted Critical
Publication of CN101340182B publication Critical patent/CN101340182B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Complex Calculations (AREA)

Abstract

The invention relates to a low complexity realization method and a device of an FIR digital filter group, the method comprises the following steps: an inner core filter of the filter group is constituted according to the common inner core of a plurality of similar filters in the filter group; the original input vector which is obtained by a tapped delay line is converted according to the corresponding linear convolution relationship between various filters in the group and the inner core filters, thereby constituting a plurality of new input vectors aiming at the inner core filter; the filter group shares a weighted sum unit of the inner core filter, and the respective filter operation of each filter in the group is realized by time-division multiplexing. The realization device of the FIR digital filter group comprises a tapped delay line unit, a vector conversion unit, a gating unit, the weighted sum unit of the inner core filter and a parallel output unit of the inner core filter. The realization method and the device multiplex the inner core filter operational structure which occupies more resources to the maximum extent, thereby reducing the complexity of the hardware realization of the FIR digital filter group and improving the flexibility and the applicability of a system.

Description

The low complex degree implementation method and the device of Finite Impulse Response filter group
Technical field
The invention belongs to digital signal processing technique field, relate to the realization of digital filter, relate in particular to the low complex degree implementation method and the device of Finite Impulse Response filter group.
Background technology
FIR (Finite Impulse Response, finite impulse response) digital filter is a basic module commonly used in digital communication and the signal processing system, occupies important status.Finite Impulse Response filter is widely used in audio-video signal processing and transmission, base band shaping filtering, goes multiple occasions such as mirror image filtering and matched filtering.In the middle of the application of reality, Finite Impulse Response filter often need be realized in various programmable logic devices, digital signal processor spare or application-specific integrated circuit (ASIC).
In recent years, along with expansion and the raising of communication with the signal processing demand, Finite Impulse Response filter more and more is subjected to extensive concern.Regular meeting runs into the situation of having to design higher order filter because of the filtering requirements raising during researcher.For example in China Digital TV terrestrial broadcasting transmission national standard (GB20600-2006), it is 0.05 that the rolloff-factor of the Square-root Raised Cosine FIR filter of modulated terminal shaping filter is required, its transition band is very narrow, the availability of frequency spectrum is very high, this has high requirements to filter order, for hardware is realized, directly caused very high complexity.Therefore the hard-wired optimization of Finite Impulse Response filter is a very real problem.
Finite Impulse Response filter can determine that establishing filter order is finite value N by system function is unique, and its system function Z-transformation can be expressed as:
H ( z ) = Σ n = 0 N h [ n ] · z - n - - - ( 1 )
H[n wherein] be the filter finite impulse response, i.e. the coefficient of Finite Impulse Response filter.Then the filter input/output relation is
y [ n ] = x [ n ] * h [ n ] = Σ m = 0 N x [ n - m ] · h [ m ] , - - - ( 2 )
Wherein * represents linear convolution, x[n] be input signal, y[n] be output signal.
The direct type implementation structure of Finite Impulse Response filter as shown in Figure 1, its part mainly comprises tapped delay line, multiplication unit and sum unit three parts, the main composition of other implementation structure (as the transposition type) is with directly type is similar.Directly the computing flow process of type structure as shown in Figure 2, length is the different delayed time that the tapped delay line of N obtains input signal, delay input signal x[n-m] with FIR filter coefficient h[m] multiply each other, whole multiplied result summations are obtained filtering output.
For the linear phase fir digital filter, according to the intrinsic symmetric properties of its filter coefficient, remake multiplication after can merging the tap of coefficient symmetry (as identical or opposite), its implementation structure is as shown in Figure 3.Make like this quantity of required multiplying almost reduce to original half, saved the resource occupation amount during filter hardware is realized greatly.
Directly merge symmetrical taps and realize it being effective processing method for the hardware of an independent linear phase fir digital filter.But this does not utilize symmetry to reduce the complexity of tapped delay line, comprises and reduces delay line length and tap number.Particularly for the realization up to the Large N FIR digital filter on hundreds of rank, the complexity that reduces tapped delay line is very necessary behave.Utilize self similitude (as symmetry) or according to the demand of practical application, a higher order filter can resolve into the subfilter of two or more lower-orders, forms bank of filters.In the digital communication or signal processing system of reality; Shi Changhui runs into this Large N FIR digital filter and decomposes back composition bank of filters is carried out filtering simultaneously to same input signal situation; such as the interpolation or the decimation filter group that adopt heterogeneous (Polyphase) structure, interpolation multiphase filter group as shown in Figure 4.Be not difficult to find that each subfilter in the group may nonlinear phase filter, may there be symmetry in himself coefficient, and the method that merges symmetrical taps just may lose efficacy.
For realizing the Finite Impulse Response filter group, traditional prioritization scheme can only be shared the tapped delay line of input signal, the delay input signal x[n-m of gained] respectively with a plurality of subfilters coefficient h separately i[m] is weighted summation operation and obtains a plurality of filtering output y i[n].Above-mentioned traditional implementation structure as shown in Figure 5, there are the following problems for it:
1. hard-wired complexity is still higher, and each subfilter works alone fully, does not make cooperate optimization.Need particularly to realize that the demand of its hardware resource is more huge under the situation of Large N FIR digital filter;
2. if the whole filter group need be adjusted or change, then wherein each subfilter all needs to realize again separately and optimize that system flexibility and adaptability are not high.
Summary of the invention
At the particularly big and flexibility of the complexity height, the hardware resource occupancy that exist of the traditional implementation method of higher order filter group in various programmable logic devices, digital signal processor spare or application-specific integrated circuit (ASIC) of Finite Impulse Response filter group, the not high many weak points of adaptability, the purpose of this invention is to provide the low complex degree implementation method and the implement device of Finite Impulse Response filter group.
For reaching above-mentioned purpose, on the one hand, technical scheme of the present invention provides the low complex degree implement device of Finite Impulse Response filter group, and this Finite Impulse Response filter group comprises M N rank Finite Impulse Response filter, and wherein M is more than or equal to 2, and this device comprises:
The clock source unit provides the clock of f and two kinds of frequencies of Mf;
The tapped delay line unit in the frequency that the clock source unit provides is under the clock driving of f, with input signal x[n] delaying time obtains input time delay signal x[n-m], be expressed as
Figure G2008101190648D00031
Input vector;
Vector transduced cell is used for the input vector that the tapped delay line unit is obtained
Figure G2008101190648D00032
Change, by changing
Figure G2008101190648D00033
In element ordering and sign and shift left accordingly, obtain the new input vector of M
Figure G2008101190648D00034
Make
Figure G2008101190648D00035
With
Figure G2008101190648D00036
Corresponding relation and similar filter coefficient vector
Figure G2008101190648D00037
With the kernel filters coefficient vector
Figure G2008101190648D00038
Corresponding relation mate fully, promptly
Figure G2008101190648D00039
1≤i≤M;
Gating unit is under the Mf that the clock source unit provides leads, to M new input vector The timesharing gating;
The kernel filters weighted sum unit under the Mf frequency that the clock source unit provides, is exported from gating unit
Figure G2008101190648D000311
Timesharing is carried out and the kernel filters coefficient vector
Figure G2008101190648D000312
The weighted sum computing, obtain M group filtering operation result;
The kernel filters output unit that walks abreast, under the Mf frequency that the clock source unit provides, the M group filtering operation result that kernel filter weight sum unit is obtained is with 0 to M-1 clock cycle of not waiting of frequency Mf time-delay, and the output M that runs simultaneously organizes the filtering operation result.
In addition, this device also comprises:
The error-filter weighted sum unit is under the f frequency that the clock source unit provides, to input vector
Figure G2008101190648D00041
Carry out and the error-filter coefficient vector The weighted sum computing;
The summation output unit with M group filtering operation result and its each self-corresponding error-filter operation result addition that the parallel output unit of kernel filters is run simultaneously and exported, obtains the filtering operation result and the output of each Finite Impulse Response filter.
On the other hand, technical scheme of the present invention also provides the low complex degree implementation method of Finite Impulse Response filter group, and this Finite Impulse Response filter group comprises M N rank Finite Impulse Response filter, and wherein M may further comprise the steps more than or equal to 2:
The input signal x[n of described Finite Impulse Response filter group] obtain input vector through tapped delay line
Figure G2008101190648D00043
The coefficient vector of each Finite Impulse Response filter is respectively in the group
Figure G2008101190648D00044
1≤i≤M;
The kernel filters constitution step, similar characteristic according to each Finite Impulse Response filter in the group, extract the same section of each Finite Impulse Response filter kernel, construct one with the kernel of the error minimum of each Finite Impulse Response filter kernel kernel as the Finite Impulse Response filter group, and the formation kernel filters, its coefficient vector is
The Finite Impulse Response filter decomposition step, from the set that the similar filter of all kernels similar to above-mentioned kernel filters constitutes, according to the filter error criterion, by search, obtain respectively and the similar filter of kernel of organizing interior each Finite Impulse Response filter coefficient error minimum, constitute coefficient vector and be respectively The similar bank of filters of kernel, 1≤i≤M, and each Finite Impulse Response filter in former group is decomposed into the similar filter of kernel of coefficient error minimum with it, and (coefficient vector is
Figure G2008101190648D00047
) (coefficient vector is with error-filter corresponding to the similar filter of this kernel ) two parts, the filtering operation of former group of interior each Finite Impulse Response filter is y as a result i[n] also correspondingly is decomposed into the similar filter operation result of kernel y Kerneli[n] and corresponding error-filter operation result y Errori[n] two parts;
The vector switch process is to input vector
Figure G2008101190648D00051
Change, by changing
Figure G2008101190648D00052
In element ordering and sign and shift left accordingly, structure is at the M of the similar filter of kernel new input vector
Figure G2008101190648D00053
Make
Figure G2008101190648D00054
With
Figure G2008101190648D00055
Corresponding relation mate fully With
Figure G2008101190648D00057
Corresponding relation, promptly
Figure G2008101190648D00058
1≤i≤M;
The multiplexing step of kernel filters weighted sum, the similar bank of filters of kernel is multiplexing kernel filters weighted sum unit in the different time periods, realizes the filtering operation of the similar filter of each kernel in the similar bank of filters of kernel;
Error-filter weighted sum step is weighted summation to the corresponding error-filter of each Finite Impulse Response filter, realizes the filtering operation of each error-filter
Figure G2008101190648D00059
1≤i≤M;
Summation output step, the error-filter that the output of the similar filter of each kernel in the similar bank of filters of above-mentioned kernel is corresponding with it is exported addition, thereby obtains the output of each Finite Impulse Response filter in the Finite Impulse Response filter group:
y i [ n ] = y kerneli [ n ] + y errori [ n ] = x → · h → kerneli + x → · h → errori = x i ′ → · h → kernel + x → · h → errori , 1≤i≤M。
Wherein, the low complex degree implementation method of above-mentioned Finite Impulse Response filter group, described filter error criterion are meant the mean square error minimum criteria of filter time domain impulse response or the mean square error minimum criteria of filter frequency domain response.
On the one hand, technical scheme of the present invention also provides the low complex degree implementation method of Finite Impulse Response filter group again, and this Finite Impulse Response filter group comprises M N rank FIR similarity number character filter, and wherein M may further comprise the steps more than or equal to 2:
The input signal x[n of described Finite Impulse Response filter group] obtain input vector through tapped delay line
Figure G2008101190648D000511
The coefficient vector of each Finite Impulse Response filter in the Finite Impulse Response filter group is respectively 1≤i≤M;
The kernel filters constitution step, the kernel total according to similar filter constitutes kernel filters, establishes its coefficient vector and is
The vector switch process is to input vector
Figure G2008101190648D000514
Change, by changing In element ordering and sign and shift left accordingly, structure is at the M of kernel filters new input vector Make
Figure G2008101190648D00062
With
Figure G2008101190648D00063
Corresponding relation mate fully
Figure G2008101190648D00064
With
Figure G2008101190648D00065
Corresponding relation, promptly
Figure G2008101190648D00066
Here 1≤i≤M;
The multiplexing step of kernel filters weighted sum, described Finite Impulse Response filter group is multiplexing kernel filters weighted sum unit in the different time periods, realizes the filtering operation of each filter in the described Finite Impulse Response filter group.
The beneficial effect of technical scheme provided by the invention is: the present invention makes full use of the similar characteristic between each Finite Impulse Response filter in the Finite Impulse Response filter group, deeply excavate optimization potentiality wherein, the multiplexing to greatest extent more kernel filtering operation structure of resource that takies is optimized and is realized particularly higher order filter group of Finite Impulse Response filter group.Compare with traditional implementation method, effectively reduce hard-wired complexity, reduced hardware resource and taken, improved the flexibility and the applicability of system, and be convenient to expand, realize the bank of filters of a plurality of Finite Impulse Response filters compositions.
Description of drawings
Fig. 1 is the direct type structural representation of existing Finite Impulse Response filter;
Fig. 2 is the computing schematic flow sheet of the direct type structure of existing Finite Impulse Response filter;
Fig. 3 is the direct type structural representation of existing linear phase fir digital filter (supposition N is an even number);
Fig. 4 is the general structural representation of existing interpolation multiphase filter group;
Fig. 5 is traditional implementation structure schematic diagram of existing Finite Impulse Response filter group;
Fig. 6 is the optimization multiplexing structure schematic diagram of similar bank of filters provided by the invention;
Fig. 7 is the operating process schematic diagram of the Finite Impulse Response filter group low complex degree implementation method of being made up of similar filter provided by the invention;
Fig. 8 is the structural representation (M=2) of the Finite Impulse Response filter group low complex degree implement device of being made up of similar filter provided by the invention;
Fig. 9 is the principle schematic (M=2) of kernel filters weighted sum unit provided by the invention;
Figure 10 is the principle schematic (M=2) of the parallel output unit of kernel filters provided by the invention;
Figure 11 is the structural representation (M=2) of the Finite Impulse Response filter group low complex degree implement device of being made up of the incomplete similarity filter provided by the invention;
Figure 12 is the structural representation (M>2) of the Finite Impulse Response filter group low complex degree implement device of being made up of the incomplete similarity filter provided by the invention.
Embodiment
Following examples are used to illustrate the present invention, but are not used for limiting the scope of the invention.
For making purpose of the present invention, technical method and advantage clearer, embodiment of the present invention is described further in detail below in conjunction with accompanying drawing.
The invention provides the low complex degree implementation method and the implement device of Finite Impulse Response filter group.Realize optimizing at the hardware of Finite Impulse Response filter group.With reference to Fig. 1, Fig. 2, represent direct type structure of Finite Impulse Response filter and computing schematic flow sheet thereof respectively; Fig. 3 represents the direct type structural representation of linear phase fir digital filter (supposition N is an even number), is characterized in that the tap that at first merges symmetry under the prerequisite of coefficient symmetry remakes multiplication; Fig. 4 represents the general structural representation (L is an interpolation factor) of interpolation multiphase filter group; Fig. 5, Fig. 6 represent the traditional implementation structure schematic diagram of bank of filters and the optimization multiplexing structure schematic diagram of similar bank of filters provided by the invention respectively.Below abbreviate Finite Impulse Response filter as filter.
In fact, above-mentioned bank of filters is made up of a plurality of filters, its input signal is identical, and general certain inner link that exists between a plurality of filters, the present invention is referred to as the similar characteristic between a plurality of filters in the bank of filters with it, make a general reference have between the coefficient of each filter absolute value identical or similarly (become 2 nRelation doubly, n is an integer, down with) number, but the difference that puts in order.So, after the tapped delay line of sharing input signal, at the weighted sum computing
Figure G2008101190648D00071
If can make full use of the similar characteristic between the filter, just can be by the weighted sum computing that repeats between multiplexing each filter, significantly reduce the operand of the inner repeated and redundant of bank of filters, thereby reduce its hard-wired complexity, save valuable hardware resource.For this reason, be necessary to propose a kind of new general implementation method and optimize the above-mentioned Finite Impulse Response filter group of realization, the independent realization of a plurality of filters in the replacement group.
The present invention is referred to as the low complex degree implementation method of Finite Impulse Response filter group, as shown in Figure 6, optimize the filtering operation that multiplexer is finished a plurality of filters, reduced hard-wired complexity with a similar filter, save hardware resource, improved the flexibility and the adaptability of system simultaneously.
For convenience of description, the bank of filters implementation method of the present invention's proposition need at first be carried out description below:
(1) radix: the coefficient of a Finite Impulse Response filter is h[n], 0≤n≤N, N are the exponent number of Finite Impulse Response filter, wherein the absolute value of arbitrary coefficient
Figure G2008101190648D00081
I ∈ n, A iBe odd number, m iBe nonnegative integer, here A iRadix for this coefficient h [i].
(2) kernel: the set of the radix of all coefficients of a Finite Impulse Response filter is the kernel of this filter, and the same cardinality of different coefficients is defined as the element of inequality in this set.
(3) similar filter: the similar each other filter of a plurality of Finite Impulse Response filters that kernel is identical.
(4) kernel filters: any filter that constitutes its coefficient that sorts of the radix by kernel is a kernel filters.
(5) the similar filter of kernel: the identical a plurality of kernel filters of kernel are the similar filter of kernel each other.
Below describe the concrete enforcement of Finite Impulse Response filter group implementation method and device thereof in detail, the Finite Impulse Response filter group of forming with two similar filters is that example is set forth technical scheme of the present invention.
If the input signal of Finite Impulse Response filter group is x[n], the coefficient of its median filter I is h 1[n], output signal is y 1[n], then its time domain linear filtering expression formula is:
y 1 [ n ] = x [ n ] * h 1 [ n ] = Σ m = 0 N x [ n - m ] · h 1 [ m ] - - - ( 3 )
* represents linear convolution in the following formula, x[n-m] be delay input signal, be expressed as input vector
Figure G2008101190648D00083
Filter coefficient is expressed as coefficient vector
Figure G2008101190648D00084
According to the vector multiplication rule, (3) formula can be converted into:
y 1 [ n ] = x → · h 1 → - - - ( 4 )
For the sake of simplicity, establish the kernel filters coefficient h of this bank of filters KernelThe ordering of [n] just in time is to construct according to the ordering of the coefficient of correspondence of filter I, and vector representation is
Figure G2008101190648D00091
H then 1[n] and h KernelThe pass of [n] is:
h 1 → = ( h 1 [ 0 ] , h 1 [ 1 ] , . . . , h 1 [ N ] ) = ( ± h kernel [ 0 ] × 2 p 0 , ± h kernel [ 1 ] × 2 p 1 , . . . , ± h kernel [ N ] × 2 p N ) - - - ( 5 )
" ± " number expresses possibility and to get one of positive negative value, p in the following formula 0, p 1..., p NBe nonnegative integer.According to (5) formula, construct new input vector
Figure G2008101190648D00093
Its " ± " number and
Figure G2008101190648D00094
With (5) formula
Figure G2008101190648D00095
Corresponding fully, then be not difficult to find out by the vector multiplication rule:
x 1 ′ → · h → kernel = x → · h → 1 = y 1 [ n ] - - - ( 6 )
If the coefficient of filter II is h in the group 2[n], output signal is y 2[n].Because two filter is similar filter, kernel is identical, so the coefficient of filter II can pass through the kernel filters coefficient h KernelThe rearrangement of [n], sign change and operation such as shift left is represented.Corresponding coefficient vector:
h 2 → = ( ± h kernel [ i 0 ] × 2 q 0 , ± h kernel [ i 1 ] × 2 q 1 , . . . , ± h kernel [ i N ] 2 q N ) - - - ( 7 )
" ± " number expresses possibility and to get one of positive negative value, i in the following formula 0, i 1..., i NMutually unequal and constitute 0 to N integer set, q 0, q 1..., q NBe nonnegative integer, then the output y of filter II 2[n] can be expressed as:
y 2 [ n ] = x → · h 2 → - - - ( 8 )
Simultaneously, construct new input vector
Figure G2008101190648D00099
Its " ± " number and
Figure G2008101190648D000910
With (7) formula Corresponding fully, then can get by the vector multiplication rule:
x → 2 ′ · h → kernel = x → · h 2 → = y 2 [ n ] - - - ( 9 )
By (6) formula and (9) Shi Kede, with the input vector of the similar characteristic between the similar filter coefficient " transfer " to tapped delay line
Figure G2008101190648D000913
On, according to the relation of each filter coefficient and kernel filters coefficient, change accordingly
Figure G2008101190648D000914
In element ordering and sign and shift left accordingly, structure is at the new input vector of kernel filters Then the different filtering operation of each filter in the weighted sum unit realization group of reusable kernel filters obtains filtering result separately.What deserves to be explained is, in the middle of hardware such as programmable logic device, digital signal processor spare and application-specific integrated circuit (ASIC) are realized, above-mentioned for input vector
Figure G2008101190648D00101
Conversion, comprise that rearrangement, sign change and shift left etc., can realize that these need take few hardware resource by corresponding line and displacement, compared to because of optimizing the required amount of hardware resources that takies of the multiplexing filtering operation that saves, can ignore fully.
With reference to Fig. 7, Fig. 8, represent the operating process schematic diagram of the Finite Impulse Response filter group low complex degree implementation method of forming by similar filter provided by the invention and the structural representation that adopts the implement device (M=2) of this method respectively.Finite Impulse Response filter group low complex degree implement device is made up of modules such as the parallel output units of clock source unit, tapped delay line unit, vector transduced cell, gating unit, kernel filters weighted sum unit and kernel filters.
The clock source unit provides the clock of f and two kinds of frequencies of 2f.
Former filter all is operated under the sample frequency f separately, described tapped delay line unit, the frequency that provides in the clock source is under the clock of f drives, with input signal x[n] delaying time obtains input time delay signal x[n-m], vector representation is
Figure G2008101190648D00102
It is input vector.
Described vector transduced cell is one of core cell of Finite Impulse Response filter group low complex degree implement device.Similar characteristic in the Finite Impulse Response filter group between each filter is fully transferred to input vector
Figure G2008101190648D00103
On, realize in this unit.In described vector transduced cell, the linear convolution that is complementary according to each filter and kernel filters in the group concerns With
Figure G2008101190648D00105
The former input vector that tapped delay line is obtained
Figure G2008101190648D00106
Change, by changing
Figure G2008101190648D00107
In ordering and sign and the operation such as shift left accordingly of element, obtain new input vector respectively at kernel filters
Figure G2008101190648D00108
With
Figure G2008101190648D00109
Make
Figure G2008101190648D001010
With
Figure G2008101190648D001011
Corresponding relation mate fully
Figure G2008101190648D001012
With
Figure G2008101190648D001013
Corresponding relation, With
Figure G2008101190648D001015
Corresponding relation mate fully
Figure G2008101190648D001016
With
Figure G2008101190648D001017
Corresponding relation.
Described gating unit, the frequency that provides in the clock source are under the clock of 2f drives, to described two groups of new input vectors
Figure G2008101190648D001018
With
Figure G2008101190648D001019
Carry out the timesharing gating.For example will
Figure G2008101190648D001020
Propose the clock cycle gating output of previous frequency 2f, form the input signal of kernel filters weighted sum unit.
With reference to Fig. 9, represent the principle schematic of kernel filters weighted sum unit provided by the invention, this unit is one of core cell of Finite Impulse Response filter group low complex degree implement device.Under the sample frequency of the 2f that this cell operation provides in the clock source,, finish respectively by the weighted sum computing structure of time division multiplexing kernel filters
Figure G2008101190648D00111
With The weighted sum computing.For example, import in the clock cycle of previous frequency 2f
Figure G2008101190648D00113
With the kernel filters coefficient vector
Figure G2008101190648D00114
Be weighted summation operation, the filtering operation that obtains filter I is y as a result 11[n]; In the clock cycle of a back frequency 2f, input
Figure G2008101190648D00115
With the kernel filters coefficient vector
Figure G2008101190648D00116
Be weighted summation operation, the filtering operation that obtains filter II is y as a result 2[n].
With reference to Figure 10, represent the principle schematic of the parallel output unit of kernel filters provided by the invention, this unit is operated under the sample frequency of the 2f that the clock source provides equally.If
Figure G2008101190648D00117
Figure G2008101190648D00118
Carry the clock cycle of previous frequency 2f and finish the weighted sum computing, then this unit is with the operation result y of filter I 11The clock cycle of a frequency 2f of [n] time-delay, obtain operation result y with filter II 2The time delayed signal y that [n] is synchronous 1[n] is with y 1[n] and y 2[n] and line output, realize described bank of filters two filter and the line output demand.
With reference to Figure 11, represent the structural representation of the Finite Impulse Response filter group low complex degree implement device of forming by two incomplete similarity filters provided by the invention.Finite Impulse Response filter group implementation method provided by the invention and device thereof are equally applicable to the realization of the bank of filters be made up of the incomplete similarity filter.The bank of filters that constitutes with two filters is an example, if the kernel major part of two filters is identical, and the similar bank of filters implementation method that then can adopt the present invention to propose equally.Concrete operations can be as follows: with the kernel of filter A kernel as bank of filters, and the structure kernel filters, coefficient vector is made as
Figure G2008101190648D00119
Then can realize filter A with kernel filters as previously mentioned; Utilize the definition of kernel filters, search for the similar filter set that the similar filter of all kernels constitutes,, obtain the similar filter B of kernel to filter B coefficient error minimum according to the filter error criterion Kernel, then can realize filter B with kernel filters as previously mentioned KernelAccording to supposition, the kernel major part of two filters is identical, so filter B KernelAnd the coefficient error between the filter B is very little, can directly realize filter B according to input vector KernelAnd the error-filter B between the filter B Error, filter B KernelWith error-filter B ErrorThe output sum is exactly the output of filter B.Above-mentioned filter error criterion is meant the mean square error minimum criteria of filter time domain impulse response or the mean square error minimum criteria of filter frequency domain response.
The low complex degree implement device of Finite Impulse Response filter group also comprises in the present embodiment: error-filter weighted sum unit and summation output unit.Wherein, the error-filter weighted sum unit is under the f frequency that the clock source unit provides, to input vector
Figure G2008101190648D00121
Carry out and the error-filter coefficient vector
Figure G2008101190648D00122
The weighted sum computing, obtain error-filter B ErrorOperation result; The summation output unit will run simultaneously output kernel filters corresponding to B KernelOperation result and error-filter B ErrorThe operation result addition, obtain operation result and the output of filter B in the former Finite Impulse Response filter group.
With reference to Figure 12, represent the structural representation of the Finite Impulse Response filter group low complex degree implement device of forming by a plurality of incomplete similarity filters provided by the invention.Finite Impulse Response filter group implementation method provided by the invention and device thereof are convenient to described Finite Impulse Response filter group is expanded to plural filter very much.If bank of filters comprises the individual filter of M (M 〉=2), then need similar characteristic according to each filter in the group, extract the same section of each filter kernel, construct one with the kernel of the error minimum of each filter kernel kernel as bank of filters, and (coefficient vector is made as to constitute kernel filters
Figure G2008101190648D00123
), preceding for another example described, (coefficient vector is made as by searching for the similar filter of kernel that obtains respectively to organizing interior each filter coefficient error minimum
Figure G2008101190648D00124
) and corresponding with it error-filter (coefficient vector is made as ).Wherein the similar bank of filters of kernel realizes according to similar bank of filters implementation method and device.
The low complex degree implement device of Finite Impulse Response filter group comprises in the present embodiment: the clock source unit provides the clock of f and two kinds of frequencies of Mf; The tapped delay line unit in the frequency that the clock source unit provides is under the clock driving of f, with input signal x[n] delaying time obtains input time delay signal x[n-m], be expressed as
Figure G2008101190648D00126
Input vector; Vector transduced cell is used for the input vector that the tapped delay line unit is obtained Change, by changing
Figure G2008101190648D00128
In element ordering and sign and shift left accordingly and wait operation, obtain M input vector newly
Figure G2008101190648D00129
1≤i≤M makes
Figure G2008101190648D001210
With
Figure G2008101190648D001211
Corresponding relation mate the similar filter coefficient vector of kernel fully
Figure G2008101190648D00131
With the kernel filters coefficient vector
Figure G2008101190648D00132
Corresponding relation, also be
Figure G2008101190648D00133
Gating unit is under the Mf frequency that the clock source unit provides, to M new input vector
Figure G2008101190648D00134
The timesharing gating; The kernel filters weighted sum unit under the Mf frequency that the clock source unit provides, is exported from gating unit
Figure G2008101190648D00135
Timesharing is carried out and the kernel filters coefficient vector
Figure G2008101190648D00136
The weighted sum computing, obtain M group filtering operation result; The error-filter weighted sum unit is under the f frequency that the clock source unit provides, to input vector
Figure G2008101190648D00137
Carry out and the error-filter coefficient vector
Figure G2008101190648D00138
The weighted sum computing; The error-filter coefficient vector
Figure G2008101190648D00139
Be to obtain by the Finite Impulse Response filter decomposition step, be specially from the set that the similar filter of all kernels similar to above-mentioned kernel filters constitutes, according to the filter error criterion, by search, obtain respectively and the similar filter of kernel of organizing interior each Finite Impulse Response filter coefficient error minimum, constitute coefficient vector and be respectively
Figure G2008101190648D001310
The similar bank of filters of kernel, 1≤i≤M, and each Finite Impulse Response filter in former group is decomposed into the similar filter of kernel of coefficient error minimum with it, and (coefficient vector is
Figure G2008101190648D001311
) (coefficient vector is with error-filter corresponding to the similar filter of this kernel ) two parts, the filtering operation of former group of interior each Finite Impulse Response filter is y as a result i[n] also correspondingly is decomposed into the similar filter operation result of kernel y Kerneli[n] and corresponding error-filter operation result y Errori[n] two parts; The kernel filters output unit that walks abreast, under the Mf frequency that the clock source unit provides, the M group filtering operation result that kernel filter weight sum unit is obtained arrives M-1 the clock cycle of not waiting with frequency Mf time-delay 0, the described M group of the output of running simultaneously filtering operation result; The summation output unit with kernel filters operation result and its each self-corresponding error-filter operation result addition of the above-mentioned output of running simultaneously, obtains the filtering operation result and the output of each Finite Impulse Response filter.
The present invention makes full use of the similar characteristic between each filter in the Finite Impulse Response filter group, deeply excavate optimization potentiality wherein, the multiplexing to greatest extent more kernel filtering operation structure of resource that takies is optimized and is realized particularly higher order filter group of Finite Impulse Response filter group.Compare with traditional implementation method, effectively reduce hard-wired complexity, reduced hardware resource and taken, improved the flexibility and the applicability of system, and be convenient to expand, realize the bank of filters of a plurality of Finite Impulse Response filters compositions.
The above only is a preferred implementation of the present invention; should be pointed out that for those skilled in the art, under the prerequisite that does not break away from the technology of the present invention principle; can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.

Claims (5)

1.FIR the low complex degree implement device of digital filter bank, this Finite Impulse Response filter group comprise M N rank Finite Impulse Response filter, wherein M is characterized in that more than or equal to 2, and this device comprises:
The clock source unit provides the clock of f and two kinds of frequencies of Mf;
The tapped delay line unit in the frequency that the clock source unit provides is under the clock driving of f, with input signal x[n] delaying time obtains input time delay signal x[n-m], be expressed as
Figure F2008101190648C00011
Input vector;
Vector transduced cell is used for the input vector that the tapped delay line unit is obtained Change, by changing
Figure F2008101190648C00013
In element ordering and sign and shift left accordingly, obtain the new input vector of M
Figure F2008101190648C00014
Make
Figure F2008101190648C00015
With Corresponding relation and similar filter coefficient vector With the kernel filters coefficient vector
Figure F2008101190648C00018
Corresponding relation mate fully, promptly
Figure F2008101190648C00019
1≤i≤M;
Gating unit is under the Mf frequency that the clock source unit provides, to M new input vector
Figure F2008101190648C000110
The timesharing gating;
The kernel filters weighted sum unit is under the Mf frequency that the clock source unit provides, to what export from gating unit
Figure F2008101190648C000111
Timesharing is carried out and the kernel filters coefficient vector
Figure F2008101190648C000112
The weighted sum computing, obtain M group filtering operation result;
The kernel filters output unit that walks abreast, under the Mf frequency that the clock source unit provides, M group filtering operation result to kernel filter weight sum unit obtains with 0 to M-1 clock cycle of not waiting of frequency Mf time-delay, runs simultaneously and exports M group filtering operation result.
2. the low complex degree implement device of Finite Impulse Response filter group as claimed in claim 1 is characterized in that, this device also comprises:
The error-filter weighted sum unit is under the f frequency that the clock source unit provides, to input vector
Figure F2008101190648C000113
Carry out and the error-filter coefficient vector
Figure F2008101190648C000114
The weighted sum computing;
The summation output unit with M group filtering operation result and its each self-corresponding error-filter operation result addition that the parallel output unit of kernel filters is run simultaneously and exported, obtains the filtering operation result and the output of each Finite Impulse Response filter.
3. implementation method of utilizing the described device of claim 1 to reach Finite Impulse Response filter group low complex degree, this Finite Impulse Response filter group comprises M N rank FIR similarity number character filter, wherein M is characterized in that more than or equal to 2, may further comprise the steps:
The input signal x[n of described Finite Impulse Response filter group] obtain input vector through tapped delay line
Figure F2008101190648C00021
The coefficient vector of each Finite Impulse Response filter in the group is respectively
Figure F2008101190648C00022
1≤i≤M;
The kernel filters constitution step: the kernel total according to similar filter constitutes kernel filters, and its coefficient vector is
Vector switch process: to input vector
Figure F2008101190648C00024
Change, by changing
Figure F2008101190648C00025
In element ordering and sign and shift left accordingly, structure is at the M of kernel filters new input vector
Figure F2008101190648C00026
Make With Corresponding relation mate fully
Figure F2008101190648C00029
With
Figure F2008101190648C000210
Corresponding relation, promptly
Figure F2008101190648C000211
1≤i≤M;
The multiplexing step of kernel filters weighted sum: described Finite Impulse Response filter group is multiplexing kernel filters weighted sum unit in the different time periods, the filtering operation of each Finite Impulse Response filter in the realization group.
4. implementation method of utilizing the described device of claim 2 to reach Finite Impulse Response filter group low complex degree, this Finite Impulse Response filter group comprises M N rank Finite Impulse Response filter, wherein M is characterized in that more than or equal to 2, may further comprise the steps:
The input signal x[n of described Finite Impulse Response filter group] obtain input vector through tapped delay line
Figure F2008101190648C000212
The coefficient vector of each Finite Impulse Response filter is respectively in the group
Figure F2008101190648C000213
1≤i≤M;
Kernel filters constitution step: according to the similar characteristic of each Finite Impulse Response filter in the group, extract the same section of each Finite Impulse Response filter kernel, construct one with the kernel of the error minimum of each Finite Impulse Response filter kernel kernel as the Finite Impulse Response filter group, and the formation kernel filters, its coefficient vector is
Figure F2008101190648C000214
Finite Impulse Response filter decomposition step: from the set that the similar filter of all kernels similar to above-mentioned kernel filters constitutes, according to the filter error criterion, by search, obtain respectively and the similar filter of kernel of organizing interior each Finite Impulse Response filter coefficient error minimum, constitute coefficient vector and be respectively
Figure F2008101190648C00031
The similar bank of filters of kernel, 1≤i≤M, and each Finite Impulse Response filter in former group is decomposed into the similar filter of kernel of coefficient error minimum with it, and (coefficient vector is ) (coefficient vector is with error-filter corresponding to the similar filter of this kernel
Figure F2008101190648C00033
) two parts, the filtering operation of former group of interior each Finite Impulse Response filter is y as a result i[n] also correspondingly is decomposed into the similar filter operation result of kernel y Kerneli[n] and corresponding error-filter operation result y Errori[n] two parts;
Vector switch process: to input vector
Figure F2008101190648C00034
Change, by changing In element ordering and sign and shift left accordingly, structure is at the M of the similar filter of kernel new input vector
Figure F2008101190648C00036
Make
Figure F2008101190648C00037
With
Figure F2008101190648C00038
Corresponding relation mate fully
Figure F2008101190648C00039
With
Figure F2008101190648C000310
Corresponding relation, promptly
Figure F2008101190648C000311
1≤i≤M;
The multiplexing step of kernel filters weighted sum: the similar bank of filters of kernel is multiplexing kernel filters weighted sum unit in the different time periods, realizes the filtering operation of the similar filter of each kernel in the similar bank of filters of kernel;
Error-filter weighted sum step: the corresponding error-filter of each Finite Impulse Response filter is weighted summation, realizes the filtering operation of each error-filter
Figure F2008101190648C000312
1≤i≤M;
Summation output step: the error-filter that the output of the similar filter of each kernel in the similar bank of filters of above-mentioned kernel is corresponding with it is exported addition, thereby obtains the output of each Finite Impulse Response filter in the Finite Impulse Response filter group:
y i [ n ] = y kerneli [ n ] + y errori [ n ] = x → · h → kerneli + x → · h → errori = x i ′ → · h → kernel + x → · h → errori , 1≤i≤M。
5. the low complex degree implementation method of Finite Impulse Response filter group as claimed in claim 4 is characterized in that, described filter error criterion is meant the mean square error minimum criteria of filter time domain impulse response or the mean square error minimum criteria of filter frequency domain response.
CN2008101190648A 2008-08-28 2008-08-28 Low-complexity implementing method and apparatus for FIR digital filter group Expired - Fee Related CN101340182B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101190648A CN101340182B (en) 2008-08-28 2008-08-28 Low-complexity implementing method and apparatus for FIR digital filter group

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101190648A CN101340182B (en) 2008-08-28 2008-08-28 Low-complexity implementing method and apparatus for FIR digital filter group

Publications (2)

Publication Number Publication Date
CN101340182A CN101340182A (en) 2009-01-07
CN101340182B true CN101340182B (en) 2010-06-02

Family

ID=40214162

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101190648A Expired - Fee Related CN101340182B (en) 2008-08-28 2008-08-28 Low-complexity implementing method and apparatus for FIR digital filter group

Country Status (1)

Country Link
CN (1) CN101340182B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101661407B (en) * 2009-09-30 2013-05-08 中兴通讯股份有限公司 Finite impulse response filter with parallel structure and processing method thereof
US10135424B2 (en) * 2015-02-27 2018-11-20 Microchip Technology Germany Gmbh Digital filter with confidence input
CN105827217B (en) * 2016-03-03 2018-10-12 深圳市紫光同创电子有限公司 One kind having limit for length's impulse response filter circuit and programmable logic device
CN106301373B (en) * 2016-08-26 2018-07-27 中国科学院地质与地球物理研究所 A kind of number multimode multi-frequency section filter group and electromagnetic method receiver
CN110138358A (en) * 2019-04-30 2019-08-16 南京大学 A kind of long linear phase limited impulse response digital filter of idol
CN116318051B (en) * 2023-03-16 2024-02-27 湖南迈克森伟电子科技有限公司 Digital shaping filter method and device, digital shaping filter and electronic equipment

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5311457A (en) * 1991-04-17 1994-05-10 Lsi Logic Kabushiki Kaisha Digital filter
CN2674763Y (en) * 2003-12-02 2005-01-26 天津大学 Full phase DFT digital filter

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5311457A (en) * 1991-04-17 1994-05-10 Lsi Logic Kabushiki Kaisha Digital filter
CN2674763Y (en) * 2003-12-02 2005-01-26 天津大学 Full phase DFT digital filter

Also Published As

Publication number Publication date
CN101340182A (en) 2009-01-07

Similar Documents

Publication Publication Date Title
CN101340182B (en) Low-complexity implementing method and apparatus for FIR digital filter group
US3979701A (en) Non-recursive digital filter employing simple coefficients
CN101567701B (en) High efficient multi-path digital down converter system
CN102098004A (en) Digital downconverter with variable bandwidth and implementation method thereof
US10003324B2 (en) Fast FIR filtering technique for multirate filters
Cheng et al. Further complexity reduction of parallel FIR filters
EP1262019B1 (en) Apparatus for splitting the frequency band of an input signal
CN102064797B (en) Parallel implementation method and device for fractional sampling rate transformation
CN102403986B (en) Multi-channel CIC (Cascade Integrator Comb) decimation filter and method for realizing same
CN102025377B (en) Improved cascaded integral comb interpolation filter
US7680869B2 (en) Interpolation and decimation using newton polyphase filters
Srivastava et al. An efficient fir filter based on hardware sharing architecture using csd coefficient grouping for wireless application
CN103117730A (en) Multichannel comb filter and implement method thereof
CN105811920B (en) A kind of FRM narrow transition band filter group structure
CN102891662B (en) A kind of general rate down-conversion, up conversion device and method
CN103237219A (en) Two-dimensional discrete cosine transformation (DCT)/inverse DCT circuit and method
CN101360087A (en) Low-complexity implementing method and apparatus for base-band forming SRRC digital filter
CN101242168A (en) A realization method and device for FIR digital filter direct-connection
Mayilavelane et al. A Fast FIR filtering technique for multirate filters
CN101610077B (en) Filter of half-band filtering implementation method
CN103078606A (en) Multichannel CIC (Cascaded Integrator Comb) interpolation-filter system and implementation method thereof
US8295388B1 (en) Efficient digital upconverter
Wu et al. An efficient architecture for two-dimensional inverse discrete wavelet transform
Lavanya et al. High speed, low complexity, folded, polymorphic wavelet architecture using reconfigurable hardware
Varghese et al. Design of Reconfigurable Filter Structure Based on FRM for Wideband Channelizer

Legal Events

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

Granted publication date: 20100602

Termination date: 20160828