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

Code läuft nicht - Hinweise

 

thomas0815
Forum-Century

Forum-Century


Beiträge: 153
Anmeldedatum: 13.05.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.05.2010, 16:13     Titel: Code läuft nicht - Hinweise
  Antworten mit Zitat      
Hallo,

meine Code mit Schleife läuft nicht, kann mir jemand einen Tipp gebene?

Code:
k=4
qk=[7,9,12,10]
bk=[1,1,1,1]
x1=0
x2=3
x3=4

for kk=k-1:-1:2
   z(kk,1)=max(z(kk+1,1)+x1-bk(kk,1),z(kk+1,1)+x2-bk(kk,1),z(kk+1,1)+x3-bk(kk,1));
   end


Danke und viele Grüße,
Thomas
Private Nachricht senden Benutzer-Profile anzeigen


Vito
Forum-Guru

Forum-Guru


Beiträge: 315
Anmeldedatum: 02.11.09
Wohnort: Stuttgart
Version: ---
     Beitrag Verfasst am: 17.05.2010, 16:21     Titel:
  Antworten mit Zitat      
Hi,

du verwendest in der Schleife z, wobei z unbekannt ist.
Was willst du erreichen oder machen?
Gruß,
Vito
Private Nachricht senden Benutzer-Profile anzeigen
 
thomas0815
Themenstarter

Forum-Century

Forum-Century


Beiträge: 153
Anmeldedatum: 13.05.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.05.2010, 16:36     Titel:
  Antworten mit Zitat      
Hallo,

nunja, ich möchte den Zustand z(kk) bestimmen und dieser hängt vom Zustand z(kk+1) ab. Also von seinem Vorgänger (Rückwärtsrekursion)!

Brauchst du detaillierte Infos; evtl. eine umfassendere Beschreibung oder ist das ausreichend?

Viele Grüße,
Thomas
Private Nachricht senden Benutzer-Profile anzeigen
 
Vito
Forum-Guru

Forum-Guru


Beiträge: 315
Anmeldedatum: 02.11.09
Wohnort: Stuttgart
Version: ---
     Beitrag Verfasst am: 17.05.2010, 17:13     Titel:
  Antworten mit Zitat      
Hi,

vll doch etwas mehr Info.
Ist das sowas wie der Kalman-Filter?

Gruß,
Vito
Private Nachricht senden Benutzer-Profile anzeigen
 
thomas0815
Themenstarter

Forum-Century

Forum-Century


Beiträge: 153
Anmeldedatum: 13.05.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.05.2010, 17:40     Titel:
  Antworten mit Zitat      
Hallo Vito,

ich habe das noch nicht mit dem kalman-filter in Verbindung gebracht, aber das worauf ich hinaus möchte könnte in Matlab vllt. ähnlich abgebildet werden. Es geht hier um eine Art Markow-Prozess. Der aktuelle Zustand darf nur von seinem Vorgänger nicht aber von der ganzen Geschichte der Zustände abhängen. Ein Beispiel:

wir befinden uns im Zustand 1 z1=5 (der vorhergehende Zustand war z.B. z0=3) ausgehend von diesem Zustand z1=5 kann nur eine Veränderung von +5 erfolgen. Entsprechend ergibt sich der Zustand 2 z2=z1+5 (NICHT: z2=z0+z1+5)

Also: Der aktuelle Zustand darf nur von seinem direkten Vorgänger abhängen, nicht von weiter zurückliegenden zuständen.

Ist es jetzt nachvollziehbarer?

Viele Grüße,
Thomas
Private Nachricht senden Benutzer-Profile anzeigen
 
Vito
Forum-Guru

Forum-Guru


Beiträge: 315
Anmeldedatum: 02.11.09
Wohnort: Stuttgart
Version: ---
     Beitrag Verfasst am: 17.05.2010, 18:18     Titel:
  Antworten mit Zitat      
Hi,

jetzt kommen wir doch der Sache näher Wink . Du könntest mal hier:

http://www.gomatlab.de/buchempfehlu.....11344,highlight,vito.html

reinschauen. Kalman-Filter Bsp mit Matlab. Oder hier

http://www.mathworks.com/matlabcent.....ltinominal-option-pricing

http://www.cs.ubc.ca/~murphyk/Software/HMM/hmm.html

