THMMY.gr

Επειδή λέω να αρχίσω να ασχοληθώ με την ΨΕΣ.  :D

Αυτά που ζητάει η εργασία θα γίνουν σε ματλαμπ; (αν θυμάμαι καλά αυτό είπαν).

Ναι,κώδικας και διαγραμμάτα σε matlab και μια παρουσίαση της εργασίας σε αρχείο κειμένου.

Σήμερα στο τελευταίο μάθημε θεωρίας ο Χατζηλεοντιάδης ανακοίνωσε ότι η παράδοση των εργασιών δε θα γίνει μέχρι την ημέρα των εξετάσεων αλλά μετά το πέρας της εξεταστικής και συγκεκριμένα μέχρι την Παρασκευή , 29/2/2008

Σωστός ο Λεό! :D

Ξερει κανεις που παρακολουθουσε αν θα δοθει και 2η εργασια?

Είπαν ότι θα δοθεί.

Έχει βγει ήδη στο site της ψυχής.. ;)

Πάλι καλά που παραδίδουμε μέχρι τις 29/02, γιατί αλλιώς θα μας έβγαινε η ψυχή...  ;D

^hello^


έχει κανείς άποψη για τη δυσκολία?

     Τι να σου πω! Δεν είμαι και σίγουρος αλλά μου φαίνεται πως παλεύεται! Έχω σκεφτεί περίπου τί πρέπει να γίνει, κάτι, βέβαια, μου λέει πως στην πράξη θα κολλήσουμε!
      Πάντως η εκφώνηση είναι όλα τα λεφτά! Το ηχητικό σήμα λέει είναι προϊον υποκλοπής και πρέπει να βρούμε ποιος αριθμός καλείται! Και μόνο για αυτό αξίζει να την κάνει κανείς! Γενικά το σήμα είναι τελείως κατεστραμμένο από το θόρυβο. Αν τελικώς καταφέρω να βρώ τον αριθμό θα τον καλέσω κιόλας για να δώ ποιός θα απαντήσει, αν φυσικά είναι αριθμός τηλεφώνου!

Παιδιά please ανεβάστε κάποιος τις εργασίες (+αρχεία ήχου)..

Έχει πέσει ο server και δεν μπορώ να τις βρω..

 :)

Ευχαριστώ πολύ!

τελικά είναι μέχρι αυτήν ή την άλλη παρασκευή? Προλαβαίνω λέτε να ασχοληθώ κι εγώ?

την άλλη κορίνα την άλλη !!

29/2

Οκ, το πιασα...Στο κυριότερο δε μου απαντήσατε! Τελικά προλαβαίνω?????? :P :P

για την πρώτη νομίζω προλαβαίνεις...

χωρίς να έχω κάνει τίποτα...αλλά έτσι λεει η διαίσθηση μου...  :P :P

Χαλαρα προλαβαινεις, αμα δεν εχεις 4 μαθηματα να δωσεις αυτές τις μέρες  :P :P...

Μόνο ηλεκτρακουστική μου μένει...Αφου έστειλα και το Πεδίο 3...

πού το στείλες το καημένο το πεδίο 3 :P :P :P

στο διάολο...μετά το πεδίο ένα δεν έχω άλλες πεδιακές αντοχές.... :P

σε καταλαβαίνω... :P :P :P

είστε οφφ τοπικ. μπαχαλάκηδες!  :P


κανένας να μου απαντήσει για το μοναδιαίο: :D

http://www.thmmy.gr/smf/index.php?topic=17406.msg435930#msg435930

αφού σου απάντησαν ρε :P :P

Έπιασα να τις κάνω :D

Υπάρχει η συνάρτηση fft.Στο help έχει οδηγίες και παραμέτρους.
Μπορείς βέβαια να γράψεις δικιά σου συνάρτηση dft αφού ο αριθμός των δειγμάτων δεν είναι απαράδεκτα υψηλός αλλά λογικά με τα 1024 δείγματα μπορεί να κάνει κανένα λεπτό να τελειώσει.

Τελικά δε μου βγάζει σωστά το ποστ μου. Το έβαλα σε code, δες τώρα τι ρωτάω.

