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

[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 με προσοχή στα ονόματα των αρχείων!

[Dimitris1989]

Για το πρώτο κομμάτι της άσκησης (πριν μπλέξουμε με το sound.mat), για να φτιάξουμε τα διανύσματα που ζητάει θα ορίσουμε εμείς έναν πίνακα που θα έχει όλα τα n;

Ο τύπος του u(n) είναι σωστός; Έχουμε u(n) και στα δύο μέλη οπότε εκτελώντας τις πράξεις προκύπτει u(n-1)=0

[Nerevar]

Για το πρώτο κομμάτι της άσκησης (πριν μπλέξουμε με το sound.mat), για να φτιάξουμε τα διανύσματα που ζητάει θα ορίσουμε εμείς έναν πίνακα που θα έχει όλα τα n;

Ο τύπος του u(n) είναι σωστός; Έχουμε u(n) και στα δύο μέλη οπότε εκτελώντας τις πράξεις προκύπτει u(n-1)=0

Είναι u(n)=-0.25*u(n-1)+v(n)

[Nerevar]

Eχει βρει κανένας πίνακες R,p για αλληλοεπιβεβαίωση (sic)?

[Dimitris1989]

Δεν το πρόσεξα 

Για τους πίνακες που αναφέρεις δε θα πρέπει να έχουμε κάποιες τιμές για να τους υπολογίσουμε;
Το διάνυσμα v(n) του θορύβου θα το φτιάξουμε εμείς (το matlab δλδ) με βάση τη διακύμανση;

[Nerevar]

Δεν το πρόσεξα 

Για τους πίνακες που αναφέρεις δε θα πρέπει να έχουμε κάποιες τιμές για να τους υπολογίσουμε;
Το διάνυσμα v(n) του θορύβου θα το φτιάξουμε εμείς (το matlab δλδ) με βάση τη διακύμανση;

Δες τα scripts που έχει στο ετημμυ.

_________________________________________________ _________________________________________________ ______________________________


Ο πιτσιάνης τρολάρει επειδή του είπα οτι ο σισμάνης δεν απαντάει στα emails και δεν έχει ανεβάσει ακόμα τις ασκήσεις...

[Dimitris1989]

Σε κάποια .m παιρνει
v = sqrt(σ)*randn(n,1); v = v - mean(v);
ενώ σε άλλα
v = σ*randn(n,1); v = v - mean(v);

Ποιό από τα δύο είναι σωστό; Επίσης για να υπολογίσουμε τον πινακα u(n) δε θα πρέπει να έχουμε την αρχική τιμή u(0);

[Nerevar]

Σε κάποια .m παιρνει
v = sqrt(σ)*randn(n,1); v = v - mean(v);
ενώ σε άλλα
v = σ*randn(n,1); v = v - mean(v);

Ποιό από τα δύο είναι σωστό; Επίσης για να υπολογίσουμε τον πινακα u(n) δε θα πρέπει να έχουμε την αρχική τιμή u(0);

Είναι sqrt(variance)

[Dimitris1989]

Η αυτοσυσχέτιση πως θα βγει;
Στο .m file την έχει υπολογίσει απο πριν και βάζει τις τιμές σε ενα 2x2 πίνακα.

Με την autocorr του matlab προκύπτει ένα μονοδιάστατο διάνυσμα

Το pdf με τα φιλτρα wiener δίνει τον τύπο R=E[u(n)*u'(n)]. Κάνω τον πολλαπλασιασμό του u με τον ανάστροφό του, παίρνω το mean() και προκύπτει ένα διάνυσμα 1xn αριθμούς.

Πως καταλήγει στους 2 αριθμούς;

[Nerevar]

u=[1;2;3;4];
z=autocorr(u);
Ruu=toeplitz(z);

[Dimitris1989]

Σκέφτηκα το εξής:
Εφόσον ζητάει φίλτρο 2 συντελεστών άρα Μ=2. Άρα σύμφωνα με τις σελίδες 3-4 των διαφανειών για τα wiener θα είναι
R=[r(0) r(1);r(1) r(0)]. Όπου r(0) η αυτοσυσχέτηση όταν lag=0 (σελ 3), το οποίο βγάζει πάντα 1 και r(1) όταν lag=1.

Δηλαδή, τρέχω την autocorr(u) και κρατάω τα δύο πρώτα αποτελέσματα.
Ομοίως με την crosscorr για το P.

Είναι σωστή σαν σκέψη; Το r(1) μου το βγάζει αρνητικό

[Nerevar]

Για ποιο λόγο να μην κάνεις το πρώτο ερώτημα με πράξεις στο χέρι αφού είναι αρκετά απλό?Ο πιτσιάνης δε διαφωνεί.Εμένα μου βγήκαν:

R=[0.32 0.08        p=[0.34  
     0.08 0.32]            0]


Για τον p δεν είμαι και πολύ σίγουρος,αλλά τα όρια ευστάθειας είναι 0<μ<5,κάτι που επιβεβαιώνεται και από την προσομοίωση.

_________________________________________________ _________________________________________________ _______________________________________

Για το δεύτερο κομμάτι της εργασίας,ο πίνακας R μου βγαίνει 21x21,ο p βγαίνει 26682624x1,υποτίθεται οτι πρέπει να κρατήσουμε μόνο τα πρώτα 21 στοιχεία του?

[Dimitris1989]

Με το χέρι πως το κάνεις; Εννοείς τους τύπους στη σελίδα 4 από το optWienerFilt.pdf;

[Nerevar]

Με το χέρι πως το κάνεις; Εννοείς τους τύπους στη σελίδα 4 από το optWienerFilt.pdf;

Δες το παράδειγμα στο pdf που λες και τις λυμμένες ασκήσεις που ανέβηκαν σήμερα.

[opcode]

Για ποιο λόγο να μην κάνεις το πρώτο ερώτημα με πράξεις στο χέρι αφού είναι αρκετά απλό?Ο πιτσιάνης δε διαφωνεί.Εμένα μου βγήκαν:

R=[0.32 0.08        p=[0.34 
     0.08 0.32]            0]


Για τον p δεν είμαι και πολύ σίγουρος,αλλά τα όρια ευστάθειας είναι 0<μ<5,κάτι που επιβεβαιώνεται και από την προσομοίωση.

_________________________________________________ _________________________________________________ _______________________________________

Να σε ρωτήσω γιατί έχω σκαλώσει. Εγώ παίρνω την σχέση u(n)=-0.25*u(n-1)+v(n) και αφου πολλαπλασιάσω με u(n-k) παίρνω την μέση τιμή και στα 2 μέλη και μου προκύπτει : ruu(k)=-0.25*ruu(k-1)+rvu(k). Αυτό το ruv(k) από πού υπολογίζεις τισς τιμές τού; Γιατί εγώ έβγαλα R=[0.3613     0.17
                         0.17     0.3613]
και νομίζω πως σε αυτό το σημείο κάνουμε κάτι διαφορετικό.

[Nerevar]

Για ποιο λόγο να μην κάνεις το πρώτο ερώτημα με πράξεις στο χέρι αφού είναι αρκετά απλό?Ο πιτσιάνης δε διαφωνεί.Εμένα μου βγήκαν:

R=[0.32 0.08        p=[0.34 
     0.08 0.32]            0]


Για τον p δεν είμαι και πολύ σίγουρος,αλλά τα όρια ευστάθειας είναι 0<μ<5,κάτι που επιβεβαιώνεται και από την προσομοίωση.

_________________________________________________ _________________________________________________ _______________________________________

Να σε ρωτήσω γιατί έχω σκαλώσει. Εγώ παίρνω την σχέση u(n)=-0.25*u(n-1)+v(n) και αφου πολλαπλασιάσω με u(n-k) παίρνω την μέση τιμή και στα 2 μέλη και μου προκύπτει : ruu(k)=-0.25*ruu(k-1)+rvu(k). Αυτό το ruv(k) από πού υπολογίζεις τισς τιμές τού; Γιατί εγώ έβγαλα R=[0.3613     0.17
                         0.17     0.3613]
και νομίζω πως σε αυτό το σημείο κάνουμε κάτι διαφορετικό.

Στη σχέση που έχεις βάλε k=0,1,τότε rvu(0)=sv² και rvu(1)=0,εκτός και αν κάνω κάποιο τραγικό λάθος,βγαίνει

ruu(0)=-0.25*ruu(1)+sv²
ruu(1)=-0.25*ruu(0)