[Τεχνικές Βελτιστοποίησης] Ασκήσεις και απορίες Εαρινό 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 με προσοχή στα ονόματα των αρχείων!

[pap-x]

Σόρρυ που άργησα να απαντήσω.
Λοιπόν, το κύριο πρόβλημα είναι η εύρεση του γ. Στην αρχή το έβαλα να το κάνει κάπως έτσι:

for i=0.0001:0.0002:0.01
     if  fun(x+i*grad(f(...))) < min
         min= fun(x+i*grad(f(...)));
         g=i;
end

Απλά έτσι εκτός το ότι δεν ξέρεις ποιές τιμές και πόσες να βάλεις στο γ (από πού να αρχίζει και που να τελειώνει), μου έπαιρνε πάρα πολύ χρόνο και δεν κατέληγε σε κάτι καλό. Οπότε στη μέγιστη κάθοδο έβαλα γ σταθερό στο 0.001 και πήγε σχετικά καλά. Έκανε πολλά βήματα για να βρει κάτι αλλά το έβρισκε.
Στη Newton όμως δεν τρέχει ούτε έτσι ούτε αλλιώς. Πώς μπορεί να βγει το γ σωστά και στις 2 περιπτώσεις? 

[pap-x]

Λάθος ότι έγραψα πριν. Βρήκα τον σωστό τρόπο για να βρίσκω τα γ. Στη δεύτερη συνάρτηση βγάζει ότι νάναι γιατί πέφτει σε σημείο που η κλίση είναι 0 και δεν ξέρει που να πάει.
Τώρα στη Newton μάλλον έχω μπερδέψει τους πίνακες και παίρνω λάθος αποτέλεσμα.
Με την σειρά:

ορίζουμε την f

hess=hessian(f)
result=subs(hess, [x1 x2], [x y])
result=inv(result)

grad=gradient(f)
result=subs(grad, [x1 x2], [x y])

Οπότε καταλήγουμε σε έναν πίνακα 2*2 εσσιανό και έναν 2*1 της κλίσης. Αν τους πολλαπλασιάσουμε παίρνουμε έναν πίνακα 2*1 στον οποίο το πρώτο στοιχείο είναι αυτό που θα αφαιρεθεί από το χ1 και το δεύτερο από το χ2 για το επόμενο βήμα?

[fcdimitr]

Παίδες, δούλεψε κανείς την τρίτη εργασία? Κι αν ναι, έχετε καμιά ιδέα πως να υπολογίσουμε το βήμα γ? Εγώ ξεκίνησα θεωρώντας το βήμα ίσο με 1 και στη συνέχεια το μείωνα μέχρι να ικανοποιούνται τα κριτήρια 3 και 4.... Έχει κανείς καμιά καλύτερη ιδέα?

[edenaxas]

ετσι νομιζω πως ειναι.


καμμια παραταση θα γινει?

[Dimitris1989]

Tο dk πώς θα το υπολογίσουμε;

[τσαι-borg]

για το dk φανταζομαι πως θελουμε να κανουμε αυτο που λεει στην 138 στο βιβλιο στη σχεση 5.2.37.
οχι οτι καταλαβαινω τι λεει, αλλα δε γαμεις.
το μ_κ κατα ποσο μεγαλυτερο της μεγαλυτερης ιδιοτιμης το πηρατε?

[τσαι-borg]

μηπως ειπωθηκε κατι για παραταση στο μαθημα?

[crystal]

Ναι , θα πάρει αλλά δεν είπε μέχρι πότε

[τσαι-borg]

χελ γεα.

[Silvo the Beautiful]

μεχρι κυριακη σιγουρα!

[ectoras]

μέχρι 9 δείτε το ethmmy

[pap-x]

Το βήμα γ το ξεκινάμε από 1 και μειώνουμε σταδιακά. Εγώ για α έβαλα 0.0001 και για β=0.5 και έτρεξε κανονικά.
Το dk το υπολογίζουμε όπως πριν απλά στην θέση του εσσιανού βάζουμε το [hess(f)+μ*I].
Το μ το πήρα μεγαλύτερο κατά 1. Δεν έχει ιδιαίτερη σημασία αρκεί να είναι μεγαλύτερο.

[τσαι-borg]

επειδη μαλλον μειωνει το iq μου αυτο το μαθημα, οταν λεει με 3 διαφορετικα παραμενοντα βηματα τι στο καλο εννοει?
επισης αν καποιος μπορει να βρει καποιο λαθος σε αυτο , ειναι ευπροσδεκτος να το κανει
d'*subs(d22,[x y],x0+g*d) > b*d'*subs(d22,[x y],x0)  &&  subs(fun,[x y],x0+g*d)<= subs(fun,[x y],x0)+a*g*d'*subs(d22,[x y],x0)

d22 ειναι το gradF και d11 o εσσιανος.

[edenaxas]

το γ δεν πρεπει να το μειωνεις σιγα σιγα?

ε, δεν θα το μειωσεις και τελειως (0) θα το φτασεις σε ενα σημειο απο οπου και δεν θα μειωνεται αλλο. Ε αυτο ειναι.

πχ εαν εχουμε παραμενων 0,5 θα μπορουσε το γ να ειναι

10-8-5-3-2-1.7-0.7-0.5-0.5-0.5 κτλ

ή εαν εχεις παραμενων 0,2
10-8-5-3-2-1.7-0.7-0.5-0.2-0.2 you get the point..

[τσαι-borg]

το γ το μειωνω κατα 0.1 σε ενα for απλως δεν εκανα copy paste εκεινο το σημειο γιατι ο κωδικας κολλαει οταν παει να μπει στο if με τα κριτηρια(αυτο που εκανα copy paste).
ρωταω αν βλεπει κανεις κανενα λαθος γιατι δεν κανει ποτε break και δεν μπορω να βρω τι φταιει.

αυτο δλδ
for g=1:-0.01:0
   if d'*subs(d22,[x y],x0+g*d) > b*d'*subs(d22,[x y],x0) && subs(fun,[x y],x0+g*d)<= subs(fun,[x y],x0)+a*g*d'*subs(d22,[x y],x0)
        break;
   end
  end