Matlab βοηθεια

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

[Πρόγραμμα Χειμερινής Εξεταστικής 2024-2025]

Πρόγραμμα Χειμερινής Εξεταστικής 2024-2025

[teslaaaa]

στην αρχη τα ειχα βαλει
αλλα και χωρις αρχικοποιηση τους,μια χαρα νουμερα βγαζει παλι!
θενξ για την τρολλο-απαντηση παντως!

[Randall Flagg]

function [aa,bb] = squarepar(e)
x = [0 1 4];
k=1;
while (x(1,3)-x(1,1))>e
    y = [ eq1(x(1,1)) eq1(x(1,2)) eq1(x(1,3))];
    d=polyfit(x,y,2);
    dd=polyder(d);
    w=roots(dd);
if w>x(1,2)
    if eq1(w)>=eq1(x(1,2))
        x(1,3)=w;
    else
        x(1,1)=x(1,2);
        x(1,2)=w;
    end;
elseif w<x(1,2)
     if eq1(w)>=eq1(x(1,2))
        x(1,1)=w;
     else
         x(1,3)=x(1,2);
         x(1,2)=w;
     end;
else
     if x(1,3)-x(1,2) < x(1,2)-x(1,1)
     w=x(1,2)-e/2;
     if eq1(w)>=eq1(x(1,2))
        x(1,1)=w;
        else
        x(1,3)=x(1,2);
        x(1,2)=w;
     end;
    else
        w=x(1,2) + e/2;
        if eq1(w)>=eq1(x(1,2))
        x(1,3)=w;
       else
        x(1,1)=x(1,2);
        x(1,2)=w;
        end;
    end;
end;
k=k+1;
end;
aa=x(1,1)
bb=x(1,3)
k=k
end


μου βγάζει συνέχεια το ίδιο αποτέλεσμα έχει κανείς ιδέα τι πάει λάθος;;;

[Merlin]

Έχω το εξής πρόβλημα. Έχω αυτή τη συνάρτηση:

Code:

function [ A ,B] = dixotomos2( func,a,b,l,e)
z=1;
A(z)=a;
A(z)=b;

while (b-a)>l
    x1=(a+b)/2 -e;
    x2=(a+b)/2 +e;
    if func(x1)<func(x2)
        A(z+1)=a;
        B(z+1)=x2;
    else
        A(z+1)=x1;
        B(z+1)=b;
    end
 z=z+1;  
    
end

end

και όταν πάω να την καλέσω πχ [a,b]=dixotomos2(@f1,0,4,0.01,0.001) έχοντας φτιάξει φυσικά την f1 κολάει και δεν γίνεται τίποτα. Προφανώς έχω λάθος στην συνάρτηση,αλλα δεν μπορώ να το βρω με τίποτα .Καμιά ιδέα;

Τσέκαρε τις τιμές των a και b,την συνθήκη στο while και βάλε ένα άνω όριο επαναλήψεων που να κάνει break την while,πχ με a=0,b=4 και l=0,01 τα x1,x2 => α,β θα αλλάζουν κατά +-e=0,01 οπότε ο βρόχος θα κάνει πολλές επαναλήψεις.Βάλε πχ

n = 100
while ()
....
if (k==n)
    break  
end
...
end

[teslaaaa]

Code:

function X = newton_stathero( x0,y0,e,gk,df,d2f )
%Είσοδος:(x0,y0)το αρχικό σημείο
%e->η σταθερά τερματισμού ε
%το βήμα γκ->εδώ σταθερό
%df->to grad tis f pou thelo na elaxistopoihso->apo cmw
%d2f->o essianos tis f os jacobian(df,v) opou v=[x,y]

k=1;
syms x y;
grad=subs(df,{x,y},{x0,y0});
%ypologismos arxikis timis gia na ksekinisei o algorithmos
%df diastaseis 1x2
%d2f diastaseis 2x2



while norm(grad)>e
    disp(norm(grad));
    grad=subs(df,{x,y},{x0,y0});%se kathe bima
    ess=subs(d2f,{x,y},{x0,y0});% 2x2
    i=inv(ess);%2x2
    dk=-i.*[grad(1);grad(2)];%(2x2)*(2x1)=(2x1)
    x0=x0+gk.*dk(1);
    y0=y0+gk.*dk(2);
    k=k+1;
end;
X(1)=x0;
X(2)=y0;
X(3)=k;


end

ενω προσπαθω να τα χω καλα με τις διαστασεις μου βγαζει αυτο:
??? Error using ==> times
Matrix dimensions must agree.

Error in ==> newton_stathero at 22
    dk=-i.*[grad(1);grad(2)];%(2x2)*(2x1)=(2x1)
μπορει κανεις να βρει τι μου διαφευγει?
ευχαριστω!

εδιτ:ακυρο,η τελεια εκανε ολη τη ζημια τλκα

[T--hmmy]

teslaaaa trolling

[Δον]

Xρειάζομαι και γω βοήθεια με το matlab.

Έχω τον κώδικα :

Code:

function m = LevenbergMarquart( f,x0,y0 )

xk=[x0, y0];
gk=0.5;
e=0.01;
k=1;

gradf=jacobian(f);
hessianf=jacobian(gradf);
gradf=matlabFunction(gradf);
d=eig(hessianf);
a=inv(hessianf);
hessianf=matlabFunction(hessianf);
a=matlabFunction(a);
d=matlabFunction(d);


while norm(gradf(xk(1),xk(2)),2)>=e
    
    
    mk=max(d(xk(1),xk(2)))+10*e;
    syms x y;
    dk=[x ; y];
    sol=(hessianf(xk(1),xk(2))+mk*eye(2))*dk+gradf(xk(1),xk(2))';
    [dx ,dy]=solve(sol(1),sol(2));
    
    
    xk=xk+gk*[dx dy]
    
    
    k=k+1;
    
    
    
