[Ψηφιακά Φίλτρα]Εργασίες 2011

[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 με προσοχή στα ονόματα των αρχείων!

[Admins]

Για αλλαγή του public name σας, επικοινωνήστε με έναν από τους Admins.

[Time_Tripper]

με το λ τι ακριβώς παίζει και τους περιορισμούς του?

θετικο.

Μα... ο κώδικας της fftfilt χρησιμοποιεί/καλεί τη συνάρτηση fft...
το x[n-3k] δεν υφίσταται. x[n-3] είναι. απλά να ξες ότι πρέπει να μετατοπίσεις τα b κατα ένα δεξιότερα.
π.χ. το b[4] θα έχει τον συντελεστή του x[n-3], το b[7] θα έχει τον συντελεστή του x[n-6] και αυτό γιατί το b[1] έχει το συντελεστή του x[n] και όχι το b[0] (b[0] δεν υπάρχει στη matlab).

τα b ειναι παντου τα ιδια(σ' ολα τα y)  οποτε απλως κανεις εναν πινακα b με τα στοιχεια του στη σειρα και δεν εχεις να ανησυχεις για τιποτα...
Το x απ την αλλη αν γραψεις απλως x[n-3] τοτε το fftfilt θα παρει τα x[n-3], x[n-4] κ.ο.κ που δεν το θες.
Εγω αυτο που εκανα ειναι εναν tempx ο οποιος καθε φορα κρατα τις 8 απαραιτητες τιμες.
Αυτο που με δυσκολεψε πιο πολυ γιατι αργησα να το καταλαβω ειναι οτι η  fftfilt δεν επιστρεφει την y[n] αλλα επιστρεφει σε πινακα τα 8 "στοιχεια" του (μην την πατησει κανενας αλλος)...
Τωρα βεβαια γιατι μπαινουμε στη διαδικασια να κανουμε fft για 8 μονο στοιχεια δεν μπορω να καταλαβω...

edit: ακυρο, δεν επιστρεφει αυτο η fftfilt ρε πατριωτες τι στο δι*ολο επιστρεφει αυτη η ^&%συναρτηση γιατι ακομα δεν εχω καταλαβει?

[Faros]

Η αλλη εκδοχη ειναι να το φτιαξουμε μονοι μας αλλα να κανω ενα function που να μετατρεπει fft την εισοδο, να κανει πολλαπλασιασμο με H και να παιρνει μετα τον αντιστροφο fft για την εξοδο....εεεεε ΔΕ ΣΦΑΞΑΝΕ!!!!!

Αυτο βασικα,πέρνεις fft για x[n] και για h[n],βρίσκεις τη Y(z) και μετα τη ζητούμενη.Η συναρτηση που θα προκύψει βέβαια δε θα μοιάζει και πολυ με την αντίστοιχη του 3ου ερωτήματος μιας και χρησημοποιούμε μονο 8 όρους,άμα βαλεις παραπάνω τοτε θα πάρει την ίδια μορφη.

P.S για να ισχύει ο τύπος Y(z)=X(z)*H(z) πρέπει να πάρεις κυκλική συνέλιξη

[Time_Tripper]

[
Αυτο βασικα,πέρνεις fft για x[n] και για h[n],βρίσκεις τη Y(z) και μετα τη ζητούμενη.Η συναρτηση που θα προκύψει βέβαια δε θα μοιάζει και πολυ με την αντίστοιχη του 3ου ερωτήματος μιας και χρησημοποιούμε μονο 8 όρους,άμα βαλεις παραπάνω τοτε θα πάρει την ίδια μορφη.

P.S για να ισχύει ο τύπος Y(z)=X(z)*H(z) πρέπει να πάρεις κυκλική συνέλιξη

  πως παιρνεις την fft της h[n]? Και πως αυτο ειναι μετα FIR φιλτρο? Αφου το  FIR εχει συγκεκριμενη μορφη. Δεν πρεπει καπου να χρησμιποποιησεις τον τυπο απο το 2?

[Faros]

Για αυτό δεν είμαι σίγουρος,εγώ πήρα την εντολή  f=fft(x)  κατευθείαν ,ίσως να χρειάζεται να γράψεις ολον τον κώδικα.Όπου h[n] βάζεις οπως λες τους ορούς που βρήκες απο το 2.

[AgentCain]

Η αλλη εκδοχη ειναι να το φτιαξουμε μονοι μας αλλα να κανω ενα function που να μετατρεπει fft την εισοδο, να κανει πολλαπλασιασμο με H και να παιρνει μετα τον αντιστροφο fft για την εξοδο....εεεεε ΔΕ ΣΦΑΞΑΝΕ!!!!!

Έτσι κάντο
Εγώ αυτό έβγαλα, μοιάζει κάπως στο IIR φίλτρο, τουλαχιστον για 8 όρους

Τελικά η χαζομάρα που έβγαζα δεν ήταν χαζομάρα, εγώ είχα κανει λάθος στο plot 

edit: στο παρακάτω έχω κάνει επίσης την εξής χαζομάρα, έβαλα όρους FIR μέχρι το z^-21 και όχι μεχρι το z^-6

[edenaxas]

3rd, is up.

[vag]

το 44100 δεν είναι ακέραια δύναμη του 2
πως συνεχίζουμε για το Decimation in Time ?

[edenaxas]

zero padding vag

[vag]

cool

thanks

petyxa k afto alla sketi zalada m fainetai
http://www.ece.cmu.edu/~ee791/lectures/L18/GeneralRadixFFT.pdf

[edenaxas]

http://cnx.org/content/m12016/latest/

kai yparxei kai afto,

http://www.mathworks.com/matlabcentral/fileexchange/17778

[vag]

αν καλέσω την fft του matlab με όρισμα το zero padded σήμα βρίσκουμε τα ίδια
αν του δώσω το original σήμα καμία σχέση
σχηματικά το αποτέλεσμα είναι ίδιο αλλά συμπιεσμένο στο διάστημα 44100 και όχι στο 65536 που είναι το next power

άρα
ή
η fft του matlab δεν δουλεύει με zero padding
ή
εγώ κάνω κάτι λάθος..

ANY IDEAS?

επίσης
wtf τον παίρνει λίγο η εργασία ή μου φαίνεται?
brainfucked

[arashi]

Δεν σβηνονται με περιεργο τροπο

το τοπικ καθαριζεται

[jimmakos]

Να ρωτήσω και γω.
Για το τελευταίο ερώτημα που λέει να σχεδιάσουμε ένα ιδανικό χαμηλοπερατό φίλτρο.
η απόκριση ενός ιδανικού χαμηλοπερατού είναι h(t)=2*Β*sinc(2*B*t) ή H(f)=rect(f/2B)
Πρώτον. μαντεύω ότι πρέπει να βάλουμε μια συχνότητα στο Β γύρω στο 100-200 διαιρούμενη με το (22050/2). Γιατί;
(Hint, το είδα από εδώ: http://www.aquaphoenix.com/lecture/matlab10/page4.html στο 10.3.2.1)

Δεύτερον. Όντως θα χρησιμοποιήσουμε αυτή την εξίσωση;

Τρίτον. Αν στο 2ο η απάντηση είναι ναι, τότε για πόσο χρόνο (t) θα την κατασκευάζαμε αν δουλεύαμε με συνέλιξη;

Τελευταίο. Είναι λογικό κάτι τέτοιο στο διάγραμμα του φιλτραρισμένου(ντεμέκ) αποσπάσματος;

[crystal]

Βγήκε και η 4η..

Τι είναι το βαθυπερατό φίλτρο? 

[AgentCain]

Κακή μετάφραση του χαμηλοπερατού?