THMMY.gr

Η τεταρτη εργασία έχει ανέβει στο e-thmmy!!

...και εδώ!!  :P

Σημειώνω τα αποτελέσματα που βρήκα για τα ερωτήματα 1 και 2:

1. Συνθήκη ευστάθειας:
     -1<μ+λ<1

2.Συνθήκες ευστάθειας:
     |μ|<|1-λ|
     |λ|<|1-μ|
     |μ|+|λ|<1

Ασχολήθηκε κανένας; Είχατε τα ίδια αποτελέσματα;

Απορία: Πώς μπορώ να κάνω plot μια συνάρτηση δύο διαστάσεων στο MATLAB;

Σωστός ;)

Για τα υπόλοιπα ερωτήματα επιφυλάσσομαι να απαντήσω άλλη στιγμή, αν δεν με προλάβει κάποιος άλλος.

H mesh είναι η εντολή που ψάχνεις.  ::)

Για τα υπόλοιπα... (marauber copy paste)  :-X

Ωραία!
Έβγαλα το διάγραμμα που επισυνάπτω. Τι βγάλατε εσείς;

H mesh ξεκινάει από το 1 (καθότι και στον πίνακα y που σχηματίζουμε για τα αποτελέσματα δεν έχει y(0,0) αλλά ξεκινάει προφανώς από το y(1,1).

Μήπως η χρήση της meshz είναι καλύτερη?

Το σχήμα ούτως ή άλλως δεν φαίνεται διακριτό αλλά ενώνει τις τιμές του y.
Θα μπορούσαμε να το κάνουμε αλλιώς να φαίνονται μπαρίτσες ας πούμε?

Βασικά επειδή έβαλα Ν=100 δε φαίνεται διακριτό. Για μικρότερες τιμές φαίνεται όπως τα σχεδιαγράμματα των σημειώσεων. Αυτό που επισυνάπτω είναι για Ν=20.

Και η meshz το ίδιο κάνει από ότι κατάλαβα. Απλά δε δείχνει τις αρνητικές τιμές...
Αυτό που επισυνάπτω είναι για Ν=50. Με μπαρίτσες θα ήταν όντως καλύτερο.

Δηλαδή εσύ βρήκες τα ίδια διαγράμματα;

Σχετικά με την IIR υλοποίηση για την διδιάστατη συνάρτηση μεταφοράς (δεύτερο ερωτήμα), εκτός από την  εξίσωση διαφορών χρειάζεται να κάνουμε και σχήμα (όπως στην μονοδιάστατη μορφή που κάνουμε ευθεία μορφή ΙΙ, ή ευθεία μορφή Ι κτλ.)??
Αν ναι πως ακριβώς είναι το σχήμα, υπάρχει κανένα παράδειγμα πουθενά γιατί δεν βρήκα κάτι?
Θα μπορούσαμε να κάνουμε κάποιο σχήμα αναλογικά με την μονοδιάστατη περίπτωση....αλλά υπάρχουν διπλοί ολοκληρωτές για την περίπτωση που έχουμε z1^(-1)*z2^(-1)??

Και κάτι ακόμα για το τρίτο ερώτημα για τα m και n τιμές στο περίπου δίνουμε?? 

Η εξισωση διαφορων πιστευω ειναι αρκετη. Εγω δεν εκανα διαγραμμα! Αν ενδιαφερεσε ομως ψαξε στην παραγραφο υλοποιηση IIR φιλτρων στις σημειωσεις των 2D.

Οσο για τα m+n θα βαλεις εσυ μεχρι να εχεις μια ικανοποιητικη απεικονιση του συνημιτονου. Αν θυμαμαι καλα το πηγα μεχρι 50.

Για να σας βγει καλα το διαγραμμα κοιταξτε να το στρεψετε καταλληλα ωστε να φαινεται ωραιο.
εγω πχ με mesh εβγαλα το επισυναπτομενο.

Αν θελετε δοκιμαστε και με την stem3 μιας και το συστημα μας ειναι διακριτο αλλα εμενα δεν μου αρεσε οπως τα εμφανιζε.

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

Για το τρίτο ερώτημα....
Στην πρώτη εργασία που ήταν ανάλογη αυτής αλλά για το μονοδιάστατο χρησιμοποιήσαμε την συνάρτηση filter η οποία παίρνει σαν εισόδους τo x στο πεδίο n, τους συντελεστές b και a της διαφορικής εξίσωσης του συστήματος και δίνει σαν έξοδο το y στο πεδίο n.
Προσπαθήσαμε να βρούμε αντίστοιχη συνάρτηση για το διδιάστατο αλλά αυτό που καταφέραμε να βρούμε  είναι την filter2 η οποία είναι για FIR φίλτρα, οπότε άκυρο και την imfilter η οποία δεν περιορίζει το είδος του   
φίλτρου:
B = imfilter(A,H)
Σ' αυτήν την συνάρτηση είναι το Α ένας πίνακας mxn που έχει τις τιμές του x(m,n) και το Β ένας πίνακας mxn που δίνει τις τιμές της εξόδου y(m,n)...σωστά?Το Η όμως τι είναι?Η κρουστική απόκριση ή η συνάρτηση μεταφοράς και με τι μορφή το δίνουμε?
Βρήκε κάποιος άλλη συνάρτηση για να το κάνει??

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

Αφου εχεις την εξισωση διαφορων υλοποιεις απλα αυτην.
κανεις ενα loop που αλλαζει την εισοδο σου και δινει τιμες στον πινακα y(n1,n2).

δλδ
for m
  for n
    if (...)...
   ...
   y(n1,n2) = ...y(n1-1,n2) ... + ... -   + x(n1,n2)....
end
end

Εγώ πάντως χρησιμοποίησα συναρτήσεις του MATLAB όπου δεν μπορούσα διαφορετικά.
Να μάθω και τίποτα...
Όταν θα μου βάλει άσκηση στην εξέταση: "Να υπολογίσεις την έξοδο..."
Τί θα απαντήσω;
"Με τη συνάρτηση filtrer του MATLAB υπολογίζεται η έξοδος.";


Όντως. Αν το στρίψεις φαίνεται πιο ωραίο!!!

Πιστός σε ότι δήλωσα προχτές 8) (αν και άργησα λίγο)

Σωστός και εδώ!

Μερικές απορίες:
 
1)Στο τρίτο ερώτημα, για να υπολογίσουμε μια ορισμένη τιμή του y(n1,n2) χρειαζόμαστε και τις προηγούμενες τιμές y(n1-1,n2) και y(n1,n2-1), σωστά; Άρα για να αρχίσουμε από κάπου δεν χρειαζόμαστε αρχικές τιμές; Ποιες είναι αυτές;

2) Επίσης που μας ζητάει στην άσκηση γραφική παράσταση; Γιατί χρειάζεται να χρησιμοποιήσουμε τη mesh ή κάποια παρόμοια;

3)Γιατί το σύστημα μας είναι διακριτό; Μας λέει πουθενά στην εκφώνηση ότι τα m και n θα είναι ακέραιοι;

Δηλαδή ρε παιδιά θα παραδώσετε εργασία σήμερα? Δεν αποφασίσαμε να απέχουμε από κάθε εκπαιδευτική δραστηριότητα μέχρι την επόμενη συνέλευση? Πάντως το upload στο eTHMMY εχει μπλοκαριστεί. Έχετε και το μέιλ του Μουσ, φαντάζομαι πως θα μου πείτε...

Ρε παιδιά, ενώ έχω εξάγει τα αποτελέσματα  (δλδ τους πίνακες στήλες  X και Y) στο Matlab,
 δεν μπορώ να τα απεικονίσω... Τί βάζουμε σαν όρισμα Ζ στην mesh(X,Y,Z) ?

Ε???? 8) 8) :P :P :-X :-X :-\ :-\ :'( :'( :o :o :o :o :o :o :o :o :o :o :o :o :) >:( :o

βάλε το y.

και προτείνω την χρήση της meshz, για ακριβέστερα αποτελέσματα :)


meshz(y) λοιπόν

1)Το συστημα θεωρειται αιτιατο. αρα για n1 ή n2 < 0 y=0.
2)Ζηταει αν θυμαμαι καλα την εξοδο του φιλτρου. Ε δεν θα του παρουσιασεις ενα πινακα και θα του πεις παρτο!
Θα το κανεις ενα διαγραμματακι!!
3)Κωστη το συστημα μας ειναι διακριτο γιατι ειμαστε στο 8ο εξαμηνο και κανουμε ΨΗΦΙΑΚΑ φιλτρα ;D ;D :P :P
Οι εκφωνησεις των εργασιων φημιζονται για την λακωνικοτητα τους. Αλλα αυτα εννοουνται!

Αυτή η εργασία νόμιζα ήταν για το αναλογικό... τα ψηφιακά φίλτρα τα πέρασα...

φτου σου...

εντάξει, κρύο.

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

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

ποιο είναι τι πεδίο υποστήριξης της χ όμως?

τέλος για την ευστάθεια πήρατε το θεώρημα decarlo-strintzis?
αυτό δεν ισχύει μόνο αν η συνάρτηση μεταφοράς είναι της μορφής :
H(z1,z2)=1/A(z1,z2)?

Να ρωτήσω κάτι και εγώ, ποιος είπε ότι πρέπει να παραδοθεί αυτή η εργασία;!  :???: 8)

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

