LAB2 5.4

[giorgoskallitsis]

Καλημερα, προσπαθησα να προσεγγισω το φασμα ισψυος των 2 σηματων με τη μεθοδο welch. Αρχικα αυτα που εχω κανει ειναι τα εξης: Παιρνω δειγματα απο το αρχικο σημα μηκους 100 σημειων και τα πολλαπλασιαζω με το παραθυρο bartlet. Εν συνεχεια, βρισκω το φασμα ισχυος αυτων με τη μεθοδο του περιοδογραμματος και το αποτελεσμα το βαζω σε μια λιστα. Καθε φορα που βρισκω ενα καινουριο φασμα ισχυος, το προσθετω στην λιστα μου στις αντιστοιχες χρονικες στιγμες (100) και στο τελος διαιρω ολα τα στοιχεια της λιστας με το πληθος των παραθυρωμενων δειγματων. Δεν ειμαι βεβαιος αν κανω κατι λαθος. Για επιβεβαιωση, εφαρμοσα στο σημα μου την ετοιμη συναρτηση welch για να δω αν ειμαι κοντα στο επιθυμητο αποτελεσμα. Παρατηρω πως αν εξαιρεσουμε το γεγονος πως η δικη μου υλοποιηση καταλαμβανει ολο το ευρος συχνοτητων απο 0 εως fs (απο fs/2 εως fs ειναι πρακτικα συμμετρικη), το αποτελεσμα ειναι κοντινο. Γνωριζω πως το φασμα ισχυος εκτεινεται απο 0 εως fs/2. Θα πρεπει να μηδενισω τα δειγματα απο fs/2 εως fs ή μπορω να το αφησω να ειναι συμμετρικο και απλως να κοιταω το πρωτο μισο?

[cgaroufis]

Καλησπέρα,

1. Η αρχική σου διαδικασία είναι ορθή, εφόσον τηρείς το overlap/hop size 50 δειγμάτων ανάμεσα στα διαφορετικά παράθυρα του σήματός σου.
2. Από τη στιγμή που η χρήσιμη πληροφορία του περιοδογράμματός σου είναι στο διάστημα μέχρι fs/2, μπορείς είτε να κοιτάς μόνο το αριστερά μισό του περιοδογράμματος κρατώντας το ολόκληρο, είτε να πλοτάρεις μόνο το αριστερά μισό.

Καλή συνέχεια, Χρήστος

[giorgoskallitsis]

Καλημερα, σας ευχαριστω πολυ!