THMMY.gr

Η πρώτη εργασία ανέβηκε στο elearning, στην ενότητα " Εργασίες και Τεστ ". Προθεσμία Παράδοσης: 14/11/2018

Ας ξεκινήσω εγώ  :P

Στο θέμα 2, ως εξόδους του συστήματος θεωρούμε τα Vr και Vc;

εγω ετσι καταλαβα :P

κι προσπαθω να φερω το κυκλωμα στην μορφη:
VRdot = ...
VCdot = ...

και ομολογω οτι δεν ειχα καμια τυχη  :P  :-\

Μάλλον είναι πολύ απλό που θα ρωτήσω αλλά κάτι δεν καταλαβαίνω. Στο πρώτο θέμα τις μετρήσεις του y από πού τις παίρνουμε για να τις χρησιμοποιήσω στον τύπο: θ_ο^Τ (Φ^ΤΦ) = yΦ ;

γενικα, για να υλοιποιησεις το 1β, χρειαζεσαι πραγματα που δεν εχεις, οπως το y που αναφερεις!

ετσι λυνεις με καποιον τροπο το 1γ (πχ ode) και παιρνεις το διανυσμα εξοδου (αλλα κι αυτο της εισοδου) και λυνεις το 1β.

Απ' ότι κατάλαβα εγώ, στο 1β απλά ζητάει να κάνεις την συνάρτηση που υλοποιεί την μέθοδο των ελαχίστων τετραγώνων, δηλαδή θα παίρνει ως είσοδο τον Y και το u και θα υπολογίζει το θ₀.

ας συνεχισω εγω τις ερωτησεις :Ρ

εν τελει στο 2ο θεμα, καταληξατε σε μια διαφορικη 2ου βαθμου και μια μηδενικου βαθμου?

bonus ερωτηση:
Οι διαφορικες, ειναι πεπλεγμενες ή μη?

Εγώ έβγαλα δύο διαφορικές 2ου βαθμού μια για το Vr και μια για το Vc. Αυτή για το Vr περιέχει μόνο παραγωγους του Vr και τις εισόδους με τις παραγώγους τους.Ομοίως και για το Vc.

αρα πρεπει να εχεις κι u δευτερη παραγωγος, στην εξισωση για το Vr?

επισης δεν μας ενοχλει που το συστημα πρεπει να βγει 2ου βαθμου?

Ναι, έχει και δεύτερες παραγώγους.

Δε ξέρω. :P

εε αυτο θελω να αποφύγω... γιατι μετα εχω θεμα αν κανω τα φίλτρα με ode, εμφανιζεται το udot

καταλαβα...ολοι στην ιδια φαση ειμαστε

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

Μια ερώτηση για την 1η άσκηση.. Η διαδικασία που ακολουθούμε για να φτάσουμε στο γραμμικά παραμετροποιήμενο μοντέλο πρέπει να υλοποιηθεί και αλγοριθμικά; Η αρκεί να την υλοποιήσουμε στο χαρτί και απλά να χρησιμοποιήσουμε το τελικό μοντέλο για το ερώτημα 1β;

Μετά απο ερωτηση στο μάθημα της Παρασκευής ειπε οτι θεωρουμε γνωστα μονο u και y, οχι τις παραγωγους

Γίνεται να βγει όμως η διαφορική κάπως και να μην έχει τις παραγώγους των εισόδων;

κοιτα να δεις :Ρ
εγω γενικα εχω στο μυαλο μου οτι επειδη εχω L,C πρεπει να μου βγει 2ης ταξης.

Ετσι κατεληξα μεσω αναλυσης (κυκλωμα στο λαπλας, βρογχος κλπ) σε μια εξισωση 2ου βαθμου ως προς Vc και εχω μια μηδενικου βθαμου που συνδεει Vr με Vc.

Για το αν ειναι σωστο θα σε γελασω :P αλλα μου φαινεται σωστο :D

ΥΓ. στο Vc 2ου βαθμου, εμφανιζεται η παραγωγος της εισοδου, ωστοσο οταν πας να κανεις τα φιλτρα εμφανιζεται μονο η εισοδος:

---

νομιζω μονο στο χαρτι...

Τωρα αν εσυ θες να κανεις τον κωδικα σου γενικο και να αντιλαμβανεται μονος του την ταξη του συστηματος και φτιαχνει το ζ, be my guest :P

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

μαλλον plotαρουμε διαφορετικα την εξοδο :Ρ

