RC Tiefpass 2.Ordnung

OP #4583095
Lesenswert?

Hallo zusammen,
Ich beschäftige mich momentan mit einem Tiefpass 2.Ordnung.
Ich habe zwei RC Tiefpässe hintereinander zu einem Tiefpass 2.Ordnung 
geschaltet. Auf diesen Tiefpass haben wir einen sehr schmalen Impuls 
(ich denke es ist eher Rechteck gemeint) gegeben mit einer niedrigen 
Wiederholfrequenz.
Dann haben wir das ganze mit einem FFT-Analysator gemessen und dadurch 
eben den Frequenzgang dargestellt. Der FFT-Analysator misst ja 
eigentlich das Spektrum, aber da hier auf den Eingang ein sehr schmaler 
Impuls gegeben wird ist das Spektrum auch der Frequenzgang.
So ich habe jetzt mal die Übertagungsfunktion durch die Multiplikation 
von zwei Tiefpässen berechnet und dann in Matlab geplottet
1 1 1
Gs = ---------- * ------------ = ------------------------------
sR1C1+1 sR2C2+1 s^2R1R2C1C2+s(R1C1+R2C2)+1
Ich habe dann in Matlab R2 auf 5kOhm und C1=C2 auf 4uF eingestellt.
R1 habe ich wie auch im Plot zu sehen auf 0,5kOhm und 10kOhm 
eingestellt.
So weit so gut, jetzt habe ich die Messung gegenübergestellt und dort 
sieht es fast so aus, als wäre die Grenzfrequenz für die verschiedenen 
Einstellungen von R die gleiche. Wir haben bei der Messung R2 fest 
gehabt und C1=C2, R1 haben wir mit einem Poti auf die minimale mittlere 
und maximale Position eingestellt. Die x-Achse läuft von 0-25kHz
Nun zu meiner Frage, in Matlab sieht man ja deutlich andere 
Grenzfrequenzen vor allem für R->0, lese ich das jetzt aus dem Bild der 
Messung nur falsch ab, oder gibt es da eine andere Erklärung dafür?
Angehängte Dateien:
Gast #4583112
Lesenswert?

Peter M. schrieb:
> lese ich das jetzt aus dem Bild der
> Messung nur falsch ab,

In der Simulation ist das Verhältnis Ausgang zu Eingang dargestellt.
Die Messung zeigt nur den Ausgang des Filters.
Du kannst nicht davon ausgehen daß der Rechteck-Impuls alle Oberwellen 
mit der gleichen Amplitude erzeugt.

Gruß Anja
OP #4583797
Lesenswert?

Schaltung ist aus zwei RC Tiefpässen aufgebaut
R1-----|--R2--------|
  @@        | @@@@@@@     |
          C1.        @@@@   C2

Die Gleichung sollte eigentlich so aussehen, da hat wohl irgendwas meine 
Formatierung zerstört:
@@@@@@@@@@@1@@@@@1
 @@@@@
Gs = ---------------- * -----------------
     @@@@@@@@@       sR1C1+1   @sR2C2

         @@@@@@@@@@@@@@@@
@@@                     1
   @   = --------------------------------------------------------
      @@@     s^2(R1R2C1C2)+s(R1C1+R2C2)+1

Ja das mit der Skalierung der Frequenzachse kann natürlich sein, dass es 
nur durch die lineare Frequenzachse so aussieht.

Irgendwie siebte in der Vorschau immer gut aus und wenn ich es absende 
dann verzieht sich alles....
Ich machs dann nachher nochmal am pc.
Die @ Zeichen sind übrigens einfach nur Platzhalter, da Leerzeichen wohl 
gerne mal ignoriert werden
OP #4583961
Lesenswert?

Ja ich habe mich verschrieben vorhin, ich wollte eigentlich schreiben, 
in der Vorschau sieht immer alles gut aus und nicht siebte, naja ich 
hätte das dann einfach so wie Helmut schreiben sollen, ist der 
einfachste Weg.
Ich habe jetzt mal von der Gleichung der Übertragungsfunktion ein Bild 
gemacht.
Wie kommst du auf den zusätzlichen Term, s*R1*C2, Helmut?

Aber davon unabhängig welche Übertragungsfunktion ich jetzt für die 
Simulation nutze, man sieht auf beiden Plots, dass die Grenzfrequenz für 
den Fall von R1->0 eine deutlich andere Grenzfrequenz hat als für die 
anderen beiden Einstellungen für R1-> Rmax/2 und R1-> Rmax
Ist das jetzt wirklich nur dadurch, dass bei der Messung mit dem Spec. 
Analysator eine lineare Frequenzachse vorliegt, oder gibt es noch 
weitere Gründe?
Angehängte Dateien:
OP #4584039
Lesenswert?

hmm, könntest du dann noch eine kurze Herleitung angeben, weil ich weiß 
nicht so recht wie ich auf den zusätzlichen Term kommen soll.

Ist eigentlich die lineare Frequenzachse oder fehlende Pufferstuffe 
maßgeblich für die Unterschiede zwischen der Messung und Simulation 
verantwortlich?
#4584745
Lesenswert?

Hallo,

Ich glaube du hast da ein grundlegendes Verständnisproblem.

Deine Übertragungsfunktion ist nicht gleich deine Sprungantwort!

Was du misst ist aber deine Sprungantwort (Ableitung vom Dirac Impuls).

