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

Messdaten aufteilen

 

ThKo26
Forum-Century

Forum-Century


Beiträge: 184
Anmeldedatum: 21.09.17
Wohnort: ---
Version: 2015b
     Beitrag Verfasst am: 25.10.2019, 10:51     Titel: Messdaten aufteilen
  Antworten mit Zitat      
Hallo zusammen,

ich benötige wieder mal eure hilfe.

es geht um splitten der Messdaten.

Bild 1 zeigt eine Übersicht der Messdaten, diese wurden schon in data und warmup getrennt.

Im Bild 2 erkennt man, dass es zwei Sweeps gibt, die abwechselnd durchgeführt werden.

Habe on Top die Daten beigefügt.

Hat jemand einen Vorschlag, wie ich dies trennen kann?

Danke für die Hilfe.

Grüße Thomas

data.mat
 Beschreibung:

Download
 Dateiname:  data.mat
 Dateigröße:  41.27 MB
 Heruntergeladen:  272 mal
Messdaten_Ziel.PNG
 Beschreibung:

Download
 Dateiname:  Messdaten_Ziel.PNG
 Dateigröße:  10.9 KB
 Heruntergeladen:  219 mal
Messdaten_komplett.PNG
 Beschreibung:

Download
 Dateiname:  Messdaten_komplett.PNG
 Dateigröße:  11.34 KB
 Heruntergeladen:  211 mal
Private Nachricht senden Benutzer-Profile anzeigen


ThKo26
Themenstarter

Forum-Century

Forum-Century


Beiträge: 184
Anmeldedatum: 21.09.17
Wohnort: ---
Version: 2015b
     Beitrag Verfasst am: 28.10.2019, 18:00     Titel:
  Antworten mit Zitat      
Hallo zusammen,

ich habe bisschen versucht, komme aber nicht zum gewünschten Ziel.

Vll kann mir jemand weiterhelfen.

DAnke.


Code:


% Variante 1
peaks = logical(data.SA(:,1) < 6 & data.SA(:,1) > -6);
pospeaks = find(double(zero)==1);
deltapeaks = diff(pospeaks);
posdiffpeaks = (find(deltapeaks>4));

% Variante 2
pks = findpeaks(data.SA(:,1));
[pks,locs] = findpeaks(data.SA);
dpeaks = diff(locs);
positiondiffpeaks = (find(dpeaks>100));

figure

plot(data.ET,data.SA,data.ET(locs),pks,'o')
 


Grüße
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.448
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 28.10.2019, 18:23     Titel:
  Antworten mit Zitat      
Hallo,

Zitat:
komme aber nicht zum gewünschten Ziel.

Und wie genau weicht das, was du bekommst, von dem ab, was du willst?

Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
ThKo26
Themenstarter

Forum-Century

Forum-Century


Beiträge: 184
Anmeldedatum: 21.09.17
Wohnort: ---
Version: 2015b
     Beitrag Verfasst am: 30.10.2019, 14:19     Titel:
  Antworten mit Zitat      
Hallo Harald,

ich möchte den Zeitbereich herausbekommen, an dem der Sweep 0--> 12 --> 0 --> -12 -->0 stattfindet. (siehe Bild)

ich habe mit den Varianten versucht erst mal die Peaks zu erhalten, um von dort aus zum Nullpunkt zukommen. Dies war bis jetzt sehr umständlich.

ich habe mit folgenden Code
Code:

[pksmax,locsmax] = findpeaks(data.SA,data.ET,'MinPeakProminence',6,'minPeakheight',5,'Annotate','extents');
 

(siehe Bild wie vorher)
ein erstes Resultat, das eindeutig und robust ist.
Nun möchte ich aber gleichzeitig, es für den negativen Anteil haben.

weißt du, wie ich die negativen erhalte ?
Code:

findpeaks(data.SA,data.ET,'MinPeakProminence',5,'maxPeakheight',-5,'Annotate','extents')
 

Code:
Error using findpeaks>parse_inputs (line 270)
'maxPeakheight' is not a recognized parameter. For a list of valid name-value pair arguments, see the documentation for this function.


Danke und Grüße
Thomas

Sweep.PNG
 Beschreibung:

Download
 Dateiname:  Sweep.PNG
 Dateigröße:  49.46 KB
 Heruntergeladen:  190 mal
Private Nachricht senden Benutzer-Profile anzeigen
 
ThKo26
Themenstarter

Forum-Century

Forum-Century


Beiträge: 184
Anmeldedatum: 21.09.17
Wohnort: ---
Version: 2015b
     Beitrag Verfasst am: 30.10.2019, 14:34     Titel:
  Antworten mit Zitat      
ThKo26 hat Folgendes geschrieben:


weißt du, wie ich die negativen erhalte ?
Code:

findpeaks(data.SA,data.ET,'MinPeakProminence',5,'maxPeakheight',-5,'Annotate','extents')
 

Code:
Error using findpeaks>parse_inputs (line 270)
'maxPeakheight' is not a recognized parameter. For a list of valid name-value pair arguments, see the documentation for this function.





die Lösung:
Code:
[pksmin,locsmin] = findpeaks(-data.SA,data.ET,'MinPeakProminence',8,'minPeakheight',-5,'Annotate','extents');




Grüße
Thoma
Private Nachricht senden Benutzer-Profile anzeigen
 
ThKo26
Themenstarter

Forum-Century

Forum-Century


Beiträge: 184
Anmeldedatum: 21.09.17
Wohnort: ---
Version: 2015b
     Beitrag Verfasst am: 30.10.2019, 16:47     Titel:
  Antworten mit Zitat      
Hallo zusammen,

ich habe das Problem gelöst.

wenn es jemanden interessiert:

Code:


clc
clear
close all
load('C:\Users\P369309\Desktop\data.mat');

figpeaks = figure('Name', 'Find Peaks Overview');
[pksmax,locsmax] = findpeaks(data.SA,'MinPeakProminence',8,'minPeakheight',5,'MinPeakDistance',400,'Annotate','extents');
[pksmin,locsmin] = findpeaks(-data.SA,'MinPeakProminence',8,'minPeakheight',5,'MinPeakDistance',400,'Annotate','extents');
plot(data.ET,data.SA,data.ET(locsmax),pksmax,'o')
hold all, grid on
plot(data.ET,data.SA,data.ET(locsmin),-pksmin,'*')


for igh = 1:length(locsmax)
    startpeakrangeSA = data.SA(1:locsmax(igh));
    posstartpeakSA = find(startpeakrangeSA<0,1,'last');
    endepeakrangeSA = data.SA(locsmin(igh):end);
    posendepeakSA = find(endepeakrangeSA>0,1,'first')+locsmin(igh);
   
    sweepET = data.ET(posstartpeakSA:posendepeakSA);
    sweepSA = data.SA(posstartpeakSA:posendepeakSA);
   
    % infos on top, but not necessary;
    startpeakET = data.ET(locsmax(igh));
    startpeakrangeET = data.ET(1:locsmax(igh));
    posstartpeakET = data.ET(posstartpeakSA);
   
    endepeakET = data.ET(locsmin(igh));
    endepeakrangeET = data.ET(locsmin(igh):end);
    posendepeakET = data.ET(posendepeakSA);

    figure
    plot(sweepET,sweepSA);
   
end

 


Grüße
Thomas
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



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.