For-Schelife mit find() zu langsam

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Winkow

Moderator


	Beiträge: 3.842

	Anmeldedatum: 04.11.11

	Wohnort: Dresden

	Version: R2014a 2015a

Verfasst am: 25.05.2015, 19:22 Titel:

Zitat:

1. Wie kriege ich denn in dem Beispiel FP mehrdimensional?

eine möglichkeit wäre es zb mittels einer schleife über die kleinere dimmension zu arbeiten. das wäre in deinem fall wieder die 31. dann hast du zwar eine schleife aber es wird ja trotzdem nicht 86400 mal find aufgerufen.

Code:

Preise=randi(10,86400,31);
[Ascend,Iasc] = sort(Preise,1,'ascend');
Descend=flipud(Ascend);
B=Descend-Ascend;
for k=1:31;
FB(Iasc(B(:,k)>0,k),k)=1;
FB(Iasc(B(:,k)<0,k),k)=-1;
end

Funktion ohne Link?

Zitat:

2. In dem Fall, dass durch die Wirkungsgradverluste Ascend > Descend wird, soll eine 0 in den Fahrplan.

kannst du das nochmal genauer erklären ?
_________________

richtig Fragen

Zuletzt bearbeitet von Winkow am 25.05.2015, 19:33, insgesamt einmal bearbeitet

Apocalypto

Themenstarter

Forum-Anfänger


	Beiträge: 23

	Anmeldedatum: 25.03.15

	Wohnort: ---

	Version: ---

Verfasst am: 25.05.2015, 19:29 Titel:

Nehmen wir mal ein Beispielpreispaar von vorhin:

9,8 und 8,5

in diesem Fall ist es ja so, dass

(9,8 * 0,9* 0,97) - (8,5 * (1/(0,9*0,97)) = -1,18

In diesem Fall würde ich keinen Gewinn mit diesem Preispaar machen, weshalb für den Index des jeweiligen Preises dann eine 0 im Fahrplan auftauchen muss

Winkow

Moderator


	Beiträge: 3.842

	Anmeldedatum: 04.11.11

	Wohnort: Dresden

	Version: R2014a 2015a

Verfasst am: 25.05.2015, 19:38 Titel:

hmm so?

Code:

Preise=40*rand(86400,31);
[Ascend,Iasc] = sort(Preise,1,'ascend');
Descend=flipud(Ascend);
B=Descend-Ascend;
eta_bat_B = 0.9;
eta_wr_B = 0.97;
P = Descend*eta_bat_B*eta_wr_B - Ascend*(1/(eta_bat_B*eta_wr_B));
for k=1:31;
FB(Iasc(B(:,k)>0,k),k)=1;
FB(Iasc(B(:,k)<0,k),k)=-1;
FB(Iasc(P(:,k)<0,k),k)=0;
end

Funktion ohne Link?

?
_________________

richtig Fragen

Apocalypto

Themenstarter

Forum-Anfänger


	Beiträge: 23

	Anmeldedatum: 25.03.15

	Wohnort: ---

	Version: ---

Verfasst am: 25.05.2015, 20:08 Titel:

Anhanhme Wirkungsgrad = 0,9

Ich versuch es nochmal an einem einfachen Beispiel:

Stunde 1 - 8,5
Stunde 2 - 9,3
Stunde 3 - 10,8
Stunde 4 - 6,5
Stunde 5 - 7,5
Stunde 6 - 12

Aufsteigend sortiert:

6,5
7,5
8,5
9,3
10,8
12

Absteigend sortiert:

12
10,8
9,3
8,5
7,5
6,5

Preispaare bilden (gegenüberstellen von absteigend und aufsteigend)

12*0,9 - 6,5*(1/0,9) = 3,5
10,8*0,9 - 7,5*(1/0,9) = 0,6
9,3*0,9 - 8,5*(1/0,9) = -1,07

Folglich würde in den Stunden 1 und 2 nicht passieren, Stunde 3 und 6 entladen und Stunde 4 und 5 geladen.

Ich benötige diese Ruhephasen der Batterie um diese nicht unnötig einer Belastung auszusetzen.

Die Variable B aus deinem Beispiel Frage ich deshalb so garnicht ab.

Winkow

Moderator


	Beiträge: 3.842

	Anmeldedatum: 04.11.11

	Wohnort: Dresden

	Version: R2014a 2015a

Verfasst am: 25.05.2015, 20:34 Titel:

Code:

Preise=[8.5,9.3,10.8,6.5,7.5,12]';
[Ascend,Iasc] = sort(Preise,1,'ascend');
[Descend,Ides] = sort(Preise,1,'descend');
eta_bat_B = 0.9;
eta_wr_B = 0.97;
P = Descend*eta_bat_B*eta_wr_B - Ascend*(1/(eta_bat_B*eta_wr_B));
FB=zeros(size(Preise));
for k=1:size(Preise,2);
FB(Iasc(P(:,k)>0,k),k)=1;
FB(Ides(P(:,k)>0,k),k)=-1;
end

Funktion ohne Link?

so?
_________________

richtig Fragen

Apocalypto

Themenstarter

Forum-Anfänger


	Beiträge: 23

	Anmeldedatum: 25.03.15

	Wohnort: ---

	Version: ---

Verfasst am: 25.05.2015, 21:01 Titel:

Genius!

Elapsed time is 0.171381 seconds.

Ich glaube es hackt! Ich danke dir vielmals!

Ich sollte dir auch mal meine anderen Codes zur verfügung stellen Very Happy

Danke Winkow! Hast mir sehr geholfen!

Gruß

Apo

Winkow

Moderator


	Beiträge: 3.842

	Anmeldedatum: 04.11.11

	Wohnort: Dresden

	Version: R2014a 2015a

Verfasst am: 25.05.2015, 22:04 Titel:

wenn du noch den tip von harald mit flipud einbaust kannst du nochmal ein bisschen spaaren.

Code:

[Ascend,Iasc] = sort(Preise,1,'ascend');
Descend = flipud(Ascend);
Ides=flipud(Iasc);
eta_bat_B = 0.9;
eta_wr_B = 0.97;
P = Descend*eta_bat_B*eta_wr_B - Ascend*(1/(eta_bat_B*eta_wr_B));
FB=zeros(size(Preise));
for k=1:size(Preise,2);
FB(Iasc(P(:,k)>0,k),k)=1;
FB(Ides(P(:,k)>0,k),k)=-1;
end

Funktion ohne Link?

grüße
edit winkow: da hatte sich der fehler teufel eingeschlichen.
_________________

richtig Fragen

Zuletzt bearbeitet von Winkow am 25.05.2015, 22:10, insgesamt einmal bearbeitet

Apocalypto

Themenstarter

Forum-Anfänger


	Beiträge: 23

	Anmeldedatum: 25.03.15

	Wohnort: ---

	Version: ---

Verfasst am: 25.05.2015, 22:08 Titel:

Das werde ich machen, danke!

	Gehe zu Seite Zurück 1, 2
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.