[Τεχνικές Βελτιστοποίησης] Ασκήσεις και απορίες Εαρινό 2012

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

Για οποιοδήποτε πρόβλημα με register/login, στείλτε email στο contact@thmmy.gr.

[edenaxas]

ποιο Matlab εχεις?

πατα σε command
>>version
και κανε paste οτι σου πει.

Σε εμενα εβγαλε κανονικα αποτελεσμα 7 με matlab R2010a

εαν εχεις παλιοτερο κανε ενα κοπο και εγκατεστησε καινουριο

[τσαι-borg]

το R2012a εχω γαμω το κερατο του γαμω.
θα δοκιμασω να κατεβασω και το 2010 μπας και δουλεψει.

[pap-x]

Μπα μην το κάνεις...Έχω και τα 2 και τρέχει και στα 2. Μήπως προσπαθείς να χρησιμοποιήσεις hessian, diff ή gradient? Γιατί σε αυτά δεν τρέχει έτσι και βγάζει τέτοιο error.

[τσαι-borg]

ουτε προσθεση δεν εκανε σου λεω.
δεν ξερω τι εφταιγε,και το path κ τα παντα κοιταξα, το εκανα reinstall τελικα και τωρα ειναι οκ.

[pap-x]

Έχω κολλήσει στο τελευταίο (Newton). Κάνω το εξής για να βρω το hessian:

s=hessian(f);
ihess=inv(s);
A=subs(ihess, [x1 x2], [0 -1]);

gradient:

g=gradient(f);
B=subs(g, [x1 x2], [0 -1]);

Και μετά για το βήμα 4 κάνω:

x=x-a*A(1, : )*B;

αντίστοιχα για το y A(2, : ) κτλ...

Κάνω κάτι λάθος? Γιατί μου βγάζει ότι νάναι.

[foul frog]

φιλε, δοκιμασε να βρεις πρωτα το dk ξεχωριστα μονο του και μετα προχωρα.

εχω ενα αλλο προβλημα, δεν αποτυγχανει να βρει το ελαχιστο σε κανενα σημειο για την f. ειναι οκ αυτο?

[pap-x]

Μα έχει διαφορά το
dk=A(1, : )*B;
x=x-a*dk

από το
x=x-a*A(1, : )*B;

?
Για αυτό που ρώτησες νομίζω είναι οκ στην Newton αλλά όχι στη μέγιστη κάθοδο.

[edenaxas]

ξερεις τι γινεται?

εχεις πχ την δυδιαστατη f = x*y (δεν εχουμε αυτη αλλα ειναι μακαρονι η αλλη)

Στις μεθοδους σου λεει να βρεις το γ που να ελαχιστοποιει την f(xk + γ * dk) αρα πρεπει να βρεις το dk, εχεις το xk και εντελει καταληγεις να εχεις 1) μια παραμετρο που αθροιζεται με τη μια διασταση του διανυσματος - το xk(1) 2) μια παραμετρο που αθροιζεται με τη δευτερη διάσταση του διανύσματος - το xk(2) 3) την 1η διάσταση του dk και 4) τη 2η διάσταση του dk. Αυτα σου δινουν στο τελος ενα γραμμικο μετασχηματισμο με μοναδικη μεταβλητη το γ. δηλαδη καταληγεις να εχεις να ελαχιστοποιησεις οχι την f(x,y) = x*y αλλα την f(γ) = (α+γ*β) * (ε+γ*δ) που ειναι μια 1D τελειως διαφορετικη συναρτηση απο την f.
Το γ που θα βρεις τοτε το κανεις Χ(κ+1) = Χ(κ) - γ * dk και τελος.

δν ξερω εαν το εχεις κανει ετσι, απλα δεν μου φαινεται straightforward στο βιβλιο οτι στην ουσια πρεπει να κανεις για καθε συναρτηση και μια ακομη που να παιρνει και 4 ακομη εισοδους και να ειναι 1D.

[pap-x]

Έτσι είναι? Δεν το φαντάστηκα καθόλου. Στη μέγιστη κάθοδο το πήρα σταθερό το γ γιατί όταν το έβαζα να ψάχνει το min έβγαινε χειρότερο, εκτός του ότι έπαιρνε 1 ώρα να τελειώσει. Θα το κάνω όπως το λες στο Newton να δω αν θα καταφέρω κάτι. Thanks!

[Silvo the Beautiful]

για τα διαγραμματα τι πρεπει να παρω? Εχω μπερδευτει λιγο και θελω να ειμαι σιγουρος

[pap-x]

Για τα διαγράμματα παίρνεις τα α και β σε κάθε βήμα και τα εμφανίζεις σε άξονες. Για να καταλάβεις δες τι έβγαλα εγώ στη μέθοδο του χρυσού τομέα:

http://i44.tinypic.com/zitv1s.jpg

Για το Newton δεν κατάφερα τίποτα...Πρέπει να ελαχιστοποιήσω την f(xk-γκ*dk) αλλά εκτός του ότι καθυστερεί πολύ, παίρνει πάντα την μικρότερη τιμή του γκ (από εκεί που αρχίζουν οι επαναλήψεις (for g=0.0001:0.0001:0.01)) σαν ελάχιστη. Αν κάποιος το έχει κάνει και τρέχει ας γράψει περισσότερες λεπτομέρειες γιατί έχω χάσει τις μισές διακοπές του Πάσχα για αυτή την εργασία...

[Dimitris1989]

Καλησπέρα και χρονια πολλα,
βάζω στο matlab την εξίσωση της g(x)=power(x+4,0.5)+(power(x,2)-6)*log(x+1) και μου πετάει error "Inner matrix dimensions must agree."

καμιά ιδέα;

[τσαι-borg]

αν και δεν εχω ιδεα γιατι το γραφεις με power γενικα πρεπει να χρησιμοποεις τις πραξεις με την τελεια.
δηλαδη να γραψεις g=(x+4).^1/2 κλπ  κλπ

[Silvo the Beautiful]

Ερωτηση.... Στο ερωτημα 6. Τι πρεπει αν αλλαζω για να παρω τιμες? Αλλαξα το l και με αυτα τα xi που διαλεξα για σημεια με εδιναν το ιδιο αποτελεσμα. Τι θα παιρνουμε καθε φορα τυχαια xi? εχει νοημα αυτο?

Α και κατι ακομα... Πως βαζω απειρο πεδιο ορισμου. Στο ερωτημα 8 δεν μπορω να βγαλω τα 3d

[pap-x]

Κανονικά για διαφορετικά l έπρεπε να σου δίνει διαφορετικἐς τιμές. Τα xi δεν τα αλλάζεις είναι τα ίδια με τα προηγούμενα. Το θέμα είναι να πάρεις σωστό x2 στην αρχή.

Δεν χρειάζεται να τις δεις σε όλο το πεδίο ορισμού βάλε από -5 έως 5. Για βοήθεια δοκίμασε να τις κάνεις στο wolframalpha.com και να δεις που παίζουν περίπου.