THMMY.gr

Ότι απορίες έχετε ???
πάνω στο
δημοφιλέστερο και πληρέστερο
πακέτο Μηχανικών
εδώ είμαστε να σας τις λύσουμε 8)

---

Edited by moderator: Αλλαγή τίτλου

Πού θα βρω ένα καλό βιβλίο για το matlab
ΤΖΑΜΠΑ; Ή τουλάχιστον φτηνά; ::)

το καλύτερο βιβλίο του MATLAΒ είναι τσάμπα
και είναι ήδη μέσα στον Η/Υ σου  ;)
Πρόκειται γισ το HELP .
Μπορώ να σε πω ότι έχω διαβάσει 3 βιβλία πάνω στο MATLAB
αλλά έμαθα από το Help με τη γνωστή μέθοδο του
by trial and error
Πάντως ένα καλό βιβλίο για αρχή
με όλα τα βασικά που χρειάζεσαι βρίσκεται
στο βιβλιαράκι του κ. Μάργαρη εργαστήριο Ηλεκτρικών Κυκλωμάτων
που μας διανέμεται μαζί με το βασικό βιβλίο.
Εάν είσαι αρχάριος τότε
όπου βλέπεις παχύ βιβλίο... μακριά!
Ξεκίνα και ότι απορία έχεις εδώ είμαστε

Χαίρεται.
Και γω δεν χρησιμοποιούσα το ΜΑΤLΑΒ μέχρι προσφάτως, οπότε και ο κύριος Μάργαρης μας ανέθεσε την εργασία στα Κυκλώματα 2.Η εργασία αυτη απ'ότι έχω διαπιστώσει δεν μπορεί να λυθεί χωρίς την χρήση του προγράμματος.Η απορία που έχω αυτή τη στιγμή είναι η εξής : Έχω βρει μια εξίσωση που περιέχει και πραγματικά και μιγαδικά στοιχεία.Αυτή η εξίσωση περιέχει ακόμα τα σύμβολα L,C,R  για τα οποία στο πρόβλημα δεν δίνεται κάποιο συγκεκριμένο νούμερο.Τέλος περιέχονται και τα σύμβολα u,Q,όπου το Q προκύπτει από τα L,C.Αύτό που χρειάζομαι αυτή τη στιγμή είναι ένας τρόπος ώστε να  δημιουργήσω το διάγραμα Η(u),όπου Η η συνάρτηση μου.Μπορεί να γίνει αυτό χωρίς να γνωρίζω τα L,C,R και αν ναι πώς? Δεν έχω ιδέα από Μatlab και οι προσπάθειες μέχρι τώρα δεν καρποφόρησαν....
Please όποιος μπορεί ας βοηθήσει γιατί το θέμα επείγη! Ακόμα,ένας γενικός τρόπος απεικόνισης φανταστικών συναρτήσεων θα ήταν καθ'όλα καλοδεχούμενος.

Ευχαριστώ για την βοήθεια,

Help the n00bs!    

Στην περίπτωσή σου αρμόζει η εντολή syms
Η εντολή αυτή σε επιτρέπει να δουλέψεις με σύμβολα
π.χ.
syms u;
H=u+1;

Αυτή τη στιγμή δεν θυμάμαι εάν λειτουργεί όμως με την plot.
Από ότι θυμάμαι όμως όχι.
ʼλλωστε δεν σου ζητά κάτι τέτοιο η εργασία.
Αυτό που θα κάνεις είναι να κάνεις το διάγραμμα, με παράμετρο το Q και να δεις πώς επιδρα στη συνάρτηση μεταφοράς.
Οπότε σε συνιστώ να γράψεις ένα m-file που να ξεκινά ενδεικτικά ως εξής:

Q=10;
L=Q+1
C=Q+2

Όπου βάλε αντίστοιχα το πώς ορίζονται τα C,L συνατρήσει του Q.
Δίνεις διάφορες τιμές στο Q και αμέσως σου επιστρέφει και τα αποτελε΄σματα των L,C.

Όσον αφορά τώρα τις μιγαδικές, γίνε λίγο πιο συγκεκριμένος ως προς τι ζητάς.
Το MATLAB προσφέρει πολλούς τρόπους απεικόνισης.
Εσύ τι ακριβώς θες

Χαίρεται και πάλι και ευχαριστώ για την άμεση απάντηση.

Δυστυχώς το πρόβλημα είναι πιο σύνθετο και φοβάμαι πως δεν μπορεί να λυθεί με το απλό Μ-file που περιγράφεις.... Βλέπεις,όταν η εξίσωση του Q είναι Q=1/R*(sqrtL/C), τί νόημα έχει να λύσεις ως προς L η ως προς C? Θα βγεί κατι της μορφής L=Q^2*R^2*C. Ό,τι και να κάνω με το Q, αν δεν δώσω τιμές στα  R και L και C δεν βγαίνει τίποτα...

Αλήθεια,μήπως την έχεις δεί την άσκηση? Γιατί  απ'οτι έχω καταλάβει ζητάει το Η(ju),το οποίο με κάποιο τρόπο αν μπεί στο Ματλαμπ,σαν συνάρτηση του u, εμφανίζει μια γραφική παράσταση με 2 μέγιστα τα οποία είναι και το ζητούμενο της άσκησης.Γι'αυτό προβληματίζομαι λίγο που λες να το λύσω με παράμετρο το Q.Δηλαδή να το πάω σαν Η(Q)?Το ίδιο δεν είναι?Να διευκρινίσω ότι το πρόβλημα έγκειται κυρίως στο γεγονός ότι η συνάρτηση μεταφοράς τάσης Η(ju) περιέχει εκτώς των j,u και τις μεταβλητές R,C,L οι οποίες είναι τελείως αδιευκρίνιστες και δεν  φαίνεται να υπάρχει τρόπος απαλοιφής τους από την συνάρτηση.Γι'αυτό δεν καταλαβαίνω πως μπορεί να γίνει γραφική παράσταση της συνάρτησης(που έχει και μιγαδικούς),χωρίς να δωθούν τιμές στα R,C,L,αλλά μόνη μεταβλητή να είναι το u το οποίο είναι γνωστό ότι μεταβάλλεται από 0-1.....

Όποιος μπορεί να βοηθήσει,θα 'ταν μεγάλη βοήθεια.

Βρε cyberwizard,μια και ασχοληθηκες με το θεμα απ'οτι καταλαβαινω,το Q ειναι σιγoυρο οτι ειναι 1/R*(sqrtL/C);Αν δεν κανω λαθος η ασκηση δινει τα ωn και Q οταν δεν υπαρχει μαγνητικη συζευξη.Στο κυκλωμα ομως υπαρχει.

Φίλε netfun...
Έχω ψάξει όλο το βιβλίο των κυκλωμάτων,και των ασκήσεων και προπαντώς τα κεφάλαια 7,10 για στοιχεία πάνω στο θέμα που αναφέρεις,αλλά και στην συμπεριφορά κυκλωμάτων συντονισμού σε σύζευξη.Λοιπόν,ΠΟΥΘΕΝΑ δεν υπάρχει κάποιο στοιχείο που να αναφέρει ότι το Q αλλάζει αν υπάρχει σύζευξη.Το πρόβλημα είναι ότι πουθενά δεν υπάρχει κάτι αντίστοιχο στο βιβλίο και απότι μου είπαν κάποιοι,ούτε σε κάποιο βιβλίο της βιβλιοθήκης υπάρχει κάτι που να βοηθάει...Ήθελα να 'ξερα πώς υποτίθεται ότι πρέπει να λύσουμε κάτι τέτοιο....Εγώ σου λέω,ότι όσο κι αν το χω ψάξει,αν δεν βρεθεί κάποιος να μου δώσει τα φώτα του στην Matlab(την οποία απιτούν να γνωρίζουμε,με ένα εργαστήριο 2 ωρών που κάναμε στο πρώτο εξάμηνο....),δεν πρόκειται να βγάλω κανένα αποτέλεσμα.Προς το παρόν,μέχρι να αποδειχτεί ότι το Q αλλάζει,θα το κρατήσω ως έχει.Ας έχουμε ένα δίλημμα λιγότερο γιατί αν και δώ που φτάσαμε είναι λάθος,θα τα παρατήσω.

Σορρυ που δεν απάντησα γρήγορα αλλά έχω ένα πρόβλημα με τη σύνδεσή μου.
ΟΧΙ συνάδελφε το χειρότερο που έχεις να κάνεις είναι να τα παρατήσεις.
Λοιπόν, θα κάνω μια προσπάθεια να απαντήσω ελπίζοντας να έχω καταλάβει τι ζητάς.
1. Το Η(ju) ΔΕΝ απεικονίζεται.
Ουσιαστικά απεικονίζουμε το |Η(ju)| που είναι η απόκριση συχνότητας.
Αυτό γίνεται είτε με την ενοτλή bode() είτε πιο σύνθετα με την ενολή freqresp πιο σ'υνθετα.
Έχουμε λοιπόν για μια συνάρτηση
%Η=(s^2+3s+1)/(7s^3+s+9)
%όπου s=jw,ju ή οτιδήποτε σχετικό

num=[1 3 1];
den=[7 0 1 9];
sys=tf(num,den);
bode(sys,u)
όπου u η μεταβλητή σου.
Τώρα όσον αφορά το Q,αυτό είναι παράμετρος και όχι μεταβλητή.
Έτσι λοιπόν θα κάνεις την ίδια δουλειά για διάφορες τιμές του Q.
Συνήθως παίρνουμε Qmin,Qmax,Qmin+Qmax/2
Για να έχουμε μια σφαιρική εικόνα.
ΕΠίσης για τα R,L,C αυτά καθορίζονται τόσο από το Q όσο και από διάφορες άλλες προδιαγραφές.
Να προσέξεις επίσης το εύρος τιμών.
Π.χ. στα Κυκλώματα για θεωρητικούς/διδακτικούς λόγους αναφέρεται τιμή C=1F κάτι που δεν υπάρχει στην πραγματικότητα.
Οι πυκνωτές είναι της τάξης των μF ή και μικρότερα.

Επομένως με τα δεδομένα της άσκησης δεν υπάρχουν κριτήρια για τον καθορισμό των R,L,C.
Αυτό αφήνεται στην κρίση του σχεσιαστή, στην προκειμένη περίπτωση σε σένα  ;)

Έχω συνεργαστεί με τον cyberwizard για την επίλυση της άσκησης και τελικά η τελική εξίσωση που προκύπτει δεν περιέχει R,L,C, αλλά είναι η Hv=(k*q^2)/((1+j*q*(u^2-1)/u)^2+k^2*u^2*q^2). Όπου το k είναι από 0-1 και το q ας πούμε ότι είναι 1. Το πρόβλημα τώρα είναι πως γίνονται τα διαγράμματα πλάτους και φάσης στο MATLAB; Απ'ότι ξέρω τα διαγράμματα βγαίνουν όπως στο παράδειγμα που βρίσκεται στο βιβλίο των ασκήσεων σελ 499 "Απόκριση συχνότητας", αλλά επειδή δεν είμαι και πολύ σχετικός με το MATLAB θα ήθελα τη βοήθεια σας.

Θέτοντας λοιπόν το Q=1
έχουμε Η=f(u) με παράμετρο κ.
Δίνοντας λοιπόν διάφορες τιμές στο κ
έχουμε μια μιγαδική συνάρτηση
η οποία είναι της τυπικής μορφής  H(jω).
Όπως είπα και προηγούμενα δεν απεικονίζουμε ποτέ
την συνάρτηση σε αυτή τη μορφή
αλλά απεικονίζουμε πλάτος και φάση.
Αυτό γίνεται με 2 τρόπους.
Αντικαθιστούμε το jω με s και χρησιμοποιούμε την εντολή bode()
όπως περιγράφεται παραπάνω.
Θυμηθείτε ότι ω^2= - (jω)^2 =...= - s^2

Ο δεύτερος τρόπος είναι με την εντολή freqresp()
ο οποίος όμως δεν συνίσταται τόσο σε αρχάριους
αλλά σε επιτρέπει να πάρεις τον έλεγχο στα χέρια σου.

Ψάξτε και στο help για τις 2 αυτές συναρτήσεις
και εδώ είμαστε για νεότερες απορίες.
Σε κάθε περίπτωση συνίσταται να χρησιμοποιείται το
έκδοση του MATLAB 6.5 ή νεότερη.

Μήπως θα μπορούσες να μου πεις που ειναι το λάθος πιο κάτω;

clear all
u=[0.1:0.01:10];
for i=1:1000
    s=j*u(i);
    k=0:1;
    q=1;
num(i)=(q^2)*(1-k);
den(i)=(2*s*(1-k)*q+2+2*q/s)-(s^2)*(q^2)*(1-k)^2;
H(i)=num(i)/den(i);
mag(i)=abs(H(i));
phase(i)=angle(H(i))*180/pi;
end
subplot(2,1,1)
loglog(k,mag,'r.')
grid
subplot(2,1,2)
semilogx(w,phase,'r.')
hold on
grid on
end

Ας τα πάρουμε με τη σειρά.
Καταρχάς δεν χρειάζεται το πρώτο for-loop
δες την εντολή .*
Πάντως κι έτσι σωστό είναι.

Πάμε τώρα στα λάθη.
Σωστό:  s(i)=j*u(i).
Κι αυτό διότι με τον δικό σου κώδικα το s ΔΕΝ θα γίνει πίνακας αλλά θα πάρει την τελευταία τιμή το u(i) δλδ thn u(1000).
Επίσης k=0:1 θα φτιάξει έναν πίνακα
κ= [0 1] καθότι εάν δε βάλεις βήμα το MATLAB έχει προκαθορισμένη τι΄μη την 1.
Και δεν χρειάζεται να είναι μέσα στο loop.

num(i)=(q^2)*(1-k);
Αυτό δυστυχώς είναι τελείως λάθος.
Αφού έχεις i στα αριστερά χρειάζεσαι και στα δεξιά.
Μάλλον υποθέτω k(i)?
Ομοίως και στο den(i)
Γιατί όπως είπα το k είναι ένας πίνακας 1*2
οπότε δεν μπορείς να ορίσεις ένα στοιχείο ( num(i) )  να είναι ίσο με πίνακα.

Λοιπόν ο τρόπος που διάλεξες σε γενικές γραμμές είναι σαν να πηγαίνεις Καβάλα μέσω Λονδίνου :(

Όταν έλεγα αντικατέστησε το jω με το σ εννοούσα απευθείας αντικατάσταση.
Δοκίμασε αυτό

num=[q^2)*(1-k)];
den=[A B C] όπου Α,Β,Γ οι συντελεστές του s^2, s, σταθερός όρος
                     π.χ. για το s-> B=2*(1-k)*q  εάν κατάλαβα καλά από
                     αυτά που γράφεις.
sys=tf(num,den)
bode(sys)

όπου το K έχει συγκεκριμένη τιμή.
Εάν θέλεις για διάφορες τιμές τότε περιέκλεισε το παραπάνω πρόγραμμα στον κώδικά

for i=0:0.1:1
  k=i;
  figure
................
.................
...................
end

Έτσι το κ έχει μια τιμή κάθε φορά.
Εάν για κάποιο λόγο θες να κρατήσει όλες τις τιμές τότε

k=0:0.1:1
for i=1:10
num=[???k(i)??]
den=[?k(i)???]
.............
end

Για περισσότερες πληροφορίες ανέτρεξε στο help στις συναρτήσεις tf(),bode()

to h mou prokeiptei ws ekseis:
h=Vo/Vs
h=(1/j*w)/(r+1/j*w+j*L*w+k^2*L^2*w^2/(r+j*L*w+1/(j*w*c))
kai antikathistoume me q kai u-exw lathos?????????
EPEIGON

Εάν μπορούσες να γίνεςι πιο σαφής
θα μπορούσαμε να σε βοηθήσουμε  ;)
Δεν γνωρίζω την άσκηση και οπότε
δεν μπορώ να σ'απαντήσω εάν η έκφρασή σου είναι σωστή.
Ένα λάθος όμως που βλέπω στον κώδικα είναι
το .....r+1/j*w+.......
αυτό το MATLAB το αντιλαμβάνεται σαν (1/j)*(w)
ενώ εσύ προφανώς θέλεις το 1/(j*w) το οποίο
γραφεται έτσι ή αν θες να αποφύγεις τις παρενθέσεις
1/j/w
Ο τρόπος των παρενθέσεων προτιμάται πάντα γιατί έτσι αποφεύγεται η σύγχυση

Υ.Γ: Και παράκληση μην χρησιμοποιείς greeklish.

Giannis Papadopoulos wrote:

 Απορία:
 Πώς γίνεται να χρησιμοποιήσω σε μια έκφραση μία σταθερά χωρίς να της δώσω τιμή?
 πχ. Έχω ενα μιγαδικό σε μια μορφή που το φανταστικό και πραγματικό μέρος δεν είναι χωρισμένα και θέλω μεσω του Matlab να μου δώσει το Imag και το real.Ο μιγαδικός περιέχει και δύο σταθερές οι οποίες θέλω να παραμείνουν και να μην τις αντικαταστήσω με συγκεκριμένες τιμές.
Ευχαριστώ

Εάν θέλεις να έχεις σταθερές σε μια παράσταση τότε πρέπει να χρησιμοποιήσεις την εντολή syms.
Κάποιες πληροφορίες για την εντολή θα βρεις εδώ (http://www.electronics-lab.com/forum_ee/index.php?board=27;action=display;threadid=721;start=msg2778#msg2778)
αλλά καλύτερα να ανατρέξεις και στο help.

Για να βρεις το real/imag ενός μιγαδικού αριθμού z δίνεις αντίστοιχα τις εντολές.

r=real(z)
im=imag(z)

για να βρεις πλα΄τος και φάση δίνεις

amp=abs(z)
p=angle(z) όπου η γωνία είναι σε ακτίνια

για μετατροπή χρησιμοποιείς τισ συναρτήσεις
rad2deg() & deg2rad()

εχω πελαγωσει-αν μπορεισ ριξε μια ματια στην εργασια
με ποια εντολη βρισκουμε μερικη παραγωγο στο matlab?

ode45()

Αλλά τι να την κάνεις την μερική παράγωγο;  ??? ???
Γενικά το αποτέλεσμα της άσκησης θα είναι μια συνάρτηση του u
Αντικατέστησε το ju->s και χρησιμοποίησε την bode()
Δες εδώ (http://www.electronics-lab.com/forum_ee/index.php?board=27;action=display;threadid=721;start=msg2849#msg2849) για αυτόν τον τρόπο.
Εάν συνεχίζεις να έχεις απορία του πώς γίνεται
απλά στείλε μου τη συνάρτηση H(ju,Q,k)
για να σε γράψω ακριβώς το πρόγραμμα που παράγει τα γραφήματα.

h=q/(j*u+q+j*q*u^2+k^2*q^2*j*u^4/(u+j*q*(u^2-1)))

christina δεν ξέρω αν έχεις κάνει κάπου λάθος αλλά εγώ έχω βρει ότι Hv=(k*q^2)/(1+j*q(u^2-1)/u)^2+k^2*u^2*q^2)  


Ίσως η απάντηση που έχεις βρει να είναι η ίδια με τη δική μου, όμως με έχει μπερδέψει ο τρόπος που χρησιμοποιείς τις παρενθέσεις

Απαντώ σε svistos και για σένα Χριστίνα, χρησιμοποιώντας την συνάρτηση της Χριστίνας.

Επειδή δεν μπορώ να απομονώσω το ju δεν γίνεται να χρησιμοποιήσουμε την bode. Οπότε πάμε με τον ορισμό.

q=1;
k=0.5;
u=[0.01:0.01:100]; % το εύρος συχνοτήτων
h=q./(j*u+q+j*q*u.^2+k^2*q^2*j*u.^4./(u+j*q*(u.^2-1)))
%προσέξτε τον τελεστή ".*" και "./"  είναι απαραίτητος για να κάνουμε πράξεις
%στοιχείο με στοιχείο.
amp=abs(h);
phase=angle(h);
plot(u,amp)
figure
plot(u,phase)

Βγαίνει χαμηλοπερατό φίλτρο.
Κάντε και zoom στις χαμηλές συχνότηες να φανεί καλύτερα ή
βάλτε u=[0.01:0.01:10]

Giannis Papadopoulos wrote:

Επισυνάπτω και τα αντίστοιχο m-file η εντολή φαίνεται να μη δουλεύει.Γιατί???

syms Q;
syms K;
u=[0:0.01:4]
z=K*(Q^2)*((1+j*Q*(u-(u.^(-1)))).^2-(Q*u.*K).^2).^(-1)
y=abs(z)
plot(u,y)
hold on
x=real(z)
y=imag(z)



Το MATLAB μπορεί να κάνει πράξεις με σύμβολα (syms) αλλά προφανώς δεν μπορεί να σου παραστήσει μια συνάρτηση που έχει σύμβολα (για ποιες τιμές να το κάνει??? ).
Βάλε τιμες στα K,Q και τρέξτο.
Η φιλοσοφία της άσκησης είναι να δεις την απόκριση για διάφορες τιμές των παραμέτρων. Οπότε δοκίμασε για π.χ. k=0,0.5,1 Q=1,10...

Να'σαι καλά fugitive μου!!! Αν δεν είχαμε
κι εσένα να μας λύνεις τις απορίες τι θα
κάναμε;! :P

Να κάνω κι εγώ μια ερώτηση:
επειδή η άσκηση ζητάει και μέγιστα της
συνάρτησης, μήπως να την παραγωγίζαμε
στο Matlab; Γίνεται; Η' μήπως κάνουμε
τίποτα κόλπα επίσης στο Matlab και τα
βρίσκουμε κατευθείαν;!

Έλα ρε Κίττεν χαλάρωσε...  8)
όχι και να βρεις παράγωγο για να δεις τα μεγιστα.
Έλεος!
το mATLAB είναι εδώ για σένα.
Μέγιστα βρίσκεις με την εντολή max()
ελάχιστα με το min()
Αναλυτική περιγραφή θα βρεις στο help
Επίσης για όσους έχουν.... καλά μάτια
μπορείς να το κάνεις και χειροκίνητα από το figure
κάνοντας κλικ στο μέγιστο/ελάχιστο με το ποντίκι

