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

Performanceverbesserung Matrizen mit Bedingung
 

R0b1
Forum-Anfänger

Forum-Anfänger

Beiträge: 34
Anmeldedatum: 11.01.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 07.10.2014, 07:40     Titel: Performanceverbesserung Matrizen mit Bedingung
  Antworten mit Zitat      
Hallo,

ich möchte gern meine Berechnungsgeschwindigkeit verbessern.

Ausgangssituation: 3 Matrizen á <double 450x450>
- M <450x450> --> soll so groß sein wie Msoll
- U <450x450> --> soll kleiner sein als Umax
- I <450x450> --> soll kleiner sein als Imax und gesucht ist das Minimum von I (inkl. der Ort in der Matrix, wo der Wert steht)

Hier mal meine bisherigen Vorschläge:
Option 1:
Code:

% Durchsuchen jedes einzelnen Matrix-Ortes
for i=1:1:450
   for j=1:1:450
   
        % wenn alle Bedingungen erfüllt, dann Wert von I abspeichern
   if M(i,j) == Msoll && U(i,j) < Umax && I(i,j) < Imax
       I_merker(i,j) = I(i,j);
        % sonst NaN setzen
        else
            I_merker(i,j) = NaN;
   end
     
   end
end

% Abschließend das Minimum von I suchen. index_i und index_j geben den Matrix-Ort an
[I_min index_j] = min( min(I_merker) );
[I_min index_i] = min( I_merker(:,index_j) );
 


Option2:
Code:

% Matrix-Ort abspeichern, bei der die jeweilige Bedingung erfüllt ist
[ind_M]     = find( M == Msoll );
[ind_U]     = find( U < Umax );
% Matrix-Orte, bei denen Bedindungen von M und U erfüllt sind
ind_M_U     = intersect(ind_M, ind_U);

[ind_I]     = find( I < Imax );
% Matrix-Orte, bei denen Bedindungen von M, U und I erfüllt sind
ind_M_U_I   = intersect(ind_M_U, ind_I);

% Abschließend das Minimum von I suchen. Is_min_ort gibt den Matrix-Ort an
[I_min index]  = min( I(ind_M_U_I) );
Is_min_ort     = ind_M_U_I(index);
 


Habt ihr schnellere Lösungen für mich?

Grüße
Private Nachricht senden Benutzer-Profile anzeigen


R0b1
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 34
Anmeldedatum: 11.01.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 08.10.2014, 08:38     Titel:
  Antworten mit Zitat      
Zu kompliziert beschrieben oder einfach keine Idee?

Question Question Question
Private Nachricht senden Benutzer-Profile anzeigen
 
Andreas Goser
Forum-Meister

Forum-Meister

Beiträge: 3.654
Anmeldedatum: 04.12.08
Wohnort: Ismaning
Version: 1.0
     Beitrag Verfasst am: 08.10.2014, 09:21     Titel:
  Antworten mit Zitat      
Ich denke es liegt hauptsächlich daran, dass der Code ohne Daten nicht ausführbar ist und deshalb zuviel geraten werden muss.

Ich habe mal folgendes Beispiel geschrieben:

Code:

n=7;
M=rand(n, n)
M_soll=0.8*ones(n, n)
M>M_soll
 


Der Vergleich M>M_soll spuckt dann eine Matrix aus, die man natürlich auch weiter verarbeiten kann.

Andreas
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden
 
Harald
Forum-Meister

Forum-Meister

Beiträge: 24.499
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 08.10.2014, 09:32     Titel:
  Antworten mit Zitat      
Hallo,

wenn ich es richtig sehe, sollte folgendes für den ersten Block reichen:
Code:
I(M ~= Msoll & U < Umax & I < Imax) = NaN;


Die von Andreas und mir vorgeschlagene Technik heißt logische Indizierung (logical indexing).

Je nachdem wie oft die Bedingung M == Msoll erfüllt ist, kann es sein, dass deine Option 1 sogar schneller ist, wenn du I_merker auf NaN vorbelegst und den else-Teil weglässt.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
sara736
Forum-Newbie

Forum-Newbie

Beiträge: 1
Anmeldedatum: 09.10.14
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 09.10.2014, 08:51     Titel:
  Antworten mit Zitat      
leider kann ich dir ohne genauere Beschreibung (am besten in Formeln) des Problems nicht sagen, wie du vorgehen musst.
Grundsätzlich ist es für ode45 aber egal, ob die Größen durch Auswertung oder Interpolation gewonnen werden. Man muss sich lediglich darüber im klaren sein, dass der Interpolationsfehler bei den Toleranzen nicht berücksichtigt wird.
_________________

http://www.testking.us/70-332.htm
http://www.stlcop.edu/
http://www.hp.com/
http://www.cam.ac.uk/
Private Nachricht senden Benutzer-Profile anzeigen
 
R0b1
Themenstarter

 Forum-Anfänger

Forum-Anfänger

Beiträge: 34
Anmeldedatum: 11.01.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 09.10.2014, 10:40     Titel:
  Antworten mit Zitat      
Danke, die logische Indizierung kannte ich noch nicht.

War sehr hilfreich!
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 - 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.