Transcript ppt/651kB

Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Τυχαία συνάρτηση

Μία τυχαία συνάρτηση (ΤΣ)
είναι ένας κανόνας με τον
οποίο σε κάθε αποτέλεσμα ζ
ενός
πειράματος
τύχης
αντιστοιχείται μία συνάρτηση
x(t,ζ) όπου t=(x, y, … , t) ένα
διάνυσμα συντεταγμένων στο
χώρο των n διαστάσεων, η
οποία ονομάζεται υλοποίησή
της. Δηλαδή μία ΤΣ είναι μία
οικογένεια συναρτήσεων με
παράμετρο το ζ.
f(t)
t
t1
tk
tk+1
1
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Υλοποιήσεις


Στο πρώτο σχήμα φαίνονται
πολλές υλοποιήσεις της ΤΣ
Χ(t) οι οποίες προκύπτουν
από ισάριθμες επαναλήψεις
του πειράματος τύχης.
Στο
δεύτερο
σχήμα
φαίνονται δύο διαφορετικές
υλοποιήσεις μίας ΤΣ σε δύο
διαστάσεις. Είναι εμφανής η
ομοιότητα μεταξύ των δύο
υλοποιήσεων,
παρά
το
γεγονός
ότι
υπάρχουν
τυχαίες διαφοροποιήσεις.
2
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Περιγραφή μίας ΤΣ

Για την πλήρη περιγραφή των στατιστικών
ιδιοτήτων μίας ΤΣ απαιτείται η γνώση της
συνάρτησης F(x1,… ,xn; t1,…, tn) για κάθε xi, ti
και n.
3
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Η πράξη


Στην
πράξη
είναι
αρκετή
η
χρησιμοποίηση των ροπών πρώτης και
δεύτερης τάξης της Χ(t) :
Συνάρτηση μέσης τιμής m(t) της Χ(t)
είναι η μέση τιμή της ΤΜ Χ(ti) για κάθε
t i:

m( t )  E{X( t )} 
m(t)=c
t
y
m(t)=at2
 xf (x, t )dx


x
Συνάρτηση συνδιασποράς C(t1,t2) της
ΤΣ X(t) είναι η συνδιασπορά των ΤΜ X1
και Χ2 για κάθε συνδυασμό t1 και t2:
C(t1,t2)=Cov(X(t1),X(t2))
t
4
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Ένα παράδειγμα

Χρησιμοποιώντας την συνάρτηση – γεννήτρια
τυχαίων αριθμών του Excel δημιουργήστε μία
στήλη με 20 κελιά ως μία υλοποίηση της
τυχαίας συνάρτησης Χ(t). Ποιά είναι η μέση
τιμή της Χ(t); Κάντε το ίδιο και στην διπλανή
στήλη. Τί παρατηρείτε στην πρώτη στήλη μετά
την δημιουργία της δεύτερης; Τι αποτελεί η
δεύτερη στήλη σε σχέση με την Χ(t);
5
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Δομή και συνδιασπορά


Σύμφωνα με τα ανωτέρω, μία ΤΣ,
η οποία δεν είναι παρά μία
ακολουθία ΤΜ με παράμετρο το
t, θα έχει στον χώρο Hilbert την
μορφή που δίνεται στο σχήμα,
όπου, βέβαια, ο περιορισμός σε
τρείς διαστάσεις είναι τεχνικός.
Τα μήκη και οι σχετικές θέσεις
των n διανυσμάτων Χ(ti), ι=1, …,
n καθορίζονται από τον εκάστοτε
μη αρνητικά ορισμένο πίνακα
συνδιασπορών Cn .
cos  
v
X(t1)
X(t2)
CovX ( t1 ), X ( t 2 ) 
Var X ( t1 ) Var X ( t 2 ) 
w
φ
X(tn)
u
6
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Απλουστευμένη εικόνα

Στο σχήμα γίνεται μία
προσπάθεια
περαιτέρω
απλούστευσης
της
παράστασης μίας ΤΣ μέσω
ισομορφισμού του χώρου
Hilbert προς τον χώρο του
κοιτάσματος. Αυτό βέβαια
είναι εφικτό μόνον στην
περίπτωση που ο χώρος
Hilbert έχει μόνον 3
διαστάσεις.
X(tn)
z
X(t2)
y
t2(x2,y2,z2)
tn(xn,yn,zn)
X(t1)
t1(x1,y1,z1)
x
7
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Μία σημαντική ιδιότητα


Επειδή η συνάρτηση συνδιασποράς C(ti,tj) προσδιορίζει
τις τιμές των στοιχείων του Cn, έπεται ότι είναι επίσης μία
μη αρνητικά ορισμένη συνάρτηση.
Η ιδιότητα αυτή έχει πολύ μεγάλη πρακτική σημασία
εφόσον περιορίζει την κλάση των συναρτήσεων που
μπορούν να θεωρηθούν ως μοντέλα συναρτήσεων
συνδιασποράς. Η σημασία της θα φανεί στο κεφάλαιο 7. Ο
έλεγχος του κατά πόσο μία συνάρτηση είναι μη αρνητικά
ορισμένη γίνεται μέσω του μετασχηματισμού Fourier της,
δηλαδή του «φάσματός» της, το οποίο πρέπει να είναι μη
αρνητικό.
8
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Μοντέλα ΤΣ: ΟΣΕ