Οσο για το θεωρημα Στριντζη μαλλον δεν ειδες τι ελεγε πριν. Λεει οτι ο αριθμητης επηρεαζει ΣΠΑΝΙΑ την ευσταθεια και γι αυτο δεν θα ασχολουμαστε κοινως μ αυτον! Οποτε μελετας τον παρονομαστη (ή αν θες ενα συστημα 1/Α) και επεκτεινεις τα αποτελεσματα ως εχουν για το Β/Α !
Κοινως ΞΕΧΝΑ τον αριθμητη.

Προσωπικά την κάνω για να την έχω έτοιμη...


Για να ενισχύσω τη φωνή του Jason, δηλώνω ότι με βρίσκει σύμφωνο σε όλα!!


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

Ρε παιδιά μόλις χθές έμαθα ότι τελικά το eTHMMY ξαναανέβηκε την Κυριακή και ότι ξαναβγήκε ημερομηνία παράδοσης για την Τρίτη.
ΕΛΕΟΣ!! Πρέπει να είμαι κολλημένος στο eTHMMY?? Κατάληψη έχουμε ρε γαμώτο!!
Πολλοί που ρώτησα δεν το είχαν πάρει χαμπάρι επισης.
Δεν γίνεται κάποιο να είναι μέσα στην κατάληψη και να τρέχου  και άλλοι να στέλνουν εργασίες στο eTHHMY.
Μήπως πρέπει να γίνει κάτι με αυτό?
Λέω να στείλω μαιλ στον Μουστάκα.

Φυσικά και να του στείλεις και να πεις το λογικό, ότι δεν ήσουν υποχρεωμένος να ελέγχεις κάθε τόσο το eTHMMY μπας και, εν καιρω κατάληψης..
Εγώ το έμαθα στις 7 η ώρα την Τρίτη και τελευταια στιγμή πλακώθηκα να τελειώσω το τμήμα που μου αντιστοιχούσε! :Ρ Εχασα και την συζητηση στη σχολή για την αξιολόγηση για την παλιοεργασια :(
Πάντως ρε παιδιά του φόρουμ ωραίοι είστε, γράφετε κανονικά στο τόπικ για το πώς θα λυθει η άσκηση, αλλά δεν προθυμοποιήθηκε κανένας να μεταφέρει την ανακοίνωση...

Έπρεπε να ανεβάσω την εργασία για όσους δεν μπόρεσαν να την κάνουν επειδή ασχολούνται με την κατάληψη. Αλλά ο βλάκας δεν το σκέφτηκα.

Έστειλα e-mail στον καθηγητή και στον μεταπτυχιακό και ζήτησα η τελευταία εργασία να μη ληφθεί υπόψην. Δεδομένου ότι την έστειλα (την εργασία) ελπίζω να με πάρουν πιο σοβαρά.

miss_elec σε παρακαλώ να κάνεις το ίδιο και όποιος άλλος ανέβασε την εργασία. Είναι άδικο για όσους τρέχουν για την κατάληψη.

Έτσι για την ιστορία. Εγώ δεν είχα σκοπό να παραδώσω την εργασία τώρα, γι αυτό και ενώ η φίλη μου είχε κάνει το δικό της τμήμα εγώ ούτε είχα ακουμπήσει το δικό μου.
Τρίτη απόγευμα μου τηλεφωνεί η φίλη μου και μου λεει ότι η προθεσμία είναι για σήμερα και ότι σχεδόν όλοι έχουν παραδώσει εργασία  :o Πλακώθηκα λοιπόν να την στείλω κι εγώ εμπρόθεσμα.

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

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

Αυτά τα εκ των υστέρων να στείλουμε μέιλ εγώ τα ακούω βερεσέ. Οι κινήσεις πρέπει να γίνονται σε συνεννόηση με τους υπολοίπους συμφοιτητές. Ή όλοι ή κανένας. Στην δικιά  μας περίπτωση βέβαια, ήταν όλοι εκτός αυτών που φαίνεται πως την πάτησαν...

Βέβαια, ο Στρίντζης πιστευω να μην έχει πρόβλημα τελικά με τους εκπρόθεσμους και να κατανοήσει.

Είναι όμως η ρημάδα η νοοτροπία παιδιά...

Ο μεταπτυχιακός μου απάντησε ότι:

Η εργασία μπορεί να παραδωθεί μέσω e-mail (digitalfilters@olympus.ee.auth.gr) μέχρι και την Τρίτη 6/6.


Επισυνάπτω τις λύσεις των ερωτημάτων της εργασίας.

Σας παρακαλώ:
1.Μην απαξιώσεται την κίνησή μου κατεβάζοντας τις λύσεις επειδή βαριέστε να κάνετε την εργασία.
2.Όσοι ασχολείστε με την κατάληψη να τις κατεβάσετε και να παραδώσεται την εργασία.

Ευχαριστούμε Φαίδων.