[Ψηφιακά Φίλτρα] 3η Εργασία 2015

[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.

[Exomag]

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

Η προθεσμία λήγει την Πέμπτη 4 Ιουνίου.

[opcode]

Στο ερώτημα (β) πως υπολογίζεται ο πίνακας R και το διάνυσμα ετεροσυσχέτισης p 100 συντελεστών; Σίγουρα γίνεται μέσω Matlab καθώς με το χέρι δεν παίζει να βγαίνει .

[Silvo the Beautiful]

ποιες διαφάνειες χρειάζονται για διάβασμα για αυτήν την εργασία;

[Silvo the Beautiful]

κανείς;

[Kodi]

ποιες διαφάνειες χρειάζονται για διάβασμα για αυτήν την εργασία;

αφού δε θα τις διαβάσεις... 

[Δον]

κανείς;

Πρέπει να διαβάσεις την 7η (αν και ντάξει, δεν βοηθάει πάρα πολύ γιατί τα λέει πολύ περιληπτικά). Βρήκα αυτό http://www.slideshare.net/UmaRajaram/linear-prediction-18659179 το slide, το οποίο βασίζεται και αυτό στο βιβλίο του Haykin, αλλά έχει και δυο σχήματα με τα οποία μπορείς να καταλάβεις πιο εύκολα τι γίνεται.

Επίσης τα σκριπτάκια του Σισμάνη αν τα καταναοήσεις (για το πως να καλέσεις levinson durbin) θα σε βοηθήσει ιδιαίτερα.

[Silvo the Beautiful]

thanks για το site! είσαι ωραίος!

[Επικο Burger]

Παιζει να ανεβασει κανεις την εκφωνηση γιατι το ετημμυ ειναι ακομα down?

[Silvo the Beautiful]

εδώ!

[TrueForce]

αυτο για τις 3 σελιδες της αναφορας το τήρησε ποτε κανεις; :Ρ

[opcode]

Εφόσον χρησιμοποιούμε έναν γραμμικό προβλέπτη που παίρνει σαν είσοδο μια καθυστερημένη κατά Δ έκδοση του d(n) παράγει έξοδο y(n) και την συγκρίνει με το ίδιο το d(n) δεν θα πρέπει η έξοδος του φίλτρου να ακολουθεί το d(n)? Πως γίνεται να καθαρίσει τον θόρυβο;

[Δον]

Εφόσον χρησιμοποιούμε έναν γραμμικό προβλέπτη που παίρνει σαν είσοδο μια καθυστερημένη κατά Δ έκδοση του d(n) παράγει έξοδο y(n) και την συγκρίνει με το ίδιο το d(n) δεν θα πρέπει η έξοδος του φίλτρου να ακολουθεί το d(n)? Πως γίνεται να καθαρίσει τον θόρυβο;

Πρέπει να προσέξεις λίγο. Δεν έχεις θόρυβο. Έχεις περιοδική παρεμβολή. Που σημαίνει ότι το ανεπιθύμητο σήμα επαναλαμβάνεται περιοδικά. Έτσι, κάνοντας πρόβλεψη/εκτίμηση του σήματος, αυτό που παίρνεις μετά την εκτήμηση έχει ουσιαστικά δύο συνιστώσες, το σήμα σου και την (σχεδόν) ίδια περιοδική παρεμβολή.
Παίρνωντας το residual τώρα ( e = d - y ) έχεις αφαιρέσει σίγουρα την επίδραση της παρεμβολής από το σήμα. Είναι όπως ο διαφορικός ενισχυτής αφαιρεί την common mode τάση.

Toυλάχιστον κάπως έτσι το αντιλαμβάνομαι εγώ.

[opcode]

Εφόσον χρησιμοποιούμε έναν γραμμικό προβλέπτη που παίρνει σαν είσοδο μια καθυστερημένη κατά Δ έκδοση του d(n) παράγει έξοδο y(n) και την συγκρίνει με το ίδιο το d(n) δεν θα πρέπει η έξοδος του φίλτρου να ακολουθεί το d(n)? Πως γίνεται να καθαρίσει τον θόρυβο;

Πρέπει να προσέξεις λίγο. Δεν έχεις θόρυβο. Έχεις περιοδική παρεμβολή. Που σημαίνει ότι το ανεπιθύμητο σήμα επαναλαμβάνεται περιοδικά. Έτσι, κάνοντας πρόβλεψη/εκτίμηση του σήματος, αυτό που παίρνεις μετά την εκτήμηση έχει ουσιαστικά δύο συνιστώσες, το σήμα σου και την (σχεδόν) ίδια περιοδική παρεμβολή.
Παίρνωντας το residual τώρα ( e = d - y ) έχεις αφαιρέσει σίγουρα την επίδραση της παρεμβολής από το σήμα. Είναι όπως ο διαφορικός ενισχυτής αφαιρεί την common mode τάση.

Toυλάχιστον κάπως έτσι το αντιλαμβάνομαι εγώ.

Κοίτα και εγώ αυτό που κατάλαβα από την εκφώνηση είναι ότι πρέπει να αφαιρέσουμε την επίδραση του x(n), που λογικά είναι η περιοδική παρεμβολή αφού το άλλο σήμα (v(n)) είναι λευκός θόρυβος και δεν εμφανίζει περιοδικότητα. Στη συνέχεια όμως ορίζει σαν επιθυμητό σήμα το d(n) που είναι άθροισμα του σήματος "πληροφορίας" και της περιοδικής παρεμβολής. Αυτό το σήμα τώρα δίνεται σαν είσοδο και στον προβλέπτη αφού υποστεί καθυστέρηση κατά Δ δείγματα (δηλαδή στο MATLAB η είσοδος του φίλτρου είναι απλά u = [zeros(D,1);d(1:1:end-D] ;όπου D=10).

Το πρόβλημά μου είναι πως μάλλον δεν υπολογίζω σωστά τους βέλτιστους συντελεστές Wiener του φίλτρου καθώς μετά την εφαρμογή του Wiener φίλτρου στο σήμα εισόδου, το σήμα εξόδου του φίλτρου ακολουθεί τις ίδιες μςταβολές με το x (ή το d) αλλά διαφέρει σημαντικά στο πλάτος. Μάλιστα μου φαίνεται πως μοιάζει περισσότερο στο d παρά στο x.

Οι εντολές για υπολογισμό των πινάκων R και p που χρησιμοποίησα είναι οι εξής:
p=xcorr(u,d);
p=p(length(u):(length(u)+M));
p=p/length(u);

r=xcorr(u,u);
r=r(length(u):(length(u)+M));
r=r/length(u);
R=toeplitz(r,r);
wo=R\p;  

Τώρα ή κάνω κάποιο λάθος στον υπολογισμό του διανύσματος εισόδου u(n) ή στον υπολογισμό των R,p καθώς ο LD αλγόριθμος μου δίνει φίλτρο που ουσιαστικά σχεδόν ταυτίζεται με τη βέλτιστη απόκριση του φίλτρου Wiener. Δηλαδή, ακολουθεί την βέλτιστη απόκριση Wiener όπως θα έπρεπε άλλωστε. Το πρόβλημα είναι πως η απόκριση Wiener δεν είναι καλή. Κάνω καμία πατάτα στις παραπάνω εντολές που δεν τη βλέπω;      

[Σα τανυστής]

αυτο για τις 3 σελιδες της αναφορας το τήρησε ποτε κανεις; :Ρ

αμα εισαι gangster και γραφεις latex 2 στηλες με μικρη γραμματοσειρα μπορει και να το πετυχεις.

μεχρι τοτε εγω θα συνεχισω να στελνω 4-7 σελιδες 

[xaris]

Εφόσον χρησιμοποιούμε έναν γραμμικό προβλέπτη που παίρνει σαν είσοδο μια καθυστερημένη κατά Δ έκδοση του d(n) παράγει έξοδο y(n) και την συγκρίνει με το ίδιο το d(n) δεν θα πρέπει η έξοδος του φίλτρου να ακολουθεί το d(n)? Πως γίνεται να καθαρίσει τον θόρυβο;

Πρέπει να προσέξεις λίγο. Δεν έχεις θόρυβο. Έχεις περιοδική παρεμβολή. Που σημαίνει ότι το ανεπιθύμητο σήμα επαναλαμβάνεται περιοδικά. Έτσι, κάνοντας πρόβλεψη/εκτίμηση του σήματος, αυτό που παίρνεις μετά την εκτήμηση έχει ουσιαστικά δύο συνιστώσες, το σήμα σου και την (σχεδόν) ίδια περιοδική παρεμβολή.
Παίρνωντας το residual τώρα ( e = d - y ) έχεις αφαιρέσει σίγουρα την επίδραση της παρεμβολής από το σήμα. Είναι όπως ο διαφορικός ενισχυτής αφαιρεί την common mode τάση.

Toυλάχιστον κάπως έτσι το αντιλαμβάνομαι εγώ.

Κοίτα και εγώ αυτό που κατάλαβα από την εκφώνηση είναι ότι πρέπει να αφαιρέσουμε την επίδραση του x(n), που λογικά είναι η περιοδική παρεμβολή αφού το άλλο σήμα (v(n)) είναι λευκός θόρυβος και δεν εμφανίζει περιοδικότητα. Στη συνέχεια όμως ορίζει σαν επιθυμητό σήμα το d(n) που είναι άθροισμα του σήματος "πληροφορίας" και της περιοδικής παρεμβολής. Αυτό το σήμα τώρα δίνεται σαν είσοδο και στον προβλέπτη αφού υποστεί καθυστέρηση κατά Δ δείγματα (δηλαδή στο MATLAB η είσοδος του φίλτρου είναι απλά u = [zeros(D,1);d(1:1:end-D] ;όπου D=10).

Το πρόβλημά μου είναι πως μάλλον δεν υπολογίζω σωστά τους βέλτιστους συντελεστές Wiener του φίλτρου καθώς μετά την εφαρμογή του Wiener φίλτρου στο σήμα εισόδου, το σήμα εξόδου του φίλτρου ακολουθεί τις ίδιες μςταβολές με το x (ή το d) αλλά διαφέρει σημαντικά στο πλάτος. Μάλιστα μου φαίνεται πως μοιάζει περισσότερο στο d παρά στο x.

Οι εντολές για υπολογισμό των πινάκων R και p που χρησιμοποίησα είναι οι εξής:
p=xcorr(u,d);
p=p(length(u):(length(u)+M));
p=p/length(u);

r=xcorr(u,u);
r=r(length(u):(length(u)+M));
r=r/length(u);
R=toeplitz(r,r);
wo=R\p;  

Τώρα ή κάνω κάποιο λάθος στον υπολογισμό του διανύσματος εισόδου u(n) ή στον υπολογισμό των R,p καθώς ο LD αλγόριθμος μου δίνει φίλτρο που ουσιαστικά σχεδόν ταυτίζεται με τη βέλτιστη απόκριση του φίλτρου Wiener. Δηλαδή, ακολουθεί την βέλτιστη απόκριση Wiener όπως θα έπρεπε άλλωστε. Το πρόβλημα είναι πως η απόκριση Wiener δεν είναι καλή. Κάνω καμία πατάτα στις παραπάνω εντολές που δεν τη βλέπω;      

μπορεί να είναι χαζή η ερώτηση αλλά για ποιόν λόγο διαιρείς με το length(u) ??