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.

[Karaμazoβ]

Βασικά αυτό με το 0 δεν το ήξερα απλά στην εκφώνηση λέει απο 0 ως μια τιμή..πως θα το κάνω αυτό?

τα ερρορ είναι 2 :

Error in ==> Askisi1_1 at 18
    G(f)=(Yo-yin)/(Yo+yin);

Error in ==> run at 57
          evalin('caller', [s ';']);

αυτό το f τι είναι??

[Lampros]

Μη κoλλάς σε λεπτομέρειες ρε συ...αντί για i έχω f...

[Lampros]

Λοιπόν για να μην μπερδευόμαστε ο κώδικας είναι αυτός

Code:

Zo=50;
f0=1;
Yo=1/Zo;
R=100;
C=2*10^(-3);
N=201;
GL=1/R;
for f=0:(4*f0/N):(4*f0)
    BL=C*2*pi*f;
    yL=GL/Yo+1i*BL/Yo;
    y1in=Yo*((yL+1i*Yo*tan(2*pi*0.2))/(Yo+1i*yL*tan(2*pi*0.2)));
    y2=Yo*1i*tan(2*pi*0.13);
    y2in=y1in+y2;
    y3in=Yo*((y2in+1i*Yo*tan(2*pi*0.1))/(Yo+1i*y2in*tan(2*pi*0.1)));
    B2=2.7*10^(-3)*2*pi*f;
    b2=B2/Yo;
    yin=y3in+1i*b2;
    G(f)=(Yo-yin)/(Yo+yin);
end

[Grecs]

Aφου εχεις f απο 0 μας δουλευεις?

βαλε G(f+1)

[Karaμazoβ]

Μη κoλλάς σε λεπτομέρειες ρε συ...αντί για i έχω f...

οκ ρε, γιατι θύμησε την συχνοτητα f, και είπα μηπως.

Λοιπον, μπορει να ειναι αυτονόητο και να το χεις κανει, αλλα το λεω έτσι κι αλλιως.
Εστω πως σου λεει απο 0 ως 10 GHz, και για καλή (?) ακρίβεια παίρνεις 100 δείγματα.

αρα Δf = 10/100 = 0.1 GHz το βημα της δειγματοληψίας. Αρα ο πινακας θα εχει n=100 στοιχεία με f(1)=0;, f(2)=0.1; f(3)=0.3; ... f(100)=10;


επομενως ετσι βγαζεις τους συντελεστες ανακλασης απο 0 ως 10...αρκει φυσικα να υπολογίζεις την σωστή Yin για καθε δείγμα βεβαια!

[bakeneko]

 f=0:(4*f0/N):(4*f0)

Εκτός από το 0 που έχεις, το πρόβλημα για το μάτλαμπ είναι αυτά τα άλλα f δε βγαίνουν ακέραιοι, δηλαδή δεν μπορείς να βάλεις G(1,5). Θα πρέπει να βρεις έναν άλλο τρόπο να το κάνεις, π.χ. να βάλεις οι κάθε 10 στήλες του G να αντιπροσωπεύουν 1 μονάδα του οριζόντιου άξονά σου! (δεν ξέρω ακριβώς τι κάνεις οπότε δεν μπορώ να σου πω και ακριβώς, απλά πετάω μια ιδέα )

[Lampros]

f=0:(4*f0/N):(4*f0)

Εκτός από το 0 που έχεις, το πρόβλημα για το μάτλαμπ είναι αυτά τα άλλα f δε βγαίνουν ακέραιοι, δηλαδή δεν μπορείς να βάλεις G(1,5). Θα πρέπει να βρεις έναν άλλο τρόπο να το κάνεις, π.χ. να βάλεις οι κάθε 10 στήλες του G να αντιπροσωπεύουν 1 μονάδα του οριζόντιου άξονά σου! (δεν ξέρω ακριβώς τι κάνεις οπότε δεν μπορώ να σου πω και ακριβώς, απλά πετάω μια ιδέα )