καταλαβα ολοι διαφορετικη υλοποιηση εχουμε  :D

Επίσης στο διάνυσμα ζ που προκύπτουν κάποιες πράξεις μεταξύ της συνάρτησης μεταφοράς του σχεδιασμένου φίλτρου(πεδίο laplace) και της εισόδου ή της εξόδου(πεδίο χρόνου), πως θα πραγματοποιηθούν; Θα πάρουμε τον αντίστροφο μετασχηματισμό της συνάρτησης μεταφοράς για να είμαστε στο πεδίο του χρόνου; Η κάτι αλλο;

Μπορείς να ορίσεις την συνάρτηση μεταφοράς του φίλτρου με την εντολή tf και στη συνέχεια να κάνεις lsim με την συνάρτηση μεταφοράς και την είσοδο/έξοδο.

Παιδιά εγώ δεν μπορώ να εγγραφώ στο μάθημα στο elearning έχει κάποιος το ίδιο πρόβλημα ? Επιπλέον ποιος είναι ο κωδικός(αν έχει) ?

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

εδιτ: παρολαυτα η εργασια ειναι αυτη:

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

Το αρχείο v.p  με τι το ανοίγουμε και τι είναι ? τα δεδομένα ?

Συνάρτηση στο matlab είναι, απλά δε μπορείς να δεις τον κώδικα. Την καλείς όπως λέει στο pdf της εκφώνησης, απλά πρόσεχε γιατί είναι:
Vr = Vout(2), Vc = Vout(1)

Με lsim() μπορούμε να βρούμε τις τιμές της εξόδου y με βάση τις τιμές της εισόδου u ;

Η lsim προσομοιώνει την απόκριση ενός συστήματος σε μια είσοδο, οπότε αν της δώσεις όρισμα τη συνάρτηση μεταφοράς του φίλτρου και είτε την είσοδο είτε την έξοδο, θα φιλτράρει την είσοδο/έξοδο αντίστοιχα.

στο δεύτερο θέμα εχω την εντύπωση ότι θέλει δειγματοληψία σε time intervals μικρότερα από 0.0004sec (0.4 χιλιοστα δευτερολεπτου) λόγω της περιόδου της μικρής ταλάντωσης λόγω του LC, που είναι 0.001. Είναι κανείς σε ίδια φάση?

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

η αλήθεια είναι ότι γενικά που το τσέκαρα με μερικά άτομα μπορείς να φτάσεις 99%++ ακρίβεια εύκολα οπότε κάτι κάνεις λάθος :/

Το να βγουν οι εκτιμήσει κατά πολύ ελάχιστο πάνω από αυτά που μας δίνονται είναι οκ ή κάπου έγινε λάθος ;

ε κοίτα δεν βαθμολογώ εγώ αλλά σίγουρα πετυχαίνεται 99% ακρίβεια και πάνω οπότε πάντα θα μπορούσες να κάνεις "σωστότερη" υλοποιηση ισως αλλα μπορεί να μην αξίζει τον κόπο εν τέλει

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

Τελικά έβγαλε κανείς άκρη σχετικά με το πως προκύπτουν οι ΔΕ του κυκλώματος στην άσκηση 2;

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

εγω (και μερικοι ακομα :D) αναλυσαμε ως προς Vcdotdot(2η παραγωγος δηλαδη) και απλα μετα αφου ειχαμε την VC λυσαμε ως προς VR. (u1 + u2 = VC +VR)

Σε αυτό το σημείο βρίσκομαι και εγώ μετά απο την ανάλυση του κυκλώματος ..εχω μια ΔΕ δευτέρου βαθμού ως προς VC η οποία όπως περιλαμβάνει την 1η παράγωγο της VR  καθώς και την VR και V1. Από εδώ πως το προχώρησες; Δεν σε  έπιασα

Στο ημικύκλωμα με το Vc κάνεις το τύπο που αν είχες αντιστάσεις θα ήταν V_R2 = V * R₂ / (R₁ + R₂)

αξιοποιείς την αλλη εξισωση που εχεις

Ευχαριστώ για την βοήθεια παιδιά.. και κάτι τελευταίο  :P

Την γραμμική παραμετροποίηση στο δεύτερο σύστημα πως την κάνατε; Πρακτικά, είναι ένα σύστημα 2 εισόδων-2 εξόδων. Δεν διαφοροποιείται αυτή η περίπτωση σε σχέση με τα βήματα που ακουλοθήσαμε στην 1η περίπτωση(σύστημα 1 εισόδου-εξόδου);

