THMMY.gr

Δεύτερη εργασία στο ethmmy
http://alexander.ee.auth.gr:8083/eTHMMY/archive/121/downloadFile/6378/AF-assignment2-2016.pdf

Προθεσμία :24 Απριλίου 8 Μαϊου

Στο γ ερώτημα θέλει M = 3000(άρα και L = 3000) σωστά?
Έχει βγάλει κανείς καμιά άκρη με τις τιμές του μ?

Μπορεί κάποιος να εξηγήσει τι ακριβώς θέλει στο 1ο ερώτημα?

για μ >=0.01 ( και μ<(2/λmax) ) οι συντελεστές μου παίρνουν πολύ μεγάλες τιμές... ίσως εγώ κάνω λάθος κάτι

τελικά κάπου είδα ότι για μεγάλες τιμές του L πρέπει μ < 2/(L * λmax) και πράγματι για τόσο μικρές τιμές, το error πέφτει στο 10^-20 τουλάχιστον

α τώρα μάλιστα φίλε..

Για να κάνουμε τα learning curves, το J θα το υπολογίζουμε μια φορά για κάθε μπλοκ, ή σε κάθε επανάληψη μέσα στο μπλοκ όπου υπολογίζεται και το e(n) ?

το δεύτερο

θέλει να μετατρέψεις σε κώδικα matlab την σχέση που έχει στην εκφωνηση(αυτή που είναι σε μορφή πινάκων, όπου τα στοιχεία του πίνακα υπάρχουν στην σχέση ακριβώς από πάνω)

προσοχή! μην κάνετε την αναστροφή των διανυσμάτων με τη χρήση του " ' " ,καθώς αυτό δημιουργεί το ανάστροφο συζυγές διάνυσμα. δεν το ήξερα και μου πήρε 1 μέρα να βρω τι είναι.

στο β που θέλει να υπολογίσουμε το υπολογιστικό κόστος δεν καταλαβαίνω ποιοι έιναι οι δύο αλγόριθμοι, το κόστος των οποίων θα υπολογίσουμε. αν έχει καταλάβει κάποιος ας μας πει.

ο block lms που ζητάει υπάρχει κάπου ή τον φτιάχνουμε εμείς?

για να υπολογίσουμε το μ πρέπει να βρούμε τις ιδιοτιμές του πίνακα αυτοσυσχετισης, σωστά? όμως επειδή τα δεδομένα είναι πολλά ο πίνακας αυτοσυσχέτισης βγάινει τεράστιος και το matlab δεν μπορεί να τον υπολογίσει. τι άλλο μπορώ να κάνω?

Η εργασία θα πάρει παράταση δυο βδομάδες από οτι ειπώθηκε,αναμένετε ανακοίνωση

Ρε παιδιά στο δεύτερο ερώτημα πώς θα τα βρούμε τα flops?  Ξέρουμε τις πολυπλοκότητες των FFT και DFT (και απ΄την σχέση που μας δίνεται βρίσκουμε την nlogn του FFT), αλλά τα flops πώς θα τα βρούμε?

Σύμφωνα με αυτό το site:  http://www.fftw.org/speed/method.html

thoughts?

Αυτό είναι σίγουρο?

Όχι λέω ψέματα για να την κάνω μονο εγώ

πως με κατάλαβες οτι το ρώτησα γι' αυτό!!  :-X
απλά δεν έχουν βγάλει ανακοίνωση για αυτό λέω τι γίνεται...

Ο Πιτσιανης είπε ανεπίσημα έχει πάρει σε κάποια φάση θα το ανακοινώσουν και επίσημα.Μην τον πιεζουμε κιόλας τόσο πράγματα έχει να κάνει

Τί μένει να κάνουμε στο 1ο ερώτημα στην fftproof; Τά'χει κάνει όλα ήδη αυτός.

την δεύτερη γραμμή της εξίσωσης.Ουσιαστικά αυτό που κάνει στο τελευταίο απλά εσύ θα βρεις τα fe και fo με τον τρόπο στη δεύτερη γραμμή της εξίσωσης.

