[anagnorisi] Re: Ðñþôåò áðïñßåò ãéáôç Â' Åñãáóßá...
- From: "Anastasios I Zafiropoulos" <zafeiropoulos@xxxxxxxxxxx>
- To: <anagnorisi@xxxxxxxxxxxxx>
- Date: Mon, 3 Apr 2006 23:07:55 +0300
Κωνσταντίνε συμφωνούμε στα αποτελέσματα!
Περιμένουμε νέα σου για το 2ο κομμάτι της άσκησης εκτός τσε αν αποφάσισες να
το αφήσεις σε αέναο βρόχο.
10^10 loops means nothing to me.
I run an endless loop in a couple of hours DOT
Mpimperis ©
_____
From: anagnorisi-bounce@xxxxxxxxxxxxx
[mailto:anagnorisi-bounce@xxxxxxxxxxxxx] On Behalf Of Konstantinos
Iliopoulos
Sent: Saturday, April 01, 2006 7:28 PM
To: anagnorisi@xxxxxxxxxxxxx
Subject: [anagnorisi] Πρώτες απορίες για τη Β' Εργασία...
Χαιρετώ άπαντες :-). Καλό σας μήνα!
Λοιπόν, νομίζω ότι ήρθε η ώρα να τεθεί σε εφαρμογή η ο σκοπός της λίστας...
Κατ'αρχάς παραθέτω τα αποτελέσματα που έχω από τα ερωτήματα της Α' άσκησης
(πρσχ: δεν κανονικοποίησα τα δείγματα....)
Α' Άσκηση
---------------
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:
