THMMY.gr

Ανέβηκε  :D :D
Καταληκτική Ημερομηνία 24-3-2017

Τα αρχεία sound.mat και noise.mat που αναφέρει στη τελευταία παράγραφο, που βρίσκονται;;;; :???:

Στην εκφώνηση λεει ότι η παράδοση είναι μέχρις την παρασκευή 24/3, ενώ στο elearning μέχρι την κυριακή 26/3; Θεωρούμε ότι μας συμφέρει;

εχει παει κανεις να κανει το quiz, εγω οταν παταω υποβολη δε μου βγαζει βαθμο, συμβαινει κ σε εσας?

μπα τίποτα δεν βγάζει... βασικά αυτό το quiz τι ακριβώς είναι γιατί κάπου χάθηκα έχει κάποια σχέση με την εργασία;

να σου πω την αληθεια δε ξερω τι ακριβως ρολο βαρανε, ισως ειναι πχ για δικη μας εξασκηση??

το quiz ειναι καθαρα για δικη μας εξασκηση και παιρνεις προσθετικο βαθμο
δεν ειναι υποχρεωτικο

Νομίζω ο Σισμάνης δεν είπε ότι είναι για προσθετικό βαθμό απαραίτητα...Είπε απλά ότι δε ξέρει και θα ρωτήσει τον Πιτσιάνη...

τι φάση εργασία ;; Πόσες θα έχουμε γενικά ; και τι χρόνο θέλουν ;;

Αν κατάλαβα καλά η εργασίες θα είναι 3, προσθετικές στο βαθμό από τις τελευταίες εξετάσεις
Δεν φάνηκε η πρώτη τουλάχιστον δύσκολη, απλά ένα διάβασμα θεωρία μέχρι την 3η σειρά διαφανειών, ασκήσεις και μια ματιά σε κάτι κώδικες που έχει ανεβασμένους στο ethmmy ... γενικά βγαίνει πιστεύω εύκολα

να ρωτησω κατι ισως χαζο

στο δευτερο ερωτημα που θελει να υπολογισουμε το ελαχιστο τετραγωνικο σφαλμα, θα παρουμε τον τυπο 18 απο τις διαφανειες - Lecture 2??

Εγώ από εκεί το πήρα

οκ ευχαριστωωω :D :D

 Αυτό που λέει να υλοποιήσουμε την διάταξη του σχήματος σε MAtlab , εννοεί μέσω του simulink ?

Όχι. Απλός κώδικας που υλοποιεί το σύστημα.

Το e(n) πρέπει να προσεγγίζει το x(n) ;

Εγώ αυτό καταλαβαίνω αφαιρούμε από το d(n) το θόρυβο οπότε θα πρέπει να προσεγγίζουμε το αρχικό σήμα χωρίς την επίδραση του θορύβου

To p και στο πρώτο και στο τελευταίο ερώτημα δεν είναι P=[σ^2 0 0]?

αυτό βρήκα αλλά στ διαγράμματα το e(n)δεν προσεγγίζει το x(n).Κάτι πρέπει να κάνουμε λάθος

Εμένα βγαίνει πάρα πολύ καλό,ταυτίζεται με το χ(ν),επομένως θα μπορούσα να πω ότι όσο αφορά το πίνακα δεν υπάρχει θέμα,αν έχεις χρησιμοποιήσει το s[0 0 u]  θα σου πρότεινα να το βγάλεις και να βάλεις απλά στην θέση του το u γιατί και εμένα με αυτό τον τρόπο έβγαινε άλλα ντ αλλων.

Όντως έφτιαξε! thanks

πώς μπορώ να ακούσω κάτι από το matlab??

sound(σήμα,συχνότητα)

https://www.mathworks.com/help/matlab/ref/sound.html (https://www.mathworks.com/help/matlab/ref/sound.html)

το δοκίμασα, αλλά μου βγάζει το εξής έρρορ:
Subscript indices must either be real positive integers or logicals.

