[anagnorisi] Πρώτες απορίες για τη Β' Εργασία...

  • From: "Konstantinos Iliopoulos" <digitaldiaphony@xxxxxxxxx>
  • To: anagnorisi@xxxxxxxxxxxxx
  • Date: Sat, 1 Apr 2006 20:27:54 +0300

Χαιρετώ άπαντες :-). Καλό σας μήνα!

Λοιπόν, νομίζω ότι ήρθε η ώρα να τεθεί σε εφαρμογή η ο σκοπός της λίστας...
Κατ'αρχάς παραθέτω τα αποτελέσματα που έχω από τα ερωτήματα της Α' άσκησης
(πρσχ: δεν κανονικοποίησα τα δείγματα....)

Α' Άσκηση
---------------

1)
- λ1 = 8.6778e+005
- u1 = 0.9955    0.0121    0.0223    0.0162    0.0895
- Θεωρητικό σφάλμα = 2.6317e-004
- Μέσο σφάλμα =  0.0149

- λ2 = 6.6192e+005
- u2 = 0.9822    0.0099   -0.0000    0.0423    0.1828
- Θεωρητικό σφάλμα =  0.0026
- Μέσο σφάλμα = 0.0446

- λ3 = 1.0721e+006
- u3 = 0.9828    0.0953    0.1550    0.0077    0.0315
- Θεωρητικό σφάλμα = 0.0206
- Μέσο σφάλμα = 0.1298

2)
- Σφάλμα (%) για τα στοιχεία του Τ1 : 0
- Σφάλμα (%) για τα υπόλοιπα στοιχεία της ομάδας 1 : 0

- Σφάλμα (%) για τα στοιχεία του Τ2 : 10
- Σφάλμα (%) για τα υπόλοιπα στοιχεία της ομάδας 2 : 1.5625

- Σφάλμα (%) για τα στοιχεία του Τ3 : 22.8571
- Σφάλμα (%) για τα υπόλοιπα στοιχεία της ομάδας 3 : 31.6667

3,4) (Αν θέλει κανείς να δεί σχήματα ας μου πει να τα επισυνάψω στο επόμενο
e-mail)

5) (Εδώ με παραξενεύουν τα αποτελέσματα)
- γωνία μεταξύ u1, u2: 5.7517 (σε μοίρες)
- γωνία μεταξύ u1, u3: 9.6185
- γωνία μεταξύ u2, u3: 13.5177

Επίσης μια απορία : Με το [V,D] = eig(A, 'nobalance') δεν έπρεπε να μου
κανονικοποιεί τα ιδιοδιανύσματα, όμως το Matlab παρουσιάζεται λιγάκι
αυταρχικό σε αυτό τον τομέα :-). Έχει κανείς κάποια άλλη λύση?



Όσον αφορά τώρα τη δεύτερη άσκηση, συγκεκριμένα στο Β1, συμβαίνει το εξής:
έχω τον κώδικα και τρέχει εδώ και κανένα 20λεπτο και αποτελέσματα δεν βλέπω.
Αυτό σημαίνει ότι ο αλγόριθμος δεν καταλήγει σε αποτέλεσμα με τίποτε!
(δηλαδή σε καμία φάση δεν είναι όλα τα στοιχεία του e(n) όλα θετικά ή όλα
αρνητικά. Αν ο δικός σας συγκλίνει, τότε για check-άρετε λιγάκι το σημείο
του κώδικα που αφορά τη σύγκλιση...

Β' Άσκηση
--------------

% mexri edw exw ypologisei ton pinaka Y, toy opoioy oi grammes
% 1-40 periexoun ta 40 deigmata tou T1, oi grammes 41-80 ta antitheta twn
% deigmatwn tou T2. Genikws ta x(i) pou periexei o pinakas einai
% epay3hmena me th monada (dhladh h sthlh 14 gia tis prwtes 40 grammes
% periexei "1" kai % h idia sthlh gia tis 40 teleutaies grammes periexei
"-1".
% Genikws oson afora ta bhmata pou akolouthw, auta einai parmena apo
% askhsh pou exei kanei o Iordanhs sthn ta3h stis 08.03.2006. Loipon...

YPseydoInverse = inv((Y' * Y)) * Y';
B = ones(80,1);
E = zeros(80,1);
doMore = true;

while doMore == true
    B = B + 1*(E + abs(E));            % c = 1
    W12 = YPseydoInverse * B;
    E = Y*W12 - B;
    if abs(E) == E
        doMore = false;
        result = 1;                         %  sygklish
    elseif abs(E) == -E
        doMore = false
        result = 0;                         % mh syglkish
    end

SoS! Δεν τελειώνει ποτέ αυτή η loop!

Αυτά για την ώρα... Θα τα πούμε σύντομα. Ελπίζω να βοήθησα περισσότερο
απ'ότι κούρασα!


P.S: Ρε Στέφανε, γιατί κάθε λέξη σου αρχίζει με κεφαλαίο γράμμα? :-)
P.S2: Για να απαντήσω σε μια προηγούμενη πρόταση, αν πιστεύετε ότι θα
βοηθήσει μια λίστα και στα ψηφιακά, εγώ τουλάχιστον είμαι μέσα. Υποθέτω
μάλιστα ότι εκεί θα υπάρξει αρκετά μεγαλύτερη συμμετοχή.
P.S3: Καλή συνέχεια σε όλους.

Other related posts: