Größe (Magnitude) eines Signals mit Verwendung eines 2-Ordungsfilter

Preg M. OP #6556937
Lesenswert?
• ▲
▼ 
Hallo,

Ich beabsichtige, die Funktion freqz() zu verwenden, um die Größen 
(Magnitude) bei verschiedene Grenzfrequenzen zu zeichnen. Aber zuerst 
muss ich die Grenzfrequenzen bestimmen.

Alles, was ich habe ist:

 - Ein Bandpassfilter reicht von ungefähr 100 Hz (f_p) für die apikale 
Elektrode bis 8 kHz (f_s) für die basale Elektrode. (22 Elektroden 
insgesamt)
1
b=log10(8000);

2
yc1=logspace(2,b,23);

3
yd=stem(yc1);

4
grid on

5
xlabel('Borders','FontSize',12);

6
ylabel('F (Hz)','FontSize',12);

7
set(gca,'YScale', 'log')


Jetzt muss ich eine Bandpassfilterbank zweiter Ordnung mit den 
gefundenen Grenzfrequenzen unter Verwendung eines Butterworth-Filters 
(-3 dB bei den Grenzfrequenzen) implementieren.
Ich weiß nur, dass freqz () Grenzfrequenzen im Bereich (0 .. 1) 
benötigt. Dies bedeutet, dass die Abtastfrequenz, die hier benötigt 
wird, größer als f_s sein muss.
Antwort Zitieren Bearbeiten Thread verschieben Thread sperren Anmeldepflicht aktivieren Anpinnen Thread löschen Thread mit anderem zusammenführen Beitrag melden
Dergute W. #6557442
Lesenswert?
• ▲
▼ 
Moin,

Preg M. schrieb:
> Dies bedeutet, dass die Abtastfrequenz, die hier benötigt
> wird, größer als f_s sein muss.

Die sollte sogar mindestens doppelt so gross sein. Wenns tatsaechlich 
auch in echt mit real aufzubauendem Antialiasing-TP vor dem ADC laufen 
soll, dann noch groesser.
Also solltest du dir da erstmal eine Samplefrequenz aussuchen/erfragen. 
Oft gibt die HW da schon was vor.
Wenn du die Samplefrequenz als Wert hast, normierst du deine 
Filterfrequenzen auf die Haelfte davon.
Beispiel: du nimmst 48kHz Samplefrequenz, dann entsprechen deine 8kHz 
dem Wert: 8kHz/(48Khz/2) = 1/3
Diese normierten Frequenzwerte kannst du dann z.b. in die "butter" 
Funktion einsetzen, die dir die Filterkoeffizienten fuer eines deiner 
Bandfilter zurueckliefert.

Edit: Kleine Betriebsblindheit meinerseits: Wenn du die Filterei eher im 
analogen machen willst, dann brauchst du keine Abtastfrequenz, aber dann 
wirds eher freqs() statt freqz() und hinten im butter() wird dann auch 
noch 's' stehen.

Gruss
WK
Antwort Zitieren Bearbeiten Löschen Beitrag melden
Thread beobachten

Antwort schreiben

Hinweis: der ursprüngliche Beitrag ist mehr als 6 Monate alt. Bitte hier nur auf die ursprüngliche Frage antworten, für neue Fragen einen neuen Beitrag erstellen.

Formatierung
  • [c]C-Code[/c]
  • [code]Code in anderen Sprachen, ASCII-Zeichnungen[/code]
  • [math]Formel in LaTeX-Syntax[/math]
  • [[Titel]] – Link zu einem Artikel
  • Verweis auf einen anderen Beitrag: Rechtsklick auf dessen Beitrags-ID (#8078628), „Adresse kopieren“, in den Text einfügen

Mehr zur Formatierung

Längeren Sourcecode als Anhang, nicht in den Text. JPG nur für Fotos und Scans; Zeichnungen und Screenshots als PNG oder GIF. Siehe Bildformate.

Dein Name und deine E-Mail-Adresse

Schon registriert? Anmelden

Wird nicht angezeigt. Wir schicken dir einen Link; sobald du damit bestätigst, erscheint dein Beitrag.

Mit dem Abschicken legst du einen kostenlosen Account an und erkennst die Nutzungsbedingungen und die Datenschutzerklärung an.