Re: Extrahierung von Polen und Nullstellen aus Messdaten

Gast #6356198 
Hallo zusammen,

Ich habe die Messdaten eines Bodediagramms  gespeichert und wollte bitte 
aus diesen Messdaten Pole und Nullstellen in Matlab extrahieren. Aus 
diesen Messdaten:
-Freq ist die Frequenz
-cs entspricht der Amplitude in dB
-phis entspricht der Phase.
-CTF5 ist der Dateiname

Ich würde mich auf Ihre Hilfe sehr freuen

Viele Grüße
Thierry Assop
Angehängte Dateien:
Antwort Zitieren Bearbeiten Thread verschieben Thread sperren Langsamer Modus Anmeldepflicht aktivieren Anpinnen Thread löschen Thread mit anderem zusammenführen Melden Danke
Gast 
Ich will aus meine Messdaten (CTF5.txt) die pole und Nullstellen mit
Matlab berechnen lassen.
Welche Befehle soll ich benutzen, um die Pole und Nullstellen aus meine 
Messdaten zu extrahieren?
Antwort Zitieren Bearbeiten Löschen Melden Danke
Gast 
Assop schrieb:
> Ich will aus meine Messdaten (CTF5.txt) die pole und Nullstellen mit
> Matlab berechnen lassen.

Falls dein Name nicht eigentlich mit 2 a geschrieben wird: wo sie 
liegen, kannst Du ja sehen, wenn Du die Daten plottest. Und was kannst / 
machst Du mit Matlab?
Antwort Zitieren Bearbeiten Löschen Melden Danke
Gast 
Hallo zusammen,
danke für eure Rücckmeldung.
@Joggel, ich habe versucht, so
- Null(CTF5.txt)
- Pole(CTF5.txt)
zu machen, es kommt aber ein Fehler. Hast du es auch zufällig mit meinen 
Messdaten versucht? Bei kommt immer einen Fehler
Antwort Zitieren Bearbeiten Löschen Melden Danke
Gast 
Um die Null- und Polstellen zu bestimmen, brauchst Du die 
Übertragunsfunktion auf der gesamten komplexen Ebene. Das Bodediagramm 
bildet aber nur die imaginäre Achse ab. Du musst also versuchen, die 
Funktion H(j*omega) analytisch auf die gesamte C- Ebene fortzusetzen. 
Wie das in Matlab geht, kann ich nicht genau sagen, aber schau Dir mal 
die Matlab- Funktion rationalfit an.
Antwort Zitieren Bearbeiten Löschen Melden Danke
Gast 
Hi @Mike,

danke für deine Hilfe. Es ist erledigt worden.

load CTF5.txt;
Freq=CTF5(:,1);
cs=CTF5(:,2);
phis=CTF5(:,3);
PHA = phis; %phase of experimental result
AMP = cs; %gain of experimental result
W = Freq*2*pi; %transforming frequency in Hz to rad/s
Ts = 0;
gain=10.^(AMP/20);
zfr = gain.*exp(1i*PHA*pi/180);
gfr = frd(zfr,W,Ts);
sys=tfest(gfr,2);
[mag,ph]=bode(sys,w);
mag=squeeze(mag);
ph=squeeze(ph);
Antwort Zitieren Bearbeiten Löschen Melden Danke
Thread beobachten

Antwort schreiben

Eine eigene Frage? Dieser Thread ist über 6 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.

Dein Name und deine E-Mail-Adresse

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.