THMMY.gr

Το τόπικ αυτό προορίζεται αποκλειστικά για ερωτήσεις και ανακοινώσεις που αφορούν τις εργασίες του μαθήματος.

Οποιοδήποτε άσχετο ποστ θα διαγράφεται.

Στο πρώτο ερώτημα που ζητάει φασματικό περιεχόμενο μέσω DFT, θα έχουμε διάγραμμα με περιεχόμενο μόνο στις συχνότητες 0 και 500; (Λόγω συχνότητα δειγματοληψίας στα 500Hz)

εμένα έτσι μου  βγήκε..αν και νομίζω είναι περιοδικός ο DFT,και γύρω από τα 250 βγαίνει επίσης συμμετρία απ ότι νομίζω

Αυτό έβγαλα εγώ. Μου φάνεται κάπως... το διάγραμμα

αν δεν έχεις περιεχόμενο στα 50 ξανακοιτά το

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

Επίσης για την εξασθένιση ενός φίλτρου, υπάρχει συγκεκριμένη διαδικασία για να βρούμε την ιδανική ή τη βρίσκουμε πειραματικά;

Να κάνω μία ερώτηση; Αφού βρείτε ας πούμε τον μετασχηματισμό z του φίλτρου ή τους συντελεστές των πολυωνύμων τέλος πάντων μετά πως φιλτράρεται το σήμα; Εγώ χρησιμοποιώ την filter() αλλά για τάξη πάνω από 7 το φίλτρο που βγαίνει είναι κάπως προβληματικό. Για χαμηλότερες τάξεις κάτι γίνεται αλλά το σήμα που βγαίνει δεν μοιάζει και πολύ με το αναμενόμενο.Έχει δηλαδή μία κυμάτωση και αυτό όχι σαν το αρχικό, πιο υψίσυχνη και χαμηλότερου πλάτους αλλά το αποτέλεσμα δεν είναι καλό. Εσείς πως το αντιμετωπίζετε;

Δοκίμασε και την filtfilt..Μπορεί να κάνει για τη δουλειά που θες. Αν και αν είσαι στο πρώτο φιλτράρισμα (στο υψιπερατό δηλαδή), όντως χαλούσε το αρχικό σήμα από μια τάξη και μετά

Πράγματι για το Highpass ρωτούσα. Η filtfilt τελικά λειτούργησε πολύ καλύτερα από την filt για το συγκεκριμένο με την ίδια διαδικασία υπολογισμού...ευχαριστώ!

'H απλά φτιάχνετε απευθείας ψηφιακό φίλτρο με το fdatool 8))

Προτείνω αν θέλετε να κάνετε το ερώτημα 3 από μέρος ΙΙ να μπείτε εδώ http://www.physionet.org/cgi-bin/atm/ATM και να διαλέξετε πάνω δεξιά export σε .mat γιατί εγώ έφαγα αρκετή ώρα να το βρώ

ή αυτό
http://www.physionet.org/physiotools/matlab/rddata.m

Τα annotation των γιατρών ξερει κανείς που τα βρίσκουμε για να μπορεσουμε να τα συγκρίνουμε με τα αποτελεσματα μας?

Νομίζω ότι είναι τα αρχεία .hea της δεξιάς στήλης από εδώ (http://www.physionet.org/physiobank/database/mitdb/)

Με την αυτοσυσχέτιση ξέρει κανείς τι γίνεται και γιατί πρέπει να την ενσωματώσουμε στον αλγόριθμο;

την χρησιμοποιεις για να βρεις  την περιοδικοτητα του σηματος εγω αυτο καταλαβα

δηλαδή σημαίνουν κάτι αυτοί οι πεταμένοι αριθμοί  :-\

Έτσι κατάλαβα
Αν πχ έχεις αυτό
στην πρώτη γραμμή πχ 100 είναι το όνομα του αρχείου, 360 η συχνότητα δειγματοληψίας, 650000 ο αριθμός των δειγμάτων
Στην πρώτη #, 69 πρέπει να είναι ο καρδιακός ρυθμός
Για τα άλλα δεν ξέρω..

Επίσης δείτε και εδώ : http://physionet.org/physiobank/database/html/mitdbdir/mitdbdir.htm
Παίζει να έχει και άλλες πληροφορίες για κάθε εγγραφή

αυτο βρήκα στο ιντερνετ, λεει και  πιο αρχειο ειναι αυτο που θες  :)

