DE10317698B4 - Verfahren zum Entwerfen von Polynomen zur Steuerung des Veränderns von anpassungsfähigen Digitalfiltern - Google Patents

Verfahren zum Entwerfen von Polynomen zur Steuerung des Veränderns von anpassungsfähigen Digitalfiltern Download PDF

Info

Publication number
DE10317698B4
DE10317698B4 DE10317698A DE10317698A DE10317698B4 DE 10317698 B4 DE10317698 B4 DE 10317698B4 DE 10317698 A DE10317698 A DE 10317698A DE 10317698 A DE10317698 A DE 10317698A DE 10317698 B4 DE10317698 B4 DE 10317698B4
Authority
DE
Germany
Prior art keywords
filter
coefficients
polynomial
cutoff frequency
characteristic
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
DE10317698A
Other languages
English (en)
Other versions
DE10317698A1 (de
Inventor
Mark W. Trenton Corless
Sunil Dearborn Shukla
J. William Novi Whikehart
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.)
Visteon Global Technologies Inc
Original Assignee
Visteon Global Technologies Inc
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 Visteon Global Technologies Inc filed Critical Visteon Global Technologies Inc
Publication of DE10317698A1 publication Critical patent/DE10317698A1/de
Application granted granted Critical
Publication of DE10317698B4 publication Critical patent/DE10317698B4/de
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0294Variable filters; Programmable filters
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0248Filters characterised by a particular frequency response or filtering method
    • H03H17/028Polynomial filters
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/04Recursive filters
    • H03H17/0416Recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0294Variable filters; Programmable filters
    • H03H2017/0295Changing between two filter characteristics
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H2218/00Indexing scheme relating to details of digital filters
    • H03H2218/04In-phase and quadrature [I/Q] signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Filters That Use Time-Delay Elements (AREA)

Abstract

Verfahren zur Signalverarbeitung mit Digitalfiltern (37, 38; 108, 109) in einer elektronischen Einheit mit einem Digitalsignalprozessor (25), wobei das Digitalfilter (37, 38; 108, 109) zum Erreichen einer zwischen einer ersten Charakteristik mit einer ersten Grenzfrequenz (81) und einer zweiten Charakteristik mit einer zweiten Grenzfrequenz (80) liegenden Bandweitencharakteristik mit einer dritten Grenzfrequenz eine Vielzahl von konfigurierbaren Filtersektionen (46, 47, 48) besitzt, wobei den Filtersektionen (46, 47, 48) Filterkoeffizienten (a1, a2, b0, b1, b2) zugeordnet sind, wobei die Filterkoeffizienten (a1, a2) mit folgenden Polynomen a1 = p0,1 + (p1,1·i) + (p2,1·i2) (I) a2 = p0,2 + (p1,2·i) + (P2,2·i2) (II)bestimmt werden,
wobei i der Bandweitenpositionsindex als Variable in einem Speicher (105) einer Steuerung (104) des Digitalsignalprozessors (25) einem Abtastfrequenzbereich (fS) zwischen der ersten Grenzfrequenz (81, i = 1) und der zweiten Grenzfrequenz (80, i = 0) mit den Indexwerten 0 ≤ i ≤ 1 entspricht und wobei der dritten Bandweitencharakteristik mit der dritten Grenzfrequenz...

