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

[Admins]

Για αλλαγή του public name σας, επικοινωνήστε με έναν από τους Admins.

[edenaxas]

Προσπαθεις στην λαθος κατευθυνση νομιζω.

Δεν χρειαζεται να κανεις explicitly ενα πινακα τιμων και να λυσεις με το χερι.

πχ εχεις την g(x) = x2+4x+2  εγω θα εκανα

syms x;
syms g;

g = x^2 + 4*x + 2;
y = solve(g);

και εχεις στη μεταβλητη y τις ριζες σε συμβολικη μορφη.

[pap-x]

Ευχαριστώ πολύ για την βοήθεια!
Εξηγώ πώς έλυσα την 7η άσκηση για να μου πείτε αν το παίρνω λάθος.
Για να βρώ το ελάχιστο πήρα την συνάρτηση σαν fx=@(x)3*x^4...... και έκανα fminsearch.
Για το μέγιστο έκανα *(-1) την συνάρτηση (αναποδογύρισα) και έκανα το ίδιο (δεν υπάρχει fmaxsearch).
Για τα σημεία καμπής έκανα syms f, μετά diff(f) και μετά βρήκα τις ρίζες μέσω solve (ρίζες παραγώγου).
Για σημεία που τέμνει τον χ έκανα solve την f.
Και για plot έκανα x=-100:0.1:100 και f=3*x^4...

Σωστή είναι η διαδικασία?

[pap-x]

Επίσης δεν μου τρέχει τα παραδείγματα που έχει στο gradient και hessian και βγάζει error. Δεν θα έπρεπε να τα τρέχει κανονικά? Μήπως θέλει κάποιο επιπλέον πακέτο?

Εrror στο gradient:

syms x y z
f = 2*y*z*sin(x) + 3*x*sin(z)*cos(y);
gradient(f, [x, y, z])
??? Error using ==> zeros
Trailing string input must be a valid numeric class name.

Error in ==> gradient at 64
   g  = zeros(size(f),class(f)); % case of singleton dimension

Error στο hessian:
syms x y z
f = x*y + 2*z*x;
hessian(f)
??? Undefined function or method 'hessian' for input arguments of type 'sym'.

[Emfanever]

Επίσης δεν μου τρέχει τα παραδείγματα που έχει στο gradient και hessian και βγάζει error. Δεν θα έπρεπε να τα τρέχει κανονικά? Μήπως θέλει κάποιο επιπλέον πακέτο?

Εrror στο gradient:

syms x y z
f = 2*y*z*sin(x) + 3*x*sin(z)*cos(y);
gradient(f, [x, y, z])
??? Error using ==> zeros
Trailing string input must be a valid numeric class name.

Error in ==> gradient at 64
   g  = zeros(size(f),class(f)); % case of singleton dimension

Error στο hessian:
syms x y z
f = x*y + 2*z*x;
hessian(f)
??? Undefined function or method 'hessian' for input arguments of type 'sym'.

Δε λειτουργούν έτσι οι συναρτήσεις.

Για το gradient θα το βρεις manually

grad = [diff(f,x) diff(f,y) diff(f,z) ]

Το hessian το βρίσκεις με πατέντα

hess = jacobian(jacobian(f))

edit: Βασικά τώρα που είδα το documentation, έτσι το έχει όπως τα έγραψες, αλλά ίσως δουλεύει μόνο σε matlab 2012. Από r2011b και πάνω τα έχει.

[crystal]

Eίπαν τίποτα για τη 2η εργασία στο μάθημα?

[edenaxas]

πηγε κανενας?

[Room]

πηγε κανενας?

[τσαι-borg]

Η Κωσταριγκα ειναι "υπευθυνη" για τις εργασιες οποτε αν ειπωθει κατι θα γινει σημερα.

[Silvo the Beautiful]

Η Κωσταριγκα ειναι "υπευθυνη" για τις εργασιες οποτε αν ειπωθει κατι θα γινει σημερα.

ααααχ 

[edenaxas]

απο ethmmy

Άσκηση 2

Παράδοση εργασίας μέχρι Τετάρτη 25/04/2012. Η αναφορά (αρχείο .word, .tex) να έχει όνομα με λατινικούς χαρακτήρες και να είναι σε μορφή .pfd κατά προτίμηση.