Σωστός ρε..λολ

Απλά πως θα το κάνω αυτό.Να βάλω μια δεύτερη λούπα?

[Karaμazoβ]

f=0:(4*f0/N):(4*f0)

Εκτός από το 0 που έχεις, το πρόβλημα για το μάτλαμπ είναι αυτά τα άλλα f δε βγαίνουν ακέραιοι, δηλαδή δεν μπορείς να βάλεις G(1,5). Θα πρέπει να βρεις έναν άλλο τρόπο να το κάνεις, π.χ. να βάλεις οι κάθε 10 στήλες του G να αντιπροσωπεύουν 1 μονάδα του οριζόντιου άξονά σου! (δεν ξέρω ακριβώς τι κάνεις οπότε δεν μπορώ να σου πω και ακριβώς, απλά πετάω μια ιδέα )

Σωστός ρε..λολ

Απλά πως θα το κάνω αυτό.Να βάλω μια δεύτερη λούπα?

αποσυμπλεξε τον υπολογισμο του f απο το loop.
Βαλε i=1:n και μεσα στο loop βαλε f(i) = i*df; , οπου df το βημα ((fmax-fmin)/n Ghz). Μετα οπου f, αντικατεστησε με f(i).


στο τελος plot(f,G)

[Grecs]

f=0:(4*f0/N):(4*f0)

Εκτός από το 0 που έχεις, το πρόβλημα για το μάτλαμπ είναι αυτά τα άλλα f δε βγαίνουν ακέραιοι, δηλαδή δεν μπορείς να βάλεις G(1,5). Θα πρέπει να βρεις έναν άλλο τρόπο να το κάνεις, π.χ. να βάλεις οι κάθε 10 στήλες του G να αντιπροσωπεύουν 1 μονάδα του οριζόντιου άξονά σου! (δεν ξέρω ακριβώς τι κάνεις οπότε δεν μπορώ να σου πω και ακριβώς, απλά πετάω μια ιδέα )

Σωστός ρε..λολ

Απλά πως θα το κάνω αυτό.Να βάλω μια δεύτερη λούπα?

Βαλε μια αλλη μεταβλητη στο G(a) px πριν την λουπα a=1 και μεσα στην λουπα a++

Αμα θες και τα f βαλτα σε μια γραμμη του G απο διπλα.

[Lampros]

f=0:(4*f0/N):(4*f0)

Εκτός από το 0 που έχεις, το πρόβλημα για το μάτλαμπ είναι αυτά τα άλλα f δε βγαίνουν ακέραιοι, δηλαδή δεν μπορείς να βάλεις G(1,5). Θα πρέπει να βρεις έναν άλλο τρόπο να το κάνεις, π.χ. να βάλεις οι κάθε 10 στήλες του G να αντιπροσωπεύουν 1 μονάδα του οριζόντιου άξονά σου! (δεν ξέρω ακριβώς τι κάνεις οπότε δεν μπορώ να σου πω και ακριβώς, απλά πετάω μια ιδέα )

Σωστός ρε..λολ

Απλά πως θα το κάνω αυτό.Να βάλω μια δεύτερη λούπα?

αποσυμπλεξε τον υπολογισμο του f απο το loop.
Βαλε i=1:n και μεσα στο loop βαλε f(i) = i*df; , οπου df το βημα ((fmax-fmin)/n Ghz). Μετα οπου f, αντικατεστησε με f(i).


στο τελος plot(f,G)

έτσι όμως η λούπα θα κάνει μόνο 4 επαναλήψεις(στη δική μου περίπτωση)...

[bakeneko]

