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

Bewertung von Signalverläufen

 

Peppone
Forum-Anfänger

Forum-Anfänger


Beiträge: 12
Anmeldedatum: 12.03.15
Wohnort: ---
Version: R2013a und R2014b
     Beitrag Verfasst am: 10.06.2015, 10:12     Titel: Bewertung von Signalverläufen
  Antworten mit Zitat      
Hallo zusammen,

ich möchte einen Signalverlauf bewerten. Siehe hierzu den geladenen Plot.

Mir ist durch das Soll Signal der Sprung im Signal bekannt.
Das Ist Signal näher sich diesem Sprung an.

Kriterien wären zum Beispiel leichte Steigung am Anfang, dann stark (streng) monotone Steigung bis kurz vor erreichen Sollsignal und dann eben ein einlaufen auf Zielniveau.

Ich hab schonmal an eine Funktion / ein Skript aus der Regelungstechnik gedacht. Anregelzeit usw. Hab da aber nichts gefunden.
Ist hier jemand 'Fit' auf diesem Gebiet?

Als Anhang mal meine Signale.

Grüße Stefanie

und Dankeschön Wink

Messung_Ist.txt
 Beschreibung:
Ist

Download
 Dateiname:  Messung_Ist.txt
 Dateigröße:  2.08 KB
 Heruntergeladen:  511 mal
Messung_Soll.txt
 Beschreibung:
Soll

Download
 Dateiname:  Messung_Soll.txt
 Dateigröße:  638 Bytes
 Heruntergeladen:  524 mal
Signalverlauf.png
 Beschreibung:
Plot des Signales

Download
 Dateiname:  Signalverlauf.png
 Dateigröße:  12.47 KB
 Heruntergeladen:  549 mal
Private Nachricht senden Benutzer-Profile anzeigen


Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 23.06.2015, 11:54     Titel: Re: Bewertung von Signalverläufen
  Antworten mit Zitat      
Hallo Peppone,

Zitat:
ich möchte einen Signalverlauf bewerten. Siehe hierzu den geladenen Plot.

Ich sehe keinen geladenen Plot.

Gruß, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
Peppone
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 12
Anmeldedatum: 12.03.15
Wohnort: ---
Version: R2013a und R2014b
     Beitrag Verfasst am: 23.06.2015, 12:58     Titel:
  Antworten mit Zitat      
Hab ich nachgeholt.

MFG

Plot.fig
 Beschreibung:

Download
 Dateiname:  Plot.fig
 Dateigröße:  22.08 KB
 Heruntergeladen:  523 mal
Private Nachricht senden Benutzer-Profile anzeigen
 
DSP
Forum-Meister

Forum-Meister



Beiträge: 2.117
Anmeldedatum: 28.02.11
Wohnort: ---
Version: R2014b
     Beitrag Verfasst am: 23.06.2015, 19:04     Titel:
  Antworten mit Zitat      
Hallo,

kennst du die Abtastfrequenz deiner Signale? Ist das Sollsignal, also der Sprung das Eingangssignal auf dein System und das Istsignal also der Ausgang?

Wenn ja...würde ich das Signal zunächst mal mit einem Tiefpass filtern. Dann versuchen mit einem ARX oder ARMAX Modell das System zu schätzen. Ich würde mal bei einem Modell 2. Ordnung beginnen.
Private Nachricht senden Benutzer-Profile anzeigen
 
Peppone
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 12
Anmeldedatum: 12.03.15
Wohnort: ---
Version: R2013a und R2014b
     Beitrag Verfasst am: 24.06.2015, 07:32     Titel:
  Antworten mit Zitat      
Hi,

Abtastung in 10 ms also 100 Hz.

Nicht ganz richtig. Ich messe nur den Ausgang. Aus meinem Ausgang gebe ich mir dann den Sollsprung vor.
Mit:
Code:


Bei dem Signalverlauf habe ich bereits alles rausgefiltert was ich nicht brauchen kann.

Code:
TPF_Ist = designfilt('lowpassfir','FilterOrder',150,'CutoffFrequency',12,'SampleRate',100);




Ich habe bereits Modelle versucht zu erstellen. Problem dabei ist, dass der Signalverlauf auch mal deutliche Dellen (Signal_2.fig) aufweisen kann.
Im Endeffekt finde ich hier keine saubere Lösung um einen reproduzierbare Auswertung machen zu können.

Signal_2.fig
 Beschreibung:

Download
 Dateiname:  Signal_2.fig
 Dateigröße:  14.66 KB
 Heruntergeladen:  499 mal
Private Nachricht senden Benutzer-Profile anzeigen
 
DSP
Forum-Meister

Forum-Meister



Beiträge: 2.117
Anmeldedatum: 28.02.11
Wohnort: ---
Version: R2014b
     Beitrag Verfasst am: 24.06.2015, 15:55     Titel:
  Antworten mit Zitat      
Die angehängten Txt-files mit den Messdaten sind dann schon mit dem FIR gefiltert?

Das Sollsignal legst du nach was genau fest? Woher weißt, dass der Sprung genau zu dem Zeitpunkt t stattfindest wie du ihn gelegt hast und woher kennst du die exakte Sprunghöhe?