Description

  • Querverweis auf entsprechende Patentanmeldungen
  • Diese Anmeldung bezieht sich auf die ebenfalls anstehende US-Patentanmeldung mit der Seriennummer US 09/247,214 , nunmehr US 6,445,735 , eingereicht am 8. Februar 1999 unter dem Titel „Digitalfilter mit geschalteter Bandweite mit reduzierten Übergängen während des Schalten” und die ebenfalls anstehende US-Patentanmeldung mit der Seriennummer US 09/382,727, nunmehr US 6,393,450, eingereicht am 25. August 1999 unter dem Titel „Digitalfilter mit geschalteter Bandweite mit Schwenkung”.
  • Diese Anmeldung bezieht sich außerdem auf die gleichzeitig unter dem Titel „Erzeugung von Laufzeitkoeffizienten für Digitalfilter mit schwenkbarer Bandweite” eingereichten US-Patentanmeldung mit der Seriennummer 10/122894, nunmehr US 7,027,502 und die gleichzeitig unter dem Titel „Verfahren zum Schwenken eines Digitalfilters, das Filtersektionen mit passender Verstärkung bereitstellt” eingereichten US-Patentanmeldung mit der Seriennummer 10/122529, nunmehr US 2002/0136289 A1.
  • Hintergrund der Erfindung
  • Die Erfindung bezieht sich allgemein auf einstellbare, in der Digitalsignalverarbeitung verwendete Digitalbandfilter und insbesondere auf das Entwerfen von Polynomen zum Erzeugen von Filterkoeffizientenwerten bei verschiedenen Filtergrenzfrequenzen, die zum Verändern der Bandbreite eines Digitalfilters verwendet werden.
  • Digitalfilter manipulieren diskrete Stichproben eines Eingangssignals zur Erzeugung eines gefilterten Ausgangssignals. Im Fachgebiet sind verschiedene Filterstrukturen bekannt, wie z. B. solche von finiten Transversalfiltern (FIR) und infiniten Transversalfiltern (IIR). IIR-Filter höherer Ordnung (sorgen für größere Wahlmöglichkeit) sind typischerweise mithilfe einer Vielzahl von in Kaskade geschalteten Filtern niedrigerer Ordnung implementiert.
  • Während der Verarbeitung eines Signals kann eine Änderung beim Filtern des Signals erforderlich werden (d. h. eine Änderung der Bandbreite, Durchlasscharakteristik, Gruppenlaufzeit oder anderer Filterparameter). Zur Minimierung von Anforderungen an Hardware und/oder Software ist es wünschenswert, dass durch Änderung lediglich der Filterkoeffizienten vor und nach der Änderung dieselbe Filterstruktur verwendet werden kann.
  • In einem digitalsignalverarbeitenden (DSP) Radioempfänger wird zum Beispiel ein Digitalkanalfilter auf ein Zwischenfrequenz-(IR-)Signal zum Auswählen des gewünschten Signals und Sperren anderer Rundfunksignale angewendet. Für das Kanalfilter kann in Abhängigkeit des Vorhandenseins von störenden benachbarten oder wechselnden Kanalrundfunksignalen ein breiter oder schmaler Durchlassbereich verwendet werden. Beim Umschalten zwischen den zwei Bandbreiten durch Umschalten zwischen zwei Koeffizientensätzen in einem DSP-Filter kann die plötzliche Änderung der Bandbreite vom Hörer bemerkt und als Qualitätsmangel wahrgenommen werden. Audiosignalverarbeitungselemente in der Signalkette, wie z. B. eine automatische Verstärkungssteuerungs-(AGC-)Schleife oder ein Stereodekoder, können in den Übergängen versagen und zu einer weiteren Qualitätsverschlechterung des wiedergegebenen Audiosignals führen.
  • Wie in der auch anstehenden US-Patentanmeldung mit der Seriennummer US 09/382,727 beschrieben, kann eine allmählichere Änderung der Bandbreite (und damit eine Reduzierung der Schalteinflüsse) durch Verändern der Filterkoeffizienten in einer Reihe von Schritten anstelle des abrupten, direkten Schaltens zwischen Breit- und Schmalband realisiert werden. Die Berechnungs- oder Speicheranforderungen zur Erzeugung einer Anzahl von verschiedenen Filtern mit den gewünschten Charakteristiken sind jedoch unerwünscht groß.
  • Zusammenfassung der Erfindung
  • Die Erfindung hat den Vorteil der Erhaltung einer gewünschten Filtercharakteristik während des Veränderns, wobei Filterkoeffizienten auf rechentechnisch wirkungsvolle Art erhalten werden.
  • In einer Ausgestaltung der Erfindung wird ein Verfahren zum Herleiten eines Polynoms für die Erzeugung von Filterkoeffizienten bereitgestellt, die in einem Digitalfilter einer elektronischen Einheit mit einem Digitalsignalprozessor verwendet werden. Das Digitalfilter besitzt zum Erreichen einer gewünschten, zwischen einer ersten Charakteristik mit einer ersten Grenzfrequenz und einer zweiten Charakteristik mit einer zweiten Grenzfrequenz liegenden Bandbreitencharakteristik eine Vielzahl von konfigurierbaren Filtersektionen. Zum Bereitstellen der ersten Charakteristik wird ein erster Satz Filterkoeffizienten bestimmt. Zum Bereitstellen der zweiten Charakteristik wird ein zweiter Satz Filterkoeffizienten bestimmt. Zum Bereitstellen einer dritten Charakteristik mit einer dritten, zwischen erster und zweiter Grenzfrequenz liegenden Grenzfrequenz wird ein dritter Satz Filterkoeffizienten bestimmt. Es wird ein Polynom bestimmt, das größtenteils für den ersten, zweiten und dritten Satz von Filterkoeffizienten passende Lösungen in Bezug auf einen von einem ersten Indexwert zu einem zweiten Indexwert, entsprechend der ersten bzw. zweiten Grenzfrequenz, abweichenden Index besitzt.
  • Kurzbeschreibung der Zeichnungen
  • 1 zeigt Kurven der Empfangsfeldstärke in einem lokalen Empfangsgebiet, in dem der gewünschte, interessierende Radiokanal durch einen benachbarten Kanal gestört ist.
  • 2 ist ein Blockdiagramm, das Teile eines DSP-Radioempfängers zeigt.
  • 3 ist ein Blockdiagramm, das die Digitalsignalverarbeitung einer in der Erfindung verwendeten Zwischenfrequenz darstellt.
  • 4 ist ein Blockdiagramm, das die innere Struktur eines Kanalfilters mit einer Kaskade von Filtersektionen zweiter Ordnung zeigt.
  • Die 5a und 5b sind Kurven, die Zwischenfiltercharakteristiken unter Verwendung des linearen Schwenken der Koeffizienten darstellen.
  • Die 6a und 6b zeigen für ideale Filterkoeffizienten passende Näherungen von Polynomen zweiter Ordnung.
  • Die 7a und 7b sind Kurven, die Zwischenfiltercharakteristiken unter Verwendung von Näherungen von Polynomen zweiter Ordnung zum Verändern der Koeffizienten darstellen.
  • 8 ist ein Flussdiagramm, das eine Vorzugsausgestaltung des Betriebs eines DSP-Radioempfängers unter Verwendung von Näherungen von Polynomen zum Erhalten von Filterkoeffizientensätzen zeigt.
  • 9 ist eine Kurve, die alternative Ausgestaltungen zum Verarbeiten eines Bandbreitenindexes darstellt, der in die Filterkoeffizientenpolynome eingegeben wird.
  • 10 ist ein Blockdiagramm, das einen DSP-Radioempfänger entsprechend einer Ausgestaltung der Erfindung darstellt.
  • 11 ist ein Flussdiagramm, das eine Vorzugsausgestaltung zum Herleiten von in einem Gerät zu speichernden Polynomen unter Verwendung der Erfindung darstellt.
  • Ausführliche Beschreibung der Vorzugsausgestaltungen
  • 1 zeigt ein Frequenzspektrum 10 einer gewünschten Rundfunksendung mit einer Zentralfrequenz 11 in einem gekennzeichneten Kanal fd zwischen einer unteren Fre quenz 12 und einer oberen Frequenz 13. Es ist ein oberer benachbarter Kanal fu dargestellt, der ein Rundfunksignal 14 mit größtenteils keinem in den gewünschten Frequenzkanal hineinreichenden Signalgehalt enthält, wodurch der obere Frequenzkanal keine Störung verursacht. Es ist jedoch ein unterer benachbarter Kanal fl einer Rundfunksendung dargestellt, der ein Frequenzspektrum 15 mit einem wesentlichen Signalgehalt oberhalb der unteren Frequenz 12 des gewünschten Kanals besitzt. Die daraus resultierende Störung verschlechtert den Empfang der gewünschten Radiosendung.
  • Die Störung eines benachbarten Kanals kann durch Schmälerung des Durchgangsbereichs eines Durchgangsfilters im Empfänger durch Reduzierung des Signalgehalts des den Empfänger passierenden benachbarten Kanals vermindert werden. Damit zeigt 1 eine Schmalbandbreite 16, die in den Signalweg der Zwischenfrequenz zwecks Abschwächung einer Störung durch einen benachbarten Kanal geschaltet werden kann. Ist keine Störung durch einen benachbarten Kanal vorhanden, wird zwecks Maximierung der Qualität des empfangenen gewünschten Signals eine Breitbandbreite 17 verwendet. Im Empfänger wird die Zentralfrequenz 11 in eine Zwischenfrequenz übersetzt, die eine Nullpunkt-Zwischenfrequenz sein kann. In diesem Fall ist das Filter ein Tiefpassfilter.
  • 2 ist ein Blockdiagramm, das einen Radioempfänger mit Digitalsignalverarbeitung darstellt. Eine an einen HF-Verstärker 21 gekoppelte Antenne 20 empfängt Rundfunk-HF-Signale. Die verstärkten HF-Signale werden an einem Eingang einer Mischeinrichtung 22 bereitgestellt. Ein lokaler Oszillator 23 stellt ein Mischsignal an einem zweiten Eingang der Mischeinrichtung 22 bereit, wobei das Mischsignal eine Frequenz unter Kontrolle einer Abstimmsteuerschaltung (nicht dargestellt) hat. Ein Trägersignal in Form eines Zwischenfrequenz-(IF-)Signals wird von der Mischeinrichtung 22 an einen Eingang eines Analog-zu-Digital-(A/D-)Wandlers 24 geliefert. Ein digitalisiertes IF-Signal wird zum Filtern, Demodulieren und zur weiteren Verarbeitung des resultierenden Audiosignals an einen Digitalsignalprozessor-(DSP-)Block 25 geliefert. Der Block 25 enthält Datenspeicher und Programmspeicher zur Durchführung dieser Funktionen. Ein Endaudiosignal wird vom DSP 25 an den Eingang des Digital-zu-Analog-(D/A-)Wandlers 26 geliefert, der ein Lautsprechersystem 27 mit analogen Audiosignalen versorgt.
  • Die Verarbeitung des digitalisierten IF-Signals im DSP 25 ist in 3 detaillierter dargestellt. Die der 3 entsprechenden Ausgestaltung ist besonders zum Verarbeiten von AM-Signalen ausgelegt. Das digitalisierte IF-Signal wird zur Erzeugung von Gleichtakt-(I-) und Quadraturphasen-(Q-)Signalen an den Eingang einer komplexen Mischeinrichtung 30 geliefert. Ein Oszillator 31 erzeugt ein Zuführungssignal fif, das nominell gleich der Zwischenfrequenz des IF-Signals ist, so dass das IF-Signal zu einer neuen IF-Frequenz mit ungefähr null Hertz gemischt wird. Das Zuführungssignal wird direkt an einen Eingang einer ersten Mischeinrichtung 32 und durch einen 90°-Phasenverschiebungsblock 33 an einen Eingang einer zweiten Mischeinrichtung 34 geleitet. Das digitalisierte IF-Signal wird zum Erzeugen der I- und Q-Signale an den entsprechenden Eingängen der Mischeinrichtungen 32 und 34 bereitgestellt. Die I- und Q-Signale werden in den Dezimierblöcken 35 bzw. 36 dezimiert, um als abtastfrequenzreduzierte Signale an den Eingängen der Kanalfilter 37 und 38 bereitgestellt zu werden.
  • In der Erfindung können andere Nichtnullpunkt-IF-Frequenzen oder nichtkomplexe Signaldarstellungen verwendet werden. Eine Nullpunkt-IF-Komplexdarstellung hat in der Digitalsignalverarbeitung viele Vorteile, wie z. B. kompakte Kodegröße, minimierte Chipfläche und effiziente Datenbehandlung.
  • Die Kanalfilter 37 und 38 können mit einem Koeffizientensatz aus einer Vielzahl von Koeffizientensätzen (d. h. einem Koeffizientensatz #1 oder einem Koeffizientensatz #2) über einen Multiplexer 40 unter Kontrolle eines Nachbarkanalerfassungsblocks 41 geladen werden. Der Koeffizientensatz #1 liefert zum Beispiel eine Breitbandbreite, während der Koeffizientensatz #2 eine Schmalbandbreite liefert. Beim Verändern zwischen einer Breitbandbreite und einer Schmalbandbreite werden Zwischenkoeffizientensätze geladen. In Abhängigkeit des Vorhandenseins benachbarter Kanalstörungen wird ein zu jeder Zeit zu verwendender Koeffizientensatz ausgewählt. Unter Verwendung einer Zwischenfrequenz von null Hertz werden die Kanalfilter 37 und 38 als Tiefpassfilter im schmaleren Filter implementiert, wobei das schmalere Filter eine obere Grenzfrequenz besitzt, die niedriger als die des breiteren Filters ist. Das Vorhandensein benachbarter Kanalstörungen kann unter Verwendung jedes üblichen im Fachgebiet bekannten Verfahrens erfasst werden. Die gefilterten Outputs der Kanalfilter 37 und 38 werden an einen Signaldetektor 42 zum Erzeugen eines Audiooutputs geliefert, das zum Beispiel linke und rechte Stereosignale enthalten kann.
  • Die Speicherblöcke für die Koeffizientensätze #1 und #2 können außerdem Koeffizientenwerte enthalten, die beim Verändern einer aus den Filtersektionen ausgewählten Filtersektion zu verwenden sind. Alternativ können Informationen gespeichert werden, die bei Bedarf die Berechnung dieser Veränderungskoeffizientenwerte ermöglichen.
  • Die Kanalfilter können unter Verwendung verschiedener Filterstrukturen und – typen implementiert werden. Hier wird ein infinites Transversalfilter (IIR), wie es in 4 gezeigt ist, als Beispiel beschrieben, da seine Verwendung wegen der Vorteile seiner kompakten Größe wünschenswert ist.
  • 4 zeigt eine typische Architektur eines IIR-Filters, das aus drei zwischen einem Eingang x(n) und einem Ausgang y(n) in Kaskade geschalteten Sektionen zweiter Ordnung besteht. Ein Filter kann einen an einem Eingang eines Multiplizierers 45 eingebrachten Gesamtverstärkungsterm G0 enthalten, wobei der Multiplizierer an seinem anderen Eingang ein Input x(n) empfängt. Alternativ kann, wie im Fachgebiet bekannt, der Verstärkungsterm auf die einzelnen Sektionen verteilt werden.
  • Die Sektionen zweiter Ordnung 46, 47 und 48 sind zum Erzeugen eines Filters sechster Ordnung in Reihe geschaltet. Die erste Sektion 46 enthält einen Multiplizierer 50 zur Multiplikation des Inputs zur Sektion 46 mit einem Koeffizienten b0,1. Das resultierende Produkt wird an einen Eingang eines Summierers 51 geliefert. Der Ausgang des Summierers 51 bildet den Ausgangsknoten der Sektion 46 und ist außerdem ein innerer Knoten des gesamten Filters.
  • Das Input zu Sektion 46 wird durch eine Abtastperiode in einem Verzögerungsblock 52 und anschließender Eingabe in einen Multiplizierer 53 verzögert. Der Koeffizient b1,1 wird in einem zweiten Eingang des Multiplizierers 53 eingebracht und das Output an den Summierer 51 geliefert. Das Output des Verzögerungsblocks 52 wird nach einer weiteren Verzögerung im Verzögerungsblock 54 zunächst in einen Multiplizierer 55 zur Multiplikation mit einem Koeffizienten b2,1 geleitet. Das Output des Multiplizierers 55 wird an einen noch anderen Eingang des Summierers 51 geliefert. Die b-Koeffizienten stellen die Mitkopplungsterms von Sektion 46 bereit. Außerdem enthält Sektion 46 Rückkopplungsterms, wobei das Output des Summierers 51 in einem Verzögerungsblock 56 verzögert wird. Das verzögerte Output wird an einen Multiplizierer 57 übertragen, der außerdem einen Koeffizienten a1,1 empfängt. Das Output des Multiplizierers 57 wird an einen anderen Eingang des Summierers 51 übertragen. Das verzögerte Output vom Verzögerungsblock 56 wird durch einen Verzögerungsblock 58 und anschließend an einen Eingang eines Multiplizierers 59 geleitet. Koeffizient a2,1 wird an einen anderen Eingang des Multiplizierers 59 geliefert und das resultierende Produkt zum Summierer 51 gesendet. Die a-Koeffizienten stellen die Rückkopplungsterms von Sektion 46 bereit. Der den inneren Knoten von Sektion 46 umfassende Ausgang vom Summierer 51 ist mit dem Eingang des zweiten Bereichs 47 gekoppelt. Die Sektionen 46 und 47 sind überlappt dargestellt, da die Verzögerungsblöcke 56 und 57 von beiden Sektionen zwecks Minimieren der Hardwareanforderungen gemeinsam genutzt werden.
  • Das Input von Sektion 47 (von Summierer 51) wird in einen Multiplizierer 60 eingebracht, der außerdem Koeffizient b0,2 empfängt. Zusätzlich werden die b-Koeffizienten b1,2 und b2,2 in den Multiplizierer 62 bzw. 63 eingebracht, und die resultierenden Produkte werden im Summierer 61 addiert. Verzögerungsblock 64, Multiplizierer 65, Verzögerungsblock 66 und Multiplizierer 67 stellen wie in der vorherigen Sektion Rückkopplungsterms unter Verwendung der Koeffizienten a1,2 und a2,2 bereit. Sektion 48 arbeitet in derselben Weise mit den für die dritte Sektion verwendeten b-Koeffizienten b0,3, b1,3 und b2,3 und a-Koeffizienten a1,3 und a2,3. Zum Bereitstellen einer Endfilterung höherer Ordnung können nach Sektion 48 weitere Sektionen zweiter Ordnung in Kaskade geschaltet werden.
  • Zur Bestimmung entsprechender Werte für die a- und b-Koeffizienten für eine gewünschte Bandbreitencharakteristik sind verschiedene Verfahren unter Verwendung unterschiedlicher Typen von Filtern (z. B. Butterworth, Chebyshev oder Elliptic) allgemein bekannt. In der DSP-Hardware implementierte Filtersektionen können unterschiedliche ausgewählte Bandbreitencharakteristiken durch Laden unterschiedlicher a- und b-Koeffizienten bereitstellen, Zum Beispiel sind erster und zweiter Koeffizientensatz für das Bereitstellen einer ersten Bandbreitencharakteristik bzw. einer zweiten Bandbreitencharakteristik zwecks Schalten zwischen den Bandbreiten bei Verwendung derselben wie in der 3 gezeigten Filterhardware und -struktur einzeln hergeleitet worden (im Voraus durch die Filterkonstrukteure). Durch Austauschkoeffizientenwerte zwischen den beiden Koeffizientensätzen kann während des Betriebs des Radioempfängers entweder das Breitband- oder Schmalbandkanalfilter ausgewählt werden. Wie hier verwendet, bezieht sich ein Koeffizientensatz auf alle a- und b-Koeffizienten für alle Filtersektionen zum Bereitstellen einer Bandbreitencharakteristik.
  • Bei Verwendung der mit einem konventionellen Prozess erhaltenen Koeffizientensätze können wesentliche hörbare Wirkungen in einem Audiooutput auftreten, wenn die Filterkoeffizienten gewechselt werden. Die ebenfalls anstehende Patentanmeldung mit der Seriennummer US 09/382,727 lehrt zwecks Verringerung der hörbaren Änderung im Frequenzgang und anderer Auswirkungen eine allmähliche Änderung in der Gesamtfiltercharakteristik.
  • Würde die Bandbreite bei allen gewünschten Zwischenbandbreiten zwischen der breiten und schmalen Einstellung in diskreten Schritten unter Verwendung konventioneller Koeffizientenentwurfsverfahren geschwenkt werden, wären ein hoher Speicherbedarf und/oder umfangreiche rechentechnische Ressourcen erforderlich. Wenn 50 Zwischenschritte gebraucht würden, müssten 50 verschiedene komplette Koeffizientensätze gespeichert werden, was hinsichtlich Speicherbelegung und zusätzlicher Prozessorzeit für den Austausch der gespeicherten Koeffizienten innerhalb und außerhalb des Prozessors teuer wäre. Die ebenfalls anstehende Patentanmeldung mit der Seriennummer US 09/382,727 verändert Koeffizienten einer einzelnen Filtersektion, wobei die veränderten Koeffizientenwerte zwischen ihren An fangs- und Endwerten linear abgestuft werden. 5a zeigt eine durch lineares Verändern von Filterkoeffizienten erreichte Sperrbereichsdämpfung für eine Filtersektion, wie in der anstehenden Patentanmeldung beschrieben. Eine Kurve 80 zeigt die Charakteristik eines Schmalbandfilters und Kurve 81 die eines Breitbandfilters. Die Zwischenkurven 82 zeigen die durch lineares Verändern erhaltenen Charakteristiken von Zwischenfiltern, die im Wesentlichen aus den idealen Formen der Kurven 80 und 81 hergeleitet worden sind. 5b zeigt eine Durchgangsbereichsdämpfung für dasselbe wie in 5a gezeigte Verändern. Es ist ersichtlich, dass nach dem Verändern des Filters durch die interpolierten Koeffizientenwerte und anschließendem Laden der für die Endfilter (d. h. schmal oder breit) hergeleiteten Idealwerte in die Filtersektion eine abrupte Änderung der Filtercharakteristik bewirkt wird.
  • Die Erfindung stellt während des Bandbreitenveränderns eine idealere Filtercharakteristik bereit, wodurch mit nichtidealen Zwischencharakteristiken verbundene Nebenwirkungen vermieden werden. Außerdem stellt die Erfindung Zwischenbandbreitencharakteristiken bereit, die Leistungsspezifikationen erfüllen, durch die für Langzeitfilterung (d. h. nicht nur für kurze Perioden während des Veränderns) verwendete Bandbreitencharakteristiken erreicht werden können. Damit kann die Erfindung verwendet werden, um ein Digitalfilter mit variabler Bandbreite, der geringe rechentechnische und Speicheranforderungen stellt, aufzubauen.
  • Die Erfindung bildet die Varianz einzelner Koeffizienten einer oder mehrerer Sektionen eines Mehrsektionendigitalfilters mittels Polynomen zweiter Ordnung, wie in den 6a und 6b gezeigt werden. Die Kurve 83 repräsentiert eine idealisierte Variabilität von Werten der a1- und a2-Terme in einer bestimmten Sektion eines Tiefpassfilters über einen Bereich von Grenzfrequenzen des Filters. Während nur einige wenige dieser kontinuierlichen, idealen Filter tatsächlich entworfen werden, werden durch die Erfindung die resultierenden Koeffizientenwerte für jeden Koeffizienten an ein Polynom zweiter Ordnung angepasst. Zum Beispiel wird ein Polynom in 6a als a1 = –0,012266 + 2,317i + (–0,48752i2) gezeigt, wobei i ein Index entsprechend der Grenzfrequenz des gewünschten Filters (d. h. anteilige Abtastfrequenz des Digitalsystems) ist. 6b führt zu einer Polynomnäherung a2 = 0,67324 + 0,0035137i + 0,21376i2. Als Reaktion auf jede gewünschte Bandbreite (d. h. Grenzfrequenz) kann der Index i in entsprechende Polynome zur Erzeugung einer Polynomnäherung (wie durch die Punkte 84 gezeigt) der a-Terme unter Verwendung sehr geringer rechentechnischer und Speicherressourcen eingegeben werden. Wie zuvor beschrieben, können die b-Terme jeder Filtersektion als Reaktion auf die Werte der a-Terme berechnet werden. Die daraus resultierende Verbesserung in den Charakteristiken der Zwischenfilter wird in den 7a und 7b gezeigt, in denen die Kurven des Sperrbereichs bzw. des Durchlassbereichs dargestellt sind.
  • Anhand eines Beispiels eines AM-Kanalfilters achter Ordnung in Form eines als Direktformsektionen zweiter Ordnung in Vierfach-Kaskade implementierten IIR-Filters wird die Übertragungsfunktion für jede Sektion bestimmt als
    Figure 00090001
    wobei ein Sektionsindex X im Bereich von 0 bis 3 liegt und die vier Filtersektionen zweiter Ordnung repräsentiert. Die Filterkoeffizienten werden in der Laufzeit mithilfe eines Satzes von Polynomen und einer- Verstärkungsberechnung berechnet, die eine DC-Verstärkung Eins sicherstellt. Die a1- und a2-Terme jeder einzelnen Sektion werden aus einem entsprechenden Satz Polynome zweiter Ordnung basierend auf dem Bandbreitenpositionsindex i wie folgt berechnet: a1 = p0,1 + (p1,1·i) + (p2,1·i2), und a2 = p0,2 + (p1,2·i) + (p2,2·i2).
  • Zur Sicherstellung einer für jede Sektion passenden Verstärkung kann die Verstärkung jeder Sektion veranlasst werden, in allen Sektionen gleich einer bestimmten Konstante, wie z. B. Verstärkung Eins (d. h. eine Verstärkung von 0 dB), zu sein. Alternativ können andere Verstärkungswerte oder unterschiedliche Verstärkungswerte für verschiedene Filtersektionen verwendet werden. Verstärkung Eins sorgt für eine gute Leistung und ermöglicht eine effiziente, auf den a-Termen in einer bestimmten Filtersektion basierende Berechnung der b-Terme.
  • Figure 00090002
  • In Abhängigkeit von der Hardware können die b-Terme auch wie folgt berechnet werden: b1 = 0,25 + 0,5a1 + 0,5a2; b0 = b1 ÷ 2; und b2 = b0.
  • Der Laufzeitbetrieb eines diese Erfindung verwendenden Radioempfängers wird in 8 gezeigt. In Schritt 86 werden die Bandbreite des Kanalfilters und der entsprechende Index initialisiert (Der Empfänger kann z. B. zu Beginn das Schmalbandfilter verwenden, bis gute Empfangsbedingungen erkannt werden.). In Schritt 87 wird festgestellt, ob die Bandbreite erhöht (oder reduziert oder nicht modifiziert) werden sollte. Hat eine Erhöhung zu erfolgen, wird der Bandbreitenpositionsindex i in Schritt 88 um eine Vergrößerungsschrittweite erhöht. Hat eine Reduzierung zu erfolgen, wird der Bandbreitenpositionsindex i in Schritt 89 um eine Verkleinerungsschrittweite herabgesetzt. In einer Vorzugsausgestaltung variiert der Index i zwischen 0 und 1, wobei 0 das schmalste Filter (d. h. die niedrigste Grenzfrequenz) und 1 das breiteste Filter repräsentiert. Die Filterauflösung (d. h. die Anzahl der Zwischenfilter, die zwischen den Endpunkten des breiten und des schmalen Filters erreicht werden kann) kann gleich der numerischen Auflösung des DSP sein, obwohl bei Anwendung des Veränderns eine viel niedrigere Auflösung verwendet werden würde. Die Vergrößerungsschrittweite muss nicht gleich der Verkleinerungsschrittweite sein, da ein schnelleres Verändern bei Schmälerung der Bandbreite vorteilhaft sein kann. Wenn das Erhöhen oder Herabsetzen des Indexes zu einem Index außerhalb des Bereichs führt, wird der Index auf 0 bzw. 1 gesetzt.
  • Die Schrittweiten während einer Veränderung in eine Richtung müssen nicht alle gleich sein. Es kann eine Änderung der Progressionsgeschwindigkeit des Veränderns erwünscht sein. Deshalb kann ein optionaler Schritt 90 verwendet werden, um den Anfangsindex bis zu einem gewandelten Index zwecks Bereitstellung einer nichtlinearen Teilung der Zwischenfilter über einen Bereich von beim Verändern verwendeter Grenzfrequenzen abzubilden. 9 zeigt eine Kurve 96 eines nicht gewandelten Index, der die Bandbreite oder Grenzfrequenz linear variiert. Die Kurven 97 und 98 stellen Beispiele einer nichtlinearen Beziehung zwischen dem Index und der Zwischenfrequenzbandweite dar. Diese nichtlinearen Beziehungen können zum Beispiel durch Verwendung entsprechender Polynome zum Wandeln des Indexes erhalten werden. Damit kann der Anfangsindexwert schrittweise in gleichen Schrittweiten geändert werden, wird aber anschließend unter Verwendung eines Polynoms gewandelt, so dass eine nichtlineare Änderung der Schrittweite erhalten wird, wie durch die Kurven 97 oder 98 gezeigt. Diese Wandlung kann außerdem während des Betriebs des Emp fängers in Abhängigkeit verschiedener Faktoren, wie z. B. der Schwierigkeit einer Störungsbedingung oder der Schwäche eines empfangenen Signals, modifiziert werden.
  • Auf der Basis des Endindexes (gewandelt oder nicht gewandelt) wird, beginnend bei Schritt 91, ein Zwischenbandfilter erhalten, wobei für jede Filtersektion Filterkoeffizienten berechnet werden. Die beste Leistung wird durch Veränderungskoeffizienten in jeder Sektion erreicht, wobei die Erfindung jedoch ein Verändern nicht aller Sektionen umfasst. Für die momentane Sektion werden die a1- und a2-Terme unter Verwendung der zuvor in Schritt 92 hergeleiteten und gespeicherten entsprechenden Polynome durch Polynome zweiter Ordnung angenähert. In Schritt 93 werden die b0-, b1- und b2-Terme für die momentane Sektion aus den Werten der a1- und a2-Terme berechnet. Es erfolgt ein Rücksprung auf Schritt 91 und die Schritte 92 und 93 werden wiederholt, bis ein Filterkoeffizient bestimmt worden ist.
  • Wenn alle Koeffizientenwerte bestimmt sind, werden sie in Schritt 94 in die Filtersektionen geladen, und in Schritt 95 wird eine Kanalfilterung unter Verwendung der neuen Zwischenfilter durchgeführt. Es erfolgt ein Rücksprung auf Schritt 87, um die Richtung zu bestimmen, die eventuell für weiteres erforderlich ist.
  • Wenn gewünscht, können weitere Schritte durchgeführt werden, um zu bestimmen, warm der Endpunkt (d. h. maximale oder minimale Grenzfrequenz) erreicht ist, und in die Filtersektionen kann anstelle einer Polynomnäherung ein separat gespeicherter Satz von idealen Filterkoeffizienten geladen werden. Umgekehrt können das untere und obere Bandbreitenextrem eines bestimmten Filters durch Aufzwingen minimaler und maximaler Werte auf den Bandbreitenindex gesteuert werden, so dass ein Verändern nur über einen Teil eines gesamten Bereichs durchgeführt wird, für den das Filter und die Polynome ursprünglich ausgelegt wurden (wobei die Polynomnäherungen immer angewendet werden).
  • In 10 ist ein Apparat zur Durchführung des Verfahrens von 8 dargestellt. An einem Qualitätsdetektor 100 und an einer Quadraturmischeinrichtung 107 wird ein Digitalsignal (z. B. ein AM-Rundfunkübertragungssignal) bereitgestellt. Der Qualitätsdetektor 100 kann ein Störungsdetektor (z. B. zum Bestimmen der Starke unerwünschter Signale bei benachbarten Kanalfrequenzen, so dass durch Schmälerung der Kanalfilterbandbreite die störenden Signale abgedeckt werden können) oder ein Signalstärkedetektor sein (z. B. zum Feststellen, ob das gewünschte Signal schwach ist, wobei das Gesamtsignal im Verhältnis zum Rauschen durch Schmälerung des Kanalfilters verbessert werden kann). Detektor 100 stellt ein eine niedrige Qualität anzeigendes Signal an einem Zähler 101 bereit, der während des Signals niedriger Qualität mit einer durch eine Uhr 102 bestimmten Zählfrequenz vorwärts zählt. Nach einer festgelegten Anzahl von Zählimpulsen wird eine Gesamtzahl im Zähler 101 mit einem Schwellenwert in einem Komparator 103 verglichen. Damit müssen die Bedingungen niedriger Qualität (Störung oder schwaches Signal) für mindestens einen festgelegten Zeitmittelwert erfasst werden, bevor ein Umschalten in das Schmalbandfilter durch den Komparator 103 angezeigt wird.
  • Eine Steuerung 104 empfängt ein Output vom Komparator 103 und erneuert den am Speicherort 105 im Digitalsignalprozessor gespeicherten Bandbreitenindex i. Die auch als ROM-Daten im Digitalsignalprozessor speicherbaren Polynome zum Berechnen der Filterkoeffizientenwerte sind in einem Speicher 106 gespeichert. Die Steuerung 104 bestimmt die Filterkoeffizientenwerte entsprechend dem momentanen Wert des Index i und lädt sie für die I- und Q-Quadratursignale in die Kanalfilter 108 bzw. 109. Anschließend werden die Quadratursignale durch einen Detektor 110 erfasst (d. h. demoduliert).
  • 11 zeigt einen Filterentwurfsprozess detaillierter. In Schritt 120 können konventionelle Digitalfilterentwurfshilfsmittel zum Berechnen von Filterkoeffizienten für mindestens das Schmalband-, Breitband- oder eines der Zwischenbandfilter verwendet werden. Vorzugsweise hat das Zwischenbandfilter eine Grenzfrequenz in der Mitte zwischen den Grenzfrequenzen des Schmal- und Breitbandfilters. Obwohl hier allgemein Tiefpassfilter dargestellt sind, könnten die Filter alternativ aus Hochpassfiltern bestehen.
  • Obwohl drei ideale Koeffizientensätze zum Bestimmen eines Polynoms ausreichen können, wird eine höhere Genauigkeit durch das Entwerfen weiterer für das Bestimmen der Polynomkoeffizienten zu verwendenden Zwischenfilter erreicht. Zum Beispiel sind in einem Tiefpassfilter mit einer Grenzfrequenz von –3 dB für das Schmalbandfilter von etwa 3 kHz und einer Grenzfrequenz von –3 dB für das Breitbandfilter von etwa 7 kHz bis zu etwa 20 ideal entworfene Filter mit Grenzfrequenzen in größtenteils gleichmäßiger Teilung von 3 kHz bis 7 kHz verwendet worden. Alternativ müssen die Teilungen der gewünschten Filter nicht linear sein. Eine engere Teilung kann in Teilen des Bereichs bereitgestellt werden, wo eine höhere Genauigkeit der Filterkoeffizientensätze erwünscht ist.
  • In Schritt 121 werden die Filtersektionen der einzeln manuell entworfenen Filter zum Bereitstellen des bestmöglichen Zusammenpassens von Q-Faktoren entsprechend gruppiert. Das bedeutet ein unabhängiges Neuordnen von Sektionen der Filter. In Schritt 122 werden die zusammenpassenden Sektionen zur Maximierung des dynamischen Bereichs und der Rauschleistung erneut zueinander geordnet, wie es den mit dem Fachgebiet vertrauten Personen bekannt ist. In Schritt 123 werden Polynome den hergeleiteten Koeffizientenwertsätzen der entsprechenden Filtersektionen angepasst. Alternativ könnte die Polynoman passung selbstverständlich vor der Gruppierung und der Neuordnung durchgeführt werden. In Schritt 124 werden die Polynomkoeffizientensätze im DSP-Speicher gespeichert.
  • Zur weiteren Darstellung des Entwurfsprozesses wird ein detaillierteres Beispiel gezeigt. Das Kanalfilter besteht aus vier in Kaskade geschalteten Sektionen zweiter Ordnung (Biquadrate). Die Übertragungsfunktion jedes Biquadrats ist nachfolgend dargestellt.
  • Figure 00130001
  • Die Übertragungsfunktion zeigt, dass zum Verändern der Bandbreite der Filtersektion ein Modifizieren der Stelle von Pol (p1) verwendet werden kann. Das erklärt, warum Polynomnäherungen nur für die a1- und a2-Terme notwendig sind. Die zuvor gezeigte Tiefpassfilterstruktur vorausgesetzt, kann eine Polynomnäherung zweiter Ordnung zum Darstellen von a1 und a2 verwendet werden, wenn die –3 dB-Grenzfrequenz geringer als ein Viertel der Abtastfrequenz ist. Die b-Terme können dann erneut aus den neuen a-Termen berechnet werden, damit die DC-Verstärkung Eins bleibt. Ein Satz Polynome wird zum Berechnen beider a-Terme in allen Stufen (z. B. 8 Polynome bei Verwendung von 4 in Kaskade geschalteten IIR-Filtersektionen zweiter Ordnung) auf der Basis eines Positionsindexes einer Kanalfilterbandbreite wie folgt hergeleitet: a1,X = p0,1,X + (p1,1,X·i) + (P2,1,X·i2), und a2,X = P0,2,X + (p1,2,X·i) + (p2,2,X·i2),wobei p Polynomkoeffizienten sind und i der Bandbreitenindex ist. Da sich die Nullpunkte von HLP auf der Nyquist-Abtastfrequenz befinden, werden die b0-, b1- und b2-Koeffizienten jeder Biquadratsektion aus den a1- und a2-Koeffizienten effizient berechnet, so dass die DC-Verstärkung jeder Sektion Eins ist. Da bei DC z = e = ej0 = 1, gilt:
    Figure 00130002
  • Unter Verwendung konventioneller Verfahren leitet ein Konstrukteur drei verschiedene Filter, und zwar das Schmalband-, das Mittelband- und das Breitbandfilter her. Vorzugsweise besitzt das Mittelband eine Bandbreite dicht am Mittelwert der Schmalband- und Breitbandbreiten. Vorzugsweise ist jedes Filter ein IIR-Filter achter Ordnung mit der Einschränkung, dass sich alle Nullpunkte auf der Nyquist-Frequenz befinden. Dieses Kriterium wird vom Entwurf jedes Filters als Butterworth-Sektion oder Zusammenstellung von Butterworth-Sektionen erfüllt.
  • Nachdem die konventionellen (d. h. idealen) Filter entworfen worden sind, werden ihre Übertragungsfunktionen in Sektionen zweiter Ordnung (Biquadrate) überführt. Vorzugsweise kann jede Sektion die DC-Verstärkung Eins und einen Nullpunkt auf der Nyquist-Frequenz haben.
  • Die Biquadrat-Sektionen der Schmalband-, Mittelband- und Breitbandfilter werden auf der Basis der allgemeinen Form der Größenwiedergabe jeder Sektion gruppiert. Zum Beispiel sollte die Sektion mit dem höchsten Q im Schmalbandfilter mit den höchsten Q-Sektionen der Mittelband- und Breitbandfilter usw. gruppiert werden. Eine direkte Lösung besteht darin, die Biquadratsektionen so zu ordnen, dass die Reihenfolge der Pole absteigend ist. Die Gruppierung der Sektionen hilft beim Umschalten zwischen den Filtern, die Sprungantwort zu minimieren.
  • Die Gruppen der Biquadratsektionen werden danach zur Berücksichtigung einer dynamischen Lautstärke und zur Minimalisierung eines übermäßigen Rauschens neu geordnet. Bei Berücksichtigung einer dynamischen Lautstärke ist es vorteilhaft abzusichern, dass maximale Signalpegel am Ausgang jeder Zwischensektion nicht überschritten werden. Übermäßiges Rauschen bezieht sich auf die Form von Abschneidungsrauschen, wie es sich durch darauffolgende kaskadengeschaltete Filtersektionen ausbreitet. Zum Minimieren des übermäßigen Rauschens wird das Ordnen der Sektionen gewählt. Ein einfaches Verfahren ist das Ordnen der Sektionen der Mittelbandfilter zur Berücksichtigung einer dynamischen Lautstärke und bei übermäßigem Rauschen und anschließendes Neuordnen der Schmalband- und Breitbandfilter, um dieselben Gruppierungen, wie zuvor hergeleitet, beizubehalten. Ein alternatives Verfahren ist der Versuch, das Rauschen basierend auf allen drei Filtern zu minimieren.
  • Zur Näherungsberechnung der a1- und a2-Koeffizienten werden Polynome zweiter Ordnung für die Modellierung der Filterkoeffizientenwerte gebildet. Es sind nur die a1- und a2-Koeffizienten erforderlich, da b0, b1 und b2 dynamisch auf der Basis von a1 und a2 berechnet werden. Es kann in einigen DSP-Implementierungen wünschenswert sein, alle a1- und a2-Koeffizienten mit einem Faktor 2 abwärts zu skalieren, um die DSP-Implementierung widerzuspiegeln. Die Polynomnäherung würde also die momentanen, durch 2 dividierten a-Koeffizienten repräsentieren.
  • Vorzugsweise werden für jeden der skalierten a1- und a2-Terme die die Koeffizientensätze jeder Filtersektion am besten anpassenden Polynome zweiter Ordnung unter Verwendung der Näherung kleinster Quadrate (d. h. eine kleinste-Quadrate-Anpassung der Koeffizienten, die das Polynom für den ersten, zweiten und dritten Filterkoeffizientensatz definiert) hergeleitet. Durch dieses Verfahren wird der Polynomkoeffizientenvektor p bestimmt, der den quadratischen Fehler ||e||2 minimiert, wie nachfolgend dargestellt:
    Figure 00150001
    wobei ind die Matrix der Indizes für das Polynom und a der Vektor der Filterkoeffizienten für das Schmalband-, Mittelband- und Breitbandfilter (oder mehr Sätze, wenn mehr ideale Filter entworfen werden) ist. Eine Lösung, die durch Einführen des Orthogonalitätsprinzips vereinfacht wird, ist explizit nachfolgend dargestellt:
    Figure 00150002
  • Da der Bandbreitenindex von 0 bis 1 reicht, sind die Werte für indschmal, indmittel und indbreit 0, 0,5 bzw. 1. Deshalb wird die vorherige Lösung zu
    Figure 00150003
  • Zum Bestimmen der a1- und a2-Koeffizientensätze für jede der Biquadratsektionen wird eine Gesamtgleichung der 8 Gleichungen der Polynomsätze erzeugt. Die resultierenden Polynomkoeffizienten werden dann in den Prozessorspeicher geschrieben.

Claims (10)

  1. Verfahren zur Signalverarbeitung mit Digitalfiltern (37, 38; 108, 109) in einer elektronischen Einheit mit einem Digitalsignalprozessor (25), wobei das Digitalfilter (37, 38; 108, 109) zum Erreichen einer zwischen einer ersten Charakteristik mit einer ersten Grenzfrequenz (81) und einer zweiten Charakteristik mit einer zweiten Grenzfrequenz (80) liegenden Bandweitencharakteristik mit einer dritten Grenzfrequenz eine Vielzahl von konfigurierbaren Filtersektionen (46, 47, 48) besitzt, wobei den Filtersektionen (46, 47, 48) Filterkoeffizienten (a1, a2, b0, b1, b2) zugeordnet sind, wobei die Filterkoeffizienten (a1, a2) mit folgenden Polynomen a1 = p0,1 + (p1,1·i) + (p2,1·i2) (I) a2 = p0,2 + (p1,2·i) + (P2,2·i2) (II)bestimmt werden, wobei i der Bandweitenpositionsindex als Variable in einem Speicher (105) einer Steuerung (104) des Digitalsignalprozessors (25) einem Abtastfrequenzbereich (fS) zwischen der ersten Grenzfrequenz (81, i = 1) und der zweiten Grenzfrequenz (80, i = 0) mit den Indexwerten 0 ≤ i ≤ 1 entspricht und wobei der dritten Bandweitencharakteristik mit der dritten Grenzfrequenz ein Indexwert i innerhalb 0 ≤ i ≤ 1 zugeordnet wird, und wobei die Polynomkoeffizienten (p) modifizierbare Polstellen der Übertragungsfunktion HLP(z) der Filtersektionen zweiter Ordnung mit der Übertragungsfunktion
    Figure 00160001
    und wobei z die Verzögerungen darstellen, und wobei die Berechnung (93) der Koeffizientenwerte für die b0, b1, b2-Koeffizienten in Verbindung mit den Koeffizientenwerten der a1, a2-Koeffizienten in der Steuerung (104) derart erfolgt, dass der Berechnung der b0, b1, b2-Koeffizientenwerte eine angepasste Verstärkung (G) der jeweiligen Filtersektion (46, 47, 48) zugrunde liegt, so dass ein dritter Satz von Filterkoeffizienten (a1, a2, b0, b1, b2) bestimmt wird, der der dritten Grenzfrequenz zugeordnet wird, und wobei folgende nachgeordnete Schritte durchgeführt werden: – Berechnung (120) von Filterkoeffizienten (a1, a2, b0, b1, b2) für mindestens drei ausgebildete Charakteristiken – zumindest für Schmalband-, Breitband- und Zwischenbandfilter – zum Herleiten der Polynomkoeffizienten (p0,1, p1,1, p2,1, p0,2, p1,2, p2,2) der Polynome zur Einspeicherung in einen Speicher des Digitalsignalprozessors (25), wobei die Polynomkoeffizienten (p0,1, p1,1, p2,1, p0,2, p1,2, p2,2) nach folgender Matrixgleichung aus den Polynomen (IV) ermittelt werden:
    Figure 00170001
    wobei ind die Matrix der Indizes i für das Polynom und a der Vektor der Filterkoeffizienten für das Schmalbandfilter, Zwischenbandfilter und Breitbandfilter sind, – Einschreiben der berechneten, resultierenden Polynomkoeffizienten (p0,1, p1,1, p2,1, p0,2, p1,2, p2,2) in den Speicher (106) des Digitalsignaiprozessors (25).
  2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass das Herleiten von Polynomen mittels drei vorgegebener Sätze von Filterkoeffizienten als Funktionen des Bandweitenpositionsindexes (i) der drei bereitgestellten Charakteristiken mit einer dritten, zwischen erster und zweiter Grenzfrequenz liegenden Grenzfrequenz erfolgt, wobei mit dem Polynom der erste, zweite und dritte Satz von Filterkoeffizienten in Bezug auf einen sich von einem ersten Indexwert zu einem zweiten Indexwert entsprechend der ersten bzw. der zweiten Grenzfrequenz verändernden Index angepasst werden.
  3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass bei einem Bandweitenpositionsindex (i) Werte von 0 bis 1 vorhanden sind, wobei die Werte indSchmal = 0, indZwischen = 0,5 und indBreit = 1 sind, so dass folgende Polynomkoeffizienten-Matrixgleichung eingesetzt wird:
    Figure 00170002
    wobei für jede der Filtersektionen (46, 47, 48) zum Bestimmen der a1, a2-Koeffizientensätze eine Gesamtgleichung von acht Gleichungen der Polynomsätze erzeugt werden.
  4. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass die Berechnung einer weiteren Vielzahl von Filterkoeffizientensätzen zum Bereitstellen weiterer Charakteristiken mit entsprechenden Grenzfrequenzen zwischen der ersten und der zweiten Grenzfrequenz errfolgt, wobei das Polynom Lösungen besitzt, die die weiteren Filterkoeffizientensätze an das vorgesehene Digitalfilter anpassen.
  5. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass mindestens ein entsprechendes Polynom für jede Filtersektion (46, 47, 48) hergeleitet wird.
  6. Verfahren nach den Ansprüchen 1 bis 5, dadurch gekennzeichnet, dass eine Gruppierung der Koeffizienten der jeweiligen Filtersektionen (46, 47, 48) für die ersten, zweiten und dritten Koeffizientensätze herbeigeführt wird.
  7. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass die Herleitung mindestens eines Polynoms mit den ersten, zweiten und dritten Koeffizientensätzen anpassenden Lösungen aus einer mindestens quadratischen Anpassung von Polynomkoeffizienten (p) für das Polynom an den ersten, zweiten und dritten Filterkoeffizienten besteht.
  8. Verfahren nach den Ansprüchen 1 bis 7, gekennzeichnet durch die Schritte zur Bereitstellung von Filterkoeffizienten (a1, a2, b0, b1, b2) in den Filtersektionen (46, 47, 48): – Vorgabe eines ersten Satzes Filterkoeffizienten (a1, a2, b0, b1, b2) zur Bereitstellung der ersten Charakteristik, – Vorgabe eines zweiten Satzes Filterkoeffizienten (a1, a2, b0, b1, b2) zur Bereitstellung der zweiten Charakteristik, – Vorgabe eines dritten Satzes Filterkoeffizienten (a1, a2, b0, b1, b2) zur Bereitstellung der dritten Charakteristik mit einer dritten, zwischen erster und zweiter Grenzfrequenz liegenden Grenzfrequenz, – Herleiten eines Polynoms
    Figure 00190001
    das für den ersten, zweiten und dritten Satz von Filterkoeffizienten (a1, a2, b0, b1, b2) passende Polynomkoeffizienten (p) in Bezug auf einen von einem ersten Indexwert indBreit zu einem zweiten Indexwert indSchmal entsprechend der ersten bzw. zweiten Grenzfrequenz abweichenden Index indZwichen besitzt, – Speichern des ersten und zweiten Filterkoeffizientensatzes und des Polynoms in der elektronischen Einheit, – Berechnung einer gewünschten Bandbreitencharakteristik während des Betriebs der elektronischen Einheit, – Berechnung eines Bandweitenpositionsindexwertes (i) entsprechend der gewünschten Bandbreitencharakteristik, – Berechnen eines Satzes von Filterkoeffizienten (a1, a2, b0, b1, b2) als Reaktion auf den Indexwert indZwichen und auf das gespeicherte Polynom, – Laden der berechneten Filterkoeffizienten (a1, a2, b0, b1, b2) in die Filtersektionen (46, 47, 48).
  9. Einrichtung zur Durchführung des Verfahrens in einem Digitalsignalprozessor (25) mit einem Digitalfilter (37, 38; 108, 109) nach den Verfahrens-Ansprüchen 1 bis 8, dadurch gekennzeichnet, dass ein bereitgestelltes Zwischenfrequenz-Digitalsignal (IF) in einer Leitungsteilung sowohl an einen Qualitätsdetektor (100) als auch an eine Quadraturmischeinrichtung (107) geführt ist, wobei der Qualitätsdetektor (100) mit einem Zähler (101) verbunden ist, der eingangsseitig mit einer bestimmten Zählfrequenz vorwärts zählenden Uhr (102) verbunden ist und ausgangsseitig mit einem Komparator (103), an den zum Vergleich ein Schwellenwert geführt ist, in Verbindung steht, wobei der Komparator (103) ausgangsseitig mit einer Steuerung (104) verbunden ist, die in einem internen Speicherort (105) den Bandweitenpositionsindex (i) spei chert und die eingangsseitig mit einem externen Speicher (106) zur Speicherung der Polynome a1 = p0,1 + (p1,1·i) + (p2,1·i2) (I) a2 = p0,2 + (p1,2·i) + (P2,2·i2) (II)zum Berechnen der Filterkoeffizientenwerte (a1, a2, b0, b1, b2) und zur Berechnung der Polynomkoeffizienten (p0,1, p1,1, p2,1, p0,2, p1,2, p2,2) zur Ausbildung von Sätzen (#1, #2) von Filterkoeffizienten verbunden ist, wobei die Steuerung (104) die Filterkoeffizientenwerte (a1, a2, b0, b1, b2) entsprechend dem momentanen Wert des Bandweitenpositionsindexes i mit 0 ≤ i ≤ 1 für die gewünschte Bandweitencharakteristik mit der gewünschten Grenzfrequenz bestimmt, wobei dem Indexwert i = 1 die erste Bandweitencharakteristik mit der ersten Grenzfrequenz (81) und dem Indexwert i = 0 die zweite Bandweitencharakteristik mit der zweiten Grenzfrequenz (80) zugeordnet sind und wobei ein Wert des Bandweitenpositionsindexes (i) jeweils einer Abtastfrequenz der zwischen den Grenzfrequenzen (81, 80) befindlichen Abtastfrequenzen (fS) zugeordnet ist, wobei die Steuerung (104) ausgangsseitig mit einem Digitalfilter aus mindestens zwei Kanalfiltern (37, 38; 108, 109) zur Eingabe der Sätze (#1, #2) von Filterkoeffizienten (a1, a2, b0, b1, b2) verbunden ist, wobei einer der Kanalfilter (108) mit einer I-Signalleitung (I) und der andere Kanalfilter (109) mit einer Q-Signalleitung (Q) jeweils eingangsseitig mit der Quadraturmischeinrichtung (107) in Verbindung steht, und die Ausgangsleitungen (R, L) der Kanalfilter (37, 38; 108, 109) zu einem Detektor (42; 110) zur Ausbildung von Audiosignalen führen.
  10. Einrichtung nach Anspruch 9, dadurch gekennzeichnet, dass der Qualitätsdetektor (100) entweder ein Störungsdetektor oder ein Signalstärkedetektor ist.
DE10317698A 2002-04-15 2003-04-11 Verfahren zum Entwerfen von Polynomen zur Steuerung des Veränderns von anpassungsfähigen Digitalfiltern Expired - Fee Related DE10317698B4 (de)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US10/122,776 2002-04-15
US10/122,776 US6988116B2 (en) 2002-04-15 2002-04-15 Method of designing polynomials for controlling the slewing of adaptive digital films
US10/122776 2002-04-15

Publications (2)

Publication Number Publication Date
DE10317698A1 DE10317698A1 (de) 2004-05-19
DE10317698B4 true DE10317698B4 (de) 2010-10-14

Family

ID=22404702

Family Applications (1)

Application Number Title Priority Date Filing Date
DE10317698A Expired - Fee Related DE10317698B4 (de) 2002-04-15 2003-04-11 Verfahren zum Entwerfen von Polynomen zur Steuerung des Veränderns von anpassungsfähigen Digitalfiltern

Country Status (3)

Country Link
US (1) US6988116B2 (de)
DE (1) DE10317698B4 (de)
GB (1) GB2391122B (de)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4074823B2 (ja) * 2003-03-05 2008-04-16 株式会社デンソー A/d変換出力データの非直線性補正方法及び非直線性補正装置
US7221975B2 (en) * 2003-12-04 2007-05-22 Maquet Critical Care Ab Signal filtering using orthogonal polynomials and removal of edge effects
US7742914B2 (en) * 2005-03-07 2010-06-22 Daniel A. Kosek Audio spectral noise reduction method and apparatus
US7797065B2 (en) * 2005-05-02 2010-09-14 Texas Instruments Incorporated Automute detection in digital audio amplifiers
US7205914B1 (en) * 2005-10-28 2007-04-17 Motorola, Inc. System and method for efficient upsampled signal processing
FI20075931L (fi) * 2007-12-18 2009-06-19 Abb Oy Menetelmä ja laite signaalidatan siirtämiseksi
FR2935850B1 (fr) * 2008-09-05 2011-04-01 Commissariat Energie Atomique Dispositif de filtrage a structure hierarchique et dispositif de filtrage reconfigurable
US20120042193A1 (en) * 2010-08-13 2012-02-16 Apple Inc. Clock frequency selection scheme
US9148328B2 (en) * 2010-10-29 2015-09-29 Samsung Electronics Co., Ltd. Digital I/Q imbalance compensation in a quadrature receiver
US8793298B2 (en) * 2010-11-01 2014-07-29 Blackberry Limited Reconfigurable digital signal filter processor
KR101156667B1 (ko) 2011-12-06 2012-06-14 주식회사 에이디알에프코리아 통신 시스템의 필터 계수 설정 방법
US8750364B2 (en) 2011-12-31 2014-06-10 St-Ericsson Sa Interpolation of filter coefficients
US8924910B2 (en) * 2012-11-09 2014-12-30 Analog Devices, Inc. Filter design tool
KR101920719B1 (ko) 2012-11-19 2019-02-13 삼성전자주식회사 논리 장치, 논리 장치를 포함하는 디지털 필터 및 논리 장치를 제어하는 방법
CN111669149B (zh) * 2020-06-21 2022-11-22 陕西航空电气有限责任公司 一种自适应巴特沃斯低通数字滤波器设计方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5502663A (en) * 1992-12-14 1996-03-26 Apple Computer, Inc. Digital filter having independent damping and frequency parameters
US5856934A (en) * 1996-08-30 1999-01-05 Icom Incorporated Digital filter system
US20020042802A1 (en) * 2000-06-07 2002-04-11 Yukihiko Mogi FIR filter and setting method of coefficient thereof
US6393450B1 (en) * 1999-08-25 2002-05-21 Visteon Global Technologies, Inc. Switched bandwidth digital filters with slewing
US6445735B1 (en) * 1999-02-08 2002-09-03 Visteon Global Technologies, Inc. Switched bandwidth digital filters with reduced transients during switching
US20020136289A1 (en) * 1999-02-08 2002-09-26 Sunil Shukla Method of slewing a digital filter providing filter sections with matched gain
US7027502B2 (en) * 2002-04-15 2006-04-11 Visteon Global Technologies, Inc. Run-time coefficient generation for digital filter with slewing bandwidth

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US102991A (en) * 1870-05-10 Improvement in wooden street pavements
US5594612A (en) * 1994-08-24 1997-01-14 Crystal Semiconductor Corporation Analog-to-digital converter with digital linearity correction
US5646569A (en) * 1995-08-30 1997-07-08 Hewlett-Packard Company Method and apparatus for AC coupling
US6154547A (en) 1998-05-07 2000-11-28 Visteon Global Technologies, Inc. Adaptive noise reduction filter with continuously variable sliding bandwidth
FR2832568A1 (fr) 2001-11-20 2003-05-23 Koninkl Philips Electronics Nv Convertisseur numerique de frequence d'echantillonnage

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5502663A (en) * 1992-12-14 1996-03-26 Apple Computer, Inc. Digital filter having independent damping and frequency parameters
US5856934A (en) * 1996-08-30 1999-01-05 Icom Incorporated Digital filter system
US6445735B1 (en) * 1999-02-08 2002-09-03 Visteon Global Technologies, Inc. Switched bandwidth digital filters with reduced transients during switching
US20020136289A1 (en) * 1999-02-08 2002-09-26 Sunil Shukla Method of slewing a digital filter providing filter sections with matched gain
US6393450B1 (en) * 1999-08-25 2002-05-21 Visteon Global Technologies, Inc. Switched bandwidth digital filters with slewing
US20020042802A1 (en) * 2000-06-07 2002-04-11 Yukihiko Mogi FIR filter and setting method of coefficient thereof
US7027502B2 (en) * 2002-04-15 2006-04-11 Visteon Global Technologies, Inc. Run-time coefficient generation for digital filter with slewing bandwidth

Also Published As

Publication number Publication date
GB2391122A (en) 2004-01-28
US6988116B2 (en) 2006-01-17
US20030195910A1 (en) 2003-10-16
GB0304722D0 (en) 2003-04-02
DE10317698A1 (de) 2004-05-19
GB2391122B (en) 2004-06-30

Similar Documents

Publication Publication Date Title
DE10317698B4 (de) Verfahren zum Entwerfen von Polynomen zur Steuerung des Veränderns von anpassungsfähigen Digitalfiltern
EP0052847B1 (de) Verfahren und Schaltungsanordnung zur Umsetzung der Abtastfrequenz einer Abtastfolge unter Umgehung der Konversion in ein kontinuierliches Signal
DE69901551T2 (de) Flexibler programmierbarer sigma-delta Analogwandler
EP0080014B1 (de) Digitaler Demodulator frequenzmodulierter Signale
DE3856269T2 (de) Digitales Entzerrungsgerät mit Möglichkeit zur getrennten Modifikation der Phasen- und Amplitudencharakteristik
DE102006047197B3 (de) Vorrichtung und Verfahren zum Verarbeiten eines reellen Subband-Signals zur Reduktion von Aliasing-Effekten
EP0758817B1 (de) Equalizer für digitalisierte Signale
DE60018624T2 (de) Digitales filter mit umschaltbarer bandbreite mit verringerten transienten während des schaltens
EP0696848B1 (de) Verfahren zur digitalen Interpolation von Signalen
DE2707936C3 (de) Einseitenband-FrequenzmultiplexÜbertragungssystem
DE10317701B4 (de) Verfahren und Digitalsignalverarbeitungseinheit zur Erzeugung von Filterkoeffizienten für Digitalfilter mit veränderlicher Bandbreite
DE2839229C2 (de) Frequenzweiche mit einem Transversalfilter
EP0135066B1 (de) Übertragungsanordnung für digitale Signale
DE60026943T2 (de) Digitales anstiegsfilter mit umschaltbarer bandbreite
DE3922469C2 (de)
US20020136289A1 (en) Method of slewing a digital filter providing filter sections with matched gain
US6954774B1 (en) Bandpass filter with ability to tune a pass frequency, passband gain and filter Q
DE10032520A1 (de) Interpolationsfilter und Verfahren zur digitalen Interpolation eines digitalen Signals
CH655213A5 (de) Sprachverarbeitungsgeraet fuer teilnehmerleitungen.
EP1440517A1 (de) Adaptives Kanalfilter für Mobilfunkempfänger und Verfahren zur Adaptiven Kanalfilterung
DE60028769T2 (de) Anti-aliasierte begrenzung bei begrenzter modulation mit stufen-funktion
DE3408384C2 (de) Schaltkreis zur Nachbildung einer Impedanz
EP1775833A1 (de) Digitales Filter und Verfahren zum Entwurf digitaler Filter mittels Integrations- und Löschfilter
EP0169961B1 (de) Digitales Filter
EP0326672A1 (de) Verfahren und Übertragungsanordung für digitale Signale

Legal Events

Date Code Title Description
OP8 Request for examination as to paragraph 44 patent law
8364 No opposition during term of opposition
R082 Change of representative

Representative=s name: MERH-IP MATIAS ERNY REICHL HOFFMANN PATENTANWA, DE

R119 Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee