THMMY.gr

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

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

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

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

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

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

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

_________________________________________________ _________________________________________________ ______________________________


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

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

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

Είναι sqrt(variance)

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

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

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

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

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

Σκέφτηκα το εξής:
Εφόσον ζητάει φίλτρο 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) μου το βγάζει αρνητικό

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

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


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

_________________________________________________ _________________________________________________ _______________________________________

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

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

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

Να σε ρωτήσω γιατί έχω σκαλώσει. Εγώ παίρνω την σχέση 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)

Ναι εγώ αρχικά το χα υπολογίσει έτσι όπως περιγράφεις και μου έδωσε αποτελέσματα κοντά στα δικά σου (για την ακρίβεια r0=0.34 και r1=0.085). Απλά δεν είμαι σίγουρος κατα πόσο ισχύει η σχέση αυτή, οτι δηλαδή ορίζεται η ετεροσυσχέτιση μετάξύ των v και u μόνο για μηδενική καθυστέρηση (δεν εβγαλα άκρη από το δίδυμο Πιτσιάνη Σισμάνη  :P ) και αποφάσισα να πάω και να υπολογίσω το rvu μέσω του rvv οπου για την αυτοσυσχέτιση ισχύει σίγουρα καθώς είναι λευκός θόρυβος.

Το αναφέρω αυτό με κάθε επιφύλαξη καθώς ο Σισμάνης με μπέρδεψε. Στην αρχή έλεγε πως ο θόρυβος είναι ασυσχέτιστος με κάθε σήμα και στην συνέχεια εμφάνιζε περιπτώσεις οπου υπήρχε συσχέτιση.  :D
Επειδή άκρη από δαύτους δεν πρόκειται να βγάλω και δεν έχω χρόνο να ψάξω όλο το βιβλίο του Haykin, ειμαι σε φάση που δεν ξέρω τι είναι σωστο και τι όχι.  :D

Αν βάλεις k=0 τότε είναι ruu(0)=-0.25*ruu(-1)+sv². (σου έφυγε ενα - )
Το ruu(-1) είναι 0;

Μπορεί αυτό να είναι πιο σωστό,ποιος ξέρει :D


Η αυτοσυσχέτιση είναι άρτια συνάρτηση,ισχύει Ruu(-k)=Ruu(k).Για το άλλο μάλλον μπέρδεψες πάλι τους δείκτες.Ας ρωτήσει κάποιος τον πιτσιάνη σήμερα (αφού ο σισμάνης δεν απαντάει ποτέ..)

Εγώ γενικά δεν έχω καταλάβει τον παλμό της εργασίας .Πως πρέπει να ξεκινήσουμε? Φτιάχνουμε μέσω matlab ένα διάνυσμα για τον θόρυβο όπως κάνει αυτός στα  scriptakia? φτιάχνουμε και το διάνυσμα u που μπαίνει στο φίλτρο και κάνουμε autocorr ? και αν ναι πόσο πρέπει να είναι το n,δηλαδή τι μέγεθος είναι το διάνυσμα του θορύβου και της εισόδου στο φίλτρο αντίστοιχα?

από αυτούς που ήδη παλεύουν ή πάλεψαν με την εργασία , πόσο χρόνο "παίρνει"  ξέρουμε;

Δεν την εχω κανει ακόμη, αλλά αφού ο Πιτσιάνης είπε ότι είναι "εύκολη", υπολόγισε κανένα 3ήμερο…

Πάλι τα νούμερα δε μου βγαίνουν.Βασικά μου βγαίνει το ruu(1) αρνητικό. Λύνω το παραπάνω σύστημα ως εξής:
Αντικαθιστώ στην πρώτη εξίσωση το ruu(1) και λύνοντας βγαίνει ruu(0)=0.36. Στη συνέχεια βαζω το 0.36 στη δεύτερη εξίσωση οπότε βγαίνει ruu(1)= -0.09

Είπε σήμερα τίποτα βοηθητικό σχετικά με την εργασία;

Εγώ αυτό που έκανα με την ετεροσυσχέτιση rvu είναι rvu(k)=E[v(n)*u(n-k)] οπότε βγαίνει rvu(k)=-0.25rvu(k+1)+σ^2*δ(k) . Μου φαίνεται το πιο λογικο καθως το v(n) εξαρτάται από το u(n)

Εγώ το έκανα και με τους 2 τρόπους, δηλαδή και θεωρώντας ότι u και v είναι ασυσχέτιστα μεταξύ τους για lag>0 και ότι συσχετίζονται και βρήκα τα ίδια αποτελέσματα για τους βέλτιστους συντελεστές του φίλτρου Wiener. Βέβαια ο πίνακας αυτοσυσχέτισης και η στήλη ετεροσυσχέτισης μου βγήκαν διαφορετικοί σε κάθε περίπτωση.

