Programmiertipps für eigene Plotfunktion

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

DSP

Forum-Meister


	Beiträge: 2.117

	Anmeldedatum: 28.02.11

	Wohnort: ---

	Version: R2014b

Verfasst am: 03.10.2015, 11:15 Titel: Programmiertipps für eigene Plotfunktion

Hallo goMatlab Forum,

ich bin gerade dabei einen eigenen SignalViewer ähnlich dem Matlab Spectrum Viewer zu programmieren. Folgende Plots sollen dabei dargestellt werden können:

- Zeitbereich
- FFT
- Spektrogramm
- Leistungsdichtespektrum nach Welch

Es gibt insgesamt 4 Signale:

- Referenzsignal mit/ohne Fensterfunktion
- Compare Signal mit/ohne Fensterfunktion
- Testsignal mit/ohne Fensterfunktion
- Fensterfunktion

Gleichzeitig sollen nur Referenz- und Comparesignal dargestellt werden. Dazu benötige ich natürlich eine eigene Plot Funktion, welche ich momentan so angefangen habe:

Code:

function Plot_function(axh, PlotTyp, SignalTyp, Data, MarkersState, handlesMarkerStrings)
%
% Function Inputs:
%
% axh = axes handle
% PlotTyp = Darstellungsform -> Zeit-/Frequenzbereich
% SignalTyp = [xxxxxx] -> x = 0/1
% Data = Signals, Ts, Fs, etc
% MarkerState = on/off
% handlesMarkerStrings = handles for text fields at marker buttongroup
%
% -------------------------------------------------------------------------

switch(PlotTyp)
case 'time'
% load right signal
switch(SignalTyp)
case 000001
% Signal - SignalTyp: 000001
signal(:,1) = Data.Signal;
case 000011
% Signal * win - SignalTyp: 000011
signal(:,1) = Data.Signal .* Data.WindowVal;
case 000100
% Signal + CompareSignal - SignalTyp: 000100
signal(:,1) = Data.Signal;
signal(:,2) = Data.CompareSignal;
case 001100
% Signal * win + CompareSignal * win - SignalTyp: 001100
signal(:,1) = Data.Signal .* Data.WindowVal;
signal(:,2) = Data.CompareSignal .* Data.WindowVal;
case 010000
% TestSignal - SignalTyp: 010000
signal(:,1) = Data.TestSignal;
case 110000
% TestSignal * win - SignalTyp: 110000
signal(:,1) = Data.TestSignal .* Data.WindowVal;
end

% Abtastrate vorhanden?
if isnan(Data.Ts)
xVal = 1:Data.NumVal;
if size(signal,2) == 1
hplot = plot(axh,xVal,signal);
else
hplot = plot(axh,xVal,signal(:,1),xVal,signal(:,2));
end
xlim([1 Data.NumVal])
xTicks = round(linspace(1,Data.NumVal,10),2);
set(axh,'XTick',xTicks);
xlabel('Samples','FontSize',12,'FontWeight','bold','Color','k')
else
t = 0:Data.Ts:(Data.NumVal - 1) * Data.Ts;
if size(signal,2) == 1
hplot = plot(axh,t,signal);
else
hplot = plot(axh,t,signal(:,1),t,signal(:,2));
end
xlim([0 t(end)])
xTicks = round(linspace(0,t(end),10),2);
set(axh,'XTick',xTicks);
xlabel('Zeit [s]','FontSize',12,'FontWeight','bold','Color','k')
end
set(hplot,'LineWidth',1.2,...
'LineStyle','-',...
'Color','b',...
'Marker','.',...
'MarkerEdgeColor','k',...
'MarkerFaceColor',[0,0,0],...
'MarkerSize',4);
title('Zeitbereich','FontSize',12,'FontWeight','bold','Color','k')
ylabel('Amplitude','FontSize',12,'FontWeight','bold','Color','k')

case 'fft'
% Magnitude:
% Signal - SignalTyp: 000001
% Signal + CompareSignal - SignalTyp: 000011
plotyy(axh)

% Phase:
% Signal - SignalTyp: 000100
% Signal + CompareSignal - SignalTyp: 001100
plotyy(axh)

% Magnitude + Phase:
% Signal - SignalTyp: 010000
% Signal + CompareSignal - SignalTyp: 110000
plotyy(axh)
case 'pwelch'
case 'spec'
% 2D
% 3D
case 'win'
if isnan(Data.Ts)
hplot = plot(axh,1:Data.NumVal,Data.WindowVal);
xlabel('Samples','FontSize',12,'FontWeight','bold','Color','k')
xlim([1 Data.NumVal])
% ylim([round(min(Data.WindowVal),1)*1.1 round(max(Data.WindowVal)*1.1,1)])
xTicks = round(linspace(1,Data.NumVal,10),2);
% %xticklabel('')
% %yTicks = round(linspace(0,max(Data.WindowVal)*1.1,10),1);
% %set(axh,'XTick',xTicks,'YTick',yTicks);
set(axh,'XTick',xTicks);
% %yticklabel('')
else
t = 0:Data.Ts:(Data.NumVal - 1) * Data.Ts;
hplot = plot(axh,t,Data.WindowVal);
xlabel('Zeit [s]')
xlim([0 t(end)])
% ylim([round(min(Data.WindowVal),1)*1.1 round(max(Data.WindowVal)*1.1,1)])
xTicks = round(linspace(0,t(end),10),2);
% %xticklabel('')
% %yTicks = round(linspace(0,max(Data.WindowVal)*1.1,10),1);
% %set(axh,'XTick',xTicks,'YTick',yTicks);
set(axh,'XTick',xTicks);
%yticklabel('')
end
set(hplot,'LineWidth',1,...
'LineStyle','-',...
'Color','b',...
'Marker','.',...
'MarkerEdgeColor','k',...
'MarkerFaceColor',[0,0,0],...
'MarkerSize',4);
title('Zeitbereich','FontSize',12,'FontWeight','bold','Color','k')
ylabel('Amplitude','FontSize',12,'FontWeight','bold','Color','k')

otherwise
% empty plot
cla(axh,'reset')
dualcursor off
end

% load dualcursors new?
if strcmp(MarkersState,'on')
dualcursor off
dualcursor('on',[], [], [], axh, handlesMarkerStrings);
end
end

Funktion ohne Link?

In der GUI wird nun festgelegt welcher PlotTyp und welche(s) Signal(e) dargestellt werden sollen. Evtl. hat ja jemand von euch eine Idee für eine geeignetere Umsetzung. Ich befürchte nämlich, dass die Funktion mit steigender Anzahl von Fällen immer unübersichtlicher wird.

Ich freue mich über eure Ratschläge!

Gruß DSP


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 |

RSS

Hosted by:

Copyright © 2007 - 2025 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.