Κι άλλη μια ερώτηση: έκανα αυτό

q=1;
k=0.5;
u=[0.01:0.01:100];
h=(k*q^2)./((1+j*q(u.^2-1)./u.)^2+k^2*u.^2*q^2);
amp=abs(h);
phase=angle(h);
plot(u,amp)
figure
plot(u,phase)

και μου λέει ότι στο h=... έχω λάθος(;!)
"identifier" expected, ")" found.
Αφού τα έχω σωστά τα ")"! :P

Το έκανα με τη μορφή:

h=k*(q^2)*((1+j*q*(u-(u.^(-1)))).^2+(q*u.*k).^2).^(-1)

και βγήκε, αλλά η απόκριση πλάτους πρέπει να έχει
2 μέγιστα κι εμένα μου βγάζει μόνο ένα!!!!!!!! :-\

Ο μάργαρης στη σελίδα 201 του βιβλίου του εχει μια γραφική παράστση με ένα μέγιστο αλλα παρακάτω βρίσκει δύο μέγιστες τιμές εξισώνοντας το φανταστικό και το πραγματικό μέρος του μιγαδικόυ που βρήκε.Η διαφορα είναι οτι ο μιγαδικος στο βιβλίο εχει πολυ απλουστερη μορφη από αυτόν που βγαζουμε εμείς.Γι' αυτό ψάχνω τρόπο για να βρω το Im και το real του μιγαδικου συναρτησει πάντα του  Q και του k ώστε εξισώνοντας να προκύψουν του μέγιστα συναρτήσει του k όπως ακριβώς ζητά η εκφώνηση,Στη συνέχεια θα δώσουμε τιμές και στα δυο και κάνουμε τις γραφικές παραστάσεις.
Σωστά???

To MATLAB δεν μπορεί να βρει το φανταστικό και πραγματικό μέρος συμβολικών αριθμών.
Αντ'αυτού δώσε τιμές στα k,Q και μετά χρησιμοποίησε τις συναρτήσεις imag(), real()
Θα πάρεις φανταστικό και πραγματικό για την εκάστοτε τιμή.

Πρέπει να γίνει κατανοητό ότι το MATLAB ΔΕΝ φτιάχτηκε με σκοπό να δουλεύει με συμβολικές παραστάσεις. Αυτές γίνονται από τον άνθρωπο.Το MATLAB προσομοιώνει τις συναρτήσεις/εκφράσεις/παραστάσεις με αριθμητικά δεδομένα.

με ποια εντολη βρισκουμε τοπικα μεγιστα?

Η εντολή που βρίσκει τα μέγιστα
καθώς και ο κώδικας που ζητάς έχει ήδη γραφτεί
στην προηγούμενη σελίδα
"Απάντηση 20"

Loipon afou dothike h eykairia as poume merika genika pragmata gia tis ergasies gia na enhmerwneste ligo oi neoteroi. Katarxhn h proswpikh mou apopsh einai oti oi ergasies kai h diplwmatikh einai oti kalytero exei na prosferei h sxolh mias kai sou dinei th dynatothta na anaptykseis ayta pou ksereis (thewria) se praksh. Gia na kaneis mia swsth ergasia, eyparousiasth, swsth episthmonika kai tekmhriwmenh xreiazesai 3 basika pragmata. Programmatismo, organwsh kai methodikothta. Me tsapatsoulies kai arpakollatzidikes katastaseis ths theleytaias stigmhs afenws den mathaineis (basikos skopos ths ergasias) kai afeterou spane ta neyra sou kai den exeis ayto pou leme xara ths dhmiourgias.
Etsi loipon molis paroume thn ergsia sta xeria mas, kanoume programma. Efoson kseroume thn hmeromhnia paradoshs, programmatizoume etsi tis meres mas wste na brethei xronos gia ta akoloutha:
1)Anazhthsh bibliografias se bibliothikes, internet kai opou allou xwraei o nous sas.
2) Diabasma kai thewrhtikh meleth ths ergasias gia na kseroume ti kanoume, pou badizoume. Oxi ainte na piasoume to matlab kai opou mas bgalei h akrh. Elege kapoios oti "otan kserw ti kanw to kanw opws thelw". Ayto isxyei kai me tis ergasies. Otan kserw ti kanw se theorhtikh bash to ylopoiw meta opws thelw (me MATLAB, me Mathematica, me Maple, me C, me oti thelw egw). Proswpika nomizw oti ayto einai to kyriotero kommati ths ergasias. Afierwste xrono arketoutsiko edw kai tha deite pws ta pragmata ginontai eykolotera sta epomena stadia.
3) Epilegw thn glwssa (Matlab,C,Java klp) me thn opoia tha montelopoihsw (programmatisw)ayto pou exw brei theorhtika kai sxediazw to montelo. Den grafw akoma kwdika. Protimotero einai na mathw na douleyw poly kala se mia glwssa (p.x. MATLAB) kai meta na epektathw se ypoloipes. Skopos den einai na mathw oles tis glwsses tou kosmou, skopos einai na kanw th douleia mou. An ayto mporw na to petyxw mono me mia glwssa tote ayth mou arkei.
4) Syntassw kwdika kai trexw thn prosomoiwsh
5) kritikarw ta apotelesmata ths prosomoiwshs kai sygkrinw me ta theorhtika mou
6) Sth diarkeia olwn twn parapanw stadiwn exw krathsei shmeiwseis tou ti ekana, ti apotelesmata phra kai ola ta synafh. Twra sygkentrwnw olo ayto to yliko apo ta prohgoumena stadia kai syntassw mia eyparousiasth anafora(na mathete Microsoft VISIO osoi den kserete, boithaei edw).

