[Ψηφιακά Φίλτρα] - Επικαιρότητα, απορίες, ασκήσεις 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.

[tomshare3]

ανέβηκε η 2η εργασία στο ethmmy

[***]

Παίδες έχει καταλάβει κανείς ακριβώς πώς κάνουμε την επαλήθευση της σχέσης πειραματικά; Η επαλήθευση πρέπει να γίνει μόνο για τον αλγόριθμο με μέγεθος block 1?

[BluE_NoTe]

Παίδες έχει καταλάβει κανείς ακριβώς πώς κάνουμε την επαλήθευση της σχέσης πειραματικά; Η επαλήθευση πρέπει να γίνει μόνο για τον αλγόριθμο με μέγεθος block 1?

Χωρίς να θέλω να σε πάρω στο λαιμό μου... νομίζω ότι εννοεί ότι το χ (τα τέσσερα ημίτονα) μαζί με το θόρυβο είναι η είσοδος. Ας το πούμε a.  Εσύ θέλεις με το blocked LMS (γιατί δεν πρόκειται για το κλασικό LMS που είχε στην εργασία των ΨΕΣ) να αφαιρέσεις τον θόρυβο οποτε το y είναι αύτο που προκύπτει αν εφαρμόσεις τους συντελεστές w στο a. Το επιθυμητό σήμα που θέλεις να πάρεις είναι το a καθαρό, άρα το χ ! Έτσι βάζεις όπου x[n] -> a, y[n] ->a*w (έτσι οπως ορίζεται στην εκφώνηση, μπορεί να μην το συμβολίζω καλά) και d[n] -> x.
Δεν θέλω να σε πάρω στο λαιμό μου όμως γιατί και εγώ ακόμα ψάχνομαι και ακόμα το λmax ψάχνω...

Ααααα το λmax μου βγήκε 2 και κάτι. Σωστό είναι ???

[christinette]

νομίζω το παιδί ρώτησε πως στο πρώτο ερώτημα επαληθεύουμε τον αλγόριθμο...

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

[4Dcube]

αρκετά μυρισοδαχτυλική εργασία

1) ο αλγόριθμος δεν εξηγείται καλά και

2)
το σήμα με awgn προκύπτει
f+N1*f+N2*f
?

όπου f η συνάρτηση με τα ημίτονα, Ν1 η κανονική κατανομή τυχαίων με μέσο 0 και διακύμανση 0.1 και Ν2 με μέσο 0 και διακύμανση 0.3
?
αυτά από ψάξιμο στο νετ

[4Dcube]

άλλη μια ερώτηση, (έχει γίνει και σχετική άσκηση αλλά δεν την έχω)
ποιο είναι το x(n) του πίνακα αυτοσυσχέτισης Ε(x xT)? το x είναι όλο το σήμα εισόδου ή μόνο αυτό που επιλέγει ή θα επιλέξει (σύνολο) ο αλγόριθμος?

[Petros Nyfantis]

To x(n) είναι όλο το σήμα (ημίτονα  μαζί με το θόρυβο).
Χρειάζονται δύο περιπτώσεις για το θόρυβο:
μία με τυπική απόκλιση 0.1 και μια με 0.3.
Νομίζω οτι οι εξισώσεις είναι επαρκείς για την
υλοποίηση.
Τελικά αυτό που πρέπει να γίνει είναι:
να υλοποιηθούν οι εξισώσεις
και να δημιουργηθεί ένα σημα αθροίζοντας τα τέσσερα ημίτονα.
Μετά για τις δύο περιπτώσεις θορύβου,
να ελέγξετε τη συνθήκη ευστάθειας.
Να πάρετε τρείς τιμές: Μια κοντά στη μέγιστη ιδιοτιμή,
μια αρκετά μικρότερη και μια μεγαλύτερη και για τις
δυο περιπτώσεις και να σχολιάσετε το αποτέλεσμα.

Η επαλήθευση του αλγορίθμου θα είναι η ελαχιστοποίηση
του τετραγωνικού σφάλματος καθώς ο χρόνος προχωρά.

[***]

έχω τα εξής ερωτήματα:

1. το 1/λmax είναι ένα αρκετά ελαστικό φράγμα. αυτό σημαίνει ότι για αρκετά μικρότερες τιμές του μ o αλγόριθμος θα συγκλίνει. μπορεί όμως να μη συγκλίνει για τιμή του μ μικρότερη από το λmax ή να συγκλίνει για τιμές λίγο μεγαλύτερες(αν δεν κάνω λάθος). άρα μου φαίνεται και αρκετά ελαστική η επαλήθευση της συνθήκης.

2. Η συνθήκη αυτή που πρέπει να επαληθεύσουμε(αν δεν κάνω λάθος) ισχύει για τον κλασικό lms δηλαδή με μέγεθος  block ίσο με 1. άρα η επαλήθευση γίνεται για αυτήν την περίπτωση;

3. αν κατάλαβα καλά από το μάθημα, ο πίνακας αυτοσυσχέτισης της εισόδου  x (Νx1) είναι ο toeplitz πίνακας(NxN) που προκύπτει από το διάνυσμα αυτοσυσχέτισης r(lag) {Nx1}, θεωρώντας ότι η είσοδος είναι στάσιμη υπό την ευρεία έννοια στοχαστική διαδικασία. με αυτό το σκεπτικό. είναι σωστό;

