THMMY.gr

Παιδιά επειδή έλειπα στα πρώτα μαθήματα και επειδή ο Στρίντζης συνήθως μοιράζει ένα φυλλάδιο που λέει τι θα έχει το μάθημα
και τις προθεσμίες των ασκήσεων,μήπως μπορεί κανείς να μου πει πότε είναι η προθεσμία της πρώτης άσκησης ή ακόμα καλύτερα να σκανάρει το φυλλάδιο?

Στις 6/3 έλεγε (Καθαρά Δευτέρα). Ο μεταπτυχιακός μας είπε ότι αυτή ήταν μία ενδεικτική ημερομηνία, άρα θα δούμε...

Μόλις ήρθε εμαιλ ότι μπήκε η πρώτη εργασία στα φίλτρα...βρίσκεται στο εΤΗΜΜΥ...

και ερωτώ ...πως θα την κάνουμε όταν δεν έχουμε καν τα βιβλία..???οεο??? :D :-\

 :-[ :-[
κιολας??

εχει παραδοση μεχρι 13/3 και μαλλον θα ειναι κατι απλο οποτε βλεπουμε...!

Μου είπανε ότι στο μάθημα της Πέμπτης ο μεταπτυχιακός που κάνει τις ασκήσεις είπε ότι τελικά η πρώτη εργασία μπορεί να παραδοθεί μέχρι την Παρασκευή αντι για τη Δευτέρα.
Ισχύει αυτό;  ^que^  Γιατι μπήκα στο eTHMMY και δεν έχει αλλάξει την ημερομηνία μέχρι στιγμής

Ναι, είπε ότι η διορία θα είναι μέχρι την Παρασκευή 17/3.

Παιδια με τα FIR τι κανατε? βρηκατε την h[n]? εμεις το καναμε ετσι αλλα βγαινει τεραστιο... πηραμε αντιστροφο Ζ του πρωτου κλασματος αντιστροφο του δευτερου και ειπαμε οτι θελουν συνελιξη αλλα μου φαινεται too much :)

Εγώ βρήκα την h(n) στην 1η με τη μέθοδο των υπολοίπων και στη 2η με επέκταση σε δυναμοσειρά.
στην 1η μου βγήκε λίγο τεραστιο (5 υπόλοιπα) αλλα απο την συνέλιξη νομίζω ειναι ευκολότερο.

Στην τρίτη άσκηση που λεει να κανουμε την υλοποίηση χρησιμοποιώντας τις υλοποιήσεις των H1 & H2
τι εννοει?
Να μην πολ/σουμε τις H1(z) H2(z) ?

Η Η3=Η1*Η2. Αρα ειναι τα δυο φιλτρα που υπολογισες σε σειρα!

thx

Στην υλοποίηση με τον υπολογιστή τι ακριβώς πρέπει να κάνουμε? Να φτιάξουμε πρόγραμμα με βάση τις IIR υλοποιήσεις που θα του δίνουμε n και μας δίνει την έξοδο?

με την μέθοδο των υπολοίπων δεν είναι λίγο περίπλοκο εως πολύ περίπλοκο...άρχισα να το κάνω και από το πρώτο υπόλοιπο κιόλας απελπίστηκα από τις πράξεις που πρέπει να κάνω... Οι πόλοι που βγαίνουν δεν είναι οι:
1.τρίτη ρίζα του μ (πολλαπλότητας 3)
2.sqrt(λ) και
3.-sqrt(λ)
???
Όταν λες 5 υπόλοιπα ότι σου προέκυψαν....?

υπαρχουν και μιγαδικοι πολοι με παραμετρο το μ!

Για το FIR, το πρώτο που έχει και παραμέτρους, με τη μέθοδο της συνεχούς διαίρεσης, δεν λύνεται;

Εμείς έτσι κάναμε και το κάναμε μέχρι να φτάσουμε στο z^19, είναι άκυρο αυτό; γιατί να βρούμε πόλους με τις παραμέτρους μας φάνηκε κουλό.

Ναι γίνεται και με διαίρεση αλλά πρέπει να είναι ατελειίωτη...
Γίνεται ή τουλαχ έτσι ελπίζω και με επέκταση σε δυναμοσειρά και βγαίνει αμέσως αλλά αυτό σχετικό είναι - πάλι 20 όρους θες..

Mε την δυναμοσειρά μέχρι το h(5) δεν βρίσκουμε?

πώς βρίσκουμε τα επόμενα?

ή δεν καταλαβαμε κάτι.....

με την δυναμοσειρα παλι μπορεις να βγαλεις απειρους ορους για το πρωτο κλασμα απειρους για το δευτερο και να τους συνελιξεις.
εμεις κρατησαμε 5 και 5!

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

Έστειλα e-mail στο μεταπτυχιακό για το τι θα πρέπει να περιέχει το συμπιεσμένο αρχείο που θα στείλουμε και μου είπε:
Θα περιέχει ένα .doc με την αναφορά και τον κώδικα σε .m ή .cpp μορφή.

Ποιό πρόγραμμα βγάζει έξοδο σε αυτό το format; Matlab και C++;

Και αν το κάνουμε σε Java θα του στείλουμε .java αρχείο;

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

Ας το ελέγξει όποιος το έχει κάνει για διασταύρωση...

ΑΑΑΑΑΑΑΑΑΑΑΑΑΑΑΑΑΑΑΑΑ, με τις δυναμοσειρές απλά πέρνεις για Pi = Qi = 0 για ν>5 και είσαι τζετ?


κατάρα.. προφανώς... εεε? δεν έχω πρόχειρα και τα νούμερα μου.

Μέγκαβαττ αν το δεις.... κανόνισε.

ohh yes. χαρηκα τωρα!  ;D
εχουμε περιπου τα ιδια...απλα εμεις πηραμε λιγοτερους ορους και εχουμε χειροτερη προσεγγιση. κατι ορους δλδ λ^6 κλπ δεν τους εχουμε αλλα γενικα ειμαστε στα ιδια!

Το είδα... και πιάνω αμέσως δουλειά!  ;)

faidon σου ξέφυγε ενα λ στο h(12) και πηρε η μπάλα και τα παρακάτω συνεπώς.....

είναι h(12)=λ^6+....κοιταξε το.. ;)
τα πρώτα είναι όλα οκ

Για ποιο λογο να κανουμε τη συνελιξη; Θεωρουμε πως εχουμε 2 φιλτρα το ενα με τα μ το αλλο με τα λ, υπολογιζουμε τα FIR για το καθενα και το τελικο φιλτρο ειναι η σειρα αυτων των δυο. Χρειαζεται να υπολογιστει και τυπος;

Επισυνάπτω τα αποτελέσματα της επέκτασης σε δυναμοσειρά τουπρώτου
φίλτρου αφού τα διόρθωσα (ευχαριστώ simone) και τα αποτελέσματα επέκτασης
σε δυναμοσειρά του δεύτερου φίλτρου.
Τα δεύτερα θα είναι λογικά σωστά αφού τα έκανα με java...

Ας το ελέγξει όποιος το έχει κάνει για διασταύρωση...

Στο δεύτερο φίλτρο η IIR υλοποίηση θα γίνει κλασικά με την εξίσωση διαφορών
παρόλο που τα μηδενικά είναι περισσότερα από τους πόλους;



Έτσι όπως είναι η εκφώνιση μου φαίνεται ότι ζητά μία συνάρτηση μεταφοράς που
θα προκύπτει από ένα άθροισμα των h(i)x(n-i) και όχι υλοποίηση του συνδιασμού
σε σειρά.

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

Κάτσε ρε Jason γιατί αυτό που λες αν και ακούγεται λογικό με μπέρδεψε...
Δηλαδή στο 3το ερώτημα το τελικό φίλτρο (σε σχηματική μορφη με τα z^-1)
δε θα είναι ο σε σειρά συνδυασμός των δύο προηγούμενων IIR φίλτρων που
υλοποιήσαμε;

Ναι!
Αλλα αλλο αυτο...εγω λεω τοση ωρα για το FIR φιλτρο που μας δυσκολεψε αποτι φαινεται περισσοτερο.

Ναι, για τα FIR:
Αν θεωρήσεις ότι είναι γινόμενο δύο συναρτήσεων μεταφοράς μία με λ και μία με μ.
Στη συνέχεια υπολογίζεις το FIR της κάθε συνάρτησης μεταφοράς ξεχωριστά.
Tελικά υλοποιείς το φίλτρο με τον σε σειρά συνδυασμό των δύο FIR.
Θα καταλήξεις δηλαδή σε μία υλοποίηση όπως του σχήματος 5.8.2 μόνο που θα έχεις
δύο πολύ μεγαλύτερα φίλτρα και όχι πολλά μικρουτσικα.
Σύμφωνα με αυτό θα μπορούσες να χωρίσεις το φίλτρο και σε μικρότερους παράγοντες
και να πάρεις τον σε σειρά συνδυασμό των επιμέρους υλοποιήσεων.
Σωστά δεν τα λέω;

Απλά εγώ θεώρησα ότι στα δύο πρώτα ερωτήματα οι υλοποιήσεις έπρεπε να γίνουν με
ένα μόνο φίλτρο.

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

Η παράδοση θα γίνει στο eTHMMY αφού δηλωθείς στις εργασίες.

το ethmmy ειναι Down!!!!

Ναι φυσικα αμα καταφερες να το βγαλεις ετσι μια χαρα μου ακουγεται...
που βρηκες ομως αλλους παραγοντες να χωρισεις?

Μπορείς να το χωρίσεις σε τέσσερις παράγοντες ως εξής:

1/(1-λz^-2)(1-μz^-3) =

= 1/(1-sqrt(λ)z^-1)(1+sqrt(λ)z^-1)(1-cubert(μ)z^-1)(1+cubert(μ)z^-1+cubert(μ²)z^-2)

Ο τελευταίος χωρίζεται και σε ακόμα δύο αλλά προκύπτουν μιγαδικοί συντελεστές.
Όπως και να έχει νομίζω ότι για τα δύο πρώτα ερωτήματα οι υλοποιήσεις πρέπει να γίνουν με ένα φίλτρο.

Μπορεί κάποιος να με βοηθήσει με το τελευταίο ερώτημα που πρέπει να βρούμε την έξοδο
για x[n]=cos(αn); Δεν έχω ιδέα τι πρέπει να κάνω...

Εμείς πήραμε την εξίσωση διαφορων και υπολογίσαμε τις 100 τιμες της εξοδου θεωρωντας y[n]=0,n<0.
ορίσαμε 2 πίνακες για να καταχωρουμε τα y[n-i] και τα x[n-i].και μετά τα βαλαμε στην εξίσωση διαφορων και με ένα loop υπολόγιζε τις 100 τιμες.θέλει και γραφική παρασταση.

λέει οτι το συστημα είναι αιτιατό.αυτο δε σημαινει οτι y[n]=0 , n<0?

Έχεις απόλυτο δίκαιο. Το βρήκα και στα διαγράμματα σελίδες 103-4.
Όταν ένα σύστημα είναι αιτιατό τότε y[n]=0 για n<0.

Ευχαριστώ

Στην είσοδο x(n)=cos(αn) το (αn) είναι σε rads ή μοίρες;

Εγώ το δοκίμασα και στα δύο. Σε μοίρες μου βγαίνουν φθίνουσες καμπύλες, ενώ σε rad μου βγαίνουν
ημιτονοειδή καμπύλες. Έτσι το θεώρησα rad...

To eTHMMY σηκώθηκε.

Το σωστό είναι σε rad.

Ανακοινώθηκε η δεύτερη εργασία!
Πληροφορίες στο εΤΗΜΜΥ

2η Εργασία

1. Δεν κατάλαβα πώς από την H_d(ω) θα βγάλουμε την h_d(n) ή την H_d(e^jΩ);

2. Στο παράθυρο Kaiser τι τιμή θα βάλουμε στο Ω_α;


Πάντως Βρήκα το τροποποιημένο πολυώνυμο Bessel:

I₀(x)  =  1 + x²/2² + x⁴/2²4² + x⁶/2²4²6² + ...

Υπάρχει στο τυπολόγιο που είχαμε αγοράσει για τα Σήματα στη σελίδα 161-Τύπος 27.35.
Στη σελίδα 261 έχει και πίνακα με τιμές του πολυωνύμου.

την h_d(n) την βρίσκεις με αντίστροφο Fourier

Όμως ο αντίστροφος μετασχηματισμός Fourier είναι άθροισμα για ω=0 έως ω=Ν-1...
Δηλαδή, δε θα λάβουμε υπόψην το τμήμα της H_d(ω) για τα αρνητικά ω;;

οχι με DFT αλλα με IFT διακριτου χρονου(είναι διαφορετικο,δεν ξερω τη διαφορα αλλα ετσι ειναι)
δες τις σημειωσεις του στριντζη :
DFT:σελ 19
IFT:σελ 17 σχέση 3.1.3
οπότε δεν είναι Σ αλλα ολοκλήρωμα απο -π/4 ως π/4.

ΥΓ:βρες το ξεχωριστα για n=0 και για n!=0

Δηλαδή θα το κάνουμε με IFT...
Τότε όμως λέει Ω=ω*Τ. Τι τιμή θα δώσουμε στο Τ;;

Βρήκα στο βιβλίο του Πανά (σελ. 52) ότι η περίοδος θεωρείται 2π και τότε για τον αντίστροφο μετασχηματισμό Fourier ισχύει ο τύπος 3.1.23. Απάντησα έτσι στη δική μου ερώτηση, αλλά το αφήνω έτσι μήπως έχει κανένας την ίδια απορία.

Τελοσπάντων έλυσα το ολοκλήρωμα και βρήκα:

h_d(n) = (1/nπ) * j * sin(nπ/4)      για n!=0
και
h_d(0) = 1/4

Σωστό είναι αυτό;; Θα πάρουμε για h_d(n) το μέτρο του αποτελέσματος ώστε να φύγει το j;;
Σας παρακαλώ!!! Βοήθεια!!

έκανες λάθος στην ολοκλήρωση.το σωστο ειναι:
h_d(n) = (1/nπ)  sin(nπ/4)      για n!=0
και
h_d(0) = 1/4
το είχε κανει ο μεταπτυχιακος οποτε ειναι σιγουρα αυτο.

Ναι λογικά έχεις δίκαιο.
Για να λύσεις το ολοκλήρωμα γράφεις τον exp(jωn) = cos(ωn) + j sin(ωn);

Και τελικά με το Ω_α του παραθύρου Kaizer τι γίνεται; Τι τιμή παίρνει;

Στην σελ. 194 στο βιβλίο του κ.ΠΑΝΑ γράφει οτι για το  Ωα του παραθύρου Kaizer
παίρνουμε τιμεσ στην περιοχή:       4<Ω(Ν-1)/2<9

Ευχαριστώ MIMIRU!!

Το βιβλίο του Πανά τα σπάει σε σχέση με τις σημειώσεις του Στριντζή!...

