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

[Πληροφορίες Καθηγητών]

Πληροφορίες Καθηγητών

[fugiFOX]

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

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

[Gostas]

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

[fugiFOX]

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

[cyberwizard]

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

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

Help the n00bs!    

[fugiFOX]

Στην περίπτωσή σου αρμόζει η εντολή 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 προσφέρει πολλούς τρόπους απεικόνισης.
Εσύ τι ακριβώς θες

[cyberwizard]

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

Δυστυχώς το πρόβλημα είναι πιο σύνθετο και φοβάμαι πως δεν μπορεί να λυθεί με το απλό Μ-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.....

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

[Netfun]

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

[cyberwizard]

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

[fugiFOX]

Σορρυ που δεν απάντησα γρήγορα αλλά έχω ένα πρόβλημα με τη σύνδεσή μου.
ΟΧΙ συνάδελφε το χειρότερο που έχεις να κάνεις είναι να τα παρατήσεις.
Λοιπόν, θα κάνω μια προσπάθεια να απαντήσω ελπίζοντας να έχω καταλάβει τι ζητάς.
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.
Αυτό αφήνεται στην κρίση του σχεσιαστή, στην προκειμένη περίπτωση σε σένα  

[svistos]

Έχω συνεργαστεί με τον 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 θα ήθελα τη βοήθεια σας.

[fugiFOX]

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

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

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

[svistos]

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

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

[fugiFOX]

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

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()

[christina]

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

[fugiFOX]

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.