THMMY.gr

Το μαθημα διδασκεται φετος?

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

ειναι ευκολο?

γιατι μου ειπανε οτι και η ιδια ρωταει για να καταλαβει....

ρωτοντας πας στην πολη

σαε 3 διδασκεις ομως??

η πλακα ειναι οτι το δηλωσα χωρις να ξερω....

λολ αράσι


Η δουλγέρη έχει διαφορετικό τρόπο σκέψης από μένα και είναι καλό για μένα να τον παρακολουθώ. Γενικά
την πρώτη φορά σου φαίνεται δυσνόητη η νοοτροπία της αλλά μετά γίνεται challenge.

sos matlab!

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

το παράδειγμα του παραρτήματος είναι το εξής:

function dx=odefun(t,x)
mi=1;
dx=[x(2);mi*(1-x(1)^ 2)*x(2)-x(1)];
[t,x]=ode45(odefun,[0 20],[2;0]);

γράφω σε διαφορετικό σκριπτ την function και καλώ την ode αλλα τιποτες. τα βάζω όλα σε ένα πάλι τίποτες. άλλες φορές λέει ότι το χ ειναι αγνωστη αλλες φορες το t ειναι άγνωστο. γενικά πως λειτουργεί το σύστημα? αν γνωρίζει κάποιος από ματλαμπ και όχι απαραίτητα σαε ας απαντήσει  :-\

βασικά δοκίμασε το εξής παράδειγμα στη σελίδα 420-421 του βιβλίου ασκήσεων των Κυκλωμάτων 3


Γράψε ένα script στο MATLAB με όνομα equation.m

function xprim=equation(~,x)
xprim=-x^2;

Και άλλο 1 με όνομα cad10.m

[t,x]=ode45('equation',0,10,1);
plot(t,x,t,x,'o')
xlabel('time(sec)')
ylabel('x')
grid

Mετά τρέξε το cad10.m

p.s.έχω διορθώσει το 'equ1' me 'equation' γιατί αλλιώς δεν τρέχει!!!(έχει λάθος!)
Φαντάζομαι αν τρέξεις αυτό θα βρεις λύση κ στο πρόβλημα σου! :)

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

Πάντως το παραπάνω παράδειγμα του Μάργαρη μόλις το δοκίμασα σε Matlab 2011a k δούλεψε!

Μία είναι η υπερλύση για όλα τα προβλήματα στην εργασία σας

http://www.google.de/search?q=pplane8&ie=utf-8&oe=utf-8&aq

η συναρτηση προσήμου ορίζεται καπως στη matlab? εκανα search χθες και ειδα ότι υπαρχει μια συναρτηση sign αλλα επιστρεφει πινακα με  στοιχεια -1 0 +1 αν τα στοιχεια των δυο πινακων που συγκρινονται ειναι<= = >=. στο β μερος της εργασιας το πορτραιτο και οι χρονικες αποκρισεις βγαινουν ότι να'ναι παρολο που η δομή του κωδικα μου δεν διαφερει πολυ απο εκεινη του γραμμικου (αν εξαιρεσεις μια συνθηκη if else if για να ελεγξεις σε ποια διαστηματα τιμων κινειται το ε).... υπαρχει περιπτωση καποιος από τα ελαχιστα ατομα που ειμαστε στα σαε να το εκανε?  :)

Εδώ έχουμε ένα μικρό τεχνικό θέμα.
Ζητάει ρητά να χρησιμοποιήσουμε τις odefun Και ode45  :-\

Eμένα μου τρέχει πάντως η ode με τον τρόπο
τώρα που το σκέφτομαι, ελπίζω χριστίνα να μην έβαλες τη γραμμή
[t,x]=ode45(odefun,[0 20],[2;0]);
μέσα στο odefun.m ? (επίσης προσοχή με τον πίνακα αρχικών τιμών στις ode, μπορεί να υπάρχει πρόβλημα με τις διαστάσεις και θέλει ανάστροφο ή όχι ανάλογα πώς το έχετε ορίσει)

ναι, αυτή είναι η συνάρτηση προσήμου
αν και δεν ξέρω αν τη χρειάζεσαι ακόμα

η δουλγέρη είπε ότι δεν έχει κανένα πρόβλημα με την πλατφόρμα, αρκεί να το υλοποιήσεις το πρόβλημα

προτείνει αυτές τις συγκεκριμένες συναρτήσεις γιατί η μέθοδος της αριθμητικής προσέγγισής τους είναι ακριβής για τις διαφορικές που αντιμετωπίζουμε
αν θέλεις περισσότερα περί αριθμητικών μεθόδων στις διαφορικές εξισώσεις, έχει λεπτομέρειες σε ένα pdf για το matlab στο ethmmy - σαε3

