diff --git a/tex/main.tex b/tex/main.tex index 2873998..9e473d2 100644 --- a/tex/main.tex +++ b/tex/main.tex @@ -52,7 +52,7 @@ Στο πρώτο επίπεδο υλοποιήσαμε τις βαθμίδες Sequence Segmentation Control και Filterbank. Επίσης, υλοποιήθηκαν και οι συναρτήσεις AACoder1, η οποία εκτελεί τις προηγούμενες ρουτίνες και αποθηκεύει τα αποτελέσματα σε ένα struct, και η αντίστροφή της iAACoder. - \subsection{\textit{function} frameType = SSC(frameT, nextFrameT, prevFrameType)} + \subsection{\texttt{function frameType = SSC(frameT, nextFrameT, prevFrameType)}} Η πρώτη συνάρτηση που υλοποιήθηκε αφορά την βαθμίδα Sequence Segmentation Control, και σκοπός της είναι να καθορίσει τον τύπο ενός frame. Το πρώτο όρισμά της είναι ένας πίνακας 2048x2 , που περιέχει 2 κανάλια ήχου στο χρόνο, του frame του οποίου τον τύπο θέλουμε να διαπιστώσουμε. Αντίστοιχα, το δεύτερο όρισμα περιέχει τα κανάλια του του αμέσως επόμενου frame ενώ το τρίτο, αφορά τον τύπο του αμέσως προηγούμενου frame, ο οποίος μπορεί να πάρει 1 από τις αλφαριθμητικές τιμές ``OLS'' (ONLY\_LONG\_SEQUENCE), ``LSS'' (LONG\_START\_SEQUENCE), ``ESH'' (EIGHT\_SHORT\_SEQUENCE) και ``LPS'' (LONG\_STOP\_SEQUENCE) , ανάλογα το περιεχόμενο του. @@ -74,7 +74,7 @@ Tα frame τα οποία είναι στατικά χαρακτηρίζοντα % Ανάλογα με την ισχύ ή όχι των παραπάνω συνθηκών και δεδομένου του τύπου του προηγούμενου frame, είναι δυνατή η απόφαση του τύπου του frame i, με τη λογική που παρουσιάζεται στο πρότυπο και την εκφώνηση της εργασίας. Για αυτό τον λόγο, στη συνέχεια, συγκρίνονται οι τύποι του επόμενου και του προηγούμενου frame και λαμβάνεται η απόφαση για τον τύπου του κάθε καναλιού. Για την απόφαση του τύπου ολόκληρου του frame, είναι απαραίτητος ο συνδυασμός των αποφάσεων των τύπων των δύο καναλιών σύμφωνα και πάλι με δοθείσα λογική. - \subsection{\textit{function} frameF = filterbank(frameT, frameType, winType)} + \subsection{\texttt{function frameF = filterbank(frameT, frameType, winType)}} Η συνάρτηση αυτή αφορά την βαθμίδα filterbank και έχει ως ορίσματα ένα frame ήχου 2048x2, τον τύπο του καθώς και τον τύπο παραθύρου που θα χρησιμοποιηθεί. Η δε έξοδός του, είναι η αναπαράσταση του ίδιου του frame στο πεδίο της συχνότητας, σε όρους MDCT. Στην περίπτωση που το frame είναι τύπου ``ESH'', η συνάρτηση θα επιστρέψει 8 υποπίνακες 128x2 με τους συντελεστές του MDCT, έναν για κάθε subframe ενώ σε άλλη περίπτωση, θα επιστρέψει έναν πίνακα 1024x2. @@ -82,13 +82,13 @@ Tα frame τα οποία είναι στατικά χαρακτηρίζοντα Στην περίπτωση που ο τύπος του frame είναι ``ESH'', πρώτα το χωρίζουμε σε 8 subframes και έπειτα εφαρμόζουμε την παραθυροποίηση και τον MDCT. -\subsection{\textit{function} frameT = iFilterbank(frameF, frameType, winType)} +\subsection{\texttt{function frameT = iFilterbank(frameF, frameType, winType)}} Η συνάρτηση της iFilterbank, αποτελεί την αντίστροφη της filterbank. Δέχεται, δηλαδή, σαν όρισμα τους συντελεστές MDCT ενός frame, τον τύπο του και τον τύπο του παραθύρου που θα χρησιμοποιηθεί ενώ ως έξοδο, επιστρέφει έναν πίνακα 2048x2 που αντιστοιχεί στο frame στο πεδίο του χρόνου. Αντίστοιχα με την filterbank, πρώτα ορίζονται τα παράθυρα ``KBD'' και ``SIN''. Πριν εφαρμοστούν, όμως, κατάλληλα στο frame, εφαρμόζουμε την συνάρτηση imdct4 για να επιστρέψουμε στο πεδίο του χρόνου. Στην περίπτωση που ο τύπος τους είναι ``ESH'', πρώτα εφαρμόζεται η συνάρτηση imdct4 σε όλα τα subframes, έπειτα γίνεται η παραθυροποίησή τους και στο τέλος αθροίζονται με κατάλληλο τρόπο ώστε στην έξοδο να επιστραφεί ολόκληρο το frame. -\subsection{\textit{function} AACSeq1 = AACoder1(fNameIn)} +\subsection{\texttt{function AACSeq1 = AACoder1(fNameIn)}} Για να εφαρμοστούν οι παραπάνω βαθμίδες, υλοποιήσαμε τον κωδικοποιητή AACoder1, ο οποίος παίρνει σαν μοναδικό όρισμα το path του WAV αρχείου που επιθυμούμε να κωδικοποιήσουμε, το οποίο πρέπει να είναι δειγματοληπτημένο στα 48KHz και να έχει 2 κανάλια, ενώ σαν έξοδο επιστρέφει ένα struct με διάσταση k x 1, όπου k, ο αριθμός των frames του αρχείου. @@ -96,13 +96,13 @@ Tα frame τα οποία είναι στατικά χαρακτηρίζοντα Στη συνέχεια, εφαρμόζουμε την filterbank σε όλα τα frames και αποθηκεύουμε για το κάθε ένα, τον τύπο του, τον τύπο του παραθύρου που χρησιμοποιήθηκε καθώς και τους συντελεστές MDCT για κάθε κανάλι. Όλα αυτά αποθηκεύονται στο struct AACSeq1. -\subsection{\textit{function} x = iAACoder1(AACSeq1, fNameOut)} +\subsection{\texttt{function x = iAACoder1(AACSeq1, fNameOut)}} Αποτελεί την αντίστροφη συνάρτηση της AACoder1, υλοποιεί δηλαδή έναν αποκωδικοποιητή. Δέχεται ως όρισμα, τον πίνακα με structs που επιστρέφεται από τον AACoder1 και το όνομα που θα δοθεί στο αποκωδικοποιημένο αρχείο μαζί με το path του, ενώ η έξοδός του είναι ένας πίνακας που περιέχει τα δείγματα του αποκωδικοποιημένου αρχείου ήχου. Κατά την αποκωδικοποίηση, τροφοδοτούνται οι συντελεστές MDCT των καναλιών στην iFilterbank, μαζί με τον τύπο του κάθε frame και τον τύπο του παραθύρου και αφού προστεθούν όλα τα frames, επιστρέφεται το αρχείο ολόκληρο στο πεδίο του χρόνου και αποθηκεύεται. - \subsection{\textit{function} SNR = demoAAC1(fNameIn, fNameOut)} + \subsection{\texttt{function SNR = demoAAC1(fNameIn, fNameOut)}} Στην συνάρτηση αυτή, επιδεικνύεται η λειτουργία του κωδικοποιητή και του αποκωδικοποιητή και επιστρέφεται ο σηματοθορυβικός λόγος μεταξύ του αποκωδικοποιημένου και του αρχικού αρχείου ήχου. Δέχεται ως ορίσματα το path του αρχείου που επιθυμούμε να κωδικοποιήσουμε και το όνομα που θα του δοθεί μετά την αποκωδικοποίηση. Για την εύρεση του σηματοθορυβικού λόγου χρησιμοποιείται η συνάρτηση του MATLAB snr. @@ -139,28 +139,28 @@ Tα frame τα οποία είναι στατικά χαρακτηρίζοντα Στο 2ο επίπεδο, υλοποιείται η βαθμίδα Temporal Noise Shaping (TNS), ένας κωδικοποιητής και οι αντίστροφες συναρτήσεις τους. -\subsection{\textit{function} [frameFout, TNScoeffs] = TNS(frameFin, frameType)} +\subsection{\texttt{function [frameFout, TNScoeffs] = TNS(frameFin, frameType)}} Η συνάρτηση TNS υλοποιεί τη βαθμίδα Temporal Noise Shaping για ένα κανάλι. Δέχεται ως εισόδους τους κβαντισμένους συντελεστές MDCT των 2 καναλιών ενός frame και τον τύπο του ενώ επιστρέφει σαν έξοδο το frame στο πεδίο συχνότητας, σε όρους MDCT συντελεστών έπειτα από την εφαρμογή του TNS καθώς και τους κβαντισμένους συντελεστές του TNS. Σαν πρώτο βήμα, αρχικοποιούμε τον πίνακα των μπαντών του ψυχοακουστικού μοντέλου όπως δίνεται στο πρότυπο. Στη συνέχεια, κατατάσσουμε του συντελεστές MDCT στις μπάντες και βρίσκουμε την ενέργεια κάθε μπάντας. Έπειτα, υπολογίζουμε τους συντελεστές κανονικοποίησης, τους εξομαλύνουμε όπως υποδεικνύεται και κανονικοποιούμε τους συντελεστές MDCT ανάλογα με την ενέργεια της μπάντας στην οποία ανήκουν. Με την βοήθεια της συνάρτησης lpc του MATLAB, βρίσκουμε τους συντελεστές γραμμικής πρόβλεψης για κάθε frame (ή subframe αν ο τύπος του είναι "ESH") και τους κβαντίζουμε με τη βοήθεια της συνάρτησης quantiz για ομοιόμορφο κβαντιστή βήματος 0.1. Εφαρμόζουμε το FIR φίλτρο $H_{TNS}(z) = 1 - a_{1}z^{-1} - a_{2}z^{2} - a_{3}z^{-3} - a_{4}z^{-4}$ στους αρχικούς συντελεστές MDCT όπου $a_{i}$ οι κβαντισμένοι συντελεστές γραμμικής πρόβλεψης. -\subsection{\textit{function} frameFout = iTNS(frameFin, frameType, TNScoeffs)} +\subsection{\texttt{function frameFout = iTNS(frameFin, frameType, TNScoeffs)}} Αποτελεί την αντίστροφη συνάρτηση της TNS. Δέχεται ως όρισμα το frame στο πεδίο του χρόνου σε όρους συντελεστών MDCT μετά το TNS, τον τύπο του και τους συντελεστές TNS, ενώ ως έξοδο επιστρέφει τους συντελεστές MDCT πριν το TNS. Για την εύρεση τους χρησιμοποιεί το αντίστροφο φίλτρο $H_{TNS}(z) = \frac{1}{1 - a_{1}z^{-1} - a_{2}z^{2} - a_{3}z^{-3} - a_{4}z^{-4}}$ -\subsection{\textit{function} AACSeq2 = AACoder2(fNameIn)} +\subsection{\texttt{function AACSeq2 = AACoder2(fNameIn)}} Υλοποιεί έναν κωδικοποιητή ο οποίος ουσιαστικά εκτελεί την συνάρτηση TNS. Δέχεται ως όρισμα το path του δείγματος ήχου που επιθυμούμε να κωδικοποίησουμε και επιστρέφει ένα struct μεγέθους K x 1, όπου Κ το πλήθος των frames που έχουν κωδικοποιηθεί και παριέχει για το κάθε ένα τον τύπο του, τον τύπο παραθύρου που χρησιμοποιήθηκε κατά την κωδικοποίηση, τους συντελεστές MDCT και τους κβαντισμένους συντελεστές TNS για κάθε κανάλι. Όπως και στον παραπάνω κωδικοποιητή, αρχικά διαβάζεται το αρχείο, χωρίζεται σε frames και με την βοήθεια της συνάρτησης SSC, καθορίζεται ο τύπος τους με την ίδια αντιμετώπιση για τα οριακά frames. Στη συνέχεια, τροφοδοτούνται τα frames ως είσοδοι στην συνάρτηση filterbank και στην συνάρτηση TNS και αποθηκεύονται τα αποτελέσματα στο struct. -\subsection{\textit{function} x = iAACoder2(AACSeq2, fNameOut)} +\subsection{\texttt{function x = iAACoder2(AACSeq2, fNameOut)}} Υλοποιεί την αντίστροφη του κωδικοποιητή AACoder1, δηλαδή έναν αποκωδικοποιητή. Δέχεται ως εισόδους το K x 1 struct του κωδικοποιημένου αρχείου καθώς και το όνομα που πρόκειται να πάρει μετά την αποκωδικοποίησή του. Ως έξοδο, επιστρέφει τα αποκωδικοποιημένα δείγματα του αρχείου και τα αποθηκεύει σε μορφή WAV. Για να το επιτύχει αυτό, χρησιμοποιεί τηις συναρτήσεις iTNS και iFilterbank. -\subsection{\textit{function} SNR = demoAAC2(fNameIn, fNameOut)} +\subsection{\texttt{function SNR = demoAAC2(fNameIn, fNameOut)}} H συνάρτηση αυτή επιδεικνύει την κωδικοποίηση του 2ου επιπέδου. Δέχεται ως ορίσματα το path του αρχείου ήχου που πρίκειται να κωδικοποιηθεί και το όνομα που πρόκειται να πάρει μετά την αποκωδικοποίησή του ενώ επιστρέφει τον συνολικό σηματοθορυβικό λόγο σε dB. Χρησιμοποιεί τις συναρτήσεις AACoder2 και iAACoder2 για την κωδικοποίηση και αποκωδικοποίηση του αρχείου ενώ για την εύρεση του σηματοθορυβικού λόγου χρησιμοποιεί την συνάρτηση snr του MATLAB. @@ -179,8 +179,8 @@ H συνάρτηση αυτή επιδεικνύει την κωδικοποίη Στο 3ο επίπεδο υλοποιείται το ψυχοακουστικό μοντέλο, οι βαθμίδες του κβαντιστή και του Huffman ένας κωδικοποιητής και οι αντίστροφές τους συναρτήσεις. -\subsection{\textit{function} SMR = psycho(frameT, frameType, frameTprev1, frameTprev2) -} +\subsection{\texttt{functionSMR = psycho(frameT, frameType, frameTprev1, frameTprev2) +} } Υλοποιεί την βαθμίδα του ψυχοακουστικού μοντέλου για ένα κανάλι. Δέχεται ως είσοδο το frame στο πεδίο του χρόνου, τον τύπο του και τα 2 προηγούμενα του frames στο ίδιο κανάλι. Σαν έξοδο επιστρέφει τον λόγο σήματος προς μάσκα (Signal to Mask Ratio - SMR). @@ -207,7 +207,7 @@ H συνάρτηση αυτή επιδεικνύει την κωδικοποίη -\subsection{\textit{function} [S, sfc, G] = AACquantizer(frameF, frameType, SMR)} +\subsection{\texttt{function [S, sfc, G] = AACquantizer(frameF, frameType, SMR)}} Υπολογίζει εσωτερικά το κατώφλι ακουστότητας T(b) και υλοποιεί τη βαθμίδα του κβαντιστή για κάθε κανάλι. Δέχεται ως εισόδους το frame στο πεδίο της συχνότητας σε όρους συντελεστών MDCT, τον τύπο του και το υπολογισμένο SMR. Σαν έξοδο επιστρέφει τον πίνακα s που περιέχει τα σύμβολα κβάντισης των συντελεστών ΜDCT του frame, τον πίνακα sfc με τις τιμές των Scalefactor για κάθε Scalefactor band και το G που είναι το global gain του frame. @@ -236,11 +236,12 @@ H συνάρτηση αυτή επιδεικνύει την κωδικοποίη \label{ref:5} \end{figure} -\subsection{\textit{function} frameF = iAACquantizer(S, sfc, G, frameType)} +\subsection{\texttt{function frameF = iAACquantizer(S, sfc, G,} frameType)} Υλοποιεί την αντίστροφη συνάρτηση της AACquantizer. Δέχεται ως ορίσματα τα σύμβολα κβάντισης των συντελεστών MDCT, τον πίνακα με τα scalefactors κάθε μπάντας και το global gain. Σαν έξοδο, με τους συντελεστές και το scalefactor προσπαθεί να επανακατασκαευάσει το frame στο πεδίο της συχνότητας σε όρους συντελεστών MDCT, το οποίο και επιστρέφει. -\subsection{\textit{function} [huffSec, huffCodebook] = encodeHuff(coeffSec, huffLUT, forcedCodebook)} +\subsection{\texttt{function [huffSec, huffCodebook] = encodeHuff(coeffSec, huffLUT, forcedCodebook)}} + Υλοποιεί τη βαθμίδα κωδικοποίσης Huffman. Δέχεται ως εισόδους κβαντισμένους συντελεστές, τους Huffman look-up tables και τον δείκτη Codebook, όπως αυτός παρουσιάζεται στο πρότυπο. @@ -248,23 +249,23 @@ H συνάρτηση αυτή επιδεικνύει την κωδικοποίη Για την υλοποίησή της ακολουθήθηκαν οι οδηγίες από το πρότυπο ενώ για να χρησιμοποιηθούν οι look-up tables κλήθηκε η συνάρτηση loadLUT.m. -\subsection{\textit{function} decCoeffs = decodeHuff(huffSec, huffCodebook, huffLUT)} +\subsection{\texttt{function decCoeffs = decodeHuff(huffSec, huffCodebook, huffLUT)}} Υλοποιεί την αντίστροφη συνάρτηση της encodeHuff, δηλαδή αποτελεί έναν αποκωδικοποιητή Huffman. Δέχεται ως είσοδο το string με '1' και '0' που αντιστοιχεί στην κωδικοποίηση Huffman, τους δείκτες codebook που χρησιμοποιήθηκαν και τους Huffman look-up tables. Σαν έξοδο, επιστρέφει τους αποκωδικοποιημένους κβαντισμένους συντελεστές. Η διαδικασία είναι και πάλι αυτή που παρουσιάζεται στο πρότυπο ενώ λόγο του zero padding που είναι απαραίτητο, το μήκος της εξόδου, δηλαδή των κβαντισμένων συντελεστών, ίσως είναι μεγαλύτερο από αυτό που θα έπρεπε. Σε αυτήν την περίπτωση οι τιμές που είναι εκτός του εύρους του αναμενόμενου μήκους, θα πρέπει να αγνοηθούν ενώ αναμένεται να είναι και μηδενικοί -\subsection{\textit{function} AACSeq3 = AACoder3(fNameIn, fnameAACoded)} +\subsection{\texttt{function AACSeq3 = AACoder3(fNameIn, fnameAACoded)}} Υλοποιεί έναν κωδικοποιητή ο οποίος δέχεται σαν είσοδο το path του αρχείου το οποίο πρόκειται να κωδικοποιηθεί και το path του .mat αρχείου που πρόκειται να εγγραφεί μετά την κωδικοποίηση. Ως έξοδο, επιστρέφει ένα struct με διάσταση k x 1, όπου k, ο αριθμός των frames του αρχείου, το οποίο περιέχει για κάθε frame, όπως και πριν τον τύπο του και το παράθυρο που χρησιμοποιήθηκε,ενώ για κάθε κανάλι επιστρέφει τους κβαντισμένους συντελεστές TNS, τα κατώφλια ακουστότητας λόγω του ψυχοακουστικού μοντέλου, τα κβαντισμένα global gains, τους κωδικοποιημένους με Huffman κβαντισμένους συντελεστές MDCT και scalefactors και το Huffman codebook. Φορτώνει αρχικά τους πίνακες TNS, διαβάζει το αρχείο προς κωδικοποίηση και αφού το χωρίσει σε frames χρησιμοποεί στη συνέχεια τις συναρτήσεις SSC, filterbank, TNS, psycho, AACquantizer και encodeHuff για να το κωδικοποίησει. Τα αποτελέσματα φαίνονται παρακάτω: -\subsection{\textit{function} x = iAACoder3(AACSeq3, fNameOut)} +\subsection{\texttt{function x = iAACoder3(AACSeq3, fNameOut)}} Υλοποιεί τον αντίστοιχο αποκωδικοποιητή. Δέχεται ως εισόδους, το struct που προκύπτει από τον κωδικοποιητή AACoder3 και το όνομα μαζί με το path που πρόκειται να πάρει το αποκωδικοποιημένο αρχείο ενώ ως έξοδο επιστρέφει έναν πίνακα με το αποκωδικοποιημένο σήμα. Για την υλοποίηση του χρησιμοποιούνται οι συναρτήσεις decodeHuff, iAACquantizer, iTNS και iFilterbank. -\subsection{\textit{function} [SNR, bitrate, compression] = demoAAC3(fNameIn, fNameOut, frameAACoded)} +\subsection{\texttt{function [SNR, bitrate, compression] = demoAAC3(fNameIn, fNameOut, frameAACoded)}} Ουσιαστικά η συνάρτηση αυτή παρουσιάζει την λειτουργία των συναρτήσεων του επιπέδου 3. Δέχεται ως εισόδους το path του αρχείου που πρόκειται να κωδικοποιηθεί και το όνομα του αρχείου που πρόκειται να επιστραφεί μετά την αποκωδικοποίηση. Σαν έξοδο, επιστρέφει τον σηματοθορυβικό λόγο του σήματος μετά την κωδικοποίηση και αποκωδικοποίησή του, το bitrate και το ποσοστό συμπίεσης.