Filter 2-ter Ordnung

Walter L. Gast #7189949
Lesenswert?
• ▲
▼ 
Hallo,

ich möchte die Erg. der Implementierung in C eines "Filters" 2-ter 
Ordnung vergleichen (realisiert über den Zustandsraum und über die 
Z-Transformation).

Ausgangspunkt ist die Dgl. eines Reihenschwingkreises mit L=1H, C=1F, 
R=0,2Ohm. y ist die Spg. über C, u die Spg. über der Reihenschaltung von 
L,C,R.

y(zwei_punkt)+0,2*y(punkt)+y=u;

"punkt" ist die Ableitung nach der Zeit,
y die Ausgangsgröße, u die Eingangsgröße.
Die Eingangsgröße springt von 400 auf 600 alle Sek. hin und her.
y gebe ich auf einen DAC und zeige das Erg auf einem Osz.

Die Zustandsraumrealisierung funktioniert und liegt als C-Code vor.

Hat jemand eine Realisierung über die Z-Transformation in C?
Wichtig ist, dass die Abtastzeit „Tab“ im C-Code eingebbar ist, und in 
der Initphase entsprechend mit den Koeffizienten verrechnet wird.

Ein Austausch des Quellcodes ist selbstverständlich.

VG Walter
Antwort Zitieren Bearbeiten Thread verschieben Thread sperren Anmeldepflicht aktivieren Anpinnen Thread löschen Thread mit anderem zusammenführen Beitrag melden
Rezy Gast #7190147
Lesenswert?
• ▲
▼ 
Walter L. schrieb:
> Rezy->Leeeesen!!!!
> L,C,R sind gegeben.
> Zeige den den rekursiven c-code.

Du verlangst von mir (in unfreundlichem Ton), den Code zu zeigen und die 
Koeffizienten des Filters (mittels einer geeigneten Transformation) 
überhaupt erst einmal korrekt zu bestimmen? Und das auch noch 
unentgeltlich?

Nein, danke. Ich habe schon genügend Hinweise gegeben. Niemand kann dir 
hier den fertigen Code geben, OHNE das Filter zuvor zu digitalisieren 
(mit welcher Methode auch immer).

Und damit bin ich hier raus ;)
Antwort Zitieren Bearbeiten Löschen Beitrag melden
Walter L. Gast #7190323
Lesenswert?
• ▲
▼ 
Zu dem Gast: Es ist das Gleiche, als wenn ich sage, Zustandsraumlösung 
ist doch ganz einfach, schau ins Internet, brauchst nur die die 
Transitionsmatrix über das C.-H-Theorem lösen, musst nur aufpassen, wenn 
der Wert unter der Wurzel neg. wird, dann A hoch minus 1 und den 
Steuervektor noch für die parti. Lösung einbinden, und schon ist alles 
zusammen.

Über die Differenzen-Methode

y(zwei punkt) ~=(y(kT)-2*y(k(T-1))+y(k(T-2)))/(Tab*Tab) und

y(punkt) ~=(y(kT)-y(k(T-1)))/Tab

mit Tab = 0,1 funktioniert dies rekursiv jetzt auch. Auch hier gerne im 
Tausch gegen den rekursiven C-code gelöst über die Z-Transformation.

VG Walter
Antwort Zitieren Bearbeiten Löschen Beitrag melden
Thread beobachten

Antwort schreiben

Bitte melde dich an, um einen Beitrag zu schreiben.

Bestehender Account

oder

Mit Google-Account einloggen

Noch kein Account?

Die Registrierung ist kostenlos und dauert nur eine Minute.

Jetzt registrieren