THMMY.gr

Παιδια καλησπερα, I'm back!!

Εχει ασχοληθει κανεις με τις εργασιες?Επειδη πηγα απ τη γραμματεια του τομεα και δεν ειχε συγγραμμα για το μαθημα,μπορει να μας πει καποιος αν βοηθανε αρκετα και τι παιζει?Υπαρχει καποιος που το εχει περασει και δεν χρειαζεται το βιβλιο να μου στειλει ενα post  ? :???: Γιατι καπου διαβασα οτι θα ξαναμοιραστουν απο τελη Απριλιου ενω η παραδοση της εργασιας ειναι 28/3 !!
Σημειωτεον εχω χασει μαθηματα και δεν εχω σημειωσεις δικες μου...
Αυτα τα site που εντοπισα σε αλλο topic απ το πανεπιστημιο της Πατρας για τα FIR IIR και γω δεν ξερω τι αλλο, αρκουν?

Ευχαριστω!!

Δίαβασε απο το βιβλίο των ΨΕΣ, τα λέει καλύτερα. Ούτως οι άλλως οι σημειώσεις για τα φίλτρα είναι χάλια

Ho Ho Ho, j j !!! Το βιβλιο των ΨΕΣ μου το λες για καλυτερο????!!!
Αυτο δεν ειναι οπως του αναλογικου και του στοχαστικου ? Φαντασου δηλαδη πως ειναι τα φιλτρα !!  :P
Το οτι δεν εχω περασει ΨΕΣ, μηπως, λεω μηπως, πειραζει λιγουλακι? εεεεεεε?  ^superinnocent^  (σφυριζω αμεριμνα)

   Κατά την γνώμη μου εάν κάποιος έχει ασχοληθεί με τα ΨΕΣ τότε δεν θα έχει πρόβλημα να βγάλει και της 4 νομίζω Εργασίες. Εμείς πάντως το 2007 καταφέραμε να πάρουμε της 2 μονάδες. Θέλει όμως δουλειά και χρόνο.Εκτός από το Βιβλίο των ΨΕΣ του Πανά, μπορείς να δεις τα ΨΕΣ του Monson Hayes (TZ), και το ξενόγλωσσο DSP των Proakis -Manolakis.
  Για το Matlab (αν το επιλέξεις) μπορείς να πειραματιστείς μόνη σου. Δεν έχει ιδιαίτερη δυσκολία.
Ένα Πρόβλημα που είχαμε πέρισυ είναι το θέμα με της απορίες..... Στο Μάθημα δεν προλαβαίνεις και δεν υπήρχε κ κάποια ώρα γραφείου να ρωτήσεις της απορίες σου :(. Αν λυθεί αυτό φέτος, όλα θα πάνε  καλά.

Ελπίζω να βοήθησα  8) 

Εγώ τις κάνω σε C γιατί είπε ότι θα έχει μια μονάδα παραπάνω... βασικά εκεί που κολλάω είναι άμα σου δίνει την ακολουθία με τη μορφή συνάρτησης τότε πώς θα βγάζεις τα στοιχεία της ακολουθίας για να τα βάλεις στο φίλτρο. Αυτό είναι math parsing και είναι πολύ δύσκολο... δε νομίζω να πρέπει να το κάνουμε έτσι!

