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

Effizienter Ersatz von Schleife
 

Bahnwärter
Forum-Newbie

Forum-Newbie

Beiträge: 2
Anmeldedatum: 03.06.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 03.06.2015, 09:17     Titel: Effizienter Ersatz von Schleife
  Antworten mit Zitat      
Hallo,
ich lese schon seit einigen Jahren intensiv im Forum und habe hierbei bisher immer viele hilfreiche Anregungen mitnehmen könnnen. Nun habe ich mich aktiv registriert und hoffe in Zukunft auch anderen etwas weiterhelfen zu können

Aktuell grüble ich an der Optimierung eines Codabschnitts und komme dabei nicht so recht weiter, vielleicht hat jemand ja eine Idee wie man dies Optimieren könnte


Code:
for i=1:n
    P(i)=CD(:,i)'*R*CD(:,i)./(CD(:,i)'*CD(:,i));
end



Als ergänzende Erklärung:

CD ist eine Matrix der Dimension n x m
R ist eine Matrix der Dimension n x n

Meine bisherigen Versuche über äußere Produkte sind gescheitert.

Würde mich über Anregungen freuen..

Grüße Bahnwärter
edit winkow: bbcode mus aktiviert sein damit die code umgebung funktioniert und der stich muss richtig rum :).
Private Nachricht senden Benutzer-Profile anzeigen


Harald
Forum-Meister

Forum-Meister

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

die Schleife sollte bis m laufen, oder?

Wie groß sind die Matrizen denn größenordnungsmäßig? Insbesondere ist interessant, wie das Größenverhältnis von m und n zueinander ist.

Code:
P2 = diag( CD'*R*CD./(CD'*CD) );

war in einem ersten Test schneller, obwohl viele überflüssige Berechnungen stattfinden.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Bahnwärter
Themenstarter

 Forum-Newbie

Forum-Newbie

Beiträge: 2
Anmeldedatum: 03.06.15
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 03.06.2015, 13:01     Titel:
  Antworten mit Zitat      
Hallo,
danke für die schnelle Antwort.
Richtig, die Schleife läuft bis m.

Von der Größenordnung liegt n im Bereich von 100-500 und m im Bereich von 10k>m>50k.

Bei meinem gegebene Zahlen lande ich bei deiner Implemtierung bei:

Out of memory. Type HELP MEMORY for your options.

Grüße
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister

Beiträge: 24.501
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 03.06.2015, 13:25     Titel:
  Antworten mit Zitat      
Hallo,

ich hab mal damit rumgespielt.

Code:
n = 500;
m = 50000;
CD = rand(n, m);
R = rand(n,n);


tic
P = zeros(1,m);
for i=1:m
P(i)=CD(:,i)'*R*CD(:,i)./(CD(:,i)'*CD(:,i));
end
toc

% 1. Verbesserungsversuch
tic
CDT = CD';
CDR = CDT*R;
P2 = sum(CDR .* CDT, 2);
DEN = sum(CDT .* CDT, 2);
P2 = (P2 ./ DEN)';
toc
max(abs(P - P2))

% Weiterer Versuch
tic
CDR = R'*CD;
P3 = sum(CDR .* CD);
DEN = sum(CD .* CD);
P3 = P3 ./ DEN;
toc
max(abs(P - P3))


Grüße,
Harald
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.