λοιπον να ρωτησω

για να υπολογισετε την αυτοσυσχετιση του u στο ματλαβ, για το πρωτο κομματι τι κανετε?


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

κανεις καμια προταση?

δεν ξέρω τι κάνεις αλλά πήγα να το τρέξω 3 φορές και τις 3 κόλλησε το laptop  :P
Μπόρεσε κανείς να ακούσει τον ήχο;
Παίρνουμε το R_u και p που υπολογίσαμε από πριν (με διαφορετικό σv);

εκει π εχει nCoeff βάζω το 3 (θελουμε φιλτρο 3 συντελεστων)
εσυ την αυτοσυσχετιση πως την υπολογισες?

βρήκα στο χαρτί τις εξισώσεις για τα r0 r1 r2 και στην ματλαμπ έλυσα το γραμμικό σύστημα

αα ενταξει ετσι, εγω πηγα να το κανω απευθειας στο ματλαμπ
υπαρχει καποιος που υπολογισε την αυτοσυσχετιση κατευθειαν σε ματλαμπ?

@lady_of_winter
υπάρχει η autocorr που κάνει αυτή τη δουλειά
η σύνταξη είναι r = autocorr( signal, max_lag);

ένα σημείο που θέλει προσοχή είναι πως ο τύπος Ε( x(n)*x(n-k) ) δεν δίνει αυτοσυσχέτιση αλλά αυτοδιασπορά. η αυτοδιασπορά είναι ίση με την αυτοσυσχέτιση όταν η διακύμανση είναι 1 (δεν ισχύει για την εργασία). έτσι, η autocorr θα δίνει πάντα r(0) = 1 γιατί τόσο κάνει εξ' ορισμού. για περισσότερες πληροφορίες μπορεί κανείς να κοιτάξει τις διαφάνειες του κουγιουμτζή για τις χρονοσειρές ( http://users.auth.gr/dkugiu/Teach/TimeSeriesTHMMY/TimeSeriesNotes.pdf )
άμα βρίσκεις διαφορετικά αποτελέσματα έλενξε εάν το ένα ισούται (στο περίπου) με το άλλο πολλαπλασιασμένο με μία σταθερά

για την εργασία πιστεύω θέλουμε αυτοδιασπορά. παίζει να είναι λίγο φάουλ οι διαφάνειες

ωω ευχαριστω, ειχα δοκιμασει και το autocorr αλλα μου βγαζε λαθος αποτελεσμα και τωρα καταλαβα γιατι

θα το κοιταξω :) :)

Ερώτηση:
Το p είναι η ετεροσυσχέτιση της εξόδους του Wiener με την επιθυμητή έξοδο (u(n) με υ(n)) ή ετεροσυσχέτιση u(n) με d(n) απλώς? Και αν είναι το δεύτερο, πως βρίσκουμε την ετεροσυσχέτιση του x(n) με το u(n)?  :???:

+
Γενικά το desired θεωρούμε ότι είναι το d(n) ή η επιθυμητή έξοδος του Wiener, δηλαδή το υ(n)?

αν δε κανω λαθος ειναι το δευτερο που λες, ο πινακας θα εχει τιμες για p(0), p(-1) και p(-2) και απο αυτά μόνο το πρώτο νομίζω είναι μη μηδενικό και ισούται με την διακύμανση του λευκού θορύβου, αν δε ειναι σωστο ας με διορθωσει καποιος

Και στον υπολογισμό του, εκεί που πρέπει να υπολογίσεις την συσχέτιση μεταξύ u(n) και x(n) τι κάνεις? (r_ux(-k)+r_uυ(-k))

κάποιος για επαληθευση αν μπορει να μου πει ποσο υπολογισε με το χερι τα r(0),r(1),r(2) :)

εγω τα χω βγαλει 0,336 0,075 και 0,02158
αλλα δε ειμαι σιγουρη γιατι δε τα επαληθευσα

