[Ψηφιακά Φίλτρα] 3η-4η εργασία 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.

[greekoo]

1)E[x(n)x(n)] = σ_x² = r_x(0) , γενικά μπορείς να το βρείς
2)E[x(n)v(n)] = 0 υποθέτω , επειδή είναι σήμα πληροφορίας με θόρυβο (στο matlab βγαινει 0.08 * 10^-3 οπότε μάλλον όντως 0 θα βγαίνει )
3)E[(v(n-10)x(n))] = 0 για τον ίδιο λόγο με πάνω
4)E[v(n-10)v(n)] = 0 επειδή ο θόρυβος δεν έχει συσχέτιση με προηγούμενες τιμές του

Άρα νομίζω μπορείς να πείς ότι τελικώς  p=E[u(n)d(n)]= E[x(n)x(n)] = σ_x² = r_x(0)

συμφωνώ φίλε μου! Αλλά για το κομμάτι ήχου, εμείς δεν ξέρουμε το x(n). Ξέρουμε μόνον το s(n) και άρα και το u(n). Οπότε πώς θα βρούμε το p?

Ή υποννοείται πώς σε όλα τα ερωτήματα του πρώτου μέρους, το x(n) είναι αυτό που δίνεται στον τύπο;

[c0ndemn3d]

Για το κομμάτι του ήχου προφανώς δεν ξέρεις το χ. Βρίσκεις τις αυτοσυσχετίσεις με matlab. Το d ειναι το s και με το u ειναι παλ το s με καποια καθυστερηση. Αρα τις αυτοσυσχετισεις μποιρεις να τις βρεις

[Sf(x)dx]

συμφωνώ φίλε μου! Αλλά για το κομμάτι ήχου, εμείς δεν ξέρουμε το x(n). Ξέρουμε μόνον το s(n) και άρα και το u(n). Οπότε πώς θα βρούμε το p?

Ή υποννοείται πώς σε όλα τα ερωτήματα του πρώτου μέρους, το x(n) είναι αυτό που δίνεται στον τύπο;

To x(n) σου το δίνει από την εκφώνηση, τουλάχιστον μέχρι το ερώτημα 3). Το θέμα είναι οτι ενώ σου το δίνει, κανονικά οχι, δέν το ξέρεις. Τα μόνα γνωστά είναι τα σήματα d(n) και u(n) που στην συγκεκριμένη περίπτωση τυχαίνει να είναι και τα ίδια, δλδ : d(n) = u(n) = x(n)+v(n)

[greekoo]

για την ετεροσυσχέτιση για 100 συντελεστές αρκεί να γράψω αυτό;

P = xcorr(u,d)
P = P(1:100);

[Sf(x)dx]

Όχι, αυτό δεν είναι σωστό.

1) Αρχικά, έχουμε κανονικοποιημένη συσχέτιση, δηλαδή για παράδειγμα:

Code:

rx(0) = E [x(n) * x(n)] = [Σ x(i) * x(i)]/n  , i = 1, ..., n     και όχι
rx(0) = Σ x(i) * x(i)  , i = 1, ..., n  

(σόρρυ για το παραπάνω αλλά δεν ξέρω Latex για να το γράψω σαν άνθρωπος)

Η xcorr by deafault κάνει το 2ο, αλλά μπορεί να κάνει και το 1ο αν προσθέσεις την επιλογή 'unbiased'

2) Η xcorr() σου επιστρέφει κανονικά και ένα άλλο διάνυσμα που λέγεται lags, δηλαδή συνήθως την καλείς έτσι :  [p,lags] = xcorr(d,u, ... )
Αν χρησιμοποιήσεις σωστά αυτό το διάνυσμα, μπορείς να κόψεις από το συνολικό διάνυσμα που σου επιστρέφει η xcorr τις τιμές που σε ενδιαφέρουν, και να σχηματήσεις έτσι το διάνυσμα ετεροσυσχέτισης p από την θεωρία.