Εφόσον έχουμε τη διαφορική εξίσωση δεύτερου βαθμού ως προς Vc, κάνουμε γραμμική παραμετροποίηση χρησιμοποιώντας αυτήν ;

Ναι, ακριβώς.

Στο δεύτερο θέμα ποιο χρονικό διάστημα θα πάρουμε για να χρησιμοποιήσουμε για παράδειγμα την lsim(..) ή μια ode(..) ; Παίρνουμε κάποιο τυχαίο που επιλέγουμε εμείς ;

Ρε παιδιά επειδή δεν την παλεύω άλλο,  έστειλα και μαιλ αλλα δεν πήρα σαφή απάντηση.. Μόλις φτάσουμε στην τελική ΔΕ του 2ου συστήματος, πως κάνουμε την γραμμική παραμετροποίηση; Αφού στην τελική εξίσωση υπάρχουν και οι 2 είσοδοι.;Επισυνάπτω και φωτο με την τελική ΔΕ

παραμετροποιησε γραμμικα το VC ... λυστο και θα σου βγει ενας πινακας θ* με 6 στοιχεια, που αφορουν τα  [VCdot VC u1dot u1 u2dot u2]

οταν βρεις το VC, χρησιμοποιησε την σχεση με τις εισοδου και βρες το VR

---

δεν ξερω αν χθες ησουν στο μαθημα του μεταπτυχιακου, αλλα μας εδωσε ενα hint και μας ειπε οτι το βημα πρεπει να ειναι καταλληλο για να πετυχουμε την ταλαντωση που κανει το VC κοντα στο 0. Ουσιαστικα try and error ειναι η φαση :Ρ

Τι εννοείς να βρω το Vc; Πρακτικά θέλω να φέρω το σύστημα σε μια μορφή y = θλ * ζ. Να το περάσω στον αλγόριθμο για να βρω τις παραμέτρους;
Το θ που καταλήγω είναι: θ = [ 1/RC 1/LC 1/RC 0 1/RC 1/LC]..

Βρες τις παραμέτρους με τον αλγόριθμο και κάνε προσομοίωση για να έχεις το Vc του μοντέλου. Μετά το Vr θα βγεί απ΄την άλλη εξίσωση απευθείας.

Παιδια τελικα εχει θεμα το αρχειο; Το εκανε κανεις το 2ο θεμα;

Ναι το έχω αντιληφθεί αυτό..πρακτικά, το πρόβλημα μου είναι όταν θα πάω να τρέξω τον αλγόριθμο για να βρώ τις παραμέτρους, το διάνυσμα ζ το οποίο πρέπει να σχηματίσω. Η είσοδος u που θα υπάρχει στο διάνυσμα ποια θα είναι; Το άθροισμα των δυο είσοδων u1+u2;

Το ζ θα έχει μια μορφή ζ = [ -s/L(s) *Vc -1/L(s)*Vc  s/L(s) *u 1/L(s)* u ]

Και τις δύο εισόδους θα έχει το διάνυσμα.

ζ = [ -s/L(s) *Vc -1/L(s)*Vc  s/L(s) *u1 1/L(s) *u1 s/L(s)*u2 1/L(s) * u2]

Ωραία ευχαριστώ πολύ..και φαντάζομαι προκύπτει ένα διάνυσμα 6 παραμέτρων , και παίρνουμε τις παραμέτρους για να λύσουμε την ΔΕ ως προς Vc; Μετά βρίσκουμε την Vr από την 2η σχέση και κάνουμε σύγκριση των τιμών των εξόδων που πρόεκυψαν με τις πραγματικές εξόδου για να δούμε πόσο κοντά στο πραγματικό σύστημα; Η κάτι κατάλαβα λάθος;

Ναι, αυτό που λες κάνεις.

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

Ρε παιδιά όταν πάμε να βρούμε τις τιμές του ζ πίνακα, για παράδειγμα την -s*y/Λ(s), με την lsim(..) δεν πρέπει να ορίσουμε ένα χρονικό διάστημα με ένα βήμα ; Παίρνουμε κάποιο συγκεκριμένο ή επιλέγουμε εμείς όποιο θέλουμε ;

Χρησιμοποιείς τον ίδιο χρόνο για τον οποίο πήρες τις μετρήσεις.