|
Anzeige
|
Hilfestellung zu FIR-Filter
Gast
#764895
Für die Filterung soll folgende Differenzengleichung verwendet werden: y(k)=x(k)+3x(k-1)+x(k-2)+0.231y(k-1)-0.123y(k-2) Nun soll man eine C-Funktion schreiben "filter_FIR", die als Parameter einen Vector mit Eingangswerten und einen Vektor mit Ausgangswerten sowie die Länge der beiden Vektoren bekommt und die Filterung durchführt.
Stimmt die Umsetzung der Differenzengleichung in die C-Funktion so? Abgesehen davon, dass du dem Compiler nicht sagst, dass er vor der Klammer ne Multiplikation machen soll, sieht es schon ganz gut aus. Evtl könntest du noch den x, bzw y Teil von k-1 zusammenfassen.Ist aber nicht nötig, könnte aber eine Multiplikation sparen.
Gast
#764920
Hallo Marcel, denke so einfach geht das nicht. k in der beschreibenden Gleichung ist doch sicher sowas wie die Zeit (zu festen Abtastwerten Delta-T). Das heißt für y(k)=x(k)+3x(k-1)+x(k-2) ... dass der Ausgangswert y zur Zeit k durch eine Summe von x- und y-Werten zur Zeit k, zur Zeit (k-1) (also das Sample davor) und (k-2), also das zweitletzte Sample gebildet wird. Hierzu musst du mindestens die 2 letzten Zeitwerte von x und y zwischenspeichern, um sie in der aktuellen Summe parat zu haben (und sie in jedem neuen Zeitschritt auch "durchshiften"). Weiteres Problem: Wie startet der Algorithmus. Für das erste Zeitsample sind die Felder für (k-1) und (k-2) ja noch leer. Nicht dass das falsch ist. Solltest dir aber dessen bewusst sein. Gruß Gerhard
Gast
#764925
Sieht eigentlich OK aus
void filter_FIR(double *x, double *y, int rec)
{
int n;
for(n=0;n<rec;n++)
y[n]=x[n]+3*x[n-1]+x[n-2]+0.231*y[n-1]-0.123*y[n-2];
}
Bei Arrays wird der Index in [] Klammern gesetzt
Beim Aufruf der Funktion musst du allerdings beachten
das du als Addresse nicht x oder y angibst sondern das ende des Arrays
also x[2] oder y[2] weil du ja in deiner Funktion negative Array werte
bekommst.
void Test(void)
{
double x[3];
double y[3];
filter_FIR(&x[2],&y[2],3);
}
Gruss Helmi
Gast
#764993
Recht herzlichen Dank für eure Hilfe. Noch ein Problem hab ich noch. Und zwar soll nun ein weitere C-Funktion mit dem Namen "filter_Forever" implementiert werden, der eine andauernde Filterung mit Hilfe der Funktion "filter_FIR" durchführen soll. Wie kann ich dies machen?
Gast
#765076
Dies könnte ich doch dann im Hauptprogramm in einer WHILE Schleife tun. Bin mir nicht ganz sicher.
Gast
#765086
Könntest du so machen. Allerdings werden da nirgendwo neue Messdaten eingelesen oder ergebnisse ausgegeben. Gruss Helmi
Gast
#765492
Wie würde die C-Funktion "filter_FIR" aussehen, wenn man das ganze mit Gedächtnis machen will? Mit Gedächtnis meine ich, dass man nur zum Beispiel wenn ein Block 100 lange ist, das nur die ersten 80 Werte benutzt werden und die anderen 20 erst beim nächsten durchlauf verwendet werden.
Gast
#765514
Wie wär's mal mit zuhören in der Vorlesung oder den Dozenten fragen, wenn man was nicht versteht? Herr Q. ist äusserst hilfsbereit und erklärt's auch der letzten Torfnase noch einmal ausführlich. Aber einen Tag vor der Klausur natürlich nicht mehr... Stürzt euch das Programm eigentlich nicht ab, da ihr auf ein negativen Array Index zugreift? Beim ersten Durchlauf der Schleife ist n = 0 und ihr addressiert n-1, n-2 !
Gast
#765529
@Nils H. (iaby) Das kann nicht abstuerzen weil beim aufruf der Funktion die Addresse vom 3. Element uebergeben wird. filter_FIR(&x[2],&y[2],3); Gruss Helmi Ok stimmt, mein Fehler! Wenn ich nun aber das 2te Element übergebe und n läuft bis 2, dann greif ich ja auf Index 4 zu, das Array hat aber nur 3 Elemente, also Index 0-2. Oder steh ich schon wieder auf dem Schlauch?!?
Gast
#765585
@Nils H. (iaby) >Wenn ich nun aber das 2te Element übergebe und n läuft bis 2, dann greif >ich ja auf Index 4 zu, das Array hat aber nur 3 Elemente, also Index >0-2. >Oder steh ich schon wieder auf dem Schlauch?!? Nein du stehst nicht auf dem Schlauch du hast recht das Array muss groesser sein. Gruss Helmi Antwort schreiben |
Anzeige
|