ΟΚ, με την υπόθεση ότι παίρνω Ν στοιχεία της ακολουθίας Χ (τώρα ξαναείδα τη θεωρία), θα πρέπει να αντιγράψω 32+64+128+256+1024 νούμερα στο Word :-X

Μπορείς να πάρεις κατευθείαν τον πίνακα που σου δημιουργεί το matlab με τον fft και είτε να τον εισάγεις με κάποιο τρόπο που δεν έχω ψάξει στο word,πάντως να γράψεις με το χέρι όλα αυτά τα νούμερα δεν νομίζω ότι ενδείκνυται..

Δεν ήξερα ότι έχει έτοιμο FFT και είχα κάνει δικό μου πρόγραμμα. Επειδή πιθανότατα είναι λάθος και θα ήταν λιγάκι άδικο να την πατήσω σε κάτι που έτσι κι αλλιώς το έχει έτοιμο :-[ ας του ρίξει μια ματιά όποιος θέλει.

Αν n είναι ο αριθμός των δειγμάτων, η x(k) υπολογίζεται από την

x=0;
for h=0:1:n-1
   x=x+2*exp(-h/128-2*pi*k*h/n*i)*cos(pi*h/2.56);
end;

Καλά το έχω;

    Δεν ξέρω αν είναι σωστό αυτό που κάνεις! Πού το βρήκες;
    Ο πιο απλός τρόπος για να το ελέγξεις είναι να δεις αν παράγει τα ίδια αποτελέσματα με τον FFT.

Αντικατέστησα στον τύπο ορισμού του DFT.

έχει έτοιμο κώδικα για FFT, ψάξε στο google

Ο ΦΦΤ κάνει επιπλεόν κανονικοποιήσεις οπότε θα είναι λίγο δύσκολη η σύκγριση.Ισως θα πρέπει να κανονικοποιήσεις και τους 2 πίνακες ως προς ένα κοινό μέγιστο.Τότε θα πρέπει να ταυτίζονται οι 2 πίνακές σου.
Η αλλιώς έλεγξε την σε κανένα παράδειγμα που μας έκανε αυτός οπότε έχεις τα αποτελέσματα ;)

:???:

Νομίζω μας είχε κάνει μια άσκηση με 5 δείγματα οπότε έχεις το αποτέλεσμα έτοιμο.Αν τρέξεις και τον αλγόριθμό σου θα μπορείς να συγκρίνεις!
Αλλά ακόμη και αν δεν την βρεις αυτή την άσκηση πάρε μόνη σου 5,6 δείγματα και κάνε dft στο χέρι και δες αν βγάζεις τα ίδια με το matlab!

