Aus Sweep per FFT Frequenzgang < Matlab < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) beantwortet  | | Datum: | 08:16 Fr 24.06.2005 | | Autor: | tronny |
Ersteinmal finde ich es Klasse, dass es dieses Forum gibt und hier soviele Infos zu finden sind.
Nun zu meinem Problem: Ich mache momentan eine Regelstrecken Untersuchung und habe zur Anregung einen Sinus Sweep in dem interessierenden Frequenzbereich zur Anregung auf die Strecke gegeben. Das Ausgangssignal wurde in Form von diskreten Werten abgetastet. Nun wollte ich dieses Ausgangssignal mit Matlab per FFT auswerten und mir daraus das Bodediagramm erstellen. Ich bekomme zwar mit fft(out) und einem anschliessenden Betrag der Transformation schon einen ungefähren Verlauf heraus, der sich aber nicht in ein Bodediagramm umsetzen lässt.
Hat irgendjemand Tips zum Vorgehen, oder einen Link, denn dieses Problem der Untersuchung sollte doch häufiger vorkommen?
Vielen Dank schon mal im voraus!
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
| |
|
hi,
das ganze geht so ziemlich in einer zeile
bode diagramm ist ja nicht anderes als die übertragungs funktion in db dargestellt
also
plot(log(fftshift(abs(fft(in))./abs(fft(out)))));
wobei gesagt sei das das so keine echten db einheiten sind und du nur eine hälft des plots brauchst.
es gibt auch den befehl db um direkt db einheiten zu erzeugen.
plot(db(fftshift(abs(fft(in))./abs(fft(out)))));
aber da les am besten nochmal die hilfe genau durch.
für die phase gilt ähnliches
plot(fftshift(angle(fft(in))-angle(fft(out))));
hoffe es hilft dir weiter.
achja. google ist eine so wunderbare suchmaschine *g* 
grüße Physikus
|
|
|
| |
|
| Status: |
(Frage) beantwortet | | Datum: | 21:46 So 03.07.2005 | | Autor: | tronny |
Ersteinmal Danke für die Antwort. Mein Problem ist nur, dass der Bode Plot nicht meinen Erwartungen entspricht. Mich interessiert bei der Struktur nur ein recht kleines Frequenzfenster bis 50Hz. Ich bekomme aber bei dieser Art der fft im Frequenzbereich kein Maximum zu sehen, obwohl ich im Zeitbereich eindeutig ein Maximum ausmachen kann, welches sich aber in der fft nicht wiederfinden lässt.
Ich bin etwas ratlos. Vielleicht ist die Anregung mit einem Sinus Sweep auch nicht die Beste Wahl?
|
|
|
| |
|
servus tronny,
ich denke mal das proble wird im sweep liegen.
du zeichest ein siganl [mm] $f(\omega(t) [/mm] ,t)$ auf. So ein siganl kannst du nicht ohne weiteres fourier transformieren. benutze ein rechteck impuls, da hast du zwar alle freuqenzen drin aber du kannst dir ja den bereich nachher raussuchen der dich interressiert. ansonsten müsstest du dir die übertragungs funktion für jede einzelne freuqenz anschauen. kein sweep!!!
grüße physikus
|
|
|
|