Ό,τι και να κάνω... το αρχείο εξόδου μένει κενό... τι διάλο κάνει :(

Δε με θέλει σήμερα... ούτε η γραφική, ούτε τα φίλτρα... :(

Συνάδελφοι, πρέπει να κάνουμε στο Matlab (ή στη C κλπ) και τις 3 ασκήσεις της πρώτης εργασίας ή μόνο την τρίτη (που χρειάζεται ούτως ή άλλως) ?????

Λογικά και τις τρεις (αλλιώς εγώ γιατί παιδεύομαι; :D)

Τώρα πλάκα μου κάνετε?
Μόνο την τρίτη θέλει.
Ούτως η άλλως δε σου δίνει δεδομένα για να υλοποιήσεις τις 2 πρώτες

@svistos: Δες τι μου απάντησαν όταν ρώτησα σε άλλο thread... πρέπει να κάνεις ένα πρόγραμμα που να του δίνεις τη x και να σου βγάζει την y.

Απορία: Όταν βάζετε τη συνάρτηση που μας δίνει [x(n)=cos(0.01n)] στο πρώτο φίλτρο (H1 με IIR) σας βγαίνει άπειρο σε κάποια σημεία; Κανονικά δε θα έπρεπε, αφού είναι ευσταθές το φίλτρο. Μήπως μπορεί κάποιος να ανεβάσει ενδεικτικά αποτελέσματα για να επαληθεύσω αν έχω κάνει σωστά το πρόγραμμα;

Πέρσυ εμείς κάναμε σε ματλαπ μόνο το 3 κ μας έδωσε όλες τις μονάδες.
Τώρα αν δε βαριέσαι κάνε κ τα υπόλοιπα...

 ^dontknow^

Πάντως άκυρη η απορία μου πριν, το πρώτο φίλτρο δουλεύει.

Μένει τώρα να βρω τα bugs στο δεύτερο...

Δεν αντέχω άλλοοοοοο ^crybaby^

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

Απαντώ στον εαυτό μου: Δε διαβάζει τα δεδομένα. ^aytoxeir^

Γι'αυτό τη μισώ τη C... είναι τόσο εύκολο να την πατήσεις με μ^%$#ες και τρομερά δύσκολο να βρεις και να διορθώσεις αυτό που φταίει! ^tickedOff^

Θες μια συμβουλή ?   ::)

                    Μ Α Τ L A B          


                   ^wav^

Για το δεύτερο ερώτημα πως κάνουμε τη διαίρεση στο matlab για να υλοποιήσουμε το FIR; Γιατί με το χέρι βγαίνουν πολλές πράξεις!!!   ^dontknow^

Ευχαριστω  schelo!
Better late than ever!

Προς το παρων το παλευω χωρις τις σημειωσεις...(ακομα!!  ^tickedOff^ )

edit: never...

Να ρωτηξω και μια ερωτηξη...Στο πρωτο και δευτερο κομματι  θελει το σχεδιαγραμμα στο matlab για την αποκριση ή το κανουμε μονο στο τριτο κομματι της εργασιας?χεμ?

  Ορίστε η απάντηξης !!!  ;)

  Το 2007 κάναμε μόνο το 3ο και πήραμε όλες της μονάδες.    ^peace^
 
  Μην σε πάρω και στο λαιμό μου όμως......       ^kremala^

 

Ναι σωστα,μετα ειδα οτι κατι γραφτηκε στην προηγουμενη σελ,σορρυ..
Ευχαριστω παντως.Ξερω αλλους που εκαναν περυσι μονο το τριτο και αλλους που εκαναν σχεδιαγραμμα και για τα 3, γι αυτο.
Anyway.

Να ρωτήσω κάτι λίγο off topic.
Αν έχει παραδώσει κάποιος πέρυσι εργασία, αλλα δεν πέρασε το μάθημα, ισχύει και για φέτος  ή πρέπει να την ξανακάνει;  :)

Και εγώ σε αυτή την κατηγορία ανήκω...Ρώτησα τον Μουστάκα και είπε ότι πρέπει να τις ξανακάνουμε...

