THMMY.gr

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

Οτιδήποτε δεν έχει σχέση με τα παραπάνω θα μεταχειρίζεται αναλόγως

ΨΗΦΙΑΚΗ ΕΠΕΞΕΡΓΑΣΙΑ ΣΗΜΑΤΟΣ    
Παρασκευή, 07 Οκτώβριος 2011 00:00

Η διδασκαλία του μαθήματος ΨΗΦΙΑΚΗ ΕΠΕΞΕΡΓΑΣΙΑ ΣΗΜΑΤΟΣ (θεωρία) θα ξεκινήσει τη Δευτέρα 10-10 (13:00-15:00).

 

Οι διδάσκοντες

Στην 7? Γιατί πηγα 1.20 αι βρήκα κλειστή πόρτα :P

Επειδή ο κ.Χ#Λεοντιάδης δεν είχε τα κλειδιά της 7, μεταφερθήκαμε στην 5 ..δεν ξέρω αν θα είναι μόνιμη η αλλαγή όμως.

έκανε τίποτα σήμερα ή απλά είπε εισαγωγικά πράγματα???

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

Μπορεί να ενημερώσει κάποιος τι έκανε χθές θεωρία;

χτες είπε για ιδιότητες DFT, κυκλική μετατόπιση και κυκλική συνέλιξη

Σήμερα έκανε ο Πανάς; Γενικά τι ειπώθηκε σμρ από ύλη;

Thanks

ακουσα οτι δεν ηρθε ο πανας αλλα εκανε μαθημα ο λεοντιος. τωρα απο υλη... ρωτα κανα τηλεπ (μπρρρρ...  :P :P )

Ο Λεόντιος μας έκανε μάθημα σήμερα όντως και μιλήσαμε για την σύνδεση μετασχηματισμών Ζ με φίλτρα FIR και IIR. Α και έλυσε και ένα θέμα εξετάσεων!

Φετος εργασίες θα έχει?

Είχε πει ότι θα έδινε προαιρετικές που δεν θα δίνουν εξτρά βαθμό αλλά θα βοηθούν στην κατανόηση, αν θυμάμαι καλά..

Θα εχει εργασια που θα πιανει 2 μοναδες. Ετσι ειπε τουλαχιστον.

@wanderer: είχα στείλει μαιλ στον Χ#λεοντιάδη και τον είχα ρωτήσει το ίδιο πράγμα και μου είπε πως είχε σκοπό να βάλει νέα homework, αλλά ακόμα και αυτά να αφήσει, είναι εξίσου καλά για την δική μας εξάσκηση, δεδομένου πως αυτά δεν είναι για βαθμό και απλά μπορούμε να τα κάνουμε και να ρωτήσουμε απορίες!

Όντως.

(short response :P )

Tην Παρασκευη 23 δεν θα γινει μαθημα.

ποιο θέμα έλυσε? μπορείς να ανεβάσεις τη λύση?   ^beg^

βγηκε η εργασια του μαθηματος!

http://psyche.ee.auth.gr/index.php?option=com_content&view=category&layout=blog&id=5&Itemid=6&lang=el (http://psyche.ee.auth.gr/index.php?option=com_content&view=category&layout=blog&id=5&Itemid=6&lang=el)

κάποιος που να το έχει διαβάσει το μάθημα μπορεί να προτείνει απο που να διαβάσουμε και με πόσες μέρες διάβασμα μπορούμε να το περάσουμε?? :)

από εδώ: http://psyche.ee.auth.gr/index.php?option=com_content&view=category&layout=blog&id=5&Itemid=6&lang=el

Όσοι ενδιαφέρονται να εκπονήσουν την εργασία στα πλαίσια του μαθήματος σε ομάδα (μέχρι 3 άτομα) θα πρέπει να στείλουν τα ονόματα των μελών της κάθε ομάδας μέσω της φόρμας επικοινωνίας (ένα e-mail ανά ομάδα) το αργότερο μέχρι την Παρασκευή 21/01.