λοιπόν στο ετημμυ αναρτήθηκε ολόκληρη η εργασία. από μια ματιά που της έριξα είναι λίγο μεγαλούτσικη (από δυσκολία ακόμη δεν μπορώ να πω κάτι). την τρίτη σκοπεύω να πω στο μάθημα την κ.Δουλγέρη αν μπορεί η παράδοσή τους να γίνει αμέσως μετά τις διακοπές των Χριστουγέννων. Δηλαδή κατά τις 8-10/1. όσοι είμαστε στο 7ο εξάμηνο για τις επόμενες τρεις βδομάδες θα έχουμε εργαστήριο μίκρο1 στα οποία πρέπει να παραδώσουμε εργασίες. οπότε ο χρόνος μας πιέζει πολύ για να ασχοληθούμε και με τα σαε. όσοι συμφωνείτε με την πρότασή μου θα σας παρακαλούσα την τρίτη να την στηρίξετε για να μην δοθεί η εντύπωση ότι απλά το προτείνω γιατί βαριόμαστε να ασχοληθούμε (του εναντίων μάλιστα μακάρι να ασχολούμασταν μόνο με σαε).

Πρέπει να παραδώσουμε την εργασία όλη στο τέλος??
Δε γίνεται να την παραδίδουμε σε κομμάτια? :/
πχ. το πρώτο τώρα, το δευτερο το Δεκεμβρη, το τρίτο μετα τα Χριστουγεννα??
Όχι τιποτε αλλο, θα βοηθησει στην οργάνωση...

Ζήτησα σήμερα από την κ.Δουλγέρη αν γίνεται να παραδώσουμε τις εργασίες μετά τα χριστούγεννα. Εκείνη δέχτηκε και μας όρισε τελική ημερομηνία παράδοσης στις 16/1/2012.

Επίσης την επόμενη βδομάδα θα λείπει και θα κάνουμε με τον Λεωνίδα 3 ώρες μόνο ασκήσεις πάνω στην ύλη σχετικά με τη γραμμικοποίηση συστημάτων αλλά και πιο σύνθετες που θα συνδυάζουν και προηγούμενα κεφάλαια της ύλης.

γεια σου ρε κριστινέτ

Έχει καταλάβει κάποιος τη λύση της άσκησης 23 όπως παρουσιάστηκε στο μάθημα στις 22/11/2011 ???   :-\

Κάπου το χάνω εκεί που κάνει τις αντικαταστάσεις του Α στη V' 

Έστω ότι έχω ένα διάνυσμα θέσης q=[q1 q2]^T και ένα διάνυσμα ταχύτητας  q'=[q1' q2']^T. αυτά συνδέονται με μια διαφορική εξίσωση της μορφής Μ(q)*q''+C(q, q')*q'=u όπου Μ, C πίνακες 2χ2 και u η είσοδος του συστήματος.

γενικότερα για να λύσω συστήματα διαφορικών εξισώσεων χρησιμοποιώ την ode45 αλλά εκεί απλά όριζα δύο παραμέτρους χ1=κάτι και χ2=κάτι άλλο (ότι έδινε τελος πάντων το σύστημα). Αν τώρα ορίσω χ1=q και χ2=q' το διάνυσμα χ από εκεί που ήταν 2χ1 όπως συνηθίσαμε στις προηγουμενες εργασίες γίνεται 4χ1! πως στο δία γίνεται η εισαγωγή των q, q' στο χ έτσι ώστε να βρούμε τα ζητούμενα?

εννοείς ότι το πρόβλημα είναι με πίνακες

δεν έχει διαφορά με πριν, γράψε μία μια τις εξισώσεις και λύσε τες με ode45

να κανω ας πουμε χ1=q1, χ2=q2, x3=q1', x4=q2' ?

δεν το γραψα καλά
και δεν έχω κάνει το ερώτημα αλλά
αν έχω σύστημα π.χ. 4 δ.ε. (4 μεταβλητές κατάστασης - είναι φασικές, έτσι; ) φτιάχνω αρχείο στυλ

function df=odefun3(t,x)
df=[μπλα μπλα 1; %x1dot=q1
μπλα μπλα 2; %x2dot=q1'
μπλα μπλα 3; %x3dot=q2
μπλα μπλα 4; %x4dot=q2'
];

και μετά γράφω την ode μου κανονικότατα
[t,x]=ode45(@odefun3,[0 10],x0);
με
x0=[α;β;γ;δ];
και ό τι επιπλέον αρχικές τιμές έχει τις γράφω στην ίδια σειρά

Αναρτήθηκαν οι βαθμοί των εργασιών στο  eTHMMY !!!

παιδιά ξέρουμε αν δίνουμε με ανοιχτές τις σημειώσεις της δουλγέρη?

Μόνο το τυπολόγιο επιτρέπεται!!!

ποιο τυπολόγιο?


edit: α, τώρα το είδα
στο ετημμυ, υλικό μαθήματος

το μάθημα μπορούν να το δηλώσουν ενεργειακοί σαν επιλεγόμενο?

σαν ελεύθερης νομίζω ναι.
τον μπελά σου ψάχνεις...?  :???:

δεν αντεχω χωρις σαε λαιμαι ;D

κλαμα ελευθερης επιλογης

μεγαλες στιγμες