[~GiA~]

μαρεσε το καλο πασχα που εχει στο τελος...
στα διαλα  μας εχουν φαει τις διακοπες παλι .. 

[pap-x]

Ρε παιδιά εγώ δεν βγάζω τίποτα από την εκφώνηση. Είναι πολύ γενική.
1. Για να λειτουργεί η συνάρτηση που θα φτιάξουμε για κάθε f(x) πρέπει να παίρνει σαν όρισμα την συνάρτηση που θα χρησιμοποιεί. Πώς γίνεται στο matlab να κάνουμε συνάρτηση που θα παίρνει σαν όρισμα άλλη συνάρτηση?
2. Η συνάρτηση τι θα παίρνει και τι θα δίνει? Στο πρώτο ερώτημα λέει να δίνει α και β και κ για διάφορα ε, ενώ στο δεύτερο λέει να δίνει τα α, β για διάφορα κ.
3. Λέει να επιλέξουμε τις βέλτιστες περιπτώσεις αλλά με βάση τι? Γιατί έχουμε σαν παράγοντες και το α,β αλλά και το κ. Μας ενδιαφέρουν τα βήματα ή το τελικό διάστημα?

[edenaxas]

Ρε παιδιά εγώ δεν βγάζω τίποτα από την εκφώνηση. Είναι πολύ γενική.
1. Για να λειτουργεί η συνάρτηση που θα φτιάξουμε για κάθε f(x) πρέπει να παίρνει σαν όρισμα την συνάρτηση που θα χρησιμοποιεί. Πώς γίνεται στο matlab να κάνουμε συνάρτηση που θα παίρνει σαν όρισμα άλλη συνάρτηση?
2. Η συνάρτηση τι θα παίρνει και τι θα δίνει? Στο πρώτο ερώτημα λέει να δίνει α και β και κ για διάφορα ε, ενώ στο δεύτερο λέει να δίνει τα α, β για διάφορα κ.
3. Λέει να επιλέξουμε τις βέλτιστες περιπτώσεις αλλά με βάση τι? Γιατί έχουμε σαν παράγοντες και το α,β αλλά και το κ. Μας ενδιαφέρουν τα βήματα ή το τελικό διάστημα?

1. χρησιμοποιεις το παπακι @func
2. ακριβως οπως το λες, σε μια for, και αποθηκευεις οτι σου δινει
3. το κ παιρνεις οτι σου βγει, βαζεις ενα αρχικο μεγαλο και καταγραφεις αυτο που θα παρεις σαν κ.

[pap-x]

Ευχαριστώ για την βοήθεια!
Στο 3 δεν κατάλαβα τι εννοείς. Εμείς αλλάζουμε τα ε και l. Παίρνουμε σαν αποτέλεσμα το τελικό διάστημα [α, β] και τα συνολικά βήματα κ. Άλλες φορές έχουμε καλή τελική προσέγγιση διαστήματος αλλά πολλά βήματα και άλλες λίγα βήματα και κακή προσέγγιση. Η επιλογή του ε και l θα γίνει βάσει ποιανού παράγοντα? Τελικού διαστήματος ή συνολικών βημάτων?

[edenaxas]

κοιτα, εγω εβαλα τους αλγοριθμους να τρεξουν με κ σταθερο = 100. Εαν δωσει λιγοτερα απο 100, παει να πει οτι τελειωσε για καποιο αλλο λογο, αρα τελειωσε επειδη εφτασε την επιθυμητη ακριβεια, αρα ειναι καλη λυση
Εαν δωσει == 100 παει να πει οτι δεν τελειωσε για καποιο αλλο λογο αρα δεν εχει καλη ακριβεια αρα δεν ειναι καλη λυση.

Αρα ψαχνεις να βρεις τη λυση με το μικροτερο λ που εχει τα μικροτερα βηματα. πχ λ=0.001 και 13 βηματα ειναι καλα ενω λ=0.00001 και κ= 100 βηματα δεν ειναι καλο μιας και δεν τελειωσε λογω ακριβειας αρα λογω βηματων