THMMY.gr

Ανέβηκε στο Υλικό Μαθήματος η 4^η εργασία.

Η προθεσμία λήγει στο τέλος της εξεταστικής.

Συμβαίνουν 2 πράγματα που μου φαίνονται λίγο περίεργα.
1. Στο ερώτημα (β) το τετραγωνικό σφάλμα μεταξύ των παραμέτρων wo και αυτών που βρίσκω με τις διαφορες μεθόδους για φίλτρο 4 συντελεστών είναι γυρω στο 10^-2 με 10^-3
2. Στο ερώτημα (γ) οι LMS και NLMS δε δινουν ικανοποιητικό αποτέλεσμα, δηλαδή εντάξει πλέον ακούγεται πιο καθαρά η φωνή του ομιλητή Β, αλλά ακούγεται αρκετά ενοχλητικά και η φωνή του ομιλητή Α. Ωστόσο, με τον RLS δεν έχω καταφέρει να πάρω αποτέλεσμα γιατί τρέχει τρομερά αργα (χρησιμοποιω αυτόν τον κωδικα που έχει ανεβασει ο Σισμάνης και φίλτρο 6700 συντελεστών)

Κανείς που να του έχουν δουλέψει? Κάνω κάτι λάθος?  :-[

Στο τετραγωνικό σφάλμα παίρνεις για ένα trial ή παραπάνω; Πάντως λογικά νούμερα βγάζεις καθώς και εγώ τα ίδια βγάζω για ένα trial.

Όσον αφορά το ερώτημα (γ) αντιμετωπίζω και εγώ το ίδιο πρόβλημα. Με τον LMS μειώνεται η φωνή του Σισμάνη αισθητά αλλά και πάλι ακούγεται περισσότερο απ' ότι την αντίστοιχη έξοδο με Wiener φίλτρο. Τι βήμα σύγκλισης πήρες για το LMS και NLMS?

Για τον RLS με 6600 συντελεστές και εγώ δεν κατάφερα να το τρέξω μέχρι το τέλος για να δω αποτέλεσμα. Έβγαλες άκρη τι πρέπει να κάνουμε; Καθώς αν βγάζει σωστά ο RLS σημαίνει πως δεν έχουμε κάνει πατάτα και πως όντως οι άλλοι δύο αλγόριθμοι απλά δεν δίνουν καλά αποτελέσματα σε αυτήν την εφαρμογή.

( BTW έτρεξα τους αλγορίθμους με είσοδο την μουσική και τον θόρυβο της 1ης εργασίας και καθάρισαν πλήρως το κομμάτι από τον θόρυβο. Άρα, παίζει να τους υλοποιήσαμε σωστά και απλά να μην μπορούν να δώσουν καλά αποτελέσματα σε αυτήν την εφαρμογή)

Για την τάξη του φίλτρου έχουμε κάποια πληροφορία;

edit: Μάλλον θα πάρω απλά Μ=11 όπως έχει στο πρότυπο του adaptRLS

Έχει καταλάβει κανείς πως θα χρησιμοποιήσουμε την πληροφορία για την καθυστέρηση των 6600 χρονικών στιγμών που λέει στο ερώτημα γ?

Παίρνεις φίλτρα 6600 συντελεστών.

Τι εννοείς για ένα trial? Παίρνω για 25000 χρονικές στιγμές π.χ.
Για τον RLS δεν έχω βγάλει ακόμα ακρη. Πέρασα και από το Σισμάνη ο οποίος έδειξε να εκπλήσσεται με το πόσο αργά τρέχει, αλλά δεν είχε να προτείνει κάτι.

O RLS για φίλτρο 6600 συντελεστών και για 1 trial δεν τελειώνει ποτέ (τον άφησα πάνω από 6 ώρες και αποτέλεσμα δεν πήρα!!). Γι αυτό ρωτάω μήπως οι 6600 χρονικές στιγμές πρέπει να χρησιμοποιηθούν κάπως αλλιώς..

Εννοώ πως για να βγάλεις το μέσο τετραγωνικό σφάλμα στον LMS πρέπει να το τρέξεις καμία 50 φορές και να πάρεις την μέση τιμή αλλιώς δεν έχεις το μέσο τετραγωνικό σφάλμα αλλά το τετραγωνικό. Για το αποτέλεσμα του LMS τι είπε; Έπρεπε να βγαίνει καλό;

Τι να σου πω δεν ξέρω. Στο ίντερνετ το βρήκα αυτό που σου είπα.Δύο τα ενδεχόμενα:
 ή κάνουμε πατάτα με τους συντελεστές
 ή ο Σισμάνης είπε δοκιμάστε τον RLS για να παραλλάξει την εργασία από την περσινή, χωρίς να το τρέξει καν για να δει πως πάει.  :P

ααα... ναι για 1 trial το έτρεξα, αλλά το έχω τρέξει πολλές φορές και δεν υπάρχει διαφορά.
Είπε ότι πρέπει να βγαίνει χειρότερο από του RLS και ότι με τον RLS πρέπει να καθαρίζει πλήρως.

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

Δεν σου είπε δηλαδή "εγώ το έτρεξα και ήταν ΟΚ"...

Μου είπε επί λέξη "Με τον RLS πρέπει να καθαρίζει τελείως. Πάντως μου κάνει εντύπωση που τρέχει τόσο αργά. Τι υπολογιστή έχεις? (συζητηση με την οποία καταλήξαμε ότι δε φταίει ο υπολογιστής) Εμένα μου έτρεχε σχετικά γρήγορα, αλλα μπορεί να είχα βάλει... Ή όχι... Δε θυμάμαι... Δεν ξέρω... Τι να σου πω?"

Απλά θεός ...  8))