Στο 3ο ερώτημα με τα σχήματα τι παρατηρείται?  :(

Έλα ντε! Τι θέλει να απαντήσουμε; :???:

Πραγματικά όμως... βλέπετε κάτι που δεν βλέπω;;

"Τι διαφορές παρατηρείτε;"

Μισο ...Δηλαδη εσεις δεν παρατηρειτε διαφορες, βγαζετε ιδια αποκριση?Το δευτερο σημα δεν ειναι πιο γρηγορο?Αφου μεχρι να ολοκληρωσει μια κυμματωση το 1ο , το cos^2 το τρωει και το περναει μακραν....Βεβαια δεν ξερω αν θελει να σχολιασουμε κατι για το φιλτρο πιο συγκεκριμενα....

Ναι, φυσικά είναι διαφορετικές οι αποκρίσεις... απλά δεν καταλαβαίνω ακριβώς τι θέλει να σχολιάσουμε.

Τελικα @Nessa δεν κατεληξα σε κατι στανταρ...Απλα ό,τι εβλεπα το σχολιασα...Μια ερωτηση εχω....Οσον αφορα τις σημειωσεις, μηπως ξαναμοιραστηκαν?Πηρε κανενας προσφατα,δλδ  τη βδομαδα που μας περασε? Οι σημειωσεις ειναι οι ιδιες με τις περσυνες, αυτες του Πητα? Ακομα να τις παρω...Να ναι καλα ο @ Spacetraveler ...Και το πανεπιστημιο της Πατρας...

Εγώ τη Δευτέρα πήρα... δεν ξέρω πώς ήταν οι παλιές.

Ωωωω!!! Ελπιζω να χει μεινει τπτ αυτη τη Δευτερα...2η φορα που θα παω και να φυγω απραγη, δε λεει...Κοιτα τωρα που θα αγχωθουμε και για τις σημειωσεις !!

Νομίζω βγήκε η 2η σειρά ασκήσεων.
Αυτή χρειάζεται τα κεφάλαια 4-5 ή θέλει και παρακάτω τίποτα

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

Επίσης, πώς θα βρούμε το α του Kaizer, αφού δε μας δίνονται προδιαγραφές; Υποτίθεται ότι το α επηρεάζει την ακρίβεια του φίλτρου, σωστά;

Κάπου έχω χαθεί. Ποιο βιβλίο βοηθάει?
Από τα ΨΕΣ δεν θυμάμαι να έχει τρόπους σχεδίασης με υπολογιστή

Προσωπικά διαβάζω μόνο από τις σημειώσεις που κρατάω από τις παρακολουθήσεις και από το βιβλίο των ΨΕΣ. Το κακό είναι ότι πολλές φορές δε βγάζω τα γράμματά μου :-X Με αυτό το δεδομένο, νομίζω ότι αν πάρεις σημειώσεις από το Α4 θα είσαι καλά.

Παίρνω τα παράθυρα ως έχουν, αλλά να ξεκινούν από το μέγιστο που βρίσκεται στο 0 (μου φαίνεται και το πιο λογικό και σε ό,τι αφορά τη λειτουργία που επιτελούν). Αν ξέρει κάποιος κάτι διαφορετικό, ας το πει :)

Μπορεί να λύνεται αυτό με μια απλή ολίσθηση στο χρόνο μετά την εφαρμογή παραθύρου.
Άσχετο αλλά κάτι εισαγωγικό για το matlab υπάρχει χωρίς πολλά πολλά μόνο τα βασικά?

Εγω παλι δε καταλαβαινω τι κανω λαθος....βρισκω την hd(n) μεσω αντιστροφου Ζ, και στη συνεχεια φτιαχνω για καθε n το γινομενου του hd(n)*w(n)....οι τιμες που παιρνω ειναι κατι μιζερες μικρες....μετα κανω fft το αποτελεσμα και προσπαθω να το τυπωσω, αλλα τπτ.....

που το στριβω λαθος?....

Έτσι όπως μας τα λες, πουθενά! ^dontknow^

