THMMY.gr

(Αν ασχολείται κάποιος) Τι παίζει με τα σφάλματα στο proof.m ?
DFT : 1.244156e-014
split output top bottom : 1.244156e-014
split input even odd : 1.288815e-014
done : 1.256074e-015

πρέπει και στο done να βγει σαν τα υπόλοιπα ?

Λογικά πρέπει να φτιάξεις τη σχέση που λέει στην εκφώνηση και να τη συγκρίνεις με το ys όπως στο script (fprintf). Το σφάλμα θα πρέπει να είναι στην ίδια τάξη (?) με το done.

Ναι το ξεδιάλυνα κι εγώ. Αν το σφάλμα είναι της τάξης του 10^-13 και κάτω είσαι μια χαρά. Εμένα κάπου εκεί μου βγήκε, όχι ίδιο πάντως με το "done".

Είμαι ο μόνος που η συνάρτηση plant μου επιστρέφει ό,τι της δίνω ως είσοδο;  :D

Κι εμένα το ίδιο :D

x = cos(1:0.1:100);
d = plant(x);
plot([x' d'])

Δε βγάζει την είσοδο..

Έχεις δίκιο. Της έδινα εισόδους μικρού μήκους γι αυτό έπαιρνα έξοδο ίδια με την είσοδο. Ανακαλώ!   8))

Στον block LMS πόσο παίρνουμε το kmax και το L; Κωλομάθημα σιχτίρ.

Στη σελίδα 10 του BlockLMS.pdf που υπάρχει στο ethmmy στο υλικό του μαθήματος στον φάκελο "Διαλέξεις"
αναφέρει :

"Choice of block size
In most application the block size is selected to be equal to the filter length L = M . It is a tradeoff of
the following drawbacks:
– For L > M the gradient is estimated using more data than the filter itself.
– For L < M the data in the current block is not enough to feed the whole tap vector, and consequently
some weights are not used."

Σε ότι αφορά την επιλογή του kmax είναι δικό σου θέμα. Δοκίμασε διάφορες τιμές και κράτησε κάποια για την οποία ο
αλγόριθμός σου θα προλαβαίνει να εντοπίσει τις τιμές των παραμέτρων. Για παράδειγμα, εάν επιλέξεις kmax = 5, κατά πάσα
πιθανότητα ο αλγόριθμος δε θα προλάβει να εντοπίσει τις σωστές τιμές των παραμέτρων.

Ναι το είδα αυτό, οπότε L = 32; Ο αριθμός των w θεωρείται το filter length?

Επίσης δε θα πρέπει το kmax*L να είναι ίσο με το n (time steps) ?

Υλοποιησε κανεις τον block lms ? Εχετε συγκλιση για μικρες τιμες του μ ? γιατι σε μενα συγκλινει για μ = 1000 πχ   :D

Μόνο εγώ δεν κατάλαβα το πρώτο ερώτημα;

Στο proof.m υπολογίζει το DFT με 3 τρόπους. 1-με την fft, 2-συμφωνα με το 3ο μέρος της εξίσωσης που δίνεται, 3-σύμφωνα με το 4ο μέρος της εξίωσης που δίνεται
Τα αποτελέσματα των παραπάνω έχουν μικρές αποκλισεις.

Εμάς τι ακριβώς μας ζητάει να συμπληρώσουμε;

Πρεπει να υλοποιησεις τα αθροισματα της 2ης γραμμης απο τον τυπο που σου δινει. Αν παρατηρησεις το done ειναι το ιδιο αλλα χρησιμοποιει την fft() της matlab. Ουσιαστικα θες σφαλμα της ταξης του 10^-13 η μικροτερο.

Έχει καταλάβει κανείς στο β) ερώτημα η αναδρομική συνάρτηση που θα καλείται κάθε φορά από τον εαυτό της τι θα υπολογίζει ακριβώς; Θα χωρίζει το αρχικό διάνυσμα σε 2, 4, 8, ... διανύσματα και θα υπολογίζει κάθε ένα ξεχωριστά, όπως και στο αρχείο fftproof.m που χωρίζει σε πάνω και κάτω μέρος το διάνυσμα;

Υπολογίζω τα αθροίσματα όπως τα ζητάει. Αποθηκεύω ένα προηγούμενο y σε μια temp.
Συγκρίνω το y που υπολογίσα, με τα προηγούμενα και ενώ είναι ίδια (τουλάχιστον τα 4 δεκαδικα που δείχνει το matlab)
Αν κάνω y-temp έχει κάποιες διαφορές σε δεκαδικά που δε φαίνονται και το τελικό αποτέλεσμα προφανώς βγαίνει λάθος

