THMMY.gr

Δύο σε ένα νοικοκυρεμένα στο ethmmy.

Προθεσμία: 15 Ιουνίου

Τζιζους φακινγκ κραιστ

Guys στην εκφώνηση να υποθέσω εκεί που λέει :
Μήπως έχει ξεχάσει στο τελευταίο συνημίτονο ένα f ? Έχει ρωτήσει κανείς τον πιτσιάνη ?

Βρήκε κανείς κάτι αξιόλογο στο πρώτο τμήμα;

define αξιόλογο

Βγάζω ό,τι να 'ναι φίλτρο. Λέω αν σε κάποιον έκανε καλή δουλειά.

Εμένα για το πρώτο μέρος μου βγήκε πάντως πως το σφάλμα e(n) θα πρέπει να συγκλίνει στο τυχαίο σήμα v(n) και συνεπώς η έξοδος του φίλτρου y(n) είναι η περιοδική παρεμβολή x(n). Δεν το έχω αποδείξει ακόμα βέβαια, αλλά με το script αμα το προσομοιώσεις αυτό βγαίνει πάντως.

Ναι, τυπογραφικό λάθος είναι.


Έτσι είναι, το σήμα y(n) προσεγγίζει το x(n) και το σήμα e(n) προσεγγίζει το σήμα v(n).

Οι συντελεστές γ στο πρώτο μέρος στο γ ερώτημα ποιοι είναι;

EDIT: Πρέπει να είναι οι g0 απ'τις διαφάνειες

δες το pdf lattice-filters σελ.8

χμμ, δεν το είχα δεί αυτό. Άρα τελικά το y θα είναι οι συντελεστές γ πολλαπλασιασμένοι με το u (με τα σωστά lags) ?

Την 4η εργασία που την είδατε; Είναι μεγάλη η 3η και δε θα μπει 4η φέτος;

Είναι 2 τμήματα. Δες από κάτω

Αρχικά, ισχύει αυτό;
p=E[u(n)d(n)]=E[(x(n)+v(n-10))(x(n)+v(n))]=E[x(n)x(n)]+E[x(n)v(n)]+E[(v(n-10)x(n))]+E[v(n-10)v(n)]

Αν ναι, πως λύνεται; Λογικά, το E[v(n-10)v(n)] είναι 0. Αλλά πως υπολογίζονται τα υπόλοιπα;

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) είναι αυτό που δίνεται στον τύπο;

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

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

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

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

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

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

(σόρρυ για το παραπάνω αλλά δεν ξέρω 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) , τα εξηγεί καλά.

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

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

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

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

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

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

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

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

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

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

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

του έστειλα και εγώ και δεν δίνει. παιδιά στείλτε και άλλοι μπας και δώσει!

Τελειωνοντας την 3 εργασια εχω να πω μονο
https://www.youtube.com/watch?v=pCXVxE_YeP4

Γερό κομμάτι, ο Σισμάνης ξέρει

Ξέρετε πόσο μετράνε οι εργασίες;

0.5 η πρώτη
0.5 η δεύτερη

για την τρίτη δεν είπανε ακριβώς, αλλά είπαν ότι θα είναι παραπάνω από 0.5 (λογικα 0.8 με 1)

Δεν πιστεύετε ότι είναι λίγο, σε σχέση με άλλα μαθήματα, για τον κόπο που απαιτούν;

Στο Α μέρος, στο δ), πως βρίσκετε τον πίνακα R, και το p;
Εγώ κάνω αυτά:

u = [zeros(Delta,1);s(1:1:end-Delta)] ;

[r, lags] = xcorr(s, M-1, 'unbiased');
r=r(lags >= 0);
R = toeplitz(r);

[p, lags] = xcorr(s, u, M-1, 'unbiased');
p = p(lags >= 0);

wo = R \ p;


Όμως, το τελικό σήμα (φιλτραρισμένο με συντελεστές Wiener) μοιάζει πιο πολύ στο αρχικό, παρά σε φιλτραρισμένο.

Kαλησπέρα , έχουμε κανά νέο για παράταση μήπως;

Γράφει πιο πάνω ότι δεν θα δοθεί. Μπορείς να ρισκάρεις και να τις δώσεις πιο μετά σε στυλ παράλληλα που η πραγματική προθεσμία είναι λίγο πριν βγουν οι βαθμοί.

Είναι δηλαδή όποτε μαντέψεις ότι δεν έχουν ξεκινήσει να διορθώνουν ! χαχαχ  ;D ;D ;D

Πώς στο @@ υπολογίζετε τους συντελεστές γ στο Α μέρος? ? >:( >:(

Που κολλάς; Και στις διαφάνειες και στο Haykin λέει ξεκάθαρα πως βγαίνουν

u(x) = x(n-10) + v(n-10)
άρα μήπως είναι p=E[u(n)d(n)]= E[x(n)x(n-10)] = σ_x² = r_x(10) ;

Ούτε εγώ κατάλαβα πως βγαίνουν τα γ για να είμαι ειλικρινής . Στη διάλεξη στα lattice εννοεις οτι το λεει ξεκάθαρα?

ναι, όπως λεει και στις διαφάνειες

Οταν λεει στο Β μερος στο γ) μπορει να ειναι μεχρι 6600 χρονικες στιγμες ,εννοει να κανουμε φιλτρο με 6600 συντελεστες?

Στο Α.γ που λέει να επαληθεύσουμε τα αποτελέσματα με τις έτοιμες συναρτήσεις του Matlab, ποιες συναρτήσεις εννοεί;

levinson()

Σύμφωνοι. Αλλά αυτό θα μας δώσει το gamma σε διαστάσεις 101x1.
Όμως το βέλτιστο w0 είναι 100x1.
Και η σχέση που συνδέει αυτά τα δύο είναι L' * γ = w0 , κάτι που δεν ταιριάζει σε διαστάσεις!!!

Εκεί κολλάω!

εχω την ιδια απορια! μπορει να απαντησει καποιος αγαπητος συναδελφος??  ::)

Στα πρώτα ερωτήματα, η περίοδος του x που δίνει είναι 8, σωστά; Το φίλτρο με βγάζει τελεια αποτελεσματα για Δ=8 ενώ για Δ=10 που λεει στην εκφωνηση τα κανει χαλια. Εχει κανεις παρομοιο προβλημα;

Ρε παιδια μονο σε εμενα βγαινει το τραγουδι στην εξοδο του φίλτρου;;

Θα πάρει καμιά παράταση η προθεσμία μήπως;

Πέτυχα τον Φλώρο προχθές και τον ρώτησα και μου είπε όχι... Λέω να την κάνω πάραυτα και να την στείλω εκπρόθεσμα....