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 με προσοχή στα ονόματα των αρχείων!

[Αντιστοίχηση Μαθημάτων ΝΠΣ με ΠΠΣ]

Επίσημη ενημέρωση για Αντιστοίχηση Μαθημάτων ΝΠΣ με ΠΠΣ και η συζήτηση στο forum.

[png]

λοιπόν:
θέλω να χρησιμοποίσω την fzero για να βρω νουμερικώς τις ριζες μιάς τριγωνομετρικής εξίσωσης η οποία δε λύνεται αναλυτικά.
η fzero θελει function handler. εγώ έχω αυτό:

Code:

syms theta;
phi=pi/6;
for h=[0 0.1 0.2 0.5 1 2 5]
        y=(sin(theta))^2+ 2*h*cos(theta-phi);
        y1 = sin(theta)*cos(theta)+h*sin(theta-phi);
        y2 = cos(2*theta)+h*cos(theta-phi);
        ezplot(y,[0 2*pi]);

//κάπου εδω θέλω να βρω της ριζες της y1.

        %ezplot(y1,[0 2*pi]);
        %ezplot(y2,[0 2*pi]);
        %hold all;
end

τα y, y1, y2 είναι symbolic variables (απ'όσο κατάλαβα). τι πρέπει να τα κάνω ώστε μια fzero μεσα στη λούπα να μη μου γυρνάει error;
(πλιζ, όχι references στο documentation της mathworks, γιατι ειναι καλό μόνο άμα έχεις χρόνο και υπομονή)

και γενικά έχω μπλεχτεί άπειρα.

[Karaμazoβ]

να σε ρωτησω, γιατι δεν χρησιμοποιείς την solve(y,theta) ??

[png]

γιατι δε μπορεί να τη λύσει.

μου βγαζει ans= arg(Re(z1)+i*Im(z2))+2*pi*k το οποίο δε βοηθάει ιδιαίτερα.

[Karaμazoβ]

γιατι δε μπορεί να τη λύσει.

μου βγαζει ans= arg(Re(z1)+i*Im(z2))+2*pi*k το οποίο δε βοηθάει ιδιαίτερα.

εισαι σιγουρη?? απο το MatLab 2011b:

Code:

syms theta;
phi=pi/6;
h = 0.1; %paradeigma

        y=(sin(theta))^2+ 2*h*cos(theta-phi);
        y1 = sin(theta)*cos(theta)+h*sin(theta-phi);
        y2 = cos(2*theta)+h*cos(theta-phi);
        ezplot(y,[0 2*pi]);
>> 
>> 
>> 
>> solve(y1,theta)
 
ans =
 
0.046027347431528771078536140575958 + 1.8436655436388304855647095976386e-34*i
 
>> solve(y2,theta)
 
ans =
 
-0.7977431803860713605522195021279
 
>> solve(y,theta)
 
ans =
 
- 0.049895392467347344096156514398011 - 0.41934903901100004962886897809882*i
 
>> 

[png]

εχω το 2010α... λες να φταιει αυτο;

(επισης θελω οι τιμες του solve να ειναι πραγματικες)

[Eru lluvatar]

δοκίμασε fsolve.
πραγματική είναι η λύση του Karaμazoβ
το φανταστικό μέρος είναι σφάλμα μηχανής 1e-34

[Karaμazoβ]

δοκίμασε fsolve.
πραγματική είναι η λύση του Karaμazoβ
το φανταστικό μέρος είναι σφάλμα μηχανής 1e-34

το μυστηριώδες είναι οτι σκαλώνει στην y(theta) και βγαζει μιγαδικό, ενω στο διαγραμμα φαινεται οτι εχει δυο ριζες

[png]

*sigh*

σας ευχαριστώ.

[Eru lluvatar]

δοκίμασε fsolve.
πραγματική είναι η λύση του Karaμazoβ
το φανταστικό μέρος είναι σφάλμα μηχανής 1e-34

το μυστηριώδες είναι οτι σκαλώνει στην y(theta) και βγαζει μιγαδικό, ενω στο διαγραμμα φαινεται οτι εχει δυο ριζες

Το y(theta) δε το δα
Η solve είναι για πολυώνυμα αν δεν κάνω λάθος.

Code:

yy =@(theta) sin(theta).^2+ 2*h*cos(theta-phi);

fsolve(yy,[2 4])

Equation solved.

fsolve completed because the vector of function values is near zero
as measured by the default value of the function tolerance, and
the problem appears regular as measured by the gradient.

<stopping criteria details>


ans =

    2.7782    3.6048

[Karaμazoβ]

γαμηστερός ρε Eru !

η solve δεν ειναι μονο για πολυώνυμα, εχω λύσει και λογαριθμικες , αλλα επειδη ειναι περιοδική μαλλον κομπλαρει, και δε μπορείς να βαλεις πεδίο ορισμου να ψαξει...

[Eru lluvatar]

δε ξέρω δεν την έχω δουλέψει την solve γι αυτό είπα νομίζω 

Γενικά fsolve χρησιμοποιώ στην οποία μπορείς να περάσεις και πολλές παραμέτρους.
Αργεί καμιά φορά αλλά την κάνει την δουλεία της.

[Karaμazoβ]

δε ξέρω δεν την έχω δουλέψει την solve γι αυτό είπα νομίζω 

Γενικά fsolve χρησιμοποιώ στην οποία μπορείς να περάσεις και πολλές παραμέτρους.
Αργεί καμιά φορά αλλά την κάνει την δουλεία της.

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

[Decadence]

a = fi(v,s,w,F) returns a fixed-point object with value v, Signed property value s, word length w, best-precision fraction length, and embedded.fimath F.

μπορεί να μου εξηγήσεις κάποιος τι παριστάνει το "embedded.fimath F"

το έψαξα στο documentation, αλλά δεν έβγαλα άκρη..

[Decadence]

νομίζω έβγαλα μια άκρη.. μάλλον δεν έπρεπε να ασχοληθώ καν..

"Specifying a Fixed-Point Input Parameter by Example

To specify fixed-point inputs, you must install Fixed-Point Toolbox software.
To specify a signed fixed-point type with a word length of 8 bits, and a fraction
length of 3 bits, in the Define by Example dialog box, enter:
fi(10, 1, 8, 3)
HDL Coder sets the type of input to embedded.fi(1x1). By default, if you
have not specified a local fimath, HDL Coder uses the global fimath.
"

κάτι τέτοιο πρέπει να παριστάνει το "embedded.fimath F"

το συγκεκριμένο παράδειγμα αναφέρεται σε αυτήν την μορφή :

a = fi(v,s,w,f) returns a fixed-point object with value v, Signed property value s, word length w, and fraction length f.

αλλά λέει ότι :  HDL Coder sets the type of input to embedded.fi(1x1).

[Elade]

η συνάρτηση δέλτα στο matlab είναι dirac έτσι? γιατί προσπαθώ νατην χρησιμοποιήσω και μου λέει ότι είναι undefined