Αν θέλεις να γίνεις λίγο μπακάλης πάντως μπορείς απλά αντί για "f" (στο G(f)) να βάλεις round(100*f), όπου round κάνει στρογγυλοποίηση στον κοντινότερο ακέραιο, και αρκεί βέβαια το 100 να σου δίνει την επιθυμιτή ανάλυση στον άξονα (π.χ. μη σου βγει 2 φορές το ίδιο νούμερο με τη στρογγυλοποίηση), αλλιώς μπορείς να βάλεις και 1000 κτλ, και απλά μετά να το έχεις υπόψη σου αυτό στους υπόλοιπους υπολογισμούς, και αν εμφανίζονται νούμερα στο plot να αναφέρεις την κλίμακα!

[Karaμazoβ]

f=0:(4*f0/N):(4*f0)

Εκτός από το 0 που έχεις, το πρόβλημα για το μάτλαμπ είναι αυτά τα άλλα f δε βγαίνουν ακέραιοι, δηλαδή δεν μπορείς να βάλεις G(1,5). Θα πρέπει να βρεις έναν άλλο τρόπο να το κάνεις, π.χ. να βάλεις οι κάθε 10 στήλες του G να αντιπροσωπεύουν 1 μονάδα του οριζόντιου άξονά σου! (δεν ξέρω ακριβώς τι κάνεις οπότε δεν μπορώ να σου πω και ακριβώς, απλά πετάω μια ιδέα )

Σωστός ρε..λολ

Απλά πως θα το κάνω αυτό.Να βάλω μια δεύτερη λούπα?

αποσυμπλεξε τον υπολογισμο του f απο το loop.
Βαλε i=1:n και μεσα στο loop βαλε f(i) = i*df; , οπου df το βημα ((fmax-fmin)/n Ghz). Μετα οπου f, αντικατεστησε με f(i).


στο τελος plot(f,G)

έτσι όμως η λούπα θα κάνει μόνο 4 επαναλήψεις(στη δική μου περίπτωση)...

γιατι 4 ρε? βαλε εκει 100-500 να δεις την υγεια σου

[Lampros]

λοιπόν έκανα τις εξής αλλαγές και μοιάζει να δουλεύει..

Code:

Zo=50;
f0=1;
Yo=1/Zo;
R=100;
C=2*10^(-3);
N=201;
GL=1/R;
j=1;
for f=0:(4*f0/N):(4*f0)
    F(j)=f;
    BL=C*2*pi*f;
    yL=GL/Yo+1i*BL/Yo;
    y1in=Yo*((yL+1i*Yo*tan(2*pi*0.2))/(Yo+1i*yL*tan(2*pi*0.2)));
    y2=Yo*1i*tan(2*pi*0.13);
    y2in=y1in+y2;
    y3in=Yo*((y2in+1i*Yo*tan(2*pi*0.1))/(Yo+1i*y2in*tan(2*pi*0.1)));
    B2=2.7*10^(-3)*2*pi*f;
    b2=B2/Yo;
    yin=y3in+1i*b2;
    G(j)=10*log((Yo-yin)/(Yo+yin));
    if (G(j)<-40)
        G(j)=-40;
    end
    j=j+1;
end

το τελευταίο κομμάτι είναι γιατί θέλω τις τιμές σε dB..μου φαίνεται όμως πως έχω κάποιο λογικό λάθος.
Τέλος πάντων ευχαριστώ για τις απαντήσεις με βοηθήσατε!

[Lampros]

Για να πάρω το μέτρο ενός μιγαδικού αριθμού, θα χρησιμοποιήσω την συνάρτηση abs(value) ?

(Το ρωτάω γιατί μου βγάζει περίεργα αποτελέσματα και ενώ έχω ψάξει το πρόγραμμά μου πολύ δεν βρίσκω κάποιο άλλο λογικό λάθος)

[il capitano]

Θεωρητικά ναι, αλλά νομίζω και μένα μου παρουσίαζε προβλήματα κατά καιρόυς, και τότε έκανα το παλιό καλό

Code:

sqrt((real(x)).^2+(imag(x)).^2)