Προγραμματισμός με Εφαρμογές στην Επιστήμη του Μηχανικού

[alexk]

Στο θέμα 4 πως βρίσκουμε σημεία ισορροπίας που ζητάει? Εμένα και στα πέντε σημεία που υπολόγισε, η παράγωγος είναι <1.. είναι ευσταθή σημεία?? Φτάνει αυτό??

εεεε..??

Φαντάζομαι αν κανεις plot τα i - Xi σου βγαζει μια γ.π η οποια για y=0.948 περίπου είναι το σημειο ισορροπιας το οποιο ειναι ευσταθες γιατι η απολυτη τιμη της παραγωγου της sin (r * x) είναι μικροτερη της μοναδος..


Επίσης εγώ έχω ερώτηση πως βρίσκουμε περίοδο πλάτος μεγιστη ταχυτητα στο θεμα 2γ
και τι ακριβως εννοει με τον ορο συγκλιση στο θεμα 5, β ερωτημα

[vampinthe]

..μα σε όλα τα σημεία του γραφήματος είναι <1 η απόλυτη τιμή της παραγώγου.. δεν είναι??

όσο για την περίοδο στο δευτερο θέμα, μετράς την οριζόντια απόσταση δυο κόμβων ή δυο κοιλίων στο x(t) και για το πλάτος, την μέγιστη απόσταση από τον άξονα x . Πάτα set grid στο gnuplot πριν πλοτάρεις για να το δεις πιο καθαρά. Την μεγιστη ταχύτητα απτο v(t).

για το 5ο θέμα αϊ ρο νόου..

[jimmakos]

ego thelo na rotiso ton vampinthe an ta kserei kala ta themata gt kati mou leei oti dne exei diavasei kala. ta pine tin paraskeui!
XAXA
an me vgalete off-topic tha plakothoume. einai sxetiko me to mathima to minima mou!

[Tsakalos]

Λοιπον, απο οτι καταλαβα σχετικα με το δευτερο θεμα που λεγαμε πριν.. ο citizen και ο chris κανουν το ιδιο πραγμα.. μια αλλαγη μονο στον ηδη υπαρχοντα κωδικα της euler. To διαγραμμα που προκυπτε για τα δεδομενα της ασκησης ειναι καπως ετσι http://img195.imageshack.us/f/euler.png/

Αυτο που προτεινει ο alexk νομιζω οτι ειναι αυτο που μας εκανε και ο καθηγητης στο εργαστηριο και προκυπτει αυτο το διαγραμμα http://img580.imageshack.us/f/alloskodikas.png/

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

EDIT: τελικα βρηκα τι φταιει.. η μεθοδος euler επειδη δεν ειναι ακριβης χρειαζεται πολλα βηματα για να βρει τη λυση. Το διαγραμμα της euler αν αυξηθουν αρκετα τα βηματα γινεται τελικα ιδιο με το αλλο διαγραμμα...

[Chris]

Εμένα μου βγαίνει το ίδιο με το alloskodikas.jpg

Προφανώς και είναι σωστό αφού το έκανα ΕΓΩ.

EDIT: τελικα βρηκα τι φταιει.. η μεθοδος euler επειδη δεν ειναι ακριβης χρειαζεται πολλα βηματα για να βρει τη λυση. Το διαγραμμα της euler αν αυξηθουν αρκετα τα βηματα γινεται τελικα ιδιο με το αλλο διαγραμμα...

Για αυτό λέει στην εκφώνηση να παίξουμε με τα steps.

[apolski]

Βγηκαν οι κατανομες http://mycourses.ntua.gr/courses/SEMFE1 ... ne2010.pdf

Ομως δεν εχει το ονομα μου ενω γραφτηκα... Τι να κανω?

[alexk]

Λοιπον, απο οτι καταλαβα σχετικα με το δευτερο θεμα που λεγαμε πριν.. ο citizen και ο chris κανουν το ιδιο πραγμα.. μια αλλαγη μονο στον ηδη υπαρχοντα κωδικα της euler. To διαγραμμα που προκυπτε για τα δεδομενα της ασκησης ειναι καπως ετσι http://img195.imageshack.us/f/euler.png/

Αυτο που προτεινει ο alexk νομιζω οτι ειναι αυτο που μας εκανε και ο καθηγητης στο εργαστηριο και προκυπτει αυτο το διαγραμμα http://img580.imageshack.us/f/alloskodikas.png/

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

EDIT: τελικα βρηκα τι φταιει.. η μεθοδος euler επειδη δεν ειναι ακριβης χρειαζεται πολλα βηματα για να βρει τη λυση. Το διαγραμμα της euler αν αυξηθουν αρκετα τα βηματα γινεται τελικα ιδιο με το αλλο διαγραμμα...

