1-port s-parameter nach rlc

OP #2959765
Hi,

kennt jemand eine Lösung um aus S11 (50 Ohm System) möglichst einfach 
auf ein Reihen bzw. Serienersatzschaltbild für einen RLC Schwingkreis zu 
kommen? (best fit)

Der DG8SAQ VNWA hat z.B. für Quarze ein Modul das die Quarz-Paramater 
aus S11 ermittelt. Scheint aber für ein einfaches RLC nicht zu 
funktionieren da er wohl die zweite Resonanz sucht und nicht findet.
Würde das dann gerne im lt-spice weiterverwenden.

Subcircuit nach dem Schema
.SUBCKT meinrlc 1 2
R1N 1 10 -5.000000e+001
R1P 10 11 1.000000e+002

*S11 FREQ REAL IMAG
E11 11 2 FREQ {V(10,2)}= R_I
+( 1.000000e+005, 9.667665e-001, 2.443275e-001)
...

Hab ich schon hinbekommen, taugt aber vermutlich kaum für .tran analysen 
da ich ja nur einen Ausschnitt aus der Ü-Funktion habe.
OP
Kap 5.15 geht in die Richtung, aber eigentlich wollte ich die Software 
zum Parameterfit nicht selber schreiben... Und vorsichtshalber sind die 
benutzten Gleichungen auch nicht abgedruckt.

Wird wohl auf so ein iteratives Verfahren mit kleinsten Fehlerquadraten 
rauslaufen...

Für's erste hab ich das jetzt mal für einen par. rlc Kreis von Hand im 
Ltspice hingefrickelt.
R(@f0), f0, und L(100 Hz) waren bekannt. Der Rest ist ein wenig Hin und 
Hergeschiebe zwischen L / C
Nicht ideal, aber für Sonntagabend ok.
OP
@EMU:
Mal sehen, ob ich die Woche abends dazukomme.
Muss noch Frohnarbeit tun, Weihnachtsfeiern besuchen, Baum kaufen, 
Geschenke usw. - Die "stade Zeit" halt ;-)

Btw. zu welchem Aufsatz gehört denn der verlinkte Anhang? Da steht was 
von Chapter 5 mit den Hintergründen.
Den Link zurückverfolgen.
http://www.brl.uiuc.edu/Downloads/
Dann nach bigelow suchen.


Serienschwingkreis mit R1, L1 und C1
Ich habe mal selber Formeln aufgestellt um aus S11 die Werte für R1, L1 
und C1 zu bestimmen. Zumindest für Startwerte sollte es gut genug sein.

Zin=Zo*(1+S11)/(1-S11) ; Zin=Z11, make a table or plot
R1 = min(Zin)
f0 frequency at min(Zin)
B = f2(Z==sqrt(2)*R1)-f1(Z==sqrt(2)*R1) ; from table or plot of Zin(f)
Q = f0/B
C1 = 1/(w0*Q*R1)
L1 = 1/(w0*w0*C1)
OP
octave ist noch genau so, wie ich es in Erinnerung habe...
Immerhin muss ich mir nicht erst noch cygwin antun.

Insgesammt klappt das Script (bei mir) nicht (direkt).
j ist z.B. predefined, das Fileformat das read_s11 erwartet ist wohl 
kein touchstone. Funktionen müssen scheinbar  erst definiert sein bevor 
man sie benutzen darf.


Muss ich erst hinfrickeln und meine octave-Kenntnisse würden nicht mal 
für eine Pizzabestellung reichen. Also wird's wohl etwas dauern.

Antwort schreiben

Eine eigene Frage? Dieser Thread ist über 13 Jahre alt. Stell sie in einem eigenen Thread – dort lesen mehr Leute mit, und sie wird öfter und schneller beantwortet.

Neuen Thread starten

Eine Lösung, Ergänzung oder Rückmeldung zu diesem Thread? Dann antworte hier.

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.

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.