ειναι λάθος :D :D :D

Προλαβαινω αμα αρχισω σημερα? Εχει πολλη δουλεια η εργασια?

νομιζω βγαινει

Υπάρχει περίπτωση για παράταση της εργασίας μερικές μέρες ; Έχετε ακούσει τίποτα ; είναι κανείς που θα ήθελε και αυτός να ζητηθεί μια παράταση ;

Παιδες μπορει μήπως κάνεις να ανεβάσεi το Noise.mat και sound.mat με link απο dropbox γιατι για κάποιο λόγο δεν μου τα κατεβάζει απο e-learning...

voila
https://www.dropbox.com/sh/fx40ttlxwic7q9q/AACrPjNMiDqG9FlVvUK2pvsSa?dl=0

Παιδιά τελικά η διορία είναι την Παρασκευή ή την Κυριακή?

στο ελεαρνινγκ λεει κυριακη 26

Λόγω της διαφοράς στην εκφώνηση της εργασίας που λέει Παρασκευή ρώτησα. Μπας και κάποιον τον είχε ρωτήσει προσωπικά. Αλλά ευχαριστώ!  :)

Lady_of_winter Cheers! Αν και τα βρήκα και από το περσινό Topic... δεν ξερω κατα ποσο βεβαια είναι ίδια

Ας τον ρωτήσει κάποιος στην τάξη και να μας πει αν είναι 24 η 26 στα σιγουρα

διαφορετικα ειναι, περσι εβγαινε αλλο κομματι

το πρώτο ερώτημα θα πρέπει ουσιαστικά να γίνει στο χέρι ;

 – Υπολογίστε:
* Τον πίνακα αυτοσυσχέτισης R του σήματος u(n), το διάνυσμα ετεροσυσχέτισης p του u(n)
και του d(n) και τους βέλτιστους συντελεστές του φίλτρου Wiener (wo), κάνοντας χρήση
της εξίσωσης Wiener-Hopf.
* Το ελάχιστο μέσο τετραγωνικό σφάλμα Jmin, για το φίλτρο Wiener του προηγούμενου
υποερωτήματος.
* Το πεδίο τιμών της παραμέτρου µ για το οποίο ο steepest descent συγκλίνει προς τη βέλτι-
στη λύση Wiener (wo).

γιατί;; είχαμε περίπου ίδια. Τι βρήκες τελικά;

0.3417, 0.0763, -0.0219

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

Καμία ιδέα για το τί κάνουμε με τον ήχο; Δεν μπορώ να κάνω την αυτοσυσχέτιση

τελικα 24 η 26 η προθεσμια?

Συμφωνώ! Για το denoising βγήκανε:
 r0 = 0.768826;
 r1 = 0.171613;
 r2 = -0.0493559;
?

Ναι, αλλά για τη λύση των γραμμικών εξισώσεων που θα προκύψουν, για να κάνεις εύκολη τη ζωή σου χρησιμοποίησε matlab ή wolframalpha...

Ρώτησα σήμερα τον Πιτσιάνη και είπε 26 είναι η προθεσμία.

Παίδες έτσι δεν υπολογίζουμε το διάνυσμα ετεροσυσχέτισης??? δηλαδή επείδη η Α ειναι Σ.Α. απο το ημιτονο με μέση τιμή 0 μηδενίζει όλον εκείνον τον ορο σωστα?

Πάντως στο σύστημα γράφει προθεσμία Κυριακή 26 Μαρτίου 12:00 π.μ. που πρακτικά είναι κάτι ενδιάμεσο. Απλά ενημερώνω, δεν έχω πάει στο μάθημα, προφανώς εσείς ξέρετε καλύτερα.
Εγώ, πάντως λέω να μην το ρισκάρω.

Παιδια κατάφερε κανείς να ακούσει το τραγούδι χωρίς ΚΑΘΟΛΟΥ χιόνια??

