Αριθμητική Ανάλυση II και Εργαστήριο

[1/2rizax]

Και οι 2 κώδικες φαίνονται σωστοί, απλά ο 2ος θα μπορούσε να ναι συντομότερος. πχ το tout1v, που δεν είναι άλλο από το tout1 χωρίς το τελευταίο του σημείο, το ορίζει, αλλά δεν το χρησιμοποιεί πουθενά. Επίσης το while k<length(err1)-1 δεν αφαιρεί 2 τιμές αντί για μία; Ίσως να θέλει <=

Βγήκαν οι ημερομηνίες για τις εργασίες; Πρέπει να τις παραδώσουμε τισ συγκεκριμένες ημερομηνίες ή μπορούμε και από τώρα για να ξεμπερδεύουμε; Έχουμε και εξαταστική...

[armaos]

παιδια ευχαριστω πολυ για την απαντηση.
Μπορειτε να μου εξηγήσετε κατι;
1)το OSF=[]; ειναι δηλωση διανυσματος χωρις να δηλωνει το μεγεθος του;
2)Η γραμμη OS=[OS;abs(f0true(tout(length(tout)))-yout(length(yout)))]; τι εκφραζει; μονο το abs καταλαβαινω! (ξερω οτι γινομαι κουραστικος, αν θελετε μου απαντατε)

Ημερομηνια για παραδοση των εργασιων δεν εχει βγει.Οταν ειχα μιλησει με την κοπελα που μας κανει το εργαστηριο μου ειχε πει παντως οτι και πριν μπορουμε να τις παραδοσουμε.Κανονικα παντως ειτε πεμπτη ειτε παρασκευη θα ειναι η εξεταση.Λογικα παρασκευη

[1/2rizax]