Μία ΤΣ Χ(s) ονομάζεται ομοιογενής υπό την
στενή έννοια (ΟΣΕ), εάν οι στατιστικές της
ιδιότητες
παραμένουν
αμετάβλητες
σε
οποιαδήποτε μετατόπιση της αρχής των αξόνων.
Αυτό σημαίνει ότι οι ΤΣ Χ(s) και Χ(s+c) έχουν τις
ίδιες στατιστικές ιδιότητες για οποιοδήποτε
διάνυσμα c του ευκλείδειου χώρου Rd.
f(x;s) = f(x)
f(x1,x2 ; s1,s2) = f (x1,x2 ; h) όπου h = s1-s2
9
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Μοντέλα ΤΣ: ΟΕΕ


Μία ΤΣ Χ(s) ονομάζεται ομοιογενής υπό την
ευρεία έννοια (ΟΕΕ) εάν η μέση τιμή της
παραμένει σταθερή
Ε{Χ(s)}=m
και η συνάρτηση συνδιασποράς της εξαρτάται
μόνον από το διάνυσμα της απόστασης h=s1-s2
C(s1,s2)=Cov[X(s+h),X(s)]=C(h)
10
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Μοντέλα ΤΣ: ΙΕΕ

Μία ΤΣ Χ(s) ονομάζεται ισοτροπική με την
ευρεία έννοια (ΙΕΕ), εάν εκτός των συνθηκών
της ΟΕΕ ισχύει:
C(si,sj) = C(|si-sj| = r) = C(|h|)
11
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Γιατί ΟΕΕ;



Η υπόθεση της ΟΕΕ είναι βασική για την εφαρμογή των
συμπερασμάτων της θεωρίας των ΤΣ, στις περιπτώσεις
που διατίθεται μόνον μία υλοποίηση της ΤΣ, περίπτωση
στην οποία εμπίπτουν και σχεδόν όλες οι εφαρμογές της
γεωπληροφορικής.
Είναι απαραίτητος ο υπολογισμός μέσων τιμών του τύπου
Ε{Χ(s)} για συγκεκριμένο σημείο s, πράγμα αδύνατο
εφόσον διατίθεται μόνον μία υλοποίηση της ΤΣ,
Γίνεται η παραδοχή ότι δύο δείγματα στις θέσεις s1 και s2
προέρχονται από την ίδια θέση s αλλά από διαφορετικές
υλοποιήσεις.
12
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Μοντέλα ΤΣ: ΟΤΣ
x
OEE
y
m(t)=c
t


OTΣ
m(t)=c
t
Εάν δοθεί μία ΤΣ Χ(s), ορίζεται ως ΤΣ των διαφορών
της η
ΧΔ(s;h)≡ ΧΔ(h)=X(s+h)- Χ(s)
όπου το h είναι
παράμετρος
Μία ΤΣ Χ(s) ονομάζεται ομοιόμορφη (ΟΤΣ) όταν η
συνάρτηση διαφορών (~παράγωγος) της είναι ΟΕΕ
με μέση τιμή 0.
13
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Η συνάρτηση ημιβαριογράμματος

Η συνάρτηση αυτή ορίζεται ως

 (h )  E X(s  h)  X(s)2
1
2


Αποδεικνύεται ότι: C(0) - C(h) = γ(h)
Εφόσον το C(0) είναι σταθερά, συμπεραίνεται ότι
κάτω από την υπόθεση ΟΕΕ, η συνδιασπορά και
το ημιβαριόγραμμα είναι δύο ισοδύναμα εργαλεία
για τον χαρακτηρισμό των συσχετίσεων της ΤΣ
14
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Τυχαίος βηματισμός
X(t)
0

T
κ κ κ γ
W(t)
t
γ
γ
κ γ
γ
0
t
Έστω το σύνολο  των ενδεχομένων ενός πειράματος
ρίψης νομίσματος άπειρες φορές. Οι ρίψεις γίνονται κάθε
Τ δευτερόλεπτα, και μετά από κάθε ρίψη γίνεται ένα βήμα
μήκους s με κατεύθυνση προς τα αριστερά εάν έρθει
κορώνα και προς τα δεξιά αν έρθουν γράμματα.
15
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Γιατί (ημι)βαριόγραμμα;

Ο υπολογισμός της συνάρτησης συνδιασποράς
μίας ΟΤΣ δεν είναι πάντοτε εφικτός επειδή δεν
εξασφαλίζεται η ύπαρξη πεπερασμένης C(0).
Πράγματι, υπάρχουν ΟΤΣ οι οποίες παρουσιάζουν
άπειρη δυνατότητα διασποράς όπως αυτή του
Τυχαίου Βηματισμού. Η μελέτη μίας ΟΤΣ γίνεται
υποχρεωτικά
βάσει
της
συνάρτησης
ημιβαριογράμματος, η οποία αντιθέτως μπορεί
πάντα να ορισθεί.
16
Κ. Μόδη: Γεωστατιστική και Εφαρμογές της
(Κεφάλαιο 5)
Άσκηση



Ξεκινώντας
από
το
0,
δημιουργήστε τις 9 υλοποιήσεις
της ΤΣ τυχαίου βηματισμού
όπως στο σχήμα.
Στη συνέχεια, χρησιμοποιώντας
τις βοηθητικές στήλες Μ έως U,
υπολογίστε την διασπορά των
διαφορών Var{Χ(n2T)-X(n1T)}
για n1=0 και παραστήστε την
γραφικά.
Συμφωνούν τα αποτελέσματα
με όσα είπαμε για τις ΟΤΣ;
17