[Ψηφιακά Φίλτρα] Απορίες σε ασκήσεις 2019

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

Για οποιοδήποτε πρόβλημα με register/login, στείλτε email στο contact@thmmy.gr.

[chamezos]

σωστο..για να υπολογισεις το R = 100x100 δεν χρειαζεσαι ενα σημα u = 100x1 το οποίο θα δημιουργηθεί απο μια AR διαδικασια?

ΑΝ τα θυμάμαι σωστα, με μαθηματικά καφενείου ισχύει το εξής: Το σήμα u είναι μια AR διαδικασία. Συνεπώς το u(n) θα συσχετίζεται με τον εαυτό του (δηλαδη θα αυτοσυσχετίζεται) με το u(παρελθοντικής στιγμής). Εφόσον απο τον τύπο της εκφώνησης βλέπουμε οτι υπάρχει συσχέτιση μέχρι και 3 χρονικές στιγμές πρίν, θα ορίζεται αυτοσυσχέτιση μέχρι 3 χρονικές στιγμές πρίν. Από 4 χρονικές στιγμές και πίσω το u(n) δεν επηρεάζει την τιμή του, άρα και η αυτοσυσχέτιση είναι 0

[delta]

ΑΝ τα θυμάμαι σωστα, με μαθηματικά καφενείου ισχύει το εξής: Το σήμα u είναι μια AR διαδικασία. Συνεπώς το u(n) θα συσχετίζεται με τον εαυτό του (δηλαδη θα αυτοσυσχετίζεται) με το u(παρελθοντικής στιγμής). Εφόσον απο τον τύπο της εκφώνησης βλέπουμε οτι υπάρχει συσχέτιση μέχρι και 3 χρονικές στιγμές πρίν, θα ορίζεται αυτοσυσχέτιση μέχρι 3 χρονικές στιγμές πρίν. Από 4 χρονικές στιγμές και πίσω το u(n) δεν επηρεάζει την τιμή του, άρα και η αυτοσυσχέτιση είναι 0

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

[A Caster]

ΑΝ τα θυμάμαι σωστα, με μαθηματικά καφενείου ισχύει το εξής: Το σήμα u είναι μια AR διαδικασία. Συνεπώς το u(n) θα συσχετίζεται με τον εαυτό του (δηλαδη θα αυτοσυσχετίζεται) με το u(παρελθοντικής στιγμής). Εφόσον απο τον τύπο της εκφώνησης βλέπουμε οτι υπάρχει συσχέτιση μέχρι και 3 χρονικές στιγμές πρίν, θα ορίζεται αυτοσυσχέτιση μέχρι 3 χρονικές στιγμές πρίν. Από 4 χρονικές στιγμές και πίσω το u(n) δεν επηρεάζει την τιμή του, άρα και η αυτοσυσχέτιση είναι 0

αυτο που λες ειναι μερική αυτοσυσχέτιση. Μπορείς να δεις κάποια αναλυτικότερα πράγματα εδώ :
http://users.auth.gr/dkugiu/Teach/TimeSeriesTHMMY/TimeSeriesNotes.pdf
Στα κεφάλαια 3.2 και 3.3 (τα σχηματάκια φτάνουν).

Τωρα στο προτυπο οντως τα εχει βαλει μαλλον με το χερι 0. Εγω για να τα υπολογίσω έφτιαξα σε 2 λουπες τον πίνακα με τα στοιχεία του a ο οποιος πολλαπλασιάζεται με το διάνυσμα των αυτοσυσχετισεων για να βγαλει το δεξι κομματι (που ειναι: το πρωτο στοιχειο το s_v1^2 kαι τα υπολοιπα μηδενικα)

Δηλαδή έκανα το A*r = [S^2; zeros(99,1)]  ( Yule - Walker Equations)

και βρήκα το r = A \ [S^2; zeros(99,1)] , R = toeplitz(r).

Προφανως το πηρε λαθος αλλα αυτός ειναι ο σωστός (μόνος?) τρόπος να υπολογίζεις αυτοσυσχετίσεις εως αυθαίρετα μεγάλο αριθμό νομίζω.

[JohnFli]

Παιδιά είναι η  άσκηση 1.1 από το φυλλάδιο των ασκήσεων.

Αρκεί απλά να κάνετε r = [r ; zeros(96,1)]; με r = [ r(1); r(2); r(3); r(4) ] όπως τα υπολογίσατε από τις Yule-Walker και φυσικά R = toeplitz(r);

[thankara]

Τελικά έχουμε κολλήσει στο 0.63/1.00 έτσι? Υπάρχει κάποιος που πήρε παραπάνω? Έχω συγκρίνει όλους του αλγορίθμους με την leastMeanSquareGold και παίρνω τα ίδια αποτελέσματα (με ένα σφάλμα 10^-4 που λογικά εξαρτάται από τις παραμέτρους εισόδου (alpha delta mu κλπ) που δίνω εγώ).

[leukosaraphs!]

Τελικά έχουμε κολλήσει στο 0.63/1.00 έτσι? Υπάρχει κάποιος που πήρε παραπάνω? Έχω συγκρίνει όλους του αλγορίθμους με την leastMeanSquareGold και παίρνω τα ίδια αποτελέσματα (με ένα σφάλμα 10^-4 που λογικά εξαρτάται από τις παραμέτρους εισόδου (alpha delta mu κλπ) που δίνω εγώ).

εδω και 3 μερες 
και ο φλωρος δεν απανταει στα μαιλ μου

[achariso]

Παιδιά, επειδή δεν έκανα το 2η εργασία ( τις έπιασα από την 3η και τώρα την τέταρτη ), τί είναι τα u_stream; Δηλ. τι θά έχουν τα args της συνάρτησης;

π.χ. για τον RLS, εάν τα u και τα d ήταν ολόκληρα τα σήματα ( όχι ως streams ), ο παρακάτω κώδικας έχω τεστάρει ότι δουλεύει και προσεγγίζει καλά τους συντ/τές Wiener:

Code:

function [y,w_out] = rls(u, d, m, lambda, delta)
% RLS adaptation
%   
  
  %persistent w
  %persistent u_stream
  %persistent P
  
    u = u(:);
    d = d(:);
    n = length(u);

    P = delta * eye(m);
    w_out = zeros(m, 1);
    y = zeros(size(u));
    e = zeros(size(u));

    for i = (m + 1):n
        ubar = u(i:-1:i - m + 1);
        y(i) = w_out.' * ubar;
        k = 1 / lambda * P * ubar / (1 + 1 / lambda * ubar' * P * ubar);
        e(i) = d(i) - y(i);
        w_out = w_out + k * e(i);
        P = 1 / lambda * P - 1 / lambda * k * ubar.' * P;
    end
  
end

Ωστόσο, τώρα μου το απορρίπτει, ( ενώ δέχεται την συνάρτηση wiener(), δηλ. μου δίνει 0.25/1.00 ).
Btw, ούτε εγώ χρησιμοποιώ το c και sigma_v2 στην wiener():

Code:

function [R, p, w_o] = wiener( a, b, ~, sigma_v1, ~ )

    %% Reformat args
    a = [1; a(:)];


    %% Υπολογισμός πίνακα αυτοσυσχέτισης
    A = [...
        a(1) a(2) a(3) a(4);...
        a(2) a(1)+a(3) a(4) 0;...
        a(3) a(2)+a(4) a(1) 0;...
        a(4) a(3) a(2) a(1);...
    ];
    
    B = [sigma_v1; 0; 0; 0];
    
    r = A \ B;
    R = toeplitz(r);    % πίνακας αυτοσυσχέτισης
    
    %% Υπολογισμός διανύσματος ετεροσυσχέτισης
    p = R * (-b');
    
    %% Υπολογισμός Wiener Coefficients
    w_o = R \ p;
  
end

[MrRobot]

Σου δίνει σαν όρισμα το u(n), δηλαδή το u σε κάθε στιγμή. Αυτό το περνάς στο u_stream. Στη προηγούμενη σελίδα έχω ανεβάσει ένα κομμάτι κώδικα για την LMS. Θα καταλάβεις τι πρέπει να κάνεις.

[dinis]

Παιδιά υπάρχει περίπτωση να έχει κάνει κάποιο λάθος με τον έλεγχο στο e-learning. Δεν γίνεται να έχουμε κολλήσει όλοι στο 0.63.

[A Caster]

Παιδιά υπάρχει περίπτωση να έχει κάνει κάποιο λάθος με τον έλεγχο στο e-learning. Δεν γίνεται να έχουμε κολλήσει όλοι στο 0.63.

ε εκει εχουμε καταλήξει όλοι, ότι μάλλον υπάρχει κάποιο πρόβλημα στον έλεγχο. Από ότι φαίνεται γίνονται 2 δοκιμές σε κάθε ενα από τα rls / lms / nlms και περνάει μόνο τον πρώτο, αρα παίρνουμε όλοι 0.13 από κάθε ένα, συν 0.25 από το wiener. 

[chamezos]

Παιδιά, επειδή δεν έκανα το 2η εργασία ( τις έπιασα από την 3η και τώρα την τέταρτη ), τί είναι τα u_stream; Δηλ. τι θά έχουν τα args της συνάρτησης;

π.χ. για τον RLS, εάν τα u και τα d ήταν ολόκληρα τα σήματα ( όχι ως streams ), ο παρακάτω κώδικας έχω τεστάρει ότι δουλεύει και προσεγγίζει καλά τους συντ/τές Wiener:

Code:

function [y,w_out] = rls(u, d, m, lambda, delta)
% RLS adaptation
%   
  
  %persistent w
  %persistent u_stream
  %persistent P
  
    u = u(:);
    d = d(:);
    n = length(u);

    P = delta * eye(m);
    w_out = zeros(m, 1);
    y = zeros(size(u));
    e = zeros(size(u));

    for i = (m + 1):n
        ubar = u(i:-1:i - m + 1);
        y(i) = w_out.' * ubar;
        k = 1 / lambda * P * ubar / (1 + 1 / lambda * ubar' * P * ubar);
        e(i) = d(i) - y(i);
        w_out = w_out + k * e(i);
        P = 1 / lambda * P - 1 / lambda * k * ubar.' * P;
    end
  
end

Ωστόσο, τώρα μου το απορρίπτει, ( ενώ δέχεται την συνάρτηση wiener(), δηλ. μου δίνει 0.25/1.00 ).
Btw, ούτε εγώ χρησιμοποιώ το c και sigma_v2 στην wiener():

Code:

function [R, p, w_o] = wiener( a, b, ~, sigma_v1, ~ )

    %% Reformat args
    a = [1; a(:)];


    %% Υπολογισμός πίνακα αυτοσυσχέτισης
    A = [...
        a(1) a(2) a(3) a(4);...
        a(2) a(1)+a(3) a(4) 0;...
        a(3) a(2)+a(4) a(1) 0;...
        a(4) a(3) a(2) a(1);...
    ];
    
    B = [sigma_v1; 0; 0; 0];
    
    r = A \ B;
    R = toeplitz(r);    % πίνακας αυτοσυσχέτισης
    
    %% Υπολογισμός διανύσματος ετεροσυσχέτισης
    p = R * (-b');
    
    %% Υπολογισμός Wiener Coefficients
    w_o = R \ p;
  
end

Τα ορίσματα των συναρτήσεων lms, rls, nlms είναι όλα scalar. Βασισμένος στη λύση της δεύτερης εργασίας, το u_stream είναι ένας vector κοινός για τις συναρτήσεις όπου στην αρχή κάθε συνάρτησης μπορείς να βάλεις το u(n) που έρχεται σαν όρισμα κάπως έτσι: u_stream = [u ; u_stream] και να έχεις ουσιαστικά το σήμα u = [ u(n) ; u(n-1) ;.... u(0)]. Δες τον κώδικα που δίνει αυτός σα λύση για την δεύτερη εργασία

[annaTZ]

Μπορει κανεις να μου εξηγησει τι κανω λαθος στην 4η ασκηση? Παιρνω 0.25, δηλαδη ο wiener βγαινει σωστος αλλα τα υπολοιπα οχι...
Δεν μπορω να καταλαβω γιατι, τα κανω οπως ο φλωρος στην απαντηση που ειχε δωσει σε προηγουμενη ασκηση

Code:

function [y,w_out] = lms(u, d, m, mu)
% LMS - Least Mean Square adaptation
%   
 
  persistent w
  persistent u_stream
  
  % initialize w
  if isempty(w)
    w = zeros( m, 1 );
  end
  
  u_stream(end+1, 1) = u;
 
  if length( u_stream ) < m                
    y = 0;   
  else                                  
    % compute output
    y = w' * u_stream(end:-1:end-m+1);  
    % get error
    e = d - y;
    
    if length( u_stream ) > m    
      % update weights
      w = w + mu * e * u_stream( end:-1:end-m+1 );      
    end    
  end
  
  if nargout > 1
    w_out = w;
  end

end

[BIC]

Code:

[quote author=annaTZ link=topic=74039.msg1297525#msg1297525 date=1560251202]
Μπορει κανεις να μου εξηγησει τι κανω λαθος στην 4η ασκηση? Παιρνω 0.25, δηλαδη ο wiener βγαινει σωστος αλλα τα υπολοιπα οχι...
Δεν μπορω να καταλαβω γιατι, τα κανω οπως ο φλωρος στην απαντηση που ειχε δωσει σε προηγουμενη ασκηση

[/quote]

Και εγώ ακριβώς το ίδιο θέμα έχω, ξέρει κανείς τι παίζει??  

[BIC]

Η διαφορά με τη 2η εργασία είναι ότι ενώ στη 2η όσο δεν είχες αρκετά δείγματα εισόδου (length(u_stream)<M) δεν υπολόγιζες φίλτρο και έδεινες έξοδο y=u, τώρα λέει από την αρχή να κάνεις την κανονική διαδικασία.

Δηλαδή αν βγάλουμε το if και αφήσουμε μόνο τον κώδικα που ήταν μέσα στο else θα τρέξει? (ρωτάω γιατί δεν έχω άλλες προσπάθειες να δοκιμάσω)

[annaTZ]

Δηλαδή αν βγάλουμε το if και αφήσουμε μόνο τον κώδικα που ήταν μέσα στο else θα τρέξει? (ρωτάω γιατί δεν έχω άλλες προσπάθειες να δοκιμάσω)

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