Παρουσίαση 5η: Εισαγωγή στα σήματα · Περιεχόμενα του...

Εφαρμογές Ανάλυσης Σήματος στη Γεωδαισία

Πρόγραμμα Μεταπτυχιακών Σπουδών: Γεωχωρικές τεχνολογίες

Βασίλειος Δ. Ανδριτσάνος Αναπληρωτής Καθηγητής

Τμήμα Πολιτικών Μηχανικών ΤΕ και Μηχανικών Τοπογραφίας και Γεωπληροφορικής ΤΕ Τεχνολογικό Εκπαιδευτικό Ίδρυμα Αθήνας

Παρουσίαση 5η: Εισαγωγή στα σήματα

Περιεχόμενα του μαθήματος (1)

• ΕΝΟΤΗΤΑ 2η Σήματα και ανάλυση Fourier (ΕΡΓΑΣΙΑ 2η)

• Εισαγωγή στα σήματα (Ορισμοί, κατηγορίες σημάτων, βασικά σήματα

συνεχή και διακριτά, κατηγορίες συστημάτων)

• Εισαγωγή στην ανάλυση Fourier (Ολοκληρώματα, σειρές, ιδιότητες,

μιγαδικές εκφράσεις, παραδείγματα, ανάλυση σε συνιστώσες συχνοτήτων,

σειρές στο τετράγωνο και στον κύκλο, παραδείγματα υπολογισμού)

• Μετασχηματισμοί Fourier (Από τα ολοκληρώματα – σειρές στους

μετασχηματισμούς, παραδείγματα, χαρακτηριστικοί μετασχηματισμοί,

ιδιότητες και αποδείξεις)

Περιεχόμενα του μαθήματος (2)

• Διακριτός μετασχηματισμός Fourier (Διαφορές από το συνεχή

μετασχηματισμό, θεώρημα δειγματοληψίας, συχνότητα Nyquist, ιδιότητες,

υπολογισμοί, προβλήματα, φασματική διαρροή, παραποίηση, ταχύς

μετασχηματισμός Fourier – FFT, παραδείγματα)

Βιβλιογραφία

• ΕΝΟΤΗΤΑ 2η

• Hsu, H. P. (1995): Signals and Systems, Schaum’s Outlines eds.

• Proakis, J.G. and D.G. Manolakis (2006): Digital Signal Processing, Fourth ed.,

Pearson, Prentice Hall eds.

• Spiegel, M.R. (1974): Ανάλυση Fourier. Schaum’s Outline Series. McGraw-Hill,

ΕΣΠΙ Αθήνα.

• Brigham, E.O. (1988): The Fast Fourier Transform and its Applications. Prentice

Hall eds.

• Bracewell, R.N. (1978): The Fourier Transform and its applications. McGraw-Hill

eds.

Περιεχόμενα παρουσίασης • Σήματα και γεωεπιστήμες

• Ορισμοί – κατηγορίες σημάτων

• Βασικά συνεχή και διακριτά σήματα και τρόποι κατασκευής τους

• Συστήματα και κατηγορίες συστημάτων

• Η αξιοποίηση των φασματικών τεχνικών στις γεωεπιστήμες

Σήματα και γεωεπιστήμες • Σήματα στις γεωεπιστήμες τελευταία 30ετία

• Οι μαθηματικές αναλύσεις της επιστήμης των ηλεκτρονικών εφαρμόζονται

και στην περιγραφή φυσικών συστημάτων

• Η αντιμετώπιση των παρατηρήσεων ως σήματα δίνει τη δυνατότητα

απλούστερων μαθηματικών υπολογισμών ταχύτεροι υπολογισμοί και

διαχείριση μεγάλου όγκου δεδομένων

• Τα πεδία των σημάτων μπορεί να είναι από απλές χρονοσειρές έως

τρισδιάστατα πεδία σημάτων περιγραφής του χώρου

Σήματα και γεωεπιστήμες

Electromagnetic Distance Measurement

Satellite Laser Ranging

Satellite Altimetry

Global Navigation Satellite Systems

Very Long Baseline Interferometry

Σήματα και γεωεπιστήμες

Geophysics and Geodynamics Geoid estimation

Ορισμοί

• Η έννοια του σήματος και του συστήματος είναι απαραίτητη σχεδόν σε

κάθε κλάδο της επιστήμης του ηλεκτρονικού μηχανικού, αλλά και σε

σύγχρονες εφαρμογές της επιστήμης της Γεωπληροφορικής

• Σήμα συνάρτηση που αναπαριστά μία φυσική ποσότητα ή μεταβλητή και

τυπικά περιέχει πληροφορίες σχετικά με τη συμπεριφορά ή τη φύση του

φαινομένου

Ανεξάρτητη μεταβλητή συνήθως χρόνος (αλλά και συντεταγμένες χώρου – μήκος)

Κατηγορίες σημάτων

• Ανάλογα με τη μορφή της ανεξάρτητης μεταβλητής t