Εγώ πάλι θα ρωτήσω για το IIR... για να σιγουρευτώ: Από το πρότυπο Butterworth για να πάμε στο υψιπερατό που θέλουμε θέτουμε όπου s το ω_c/s ή το tan(ω_c/2)/s; Γιατί αυτό που γίνεται με τη συχνότητα αποκοπής δεν το έχω πολυκαταλάβει :-[

Ακολουθώ κατα γράμμα τις σημειώσεις, αλλά η h_D μου βγαίνει χαμηλοπερατό φίλτρο αντί για υψιπερατό, όπως φαίνεται από τη γραφική παράσταση της |H_D|.

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

Πριν λίγο βρήκα αυτήν τη σελίδα (http://www.phon.ucl.ac.uk/courses/spsci/dsp/filter.html) όπου εξηγείται σχετικά καλά η υλοποίηση φίλτρων και μάλιστα με παραδείγματα σε C (ωραία, έτυχα στην ίδια γλώσσα που χρησιμοποιώ :D αν και έχω ήδη φτιάξει το πρόγραμμα...).

Δείτε και αυτό το link για το πώς να μετατρέψεις ένα low-pass σε high-pass:
http://www.dspguide.com/ch14/5.htm

Τη γραφική παράσταση |HD| στο πεδίο συχνοτήτων πώς τη κάνετε; Με έτοιμη συνάρτηση ή αλλιώς;

Παιδιά στο highpass filter επειδή λογικά θα κάνουμε προσέγγιση θα πάρουμε π/4 σαν -3db ή στο π/4 θα είναι το τέλος της ζώνης μετάβασης δλδ 0db.
PS: Τα διαγράμματα σε db τα θέλει ή όχι, γιατί μόνο με db έχω βρεί, πως αλλάζει αυτό στο matlab?

Επειδή έχω μπερδευτεί λίγο με τα FIR ίσως να είναι λίγο χαζό αυτό που θα ρωτήσω. Αν έχουμε ένα low-pass που δίνεται από τον τύπο: sinc(n*omega)*omega/pi; τότε το φίλτρο είναι για n από 0 ως N-1 ή -(Ν-1) ως Ν-1 δηλαδή συμμετρικό; Και το παράθυρο hamming πώς πρέπει να πέσει αντίστοιχα; Εγώ δοκίμασα από 0 ως N-1 με μισό παράθυρο hamming όπως φαίνεται στο διάγραμμα. Σωστό είναι;

Λοιπον...Δεν ξερω απο που στο καλο να διαβασω....Βλεπω το βιβλιο ΨΕΣ, τις σημειωσεις του Πητα που πηρα απο ενα παιδι ( να ναι καλα!) που το χει περασει,και σημειωσεις απο αλλα πανεπιστημια...Υπαρχει ενα σημειο που να λεει καθαρα πως γινεται η υλοποιηση του highpass φιλτρου με IIR ? Δεν βγαζω ακρη...Μαλλον ξεκιναμε επ το κατωδιαβατο και μετα το μετατρεπουμε, αλλα που υπαρχει ενα παραδειγμα? Αυτες οι σημειωσεις που ειδα εγω, συνηθως δινουν προδιαγραφες, δλδ Α, fs fp κτλ.. αλλα στο συγκεκριμενο παραδειγμα που δε δινει  τι κανουμε?

ρίξε μια ματιά στο βιβλίο των Σκόρδα Αναστασόπουλου στην Αλεξάνδρεια μπορεί να σε βοηθήσει

Παιδιά το hamming παράθυρο βγαίνει κανονικά. Το ίδιο ακριβώς με το kaiser ττην κάνω τη μέθοδο μόνο ο τύπος του παραθύρου αλλάζει και ενώ kaiser έχω κανονικά διάγραμμα στο hamming έχω σωστή μορφή αλλά μέχρι 0.25 πηγαίνει max η |H(w)|!!!
Επίσης αν ξέρει κάποιος με τη freqz μου βγαίνει μέχρι π η γραφική παράσταση, είναι σωστό ή θέλει μέχρι 2π να φανεί όλο το πεδίο?

@alekos, δε νομίζω να το έχει το βιβλίο που λες (έψαξα και δεν το βρήκα) :???:

Πάντως έχει ένα αρκετά καλό του Mohamed Najim στα αγγλικά.

Για τη μετατροπή, τελικά ήταν απλό: H_dhpf=1-H_dlpf άρα με IDFT h_dhpf=δ-h_dlpf. Η συνάρτηση δ είναι φυσικά η διακριτή, δηλ. δ(0)=1 και τα υπόλοιπα 0.

Άρα h_dhpf(0)=1-h_dlpf(0) και h_dhpf=-h_dlpf για τα υπόλοιπα.

Εξακολουθώ να έχω αυτήν την απορία: :(

Το πρότυπο butterworth τι το θέλουμε? Για το iir με τη yulewalk βγαίνει κατευθείαν.
Και επειδή το help δεν βοηθάει είναι πολύ πιο καλό την yulewalk να τη δείτε από Help-->Product Help έχει ένα πολύ καλό παράδειγμα

Τι είναι το yulewalk;

Συναρτηση του Μatlab για να υλοποιησεις το IIR... Παραπανω δε λεω γιατι θα μπερδεψω κοσμο...Οταν το μελετησω καλυτερα θα ξαναχτυπησω!  :)

Για τη yulewalk δείτε αυτό το παράδειγμα που λέω παραπάνω είναι πολύ εύκολη συνάρτηση.
Στη γραπτή αναφορά τι γράφουμε τελικά? Έχει καταλάβει κανείς? Να ορίσουμε φίλτρο δηλαδή τι θέλει, ευθεία μορφή κτλ ?

Αν δεν προσεχεις τα παραπανω posts αυτα παθαινεις...επαναλαμβανεις κατι που εχει πει ηδη ο kekropas που τα ξερει καλυτερα!  :D

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

Μαλλον την υλοποιηση οπως στο Α μερος της εργασιας.... Μηπως?
Λεω μηπως ?

χμ θα κοιτάξω τι έλεγε η πρώτη άσκηση γιατί δεν την έχω κάνει.
Αλήθεια το hamming FIR σας βγαίνει κανονικά? Γιατί εμένα έχει max απόκριση |H(w)|=0.25 περίπου.

Ναι, κανονικά.

Ρε παιδιά όλοι εσείς που λέτε ότι δε σας βγαίνει δώστε τον κώδικά σας ή έστω κάποιες σχέσεις που βγάλατε για να δούμε πού κάνετε λάθος. Αλλιώς τι θα κάνουμε, θα μαντέψουμε;

 Κι εγώ το ΙIR φίλτρο το έκανα παίρνοντας Butterworth και στη συνέχεια έθεσα όπου s->ω/s. Ωστόσο στο σχήμα που υπολογίζεται στο matlab το φίλτρο βγαίνει όπως να'ναι... Το ότι θέτουμε s->ω/s το είδα από την σελίδα του μαθήματος ΨΕΣ που υπάρχει στο internet. Γιατί όμως τελικά το σχήμα βγαίνει λάθος; Δεν γίνεται με αυτόν τον τρόπο;  :-\ :-\

Δοκίμασε όπου s να βάλεις tan(ω/2)/s.

Γιατί; Από που το βρήκες αυτό; Εσένα δούλεψε έτσι;

Όχι, το δικό μου δε δουλεύει... συγκεκριμένα βγάζει στο |H(k)| τη γραφική παράσταση που επισυνάπτω. Ακόμα ψάχνω να βρω τι δεν έκανα καλά :-\

Ο τύπος που σου δίνω όμως είναι ένας που μας είχε πει στο μάθημα.

Έκανα την αντικατάσταση και κάπως έτσι βγήκε και σε μένα... Οπότε κάτι είναι λάθος...  :(

Εμενα παλι οταν βαζω το tan που αναφερατε μου βγαζει ενα χαζο-υψιπερατο.....οταν κανω αντικατασταση απο βαθυπερατο σε απορριψη ζωνης συμφωνα με τις σημειωσεις, εχω κατι καλυτερο....τωρα αν το γραφημα ειναι σωστο ιδεα δεν εχω...

Ποια τελική συνάρτηση βγάζεις; Εμένα μου βγαίνει:   y(n)=1.4330*y(n-1)-0.5615*y(n-2)+0.7486*x(n)-1.4972*x(n-1)+0.7486*x(n-2);

y(n)=1.45904y(n-1)-0.91038y(n-2)+0.19783y(n-3)+0.4459(x(n)-x(n-1))

Παίρνεις Butterworth 2ης τάξης ε;  Και μετά αφού κάνεις τις πράξεις χρησιμοποιείς την συνάρτηση bilinear του matlab για να βρεις την συνάρτηση με παράμετρο την z. Έτσι το έκανες;

Όχι βασικά εγώ το έκανα με τον κλασικό μαζοχιστικό τρόπο πράξεις-πράξεις-πράξεις :P

3ης τάξης φίλτρο πήρα.

 :D :D

Ε γι'αυτό προφανώς έχει διαφορά. Πάντως το σχήμα δεν είναι και πολύ διαφορετικό οπότε τζάμπα συζητάμε το λάθος... :P
Εσύ l0v3null με τί τρόπο το έκανες;

Εγω το εκανα με τη μεθοδο "βαριεμαι να τα μεταφερω στο πεδιο της συχνοτητας και τυπωνω το γραφημα στον s χωρις αντικατασταση", κι αφου η το γραφημα δε με χαλασε πολυ, το αφησα ετσι....θα μου πειτε για να βρω εξισωση διαφορων κλπ πρεπει να κουνησω και λιγο το δαχτυλακι μου,....μαλλον διγραμμικο το κοβω και καμια πραξη.....

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

Οι εργασίες είναι φέτος υποχρεωτικές ή προαιρετικές? Τί ισχύει πάνω κάτω?

Θέλω να ξέρω επειδή τώρα είμαι Βερολίνο και θα το δώσω Σεπτέμβρη... :D

Ευχαριστώ

Προαιρετικές είναι, μην ανησυχείς. :)

Dankeschön!

Αυτό που μου τη σπάει με τις εργασίες στα Φίλτρα είναι ότι αφενός δεν έχουμε καθόλου feedback ώστε να ξέρουμε τι δεν κάναμε καλά, αφετέρου οι προηγούμενες εργασίες χρειάζονται για τις επόμενες.

Υπάρχει κάποιος που να έχει κάνει σίγουρα σωστά τη δεύτερη εργασία, για να μας πει τι έκανε; Όχι μόνο για να δούμε αν την κάναμε σωστά, αλλά και επειδή χρειάζεται για αυτήν που κάνουμε τώρα :(

Αν τα σχήματα σου είχαν βγεί λογικά σωστά θα είναι. Εξάλλου τώρα φαντάζομαι μπορούμε να πάρουμε και έτοιμες συναρτήσεις για φίλτρα.
Κάτι άλλο όταν έχουμε 44100 δείγματα της χ(n) και κάνουμε fft έχουμε πάλι 44100 δείγματα της χ(κ) τώρα. Αυτά έχουν κάποια περιοδικότητα? Την fs ή την fs/2 και γιατί?

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

Δοκίμασε και με τη fir1 νομίζω να πάρεις απευθείας το φίλτρο. Φαντάζομαι λες για το 2ο μέρος, δεν έχω φτάσει ακόμα εκεί και για αυτό δεν ξέρω αν όντως υπάρχει κάποιο πρόβλημα γενικά.
Για το άλλο θέμα κανείς? Και κάτι ακόμα τα φάσματα θα είναι από -fs/2 ως fs/2 ή από 0 ως fs/2?

Λογικά είναι από 0 μέχρι και 2π.

Nessa το έβγαλες αυτό που έλεγες? Και εμένα ένας πανικός γινότανε κάπου το έχανα στην κανονικοποίηση, πήρα συνέλιξη στο χρόνο για να εφαρμόσω το φίλτρο και όλα οκ. Η συνέλιξη γίνεται με την w=conv(u,v)

Βασικά το έκανα λάθος... έπαιρνα το μισό φίλτρο! Κανονικά ήθελε shift κατά Ν/2.

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

Για το παράθυρο λες? Υπάρχουν και τύποι από 0 εως Ν και υπάρχουν και τύποι από -Ν/2 εως N/2 για τον ίδιο τύπο παραθυρου. Ανάλογα ποιον θα χρησιμοποιήσεις είναι. Εγώ σε αυτή την εργασία ήθελα να αποφύγω πάση θυσία το shift και δεν έκανα με τους τύπους από 0-Ν.

Αυτό που δεν έχω καταλάβει ακόμα είναι αν το παράθυρο είναι Ν=50 τότε το μη ιδανικό φίλτρο μας θα είναι ακολουθία h(n) με 50 όρους. Το σήμα εισόδου μας είναι το wav αρχείο ήχου μια ακολουθία μετά τη δειματολψία με παράδειγμα 44100 όρους. Πως θα γίνει στο πεδίο της συχνότητας ο πολλαπλασιασμός h(k)x(k) ?

edit:Καλά η φωνή του τύπου είναι σε όλο το φάσμα!

Αν κάνεις συνέλιξη στο χρόνο, δε χρειάζεται να κάνεις και πολ/μό στη συχνότητα, το ένα από τα δύο θα κάνεις... εκτός αν δεν κατάλαβα τι εννοείς :???:

BTW αυτό με το shift μας το είχε πει τελικά αλλά εγώ δεν το είχα καταλάβει. Κάπου μας είχε δώσει τον τύπο H'(z)=z^-(N-1)/2H(z) (στα βήματα σχεδίασης φίλτρου Kaiser) αλλά δε μού 'κοψε ότι είναι μετατόπιση στο χρόνο ^tomato^

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

Πριν κάνεις Fourier στην h, συμπλήρωσέ την με μηδενικά. Νομίζω έτσι γίνεται.

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

Καλα τι παιζει με την υποβολη....παω και δεν εμφανιζει τπτ.....τα ιδια και σε εσας?...

Υ.Γ Αντε εχουμε και τελικο αυριο.....

......ΕΙΣΑΙ ΑΡΡΩΣΤΕΙΑ ΜΟΥ ΜΕΓΑΛΗ,...ΠΟΥΘΕΝΑ ΔΕΝ ΕΧΕΙ,...ΘΑ ΤΡΕΛΑΘΩ ΑΝ ΔΕΝ ΣΕ ΔΩ ΘΑ ΤΡΕΛΑΘΩΩΩΩ.....Σ'ΑΓΑΠΩ ΚΑΙ ΘΑ ΣΤΟ ΛΕΩ.....ΑΡΗ ΟΛΕΕΕΕ ΟΛΕ ΟΛΕΟοοο.....ΔΥΝΑΤΑ ΠΙΟ ΔΥΝΑΤΑ, ΠΙΟ ΔΥΝΑΤΑΑΑΑ......

Ελπίζω να πρόσεξες ότι η προθεσμία αυτήν τη φορά λήγει μεσάνυχτα Κυριακής και όχι Παρασκευής όπως άλλες φορές. Μη στεναχωριέσαι τσάμπα δλδ :)

Εσείς μου τρομάζετε το σκύλο κάθε φορά ρε; >:(

Μωρε ας το σηκωσουμε σημερα και δε τπτ δε με χαλαει.....ελα κι εσυ Αρειανακι εισαι,....αστο ελευθερο μη το καταπιεζεις και θα νιωσεις καλυτερα....:D

Υ.Γ Κυριακη βραδυ ε?....καταλαβα,...θα εχει προβλημα μεχρι τις 11:58,...και θα το ανοιξει για 2'.....

Υ.Γ οοοΟΟΟΟΟΟοοοοο ΑΡΗ ΜΟΥ ΣΕ ΑΓΑΠΩΩΩΩ....ΦΕΡΕ ΜΑΣ ΤΟ ΚΥΠΕΛΟΟΟΟοοοο.....ΟΟΟΟΟΟοοοοΟΟΟΟ!!!!!!!!!

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

Εμείς τη στείλαμε με email καλού κακού...

Σειρά 4, στο Β, πώς βγαίνει η h; Έχει καταφέρει κανείς τον αντίστροφο μετασχηματισμό; :???:

Να ρωτήσω κάτι γενικό για τις εργασίες.Υπάρχει πρόβλημα αν δουλεύουμε με το toolbox του Matlab για τις διάφορες γραφικές παραστάσεις(αποκρίσεις συχνότητας κτλ);

Στο Α3 τις οριακές συνθήκες πως τις παίρνουμε? Για m<0 ή n<0 όλα 0 ή θέλει ψάξιμο?

Μας είχε πει στο μάθημα ότι όταν δε δίνονται αρχικές συνθήκες θα τις παίρνουμε 0.

Απλά το λέω γιατί έχει ένα καφάλαιο στο βιβλίο που μηδενικές αρχικές συνθήκες δεν παίρνει για n1,n2<0 αλλά κάνει κάτι με τα πεδία υποστήριξης των h και x υπολογίζει το αντίστοιχο πεδίο της y με συνέλιξη κτλ.
Αυτή τη διαδικασία πότε την κάνουμε?

Στις εξετάσεις :P

+1
Τι εννοεί υλοποίηση? Τι ζητάει να κάνουμε? Την h(n1,n2) θέλει?

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

Hint: Βρες πρώτα την εξίσωση διαφορών. Θα παρατηρήσεις κάτι που απλοποιεί τα πράγματα ;)

Δεν το βλέπω αυτό που λες νομίζω ότι βρήκα τρόπο για τον αντίστροφο μετασχηματισμό Ζ όμως

Να ρωτησω επειδη θα φαω τα νιατα μου, υπαρχει συναρτηση στο Matlab αντιστοιχη με το difilt για δισδιαστατα φιλτρα? Για την εφαρμογη φιλτρου σε εισοδο υπαρχει η filter2 αλλα δε βρικσω σε τι φιλτρο την εφαρμοζει.