THMMY.gr

Παρέμβαση συντονιστή

Αλλάχτηκε ο τίτλος έτσι ώστε οι συζητήσεις γύρω από την εργασία να συγκεντρωθούν σε αυτό το τόπικ όπως συνηθίζεται και βολεύει. Καλή δουλειά σε όλους!!!

----------------------------------------------------------------------------------------------------

ξερει καποιος αν εχει ξεκινησει η εργασια; υπαρχουν καπου τα θεματα;

 Εργασία 2008-9

19 Ιαν 2009 9:29 μμ
Ντελόπουλος

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

Μιας και έχουμε παράταση μέχρι τις 15 Μαρτίου...είπαμε να ασχοληθούμε με την εργασία. :)

Έχουμε κολλήσει λιγάκι στα διαγράμματα, μας βγαίνουν λίγο περίεργα. Χρησιμοποιούμε τη freqz, αλλά δεν μας αποκρίνεται καλά. :o :(

καμιά άλλη ιδέα?? 8)

Τι περιεργο εχουν δηλαδη;
βάλ' τα λιγο εδώ αν μπορείς.

Λέτε για τα διαγράμματα στην πρώτη ενότητα; (Βαρέθηκα να διαβάσω τις υπόλοιπες. :D)

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

Ναι, για τα διαγράμματα στην πρώτη ενότητα λέγαμε. Κάτι βρήκαμε, αλλά ακόμη το παλεύουμε.

Και οι απορίες συνεχίζονται: Στην πρώτη ενότητα και πάλι, στο 4.(α) θέλουμε σε κάθε επανάληψη, να διαβάζουμε 1152 δείγματα του αρχικού ηχητικού σήματος, σωστά?? Θα κάνουμε δηλαδή 514944/1152 επαναλήψεις για να καλύψουμε όλα τα δείγματα?

Μετά πρέπει να δημιουργήσουμε το xbuff, ώστε να έχει 1632 δείγματα? ΠΩΣΣΣ?? :o :-[

Για τα διαγράμματα νομίζω ότι πρέπει να κάνεις FFT σειρά-σειρά και στο plot για τον άξονα των X να πολλαπλασιάσεις τον αύξοντα αριθμό με τη συχνότητα δειγματοληψίας.

Στην άλλη σου ερώτηση πρέπει να γίνεις λίγο πιο συγκεκριμένη...

Για τη δεύτερη ερώτηση, αυτό που καταλάβαμε (διαβάζοντας τον κώδικα της συνάρτησης Y=frame_sub_analysis(xbuff,H,q) ) είναι ότι το xbuff θα πρέπει να έχει πλήθος στοιχείων μεγαλύτερο ή ίσο του 1632 γιατί στο κώδικα εκτελεί την εντολή: X = xbuff(ind) και υπολογίζοντας τον πίνακα-δείκτη ind είδαμε ότι το μεγαλύτερο στοιχείο του έχει τιμή 1632. Σωστά μέχρι εδώ?

Χρησιμοποιήσαμε την εντολή buffer για να δημιουργήσουμε στήλες των 1632 δειγμάτων του αρχικού σήματος και μετά κάναμε κλήση της Y=frame_sub_analysis(xbuff,H,q)  όπου xbuff μια στήλη δειγμάτων(από τον νέο πίνακα που δημιουργήσαμε), Η όπως δίνεται  και q=N...
Αλλά εμφανίζει πρόβλημα στο σημείο X = xbuff(ind) ...πιο συγκεκριμμένα: In an assignment  A(I) = B, the number of elements in B and I must be the same. 

Τι καταλαβαίνουμε λάθος? Και γιατί έχουμε αυτό το warning?  ::) :???:

Ρε παιδιά δεν ασχολείται κανείς με αυτή την εργασία;Αν είναι να την παρατήσουμε και εμείς...  :'(
Καμιά συμβουλή;;;;;;;;;;;;;;;;; :-\

απο ετημμυ:

Συνάντηση για την εργασία.

9 Μαρ 2009 3:48 μμ
Ντελόπουλος