Bebaia h seira twn bhmatwn den oritheteitai apolyta. Mporei dhladh meta thn prosomoiwsh na xreiastei na gyrisw pali gia ligo sta thewrhtika. Alla kalo einai na exw proetoimastei arketa kala (apo ta prwta bhmata) wste ayta ta piswgyrismata na elaxistopoiountai.

Skopos mou einai na teleiwsw olh ayth th diadikasia mia mera prin thn prothesmia paradoshs. Ayto giati sigoura tha exw ksexasei kati, sigoura kati den tha mou bgainei kalo, mporei o ektypwths mou na xalase kai na mhn mporw na typwsw kai xilia dyo alla tis teleytaias stigmhs. Gia ayto mia mera prin kai ola tha einai kala.

Anaferw endeiktika pws katanemw to xrono egw toulaxiston sta parapanw (p.x. gia prothesmia 15 hmerwn)
1) 1-1.5 meres
2) 3-4 meres
3) 0.5-1 meres
4) 4-5 meres
5) 0.5-1 meres
6) 1-1.5 meres
To synolo kanei sth xeiroterh 14 meres syn mia mera opws eipame gia pan endexomeno

Egw proteinw ayto to montelo. Meta apo arketh douleia panw se ergasies diapistwsa oti douleyei arista.

Dokimaste to an thelete kai peite mou entypwseis.

Flo

Na prosthesw oti se ergasies oi opoies perilambanoun peiramatikes metrhseis, aytes ginontai panta meta ti thewrhtikh meleth, wste na ginoun swsta. Prospathoume na eimaste oso akriveis ginetai. Oxi mageirema kai alatopipero sta apotelesmata, gia na bgoun opws ta theloume emeis. Kalytera einai na diabasoume apo prin kai na paroume thn wra ths metrhshs swsta apotelesmata para na psaxnomaste meta me mageiremata.

Epishs ergasies oi opoies perilambanoun ylopoihsh se hardware kanoume ola ta parapanw kai eidika thn prosomoiwsh sthn opoia dinoume megalh bash kai meta ksekiname siga siga thn kataskeyh me pollh prosoxh. Gia kataskeyastikes ergasies perissotera tips sto mellon

Flo

Καλό παίδεμα οι χρονοσειρές..  :-\
Thanks για την πληροφορια. Και στο χέρι εύκολα λύνεται, αλλά ας μάθουμε και κάτι από Matlab. Α, ο Στρίντζης είπε σήμερα ότι μας δίνει προθεσμία μέχρι την άλλη τετάρτη σε περίπτωση που δεν ήσουν στο μάθημα.  ;)

Λοιπόν αν και δεν παρέδωσα εργασία παιδεύτικα λίγο ξανά. Βάζω την συνάρτηση

h(u)=K/((1/Q^2)+((K^2)*(u^2))-((((u^2)-1)/u)^2)+j*(2/Q)*(((u^2)-1)/u)) με Q=1, Κ=0.5 και υ=[0.01:0.01:100]; (το πήρα από άπάντηση στο φόρουμ) και βγάζει το εξής μήνυμα: Error using ==> mpower
Matrix must be square. Γιατί;

Λοιπόν, αντί να σου εξηγήσω εγώ με κίνδυνο να μην καταλάβεις τίποτα, πήγαινε σελ. 483 στο βιβλίο των ασκήσεων των ηλεκτρικών κυκλωμάτων :)

Πήρες το μήνυμα από το φόρουμ αλλά το πήρες από λάθος ποστ!
Πρόσεξε λίγο καλύτερα.
Ο τελεστής που πρέπει να χρησιμοποιήσεις σε element-wise πράξεςι είναι ο "./" ή ".*" και όχι απλό "/" ή "."
Ψάξε την "Απάντηση 20"

Μεταφέρθηκαν.
Στο Κυκλώματα ΙΙ-> Εργασία

έφτιαξα την συνάρτηση:

H(u)=K./(1/Q^2)+((K^2.*(u.^2))-((((u.^2)-1)./u).^2)+j*(2/Q).*(((u.^2)-1)./u))
τη βάζω σύμφωνα με την απάντηση 20 αλλά δεν μου βγάζει τα μέγιστα και τα διαγράμματα που ζητάει!

Πρόσεξε ότι δεν βάζουμε H(u)= αλλά σκέτο Η=

Σου επισυνάπτω το γράφημα ου έβγαλα εγώ με 2 μέγιστα και τον κώδικα που χρησιμοποίησα

» plot(u,abs(H))
» K=0.5;
» Q=8;
» u=[0.1:0.01:10];
» H=K./(1/Q^2)+((K^2.*(u.^2))-((((u.^2)-1)./u).^2)+j*(2/Q).*(((u.^2)-1)./u));
» plot(u,abs(H))