Wenn der Eingang nicht mit gemessen wird, werden solche Modelle vor allem mit derartigem welligem Verlauf ungenau. Es sind dann nur noch geschätzte Modelle auf Grund des Ausgangs, also AR oder ARMA Modelle. Einfacher wäre es eben wenn man das Eingangssignal auch hätte oder es zumindest annäherend exakt nachbilden kann, in dem der Sprungzeitpunkt und die Höhe bekannt sind.

Ich werde aber später mal was versuchen.
Private Nachricht senden Benutzer-Profile anzeigen
 
Peppone
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 12
Anmeldedatum: 12.03.15
Wohnort: ---
Version: R2013a und R2014b
     Beitrag Verfasst am: 25.06.2015, 08:54     Titel:
  Antworten mit Zitat      
Guten Morgen,

den Zeitpunkt des Sprunges kenn ich, da ich dort eben das "Knöpfchen" drücke -> Trigger.

Mit:
Code:

lege ich die Sprunghöhe nach der Histogrammmethode fest.


Als mit den Modellen hab ich keine brauchbaren Ergebnisse erzielen können. Ich kann zwar für eine ausgewählte Messung alles schön hinfeilen. Passt dann aber für eine andere Messung nicht mehr.
Private Nachricht senden Benutzer-Profile anzeigen
 
DSP
Forum-Meister

Forum-Meister



Beiträge: 2.117
Anmeldedatum: 28.02.11
Wohnort: ---
Version: R2014b
     Beitrag Verfasst am: 25.06.2015, 15:49     Titel:
  Antworten mit Zitat      
Welche Genauigkeit (Fit) erhälst du denn für einen Datensatz?
Evtl. repräsentiert ein Datensatz eben das System nicht gut genug oder aber die Störungen sind immer wieder unterschiedlich.

Hilft evtl. eine Mittelung mehrerer Datensätze und dann eine Schätzung des gemittelten Datensatzes?
Private Nachricht senden Benutzer-Profile anzeigen
 
Peppone
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 12
Anmeldedatum: 12.03.15
Wohnort: ---
Version: R2013a und R2014b
     Beitrag Verfasst am: 26.06.2015, 11:54     Titel:
  Antworten mit Zitat      
78 %

mh... das mit der Mitterlung wird glaub ich keine entsprechende Trennschärfe mehr liefern.
Private Nachricht senden Benutzer-Profile anzeigen
 
DSP
Forum-Meister

Forum-Meister



Beiträge: 2.117
Anmeldedatum: 28.02.11
Wohnort: ---
Version: R2014b
     Beitrag Verfasst am: 27.06.2015, 07:54     Titel:
  Antworten mit Zitat      
Ich komme mit einem linearen ARX Modell auch nur auf ca. 81% und mit einem nichtlinearen auf knapp 91%.

Code:

clear;
Fs = 100;
Ts = 1/Fs;
load('sprungantowrt_laupl.mat')
y = MessungIst;
u = MessungSoll;

dataset = iddata(y,u,Ts); % Datensatz erstellen

% Modell Ordnung festlegen, nk = Verzögerung um x Werte
% lineares ARX Modell
na=3; nb=1; nk=0;
sysid = arx(dataset,[na nb nk]);
% nichtlineares ARX Modell
na=3; nb=3; nk=0;
sysid_nl = nlarx(dataset,[na nb nk]);

sys = tf(d2d(sysid,Ts)) % Übertragungsfunktion erstellen
fpe1 = fpe(sysid)
[yh,fit] = compare(dataset,sysid,sysid_nl);
disp('fit ARX Modell'); disp(fit);
compare(dataset,sysid,'b.-',sysid_nl,'r');
T = 0:Ts:((length(y)-1)*Ts) ;
plot(T,y,'b.--');
hold on;
plot(T(1:length(u)),u,'b');
hold on;
% Sprungantwort des Modells auf Eingang -> MessungSoll
[y_sim,t_sim] = lsim(sys,u,T);
plot(t_sim,y_sim,'r');
ylabel('y(t)');
xlabel('Zeit in s');
title('Vergleich Modell und reales System');
legend('MessungIst', 'MessungSoll','estimated y')
grid on;
 


Ein ARMAX Modell liefert durchweg deutlich schlechtere Ergebnisse. Andere Modelle habe ich noch nicht getestet.

Aber mir stellt sich mittlerweile die Frage was du mit diesem Modell überhaupt anfangen willst, selbst wenn dessen Genauigkeit nahe 100% wäre. Auf das Soll Signal kommst du damit auch nicht. Dazu bedarf es einer Regelung, welche aber nur bei einem linearen Modell möglich ist. Das 2 Signal hat ja auch noch eine derart große Abweichung, in dessen Verlauf sich kein Einschwingen auf einen Sollwert erkennen lässt.

Wozu also die Bewertung der realen Signale?
Private Nachricht senden Benutzer-Profile anzeigen
 
Peppone
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 12
Anmeldedatum: 12.03.15
Wohnort: ---
Version: R2013a und R2014b
     Beitrag Verfasst am: 28.06.2015, 12:49     Titel:
  Antworten mit Zitat      
Es geht um Druck in einer Presse. Soll so schnell und schwingungsfrei wie möglich sein.
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.