WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Leistungsdichte, pwelch

 

flyingernst
Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.07.2010, 10:13     Titel: Leistungsdichte, pwelch
  Antworten mit Zitat      
Hallo,

ich komme mit den vielen Eingabeparametern nicht zurecht, vielleicht kann mir einer helfen

Quellsignal: Kraft-Zeitverlauf, Kraft und Zeit jeweils in einem Vektor
Fs=250Hz

mit pwelch(Kraft) bekomme ich ein Bild, welches mir hoffentlich zumindest Qualitativ sagen kann ob es bei einer bestimmten Frequenz zu extremen Leistungsumsätzen kommt, siehe Anhang.
Wie bekomme ich hin dass die Skalierung der X-Achse nicht die normierte Frequenz sondern die tatsächliche Frequenz ist? Welche Parameter wären beim Aufruf des Befehls noch sinnvoll?!

Danke, Gruß Michael

untitled.tif
 Beschreibung:

Download
 Dateiname:  untitled.tif
 Dateigröße:  46.29 KB
 Heruntergeladen:  1895 mal
Private Nachricht senden Benutzer-Profile anzeigen


Fourier
Forum-Guru

Forum-Guru


Beiträge: 329
Anmeldedatum: 05.04.10
Wohnort: Karlsruhe
Version: 2009b
     Beitrag Verfasst am: 07.07.2010, 10:25     Titel:
  Antworten mit Zitat      
Hallo,

deine normalized frequenzen müssen mit Fs/2 multipliziert werden.
vlt so:
Code:

set(gca,'XTickLabel',(Fs/2)*get(gca,'XTick'))
 


wenn es klappen sollte musst du vlt noch Xlabel ändern.

ABER ich glaube bin mir aber nicht sicher wenn du bei dem befehl pwelch fs als Argument mitgibst, dann sind deine frequenzen nicht normiert.
Code:

pwelch(x,.....,....,[],Fs,'onesided')
 


gruss
_________________

Ein schlechter General ist besser als zwei gute.
Private Nachricht senden Benutzer-Profile anzeigen
 
flyingernst
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.07.2010, 10:40     Titel:
  Antworten mit Zitat      
das klappt, prima danke.

Ist es möglich den pwelch eines 2. Kanals in den gleichen Graphen zu setzen? Also nicht mit subplot sondern im gleichen Plot, weil dieser in meinem eigentlich Programm nämlich schon ein subplot ist
Private Nachricht senden Benutzer-Profile anzeigen
 
Fourier
Forum-Guru

Forum-Guru


Beiträge: 329
Anmeldedatum: 05.04.10
Wohnort: Karlsruhe
Version: 2009b
     Beitrag Verfasst am: 07.07.2010, 10:48     Titel:
  Antworten mit Zitat      
Code:


damit kannst du im gleichen graphen ploten, ohne dene rsten zu überschreiben.

es interessiert mich mit welcher Variante du dein erstes problem gelöst hast.
Code:

set(gca,'XTickLabel',(Fs/2)*get(gca,'XTick'))
 


oder
Code:

pwelch(x,.....,....,[],Fs,'onesided')
 

_________________

Ein schlechter General ist besser als zwei gute.
Private Nachricht senden Benutzer-Profile anzeigen
 
flyingernst
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.07.2010, 11:22     Titel:
  Antworten mit Zitat      
mit beidem also mit


Code:
PWELCH(test(:,2),[],[],[],250)


und

Code:
>> pwelch(test(:,2))
>> set(gca,'XTickLabel',(Fs/2)*get(gca,'XTick'))
>> xlabel('Frequenz')


klappt es.
das mit hold on geht auch, allerdings haben da beide Datensätze die gleiche Farbe. Ich kann zwar mit der Y-Skalierung auch nichts anfangen, aber es gibt mit zumindest mal einen überblick ob es hochfrequent hohe Leistungen gibt, die ich nicht vermutet habe

Vielen Dank, Gruß Michael
Private Nachricht senden Benutzer-Profile anzeigen
 
Fourier
Forum-Guru

Forum-Guru


Beiträge: 329
Anmeldedatum: 05.04.10
Wohnort: Karlsruhe
Version: 2009b
     Beitrag Verfasst am: 07.07.2010, 11:35     Titel:
  Antworten mit Zitat      
ich kanns leider gerade nicht testen, aber probieren geht vor studieren:
Code:

pwelch(.......)
.
.
hold on
h=pwelch(........)            % pwelch des 2.kanal

set(h,'Color','r')
 


wenn der 2.pwelch nichts plotet, dann muss man was anderes probieren (was ich leider gerade nicht kann)

oder vlt so:

Code:

pwelch(......,Fs,'Color','r')    % pwelch des 2.kanal
 

_________________

Ein schlechter General ist besser als zwei gute.
Private Nachricht senden Benutzer-Profile anzeigen
 
flyingernst
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.07.2010, 11:38     Titel:
  Antworten mit Zitat      
ja der plottet nicht sobald man etwas zurück geben lässt
Private Nachricht senden Benutzer-Profile anzeigen
 
