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

Fehler beim plotten von Diagramm

 

kwd46056
Forum-Anfänger

Forum-Anfänger


Beiträge: 39
Anmeldedatum: 27.10.17
Wohnort: Regensburg
Version: ---
     Beitrag Verfasst am: 03.11.2017, 14:51     Titel: Fehler beim plotten von Diagramm
  Antworten mit Zitat      
Hey,

ich bins mal wieder.
Folgendes Problem: Möchte für verschiedene Frequenzen
Diagramme in 3 Spalten plotten lassen.
Hierfür habe ich zwei Funktionen geschrieben.
evalplotMessSim soll die Diagramme plotten
und die Funktion Formeln gibt die Formeln für die Berechnung
der Kennwerte vor.

Code:
function  [result, b] = evalPlotMessSim(...
                      strHall,...
                      parvars,...
                      freqs,...
                      file,...
                      simfile)
%Daten von Messungen einlesen
measData = readMessung(file);
%Daten von Simulation einlesen
simDatas  = readSimulation(simfile, freqs);
%für alle Frequenzen
allFreqs = sort(unique([measData.f, simDatas.f]));                   %sortieren und einmalig anzeigen aus Funktionen
       
%% Kraft-Weg (x-F) Plots mit Einrichtung der Seite
strParvar = strtrim(sprintf('% s', parvars{:}));
b(2) = figure(...
    'Name', sprintf('%s x-F %s', strHall, strParvar),...  
    'PaperType',    'A4', ...
                    'Units','centimeters', ...
                    'PaperUnits','centimeters');
pos   = get(b(2), 'Position');
psize = get(b(2), 'Papersize');
set(b(2), 'Position', [pos(1:2)+pos(3:4)-psize psize]);
set(0, 'defaultAxesFontSize', 6);

%Plots in drei Spalten verteilen
ncols = 3;                                                                 %Seite wird in 3 Spalten aufgteilt
nrows = ceil(length(allFreqs)/ncols);                                      %Anzahl der Reihen auf der Seite mit Obergrenze

icol  = 0;
irow   = 1;
result    = Formeln();

for iplot = 1:length(allFreqs)
    freq = allFreqs(iplot);
    icol = icol + 1;
    if icol>ncols, icol=1; irow=irow+1;
    end
   
    ax = subplot(nrows, ncols, iplot);
    hold (ax, 'on');
    grid (ax, 'on');
    set  (ax, 'XMinorGrid', 'on', 'YMinorGrid', 'on');
    claLegend = {};
% Messdaten plotten
        indData = find(freq == [measData.f]);
        if ~isempty(indData) &&  ~isempty(measData(indData).x)
            plot(ax, 1000*measData(indData).x, 0.001*measData(indData).F, 'x-');
            claLegend{end+1}  = 'meas';
            result(end+1) = Formeln(measData(indData), 'meas');
        end
%Simulationsdaten plotten
  for iSim = 1:length(parvars)
            indData   = find(freq == [simDatas(iSim).f]);
            if ~isempty(indData) &&  ~isempty(simDatas(iSim).x)
                plot(ax, 1000*simDatas(iSim).x, 0.001*simDatas(iSim).F);
                claLegend{end+1}    = parvars{iSim};
                result(end+1)       = Formeln(simDatas(iSim), parvars{iSim});
            end
         
 end


Code:
function chars = Formeln(data, ID)
%   Detailed explanation goes here
if nargin == 0
    chars = struct( ...
        'ID',       {}, ...                                                % gibt die Anzahl der im Aufruf
        'freq',     {}, ...                                                % übergebenen Eingabeargumente an die
        'xamp',   {}, ...                                                  % aktuelle ausgfeührte Funktion zurück
        'Famp',   {}, ...
        'dFx0',     {}, ...
        'dFend',    {}, ...
        'Aw',        {}, ...
        'delta',    {}, ...
        'cdyn',     {});
    return
end

chars       = {};                                                            
chars.ID    = ID;
chars.freq  = data.f;

%%Formeln zur Berechnung der Werte

%max. Verschiebung
chars.xamp = 0.5*(max(data.x) - min(data.x));

%max. Kraft
chars.Famp = 0.5*(max(data.F) - min(data.F));

%max. Kraft bei x=0
idxX0         = (data.x(:).*circshift(data.x(:), [-1 0]) <= 0);
chars.dFx0      = 0.5*abs(diff(data.F(idxX0)));

%max. Kraft bei min und max. x
    [~, indXmin]    = min(data.x);
    [~, indXmax]    = max(data.x);
    chars.dFend     = 0.5*abs(diff(data.F([indXmin, indXmax])));

%Verlustarbeit
chars.Aw = trapz(data.x, data.F);

%Verlustwinkel
chars.delta = asin(chars.Aw/(pi*chars.Famp*chars.xamp));

%dynamische Steifikeit
chars.cdyn = chars.Famp / chars.xamp;





end

 


Matlab gibt mir beim ausführen des Hauptprogramms folgende Fehlermeldung raus.

Error using permute
ORDER contains an invalid permutation index.

Error in trapz (line 42)
y = permute(y,perm);

Error in Formeln (line 39)
chars.Aw = trapz(data.x, data.F);

Error in evalPlotMessSim (line 50)
result(end+1) = Formeln(measData(indData), 'meas');

Versteh nicht ganz was der Fehler in line42 ist.
Hat jemand einen Tipp?

Gruss
kwd46056
Private Nachricht senden Benutzer-Profile anzeigen


SigiSorglos
Forum-Century

Forum-Century


Beiträge: 147
Anmeldedatum: 24.09.14
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 03.11.2017, 15:18     Titel:
  Antworten mit Zitat      
wenn man sich die Hilfe zu "permute(A,order)" anschaut, wird schnell klar, dass etwas mit der Variable "perm" nicht stimmt. Bzw. diese einen nicht gültigen Wert enthält. Da die Funktion von "trapz()" aufgerufen wird, solltest du dir die Eingabewerte von "trapz()" nochmal anschauen.

Wahrscheinlich stimmt was bei "data.x" oder "data.F" nicht.
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.