Για την άσκηση 1 αφού βρω τις τιμές για τα w(n) και hd(n), βρίσκω το h(n)=hd(n)w(n), οπότε έχω Ν τιμές για το h(n). Πως όμως μετά θα βρω το H(ω) για να σχεδιάσω και το μέτρο της?Υπάρχει κάποιος τρόπος να το υπολογίσω το H(ω) από το h(n) μέσω matlab? Βρήκα μία συνάρτηση, την freqz η οποία μάλλον παίρνει σαν ορίσμα την h(n) και δίνει την H(ω),αλλά από το διάγραμμα δεν μου φάίνεται ότι είναι το σωστό αποτέλεσμα.
 Ξέρει κανείς αν είναι η συνάρτηση που χρησιμοποιείται γενικά γι' αυτό το θέμα?   

Όποιος έχει ασχοληθεί με την άσκηση 2...ας δώσει μια μικρή βοήθεια σχετικά με την yulewalk...
Αυτό που βρίσκει η yulewalk από ότι είδα στο matlab είναι συντελεστές από πηλίκο δύο πολυωνύπων 30ης τάξης το καθένα.. Αυτό όμως που ζητάει αυτός είναι συντελεστές από 30 πηλίκα πολυωνύμων 1ης τάξης το καθένα...Πως θα τα συνδυάσουμε αυτά τα δύο..??Με επιμεριστική αν το σκεφτώ προκύπτουν για αριθμητή και παρανομαστή ξεχωριστά 30 εξισώσεις με 60 αγνώστους συνολικά...!!!Άρα δεν μπορούμε να βρούμε τις τιμές των συντελεστών που ζητάει στην εκφώνηση από τους συντελεστές που υπολογίζουμε με την yulewalk?
Σκέφτομαι κάτι λάθος και αν όχι τι άλλος τρόπος υπάρχει?? ^confused^

Για την άσκηση 1 εγώ έβαλα το MATLAB να υπολογίζει τιμές του μετασχηματισμού Fourier της h_d(n) και μετά τα απεικόνισα...

Για το Hamming (N=7): Έβγαλα ένα χαμηλοπερατό φίλτρο που έχει ενίσχυση 0.27 για κάτω από π/4 μετά πέφτει και φτάνει στα 0.0075 μετά τα 3π/4.

Για το Kaiser (N=8, Ω_α=1.857: Έβγαλα ένα χαμηλοπερατό φίλτρο που έχει ενίσχυση 2.04 για κάτω από π/4 μετά πέφτει και φτάνει στα 0.002 μετά τα 3π/4.

Δε θα έπρεπε να έχουν ενίσχυση 1 για κάτω από π/4; Εσείς τι βγάλατε;

1)Πειτε καμμια εντολη για να υπολογισουμε το fourier και να το απεικονισουμε? με την fft το κανετε?

2)Σχεδιάσατε το ιδανικό φίλτρο για να το συγκρίνετε (ασκηση 1)? Εμεις το σχεδιάσαμε απο το πεδιο του χρονου (ξέροντας την hd[n]) και μεσω της fvtool  στο πεδιο της συχνοτητας αλλα δεν βγαινει ακριβως ιδανικο. εχει gibbs!
εσεις τι κανατε?

3)Για το Καιζερ τι Ω θα επιλεξουμε? αφου δεν ξερουμε την συχνοτητα δειγματοληψιας? θα την επιλεξουμε εμεις?

thanks

1)εγω βρήκα τον τυπο με το "χερι"(0,02+Σ...απο 0 εως Ν-1)και απλά υπολόγισα με ενα loop το άθροισμα.μετα plot του μετρου της συνάρτησης που υπολόγισε.

2)ουσιαστικά θέλουμε να βάλουμε την u(w+p/4)-u(w-[/4)
δεν ξερω πως γίνεται αυτο αλλα μου είπε ο μεταπτυχιακός οτι γινεται με εισαγωγη πίνακα.
θα το δω και αν το κάνω θα το κανω post

3)οτι δεν ειναι σαφες αυθαίρετα αποτι φανταζομαι

thanks αλλα τελικα το εκανα με την freqz για οποιον ενδιαφερεται.

Πως το έκανες με τη freqz;

Εγώ το Kaizer το έκανα όπως η elisabeth αλλά δεν είμαι καθόλου σίγουρος για την ορθότητα της λύσης.
π.χ. για Ν=23 βγάζω το διάγραμμα που επισυνάπτω. Μπορεί να το δεί κάποιος να μου πει αν είναι σωστό;;

Για επιβεβαίωση και για όσους δεν είναι γραμμένοι στην electoauth..
Μόλις διάβασα στη λίστα ότι η εργασία μεταφέρεται για Μ.Δευτέρα.
Στάνταρ αυτό???

ναι το είπε ο μεταπτυχιακος στην τάξη

Δεν θελω να σε απογοητευσω αλλα δεν μου φαινεται σωστο.
Εγω εβγαλα ενα που παει αρκετα καλα προς το ιδανικο με τις κυματωσεις του κλπ
Για να δουλεψεις την freqz βρισκεις την κρουστικη αποκριση του φιλτρου σου στο πεδιο του χρονου (πολλαπλασιαζοντας την ιδανικη hd[n] με το παραθυρο w[n] για Ν ορους).
Το βαζεις στην freqz και στο παει στο πεδιο της συχνοτητας και επειτα το σχεδιαζεις.