Κοιταξε με mail που εστειλα στον Αναγνωστοπουλο μου απαντησε την εντελως απλη μορφη. Αυτη που προτεινα.
Μου το παρουσιασε σαν πιο απλη λύση. Βεβαια πιο απλο νομιζω οτι ειναι να πειραξεις μονο την επιταχυνση στον κωδικα της euler.f
Θεωρω οτι και τα δυο σωστα θα τα πιασει αν κρινω απο την εκφωνηση η οποια λεει "Να μειώσετε το βημα ολοκληρωσης dt μεχρι τα αποτελσματα σας να σταθεροποιηθουν (εναλλακτικα αυξηστε τον αριθμο βηματων steps"
Βασικα τωρα που το ξαναβλεπω σίγουρα θα πιασει και τις δυο απαντησεις σωστες!

[citizen]

Ο πασιπουλαρίδης πάντως στο μάθημα άλλαξε μόνο την συνάρτηση αccel στο euler.f

[zoe]

Ξέρει κανείς αν αύριο θα μπορούμε να έχουμε μαζί μας τα βιβλία που μας έχουν δώσει??

[alexk]

Ξέρει κανείς αν αύριο θα μπορούμε να έχουμε μαζί μας τα βιβλία που μας έχουν δώσει??

Ο,τι θέλουμε μπορούμε να έχουμε..

[diosak]

ποσα θεματα πρεπει να λυσουμε??

[dj manos]

μας δινουν 1 θεμα & το λυνουμε.

[alexk]

Έχει λύσει κανείς το θέμα 6?

[sunday]

Εγώ αυτό έγραψα για το 6 τώρα το αν ειναι σωστό παίζεται..Το έκανα βάσει του θ.4

Κώδικας: Επιλογή όλων

 program Newton_Raphson
     
      implicit none
      integer P,i
      parameter(P=50)
      real  x(P)
      real x0,e,ei
      real f,fd
       

      print *,"# Enter x0:"
      read(5,*)x0
      print *,"# Enter the error:"
      read(5,*)e
       
      if (e .le. 0.0) stop 'illegal value of error'
      i=2
      x(1)=x0 

      open(unit=20,file="newton_r2.dat")
      write(20,*)1,x(1)

      do while(i .le. P)
      x(i)=x(i-1)-f(x(i-1))/fd(x(i-1))
      ei=abs(x(i)-x(i-1))/x(i)
      print *,i,x(i)
      write(20,*)i,x(i),ei
      if (ei .le. e) stop 'Procedure completed successfully'
      i=i+1
      enddo
      print *,'# Method failed after P=',P,'iterations'
      close(20)
      end
     
      real function f(x)
      implicit none
      real x
      f=sqrt(2.0)*x*(1-x)-x
      return
      end
     
      real function fd(x)
      implicit none
      real x
      fd=sqrt(2.0)+2*sqrt(2.0)*x-1
      return
      end
      

[Tsakalos]

Εκανα το πρωτο ερωτημα απο το θεμα 7 και μου φαινεται οτι παιζει σωστα. To εγραψα απο κατω αν θελετε να το τσεκαρετε.
Στην τροχια μου βγαζει τα γνωστα σχηματακια λισαζου!
Ξερετε πως θα υπολογισουμε την ενεργεια? Λογικα για την κινητικη ενεργεια θα βρουμε τη συνισταμενη ταχυτητα αφου οι ταχυτητες ειναι καθετες v=ριζα(vx^2+vy^2) σωστα?
Η δυναμικη ενεργεια ομως πως βγαινει??
Τωρα θα δοκιμασω και το θεμα 6 ενω ακομα δεν εχω κανει καμια ασκηση με τη βλακεια τη λογιστικη εξισωση

Κώδικας: Επιλογή όλων

      program kinisi
      implicit none
      real*8 t,tf,dt,x0,y0,w1,w2,vx,vy,ax,ay,x,y
c------INPUT
      print *, 'enter x0, y0, w1, w2, tf, dt'
      read(5,*) x0, y0, w1, w2,tf,dt
      print *, 'ok malaka'
c------INITIALIZE
      x=x0
      y=y0
      t=0
      open(unit=11,file='kinisi.dat')
c------COMPUTE
      do while (t .le. tf)
         x = x0*cos(w1*t)
         vx=(x0/w1)*sin(w1*t)
         ax=(-x0/(w1*w1))*cos(w1*t)
         y=y0*sin(w2*t)
         vy=(-y0/w2)*cos(w2*t)
         ay=(-y0/(w2*w2))*sin(w2*t)                        
         t=t+dt
         write(11,*)t,x,y,vx,vy,ax,ay
      enddo
      close(11)
      end