ΨΗΦΙΑΚΗ ΕΠΕΞΕΡΓΑΣΙΑ ΣΗΜΑΤΟΣ    
Τετάρτη, 18 Ιανουάριος 2012 00:00

Το μάθημα των ασκήσεων της Παρασκευής 20/01 δε θα γίνει.

φετινές σημειώσεις κανείς;

τελικα στην ύλη των φιλτρων τι ειναι μεσα και τι εκτος?

Σχεδίαση φίλτρων ΙΙR με διγραμμικό μετασχηματισμό, ναι. Όλα τα άλλα, πάνω κάτω όχι.

άκουσα ότι βγήκαν εκτός τα chebyshev και κάτι άλλα... μήπως θα μπορούσες φίλτατη εσύ ή κάποιος άλλος να μας πείτε με σελιδες ή ενότητες τι είναι μέσα από το κεφάλαιο 6?  :)

Εγώ πιστεύω-σύμφωνα με όσα ειπώθηκαν στην θεωρία- πως μέσα είναι:
6.1, 6.2 , 6.3 , 6.4 ,  6.6 , 6.8 , 6.13

παντού ασχοληθήκαμε μόνο με Butterworth, σε Chebyshev + ελλειπτικά δεν δόθηκε σημασία γενικά. Απλή αναφορά τους έγινε μόνο..

merci beaucoup  ^miss^

επειδη κατι νομιζω ειχε ειπωθει παλιοτερα, δινουμε με ανοιχτο τον Hayes?

ναι
και με τον Πανα

ευχαριστω (ελπιζω ο Hayes να μου ειναι αρκετος :) )

το 6.5 είναι μέσα??? επίσης από το 6.4 είναι το Α μόνο? γενικα πάντως εγώ την ύλη για το 6ο την ξέρω ότι πάντα ήταν

6.1- 6.2-6.3-6.4-6.6-6.8-6.9 κ 6.13

τουλάχιστον όταν ήταν η σειρά μου! φέτος νάσια είσαι σίγουρη ότι είναι αυτή η ύλη?

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

Φέτος δεν κάναμε ούτε chebychev πολύ πολύ-τα αναφέραμε απλά- ούτε και ελλειπτικά.. ενιγουει.


edit: Ανακάλεσα για το 6.13 και τελικά και εγώ πιστεύω πως το 6.5 δεν είναι μες στην ύλη.
(από το 6.9 έχουμε κάνει μόνο μια άσκηση, δεν θυμάμαι καν να τα είπαμε στην θεωρία, anyway..)

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

είναι κανονικα στις 29 έτσι?

Ναι κανονικα για τις 29.

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

έλα ηλία, να σου πω ρε συ... μίλησες όντως εσύ με τον λεόντιο ή κάποιος άλλος?

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

κρίμα είναι!

 

καλε απαντηστε κάποιος μην ντρέπεστε

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

εγώ ο ίδιος μίλησα
ναι, φαίνεται ελαστική η προθεσμία

ναι εστειλα ! περιμένω απαντηση αν απαντήση θα ενημερώσω!

ΘΑ ΔΟΘΕΙ ΠΑΡΑΤΑΣΗ ΜΕΧΡΙ 7/2

μόλις μου απάντησε η διδακτορικός!!!

Από τις Ανακοινώσεις στο psyche:

"Η εργασία στην Ψηφιακή Επεξεργασία Σήματος 2011/12 παίρνει παράταση ως τις 07/02"

τα αποτελεσματα δηλαδη για τα γραπτα θα βγουν μετα?μαζι με τις εργασιες??

 ναι

παιδιά, για πείτε γιατί μάλλον δεν έχω καταλάβει το νόημα
στο ΙΙ της εργασίας

φιλτράρουμε το y (σήμα με ηχώ) με το adaptive filter που φτιάξαμε και τι προκύπτει;
η ηχώ;
το y χωρίς ηχώ;

γιατί το σήμα-αποτέλεσμα μου ακούγεται σα θόρυβος;
εσάς σας βγαίνει νορμάλ;