Δεν το βρίσκω με τίποτα... Και με τη freqz τα ίδια αποτελέσματα έχω σαν αυτά που επισύναψα πιο πάνω...
Κάπου αλλού έχω λάθος...

Jason μπορείς να επισυνάψεις ένα διάγραμμα; Να δω τι ψάχνω ο καημένος...

κατι τετοιο...

Έχω απογοητευτεί, δε μου βγαίνει με τίποτα... 5 μέρες το ψάχνω...
Επισυνάπτω τον κώδικα MATLAB για το Hamming. Σας παρακαλώ κοιτάξτε το και πείτε μου που έχω λάθος;;;;

Δεν ξερω κατα ποσο θα παιξει ρολο αυτο αλλα δοκιμασε να κανεις το Loop σου 
for m=-(N-1)/2:1:N/2 να μπαινουν και αρνητικες τιμες δλδ.
Για την δεικτοδοτηση ομως των πινακων (h(n) κλπ) θα χρησιμοποιησεις αλλη μεταβλητη.

 ^banghead^ Ότι νάναι βγάζει πάλι...

1. Δηλαδή άλλο λάθος δε βρήκες; Τι στο διάκαλο φταίει τότε;
2. Εσύ τι όρια έβαλες στα αθροίσματα;
3. elisabeth εσύ πως το έκανες;;

 ^banghead^

Κάντε ένα διαλλειμματάκι από την εργασία και τσεκάρετε τα java applets αυτής της σελίδας
http://www.dsptutor.freeuk.com/index.htm.

 Digital filters are fun ! 8)

Jason σίγουρα σωστό είναι αυτό το διάγραμμα που έδωσες πριν...???εμένα μου έβγαλε κάτι τέτοιο για διάφορα n πάλι.....
νομίζω ότι πρέπει να τέμνει τον άξονα της συχνότητας στο ω=π/4=~0,78..όπως στο δικό μου....

τι λες??δεν είμαι και απολύτως σίγουρη.. :-\

Εκει τον τεμνει και το δικο μου!
Αν προσεξεις ο οριζοντιος αξονας λεει στις μοναδες
Normalized frequency ( x π rad/s)
ειναι δηλαδη κανονικοποιημενος με 1 να δειχνει την συχνοτητα αναδιπλωσης (ή μισο της συχν δειγματοληψιας).
αρα 0,25*Pi !

εγω εκανα ενα Loop οπως σου ειπα απο -(Ν-1)/2 εως Ν/2 (δεν ξερω αν ειναι σωστο)

Η ιδανικη σου ειναι η hd(n)=((sin(m*wc))/(m*pi));
για n=0 ομως ειναι = 2Fc

Αυτο θα το πολλαπλασιασεις με το παραθυρο που θελεις!

ναι Jason το είδα μετά ότι το έχετε κανονικοποιήσει.... 8)

Τί είναι το wc;; Εγώ έχω pi/4...


Βασικά πήγα στο site που λέει η (Συγνώμη) miss_elec  και αντέγραψα τα αποτελέσματα για το h(n) στο MATLAB και έβγαζα καλές αποκρίσεις. Άρα το λάθος μου είναι στον υπολογισμό των h(n)...

Εσύ το βιβλίο του Πανά χρησιμοποίησες ή τις σημειώσεις του Στρίντζη;

Αααααα!!!!  :o Όχι και "ο" miss_elec!! Miss..Υou know, κoπέλα, δεσποινίς!! :) ;D :)

ναι το wc το εχω pi/4

δοκιμασες να το Plotρεις μονο του?

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

πρωτα έχεις λαθος το φίλτρο(βασικα ο στρίντζης το έχει λαθος)
εγω έκανα:


w=[-4*pi:0.01:4*pi];
h=0,25
for n=1:(N-1)
    h(n)=sin(n*pi/4)*(0.54+0.46*cos(2*n*pi/(N-1)))/(n*pi);
end
y1=abs(h)
plot(w,y1)

και ενω μου φαίνονται σωστα μου βγαινει περιοδικο.
επισυναπτω και το διαγραμμα.

Αν εννοείς αυτό: πήρα την h(n)=sin(n*pi/4)/(n*pi) , δηλαδή χωρίς παράθυρο.
Η απόκριση που έιχα ήταν καλή μόνο για άθροισμα από -Ν/2 έως Ν/2 και είναι αυτή που επισυνάπτω.
Ίσως τελικά το πρόβλημα να είναι στην w(n). Τέσπα θα κάνω μία επανασχεδίαση σύμφωνα με το βιβλίο του Πανά και θα δω.

...και ευχαριστώ για το όλο ενδιαφέρον.