Gruß,
Vito
Private Nachricht senden Benutzer-Profile anzeigen
 
denny
Supporter

Supporter



Beiträge: 3.853
Anmeldedatum: 14.02.08
Wohnort: Ulm
Version: R2012b
     Beitrag Verfasst am: 17.05.2010, 18:22     Titel:
  Antworten mit Zitat      
Hallo

wenn ich Kalmanfilter jetzt vernachlässige, und nehme an dass z einfach ein Array ist.
Dann sehe ich schon folgende Fehler:

Fehler 1: du erwartest in der FOR-Schleife immer ein Spaltenvektor.
In Beispiel hast du aber einen Zeilenvektor.

Fehler2: du übergibst an die MAX-Funktion 2 Werte, wird aber nur ein Wert akzeptiert.

nicht max([3 5 6], 2) sondern max([3 5 6 2])


Code:

k=4
qk=[7,9,12,10]'
bk=[1,1,1,1]'
x1=0
x2=3
x3=4

z= round(rand(k,1)*10)

for kk=k-1:-1:2
   z(kk,1)=max([z(kk+1,1)+x1-bk(kk,1),z(kk+1,1)+x2-bk(kk,1),z(kk+1,1)+x3-bk(kk,1)]);
end
 
Private Nachricht senden Benutzer-Profile anzeigen
 
thomas0815
Themenstarter

Forum-Century

Forum-Century


Beiträge: 153
Anmeldedatum: 13.05.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 17.05.2010, 19:14     Titel:
  Antworten mit Zitat      
Hallo Denny, Hallo Vito,

danke für die vielen nutzlichen Hinweise.

@denny ... welche Funktion erfüllt das round, da bin ich jetzt noch nicht so ganz dahinter gestiegen.

Gibt es jetzt auch noch die Möglichkeit mir ausgeben zu lassen, ob x1, x2 oder x3 die Funktion maximiert, um mit diesem Wert dann weiter zu rechnen?


Soweit erst mal vielen Dank, das hat mich ein gutes Stück nach vorn gebracht...

Viele Grüße,
Thomas
Private Nachricht senden Benutzer-Profile anzeigen
 
denny
Supporter

Supporter



Beiträge: 3.853
Anmeldedatum: 14.02.08
Wohnort: Ulm
Version: R2012b
     Beitrag Verfasst am: 17.05.2010, 20:17     Titel:
  Antworten mit Zitat      
Code:

z= round(rand(k,1)*10)
 

1) Das RAND erzeugt eine Zufallsvektor zwischen 0 und 1 (Kommazahlen)
2) Nach Multiplikation mit 10 habe ich Werte zwischen 0 und 10, aber immer noch Kommazahlen
3) Jetzt kommt ROUND ins Spiel und rundet die Werte.
Private Nachricht senden Benutzer-Profile anzeigen
 
denny
Supporter

Supporter



Beiträge: 3.853
Anmeldedatum: 14.02.08
Wohnort: Ulm
Version: R2012b
     Beitrag Verfasst am: 17.05.2010, 20:24     Titel:
  Antworten mit Zitat      
thomas0815 hat Folgendes geschrieben:


Gibt es jetzt auch noch die Möglichkeit mir ausgeben zu lassen, ob x1, x2 oder x3 die Funktion maximiert, um mit diesem Wert dann weiter zu rechnen?



Weiß ich nicht, ob ich richtig verstehe, aber MAX kann auch die Position ausgeben, wo Maximumwert aufgetreten ist

Code:

for kk=k-1:-1:2
   [z(kk,1), pos_max_wert]=max([z(kk+1,1)+x1-bk(kk,1),z(kk+1,1)+x2-bk(kk,1),z(kk+1,1)+x3-bk(kk,1)]);
end
 
Private Nachricht senden Benutzer-Profile anzeigen
 
thomas0815
Themenstarter

Forum-Century

Forum-Century


Beiträge: 153
Anmeldedatum: 13.05.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.05.2010, 11:33     Titel:
  Antworten mit Zitat      
Hallo Denny,

danke für die Hinweise. Ich werde jetzt damit erst einmal arbeiten und sehen, wann ich an die Grenzen stoße und dann ggf. wieder posten.