Guys τώρα έχω πρόβλημα όχι με τον DFT αλλά με το σήμα... το plot βγαίνει σα να είναι όχι ένα σήμα, αλλά πολλά! Ας μου εξηγήσει κάποιος τι συμβαίνει :-[

Επίσης... ζητάει το σήμα που προκύπτει από δειγματοληψία του X(f) (του αναλογικού)... αλλά αυτό έχει και δέλτα μέσα (νομίζω ότι βγαίνει 10[δ(f-52)+δ(f+48)])... οπότε τι κάνουμε; :-[

Την κάνει κανένας άλλος; :(

Αν χρησιμοποιήσουμε τη συνάρτηση fft πως βγάζουμε τον πίνακα με τα αποτελέσματα; :D

Anonymous καλέ μου άνθρωπε! Λίγη βοήθεια στην άσχετη! :-[

Βασικά νομίζω ότι κάτι δεν κάνω καλά στον κανονικό FT, γιατί τότε στο τελευταίο ερώτημα θα είναι Y(f)=X(f) :-X

Τι εννοείς?
Να βγάλεις τον πίνακα και να τον πας πού?
Για να το βάλεις στο word ως έναν πίνακα 1024 νουμερών?
Ετσι θα το κάνεις εσύ?
Ασχετο,σίγουρα σε χαρτί τη θέλει ή μπορούμε και σε ηλεκτρονική μορφή?
Αν είναι έτσι του δίνουμε τον πίνακα από το ματλΆμπ.

Δεν ξέρω πως θα το κάνω. Ο πίνακας από το ματλαμπ πως βγαίνει;
(αν το κάνω σωστά :D ) :P

Όταν σου δίνει σήμα χρόνου και ζητάει τον DFT, στην ουσία δε σου ζητάει τους συντελεστές της ακολουθίας;

Άρα τον πίνακα... Σωστά;

Τελικά βρήκα πως βγαίνει ο πίνακας... δεν έπρεπε να βάλω ; στο fft(x).

Nessa τι πρόβλημα έχεις;
Χρησιμοποίησες την έτοιμη συνάρτηση του matlab;

Η σωστή ερώτηση είναι τι προβλήματα έχω ^bangwall^

1) Όταν τον βάζω να μου βγάλει σε γράφημα το σήμα στο Β μου βγάζει άλλα αντί άλλων, όπως βλέπεις στο σχήμα. Το σήμα που τον βάζω να απεικονίσει είναι το x(n/F).

2) Δε μπορώ να κάνω δειγματοληψία του FT γιατί (νομίζω ότι;) είναι αυτός που έγραψα παραπάνω με τα δέλτα.

3) Το Y(f) ως εκ τούτου βγαίνει ίδιο με το X(f).

Επίσης, δεν είμαι βέβαιη ότι ο DFT είναι αυτός που έγραψα παραπάνω.

Γενικά δηλαδή είναι θεωρητικά τα προβλήματά μου και όχι του Matlab... :-[

Το phone.mat από τη δεύτερη εργασία με τι ανοίγει; Γιατί το matlab δεν το ανοίγει :???:

Αρχείο matlab είναι σε μορφή σαν struct.Το κάνεις import από το μενού ή load στον κώδικά σου στο m-file.

Ευχαριστώ...

Βρε ikoufis, εσύ στο Γ τι έκανες με τη δειγματοληψία; Τι σήμα βγαίνει;

"What column separator does your data use?"

Και το FileInterpretError τι είναι;

Εγώ FFT έκανα και βγαίνει μια σαν δ αλλά πλατιάζει όσο κατεβαίνει,σε 2 συχνότητες,στα -50 και στα +50 Hz.
Εχω δηλαδή 2 "υψώματα" στα -50 και στα +50 Hz και πολύ χαμηλά στο ενδιάμεσο.

Για το άλλο,πας file import data βάζεις το phone.mat και πατάς finish.Το εισάγει στο πεδίο που σου βγάζει τις τρέχουσες μεταβλητές σου.

Λογικά έτσι θα ήταν, άρα σωστά το έχω και δεν έκανα πατάτα :) απλά για το άλλο που λέει ότι προκύπτει από δειγματοληψία του "κανονικού" FT...

Γιατί έχουν μέγεθος 0x0;

Βασικά εδώ τι πρέπει να κάνω; :???: :???: :???:

εγώ δεν έχω αυτό το μενού.Μάλλον λόγω άλλης έκδοσης ματλαμπ.

Δε μπορώ με καμία Παναγία να το κάνω να δει το σήμα! ^banghead^

Μήπως μπορεί κάποιος να μου το δώσει με άλλη μορφή; :-\

Ρε παιδιά κανένας; :o

Δε γίνεται να βάλετε τα νούμερα στο notepad και να μου τα στείλετε σκέτα;

Δοκίμασα να τα κάνω copy αλλά μου κόλλησε, είναι πολλά τα νούμερα... :o

Δοκίμασε να γράψεις στο Command Window την εντολή load Phone.mat, με την προϋπόθεση βέβαια ότι το Phone.mat βρίσκεται στο current directory, και λογικά στο workspace θα δημιουργηθεί μια δομή με το όνομα phone και αν κάνεις double-click πάνω της θα σου βγάλει δίπλα τα νούμερα...

??? Error using ==> load
Unable to read MAT file C:\MATLAB6p5\work\phone.mat

File may be corrupt.

Εν τω μεταξύ το έχω ήδη κατεβάσει δυο φορές, δε γίνεται να είναι κάθε φορά corrupt :(

Nessa μήπως να έβαζες μια πιο καινούρια έκδοση του ματλαμπ;

Ο DFT υπολογίζεται με τη fft του matlab. Ο πραγματικός μετασχηματισμός Fourier πως; :???:

Πού να τη βρω; :(

Με το χέρι θα το κάνεις. Στη δειγματοληψία κολλάω εγώ, γιατί έχει δέλτα μέσα.

Το αρχείο phone.mat κατάφερε να το ακούσει κανείς;

Νεσσα δοκιμασε να φορτώσεις το παρακάτω αρχείο phone στο ματλαμπ.

Είναι το phone.mat αποθηκευμένο για παλαιότερη έκδοση του ματλαμπ.

Επιτέλους!!!!! Το ανοίγει!

δεν εχεις δελτα...δε ειναι σκετο συνημιτονο εχεις κ ενα εκθετικο μπροστα....

Πως θα σχεδιάσουμε τα φίλτρα στην 2η εργ. όταν δεν μας δίνει συχνότητα και πλάτος ζώνης διέλευσης;

Αναλογικό σήμα,στους πίνακες με το Fourier έχει αντίστοιχη συνάρτηση.
Με την εντολή wavwite το κάνεις .wav και το ακούς.
Εχει τέτοιο αυτί που θα βρεις από εκεί τις συχνότητες? :D
Μάλλον τα διαλέγεις εσύ αυθαίρετα με την κρίση του μηχανικού που σου έχει αναπτύξει το τμήμα.Αν και βέβαια εσείς δεν έχετε ασύρματο τηλεπικοινωνία οπότε θα δυσκολευτείτε λίγο,ενώ εμείς :D :D

Προσωπικά έχω... αλλά κάτι μου λέει ότι αυτόν τον τρόπο δε θα μας τον πάρει σωστό ;D

Για επαλήθευση λέει όμως :P

Πες ένα πλάτος ενδεικτικά :)

Όταν κάνουμε το αναλογικό πρότυπο, τι τάξης φίλτρο θέλει; Ή πρέπει να το μαντέψουμε κι αυτό :(


Εννοώ εκεί που λέει δειγματοληψία του FT, στο Γ ερώτημα. Ακόμα είμαι κολλημένη εκεί :-\

-3 db στη συχνότητα αποκοπής και -20 db στη συχνότητα που είναι χαμηλή.Οπως τιας ασκήσεις που κάναμε.Αν τα ορίσεις αυτά και το εύρος ζώνης του bandstop φίλτρου σου τότε το n πρκύπτει από τον τύπο.

thanks



Δεν θέλω να το ακούσω εγώ... Ο Μάκης θέλει (που μπορεί)


Τρελή κρίση μηχανικού!
Και σκεφτόμουν να πάρω ασύρματο Ι  :P

Τα -3dB τα σκέφτηκα, το άλλο το παίρνουμε αυθαίρετα έτσι;
Καλή φάση.

copy paste δεν παίζει;


Οι εργασίες μπορούν να παραδοθούν και σε ηλεκτρονική μορφή στο e-mail του Χατζηλεοντιάδη.

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

Απλά γράφεις fft(x) (ή όποια είναι η συνάρτηση τελοσπάντων) χωρίς το ; στο τέλος (εγώ τα έκανα μισά μισά, γιατί δε χωρούσανε και τα 1024 στην οθόνη, έφευγαν επάνω) και σου τα βγάζει ans=τόσο ans=τόσο κλπ. Θα τα κάνεις αυτά copy-paste σε ένα txt αρχείο και μετά θα βγάλεις με κάποιον αυτόματο τρόπο όλα τα ans= για να το κάνεις copy-paste στο Word.

Κάτι δεν πάει καλά με τη φάση... πώς πρέπει να βγαίνει η φάση του μετασχηματισμένου;

Λοιπόν τα επιλέγεις όλα και κάνεις αντιγραφή επικόλληση
μετά αν έχεις όρεξη τα φτιάχνεις και ποιο όμορφα :D ακριβώς όπως τα λέει η Νεσσα.

Εγώ τουλάχιστον δεν έχω σκεφτεί άλλο τρόπο!


Η φάση είναι "unwrap(angle(y))"

Όπα, τι κάνει αυτό το unwrap;

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

Γειά και πάλι..Να σε ρωτήσω κατι επειδή δεν έχω ιδέα απο το matlab η πρώτη εργασία γίνεται στο χέρι η να μην το παλέψω καθόλου? :-\

Στο χέρι θα ήθελε χιλιάδες πράξεις. Χρειάζεσαι κάποια γλώσσα προγραμματισμού.

Στον διακριτό ΜΣ φουριέρ που ζητάει, πρέπει να βρεις 1024 νούμερα... με το χερί είναι λίγο ζόρικο :D

Α, κι επίσης, η δεύτερη γίνεται μόνο σε Matlab γιατί το σήμα δίνεται σε μορφή που ανοίγει με Matlab :-\

Προτείνετε κάτι για bandwidth στα φίλτρα...

Ακόμα, όταν λέμε -3db, για ποια συχνότητα λέμε; Σόρι αν είμαι εκτός τόπου και χρόνου, εγώ βλέπω τη σχέση 6.4.5 από το βιβλίο που λέει x db για ω συχνότητα... και μετά αυτό το μετατρέπουμε σε ζωνοφρακτικό με τους τύπους στη σελ. 156, σωστά;

Έχω κολλήσει στη δειγματοληψία...

Παίρνουμε το x(t) το κανουμε δειγματοληψία για Ν=..., βρίσκουμε τα δείγματα (που είναι αριθμοί) και μετά; :???:


ή παίρνουμε με το χέρι τον FT της x(t) και μετά;

Ε, FFT θα κάνεις (στο διάνυσμα που βρήκες). Ή το προγραμματάκι που έγραψα εγώ πιο πάνω.


Καταρχάς πες μου πόσο τον έβγαλες να επιβεβαιώσω αν τον έκανα σωστά...

Μετά κάνεις δειγματοληψία σ'αυτό που βρήκες.

Α) τελικά ναι το βρήκα.

β) δεν έβγαλα κάτι. Απλά έκανα φουριερ την x(t) όπως ήτανε.

Ναι ρε ανόνυμους :D αυτό σου λέω, βγαίνει ένας τύπος με δέλτα;

Λίγο βοήθεια στη δεύτερη... έχω φτιάξει ένα πρότυπο αναλογικό κατωδιαβατό φίλτρο παίρνοντας αυθαίρετα την τάξη (n=10), και θέλω να το μετατρέψω σε ζωνοφρακτικό με εύρος ζώνης 6Hz (πάλι αυθαίρετα) γύρω από την πρώτη ανεπιθύμητη συχνότητα. Παίρνω τον τύπο μετατροπής από τη σελίδα 156, παίρνω την H(s) του κατωδιαβατού, αντικαθιστώ όπου s εκείνο το κλάσμα και αυτό είναι όλο; Just checking...

Βασικά δε μπορώ να το κάνω... δεν ξέρω να χειρίζομαι αντικείμενα tf :( (κάποιος δεν πρόσεχε αρκετά στο εργαστήριο των ΣΑΕ :P)

Ρε παιδιά, κατάφερα να μειώσω το θόρυβο και ακούγονται 10 πληκτρολογήσεις αριθμών στα πρωτα 2 sec και μετά άλλες 7 από τα 5sec περίπου μέχρι το τέλος. Πόσα είναι τα ψηφία τελικά ξέρουμε?

Αφου βρουμε την H(s) πρεπει να βρουμε την αντίστοιχη H(z). Σωστά;
Το matlab έχει την συναρτηση bilinear και το αποτελεσμα είναι δυο πολυώνυμα
για τον αριθμητη και τον παρωνομαστή. Αυτα τα πολυώνυμα ειναι συντελεστες
σε θετικές ή αρνητικές δυναμεις του z; 
Επίσης την εξοδο την βρίσκουμε με αναδρομική σχέση  η υπάρχει
και καποιος άλλος τρόπός; 

Πες μου λίγο μιας και βρήκες την H(s), πώς έκανες το μετασχηματισμό σε ζωνοφρακτικό;

Εκανα ένα κατωδιαβατό και ένα ανωδιαβατό που σε σειρά δημιουργουν ενα
ζωνοδιαβατό. Το ζωνοφρακτικό γιατι το θες;

Για να αποκόψω την ανεπιθύμητη συχνότητα.

Τις 110Hz και 5 kHz εννοείς; Εγώ αφησα απλά να περνανε οι συχνότητες 600 - 1500 περίπου
 που αντιστοιχουν στο DTMF. Υποθετω οτι ειναι σωστό το αναλογικό φίλτρο.
 

Καλή ιδέα... απλά έλεγα όταν έφτιαξες το κατωδιαβατό, μετά πώς έβαλες τον τύπο μετασχηματισμού για να πάρεις το άλλο. Σίγουρα θα υπάρχει κάποιος κομψός τρόπος.

      Τον μετασχηματισμό z του σήματος που δίνεται πώς τον κάνατε βρε παιδιά; Δεν γίνεται με τίποτα; Η εντολή ztrans δεν λειτουγεί για ακολουθία αριθμών αλλά μόνο για αναλυτική μορφή!

Δε θα πάρουμε τον τύπο του διγραμμικού μετασχηματισμού από το βιβλίο;

Ναι αλλά τις τιμές του phone; Δεν υπάρχει εντολή για Z-transform των τιμών του phone.

Για τα φίλτρα θέλουμε z-transform, για το phone δε χρειάζεται.

Α ναι, για τη συνέλιξη :D πάρε τον τύπο ορισμού :???:

    Για τη συνέλιξη υπάρχει η εντολή conv. Το κατάφερε κανένας με συνέλιξη στο χρόνο; Έχω πολύ συγκεκριμένο λόγο που ρωτάω...
    Edit: Η εντολή αυτή θέλει υπερ-υπολογιστή για να τρέξει!

τελικά ξέρει κανένας πως παραδίδουμε?
ηλεκτρονικά,ή παμε εκει?δεν εχει βγαλει καμμια ανακοινωση...

Χρησιμοποιήσα τη συναρτηση [numd,dend] = bilinear(num,den,fs) του matlab και βρηκα τους συντελεστες της H(z).
Μετά τους χρησιμοποιήσα στον τύπο y(n)=b0x(n)+b1x(n-1)+... - a1y(n-1) - a2y(n-2) - ...
Απλά στον υπολογιστή μ καθυστερεί παρα πολύ να υπολογίσει για n 441000 δείγματα.

Είτε στο γραφείο του είτε με ε-μαιλ η παράδοση.


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

Γειαααα...λοιπόν σχετικά με τον υπολογισμό του πραγματικού μετασχηματισμού Fourier στην πρώτη άσκηση γιατί δεν δουλεύει η συνάρτηση fourier του Matlab? την βάζω και μου βγάζει "??? Undefined function or method 'fourier' for input
arguments of type 'double'."
Έχει κανείς ιδέα αν γίνεται να την κάνουμε να δουλέψει ?
Διαφορετικά για να τον κανουμε με το χέρι παίρνουμε τον τύπο στη σελίδα 73 του βιβλίου και βάζουμε όπου z=exp(jωΔτ)??Και όπου Δτ τι βάζουμε?το Ts=1/256???
Thanks in advance!

Για να κάνεις Fourier στο matlab η εντολή είναι fft(x)

Αν δε θέλεις να το κάνεις με τον έτοιμο FFT, πάρε τον τύπο στη σελ. 36. Ενδεικτικά έγραψα σε κάποια από τις προηγούμενες σελίδες μια συνάρτηση που υπολογίζει το DFT του συγκεκριμένου σήματος. Δεν είμαι 100% σίγουρη ότι δεν έχει λάθος, αλλά τα αποτελέσματα που βγάζει φαίνονται πολύ λογικά.

Η fft(x) όμως είναι για να βρούμε τον dft του σήματος.Αυτο που ρωτώ είναι για το δευτερο κομμάτι που ζητά τον υπολογισμό του πραγματικού μετασχηματισμού fourier.

Αυτόν, με το χέρι ^dontknow^

Για 1024 δέιγματα θα χρειαστούμε γύρω στις 10 σελίδες στο Word. Σίγουρα αυτό θέλει ? Έτσι κι αλλίώς τα δείγματα μπορεί  να τα δεί και από το mat file.

Δηλαδή να πάρουμε την εξίσωση 3.1.1 στη σελίδα 47 και να κάνουμε πράξεις? Και στα διαγράμματα, που βγάλαμε πριν με τον dft, πως θα βάλουμε τα αποτελέσματα? Αυτό που καταλάβαμε εμείς είναι να χρησιμοποιήσουμε τον τελευταίο τύπο της σελίδας 73 αντικαθιστωντας τα z και δτ όπως είπα και πριν, και να το κάνουμε αυτό μέσω matlab για τα διάφορα Ν που δίνονται , για να βγάλουμε τα διαγράμματα και να τα παρουσιάσουμε σε κοινό διάγραμμα με αυτά που βγάλαμε με την εντολή fft(x).

Τις κορυφές στον DFT σε ποιες συχνότητες τις βγάζετε;

στα 50 και στα 200 περίπου

Ωραία, κανονικοποίηση δε χρειάζεται έτσι;

αν πηγες με fft χρειαζεται να πολαπλασιασεις με ΔΤ

Στα πολύ αρχικά αν μπορεί κάποιος να βοηθήσει:

Στην 1η άσκηση, στο Β ερώτημα, μας ζητάει ουσιαστικά να σχεδιάσουμε το σχήμα γ) της σελίδας 49;

Στην 1η άσκηση, στο Α ερώτημα, έχετε κάνει plot στα DFT που προκύπτουν; Μου βγαίνουν περίεργα σχήματα. Στα attachments φαίνονται τα σχήματα μου για N=1024 της x(n) και της Χ(k).

το xn έτσι βγαίνει και σε μας, αλλα το χκ τι είναι?πως το έβγαλες?

ποιο; τη μπάντα με τις συχνότητες;

πολαπλασιαζεις το fft που βρήκες με το ΔΤ=1/256.. ετσι το κανονικοποιεις

ioannisnirvana είναι το plot που έκανα για το σήμα που έχει υποστεί fft. Δεν χρειάζεται κάπου αν έχω καταλάβει καλά αλλά το έκανα για να δω πως βγαίνει. Το x(n) που βγάζω είναι αυτό που ζητάει στο β ερώτημα;

ναι, αυτο που βγάζεις είναι για το β.

ξέρει κανείς τι είναι η μονότονη απόκριση συχνότητας? Τα φίλτρα στη δεύτερη άσκηση τα παίρνουμε χαμηλοπερατά?

αυτό στον διακριτό, στον πραγματικό ή και στους 2;

(thanks!)

τον πραγματικο τον αφήνεις όπως είναι..ουσιαστικά αυτόν ψάχνεις

Τα plots στον DFT και στον πραγματικό βγαίνουν τα ίδια περίπου; ή εκανα μαλακία; :D

Τι να σου πω, εγώ τη φάση στον DFT τη βγάζω γραμμική αύξουσα και στον FT μηδέν :(

Στο πλάτος;

Κοίτα, αν ο FT είναι όντως αυτός που έγραψα παραπάνω (?) και αν απεικονίσεις τα δέλτα σαν ψηφιακά δέλτα (?) τότε θα δεις ότι εκεί που έχει μονάδες στον ψηφιακό έχει κάτι σα δέλτα στον DFT και αλλού ο DFT πλησιάζει το 0 (ο FT είναι 0).

Μπουχουχου... δεν καταλαβαίνω τίποτα... όταν κάνω τον αντίστροφο μετασχηματισμό αυτού που είχα μετασχηματίσει δεν έχει καμία σχέση με το αρχικό... ουααααααααα :'( :'( :'( :'(

αν λες για αυτο που ζητάει στο γ) το χρησιμοποιεις γιατι ουσιαστικα θελει να δειγματοληπτήσεις το FT με Ν δειγματα καθε φορα

Στην πρώτη άσκηση στο γ ερώτημα, πόσο βγάλατε τα πλάτη του DFT και του FT; Εμάς τα peaks στα 50 Hz μας βγαίνουν με μεγάλη διαφορά μεταξύ τους. Ειδικά για N=1024.

Ευχαριστούμε...

Το e^-2t δε μετακινεί παρακάτω το peak του FT;

Για N=1024 στον FT βγάζουμε πλάτος περίπου 5 και στον DFT περίπου 1280. Δηλαδή περισσότερο στο DFT. Γιατί άραγε;

Και εγω το ιδιο προβλημα εχω. Μαλλον καποια κανονικοποιηση θελει ο FFT. Αν καταλαβα καλα το διβλιο σελ 54 λεει να πολλαπλασιαστει με ΔΤ διαστημα δειγματοληψιας αλλα παλι εχουν μεγαλη διαφορα τα 2 πλατη

Κανένας άλλος κάποια ιδέα; Γιατί έχουν τόση διαφορά τα πλάτη;

Παιδιά τι μορφή έχει ο μετασχηματισμός που βγάζετε; Γιατί εγώ μόνο δυο δέλτα βγάζω.

Ο μετασχηματισμός μας έχει για τον FT (από τον τύπο του βιβλίου του αναλογικού στην σελ. 118) ένα peak και για τον DFT (με την εντολή fft του matlab) 2 peaks.

Shuffle o πολλαπλασιασμός με το ΔΤ δίνει λογικά αποτελέσματα. Thanks.

Είναι λογικό ο DFT να έχει 2 peaks; Λογικά θα πρέπει να μοιάζει με τον FT οπότε να έχει ένα peak. Ξέρει κανείς γιατί συμβαίνει αυτό;

Τώρα πήγα να κάνω μετασχηματισμό z με τον τύπο του βιβλίου και symbolic object το z και κάνει πόση ώρα... και έχω και υπολογιστή τελευταίας τεχνολογίας :o

Ναι 2 peaks δινει

Στο τελευταιο ερωτημα της πρωτης εργασιας ο ifft σας δινει σωστο αποτελεσμα? Καποιο λαθος υποψιαζομαι με το πολλαπλασιασμο με τη u[f ] γιατι αλλιως δουλευει κανονικα.

Σε ποιο email στέλνουμε την εργασία;
Η παραδοση ειναι μεχρι σημερα η μπορουμε να τη στειλουμε και αυριο;

Ναι, δώστε το email... υπάρχει κάπου;

Δε βαριέσαι, θα τις στείλω ημιτελείς... :( :( :(

leontios at auth.gr

Δείτε και στο psyche.ee.auth.gr

Μέχρι σήμερα είναι η παράδοση.

εγω την εστειλα στη βικυ και στο μαιλ του μαθηματος...

Μόλις ενημερώθηκα ότι μπορούμε να παραδώσουμε τις εργασίες του μαθήματος και τώρα τον Σεπτέμβρη, την ημέρα και ώρα της εξέτασης. Για όποιον ενδιαφέρεται... ;)

Που είναι οι εργασίες?
το site?

http://psyche.ee.auth.gr/hcourses_fr.htm

Παιδιά να ρωτήσω τελικά.. Στην πρώτη άσκηση.. Στο α ερώτημα που λεει να βρουμε τον μετασχηματισμό Φουριε, στην ουσια χρησιμοποιούμε τη συναρτηση fft(x), μετά τι θέλει να βάλουμε στο word?? ολα αυτα τα νουμερα που προκύπτουν από τον πίνακα?? τον έχει ρωτήσει κάποιος μήπως?

Υποθέτω ότι σας ζητάει τα διαγράμματα μέτρου και φάσης από το Matlab.

Απ' ό,τι είδα στις προηγούμενες σελίδες όλα αυτά τα νούμερα θέλει.. :(
Tα διαγράμματα τα ζητάει παρακάτω..

Παιδιά τώρα πήρα χαμπάρι ότι μπορούμε να δώσουμε εργασίες και Σεπτέμβρη.Πόσο χρόνο περίπου χρειάζονται????

Η πρώτη πιστεύω σε 2 μέρες βγαίνει άνετα, για τη 2η δεν ξέρω, δεν την άγγιξα ακόμα.. :P

Στη 2η εργασία, αφού λέει ότι το σήμα είναι διάρκειας 10 sec, γιατί όταν το μετατρέπω σε wav είναι 55 sec? :???:

Άκυρο, δεν είχα βάλει τη σωστή συχνότητα δειγματοληψίας.. ::)