Σχετικά με την yulewalk και πιο συγκεκριμένα με το όρισμα F που παίρνει...
Λέει στο matlab ότι είναι ένας πίνακας με πρώτη τιμή 0 και τελευταία 1...και το 1 αντιστοιχεί στο μισό της συχνότητας δειγματοληψίας...Το μισό της συχνότητας δειγματοληψίας είναι pi/4..σωστά?
Για ενδιάμεσες τιμές σε αυτόν τον πίνακα τι θα πάρουμε??
Νόμιζα πως ενδιάμεσες τιμές θα πέρναμε ίσες με pi/4... :-\

Ο πινακας αυτος παει σε συνδυασμο με τον πινακα Μ που λεει και στην ουσια σχεδιαζουν την ιδανικη σου αποκριση!

Ο F ξεκιναει με 0 και καταληγει σε 1 (ο αξονας ειναι κανονικοποιημενος - αρα *π).
Ενδιαμεσες τιμες ειναι η συχνοτητα αποκοπης δλδ 0,25

Ω θεοί!!! Τα κατάφερα επιτέλους!!...

Το σωστό παράθυρο Hamming είναι:

0.54+0.46*cos(2*n*pi/N) για -Ν/2<=n<=Ν/2, όπως λέει το βιβλίο του Πανά.

Οι σημειώσεις είναι για τον πούτσο!! Και ας μου επιτραπεί η λέξη. Κάπου πρέπει να ξεδώσω. Μία βδομάδα βγάζω λάθος αποκρίσεις γιατί οι σημειώσεις είναι για κάψιμο... Ζήτω το βιβλίο του Πανά!!


Όσον αφορά το παράθυρο Kaizer (ισχύουν τα όρια -Ν/2<=n<=Ν/2) ακόμα δεν έχω καταλάβει τι τιμή πρέπει να δώσουμε στο Ωα... Αφού έχω τόσο καλύτερη απόκριση όσο μικρότερο είναι το Ωα. Συγκεκριμένα παίρνω την καλύτερη απόκριση για Ωα=0. Αλλά αυτό δεν έχει νόημα, αφού τότε προκύπτει το τετραγωνικό παράθυρο...
Όσο για τις τιμές 8/(Ν-1)<Ωα<18/(Ν-1) , που προτείναι το βιβλίο του Πανά οι αποκρίσεις είναι χάλια, όπως αυτή που επισυνάπτω...

Υ.Γ. Εντάξει και το βιβλίο του Πανά δεν είναι τέλειο...

Ακόμη και στο Ματλαβ, το παράθυρο αυτό το έχει όπως στις σημειώσεις μας.
Τώρα γιατί στο βιβλίο του Πανά το έχει διαφορετικά, δεν ξέρω....

Γιατί παιδεύεστε με τυν υλοποίηση παραθύρων? Νομίζω ότι υπάρχουν και έτοιμα...


HAMMING   Hamming window.
    HAMMING(N) returns the N-point symmetric Hamming window in a column vector.
 
    HAMMING(N,SFLAG) generates the N-point Hamming window using SFLAG window
    sampling. SFLAG may be either 'symmetric' or 'periodic'. By default, a
    symmetric window is returned.

Κάτι χαζό.

Γνωρίζουμε ότι Ω=ωΤ.

Η Wn στο matlab τι παριστάνει? Τη συχνότητα αποκοπής φυσικά κανονικοποιημένη ως προς τη Nyquist/2.

Ποια είναι όμως η Nyquist? Η συχνότητα δειγματοληψίας?
Μπορώ να βάλω αυθαίρετα μία τιμή μεγαλύτερη του pi/2?
Μπορεί κάποιος να βοηθήσει?

Ευχαριστω!

Η yulewalk μου βγάζει (Ν+1) τιμές ως αποτέλεσμα!... Τι είναι αυτές οι τιμές;;
Εγώ περίμενα 2(Ν+1) για αριθμητή και παρανομαστή...

δες καλα το help!
η yulewalk συντασεται ως εξης
[β, α]= yulewalk...
αρα στο β εχεις Ν (αριθμητη) και στο α αλλες Ν (παρονομαστη).

Τασο νομιζω το βρηκα!
εχει ξεχασει ενα Τ εκει που παει να λυσει τον τυπο βεβαια.

αυτο δεν ρωτησες?

Στη yulewalk το Ν=30 αν δεν κάνω λάθος.

Στην εκφώνιση της εργασίας λέει:

Διαλέξτε τα a_k, b_k, c_k, d_k ώστε να ελαχιστοποιείται το μέσο τετραγωνικό σφάλμα της προσέγγισης.

Αυτό όμως δεν προκύπτει από τη yulewalk αφού δίνει τους συντελεστές της συναρτησης μεταφοράς:

          B(z)   b(1) + b(2)z^-1 + .... + b(n)z^-(ν-1)
          ---- = ---------------------------------------
          A(z)    1   + a(1)z^-1 + .... + a(n)z^-(ν-1)
 
Τι θα δώσουμε ως απάντηση στο δεύτερο ερώτημα; Την H(z); Το διάγραμμα της απόκρισης της Η(ω);
Τα εν σειρά φίλτρα που να τα βρούμε με Ν=30;...

Εσείς τι απαντήσατε;