Την Πέμπτη 12.2.2009 και ώρα 11:00 θα είμαι διαθέσιμος για απορίες σχετικά με την εργασία που είναι σε εξέλιξη. Κάπου στις αίθουσες...

Το είδαμε αυτό! αλλά το αύριο μοιάζει λίγο αργά... ::) και δεν μπορούμε να προχωρήσουμε!! (βλέπε απορία από Anastasia_ πιο πάνω)καμιά άλλη ιδέα;;;;;

Βεβαιωθείτε ότι σε κάθε επανάληψη παράγονται Μ*Ν δείγματα.

Δηλαδή το xbuff(που δίνουμε σαν όρισμα στη frame_sub_analysis) σε κάθε επανάληψη να είναι ένα διάνυσμα μήκους 1152???

Σύμφωνα όμως με τις επισημάνσεις το buffer εισόδου πρέπει να είναι διάστασης (N-1)M+L .Αυτό δεν βγαίνει 1632?Αυτό πού κολλάει?

Όχι ο buffer, τα δείγματα που υπολογίζετε κάθε φορά με την frame_sub_synthesis πρέπει να είναι M*N. Δηλαδή η εντολή πρέπει να είναι xhat(i:i+M*N-1)=frame_sub_synthesis(μπλαμπλα) όπου i βάζετε ό,τι σας βολεύει ανάλογα με τη δομή της codec0.

Πόσες επαναλήψεις πρέπει να γίνουν στο ερώτημα 1.4 και πως προκύπτει?

Το xbuff πως το υπολογίζουμε?
Η σχέση (Ν-1)Μ + L =1632  πως προκύπτει αφού αρχικά λέει ότι σε κάθε επανάληψη παίρνουμε ΝxΜ=1152 δείγματα?
Κάποιος που έχει κάνει το πρώτο ερώτημα ...ας βοηθήσει plzzzzzzz !

Ο αριθμός των επαναλήψεων προκύπτει διαιρώντας το συνολικό πλήθος των δειγμάτων με M*N.

Το μήκος του buffer είναι μεγαλύτερο από τα δείγματα που θα πάρεις, γιατί κρατάει και τα L προηγούμενα, που θα χρησιμοποιηθούν στον υπολογισμό.

Στις συναρτήσεις RLE και IRLE το Κ τι χρειάζεται; Δεν είναι το μήκος του symb_index ή τα έχω καταλάβει τελείως λάθος;

Ναι κι εγώ αυτό κατάλαβα και μάλλον δεν χρειάζεται, αλλά νταξ, ας το βάλουμε.. :D

Περαιτέρω απορίες:

Στον iframeDCT δεν πρέπει να είναι είσοδος και το Μ ή το Ν; Αλλιώς πώς θα ξέρουμε σε πόσες στήλες θα το χωρίσουμε;

Στην κωδικοποίηση Huffman τις πιθανότητες πώς θα τις βγάλουμε; Άκυρο, ξαναδιάβασα την εκφώνηση :P

Μπορούμε να χρησιμοποιήσουμε έτοιμες συναρτήσεις για το Huffman;

Βασικά οι έτοιμες συναρτήσεις δέχονται ως σύμβολα μόνο νούμερα και strings, ενώ εμείς έχουμε δυάδες αριθμών. Το έχει λύσει κανένας αυτό;

Ναι, αλλά προφανώς τα βάζουμε μόνοι μας.


Ναι.


Αντιστοιχείς τα καινούρια σύμβολα στις δυάδες που έχεις...

Απ' το 3.3 το ερώτημα 2ο...έχοντας μόνο το Κμαξ πώς χωρίζουμε σ'αυτά τα διαστήματα τους συντελεστές???? Βασικά νομίζω ότι δεν πολυκατάλαβα αυτή τη σχέση που μας δίνει για τα Δκ!! Κάποιος λίγο να με διαφωτίσει!!!

Φτιάχνεις έναν πίνακα (Kmax)x(Kmax). Οι γραμμές αντιχτοιχούν στους συντελεστές k και οι στήλες στις συχνοττητες. Άρα π.χ. για k=500, στη γραμμή 500 του πίνακα, θα έχεις μονάδες στις στήλες 2 μέχρι 13.

Ποια είναι τα καινούρια σύμβολα; :???:

Αυτά που προκύπτουν απο το huffmandict..

Και πώς θα τρέξω το Huffmandict αφού δεν έχω αρχικά σύμβολα;

Ε τα δημιουργείς εσύ. Η πρώτη δυάδα θα είναι το σύμβολο 1, η 2η το 2, κτλ.
Δες λίγο τα παραδείγματα στο help, βοηθούν πολύ..

Ευχαριστώ... αν και βασικά το έχω ψιλοπαρατήσει το Huffman, παιδεύομαι με κάτι άλλο τώρα...

Το Tq έχει κάποιες τιμές στο άπειρο, οπότε... θέλουμε άπειρο SNR; Πώς στην ευχή θα το κάνουμε αυτό; ^seestars^

Έλα ντε, κι εγώ αυτό αναρωτιέμαι.. Έτσι αυξάνει συνέχεια τα bits... :(

Btw, στο μάθημα που έγινε την Πέμπτη μας είπε ότι αυτό το SNR είναι λάθος, πρεπει να το πάρουμε: 10log(eb(i)^2).
Επίσης, μπορούμε να την παραδώσουμε μέχρι και τις 25/3 περίπου. Είπε θα βγάλει ανακοίνωση, αλλά ακόμα τίποτα.

Ποιο μάθημα; Αφού άρχισε το θερινό :???:

Α :D

Κάτσε λίγο όμως, αυτή η παράσταση για το SNR αυξάνεται με την ισχύ του σφάλματος αντί να μειώνεται... είσαι σίγουρη ότι αυτό σας είπε;

Ναι γιατί μας είπε στην ουσία είναι η ισχύς του σφάλματος όχι το SNR.

Τότε αν είναι έτσι πρέπει να είναι μικρότερο από το Tg, άρα δεν υπάρχει πρόβλημα με το άπειρο! ^super^

Έχω μπερδευτεί λιγάκι με το DCT (ναι, αυτό με τις 0 μονάδες :P). Κάνουμε DCT κάθε Ν-άδα χωριστά;

Ναι κάνεις DCT σε κάθε στήλη και μετά reshape.

Ξανά μπέρδεμα με το Tg... η ισχύς του λάθους βγαίνει ΠΑΝΤΑ μικρότερη με αποτέλεσμα όλα τα B να ισούνται με 1! :-[

Ναι, αυτό ακόμα δεν το έλυσα ούτε εγώ. Είχε πει ο Ντελόπουλος ότι θα βγάλει διευκρινίσεις στο ethmmy... ^dontknow^

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

Βασικά μας είπε ότι ούτως ή άλλως δεν πρόκειται να τις διορθώσει πριν τις 25/3.. :D
Μπορείς να του στείλεις ένα μαιλ να τον ρωτήσεις για πιο σίγουρα.. ;)

Αν την ανεβάσω σήμερα, και μετά τη βελτιώσω, μπορώ μετά να αντικαταστήσω το αρχείο;

Ναι είναι αλήθεια η προθεσμία για την παράδοση τη εργασίας μετατίθεται στις 29.3.09 (29 όχι 25) και επίσης ετοιμάζω κάποιες διευκρινίσεις-διορθώσεις.
Όμως το eTHMMY με δυσκολεύει...

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

Σχήμα 2 σελίδα 6 έχει τις critical bands (λείπει το νούμερο 10, αλλά φαίνεται από τα άλλα από πού μέχρι πού είναι).

Πριν λίγο βγήκαν διευκρινίσεις.

Όποιος κατάφερε να μετασχηματίσει το Tq με τρόπο που να βγάζει λογικά αποτελέσματα, ας δώσει ένα hint...

Λοιπόν γυρίζοντας προς τα πίσω έχω κάνει λάθος τη συνάρτηση Dksparse και την STinit!!!Έβαλα τον πίνακα D να είναι: από την 3 μέχρι τη 281 γραμμή να έχει άσσο στη δεύτερη στήλη, από τη 282 γραμμή μέχρι την 569 να έχει άσσους στις στήλες 2-13 και από εκεί μέχρι τη 1151 γραμμή άσσους από το 2 μέχρι το 27.Σωστά μέχρι εδώ???
Μετά όμως στην STinit, όταν πρέπει να ελέγξω τη γειτονιά, πρέπει να βγω εκτός ορίων. Όταν  είμαι κοντά στο 1152 πρέπει να παίρνω μόνο αυτά που είναι P(i-Dk) κι όχι τα P(i + Dk).Μπερδεύτηκα :o :o


Δοκίμασε αντί για το Tq το (Tq - max(Tq))!

Ναι λογικά, κοντά στα όρια, θα παίρνεις μόνο τις ισχείς που έχουν νόημα... ;)

