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

For Schleife mit if-Bedingung im inneren, ersetzen?!

 

studi13

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 06.06.2013, 10:03     Titel: For Schleife mit if-Bedingung im inneren, ersetzen?!
  Antworten mit Zitat      
Hallo zusammen,

ich würde gerne eine for Schleife die im inneren eine if-Abfrage macht, ersetzen um das Program schneller zu machen.
Code:

for i=1:n
  for j=1:m
    if (B(i,j)>0)
      ...
    else
      ..
    end
  end
end
 


Habt ihr eine Idee, wie man die For Schleifen eliminieren kenn?


Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 06.06.2013, 10:07     Titel: Re: For Schleife mit if-Bedingung im inneren, ersetzen?!
  Antworten mit Zitat      
Hallo studi13,

Das kommt ganz darauf an, was sich hinter den "..." verbirgt. Kannst Du Details dazu posten?

Wahrscheinlich funktioniert dies dann gemäß:
Code:
index = (B > 0);
Result(index) = rand;
Result(~index) = -rand;


Gruß, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
studi13

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 06.06.2013, 10:13     Titel:
  Antworten mit Zitat      
Code:

      for i=1:m
         for j=1:n
            if (J_roh(2*i,2*j-1)<10^-14)
               w1(i,j,1) = 1;
                                        ...
            elseif(abs(J_roh(2*i,2*j)-J_roh(2*i-1,2*j-1))<10^-14)
               w1(i,j,1) = 1;
                                   ...
            else
               w1(i,j,1)= 2*J_roh(2*i-1,2*j)/sqrt((2*mu1(i,j)-2*J_roh(2*i-1,2*j-1))^2+4*J_roh(2*i-1,2*j).^2);
               ...
            end
         end
      end
 
 
Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 06.06.2013, 10:37     Titel:
  Antworten mit Zitat      
Hallo studi13,

Ja, nun wird es klarer und die Details sind für die Antwort wichtig.

10^-14 ist eine teurer Berechnung die in jeder Iteration ausgeführt wird. 1e-14 benötigt dagegen keine Kalkulationen.

Code:
w(:,:) = 2*J_roh(1:2:2*m - 1, 2:2:2*n) / ...
      sqrt((2*mu1 - 2*J_roh(1:2:2*m-1, 1:2:2*n-1)) .^ 2 + ...
      4 * J_roh(1:2:2*m-1, 2:2:2*n) .^ 2);

index = J_roh(2:2:2*m, 1:2:2*j-1) < 1e-14;
index = index | abs(J_roh(2:2:2*m, 2:2:2*n) - J_roh(1:2:2*m-1, 1:2:2*n-1)) < 1e-14;
w(index) = 1:
...

Das Prinzip sollte klar sein. Es würde mich nicht wundern, wenn ich bei dieser Index-Schlacht einen Tippfehler gemacht habe. Bitte prüfe das nochmal nach.

Gruß, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
studi13

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 06.06.2013, 14:39     Titel:
  Antworten mit Zitat      
Super, vielen vielen Dank!
 
studi13

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 14.06.2013, 12:03     Titel:
  Antworten mit Zitat      
Hallo zusammen,

ich habe nochmal eine Frage:
Ist es auch hier möglich die for-Schleife zu ersetzen:

Code:

for i=1:m
  for j=1:n
    A(i*(1:2),j*(1:2) = imfilter(A(i*(1:2),j*(1:2), maske, sigma)
  end
end
 


Dankeschön!
 
Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 17.06.2013, 00:38     Titel:
  Antworten mit Zitat      
Hallo studi13,

Das ist ein merkwürdiges Programm, da einerseits ein paar Klammern feheln, andererseits wird z.B. A(2,2) mehrfach berechnet und auch überschrieben. effizient ist das also nicht.

Gruß, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
studi13

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.06.2013, 07:22     Titel:
  Antworten mit Zitat      
Ja ...., war auch kompletter schwachsinn :/ Sorry. So jetzt nochmal richtig:
Code:

for i=1:m
  for j=1:n
    A((2*i-1)+(0:1),(2*j-1)+(0:1)) = imfilter((2*i-1)+(0:1),(2*j-1)+(0:1), maske, sigma)
  end
end
 


Ich hoffe das macht jetzt mehr Sinn. Aber sicherheitshalber nochmal in Worten. Die Matrix besteht aus m*n 2x2-Matrizen. Und diese sollen mit 'imfilter' bearbeitet und überschrieben werden.
 
studi13

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 19.06.2013, 10:07     Titel:
  Antworten mit Zitat      
Nicht möglich, oder unversändlich was ich möchte?
 
Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 19.06.2013, 23:17     Titel:
  Antworten mit Zitat      
Hallo studi13,

Ich zumindest verstehe nicht, was Du erreichen möchtest. Ich kenne IMFILTER mit den Inputs IMFILTER(A,H) und noch den Options. Du rufst es mit zwei Vektoren auf und die Werte von "maske" und "sigma" sind wahrscheinlich wichtig.

Zitat:
Die Matrix besteht aus m*n 2x2-Matrizen.

Wenn Du m*n 2x2 Matrizen hast, wäre ein [m x n x 2 x 2] Array nicht sinnvoller?

Gruß, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
studi13

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 20.06.2013, 14:44     Titel:
  Antworten mit Zitat      
Ja -.- mein Code war auch wieder falsch...

Code:

for i=1:m
  for j=1:n
    A((2*i-1)+(0:1),(2*j-1)+(0:1)) = imfilter(A((2*i-1)+(0:1),(2*j-1)+(0:1)), maske, sigma)
  end
end
 
 


Also A((2*i-1)+(0:1),(2*j-1)+(0:1))  \in \mathbb R ^2
Und für diese 2x2 Matrizen werden gefaltet.

Imoment ist mein A  \in \mathbb R ^{2m \times 2n}  .

Könnte es auch als [m x n x 2 x 2] Array schreiben. Wäre das ein Vorteil.

Ist das jetzt Verständlicher oder soll ich alles nochmal strukturiert erklären?!
 
studi13

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 20.06.2013, 15:53     Titel:
  Antworten mit Zitat      
Okay, klar ist es sinvoller es als [m x n x 2 x 2] array zu machen.
Weil dann kann ich ja wieder über die Indizes laufen lassen, z.B.
A(:,:,(1:2),(1:2)) ...

Ich glaube der Tipp hat schonmal geholfen. Ich werde mal rumprobieren!
Vielen Dank nochmal!
 
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.