THMMY.gr

Οποιαδήποτε απορία έχετε για τις εργασίες του μαθήματος.

Ανέβηκε η 1η Εργασία. Λήξη προθεσμίας υποβολής Δευτέρα 5 Απριλίου 2021 (23:59).

Μπορεί να ειναι χαζό, αλλά για την εύρεση του θο, έχουμε ένα άθροισμα γινομένων πινάκων, του οποίου ψάχνουμε τον αντίστροφο ή αλλιώς το (Σ(φ•φ^Τ))^-1. Αυτό το φ είναι το διάνυσμα οπισθοδρόμησης που έχει βρεθεί από το 1α, δηλαδή τα 3 φιλτραρισμένα σήματα από τις εισόδους-εξόδους;

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

Άρα για να καταλάβω, τελικά ο προς αντίστροφη πίνακας είναι οι Ν μετρήσεις (άρα Ν γραμμές) επί το εκάστοτε φίλτρο (3 στήλες)
Σωστά;

Τώρα δε θυμάμαι ακριβώς ποιος αντιστρέφεται, αλλά με βάση τον τύπο που έγραψες (Σ(φ•φ^Τ))^-1 το φ έχει τις διαστάσεις που λες ναι, άρα το γινόμενο είναι Νx3 • 3xN = NxN

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

Έχει κανεις κανα λινκ με θεωρια διαφορετικη απο τις σημειωσεις του για αυτην την γραμμικη παραμετροποιηση που κανει με τα φιλτρα, γτ αυτος τα δεχνει λιγο περιεργα εκει με τα Δ

Θέμα 2, πρώτο υποερώτημα, χρειάζεται και η θεωρητική μελέτη, ή μπορούμε να το κάνουμε κατευθείαν στο matlab?

Στη δεύτερη πως γίνεται να πάρω τιμές της vp για πολλές χρονικές στιγμές χωρίς να τις καλώ μία μία; Γιατί στην v(:,1) βγάζει σφάλμα όπως επίσης και όταν ορίζω τον χρόνο σαν διάνυσμα.

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


Αν το διάνυσμα των χρονικών στιγμών σου είναι π.χ. t, δοκίμασε να το γράψεις μέσω μιας for σε v(t(i)).

Συγνώμη, να κάνω μια χαζή ερώτηση. Ωραία, το φ είναι το διάνυσμα οπισθοδρόμησης. Κάθε στοιχείο του μέσα αποτελείται από ένα κλάσμα με s , το οποίο είναι είτε πολλαπλασιασμένο με y ή με u. Όταν εγώ θα πάω για να βρω το διάνυσμα φ για μια συγκεκριμένη χρονική στιγμή , το κλάσμα αυτό με τα s θα το έχω πάει στο πεδίο του χρόνου με αντίστροφο laplace  και τότε θα βάλω όπου t την χρονική στιγμή μου ή ολο το φ θα πρέπει να το έχω βάλει σε αντίστροφο Laplace και άρα κάθε στοιχειο του φ πρακτικά θα αποτελείται από συνελίξεις;

Κάθε στοιχείο του Φ θα αποτελείται από τα στοιχεία του ζ άρα πχ 4 γραμμές και οι στήλες θα είναι όσες και οι χρονικές στιγμές που θες. Για να το κάνεις αυτό πρέπει το κάθε φi να είναι της μορφής lsim(tf(συνάρτηση μεταφορας),y,t) όπου y βάζεις είτε είσοδο είτε έξοδο ανάλογα με το τι έχεις στο ζ. 

Επειδή πολλοί φάγαμε πολύ χρόνο μέχρι να το καταλάβουμε, το γράφω εδώ να υπάρχει:

Η συνάρτηση v που μας δίνει έχει σαν πρώτη συνιστώσα (v(1)) την τάση στον Πυκνωτή και όχι σαν δεύτερη, όπως γράφει στην εκφώνηση.

Αν την πλοτάρετε θα φανεί αμέσως!

Στο δεύτερο θέμα που θέλει εκτίμηση του πίνακα μεταφοράς του κυκλώματος σε τι αναφέρεται; Στους συντελεστές απλά της διαφορικής;

Ουσιαστικά ναι, αφού εκτιμήσεις τις παραμέτρους θα μπορείς να υπολογίσεις και τα στοιχεία του πίνακα μεταφοράς G(s) = Y(s)/U(s).

Ερώτηση. τους κώδικες τους θέλει και σαν αρχεία matlab ή μόνο σαν text στο pdf της αναφοράς;

εγω περασα την αναφορα μου και τα αρχεια .m

σε .rar να φανταστώ

You can never go wrong εάν τα συμπιέσεις. Αν και καλύτερα σε zip θα έλεγα εγώ, αλλά ό,τι σε βολεύει.

Ίσως είναι χαζό αλλά στο πρώτο ερώτημα έχω το εξής πρόβλημα:

Στο πρώτο ερώτημα βρίσκω το θ καπελλάκι, το οποίο είναι 2x1 έχω και το διάνυσμα φ, που κι αυτό είναι 2x1.

Αν θέλω να βρω το x_predicted = θ' * φ το x_predicted θα είναι 1x1.

Αν τα ανάστροφα πάνε ανάποδα, το x_predicted θα είναι 2x2.

Για να το κάνω 2x1, όπως θέλω να είναι, πρέπει να το κάνω ξεχωριστά, θ(1) * φ(1) και θ(2) * φ(2);;

edit: όντως χαζό είναι, η έξοδος που θέλουμε είναι βαθμωτό μέγεθος

Στα thread για την αντίστοιχη εργασία των προηγούμενων ετών λένε ότι οι σημειώσεις έχουν λαθάκια ως προς την γραφή των σχέσεων.
Εγώ το έκανα με το ανάστροφο.

Τα x θα τα πάρουμε με λύση της διαφορικής για τις τιμές των παραμέτρων που δίνονται και μετά θα τα θεωρήσουμε γνωστά στην παρακάτω ανάλυση; Επίσης εφόσον προκύπτει σύστημα διαφορικών με άγνωστα x_cup,θ1cupdot, θ2cupdot πώς μπορεί να λυθεί; Υπάρχει κάπου κάποιο αντίστοιχο παράδειγμα που να έχει γίνει στο μάθημα;

Να ρωτήσω μια χαζή ερώτηση
Για το δεύτερο ερώτημα
Δεν ισχύει για Π ότι θ1_hat_dot = -e*x_hat και
για Μ ότι θ1_hat_dot = -e*x ;

Ενώ για θ2_hat_dot = e*u ?
Εμένα τα α και β πάνε να συγκλίνουνε και στο Π και στο Μ και το Π ενώ πάει να τείνει στις πραγματικές τιμές, μου ξεφεύγει λίγο και ξαναγυρνά (ταλαντώνεται αρκετούτσικα) ενώ το Μ ενώ πάει να βρει τις πραγματικές τιμές , εκεί που πάει να τις βρει, εκτινάσσεται στο σύμπαν .
Δεν έχω ιδέα τι στο καλό πάει στραβά, έχει κανείς καμιά ιδέα;
Μιλάω για πριν εντάξω θόρυβο

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


Το Π όντως πρέπει να κάνει αυτό που λες, αλλά το Μ σίγουρα όχι. Ή θα έχεις κάνει κάποιο λάθος στις υπόλοιπες εξισώσεις (πχ του φ) ή θα έχεις gains (τα γ) που οδηγούν σε αστάθεια

Από e-learning: