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

Schleife für Volumenberechnung

 

Lord
Forum-Anfänger

Forum-Anfänger


Beiträge: 16
Anmeldedatum: 20.10.08
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 20.10.2008, 12:23     Titel: Schleife für Volumenberechnung
  Antworten mit Zitat      
Hallo,
Ich hoffe es ist noch kein ähnliches Problem aufgetreten…
Folgende Vektoren habe ich:
x (34x1)
y (70x1)
z (70x1)

Jetzt will ich wie folgt rechnen: s = y – xi …also für jedes y – x1, y – x2 usw..
Nach jedem Schleifendurchlauf soll das Ergebnis s mit z und einer weiteren zahl (z.B. 24) multipliziert und anschließend der Vektor aufsummiert werden.
z.b. r = sum (s1 * z1 * 24, s2*z2*24…send*zend*24)…danach die nächste Schleife mit Subtraktion des nächsten x…
Der Ergbnisvektor sollte dann so groß wie x sein.

Hoffentlich ist das verständlich ausgedrückt. Tut mir leid, wenn es kein großes Problem ist, aber ich bin nicht der große Programmierkünstler und auch Matlab Anfänger.

Danke schonmal im vorraus!
Private Nachricht senden Benutzer-Profile anzeigen


Schrank
Forum-Century

Forum-Century


Beiträge: 203
Anmeldedatum: 10.07.08
Wohnort: Berlin
Version: ---
     Beitrag Verfasst am: 21.10.2008, 15:06     Titel:
  Antworten mit Zitat      
Hallo,
leider nicht verständlich ausgedrückt. Das klingt doch sehr nach einer Hausaufgabe... wäre es da nicht einfacher die Aufgabenstellung online zu stellen?
Gruß Schrank
Private Nachricht senden Benutzer-Profile anzeigen
 
Lord
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 16
Anmeldedatum: 20.10.08
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 21.10.2008, 19:36     Titel:
  Antworten mit Zitat      
Nein ist keine Hausaufgabe...
Ich dachte es ist verständlicher, wenn ich es so beschreiben.
Was ich aber im Grunde machen muss, ist das Volumen eines Sees zu berechnen.
Ich hab ein DGM (Rasterkarte mit Höhen). Außerdem eine Messreihe der Wasserstandsmessung bei einer bestimmten topologischen Höhe. Daraus soll ich jetzt das Volumen berechnen.
Das heißt es ergibt sich eine Differenz von der topologischen Höhe und des gemessenen Wasserstand bei der Wasserstandsmessung. Diese muss ich bei dem Rest der Höhen im Raster abziehen.
Anschließend muss ich, um auf das Volumen zu kommen diese einzelnen Höhe der Rasterfläche multiplizieren und die Summe berechnen

Die Formeln würden nach der Beschreibung meines ersten Posts so aussehen:
1. s1 = y1(topol. höhe) - x1 (differenzbetrag aus topolog. und wasserstand am ort der wasserstandmessung), s2 = y1 - x2, ...
2. r1 = s1 * 24(rasterfläche), r2 = s2*24
3. Summe aus r1, r2 usw

diese Schleife muss nochmals durchlaufen für den nächsten Differenzsbetrag aus der topolog. und wasserstand am ort der wasserstandmessung)

jetzt besser verständlich?
danke nochmals für die, die sich meinem Problem widmen.
Private Nachricht senden Benutzer-Profile anzeigen
 
Helmert
Forum-Century

Forum-Century


Beiträge: 121
Anmeldedatum: 04.09.08
Wohnort: Karlsruhe
Version: 2011a
     Beitrag Verfasst am: 22.10.2008, 08:04     Titel:
  Antworten mit Zitat      
Das Problem ließe sich wohl (wenn ich alles richtig verstanden habe) in einer Zeile lösen:
Code:

sum(y*x'*24)'
 

y*x' ist eine Matrix, deren Spalten die Vektoren s_i enthalten. Diese werden noch mit 24 multipliziert. Anschließend werden mit sum() die Spaltensummen gebiltet.

DIESE LÖSUNG IST FALSCH! (siehe unten)

Zuletzt bearbeitet von Helmert am 22.10.2008, 08:23, insgesamt einmal bearbeitet
Private Nachricht senden Benutzer-Profile anzeigen
 
Schrank
Forum-Century

Forum-Century


Beiträge: 203
Anmeldedatum: 10.07.08
Wohnort: Berlin
Version: ---
     Beitrag Verfasst am: 22.10.2008, 08:17     Titel:
  Antworten mit Zitat      
Hallo,
ich hab es jetzt so verstanden:
Code:

for ZaehlerX = 1: length(x)
    s(:,ZaehlerX) = y1(:) -x(ZaehlerX);
    r(:,,ZaehlerX) = s(:,,ZaehlerX) *24;
    Erg(,ZaehlerX) = sum(r(:,,ZaehlerX));
end
bzw. halt
Code:
for ZaehlerX = 1: length(x)
    Erg(,ZaehlerX) = sum((y1(:) -x(ZaehlerX))*24);
end


Wenn man den Kopf anschaltet geht das bestimmt auch ohne Schleife Smile

Gruß Schrank
Private Nachricht senden Benutzer-Profile anzeigen
 
Helmert
Forum-Century

Forum-Century


Beiträge: 121
Anmeldedatum: 04.09.08
Wohnort: Karlsruhe
Version: 2011a
     Beitrag Verfasst am: 22.10.2008, 08:23     Titel:
  Antworten mit Zitat      
stimmt, da habe ich in meinem obigen Beitrag einen Denkfehler gemacht. Hier meine Version ohne Schleife:
Code:

sum((kron(y,ones(1,34))-kron(x',ones(70,1)))*24)'
 
Private Nachricht senden Benutzer-Profile anzeigen
 
Schrank
Forum-Century

Forum-Century


Beiträge: 203
Anmeldedatum: 10.07.08
Wohnort: Berlin
Version: ---
     Beitrag Verfasst am: 22.10.2008, 08:32     Titel:
  Antworten mit Zitat      
ui,
Code:
kannte ich noch nicht.

Ich frag mich nur, was z aus dem ersten Beitrag ist?

@Lord: aussagekräftigen Variablennamen helfen beim verstehen von Programmen ungemein Smile
Private Nachricht senden Benutzer-Profile anzeigen
 
Lord
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 16
Anmeldedatum: 20.10.08
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 22.10.2008, 11:06     Titel:
  Antworten mit Zitat      
Vielen Dank schonmal für die Vorschläge! Werde ich dann direkt mal ausprobieren...

Ja Pardon, werde beim nächsten mal die Variablen besser oder überhaupt deklarieren.
z sind die Anzahl der Rasterflächen für eine topologische Höhe. Das hatte ich in meinem 2ten Post vergessen. Damit muss das berechnete Volumen einfach nochmal multipliziert werden. Also nicht nur mit 24 sonderen mit dem Vektor z auch noch.
Anschließend kann die Summe berechnet werden.
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.