φιλτράρεις το y και παίρνεις το σήμα χωρίς ηχώ.

Το κακό είναι ότι η καθυστέρηση των 200 δειγμάτων αντιστοιχεί σε χρονική καθυστέρηση 200/8000 =0.025sec δηλ το σήμα με ηχώ σχεδόν δε διακρίνεται (ακουστικά) από αυτό χωρίς ηχώ.

Βέβαια, αυτό δεν είναι προβλημα στην πε΄ριπτωση του σήματος εκπαίδευσης, γιατί απλά με ένα mean square error μπορείς να ξέρεις την απόδοση του φίλτρου.

στο πρώτο μέρος ο ήχος που κόψατε σας βγήκε ακέραιος??? δλδ ακούγεται λίγο και ο άλλος? πολύ λίγο ομως

ναι θα ακουγεται και ο αλλος λιγο δεν μπορεις να κανεις κατι.. καθε ατμομηχανη περιέχει συχνότητες χαμηλού πλάτους και στην περιοχή συχνοτήτων της αλλης.

ναι ακριβως! απλα λεω μήπως! αλλα είναι αναποφεκτο

ρε παιδια χρειαζομαι λιγο βοηθεια στην εργασια

Τον ηχο πως θα τον ξεχωρισω σε 2 σηματα ?? εψαξα αλλα δεν βρηκα καποιον τροπο να κανω την αντιστροφη διαδικασια της συνελιξης ωστε να βρω τα 2 σηματα. Και επιπλέον δεν ξερω τη χρονικη διαρκεια ξεχωριστα των δυο κυματομορφων.

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

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

μπορείς να ψάξεις τον δεύτερο ήχο αν θέλεις, αλλά είναι πιο εύκολο με τον πρώτο ήχο





γενικά στο δεύτερο μέρος, πώς σας βγαίνει ο ήχος στο (γ), εμένα χάλια, σα μουρμούρα :P

αυτό σημαίνει ότι το φιλτραρισμένο σήμα πρέπει να ακούγεται σαν το ιδεατό, ε;
ακόμα
τι mse, αφού στο φίλτρο που υλοποιείς στο (α) και χρησιμοποιείς στο (γ) οι παράμετροι είναι στάνταρ για τον LMS
επίσης
το μοντέλο της ηχούς λέει ότι έχουμε δυο σήματα ηχούς, ένα με καθυστέρηση nd και ένα με 2*nd
ισχύει αυτό, ή θεωρητικά το έγραψε στην εκφώνηση;

qube thanx  και κατι ακόμα.....
Όταν ζηταει τα φασματικα χαρακτηριστικα ζητάει τη φασματικη πυκνότητα ισχύος η το εύρος ζώνης??
(απτο διάγραμμα του fft πάντως δεν βγαινει και μεγάλη άκρη)

δες στο χελπ του ματλαμπ πως υλοποιει το fft.. πρεπει να παρεις την απολυτη τιμη και να ορισεις σωστα τον οριζοντιο αξονα..

Μερικές απορίες για το 2ο μερος.

Πρωτον λεει οτι η Ηec θα έχει συνολικά 600 δείγματα; Κάνοντας τη συνέλιξη για να πάρω το x=conv(hec,y) προκύπτει ένας πινακας με 19043x1. Για να πάρει το σφάλμα παίρνει τη διαφορά αυτού που υπολογίσαμε με το x[n]. Το θέμα ειναι οτι οι πινακες έχουν διαφορετικές διαστάσεις, άρα πώς θα πάρουμε τη διαφορα; Επίσης εκει που παίρνει τη δεύτερη παραγωγο, ουσιαστικα θα δημιουργηθεί δισδιάστατος πινακας ως προς n και j?
Thnx