end

m=xk;

end

Και μου βγάζει το error :

??? Undefined function or method 'norm' for input arguments of type 'sym'.

Error in ==> LevenbergMarquart at 18
while norm(gradf(xk(1),xk(2)),2)>=e
 
Ενώ το norm() στην 1η επανάληψη υπολογίζεται κανονικά, στη 2η με πετάει από τη συνάρτηση.

Εν τω μεταξύ δοκίμασα και "με το χέρι" στο command line του matlab να τρέξω τις ίδιες ακριβώς εντολές, με την ίδια ακριβώς σειρά και μου τρέχει κανονικότατα.

Έχει κανείς καμιά ιδέα τι παίζει  ???

[edenaxas]

γιατι δεν το τρεχεις σε debug να δεις τι παιζει αναλυτικα?

[Merlin]

Xρειάζομαι και γω βοήθεια με το matlab.

Έχω τον κώδικα :

Code:

function m = LevenbergMarquart( f,x0,y0 )

xk=[x0, y0];
gk=0.5;
e=0.01;
k=1;

gradf=jacobian(f);
hessianf=jacobian(gradf);
gradf=matlabFunction(gradf);
d=eig(hessianf);
a=inv(hessianf);
hessianf=matlabFunction(hessianf);
a=matlabFunction(a);
d=matlabFunction(d);


while norm(gradf(xk(1),xk(2)),2)>=e
    
    
    mk=max(d(xk(1),xk(2)))+10*e;
    syms x y;
    dk=[x ; y];
    sol=(hessianf(xk(1),xk(2))+mk*eye(2))*dk+gradf(xk(1),xk(2))';
    [dx ,dy]=solve(sol(1),sol(2));
    
    
    xk=xk+gk*[dx dy]
    
    
    k=k+1;
    
    
    
end

m=xk;

end

Και μου βγάζει το error :

??? Undefined function or method 'norm' for input arguments of type 'sym'.

Error in ==> LevenbergMarquart at 18
while norm(gradf(xk(1),xk(2)),2)>=e
 
Ενώ το norm() στην 1η επανάληψη υπολογίζεται κανονικά, στη 2η με πετάει από τη συνάρτηση.

Εν τω μεταξύ δοκίμασα και "με το χέρι" στο command line του matlab να τρέξω τις ίδιες ακριβώς εντολές, με την ίδια ακριβώς σειρά και μου τρέχει κανονικότατα.

Έχει κανείς καμιά ιδέα τι παίζει  ???

Η συνάρτησή σου f, είναι sym, δηλαδή expression, ενώ η Norm δέχεται σαν ορίσματα διανύσματα.Πρέπει με κάποιο τρόπο να περάσεις την κανονική συνάρτηση f μέσα στην συνάρτησή σου.Δηλαδή με κανονικό πεδίο ορισμού,ας πούμε
x = linspace(-1:1:400), αντί για syms x

[Δον]

Η συνάρτησή σου f, είναι sym, δηλαδή expression, ενώ η Norm δέχεται σαν ορίσματα διανύσματα.Πρέπει με κάποιο τρόπο να περάσεις την κανονική συνάρτηση f μέσα στην συνάρτησή σου.Δηλαδή με κανονικό πεδίο ορισμού,ας πούμε
x = linspace(-1:1:400), αντί για syms x

Μααα, μέσα στη norm δεν έχω την f αλλά την gradf συνάρτηση η οποία δεν είναι sym :S

[fugiFOX]

Δεν την έχω δουλέψει αλλά έτσι πως βλέπω τον κώδικα νομίζω και η gradf symbolic ειναι

[Δον]

Δεν την έχω δουλέψει αλλά έτσι πως βλέπω τον κώδικα νομίζω και η gradf symbolic ειναι

Έστω ότι είναι symbolic.
 
Στην 1η επανάληψη της while τρέχει κανονικά! Αυτό είναι που με τρελαίνει. 

Έχω αφήσει μια εντολή χωρίς ερωτηματικό για να βλέπω τι υπολογίζει και στην 1η επανάληψη την υπολογίζει .
Μετά μόλις πάει να υπολογίσει το νέο norm κρασάρει.
Και το δοκίμασα με τον ευκλείδιο τύπο sqrt(χ^2 + y^2) και μου λέει ότι το ge ( >= ) δεν υποστηρίζει και αυτό symbolic. Μα αφού βάζω τιμές(?)! 

(Ευχαριστώ για την ανταπόκριση πάντως )

[fugiFOX]

κάνε debugging by step και δες κάθε φορά τι input μπαίνει.
έτσι θα το βρεις το σφάλμα.

π.χ. όταν ο δείκτης του debugging είναι 1 βήμα πριν το while
τρέξε στο command window to gradf(xk(1),xk(2)),2
και δες τι δίνει ως αποτέλεσμα
(αυτό θα ειναι το input του norm() )

[edenaxas]

γιατι δεν το τρεχεις σε debug να δεις τι παιζει αναλυτικα?

[Δον]

γιατι δεν το τρεχεις σε debug να δεις τι παιζει αναλυτικα?

A καλά τώρα το είδα ( σόρρυ    )

Θα το δοκιμάσω

[christineL]

Πώς γίνεται όταν μια συνάρτηση επιστρέφει ως αποτέλεσμα ένα τεράστιο κλάσμα αυτό να το κάνω έναν απλό αριθμό; Δοκίμασα με FORMAT αλλά δεν γίνεται τίποτα,εκτός και αν το χρησιμοποιώ λάθος.