Um zu sehen ob dein Ausgangsspektrum mit deiner Übertragungsfunktion 
übereinstimmt, musst du die Sprungfunktion (Sigma) Laplace 
Transformieren -> 1/s und mit deiner Übertragungsfunktion 
multiplizieren. Dann machst du eine Rückwandlung in den Zeitbereich 
(Matlab kann das) und von dort aus eine FFT. Das Spektrum vergleichst du 
mit deinem Gemessenen.

Wenn du die Übertragungsfunktion messtechnisch ermitteln willst ist es 
besser einen Frequenzsweep mit einem Funktionsgenerator zu machen und 
konstanter Amplitude. Dann nimmst du die Amplitude auf (rechnest Ein und 
Ausgang in RMS um) und machst in Matlab 20*log(Vout/Vin). Vergleichst 
Theorie mit Praxis und bist fertig.

Man kann sich den sweep auch ersparen und nimmt ein weißes Rauschen.

Gruß.
#4585381
Lesenswert?

Man schaut da mal kurz in die Help und schon sieht man, dass man da kein 
Eingangssignal anlegen muss/kann.

BODE(SYS,{WMIN,WMAX}) draws the Bode plot for frequencies between WMIN 
and WMAX (in radians/second).

Allerdings hilft auch Matlab nicht, wenn amn die falsche 
Übertragungsfunktion eingibt. Nimm die Funktion die ich angegeben habe, 
wenn du es richtig machen willst.
OP #4585397
Lesenswert?

Ich nehme auch parallel noch die Übertragungsfunktion von dir, es 
verschiebt ein wenig die Grenzfrequenz aber ansonsten ist es von der 
Form her gleich.
Ja ich hab da ein wenig durcheinander gekommen, weil wir bei früheren 
Versuchen Eingangs- und Ausgangsspannungen gemessen haben für 
verschiedene Frequenzen und damit dann das Bode Diagramm erstellt haben, 
aber ich habe eben auch nochmal nachgelesen, das Bode-Diagramm 
beschreibt wirklich nur rein die Übertragungsfunktion.
OP #4585517
Lesenswert?

ja das ist ein guter Hinweis, danke.

Kann mir eventuell noch jemand bei dieser Schaltung hier helfen?

Ich habe jetzt zwar einige Beispiele für Widerstände in Reihe gefunden, 
aber keine mit welchen parallel.
Insbesondere die Übertragungsfunktion wäre wieder interessant, damit ich 
das ganze wieder simulieren kann. In der Messung habe ich festgestellt 
das für ein sehr kleines R, der Ausgang gegen 0 geht und für ein großes 
R der Ausgang gegen die Eingangspannung geht.
Angehängte Dateien:
OP #4587103
Lesenswert?

Mit der Formel

Ua/Ue = (R2/(R1+R2))*(1+s*R1*C1)/(1+s*R2*C2)

bekomme ich in Bode irgendwie keine horizontale Verläufe, so wie ich sie 
in der Messung hatte.
Also ich habe jetzt in Matlab C1 gleich C2 gesetzt und R1=10kOhm und 
R2=1 Ohm, dann habe ich R2=5 kOhm und R1=5k Ohm und als letztes R2=10 
kOhm und R1=10 Ohm
In der Messung hatte ich Sprünge auf die maximale Spannung vom Sprung 
mit R2->0hm, auf die halbe Spannung mit R2 auf Rmax/2 und für R2->Rmax 
ging die Spannung gegen 0.
Gast #4587122
Lesenswert?

Nico B. schrieb:
> Deine Übertragungsfunktion ist nicht gleich deine Sprungantwort!
>
> Was du misst ist aber deine Sprungantwort (Ableitung vom Dirac Impuls).

Als Sprungantwort wird üblicherweise ein Sprung von einer Spannung U1 
auf U2 verstanden und nicht die Antwort auf einen kurzen Impuls (Dirac 
Impuls, idealerweise unendlich kurz und unendlich hoch).

Peter M. schrieb:
> Hmm seltsam, bei uns im Labor haben sie gemeint, es ist näherungsweise
> der Frequenzgang.

Die Fouriertransformierte der Impulsantwort ist die 
Übertragungsfunktion.
#4590448
Lesenswert?

@ Wolfgang


Wolfgang schrieb:
> Als Sprungantwort wird üblicherweise ein Sprung von einer Spannung U1
> auf U2 verstanden und nicht die Antwort auf einen kurzen Impuls (Dirac
> Impuls, idealerweise unendlich kurz und unendlich hoch).


Auch wenn Wikipedia keine gute Quellenangabe ist:

https://de.wikipedia.org/wiki/Sprungantwort

Ich habe explizit von der Sprungantwort geredet, da die Sprungfunktion 
leichter zu erzeugen ist als ein Dirac Impuls. Ableitung vom Dirac 
Impuls -> Sprungfunktion.

Wolfgang schrieb:
> Die Fouriertransformierte der Impulsantwort ist die
> Übertragungsfunktion.

Absolut richtig. Gedankenfehler von mir.

@ Peter

Du hast deinen errechneten Frequenzgang mit doppel log dargestellt. Auf 
dem Analysator ist die Darstellung (zumindest der x-Achse) linear.

Antwort schreiben

Bitte melde dich an, um einen Beitrag zu schreiben.

oder

Mit Google-Account einloggen

Die Registrierung ist kostenlos und dauert nur eine Minute.

Jetzt registrieren