μη χρησιμοποιήσεις το conv, κάντο χειροκίνητα
θα πάρεις 600 δείγματα από το y (που μετακινούνται αρνητικά, όπως θέλει η συνέλιξη) και θα πολλαπλασιάσεις το hec με τα επιλεγμένα y για να πάρεις μια τιμή του xhat (θα βγει scalar)
και θα κάνεις x-xhat, μην κάνεις το xhat πίνακα και μπλεχτείς

τι δεύτερη παράγωγο...?

να ρωτησω λιγο... στο δευτερο μέρος πρώτο ερώτημα για να υπολογίσω την h[n] χρησιμοποιώ τα x και y που μου δίνει ή μόνο τις παραμέτρους Μ και Δ?

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

Το mse που λέω είναι το mse(x_train_orig - x_hat) όπου x_hat = conv(hEC, y_train_echoed) και h_EC είναι η τελική εκτίμηση που προκύπτει από τον LMS. [Βέβαια υπάρχει και το x_hat που προκύπτει στη διάρκεια εκτέλεσης του LMS]

Όσον αφορά την καθυστέρηση των 0.025sec, αν δοκιμάσει κανείς sound(x_train_orig,8000) και sound(y_train_echo,8000) δε θα καταλάβει διαφορά εκτός αν έχει το supreme αυτί  :D :D

Κάτι που είπες παραπάνω να μη χρησιμοποιηθεί η conv αλλά να γίνει χειροκίνητα η συνελιξη, γιαυτό είσαι σίγουρος;

εννοείς πχ y = [y(1000), y(999), ..., y(401)] σωστά;;

τι να σου πω, άμα την παλεύεις με όλα τα μηδενικά που βγάζει το conv στο x_hat (πίνακας πλέον), είσαι ωραίος

εμένα με μπέρδεψε



στο δεύτερο ποστ σου, ναι, αυτό λέω

Που βγάζει μηδενικά;; Κανονικά είτε γίνει χειροκίνητα είτε όχι πρέπει να δίνει το ίδιο αποτέλεσμα... Μήπως αναφέρεσαι στη συνέλιξη που γίνεται κατά τη διάρκεια εκτέλεσης του LMS (στην εκπαίδευση δηλαδή);;;;

[Αναφερόμουν στη συνέλιξη στο (γ) ερώτημα, στην υλοποίηση του LMS όντως δεν ενδείκνυται να γίνει με conv]

η παράμετρος Μ είναι ουσιαστικά ο αριθμός των συντελεστών του πολυωνύμου που είναι το φίλτρο FIR, Δ είναι το βήμα του αλγόριθμου που κατασκευάζει το adaptive filter

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

εδώ συντελεστές του συστήματος είναι οι συντελεστές του φίλτρου FIR που είναι Μ σε αριθμό




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

και ναι, στη συνέλιξη μέσα στον LMS αναφέρομαι, υπάρχει λόγος να κάνεις συνέλιξη έξω;


edit: αχα

Στο (γ) ερώτημα θέλει συνέλιξη (αναγκαστικά, δε γίνεται αλλιώς)

ναι οκ, εγώ το κανα με την filter οπότε υπονοείται συνέλιξη

Για το γ ερώτημα (2ο μερος) εγώ έκανα το εξής :

>>load x_echo.mat
>> N3 = length(x_echo)

N3 =

       37087

>> for n = 600:N3
x_orig2(n)=hce*x_echo(n:-1:n-600+1);
end

Επειδή στο σφάλμα προσεγγίζαμε την x_echo_trainset_orig με αντίστοιχο γινόμενο, σκέφτηκα πως για να προσεγγισω την x_orig2 με Input signal την x_echo και τους ίδιους συντελεστές θα κάνω το ίδιο...

Είναι εντελώς λάθος ??? Πάντως οι κυματομορφές που παίρνω δεν ειναι κακές...

στις εξετάσεις πηγαίνουμε με σημειώσεις, με βιβλίο ???

επίσης η εργασία είναι προσθετική? κι αν ναι πόσο?

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

ευχαριστώ! hayes ή του πανά το βιβλίο.. ή δν εχει σημασία?  :P

Και τα 2 μπορείς να έχεις ανοιχτά!  ;)