Μπορεί κάποιος να δώσει μία κατατοπιστική απάντηση περί του ζητήματος?

θα μείνουμε με την απορία φίλε :Ρ.

Επίσης, στο δ) ερώτημα του BLOCKLMS στην περίπτωση του  unconstrained method, κατάφερε κανείς να το κάνει; Το βήμα μας θα είναι το alpha ή το mu ξεχωριστό για κάθε βάρος; Και την αρχική ισχύς πόσο να την πάρουμε; (γαμω τα φιλτρα  ;D ;D)

Μπορείς να το κάνεις είτε με τον έναν τρόπο είτε με τον άλλο και το ίδιο μπορείς να κάνεις και για το γ ερώτημα. Αν δεν υλοποιήσεις την παραλλαγή του αλγορίθμου με το alpha ως παράμετρο, τότε το mu παραμένει το ίδιο για όλα τα βάρη. Η διαφορά ανάμεσα στον constrained block LMS και στον αντίστοιχο unconstrained είναι ότι ο πρώτος περιλαμβάνει δύο επιπλέον FFTs.

Άσχετο, οι υποβολές για την εργασία θα ανοίξουν ποτέ;

Σας βγαίνει το d ίδιο με το u?  :o

Αν εννοείς πως το plant βγαζει την ίδια είσοδο με την έξοδο, τότε ΝΑΙ.

Δεν είναι περίεργο;;

Τελικά τα πρώτα 3000 είναι ίδια.

Μπορεί να μου εξηγήσει κάποιος πώς βγαίνει ο αναδρομικός τύπος για την πολυπλοκότητα; Εγώ όπως το βλέπω, έχουμε n/2 πολλαπλασιασμούς, οπότε το n/2 * 6 flops που έχει στον τύπο και μετά n προσθέσεις. Εδώ τώρα αυτός έχει n/2 * T(2), αυτό από πού βγαίνει; Και πόσο δηλαδή είναι το Τ(2);

Η αναδρομή είναι ότι να 'ναι παρ' όλα αυτα αν την λύσεις βγαζει σωστό αποτέλεσμα και ισοδυναμο με αυτό που λες. Επίσης θεωρώ ότι Τ(2)=4

Εντωμεταξύ το σύστημα που μάς δίνεται συμπεριφέρεται πολύ περίεργα στις πρώτες 3000 εισόδους, καθώς βγάζει έξοδο ίδια με την είσοδο, ενώ πιο μετά φαίνεται να παίρνονται υπ' όψιν μόνο πολύ παρελθοντικές τιμές, άρα δεν είναι γραμμικό. Λέω βλακεία; Το φίλτρο που βγάζω έχει άριστη συμπεριφορά για n > 3000

Tα λες εξαιρετικά!

Η είσοδος είναι ίδια με την έξοδο για τις πρώτες 3000 εισόδους γιατί το σύστημα είναι ένα φίλτρο με 3000 taps, οπότε για αρχικές τιμές βάζει την είσοδο. Αυτή είναι και η default συμπεριφορά της filter() πχ. Γιαυτό ζητάει και φίλτρο τόσων συντελεστών άλλωστε.

μπορει καποιος να δειξει μια καμπυλη εκμαθησης του συστηματος του; γιατι ειμαι τερμα μπερδεμενος στο αν βγαζω σωστα αποτελεσματα. Επισης ποια ειναι μια λογικη τιμη του μ για να βάλουμε; Ξερω για το 2/λmax αλλα νομιζω μου βγαινει πολυ μικρο.

Καμπύλες εκμάθησης δεν έχω φτιάξει αλλά (αν δεν έχω κάνει μαλακία) θέλει ελεεινά μικρό μ για να συγκλίνει. Χρησιμοποίησα mu=0.0001/max(eig(R));

Α οκ και εγω ελεεινα μικρο. Ομως ετσι οπως το υπολογισα μου εβγαινε το 2/max(eig(R)) = 2.κατιe-04 και εβαλα mu=2.0e-4 και ηταν ψιλοκομπλε