Θα το δοκιμάσω αργότερα και θα σου πω. Ευχαριστώ πολύ!

Μπορείς να εξηγήσεις πωσ βρίσκουμε τα μέγιστα και ελάχιστα μιας συνάρτησης; Στη βοήθεια δεν λέει τίποτα για το ΜΑΧ (λέει μόνο για εύρεση του μεγαλύτερου από ένα πίνακα)! επείγον!

Ta panta sto MATLAB einai pinakes. Ara h synarthsh soy einai sthn ousia enas pinakas me times. An kaleseis th MAX tha sou brei to megisto tou pinaka dld to megisto ths synarthshs.
So simple...
Flo

H tanx=5*x πως λύνεται με το Matlab;
Βασικά δοκίμασα με τη solve αλλά χρειάζομαι μη μηδενική λύση και...δεν οίδα. ::)

δοκίμασε το εξής

f=tan(x)-5*x
solve(f,x)

Θα σε δώσει όλες τις λύσεις μηδενικές και μη

atan(f) δοκίμασες;

Δε βοηθάει. :'(
Δοκίμασα την solve('tanx-5*x=0') αλλά μου δίνει μόνο x=0.
Any other thoughts?

Ναι ρε συ έχει δίκιο το MATLAB.

Σου δίνει τη βασική γωνία του 1ου τεταρτημορίου ως λύση.
ΑΠό κει και περα προφανως υπάρχουν ακόμη άπειρες γωνίες που την επαληθεέυουν δεν μπορεί να σου δώσει άπειρες λύσεις.
Χρησιμοποιείς μετά την τριγωνομετρική ταυτότητα για την εφαπρομένη.
Νομίζω είναι +-2κπ

χαλαρά

Αν έχουμε μια συνάρτηση που έχει μέσα τον νεπέριο αριθμό ή το π (3,14)
Στο matlab υπάρχουν έτοιμοι αυτοί οι αριθμοί ή θα πρέπει να τους ορίσουμε;

πχ η συνάρτηση e^100πτ, μπορώ να την βάλω αυτούσια στον κώδικα με κάποιο τρόπο ή πρέπει να ορίσω το e και το π;
(αν και τώρα που το σκέφτομαι τo π το έχει ενσωματωμένο το matlab)

η παράσταση e^kati είναι η συνάρτηση exp(kati)

Μου φαινεται το π ειναι η μεταβλητη pi. Αλλα δεν ειμαι σιγουρος. Το e ειναι το exp(1).

Έτσι είναι.

Ναι το pi είναι, αλλά το matlab σ'αυτή τη μεταβλητή κρατάει μόνο 4 δεκαδικά. Αν θέλεις να έχεις περισσότερα βάλτο σε δικιά σου μεταβλητή... :)

Ή όπως λέει η τεκμηρίωση του matlab
pi=4*atan(1)
=imag(log(-1))
 :P

Αυτό είναι λάθος!Το pi στο μάτλαμπ έχει την μέγιστη δυνατή ακρίβεια (double) όπως μπορείς να διαπιστώσεις αν γράψεις:

format long
disp(pi)

Απλά εμφανίζει τα 4 πρώτα δεκαδικά για ευκολία(μπορείς να το αλλάξεις αυτό).Θα ήταν παράλογο εξάλλου να μην είχε ακριβή τιμή το π σε ένα επαγγελματικό πακέτο μαθηματικών.

Μάλιστααα! Δεν το ήξερα αυτό... Και μου φαινόταν παράλογο όντως! :)

ok παιδιά thanks!

  Παιδιά, γράφω εδώ για να μην ανοίγω καινούργιο νήμα ;). Ξέρετε τι γίνεται με το matlab και τα Linux; Έβαλα πριν λίγο καιρό ubuntu δοκιμαστικά μπας και μείνω.
  Ας, κάνω λίγο πιο συγκεκριμένη την ερωτησή μου. Εννοώ, θυμάμαι οτι συζητιόταν πριν κανα χρόνο πως γίνεται η εγκατάσταση σε λινουκς, το έχει κάνει κανείς;

Η u(t) δεν υπάρχει έτσι; απλα ορίζουμε τον χρόνο απο 0...όσο θέλουμε. σωστά;




Υπάρχει ειδική έκδοση για το λινουξ (δεν ξέρω περρισότερα) αυτή έχεις;

Εχω το Ματλαμπ για Linux, ευγενικη χορηγια του chggr. Εχει επακριβεις οδηγιες για το πως εγκαθισταται. Για περισσοτερες πληροφοριες πμ.

Η u(t) ειναι απλα μοναδα για t>0.

Αρα βαζεις

fτιαχνεις ενα m file που να εχει

function ret = u

if t >= 0
  ret = 1 ;
else
  ret =0 ;
return u;

Μπορει να υπαρχει και ετοιμη. Γραψε step function

To matlab για λίνουξ υπάρχει και λειτουργεί άψογα.Εγώ χρησιμοποιώ την τελευταία έκδοση για 64 bit(R2007b).

Τα ίδια ισχύουν και για το Mathematica(32 και 64 bits).Μάλιστα από την έκδοση 6 και έπειτα,το γραφικό περιβάλλον του τελευταίου(frontend) έχει μεταγραφεί από Μotif σε QT, οπότε είναι πλέον το ίδιο όμορφο και λειτουργικό με το αντίστοιχο των windows!

Εγω ειχα το Matlab 7, αλλα δεν δουλευε αψογα. Οταν πηγαινες να κανεις ltiview κολλουσε το σμπαν και δεν εκανε τιποτα. Επανειλημμενα.