και μετα πως το λυνεις αυτό?

Κι εμένα τόσο μου βγαίνουν τα νούμερα.. Υποθέτοντας ότι δεν ορίζεται ετεροσυσχέτιση μεταξύ των v και u
για lag>0 και ότι η αυτοσυσχέτιση του u είναι άρτια. Άρα κάτι δεν πάει καλά... Powerslave δε γίνεται να βγάζεις
ruu(1) = -0.25ruu(0) και να σου βγαινουν και τα 2 θετικα..

Καλά όλα αυτά, με το μουσικό κομμάτι έχει βγάλει κανείς καμιά άκρη;

Καμία απολύτως. Δε θα πρέπει να παίξουμε το d-y?

Η λογική αυτό λέει. Εγώ πάντως δεν είχα αποτελέσματα ακόμα.

Κι εγώ το παλεύω με το τραγούδι.
Όταν εφαρμόζετε Steepest Descend για το σήμα x με τα ημίτονα
σας βγαίνει όντως το e(n) = x(n) όταν επιτευχθεί σύγκλιση ?
Γιατί στην εκφώνηση λέει ότι πρέπει να παίρνουμε το καθαρό σήμα σαν το σφάλμα,
άρα το ίδιο πρέπει να συμβαίνει και για να πάρουμε το καθαρό από θόρυβο τραγούδι.

Εμένα μου βγαίνει x(n) με κάποια μικρή παραμόρφωση. Αλλά με τον ήχο καμία τύχη.

Και σε μένα το ίδιο είναι. Απλά ο λευκός θόρυβος που προστίθεται στο x είναι έτσι κι αλλιώς πολύ μικρής τιμής.
Πώς ξέρουμε αν παίρνουμε μια καλή προσέγγιση της x (εφόσον διαφέρει όντως) ή
απλά παίρνουμε ότι να ναι ? :/

Φαντάζομαι μπορείς να δεις αν το e(n) είναι πιο κοντά στο x(n), απ ότι το το d(n).

Για το x(n) πόσα δείγματα βάζετε; Γιατί επηρεάζει αρκετά το αποτέλεσμα

Πείτε μου ότι δεν το κάνουμε στο χέρι!!

Σπαμ για παράταση ώστε να έχουμε και το Σ/Κ?

Μπορούμε να το ζητήσουμε τ Τεταρτη...μπορεί ν πει κάποιος τ είπε τ Δευτέρα για την εργασία ???ρωτήσατε απορίες...δυστυχως εγώ ήρθα όταν είχαν τελειώσει :S

Εγώ έβαλα n=1000. Την αυτοσυσχέτιση για το πρώτο ερώτημα θα την κάνεις με το χέρι σίγουρα,
για το τραγούδι αριθμητικά όπως λέει και η εκφώνηση (με Matlab δλδ).
Μπορεί να επιβεβαιώσει κάποιος τον πίνακα αυτοσυσχέτισης για την περίπτωση του
τραγουδιού ? Βρήκα Ru = [ 0.4951 0.1930; 0.1930 0.4951]. Το τραγούδι δεν παίζει βέβαια...

Γιατί μόνο 2x2 ο πίνακας αυτοσυσχέτισης?Mε

z=autocorr(u);
Ruu=toeplitz(z);

όπου u το διάνυσμα που δίνεται από το ethmmy σε .mat extension,εμένα μου βγαίνει 21x21 και ο p nx1 (όπου n το μήκος των u,d)

(δεν μπορώ να κάνω πράξεις μεταξύ των R,p,αλλά αυτό μάλλον είναι άλλο θέμα)

Εφόσον έχεις 2 συντελεστές w ο R πρέπει να είναι 2x2 και ο p το ίδιο (2x1).
Με την autocorr βγαίνει όπως το λες, αλλά είναι κανονικοποιημένος.
Για να μαι σίγουρος το έκανα και με πολλαπλασιασμούς, όπως μας το πε στην τάξη, αλλά φαντάζομαι
αν "ξεκανονικοποιήσεις" τον R που βγάζει η autocorr και πάρεις μόνο τα 2 πρώτα στοιχεία με toeplitz
θα βρεις το ίδιο..

 Και εγω ιδιο με τον Aloko το βγαζω τον πινακα αυτοσυσχετισης αλλα τραγουδι δεν ακουω και εχω δοκιμασει ποσα διαφορετικα μ ..

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

