The invention relates to a multirate filter as well as a display system and a mobile phone comprising a multirate filter.
Digital filters find widespread use in audio and video processing, for example in mobile phones, set top boxes, digital television sets, and other consumer, or professional products. Symmetric filters form an important class, because of their linear phase property and the possibility to exploit this symmetry to simplify the architecture of the filter and therewith reduce costs.
In particular, multirate filters are used in applications where the output signal and the input signal of the filter should have mutually different sample rates. Such filters are applied for example in image processing to effect a scaling of a digitally encoded image. One of the most important concepts in multirate filtering is the polyphase decomposition and the closely related polyphase structure. This concept allows for very efficient implementations both in hardware and software of interpolating and decimating filters.
A straightforward application of the polyphase decomposition however usually introduces asymmetric polyphase components, which substantially reduces the efficiency of the implementation
It is a purpose of the invention to provide an architecture for a multirate filter wherein the symmetry of its components is recovered. For that purpose the multirate filter according to the invention has a construction as described herein. It has been recognized by the inventors that a symmetric multirate may be constructed from an input unit, a filter unit and an output unit, wherein the filter unit has symmetric modules derived from the polyphase components of the multirate filter provided that the modules are provided in pairs having transfer functions H0(z) and H1(z), which are derived from a basic transfer function HB(z) as follows. The first one H0(z) of the transfer functions is based on the sum of the basic transfer function HB(z) and its mirrored version:
H0(z)=c0(HB(z)+Mα,ψHB(z)), and
the second one H1(z) is based on the difference of the basic transfer function HB(z) and its mirrored version:
H1(z)=c1(HB(z)−Mα,ψHB(z)).
The mirror operation Mα,ψon the basic transfer function HB(z) is defined as:
In one embodiment, the basic transfer function HB(z) is a polyphase component of the multirate filter. In another embodiment, the input unit comprises a combination unit and in another embodiment, the output unit comprises such a combination unit. In yet another embodiment, the basic transfer function HB(z) is a linear function of two polyphase components is of the multirate filter.
The invention further relates to a method of designing a multirate filter according to the invention. Such a method could either be part of a design tool, but could alternatively be part of a device having an adaptive filter. In that case the method enables that device to replace the polyphase components which it has calculated by and equivalent set of symmetrical modules.
These and other aspects of the invention are described in more detail with reference to the drawings. Therein:
Likewise the multirate filter according to the invention is applicable in the sending part of a mobile telephone. The architecture of the sending part (not shown) of a mobile telephone is globally the inverse of the receiving part. Hence, such a mobile telephone comprises an analog to digital converter for converting an analog speech signal into a digital speech signal. A speech encoder subsequently compresses the digital speech signal and provides the compressed signal. A signal encoder performs a channel encoding operation at the compressed signal and provides a channel signal. A modulator modulates the channel signal. The channel signal is filtered by a multirate filter according to the invention as described in more detail below. The filtered signal is then converted into an analog signal by a digital to analog converter. A transmitter transmits the analog signal.
In the present description a filter h[n] having z-transform H(z) is defined as (α,ψ)-symmetric if it fulfills the following relationship:
H(z)=αz−2ψH
which corresponds to
A special case is (α=1,ψ), h[n] being a real-valued function. For these filters h[n] the following relation holds:
Another special case is (α=−1,ψ). For these filters either h[2ψ−i]=−h[i] wherein i is an integer,
These filters are also referred to as anti-symmetric filters, not to be confused with asymmetric filters.
Examples are:
The sign of symmetry is however not limited to real values, as is illustrated by the following examples:
Filters having (α,ψ)-symmetry as defined above can be implemented efficiently. It is known to exploit their symmetry by reducing the number of components. This illustrated in
Analogously
Another known technique in filter design is the decomposition into polyphase components. This is in particular relevant for decimating or interpolating filters, or filters using a combination of decimation and interpolation.
Any filter H(z) can be decomposed in R polyphase components HR:r(z), which are related to H(z) as:
wherein the polyphase components relate to the transfer function h( ) as:
Such a polyphase decomposition of the filter has the advantage that the polyphase filter components can operate at a lower clock speed than a filter which is not decomposed.
This is illustrated in
Unfortunately, even a symmetric filter may have asymmetric polyphase components, as is illustrated by the following example:
Consider the filter h[n] having symmetric impulse response <a,b,c,d,d,c,b,a>. In case of multirate factor R=2, both its polyphase components h2:0[n] and h2:1[n] are non-symmetric:
It has been found by the inventors that the polyphase components of the multirate filter can be combined in (α,ψ) symmetric filter modules. The general structure of such a multirate filter according to the invention is illustrated in
The multirate filter according to the invention shown therein has:
The filter unit 20 comprises at least a first and a second filter module 21, 22, having a transfer function H0(z) and a transfer function H1(z) respectively, which are mutually related according to the relations
H0(z)=c0(HB(z)+Mα,ψHB(z)) and
H1(z)=c1(HB(z)−Mα,ψHB(z)).
Therein Mα,ψHB(z)=αz−2ψHB
By combining a non-symmetric filter HB(z) with its mirrored version Mα,ψHB(z), either a symmetric filter (in case of combining by addition) or an anti-symmetric filter (in case of combining by subtraction) is obtained. This is illustrated by the following example. When adding the two non-symmetric polyphase components
The multirate filter comprises a combination unit 11 coupled to the filter modules 21, 22 for generating a first combination signal Ssum and a second combination signal Sdiff. The combination unit 11 may be a part of the input unit 10, or a part of the output unit 30. Alternatively both the input unit 10 and the output unit 30 may comprise one or more combination units.
In an embodiment of the invention the basic transfer function HB(z) is a polyphase component i.e.
As is clear from the above, the transfer functions of the filter modules 21 and 22 are symmetric. The transfer function Ho(z) of the first filter module is the sum of the basic transfer function HB(z) and its mirrored counterpart, and the transfer function H1(z) of the second filter module is the difference of the basic transfer function HB(z) and its mirrored counterpart. Hence, both transfer functions H0(z) and H1(z) allow an efficient implementation, having a relatively small number of multipliers. On the one hand the filter modules can operate relatively slowly as they process the signal Ssum and Sdiff, which are derived from the decimated signals IS1 and IS2. The optional interpolation stage only takes place after the filter unit 20.
The transfer functions H0(z) and H1(z) of the first and the second filter module are the sum of a basic transfer function HB(z) and its mirrored version Mα,ψHB(z). In the embodiments shown in
HB(z)=HR:r
The combination unit 211 is comprised in the input unit 210, and generates the first and the second combination signal Ssum, Sdiff from a pair of the said intermediate signals IS1, IS2. The first filter module 221 filters the first combination signal Ssum and the second filter module 222 filters the second combination signal Sdiff. The filter unit 220 further comprises a third filter module 223 with a transfer function H2(z) and a fourth filter module 224 with a transfer function H3(z) which are mutually related according to the relations
H2(z)=c2(H′B(z)−MH′B(z))
H3(z)=c3(H′B(z)+MH′B(z)).
Also the further basic transfer function HB′(z) is derived from a combination of two polyphase components:
H′B(z)=HR:r
The third 223 and fourth filter module 224 each filter a respective one of the first combination signal Ssum and the second combination signal Sdiff. The filter comprises a first further combination unit 231 for generating a first auxiliary signal AS1 representing the sum of the output signals of the first filter module 221 and the second filter module 222, and a second auxiliary signal AS2 representing the sum of the output signals of the third filter module 222 and the fourth filter module 223. The filter further comprises a second further combination unit 232 for generating a further sum signal Ssum″ and a further difference signal Sdiff′ from the first and the second auxiliary signal AS1, AS2. These signals are recombined into an output signal Sout after an optional interpolation with a factor I. The signals Ssum″ and Sdiff″ may further be delayed with an additional delay z−r
As in the embodiments shown in
The transfer functions for the second group are
It is noted that the skilled person will consider many alternatives within the scope of the invention as defined by the claims. For example the delay functions z2, z0 and z1 in the input unit may be replaced by delay functions z2−k, z−k and z1−k in combination with a delay function zk in the signal line providing the input signal X(z). It is noted that merely the relative delays generated between the signals provided by the delay functions in the input unit 310 is of importance.
It is also clear to the skilled person that several functions may be interchanged. For example, the decimate units indicated by the downpointing arrows with label 3 may be located at the outputs of the combination unit 311. Furthermore the delay units indicated with z0 and z1 may be included in the combination unit 311. In that case a first combination element of the combination unit 311 (the upper one in the drawing) generates the signal X(z)+z1X(z)=X(z)(1+z). The other combination element of said combination unit generates X(z)(1−z).
The skilled person will further realized that other operations can be interchanged for example the sequence of operations:
interpolating with a factor L
delaying with za,
decimating with a factor D
may be replaced by the sequence
delaying with zaq,
decimating with a factor D
interpolating with a factor I,
delaying with zap.
Provided that D and I are coprime.
Other examples are that multiplication operations can be interchanged with delay operations and with decimation or interpolation operations.
A multirate filter according to the invention can be designed by carrying out the following steps:
a first module with transfer function H0=HR:r1+HR:r2,
a second module with transfer function H1=HR:r1−HR:r2, wherein the first and the second module sharing the common input, and by
a combination unit for generating a first and a second combination signal from the output signals of the first and the second module.
a first module with transfer function H0=HR:r1+HR:r2,
a second module with transfer function H1=HR:r1−RR:r2, which modules share the common output, and by
a combination unit for generating a first combination signal from a first and a
second intermediate signal, and providing said combination signal to the third module, and for generating a second combination signal from a first and a second intermediate signal, and providing said combination signal to the fourth module.
a first pair HR:r1, HR:r2 which are related by HR:r2=Mα,ψ(HR:r1)
a second pair HR:r3, HR:r4 which are related by HR:r4=Mα,ψ(HR:r3) wherein,
the components HR:r1 and HR:r3 share a first common input,
the components HR:r2 and HR:r4 share a second common input,
the components HR:r1 and HR:r2 share a first common output, and
the components HR:r3 and HR:r4 share a second common output.
Each such a quadruplet can be replaced by a quadruplet of symmetrical filter modules and three combination units. The quadruplet of filter modules comprises
a first module having transfer function H0=HR:r0+HR:r1+HR:r2+HR:r3,
a second module having transfer function H1=HR:r0−HR:r1−HR:r2+HR:r3,
a third module having transfer function H2=HR:r0+HR:r1−HR:r2−HR:r3,
a fourth module having transfer function H3=HR:r0−HR:r1+HR:r2−HR:r3.
A first combination unit generates a first combination signal from input signals received at the first and the second common input, and provides said first combination signal to the first and the second module. It also generates a second combination signal from those input signals, and provides that to the third and the fourth module.
A second combination unit generates a first auxiliary signal from output signals generated by the first and the third unit. It generates a second auxiliary signal from output signals generated by the second and the fourth unit.
A third combination unit generates a first and a second output signal from the first and the second auxiliary signal.
It is noted that the situations b, c and d may occur in combinations.
Number | Date | Country | Kind |
---|---|---|---|
02080103 | Dec 2002 | EP | regional |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB03/05647 | 12/4/2003 | WO | 00 | 6/1/2005 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2004/054103 | 6/24/2004 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5103310 | Gibson et al. | Apr 1992 | A |
5881107 | Termerinac et al. | Mar 1999 | A |
6279019 | Oh et al. | Aug 2001 | B1 |
6834292 | Jiang et al. | Dec 2004 | B2 |
7161979 | Wildhagen | Jan 2007 | B2 |
Number | Date | Country |
---|---|---|
44 02 632 | Aug 1995 | DE |
Number | Date | Country | |
---|---|---|---|
20060056553 A1 | Mar 2006 | US |