Θα συμφωνήσω με Ζαρατούστρα. Το MATLAB στο Linux τρέχει πολύ πιο γρήγορα και ελαφριά. Σε μένα τουλάχιστον δεν έχει βγάλει κανένα πρόβλημα.

Για χρησιμοποιησε το ltiview και θα σου πω. Γενικα, καθε φορα που το ετρεχα ελεγε οτι ειχε προβλημα στην ltiview. Μπορει να ειναι θεμα του αρχειου εγκαταστασης.

Δε το έχω χρησιμοποιήσει ποτέ και δεν ξέρω. Δώσε ένα δείγμα κώδικα να σου πω αν κολλάει και σε μένα ή όχι...

Βρε γραψε απλα αυτο
(αν θυμαμαι καλα)

fun = tf([1 1], [1 1 1]) ;

ltiview(fun) ;

και δες αν δουλευει.

Αν θυμαμαι καλα, η tf με αυτον τον τροπο φτιαχνει μια συναρτηση μεταφορας και η ltiview σου βγαζει τα παντα, φασμα κτλ...

Όχι ότι θέλω να σου μπω στο μάτι αλλά δεν κολλάει σε μένα :P :P

Το ltiview δεν παρουσιάζει κανένα πρόβλημα εδώ ( 7.5.0.338(R2007b) )

Επίσης,από διάφορα benchmarks που έχω δει προκύπτει ότι το matlab στο linux έχει καλύτερη απόδοση από ότι στα windows,πράγμα διόλου περίεργο αν λάβει κανείς υπόψη πως η φιλοσοφία του πακέτου είναι καθαρά unixοειδής.Είμαι σίγουρος πως οι πρώτες εκδόσεις έτρεχαν σε κάποιο unix.

Ναι δεν λεω πως σε ολα εχει προβλημα. Λεω οτι αυτο που ειχα εγω ειχε προβλημα. Αρα δεν σημαινει οτι υπαρχει σε ολους.

Δεν είναι ακριβώς έτσι, γιατί η συνάρτηση u(t) είναι ασυνεχής και ιδανικά
πετάγεται ακαριαία από το μηδέν στο 1 μόλις γίνει θετικός ο χρόνος,

Για να το προσομοιώσεις αυτό στο MATLAB θα πρέπει να πάρεις πολύ
μικρά χρονικά διαστήματα, το πόσο μικρά θα είναι αυτά εξαρτάται από την εφαρμογή σου.

Τέλος υπάρχει συνάρτηση step() αλλά από όσο θυμάμαι,
χρησιμεύει για να ορίσεις τον τ΄ροπο διέγερσης σε σύστημα αυτομάτου ελέγχου,
δεν μπορεί δλδ να χρησιμοποηθεί ως χρονική συνάρτηση όπως π.χ. η cos(t)

Με το που βαζω το cd για να κανω εγκατασταση, μου βγαζει αυτο

"There is an invalid character (!) in Your temp directory. Please remove the invalid character from your TEMP (or TMP) environment variable and run the installer again"

Πως το διορθωνω???Εψαξα λιγο στο google αλλα δεν καταλαβα :(

παραθύρια ή πιγκουίνο;

παραθυρια

θα μιλήσω για ΧΡ γιατί τα σΒΗΣΤΑ δεν τα έχω δουλέψει και δεν ξέρω.
Καταρχάς άνοιξε το run και γράψε %Temp%
αυτό θα σε πάει κατευθείαν στο τεμπ folder.
Σβήσε ότι υπάρχει εκεί μέσα (κάποια αρχεία δεν θα μπορείς γιατί είναι σε χρήση).

Μετά πήγαινε στο Control Panel ->ADvanced -> Environment variables
και δες τη μεταβλητή για το Temp πως είναι το path
Θα πρέπει να είναι κάτι τέτοιο:
%USERPROFILE%\Local Settings\Temp

Εάν όλα αυτά δείχνουν ΟΚ δεν μπορώ να πώ κάτι άλλο,
μήπως φταίει και το αρχείο εγκατάστασης

Eχει παρα πολλα αρχεια, ολα ολα????? Το αρχειο εγκαταστασης δεν πρεπει να εχει προβλημα σε αλλους το ιδιο cd ειναι οκ :(

# Κάντε κλικ με το δεξιό κουμπί του ποντικού στο εικονίδιο Ο Υπολογιστής μου (My Computer) και, στη συνέχεια, κάντε κλικ στην εντολή Ιδιότητες (Properties).
# Κάντε κλικ στην καρτέλα Για προχωρημένους (Advanced).
# Κάντε κλικ στην επιλογή Μεταβλητές περιβάλλοντος (Environment Variables).



γράψε τι λεει στα ΤΕMP TMP,  στις μεταβλητές χρήστη και μεταβλητές συστάματος

Αυτα τα εψαξα και χθες γραφει οτι λεει ο fugifox..

Σβηνω δηλαδη οτι εχει στον Temp και ξαναπροσπαθω για εγκατασταση??

ε γράψε τι λέει και στο userprofile

Μπορεί να μου πει κανείς με ποιά λογική ο τελεστής ' (δηλαδή transpose) δεν είναι μόνο transpose αλλά είναι transpose conjugate;;;
Μου πήρε καμιά ώρα για να βρω ότι μόνο αυτό έφταιγε που δε δούλευε σωστά όλος ο κώδικάς μου!!!
By the way το απλό transpose για πίνακες μιγαδικών είναι: .'

Γκρίνιααααα!!!