Vielen Dank und viele Grüße,
Thomas
Private Nachricht senden Benutzer-Profile anzeigen
 
thomas0815
Themenstarter

Forum-Century

Forum-Century


Beiträge: 153
Anmeldedatum: 13.05.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.05.2010, 12:05     Titel: schleife mit If-Anweisung läuft nicht
  Antworten mit Zitat      
... so, nun habe ich das Ganze mit einer If-Anweisung verknüpft und schon wieder läuft der Code nicht Sad bzw. gibt mir falsche Ergebnisse aus!

Code:
clear all
k=4
x1=0
x2=2
x3=-2

z=zeros(k+1,1)


for kk=k:-1:1
   if z(kk,1)>=5
   z(kk,1)=max([z(kk+1,1)+x1,z(kk+1,1)+x3])
   elseif z(kk,1)<=0
   z(kk,1)=max([z(kk+1,1)+x1,z(kk+1,1)+x2])
   else
   z(kk,1)=max([z(kk+1,1)+x1,z(kk+1,1)+x3,z(kk+1,1)+x2])
   end
   
 end


Aber wo liegt der Fehler, ich meine die Bedingen richtig formuliert zu haben. Ich will damit zum Ausdruck bringen, dass wenn z(kk,1) einen Wert von größer 5 annimmt, eine andere Berechnungsvorschrift für z(kk,1) gelten soll, damit im Ergebnis verhindert wird, das z(kk,1) größer 5 ist. Selbes gilt analog für den Umstand z(kk,1)<0.

Bin für Hilfe und Hinweise wie immer dankbar.

Viele Grüße,
Thomas
Private Nachricht senden Benutzer-Profile anzeigen
 
denny
Supporter

Supporter



Beiträge: 3.853
Anmeldedatum: 14.02.08
Wohnort: Ulm
Version: R2012b
     Beitrag Verfasst am: 18.05.2010, 16:18     Titel:
  Antworten mit Zitat      
Hallo

weil vielleicht du nicht auf aktuellen, sondern auf vorherigen Zustand beziehen sollst.

Code:

clear all
k  = 4
x1 = 0
x2 = 2
x3 =-2

z  = zeros(k+1,1)


for kk=k:-1:1
   if z(kk+1,1)>=5
      z(kk,1) = max([z(kk+1,1)+x1,z(kk+1,1)+x3])
   elseif z(kk+1,1)<=0
      z(kk,1) = max([z(kk+1,1)+x1,z(kk+1,1)+x2])
   else
      z(kk,1) = max([z(kk+1,1)+x1,z(kk+1,1)+x3,z(kk+1,1)+x2])
   end  
end
 
Private Nachricht senden Benutzer-Profile anzeigen
 
thomas0815
Themenstarter

Forum-Century

Forum-Century


Beiträge: 153
Anmeldedatum: 13.05.09
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 19.05.2010, 18:25     Titel:
  Antworten mit Zitat      
ja, danke Smile

soweit so gut, funktioniert und wieder bin ich einen Schritt weiter. Jetzt habe ich versucht, mir aus der Schleife den x-Wert geben zu lassen, der das ganze maximiert und die if Anweisung erfüllt. Aber das funktioniert ja nun nicht mehr so, wie es sich am Anfang ergeben hat. Und wenn ich in die for-Schleife am Ende [z(kk,1), pos_max_wert] eingebe, dann gibt mir Matlab einen Fehler aus Sad

Woran liegt das?
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: 19.05.2010, 20:54     Titel:
  Antworten mit Zitat      
Hallo Thomas,

ich klinke mich mal ein, da ich auch in einem anderen Thread von dir geschrieben habe.

[z(kk,1), pos_max_wert] wird die Fehlermeldung geben, dass pos_max_wert nicht definiert wurde. Das sollte nicht weiter verwundern, weil das ja auch davor nie vorkommt.

Zwei Bitten / Vorschläge:
- Tu dir eine systematische Einführung in MATLAB (Buch / Vorlesung / Kurs) an. Es besteht sonst die Gefahr, dass man Stückwerk produziert und sich nur von einem Fehler / Problem zum nächsten hangelt.
- Wenn du einen Fehler bekommst und Hilfe dazu möchtest, kopiere bitte den Fehler hier rein!

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.