διαγραμμα της αποκρισης θα βαλεις και την Hz Παραγοντοποιημενη στη μορφη που δινει η εκφωνηση.
βεβαια λιγο δυσκολο να το κανεις για 30...
καντο για 3-4 και μετα για 30 απλα γραψε τα αποτελεσματα!

Είναι σίγουρα αυτό? Σε ένα ebook που κατέβασα το έχει όπως στις σημειώσεις του Στρίντζη. Το κακό είναι ότι και το βιβλίο του Πανα και οι σημειώσεις του Στρίντζη φημίζονται τα λάθη τους!

Εγώ πάντως το έκανα και με τους δύο τρόπους. Με τις σημειώσεις είχα άθλια αποτελέσματα αντίθετα με το βιβλίο του Πανά...

Για να φέρω τη Η(z) στη μορφή που ζητάει η εκφώνιση χρησιμοποίησα τη συνάρτηση roots του MATLAB. Το πρόβλημα είναι ότι μου βγάζει και μιγαδικές ρίζες... Πρέπει οι παράγοντες να είναι πρώτης τάξης ή μπορεί να είναι και δεύτερης; (όπως στο βιβλίο του Πανά, σελίδα  109, τύπος 5.7.6)

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

Επισυνάπτω τον κώδικα (χρησιμοποίησα την freqz) για το πρώτο ερώυημα και το διάγραμμα .Υπάρχει κάποιο λάθος και δε βγαίνει σωστή η ενίσχυση.

Ποιο είναι το λάθος; Faidon μήπως ξέρεις;Η κανένας άλλος;Ρε παιδιά, ΣΟΣ..κοντεύει η  Μεγάλη Δευτέρα ^rotate^

Ξέρει κανείς αν χρειάζεται στην αναφορά να βάλουμε τίποτα ενδιάμεσα αποτελέσματα όπως συντελεστές παραθύρων,αποτελέσματα της yulewalk κ.τ.λ;Κάτι είχε πει μέσα στο μάθημα να βάλουμε κάποια αποτελέσματα μέσα σε txt αλλά δεν είχα ακούσει τι ακριβώς.Θυμάται κανένας τι ακριβώς είχε πει;Λέω μήπως είναι τίποτα αναγκαίο για την αναφορά

Αυτο ειναι λιγο υποκειμενικο. Οτι κρινεις απαραιτητο για να παρουσιασεις τον τροπο που εφτασες στη λυση σου το βαζεις!

Με απογοητεύεις Megawatt. Δεν κοίταξες τα προηγούμενα ποστ; Υπάρχουν όλες οι απαντήσεις.
Τέλως πάντων επειδή βλέπω ότι είσαι ακόμα στην αρχή...

1.Τα αθροίσματα στα παράθυρα να τα κάνεις από -Ν/2 έως Ν/2.

2.Το παράθυρο Hamming είναι 0.54+0.46*cos(2*n*pi/N) για -Ν/2<=n<=Ν/2.

3.Στο παράθυρο Kaiser πρέπει να βάλεις Ωα=1/(Ν-1).

4.Μπορείς να χρησιμοποιήσεις έτοιμες συναρτήσεις του MATLAB. (π.χ. hamming(N) ή kaizer(N,0.5), δες τα help)

5.Το ερώτημα 2 κάντο με τη yulewalk (δες το help). Και μετά χρησιμοποίησε τη roots (δες το help) για να παραγοντοποιήσεις.

Τελικά και οι 2 σωστά πρέπει να το έχουν:
Στον Πανά :   ω(n)=0,54 + 0.46cos(2πn/N)  για  -N/2<= n <= N/2
Στον Στρίντζη :  ω(n) = 0,54 - 0.46cos(2πn / N-1)  για  0<= n <=N-1 

Ο τύπος αλλάζει ανάλογα με το τι όρια δίνεις στο n , αλλά το τελίκο αποτέλεσμα είναι το ίδο . Το έλεξγα με τα όρια του n και
στις 2 περιπτώσεις και μια τυχαία ενδιάμεση τιμή.
Δοκίμασε να βάλεις διάφορες τιμές του n με το χέρι να δεις ότι βγαίνουν τα ίδια.

Στο Matlab ακόμα δεν το δοκίμασα  :)

Δεν έχεις τι να κάνεις μου φαίνεται...
Αν το δοκιμάσεις στο MATLAB πες τι έβγαλες γιατί έχω και εγώ απορία! :)

Νομίζω είπε να βάλουμε τους συντελεστές a,b της yulewalk σε ξεχωριστό αρχείο.
Ξέρει κανείς πως από τις ρίζες το παραγοντοποιούμε σε 1ης και 2ης τάξης;να βρούμε α_k,β_k,c_kκτλ...

Ρε συ σούπερνόβα.... η yulewalk βρίσκει τα καλύτερα ακ, βκ, ck,.... και η roots μετά τις ρίζες από το τέρας εκείνο.

Σωστά;

Αφού έχεις τους συντελεστές των πολυωνύμων αριθμιτή και παρανομαστη, έχεις και τις ρίζες, τι άλλο μπορεί να θέλεις από την ζωή σου;

Την προχώρησες καθόλου την ιδέα σου Τουραμ να δεις τι γίνεται?