Η προθεσμία ειναι μεχρι το τέλος ολης της εξεταστικής ή μεχρι σήμερα που δίναμε;

από eTHMMY

4η εργασία (Hμ/νία λήξης δηλώσεων: 08/07/2015)

Εξαιρετικά. Ευχαριστώ πολύ :)

Μπόρεσε κανείς να τρέξει τελικά τον RLS για το φιλτρο των 6600 συντελεστών; Αν ναι πως;

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

r = autocorr(u, M - 1);
r = r*var(u);
R = toeplitz(r);
p = xcorr(d, u, M/2, 'unbiased');
p = p(2:end);
Wo = R\p;
y = zeros(N, 1);
for n = M:N
    y(n) = u(n:-1:(n - M + 1))'*Wo;
end
e = d - y;

Μήπως έχει κάποιος καμιά ιδέα τι κάνω λάθος?

Στο ερώτημα γ, για τον υπολογισμό του σήματος y στον LMS χρησιμοποιώ τη σχέση
y(i) = wl'*u(i:-1:(i-M+1));
Από ένα σημείο και μετά προκύπτουν τιμές NaN. Τι μπορεί να πηγαίνει στραβά; (στο ερώτημα β δουλεύει κανονικά)

Ουτε εγω εχω καταφερει να καθαρισω απο ηχο.

Εγω κανω αυτο   

yl(i) = wl(:,i-1)'*u(i:-1:(i-M+1));
el(i) = d(i) - yl(i);
wl(:,i) = wl(:,i-1) + mu*el(i)*u(i:-1:(i-M+1));

καμια ιδεα?

Βγάζεις ήχο αλλά δεν έχει καθαρίσει τελείως από την ηχώ;

οχι βγαζω μια βλακεια 5 δευτερολεπτα περιπου, μονο θορυβο.

Ο normalized σου τρέχει; το wl πώς το ορίζεις; γιατί είναι δύο διαστάσεων;

Βασικα μαλακια ποσταρα, το διορθωσα αλλα εχω το ιδιο προβλημα με εσενα, βγαζω NaN τιμες

Δεν κατέληξα κάπου πάντως.. Ο normalized όμως λειτουργεί και δεδομένου ότι έχουν παρόμοια λειτουργία είναι περίεργο. Κάπου στην ανανέωση των βαρών γίνεται η μαλακία.

Ναι το διανυσμα του wl einai olo γεματο ΝaN, και εγω δεν κατεληξα καπου

Γμτ εμενα ουτε ο nLMS δουλευει...

Με τους αλγορίθμους LMS, normalized LMS και RLS προσεγγίζετε τους συντελεστές Wiener ή αυτοί αποκλίνουν;