Θα κάνω και γω μια πιο απλή ερωτησούλα γιατί από ότι βλέπω έχετε προχωρήσει αρκετά...
στο 1ο ερώτημα το snr το υπολογίζουμε σύμφωνα με τον κλασικό τύπο 10*log10(var(x)/var(noise))
όπου noise=x-xhat;
γιατί έτσι μου βγαίνει αρνητικό...ουψ!
υπάρχει κάτι που πρέπει να προσέξω; και πόσο περίπου πρέπει να βγαίνει;

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

Σας βγαίνει αρνητικό γιατί δεν παίρνετε υπόψη την καθυστέρηση που εισάγει η χρήση των buffers. Το λέει και στην εκφώνηση. Το SNR πρέπει να βγει πολύ μεγάλο στο 1ο ερώτημα αφού δεν έχει γίνει συμπίεση.

Έτσι βγαίνει άπειρο όμως.

και εγώ αυτό περίμενα...να βγει μεγάλο και παραξενεύτηκα με το αρνητικό.
το θέμα με την καθυστέρηση των buffers  πως μπορούμε να το χειριστούμε;
με κάποιο pause για παράδειγμα ή breakpoint;

Να σου πω κι εγώ δεν κατάλαβα γιατί έχει όση καθυστέρηση έχει, αλλά παρατήρησα ότι ισούται με το πλήθος των επαναλήψεων συν το μήκος Μ. Βαρέθηκα να κάτσω να δω τον κώδικα για να εξηγήσω γιατί :P

Νομίζω ότι υπάρχει γενικό πρόβλημα στον υπολογισμό του Tg.. Δεν αρκούν μετατοπίσεις για ικανοποιητικό αποτέλεσμα.. Σε εμένα τουλάχιστον..

Εγώ δοκίμασα να εφαρμόσω και απευθείας τον τύπο που έχει εδώ http://www.diracdelta.co.uk/science/source/t/h/threshold%20of%20hearing/source.html αλλά πάλι δε βγαίνει σωστά :(

Για να μη γράφουν στο τόπικ τόσον καιρό, μάλλον κανένας δεν το κατάφερε... :-\

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

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

Παίζει κανονικά;

Έβαλα απόλυτη τιμή ισχύος και δεν παίζει τίποτα :(

Τη συνάρτηση DCT_band_scale την κάνατε να βγάζει τους κανονικοποιημένους συντελεστές για κάθε μία μπάντα ξεχωριστά ή να δίνει όλους τους συντελεστές μαζεμένους και κανονικοποιημένους σε έναν 1x1152 πίνακα????Εγώ έκανα το δεύτερο αλλά κάπου μου δημιουργεί πρόβλημα!!!

Το πρόβλημα με το ότι ο αριθμός των bits βγαίνει πάντα 1, το έλυσε κανείς??? ^tafos^

Το πρώτο πρέπει να κάνεις, και θα βγάλεις και έναν συντελεστή sc για κάθε μπάντα.


Από ό,τι φαίνεται όχι... :-\

Κάτι χαζό...στον τύπο για την αποκανονικοποίηση των συντελεστών DCT, το σύμβολο μεταξύ του c και του Sc είναι απλό γινόμενο,έτσι δεν είναι????

Ναι, φυσικά.

Στην πρώτη άσκηση, τα διαγράμματα πώς περίπου πρέπει να μοιάζουν;

με μικρή καθυστέρηση, ανεβάσαμε τη δική μας εργασία.

Ναι, κι εγώ μπήκα, γι' αυτό το έσβησα τελικά. :)