• Συνεχές σήμα (continuous-time signal) ανεξάρτητη μεταβλητή

συνεχούς μορφής x(t)

• Διακριτό σήμα (discrete-time signal) ανεξάρτητη μεταβλητή t διακριτής

μορφής x[n] ή xn (n ακέραιος διακριτός αριθμός)

• Εφόσον ένα διακριτό σήμα ορίζεται σε διακριτές τιμές του t συνήθως

ταυτοποιείται ως μία σειρά αριθμών (sequence of numbers)

Κατηγορίες σημάτων

• Ένα διακριτό σήμα αναπαριστά ένα φυσικό φαινόμενο με διακριτή

ανεξάρτητη μεταβλητή t x[n], π.χ. θερμοκρασία τις στιγμές της μέτρησης

• Ένα διακριτό σήμα μπορεί να προκύψει από ένα συνεχές με μία διαδικασία

που ονομάζεται διακριτοποίηση ή δειγματοληψία (sampling)

• Οι διακριτές τιμές κατά τη μετατροπή από συνεχές σε διακριτό ονομάζονται

δειγματικές τιμές (samples) και το διάστημα των δειγμάτων διάστημα

δειγματοληψίας (sampling interval)

Κατηγορίες σημάτων • Ανάλογα με διαστάσεις της ανεξάρτητης μεταβλητής

• Μονοδιάστατα σήματα (one-dimensional signals) συνάρτηση μίας

ανεξάρτητης μεταβλητής, π.χ., χρονοσειρές παλίρροιας

• Διδιάστατα σήματα (two-dimentional signals) συνάρτηση δύο

ανεξάρτητων μεταβλητών, π.χ., εικόνα, πεδίο ανωμαλιών βαρύτητας, κ.λπ.

Παλίρροιες στερεού φλοιού Ανωμαλίες ελευθέρου αέρα

Κατηγορίες σημάτων • Ανάλογα με τις συνιστώσες της συνάρτησης του σήματος (εξαρτημένη

μεταβλητή)

• Πολυκάναλα σήματα (multichannel signals) σήματα που μπορούν να

διακριθούν σε συνιστώσες, π.χ., μελέτη σεισμικών κυμάτων μετά από

σεισμό: συνιστώσες στον άξονα ανατολής – δύσης, βορρά – νότου (S-

waves) και εγκάρσια συνιστώσα (P-wave)

Κατηγορίες σημάτων • Ανάλογα με τις τιμές που λαμβάνει στο πεδίο ορισμού του

• Ένα συνεχές σήμα x(t) λαμβάνει οποιαδήποτε τιμή στο συνεχές διάστημα (a,

b) αναλογικό σήμα (analog signal) (a και b μπορεί -∞, + ∞)

• Ένα διακριτό σήμα x[n] που μπορεί να λάβει συγκεκριμένες τιμές στο

διάστημα ορισμού του ψηφιακό σήμα (digital signal)

Κατηγορίες σημάτων • Ανάλογα με το χώρο των συναρτήσεων, στις οποίες αντιστοιχεί

• Πραγματικό σήμα (real signal) οι τιμές της συνάρτησης ανήκουν στο

χώρο των πραγματικών αριθμών

• Μιγαδικό σήμα (complex signal) οι τιμές της συνάρτησης x(t) ανήκουν

στο χώρο των μιγαδικών αριθμών (πραγματικό και φανταστικό μέρος)

( ) ( ) ( )[ ] [ ] [ ]njxnxnx

tjxtxtx

21

21

+=+= Πραγματικά

σήματα

1−=j

Κατηγορίες σημάτων • Ανάλογα με τη στατιστική συμπεριφορά τους

• Ντετερμινιστικά σήματα (deterministic signals) οι τιμές τους είναι

καθορισμένες πλήρως για οποιαδήποτε τιμή της ανεξάρτητης μεταβλητής t

• Στοχαστικά ή τυχαία σήματα (stochastic or random signals) οι τιμές

που λαμβάνουν είναι τυχαίες και είναι γνωστή μόνο η στατιστική

συμπεριφορά τους

Προκύπτει από τιμές συγκεκριμένης συνάρτησης

Γνωστή στατιστική συμπεριφορά (μ, σ2)

Κατηγορίες σημάτων • Ανάλογα με την απόκρισή τους στην αλλαγή προσήμου της μεταβλητής t ή τη

συμμετρία τους ως προς τον άξονα x(t) ή x[n]

• Άρτια – συμμετρικά σήματα (even – symmetric signals) όταν ισχύει

• Περιττά – αντισυμμετρικά σήματα (odd – antisymmetric signals) όταν

ισχύει

( ) ( ) [ ] [ ]nxnxtxtx =−=− ,

( ) ( ) [ ] [ ]nxnxtxtx −=−−=− ,

Κάθε σήμα μπορεί να εκφραστεί ως το άθροισμα

ενός άρτιου και ενός περιττού σήματος