4. πώς υπολογίζετε τον πίνακα αυτοσυσχέτισης της εισόδου με συναρτήσεις του matlab? z=autocorr(x,N-1) και toeplitz(z)? πάντως αυτός ο τρόπος βγάζει διαφορετικά αποτελέσματα από αλγόριθμο που υπολογίζει τον πίνακα αυτοσυσχέτισης όπως περιγράφεται στο 3.

[BluE_NoTe]

έχω τα εξής ερωτήματα:

1. το 1/λmax είναι ένα αρκετά ελαστικό φράγμα. αυτό σημαίνει ότι για αρκετά μικρότερες τιμές του μ o αλγόριθμος θα συγκλίνει. μπορεί όμως να μη συγκλίνει για τιμή του μ μικρότερη από το λmax ή να συγκλίνει για τιμές λίγο μεγαλύτερες(αν δεν κάνω λάθος). άρα μου φαίνεται και αρκετά ελαστική η επαλήθευση της συνθήκης.

2. Η συνθήκη αυτή που πρέπει να επαληθεύσουμε(αν δεν κάνω λάθος) ισχύει για τον κλασικό lms δηλαδή με μέγεθος  block ίσο με 1. άρα η επαλήθευση γίνεται για αυτήν την περίπτωση;

3. αν κατάλαβα καλά από το μάθημα, ο πίνακας αυτοσυσχέτισης της εισόδου  x (Νx1) είναι ο toeplitz πίνακας(NxN) που προκύπτει από το διάνυσμα αυτοσυσχέτισης r(lag) {Nx1}, θεωρώντας ότι η είσοδος είναι στάσιμη υπό την ευρεία έννοια στοχαστική διαδικασία. με αυτό το σκεπτικό. είναι σωστό;

4. πώς υπολογίζετε τον πίνακα αυτοσυσχέτισης της εισόδου με συναρτήσεις του matlab? z=autocorr(x,N-1) και toeplitz(z)? πάντως αυτός ο τρόπος βγάζει διαφορετικά αποτελέσματα από αλγόριθμο που υπολογίζει τον πίνακα αυτοσυσχέτισης όπως περιγράφεται στο 3.

Για την αυτοσυσχέτιση πήρα

c = xcorr(x);
cor = toeplitz(c);
lmax = eigs(cor, 1);

και μου βγήκε λmax = 2.7076... Ας μας πει ομως κάποιος για σίγουρα γιατί εγώ δεν έχω διαβάσει την θεωρία για να ξέρω αν έχει νόημα αυτό που έκανα 

[Petros Nyfantis]

Το τελευταίο είναι το σωστό.
Αν θέλετε, επαληθευστε τη συνθήκη για για κάποιες κοντινές
τιμές του στο μ. πχ αν υπολογίσατε 2.5, παρτε
με διάστημα 0.05 τις τιμές 2.3 με 2.7 και σαν έξοδο
το τετραγωνικό σφάλμα μετα απο πχ 250 δείγματα.

[***]

Για την αυτοσυσχέτιση πήρα

c = xcorr(x);
cor = toeplitz(c);
lmax = eigs(cor, 1);

και μου βγήκε λmax = 2.7076... Ας μας πει ομως κάποιος για σίγουρα γιατί εγώ δεν έχω διαβάσει την θεωρία για να ξέρω αν έχει νόημα αυτό που έκανα 

πώς γίνεται να βγαίνει τόσο μικρή η μέγιστη ιδιοτιμή; πόσα δείγματα παίρνεις για την είσοδο. με τη μέθοδο σου
η λmax μου βγαίνει 444

[Petros Nyfantis]

Δοκιμάστε
c = xcorr(x(1:tr),'biased');

[***]

κανονικά ποια ποια πρέπει να βγαίνει η μέγιστη ιδιοτιμή του πίνακα αυτοσυσχέτισης της εισόδου αν πάρουμε 500 δείγματα?

[***]

επίσης παρατηρώ ότι για διαφορετική τάξη του φίλτρου ο αλγόριθμος συγκλίνει για διαφορετικές τιμές του μ (δηλαδή το φράγμα αλλάζει)

[BluE_NoTe]

Όταν λέει πως η τάξη του φίλτρου είναι 32 υποθέτω πως εννοεί kmax = 32 και ότι αυτό είναι το μέγεθος του διανύσματος στήλης w. Ο αριθμός k όμως είναι ο αριθμός των blocks που χωρίσαμε τον χώρο δειγμάτων μας. Δηλαδή αν χ τα δείγματα k = χ/L. Εάν θέλουμε το k σταθερό και η εκφώνηση θέλει L =1 ή L=32 τότε φαντάζομαι πως τα δείγματά μας θα είναι 32 ή 1024 αντίστοιχα. Το λ=2.7076 το βρήκα για 32 δείγματα. Ελπίζω να είναι έτσι τα πράγματα...
Για όσους έχουν λίγο χρόνο ας ρίξουν μια ματιά σε αυτο το pdf. Νομίζω πως είναι πολύ βοηθητικό :

cwww.ee.nctu.edu.tw/course/asp/ASP06.pdf