Στο πρώτο ερώτημα ,μπορεί κάποιος εξηγήσει τί ζητάει? Κάποιος παραπάνω είπε πως θέλει τη 2η γραμμή της εξίσωσης που δίνει στην εργασία(με τους πίνακες) αλλά δε το πιάνω .

πρέπει να γράψεις την εξίσωση που υπάρχει στη δεύτερη γραμμή σε κώδικα ματλαμπ για να υπολογιστεί ο FFT

Στο ερώτημα για τον BLMS όταν ζητάει υλοποίηση με 2 nested loops αυτά τα loops είναι:
1) Ενα για το k (block index),
2) Ενα για το i μεσα σε καθε block?

Ή μήπως υπάρχει και ένα τρίτο loop για τα αθροίσματα στο βήμα 2.1.1 στον ψευδοκώδικα από τις διαφάνειες? Ρωτάω γιατί αν δεν χρειάζεται το 3ο loop τότε στο β ερώτημα θα πρέπει να μείνει μόνο το iteration για το k?

ετσι το καταλαβα και εγω.

ναι έτσι είναι. Και στο β ερώτημα έχεις μόνο το loop για το k.

Το αποτελεσμα μας θα πρεπει να βγαίνει ίδιο με το δικό του -> fprintf('done : %e\n', norm(y - ys)) ??

Ναι το norm(y-ys) πρεπει να βγαινει πολυ μικρο

Το plant.p περιέχει το σήμα d ?Το κάνουμε load ?γιατι οταν το κάνω load μου πετάει error οτι δε μπορεί να το διαβάσει. Επίσης σαν αρχικές τιμές για τα W μπορούμε να χρησιμοποιήσουμε τον κώδικο απτο ethmmy του LMS?

Το plant.p το χρησιμοποιεις σαν συναρτηση.
d = plant(u);

To Μέγεθος του Κ θα είναι 3000 (όσο και του φίλτρου) (?) .
Για το μέγεθος του L αυτοσχεδιάζουμε?εσεις πόσο βάλατε?

Έχω φάει σκάλωμα, αλλά μου φαίνεται εντελώς λάθος η σχέση που μας δίνει στο α) ερώτημα

Έχουμε ενα Ω NxN που πολλαπλασιάζεται με ενα διάνυσμα x⁽²⁾_k διαστάσεων 1x N/2 ( ή N/2 x 1), κάτι το οποίο δε γίνεται.. τί φάση ?

και απο την άλλη το x γιατί έχει 2 στήλες κάθεμια (υποτίθεται) διαστάσεων 1xN ...

Καμία διευκρίνηση ?

Ναί, το plant( )  είναι συνάρτηση, και βάζεις ώς είσοδο ένα διάνυσμα γραμμή τουλάχιστον 3000 στοιχείων


Το k(block index) είναι loop variable από 0 έως k_max (ο αριθμός των συνολικών blocks, πρακτικά πόσες φορές θα ανανεώσεις την εκτίμηση που έχεις για το διάνυσμα w). Το k_max το διαλέγεις αυθαιρετα. Από την άλλη το L είναι το block size και επιλέγεται νομίζω ίσο με  M, δηλαδή 3000 (εδώ δεν είμαι 100% σίγουρος).


Ο πίνακας Ω δεν είναι NxN, είναι kxk και εκεί εννοεί ότι το κ είναι από 0 έως N/2-1, δηλαδή Ν/2 στοιχεία. Συνεπώς και ο πίνακας Ω είναι (N/2 x N/2), όπως και το διάνυσμα  x₂ το οποίο είναι (Ν/2 x 1) συνεπώς ο πολλαπλασιασμός γίνεται. Μπορείς να το κάνεις όπως το κάνει ο Πιτσιάνης στο script που μας έχει δώσει. Συγκεκριμένα σε αυτές τις γραμμές:

στο 3ο ερωτημα,  τα βήματα (n) πόσο τo βάλατε?

1.2e6

παιδιά ουσιαστικά στο πρώτο ερώτημα κάνω μία loop που βάζω την δεύτερη γραμμή και μετά επαλήθευση?

max(eig(R)) πόσο έβγαλες?