Fourier
Forum-Guru

Forum-Guru


Beiträge: 329
Anmeldedatum: 05.04.10
Wohnort: Karlsruhe
Version: 2009b
     Beitrag Verfasst am: 07.07.2010, 11:40     Titel:
  Antworten mit Zitat      
2.variante auch nicht??
_________________

Ein schlechter General ist besser als zwei gute.
Private Nachricht senden Benutzer-Profile anzeigen
 
flyingernst
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.07.2010, 11:42     Titel:
  Antworten mit Zitat      
Code:
>> PWELCH(test(:,2),[],[],[],250,'Color','r')
??? Error using ==> signal\private\welchparse
Ambiguous or invalid string option specified.


nein, mit den Parametern kommt die funktion nicht zurecht
Private Nachricht senden Benutzer-Profile anzeigen
 
flyingernst
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.07.2010, 11:47     Titel:
  Antworten mit Zitat      
Moment, ich weis jetzt wie ich das plotte wenn ich werte ausgeben lasse

[Pxx1, Fxx1] = pwelch (test(:,1), [], [], [], 250);
[Pxx2, Fxx2] = pwelch (test(:,2), [], [], [], 250);
plot (Fxx1, 10 * log10 (Pxx1),Fxx2, 10 * log10 (Pxx2));

damit klappts, danke
Private Nachricht senden Benutzer-Profile anzeigen
 
flyingernst
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.07.2010, 13:30     Titel:
  Antworten mit Zitat      
Eine Frage habe ich jetzt doch noch.

Wenn ich jetzt einen Kraftverlauf als Datenquelle für pwelch nutze: Ist das überhaupt eine Leistungsdichte, was dort dargestellt wird? Welche Physikalische Größe muss denn verarbeitet werden, so dass Leistung/Frequenz überhaupt dargestellt wird und nicht irgendwas (die Y-Bezeichnung in meinem Graph ist wahrscheinlich quatsch)?

Danke, Gruß Michael

beispiel.tif
 Beschreibung:

Download
 Dateiname:  beispiel.tif
 Dateigröße:  106.08 KB
 Heruntergeladen:  1742 mal
Private Nachricht senden Benutzer-Profile anzeigen
 
flyingernst
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 13.07.2010, 16:42     Titel:
  Antworten mit Zitat      
kann mir hier nochmal bitte jemand helfen, auch wenn es eher eine physikalische Frage ist?
Normal wird für das Leistungsdichtespektrum (so wie ich zumindest im Internet gesehen und verstanden habe) für Frequenzgänge, also akkustische Messungen genommen. Ich jedoch habe einen Kraftverlauf (mit variablen Amplituden und Frequenzen). Reicht der Graph trotzdem für eine qualitative Beurteilung wieviel Arbeit bei einer bestimmten Frequenz vom lager geleistet werden muss während einer Messung

Danke, Gruß Michael
Private Nachricht senden Benutzer-Profile anzeigen
 
_Fourier

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 14.07.2010, 10:03     Titel: a
  Antworten mit Zitat      
Hallo,

ich will dir nichts falsches erzählen.
aber ich werde für dich den Meister in signalverarbeitung fragen (heute abend email schicken). wenn er antwortet sage ich dir bescheid Wink

deine frage interessiert mich natürlich.

gruss
 
flyingernst
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 25
Anmeldedatum: 13.04.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 14.07.2010, 16:05     Titel:
  Antworten mit Zitat      
das finde ich prima, vielen Dank

Gruß michael
Private Nachricht senden Benutzer-Profile anzeigen
 
Fourier
Forum-Guru

Forum-Guru


Beiträge: 329
Anmeldedatum: 05.04.10
Wohnort: Karlsruhe
Version: 2009b
     Beitrag Verfasst am: 16.07.2010, 07:44     Titel:
  Antworten mit Zitat      
Hallo Michael,

das ist die Antwort des Meisters kopiert und eingefügt, ich hoffe du kannst damit was anfangen:

Zitat:

schön, dass Sie sich mal wieder gemeldet haben.
Wenn das Signal eine Spannung ist dann ist die Einheit der spektralen
Dichte Volt^2/Hz. Volt^2 wird hier als Leistung für einen
Referenzwiderstand von 1 Ohm betrachtet.

Wenn das Signal eine mechanische Größe ist, wie z.B. Kraft in N
(Newton), dann nennt sich N^2 Leistung, die nicht direkt mit der
bekannten mechanischen Leistung verbunden werden kann, die Nm/s ist.
Diese Leistung wird aber nur mit der spektralen Leistungsdichte
korrekt so genannt.

Ähnlich ist z.B. eine Lage in m dann in Verbindung mit der
Spektralanalyse in m^2 als Leistung betrachtet.

Schöne Grüße und viel Erfolg


von mir auch viele Grüsse und viel Erfolg ;)
_________________

Ein schlechter General ist besser als zwei gute.
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen

Gehe zu Seite 1, 2  Weiter

Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.