Για τον υπολογισμό του P, επειδή έχω μπουρδουκλωθεί σε κάτι που μάλλον είναι απλό, είμαι σε καλό δρόμο ή όχι?

Παιδιά στο quiz ξέρετε τι ακριβώς συμπληρώνουμε; Μόνο ένα τελικό νούμερο για την τιμή που ζητάει στο καθένα ή και τις εξισώσεις που ζητάει;

Μόνο το τελικό νούμερο νομίζω...

Πού υπάρχουν οι τύποι για τις λύσεις των ερωτήσεων του κουίζ ;

Τσέκαρε το pdf των ασκήσεων! Έχει ίδιες ή σχεδόν ίδιες ασκήσεις!

Το Ru που ζητάει στο ερώτημα με το τραγούδι το υπολογίζουμε χρησιμοποιώντας τον πίνακα του noise.mat? Και πώς ακριβώς; με τον τύπο ή υπάρχει συνάρτηση στο matlab?

Έχει καταφέρει κανείς σε matlab να κάνει διάγραμμα για το error surface and adaptation process για w=[w0 w1 w2], όπως κάνει ο Σισμάνης στο script του (gradWienerFilt.m) ?

ναι χρησιμοποεις τον πίνακα εκείνο και παίρνεις το u του

a = xcorr(u,u,nCoeff-1,'unbiased');
a = a(nCoeff:(2*nCoeff-1));
R = toeplitz(a);

nCoeff προφανως = 3;
και το p το πρωτο στοιχειο 0.72 (διασπορα του λευκου θορυβου) και όλα τα άλλα μηδενικα...

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

Παιδιά το σd ^ 2 πόσο σας βγαίνει και αν είναι εύκολο τον τρόπο που προκύπτει ??

Το διάγραμμα αυτό είχε πει στην τάξη να μην το κάνουμε

Α ευχαριστώ πολύ! Δεν ήξερα για την παράμετρο 'unbiased'... Είδα τώρα λίγο και τις υπόλοιπες παραμέτρους αλλά δεν κατάλαβα, γιατί επιλέγουμε την 'unbiased'? Με ποιο κριτήριο;

Παιδιά το quiz είναι μπονους ή απλά για εξάσκηση?

w=[-1;-1;-1];
mu=0.1;

wt=zeros([3,n]);
wt(:,1)=w;
y=zeros(n,1);

s=[0;0;u];
for i=2:n
  w=w+mu*(p-R*w);
  wt(:,i)=w;
  y(i)=s(i:-1:i-1)' * w;
end

παιδια τι κανω λαθος και δεν τρεχει?? :-\ :-\

Πήρε παράταση η εργασία?? Γιατί στο elearning εμένα μου λέει λήξη "Τετάρτη, 29 Μάρτιος 2017, 12:00 πμ"

Λογικά ΝΑΙ!  ^hello^  ^super^  ^beer^  :D

Είναι μπόνους και δίνει +0.125. Ενώ η εργασία +0.3. Τσέκαρε τις "πληροφορίες μαθήματος" στο elearning.

Νομίζω είναι αδύνατον

Στο quiz σας βγάζει αν είναι σωστές οι απαντήσεις σας;

όχι!

Αν ξεκινήσω αύριο έχω ελπίδες να παραδώσω έστω και κάτι αξιόλογο;;( Δεδομένου πάντα ότι δεν έχω πάει σε μάθημα :D )

Έχεις ελπίδες να παραδώσεις κάτι που θα πάρει 10/10 (αναλόγως πόσο θες να ασχοληθείς πάντα). Αύριο βγάζεις θεωρία-ασκήσεις και έχεις δυο μέρες να δουλέψεις εργασία και quiz (βγαίνουν άνετα σε τόσο).
Good luck! :)