Για την yulewalk....
Παιδιά μήπως μπορεί κάποιος να μου πει ποιες είναι οι ενδιάμεσες τιμές που πρέπει να μπουν στο διάνυσμα f [0,........,1] και γιατί??
Καταλαβαίνω ότι υπάρχει αντιστοιχία μεταξύ m και f διανύσματος.
thanks.

Από τις ρίζες παιδιά αν θεωρήσουμε ότι το πολυώνυμο παραγοντοποιείται (z^-1-r1)*(z^-1-r2) ktl
και πολ/σουμε ανα 2 τις παρενθέσεις μπορούμε να βρούμε τους συντελεστές της δευτεροβάθμιας σε συνάρτηση με τα r1,r2.

με τα m,f ουσιαστικά ορίζεις την συνάρτηση του φίλτρου που θέλεις να προσεγγίσεις.

το f είναι ο άξονας των συχνοτήτων και το m η αντίστοιχη τιμή του φίλτρου σου.
τα ζευγάρια τιμών f(i),m(i) προσδιορίζουν το φίλτρο σου. δηλαδή στην συχνότητα f(i) το φίλτρο σου έχει απόκριση m(i)
οι τιμές του f είναι κανονικοποιημένες και γι αυτό το f(i) μπορεί να έχει τιμές μεταξύ 0-1
πχ

f m
0 1
0.01 1
ktl

Ωραία BoboMastora!
Δηλαδή για το συγκεκριμένο φίλτρο το f=[0,0.25,0.25,1] και m=[1,1,0,0] είναι σωστό?

Όντως μπορούμε να τις βρούμε. Το θέμα είναι γιατί να το κάνουμε; Αφού ζητά υλοποίηση με φίλτρα της μορφής [ak+bk*exp(-jω)]/[ck+dk*exp(-jω)]. Οπότε από τη roots θα πάρουμε τα -bk και -dk ενώ τα αk και ck θα είναι πάντα 1. Υπάρχει πρόβλημα με αυτό;


Σωστό είναι, χωρίς να είναι και το μοναδικό. Εξαρτάται τι συχνότητα δειγματοληψίας θεωρείς ότι έχεις. Εδώ είναι σαν να λες ότι fs=2π.

Mε τη roots βρήκαμε τις 2 πρώτες συζυγείς  ρίζες για τον αριθμητή:
  -0.9906 + 0.1506i
  -0.9906 - 0.1506i
Και για τον παρονομαστή:
-0.5663 + 0.8195i
  -0.5663 - 0.8195i 
Μετά από πράξεις καταλήξαμε στο πρώτο φίλτρο 2ης τάξης

(0.996*z^-2+1.973*z^-1+1)
------------------------------------------------
(1.0078*z^-2+1.141*z^-1+1)

Βρήκε κανείς κάτι τέτοιο;;Φυσικά  θα το κάνουμε για κανα 2 ζευγάρια μόνο
και όχι σε ματλαβ,ε;;; :???:

(0.996*z^-2+1.973*z^-1+1)
------------------------------------------------
(1.0078*z^-2+1.141*z^-1+1)

Κάποιο λάθος θα έχεις κάνει. Αυτό που έχω υπογραμμίσει θα πρέπει λογικά να είναι μικρότερο του 1, αφού είναι το τετράγωνο του μέτρου του πόλου. Και ως γνωστόν για να είναι το φίλτρο ευσταθές πρέπει |z₂|<1 άρα και |z₂|²<1 για κάθε πόλο.

Ο συντελεστής  του z^-2 βγήκε 1/(r1*r2) οπότε  μπορεί να βγει μεγαλύτερο του 1!;
πρέπει να τις κάνουμε τις πράξεις ή απλά να αναφέρουμε πώς από τη yulewalk πάμε στη μορφή της εκφώνησης;;;;

Για μιγαδικές ρίζες a+jb και a-jb που έχουν μέτρο sqrt(a²+b²) ισχύει:

(z-a+jb)(z-a-jb) = z² - 2az + a²+b² = [1 - 2az^-1 + (a²+b²)z^-2]*z²

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

Η τρίτη εργασία στο eTHMMY!! Προθεσμία παράδοσης Δευτέρα 8-5-2006.

Κανείς?

Ακόμα χωνεύετε?
ΟΕΕΕΕΟ?

Έτοιμος για δράση ε;
Εγώ σήμερα την έπιασα. Λέω να κάνω μία FIR υλοποίηση. Ακόμα δεν είμαι σίγουρος για τη συχνότητα αποκοπής. Πάντως αν δεν κάνω λάθος η μπάσα νότα είναι η ΛΑ, για να μην κοπεί θέλει συχνότητα αποκοπής μεγαλύτερη του 440Hz.

τι θα λέγατε να μεταφέρουμε την συζήτηση της τρίτης εργασίας εδώ... (http://www.thmmy.gr/smf/index.php?topic=5091.msg62459#msg62459)

http://www.thmmy.gr/smf/index.php?topic=5091.msg62459#msg62459


για να μπορούμε να βρίσκουμε αυτά που θέλουμε ακριβώς και να μη χανόμαστε στα παλιά ποστ..