( ) ( ) )(txtxtx oe +=

Το γινόμενο δύο άρτιων ή δύο περιττών σημάτων

είναι άρτιο σήμα ενώ το γινόμενο ενός άρτιο με

ένα περιττό σήμα είναι περιττό σήμα

Κατηγορίες σημάτων • Ανάλογα με την περιοδικότητά εμφάνισής τους

• Περιοδικά σήματα (periodic signals) συνεχές σήμα x(t) καλείται

περιοδικό με περίοδο Τ όταν Τ > 0 και , για κάθε t

• Ισχύει επίσης:

• Θεμελιώδης περίοδος (fundamental period) To η μικρότερη θετική τιμή

του Τ για την οποία ισχύει

( ) ( )txTtx =+

( ) ( )txmTtx =+


• Περιοδικά σήματα (periodic signals) διακριτό σήμα x[n] καλείται

περιοδικό με περίοδο N όταν N > 0 και , για κάθε n

• Ισχύει επίσης:

• Θεμελιώδης περίοδος Νο ο μικρότερος θετικός ακέραιος Ν για τον οποίο

ισχύει

[ ] [ ]nxNnx =+

[ ] [ ]nxmNnx =+


• Ιδιότητες περιοδικότητας:

• Μία σειρά, η οποία λαμβάνεται με τη διαδικασία ομογενούς δειγματοληψίας

(uniform sampling) ενός περιοδικού συνεχούς σήματος μπορεί να μην είναι

περιοδική

• Το άθροισμα δύο συνεχών περιοδικών σημάτων μπορεί να μην είναι

περιοδικό, ενώ το άθροισμα δύο περιοδικών σειρών είναι πάντοτε περιοδικό

Κατηγορίες σημάτων • Σήματα Ενέργειας και Ισχύος για συνεχή σήματα

• Για αυθαίρετο συνεχές σήμα x(t)

• Η ομαλοποιημένη ενέργεια Ε (normalized energy) ενός σήματος

• Η ομαλοποιημένη μέση ισχύς P (normalized average power) ενός

σήματος

( )∫∞

∞−

= dttxE 2

( )∫−

∞→=

2/

2/

21limT

TT

dttxT

P

Κατηγορίες σημάτων • Σήματα Ενέργειας και Ισχύος για διακριτά σήματα (σειρές)

• Για αυθαίρετο διακριτό σήμα x[n]

• Η ομαλοποιημένη ενέργεια Ε (normalized energy) ενός σήματος

• Η ομαλοποιημένη μέση ισχύς P (normalized average power) ενός

σήματος

[ ]∑∞

−∞=

=n

nxE 2

[ ]∑−=

∞→ +=

N

NnN

nxN

P 2

121lim

Κατηγορίες σημάτων • Σήματα Ενέργειας και Ισχύος

• Με τη βοήθεια των παραπάνω ορισμών μετριέται το «μέγεθος» ενός

σήματος

• Η ενέργεια ενός σήματος δίνεται από το εμβαδόν της περιοχής που καλύπτει

το τετράγωνο της συνάρτησης

• Αρνητική τιμή της συνάρτησης δεν αφαιρεί ενέργεια για τον υπολογισμό

λαμβάνεται η απόλυτη τιμή

Κατηγορίες σημάτων • Το πρόβλημα με τον υπολογισμό της ενέργειας ενός σήματος παρουσιάζεται

όταν το σήμα δεν διακόπτεται άπειρη ενέργεια

• Στη περίπτωση αυτή η ενέργεια θεωρείται ακατάλληλη για το υπολογισμό του

«μεγέθους» χρήση της ισχύος

• Η ισχύς αντιπροσωπεύει την ενέργεια ανά μονάδα χρόνου (ή χώρου)

Κατηγορίες σημάτων • Κατηγορίες σημάτων ενέργειας και ισχύος

• Ένα σήμα ονομάζεται σήμα ενέργειας (energy signal) όταν και μόνο όταν

• Ένα σήμα ονομάζεται σήμα ισχύος (power signal) όταν και μόνο όταν

• Σημείωση: ένα περιοδικό σήμα είναι ένα σήμα ισχύος εάν η περιεχόμενη

ενέργεια ανά περίοδο είναι πεπερασμένη και η μέση ισχύς του χρειάζεται

να υπολογιστεί μόνο μέσα σε μία περίοδο

00 =⇒∞<< PE

∞=⇒∞<< EP0

Βασικά συνεχή σήματα Μοναδιαία συνάρτηση βηματισμού – Unit step function (Heaviside)

• Ορισμός

( )

<>

=0001

tt

tu ( )

<>

=−o

oo tt

ttttu

01

Βασικά συνεχή σήματα Μοναδιαία συνάρτηση βηματισμού – Unit step function (Heaviside)

• Εντολές MATLAB – octave (packages signal, control, symblink)

( )

<>

=0001

tt

tu ( )

<>

=−o

oo tt

ttttu

01

>> n=-10:0.001:10; >> y=heaviside(n); >> plot(n,y)

>> n=-10:0.001:10; >> y=heaviside(n-3); >> plot(n,y)

Βασικά συνεχή σήματα Μοναδιαία συνάρτηση παλμού – Unit impulse function (Dirac delta)

• Είναι γνωστή και ως συνάρτηση δέλτα του Dirac δ(t). Αποτελεί μία σημαντική

συνάρτηση στην ανάλυση συστημάτων

• δ(t) το όριο μίας κατάλληλα επιλεγμένης συμβατικής συνάρτησης με

εμβαδό ίσο με τη μονάδα

ε

0→ε

ε1

( )

=∞≠

=δ000

tt

t

( ) 1=δ∫ε

ε−

dtt


• Μαθηματικά ορίζεται και ως:

• Η συνάρτηση υστέρησης δέλτα (delayed delta function) ορίζεται

σύμφωνα με τα παραπάνω

( ) ( ) ( )0φ=δφ∫∞

∞−

dttt

( ) ( ) ( )00 tdtttt φ=−δφ∫∞

∞−


• Η παράγωγος της μοναδιαίας συνάρτησης βηματισμού ως προς την

μεταβλητή t είναι η συνάρτηση δέλτα του Dirac

• Κάποιες ιδιότητες της συνάρτησης δέλτα

( ) ( ) ( )dt

tdutut =′=δ

>> dirac(1) ans = 0 >> dirac(0) ans = Inf >>


( ) ( )ta

at δ=δ1

( ) ( )tt δ=−δ ( ) ( ) ( ) ( )000 tttxtttx −δ=−δ

( ) ( ) ( ) ( )txttx δ=δ 0

Βασικά συνεχή σήματα Μιγαδικό εκθετικό σήμα – Complex exponential signal

• Ορίζεται ως

• Πολύ σημαντικό σήμα στις αναλύσεις Fourier

• Χρησιμοποιώντας τον τύπο του Euler είναι δυνατό να οριστεί και ως

• Είναι ένα μιγαδικό σήμα του οποίου το πραγματικό μέρος είναι το cosω0t και

το φανταστικό μέρος το sinω0t

• Το μιγαδικό εκθετικό σήμα είναι περιοδικό. Η θεμελιώδης περίοδός του είναι

( ) tjetx 0ω=

( ) tjtetx tj00 sincos0 ω+ω== ω

00

2ωπ

=T


• Γενική μορφή: έστω s = σ + jω

• Όταν σ > 0 το σήμα ονομάζεται εκθετικά αυξανόμενο ημιτονοειδές

(συνημιτονοειδές) – exponentially increasing sinusoidal signal

• Όταν σ < 0 το σήμα ονομάζεται εκθετικά ελαττούμενο ημιτονοειδές

(συνημιτονοειδές) – exponentially decreasing sinusoidal signal

( ) ( ) ( ) tjetetjteeetx ttttjst ω+ω=ω+ω=== σσσω+σ sincossincos

Πραγματικό μέρος Φανταστικό μέρος


• Παράδειγμα

Βασικά συνεχή σήματα Πραγματικό εκθετικό σήμα – Real exponential signal

• Πραγματικό εκθετικό σήμα (s = σ πραγματικός αριθμός)

( ) tetx σ=

0>σ

0<σ


• Παράδειγμα σχεδίασης μιγαδικού εκθετικού σήματος

( ) ( ) 100010/ ≤≤= π tetx tj

>> t=0:0.1:100; >> w=2*pi/20; >> x=exp(i*w*t); >> plot(t,real(x))

>> t=0:0.1:100; >> w=2*pi/20; >> x=exp(i*w*t); >> plot(t,imag(x))

( )t10/cos π

( )t10/sin π

Βασικά συνεχή σήματα Εκθετικά αυξανόμενο ημιτονοειδές σήμα

• Παράδειγμα σχεδίασης εκθετικά αυξανόμενου ημιτονοειδούς σήματος

( ) ( ) ( )( ) 100010/20/1 ≤≤= π+ tetx tj

>> t=0:0.1:100; >> w=2*pi/20; >> x=exp((1/20+i*w)*t); >> plot(t,real(x))

>> t=0:0.1:100; >> w=2*pi/20; >> x=exp((1/20+i*w)*t); >> plot(t,imag(x))

Βασικά συνεχή σήματα Εκθετικά ελαττούμενο ημιτονοειδές σήμα

• Παράδειγμα σχεδίασης εκθετικά ελαττούμενου ημιτονοειδούς σήματος

( ) ( ) ( )( ) 100010/20/1 ≤≤= π+− tetx tj

>> t=0:0.1:100; >> w=2*pi/20; >> x=exp((-1/20+i*w)*t); >> plot(t,real(x))

>> t=0:0.1:100; >> w=2*pi/20; >> x=exp((-1/20+i*w)*t); >> plot(t,imag(x))

Βασικά συνεχή σήματα Πραγματικό εκθετικό σήμα – Real exponential signal

• Παράδειγμα σχεδίασης πραγματικού εκθετικού σήματος

( ) ( ) 100020/1 ≤≤= tetx t

>> t=0:0.1:100; >> x=exp((1/20)*t); >> plot(t,x)

>> t=0:0.1:100; >> x=exp((-1/20)*t); >> plot(t,x)

Βασικά συνεχή σήματα Ημιτονοειδές σήμα – Sinusoidal signal

• Ένα συνεχές ημιτονοειδές σήμα (sinusoid) έχει τη μορφή

• A εύρος σήματος (amplitude)

• ω0 γωνιακή συχνότητα (σε ακτίνια ανά δευτερόλεπτο)

• θ γωνία φάσης (σε ακτίνια)

( ) ( )θ+ω= tAtx 0cos

Σήμα ημιτόνου

Βασικά συνεχή σήματα Ημιτονοειδές σήμα – Sinusoidal signal

• Το συνεχές ημιτονοειδές σήμα είναι περιοδική συνάρτηση με θεμελιώδη

περίοδο

• Το αντίστροφο της θεμελιώδους περιόδου ονομάζεται θεμελιώδης συχνότητα

• Χρησιμοποιώντας την ταυτότητα του Euler

• Ισχύει αντιστοίχως

00

2ωπ

=T

)(10

0 HzhertzT

f =

( ) ( ) θ+ω⋅=θ+ω tjerealAtA 00cos

( ) ( )θ+ω=⋅ θ+ω tAeimaginaryA tj0sin0

Στόχοι φασματικής ανάλυσης

• Κυματοειδής μορφή σήματος δυνατότητα ανάλυσης σε συχνότητες

• Ο χώρος των συχνοτήτων επιτρέπει ευκολότερους υπολογισμούς

πολύπλοκες συναρτήσεις αναλύονται σε απλής μορφής διαγράμματα

συχνοτήτων

Σημαντικό να γνωρίζουμε για

κάθε μέτρηση – κύμα τις

κυρίαρχες συχνότητές της

φάσμα (spectrum) της

μέτρησης

Βασικές έννοιες κυματικής μορφής σήματος

• Από το πεδίο του χρόνου/χώρου στο πεδίο των συχνοτήτων

t

Σήμα f(t) μιας μοναδικής συχνότητας (μονοχρωματικό):



t


sin cos

Επανάληψη ημιτόνων ή συνημιτόνων



t


sin cos

Επανάληψη ημιτόνων ή συνημιτόνων

κύκλος (= 1 επανάληψη) ημιτόνου κύκλος (= 1 επανάληψη) συνημιτόνου



t


Συχνότητα (Frequency) f: αριθμός κύκλων στη μονάδα χρόνου (κύκλοι ανά δευτερόλεπτο)

Τ

Τ

Περίοδος (Period) Τ: χρόνος που χρειάζεται για να επαναληφθεί ένας κύκλος

Μήκος κύματος (wavelength) λ = cT: διάστημα που διανύει το σήμα σε μία περίοδο c : ταχύτητα φωτός (ηλεκτρομαγνητικής ακτινοβολίας) στο κενό

1 sec (f = 2.5 κύκλοι ανά δευτερόλεπτο)


χαμηλή συχνότητα

υψηλή συχνότητα


μεγάλη περίοδος

μικρή περίοδος


μεγάλο μήκος κύματος

μικρό μήκος κύματος

Από τα συνεχή στα διακριτά

• Η διακριτοποίηση συνεχούς σήματος οδηγεί σε σειρές αριθμών

(sequences of numbers)

• Η δημιουργία ενός διακριτού σήματος είναι δυνατή:

1. Με τη διαδικασία της δειγματοληψίας (sampling) από ένα συνεχές

σήμα σε συγκεκριμένες χρονικές στιγμές (ή χωρικά διαστήματα)

μετρητικά όργανα με μετρήσεις ανά συγκεκριμένες χρονικές περιόδους

ή χωρικά διαστήματα

2. Με τη συσσώρευση μίας μεταβλητής μέσα σε μία χρονική περίοδο,

π.χ. μέτρηση κυκλοφορίας ανά ώρα

Από τα διακριτά στα ψηφιακά

• Το διακριτό σήμα λαμβάνει τιμές ανά συγκεκριμένα διαστήματα (χρονικά ή

χωρικά) και συνήθως ισοδιαστήματα

• Όταν το διακριτό σήμα λαμβάνει συγκεκριμένες τιμές (εξαρτημένη μεταβλητή

συγκεκριμένων) τιμών τότε ονομάζεται ψηφιακό σήμα (digital signal)

• Για να είναι δυνατή η ψηφιακή επεξεργασία ενός σήματος πρέπει να είναι

διακριτό στο χρόνο (χώρο) και να έχει συγκεκριμένες διακριτές τιμές

• Η διαδικασία «ψηφιοποίησης» διακριτού σήματος ονομάζεται

κβαντοποίηση (quantization)

Βασικά διακριτά σήματα Μοναδιαία σειρά βηματισμού – Unit Step Sequence

• Η μοναδιαία σειρά βηματισμού u[n] ορίζεται ως:

• Σε αντίθεση με το συνεχές σήμα η σειρά ορίζεται και για n = 0

[ ]

<≥

=0001

nn

nu


• Η μετατοπισμένη μοναδιαία σειρά βηματισμού (shifted unit step

sequence) u[n – k] ορίζεται ως

[ ]

<≥

=−knkn

knu01



• Η συνάρτηση heaviside που χρησιμοποιείται στα συνεχή σήματα δε δίνει

σωστά αποτελέσματα στα διακριτά γιατί ορίζεται σύμφωνα με τα

συνεχή δε δίνει τιμή στο n = 0

>> n=-5:5; >> u=heaviside(n); >> stem(n,u) >> axis([-5 5 0 2])


• Ορισμός συνάρτησης σύμφωνα με τα προηγούμενα

>> function u=step(n) u=0*n; u(find(n>=0))=1; end

>> n=-5:5; >> u=step(n); >> stem(n,u) >> axis([-5 5 -1 2])

[ ]

<≥

=0001

nn

nu

>> n=-5:5; >> u=step(n-3); >> stem(n,u) >> axis([-5 5 0 2])

[ ]

<≥

=−3031

3nn

nu

Βασικά διακριτά σήματα Μοναδιαία σειρά παλμού – Unit Impulse Sequence

• Ορίζεται και ως μοναδιαία δειγματική σειρά (unit sample sequence)

• Σε αντίθεση με το αντίστοιχο συνεχές σήμα ορίζεται χωρίς μαθηματική

δυσκολία απουσιάζει η τιμή του απείρου

• Η σύνδεση μεταξύ των σειρών βηματισμού και παλμού δίνεται από:

[ ]

≠=

=δ0001

nn

n

[ ] [ ] [ ]

[ ] [ ]∑−∞=

δ=

−−=δn

k

knu

nunun 1

Βασικά διακριτά σήματα Μετατοπισμένη μοναδιαία σειρά παλμού – Shifted Unit Impulse Sequence

• Ορίζεται και ως μετατοπισμένη μοναδιαία δειγματική σειρά (shifted unit

sample sequence)

• Κάποιες από τις ιδιότητες του συνεχούς σήματος όπως εκφράζονται ως

διακριτές

[ ]

≠=

=−δknkn

kn01

[ ] [ ] [ ] [ ][ ] [ ] [ ] [ ]knkxknnx

nxnnx−δ=−δ

δ=δ 0 Αρχή θεωρήματος δειγματοληψίας



• Η συνάρτηση dirac που χρησιμοποιείται στα συνεχή σήματα δε δίνει

σωστά αποτελέσματα στα διακριτά γιατί ορίζεται σύμφωνα με τα

συνεχή δίνει τιμή ∞ στο n = 0

>> n=-5:5 n = -5 -4 -3 -2 -1 0 1 2 3 4 5 >> dirac(n) ans = 0 0 0 0 0 Inf 0 0 0 0 0 >>



>> t=-5:5 >> impulse = t==0; >> stem(t,impulse) >> axis[(-5 5 0 2)]

>> t=-5:5 >> unitstep = t>=0; >> stem(t,unitstep) >> axis[(-5 5 0 2)]

[ ]

≠=

=δ0001

nn

n

Εναλλακτικός τρόπος κατασκευής μοναδιαίας σειράς βηματισμού

Βασικά διακριτά σήματα Μιγαδική εκθετική σειρά – Complex Exponential Sequence

• Ορίζεται ως

• Ακολουθώντας την ταυτότητα του Euler

• Για να είναι περιοδική σειρά με περίοδο Ν > 0 πρέπει

• Η μιγαδική εκθετική σειρά δεν είναι περιοδική για κάθε Ω0. Είναι περιοδική

μόνο όταν Ω0 /2π είναι ρητός αριθμός σε αντίθεση με το συνεχές μιγαδικό

εκθετικό σήμα που είναι περιοδικό για κάθε τιμή του ω0

[ ] njenx 0Ω=

[ ] njnenx nj00 sincos0 Ω+Ω== Ω Φανταστικό μέρος Πραγματικό μέρος

Nm

=Ωπ2

0 Θετικός ακέραιος

Βασικά διακριτά σήματα Μιγαδική εκθετική σειρά – Complex Exponential Sequence

• Η θεμελιώδης περίοδος για την σειρά x[n] με την προϋπόθεση της

περιοδικότητας που παρουσιάστηκε

Ω

=0

0

2πmN

Βασικά διακριτά σήματα Μιγαδική εκθετική σειρά – Complex exponential sequence

• Παράδειγμα σχεδίασης μιγαδικού εκθετικού σήματος

[ ] ( ) 100010/ ≤≤= nenx nj π

>> n=0:1:100; >> w=2*pi/20; >> x=exp(i*w*n); >> stem(n,imagl(x))

( )n10/cos π

( )n10/sin π

>> n=0:1:100; >> w=2*pi/20; >> x=exp(i*w*n); >> stem(n,real(x))

Βασικά διακριτά σήματα Εκθετικά αυξανόμενη ημιτονοειδής σειρά

• Παράδειγμα σχεδίασης εκθετικά αυξανόμενης ημιτονοειδούς σειράς

[ ] ( ) ( )( ) 100010/20/1 ≤≤= + nenx nj π

>> n=0:1:100; >> w=2*pi/20; >> x=exp((1/20+i*w)*n); >> stem(n,real(x))

>> n=0:1:100; >> w=2*pi/20; >> x=exp((1/20+i*w)*n); >> stem(n,imag(x))

Βασικά διακριτά σήματα Εκθετικά ελαττούμενη ημιτονοειδής σειρά

• Παράδειγμα σχεδίασης εκθετικά ελαττούμενης ημιτονοειδούς σειράς

[ ] ( ) ( )( ) 100010/20/1 ≤≤= +− nenx nj π

>> n=0:1:100; >> w=2*pi/20; >> x=exp((-1/20+i*w)*n); >> stem(n,real(x))

>> n=0:1:100; >> w=2*pi/20; >> x=exp((-1/20+i*w)*n); >> stem(n,imag(x))

Βασικά διακριτά σήματα Εκθετικά ελαττούμενη ημιτονοειδής σειρά

• Τα μιγαδικά σήματα και οι σειρές είναι δυνατό να παρουσιαστούν και σε

τρισδιάστατα σχήματα

[ ] ( ) ( )( ) 100010/20/1 ≤≤= +− nenx nj π

>> n=0:1:100; >> w=2*pi/20; >> x=exp((-1/20+i*w)*n); >> stem3(n,real(x),imag(x))

>> n=0:0.1:100; >> w=2*pi/20; >> x=exp((-1/20+i*w)*n); >> plot3(n,x)

Βασικά διακριτά σήματα Πραγματική εκθετική σειρά – Real exponential sequence

• Παράδειγμα σχεδίασης πραγματικού εκθετικού σήματος

[ ] ( ) 100020/1 ≤≤= nenx n

>> n=0:1:100; >> x=exp((1/20)*n); >> stem(n,x)

>> n=0:1:100; >> x=exp((-1/20)*n); >> stem(n,x)

Βασικά διακριτά σήματα Ιδιότητες μιγαδικής εκθετικής σειράς

• Μελετώντας μία εκθετική σειρά ισχύει:

• Η μιγαδική εκθετική σειρά είναι περιοδική στο διάσημα 0 ≤ Ω0 ≤ 2π

• Μελετώντας μία τέτοια σειρά μπορούμε να περιοριστούμε στο διάστημα

0 ≤ Ω0 ≤ 2π, μέσα στο οποίο επιλέγεται η συχνότητα Ω0

• Η ιδιότητα αυτή χρησιμοποιείται στους μετασχηματισμούς, όπου περιοδικές

συναρτήσεις είναι δυνατό να περιγράψουν μη περιοδικά φαινόμενα

( ) njknjnjnkj eeee 000 22 ΩπΩπ+Ω ==

Βασικά διακριτά σήματα Ημιτονοειδής σειρά – Sinusoidal Sequence

• Ορίζεται ως:

• Αν το n είναι αδιάστατο τότε το Ω0 και το θ έχουν μονάδες ακτίνια

• Η ημιτονονειδής σειρά μπορεί να εκφραστεί σε συνάρτηση με την εκθετική

μιγαδική σειρά, όπως και στα συνεχή σήματα

• Ισχύουν οι ίδιες ιδιότητες και απαιτήσεις περιοδικότητας, όπως και στην

εκθετική μιγαδική σειρά

• Εκθετικές και ημιτονειδείς σειρές είναι οι βασικές σειρές περιγραφής των

γεωδαιτικών δεδομένων κατά τη διαδικασία των διακριτών μετασχηματισμών

Fourier

[ ] ( )θ+Ω= nAnx 0cos

( ) ( ) θ+Ω⋅=θ+Ω njerealAnA 00cos

Nm

=Ωπ2

0

Βασικά διακριτά σήματα Ημιτονοειδής σειρά – Sinusoidal Sequence

• Παράδειγμα ημιτονοειδούς σειράς

Περιοδική Θεμελιώδης περίοδος Ν0 = 12

Μη περιοδική άρρητο

Nm

≠π

=π

=π

Ω41

22/1

20

121

26/

20 =

ππ

==π

ΩNm

Συστήματα και κατηγοριοποίηση

• Σύστημα στις γεωεπιστήμες καλείται ένα μαθηματικό μοντέλο περιγραφής

ενός φυσικού φαινομένου

• Το σύστημα σχετίζει τα σήματα εισόδου ή διέγερσης (input or excitation)

με τα σήματα εξόδου ή απόκρισης (output or response)

• Ένα σύστημα μπορεί να θεωρηθεί ως μετασχηματισμός ή απεικόνιση ή

φιλτράρισμα του σήματος εισόδου στο σήμα εξόδου

Γωνίες Αποστάσεις

Μαθηματικό μοντέλο Ευκλείδειας Γεωμετρίας Φυσικό

Σύστημα

Συντεταγμένες


• Εάν θεωρηθεί ως y το σήμα εισόδου (παρατήρηση στις γεωεπιστήμες) και

ως x το σήμα εξόδου (άγνωστο στις γεωεπιστήμες) τότε ο μετασχηματισμός

yx T=Τελεστής αναπαράστασης ενός καλώς ορισμένου κανόνα μετασχηματισμού x y

Cretan gyre

Ierapetra Anticyclone

Rhodes gyre and Mid-Mediterranean Jet

Mersa-Matruh Anticyclone

Παράδειγμα συστήματος πολλαπλής εισόδου – εξόδου με θόρυβο

Εκτίμηση Δυναμικής Θαλάσσιας Τοπογραφίας από δεδομένα αλτιμετρικών δορυφόρων και γεωδυναμικά

μοντέλα

Andritsanos, V. D. and I. N. Tziavos (2016): Quasi-Stationary SST Estimation in the Eastern Mediterranean Sea using marine gravity, GOCE/GRACE

gravity information and recent altimetry missions through the Multiple Input / Multiple Output System Theory, European Space Agency (ESA) Living

Planet Symposium, Prague


• Ανάλογα με τη φύση του σήματος εισόδου και εξόδου

• Σύστημα συνεχών σημάτων χρόνου/χώρου (continuous-time system)

• Σύστημα διακριτών σημάτων χρόνου/χώρου (discrete-time system)


• Ανάλογα με την εξάρτηση από τη ανεξάρτητη μεταβλητή

• Σύστημα χωρίς μνήμη το σήμα εξόδου σε κάθε χρονική στιγμή (χωρική

συντεταγμένη) εξαρτάται μόνο από το σήμα εισόδου της ίδιας χρονικής

στιγμής (χωρικής συντεταγμένης), π.χ. το ρεύμα που εισέρχεται σε μία

αντίσταση και η τάση που εξέρχεται

• Σύστημα με μνήμη κάθε σύστημα που το σήμα εισόδου επιδρά στο

σύνολο του σήματος εξόδου


• Ανάλογα με τη γραμμικότητα του τελεστή του συστήματος

• Γραμμικός τελεστής Τ και γραμμικό σύστημα όταν ισχύουν δύο ιδιότητες 1. Προσθετικότητα: για κάθε Τy1 = x1 και Τy2 = x2 Τy1 + y2 = x1 + x2

2. Ομογένεια: Ταx=αy

• Kάθε σύστημα που δεν ικανοποιεί τις παραπάνω ιδιότητες ονομάζεται

μη γραμμικό, π.χ.

2yx =yx cos=


• Ανάλογα με την επίδραση της αλλαγής της ανεξάρτητης μεταβλητής

• Χρονικά (ή χωρικά) αμετάβλητα συστήματα (time (space) invariant)

εάν μία μεταβολή στην κλίμακα του χρόνου (χώρου) του σήματος εισόδου

προκαλεί την ίδια μεταβολή στο σήμα εξόδου

• Συνεχή χρονικά (χωρικά) αμετάβλητα συστήματα

• Διακριτά χρονικά (χωρικά) αμετάβλητα συστήματα

( ) ( )τ−=τ− txtyT

[ ] [ ]knxkny −=−T

Η αξιοποίηση των φασματικών τεχνικών

• Πλήθος δεδομένων γεωεπιστημών επίγεια, από αέρα, θαλάσσια, δορυφορικά

δεδομένα τροποποίηση μεθόδων επεξεργασίας

• Οι διαφορετικές πηγές δεδομένων αντιστοιχούν σε διαφορετικές περιοχές του

φάσματος του πεδίου που αναλύεται ανάγκη συνδυασμού ετερογενών δεδομένων,

π.χ., χαμηλές συχνότητες δορυφορικών δεδομένων

• Τροποποίηση κλασικών φασματικών τεχνικών ώστε να συμπεριλαμβάνεται στη

διαδικασία κατάλληλα κριτήρια ελαχιστοποίησης των σφαλμάτων της εκτίμησης

• Η αντιμετώπιση των μαθηματικών μοντέλων ως συστημάτων είναι δυνατό να οδηγήσει

στην κατάλληλη τροποποίηση τους με την ελαχιστοποίηση του λόγου σήμα / θόρυβος

Ανακεφαλαίωση

• Εισαγωγή στα σήματα – Ορισμοί

• Κατηγορίες σημάτων

• Παραδείγματα βασικών συνεχών και διακριτών σημάτων

• Μέθοδοι κατασκευής σήματος

• Συστήματα και κατηγοριοποίηση

• Η αξιοποίηση των φασματικών τεχνικών στις γεωεπιστήμες

Παρουσίαση 5η: Εισαγωγή στα σήματα · Περιεχόμενα του...

Documents

Transcript of Παρουσίαση 5η: Εισαγωγή στα σήματα · Περιεχόμενα του...