WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Erstellung eines Filters

 

josekamara
Forum-Meister

Forum-Meister


Beiträge: 529
Anmeldedatum: 26.04.08
Wohnort: München
Version: ---
     Beitrag Verfasst am: 26.09.2008, 15:50     Titel: Erstellung eines Filters
  Antworten mit Zitat      
Servus Allerseits,

wenn man z.B den Verlauf eines Filters bzw seiner Übertragungsfkt. kennt, d.h gegeben ist Amplitude A über die Frequenz f punktweise. Das zu filternde Signal sei X.
Zum Erstellen des Filters geht man folgendermassen vor:
Code:

[b,a]=freqz(A,f,4,5);
X_gefilter = filter(b,a,X);
 


Das Problem ist dass die bekannte Frequenzachse nicht normalisiert ist also nicht zwischen 0 und pi, z.b:
f = [0 2 .........400];

Wie kann ich es jetzt machen damit das Signal auch richtig gefiltert wird?

Ich freue mich auf jeden Vorschlag.

Jose.[/code]
_________________

Simulation
Private Nachricht senden Benutzer-Profile anzeigen


student
Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 10.03.08
Wohnort: ---
Version: 7.0
     Beitrag Verfasst am: 26.09.2008, 17:29     Titel:
  Antworten mit Zitat      
Hallo,

vieleicht habe ich eine Lösung für dich.

Die Winkelgeschwindigkeit wird allgemein als
Zitat:
omega = 2 * pi / T = 2 * pi * f

beschrieben.
Die Sampling Frequenz Fs kennst du bereits, da du sie schon eingegeben hast Wink [H,F] = FREQZ(B,A,N,Fs)
Die MAXIMALE-frequenz, die "gesehen" werden kann ist gerade die sampling frequenz
Zitat:
omega_max = 2 * pi *fs

Da aber das Shannon/Nyquist Theorem beachtet werden sollte (muss), ergibt sich sinnvolle Werte für alle Frequenzen unter
Zitat:
omega = pi * fs

Normalerweise werden die Winkelgeschwindigkeiten noch mit der Samplingfrequenz normiert um vergleichbare resultate zu bekommen.

Also, (endlich) die eigentliche Lösung für dein Problem.
Multipliziere den Frequenzvektor f = [0 2 .........400]
Zitat:
w=2*pi*f

und normiere mit
Zitat:
w_n=w/fs

Verwende NUR die resultate bis w_n<=pi!!!!!!!! (Amplituden von grösserer Frequenz sind alles nur Aliase Wink )

Viele Grüsse Student
Private Nachricht senden Benutzer-Profile anzeigen
 
josekamara
Themenstarter

Forum-Meister

Forum-Meister


Beiträge: 529
Anmeldedatum: 26.04.08
Wohnort: München
Version: ---
     Beitrag Verfasst am: 27.09.2008, 16:07     Titel:
  Antworten mit Zitat      
Kann man nicht auf die Grenzfrequenz z.b hier 400 normieren:
W_n=pi* fi/fn
_________________

Simulation
Private Nachricht senden Benutzer-Profile anzeigen
 
josekamara
Themenstarter

Forum-Meister

Forum-Meister


Beiträge: 529
Anmeldedatum: 26.04.08
Wohnort: München
Version: ---
     Beitrag Verfasst am: 06.10.2008, 14:36     Titel:
  Antworten mit Zitat      
Irgendwie läuft mit dem Filter überhaupt nicht.
Normalerweise wie ich sagte, zur Filterung eines Signal X durch beliebige Übertragungsfkt. geht es folgendermaßen:
den Verlauf der Übertragungsfunktion nach Frequenzachse ist scho gegeben: [H,F];

Zur Berechnung der Zählen und Nenner-Koeffizienten der Übertragungsfunktion:

[code]
[bb,aa]=invfreqz(H,F,n,m);

Dann braucht ich den Signal X zu filtern:

[code]
X_neu =filter(bb,aa,XX);

Aber der Filter gibt mir meisten NaN-Werte,
Bitte sobald eine Idee hat, würde ich die dringend brauchen.

Danke Euch.

Jose[/code]
_________________

Simulation
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.