Μα και στο δεύτερο λέει να υλοποιήεις φίλτρο 2 συντελεστών. Εννοείς ότι ίσως χάνουμε ακρίβεια αν πάρουμε 2 συντελεστές (άρα 2x2 πίνακα R) ? Το τραγούδι το παίρνω ως d-y. nx1 δηλαδή. Αλλά δεν έχει φιλτραριστεί ο θόρυβος.

Ντοινγκ

Ναι τώρα το είδα :D :D

Το ερώτημα με το μουσικό κομμάτι υποτίθεται πως θα το κάνουμε χρησιμοποιώντας βέλτιστους συντελεστές wiener ή τον αλγόριθμο steepest descent;

Eχεις τα σήματα u,d από το ethmmy,βρίσκεις R,p τα κολλάς στο script και υποτίθεται οτι θα ακούσεις το κομμάτι.

Τα έκανα όλα αυτά, βρήκα τα αποτελέσματα που αναφέρονται και παραπάνω και είμαι στο φάση που υποτίθεται ότι ακούω το κομμάτι.

Κάπου εκεί σκαλώσαμε όλοι :P

Ξανακατεβαστε τα αρχεια  του ηχου και θορυβου γιατι ειχαν προβλημα.

Πλάκα κάνεις...

Πήγαμε στον Σισμάνη το πρωί και κατάλαβε ότι έχει λάθος! Οπότε με τα νέα αρχεία σε εμένα τρέχει μία χαρά με τον πίνακα αυτοσυσχέτισης που έχει και ο Αλόκο .

Ηλίθιοι.

Μπράβο ρε παιδιά, φαινόταν ότι κάτι δεν πάει καλά..

Για τον πίνακα αυτοσυσχέτισης χρησιμοποιείται τις εντολές π είπατε πριν???σ μένα στην είσοδο γιατί μ βγαζει διαφορετικό αποτέλεσμα με την autocorr XD

Το ανέφερα και νωρίτερα, η autocorr κάνει κανονικοποίηση στο αποτέλεσμα ώστε r(0)=1.
Νομίζω ότι αν πολλαπλασιάσεις με την διακύμανση var(u) τον πίνακα Ru που βγάζει η autocorr θα πάρεις το σωστό αποτέλεσμα. Διαφορετικά μπορείς να υπολογίσεις τον R με πολλαπλασιασμό πινάκων.
Αυτό πρότεινε και ο Σισμάνης στο μάθημα.

Εγω πλεον με τα καινουργια αρχεια δεν μπορω να τα κανω load πια .μου βγαζει σφαλματα.καμια ιδεα για το το παιζει?

Χωρίς να πεις τι σφάλμα σου βγάζει,λίγο δύσκολο.

Ξανακατέβαστα....δε παιζει να μην μπορεις να τα βάλεις :P

Error using load
Unknown text on line number 1 of ASCII file C:\Users\Dre\Downloads\noise (1).mat
"MATLAB".

Error in uiimport/runImportdata (line 447)
                    datastruct = load('-ascii', fileAbsolutePath);

Error in uiimport/gatherFilePreviewData (line 415)
        [datastruct, textDelimiter, headerLines]= runImportdata(fileAbsolutePath, type);

Error in uiimport (line 231)
[ctorPreviewText, ctorHeaderLines, ctorDelim] = ...

αυτά τα errors μου βγαζει .....τι φαση

εχει σημασια πως ονομάζεται το αρχειο ????import το κανω και μου βγαζει αυτά τα σφαλματα

Απλά αγνόησε αυτό που έγραψα.

P πόσο βγάζετε;

ασχετο, ζητησατε παραταση?

παμ παμ παμ; παραμ; τι έγινε;

Έχω ζητήσει εδώ και 2 μέρες.Καμία απάντηση.

Σήμερα νομίζω ανέφερε ότι θα πάρει παράταση μερικές μέρες..

