Verfasst am: 18.03.2013, 11:51
Titel: Plot in Funktion übergeben
Hallo Liebe MATLAB-Gemeinde,
Ich hab ein kleines Problem.
Ich erkläre das mal kurz und poste dann meine Codes.
Also ich hab eine kleine GUI mit guide gebastelt...
In dieser GUI kann man einiges auswählen und dann kann man ein Programm starten. Alles ohne prob.
Jetzt soll folgenes passieren.
Meine aufgerufene Funktion rechnet einwenig.. auch ohne Probleme.
Meine eigene Funktion ruft jedenfalls eine weitere Funktion auf die ich im I-net gefunden habe...
In dieser Funktion wird etwas geplotet.. Ich möchte aber gerne das die plots nicht auf meinem Bildschirm auftauchen... da sie in der GUI geplotet werden sollen.
Ich will die plots aber unter einen bestimmten Namen speichern wie bekomm ich das hin, denn wenn ich den obenstehenden Code verwende und ich das bild wieder mit gcf aufrufe macht er eine neue Figure.
Ich poste mal meine Codes:
Das aus dem I-net:
Code:
functionvarargout = peakfinder(x0, sel, thresh, extrema)
%PEAKFINDER Noise tolerant fast peak finding algorithm
% INPUTS: % x0 - A real vector from the maxima will be found (required) % sel - The amount above surrounding data for a peak to be % identified (default = (max(x0)-min(x0))/4). Larger values mean % the algorithm is more selective in finding peaks. % thresh - A threshold value which peaks must be larger than to be % maxima or smaller than to be minima. % extrema - 1 if maxima are desired, -1 if minima are desired % (default = maxima, 1) % OUTPUTS: % peakLoc - The indicies of the identified peaks in x0 % peakMag - The magnitude of the identified peaks
%
% [peakLoc] = peakfinder(x0) returns the indicies of local maxima that % are at least 1/4 the range of the data above surrounding data.
%
% [peakLoc] = peakfinder(x0,sel) returns the indicies of local maxima % that are at least sel above surrounding data.
%
% [peakLoc] = peakfinder(x0,sel,thresh) returns the indicies of local % maxima that are at least sel above surrounding data and larger % (smaller) than thresh if you are finding maxima (minima).
%
% [peakLoc] = peakfinder(x0,sel,thresh,extrema) returns the maxima of the % data if extrema > 0 and the minima of the data if extrema < 0
%
% [peakLoc, peakMag] = peakfinder(x0,...) returns the indicies of the % local maxima as well as the magnitudes of those maxima
%
% If called with no output the identified maxima will be plotted along % with the input data.
%
% Note: If repeated values are found the first is identified as the peak
%
% Ex: % t = 0:.0001:10; % x = 12*sin(10*2*pi*t)-3*sin(.1*2*pi*t)+randn(1,numel(t)); % x(1250:1255) = max(x); % peakfinder(x)
%
% Copyright Nathanael C. Yoder 2011 (nyoder@gmail.com)
% Perform error checking and set defaults if not passed in error(nargchk(1,4,nargin,'struct'));
error(nargoutchk(0,2,nargout,'struct'));
s = size(x0);
flipData = s(1) < s(2);
len0 = numel(x0);
if len0 ~= s(1) && len0 ~= s(2) error('PEAKFINDER:Input','The input data must be a vector') elseifisempty(x0) varargout = {[],[]};
return;
end if ~isreal(x0) warning('PEAKFINDER:NotReal','Absolute value of data will be used')
x0 = abs(x0);
end
ifnargin < 2 || isempty(sel)
sel = (max(x0)-min(x0))/4;
elseif ~isnumeric(sel) || ~isreal(sel)
sel = (max(x0)-min(x0))/4;
warning('PEAKFINDER:InvalidSel',...
'The selectivity must be a real scalar. A selectivity of %.4g will be used',sel) elseifnumel(sel) > 1 warning('PEAKFINDER:InvalidSel',...
'The selectivity must be a scalar. The first selectivity value in the vector will be used.')
sel = sel(1);
end
ifnargin < 3 || isempty(thresh)
thresh = [];
elseif ~isnumeric(thresh) || ~isreal(thresh)
thresh = [];
warning('PEAKFINDER:InvalidThreshold',...
'The threshold must be a real scalar. No threshold will be used.') elseifnumel(thresh) > 1
thresh = thresh(1);
warning('PEAKFINDER:InvalidThreshold',...
'The threshold must be a scalar. The first threshold value in the vector will be used.') end
ifnargin < 4 || isempty(extrema)
extrema = 1;
else
extrema = sign(extrema(1)); % Should only be 1 or -1 but make sure if extrema == 0 error('PEAKFINDER:ZeroMaxima','Either 1(for maxima) or -1(for minima) must be inputfor extrema');
end end
x0 = extrema*x0(:); % Make it so we are finding maxima regardless
thresh = thresh*extrema; % Adjust threshold according to extrema.
dx0 = diff(x0); % Find derivative
dx0(dx0 == 0) = -eps; % This is so we find the first of repeated values
ind = find(dx0(1:end-1).*dx0(2:end) < 0)+1; % Find where the derivative changes sign
% Include endpoints in potential peaks and valleys
x = [x0(1);x0(ind);x0(end)];
ind = [1;ind;len0];
% x only has the peaks, valleys, and endpoints
len = numel(x);
minMag = min(x);
% Set initial parameters for loop
tempMag = minMag;
foundPeak = false;
leftMin = minMag;
% Deal with first point a little differently since tacked it on % Calculate the sign of the derivative since we taked the first point % on it does not neccessarily alternate like the rest.
signDx = sign(diff(x(1:3)));
if signDx(1) <= 0% The first point is larger or equal to the second
ii = 0;
if signDx(1) == signDx(2)% Want alternating signs
x(2) = [];
ind(2) = [];
len = len-1;
end else% First point is smaller than the second
ii = 1;
if signDx(1) == signDx(2)% Want alternating signs
x(1) = [];
ind(1) = [];
len = len-1;
end end
% Preallocate max number of maxima
maxPeaks = ceil(len/2);
peakLoc = zeros(maxPeaks,1);
peakMag = zeros(maxPeaks,1);
cInd = 1;
% Loop through extrema which should be peaks and then valleys while ii < len
ii = ii+1; % This is a peak % Reset peak finding if we had a peak and the next peak is bigger % than the last or the left min was small enough to reset. if foundPeak
tempMag = minMag;
foundPeak = false;
end
% Make sure we don't iterate past the length of our vector if ii == len
break; % We assign the last point differently out of the loop end
% Found new peak that was lager than temp mag and selectivity larger % than the minimum to its left. if x(ii) > tempMag && x(ii) > leftMin + sel
tempLoc = ii;
tempMag = x(ii);
end
ii = ii+1; % Move onto the valley % Come down at least sel from peak if ~foundPeak && tempMag > sel + x(ii)
foundPeak = true; % We have found a peak
leftMin = x(ii);
peakLoc(cInd) = tempLoc; % Add peak to index
peakMag(cInd) = tempMag;
cInd = cInd+1;
elseif x(ii) < leftMin % New left minima
leftMin = x(ii);
end end
% Check end point if x(end) > tempMag && x(end) > leftMin + sel
peakLoc(cInd) = len;
peakMag(cInd) = x(end);
cInd = cInd + 1;
elseif ~foundPeak && tempMag > minMag % Check if we still need to add the last point
peakLoc(cInd) = tempLoc;
peakMag(cInd) = tempMag;
cInd = cInd + 1;
end
% Create output
peakInds = ind(peakLoc(1:cInd-1));
peakMags = peakMag(1:cInd-1);
else% This is a monotone function where an endpoint is the only peak [peakMags,xInd] = max(x);
if peakMags > minMag + sel
peakInds = ind(xInd);
else
peakMags = [];
peakInds = [];
end end
% Apply threshold value. Since always finding maxima it will always be % larger than the thresh. if ~isempty(thresh)
m = peakMags>thresh;
peakInds = peakInds(m);
peakMags = peakMags(m);
end
% Rotate data if needed if flipData
peakMags = peakMags.';
peakInds = peakInds.';
end
% Change sign of data if was finding minima if extrema < 0
peakMags = -peakMags;
x0 = -x0;
end % Plot if no output desired ifnargout == 0 ifisempty(peakInds) disp('No significant peaks found') else figure;
plot(1:len0,x0,'.-',peakInds,peakMags,'ro','linewidth',2);
set(gcf,'Position',[0,0,1500,1300]); %Größe Einstellen
end else varargout = {peakInds,peakMags};
end
x=1;
while(1)
File=names{x}; %Dateien einzeln nehmen
comma2point(File) %extra Funktion aus dem I-net gefischt und verändert Ersetzt die Kommas durch Punkte MATLAB braucht Punkte
Einlesen=dlmread(File,'\t',8,0); %Da die Dateien Einen Text Kopf und auch die Spalten Text enthalten und MATLAB nicht viel damit anfagen kann werden sie Rausgefilter (Nicht eingelesen)
Einlesen(:,4)=Einlesen(:,4).*1000;
Einlesen(:,6)=Einlesen(:,2).*Einlesen(:,3);
Einlesen(:,7)=Einlesen(:,4).*Einlesen(:,5);
Verfasst am: 18.03.2013, 12:23
Titel: Re: Plot in Funktion übergeben
Hallo Florre,
Florre hat Folgendes geschrieben:
Ich möchte aber gerne das die plots nicht auf meinem Bildschirm auftauchen... da sie in der GUI geplotet werden sollen.
Wenn ich mich nicht täusche, musst du einfach den handle der axes in die geplottet werden soll an die Funktion plot übergeben.
Florre hat Folgendes geschrieben:
Ich will die plots aber unter einen bestimmten Namen speichern wie bekomm ich das hin, denn wenn ich den obenstehenden Code verwende und ich das bild wieder mit gcf aufrufe macht er eine neue Figure.
Wo hast du oben stehenden Code denn eingefügt? Vor oder nach dem Plot-Aufruf? Sollte davor geschehen
naja wenn ich das so lasse öffnet er ja den plot...
Einstellungen und Berechtigungen
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
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.
RSS