1)Σωστά
2)length(tout) = το μήκος του tout (ή αλλιώς η θέση του τελευταίου σημείου της διαμέρισης)
tout(length(tout)) = η τιμή του τελευταίου σημείου της διαμέρισης
Η υπόλοιπη έκφραση abs (...) είναι η γνωστή έκφραση για το σφάλμα. Έτσι υπολόγισες το σφάλμα στο τελευταίο σημείο της διαμέρισης. Επειδή η μέθοδος έχει σταθερό βήμα η έκφραση tout(length(tout)) θα είναι ίση με το δεξί άκρο του διαστήματος, οπότε μπορείς να την αντικαταστήσεις κατευθείαν με την τιμή αυτή.
Με το υπόλοιπο OS=[OS; απλά αποθηκεύεις αυτήν την τιμή του σφάλματος στο ήδη υπάρχον διάνυσμα OS προσθέτοντάς του μια καινούρια γραμμή.

[armaos]

θεος/-α!!! Ευχαριστω παρα πολυ! ναου ιτ μεικς σενς!!

βασικά πανω σε αυτο που ειπες πιο πανω. Την εργασια την εχω σχεδον τελειωσει αλλα μου βγαζει οτι τα διανυσματα δεν εχουν ιδιο μεγεθος, με αποτελεσμα να μην εχω κκαι γραφικη σφαλματος. Ο κωδικας που εγραψα ειναι παρομοιος με του παιδιου.

clf
hold off
clear all
format long

t0=1;
tfinal=4;
y1=0;
tol=1.e-4;

[tout,yout]=frk45a('p1',t0,tfinal,y1,tol)
err=abs(yout-p1true(tout))
d=diff(tout);
k=1;
v=[];
toutv=[];

while k<length(err)-1
toutv=[toutv;tout(k)]
k=k+1;
end

while k<length(err)-1
v=[v;err(k)/(d(k))^5]
k=k+1;
end

subplot(2,1,1);plot(tout,yout,'*',tout,p1true(tout),'-');
subplot(2,1,2);plot(toutv,v,'-');

Πιστευεις μπορεις να με βοηθησεις καθολου με αυτο;

[O kanenas]

Δεν είδα την ανακοίνωση, αλλά πληροφορήθηκα ότι το τμήμα της Τετάρτης 2:00-3:30 αναπληρώνει αύριο στις 11 εργαστήριο.

[Opeth]

και των 3.30-5 μετα? αν ξερει κανεις να απαντησει plz:!

[1/2rizax]

Η ανακοίνωση δεν διαχώριζε τα τμήματα. Όλοι μαζί θα εννοεί στις 11.

@αρμάος
Η μεταβλητή σου k φτάνει μεχρι τη μέγιστη τιμή και στο επόμενο λουπ ξεκινά από κει με αποτέλεσμα το διάνυσμα ν να είναι κενό! Βάλε διαφορετική μεταβλητή. πχ:

clf
hold off
clear all
format long
tol=10^(-4);

[tout,yout]=frk45('r1',1,4,0,tol)
err=abs(yout-r1true(tout))
d=diff(tout);
k=1;
m=1;
v=[];
toutv=[];

while k<length(err)-1
toutv=[toutv;tout(k)]
k=k+1;
end

while m<length(err)-1
v=[v;err(m)/(d(m))^5]
m=m+1;
end
subplot(2,1,1);plot(tout,yout,'*',tout,r1true(tout),'-');
subplot(2,1,2);plot(toutv,v,'-');

[armaos]

ok eyxaristv poly!!

[garfield]

Θα μπορουσε καποιος να ενημερωσει σχετικα με το τι ειπωθηκε σημερα στο εργαστηριο για τα ακολουθα θεματα?

->ποτε ειναι η καταληκτικη ημερομηνια παραδοσης των εργασιων?
->θα γινει η προφορικη εξεταση τους?
->αν ναι,ποτε θα συμβει αυτο ή πώς θα ενημερωθουμε?
->ποτε θα γινει η γραπτη εξεταση του εργαστηριου?πώς θα ενημερωθουμε και γι αυτην?


αληθεια,η υπευθυνη γιατι δεν απαντα στα e-mails?

[ran-tan-plan]

Μα καλα,κανεις δεν πηγε σ αυτην την αναπληρωση?Ειναι πολυ σημαντικο να απαντησει καποιος,ειναι φλεγοντα τα ερωτηματα που θετει ο garfield!!!

off-topic:μεγαλη ταλαιπωρια φετος με το εργαστηριο..

[de profundis]

Λοιπόν,στο τελευταίο εργαστήριο που τη ρωτήσαμε για τη γραπτή εξέταση του εργαστηρίου είπε πως θα γίνει μες στο φλεβάρη(κατατοπιστικότατο) και θα αναρτηθεί ανακοίνωση για την ακριβή ημερομηνία στον πίνακα ανακοινώσεων έξω απ το pc lab.Τη ρώτησα αν θα βγάλουν ανακοίνωση κάπου στο internet,γιατί είναι περίοδος εξεταστικής και δε μπορούμε να πηγαίνουμε κάθε μέρα,και μου είπε αν θέλουμε να της στείλουμε mail(ελπίζω τότε να απαντάει).Οσο για την προφορική εξέταση είναι υπεύθυνος ο Φιλλιπάκης και δεν έχει βγάλει ημερομηνίες ακόμα...Για τις εργασίες απ ό,τι άκουσα η παράδοση πρέπει να γίνει πριν τη γραπτή εξέταση του μαθήματος αλλά το λέω με επιφύλαξη γιατί το είχε πει στο προηγούμενο μάθημα που δεν ήμουν εκεί και μη σε πάρω στο λαιμό μου..Τέλος,για το μάθημα που έγινε την προηγούμενη εβδομάδα(μετά την ακύρωση της γσ) πήραμε απουσίες όσοι λείπαμε αλλα γενικά και δύο απουσίες να έχεις δεν υπάρχει πρόβλημα..

[armaos]

να ρωτήσω κάτι παιδιά, στην 2η εργασία στο 2ο ερώτημα της 5ης ασκησης μου βγάζει λαθος στο μέγεθος του διανύσματος.Το προγραμμα που εχω γραψει ειναι αυτο
clear all;
clf;
format long e;

t0=0;
tfinal=10;
x0=1;
y0=3;
tol1=10^(-1);
tol2=10^(-2);
tol3=10^(-4);
tol4=10^(-5);
tol5=10^(-6);

[tout1,yout1]=nrk45('fy',t0,tfinal,y0,tol1);
[tout2,yout2]=nrk45('fy',t0,tfinal,y0,tol2);
[tout3,yout3]=nrk45('fy',t0,tfinal,y0,tol3);
[tout4,yout4]=nrk45('fy',t0,tfinal,y0,tol4);
[tout5,yout5]=nrk45('fy',t0,tfinal,y0,tol5);

G1=exp(yout1(:,1)+2*yout1(:,2))./(yout1(:,1).*(yout1(:,2).^2));
G2=exp(yout2(:,1)+2*yout2(:,2))./(yout2(:,1).*(yout2(:,2).^2));
G3=exp(yout3(:,1)+2*yout3(:,2))./(yout3(:,1).*(yout3(:,2).^2));
G4=exp(yout4(:,1)+2*yout4(:,2))./(yout4(:,1).*(yout4(:,2).^2));
G5=exp(yout5(:,1)+2*yout5(:,2))./(yout5(:,1).*(yout5(:,2).^2));

subplot(5,2,1);plot(tout1,G1,'r')
subplot(5,2,3);plot(tout2,G2,'b')
subplot(5,2,5);plot(tout3,G3,'g')
subplot(5,2,7);plot(tout4,G4,'y')
subplot(5,2,9);plot(tout5,G5,'v')
subplot(5,2,2);plot(yout1(:,1),yout1(:,2))
subplot(5,2,4);plot(yout2(:,1),yout2(:,2))
subplot(5,2,6);plot(yout3(:,1),yout3(:,2))
subplot(5,2,8);plot(yout4(:,1),yout4(:,2))
subplot(5,2,10);plot(yout5(:,1),yout5(:,2))


Στο πρόγραμμα που εχω γράψει μου βγαζει ενα ακυρο παραθυρο γραφικης παραστασης. που οφειλεται;

[1/2rizax]

(yout1(:,1)+2*yout1(:,2)
Όλα τα yout είναι πίνακες-στήλες, οπότε δε γίνεται να έχουν και 2η στήλη. Το πρόβλημα είναι το yout1(:,2)
Οπότε δοκίμασε να αντικαταστήσεις πχ: την
G1=exp(yout1(:,1)+2*yout1(:,2))./(yout1(:,1).*(yout1(:,2).^2));
με την
G1=exp(tout1(:,1)+2*yout1(:,1))./(tout1(:,1).*(yout1(:,1).^2));
και ομοίως για τις υπόλοιπες

εντιτ:
γράψε άκυρο. πες λίγο την fy σου γιατί μάλλον δε βάλαμε την ίδια.

[armaos]

η fy Μου ειναι η
function yprime=fy(t,y);
yprime=[2*(y(1)-y(1)*y(2)); -(y(2)-y(1)*y(2))];

Μηπως επρεπε να ειναι
yprime=[2*(y(1)-y(1)*y(2)): -(y(2)-y(1)*y(2))];

[armaos]

λοιπον παιδια βγηκε ανακοινωση οτι η γραπτη εξεταση του εργαστηριου θα γινει την τριτη 9 φεβρουαριου 12 με 3. Την ιδια μερα κατα πασα πιθανοτητα θα εξεταστουν κ οι εργσιες προφορικα. Δεν ξερω για σας αλλα εγω εκεινη τη στιγμη θα δινω αναλυση παλινδρομησης, λα'ι'φ ιζ φουλ οφ πισιμπιλιτιζ, αλλα δεν νομιζω οτι ειναι λυση εκεινη τη μερα. Για να μην το λεμε λιγο πριν, εγω εστειλα μειλ κ ενημερωσα την κοπελα που μας εκανε εργαστηριο για το θεμα.Αν μαθω κατι νεοτερο θα ποσταρω