Είναι σωστός ο τρόπος ή θέλει διόρθωση;

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

Κάθε μέρος της ισότητας είναι μια απόδειξη. Αυτός χρησιμοποιεί όλους τους ορισμούς εκτός απ' τον τελευταίο με τους πίνακες όπου τον αποδεικνύει με το:


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

Για το blockLMS έχει κανείς ιδέα που σκατά κάνω λάθος;

Νομίζω δεικτοδοτείς λάθος το u. Πρεπει να γράψεις:


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

Έκανε κανείς το β από το 3 ? Στην ουσία θέλει να πετάξουμε τον εσωτερικό βρόχο και να κάνουμε την ίδια δουλειά με πίνακες ή είναι κάτι αλλο?

Στο β ερώτημα ξέρει κανείς τι πρέπει να κάνουμε?

Το plant.p πως το χειριζόμαστε? δεν μπορούμε να το ανοιξουμε,πως το χρησιμοποιούμε?

Το plant.p είναι ένα γραμμικό σύστημα. Σκέψου το σαν μια συνάρτηση μεταφοράς. Του δίνεις μια είσοδο και σου επιστρέφει μια διαφορετική έξοδο.

y = plant(x)  :)

Εκανε κανένας το Block Lms (nested loops)? τι βήμα παίρνουμε ? εμένα για οτιδήποτε πάνω από mu=0.1 δεν συγκλίνει ,παρόλο που με τις ιδιοτιμές βγάζω mu<4.4 .Για τιμές κάτω του 0.1 λειτουργεί κομπλέ

Κι εγώ το ίδιο

Και εμένα έτσι συγκλίνει. Με πολύ μικρό βήμα < 0.1 .

Πήραμε καμιά παράταση;

Μέχρι Κυριακή

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

(Απο ethmmy)

Στον block lms δεν έχουμε πρόβλημα για τις τιμές του u όταν κ=0?

       y(k*L + i) = w'*u(k*M+i:-1:k*M+i-M+1);

        e(k*L + i) = d(k*L + i) - y(k*L + i);
        fi = fi + mu*e(k*L + i)*u(k*M+i:-1:k*M+i-M+1);

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

Να ρωτήσω κάτι για το 2ο ερώτημα.. Το Τ(2) στη σχέση είπε ότι είναι 4; Γιατί εμένα δε μου βγαίνει για Τ(2)=4

Για την αναδρομική σχέση Τ(n)  αν μπορει καποιος να εξηγησει τι ειναι και πως
την χρησιμοποιουμε , γιατι δεν παρακολουθω το μαθημα

3η εργασια θα βαλει; Γιατι τοσο μεγαλη παραταση;

Λογικά 4/5 είναι το σωστό, μάλλον έγινε λάθος..

Για το 1ο ερώτημα ....
Θέλει στο κάτω μέρος τα fe, fo που υπολογίζει με fft να τα υπολογίσουμε με με x(n)W? ?....σαν αθροισμα γινομένων

Κάποιος;

Η αναδρομική σχέση που δίνει στην εκφώνηση δίνει τον αριθμό των flops του FFT σε συνάρτηση με το n (μήκος του σήματος εισόδου). Υπάρχουν διάφοροι τρόποι να υπολογίσεις το T(n) από τον κώδικα/αλγόριθμο (στο μάθημα των Αλγορίθμων τα είχαμε κάνει αλλά και με googlάρισμα είναι εύκολο να βρεις υλικό). Το ερώτημα ζητάει τη βελτίωση του FFT σε σχέση με τον DFT υπολογίζοντας τα T(n) του κάθε ενός.

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

βασικα αν δεις τα fe,fo τα υπολογιζει με την fft της matlab. εμεις πρεπει να τα υπολογισουμε απο τις σχεσεις που μας δινονται

Να φανταστώ χρησιμοποιείται το κριτήριο μ<2/λmax. Κάπου είδα οτι ειναι 0<μ<2/L*λmax


Τώρα βεβαια στο pdf του Haykin δίνει το 2/λmax αλλα τ να πω xD

Για να κάνω το ερώτημα με 1 loop στον LMS: πάλι δε χρειάζομαι 2ρη loop???πως δημιουργώ απευθειας τέτοιο πίνακα????