Σας τρέχει με τις νέες εκδόσεις ? Σε μένα υπάρχει βελτίωση, ακούγεται το τραγούδι,
αλλά και πάλι ακούω και αρκετό θόρυβο.  :(

To p είναι [0.42; 0] από την εκφώνηση του ερωτήματος. Είπε ότι θα αφήσει ανοιχτά την υποβολή εργασιών μέχρι Δευτέρα-Τρίτη.

@aloko αν ο πίνακας εξόδου του φίλτρου y έχει στην πρώτη θέση μηδέν κάντο y(1) = []; και βάλε μηδενικό στο τέλος του πίνακα.

Σε ευχαριστώ Neal, όντως αυτό ήταν το πρόβλημα !

Παιδιά έχω βγάλει κι εγώ τον ίδιο πίνακα με τον aloko, για το τραγούδι, αλλά όταν εφραμόζω τον αλγόριθμο, απλά το y(n) πλησιάζει όσο μπορεί το d(n), και πρακτικά μηδενίζει το σφάλμα. Έχετε ιδέα τι μπορεί να κάνω λάθος ; Επίσης, απ΄ότι κατάλαβα το κομμάτι είναι δικαναλικό , πρέπει πρώτα να χωρίσω τα κανάλια , και για το καθένα να βρω τα R,p και να εφαρμόσω το φίλτρο, ή σε όλο μαζί το αρχείο;

Εdit: E βασικά το έλουσα, είχα ονομάσει το σφάλμα e, που έτσι ονομάζει αυτός στον κώδικα το παραμετρικό σφάλμα, οπότε έβαζα να ακούσω το παραμετρικό σφάλμα.Να μάθω να κάνω copy paste χωρίς να τα βλέπω ^banghead^

Απο ethmmy

    

1η εργασία (Hμ/νία λήξης δηλώσεων: 15/04/2014)

 ;D  ;D  ;D  ;D

 ^rotate^

https://www.youtube.com/watch?v=r2S1I_ien6A

Για όποιον θέλει να σταματήσει το τραγούδι,η εντολή sound(music,fs) δε μπαίνει στο script.Aν γίνει save και μετά load (ως music.mat),τότε:


player = audioplayer(music,fs)

% start the playback
play(player);

% pause the playback
pause(player);

% resume the playback
resume(player)

% stop the playback
stop(player)

Άρα έχουμε παράταση;

Έχει ακούσει κανείς όλο το τραγούδι; Εγω που δοκιμασα, δεν ακουω καθολου θορυβο, αλλά σε πολλά σημεία στο background ακούγεται σαν να παίζει ταυτόχρονα άλλο κομμάτι του τραγουδιού.  :D

https://www.youtube.com/watch?v=l_tbHLYsLBM

invalid utube link του παντελο???FFS!

Αν δν ακούς αυτο καθαρα τότε λάθος εργασία έκανες....


Το κομμάτι το ακούω καθαρα παντως γω!

Παιδιά ανεβάσατε και το τραγούδι;

Για κάποιον που δεν έχει ασχοληθεί ποτέ με Matlab, τι θα προτείνατε να διαβάσει προκειμένου να διεκπεραιώσει τις εργασίες του μαθήματος;

Είναι 200ΜΒ οπότε όχι...

Επειδή μάλλον κάτι κάνω πολύ λάθος και μου κολλάει το MATLAB, εσείς πως βρίσκεται τον p?

Ομοίως για το άλλο p [0.34; 0]

Γενικά αμα μπεις στο κλίμα των διαλέξεων και των φίλτρων θα μπορείς να την κάνεις την πρώτη εργασία τουλαχιστον...άλλωστε τους κώδικες τους δίνει στις διαλέξεις αντίστοιχα....

Α μάλιστα, ευχαριστώ, θα το ψάξω  :)

Μιλάς για τα σκριπτ που έχει στο ethmmy? Τα ar, arma, ma?

Eίναι το gradWienerFilt.m

Γιατί δεν το βρίσκω ούτε στο υλικό ούτε στην εργασία...;  :P

Ψάξε καλύτερα,εκεί είναι.

οταν κανετε plot x(n),e(n) εχουν μεγαλες διαφορες?

Αρχικά ούτε καν, δηλαδή ήταν αρκετά κοντά, αλλά αν σκεφτείς ότι έτσι κι αλλιώς ο θόρυβος έχει πολύ μικρή τιμή,
πρέπει να ταυτίζονται σχεδόν τα σήματα. Απλά πρέπει να προσέξεις ότι το τελικό σήμα έχει τιμή μηδέν στην πρώτη χρονική στιγμή. Αυτή την τιμή την διαγράφεις από τον πίνακα και κάνεις μηδέν την τελευταία. Το περιγράφει και συνάδελφος σε προηγούμενο post. Μόλις το έκανα αυτό ταυτίστηκαν τα σήματα. Ουσιαστικά ήταν κάτι σαν lag.

R=[0.36 -0.09;-0.09 0.36]  p=[0.34 0] αυτα ειναι που βρηκατε ?
γτ καπου τ χυνω ,μ βγαινει τελειως  λαθος

Σωστά φαίνονται.

Κι εγώ αυτά έχω. Δες το θέμα με το μηδέν στην έξοδο.

1000 ty αυτο ηταν!..οποιος τ βρηκε ειναι σκυλι

Μιας και η ημερομηνια παραδοσης των εργασιων περασε, μπορει καποιος να μου στειλει το τελευταιο ερωτημα της 1ης (που να τρεχει/ειναι σωστο); Ή αν ειναι δυνατον και ολοκληρη την εργασια για λογους μελετης