ακόμη καλύτερα! ευχαριστώ πολυ

Γίνεται κάποιος που έχει το βιβλίο του Πανά να μου το δανείσει μέχρι την εξέταση?
Αν το χρειαζεστε απλά να φωτοτυπήσω το κεφάλαιο που χρειάζομαι..
PLZ...

Για πειτε... τι διαβαζουμε για το μαθημα?

ανεβασα τα σημερινα θεματα και τις λυσεις μου.
https://www.thmmy.gr/smf/index.php?action=tpmod;dl=item1985

πειτε και το 1β πως στο διαολο λυνεται :D

Στο 2ο θέμα έχει ίδιο πόλο με μηδενικό, αυτά δεν απλοποιούνται κανονικά?

ξερω γω; (http://shmmy.ntua.gr/forum/images/smilies/smilies/MuchCry.gif)

Τι ανεβαζεις λυσεις ρε κεφτεδα αμα δεν εισαι σιγουρος. Θες να καψεις ατομα?

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

χοχοχο, παει κι αυτο το λινκ.

τουλαχιστον πειτε μου αν περναει το γραπτο :D :D :D :D :D :D

Για το θέμα 1Β, για άκου μια πρόταση :

Χρησιμοποιώντας τον τύπο που σου υπολογίζει το h[n] από το Η(e^jω), εκείνον με το ολοκλήρωμα, βγαίνει ένα πράμα :

1⁄2π ∫e^(jω(n-60)) e^(-j25ω^3 ) dω

Και τέρμα χοντρικά, αν πούμε ότι όλα τα e^j(x) όταν ολοκληρώνονται από -π ως π πάντα θα βγάζουν κάτι μικρότερο της μονάδας, τότε θα έχουμε μάξ πλάτος 1/2π. Ακόμα αν πούμε ότι ο 1ος όρος του ολοκληρώματος είναι η βασική συχνότητα του σήματος και ο 2ος όρος είναι αυτός που προκαλεί την "απόσβεση" ως προς το n, μπορούμε να πούμε ότι εμφανίζεται μια μετακίνηση κατά 60 δείγματα.

Με αυτά τα δεδομένα, δικαιολογώ την απάντηση 3

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


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

Πέστε απόψεις.

Στο 3ο θέμα εγώ πιστεύω ότι πρέπει να βρεθεί με συνελίξεις. Δηλαδή ότι ισχύει: y[n]=x[n]*h[n] = Σx(k)h(n-k). Οπότε y(0)=Σx(k)h(-k)=x(0)h(0) (όλοι οι άλλοι όροι του αθροίσματος μηδέν προφανώς). Από εξίσωση διαφοράς έχουμε: y(0) = -b0x(0). Αρα b0 = -h(0) = -35. Μετά πάλι απο συνέλιξη: y(1) = Σx(k)h(1-k) = x(0)h(1) + x(1)h(0) = 13x(0) + 35x(1). Από εξίσωση διαφοράς y(1) = -a1y(0) -b0x(1) -b1x(0), το y(0) = 35x(0), το βάζω στη προηγούμενη, εξισώνω με το y(1) που έχω από την συνέλιξη κλπ κλπ. Άρα βρίσκω τα α0, α1, α2, β0, β1.

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

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

Στο τρίτο θέμα ήταν σύστημα 4 εξισώσεων με 4 αγνώστους, άρα λύνεται κατευθείαν  ;D

μπουρδελο ηταν οι εξισωσεις ρε, να κανει τις πραξεις μονος του και να μας πει κι εμας τις λυσεις.
αν δεν εχει μπει το 2013 εμενα να μου τρυπησεις τη μυτη

Αυτό που λέει η αιμιλία, ήταν μπουρδέλο οι πράξεις δεν βγαίνουν. Ενώ με συνελίξεις (ή δυναμοσειρες παρομοια) είναι πιο ανθρώπινο.

Το Ti 89  βρήκε τη λύση μια χαρά πάντως  ::)