In the mit bih database (i) click ATM (ii) in the input column select MIT BIH arry database (iii)select the signals in record(u have number of signals) (iv)in the signals, select any one either v5 or ml11 (v) in the toolbox column select the 'export signals as .mat (vi) then click download (vii) a matlab file will be downloaded and drag this file in the workspace. and plot it. thats it.

Από το παραπάνω site έχει κανείς ίδεα ποιο από τα 2 σήματα χρειαζόμαστε?Για κάποια από τα 100,102,104,106 και τα 2 σήματα μοιάζουν αρκετά,για άλλα όπως το 106 είναι αρκετά διαφορετικά...

κράτα όποιο σου βγάλει πιο νορμάλ καρδιακό ρυθμό σε σχέση με τα annotation files :P

Οπως λέει και ο θεοχάρης:

Του γκουντ του μπι τρου...

Βρήκα κι αυτό http://www.physionet.org/lightwave/ (http://www.physionet.org/lightwave/)

Το φίλτρο μπορούμε να το φτιάξουμε χρησιμοποιώντας τα tools που δίνει το matlab για σχεδίαση φίλτρων ή πρέπει να το κάνουμε όπως περιγράφεται στα τελευαταία κεφάλαια του βιβλίου;

Όπως θέλεις,αρκεί να κάνει τη δουλειά του.

όσον αφορά το ερώτημα ΙΙ.3, που βρίσκουμε πληροφορίες για το QRS των καταγραφών?

Στα .hea αρχεία έχει το bpm:

ναι αυτό το βρηκα, αλλα μπορούμε απο εκεί να πάρουμε πληροφορίες για το που εκτείνεται το QRS?

Νομίζω ότι το μόνο που χρειάζεσαι είναι το διάστημα R-R',όχι η θέση των Q,S

λες ε? επειδή απ αυτό που λέει "συγκρίνετε τα αποτελέσματα των αλγορίθμων σας ΙΙ-1 και ΙΙ-2" και δεδομενου οτι το ΙΙ-2 λεει  να βρούμε το συμπλεγμα QRS, φαντάζομαι οτι θα θέλει κι αυτή την πληροφορία. Ή απλα παιδεύομαι τζαμπα  :P

Αν δεις και τα διαγράμματα στο άλλο pdf (εκει που έχει τα blocks και απο διπλά τις αλλάγες στο σήμα),δε νομίζω ότι μπορείς να βρεις κάτι άλλο.

Και εγώ αυτό του λέω, θα παραήταν δύσκολο μετά ακόμα και αν υπήρχαν μεσα στα Attributes

μαλιστα, ευχαριστώ

"Ασχετo",αλλά ο τύπος για το bpm είναι ο ανάστροφος από αυτόν που έχει η εκφώνηση,δηλαδή είναι 60/(R-R') και μετά επι fs για είναι σε σωστές μονάδες (thx teslaaaaaaa και internet).

Δεν κατάλαβα γιατί πρέπει να ενσωματώσουμε αυτοσυσχέτιση για να βρούμε το bpm. Έκανα έναν αλγόριθμο που βρίσκει που μετράει το χρόνο μεταξύ 2 peaks. Γιατί θέλει υποχρεωτικά αυτοσυσχέτιση;

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

Ναι με threshold το κανω και νομιζω βγαζει σωστο νουμερο για το ecg4.
Τα φιλτρα θσ εφαρμοστουν στο κανονικοποιημενο η στο αρχκικο ecg;

Ποια εντολη πρεπει να χρησιμοποιησω γιααυτοσυσχετιση; η autocorr μου δινει μια φθινουσ καμπυλη κ η xcorr βγαζει ενα ημιτονοειδες συμμετρικο plot. Δεν προσφερεται καποιο για μετρηση

Απ οτι καταλαβαίνω, τα φιλτρα εφαρμοζονται στο κανονικοποιημενο ecg. Για αυτοσυσχέτιση την autocorr

Δηλαδή εφαρμόζω πρώτα την αυτοσυσχέτιση και στη συνέχεια τρέχω τον αλγόριθμο;

Οι αποστάσεις μεταξύ των R-R' δεν είναι σταθερές,οπότε θα πάρουμε μέσο όρο;

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

Το 69 είναι η ηλικία και το M είναι για τον άνδρα(male), άμα δείτε κι εδώ: http://physionet.org/physiobank/database/html/mitdbdir/records.htm#101

σοβαρά? ;D
σόρρυ guys :-[

 :angry4: