[Ψηφιακά Φίλτρα] 2η εργασία 2016

[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 με προσοχή στα ονόματα των αρχείων!

[Αντιστοίχηση Μαθημάτων ΝΠΣ με ΠΠΣ]

Επίσημη ενημέρωση για Αντιστοίχηση Μαθημάτων ΝΠΣ με ΠΠΣ και η συζήτηση στο forum.

[orestisf]

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

[Κάρβουνο]

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

[orestisf]

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

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

[WatchDog]

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

[themis93]

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

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

[Conflict]

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

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

[orestisf]

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

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

[Gru]

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

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

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

[DarkPassenger]

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

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

[orestisf]

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

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

[WatchDog]

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

[orestisf]

Το plant.p περιέχει το σήμα d ?Το κάνουμε load ?γιατι οταν το κάνω load μου πετάει error οτι δε μπορεί να το διαβάσει.

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

[WatchDog]

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

[Xbaremenos]

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

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

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

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

[Sf(x)dx]

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

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

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

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

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

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

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

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

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

Code:

% to complete the proof
fe = fft(x((0:2:n-1) +1));
fo = fft(x((1:2:n-1) +1));

wfo = w(n,(0:n/2-1)') .* fo; 
y = [fe + wfo; fe - wfo];