Ρίξε μια ματιά εδώ http://www.mathworks.com/help/signal/ref/xcorr.html , τα εξηγεί καλά.

[greekoo]

Σε ευχαριστώ πολύ φίλε, αυτό με τα lags ήταν πολύ βοηθητικό!

Kάτι τελευταίο. Ισχύει ότι u(n) = s(n-Δ)
άρα στο matlab μπορώ να γράψω αυτό u = [zeros(Delta,1);s(1:1:end-Delta)] ; ??

[Sf(x)dx]

yeap, αυτό θα δουλέψει για να βρείς το delayed διάνυσμα εισόδου.

[greekoo]

Ρε παιδιά στο Echo Cancellation, καταφέρατε να το τρέξετε με RLS για 6600 συντελεστές? Μου τρέχει πραγματικά για ώρες... Και βλέπω πώς και πέρυσι υπήρχε το ίδιο πρόβλημα...

[c0ndemn3d]

Λοιπόν μου 'χει μείνει για το τέλος αλλά θα πρέπει να βρεις λογικά μία συνθήκη τερματισμού ικανή να σου πει ότι έχει εκπαιδεύσει το φίλτρο σου. Είναι RLS ο αλγόριθμος άρα θα έχεις σίγουρα γρήγορα σύγκλιση. Το θέμα είναι να δεις και λίγο το u πώς μεταβάλλεται, δηλαδή στην αρχή έχει αρκετά πολλές μηδενικές τιμές άρα εκεί δεν υπάρχει μεταβολή του w. Μετά από εκείνο το στάδιο μάλλον θα πρέπει να περιμένεις μέχρι το u να παίρνει μεγάλες τιμές, ώστε αν παρατηρήσεις μετά ότι το w δεν μεταβάλλεται πολύ, να σταματάς τον αλγόριθμο και τα υπόλοιπα y να τα βρίσκεις θεωρώντας ότι έχεις αυτά τα συγκεκριμένα w που βρήκες στο τέλος.

[greekoo]

Λοιπόν μου 'χει μείνει για το τέλος αλλά θα πρέπει να βρεις λογικά μία συνθήκη τερματισμού ικανή να σου πει ότι έχει εκπαιδεύσει το φίλτρο σου. Είναι RLS ο αλγόριθμος άρα θα έχεις σίγουρα γρήγορα σύγκλιση. Το θέμα είναι να δεις και λίγο το u πώς μεταβάλλεται, δηλαδή στην αρχή έχει αρκετά πολλές μηδενικές τιμές άρα εκεί δεν υπάρχει μεταβολή του w. Μετά από εκείνο το στάδιο μάλλον θα πρέπει να περιμένεις μέχρι το u να παίρνει μεγάλες τιμές, ώστε αν παρατηρήσεις μετά ότι το w δεν μεταβάλλεται πολύ, να σταματάς τον αλγόριθμο και τα υπόλοιπα y να τα βρίσκεις θεωρώντας ότι έχεις αυτά τα συγκεκριμένα w που βρήκες στο τέλος.

Έχεις δίκιο φίλε, όντως πρέπει να έχουμε μια συνθήκη τερματισμού.
Όμως μόνο και μόνο αυτή η εντολή: P = (lambda^-1)*P - (lambda^-1)*k*u(i:-1:(i-M+1))'*P;
μου παίρνει 20-25 δευτερόλεπτα!!! 

[chrimpou]

Προθεσμία: 15 Ιουνίου
Έχει ακούσει κανείς για παράταση προθεσμίας?

[c0ndemn3d]

Όχι γιατί κανείς δεν ζήτησε, λογικά όταν θα πλησιάζει ο καιρός και ζητήσουν πολλοί, θα δοθεί.

[jiannist]

Έστειλα mail, και μου είπε ότι δε θα δώσει γιατί δεν έχει νόημα να την κάνουμε αφού δώσουμε το μάθημα.

[WatchDog]

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

[jiannist]

Ρώτησε κάποιος άλλος για παράταση;