Έχουμε ρωτήσει τον Σισμάνη, δεν θα παραδώσουμε τέτοιο διάγραμμα.
Είναι πολύ δυσκολότερο (αν όχι αδύνατον) να απεικονίσεις αυτή τη πληροφορία με παραπάνω από δύο συντελεστές Wiener.

Αυτές είναι οι τιμές για την ruu ή την ruv;

r_uu

Επειδή έχεις πλέον τρείς συντελεστές θα πρέπει η for σου να έχει αυτή τη μορφή:

for i=3:n
  w = w + mu*(p-R*w); % Adaptation steps
  wt(:,i) = w;
  y(i) = s(i:-1:i-2)' * w; % filter
end

ουπς σωστα ! ευχαριστω πολυυ   :D :D :D

Μπορεί κάποιος να με βοηθήσει με τον πίνακα αυτοσυσχέτισης του noise.mat; Οποιοσδήποτε συνδυασμός εντολών που περιλαμβάνει αυτοσυσχέτιση για τον τρέχοντα μεγάλο όγκο δεδομένων καταλήγει σε crash σε ένα pc που το restart είναι μια όμορφη περιπέτεια!
Για να τεστάρω τα xcross κλπ τα εκτέλεσα σε 200 steps σε white noise που δημιούργησα μόνος μου, αλλά κάθε φορά βγάζει διαφορετικά αποτελέσματα και καμία σχέση με τα θεωρητικά που -σωστά- έχω βγάλει.

Οποιαδήποτε ιδέα, παρατήρηση, σκέψη είναι ευπρόσδεκτη. Επί ΚΥΡΙΟΛΕΚΤΙΚΑ 5.5 ώρες ασχολούμαι με αυτό.

EDIT: Βρήκα σε παλαιότερο ποστ (2015) για την εργασία έναν αλγόριθμο να υπολογίζεις "χειροκίνητα" autocorrelation. Εφόσον δε λειτουργούσε καμία συνάρτηση κατέφυγα σε αυτό.

Το κομμάτι ακούγεται τέλεια. Ο κώδικας είναι ο εξής

Υλοποιήστε τη διάταξη του σχήματος σε MATLAB. Για την προσαρμογή των συντελεστών του
φίλτρου χρησιμοποιήστε τον αλγόριθμο steepest descent. Δοκιμάστε διαφορετικές τιμές της
παραμέτρου µ που να βρίσκονται εντός και εκτός του διαστήματος σύγκλισης.

Για το v(n) που είναι λευκός θόρυβος χρησιμοποιούμε την rand και όχι την randn ;
Και με πόσο πρέπει να την πολλαπλασιάσουμε για να έχει διακύμανση 0.32;

edit: Το βρήκα: v=sqrt(0.32)*2*sqrt(3)*rand(N, 1);

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

και γω αυτο εκανα..αλλα μετα sto e=d-y δεν εχουν ιδιες διαστασεις τα d και y...που το χανω?

Αν δεν έχουν ίδιες διαστάσεις απλά κάνε ένα transpose τον y (και σε μένα έτσι βγαίνει). Αν δεν έχουν και το ίδιο μέγεθος δεν ξέρω τι να σου πω, πρέπει να γράψεις κάτι παραπάνω για να βοηθήσουμε.

Παρεπιπτόντως πήρε παράταση η εργασία σύμφωνα με το σύστημα!

εγω εβαλα το n ισο με το μεγεθος του u και το ακουω οκ το τραγουδι? ειναι λαθος η λογικη?

Κατ'αρχάς εφόσον δουλεύει είναι σωστό! :D

Επίσης κι εγώ αυτό κάνω απλά χρησιμοποιώ τη μεταβλητή usize ως n.

νταξει οκ το ιδιο πραγμα κανουμε :D :D

Στο κουιζ στο πρωτο ερωτημα πως υπολογιζεται η διακυμανση? Υπαρχει καπου καποιος τυπος?

Τσέκαρε παλαιότερα ποστ.  ;)