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

[O kanenas]

Ρε σεις, στο πρώτο μέρος της εργασίας, στην δεύτερη ερώτηση (που λέει για τις μαθόδους Secant και Muller) περιμένει να φτιάξουμε εμείς functions που να υλοποιούν αυτές τις μεθόδους?

By the way, ανέβασα την φετινή εργασία, η οποία είναι ένα κολάζ της 2ης και 3ης εργασίας του 2003

[sunday]

Την muller μας την δίνει ,τον κώδικα της secant τον έχουμε απο προηγούμενο εργαστήριο.
Εγώ για α=1,ορίζω το πολυώνυμο αλλά δεν έχει ρίζες...

[O kanenas]

Υπάρχει κανείς/μια να μου στείλει pm να του/της δώσω το msn μου ή το e-mail, να μου στείλει την secant, γιατί δεν την έχω?

[DeLarge]

Από το δεύτερο εργαστήριο, η secant:

function [x,xlist,iter] = secant(f,x0,x1,tol,maxit)
if nargin < 3, error('Insufficient input \n'); end;
if nargin < 4,tol = eps; end;
if nargin < 5, maxit = 50; end;
f0 = feval(f,x0);f1 = feval(f,x1);
iter = 0; xdiff = inf; xlist=[x0;x1]; x=x1;
while xdiff >= tol
xold=x;
x = x - f1*(x1-x0)/(f1-f0);
xdiff = abs(x-xold)/abs(x);
xlist=[xlist;x]; iter = iter + 1;
if iter >= maxit, error('Not converged after maxit iterations.'); end
x0=x1; f0=f1; x1=x; f1 = feval(f,x);
end

[O kanenas]

Σωστό! Δεν είχα σκεφτεί ότι θα μπορούσα να μάθω κι έτσι τον κώδικα της secant...
Ευχαριστώ!

[Ulm]

Βασικά τον κώδικα της τέμνουσας (secant) το βιβλίο την έχει κάπως διαφορετικά από τον παραπάνω.
Ας το τσεκάρει κάποιος γιατί δεν είμαι και σίγουρη. (σελίδα 255)

[eloxim]

Ο Φιλιππάκης ζήτησε το έκτο ερώτημα του τέταρτου εργαστηρίου... Το έχει λύση κανείς; Αν ναι,είναι εύκολο να το ανεβάσει, ή αν δεν θέλει να ανεβάσει έστω ένα σχεδιάγραμμα του πως λύνετε; Ευχαριστώ, προκαταβολικά!!!

[jova]

Καλησπέρα σε όλους.

Στην άσκηση με την μέθοδο Muller στο ερώτημα 4 τι εννοεί λέγοντας

Επιλέγοντας ως αρχικές τιμές άκρα διαστήματος που περιέχουν την αρχή των αξόνων?

Please Help!!!!

[remilt]

Ενοεί ακριβώς αυτό! Πάρε ώς άκρα για παράδειγμα τα -1 και 1 ενώ για την Muller πάρε και ένα εσωτερικό σημείο, λογικά η Muller θα βρεί μία από τις δύο συζηγείς μιγαδικές ρίζες ενώ όχι η Secant μιάς και είναι τέτοιος ο αλγόριθμος της μεθόδου που δεν μπορεί να περάσει από πραγματικές αρχικές τιμές σε μιγαδικές. άν θέσεις μιγαδικές αρχικές τιμές των οποίων τα πραγματικά μέρη να "εγκλείουν" την αρχή των αξώνων (πχ -1+i και 1+ i) μπορείς να εντοπίσεις τις μιγαδικές ρίζες και από τις δύο μεθόδους.

[jova]

thanx..!!!

Δεν ήμουνα σίγουρος αν έπρεπε να πάρω τυχαιες δικές μου.

[O kanenas]

Επειδή χρειάζεται για την εργασία, μπορεί κάποιος/α να κάνει ένα copy-paste τον κώδικα της linls, παρακαλώ πολύ?

[Ulm]

Ο κώδικας της linls:

function p=linls(x,y);
% linear least squears
n=length(x);
sx=x'*ones(n,1);
sy=y'*ones(n,1);
sx2=(x.^2)'*ones(n,1);
sxy=(x.*y)'*ones(n,1);
par=n*sx2-sx^2;
a=(n*sxy-sx*sy)/par;
b=(sx2*sy-sxy*sx)/par;
p=[a;b];

[O kanenas]

Merci beaucoup!

[fael_adan]

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

[armaos]

Παιδιά παρακαλώ ας ενημερώσει κάποιος πότε θα εξεταστεί η 2η εργασία γιατί είμαι παλιός και δεν παρακολουθώ.