01.11.2014 Views

TEORIJA SIGNALOV - Laboratorij za obdelavo signalov in daljinska ...

TEORIJA SIGNALOV - Laboratorij za obdelavo signalov in daljinska ...

TEORIJA SIGNALOV - Laboratorij za obdelavo signalov in daljinska ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

UNIVERZA<br />

V<br />

MARIBORU<br />

Žarko ČUČEJ <strong>in</strong> Peter PLANINŠIČ<br />

<strong>TEORIJA</strong> <strong>SIGNALOV</strong><br />

uvod


CIP - Kataloški <strong>za</strong>pis o publikaciji<br />

Univerzitetna knjižnica Maribor<br />

681.51/.52(075.8)<br />

ČUČEJ, Žarko<br />

Teorija <strong>signalov</strong>: uvod / Žarko Čučej <strong>in</strong> Peter Plan<strong>in</strong>šič ;<br />

[risbe Žarko Čučej]. - Maribor: Fakulteta <strong>za</strong><br />

elektrotehniko, računalništvo <strong>in</strong> <strong>in</strong>formatiko, 2000<br />

ISBN 86-435-0330-4<br />

1. Plan<strong>in</strong>šič, Peter<br />

COBBIS-ID 44994817<br />

naslov <strong>TEORIJA</strong> <strong>SIGNALOV</strong>: uvod revizija:20010328<br />

avtor<br />

Žarko ČUČEJ<br />

soavtor<br />

Peter PLANINŠIČ<br />

recenzija<br />

izr. prof. dr. Rajko SVEČKO<br />

doc. dr. Gorazd LEŠNJAK<br />

jezik<br />

prof. Milena MILANOVIČ<br />

uredil <strong>in</strong> oblikoval Žarko ČUČEJ<br />

risbe<br />

Žarko ČUČEJ<br />

uporabljani programi MiTEX1.20e z AMS-TEX, CorelDraw 7<br />

<strong>za</strong>ložba<br />

Univer<strong>za</strong> v Mariboru<br />

Fakulteta <strong>za</strong> elektrotehniko, računalništvo <strong>in</strong> <strong>in</strong>formatiko 28. marec 2001<br />

vse pravice pridržane


Predgovor<br />

Prva knjiga je pošla. Da bi vsaj delno premostila to vrzel, sva<br />

se odločila objaviti datoteko z njeno vseb<strong>in</strong>o na domači strani<br />

laboratorija:<br />

http://SPaRC.feri.uni-mb/publikacije.<br />

Datoteka se od orig<strong>in</strong>ala razlikuje v odpravljenih napakah, ki<br />

sva jih avtorja odkrila tudi s pomočjo študentov. Med njimi se<br />

je posebej izka<strong>za</strong>l Uroš Lebar, <strong>za</strong> kar se mu posebej <strong>za</strong>hvaljujeva.<br />

Je pa tudi nekaj razlik v organi<strong>za</strong>ciji poglavij (razlaga korelacije)<br />

<strong>in</strong> oblikovnih dopolnitev. Te so nastale kot priprava<br />

knjige <strong>za</strong> drugo izdajo, ki bo dopolnila obravnavano snov <strong>in</strong> jo<br />

bo razširila še na obsežnejšo obravnavo naključnih <strong>signalov</strong>. Zanjo<br />

upava, da bo doživela natis konec letošnjega leta.<br />

marec 2001<br />

Žarko Čučej <strong>in</strong> Peter Plan<strong>in</strong>šič<br />

i


ii<br />

Vseb<strong>in</strong>a knjige<br />

Snov je podana v dveh nivojih. V osnovnem, elementarnem nivoju<br />

so opisi kar se da preprosti <strong>in</strong> dopolnjeni s številnimi zgledi.<br />

Zahtevnejšim študentom so namenjena poglavja, ki imajo oznako<br />

“<strong>za</strong>htevna snov”. V njih so strožje matematične izpeljave <strong>in</strong> dokazi<br />

nekaterih teoremov. Seveda študentom priporočava, da se<br />

ozrejo še po drugih slovenskih <strong>in</strong> tujih virih. Najpomembnejše,<br />

po katerih sva se tudi zgledovala, sva naštela v seznamu literature.<br />

Uvod V uvodu se seznanimo s pojmom signal, s področji uporabe<br />

teorije <strong>signalov</strong> <strong>in</strong> s pojmom sistem.<br />

Signali Po vpeljavi funkcij kot matematičnega opisa <strong>signalov</strong><br />

sledi opis razdelitve <strong>signalov</strong> glede na njihove značilnosti. Primerom<br />

pogostih <strong>signalov</strong> sledi matematični opis. Nadaljuje se z<br />

elementarnimi operacijami nad signali, ki obsegajo transformacije<br />

signalnega območja <strong>in</strong> signalne osi.<br />

Prostor <strong>signalov</strong> je predstavljen z l<strong>in</strong>earnim prostorom. Potrebujemo<br />

ga <strong>za</strong> def<strong>in</strong>iranje norm <strong>signalov</strong>, ki so predstavljene<br />

kot mera značilnosti <strong>signalov</strong>. Fizikalni pomen norm je poka<strong>za</strong>n<br />

z izračunom moči, energije <strong>in</strong> srednje vrednosti signala. Pri srednjih<br />

vrednostih je kot <strong>za</strong>htevna snov dodana izpeljava drugega<br />

izreka o srednji vrednosti, ki se pogosto koristi pri izpeljavi <strong>in</strong> dokazovanju<br />

lastnosti amplitudno omejenih, končno krat nezveznih<br />

<strong>signalov</strong>.<br />

Razlaga skalarnega produkta predstavi signalni prostor kot<br />

vektorski prostor. Razložen je tako <strong>za</strong> časovno diskretne signale<br />

kot <strong>za</strong> zvezne signale. Podana je tudi pomembna Schwartzova<br />

neenakost.<br />

Korelacija je podana kot mera podobnosti <strong>signalov</strong>. Utemeljena<br />

je kot količ<strong>in</strong>a pretoka energije med signali. Razlaga poteka<br />

vzporedno <strong>za</strong> časovno zvezne <strong>in</strong> diskretne signale.<br />

Poglavje <strong>za</strong>ključuje opis posplošenega signala δ(t). Elementarnemu<br />

opisu sledi <strong>za</strong>htevnejši opis, v katerem so izpeljane <strong>in</strong><br />

doka<strong>za</strong>ne lastnosti Diracovega impul<strong>za</strong>. Snov <strong>za</strong>okroži opis uporabe<br />

Diracovega impul<strong>za</strong> pri vzorčenju <strong>signalov</strong>.<br />

Izražanje <strong>signalov</strong> z baznimi funkcijami V tem poglavju<br />

je vpeljan kvantitativen opis <strong>signalov</strong>. Izražanje signala z bazno<br />

funkcijo služi <strong>za</strong> vpeljavo srednjega kvadratičnega pogreška<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


iii<br />

kot mero podobnosti med približkom <strong>in</strong> orig<strong>in</strong>alom signala. Za<br />

predstavitev signala s komb<strong>in</strong>acijo baznih funkcij se vpeljejo ortogonalne<br />

funkcije. S Parsevalovo identiteto je poudarjen pomen<br />

kvantitativnega opisovanja <strong>signalov</strong> z ortogonalnimi funkcijami.<br />

Poglavje <strong>za</strong>ključuje predstavitev signalnega prostora kot vektorskega<br />

prostora, ki ga razpenjajo bazne funkcije. V njem so<br />

signali predstavljene kot točke, ki jih določajo konice vektorjev<br />

oziroma koord<strong>in</strong>ate - koeficienti, s katerimi množimo bazne funkcije<br />

pri kvantitativnem izražanju <strong>signalov</strong>.<br />

Poglavje <strong>za</strong>ključuje pregled različnih razredov ortogonalnih<br />

funkcij. Med njimi so širše opisane le Walsheve funkcije. Dodan<br />

je še Gram-Schmidtov postopek snovanja ortogonalnih funkcij.<br />

Sistemi To poglavje vsebuje le kratek pregled sistemov. Razlaga<br />

je omejena na vhodno-izhodni opis sistemov. Opisane so<br />

osnovne lastnosti <strong>in</strong> nač<strong>in</strong>i povezovanja sistemov. Opis l<strong>in</strong>earnih<br />

sistemov je namenjen razlagi konvolucije <strong>in</strong> lastnosti konvolucijskih<br />

sistemov.<br />

Zahvala<br />

Avtorja se <strong>za</strong>hvaljujeva prof. Rajku Svečku <strong>za</strong> opravljeno recenzijo.<br />

Njegove pripombe <strong>in</strong> nasveti so v veliki meri izboljšale<br />

knjigo. Uskladil je tudi uporabljeno izrazoslovje z izrazoslovjem<br />

v Teoriji sistemov <strong>in</strong> Teoriji regulacij. Posebna <strong>za</strong>hvala velja<br />

doc. dr. Gorazdu Lešnjaku. Njegova recenzija je s številnimi<br />

predlogi <strong>in</strong> popravki zelo pripomogla k izboljšanju izpeljav <strong>in</strong><br />

matematičnega jezika ter s tem k boljši razumljivosti teksta. Zahvaljujeva<br />

se prof. Mileni Milanovič, ki je pomagala izpiliti nelahek<br />

jezik. Posebna hvala gre mag. Jože Mohorku <strong>in</strong> dr. Bojanu<br />

Gergiču <strong>za</strong> natančno čitanje ”krtačnega odtisa”, <strong>za</strong> vse pripombe<br />

<strong>in</strong> <strong>za</strong> ponoven preračun vseh primerov.<br />

Maribor, marec 2000<br />

Žarko Čučej<br />

Peter Plan<strong>in</strong>šič


iv<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


Ka<strong>za</strong>lo<br />

Uvod . . . . . . . . . . . . . . . . . . . . . ii<br />

Signali . . . . . . . . . . . . . . . . . . . . . ii<br />

Izražanje <strong>signalov</strong> z baznimi funkcijami . . ii<br />

Sistemi . . . . . . . . . . . . . . . . . . . . iii<br />

1 Uvod 1<br />

1.1 Uporaba teorije <strong>signalov</strong> . . . . . . . . . . . . . . . 1<br />

1.1.1 Odelava <strong>signalov</strong> . . . . . . . . . . . . . . . 2<br />

1.1.2 Komunikacije . . . . . . . . . . . . . . . . . 2<br />

1.1.3 Regulacije . . . . . . . . . . . . . . . . . . . 4<br />

1.2 Sistemi . . . . . . . . . . . . . . . . . . . . . . . . . 5<br />

1.3 Razlikovanje med <strong>in</strong>formacijami, signali <strong>in</strong><br />

prenosnimi mediji . . . . . . . . . . . . . . . . . . . 5<br />

2 Signali 7<br />

2.1 Vrste <strong>signalov</strong> . . . . . . . . . . . . . . . . . . . . . 8<br />

2.1.1 Periodični - aperiodični signali . . . . . . . 8<br />

2.1.2 Naključni - nenaključni signali . . . . . . . 9<br />

2.1.3 Zvezni - nezvezni (diskretni) signali . . . . . 9<br />

Zvezni signali . . . . . . . . . . . . . . . . . 9<br />

Amplitudno diskretni signali . . . . . . . . 9<br />

Časovno diskretni signali . . . . . . . . . . 9<br />

Amplitudno <strong>in</strong> časovno diskretni signali . . 9<br />

Digitalni signali . . . . . . . . . . . . . . . . 10<br />

B<strong>in</strong>arni signali . . . . . . . . . . . . . . . . 10<br />

2.1.4 Realni - kompleksni signali . . . . . . . . . 10<br />

2.2 Primeri pomembnih <strong>signalov</strong> . . . . . . . . . . . . 11<br />

Diskretni enotski impulz . . . . . . . . . . . 11<br />

Enotska stopnica . . . . . . . . . . . . . . . 11<br />

v


vi<br />

KAZALO<br />

Klanec . . . . . . . . . . . . . . . . . . . . . 11<br />

Pravokotni pulz . . . . . . . . . . . . . . . . 12<br />

Trikotni pulz . . . . . . . . . . . . . . . . . 12<br />

2.2.1 Razlika med impulzi <strong>in</strong> pulzi . . . . . . . . 12<br />

2.3 Matematični opisi <strong>signalov</strong> . . . . . . . . . . . . . . 13<br />

2.4 Elementarne operacije nad signali . . . . . . . . . . 17<br />

2.4.1 Transformacija signalnega območja . . . . . 17<br />

2.4.2 Kvanti<strong>za</strong>cija . . . . . . . . . . . . . . . . . . 18<br />

2.4.3 Transformacija signalne osi . . . . . . . . . 20<br />

2.4.4 Premik po signalni osi . . . . . . . . . . . . 20<br />

2.5 Soda <strong>in</strong> liha simetričnost signala . . . . . . . . . . 22<br />

2.6 L<strong>in</strong>earni prostor . . . . . . . . . . . . . . . . . . . . 24<br />

2.7 Norme <strong>signalov</strong> . . . . . . . . . . . . . . . . . . . . 25<br />

2.7.1 Norme v prostoru L . . . . . . . . . . . . . 27<br />

2.7.2 Fizikalni pomen norm . . . . . . . . . . . . 27<br />

2.7.3 Moč signala, efektivna <strong>in</strong><br />

povprečna vrednost signala . . . . . . . . . 28<br />

Trenutna moč signala . . . . . . . . . . . . 29<br />

Trenutna medsebojna moč dveh <strong>signalov</strong> . . 29<br />

Povprečna moč . . . . . . . . . . . . . . . . 30<br />

Povprečna moč časovno neomejenih <strong>signalov</strong> 31<br />

Efektivna vrednost signala . . . . . . . . . . 31<br />

Povprečna vrednost signala . . . . . . . . . 32<br />

2.7.4 Drugi izrek o povprečni vrednosti . . . . . . 33<br />

2.8 Skalarni produkt . . . . . . . . . . . . . . . . . . . 36<br />

2.8.1 Skalarni produkt v prostoru l . . . . . . . . 37<br />

2.8.2 Skalarni produkt v prostoru L . . . . . . . 38<br />

2.8.3 Cauchy-Schwartzova neenakost . . . . . . . 39<br />

2.8.4 Signalni prostor . . . . . . . . . . . . . . . . 40<br />

2.9 Podobnost <strong>signalov</strong> . . . . . . . . . . . . . . . . . . 40<br />

2.9.1 Koncept . . . . . . . . . . . . . . . . . . . . 41<br />

2.9.2 Korelacija pri aperiodičnih signalih . . . . . 42<br />

2.9.3 Korelacija pri periodičnih signalih . . . . . 44<br />

2.9.4 Korelacija pri kompleksnih aperiodičnih<br />

signalih . . . . . . . . . . . . . . . . . . . . 44<br />

2.9.5 Izračun energije <strong>in</strong> moči s korelacijo . . . . 45<br />

2.9.6 Robni efekt . . . . . . . . . . . . . . . . . . 45<br />

2.9.7 Korelacijski koeficient . . . . . . . . . . . . 48<br />

2.10 Posplošeni signali . . . . . . . . . . . . . . . . . . . 49<br />

2.10.1 Osnove teorije s<strong>in</strong>gularnih funkcij . . . . . . 50<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


KAZALO<br />

vii<br />

Vrednost signala v trenutku t = 0 . . . . . . 50<br />

Vrednost signala v trenutku t = a . . . . . . 51<br />

2.10.2 L<strong>in</strong>earne komb<strong>in</strong>acije delta funkcij . . . . . 52<br />

Seštevanje <strong>in</strong> množenje s skalarjem . . . . . 52<br />

Vzorčenje <strong>signalov</strong> . . . . . . . . . . . . . . 53<br />

Odvod Diracovega impul<strong>za</strong> . . . . . . . . . 53<br />

2.10.3 Pove<strong>za</strong>va med enotsko stopnico<br />

<strong>in</strong> Diracovim impulzom . . . . . . . . . . . 54<br />

3 Izražanje <strong>signalov</strong> z baznimi funkcijami 57<br />

3.1 Aproksimacija <strong>signalov</strong> z bazno funkcijo . . . . . . 58<br />

3.2 Aproksimacija <strong>signalov</strong> z l<strong>in</strong>earno komb<strong>in</strong>acijo<br />

baznih funkcij . . . . . . . . . . . . . . . . . . 61<br />

3.3 Izražanje <strong>signalov</strong> z ortogonalnimi funkcijami . . . 62<br />

3.3.1 Izražanje <strong>signalov</strong> z<br />

ortonormiranimi <strong>za</strong>poredji . . . . . . . . . . 63<br />

Določitev koeficientov c n s pomočjo<br />

ortogonalnosti . . . . . . . . . . . 63<br />

Določitev koeficientov c n po metodi<br />

najmanjšega kvadratičnega<br />

pogreška . . . . . . . . . . . . . . 64<br />

3.3.2 Parsevalova identiteta . . . . . . . . . . . . 65<br />

3.4 Ortogonalne funkcije . . . . . . . . . . . . . . . . . 66<br />

3.4.1 Pregled nekaterih ortogonalnih funkcij . . . 67<br />

3.4.2 Walsheve funkcije . . . . . . . . . . . . . . 68<br />

3.4.3 Gram - Schmidtov postopek . . . . . . . . . 71<br />

3.5 Signalni prostor (drugič) . . . . . . . . . . . . . . . 74<br />

3.5.1 Vektorski prostor . . . . . . . . . . . . . . . 74<br />

3.5.2 Signalni prostor je vektorski prostor . . . . 76<br />

3.5.3 Predstavitev <strong>signalov</strong> kot vektorjev . . . . . 76<br />

4 Sistemi 77<br />

4.1 Vhodno-izhodni sistemi . . . . . . . . . . . . . . . 77<br />

4.1.1 Primeri sistemov . . . . . . . . . . . . . . . 80<br />

4.1.2 Časovno zvezni <strong>in</strong> časovno diskretni sistemi 80<br />

4.1.3 Povezovanje sistemov . . . . . . . . . . . . . 81<br />

4.2 Sistemi s povratno <strong>za</strong>nko . . . . . . . . . . . . . . 82<br />

4.3 Lastnosti sistemov . . . . . . . . . . . . . . . . . . 84<br />

Sistem s pomnjenjem . . . . . . . . . . . . . 84<br />

Invertibilnost . . . . . . . . . . . . . . . . . 85


viii<br />

KAZALO<br />

Inverzni sistem . . . . . . . . . . . . . . . . 85<br />

Vzročnost . . . . . . . . . . . . . . . . . . . 86<br />

Stabilnost . . . . . . . . . . . . . . . . . . . 87<br />

Časovna neodvisnost . . . . . . . . . . . . . 87<br />

4.4 L<strong>in</strong>earni sistemi . . . . . . . . . . . . . . . . . . . . 88<br />

4.4.1 Odziv l<strong>in</strong>earnega časovno diskretnega<br />

sistema . . . . . . . . . . . . . . . . . . . . 88<br />

4.4.2 Odziv l<strong>in</strong>earnega časovno<br />

neodvisnega zveznega sistema . . . . . . . . 91<br />

4.4.3 Konvolucija ali pregib . . . . . . . . . . . . 93<br />

Konvolucijski sistemi . . . . . . . . . . . . . 94<br />

Pregib . . . . . . . . . . . . . . . . . . . . . 94<br />

4.4.4 Impulzni odziv . . . . . . . . . . . . . . . . 98<br />

4.5 Lastnosti konvolucije . . . . . . . . . . . . . . . . . 100<br />

Komutativnost . . . . . . . . . . . . . . . . 100<br />

Asociativnost . . . . . . . . . . . . . . . . . 100<br />

Distributivnost . . . . . . . . . . . . . . . . 100<br />

Komutativnost množenja s skalarjem . . . . 100<br />

Premik funkcije . . . . . . . . . . . . . . . . 100<br />

Lastnosti odvajanja . . . . . . . . . . . . . 100<br />

4.5.1 Konvolucija signala z Diracovim impulzom 101<br />

4.5.2 Konvolucija signala z odvodom Diracovega<br />

impul<strong>za</strong> . . . . . . . . . . . . . . . . 102<br />

4.6 Obstoj konvolucije . . . . . . . . . . . . . . . . . . 103<br />

4.7 Zaključek . . . . . . . . . . . . . . . . . . . . . . . 103<br />

Seznam oznak 105<br />

Literatura 107<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


Poglavje 1<br />

Uvod<br />

S signali pojmujemo abstraktni opis poteka nekega razvijajočega<br />

se fizikalnega procesa, kot ga lahko opazujemo neposredno s čutili<br />

ali posredno z izvajanjem meritev fizikalne količ<strong>in</strong>e, ki nas <strong>za</strong>nima.<br />

Torej nas signali ”<strong>in</strong>formirajo”, kaj se z nekim procesom<br />

dogaja oziroma lahko z njim o tem sporočamo. Povedano<br />

drugače, signali so nosilci <strong>in</strong>formacij. In kaj je <strong>in</strong>formacija? Danes<br />

govorimo, da je <strong>in</strong>formacija poleg mase <strong>in</strong> energije tretja najpomembnejša<br />

osnovna velič<strong>in</strong>a. Njen pomen z razvojem računalništva,<br />

ki omogoča avtomatsko <strong>obdelavo</strong> <strong>in</strong>formacij, nenehno<br />

raste. Današnjo dobo nekateri opisujejo tudi kot <strong>in</strong>formacijsko<br />

dobo. Z <strong>in</strong>formacijami na splošno se ukvarja teorija <strong>in</strong>formacij,<br />

v tej knjigi pa je predstavljen uvod v teorijo <strong>signalov</strong>. Teorijo<br />

<strong>signalov</strong> uporabimo, ko želimo iz <strong>signalov</strong> izluščiti <strong>in</strong>formacijo, ki<br />

jo nosijo, ali pa <strong>in</strong>formacijo vpisati, ko jo želimo z njimi prenesti.<br />

V tej knjigi so opisani uvodni pojmi iz teorije <strong>signalov</strong>. Ti<br />

obsegajo razdelitev <strong>signalov</strong> v značilne skup<strong>in</strong>e, elementarne signale,<br />

nač<strong>in</strong>e njihovega opisovanja, mere značilnih lastnosti <strong>signalov</strong>,<br />

mere podobnosti, ortogonalne funkcije <strong>in</strong> <strong>za</strong>pis <strong>signalov</strong><br />

z ortogonalnimi funkcijami. Zapis <strong>za</strong>ključuje opis sistemov. Ti<br />

so pomembni <strong>za</strong>to, ker z njimi signale obdelujemo.<br />

1.1 Uporaba teorije <strong>signalov</strong><br />

Teorija <strong>signalov</strong> <strong>in</strong> teorija sistemov se uporabljata v mnogih področjih<br />

znanosti <strong>in</strong> tehnike. V tem razdelku bomo opisali tri<br />

1


2 1. Uvod<br />

primere uporabe teorije <strong>signalov</strong>. Primeri so s področja obdelave<br />

<strong>signalov</strong>, komunikacij <strong>in</strong> regulacij.<br />

1.1.1 Odelava <strong>signalov</strong><br />

Značilni problem, ki ga rešujemo z <strong>obdelavo</strong> <strong>signalov</strong>, je odstranjevanje<br />

neželenega šuma v signalu. Tega na primer želimo odstraniti<br />

v sprejemniku telefonskega signala, vhodnem ojačevalniku<br />

magnetofona <strong>in</strong> podobno. Model, s katerim običajno opišemo<br />

problem, ima naslednjo matematično obliko:<br />

v = s + n ,<br />

kjer so v sprejeti signal, s poslani signali <strong>in</strong> n šum, ki se je dodal<br />

- prištel - k signalu.<br />

Kako razdvojiti signal od šuma? Šum skušamo odstraniti<br />

ali vsaj zelo zmanjšati s sitom (slika 1.1-1), s katerim presejemo<br />

sprejeti signal. Izhod sita je signal y, <strong>za</strong> katerega želimo, da je<br />

enak oziroma bolj podoben poslanemu signalu, kot je bil sprejeti<br />

signal.<br />

Slika 1.1-1<br />

Problem načrtovanja<br />

vhodnega sita <strong>za</strong> izločevanje<br />

šuma.<br />

signal<br />

s<br />

šum<br />

n<br />

sprejeti<br />

signal<br />

v = s+n<br />

sito<br />

izsejani<br />

signal<br />

y<br />

Pri načrtovanju sita predpostavimo, da je energija šuma koncentrirana<br />

v višjem frekvenčnem področju kot je več<strong>in</strong>a energije<br />

signala, torej da je sprem<strong>in</strong>janje vrednosti šuma mnogo hitrejše<br />

kot je sprem<strong>in</strong>janje vrednosti signala. Če iz sprejetega signala<br />

odrežemo komponente, ki so nad določeno frekvenco, signal zgladimo.<br />

S tem sicer odstranimo določeno količ<strong>in</strong>o šuma, hkrati pa<br />

spremenimo tudi signal. Močnejši je uč<strong>in</strong>ek glajenja, več šuma<br />

odstranimo, hkrati pa bolj spremenimo signal. Zato v takih primerih<br />

vedno iščemo kompromis. Pri tem si pomagamo s teorijo<br />

<strong>signalov</strong>, s katero proučujemo uč<strong>in</strong>ek sita na potek signala.<br />

1.1.2 Komunikacije<br />

Pri komunikacijah imamo opravka s prenosom sporočil, ki so vpisane<br />

v obliki na primer radijskega ali televizijskega signala, ali pa<br />

<strong>za</strong>pisana v toku podatkov, ki oblikujejo izbrani signal - nosilec<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


1.1 Uporaba teorije <strong>signalov</strong> 3<br />

<strong>in</strong>formacije. Pri tem mnogokrat <strong>za</strong>htevamo prenos <strong>signalov</strong> na<br />

večje razdalje.<br />

Ponavadi so orig<strong>in</strong>alni signali neprimerni <strong>za</strong> prenos na večje<br />

razdalje, <strong>za</strong>to jih ”nasadimo”na poseben signal, ki ga imenujemo<br />

nosilec. Nosilec izberemo tako, da je primeren tudi <strong>za</strong> prenose<br />

na zelo velike razdalje. Proces ”nasajanja”<strong>in</strong>formacije na nosilec<br />

imenujemo modulacija. Rezultat modulacije je modulirani signal,<br />

ki se potem širi po prenosnem mediju. Ta je lahko kov<strong>in</strong>ski<br />

vodnik, svetlo(bo)vodnik, prostor ali drugi medij. Na sprejemni<br />

strani poslano sporočilo izluščimo iz moduliranega signala s<br />

postopkom, ki je <strong>in</strong>verzen modulaciji, <strong>za</strong>to ga imenujemo demodulacija.<br />

Modulacija <strong>in</strong> demodulacija se izvedeta v napravah, ki jih<br />

imenujemo modulatorji oziroma demodulatorji (slika 1.1-2). Te<br />

naprave svoji funkciji izvajata tako, da so neželene posledice prenosa<br />

- to je šum, ki se prišteje signalu - odstranjene v največji<br />

možni meri. Tudi pri tem si pomagamo s teorijo <strong>signalov</strong>. Z njeno<br />

pomočjo lahko na primer načrtamo optimalni sprejemnik, ki maksimira<br />

razmerje med amplitudo signala <strong>in</strong> povprečno energijo<br />

šuma.<br />

modulacijski<br />

signal s<br />

modulator<br />

nosilec<br />

modulirani<br />

signal x<br />

prenosni<br />

medij<br />

demodulator<br />

sprejeti<br />

signal y<br />

Slika 1.1-2<br />

Prenos <strong>in</strong>formacije z<br />

modulacijo nosilca.<br />

Postopek modulacije si oglejmo na primeru amplitudne modulacije,<br />

ko <strong>in</strong>formacijo določa digitalni signal m, ki ga imenujemo<br />

tudi modulacijski signal, nosilec pa je s<strong>in</strong>usni val z visoko<br />

frekvenco (slika 1.1-3). Vidimo, da je rezultat modulacijskega<br />

postopka signal, katerega amplituda ni več konstantna, ampak<br />

je odvisna od modulacijskega signala. Lastnosti takega signala<br />

določimo s pomočjo teorije <strong>signalov</strong>.<br />

Modulacijski <strong>in</strong> demodulacijski postopki so zelo pomemben<br />

del komunikacijske teorije. Poznamo jih celo vrsto, od tako imenovanih<br />

amplitudnih, kotnih pa vse do digitalnih modulacij, katere<br />

primer smo poka<strong>za</strong>li na sliki 1.1-3. Posebna zvrst modulacij<br />

so modulacije z razprševanjem spektra.


4 1. Uvod<br />

modulacijski signal (<strong>in</strong>formacija)<br />

Slika 1.1-3<br />

Primer modulacijskega <strong>in</strong><br />

moduliranega signala. Nosilec<br />

je s<strong>in</strong>usni val.<br />

1<br />

0<br />

-1<br />

modulirani signal<br />

3<br />

0<br />

-3<br />

t<br />

t<br />

1.1.3 Regulacije<br />

Pri regulacijah želimo načrtati naprave, ki samodejno vplivajo<br />

na d<strong>in</strong>amične lastnosti danega sistema. To dosežemo s sklenitvijo<br />

povratne <strong>za</strong>nke, to je s krmiljenjem sistema z razliko med<br />

želenim <strong>in</strong> izhodnim signalom. S tem skušamo odpraviti odstopanje<br />

od želenega obnašanja oziroma lastnosti (slika 1.1-4). Izbira<br />

primernega regulatorja je osrednji problem teorije regulacij. S<br />

povratno <strong>za</strong>nko se bomo bežno srečali pri obravnavi sistemov.<br />

Slika 1.1-4<br />

Uravnavanje d<strong>in</strong>amičnih<br />

lastnosti sistema z regulacijo.<br />

referenèni<br />

signal v<br />

regulacijsko<br />

odstopanje<br />

v - y<br />

regulator<br />

regulirni<br />

signal x<br />

regulirani<br />

objekt<br />

signal<br />

regulirane<br />

koliè<strong>in</strong>e<br />

y<br />

Komunikacije <strong>in</strong> regulacije danes obravnavamo kot samostojni<br />

znanstveni discipl<strong>in</strong>i. S prvimi se predvsem ukvarjajo komunikacijski<br />

<strong>in</strong>ženirji, z drugimi pa <strong>in</strong>ženirji avtomatike. Obe discipl<strong>in</strong>i<br />

pa imata isto izhodišče: teorijo <strong>signalov</strong>. Še več, težko je najti<br />

komunikacijsko napravo ali postopek, ki ne bi vseboval elementov<br />

regulacij. Pri regulacijah pa morajo regulatorji <strong>in</strong> regulirani<br />

objekti medsebojno izmenjevati <strong>in</strong>formacije o svojem stanju, torej<br />

morajo medsebojno komunicirati. Danes obstajajo mnogi sistemi<br />

vodenja, kjer so regulacije porazdeljene, njihovi gradniki pa<br />

so medsebojno pove<strong>za</strong>ni z ustreznimi komunikacijskimi sistemi.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


1.2 Sistemi 5<br />

1.2 Sistemi<br />

S term<strong>in</strong>om sistemi opredeljujemo vsak del okolja, ki ustvarja<br />

signale oziroma se nanje odziva. Primer sistema je električno<br />

vezje. V njem so lahko signali električne napetosti <strong>in</strong> tokovi,<br />

magnetni pretoki ali električni naboji. Poleg električnih sistemov<br />

so še na primer mehanski, kjer so signali na primer hitrost<br />

gibanja objekta, hitrost vrtenja, pospešek itd. Iz teh primerov<br />

lahko <strong>za</strong>ključimo, da so signali tesno pove<strong>za</strong>ni s sistemi. Določajo<br />

jih <strong>in</strong>terni procesi sistemov. Na primer napetosti <strong>in</strong> tokovi so z<br />

električnim vezjem pove<strong>za</strong>ni z elektromagnetnimi <strong>za</strong>koni; lega,<br />

hitrost, pospešek so z mehanskim sistemom pove<strong>za</strong>ni z Newtonovimi<br />

<strong>za</strong>koni itd. Pove<strong>za</strong>ve med sistemi <strong>in</strong> signali, ki so z njimi<br />

pove<strong>za</strong>ni, imenujemo <strong>za</strong>konitosti sistema. Zakonitosti sistemov<br />

tudi določajo, kako se sistem odzove na določen signal. Odziv<br />

sistema je spet signal, <strong>za</strong>to lahko s sistemi signale preoblikujemo<br />

oziroma jih obdelujemo. V tej zvezi ni teorija obdelave <strong>signalov</strong><br />

nič drugega kot skupek <strong>za</strong>konitosti, ki povezujejo sisteme <strong>in</strong><br />

signale.<br />

1.3 Razlikovanje med <strong>in</strong>formacijami,<br />

signali <strong>in</strong> prenosnimi mediji<br />

V prvem odstavku smo signal opredelili kot abstraktni opis poteka<br />

nekega fizikalnega procesa. Dodali smo, da nas signali <strong>in</strong>formirajo,<br />

kako poteka proces ter jih <strong>za</strong>to obravnavamo kot nosilce<br />

<strong>in</strong>formacij. Še v isti sapi pa smo dodali, da nas <strong>in</strong>formacije ne<br />

<strong>za</strong>nimajo, ter da se bomo ukvarjali le s signali.<br />

Zato se seveda takoj postavi vprašanje, kako razlikovati med<br />

<strong>in</strong>formacijami <strong>in</strong> signali. Dejansko imamo opraviti z več pojmi,<br />

katere bomo skušali na kratko tako pojasniti, da jih bomo razumeli<br />

<strong>in</strong> zlahka medsebojno razlikovali.<br />

Informacija<br />

V vsakodnevnem življenju se pogosto srečujemo z<br />

”<strong>in</strong>formacijami”. Ta pojem pogosto uporabljamo<br />

zelo ohlapno, <strong>za</strong>to poglejmo, kaj o njem pravi Slovenski<br />

slovar knjižnega jezika:<br />

”Informacija je, kar se o določeni stvari pove.”<br />

Torej so <strong>in</strong>formacije nova dejstva, pojasnila, razna


6 1. Uvod<br />

Signal<br />

Nosilec<br />

signala<br />

Prenosni<br />

medij<br />

sporočila <strong>in</strong> drugo. Takoj pa poudarimo, da <strong>in</strong>formacije<br />

niso znanje ali razumevanje. Za to so potrebni<br />

miselni procesi, oziroma (naravna ali umetna) <strong>in</strong>teligenca.<br />

Na <strong>in</strong>formacije lahko gledamo z različnih vidikov, od<br />

splošnega filozofskega stališča pa do povsem tehniškega.<br />

S tehniškega stališča so <strong>in</strong>formacije množica<br />

elementov, ki predstavljajo podatke, nova dejstva,<br />

pojasnila, merilne rezultate <strong>in</strong> podobno. To množico<br />

imenujemo <strong>in</strong>formacijska množica ali tudi <strong>in</strong>formacijski<br />

prostor.<br />

Akademik prof. Ludvig Gyergyek je v svoji knjigi Signali<br />

<strong>in</strong> statistične metode [16] <strong>za</strong>pisal, da so signali<br />

kočije, ki vozijo <strong>in</strong>formacije. Z drugimi besedami,<br />

<strong>in</strong>formacije določajo potek signala. Ker <strong>za</strong> opis poteka<br />

uporabljamo funkcije, je signal kar s<strong>in</strong>onim <strong>za</strong><br />

(časovne) funkcije.<br />

Nosilec signala je fizikalna količ<strong>in</strong>a, katere potek določa<br />

signal. Zato mnogokrat govorimo o električnem,<br />

akustičnem, svetlobnem, . . . signalu.<br />

Prenosni medij je snov, po kateri se širi nosilec signala.<br />

Snov ima lahko posebno obliko, na primer<br />

kov<strong>in</strong>ski vodnik, valovod, svetlovodnik.<br />

Elektromagnetno valovanje, tako radijsko kot tudi<br />

svetloba, <strong>za</strong> svoje širjenje ne potrebujeta posebnega<br />

medija. Širi se lahko tudi po praznem prostoru. Ker<br />

sprva to znanstvenikom ni bilo razumljivo, so predpostavljali,<br />

da je ves prostor <strong>za</strong>polnjen s posebnim<br />

fluidom, po katerem se širi svetloba (oziroma elektromagnetno<br />

valovanje na splošno). Zato še danes<br />

pogosto rečemo, da radijski program oddajamo v<br />

”eter”.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


Poglavje 2<br />

Signali<br />

Kako predstaviti oziroma opisati signal? V ta namen uporabljamo<br />

matematične funkcije. Ko je potek <strong>signalov</strong> odvisen od<br />

časa - več<strong>in</strong>a <strong>signalov</strong>, ki jih bomo obravnavali, bo takšnih - so<br />

funkcije časovne funkcije. Zato bomo časovne signale več<strong>in</strong>oma<br />

označevali z<br />

x(t) .<br />

Pri tem je čas t argument, ki <strong>za</strong>jema svoje vrednosti iz def<strong>in</strong>icijskega<br />

območja signala. To območje imenujemo tudi signalna os<br />

<strong>in</strong> jo bomo v primeru časovnih <strong>signalov</strong> označevali s:<br />

T : signalna os signala x(t) .<br />

Vrednosti x(t) določajo območje signala A, ki ga imenujemo tudi<br />

amplitudni razmah signala. S tem dogovorom lahko signale predstavimo<br />

z grafi. Ti so pri pona<strong>za</strong>rjanju funkcij običajni. Na primer<br />

signal, ki ga opišemo s funkcijo x(t) = 1 + s<strong>in</strong> t, t ∈ [0, π],<br />

kaže slika 2.0-1.<br />

signalno obmoèje<br />

2<br />

1<br />

0<br />

0 p/2 p<br />

signalna os<br />

Slika 2.0-1<br />

Graf signala<br />

x(t) = 1 + s<strong>in</strong> t , t ∈ [0, π].<br />

7


8 2. Signali<br />

Sklepamo lahko, da bomo pri obdelavi <strong>signalov</strong> uporabljali<br />

teorijo funkcij oziroma da bomo pri tem izkoriščali tiste lastnosti<br />

funkcij, ki bodo vodile do želenih rezultatov obdelave. Zato takoj<br />

def<strong>in</strong>irajmo signal v matematičnem jeziku:<br />

DEFINICIJA 1<br />

Če je A množica realnih števil <strong>in</strong> T podmnožica realnih števil,<br />

potem vsaka funkcija, ki preslikuje x : T → A, določa signal s<br />

signalno osjo T <strong>in</strong> s signalnim območjem A.<br />

Množico vseh <strong>signalov</strong> s signalno osjo T <strong>in</strong> s signalnim območjem<br />

A določa množica vseh funkcij, ki preslikujejo T v A. To množico<br />

označimo z {x : T → A} ali s krajšim <strong>za</strong>pisom s potenčno<br />

množico:<br />

A T .<br />

Iz matematike nam je znano, da obstaja neskončno funkcij. Tako<br />

je <strong>za</strong> opis <strong>signalov</strong> na voljo obsežno, raznoliko, predvsem pa<br />

uč<strong>in</strong>kovito orodje, ki ga bomo s pridom uporabili tako pri analizi<br />

<strong>in</strong> s<strong>in</strong>tezi <strong>signalov</strong>, kakor tudi pri njihovem prepoznavanju<br />

(detekciji) oziroma medsebojnem razlikovanju.<br />

2.1 Vrste <strong>signalov</strong><br />

Preden se bomo podrobneje spoznali z opisom <strong>signalov</strong>, si oglejmo<br />

nekatere razdelitve <strong>signalov</strong>. Razvrstitev ali klasifikacijo <strong>signalov</strong><br />

lahko izvedemo na več nač<strong>in</strong>ov. Oglejmo si le najbolj pogoste.<br />

2.1.1 Periodični - aperiodični signali<br />

Periodični signal je tisti, pri katerem se v enakih časovnih <strong>in</strong>tervalih<br />

ponavlja <strong>za</strong>poredje njegovih vrednosti. Interval ponavljanja<br />

imenujemo perioda signala. Signal x(t) je periodičen, če obstaja<br />

neko število T , ki izpolni enačbo:<br />

x(t) = x(t + T ) = x(t + 2T ) = x(t + 3T ) = · · · (2.1-1)<br />

pri neomejeno veliki signalni osi T. Najmanjše pozitivno število<br />

T , ki ustre<strong>za</strong> enačbi (2.1-1), je osnovna perioda <strong>in</strong> določa trajanje<br />

ene periode ali cikla signala. Frekvenca periodičnega signala je<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.1 Vrste <strong>signalov</strong> 9<br />

število ponovitev signala v eni sekundi, torej:<br />

f = 1 T .<br />

Aperiodični signal je tisti, pri katerem ni možno najti nobene<br />

vrednosti T , ki bi izpolnila (2.1-1). To je zelo pomembna skup<strong>in</strong>a<br />

<strong>signalov</strong>, saj so v naravi vsi signali aperiodični. Signali, ki jih<br />

rojevajo fizikalni procesi, imajo svoj <strong>za</strong>četek (<strong>in</strong> verjetno tudi svoj<br />

konec), <strong>za</strong>to ne morejo, kljub temu, da v določenem časovnem<br />

segmentu izpolnijo pogoj (2.1-1), povsem <strong>za</strong>dostiti temu pogoju.<br />

2.1.2 Naključni - nenaključni signali<br />

Naključni ali stohastični signali so preprosto povedano signali,<br />

<strong>za</strong> katere obstaja določena stopnja ne<strong>za</strong>nesljivosti, da bo signal<br />

nastopil ali da bo <strong>za</strong>vzel določeno vrednost. Če ne moremo napovedati<br />

vrednosti signala v naslednjem trenutku, tudi če poznamo<br />

njegovo preteklost do tega trenutka, tedaj signala ni mogoče opisati<br />

v eksplicitni matematični obliki.<br />

Nenaključni ali determ<strong>in</strong>istični signali so tisti, pri katerih<br />

ni nedoločenosti o njihovem nastanku <strong>in</strong> poteku. Skoraj v vseh<br />

primerih lahko determ<strong>in</strong>istične signale predstavimo v eksplicitni<br />

matematični obliki.<br />

2.1.3 Zvezni - nezvezni (diskretni) signali<br />

Druga pomembna delitev <strong>signalov</strong> izhaja iz amplitudne oblike <strong>in</strong><br />

časovne ločljivosti <strong>signalov</strong> (tabela 2.1). Glede na to delitev so<br />

se pri obdelavi <strong>signalov</strong> uveljavile posebne oznake <strong>in</strong> poimenovanja,<br />

ki jih podajamo v nadaljnjih podrazdelkih.<br />

Zvezni signali Ti signali so po času <strong>in</strong> amplitudi zvezne funkcije.<br />

Pri analognem signalu smejo amplitude znotraj vnaprej<br />

predpisanega signalnega območja <strong>za</strong>vzeti poljubno vrednost.<br />

Amplitudno diskretni signali To so signali, katerih amplitude<br />

lahko imajo le diskretne vrednosti, vrednost amplitude pa<br />

je znana v vsakem trenutku obstoja signala.<br />

Časovno diskretni signali Pri teh signalih so vrednosti amplitud<br />

znane le v diskretnih časovnih trenutkih.<br />

Amplitudno <strong>in</strong> časovno diskretni signali Te signale pogosto<br />

imenujemo kar diskretni signali. Njihove amplitude lahko


10 2. Signali<br />

Tabela 2.1<br />

Delitev <strong>signalov</strong> po amplitudni <strong>in</strong> časovni ločljivosti.<br />

čas<br />

amplituda zvezen diskreten<br />

zvezna<br />

x( t)<br />

x( nT)<br />

x( t)<br />

t<br />

T<br />

t<br />

diskretna<br />

x t q ( )<br />

x( t)<br />

x n q ( )<br />

x q (t)<br />

t<br />

t<br />

<strong>za</strong>v<strong>za</strong>mejo le diskretne vrednosti, določene pa so v diskretnih<br />

časovnih trenutkih.<br />

Digitalni signali To so diskretni signali, ki imajo svojo vrednost<br />

podano s številko (s končnim številom cifer).<br />

B<strong>in</strong>arni signali To je posebna oblika digitalnih <strong>signalov</strong>, kjer<br />

amplituda lahko ima le dve diskretni vrednosti.<br />

2.1.4 Realni - kompleksni signali<br />

Glede na to, ali signalno območje A <strong>za</strong>jema vrednosti iz množice<br />

realnih ali množice kompleksnih števil, ločimo realne <strong>in</strong> kompleksne<br />

signale:<br />

realni signali: A = R<br />

kompleksni signali: A = C<br />

OPOMBA 2-1 V tekstu bomo pri označevanju številskih množic uporabljali<br />

naslednje uveljavljene oznake:<br />

N :<br />

naravna števila<br />

Z : cela števila → Z + : nenegativna<br />

Z − : nepozitivna<br />

R :<br />

Z :<br />

realna števila<br />

cela števila<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.2 Primeri pomembnih <strong>signalov</strong> 11<br />

2.2 Primeri pomembnih <strong>signalov</strong><br />

Med vsemi signali je skup<strong>in</strong>a <strong>signalov</strong>, ki imajo v obdelavi <strong>signalov</strong><br />

<strong>in</strong> pri obravnavi sistemov <strong>za</strong> <strong>obdelavo</strong> <strong>signalov</strong> poseben,<br />

praktični pomen. Za nekatere med njimi so se uveljavila tudi<br />

posebna imena <strong>in</strong> oznake.<br />

Diskretni enotski impulz Diskretni enotski impulz je določen<br />

z:<br />

{<br />

1 pri n = 0<br />

∆(n) =<br />

n ∈ T , (2.2-1)<br />

0 drugje<br />

kjer je lahko T poljubna diskretna časovna os, ki svoje vrednosti<br />

<strong>za</strong>jema iz množic N, Z + ali Z − . Diskretni enotski impulz imenujemo<br />

tudi Kroneckerjev impulz. Njegov graf kaže slika 2.2-1.<br />

D<br />

1<br />

0<br />

-1 0 1 2 3 t<br />

Slika 2.2-1<br />

Diskretni enotski impulz.<br />

Enotska stopnica Enotska stopnica je def<strong>in</strong>irana z:<br />

{<br />

1 t 0<br />

u(t) =<br />

t ∈ T = R , (2.2-2)<br />

0 t < 0<br />

Njen graf kaže slika 2.2-2.<br />

1<br />

u( t) Slika 2.2-2<br />

Enotska stopnica<br />

0<br />

-1 0 1 2 3 t<br />

Klanec Klanec je def<strong>in</strong>iran z:<br />

{<br />

t t 0<br />

ramp(t) =<br />

0 t < 0<br />

t ∈ T = R , (2.2-3)<br />

Njen graf kaže slika 2.2-3.


12 2. Signali<br />

Slika 2.2-3<br />

Klanec<br />

ramp<br />

1<br />

0<br />

-1 0 1 2 3 t<br />

Pravokotni pulz Pravokotni pulz je def<strong>in</strong>irana z:<br />

{<br />

1 −1/2 t < 1/2<br />

rect(t) =<br />

t ∈ T = R , (2.2-4)<br />

0 drugje<br />

Njegov graf kaže slika 2.2-4.<br />

Slika 2.2-4<br />

Pravokotni pulz<br />

rect<br />

1<br />

0<br />

-1/2 0 1/2 t<br />

Trikotni pulz Trikotni pulz je def<strong>in</strong>irana z:<br />

{<br />

1 − |t| pri |t| < 1<br />

trian(t) =<br />

t ∈ T = R , (2.2-5)<br />

0 drugje<br />

Njegov graf pa kaže slika 2.2-4.<br />

Slika 2.2-5<br />

Trikotni pulz<br />

trian<br />

1<br />

0<br />

-1 0 1 t<br />

2.2.1 Razlika med impulzi <strong>in</strong> pulzi<br />

V predstvitvi pomembnih <strong>signalov</strong> smo uporabili dva term<strong>in</strong>a:<br />

impulz <strong>in</strong> pulz. Kakšna razlika pa je med njima? Iz prikazov na<br />

slikah 2.2-1 - 2.2-5 lahko uvidimo, da je med njimi bistvena razlika<br />

v njihovem trajanju. Med tem ko je diskretni impulz trenutnega<br />

trajanja, torej ničte dolž<strong>in</strong>e, so pulzi merljive dolž<strong>in</strong>e.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.3 Matematični opisi <strong>signalov</strong> 13<br />

V naravi lahko ustvarjamo le pulze, impul<strong>za</strong> sta le izjemno<br />

pomembna matematična pripomočka. V teoriji <strong>signalov</strong> uporabljamo<br />

dva impul<strong>za</strong>: diskretni ali Kroneckerjev impulz ter Diracov<br />

ali δ impulz. Slednjega bomo spoznali kasneje pri obravnavi<br />

posplošenih funkcij. Omenimo le, da je Kroneckerjev impulz diskretna<br />

prireditev Diracovega impul<strong>za</strong>.<br />

2.3 Matematični opisi <strong>signalov</strong><br />

Pri prikazu pogosto uporabljanih <strong>signalov</strong> smo podali tudi njihove<br />

matematične opise. V tem razdelku bomo te opise razširili<br />

<strong>in</strong> sistematizirali. Poka<strong>za</strong>li bomo tehnike opisov <strong>signalov</strong>, ki jih<br />

bomo kasneje, pri analizi <strong>signalov</strong>, s pridom uporabili.<br />

ZGLED 2.3-1 Risanje <strong>signalov</strong> z znanimi grafi<br />

Narišimo graf signala:<br />

x(t) = 3u(t) + t · u(t) − (t − 1) · u(t − 1) − 5u(t − 2) . (2.3-1)<br />

Grafe členov (2.3-1) kaže slika 2.3-1a, graf signala pa slika 2.3-1b.<br />

a<br />

x( t)<br />

4<br />

2<br />

3u( t)<br />

tu( t)<br />

0<br />

-( t -1) ◊u( t -1)<br />

-2<br />

-4<br />

-5u( t - 2)<br />

b<br />

x( t)<br />

4<br />

2<br />

0<br />

1<br />

2<br />

3<br />

4<br />

t<br />

0<br />

-2<br />

0<br />

1<br />

2<br />

Slika 2.3-1<br />

Primer konstrukcije grafa signala, a: sestav<strong>in</strong>e signala, b: signal<br />

3<br />

4<br />

t


14 2. Signali<br />

Posamezne člene x(t) v (2.3-1) določajo funkcije:<br />

{<br />

3 t 0<br />

3u(t) =<br />

t · u(t) =<br />

0 t < 0<br />

{<br />

t − 1 t 1<br />

−(t − 1) · u(t − 1) =<br />

−5u(t − 2) =<br />

0 t < 1<br />

{<br />

t t 0<br />

0 t < 0<br />

{<br />

5 t 2<br />

0 t < 2<br />

Z njimi opišemo x(t) kot vsoto segmentov signala, ki jih določajo zgornje funkcije:<br />

t < 0 : x(t) = 0 + 0 − 0 − 0 = 0<br />

0 1 : x(t) = 3 + t − 0 − 0 = 3 + t<br />

1 2 : x(t) = 3 + t − (t − 1) − 0 = 4<br />

2 : x(t) = 3 + t − (t − 1) − 5 = −1<br />

Vidimo, da se dobljeni <strong>za</strong>pis signala ujema z grafom na sliki 2.3-1b.<br />

♦<br />

V primeru 2.3-1 smo poka<strong>za</strong>li konstrukcijo grafa signala z<br />

odseki premic. Premico, ki gre skozi točko P 1 (x 1 , t 1 ) <strong>in</strong> oklepa s<br />

pozitivno smerjo osi t (signalno osjo) kot ϕ (slika 2.3-2), določa<br />

x<br />

P1 ( x1, t1<br />

)<br />

Slika 2.3-2<br />

Graf premice<br />

x 1<br />

x 2<br />

P0 ( x0, t0<br />

)<br />

t 0<br />

x x<br />

, tan<br />

k <br />

t1 t0<br />

t<br />

t 1<br />

1 0<br />

enačba:<br />

x − x 1 = k(t − t 1 ) kjer je k = tan ϕ . (2.3-2)<br />

Strm<strong>in</strong>o premice izračunamo z izbiro dveh točk na premici, na<br />

primer P 1 (x 1 , t 1 ) <strong>in</strong> P 0 (x 0 , t 0 ):<br />

k = x 1 − x 0<br />

t 1 − t 0<br />

. (2.3-3)<br />

Tehniko pisanja enačb <strong>za</strong> funkcije, ki jih sestavljajo odseki premic,<br />

bomo poka<strong>za</strong>li na primeru signala s slike 2.3-3. Strm<strong>in</strong>e<br />

daljic, ki opisujejo signal v posameznih odsekih, so označene s<br />

k i . Zapis daljic, ki opisujejo potek signala v posameznih segmentih<br />

signala x(t), je preprost. Zanje uporabimo (2.3-2) <strong>in</strong> (2.3-3),<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.3 Matematični opisi <strong>signalov</strong> 15<br />

x<br />

P 1<br />

k<br />

k 2<br />

k<br />

P 3<br />

1<br />

2<br />

P 0<br />

t 0 t 1 t t<br />

2<br />

Slika 2.3-3<br />

Graf signala<br />

kjer <strong>za</strong> P i−1 (x i−1 , t i−1 ) izberemo <strong>za</strong>četno točko opazovanega segmenta,<br />

<strong>za</strong> P i (x i , t i ) pa končno točko. Za prvi segment signala<br />

x(t) s slike 2.3-3 lahko opišemo s poltrakom x 1 (t), lahko <strong>za</strong>pišemo:<br />

x(t 0 < t t 1 ) = x 1 (t) = k 1 (t − t 0 )u(t − t 0 ) , (2.3-4)<br />

kjer smo z množenjem z <strong>za</strong>maknjeno enotsko stopnico u(t − t 0 )<br />

dosegli, da je signal <strong>za</strong> t < t 0 enak nič. Drugi segment signala<br />

x(t), ki leži nad <strong>in</strong>tervalom t 1 < t t 2 , opišemo s poltrakom,<br />

ki ima izhodišče v točki s koord<strong>in</strong>atami (0, t 1 ) ter v <strong>in</strong>tervalu<br />

t 1 t < t 2 poteka vzporedno s signalom. Zanj velja:<br />

x 2 (t) = k 2 (t − t 1 )u(t − t 1 ) . (2.3-5)<br />

Hitro uvidimo, da seštevek x 1 (t)+x 2 (t) ne da poteka signala x(t)<br />

na <strong>in</strong>tervalu (t 0 , t 2 ). Da ga dobimo, moramo v drugem segmentu<br />

odšteti še strm<strong>in</strong>o signala, ki prihaja iz prvega segmenta:<br />

x(t 1 < t t 2 ) = x 1 (t) −x 1 (t − t 1 ) +x<br />

} {{ } 2 (t)<br />

izravnava x 1 (t)<br />

= k 1 (t − t 0 )u(t − t 0 ) − k 1 (t − t 1 )u(t − t 1 )<br />

+ k 2 (t − t 1 )u(t − t 1 )<br />

= k 1 (t − t 0 )u(t − t 0 ) + (k 2 − k 1 )(t − t 1 )u(t − t 1 ) .<br />

Podobno velja <strong>za</strong> tretji segment. Poltrak, ki poteka skozi točko<br />

s koord<strong>in</strong>atama (0, t 2 ) <strong>in</strong> vzporedno s potekom signala v tem<br />

segmentu, je:<br />

x 3 (t) = k 3 (t − t 2 ) · u(t − t 2 ) .<br />

Potek signala v tem segmentu pa izračunamo na enak nač<strong>in</strong> kot<br />

v drugem segmentu. Torej:<br />

x(t 2 t < t 3 ) = x 1 (t) − x 1 (t − t 1 ) + x 2 (t) − x<br />

} {{ } 2 (t − t 2 ) + x 3 (t)<br />

} {{ }<br />

x(t 1 t


16 2. Signali<br />

<strong>in</strong> po upoštevanju (2.3-4) <strong>in</strong> (2.3-5) uredimo v:<br />

x(t) = k 1 (t − t 0 )u(t − t 0 ) + (k 2 − k 1 )(t − t 1 )u(t − t 1 )−<br />

k 2 (t − t 2 )u(t − t 2 ) + k 3 (t − t 2 )u(t − t 2 )<br />

= k 1 (t − t 0 )u(t − t 0 ) + (k 2 − k 1 )(t − t 1 )u(t − t 1 )+<br />

(k 3 − k 2 )(t − t 2 )u(t − t 2 ) .<br />

Dobljeni rezultat velja na splošno <strong>za</strong> signale, ki jih lahko<br />

opišemo z daljicami. Tako vidimo, da lahko tudi <strong>za</strong>htevne oblike<br />

opišemo z vsoto členov, ki opisujejo potek posameznih segmentov<br />

poteka signala. Pri tem moramo poudariti:<br />

Prika<strong>za</strong>na tehnika matematičnega opisa <strong>signalov</strong> velja le <strong>za</strong> determ<strong>in</strong>istične<br />

zvezne signale.<br />

Postopek je precejšnjega praktičnega pomena, saj je na takšen<br />

nač<strong>in</strong> modeliranih mnogo fizikalnih <strong>signalov</strong>. Opis omogoča določitev<br />

mnogih lastnosti <strong>signalov</strong>, ki jih bomo opisali v nadaljevanju<br />

knjige.<br />

Seveda pa to ni ed<strong>in</strong>i nač<strong>in</strong> <strong>za</strong> opis teh <strong>signalov</strong>. V nadaljevanju<br />

te knjige bomo izpeljali še kvantitativni opis <strong>signalov</strong> z<br />

baznimi funkcijami, na katerem prav<strong>za</strong>prav sloni več<strong>in</strong>a teorije<br />

obdelave <strong>signalov</strong>. Njej je posvečena vsa druga knjiga Teorije<br />

<strong>signalov</strong>. V njej so opisane Fourierove vrste ter Fourierova <strong>in</strong> Laplaceova<br />

transformacija. Opisano metodo matematičnega opisa<br />

<strong>signalov</strong> utrdimo še z zgledom:<br />

ZGLED 2.3-2 Enačba žagastega signala<br />

Žagasti signal je zelo pogosta oblika signala v elektroniki. Z napetostjo žagaste<br />

oblike na primer odklanjamo elektronski curek v <strong>za</strong>slonski cevi osciloskopa. Podobno<br />

funkcijo ima tok žagaste oblike pri televizijskih ekranih. Primer žagaste<br />

napetosti kaže slika 2.3-4.<br />

Slika 2.3-4<br />

Graf žagaste napetosti<br />

V<br />

x( t)<br />

0<br />

-T<br />

0 T 2T t<br />

S slike vidimo, da je signal periodičen s periodo T .<br />

določena s:<br />

k = V T ,<br />

Strm<strong>in</strong>a signala je<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.4 Elementarne operacije nad signali 17<br />

traja pa čas periode T . Za segment med 0 <strong>in</strong> T lahko torej pišemo:<br />

x 1(t) = V t[u(t) − u(t − T )]<br />

T<br />

( )<br />

= V t −<br />

T<br />

T t rect 2<br />

.<br />

T<br />

x 1(t) je def<strong>in</strong>iran na <strong>in</strong>tervalu [0, T ]. Ker je periodični signal, lahko njegov <strong>za</strong>pis<br />

posplošimo <strong>za</strong> <strong>in</strong>terval [nT , (n + 1)T ]:<br />

x n(t) = x 1(t − nT ) = = V (t − nT )[u(t − nT ) − u(t − (n + 1)T )]<br />

T<br />

( )<br />

= V t − nT −<br />

T<br />

T (t − nT ) rect 2<br />

. (2.3-6)<br />

T<br />

Zgornja enačba velja tako pri pozitivnih kot pri negativnih n. Če seštejemo <strong>za</strong>pise<br />

<strong>za</strong> posamezne segmente, dobimo:<br />

x(t) =<br />

∞∑<br />

n=−∞<br />

x 1(t − nT ) ,<br />

kjer je x 1(t − nT ) opisan v (2.3-6).<br />

♦<br />

2.4 Elementarne operacije<br />

nad časovnimi signali<br />

V tem razdelku so opisane osnovne operacije, s katerimi preoblikujemo<br />

signale. Razlikovali bomo operacije, ki se izvajajo na<br />

(enem) signalu <strong>in</strong> operacije nad dvema signaloma. Med prve operacije<br />

spadajo na primer transformacija signalnega <strong>in</strong> časovnega<br />

območja, vzorčenje, <strong>in</strong>terpolacija <strong>in</strong> kvanti<strong>za</strong>cija, med druge pa<br />

seštevanje, odštevanje, množenje <strong>in</strong> deljenje dveh <strong>signalov</strong>.<br />

2.4.1 Transformacija signalnega območja<br />

S transformacijo signalnega območja spremenimo amplitudni razmah<br />

signala. Transformacijo signalnega območja imenujemo tudi<br />

amplitudna transformacija. Zgodi se na primer pri ojačanju signala.<br />

Nekateri ojačevalniki lahko poleg ojačanja premaknejo še<br />

signalno območje s prištetjem ali odštetjem konstantnega (enosmernega)<br />

signala, lahko pa spremenijo tudi predznak signala.


18 2. Signali<br />

DEFINICIJA 2 (Transformacija signalnega območja)<br />

Transformacija signalnega območja je preslikava T : A star →<br />

A nov , kjer je A star orig<strong>in</strong>alni amplitudni razmah signala, ki izvede<br />

skaliranje <strong>in</strong>/ali premaknitev orig<strong>in</strong>alnega signalnega območja.<br />

Transformacija T torej preslika signal x star ∈ A T star v signal<br />

x nov ∈ A T nov. Def<strong>in</strong>irana je z:<br />

x nov (t) = T [x star (t)] ,<br />

t ∈ T<br />

= A · x star + B . (2.4-1)<br />

kjer je A skalirni faktor <strong>in</strong> B premik signalnega območja.<br />

ZGLED 2.4-1 Transformacija signalnega območja<br />

Izvedimo amplitudno transformacijo signala rect(t) z izbiro: A = −2 <strong>in</strong> B = 0,<br />

<strong>in</strong> A = −2 <strong>in</strong> B = 2.<br />

1<br />

x star<br />

1<br />

x star<br />

0 0<br />

-1/2 0 1/2 t<br />

A 2<br />

-1/2 0 1/2<br />

2 2<br />

x nov<br />

x nov<br />

0 0<br />

t<br />

A 2<br />

B 2<br />

-1 -1<br />

-2 -2<br />

-1/2 0 1/2 t<br />

-1/2 0 1/2<br />

Slika 2.4-1<br />

Amplitudna transformacija: levo: A = −2 <strong>in</strong> B = 0, desno: A = −2 <strong>in</strong> B = 2.<br />

t<br />

S slike 2.4-1 vidimo, da skaliranje izvede množilnik, premik pa seštevalnik.<br />

Oba elementa sta pri analognih signalih ponavadi izvedena z ustrezno ve<strong>za</strong>vo<br />

operacijskega ojačevalnika.<br />

♦<br />

2.4.2 Kvanti<strong>za</strong>cija<br />

Pomembna transformacija amplitudnega razmaha signala je kvanti<strong>za</strong>cija.<br />

To transformacijo potrebujemo, ko želimo obdelati signal<br />

z digitalnim računalnikom. Digitalni računalniki lahko obravnavajo<br />

le signale s končnim številom možnih vrednosti amplitude.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.4 Elementarne operacije nad signali 19<br />

DEFINICIJA 3 (Kvanti<strong>za</strong>cija)<br />

Vsaka transformacija (končnega) amplitudnega razmaha signala,<br />

ki povzroči, da ima A nov končno število vrednosti, se imenuje<br />

kvanti<strong>za</strong>cija.<br />

ZGLED 2.4-2 Uniformna kvanti<strong>za</strong>cija<br />

Predpostavimo, da velja A star = R (torej transformacija <strong>za</strong>jema svoje vrednosti<br />

iz <strong>za</strong>loge realnih števil). Transformacijo T : R → Z + · q amplitudnega razmaha<br />

signala izvedemo tako, da bo nov amplitudni razmah signala določen z A nov =<br />

Z + · q = {0, q, 2q, · · · , Nq}. V tem primeru je transformacija T določena z:<br />

⎧<br />

⎪⎨ ⌊<br />

0 x < 0,<br />

x<br />

x q =<br />

q + 1 ⌋<br />

0 x < (N − 1<br />

2 )q,<br />

(2.4-2)<br />

2<br />

⎪⎩<br />

Nq x (N − 1 )q, 2<br />

kjer ⌊ ⌋ označuje funkcijo, ki <strong>za</strong>okroži realno število x na največje naravno število,<br />

ki je enako ali manjše od x; N je celo število, q pa je pozitivno število, ki ga<br />

imenujemo kvanti<strong>za</strong>cijski <strong>in</strong>terval ali tudi kvanti<strong>za</strong>cijski korak. Grafični prikaz<br />

transformacije T je na sliki 2.4-2.<br />

Nq<br />

T<br />

x q<br />

( N -1)<br />

q<br />

2q<br />

q<br />

0<br />

0 q 2q<br />

( N - 1+<br />

1 ) q<br />

2<br />

( N -1-<br />

1 ) q ( N -1)<br />

q<br />

2<br />

x<br />

x ŒR<br />

x Z ◊ q<br />

q Œ<br />

+<br />

Slika 2.4-2<br />

Transformacija amplitudnega<br />

razmaha signala pri uniformni<br />

kvanti<strong>za</strong>ciji<br />

Prika<strong>za</strong>na kvanti<strong>za</strong>cija odreže vse vrednosti signala, ki so manjše od 1 2 q<br />

<strong>in</strong> omeji maksimalno vrednost na Nq. Uč<strong>in</strong>ek uniformne kvanti<strong>za</strong>cije z osmimi<br />

kvanti<strong>za</strong>cijskimi <strong>in</strong>tervali na signal x(t) = 1/2[1 − cos(2πt)], t ∈ R kaže slika<br />

1 1<br />

signal: x( t)<br />

<br />

kvantiziran<br />

Slika 2.4-3<br />

1<br />

cos( t)<br />

signal: xq t Uniformna kvanti<strong>za</strong>cija, levo:<br />

2<br />

2.4-3. Tu je N = 8 <strong>in</strong> q = 1 . 8<br />

0 0<br />

orig<strong>in</strong>alni signal, desno:<br />

kvantizirani signal<br />

0 èas<br />

1 0 èas<br />

1


20 2. Signali<br />

To kvanti<strong>za</strong>cijo imenujemo uniformna <strong>za</strong>to, ker so vsi kvanti<strong>za</strong>cijski <strong>in</strong>tervali<br />

enaki. Poleg uniformne kvanti<strong>za</strong>cije poznamo <strong>in</strong> uporabljamo še neuniformne<br />

kvanti<strong>za</strong>cije. Pri njih je kvanti<strong>za</strong>cijski <strong>in</strong>terval odvisen od amplitude signala.<br />

Ponavadi je pri malih signalih kvanti<strong>za</strong>cijski <strong>in</strong>terval mali, pri velikih pa veliki.<br />

Neuniformna kvanti<strong>za</strong>cija se pogosto uporablja pri digitalnem prenosu govora. ♦<br />

2.4.3 Transformacija signalne osi<br />

S to transformacijo skrčimo ali podaljšamo domeno, iz katere<br />

<strong>za</strong>jemamo vrednosti argumenta funkcije, s katero opišemo signal.<br />

DEFINICIJA 4 (Transformacija signalne osi)<br />

Transformacija signalne osi T star , nad katero je def<strong>in</strong>irano območje<br />

signala A, je preslikava T : T star → T nov , kjer je T nov nova signalna<br />

os.<br />

x nov (t) = x star ( T[t] ) , <strong>za</strong> vse t ∈ T star .<br />

Transformacija je bijektivna, če obstaja <strong>in</strong>verzna transformacija<br />

T −1 : T nov → T star . Ta preslika x nov ∈ A Tnov v signal x star ∈<br />

tako, da velja:<br />

A Tstar<br />

x star (t) = x nov<br />

(<br />

T −1 [t] ) , <strong>za</strong> vse t ∈ T nov .<br />

Primer transformacije signalne osi je T(t) = 1 at, a ≠ 0, pri<br />

kateri je <strong>in</strong>verzna transformacija enaka T −1 (t) = at. Pri tem je<br />

a poljubna neničelna konstanta. V tem primeru velja:<br />

(<br />

x nov (t) = x star (a · t) <strong>in</strong> x star = x t<br />

)<br />

nov a<br />

.<br />

ZGLED 2.4-3 Razširitev, skrčenje <strong>in</strong> <strong>za</strong>suk časovne osi<br />

Grafično predstavitev uč<strong>in</strong>ka izbire konstante a pri transformaciji signalne osi pri<br />

zveznih <strong>in</strong> diskretnih signalih podaja slika 2.4-4 na strani 21.<br />

Pri skaliranju signalne osi velja x nov = x star(at), <strong>za</strong>to se pri a > 1 signal<br />

razširi, pri 0 < a < 1 pa skrči. Pri negativnih vrednostih a dobimo poleg skrčitve<br />

ja (|a| > 1, a < 0) ali razširitve (0 < |a| < 1, a < 0) še <strong>za</strong>suk signala. ♦<br />

2.4.4 Premik po signalni osi<br />

Premik signala je transformacija, pri kateri premaknemo signal<br />

po signalni osi. Inverzna transformacija vrne signal v staro lego.<br />

Pri časovnih signalih sta določeni z:<br />

T(t) = t + a t ∈ T star<br />

(prehitevanje) (2.4-3)<br />

T −1 (t) = t − a t ∈ T nov<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.4 Elementarne operacije nad signali 21<br />

x<br />

1<br />

t<br />

star( )<br />

orig<strong>in</strong>alni signal<br />

časovno zvezni signal<br />

0<br />

0 b t<br />

razširjeni signal<br />

skrèeni signal<br />

<strong>za</strong>sukani signal<br />

1<br />

1 1<br />

x ( nov<br />

t ) a 1 0 a 1<br />

x ( nov<br />

t ) x ( t)<br />

nov<br />

a 1<br />

0<br />

0<br />

a b<br />

t<br />

0 0<br />

0 a b t<br />

-b 0<br />

t<br />

orig<strong>in</strong>alni signal<br />

1<br />

xstar ( n)<br />

0<br />

0 b<br />

t<br />

t<br />

časovno diskretni signal<br />

razširjen signal skrèen signal <strong>za</strong>sukan signal<br />

1<br />

a 2<br />

1 a 0, 5 1<br />

a 1<br />

xnov ( n) xnov ( n) xnov ( n)<br />

0 0 0<br />

a b<br />

t a b<br />

t 0 t<br />

0 at<br />

0 at<br />

b t<br />

Slika 2.4-4<br />

Razširitev, skrčenje <strong>in</strong> <strong>za</strong>suk časovne osi žagastemu pulzu, glej zgled 2.4-3<br />

oziroma:<br />

T(t) = t − a t ∈ T star<br />

T −1 (t) = t + a t ∈ T nov<br />

(<strong>za</strong>kasnitev) , (2.4-4)<br />

kjer je a pozitivna realna konstanta.<br />

translacije signala.<br />

Slika 2.4-5 kaže primer


22 2. Signali<br />

orig<strong>in</strong>alni signal<br />

1<br />

xstar( t) xstar( t) xnov( t a)<br />

Slika 2.4-5<br />

Premik signala po signalni osi<br />

0<br />

0 b t<br />

premaknjeni (<strong>za</strong>kasnjeni) signal<br />

1<br />

x ( nov<br />

t )<br />

x ( t) x ( t a)<br />

nov<br />

star<br />

0<br />

0<br />

a<br />

b+a<br />

t<br />

2.5 Soda <strong>in</strong> liha simetričnost signala<br />

Simetrije v izgledu signala so v analizi <strong>in</strong> s<strong>in</strong>tezi <strong>signalov</strong> zelo pomembne.<br />

V mnogih primerih poenostavijo računanje, poleg tega<br />

pa nam te značilnosti pomagajo razumeti ali doka<strong>za</strong>ti določene<br />

lastnosti <strong>signalov</strong>.<br />

Simetrije lahko določimo le <strong>za</strong> determ<strong>in</strong>istične signale. Glede<br />

na simetrijo preko vertikalne osi koord<strong>in</strong>atnega sistema, v katerem<br />

predstavimo graf signala, ločimo dve osnovni skup<strong>in</strong>i <strong>signalov</strong>:<br />

Sode signale oziroma funkcije, <strong>za</strong> katere velja:<br />

x(t) = x(−t) (2.5-1)<br />

Iz (2.5-1) sledi, da je potek signala pri t 0 zrcalna slika<br />

poteka signala pri t 0 (slika 2.5-1).<br />

Slika 2.5-1<br />

Primer sodo-simetričnega<br />

signala<br />

1<br />

x( t)<br />

0<br />

-a<br />

0 a t<br />

Primer znane sode funkcije je cos(ωt), saj velja cos(ωt) =<br />

cos(−ωt).<br />

Lihe signale oziroma funkcije, <strong>za</strong> katere velja:<br />

x(t) = −x(−t) . (2.5-2)<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.5 Soda <strong>in</strong> liha simetričnost signala 23<br />

Lihi signali oziroma funkcije so simetrične glede na koord<strong>in</strong>atno<br />

izhodišče (slika 2.5-2).<br />

1<br />

x( t)<br />

0<br />

-1<br />

-a<br />

0 a t<br />

Slika 2.5-2<br />

Primer liho-simetričnega<br />

signala<br />

Primer znane lihe funkcije je s<strong>in</strong>(ωt), saj velja s<strong>in</strong>(ωt) = − s<strong>in</strong>(−ωt).<br />

Vsako funkcijo lahko razstavimo na lihi <strong>in</strong> sodi del:<br />

x(t) = x lih (t) + x sod (t) (2.5-3)<br />

kjer je x lih liha funkcija <strong>in</strong> x sod soda funkcija. Z <strong>za</strong>menjavo argumenta<br />

t z −t iz (2.5-3) z upoštevanjem (2.5-1) <strong>in</strong> (2.5-2), sledi:<br />

x(−t) = x lih (−t) + x sod (−t) = −x lih (t) + x sod (t) (2.5-4)<br />

Seštevek (2.5-3) <strong>in</strong> (2.5-4) da:<br />

x(t) + x(−t) = x lih (t) + x sod (t) − x lih (t) + x sod (t) = 2x sod (t)<br />

oziroma<br />

x sod (t) = 1 [x(t) + x(−t)] , (2.5-5)<br />

2<br />

iz razlike (2.5-3) <strong>in</strong> (2.5-4) pa dobimo:<br />

Razdelek <strong>za</strong>ključimo s pregledom lastnosti lihih <strong>in</strong> sodih funkcij:<br />

x lih (t) = 1 [x(t) − x(−t)] . (2.5-6)<br />

2<br />

Vsota dveh sodih <strong>signalov</strong> je sodi signal.<br />

Vsota dveh lihih <strong>signalov</strong> je lihi signal.<br />

Vsota sodega <strong>in</strong> lihega signala ni ne sodi ne lihi signal.<br />

Produkt dveh sodih <strong>signalov</strong> je sodi signal.<br />

Produkt dveh lihih <strong>signalov</strong> je lihi signal.<br />

Produkt sodega <strong>in</strong> lihega signala je lihi signal.


24 2. Signali<br />

2.6 L<strong>in</strong>earni prostor<br />

V tem podpoglavju bomo predstavili različne množice <strong>signalov</strong>,<br />

ki bodo pomembne v nadaljnji obravnavi <strong>signalov</strong>. Vse imajo<br />

strukturo l<strong>in</strong>earnega prostora. L<strong>in</strong>earni prostor, preprosto povedano,<br />

je množica z naslednjimi lastnostmi:<br />

1. Vsota dveh funkcij, ki pripadata temu razredu funkcij, je prav<br />

tako v tem razredu.<br />

2. Zmnožek skalarja <strong>in</strong> funkcije iz tega razreda funkcij je tudi pripadnik<br />

tega razreda.<br />

ZGLED 2.6-1 L<strong>in</strong>earni prostor R N <strong>in</strong> C N<br />

Zelo znani l<strong>in</strong>earni prostor je množica vseh N-teric x = (x 1, x 2, · · · , x N ) z<br />

x k ∈ R, k = 1, 2, · · · N. Če so N-terice sestavljene iz realnih števil, jo označimo<br />

z R N , če pa so sestavljene iz kompleksnih števil, pa s C N . V R N oziroma C N<br />

je vsota dveh elementov x = (x 1, x 2, · · · , x N ) <strong>in</strong> y = (y 1, y 2, · · · , y N ) enaka:<br />

x + y = (x 1 + y 1, x 2 + y 2, · · · , x N + y N )<br />

množenje elementa s skalarjem pa je def<strong>in</strong>irana z:<br />

αx = (αx 1, αx 2, · · · , αx N )<br />

Torej je množica vseh R N l<strong>in</strong>earni prostor nad množico realnih števil. Analogno<br />

je množica vseh C N l<strong>in</strong>earni prostor nad množico kompleksnih števil. ♦<br />

V obravnavi <strong>signalov</strong> se bomo pogosto sklicevali na l<strong>in</strong>earni<br />

prostor vseh časovnih <strong>za</strong>poredij (časovno diskretnih <strong>signalov</strong>) <strong>in</strong><br />

prostor vseh časovno zveznih <strong>signalov</strong>. Njuni def<strong>in</strong>iciji sta:<br />

DEFINICIJA 5 (Prostor časovnih <strong>za</strong>poredij l)<br />

Množica vseh kompleksnih <strong>za</strong>poredij<br />

x = (· · · , x(−1), x(0), x(1), · · · )<br />

je l<strong>in</strong>earni prostor l nad kompleksnimi števili, če velja:<br />

(x + y)(n) = x(n) + y(n)<br />

(αx)(n) = αx(n) .<br />

Ta l<strong>in</strong>earni prostor je razširitev prostora C.<br />

n ∈ Z<br />

DEFINICIJA 6 (Prostor zveznih časovnih funkcij L)<br />

Množica vseh kompleksnih časovnih <strong>signalov</strong> določa l<strong>in</strong>earni prostor<br />

L, če <strong>za</strong> vsako točko funkcije velja:<br />

(x + y)(t) = x(t) + y(t)<br />

(αx)(t) = αx(t) .<br />

t ∈ R<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.7 Norme <strong>signalov</strong> 25<br />

2.7 Norme <strong>signalov</strong><br />

Mnogokrat želimo signale jedrnato opisati. Katere lastnosti signala<br />

pa so tiste, ki so karakteristične <strong>za</strong> določeni signal? Med<br />

fizikalnimi količ<strong>in</strong>ami so na primer <strong>za</strong>nimive amplituda, moč <strong>in</strong><br />

energija. Te količ<strong>in</strong>e ponavadi podajamo kot parametre signala.<br />

Na splošno pa so parametri, s katerimi opišemo značilnost signala,<br />

norme.<br />

Normo signala x označimo z ‖x‖. Ker poznamo mnogo norm,<br />

bomo takrat, ko želimo poudariti, <strong>za</strong> katero normo gre, to označili<br />

z ustreznim <strong>in</strong>deksom ‖x‖ p . In kaj so norme? V splošnem <strong>za</strong><br />

norme velja:<br />

DEFINICIJA 7 (Norma)<br />

Norma ‖ · ‖ p elementa x iz l<strong>in</strong>earnega prostora X je nenegativno<br />

realno število, ki je nič le takrat <strong>in</strong> samo takrat, ko je x ničelni<br />

element.<br />

DEFINICIJA 8 (Homogenost norme)<br />

Norma je homogena mera (glede na skaliranje) takrat, ko velja:<br />

‖λx‖ = |λ| · ‖x‖ (2.7-1)<br />

<strong>za</strong> vsak skalar λ <strong>in</strong> vsak element x signalnega prostora, hkrati pa<br />

izpolni trikotniški izrek:<br />

<strong>za</strong> vsak x <strong>in</strong> y iz l<strong>in</strong>earnega prostora.<br />

‖x + y‖ ‖x‖ + ‖y‖ (2.7-2)<br />

V zgornjih def<strong>in</strong>icijah smo povedali naslednje: (i) rezultat<br />

izračuna norme je nenegativno število, ki kvantitativno določa<br />

lastnost, ki jo z izbrano normo merimo; (ii) norma je homogena<br />

(enakomerna) mera. To pomeni, da s transformacijo signala, ki<br />

<strong>za</strong>jema skaliranje amplitudnega ali signalnega območja, spremenimo<br />

tudi normo signala. Ko transformacija povzroči le premik<br />

signala po signalni osi ali <strong>za</strong>suk njegovega signalnega območja, se<br />

velikosti norme ne spremeni. Veljavnost <strong>in</strong> pomen trikotniškega<br />

izreka bomo uvideli iz def<strong>in</strong>icije norm.<br />

DEFINICIJA 9 (Norme v prostorih R N <strong>in</strong> C N )<br />

Naj bo p celo pozitivno število, 1 p ∞. Potem je p-norma


26 2. Signali<br />

elementa x določena z:<br />

⎧( ∑N<br />

) 1/p<br />

⎪⎨ k=1 |x k| p <strong>za</strong> 1 p < ∞<br />

‖x‖ p =<br />

⎪⎩ max |x k| <strong>za</strong> p = ∞<br />

1kN<br />

(2.7-3)<br />

ZGLED 2.7-1<br />

Izračunajmo norme ‖ · ‖ 1, ‖ · ‖ 2 <strong>in</strong> ‖ · ‖ ∞ <strong>za</strong> element l<strong>in</strong>earnega prostora x =<br />

(1, j), x ∈ C 1 .<br />

‖x‖ 1 = ( |1| 1 + |j| 1) 1/1<br />

= 1 + 1 = 2<br />

‖x‖ 2 = ( |1| 2 + |j| 2) 1/2<br />

=<br />

√<br />

1 + 1 =<br />

√<br />

2<br />

‖x‖ ∞ = max (|1|, |j|) = max(1, 1) = 1<br />

k<br />

Normo ‖x‖ 2 si lahko v prostoru C 1 enostavno predstavimo: določa dolž<strong>in</strong>o kompleksnega<br />

ka<strong>za</strong>lca x s koord<strong>in</strong>atami (1, j) (slika 2.7).<br />

♦<br />

Slika 2.7-1<br />

j<br />

Evklidska norma, Levo v<br />

prostoru C 2 , desno v prostoru<br />

R 3 | x 1 | 1<br />

<br />

| x 2 | 1<br />

x 1<br />

2 3<br />

1 x 2<br />

1<br />

x 3<br />

1<br />

Normo p = 2 v prostorih R N imenujemo tudi Evklidska norma.<br />

V prostorih R 2 <strong>in</strong> R 3 ima Evklidska norma predočljivo geometrijsko<br />

predstavo (slika 2.7).<br />

O Grškem matematiku Evklidu (cca 365 do cca 300<br />

let p.n.š.) je le malo znanega. Njegovo glavno delo<br />

Elementi ima logično zgradbo, ki je bila vzor <strong>za</strong> resno<br />

matematično obravnavo problemov vse do 19. stoletja.<br />

Def<strong>in</strong>icije norm v prostorih R N <strong>in</strong> C N lahko preprosto posplošimo<br />

na prostor l. Pri tem velja:<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.7 Norme <strong>signalov</strong> 27<br />

⎧( ∞<br />

) 1/p<br />

∑<br />

⎪⎨ |x n | p <strong>za</strong> 1 p < ∞ če vrsta obstaja<br />

‖x‖ p = n=−∞<br />

⎪⎩ max |x n | <strong>za</strong> p = ∞ če obstaja maksimum<br />

k<br />

(2.7-4)<br />

2.7.1 Norme v prostoru L<br />

Kot neskončno dimenzionalne vektorje si lahko predstavljamo<br />

tudi časovno zvezne funkcije z def<strong>in</strong>icijsko domeno na <strong>in</strong>tervalu<br />

(a, b) ∈ R. Če predpostavimo, da vsaka točka def<strong>in</strong>icijskega<br />

območja določa eno dimenzijo vektorja, vrednost funkcije v tej<br />

točki pa komponento vektorja, je funkcija ”neskončno dimenzionalni<br />

vektor”. Ker pa so def<strong>in</strong>icijske točke <strong>in</strong>f<strong>in</strong>itezimalno<br />

blizu, je v normah v prostoru množic zveznih <strong>signalov</strong> L operacija<br />

seštevanja <strong>za</strong>menjana z <strong>in</strong>tegriranjem. Pri tem velja:<br />

DEFINICIJA 10 (Norme zveznih funkcij)<br />

p−norma elementa x iz prostora L množice zveznih funkcij je<br />

def<strong>in</strong>irana z:<br />

⎧( ∫ )<br />

∞<br />

1/p<br />

⎪⎨ −∞ |x(t)|p <strong>za</strong> 1 p < ∞<br />

‖x‖ p =<br />

⎪⎩<br />

sup<br />

t∈R<br />

|x t | <strong>za</strong> p = ∞<br />

S ”sup” smo označili supremum, to je zgornjo mejo signala.<br />

Zgornja meja signala je enaka najmanjšemu realnemu številu α,<br />

<strong>za</strong> katerega velja:<br />

|x(n)| α n ∈ Z oziroma |x(t)| α t ∈ R .<br />

2.7.2 Fizikalni pomen norm<br />

Norme, ki jih najpogosteje uporabljamo, so ‖ · ‖ 1 , ‖ · ‖ 2 <strong>in</strong> ‖ · ‖ ∞ .<br />

Te norme imajo fizikalni pomen. Določajo naslednje:<br />

1. Norma ‖ · ‖ ∞ določi (maksimalno) amplitudo signala.<br />

2. Z normo ‖ · ‖ 2 lahko določimo energijo signala:<br />

∞∑<br />

‖x(n)‖ 2 2 = |x(n)| 2 n ∈ Z (2.7-5)<br />

i=−∞


28 2. Signali<br />

orig<strong>in</strong>alni signal<br />

1<br />

x<br />

0<br />

t<br />

oziroma<br />

‖x(t)‖ 2 2 =<br />

∫ ∞<br />

−∞<br />

|x(t)| 2 dt t ∈ R (2.7-6)<br />

Ime energija signala opravičujemo s podobnostjo obrazcev (2.7-5)<br />

<strong>in</strong> (2.7-6) obrazcu <strong>za</strong> izračun energije električnega signala (napetosti<br />

ali toka), ki se troši na bremenu 1 ohma.<br />

3. Norma ‖ · ‖ 1 določa vitalnost signala:<br />

∞∑<br />

‖x(n)‖ 1 = |x(n)| n ∈ Z (2.7-7)<br />

oziroma<br />

‖x(t)‖ 1 =<br />

i=−∞<br />

∫ ∞<br />

−∞<br />

|x(t)| dt t ∈ R . (2.7-8)<br />

Grafično predstavitev vseh treh norm kaže slika 2.7-2.<br />

norma x <br />

norma x 1<br />

norma x 2<br />

površ<strong>in</strong>a =<br />

1<br />

amplituda 1 vitalnost signala 1<br />

|x| |x| |x| 2<br />

0<br />

Slika 2.7-2<br />

Amplituda, vitalnost <strong>in</strong> energija signala<br />

t<br />

0 0<br />

t<br />

površ<strong>in</strong>a =<br />

energija signala<br />

t<br />

2.7.3 Moč signala, efektivna <strong>in</strong><br />

povprečna vrednost signala<br />

Poleg amplitude, energije <strong>in</strong> vitalnosti signala so mnogokrat <strong>za</strong>nimive<br />

tudi druge količ<strong>in</strong>e, ki jih lahko izpeljemo iz norm signala.<br />

Med njimi so posebej pomembne moč (trenutna <strong>in</strong> povprečna),<br />

efektivna vrednost signala <strong>in</strong> povprečna vrednost. Te velič<strong>in</strong>e<br />

poznamo že iz fizike <strong>in</strong> osnov elektrotehnike, <strong>za</strong>to bomo v nadaljevanju<br />

le povzeli njihove def<strong>in</strong>icije.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.7 Norme <strong>signalov</strong> 29<br />

Trenutna moč signala Opazujemo zvezni kompleksni signal<br />

x(t):<br />

x(t) = a(t) + jb(t) , (2.7-9)<br />

kjer sta a(t) <strong>in</strong> b(t) realni funkciji, ki opisujeta obliko signala.<br />

Trenutna moč zveznega signala p(t) je določena z:<br />

p(t) = [a(t)] 2 + [b(t)] 2 = |x(t)| 2 . (2.7-10)<br />

Do zgornje enačbe lahko pridemo po analogiji z električno močjo<br />

pri znanem toku, ko je upornost bremena 1 ohm (slika 2.7-3). S<br />

faktori<strong>za</strong>cijo:<br />

[a(t)] 2 + [b(t)] 2 = [a(t) + jb(t)][a(t) − jb(t)]<br />

lahko trenutno moč izrazimo s konjugirano kompleksnim parom:<br />

p(t) = x(t)x ∗ (t) = |x(t)| 2 . (2.7-11)<br />

p( t)<br />

x( t)<br />

a( t)<br />

pabtg<br />

b g<br />

p( t)<br />

= x t ¥ 1<br />

= p ( t) + p ( t)<br />

a<br />

b g<br />

jb( t) p t b t<br />

pb ( t)<br />

2<br />

b g<br />

b<br />

2 2<br />

= a t + b t<br />

b ( ) = b g 2<br />

ab g<br />

b g<br />

p t = a t<br />

2<br />

Slika 2.7-3<br />

Analogija med močjo signala<br />

<strong>in</strong> močjo električne energije<br />

Trenutna medsebojna moč dveh <strong>signalov</strong> Imejmo dva kompleksna<br />

signala x <strong>in</strong> y. Njuna trenutna medsebojna moč p xy oziroma<br />

p yx je def<strong>in</strong>irana z:<br />

Časovno zvezni signali:<br />

Časovno diskretni signali:<br />

p xy (t) = x(t)y ∗ (t)<br />

(2.7-12a)<br />

p yx (t) = x ∗ (t)y(t) , t ∈ R . (2.7-12b)<br />

p xy (n) = x(n)y ∗ (n)<br />

(2.7-13a)<br />

p yx (n) = x ∗ (n)y(n) n ∈ Z . (2.7-13b)


30 2. Signali<br />

Tako <strong>za</strong> zvezne <strong>in</strong> časovno diskretne signale velja zve<strong>za</strong>:<br />

Če sta oba signala realna seveda velja:<br />

p xy = p ∗ yx . (2.7-14)<br />

p xy = p yx = xy . (2.7-15)<br />

Povprečna moč Povprečno moč signala na <strong>in</strong>tervalu dolž<strong>in</strong>e<br />

T def<strong>in</strong>iramo kot povprečno energijo signala na tem <strong>in</strong>tervalu.<br />

Označimo jo s P .<br />

Časovno zvezni signali:<br />

Časovno diskretni signali:<br />

P (t 0 , T ) = 1 T<br />

∫ t0 +T<br />

t 0<br />

|x(t)| 2 dt = 1 T ‖x(t)‖2 2<br />

(2.7-16)<br />

P (N 0 , N) = 1 N<br />

N∑<br />

0 +N<br />

n=N 0<br />

|x(n)| 2 = 1 N ‖x(n)‖2 2<br />

(2.7-17)<br />

V zgornji def<strong>in</strong>iciji smo upoštevali časovno omejen signal s<br />

končno energijo. Z nekaj razmisleka zlahka uvidimo, da sta to<br />

pomembni omejitvi, ki ju vsi signali ne izpolnijo.<br />

OPOMBA 2-2 Zaradi krajšega pisanja bomo nadalje uporabljali naslednji<br />

oznaki. Za <strong>in</strong>tegriranje preko <strong>in</strong>tervala T :<br />

∫ t0 +T<br />

∫<br />

=<br />

t 0 T<br />

<strong>za</strong> seštevanje N diskretnih vrednosti<br />

N 0 +N<br />

∑<br />

n=N 0<br />

= ∑ N<br />

Pripadajoče moči pa bomo pisali brez argumenta.<br />

Z upoštevanjem (2.7-11) lahko (2.7-16) <strong>za</strong>pišemo tudi v naslednji<br />

obliki:<br />

P xy = 1 ∫<br />

x(t)x ∗ (t) dt . (2.7-18)<br />

T<br />

T<br />

Povprečno medsebojno moč dveh <strong>signalov</strong> dobimo z ustreznim<br />

upoštevanjem obeh <strong>signalov</strong> v zgornji enačbi:<br />

P xy = 1 ∫<br />

x(t)y ∗ (t) dt (2.7-19)<br />

T T<br />

P yx = 1 ∫<br />

x ∗ (t)y(t) dt , (2.7-20)<br />

T<br />

T<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.7 Norme <strong>signalov</strong> 31<br />

kjer velja P xy = Pyx. ∗ Pri realnih signalih se enačbi (2.7-19) <strong>in</strong><br />

(2.7-20) poenostavita v:<br />

P xy = P yx = 1 ∫<br />

x(t)y(t)dt . (2.7-21)<br />

T<br />

Povprečna moč časovno neomejenih <strong>signalov</strong> Če je vsaj<br />

eden od <strong>signalov</strong> x(t), y(t), t ∈ R prehoden, to se pravi, da je pred<br />

<strong>in</strong> po določenem času, na primer (t 1 , t 2 ), enak nič, ne moremo<br />

def<strong>in</strong>irati povprečne moči. V tem primeru deljenje končno velike<br />

energije z neskončno velikim <strong>in</strong>tervalom, nad katerim sta signala<br />

def<strong>in</strong>irana, enaka nič. Zato pri aperiodičnih <strong>in</strong> časovno omejenih<br />

signalih def<strong>in</strong>iramo le energijo signala.<br />

V primeru neprehodnih <strong>signalov</strong>, torej <strong>signalov</strong>, ki imajo <strong>za</strong>radi<br />

svojega trajanja neskončno energijo, pa lahko povprečno moč<br />

(če obstaja), def<strong>in</strong>iramo z:<br />

oziroma<br />

1<br />

P = lim<br />

T →∞ T<br />

1<br />

P = lim<br />

N→∞ N<br />

∫<br />

T<br />

T<br />

|x(t)| 2 dt t ∈ R (2.7-22)<br />

∑<br />

|x(n)| 2 n ∈ Z . (2.7-23)<br />

N<br />

Efektivna vrednost signala Efektivno vrednost signala bomo<br />

označili s x rms . Pri oznaki smo si <strong>in</strong>deks sposodili pri angleški<br />

kratici <strong>za</strong> efektivno vrednost, r oot m ean s quare, ki slikovito<br />

pove njeno def<strong>in</strong>icijo:<br />

Časovno zvezni signali:<br />

Časovno diskretni signali:<br />

x rms (t) = lim<br />

T →∞<br />

√<br />

1<br />

T<br />

∫<br />

T<br />

|x(t)| 2 dt<br />

= √ P , t ∈ R . (2.7-24)<br />

x rms (n) = lim<br />

N→∞<br />

√<br />

1<br />

N<br />

Iz def<strong>in</strong>icije sledi, da je efektivna vrednost signala tista srednja<br />

vrednost, katere kvadrat da moč signala. Efektivna vrednost<br />

signala je vedno večja od nič.<br />

∑<br />

|x(n)| 2<br />

N<br />

= √ P , n ∈ Z . (2.7-25)


32 2. Signali<br />

Povprečna vrednost signala Povprečna vrednost signala x<br />

na <strong>in</strong>tervalu T je def<strong>in</strong>irana z:<br />

∫<br />

1<br />

T<br />

x(t) dt (2.7-26)<br />

T<br />

oziroma<br />

1<br />

N<br />

∑<br />

x(n) . (2.7-27)<br />

N<br />

Če je def<strong>in</strong>icijsko območje neskončno veliko, pri izračunu povprečne<br />

vrednosti signala v zgornjih enačbah uporabimo limitni<br />

postopek, podobno kot pri izračunu povprečne moči.<br />

Matematično gledano, je povprečna vrednost zveznega signala<br />

na nekem <strong>in</strong>tervalu povprečna vrednost določenega <strong>in</strong>tegrala.<br />

To lahko ocenimo s prvim izrekom o povprečni vrednosti,<br />

ki se glasi:<br />

IZREK 1 (Prvi izrek o povprečni vrednosti)<br />

Naj bo funkcija x(t) na <strong>in</strong>tervalu [a, b] zvezna, potem znotraj<br />

<strong>in</strong>tervala obstaja vsaj eno število ξ, pri katerem velja:<br />

∫ b<br />

a<br />

x(t) dt = (b − a)x(ξ) . (2.7-28)<br />

Slika 2.7-4<br />

Povprečna vrednost signala je<br />

enaka viš<strong>in</strong>i pravokotnika z<br />

dolž<strong>in</strong>o (a, b), ki ima enako<br />

plošč<strong>in</strong>o kot jo določa <strong>in</strong>tegral<br />

Geometrijska <strong>in</strong>terpretacija je naslednja: med točkama a <strong>in</strong> b<br />

obstaja točka ξ, v kateri je vrednost funkcije x takšna, da je njen<br />

zmnožek z def<strong>in</strong>icijskim <strong>in</strong>tervalom enak plošč<strong>in</strong>i, ki jo določa<br />

x(t) nad tem <strong>in</strong>tervalom (slika 2.7-4). Z drugimi besedami, v tej<br />

točki signal <strong>za</strong>v<strong>za</strong>me svojo srednjo vrednost.<br />

2.7-28. a b<br />

x( t)<br />

x( )<br />

z<br />

b<br />

x ( t ) dt <br />

a<br />

( b a) x( )<br />

t<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.7 Norme <strong>signalov</strong> 33<br />

2.7.4 Drugi izrek o povprečni vrednosti<br />

To poglavje je namenjeno <strong>za</strong>htevnejšim študentom. Za nadaljno<br />

osnovno obravnavo teorije <strong>signalov</strong> je poznavanje tega poglavja<br />

koristno, ni pa nujno.<br />

□ □ □ □ <br />

<strong>za</strong>htevna snov<br />

V izpeljavah <strong>in</strong> dokazih, s katerimi se bomo seznanili kasneje, bomo<br />

mnogokrat uporabili drugi ali tudi posplošeni izrek o povprečni vrednosti<br />

[3]. Glasi se:<br />

∫ b<br />

a<br />

x(t)y(t) dt = x(a + 0)<br />

∫ ξ<br />

a<br />

y(t) dt + x(b − 0)<br />

∫ b<br />

ξ<br />

y(t) dt (2.7-29)<br />

pri a ξ b. x(a + 0) je desna limita funkcije x(t) v točki a, x(b − 0)<br />

pa leva limita v točki b:<br />

x(a + 0) = lim<br />

∆t→0<br />

∆t>0<br />

x(a + ∆t) , x(b − 0) = lim<br />

∆t→0<br />

∆t>0<br />

x(b − ∆t) .<br />

Obrazec bomo kasneje še <strong>za</strong>pisali v obliki izreka, kjer bomo natančno<br />

podali pogoje veljavnosti izreka. Drugi izrek o povprečni vrednosti<br />

ima poleg že <strong>za</strong>pisane splošne oblike (2.7-29) še dve posebni<br />

obliki, ki veljata pri ožjih pogojih izbire funkcije x(t). Z njima bomo<br />

doka<strong>za</strong>li veljavnost (2.7-29).<br />

IZREK 2<br />

Če je funkcija x(t) na <strong>in</strong>tervalu (a, b) omejena, monotono padajoča <strong>in</strong><br />

nenegativna, funkcija y(t) pa na istem <strong>in</strong>tervalu omejena, <strong>in</strong>tegrabilna<br />

<strong>in</strong> ima končno število menjav predznaka, potem velja:<br />

∫ b<br />

a<br />

x(t)y(t) dt = x(a)<br />

∫ ξ<br />

a<br />

y(t) dt , a ξ b . (2.7-30)<br />

Pojasnimo izrek. Ker ima funkcija y(t) lahko le končno število<br />

menjav predznaka (torej prehodov skozi absciso), lahko <strong>in</strong>terval (a, b)<br />

razmejimo s končnim številom točk a r , r = 0, . . . , p + 1:<br />

a = a 0 , a 1 , a 2 , . . . a p , a p+1 = b<br />

tako, da funkcija y(t) znotraj <strong>in</strong>tervala (a r , a r+1 ) ne spremeni predznaka.<br />

Upoštevaje to razdelitev, lahko <strong>za</strong>pišemo:<br />

∫ b<br />

a<br />

x(t)y(t) dt =<br />

p∑<br />

r=0<br />

∫ ar+1<br />

a r<br />

x(t)y(t) dt .


34 2. Signali<br />

Na osnovi prvega izreka o povprečni vrednosti (glej teorem 1 na strani<br />

32) lahko sklepamo, da obstaja število µ r , x(a r ) µ r x(a r+1 ), pri<br />

katerem velja:<br />

∫ ar+1<br />

a r<br />

x(t)y(t) dt = µ r<br />

∫ ar+1<br />

a r<br />

y(t) dt ,<br />

Integral na desni izračunamo:<br />

∫ a ar+1<br />

r+1<br />

y(t) dt = F (t)<br />

= [F (a<br />

a r<br />

∣<br />

r+1 ) − F (a r )] ,<br />

a r<br />

kjer je F (t) = ∫ t<br />

y(t) dt. Zato <strong>za</strong>nj velja:<br />

a<br />

F (a) =<br />

∫ a<br />

a<br />

y(t) dt = 0 , F (a r ) =<br />

Sedaj lahko (2.7-30) preuredimo v:<br />

∫ b<br />

a<br />

∫ ar<br />

a<br />

y(t) dt , F (a r+1 ) =<br />

x(t)y(t) dt = µ 0 [F (a 1 ) − F (a 0 )] + µ 1 [F (a 2 ) − F (a 1 )]+<br />

∫ ar+1<br />

a<br />

y(t) dt .<br />

· · · + µ r−1 [F (a r ) − F (a r−1 )] + µ r [F (a r+1 ) − F (a r )] + · · ·<br />

= −µ 0 F (a 0 ) + (µ 0 − µ 1 )F (a 1 ) + (µ 1 − µ 2 )F (a 2 )+<br />

· · · + (µ p−1 − µ p )F (a r ) + µ p F (a p+1 ) .<br />

Ker je F (a) = 0, lahko prejšnji vsoti prištejemo x(a)F (a), ne da bi jo<br />

s tem spremenili:<br />

∫ b<br />

a<br />

x(t)y(t) dt = [x(a) − µ 0 ]F (a 0 ) + [µ 0 − µ 1 ]F (a 1 ) + [µ 1 − µ 2 ]F (a 2 )+<br />

· · · + [µ r−1 − µ r ]F (a r ) + µ r F (a r+1 ) + · · · .<br />

(2.7-31)<br />

V naslednjem koraku upoštevamo, da so <strong>za</strong>radi <strong>in</strong>tegrabilnosti y vrednosti<br />

F (a r ) omejene. Zato obstaja tak M, da velja:<br />

m<strong>in</strong> F (a r ) < M < max F (a r ) .<br />

a r<br />

a r<br />

Hkrati lahko desno stran v (2.7-31) nadomestimo z:<br />

[<br />

M (x(a) − µ 0 ) + (µ 0 − µ 1 ) + (µ 1 − µ 2 ) + · · ·<br />

+ (µ r−1 − µ r ) + µ r<br />

]<br />

= Mx(a) .<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.7 Norme <strong>signalov</strong> 35<br />

Ker je funkcija F (t) <strong>za</strong>radi omejenosti y zvezna na <strong>in</strong>tervalu (a, b),<br />

obstaja na <strong>in</strong>tervalu (a, b) točka ξ, a ξ b, pri kateri velja:<br />

ter <strong>za</strong>to tudi:<br />

∫ b<br />

a<br />

M = F (ξ) =<br />

∫ ξ<br />

a<br />

x(t)y(t) dt = x(a)<br />

y(t) dt<br />

∫ ξ<br />

a<br />

y(t) dt .<br />

Drugi posebni primer predpostavlja, da je x(t) omejena, monotono<br />

naraščajoča funkcija, torej:<br />

IZREK 3<br />

Če je funkcija x(t) na <strong>in</strong>tervalu (a, b) omejena, monotono naraščajoča<br />

<strong>in</strong> nenegativna, funkcija y(t) pa ima enake lastnosti kot v izreku 2,<br />

potem velja:<br />

∫ b<br />

kjer <strong>za</strong> ξ velja a ξ b.<br />

a<br />

x(t)y(t) dt = x(b)<br />

∫ b<br />

ξ<br />

y(t) dt ,<br />

Dokaz drugega posebnega primera lahko izpeljemo po podobni<br />

poti, kot smo doka<strong>za</strong>li prvi posebni primer. Zato se napotimo kar k<br />

splošnemu izreku. Ta manj omejuje potek funkcije x(t), glasi pa se:<br />

IZREK 4 (Drugi izrek o povprečni vrednosti)<br />

Če je funkcija x(t) na <strong>in</strong>tervalu (a, b) omejena <strong>in</strong> monotona, funkcija<br />

y(t) pa ima enake lastnosti kot v izreku 2, potem velja:<br />

∫ b<br />

a<br />

y(t)x(t) dt = x(a)<br />

kjer <strong>za</strong> ξ velja a ξ b.<br />

∫ ξ<br />

a<br />

y(t) dt + x(b)<br />

∫ b<br />

ξ<br />

y(t) dt ,<br />

DOKAZ Dokažimo veljavnost <strong>za</strong> monotono padajočo funkcijo x(t).<br />

Upoštevajmo dokaz izreka 2. Tam smo poka<strong>za</strong>li, da <strong>za</strong> monotono padajočo<br />

funkcijo f(t) velja:<br />

∫ b<br />

f(t)y(t) dt = f(a)<br />

∫ ξ<br />

a<br />

a<br />

y(t) dt . (2.7-32)<br />

Če izberemo f(t) tako, da velja f(t) = x(t) − x(b), potem je tudi<br />

f(a) = x(a) − x(b) <strong>in</strong> x(t) = f(t) + x(b), <strong>za</strong>to lahko izpeljemo:<br />

∫ b<br />

a<br />

x(t)y(t) dt =<br />

=<br />

∫ b<br />

a<br />

∫ b<br />

[f(t) + x(b)] y(t) dt<br />

f(t)y(t) dt + x(b)<br />

∫ b<br />

a<br />

a<br />

y(t) dt .<br />

(2.7-33)


36 2. Signali<br />

Za prvi <strong>in</strong>tegral na desni strani (2.7-33) upoštevamo (2.7-32), <strong>za</strong>to:<br />

∫ b<br />

x(t)y(t) dt = f(a)<br />

∫ ξ<br />

a<br />

a<br />

= [x(a) − x(b)]<br />

y(t) dt + x(b)<br />

∫ ξ<br />

∫ b<br />

a<br />

y(t) dt<br />

y(t) dt + x(b)<br />

∫ b<br />

a<br />

a<br />

y(t) dt .<br />

Upoštevajmo še lastnost ∫ c<br />

a = ∫ b<br />

a + ∫ c<br />

, a b c ter drugi <strong>in</strong>tegral na<br />

b<br />

desni razcepimo:<br />

∫ b<br />

a<br />

x(t)y(t) dt = [x(a) − x(b)]<br />

= x(a)<br />

kar je iskani rezultat.<br />

+ x(b)<br />

∫ ξ<br />

∫ ξ<br />

a<br />

∫ ξ<br />

a<br />

y(t) dt<br />

y(t) dt + x(b)<br />

y(t) dt + x(b)<br />

∫ b<br />

a<br />

ξ<br />

∫ b<br />

ξ<br />

y(t) dt ,<br />

y(t) dt<br />

(2.7-34)<br />

□<br />

□ □ □ □ <br />

Do rezultata (2.7-34) pridemo tudi ob predpostavki, da je f(t)<br />

monotono naraščajoča funkcija. V tem primeru bi upoštevali drugi posebni<br />

primer. Rezultat (2.7-34) je pomemben, ker ga uporabljamo v<br />

mnogih dokazih. Mi ga bomo uporabili pri dokazovanju veljavnosti Dirichletovega<br />

<strong>in</strong>tegrala oziroma Fourierovega <strong>in</strong>tegralskega izreka v drugi<br />

knjigi Teorije <strong>signalov</strong>.<br />

2.8 Skalarni produkt<br />

Zadnja snov, ki jo bomo obdelali v okviru signalnega prostora, je<br />

skalarni produkt. Za njegovo razumevanje si ta pojem razložimo<br />

s primerom skalarnega produkta dveh vektorjev, ki ležita v ravn<strong>in</strong>i.<br />

S φ x <strong>in</strong> φ y označimo kota med vektorjema <strong>in</strong> absciso pravokotnega<br />

koord<strong>in</strong>atnega sistema (slika 2.8-1). Kos<strong>in</strong>us kota φ x −φ y ,<br />

ki je med vektorjema x <strong>in</strong> y, je:<br />

cos(φ x − φ y ) = cos φ x cos φ y + s<strong>in</strong> φ x s<strong>in</strong> φ y<br />

x 2 y<br />

= √ · √ 2<br />

+<br />

x<br />

2<br />

1 + x 2 2 y<br />

2<br />

1 + y2<br />

2<br />

=<br />

x 1 y 1 + x 2 y<br />

√ √ 2<br />

x<br />

2<br />

1 + x 2 2 y<br />

2<br />

1 + y2<br />

2<br />

x<br />

√ 1<br />

·<br />

x<br />

2<br />

1 + x 2 2<br />

y 1<br />

√<br />

y<br />

2<br />

1 + y 2 2<br />

, (2.8-1)<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.8 Skalarni produkt 37<br />

x 1<br />

y 1<br />

y 2<br />

x<br />

y<br />

x<br />

y<br />

Slika 2.8-1<br />

Skalarni produkt<br />

j<br />

i<br />

x 2<br />

kjer so x k <strong>in</strong> y k komponente vektorjev x <strong>in</strong> y. Imenovalec v<br />

izpeljanem izrazu (2.8-1) je produkt dolž<strong>in</strong> obeh vektorjev, števec<br />

pa je znan skalarni produkt vektorjev x <strong>in</strong> y. Pri tem vemo, da<br />

lahko dolž<strong>in</strong>o vektorja določimo z normo p = 2. Torej lahko<br />

(2.8-1) <strong>za</strong>pišemo kot:<br />

cos(φ x − φ y ) = x 1y 1 + x 2 y 2<br />

‖x‖ 2 · ‖y‖ 2<br />

(2.8-2)<br />

<strong>in</strong> izračunamo skalarni produkt vektorjev x <strong>in</strong> y:<br />

x 1 y 1 + x 2 y 2 = ‖x‖ 2 · ‖y‖ 2 · cos(φ x − φ y ) (2.8-3)<br />

Za dani dolž<strong>in</strong>i vektorjev x <strong>in</strong> y je skalarni produkt mera medsebojne<br />

lege dveh vektorjev. Če je skalarni produkt 0, mora biti<br />

cos(φ x − φ y ) = 0 (saj so dolž<strong>in</strong>e neničelnih vektorjev vedno večje<br />

od nič). Torej je kot med vektorjema x <strong>in</strong> y enak π/2 oziroma<br />

90 o <strong>in</strong> sta vektorja pravokotna drug na drugega. Skalarni produkt<br />

je maksimalni, ko je je cos(φ x −φ y ) = 1. Takrat sta vektorja<br />

vzporedna oziroma je kot med njima enak 0.<br />

2.8.1 Skalarni produkt v prostoru l<br />

Skalarni produkt dveh vektorjev x = (x 1 x 2 x 3 ) <strong>in</strong> y = (y 1 y 2 y 3 )<br />

v tridimenzionalnem prostoru je x 1 y 1 + x 2 y 2 + x 3 y 3 <strong>in</strong> ga lahko<br />

geometrijsko predstavimo podobno kot skalarni produkt dvodimenzionalnih<br />

vektorjev. Podobno lahko določimo skalarni produkt<br />

N-dimenzionalnih vektorjev x <strong>in</strong> y:<br />

〈x, y〉 = x 1 y 1 + x 2 y 2 + · · · x N y N =<br />

N∑<br />

x n y n . (2.8-4)<br />

n=1


38 2. Signali<br />

Veljavnost <strong>za</strong>pisa (2.8-4) ni omejena le na geometrijsko predstavljive<br />

vektorje, ampak velja tudi <strong>za</strong> neskončno dimenzionalne vektorje.<br />

V primeru časovnih <strong>za</strong>poredij x(n) <strong>in</strong> y(n), ki pripadata<br />

prostoru l, njun skalarni produkt izpeljemo iz (2.8-4):<br />

〈x(n), y(n)〉 = ∑ n∈Z<br />

x n y ∗ n , (2.8-5)<br />

kjer je y ∗ n konjugirano kompleksna vrednost y n . Skalarni produkt<br />

<strong>za</strong>poredja x(n) s seboj določa energijo signala, torej kvadrat njegove<br />

Evklidske norme:<br />

〈x(n), x(n)〉 = ∑ n∈Z<br />

x n x ∗ n = ‖x‖ 2 2 . (2.8-6)<br />

Ker lahko imata dva poljubna elementa prostora l, na primer<br />

neskončni <strong>za</strong>poredji x = y = (· · · , 1, 1, 1, · · · ), nedoločljivi (neskončni)<br />

skalarni produkt, prostor l na splošno ni prostor skalarnega<br />

produkta. Iz Cauchy-Schwartzovega izreka (njegova razlaga<br />

sledi) izhaja, da lahko le elementi prostora l, ki imajo končno<br />

energijo, (paroma) tvorijo skalarni produkt. Ti elementi spadajo<br />

v podprostor l 2 , ki je prostor s skalarnim produktom. To ugotovitev<br />

poudarimo z naslednjo def<strong>in</strong>icijo:<br />

DEFINICIJA 11 (Prostor skalarnega produkta l 2 )<br />

Vse množice <strong>za</strong>poredij iz l s končno energijo tvorijo prostor skalarnega<br />

produkta l 2 .<br />

Ta prostor imenujemo tudi Hilbertov prostor.<br />

2.8.2 Skalarni produkt v prostoru L<br />

Kot neskončno dimenzionalne vektorje si lahko predstavljamo<br />

tudi funkcije, na primer x(t) <strong>in</strong> y(t), z def<strong>in</strong>icijsko domeno na<br />

<strong>in</strong>tervalu (a, b) ∈ R. Če predpostavimo, da vsaka točka def<strong>in</strong>icijskega<br />

območja določa eno dimenzijo vektorja, vrednost funkcije<br />

v tej točki pa komponento vektorja, sta funkciji ”neskončno dimenzionalna<br />

vektorja”<strong>in</strong> je njun skalarni produkt določen z:<br />

〈x, y〉 = lim<br />

n→∞<br />

∆→0<br />

b∑<br />

k=a<br />

k=k+n∆<br />

x k y ∗ k = ∫ b<br />

a<br />

x(t)yk ∗ (t) dt . (2.8-7)<br />

Rezultat (2.8-5) je enak nič, če sta si ”vektorja”, torej funkciji<br />

x(t) <strong>in</strong> y(t) na <strong>in</strong>tervalu (a, b) medsebojno pravokotni oziroma<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.8 Skalarni produkt 39<br />

s tujko ortogonalni. Ortogonalne funkcije so v obdelavi <strong>signalov</strong><br />

zelo pomembne, <strong>za</strong>to jim je posvečen v tej knjigi poseben<br />

razdelek. Za sedaj se <strong>za</strong>dovoljimo z napovedjo, da je razred ortogonalnih<br />

funkcij zelo pomemben pri obravnavi <strong>signalov</strong> <strong>in</strong> jih<br />

bomo <strong>za</strong>to kasneje še podrobneje obravnavali.<br />

Podobno kot v prejšnjem podrazdelku, lahko tudi tu def<strong>in</strong>iramo<br />

skalarni produkt funkcije same s seboj:<br />

〈x, y〉 =<br />

∫ b<br />

a<br />

x(t)x ∗ (t) dt =<br />

∫ b<br />

a<br />

|x(t)| 2 dt . (2.8-8)<br />

Tudi tukaj hitro ugotovimo, da vsi elementi L nimajo končne<br />

norme. Zato elementi s končno energijo, ki lahko medsebojno tvorijo<br />

skalarni produkt, določajo podprostor skalarnega produkta<br />

L. Poudarimo to z def<strong>in</strong>icijo:<br />

DEFINICIJA 12 (Prostor skalarnega produkta L 2 )<br />

Vse funkcije iz L s končno energijo tvorijo prostor s skalarnim<br />

produktom L 2 .<br />

2.8.3 Cauchy-Schwartzova neenakost<br />

Za vsak par elementov, na primer x <strong>in</strong> y, prostora s skalarnim<br />

produktom velja relacija:<br />

|〈x, y〉| 2 〈x, x〉〈y, y〉 . (2.8-9)<br />

Enakost velja takrat <strong>in</strong> samo takrat, če obstaja tak skalar α, da<br />

je:<br />

x = αy ali y = αx . (2.8-10)<br />

August<strong>in</strong> Louis Cauchy (1789 – 1857) je bil po izobrazbi<br />

gradbeni <strong>in</strong>ženir. Danes je znan predvsem po<br />

svojih številnih prispevkih v razvoju matematike.<br />

Amandus Schwartz (1843 – 1921) je bil nemški matematik.<br />

Dokaz Cauchy-Schwartzove neenakosti presega namen tega<br />

učbenika, <strong>za</strong>to si oglejmo njegovo razlago v prostoru realnih števil,


40 2. Signali<br />

kjer je skalarni produkt x · x enak kvadratu Evklidske norme:<br />

‖x‖ 2 2 . V tem primeru lahko Cauchy-Schwartzovo neenakost preoblikujemo<br />

v:<br />

〈x, y〉<br />

−1 1 .<br />

‖x‖ 2 ‖y‖ 2<br />

Spomnimo se, da srednji člen v zgornji enačbi pomeni ”kos<strong>in</strong>us<br />

kota med x <strong>in</strong> y”. Iz tega sledi, da se x <strong>in</strong> y maksimalno ujemata<br />

(v isti ali nasprotni smeri), ko je njun skalarni produkt enak<br />

±‖x‖ 2 ‖y‖ 2 , <strong>in</strong> maksimalno neujemata (razlikujeta), ko je enak 0.<br />

2.8.4 Signalni prostor<br />

Iz razlage v <strong>za</strong>dnjih treh podpoglavjih lahko <strong>za</strong>ključimo, da je<br />

v teoriji <strong>signalov</strong> obravnava <strong>signalov</strong> v l<strong>in</strong>earnem prostoru, kjer<br />

je def<strong>in</strong>iran skalarni produkt, izjemnega pomena. Zato l<strong>in</strong>earni<br />

prostor tipa l 2 ali L 2 , v katerem so le signali s končno energijo<br />

<strong>in</strong> pri katerih lahko določimo skalarni produkt, imenujemo tudi<br />

signalni prostor. K term<strong>in</strong>u signalnega prostora se bomo še<br />

vrnili v <strong>za</strong>ključku naslednjega poglavja.<br />

Rezultat operacije je<br />

korelacijska funkcija<br />

2.9 Podobnost <strong>signalov</strong><br />

Pri opisu <strong>signalov</strong> nas mnogokrat <strong>za</strong>nima podobnost oziroma razlika<br />

med dvema signaloma ali zbirkama podatkov ali procesoma.<br />

Primerjavo med njimi lahko izvedemo s korelacijo.<br />

Korelacijo uporabljamo v mnogih področjih znanosti, tehnike,<br />

ekonomije itd. Na primer pri obdelavi slike v robotskem<br />

(računalniškem) vidu, dalj<strong>in</strong>skem opazovanju s satelitov, kjer primerjamo<br />

podatke med različnimi slikami; potem v radarskih <strong>in</strong><br />

sonarskih sistemih, kjer s korelacijo oddajanega <strong>in</strong> odbitega signala<br />

določamo oddaljenost <strong>in</strong> smer opazovanih objektov; pri detekciji<br />

<strong>in</strong> <strong>in</strong>dentifikacji <strong>signalov</strong> v šumu, v regulacijah, kjer lahko<br />

s korelacijo vhodnega <strong>in</strong> izhodnega signala opazujemo uč<strong>in</strong>ek sistema<br />

na signal. Uporabljamo jo tudi v <strong>in</strong>dentifikacijah sprejetih<br />

kodnih besed pri pulzno kodni modulaciji itd.<br />

Korelacijo bomo označevali s simbolom r. Ločimo:<br />

križno korelacijo: r 12 , ki primerja dva različna signala,<br />

avtokorelacijo: r 11 , ki primerja signal s samim seboj.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.9 Podobnost <strong>signalov</strong> 41<br />

2.9.1 Koncept<br />

Pojem korelacije vpeljimo z naslednjim premislekom. Opazujmo<br />

dve <strong>za</strong>poredji podatkov x 1 (n) <strong>in</strong> x 2 (n), ki sta nastajali sočasno<br />

<strong>in</strong> <strong>za</strong>jemali svoje vrednosti iz R. Njihovo medsebojno podobnost<br />

lahko merimo z vsoto produktov posameznih podatkov:<br />

r 12 = ∑ n∈N<br />

x 1 (n)x 2 (n) ,<br />

kjer je N število podatkov v <strong>za</strong>poredju. Smiselnost te mere pride<br />

do izra<strong>za</strong>, če podatki v x 1 (n) <strong>in</strong> x 2 (n) naključni <strong>in</strong> medsebojno<br />

neodvisni. V tem primeru se vsota produktov z naraščanjem<br />

števila N manjša proti nič (<strong>in</strong> s tem kaže na naključnost <strong>in</strong> neodvisnost<br />

podatkov). Obstoj pozitivne vsote kaže na podobnost<br />

med podatki, obstoj negativne vsote pa, da so si podatki podobni,<br />

vendar prevladujejo nasprotni predznaki.<br />

Zapisana def<strong>in</strong>icija korelacije da rezultat, ki je odvisen tudi<br />

od števila podatkov. To anomalijo lahko odstranimo z normaliziranjem<br />

rezultata - podelimo ga s številom podatkov:<br />

r 12 = 1 N<br />

∑<br />

x 1 (n)x 2 (n) ,<br />

n∈N<br />

S to normali<strong>za</strong>cijo smo sicer rešili eno od težav pri def<strong>in</strong>iranju<br />

korelacije, nakopali smo si pa novo. V primeru, ko N narašča<br />

čez vse meje, R 12 ne glede na medsebojno podobnost x 1 (n) <strong>in</strong><br />

x 2 (n) upada proti nič. Poleg teh težav ta def<strong>in</strong>icija ne odkrije<br />

podobnost <strong>signalov</strong>, ki jih kaže slika 2.9-1. Vidimo, da <strong>za</strong> ta<br />

signala velja x 1 (n) = x 2 (n − m), vendar je rezultat korelacije<br />

x n 2( )<br />

x n 1( )<br />

m<br />

n<br />

Slika 2.9-1<br />

n


42 2. Signali<br />

računan po zgornjih formulah enak nič. Maksimalnega pa bi<br />

dobili v primeru, če bi zgornji obrazec dopolnili tako:<br />

r 12 (m) = 1 ∑<br />

x 1 (n)x 2 (n − m)<br />

N<br />

n∈N<br />

= r 12 (−m) = 1 N<br />

∑<br />

x 1 (n)x 2 (n + m) ,<br />

n∈N<br />

kjer m predstavlja premik enega <strong>za</strong>poredja napram drugemu.<br />

Zamiki med signali ponavadi nastanejo <strong>za</strong>radi <strong>za</strong>kasnitev, ki<br />

nastanejo pri prehodu <strong>signalov</strong> skozi nek sistem. Očitno je, da<br />

lahko s sprem<strong>in</strong>janjem m določimo <strong>za</strong>kasnitev med dvema signaloma.<br />

To najdemo takrat, ko ima r 12 (m) maksimum.<br />

Seveda s tem še ni konec vseh težav z računanjem korelacije.<br />

Na primer robni efekt, pa odvisnost korelacije od velikosti<br />

amplitudnega območja. Pov<strong>za</strong>memo lahko, da moramo obrazec<br />

<strong>za</strong> izračun korelacije prilagoditi naravi signala. Zato bomo v nadaljevanju<br />

podali več def<strong>in</strong>icij korelacij, ki jih uporabljamo <strong>za</strong><br />

posamezne vrste <strong>signalov</strong>.<br />

2.9.2 Korelacija pri aperiodičnih signalih<br />

Pri signalih s končno energijo, na primer pri aperiodičnih pulzih,<br />

korelacije ne moremo def<strong>in</strong>irati kot povprečno vrednost <strong>za</strong>maknjenih<br />

produktov preko časovnega <strong>in</strong>tervala T , ker s T → ∞<br />

korelacija teži k nič. Zato <strong>za</strong> te signale uporabljamo <strong>za</strong> izračun<br />

korelacije naslednja obrazca:<br />

Časovno zvezni signali:<br />

Imamo realna signala x 1 (t) <strong>in</strong> x 2 (t), t ∈ R.<br />

Križna korelacija med njima je def<strong>in</strong>irana<br />

z naslednjim:<br />

r 12 (τ) =<br />

∫ ∞<br />

−∞<br />

x 1 (t)x 2 (t + τ) dt , (2.9-1)<br />

kjer je τ ∈ R <strong>za</strong>mik med signaloma. Križna<br />

korelacija obstaja, če sta funkciji x 1 (t) <strong>in</strong><br />

x 2 (t) absolutno <strong>in</strong>tegrabilni:<br />

∫ ∞<br />

−∞<br />

|x i (t)| dt < ∞ , i = 1, 2 . (2.9-2)<br />

Časovno diskretni signali:<br />

Imamo realna signala x 1 (n) <strong>in</strong> x 2 (n), n ∈<br />

Z. Križna korelacija med njima je def<strong>in</strong>irana<br />

z vsoto:<br />

r 12 (m) = ∑ n∈N<br />

x 1 (n)x 2 (n + m) , (2.9-3)<br />

kjer je m ∈ Z <strong>za</strong>mik med signaloma.<br />

Križna korelacija obstaja, če sta funkciji<br />

x 1 (n) <strong>in</strong> x 2 (n) absolutno konvergentni:<br />

∞∑<br />

n=−∞<br />

|x i (n)| < ∞ , i = 1, 2 . (2.9-4)<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.9 Podobnost <strong>signalov</strong> 43<br />

Pri avtokorelaciji signal primerjamo s samim seboj: x 1 = x 2 :<br />

Časovno zvezni signali:<br />

Časovno diskretni signali:<br />

r 11 (τ) =<br />

∫ ∞<br />

−∞<br />

x 1 (t)x 1 (t + τ) dt (2.9-5)<br />

r 11 (m) =<br />

∞∑<br />

n=−∞<br />

Grafično predstavitev računanja križne korelacije dveh aperiodičnih<br />

zveznih <strong>signalov</strong> vidimo na sliki 2.9-1. S slike vidimo, da<br />

se korelacija r 12 (τ) sprem<strong>in</strong>ja sorazmerno s površ<strong>in</strong>o, ki jo določa<br />

produkt obeh <strong>signalov</strong> pri njunem danem medsebojnem <strong>za</strong>miku.<br />

x 1 (n)x 1 (n + m) (2.9-6)<br />

1<br />

x t<br />

1 ( )<br />

1<br />

x ( 2<br />

t )<br />

0<br />

t<br />

0<br />

t<br />

1<br />

<br />

0<br />

t<br />

1<br />

<br />

1<br />

<br />

0<br />

t<br />

Slika 2.9-2<br />

Postopek računanja križne<br />

korelacije aperiodičnih<br />

<strong>signalov</strong><br />

0<br />

t<br />

1<br />

<br />

1<br />

12( )<br />

0<br />

t<br />

0


44 2. Signali<br />

2.9.3 Korelacija pri periodičnih signalih<br />

Iz lastnosti periodičnih <strong>signalov</strong> sklepamo, da se rezultat korelacije<br />

ponavlja s periodo <strong>signalov</strong>. Zato je smiselno pri računanju<br />

korelacije signale <strong>za</strong>mikati le na območju ene periode signala.<br />

Drugi problem, ki nastane pri periodičnih signalih, je absolutna<br />

<strong>in</strong>tegrabilnost. Hitro lahko uvidimo, da periodični signali niso<br />

absolutno <strong>in</strong>tegrabilni, saj imajo neskončno energijo. Def<strong>in</strong>icija<br />

korelacije pri periodičnih signalih se <strong>za</strong>to razlikuje od njene def<strong>in</strong>icije<br />

pri aperiodičnih signalih.<br />

Časovno zvezni signali:<br />

Križna korelacija med realnima periodičnima<br />

časovno zveznima signaloma x 1 (t) <strong>in</strong><br />

x 2 (t), t ∈ R z isto periodo T je def<strong>in</strong>irana<br />

z <strong>in</strong>tegralom:<br />

r 12 (τ) = 1 ∫<br />

x 1 (t)x 2 (t + τ) dt . (2.9-7)<br />

T<br />

T<br />

kjer je τ ∈ R <strong>za</strong>mik med signaloma. Za<br />

avtokorelacijo pa velja:<br />

r 12 (τ) = 1 ∫<br />

x 1 (t)x 1 (t + τ) dt . (2.9-8)<br />

T<br />

T<br />

Časovno diskretni signali:<br />

Križna korelacija med realnima periodičnima<br />

časovno diskretnima signaloma x 1 (n) <strong>in</strong><br />

x 2 (n), n ∈ Z z enako periodo <strong>in</strong> številom<br />

diskretnih vrednosti je def<strong>in</strong>irana z vsoto:<br />

r 12 (m) = 1 ∑<br />

x 1 (n)x 2 (n + m) , (2.9-9)<br />

N<br />

N<br />

kjer je m ∈ Z <strong>za</strong>mik med signaloma.<br />

avtokorelacijo pa velja:<br />

Za<br />

r 12 (m) = 1 ∑<br />

x 1 (t)x 1 (n + m) . (2.9-10)<br />

N<br />

N<br />

Iz def<strong>in</strong>icije sledi, da je križna korelacija povprečna vrednost<br />

<strong>in</strong>tegrala oziroma vsote produkta obeh <strong>signalov</strong> pri medsebojnem<br />

premiku τ oziroma m. Na kraju še poudarimo:<br />

Korelacijo lahko v fizikalnem pomenu razumemo kot vplivanje med<br />

signali oziroma kot pretok energije med signali. Če pretoka energije<br />

med njimi ni, še ne pomeni, da sta signala medsebojno neodvisna.<br />

2.9.4 Korelacija pri kompleksnih aperiodičnih signalih<br />

V obrazcih (2.9-1), (2.9-3) <strong>in</strong> (2.9-5), (2.9-6) smo <strong>za</strong>pisali korelacijo<br />

<strong>za</strong> realne signale. Če so signali kompleksni, moramo drugi<br />

signal upoštevati v konjugirano kompleksni obliki:<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.9 Podobnost <strong>signalov</strong> 45<br />

časovno zvezni signali:<br />

Časovno diskretni signali:<br />

r 12 (τ) =<br />

r 11 (τ) =<br />

∫ ∞<br />

−∞<br />

∫ ∞<br />

−∞<br />

x 1 (t)x ∗ 2(t + τ) dt (2.9-11)<br />

x 1 (t)x ∗ 1(t + τ) dt (2.9-12)<br />

∞∑<br />

r 12 (m) = x 1 (n)x ∗ 2(n + m) (2.9-13)<br />

n=−∞<br />

∞∑<br />

r 11 (m) = x 1 (n)x ∗ 1(n + m) (2.9-14)<br />

n=−∞<br />

2.9.5 Izračun energije <strong>in</strong> moči s korelacijo<br />

Iz def<strong>in</strong>icije korelacije pri aperiodičnih signalih sledi, da pri ničtem<br />

<strong>za</strong>miku <strong>signalov</strong> s korelacijo izračunamo energijo E, ki jo pri<br />

križni korelaciji določata dva signala, na primer x <strong>in</strong> y, pri avtokorelaciji<br />

pa energijo, ki jo vsebuje obravnavani signal:<br />

Časovno zvezni signali:<br />

Časovno diskretni signali:<br />

r 11 (0) =<br />

∫ ∞<br />

−∞<br />

x 1 (t)x ∗<br />

} {{ 1(t)<br />

}<br />

dt = E .<br />

|x 1 (t)| 2<br />

∞∑<br />

r 11 (0) =<br />

n=−∞<br />

x 1 (n)x ∗<br />

} {{ 1(n)<br />

}<br />

= E .<br />

|x 1 (n)| 2<br />

Pri periodičnih signalih <strong>in</strong> signalih z omejeno močjo pa povprečno<br />

moč signala:<br />

Časovno zvezni signali:<br />

Časovno diskretni signali:<br />

r 11 (0) = 1 T<br />

∫ ∞<br />

−∞<br />

x 1 (t)x ∗ 1(t) dt = P .<br />

r 11 (0) = 1 N<br />

∞∑<br />

n=−∞<br />

x 1 (n)x ∗ 1(n) = P .<br />

2.9.6 Robni efekt<br />

Robni efekt nastane pri medsebojnem premiku <strong>signalov</strong> omejene<br />

dolž<strong>in</strong>e. Na primer, pri <strong>za</strong>poredju diskretnih vrednosti <strong>za</strong>poredja<br />

podatkov x 1 (n) <strong>in</strong> x 2 (2), podatki nimajo več svojega primerjalnega<br />

para <strong>in</strong> so <strong>za</strong>to izpuščeni iz računanja korelacije. Posledica<br />

je, da r 12 (m) upada z naraščanjem medsebojnega premika m.


46 2. Signali<br />

Seveda je tako izračunana korelacija vprašljiva.<br />

Ena izmed možnih rešitev je podvojiti <strong>za</strong>htevano dolž<strong>in</strong>o enega<br />

od <strong>za</strong>poredij. To lahko dosežemo z <strong>za</strong>pisom večjega števila podatkov,<br />

ali pri periodičnem signalu s ponavljanjem <strong>za</strong>pisa.<br />

Druga rešitev je odstraniti robni efekt iz izračunanega rezultata.<br />

Popravek ali korekcijo naredimo na osnovi naslednjega<br />

premisleka. Predpostavimo, da imata <strong>za</strong>poredji x 1 (n) <strong>in</strong> x 2 (2)<br />

enako dolž<strong>in</strong>o ter vse svoje vrednosti enake. V tem primeru dobimo<br />

<strong>za</strong>radi robnega efekta l<strong>in</strong>earno upadanje r 12 (m) od vrednosti<br />

r 12 (m) = r 12 (0), ko upoštevamo vse pare podatkov do vrednosti<br />

r 12 (m) = r 12 (N) = 0, ko se <strong>za</strong>poredje več ne prekrivata<br />

(slika 2.9-3). Vmes, pri premiku <strong>za</strong> m, 0 < m < N, je prava<br />

r 12 (0)<br />

r 12 ( )<br />

r 12 ( )<br />

pravi<br />

Slika 2.9-3<br />

Robni efekt<br />

r 12 ( )<br />

t<br />

t<br />

vrednost korelacije r 12 (m) pravi , izračunana pa je <strong>za</strong>radi robnega<br />

efekta enaka r 12 (m). Glede na sliko 2.9-3 sledi:<br />

oziroma<br />

r 12 (m) pravi − r 12 (m)<br />

m<br />

= r 12(0)<br />

N<br />

r 12 (m) pravi = r 12 (m) + m N r 12(0) .<br />

,<br />

Sledi, da pravo vrednost korelacije preprosto dobimo tako, da<br />

izračunani vrednosti prištejemo r 12 (0)m/N.<br />

Pri periodičnih signalih ponavadi uporabimo prvo rešitev, to<br />

je s podaljšanjem <strong>za</strong>pisa podatkov, saj se ti periodično ponavljajo<br />

<strong>in</strong> pri podaljševanju <strong>za</strong>pisa s periodičnim ponavljanjem ne<br />

naredimo nobene napake. Vprašanje pa je, <strong>za</strong> koliko moramo podaljšati<br />

<strong>za</strong>pis, da bo napaka v izračunu korelacije <strong>za</strong>radi robnega<br />

efekta sprejemljivo majhna. Oceno napako <strong>za</strong> primer zveznega<br />

signala A s<strong>in</strong> ωt:<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.9 Podobnost <strong>signalov</strong> 47<br />

r 12 (τ) = lim<br />

T →∞<br />

1<br />

2T<br />

A 2<br />

= lim<br />

T →∞ 2<br />

∫ T /2<br />

T /2<br />

A s<strong>in</strong>(ωt)A s<strong>in</strong>(ωt + τ) dt<br />

[<br />

cos(ωτ) − cos(ωT )<br />

2ωT s<strong>in</strong>(ωτ) ]<br />

. (2.9-15)<br />

Pregled (2.9-15) kaže, da se drugi člen v pravokotnem oklepaju<br />

manjša proti nič, ko T narašča preko vseh mej. V njem<br />

cos(ωT ) predstavlja periodično ponavljanje pogreška v odvisnosti<br />

od T , 1/(ωT ) pa določa njegovo manjšanje. Pogrešek je periodičen<br />

tudi glede na s<strong>in</strong>(ωτ).<br />

Člen cos(ωT ) ima najmanjšo vrednost, ko velja ωT = [(2n +<br />

1)/2]π. Ker je ω = 2π/T p , T p je perioda signala, iščemo pa<br />

velikost T , kjer so to zgodi, <strong>za</strong>pišemo:<br />

T (2n + 1) T p<br />

4<br />

. (2.9-16)<br />

s<strong>in</strong>(ωτ) je najmanjši, ko velja ωτ = mπ. Iz tega sledi, da se to<br />

zgodi pri:<br />

τ = m 2 T p . (2.9-17)<br />

Predpostavimo, da je pogrešek želeno majhen, če izberemo n <br />

10. V tem primeru je T nT p /2, oziroma:<br />

τ = m 2 T p . (2.9-18)<br />

Iz (2.9-17) sledi, največji τ <strong>za</strong> najnižjo frekvenčno komponento<br />

mora biti τ < T p , oziroma z upoštevanjem (2.9-18) dobimo<br />

τ T 5<br />

.<br />

Povzemimo. Pogrešek <strong>za</strong>radi končne dolž<strong>in</strong>e signala lahko m<strong>in</strong>imiziramo,<br />

če:<br />

<strong>za</strong>gotovimo, da velja T 5T p , kjer je T p najdaljša perioda<br />

komponent, ki sestavljajo signal,<br />

<strong>za</strong>mik med signaloma ne sme presegati 20% njihove dolž<strong>in</strong>e<br />

Na primer, želimo določiti avtokorelacijo govornega signala. Ta<br />

se nahaja v frekvenčnem območju med 300 <strong>in</strong> 3000 Hz. Najnižja


48 2. Signali<br />

frekvenca določa najdaljšo periodo, torej je T p = 1/300 = 3, 3 ×<br />

10 −3 [s]. Najmanjša sprejemljiva dolž<strong>in</strong>a signala potem je 5 ×<br />

3, 3×10 −3 [s] = 16, 7 [ms], največji <strong>za</strong>mik <strong>signalov</strong> je 0, 2×16, 7 =<br />

3, 33 [ms].<br />

2.9.7 Korelacijski koeficient<br />

Vrednost korelacije, ki jo izračunamo po podanih obrazcih, je<br />

odvisna od absolutne vrednosti podatkov oziroma amplitude signala.<br />

Za primer si oglejmo signale na naslednji sliki. Vidimo,<br />

x1( n) x3( n)<br />

Slika 2.9-4<br />

Postopek računanja križne<br />

korelacije aperiodičnih<br />

<strong>signalov</strong>)<br />

x2( n) x4( n)<br />

n<br />

n<br />

n<br />

n<br />

da sta si signala x 1 (n) <strong>in</strong> x 2 (n) podobna, razlikujeta se le v<br />

amplitudi. Enako velja <strong>za</strong> signala x 3 (n) <strong>in</strong> x 4 (n). Kljub temu<br />

pa se korelaciji r 12 (m) <strong>in</strong> r 34 (m) medsebojno razlikujeta, sa velja<br />

r 12 (m) < r 34 (m). Iz njune medsebojne primerjave bi lahko<br />

napačno sklepali, da sta si signala x 3 (n) <strong>in</strong> x 4 (n) bolj ”podobna”<br />

kot signala x 1 (n) <strong>in</strong> x 2 (n). Seveda temu ni tako. Izravnavo te<br />

anomalije da normali<strong>za</strong>cija rezultata s faktorjem:<br />

[<br />

1<br />

N<br />

N−1<br />

∑<br />

n=0<br />

x 2 1(n) × 1 N<br />

N−1<br />

∑<br />

n=0<br />

= 1 N<br />

x 2 2(n)<br />

] 1/2<br />

[ N−1 ∑<br />

x 2 1(n)<br />

n=0<br />

N−1<br />

∑<br />

n=0<br />

oziroma je normalizirana korelacija enaka:<br />

ρ 12 (m) =<br />

r 12 (m)<br />

[ N−1<br />

1 ∑<br />

x 2<br />

N<br />

1(n)<br />

n=0<br />

N−1<br />

∑<br />

n=0<br />

x 2 2(n)<br />

x 2 2(n)] 1/2<br />

, (2.9-19)<br />

] 1/2<br />

. (2.9-20)<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.10 Posplošeni signali 49<br />

ρ 12 (m) imenujemo korelacijski koeficient. Njegova vrednost je<br />

med −1 <strong>in</strong> +1. Pri tem +1 pomeni 100% ujemanje v direktni<br />

smeri, −1 pa 100% ujemanje v nasprotni smeri, na primer, da<br />

sta si signala v protifazi. Vrednost 0 pripada nični korelaciji. To<br />

pomeni, da sta signala medsebojno povsem neodvisna. Primer<br />

takih <strong>signalov</strong> so naključni ali slučajni signali.<br />

2.10 Posplošeni signali<br />

Za konec uvodnega pregleda <strong>signalov</strong> smo pustili posebne signale,<br />

ki jih ne moremo opisati z običajnimi funkcijami. Za te signale<br />

se uporabljata tudi term<strong>in</strong>a s<strong>in</strong>gularne funkcije <strong>in</strong> generalizirane<br />

funkcije. S tem poudarimo, da niso prave funkcije, ampak<br />

računski pripomočki, s katerimi v fiziki <strong>in</strong> <strong>in</strong>ženirstvu uč<strong>in</strong>kovito<br />

rešujemo različne probleme.<br />

Na vlogo <strong>in</strong> pomen posplošenih funkcij lahko gledamo kot na<br />

√ −1, ki ni realno število, pa ga kljub temu s pridom uporabljamo,<br />

posebno v elektrotehniki. Podobno je s posplošenimi funkcijami.<br />

Kljub temu, da jih v naravi ni, so izjemno koristno orodje.<br />

Posplošene funkcije so tiste, ki so pove<strong>za</strong>ne z impulzi. Opazujmo<br />

pulz δ h (t) (slika 2.10-1). Zanj velja:<br />

δ h (t) =<br />

Površ<strong>in</strong>a tega pul<strong>za</strong> je enaka<br />

{<br />

1<br />

2h<br />

pri − h < t < h<br />

0 sicer<br />

2h ·<br />

1<br />

2h = 1 .<br />

Z manjšanjem h narašča njegova amplituda δ h , vendar pri tem<br />

0<br />

d h ( n)<br />

-h 0 h t<br />

1<br />

2h<br />

Slika 2.10-1<br />

Pravokotni pulz δ h (t)<br />

njegova površ<strong>in</strong>a ostaja konstantna, to je enaka 1. Površ<strong>in</strong>a pul<strong>za</strong><br />

ostane enaka 1 tudi takrat, ko v limitnem postopku manjšamo h


50 2. Signali<br />

proti nič:<br />

∫ {<br />

h<br />

1<br />

lim dt = lim<br />

h→0 −h 2h {2h · 1<br />

h→0 2h } = 1 pri h = t = 0<br />

0 sicer<br />

Pri tem vidimo, da je rezultat limitnega postopka ”pulz” s šir<strong>in</strong>o<br />

nič (<strong>za</strong>to ga imenujemo impulz) <strong>in</strong> neskončno viš<strong>in</strong>o, med tem<br />

ko je plošč<strong>in</strong>a enaka 1. Funkcijo, ki opisuje ta pulz, imenujemo<br />

Diracova delta funkcija ali na kratko Diracova funkcija oziroma<br />

Diracov impulz. Označujemo ga s δ(t).<br />

To funkcijo je vpeljal P.A.M. Dirac v svojem delu ”The<br />

pr<strong>in</strong>ciples of Quantum Mechanics“.<br />

Njena matematična razlaga je bila vprašljiva kar nekaj<br />

časa. Teorijo delta funkcij je matematično vpeljal šele<br />

Francoski matematik Laurent Schwartz (1915 -) med<br />

leti 1945-1950 z razvojem teorije distribucij.<br />

2.10.1 Osnove teorije s<strong>in</strong>gularnih funkcij<br />

Lahko je spoznati, da ima δ(t) naslednjo lastnost:<br />

∫ ∞<br />

−∞<br />

x(t)δ(t) dt = x(0) , (2.10-1)<br />

□ □ □ □ <br />

<strong>za</strong>htevna snov<br />

kjer je x(t) omejena regularna funkcija s signalno osjo T ∈ R <strong>in</strong><br />

zveznostjo v trenutku t = 0. Zaradi tega, ker <strong>in</strong>tegral (z mejama<br />

−∞, ∞) produkta Diracove funkcije s to funkcijo funkcijo ohrani<br />

samo vrednost pri t = 0, imenujemo Diracovo funkcijo tudi enotski<br />

impulz. Grafični prikaz (2.10-1) vidimo na sliki 2.10-2.<br />

V razdelkih, ki sledijo, so namenjeni <strong>za</strong>htevnejšim študentom. V<br />

njih bomo povzeli osnove teorije s<strong>in</strong>gularnih funkcij, to je distribucij.<br />

Formalno bomo izpeljali nekatere lastnosti Diracove funkcije. S<br />

formalno izpeljavo želimo ustvariti le jasno sliko pomena Diracovega<br />

impul<strong>za</strong>, ne pa podati rigoroznega matematičnega opisa.<br />

Vrednost signala v trenutku t = 0 Pričnimo z dokazom veljavnosti<br />

(2.10-1), Najprej <strong>za</strong>menjajmo δ(t) z δ h (t). Tedaj je:<br />

∫ ∞<br />

−∞<br />

x(t)δ h (t) dt =<br />

∫ h<br />

−h<br />

x(t)δ h (t) dt = 1 ∫ h<br />

x(t) dt . (2.10-2)<br />

2h −h<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.10 Posplošeni signali 51<br />

x( t)<br />

z<br />

0<br />

•<br />

x( t) d ( t)<br />

dt x( 0)<br />

-•<br />

0<br />

d ( t)<br />

Slika 2.10-2<br />

t<br />

t<br />

Otipanje zveznega signala:<br />

x(0) = ∫ ∞<br />

x(t)δ(t) dt<br />

−∞<br />

0<br />

a<br />

t<br />

Preostali <strong>in</strong>tegral uženimo s prvim izrekom o povprečni vrednosti (izrek<br />

1, stran 32):<br />

∫ h<br />

−h<br />

x(t) dt = 2hx(ξ h ) − h ξ h h . (2.10-3)<br />

Rezultat v (2.10-3) upoštevamo v (2.10-2):<br />

∫ ∞<br />

−∞<br />

x(t)δ h (t) dt = 1<br />

2h 2hx(ξ h) = x(ξ h ) − h ξ h h .<br />

Sedaj manjšamo h proti nič:<br />

∫ ∞<br />

∫ h<br />

x(t)δ(t) dt = lim x(t)δ h (t) dt = x(ξ = 0) = x(0) ????<br />

−∞<br />

h→0 −h<br />

(2.10-4)<br />

Zapisani obrazec (2.10-1) torej velja.<br />

Vrednost signala v trenutku t = a Premaknimo x(t) po signalni<br />

osi <strong>za</strong> konstanto a. V tem primeru leva stran (2.10-1) postane:<br />

∫ ∞<br />

−∞<br />

x(t + a)δ(t) dt .<br />

Z <strong>za</strong>menjavo spremenljivk t s ζ = t+a (torej velja t = ζ −a <strong>in</strong> dt = dζ)<br />

po ponovitvi izpeljave v prejšnjem podrazdelku dobimo:<br />

∫ ∞<br />

−∞<br />

x(ζ)δ(ζ − a) dζ = x(ζ − a = 0) = x(ζ = a) = x(a) . (2.10-5)<br />

Vidimo, da je <strong>in</strong>tegral z mejama −∞, ∞ produkta signala <strong>za</strong> a <strong>za</strong>kasnjenim<br />

Diracovim impulzom enak vrednosti signala x(a) (slika 2.10-3).<br />

Rezultat v (2.10-5) je posplošitev (2.10-4).<br />

vzorčevanju signala.<br />

Uporabljamo ga pri<br />

□ □ □ □


52 2. Signali<br />

x( t)<br />

Slika 2.10-3<br />

Vrednost zveznega signala v<br />

trenutku t = a:<br />

x(a) = ∫ ∞<br />

x(t)δ(t − a) dt<br />

−∞<br />

t<br />

z•<br />

-•<br />

0<br />

d ( t)<br />

d ( t - a)<br />

x( t) d ( t)<br />

dt<br />

0 a<br />

x( a)<br />

t<br />

t<br />

0<br />

a<br />

2.10.2 L<strong>in</strong>earne komb<strong>in</strong>acije delta funkcij<br />

Iz def<strong>in</strong>icij, ki smo jih <strong>za</strong>pisali, vemo, da so si s<strong>in</strong>gularne funkcije<br />

enake, če povzročijo enak izračun <strong>in</strong>tegrala, v katerem nastopajo.<br />

Kasneje, pri obravnavi sistemov bomo rekli, da so si enake, če<br />

povzročijo enak odziv sistema.<br />

Sedaj si oglejmo še nekaj osnovnih operacij, kot so l<strong>in</strong>earna<br />

komb<strong>in</strong>acija, skaliranje <strong>in</strong> diferenciranje.<br />

Seštevanje <strong>in</strong> množenje s skalarjem Za regularna signala<br />

x 1 <strong>in</strong> x 2 ter skalarja α 1 <strong>in</strong> α 2 lahko <strong>za</strong>pišemo:<br />

∫ ∞<br />

−∞<br />

[α 1 x 1 + α 2 x 2 ] ξ(t) dt<br />

∫ ∞<br />

∫ ∞<br />

= α 1 x 1 (t)ξ 1 (t) dt + α 2 x 2 (t)ξ 2 (t) dt ,<br />

−∞<br />

kar velja pri vsaki funkciji ξ, <strong>za</strong> katero obstaja <strong>za</strong>pisani <strong>in</strong>tegral.<br />

Podobno lahko <strong>za</strong>pišemo <strong>za</strong> s<strong>in</strong>gularni funkciji ɛ 1 <strong>in</strong> ɛ 2 :<br />

∫ ∞<br />

−∞<br />

[α 1 ɛ 1 + α 2 ɛ 2 ] ξ(t) dt<br />

−∞<br />

∫ ∞<br />

∫ ∞<br />

= α 1 ɛ 1 (t)ξ 1 (t) dt + α 2 ɛ 2 (t)ξ 2 (t) dt ,<br />

−∞<br />

pri čemer ξ mora biti zvezna <strong>in</strong> potrebno krat odvedljiva funkcija<br />

v točki t = 0. Če sta s<strong>in</strong>gularni funkciji delta funkciji, na primer<br />

−∞<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.10 Posplošeni signali 53<br />

δ(t) <strong>in</strong> δ(t − τ), z upoštevanjem (2.10-5) dobimo enačbo:<br />

∫ ∞<br />

−∞<br />

[α 1 δ(t) + α 2 δ(t − τ)] ξ(t) dt<br />

= α 1<br />

∫<br />

T<br />

δ(t)ξ(t) dt + α 2<br />

∫<br />

T<br />

δ(t − τ)ξ(t) dt , (2.10-6)<br />

ki velja <strong>za</strong> vsako regularno funkcijo ξ(t), ki je zvezna v trenutkih<br />

t = 0 <strong>in</strong> t = τ.<br />

Vzorčenje <strong>signalov</strong> Rezultat v (2.10-6) je zelo pomemben v<br />

obdelavi <strong>signalov</strong>, saj kaže pot, kako iz zveznih <strong>signalov</strong> generirati<br />

časovno diskretne vzorce.<br />

Vzorčenje <strong>signalov</strong> je postopek, pri katerem zvezni signal x(t)<br />

opišemo s časovnim <strong>za</strong>poredjem otipkov x(n), to je z vrednostmi<br />

signala v časovno diskretnih trenutkih nT . To časovno <strong>za</strong>poredje<br />

imenujemo vzorec signala, elemente vzorca, to je vrednosti<br />

signala v trenutkih nT pa otipke. Postopek vzorčenja imenujemo<br />

tudi tipanje signala.<br />

Pogoj, da je opis signala z x(n) ekvivalenten zveznemu signalu<br />

x(t), je spoštovanje tako imenovanega Shannonovega pravila.<br />

Razlaga Shannonovega pravila presega namen tega <strong>za</strong>pisa.<br />

Oglejmo si le, kako ustvarimo vzorec signala.<br />

Če s<strong>in</strong>gularno funkcijo δ(t − τ) <strong>za</strong>menjamo z:<br />

δ(t − nT ) , n ∈ Z ,<br />

ustvarimo neskončno množico Diracovih impulzov, ki so medsebojno<br />

razmaknjeni <strong>za</strong> T . T imenujemo <strong>in</strong>terval tipanja. Če to<br />

množico uredimo v (neskončno) <strong>za</strong>poredje impulzov, s katerim<br />

množimo signal x(t) <strong>in</strong> zmnožek <strong>in</strong>tegriramo, iz (2.10-5) sledi, da<br />

je rezultat diskretni signal x(n) (slika 2.10-4).<br />

Odvod Diracovega impul<strong>za</strong> Oglejmo si še odvajanje Diracovega<br />

impul<strong>za</strong>. Če je x(t) regularna funkcija, zvezna v trenutku<br />

t = 0, potem s pomočjo <strong>in</strong>tegracije po delih lahko izpeljemo:<br />

∫ ∞<br />

−∞<br />

dx(t)<br />

dt<br />

∫ ξ(t) dt = x(t)ξ(t) ∣ ∞ ∞<br />

−<br />

−∞<br />

−∞<br />

x(t) dξ(t)<br />

dt<br />

dt<br />

∫ ∞<br />

= − x(t) dξ(t) dt (2.10-7)<br />

−∞ dt


54 2. Signali<br />

x( t)<br />

( t)<br />

0<br />

t<br />

( t n T )<br />

( t nT<br />

)<br />

x( t) ( t nT<br />

)<br />

z<br />

-<br />

x( n)<br />

z<br />

x( n) x( t) ( t nT)<br />

dt<br />

-<br />

0 n<br />

T<br />

x( a)<br />

t<br />

0<br />

a<br />

Slika 2.10-4<br />

Vzorčenje zveznega signala x(n) = ∫ ∞<br />

x(t)δ(t − nT ) dt, n ∈ Z.<br />

−∞<br />

t<br />

<strong>za</strong> vsak ξ, <strong>za</strong> katerega obstaja <strong>in</strong>tegral <strong>in</strong> velja x(−∞)ξ(−∞) =<br />

x(∞)ξ(∞) = 0. Če je x delta funkcija, potem se desna stran<br />

(2.10-7) pri vsaki regularni funkciji odvedljivi v točki t = 0 reducira<br />

na: ∫ ∞<br />

δ(t) dξ<br />

dt dt = −ξ(1) (0) .<br />

−∞<br />

Na enak nač<strong>in</strong> lahko določimo n-ti odvod delta funkcije: δ (n) .<br />

S ponavljanjem parcialnega odvajanja je preprosto poka<strong>za</strong>ti, da<br />

velja:<br />

∫ ∞<br />

−∞<br />

δ (n) ξ(t) dt = (−1) n ξ (n) (0) , (2.10-8)<br />

ki je def<strong>in</strong>irana <strong>za</strong> vsako regularno, n-krat v točki t = 0 zvezno<br />

<strong>in</strong> odvedljivo funkcijo ξ(t).<br />

2.10.3 Pove<strong>za</strong>va med enotsko stopnico<br />

<strong>in</strong> Diracovim impulzom<br />

Enotsko stopnico smo že opisali <strong>in</strong> grafično prika<strong>za</strong>li v podpoglavju<br />

o zelo uporabnih signalih. Zanjo velja:<br />

{<br />

0 t 0<br />

u(t) =<br />

.<br />

1 sicer<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


2.10 Posplošeni signali 55<br />

Označimo odvod stopnice z ˙u(t) <strong>in</strong> po delih izračunajmo <strong>in</strong>tegral:<br />

∫ b<br />

−a<br />

x(t) ˙u(t) dt =<br />

] ∫<br />

∣<br />

b b<br />

∣x(t)u(t) − ẋ(t)u(t) dt , (2.10-9)<br />

−a −a<br />

kjer sta a <strong>in</strong> b večja od nič. Izračunajmo desno stran (2.10-9).<br />

Najprej oglati oklepaj:<br />

x(t)u(t) ∣ b ∣ ∣∣ = x(t)u(t) b<br />

= x(b) − 0 ,<br />

−a −a<br />

kjer smo upoštevali def<strong>in</strong>icijo enotske stopnice u(t), <strong>za</strong>to x(−a)u(−a) =<br />

0 <strong>in</strong> ostane le vrednost pri zgornji meji. Izračunajmo še <strong>in</strong>tegral:<br />

∫ b<br />

−a<br />

ẋ(t)u(t) dt =<br />

∫ b<br />

0<br />

ẋ(t) dt = x(b) − x(0) .<br />

Dobljene rezultate vstavimo v (2.10-9) <strong>in</strong> dobimo:<br />

∫ b<br />

−a<br />

x(t) ˙u(t) dt = [x(b) − 0] − [x(b) − x(0)] = x(0) . (2.10-10)<br />

Iz primerjave (2.10-1) <strong>in</strong> (2.10-10) sledi, da velja:<br />

δ(t) = ˙u(t) . (2.10-11)<br />

Torej je odvod enotske stopnice enak Diracovi funkciji.


56 2. Signali<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


Poglavje 3<br />

Izražanje <strong>signalov</strong> z<br />

baznimi funkcijami<br />

Za kvantitativni opis <strong>signalov</strong> je potrebno, da jih izrazimo z<br />

določenimi funkcijami. Če imajo signali časovni potek, so to<br />

časovne funkcije. Signale lahko opišemo na več nač<strong>in</strong>ov. Izbira<br />

določenega opisa je odvisna od tega, kaj želimo poudariti. Na<br />

primer, da je opis primeren <strong>za</strong> matematično <strong>obdelavo</strong>, da ga je<br />

mogoče vizualno prika<strong>za</strong>ti, ali pa, da je primeren <strong>za</strong> določeno<br />

uporabo.<br />

Opis <strong>signalov</strong> z odseki premic, podali smo ga v razdelku 2.3<br />

na strani 13 je primeren le <strong>za</strong> signale, ki imajo odsekoma l<strong>in</strong>earni<br />

potek. Ta opis omogoča izračun le nekaterih lastnosti signala,<br />

na primer njegove moči ali energije. Glavna pomanjkljivost tega<br />

<strong>za</strong>pisa je, da z eno funkcijo (premico) opišemo le segment signala,<br />

ne pa ves signal. Zaradi tega se oziramo <strong>za</strong> funkcijami,<br />

ki zmorejo opisati signal v celoti. Primernost matematične obdelave<br />

ponavadi <strong>za</strong>hteva, da signal x(t) ponazorimo z l<strong>in</strong>earno<br />

komb<strong>in</strong>acijo množice znanih elementarnih funkcij. Te funkcije<br />

izbiramo tako, da imajo želene lastnosti, na primer neodvisno<br />

izračunljivost vsake elementarne funkcije, <strong>in</strong>variantnost na odvajanje<br />

itd.<br />

Elementarne funkcije glede na to, da jih uporabljamo kot<br />

”gradnike”, s katerimi sestavimo signale, imenujemo tudi osnovne<br />

ali bazne funkcije.<br />

57


58 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

3.1 Aproksimacija <strong>signalov</strong><br />

z eno bazno funkcijo<br />

Za <strong>za</strong>četek rešimo preprost primer, ko želimo realni signal x(t)<br />

izraziti z realno bazno funkcijo φ(t). Signal <strong>in</strong> bazna funkcija sta<br />

def<strong>in</strong>irani v istem časovnem <strong>in</strong>tervalu (a, b) <strong>za</strong> vsak t ∈ (a, b).<br />

Seveda želimo, da je opis signala z bazno funkcijo tak, da bo<br />

znotraj <strong>in</strong>tervala (a,b) čim bolj veljalo:<br />

x(t) ≈ c · φ(t) <strong>za</strong> vsak t ∈ (a, b) (3.1-1)<br />

kjer je c neka konstanta, ki jo lahko določimo tako, da bo produkt<br />

ˆx(t) = cφ(t) čim bolje pona<strong>za</strong>rjal signal x(t):<br />

ˆx(t) = cφ(t) ≈ x(t) (3.1-2)<br />

Funkcijo ˆx(t) imenujemo ocena ali približek funkcije x(t). Za<br />

oceno ”bliž<strong>in</strong>e”približka <strong>in</strong> prave vrednosti signala potrebujemo<br />

neko merilo. Najpogosteje se v ta namen uporablja srednji kvadratični<br />

pogrešek ali srednje kvadratično odstopanje. Def<strong>in</strong>iramo<br />

ga z:<br />

ε 2 = 1<br />

b − a<br />

= 1<br />

b − a<br />

∫ b<br />

a<br />

∫ b<br />

a<br />

[x(t) − ˆx(t)] 2 dt<br />

[x(t) − cφ(t)] 2 dt (3.1-3)<br />

Zakaj <strong>za</strong> mero podobnosti uporabljamo kvadratični pogrešek?<br />

Zaradi njegovih lastnosti. Zato poudarimo:<br />

Srednji kvadratični pogrešek je vedno pozitiven, saj je kvadratna<br />

funkcija. Zato tudi ima m<strong>in</strong>imum. Njegova srednja vrednost daje<br />

povprečno oceno pogreška nad vsem <strong>in</strong>tervalom, nad katerim je<br />

def<strong>in</strong>iran signal.<br />

Naloga, ki si jo <strong>za</strong>stavljamo, je poiskati konstanto c tako, da<br />

bo srednji kvadratični pogrešek med signalom <strong>in</strong> njegovim približkom<br />

čim manjši. Z drugimi besedami, poiskati moramo tako<br />

konstanto c, da bo ε 2 m<strong>in</strong>imalni. Potrebni pogoj <strong>za</strong> določitev<br />

m<strong>in</strong>imuma je, da je prvi odvod pogreška enak nič. Da je ekstrem,<br />

ki ga s tem izračunamo, tudi m<strong>in</strong>imalni pogrešek, sledi<br />

iz def<strong>in</strong>icije kvadratičnega pogreška. Seveda ga lahko preverimo<br />

s testom vrednosti drugega odvoda pogreška. Ta je v tej točki<br />

vedno pozitiven.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.1 Aproksimacija <strong>signalov</strong> z bazno funkcijo 59<br />

Rezultat parcialnega odvajanja, v katerem upoštevamo l<strong>in</strong>earnost<br />

operacij odvajanja <strong>in</strong> <strong>in</strong>tegriranja, je:<br />

∂ε<br />

∂c = ∂ 1<br />

∂c b − a<br />

= 1<br />

b − a<br />

= 1<br />

b − a<br />

∫ b<br />

a<br />

∫ b<br />

a<br />

∂<br />

∂c<br />

[<br />

0 − 2<br />

[x(t) − cφ(t) dt] 2 = 0<br />

[<br />

x 2 (t) − 2x(t)cφ(t) + c 2 φ 2 (t) ] dt<br />

∫ b<br />

Torej je konstanta c določena z:<br />

a<br />

c =<br />

x(t)φ(t) dt + 2c<br />

∫ b<br />

a<br />

∫ b<br />

a<br />

x(t)φ(t) dt<br />

φ 2 (t) dt<br />

∫ b<br />

a<br />

]<br />

φ 2 (t)dt = 0 .<br />

. (3.1-4)<br />

Če smo <strong>za</strong> φ(t) izbrali realno funkcijo, potem je imenovalec zgornje<br />

enačbe realno pozitivno število, na primer K:<br />

<strong>in</strong> je koeficient c enak:<br />

c = 1 K<br />

∫ b<br />

a<br />

K =<br />

∫ b<br />

x(t)φ(t) dt oziroma cK =<br />

a<br />

φ 2 (t) dt (3.1-5)<br />

∫ a<br />

b<br />

x(t)φ(t) dt (3.1-6)<br />

Ko kvadriramo oglate oklepaje v (3.1-3) <strong>in</strong> upoštevamo (3.1-5)<br />

ter (3.1-6), dobimo:<br />

ε 2 = 1 [∫ b ∫ b<br />

∫ b<br />

]<br />

x 2 (t) dt − 2x(t)cφ(t) dt + c 2 φ 2 dt<br />

b − a a<br />

a<br />

a<br />

= 1 [ ∫ b<br />

∫ b<br />

∫ b ]<br />

x 2 (t) dt − 2c x(t)φ(t) dt +c 2 φ 2 dt<br />

b − a a<br />

a<br />

a<br />

} {{ } } {{ }<br />

=cK<br />

=K<br />

= 1 [∫ b<br />

]<br />

x 2 (t) dt − 2c 2 K + c 2 K .<br />

b − a a<br />

Srednji kvadratni pogrešek je:<br />

ε 2 = 1 [∫ b<br />

]<br />

x 2 (t) dt − c 2 K<br />

b − a a<br />

(3.1-7)


60 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

Tako vidimo, da je ε 2 realen, nenegativen <strong>in</strong> da je pri izbiri konstante<br />

c z enačbo (3.1-4) m<strong>in</strong>imalen.<br />

ZGLED 3.1-1<br />

Aproksimirajmo signal, ki ga kaže slika 3.1-1, z bazično funkcijo φ(t) = s<strong>in</strong> t<br />

tako, da bo srednji kvadratični pogrešek m<strong>in</strong>imalen!<br />

Slika 3.1-1<br />

amplituda<br />

1<br />

0 2 t<br />

-1<br />

Konstanto c izračunamo s pomočjo (3.1-4):<br />

c =<br />

∫ 2π<br />

0<br />

∫ 2π<br />

0<br />

x(t) s<strong>in</strong> t dt<br />

=<br />

s<strong>in</strong> 2 t dt<br />

∫ π<br />

(1) s<strong>in</strong> t dt +<br />

0<br />

} {{ }<br />

=2<br />

∫ 2π<br />

(−1) s<strong>in</strong> t dt<br />

π<br />

} {{ }<br />

=2<br />

( 1<br />

2 t − 1 4 s<strong>in</strong> t ) ∣ ∣ ∣∣∣<br />

2π<br />

} {{<br />

0<br />

}<br />

K=π<br />

= 4 π<br />

. (3.1-8)<br />

Ocena signala je:<br />

ˆx(t) = 4 π s<strong>in</strong> t<br />

Oceno signala vidimo na sliki 3.1-2. M<strong>in</strong>imalni srednji kvadratični pogrešek ε 2<br />

izračunamo iz (3.1-7). V obravnavanem primeru sta meji a = 0, b = 2π, iz<br />

imenovalca (3.1-8) pa sledi, da je K = π, <strong>za</strong>to velja:<br />

ε 2 = 1 [∫ 2π<br />

] [ ∫<br />

x 2 (t)dt − c 2 K = 1 2π<br />

( ) 2<br />

x 2 4<br />

(t)dt − π]<br />

2π 0<br />

2π 0<br />

π<br />

] [1 − 8 ]<br />

π 2<br />

= 1<br />

2π<br />

[<br />

2π − 16<br />

π<br />

=<br />

≈ 0, 18943<br />

oziroma v procentih znaša 18,9% energije signala.<br />

♦<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.2 Aproksimacija <strong>signalov</strong> z l<strong>in</strong>earno komb<strong>in</strong>acijo baznih funkcij 61<br />

amplituda<br />

4/<br />

1<br />

x<br />

x<br />

Slika 3.1-2<br />

0<br />

<br />

2<br />

t<br />

-1<br />

3.2 Aproksimacija <strong>signalov</strong> z l<strong>in</strong>earno komb<strong>in</strong>acijo<br />

baznih funkcij<br />

Funkcijo x(t) izrazimo z l<strong>in</strong>earno komb<strong>in</strong>acijo baznih funkcij φ i<br />

iz množice baznih funkcij Φ:<br />

Φ = {φ 0 (t), φ 1 (t), φ 2 (t), . . . , φ N (t)} ,<br />

kjer je lahko <strong>in</strong>deks N v določenih primerih tudi neskončen. Pri<br />

končnem N lahko približek ˆx(t) izrazimo z:<br />

ˆx(t) = c 0 φ 0 (t) + c 1 φ 1 (t), +c 2 φ 2 (t), . . . , c N φ N (t)<br />

=<br />

N∑<br />

c n φ n (t) . (3.2-1)<br />

n=0<br />

Koeficienti c 0 , c 1 , c 2 , . . . , c N so praviloma elementi množice kompleksnih<br />

števil C. Problem izbire najprimernejšega <strong>za</strong>poredja baznih<br />

funkcij Φ <strong>in</strong> izbire najprimernejšega <strong>za</strong>poredja koeficientov<br />

{c n } v splošnem ni rešen. Poznamo pa že veliko število rešitev <strong>za</strong><br />

funkcije s posebnimi lastnostmi. Ena takih lastnosti, ki je ponavadi<br />

zelo <strong>za</strong>želena, je neodvisnost koeficientov c n ∈ C. Ta lastnost<br />

omogoča, da lahko določimo posamezne koeficiente, čeprav ne<br />

poznamo nobenega drugega koeficienta. Povedano drugače, če<br />

lahko l<strong>in</strong>earni komb<strong>in</strong>aciji, ki določa približek signala, dodamo<br />

člene, ne da bi morali že znane sprem<strong>in</strong>jati, potem so členi v<br />

l<strong>in</strong>earni komb<strong>in</strong>aciji (medsebojno) neodvisni.


62 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

3.3 Izražanje <strong>signalov</strong> z<br />

ortogonalnimi funkcijami<br />

Neodvisnost členov dosežemo, če je množica baznih funkcij φ(t) ∈<br />

Φ ortogonalna na celotnem časovnem <strong>in</strong>tervalu, <strong>za</strong> katerega velja<br />

opis funkcije x(t) z baznimi funkcijami. S term<strong>in</strong>om ortogonalna<br />

množica povemo, da so vsi elementi množice medsebojno pravokotni<br />

ali ortogonalni: Φ : φ m ⊥ φ n , kjer znak ⊥ pomeni<br />

”pravokotno“.<br />

Dve funkciji sta medsebojno ortogonalni, ko je njun skalarni<br />

produkt enak nič. Torej je:<br />

〈φ i φ j 〉 = 0<br />

pogoj ortogonalnosti. V razdelku 2.8.2 na strani 38 smo skalarni<br />

produkt <strong>za</strong> (kompleksne) funkcije def<strong>in</strong>irali z:<br />

〈φ m φ n 〉 =<br />

∫ b<br />

a<br />

φ m (t)φ ∗ n(t) dt ,<br />

ki lahko v primeru ortogonalnih funkcij <strong>za</strong>v<strong>za</strong>me le eno od dveh<br />

možnih vrednosti:<br />

∫ {<br />

b<br />

0 m ≠ n<br />

φ m (t)φ n (t)dt =<br />

φ m , φ n ∈ Φ (3.3-1)<br />

a<br />

m = n<br />

K n<br />

<strong>in</strong> to pri vseh m <strong>in</strong> n. Če je K n = 1, so bazne funkcije normirane<br />

<strong>in</strong> <strong>za</strong>to ortonormalne. Normiranje funkcij ima podoben pomen<br />

kot normiranje vektorjev. Pri tem je normirani vektor, ki ga imenujemo<br />

tudi enotski vektor, vektor z dolž<strong>in</strong>o enako enoti. Torej<br />

velja:<br />

〈aa〉 = a 2 = 1 .<br />

Podobno lastnost imajo tudi funkcije φ n ∈ Φ. Funkcija φ n (t) je<br />

normirana, če velja:<br />

〈φ n φ n 〉 =<br />

∫ b<br />

a<br />

φ 2 n(t) dt = 1 .<br />

Iz te kratke obravnave ortogonalnih funkcij lahko pov<strong>za</strong>memo<br />

naslednje: če obstaja <strong>za</strong>poredje funkcij φ i ∈ Φ, <strong>za</strong> katere velja:<br />

∫ {<br />

b<br />

0 m ≠ n<br />

φ m (t)φ n (t)dt =<br />

φ m , φ n ∈ Φ (3.3-2)<br />

a<br />

1 m = n<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.3 Izražanje <strong>signalov</strong> z ortogonalnimi funkcijami 63<br />

je vsaka funkcija φ n (t) ∈ Φ ortonormalna k ostalim funkcijam v<br />

tem <strong>za</strong>poredju. Tako <strong>za</strong>poredje funkcij imenujemo ortonormalno<br />

ali ortonormirano <strong>za</strong>poredje v <strong>in</strong>tervalu (a, b). Meji a <strong>in</strong> b sta<br />

pri tem lahko končni ali tudi neskončni. To je odvisno od narave<br />

problema, ki ga želimo rešiti s tem <strong>za</strong>poredjem ortonormiranih<br />

funkcij.<br />

3.3.1 Izražanje <strong>signalov</strong> z<br />

ortonormiranimi <strong>za</strong>poredji<br />

Tako kot lahko vsak vektor izrazimo z l<strong>in</strong>earno komb<strong>in</strong>acijo ortogonalnih<br />

enotskih vektorjev i, j, k, . . ., na primer z r = c 1 i +<br />

c 2 j + c 3 k, lahko tudi x(t) izrazimo z l<strong>in</strong>earno komb<strong>in</strong>acijo ortonormiranih<br />

funkcij:<br />

x(t) = ∑ n<br />

c n φ n (t) . (3.3-3)<br />

Pri tem moramo rešiti dve nalogi. Prva je izbrati <strong>za</strong>poredje baznih<br />

funkcij, druga pa je izračunati koeficiente c n tako, da bo<br />

pogrešek med orig<strong>in</strong>alnim signalom <strong>in</strong> njegovim približkom, sestavljenim<br />

iz <strong>za</strong>poredja ortonirmiranih koeficientov, m<strong>in</strong>imalen.<br />

Najprej si bomo ogledali možnosti določanja koeficientov c n , o<br />

izbiri ortonormiranih funkcij pa bo govora kasneje.<br />

Določitev koeficientov c n s pomočjo ortogonalnosti Kako<br />

določiti koeficiente c n ∈ C? Recimo, da vrsta na desni strani<br />

(3.3-3) konvergira k x(t) <strong>in</strong> da so bazne funkcije ortonormirane.<br />

Obe strani zgornje enačbe pomnožimo s φ m (t) <strong>in</strong> ju nato <strong>in</strong>tegriramo:<br />

∫ b<br />

a<br />

x(t)φ m dt =<br />

∫ b<br />

a<br />

= ∑ n<br />

∑<br />

c n φ n (t)φ m (t) dt<br />

n<br />

∫ b<br />

c n φ n (t)φ m (t) dt .<br />

a<br />

(3.3-4)<br />

Zaradi lastnosti ortogonalnosti, glej (3.3-2), je desna stran (3.3-4)<br />

različna od nič le, ko je m = n. V vseh ostalih primerih je enaka<br />

nič. Še več, desna stran je v primeru m = n enaka c n . Torej je<br />

enačba <strong>za</strong> izračun koeficienta c n :<br />

c n =<br />

∫ b<br />

a<br />

x(t)φ n (t) dt . (3.3-5)


64 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

Določitev koeficientov c n po metodi najmanjšega kvadratičnega<br />

pogreška Približek ˆx(t) s končno l<strong>in</strong>earno komb<strong>in</strong>acijo<br />

ortonormiranih baznih funkcij tvorimo tako kot smo to<br />

<strong>za</strong>pisali v (3.2-1), <strong>in</strong> ga določimo tako, da bo srednji kvadratni<br />

pogrešek med signalom <strong>in</strong> njegovo aproksimacijo m<strong>in</strong>imalen.<br />

Bazne funkcije, s katerimi določimo približek, vzemimo iz<br />

<strong>za</strong>poredja ortonormiranih funkcij Φ, želimo pa m<strong>in</strong>imizirati:<br />

ε 2 = 1<br />

b − a<br />

∫ b<br />

= 1<br />

b − a<br />

[<br />

x(t) − ˆx(t)] 2<br />

dt<br />

a<br />

∫ b<br />

a<br />

[<br />

x(t) − ∑ n<br />

c n φ(t)] 2<br />

dt<br />

(3.3-6)<br />

z izbiro koeficientov c n ∈ C. Pri tem vemo, da je pri iskanju m<strong>in</strong>imuma<br />

več spremenljivk potreben (vendar ne <strong>za</strong>dosten) pogoj,<br />

da je funkcija odvedljiva po vseh koeficientih <strong>in</strong> obstaja nabor<br />

vrednosti c 0 , c 1 , c 2 , . . . c N , ki <strong>za</strong>dovoljuje sistem enačb:<br />

∂ε<br />

= ∂ε = ∂ε = · · · = ∂ε = 0 .<br />

∂c 0 ∂c 1 ∂c 2 ∂c n<br />

Izračunajmo parcialni odvod po koeficientu c k . Ker je <strong>in</strong>terval<br />

konstanten, njegova vrednost ne vpliva na parcialni odvod. Zato<br />

lahko <strong>za</strong>pišemo:<br />

∂<br />

∂c k<br />

ε =<br />

∂ ∫ b<br />

∂c k<br />

a<br />

[<br />

x(t) − ∑ n<br />

c n φ n (t)] 2<br />

dt = 0 . (3.3-7)<br />

Izračunajmo kvadrat oglatega oklepaja v (3.3-7), upoštevajmo<br />

lastnost ortonormiranosti <strong>in</strong> uredimo rezultat računanja. Pri tem<br />

dobimo:<br />

∫<br />

∂ b<br />

x 2 (t) dt −<br />

∂ ∫ b<br />

2x(t) ∑ ∂c k a<br />

∂c k a n<br />

∫ [<br />

∂ b ∑<br />

∂c k a n<br />

c n φ n (t) dt+<br />

c n φ n (t)] 2<br />

dt = 0 .<br />

Po odvajanju po c k odpadejo vsi členi, ki ne vsebujejo koeficienta<br />

c k , <strong>za</strong>to iz zgornje enačbe ostaneta le dva člena:<br />

2<br />

∫ b<br />

a<br />

x(t)φ k (t) dt + 2c k<br />

∫ b<br />

a<br />

φ 2 k (t) dt = 0 . (3.3-8)<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.3 Izražanje <strong>signalov</strong> z ortogonalnimi funkcijami 65<br />

Ker smo <strong>za</strong> bazne funkcije izbrali ortonormirane funkcije, je vrednost<br />

<strong>in</strong>tegrala v drugem členu na levi strani (3.3-8) enak ena,<br />

<strong>za</strong>to je pot do enačbe, ki določa koeficient c k , kratka:<br />

c k =<br />

∫ b<br />

a<br />

x(t)φ n (t) dt . (3.3-9)<br />

Pri tem vidimo, da je rezultat našega truda v (3.3-9) enak (3.3-5).<br />

Iz tega sledi, da so koeficienti c k , ki jih izračunamo na osnovi<br />

lastnosti ortonormiranih funkcij, optimalni glede na srednji kvadratični<br />

pogrešek.<br />

3.3.2 Parsevalova identiteta<br />

Predpostavimo, da aproksimiramo signal x(t) z l<strong>in</strong>earno komb<strong>in</strong>acijo<br />

N + 1 baznih funkcij φ n (t). Pogrešek med signalom x(t)<br />

<strong>in</strong> njegovo oceno ˆx(t) izračunajmo z (3.3-6):<br />

ε 2 = 1 ∫ [<br />

2 b<br />

N∑<br />

x(t) − c n φ n (t)]<br />

dt<br />

b − a a<br />

n=0<br />

[ ∫<br />

= 1 b<br />

N∑<br />

∫ b<br />

x 2 (t) dt − 2 c n x(t)φ n (t) dt<br />

b − a a<br />

n=0 a<br />

} {{ }<br />

[ ∫<br />

= 1 b<br />

b − a a<br />

[ ∫<br />

= 1 b<br />

b − a a<br />

x 2 (t)dt − 2<br />

x 2 (t) dt −<br />

N∑<br />

c 2 n +<br />

n=0<br />

N∑<br />

n=0<br />

c 2 n<br />

]<br />

∫ b N∑<br />

+<br />

a n=0<br />

=c n<br />

N∑<br />

∫ ]<br />

b<br />

c 2 n φ 2 n(t) dt<br />

n=0 a<br />

} {{ }<br />

=1<br />

c 2 nφ 2 n(t) dt<br />

. (3.3-10)<br />

V izračunu (3.3-10) smo upoštevali l<strong>in</strong>earnost operacij seštevanja<br />

<strong>in</strong> <strong>in</strong>tegriranja, <strong>za</strong>to smo <strong>za</strong>menjali njuno <strong>za</strong>porednost <strong>in</strong> ortogonalnost.<br />

Iz (3.3-10) jasno sledi, da se z večanjem členov v oceni<br />

signala - to je funkcije, ki ga opisuje - manjša razlika med potekom<br />

signala <strong>in</strong> njegovo aproksimacijo. V limitnem postopku, v<br />

katerem večamo N preko vseh meja, postane praviloma srednji<br />

kvadratični pogrešek enak nič:<br />

[ ∫ ]<br />

1 b<br />

N∑<br />

lim<br />

N→∞ ε2 = lim<br />

x 2 (t)dt − c 2 n = 0<br />

N→∞ b − a<br />

a<br />

n=0<br />

]


66 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

oziroma:<br />

∫ b<br />

a<br />

x 2 (t) dt =<br />

∞∑<br />

c 2 n . (3.3-11)<br />

n=0<br />

V tem primeru je funkcija x(t) izražena z neskončnim <strong>za</strong>poredjem<br />

ortonormiranih funkcij; torej med oceno funkcije <strong>in</strong> funkcijo<br />

ni pogreška. Tako <strong>za</strong>poredje imenujemo polno ali kompletno<br />

<strong>za</strong>poredje. Enačbo (3.3-11), ki povezuje polno <strong>za</strong>poredje ortonormiranih<br />

funkcij <strong>in</strong> orig<strong>in</strong>alni signal, imenujemo Parsevalova<br />

identiteta.<br />

Rečemo lahko, da je pogoj Parsevalove <strong>in</strong>dentitete obstoj polnega<br />

<strong>za</strong>poredja, <strong>za</strong>to njen pomen poudarimo z naslednjim izrekom:<br />

IZREK 5 (Polno <strong>za</strong>poredje)<br />

Zaporedje ortonormiranih baznih funkcij Φ je polno takrat <strong>in</strong><br />

samo takrat, ko <strong>za</strong> vsak x ∈ L 2 (a, b) velja:<br />

∫ b<br />

a<br />

x 2 (t) dt =<br />

∞∑<br />

c 2 n .<br />

n=0<br />

To pomeni, da ne obstaja neničelna funkcija φ(t), ki ni član Φ <strong>in</strong><br />

<strong>za</strong> katero velja pogoj ortogonalnosti.<br />

Še kratek kometar k izreku. Če bi mogli najti tako funkcijo φ(t),<br />

tedaj bi bila ta funkcija ortogonalna na vsak člen v <strong>za</strong>poredju Φ<br />

<strong>in</strong> bi tudi sama pripadala temu <strong>za</strong>poredju.<br />

Pomen Parsevalove <strong>in</strong>dentitete je zelo velik, <strong>za</strong>to se bomo še<br />

večkrat vračali k njej. Napovejmo le njeno fizikalno razlago. S<br />

primerjavo (3.3-11) z obrazci <strong>za</strong> izračun moči vidimo, da Parsevalova<br />

<strong>in</strong>dentiteta govori o ohranitvi moči. To bomo posebej<br />

poudarili v poglavjih, v katerih bomo obravnavali moč spektrov.<br />

3.4 Ortogonalne funkcije<br />

Lastnosti sistemov (<strong>za</strong>poredij) ortogonalnih funkcij smo spoznali.<br />

Sedaj moramo takšna <strong>za</strong>poredja le še poiskati. Obstaja mnogo<br />

druž<strong>in</strong> ortonormalnih funkcij. Med njimi so najbolj znane trigo-<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.4 Ortogonalne funkcije 67<br />

nometrična <strong>za</strong>poredja:<br />

eksponentna <strong>za</strong>poredja:<br />

{. . . , 1, cos ωt, cos 2ωt, cos 3ωt, . . .}<br />

{. . . , 1, s<strong>in</strong> ωt, s<strong>in</strong> 2ωt, s<strong>in</strong> 3ωt, . . .}<br />

{. . . , e −ωt , 1, e ωt , e 2ωt , e 3ωt , . . . }<br />

<strong>in</strong> mnoge druge. V naslednjih razdelkih si bomo nekatere med<br />

njimi ogledali, zgoraj naštetim pa bomo posvetili posebno knjižico.<br />

3.4.1 Pregled nekaterih ortogonalnih funkcij<br />

Legendrove funkcije<br />

√<br />

2n + 1<br />

φ n (t) = P n (t)<br />

2<br />

, −1 ≤ t ≤ 1 ,<br />

kjer je P n (t) Legendrov pol<strong>in</strong>om:<br />

d n<br />

P n (t) = 1<br />

2 n n! dt n (tn − 1) n .<br />

Bazne Legendrove funkcije računamo s pomočjo Gram-Schmidtovega<br />

postopka (glej razdelek 3.4.3 na strani 71).<br />

Laguerrove funkcije<br />

φ n (t) = 1 n! e−t/2 L n (t) , 0 ≤ t < ∞ ,<br />

kjer je L n (t) Laguerrov pol<strong>in</strong>om reda n:<br />

Hermiteove funkcije<br />

φ n (t) =<br />

L n (t) = e t dn<br />

dt n tn (e −t ) .<br />

1<br />

2 n n! √ π e−t2 /2 H n (t) , −∞ < t < ∞<br />

kjer je H n (t) Hermiteov pol<strong>in</strong>om reda n:<br />

H n (t) = (−1) n t2 dn<br />

e<br />

dt n (e−t2 )<br />

Poleg naštetih ortogonalnih funkcij, ki so vse poljubnokrat odvedljive,<br />

so pomembne še Walsheve funkcije - posvečeno jim je<br />

naslednje poglavje.


68 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

3.4.2 Walsheve funkcije<br />

Polno <strong>za</strong>poredje Walshevih funkcij tvori <strong>za</strong>poredje medsebojno<br />

ortogonalnih pravokotnih valovnih oblik v končnem časovnem <strong>in</strong>tervalu.<br />

Te funkcije so tudi ortonormalne <strong>in</strong> lahko <strong>za</strong>vzemajo le<br />

dve vrednosti: +1 <strong>in</strong> −1. Časovni <strong>in</strong>terval def<strong>in</strong>iramo kot normirani<br />

<strong>in</strong>terval (0, 1), <strong>in</strong>tervale drugih dolž<strong>in</strong> pa dobimo s preprostim<br />

skaliranjem tega <strong>in</strong>tervala.<br />

Nekaj Walshevih funkcij, ki so <strong>in</strong>deksirane po številu prehodov<br />

funkcije skozi absciso, kaže slika 3.4-1. S slike vidimo, da<br />

1 1<br />

0 ( t)<br />

0 3 ( t)<br />

0<br />

1/2<br />

1 1/2<br />

1<br />

-1 -1<br />

1 1<br />

1 ( t)<br />

0 4 ( t)<br />

0<br />

1/2 1 1/2<br />

1<br />

-1 -1<br />

1 1<br />

2 ( t)<br />

0 5 ( t)<br />

0<br />

1/2<br />

1 1/2<br />

1<br />

-1 -1<br />

Slika 3.4-1<br />

Walsheve bazne funkcije<br />

se da Walsheve funkcije enostavno generirati z digitalnimi vezji.<br />

Vidimo tudi, da velja:<br />

{<br />

1 0 ≤ t ≤ 1<br />

φ 0 (t) =<br />

.<br />

0 sicer<br />

Ostale temeljne funkcije φ i lahko najdemo eno <strong>za</strong> drugo, če uporabimo<br />

pogoje ortonormiranosti. Splošni algoritem <strong>za</strong> ta postopek<br />

je tako imenovani Gram-Schmidtov postopek, ki je opisan v<br />

naslednjem razdelku. Tukaj le preverimo, ali grafi funkcij na sliki<br />

3.4-1 res pripadajo baznim funkcijam. Za φ 1 (t) velja:<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.4 Ortogonalne funkcije 69<br />

∫<br />

T<br />

φ 0 (t)φ 1 (t) dt = 0 ,<br />

∫<br />

T<br />

φ 2 1(t) dt = 1 <strong>in</strong> φ 1 (t) =<br />

{<br />

1 0 ≤ t ≤ 1/2<br />

−1 1/2 ≤ t ≤ 1<br />

.<br />

Temeljna funkcija φ 2 (t) mora ustre<strong>za</strong>ti enačbam:<br />

∫<br />

T<br />

φ 0 (t)φ 2 (t) dt = 0 ,<br />

∫<br />

T<br />

φ 1 (t)φ 2 (t) dt = 0<br />

<strong>in</strong><br />

∫<br />

T<br />

φ 2 2(t) dt = 1 .<br />

Funkcija φ 2 (t) na sliki 3.4-1 izpolnjuje vse te pogoje. Opisano<br />

pot lahko nadaljujemo <strong>za</strong> vse ostale temeljne funkcije. Pri tem<br />

vidimo, da moramo <strong>za</strong> bazno funkcijo φ k (t) rešiti k + 1 enačb.<br />

ZGLED 3.4-1<br />

Ponazorimo signal x(t) = 6t, 0 ≤ t ≤ 1 z Walshevimi funkcijami. Signal kaže<br />

slika 3.4-2.<br />

x ( t )<br />

6<br />

3<br />

-1 0 1/2 1<br />

t<br />

Slika 3.4-2<br />

Za izračun optimalnih koeficientov c k uporabimo (3.3-9). Prvi koeficient je<br />

c 0 <strong>in</strong> je enak:<br />

Drugi koeficient je c 1:<br />

c 1 =<br />

∫ 1/2<br />

0<br />

c 0 =<br />

∫ 1<br />

0<br />

6t(1)dt .<br />

∫ 1<br />

6t(1)dt + 6t(−1)dt = − 3<br />

1/2 2<br />

Po podobnem postopku določimo še ostale koeficiente. Dobimo:<br />

c 2 = 0, c 3 = − 3 4 , c4 = 0, c5 = 0, c6 = 0, c7 = − 3 8 · · ·<br />

Približek oziroma ocena ˆx(t) funkcije x(z) izražena z osmimi <strong>za</strong>porednimi baznimi<br />

funkcijami φ k (t), je:<br />

ˆx(t) = 3φ 0(t) − 3 2 φ1(t) − 3 4 φ3 − · · · 3<br />

2 k φ 2 k −1(t) · · · (3.4-1)<br />

Zgradbo ˆx(t) kaže slika 3.4-1. Srednji kvadratični pogrešek, ki ga naredimo z<br />

opisom x(t) z ˆx(t), če v<strong>za</strong>memo le prvih osem členov, je enak:<br />

.


70 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

6<br />

5<br />

4<br />

3<br />

2<br />

x( t)<br />

3<br />

3<br />

( t) ( t)<br />

x<br />

2<br />

0 1 0<br />

3 0 ( t)<br />

1<br />

Slika 3.4-3<br />

0<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1/2<br />

1<br />

t<br />

2 3<br />

3 ( t) ( t) ( t)<br />

x<br />

3 4<br />

3<br />

0 1 3 1<br />

3<br />

( t) ( t)<br />

2<br />

0 1<br />

1<br />

0<br />

6<br />

5<br />

4<br />

3<br />

1/2<br />

1<br />

t<br />

2 3 3<br />

3 ( t) ( t) ( t)<br />

x<br />

3 4 8<br />

2 3<br />

3 0( t) 1( t) 3( t)<br />

3 4<br />

0 1 3 7 2<br />

2<br />

1<br />

0<br />

1/2<br />

1<br />

t<br />

ε 2 = 1 [∫ 1<br />

b − a 0<br />

[ ∫ 1<br />

= 1<br />

1 − 0<br />

0<br />

[<br />

] ]<br />

x(t) 2 dt − c 2 0 + c 2 1 + c 2 7<br />

(<br />

(6t) 2 dt − 3 2 + − 3 ) 2 (<br />

+ − 3 ) ] 2<br />

= 12 − 765<br />

2 8<br />

64 ≈ 0.04685 .<br />

Čeprav je na videz razlika med ˆx(t) <strong>in</strong> x(t) velika, je srednji kvadratični<br />

pogrešek že po nekaj členih <strong>za</strong>poredja Φ(t) mali. Na primer, če bi v <strong>za</strong>poredju<br />

(3.4-1) dodali še en člen, bi padel srednji kvadratični pogrešek na četrt<strong>in</strong>o prejšnje<br />

vrednosti.<br />

♦<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.4 Ortogonalne funkcije 71<br />

3.4.3 Gram - Schmidtov postopek<br />

Gram-Schmidtov postopek je splošna metoda določanja ortonormalnih<br />

osnovnih funkcij ali baz <strong>za</strong> druž<strong>in</strong>o izbranih funkcij. Baze<br />

izpeljemo kar iz oblike funkcije, <strong>za</strong> katero jih iščemo. Bodi x i (t), i =<br />

1, 2, . . . , M množica <strong>signalov</strong>, <strong>za</strong> katero iščemo množico ortonormiranih<br />

baz φ i (t), i = 1, 2, . . . , N. Postopek iskanja je naslednji:<br />

Prva bazna funkcija je enaka normalizirani vrednosti izbranega<br />

signala, na primer:<br />

φ 1 (t) = x 1(t)<br />

√<br />

E1<br />

, (3.4-2)<br />

kjer je E 1 energija signala x 1 (t):<br />

∫<br />

E 1 = x 2 1(t) dt . (3.4-3)<br />

T<br />

Signal x 1 (t) lahko sedaj predstavimo z bazno funkcijo φ 1 (t)<br />

kot:<br />

x 1 (t) = x 1,1 φ 1 (t) , (3.4-4)<br />

kjer <strong>za</strong> koeficient x 1,1 iz (3.4-2) sledi x 1,1 = √ E 1 .<br />

Drugo bazno funkcijo določimo iz x 2 (t) tako, da najprej izračunamo<br />

projekcijo φ 1 (t) na x 2 (t). Enaka je:<br />

∫<br />

x 2,1 = x 2 (t)φ 1 (t) dt . (3.4-5)<br />

T<br />

Obrazec (3.4-5) je enak obrazcu (3.3-4) (paragraf 3.3.1 na<br />

strani 63), s pri katerem smo pri določitvi optimalne vrednosti<br />

koeficienta c i izkoristili lastnost ortogonalnih funkcij. S<br />

pomočjo projekcije (3.4-5) določimo pomožno funkcijo θ 2 :<br />

θ 2 (t) = x 2 (t) − x 2,1 φ 1 (t) , (3.4-6)<br />

ki je ortogonalna na φ 1 (t). Ortonormalnost dosežemo z normali<strong>za</strong>cijo<br />

- podelimo jo z njeno energijo. tako je bazna funkcija<br />

φ 2 (t) enaka:<br />

φ 2 (t) = θ 2(t)<br />

√<br />

Eθ2<br />

. (3.4-7)<br />

kjer je E θ2 = ∫ T θ2 2 (t) dt. Signal x 2(t) lahko <strong>za</strong>pišemo kot:<br />

x 2 (t) = x 2,1 φ 1 (t) + x 2,2 φ 2 (t) . (3.4-8)


72 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

Bazni funkciji φ 1 (t) <strong>in</strong> φ 2 (t) ima <strong>za</strong>radi nač<strong>in</strong>a, kako smo ju<br />

določili, energijo enote, torej sta normirani. Preizkusimo še<br />

njuno ortogonalnost:<br />

∫<br />

φ 1 (t)φ 2 (t) dt = √ 1 ∫<br />

φ 1 (t)θ 2 (t) dt<br />

Eθ2<br />

T<br />

= √ 1 ∫<br />

Eθ2<br />

T<br />

T<br />

= 1 √<br />

Eθ2<br />

[∫<br />

T<br />

[<br />

]<br />

φ 1 (t) x 2 (t) − x 2,1 φ 1 (t) dt<br />

x 2 φ 1 (t)dt − x 2,1<br />

∫<br />

= 1 √<br />

Eθ2<br />

[x 2,1 − x 2,1 ] = 0 .<br />

T<br />

]<br />

φ 2 1(t) dt<br />

(3.4-9)<br />

V vsakem naslednjem koraku določanja baznih funkcij določimo<br />

pomožno funkcijo θ j , ki je enaka razliki signala <strong>in</strong> njegovih do<br />

takrat - to je (j − 1) izračunanih komponent določenih z že<br />

dobljenimi baznimi funkcijami. Pomožno funkcijo normiramo<br />

z njeno energijo. Torej je j-ti korak:<br />

φ j (t) = θ j(t)<br />

√<br />

Eθj<br />

(3.4-10)<br />

kjer sta<br />

<strong>in</strong><br />

j−1<br />

∑<br />

θ j (t) = x j − x i,j φ i (t) (3.4-11)<br />

∫<br />

x i,j =<br />

T<br />

i=1<br />

x j (t)φ i dt , i = 1, 2, . . . j − 1 . (3.4-12)<br />

Z Gram-Schmidtovim potopkom dobimo število baznih funkcij,<br />

ki je zmeraj manjše ali kvečjemu enako številu <strong>signalov</strong> v<br />

množici, <strong>za</strong> katero smo iskali bazne funkcije.<br />

ZGLED 3.4-2<br />

Z Gram-Schmidtovim postopkom poiščimo bazne funkcije <strong>za</strong> signale s slike 3.4-4.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.4 Ortogonalne funkcije 73<br />

2<br />

s1( t) s2( t) s3( t)<br />

2 2<br />

1 1 1<br />

0 0 0<br />

1 1 1<br />

1 2 1 2 1 2<br />

Slika 3.4-4<br />

Signali, <strong>za</strong> katere iščemo<br />

bazne funkcije.<br />

Pričnimo s s 1(t). Prva bazna funkcija je:<br />

φ 1(t) = s1(t) √<br />

E1<br />

=<br />

1 √<br />

2<br />

s 1(t)<br />

Pomožna funkcija je:<br />

θ 2(t) = s 2(t) − s 2,1φ 1(t) , s 2,1 =<br />

<strong>in</strong> druga bazna funkcija:<br />

φ 2(t) = θ2(t) √<br />

Eθ2<br />

∫ T<br />

= s2(t) √<br />

Eθ2<br />

= 1 √<br />

2<br />

s 2(t)<br />

0<br />

s 2(t)φ 1(t)dt = 0<br />

Bazni funkciji φ 1(t) <strong>in</strong> φ 2(t) kaže slika 3.4-5.<br />

1 ( t)<br />

2 ( t)<br />

1/ 2<br />

1/<br />

2<br />

0 0<br />

1/ 2<br />

1/<br />

2<br />

1 2 1 2<br />

Slika 3.4-5<br />

Bazni funkciji, ki smo ju<br />

izpeljali z Gram-Schmidtovim<br />

postopkom <strong>za</strong> funkcije s slike<br />

3.4-4.<br />

Končna pomožna funkcija je:<br />

kjer sta<br />

Sledi:<br />

s 3,1(t) =<br />

s 3,2(t) =<br />

θ 3(t) = s 3(t) − s 3,1φ 1(t) − s 3,2φ(t) ,<br />

∫ T<br />

0<br />

∫ T<br />

0<br />

s 3(t)φ 1(t) dt =<br />

s 3(t)φ 2(t) dt =<br />

∫ 2<br />

1<br />

∫ 2<br />

1<br />

√<br />

2<br />

(<br />

− 1 √<br />

2<br />

)<br />

dt = −1<br />

√<br />

2<br />

1 √2 dt = 1 .<br />

θ 3(t) = s 3(t) + φ 1(t) − φ 2(t) = 0 .<br />

Zato lahko ves signal opišemo le z dvema baznima funkcijama, ki sta prika<strong>za</strong>ni<br />

na sliki 3.4-5. Z njima lahko aproksimiramo vse tri signale s slike 3.4-4.<br />

Če bi Gram-Schmidtov postopek <strong>za</strong>čeli z drugim signalom, na primer s 2(t),<br />

bi bile bazne funkcije drugačne, vendar pa bi predstavitve <strong>signalov</strong> v signalnem<br />

prostoru ostale podobne.<br />


74 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

3.5 Signalni prostor (drugič)<br />

V tem podpoglavju bomo <strong>za</strong>okrožili razlago l<strong>in</strong>earnih prostorov.<br />

Term<strong>in</strong> signalni prostor, ki smo ga vpeljali v podpoglavju 2.5 na<br />

strani 22, bomo predstavili kot vektorski prostor. V njem bomo<br />

signale izrazili kot vektorje.<br />

Razlago bomo pričeli s kratkim pregledom prika<strong>za</strong> <strong>in</strong> lastnosti<br />

vektorjev, nadaljevali pa z izbiro koord<strong>in</strong>atnega sistema, ki ga<br />

določajo ortonormirane bazne funkcije. V tem koord<strong>in</strong>atnem sistemu<br />

bomo signale predstavili kot l<strong>in</strong>earno komb<strong>in</strong>acijo, ki določa<br />

vektor.<br />

3.5.1 Vektorski prostor<br />

Na kratko povedano, vektorski prostor je l<strong>in</strong>earni prostor. V<br />

splošnem lahko vektorske prostore delimo na končno <strong>in</strong> na neskončno<br />

dimenzionalne prostore. V nadaljevanju se bomo omejili<br />

na N-dimenzionalne prostore. V njem lahko N-dimenzionalni<br />

vektor v i = [v i1 , v i2 , · · · , v iN ] predstavimo kot l<strong>in</strong>earno komb<strong>in</strong>acijo<br />

enotskih vektorjev e m , 1 m N, torej kot:<br />

v i =<br />

N∑<br />

v im e m ,<br />

m=1<br />

kjer so enotski vektorji e m po svoji def<strong>in</strong>iciji vektorji z dolž<strong>in</strong>o 1,<br />

koeficienti v im pa razmerje dolž<strong>in</strong>e projekcije vektorja v i na m-ti<br />

enotski vektor z dolž<strong>in</strong>o enotskega vektorja (torej, kolikokrat je<br />

projekcija daljša od enotskega vektorja).<br />

Vektorski prostor, ki nas <strong>za</strong>nima, ima def<strong>in</strong>iran skalarni produkt:<br />

N∑<br />

〈v i , v j 〉 = v im v jm ,<br />

m=1<br />

kjer sta v i <strong>in</strong> v j N-dimenzionalna vektorja. Če sta si ta dva<br />

vektorja medsebojno pravokotna, torej v i ⊥ v j , potem je njun<br />

skalarni produkt enak nič.<br />

Za vektorje rečemo, da so normirani, kadar je njihova Evklidska<br />

norma (norma p = 2) enaka ena:<br />

∑<br />

‖v i ‖ = √ N vim 2 = 1 .<br />

m=1<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


3.5 Signalni prostor (drugič) 75<br />

Pri tem je zelo razširjena navada, da pri simbolnem <strong>za</strong>pisu Evklidske<br />

norme izpuščamo <strong>in</strong>deks, ki določi tip norme, torej jo<br />

označimo le z ‖ · ‖ <strong>in</strong> jo na kratko imenujemo norma. Normo<br />

lahko določimo tudi s skalarnim produktom:<br />

‖v i ‖ = 〈v i , v i 〉 1/2 .<br />

Norma (natančneje Evklidska norma) vektorja seveda ni nič drugega<br />

kot dolž<strong>in</strong>a vektorja. Normirani <strong>in</strong> ortogonalni vektorji tvorijo<br />

množico ortonormiranih vektorjev. Zanje velja:<br />

{<br />

0 i ≠ j<br />

v i · v j =<br />

.<br />

1 i = j<br />

Dva vektorja izpolnjujeta trikotniško neenakost:<br />

‖v i + v j ‖ ‖v i ‖ + ‖v j ‖ ,<br />

v kateri velja enakost le takrat, ko sta vektorja v i <strong>in</strong> v j vzporedna,<br />

torej ko velja v i = αv j , kjer je α pozitivno realno število.<br />

V razdelku 2.8.3 na strani 39 smo izpeljali Cauchy-Schwartz neenakost.<br />

Ta se pri vektorjih glasi:<br />

‖v i · v j ‖ ‖v i ‖ · ‖v j ‖ ,<br />

v kateri velja enakost, če velja v i = αv j .<br />

Na kraju še povzemimo postopek konstrukcije množice ortonormalnih<br />

vektorjev iz množice N-dimenzionalnih vektorjev, ki<br />

jo naredimo z Gram-Schmidtovim postopkom. Konstrukcijo ortonormiranih<br />

vektorjev pričnemo z izbiro poljubnega vektorja iz<br />

množice vektorjev, <strong>za</strong> katero konstruiramo ortonormalno množico.<br />

Naj bo izbrani vektor na primer vektor v 1 . Ta vektor najprej<br />

normaliziramo, torej podelimo z njegovo dolž<strong>in</strong>o. Dobimo:<br />

u 1 = v 1<br />

‖v 1 ‖<br />

Nadaljujemo z izbiro naslednjega vektorja, na primer vektorja<br />

v 2 . Od njega najprej odštejemo projekcijo v 2 na u 1 :<br />

<strong>in</strong> rezultat normaliziramo:<br />

u ′ 2 = v 2 − 〈v 2 , u 1 〉u 1 .<br />

u 2 = u′ 2<br />

‖u ′ 2 ‖ .<br />

.


76 3. Izražanje <strong>signalov</strong> z baznimi funkcijami<br />

Postopek nadaljujemo z izbiro vektorja v 3 . Najprej odštejemo<br />

njegovo projekcijo na vektorja u 1 <strong>in</strong> u 2 :<br />

u ′ 3 = v 3 − 〈v 3 , u 1 〉u 1 − 〈v 3 , u 2 〉u 2 .<br />

Po normali<strong>za</strong>ciji u ′ 3 nadaljujemo s postopkom, dokler ne izčrpamo<br />

vseh vektorjev, <strong>za</strong> katere bi radi zgradili množico baznih, ortonormalnih<br />

funkcij. Gram-Schmidtov postopek na splošno (<strong>za</strong><br />

funkcije) smo izpeljali že v prejšnjem razdelku (stran 71).<br />

3.5.2 Signalni prostor je vektorski prostor<br />

V razdelkih od 2.5 do vključno 2.8 smo izpeljali koncept signalnega<br />

prostora. Pri tem smo izhajali iz predstave, da je funkcija<br />

nad <strong>in</strong>tervalom (a, b), neskončno dimenzionalni vektor. S to predstavo<br />

smo izpeljali skalarni produkt med funkcijami, def<strong>in</strong>irali<br />

ortogonalnost <strong>in</strong> ortonormalnost med funkcijami, izmerili signal<br />

z normami itd. Zato dopolnimo gledanje na signalni prostor z<br />

naslednjo def<strong>in</strong>icijo:<br />

DEFINICIJA 13 (Signalni prostor)<br />

Signalni prostor je l<strong>in</strong>earni vektorski prostor, v katerem je def<strong>in</strong>iran<br />

skalarni produkt.<br />

To def<strong>in</strong>icijo ponavadi omejimo še z:<br />

s 3<br />

1<br />

2<br />

s 2<br />

1<br />

s 1<br />

-1<br />

1 2<br />

2<br />

Slika 3.5-1<br />

DEFINICIJA 14 (Koord<strong>in</strong>atni sistem signalnega prostora)<br />

Koord<strong>in</strong>atni sistem signalnega prostora tvorijo bazne funkcije.<br />

S tem pogledom na signalni prostor je seveda izražanje signala<br />

z ortonormalnimi baznimi funkcijami opis vektorja z njegovimi<br />

koord<strong>in</strong>atami. Povedano drugače, signal si lahko predstavljamo<br />

kot vektor, katerega izrazimo z l<strong>in</strong>earno komb<strong>in</strong>acijo<br />

koord<strong>in</strong>at - to je l<strong>in</strong>earno komb<strong>in</strong>acijo baznih funkcij.<br />

3.5.3 Predstavitev <strong>signalov</strong> kot vektorjev<br />

Opis signala kot vektorja je geometrijsko predstavljiv le pri dvo <strong>in</strong><br />

tri-dimenzionalnih signalnih prostorih. Ker je klasična predstavitev<br />

vektorja s poltrakom nerodna, ponavadi predstavimo vektor<br />

s točko na konici vektorja. Na primer, signale s slike 3.4-4 iz<br />

prejšnjega primera predstavimo v signalnem s točkami s koord<strong>in</strong>atami<br />

(slika 3.5-1):<br />

(√ )<br />

s 1 = 2, 0<br />

, s 2 =<br />

(<br />

0, √ )<br />

2<br />

, s 3 = (−1, 1) .<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


Poglavje 4<br />

Sistemi<br />

Sistem je del okolja, ki ustvarja signale, ali pa se na njih odziva.<br />

Zato s sistemi lahko preoblikujemo signale. Na splošno sisteme<br />

opisujemo na dva nač<strong>in</strong>a:<br />

s spremenljivkami stanja sistema,<br />

z vhodno-izhodnim opisom.<br />

Opis s spremenljivkami stanja se uporablja predvsem v teoriji<br />

sistemov <strong>in</strong> regulacijah, <strong>za</strong>to ga ne bomo obravnavali. Pozornost<br />

bomo namenili le vhodno-izhodnemu opisu. Te bomo imenovali<br />

kar vhodno-izhodni sistemi. Predstavili jih bomo kot operatorje,<br />

ki vhodne signale preslikajo ali transformirajo v izhodne signale.<br />

Mnogokrat bomo<br />

<strong>za</strong>radi krajšega pisanja<br />

vhodne signale<br />

imenovali kar vhod <strong>in</strong><br />

izhodne signale izhod.<br />

4.1 Vhodno-izhodni sistemi<br />

Vhodno-izhodni sistem določa transformacija T, ki priredi določeni<br />

množici vhodnih <strong>signalov</strong> (končno) množico izhodnih <strong>signalov</strong>.<br />

Povedano drugače: transformacija T je določena z množico parov,<br />

ki jih tvorijo vrednosti vhodnega signala <strong>in</strong> njej pripadajoči<br />

možni izhodni signali. Take sisteme formalno opišemo z naslednjo<br />

def<strong>in</strong>icijo:<br />

DEFINICIJA 15 (Vhodno-izhodni sistem)<br />

Vhodno-izhodni sistem določa vhodna množica <strong>signalov</strong> V, izhodna<br />

množica <strong>signalov</strong> Y <strong>in</strong> podmnožica pravil preslikave T, ki so<br />

pove<strong>za</strong>ni s produktom V × Y.<br />

77


78 4. Sistemi<br />

Vsak par (v, y) ∈ T z v ∈ V <strong>in</strong> y ∈ Y imenujemo vhodnoizhodni<br />

par, ki vhodu v priredi izhod y.<br />

Grafično predstavitev def<strong>in</strong>icije sistema vidimo na sliki 4.1-1. Za<br />

splošni sistem vidimo, da ima lahko pri danem vhodu v množico<br />

mogočih izhodnih parov (v, y) ∈ T. Kadar sistem priredi danemu<br />

vhodu le en izhod, imenujemo tak sistem enolično preslikavni<br />

sistem. V nadaljevanju se bomo ukvarjali le s sistemi z<br />

Slika 4.1-1<br />

Vhodno-izhodni sistem, levo:<br />

splošni vhodno-izhodni sistem,<br />

desno: vhodno-izhodni sistem<br />

z enolično preslikavo<br />

izhodna mnoica Y<br />

y v<br />

v<br />

vhodna mnoica V<br />

izhodna mnoica Y<br />

y<br />

v<br />

vhodna mnoica V<br />

enolično preslikavo. Zato jih bomo na kratko imenovali sistemi.<br />

Na splošno predstavimo transformacijo oziroma preslikavo T z:<br />

y(t) = T [v(t)] . (4.1-1)<br />

<strong>in</strong> prikažemo z blokovnim diagramom, kot to kaže slika 4.1-2.<br />

Sistem torej priredi vhodnemu signalu izhodni signal. Pri tem<br />

poudarimo:<br />

Transformacija ne označuje matematične funkcije, s katero bi lahko<br />

pri danem vhodu izračunali izhod. V ta namen rabimo matematični<br />

model ali, na kratko, model sistema.<br />

Slika 4.1-2<br />

Abstraktna predstavitev<br />

splošnega sistema<br />

v( t)<br />

v( n)<br />

T<br />

y( t) = T[ v( t)]<br />

y( n) = T[ v( n)]<br />

Pri časovno zveznih sistemih model ponavadi opišemo s sistemom<br />

diferencialnih enačb, pri časovno diskretnih pa s sistemom<br />

diferenčnih enačb.<br />

ZGLED 4.1-1 Opis sistema s transformacijo<br />

Opišimo sistem na sliki 4.1-3 s transformacijo <strong>in</strong> <strong>za</strong>nj poiščimo pove<strong>za</strong>vo med<br />

vhodom <strong>in</strong> izhodom.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.1 Vhodno-izhodni sistemi 79<br />

u( t)<br />

i( t)<br />

R<br />

L<br />

Slika 4.1-3<br />

Električno vezje z elementoma<br />

R <strong>in</strong> L<br />

Vhod sistema je napetost u(t), izhod pa tok i(t). Zapis tega sistema s<br />

transformacijo je:<br />

i(t) = T[u(t)] . (4.1-2)<br />

Diferencialna enačba, s katero opišemo model sistema, to je dogajanje v sistemu,<br />

dobimo iz Ohmovega <strong>in</strong> 1. Kirchoffovega <strong>za</strong>kona. Za ta sistem velja:<br />

L di(t)<br />

dt<br />

+ Ri(t) = u(t) . (4.1-3)<br />

Vidimo, da <strong>za</strong>pis v obliki transformacije v (4.1-2) simbolično podaja eksplicitno<br />

pove<strong>za</strong>vo med vhodom <strong>in</strong> izhodom sistema, ki jo določa (4.1-3), ki pa je l<strong>in</strong>earna<br />

diferencialna enačba prvega reda.<br />

♦<br />

Posvetimo se še malo razliki med fizikalnim sistemom <strong>in</strong> njegovim<br />

modelom. Vezje na sliki 4.1-3 lahko na primer predstavlja:<br />

stvarno pove<strong>za</strong>vo upora, tuljave <strong>in</strong> generatorja,<br />

analogijo fizikalnega sistema, ki ga opišemo z diferencialno<br />

enačbo, ki je podobna (4.1-3), ali<br />

nadomestno vezje fizikalnega sistema, ki vsebuje poljubno število<br />

fizikalnih elementov.<br />

V slednjem primeru govorimo, da modeliramo sistem z določeno<br />

nenatančnostjo. Poudarimo:<br />

Natančnost modela, s katerim opišemo neki fizikalni sistem, je ena<br />

najtežjih <strong>in</strong> najpomembnejših <strong>in</strong>ženirskih nalog<br />

Kot <strong>za</strong>nimivost omenimo, da so na primer <strong>za</strong> določitev modela<br />

raketnega ”boosterja“ (pomožnega <strong>za</strong>gonskega raketnega<br />

sistema) rakete Saturn V, s katerim so določali viš<strong>in</strong>o doleta, uporabili<br />

diferencialno enačbo 27 reda! Šele tako <strong>za</strong>pleteni enačbi so<br />

verjeli, da dovolj natančno opisuje dogajanje po izstrelitvi rakete<br />

v času delovanja <strong>za</strong>gonskega raketnega sistema. Zaključimo z<br />

naslednjim poudarkom:


80 4. Sistemi<br />

Ponavadi, ko govorimo o sistemu, mislimo na matematični model<br />

sistema <strong>in</strong> ne na njegovo fizikalno sestavo.<br />

Sistemi (s katerimi se bomo ukvarjali) so po dani def<strong>in</strong>iciji vzročni<br />

ali kav<strong>za</strong>lni. To pomeni, da je izhodni signal sistema reakcija sistema<br />

na vhodni signal.<br />

4.1.1 Primeri sistemov<br />

Kaj vse pa je lahko sistem? Primer fizikalnega sistema je električni<br />

bojler. Vhodni signal oziroma vhod vanj je električna<br />

napetost u(t), izhod pa temperatura vode ϑ(t) (slika 4.1-4).<br />

Slika 4.1-4<br />

Električni bojler predstavljen<br />

kot sistem.<br />

elektrièni<br />

bojler<br />

v( t) = u( t) y ( t)<br />

= J( t)<br />

napetost<br />

T<br />

temperatura<br />

vode<br />

Sistem je lahko na primer ojačevalnik napetosti. V zgledu<br />

2.4-1 na strani 18 smo ojačali signal. Signal je lahko na primer<br />

električna napetost, ki jo generira mikrofon. Da lahko predvajamo<br />

glas po zvočniku, ga moramo ustrezno ojačati, na primer<br />

<strong>za</strong> deset krat (slika 4.1-5).<br />

Slika 4.1-5<br />

Primer sistema, ki ga sestavlja<br />

idealni ojačevalnik z ojačenjem<br />

10.<br />

v( t)<br />

idealni ojaèevalnik<br />

z ojaèenjem 10<br />

T<br />

y( t) = 10v( t)<br />

4.1.2 Časovno zvezni <strong>in</strong><br />

časovno diskretni sistemi<br />

Več<strong>in</strong>oma bomo proučevali časovne signale. Zato so tudi sistemi,<br />

ki jih bomo opisali, določeni s časovnimi funkcijami. Glede na<br />

obliko vhodnega <strong>in</strong> izhodnega signala ločimo časovno zvezne sisteme<br />

<strong>in</strong> časovno diskretne sisteme. Pri časovno zveznih sistemih<br />

je signalna os zvezna, pri časovno diskretnih sistemih pa diskretna.<br />

V tem poglavju bomo spoznali le splošne lastnosti sistemov,<br />

ki veljajo tako <strong>za</strong> časovno zvezne kot tudi <strong>za</strong> časovno diskretne<br />

sisteme.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.1 Vhodno-izhodni sistemi 81<br />

4.1.3 Povezovanje sistemov<br />

Sisteme lahko gradimo iz elementov tako, da jih medsebojno povezujemo.<br />

Pri tem uporabljamo pri grafični predstavitvi pove<strong>za</strong>v<br />

različne simbole. Najpogosteje uporabljani so narisani v naslednjih<br />

slikah.<br />

v t 1( )<br />

v2( t)<br />

v3( t)<br />

v t 1( )<br />

v2( t)<br />

y( t) = v ( t) + v ( t) + v ( t)<br />

1 2 3<br />

y( t) = v ( t) ◊v ( t)<br />

1 2<br />

Slika 4.1-6<br />

Primera simbolov gradnikov<br />

sistemov; zgoraj: seštevalnik<br />

<strong>signalov</strong>, spodaj: množilnik<br />

Tudi sisteme lahko medsebojno povezujemo. Pri tem poznamo<br />

dve osnovni pove<strong>za</strong>vi: vzporedno ali paralelno ter <strong>za</strong>poredno<br />

ali serijsko ve<strong>za</strong>vo (slika 4.1-7). Najprej si oglejmo vzporedno<br />

v( t)<br />

v( t)<br />

sistem<br />

1<br />

sistem<br />

2<br />

sistem<br />

1<br />

y t 1( )<br />

y t 1( )<br />

y2( t)<br />

sistem<br />

2<br />

y( t)<br />

y( t)<br />

Slika 4.1-7<br />

Vzporedna (zgoraj) <strong>in</strong><br />

<strong>za</strong>poredna (spodaj) pove<strong>za</strong>va<br />

dveh sistemov.<br />

pove<strong>za</strong>vo. Vhod v skupni sistem je v(t), izhod pa določa vsota<br />

izhodov pove<strong>za</strong>nih sistemov:<br />

y(t) = y 1 (t) + y 2 (t) = T 1 [v(t)] + T 2 [v(t)] = T[v(t)] . (4.1-4)<br />

Pri <strong>za</strong>poredni ve<strong>za</strong>vi sistemov pa je izhod skupnega sistema določen<br />

z:<br />

y(t) = T 2 [y 1 (t)] = T 2 (T 1 [v(t)]) = T[v(t)] . (4.1-5)<br />

Enačb (4.1-4) <strong>in</strong> (4.1-5) brez poznavanja modela sistema ni mogoče<br />

nadalje poenostaviti.


82 4. Sistemi<br />

Slika 4.1-8<br />

Primer sistema. Določa ga<br />

enačba (4.1-7)<br />

v( t)<br />

sistem<br />

1<br />

sistem<br />

2<br />

y t 1( )<br />

y2( t)<br />

v3( t)<br />

sistem<br />

3<br />

sistem<br />

4<br />

y3( t)<br />

y4( t)<br />

y( t)<br />

ZGLED 4.1-2 Pove<strong>za</strong>ve sistemov<br />

Izpišimo simbolični <strong>za</strong>pis sistema, ki ga kaže slika 4.1-8.<br />

v 3(t) = y 1(t) + y 2(t) = T 1[v(t)] + T 2[v(t)]<br />

y 3(t) = T 3[v 3(t)] = T 3<br />

[<br />

T1[v(t)] + T 2[v(t)] ] (4.1-6)<br />

<strong>in</strong> izhod sestavljenega sistema je<br />

y(t) = y 3(t)y 4(t) = T 3<br />

[<br />

T1[v(t)] + T 2[v(t)] ] T 4[v(t)] . (4.1-7)<br />

Iz tega primera vidimo, da podani opis pove le to kako so sistemi medsebojno pove<strong>za</strong>ni,<br />

ne daje pa modela skupnega sistema. Če pa so modeli sistemov, iz katerih<br />

je zgrajen sistem, znani, ta opis omogoča (olajša) <strong>za</strong>pisati model sestavljenega<br />

sistema.<br />

♦<br />

4.2 Sistemi s povratno <strong>za</strong>nko<br />

Zelo pomembni sistemi so sistemi s povratno <strong>za</strong>nko. Uporabljamo<br />

jih v regulacijah. S term<strong>in</strong>om regulacije ponavadi opredeljujemo<br />

uravnavanje delovanja ali lastnosti brez posega človeka. Včasih to<br />

dejstvo poudarimo s pridevnikom avtomatsko, torej avtomatska<br />

regulacija.<br />

Primeri regulacij so številni. Prav<strong>za</strong>prav si danes le težko<br />

predstavljamo življenje brez regulacijskih sistemov, ki samodejno<br />

skrbijo <strong>za</strong> mnoge naše potrebe: na primer <strong>za</strong> vzdrževanje konstantne<br />

temperature v stanovanju ali na delovnem mestu, usmerjanje<br />

letala ali ladje z ”avtopilotom“ itd. .<br />

Osnovno zgradbo sistema s povratno <strong>za</strong>nko kaže slika 4.2-<br />

1. Regulator je fizični sistem, ki ga načrtamo tako, da bo imel<br />

skupni sistem želene lastnosti oziroma karakteristike. Med najpomembnejše<br />

lastnosti sodi stabilnost. Senzor meri signal, ki ga<br />

želimo regulirati. Signal pogreška e(t) je mera razlike med želeno<br />

<strong>in</strong> dejansko vrednostjo izhoda sistema.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.2 Sistemi s povratno <strong>za</strong>nko 83<br />

v( t)<br />

regulator<br />

e( t) sistem 1:<br />

T<br />

1<br />

y3( t)<br />

y ( t) = v ( t)<br />

1 2<br />

senzor<br />

sistem 3:<br />

T<br />

3<br />

regulirani<br />

sistem<br />

sistem 2:<br />

T<br />

2<br />

v1( t) = y( t)<br />

y( t)<br />

Slika 4.2-1<br />

Sistem s povratno <strong>za</strong>nko.<br />

Sistem s povratno <strong>za</strong>nko na sliki 4.2-1 opišemo z:<br />

e(t) = v(t) − y 3 (t) = v(t) − T 3 [y(t)]<br />

<strong>in</strong><br />

y(t) = T 2 [v 2 (t)] = T 2 [ T 1 [e(t)] ] (4.2-1)<br />

= T 2 [ T 1 [v(t) − T 3 [y(t)] ] ] . (4.2-2)<br />

Izhod sistema smo izrazili kot funkcijo vhoda <strong>in</strong> izhoda sistema.<br />

Zapisa (4.2-2) ne moremo še bolj poenostaviti, če ne poznamo<br />

matematičnih modelov regulatorja, reguliranega sistema <strong>in</strong> senzorja.<br />

Pri poznavanju teh modelov pa ga uporabimo <strong>za</strong> izpeljavo<br />

skupnega modela.<br />

ZGLED 4.2-1 Sistem s povratno <strong>za</strong>nko<br />

Za sistem na sliki 4.2-1 predpostavimo, da lahko regulator <strong>in</strong> senzor modeliramo<br />

s preprostim ojačenjem, torej skaliranjem amplitude (v teoriji sistemov <strong>in</strong> regulacij<br />

imenujemo take sisteme proporcionalni sistemi), regulirani sistem (objekt) pa z<br />

l<strong>in</strong>earno diferencialno enačbo prvega reda oblike y = dv 2<br />

. V tem primeru velja:<br />

dt<br />

y 1(t) = v 2(t) = T 1[v(t)] = K 1e(t)<br />

y 3(t) = T 3[y(t)] = K 3y(t)<br />

y(t) = T 2[v 2(t)] = dv2(t)<br />

dt<br />

kjer sta K 1 <strong>in</strong> K 2 realni konstanti (K 1, K 2 ∈ R), ki določata ojačenje sistema<br />

1 <strong>in</strong> 2. Iz regulacijskega pogreška e(t):<br />

e(t) = v(t) − K 3y(t)<br />

izračunamo regulacijski signal<br />

v 2(t) = K 1e(t) = K 1[v(t) − K 3y(t)] = K 1v(t) − K 1K 3y(t) .<br />

,<br />

Z upoštevanjem modela objekta dobimo <strong>za</strong> izhod sistema:<br />

y(t) =<br />

d dt [K1v(t) − K1K3y(t)] = K1 dv(t)<br />

dt<br />

− K 1K 3<br />

dy(t)<br />

dt<br />

.


84 4. Sistemi<br />

Zgornja enačba je model opazovanega sistema. Združimo še spremenljivke na<br />

eni strani enačaja:<br />

dy(t)<br />

dv(t)<br />

K 1K 3 − y(t) = K 1<br />

dt<br />

dt<br />

<strong>in</strong> dobimo v matematiki običajen <strong>za</strong>pis l<strong>in</strong>earne diferencialne enačbe prvega reda<br />

s konstantnimi koeficienti.<br />

♦<br />

V zgledu 4.2-1 smo večkrat uporabili lastnosti l<strong>in</strong>earnih funkcij.<br />

Te so nam že znane iz obravnave l<strong>in</strong>earnih prostorov (podpoglavje<br />

2.5 na strani 22). Pomen l<strong>in</strong>earnih funkcij pri modeliranju<br />

sistemov bomo opisali v naslednjem razdelku.<br />

4.3 Lastnosti sistemov<br />

V tem podpoglavju se bomo seznanili z osnovnimi lastnostmi<br />

sistemov. Te lastnosti bomo določili na splošno, ne glede na<br />

posamezni matematični model sistema.<br />

Sistem s pomnjenjem Sistem s pomnjenjem imenujemo d<strong>in</strong>amični<br />

sistem. Njegova def<strong>in</strong>icija je:<br />

DEFINICIJA 16 (Sistem s pomnjenjem: d<strong>in</strong>amični sistem)<br />

Če je izhod sistema y(t) odvisen od trenutne <strong>in</strong> predhodnih vrednosti<br />

vhoda v(t), je sistem s spom<strong>in</strong>om, torej d<strong>in</strong>amični sistem.<br />

V primeru, ko je trenutni izhod odvisem le od trenutnega<br />

vhoda, je sistem brez pomnjenja.<br />

Primer d<strong>in</strong>amičnega sistema je <strong>in</strong>tegrator. Def<strong>in</strong>iran je z:<br />

y(t) =<br />

∫ t<br />

−∞<br />

v(τ) dτ ,<br />

torej je izhod <strong>in</strong>tegratorja odvisen od vseh predhodnih vrednosti<br />

vhoda v(t), <strong>za</strong>to ga lahko obravnavamo kot d<strong>in</strong>amični sistem.<br />

Primer vezja s pomnjenjem, ki ga opišemo z <strong>in</strong>tegratorjem, je<br />

konden<strong>za</strong>tor s tokom i kot vhodom <strong>in</strong> napetostjo u kot izhodom:<br />

u(t) = 1 C<br />

∫ t<br />

−∞<br />

i(τ) dτ = q(t)<br />

C<br />

. (4.3-1)<br />

V (4.3-1) smo s q(t) označili naboj, ki se ob pritisnjeni napetosti<br />

na konden<strong>za</strong>tor kopiči v konden<strong>za</strong>torju. Torej je naboj v konden<strong>za</strong>torju<br />

odvisen od vsega toka, ki je do opazovanega trenutka<br />

pritekel v konden<strong>za</strong>tor.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.3 Lastnosti sistemov 85<br />

Sisteme brez pomnjenja imenujemo tudi statični sistemi. Kot<br />

smo že def<strong>in</strong>irali, je pri njih trenutni izhod odvisen le od trenutnega<br />

vhoda. Primer takega sistema je idealni ojačevalnik, ki<br />

smo ga že spoznali. Modelirali smo ga z:<br />

y(t) = Kv(t) t ∈ R .<br />

Sistem brez pomnjenja je tudi upor (spomnimo se Ohmovega <strong>za</strong>kona!)<br />

ali pa seštevalnik dveh <strong>signalov</strong>, množilnik dveh <strong>signalov</strong><br />

itd.<br />

Invertibilnost Invertibilnost je lastnost, pri kateri je med dvema<br />

spremenljivkama enolična pove<strong>za</strong>va. To pomeni, da lahko s poznavanjem<br />

ene enoumno določimo drugo. Pri sistemih to pomeni<br />

naslednje:<br />

DEFINICIJA 17 (Invertiranje)<br />

Sistem je <strong>in</strong>vertibilen, če <strong>za</strong> vsak izhod obstaja le en sam vhod.<br />

Primer <strong>in</strong>vertibilnega sistema je idealni ojačevalnik z ojačenjem<br />

K. Zanj lahko <strong>za</strong>pišemo:<br />

y(t) = Kv(t) ⇒ v(t) = 1 y(t) . (4.3-2)<br />

K<br />

Pri sistemu, kjer je izhod potenca vhoda, <strong>in</strong>vertibilnosti ni. Na<br />

primer:<br />

y(t) = x 2 (t) ⇒ v(t) = ± √ y(t) . (4.3-3)<br />

Če je v (4.3-3) signal napetost, je na primer izhodna napetost 9<br />

V, je vhodna napetost lahko 3 V ali −3 V.<br />

Inverzni sistem Če sistem opravi preslikavo v(t) → y(t), <strong>in</strong>verzni<br />

sistem opravi preslikavo y(t) → v(t), seveda le, če je preslikava<br />

v(t) → y(t) <strong>in</strong>vertibilna. Iz tega sledi, da <strong>za</strong>poredna ve<strong>za</strong>va sistema<br />

<strong>in</strong> njemu <strong>in</strong>verznega sistema tvori tako imenovani identični<br />

sistem. To je sistem, ki ima izhod enak vhodu. Primeri identičnih<br />

sistemov so:<br />

Idealen ojačevalnik z ojačenjem 1.<br />

Zaporedna ve<strong>za</strong>va idealnih ojačevalnikov z ojačenji K <strong>in</strong> 1/K.<br />

V tej ve<strong>za</strong>vi ojačevalnikov drugi ojačevalnik oslabi signal ravno<br />

<strong>za</strong> toliko, kot ga prvi ojači, <strong>za</strong>to je izhodni signal enak vhodnemu.


86 4. Sistemi<br />

itd.<br />

Povedano z drugimi besedami: <strong>in</strong>verzni sistem opravi <strong>in</strong>verzno<br />

transformacijo, <strong>za</strong>to njegovo delovanje simbolično označimo s<br />

T −1 . Za sistem <strong>in</strong> njemu <strong>in</strong>verzni sistem torej velja:<br />

y(t) = T[v(t)] ⇒ v(t) = T −1 [y(t)] . (4.3-4)<br />

Zaporedno ve<strong>za</strong>vo poljubnega sistema <strong>in</strong> njemu <strong>in</strong>verznega sistema<br />

kaže slika 4.3-1.<br />

Slika 4.3-1<br />

Inverzni sistem.<br />

v( t)<br />

<strong>in</strong>verzni<br />

sistem<br />

sistem<br />

y( t)<br />

T T -1<br />

v( t)<br />

Vzročnost Vsi fizikalni sistemi so vzročni. To na kratko pomeni,<br />

da posledica ne more nastati pred svojim vzrokom oziroma<br />

v obliki def<strong>in</strong>icije:<br />

DEFINICIJA 18 (Vzročni ali kav<strong>za</strong>lni sistem)<br />

Sistem je vzročen, ko je trenutna vrednost izhodnega signala sistema<br />

odvisna samo od trenutne <strong>in</strong> preteklih vrednosti vhodnega<br />

signala:<br />

y(t 0 ) = T[v(t t 0 )] . (4.3-5)<br />

D<strong>in</strong>amični sistemi so kav<strong>za</strong>lni sistemi.<br />

Vzročne sisteme imenujemo tudi neanticipativne sisteme. Iz<br />

4.3-5 sledi, da je sistem s preslikavo<br />

y(t) = v(t − 2)<br />

vzročen, saj je njegov trenutni izhod odvisen od vhoda pred<br />

dvema časovnima enotama. Sistem s preslikavo:<br />

y(t) = v(t + 2)<br />

pa ni vzročen, saj predvideva, da je njegov trenutni izhod enak<br />

vrednosti, ki jo bo vhod <strong>za</strong>vzel čez dve časovni enoti. Tak sistem<br />

imenujemo anticipativni sistem. Anticipativni sistem je idealni<br />

prediktor.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.3 Lastnosti sistemov 87<br />

ded Input<br />

utput<br />

Stabilnost S stabilnostjo sistemov se predvsem ukvarja teorija<br />

sistemov, <strong>za</strong>to se bomo tukaj omejili le na def<strong>in</strong>icijo tako imenovane<br />

BIBO stabilnosti:<br />

DEFINICIJA 19 (BIBO stabilni sistem)<br />

Sistem je stabilen, če je pri vsakem amplitudno omejenem vhodu<br />

tudi izhod amplitudno omejen:<br />

|v(t)| M ⇒ |y(t)| = T[v(t)] N , M, N < ∞ . (4.3-6)<br />

Grafično predstavitev pogoja BIBO stabilnosti kaže slika 4.3-2.<br />

v( t)<br />

N<br />

y( t)<br />

N<br />

-M<br />

t<br />

-M<br />

t<br />

Slika 4.3-2<br />

Amplitudno omejeni signali<br />

Primer BIBO stabilnega sistema je identični sistem. Pri njem<br />

je <strong>za</strong>radi v(t) = y(t) tudi M = N. Idealni <strong>in</strong>tegrator, ki ga določa<br />

y(t) =<br />

∫ t<br />

0<br />

v(τ) dτ ,<br />

ni stabilen sistem. Če je na primer vhod enotska stopnica u(t), ki<br />

je omejena funkcija, je izhod strm<strong>in</strong>a y(t) = t, ki pa ni omejena,<br />

saj z naraščanjem t neprestano narašča.<br />

Stabilnost je osnovna lastnost, ki jo <strong>za</strong>htevamo od vseh sistemov.<br />

Nestabilni sistemi so brez praktične vrednosti.<br />

Časovna neodvisnost Časovno neodvisnim ali tudi časovno<br />

<strong>in</strong>variantnim sistemom se lastnosti s časom ne sprem<strong>in</strong>jajo. Def<strong>in</strong>icija<br />

tega se glasi:<br />

DEFINICIJA 20 (Časovno nedvisni sistem)<br />

Sistem je časovno neodvisen, če časovni premik (<strong>za</strong>kasnitev) vhoda<br />

povzroči le časovni premik izhoda:<br />

v(t − t 0 ) → y(t − t 0 ) .


88 4. Sistemi<br />

Časovno neodvisne sisteme imenujemo tudi pomično neodvisne<br />

sisteme. Test časovne neodvisnosti ilustrira slika 4.3-3. Na<br />

Slika 4.3-3<br />

Časovno neodvisni sistem<br />

<strong>za</strong>kasnitev t 0<br />

y( t - t0)<br />

<strong>za</strong>kasnitev t 0<br />

v( t - t0)<br />

sistem<br />

y( t - t0)<br />

v( t) y( t)<br />

sistem<br />

v( t)<br />

njej vidimo, da lahko signal y(t − t 0 ) dobimo ali z <strong>za</strong>kasnitvijo<br />

izhoda <strong>za</strong> t 0 ali z <strong>za</strong>kasnitvijo vhoda <strong>za</strong> t 0 .<br />

4.4 L<strong>in</strong>earni sistemi<br />

L<strong>in</strong>earnost je ena najpomembnejših lastnosti sistemov, ki jih uporabljamo<br />

pri teoriji <strong>signalov</strong>. Uporabili smo jo že pri obravnavi<br />

<strong>signalov</strong>, <strong>za</strong>to tukaj le ponovimo:<br />

DEFINICIJA 21 (L<strong>in</strong>earnost)<br />

Sistem je l<strong>in</strong>earen, če izpolni naslednja kriterija:<br />

1. Aditivnost: če sistem preslika v 1 (t) → y 1 (t) <strong>in</strong> v 2 (t) → y 2 (t),<br />

potem preslika tudi v 1 (t) + v 2 (t) → y 1 (t) + y 2 (t)<br />

2. Homogenost: če sistem preslika v(t) → y(t), potem preslika tudi<br />

αv(t) → αy(t), kjer je α konstanta.<br />

Navedene lastnosti morajo veljati <strong>za</strong> vse v 1 (t), v 2 (t) <strong>in</strong> α.<br />

Kriterija v def<strong>in</strong>iciji vodita do pomembne lastnosti l<strong>in</strong>earnih<br />

sistemov, to je superpozicije. Če sistem preslika v n (t) → y n (t)<br />

potem tudi velja:<br />

α 1 v 1 (t) + α 2 v 2 (t) → α 1 y 1 (t) + α 2 y 2 (t) ,<br />

kjer sta α 1 <strong>in</strong> α 2 konstanti.<br />

superpozicije, l<strong>in</strong>earni sistem.<br />

Torej je sistem, ki izpolni pogoj<br />

4.4.1 Odziv l<strong>in</strong>earnega časovno diskretnega sistema<br />

Opazujmo poseben primer, ko je vhodni signal enotski impulz<br />

∆(n). Sistem, ki je časovno diskretni, l<strong>in</strong>earni, pomično neodvisni<br />

ter relaksiran, se naj nanj odzove s signalom h(n). Predpo-<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.4 L<strong>in</strong>earni sistemi 89<br />

stavimo še, da ta odziv poznamo (slika 4.4-1).<br />

D( n)<br />

1 1<br />

h( n)<br />

0 0<br />

-1 0 1 2 3 t<br />

-1 0 1 2 3<br />

T<br />

v<br />

y<br />

t<br />

= h( n)<br />

Slika 4.4-1<br />

Odziv l<strong>in</strong>earnega, časovno<br />

neodvisnega sistema na<br />

enotski impulz<br />

S term<strong>in</strong>om relaksirani sistem opisujemo sistem, ki ima prazne<br />

pomnilniške elemente:<br />

<strong>in</strong> je izhod sistema enak:<br />

T : ∆(n) → h(n)<br />

y(n) = h(n) .<br />

Zaradi naštetih lastnosti sistema je odziv sistema na <strong>za</strong>kasnjen<br />

enotski impulz enak kot prej, le da je <strong>za</strong>kasnjen <strong>za</strong> toliko,<br />

kot je <strong>za</strong>kasnjen enotski impulz (slika 4.4-2):<br />

na izhodu pa dobimo:<br />

T : ∆(n − m) → h(n − m) ,<br />

y(n) = h(n − m) .<br />

V primeru, da si na vhodu sledita dva enotska impul<strong>za</strong>, na primer<br />

∆(n) <strong>in</strong> ∆(n − m), dobimo na izhodu signal, ki je vsota odzivov<br />

na posamezni enotski impulz (slika 4.4-3):<br />

y(n) = h(n) + h(n − m) .<br />

Z <strong>in</strong>duktivnim sklepanjem ugotovimo, da je izhod obravnavanega<br />

sistema pri <strong>za</strong>poredju enotskih impulzov na vhodu, ki se prične<br />

( n m)<br />

( m 2 )<br />

h( n m)<br />

( m 2)<br />

1 1<br />

0 0<br />

-1 0 1 2 3 t<br />

T<br />

v<br />

-1 0 1 2 3 t<br />

y h( n m)<br />

Slika 4.4-2<br />

Odziv l<strong>in</strong>earnega, časovno<br />

neodvisnega sistema na<br />

<strong>za</strong>kasnjeni enotski impulz.


90 4. Sistemi<br />

D( n)<br />

D( n m)<br />

1 - h( n)<br />

1<br />

h ( n - m )<br />

Slika 4.4-3<br />

Odziv l<strong>in</strong>earnega, časovno<br />

neodvisnega sistema na dva<br />

enotska impul<strong>za</strong>.<br />

0 0<br />

-1 0 1 2 3 t<br />

T<br />

v<br />

1<br />

-1<br />

0<br />

1<br />

2<br />

3<br />

y( n) = h( n) + h( n - m)<br />

t<br />

0<br />

-1<br />

0<br />

1<br />

2<br />

3<br />

t<br />

v trenutku −M <strong>in</strong> konča v trenutku N, enak vsoti odzivov na<br />

posamezne enotske impulze:<br />

y(n) = h(n + M) + h(n + M − 1) + · · · h(n) + · · · h(n − N)<br />

=<br />

N∑<br />

m=−M<br />

h(n − m) .<br />

Pri l<strong>in</strong>earnih sistemih se s skalarjem pomnoženi vhod preslika v<br />

s skalarjem pomnoženi izhod. Torej velja:<br />

<strong>in</strong> <strong>za</strong>to tudi:<br />

α m ∆(n − m) → α m h(n − m)<br />

y(n) =<br />

N∑<br />

m=−M<br />

h(n − m) · α m . (4.4-1)<br />

Sedaj predpostavimo, da imamo na vhodu sistema signal v(m),<br />

<strong>za</strong> katerega velja:<br />

v(m) = α m · ∆(n − m) = α m . (4.4-2)<br />

V (4.4-2) smo upoštevali, da je ∆(n) enotski impulz, torej ima<br />

∆(n − m) v trenutku m vrednost 1 ( v vseh ostalih pa je nič).<br />

Zato lahko (4.4-1) z upoštevanjem (4.4-2) <strong>za</strong>pišemo v obliki:<br />

y(n) =<br />

N∑<br />

m=−M<br />

h(n − m) · v(m) . (4.4-3)<br />

Obrazec (4.4-3) imenujemo l<strong>in</strong>earna konvolucija 1 ali na kratko<br />

kar konvolucija v diskretni obliki. Iz njega sledi, da lahko pri<br />

1 Poznamo tudi krožno konvolucijo, ki pa je tukaj ne obravnavamo.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.4 L<strong>in</strong>earni sistemi 91<br />

poznavanju odziva sistema na enotski impulz (pri časovno neodvisnih<br />

<strong>in</strong> l<strong>in</strong>earnih sistemih) eksplicitno določimo izhod.<br />

Veljavnost (4.4-3) se pri BIBO stabilnih sistemih ohrani tudi<br />

pri neprehodnih signalih, to je pri signalih neskončnega trajanja.<br />

Pri njih velja −M → −∞, N → ∞, konvolucijo pa določa<br />

obrazec:<br />

∞∑<br />

y(n) = h(n − m) · v(m) . (4.4-4)<br />

m=−∞<br />

4.4.2 Odziv l<strong>in</strong>earnega časovno<br />

neodvisnega zveznega sistema<br />

Obrazec <strong>za</strong> konvolucijo pri časovno zveznih, časovno neodvisnih<br />

sistemih lahko izpeljemo z limitnim postopkom, v katerem diskretni<br />

približek zveznega signala prevedemo v zvezni signal. Poglejmo!<br />

Zvezni signal v(t) nadomestimo s približkom, ki ga tvori <strong>za</strong>poredje<br />

impulzov {ˆv(n)}, katerih plošč<strong>in</strong>a je sorazmerna plošč<strong>in</strong>i<br />

pravokotnih pulzov šir<strong>in</strong>e ∆τ <strong>in</strong> viš<strong>in</strong>e v(n∆τ) (slika 4.4-4):<br />

v( t)<br />

v( n<br />

)<br />

<br />

v( n)<br />

-2<br />

-1<br />

0<br />

t<br />

v( n ) ( t n<br />

)<br />

Slika 4.4-4<br />

Pove<strong>za</strong>va med zveznim<br />

signalom <strong>in</strong> njegovim<br />

diskretnim približkom.<br />

-2 -1 0<br />

M a<br />

( n 1)<br />

<br />

n<br />

t<br />

N b<br />

ˆv(n) = v(n∆τ)∆τ · δ(t − n∆τ) . (4.4-5)<br />

To <strong>za</strong>poredje impulzov pošiljamo v trenutkih<br />

a = −M∆τ, . . . , −∆τ, 0, ∆τ, 2∆τ, . . . , n∆τ, . . . , N∆τ = b<br />

v sistem, ki povzroči <strong>za</strong>poredje odzivov:<br />

{v(n∆τ)∆τ · h(t − n∆τ)} ,<br />

Zaradi krajšega pisanja<br />

sta v nadaljevanju <strong>za</strong><br />

l<strong>in</strong>earni časovno<br />

neodvisni zvezni sistem<br />

<strong>in</strong> l<strong>in</strong>earni časovno<br />

neodvisni diskretni<br />

sistem uporabljana<br />

term<strong>in</strong>a zvezni oziroma<br />

diskretni sistem.


92 4. Sistemi<br />

kjer je h(t−n∆τ) odziv zveznega sistema na Diracov impulz δ(t−<br />

n∆τ) (slika 4.4-5). Odziv zveznega sistema v trenutku n∆τ na<br />

Slika 4.4-5<br />

Odziv zveznega sistema na<br />

Diracov impulz.<br />

0 0<br />

-1 0 1 2 3 t<br />

-1 0<br />

T<br />

v<br />

( n m<br />

)<br />

m 2<br />

h( t m<br />

)<br />

1<br />

2<br />

3<br />

y h( n m<br />

)<br />

<br />

( m 2)<br />

t<br />

vhodno <strong>za</strong>poredje pulzov {ˆv(n)} sedaj določimo (podobno kot pri<br />

diskretnih sistemih) s seštevkom trenutnega <strong>in</strong> vseh predhodnih<br />

odzivov:<br />

y(n∆τ) =<br />

=<br />

n∑<br />

r=−M<br />

n∑<br />

r=−M<br />

y τ (n − r)<br />

v[(n − r)∆τ]∆τ · h[t − (n − r)∆τ]<br />

vpeljimo novo oznako n − r = m (iz nje sledi r = n − m):<br />

=<br />

n∑<br />

m=−M<br />

v(m∆τ)∆τ · h[t − m∆τ) . (4.4-6)<br />

Če v limitnem postopku manjšamo ∆τ proti dτ, je v tem primeru<br />

<strong>za</strong> vsako točko na signalni osi določena vrednost v(m∆τ) = v(t),<br />

<strong>za</strong>to n∆τ → t, m∆τ → τ <strong>in</strong> vsota v (4.4-6) preide v <strong>in</strong>tegral:<br />

y(n∆τ }{{}) = lim<br />

t<br />

y(t) =<br />

∆τ→0<br />

n→∞<br />

∫ t<br />

a<br />

n∆τ=t<br />

∑<br />

m∆τ=a<br />

m=m∆τ+∆τ<br />

} {{ }<br />

∫ t<br />

a<br />

v(m∆τ } {{ }) }{{} ∆τ<br />

τ dτ<br />

h(t − m∆τ } {{ })<br />

τ<br />

v(τ)h(t − τ) dτ . (4.4-7)<br />

Integral v (4.4-7) povezuje izhod zveznega sistema z zveznim<br />

vhodnim signalom v(t) <strong>in</strong> s h(t − τ). Imenujemo ga konvolucijski<br />

<strong>in</strong>tegral.<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.4 L<strong>in</strong>earni sistemi 93<br />

4.4.3 Konvolucija ali pregib<br />

Simbolično konvolucijo <strong>za</strong>pišemo z y = h∗v oziroma, kadar želimo<br />

poudariti, ali gre <strong>za</strong> časovno diskretne ali pa <strong>za</strong> časovno zvezne<br />

sisteme z:<br />

y(t) = h(t) ∗ v(t) , t ∈ R y(n) = h(n) ∗ v(n) , n ∈ Z .<br />

Pri izpeljavi konvolucijskih obrazcev smo upoštevali lastnosti l<strong>in</strong>earnih,<br />

časovno neodvisnih sistemov. Pri teh sistemih velja teorem<br />

superpozicije.<br />

Pri časovno spremenljivih sistemih se odziv s časom sprem<strong>in</strong>ja,<br />

<strong>za</strong>to <strong>za</strong>nje ne konvolucijski obrazec ne obstoja. pri zveznih<br />

sistemi namesto njega uporabljamo superpozicijski <strong>in</strong>tegral, v katerem<br />

namesto impulznega odziva h(t−τ) nastopa impulzni odziv<br />

h(t, τ):<br />

y(t) =<br />

∫ t<br />

a<br />

v(τ)h(t, τ) dτ . (4.4-8)<br />

V dosedanji obravnavi smo upoštevali, da je vhodni signal<br />

končnega trajanja. To seveda pri izračunu odziva ni omejitev.<br />

Zato poudarimo:<br />

Če se vhodni signal v(t) prične v določenem trenutku, na primer v<br />

t = 0 <strong>in</strong> še traja, je odziv zveznega sistema določen z:<br />

y(t) =<br />

∫ t<br />

0<br />

h(t − τ)v(τ) dτ . (4.4-9)<br />

Določeni <strong>in</strong>tegrali s spremenljivo zgornjo mejo so zvezne funkcije te<br />

meje. Z drugimi besedami, konvolucija je funkcija, ki opisuje potek<br />

odziva sistema.<br />

ZGLED 4.4-1 Sistem z <strong>in</strong>tegratorjem<br />

Določimo impulzni odziv <strong>in</strong>tegratorja!<br />

Iz obravnave posplošenih <strong>signalov</strong> (podpoglavje 2.10 na strani 49) vemo, da<br />

<strong>za</strong> Diracov impul<strong>za</strong> velja:<br />

{<br />

∫<br />

1 pri t = 0<br />

∞<br />

δ(t) =<br />

<strong>in</strong> x(t)δ(t) dt = x(0) . (2.10-1)<br />

0 sicer<br />

−∞<br />

Iz tega sledi:<br />

y(t) =<br />

=<br />

∫ t<br />

0<br />

∫ t<br />

0<br />

δ(t) dt (4.4-10)<br />

h(t − τ)δ(τ) dτ =<br />

∫ t<br />

0<br />

h(t)δ(τ) dτ . (4.4-11)


94 4. Sistemi<br />

Ker imamo le en vhodni pulz v trenutku nič, je τ = 0. Iz podobnosti <strong>in</strong>tegrala<br />

(4.4-10) z <strong>in</strong>tegralom (4.4-11) <strong>in</strong> iz lastnosti (2.10-1) sledi:<br />

y(t) = h(t) =<br />

∫ t<br />

0<br />

1δ(t) dt =<br />

∫ ∞<br />

−∞<br />

δ(t) dt = 1 , t 0 . (4.4-12)<br />

Torej je impulzni odziv <strong>in</strong>tegratorja enak ena, njegov odziv na Diracov impulz pa<br />

enotska stopnica:<br />

y(t) =<br />

∫ t<br />

0<br />

δ(t) dt = u(t) .<br />

Blokovno shemo <strong>in</strong>tegratorja kaže slika 4.4-6. Še mimogrede: <strong>in</strong>tegrator je zelo<br />

Slika 4.4-6<br />

Primer sistema z<br />

<strong>in</strong>tegratorjem. Zgoraj: sistem<br />

je predstavljen s simbolom <strong>za</strong><br />

<strong>in</strong>tegrator, spodaj: sistem je<br />

predstavljen z impulznim<br />

odzivom.<br />

( t)<br />

t<br />

y( t) z ( t) dt u( t)<br />

<br />

( t)<br />

y( t) = h( t)* ( t) = u( t)<br />

h( t) = 1<br />

pogost gradnik sistemov.<br />

♦<br />

Konvolucijski sistemi Dejansko je odziv sistema na poljubni<br />

vhodni signal superpozicija impulznih odzivov sistema uteženih<br />

z vhodnim signalom. Zato poudarimo:<br />

L<strong>in</strong>earni, časovno diskretni preslikavni sistem s signalno<br />

osjo T ∈ Z <strong>in</strong> realnim ali kompleksnim<br />

vhodom <strong>in</strong> izhodom je časovno neodvisni sistem<br />

takrat <strong>in</strong> samo takrat, če je konvolucijski sistem.<br />

L<strong>in</strong>earni, časovno zvezni preslikavni sistem s signalno<br />

osjo T ∈ R <strong>in</strong> realnim ali kompleksnim<br />

vhodom <strong>in</strong> izhodom je časovno neodvisni sistem<br />

takrat <strong>in</strong> samo takrat, če je konvolucijski sistem.<br />

Vsi sistemi, s katerimi se bomo ukvarjali v tej knjigi, so konvolucijski<br />

sistemi, torej l<strong>in</strong>earni, časovno neodvisni sistemi. Te<br />

sisteme se pogosto označuje z LTI (L<strong>in</strong>ear Time Invariant) sistemi.<br />

LTI sistemi so seveda lahko časovno zvezni ali časovno<br />

diskretni.<br />

Pregib V naslovu podpoglavja smo h konvoluciji <strong>za</strong>pisali besedo<br />

pregib. Čeprav se ta term<strong>in</strong> v slovenšč<strong>in</strong>i le redko uporablja<br />

kot s<strong>in</strong>onim <strong>za</strong> konvolucijo, ga bomo v tem razdelku opravičili.<br />

Pri tem bomo spoznali tudi o<strong>za</strong>dje obrazcev <strong>za</strong> konvolucijo.<br />

Beseda pregib je prevod term<strong>in</strong>ov faltung (nemško) oziroma<br />

convolution (angleško). Term<strong>in</strong>a poudarjata, da pri računanju<br />

konvolucije <strong>za</strong>sučemo, to je prepognemo signalno os. To uvidimo<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.4 L<strong>in</strong>earni sistemi 95<br />

iz naslednje izpeljave. Vpeljimo <strong>za</strong>menjavo spremenljivk:<br />

t − τ = η<br />

<strong>in</strong> predpostavimo, da se je signal pričel v trenutku nič ter da še<br />

traja. Odziv sistema, določa ga (4.4-10), je:<br />

y(t) =<br />

∫ t<br />

0<br />

h(η)v(t − η) dη . (4.4-13)<br />

Opazujmo signal v(t + η)! To je signal v(η), ki je na signalni osi<br />

<strong>za</strong> t premaknjen v levo (slika 4.4-7). Torej del tega signala, ki<br />

v( t)<br />

t<br />

0<br />

<br />

Slika 4.4-7<br />

Potek signala v(t + η).<br />

je levo od izhodišča η = 0, pripada preteklosti funkcije v(t + η).<br />

Odsek (t + η), ki pa je desno od izhodišča η = 0, je prihodnost<br />

opazovanega signala.<br />

Sedaj si oglejmo potek signala v(t − η), ki se nahaja pod <strong>in</strong>tegralom<br />

v (4.4-13)! Signal v(t − η) nastane iz signala v(t + η)<br />

tako, da η spremenimo predznak. To pomeni, da signal preganemo<br />

tako, da je njegov <strong>za</strong>četek v točki t (slika 4.4-8).<br />

v( t)<br />

<br />

0<br />

h( )<br />

t<br />

v( t)<br />

<br />

Slika 4.4-8<br />

Signal v(t − η) <strong>in</strong> odziv h(η)<br />

<strong>in</strong><br />

Glede na to, da obravnavamo vzročni sistem, velja:<br />

h(η) = 0 če je η < 0<br />

v(t − η) = 0 , η > t + a ,<br />

smemo raztegniti meji <strong>in</strong>tegrala na vso realno os. Vzporedno<br />

<strong>za</strong>pišimo konvolucijska obrazca <strong>za</strong> časovno diskretne <strong>in</strong> <strong>za</strong> časovno<br />

zvezne signale <strong>in</strong> sisteme:


96 4. Sistemi<br />

Časovno zvezni signali:<br />

∫ ∞<br />

y(t) = h(τ)v(t − τ) dτ (4.4-14)<br />

oziroma:<br />

−∞<br />

∫ ∞<br />

y(t) = h(t − τ)v(τ) dτ . (4.4-15)<br />

−∞<br />

Časovno diskretni signali:<br />

∞∑<br />

y(n) = h(n)v(n − m) (4.4-16)<br />

m=−∞<br />

oziroma:<br />

∞∑<br />

y(n) = h(n − m)v(m) . (4.4-17)<br />

m=−∞<br />

Iz <strong>in</strong>tegrala v (4.4-14) <strong>in</strong> s slike 4.4-8 lahko razberemo, da<br />

prispeva sedanji vrednosti izhoda y(t) vsa preteklost vhodnega<br />

signala v(t) uteženega z impulznim odzivom h(t), ki pripada sistemukonvolucija<br />

je<br />

Za dobro razumevanje konvolucije izpišimo še algoritem nje-<br />

komutativna operacija nega računanja pri zveznih signalih. Sledita algoritma <strong>za</strong> primer,<br />

ko <strong>za</strong>sučemo vhodni signal (na levi strani) <strong>in</strong> <strong>za</strong> primer, ko<br />

<strong>za</strong>sučemo impulzni odziv. Oba algoritma imata enak rezultat!<br />

1. preganemo vhodni signal:<br />

v(t) → v(−t),<br />

2. preganjeni signal premikamo po signalni<br />

osi <strong>za</strong> τ,<br />

3. množimo (točko <strong>za</strong> točko) vhodni signal z<br />

impulznim odzivom:<br />

h(t)v(τ − t),<br />

4. zmnožke <strong>in</strong>tegriramo, glej (4.4-14).<br />

1. preganemo impulzni odziv:<br />

h(t) → h(−t),<br />

2. preganjeni impulzni odziv premikamo po<br />

signalni osi <strong>za</strong> τ,<br />

3. množimo (točko <strong>za</strong> točko) vhodni signal z<br />

impulznim odzivom:<br />

h(τ − t)v(t),<br />

4. zmnožke <strong>in</strong>tegriramo, glej (4.4-15).<br />

Naštete korake ponovimo <strong>za</strong> vsak τ. Postopek konvolucije<br />

si lahko predstavljamo kot lokalno povprečenje vhodnega signala<br />

z <strong>za</strong>sukanim oziroma preganjenim <strong>in</strong> premikajočim se impulznim<br />

odzivom (slika 4.4-9). Ker pri računanju izhodnega signala uporabimo<br />

pregib impulznega odziva, opravičujemo ime pregib (nemško<br />

faltung, angleško convolution) <strong>za</strong> to operacijo. To ime so vpeljali<br />

nemški matematiki, ki so izpeljali konvolucijske obrazce. V nadaljevanju<br />

teksta bomo <strong>za</strong> te obrazce uporabljali udomačeno ime<br />

konvolucija.<br />

Iz obrazcev konvolucije <strong>in</strong> s slike 4.4-9 vidimo, da je konvolu-<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.4 L<strong>in</strong>earni sistemi 97<br />

h( t)<br />

h( t)<br />

0 t<br />

0<br />

t<br />

v( t)<br />

0<br />

a<br />

t<br />

<br />

0<br />

a<br />

t<br />

<br />

0<br />

a<br />

t<br />

Slika 4.4-9<br />

Postopek računanja<br />

konvolucije<br />

<br />

0<br />

a<br />

t<br />

<br />

y( t)<br />

0<br />

a<br />

t<br />

0<br />

a<br />

t<br />

cija po postopku računanja zelo podobna korelaciji. Razlika med<br />

njima je le v <strong>za</strong>suku enega od <strong>signalov</strong>, <strong>za</strong>to imata korelacija <strong>in</strong><br />

konvolucija <strong>signalov</strong>, kjer je signal h sodi, enak rezultat!


98 4. Sistemi<br />

4.4.4 Impulzni odziv<br />

Funkcijo h, ki je odziv sistema na vzbujanje z enotskim impulzom,<br />

imenujemo impulzni odziv. V diskretnih sistemih je to odziv<br />

na Kroneckerjev impulz ∆, pri zveznih pa na Diracov impulz δ.<br />

Uč<strong>in</strong>ek konvolucije na izhod sistema lahko simbolično prikažemo<br />

z:<br />

v −−−−−−−−−−−−→<br />

konvolucija s h<br />

y , ∆ ali δ −−−−−−−−−−−−→<br />

konvolucija s h h .<br />

Iz impulznega odziva lahko na primer takoj vidimo, ali je sistem<br />

vzročen ali ne. Če je vzročen, potem seveda velja že znani pogoj:<br />

h(t) = 0 , t < 0 . (4.4-18)<br />

Nadalje lahko vidimo, ali je sistem stabilen. Za BIBO stabilnost<br />

konvolucijskega sistema lahko def<strong>in</strong>iramo:<br />

DEFINICIJA 22 (BIBO stabilnost konvolucijskih sistemov)<br />

Časovno diskretni ali časovno zvezni konvolucijski sistem je BIBO<br />

stabilen, če je odziv sistema y na vsak vhod v s končno amplitudo<br />

končno velik. Torej,<br />

če je ‖v‖ ∞ < ∞, potem je ‖y‖ ∞ < ∞<br />

Ker je odziv konvolucijskega sistema določen z impulznim odzivom<br />

h, seveda ne preseneča, da lahko stabilnost takega sistema<br />

preverimo s proučitvijo impulznega odziva. Iz zgornje def<strong>in</strong>icije<br />

BIBO stabilnosti lahko pov<strong>za</strong>memo, da je konvolucijski sistem<br />

BIBO stabilen takrat (<strong>in</strong> samo takrat), ko ima impulzni odziv h<br />

končno vitalnost: ‖h‖ 1 < ∞.<br />

□ □ □ □ <br />

<strong>za</strong>htevna snov<br />

Doka<strong>za</strong>li bomo, da v primeru, ko je konvolucijski sistem BIBO<br />

stabilen, velja <strong>za</strong> odziv sistema naslednje:<br />

‖y‖ ∞ ‖h‖ 1 · ‖v‖ ∞ (4.4-19)<br />

pri vsakem vhodu s končno amplitudo. Enakost v (4.4-19) lahko vedno<br />

dosežemo. Na primer, če velja:<br />

⎧<br />

⎨ h ∗ (−t)<br />

pri h(−t) ≠ 0<br />

v(t) = |h(−t)|<br />

⎩<br />

0 drugače<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.4 L<strong>in</strong>earni sistemi 99<br />

Veljavnost (4.4-19) bomo doka<strong>za</strong>li le pri časovno diskretnih sistemih.<br />

Za časovno zvezne sisteme je izpeljava doka<strong>za</strong> ekvivalentna.<br />

DOKAZ<br />

Izhod časovno diskretnega konvolucijskega sistema y(n) je:<br />

y(n) =<br />

n∑<br />

m=−∞<br />

h(n − m)v(m) , n ∈ Z<br />

Z uporabo trikotniškega izreka <strong>za</strong> kompleksna števila lahko <strong>za</strong>pišemo:<br />

n∑<br />

|y(n)| =<br />

h(n − m)v(m)<br />

∣<br />

∣ <br />

∑ ∞<br />

|h(n − m)v(m)|<br />

m=−∞<br />

m=−∞<br />

( n<br />

)<br />

∑<br />

= |h(n − m)| ‖v‖ ∞ , n ∈ Z<br />

m=−∞<br />

Pri tem smo upoštevali, da je<br />

|v(m)| sup |v(m)| = ‖v‖ ∞ (4.4-20)<br />

m<br />

pri vseh m ∈ Z. Z <strong>za</strong>menjavo spremeljivk n − m = k dobimo:<br />

( n<br />

)<br />

∑<br />

|y(n)| |h(k)| ‖v‖ ∞ = ‖h‖ 1 · ‖v‖ ∞ , n ∈ Z<br />

m=−∞<br />

kar z upoštevanjem (4.4-20) da želeno oceno: ‖y(t)‖ ‖h‖ 1 · ‖v‖ ∞ . □<br />

Iz doka<strong>za</strong> lahko pov<strong>za</strong>memo:<br />

če je vitalnost impulznega odziva ‖h‖ 1 , končna, potem je maksimalna<br />

amplituda odziva ‖y‖ ∞ na vsak vhod s končno amplitudo<br />

‖v‖ ∞ ∞ omejena s ‖h‖ 1 · ‖v‖ ∞ .<br />

□ □ □ □ <br />

V blokovnih shemah konvolucijske sisteme ponavadi označimo<br />

s pripadajočim impulznim odzivom (slika 4.4-10).<br />

v( t) y( t) = h( t v t v( n)<br />

h( t) )* ( ) y( n) = h( n<br />

h( n)<br />

)* v( n)<br />

Slika 4.4-10<br />

Levo: zvezni sistem, desno:<br />

diskretni sistem<br />

Zaključimo s poudarkom:<br />

Impulzni odziv h je osnovni podatek o sistemu.<br />

obnaša nek sistem.<br />

Pove, kako se


100 4. Sistemi<br />

4.5 Lastnosti konvolucije<br />

Konvolucija je l<strong>in</strong>earna operacija, <strong>za</strong>to <strong>za</strong>njo veljajo vse lastnosti<br />

l<strong>in</strong>earnih funkcij. Njene lastnosti si oglejmo na primeru dveh<br />

splošnih <strong>signalov</strong> x <strong>in</strong> y.<br />

Komutativnost Če obstaja x ∗ y, potem velja tudi:<br />

x ∗ y = y ∗ x .<br />

Asociativnost<br />

Če obstaja (x ∗ y) ∗ z, velja tudi:<br />

(x ∗ y) ∗ z = y ∗ (x ∗ z) .<br />

Distributivnost<br />

Če obstajata x ∗ y <strong>in</strong> x ∗ z, velja tudi:<br />

x ∗ (y + z) = x ∗ y + x ∗ z .<br />

Če obstaja x∗y, po-<br />

Komutativnost množenja s skalarjem<br />

tem <strong>za</strong> vsak skalar α ∈ C velja:<br />

α(x ∗ y) = (αx) ∗ y = x ∗ (αy) .<br />

Premik funkcije σ t naj bo operator pomika funkcije po časovni<br />

osi (signala po signalni osi) <strong>za</strong> čas t, kjer je n ∈ Z pri diskretnih<br />

funkcijah <strong>in</strong> t ∈ R pri zveznih funkcijah. Če obstaja konvolucija<br />

x ∗ y, potem obstaja tudi <strong>za</strong> premaknjene funkcije:<br />

σ t (x ∗ y) = (σ t x) ∗ y = x ∗ (σ t y) .<br />

Lastnosti odvajanja<br />

D naj bo operator odvajanja:<br />

Dz(t) = dz(t)<br />

dt<br />

,<br />

z(t) je zvezna, odvedljiva funkcija, potem pri konvoluciji zveznih<br />

<strong>signalov</strong>, od katerih mora biti vsaj eden odvedljiv, velja:<br />

če je odvedljiv signal x, <strong>in</strong><br />

če je odvedljiv signal y.<br />

D(x ∗ y) = (Dx) ∗ y ,<br />

D(x ∗ y) = x ∗ (Dy) ,<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.5 Lastnosti konvolucije 101<br />

4.5.1 Konvolucija signala z Diracovim impulzom<br />

Ponavadi se v učbenikih <strong>za</strong> teorijo sistemov <strong>in</strong> tudi <strong>za</strong> <strong>obdelavo</strong><br />

<strong>signalov</strong> konvolucija časovno zveznih <strong>signalov</strong> izpelje s pomočjo<br />

lastnosti, ki jo ima Diracov impulz. Mi smo v tej knjigi ubrali<br />

nekoliko drugačno pot, <strong>za</strong>to si posebej poglejmo konvolucijo signala<br />

z Diracovim impulzom. Ta pristop nam bo koristil kasneje,<br />

ko bomo v nadaljnjih knjižicah o obdelavi signala obravnavali<br />

transformacije <strong>signalov</strong> <strong>in</strong> operacije nad njimi.<br />

Iz razdelka 2.10.1 na strani 51 vemo, da <strong>in</strong>tegral produkta<br />

časovno zveznega signala x(t) z Diracovim impulzom δ(t − τ) da<br />

vrednost signala v trenutku τ:<br />

∫ ∞<br />

−∞<br />

x(t) · δ(t − τ) dt = x(τ) (4.5-1)<br />

Ta lastnost izhaja iz def<strong>in</strong>icije δ(t) <strong>in</strong> jo uporabljamo pri vzorčenju<br />

<strong>signalov</strong> (o tem bo govora v knjižici o digitalni obdelavi <strong>signalov</strong>)<br />

<strong>in</strong> pri računanju konvolucije. Če v konvolucijskem <strong>in</strong>tegralu vhodni<br />

signal <strong>za</strong>menjamo z Diracovim impulzom, iz def<strong>in</strong>icije konvolucije<br />

sledi:<br />

∫ t<br />

−∞<br />

kar simbolično <strong>za</strong>pišemo z:<br />

x(t − τ)δ(τ) dτ = x(t) ,<br />

x(t) ∗ δ(t) = x(t) ,<br />

oziroma pri upoštevanju (4.5-1) dobimo (slika 4.5-1):<br />

∫ t<br />

−∞<br />

x(t − τ)δ(t − τ) dτ = x(t − τ) ,<br />

x( t)<br />

- <br />

0<br />

( ) ( t - )<br />

t<br />

kar simbolično <strong>za</strong>pišemo z:<br />

x(t) ∗ δ(t − τ) = x(t − τ) .<br />

y( t)<br />

0<br />

<br />

y( t) = x( t - )<br />

t<br />

Zgornje enačbe so zgovorne same po sebi, <strong>za</strong>to <strong>za</strong>ključimo samo<br />

z naslednjim poudarkom:<br />

Konvolucija signala z Diracovim impulzom da signal.<br />

Konvolucija signala s časovno premaknjenim Diracovim impulzom<br />

premakne signal <strong>za</strong> <strong>za</strong>mik Diracovega impul<strong>za</strong>.<br />

0<br />

Slika 4.5-1<br />

t


102 4. Sistemi<br />

4.5.2 Konvolucija signala z odvodom Diracovega<br />

impul<strong>za</strong><br />

Nadaljujmo s pregledom lastnosti konvolucije z odvodom Diracovega<br />

impul<strong>za</strong>. Predpostavimo, da je x n-krat odvedljiva regularna<br />

funkcija z n ∈ Z + . Z uporabo lastnosti:<br />

∫ ∞<br />

−∞<br />

δ (n) (t)φ(t) dt = (−1) n φ (n) (0) ,<br />

kjer je n ∈ Z + <strong>in</strong> φ najmanj n-krat v točki 0 (zvezno) odvedljiva<br />

funkcija, dobimo:<br />

(<br />

x ∗ δ (n)) ∫ ∞<br />

= x(t − τ)δ (n) (τ) dτ (4.5-2)<br />

−∞<br />

kjer je x (n) n-ti odvod funkcije x. Sledi<br />

∣<br />

(n) d(n)<br />

= (−1)<br />

dτ n x(t − τ) ∣∣τ=0<br />

(4.5-3)<br />

= x (n) (t) ; t ∈ R , (4.5-4)<br />

x ∗ δ (n) = x (n) , m, n ∈ Z . (4.5-5)<br />

Lastnosti konvolucije z δ impulzom so povzete v tabeli 4.1.<br />

Tabela 4.1<br />

Lastnosti konvolucije z Diracovim δ impulzom<br />

lastnost<br />

pogoj<br />

1 x ∗ δ (n) = x (n) x je katerakoli regularna ali<br />

generalizirana funkcija, n ∈ Z +<br />

2 δ (m) ∗ δ (n) = δ (m+n) m, n ∈ Z +<br />

ZGLED 4.5-1 Sistem z diferenciatorjem<br />

Določimo impulzni odziv diferenciatorja!<br />

Diferenciator je zvezni vhodno-izhodni sistem s pove<strong>za</strong>vo y = Dv, D je<br />

operator odvajanja. Polni <strong>za</strong>pis odvajanja je:<br />

y(t) = dv(t) , t ∈ R .<br />

dt<br />

Preprosto lahko preverimo, da je ta sistem l<strong>in</strong>earen <strong>in</strong> časovno neodvisen, torej<br />

je konvolucijski sistem. Njegov impulzni odziv poiščimo tako, da damo na vhod<br />

Diraciv impulz δ(t). Vemo, da je v tem primeru odziv sistema y(t) = h(t):<br />

h(t) = Dδ(t) = δ (1) (t) .<br />

Torej je impulzni odziv diferenciatorja enak δ (1) . Ker je δ (1) generalizirana funkcija<br />

oziroma distribucija, mnogo avtorji pravijo, da diferenciator kot sistem ne<br />

obstaja, saj v praksi tako δ kot tudi δ (1) ne znamo realizirati.<br />

♦<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


4.6 Obstoj konvolucije 103<br />

4.6 Obstoj konvolucije<br />

Konvolucijo smo v (4.4-14) <strong>in</strong> (4.4-15), oziroma v (4.4-16) <strong>in</strong><br />

(4.4-17) def<strong>in</strong>irani z neskončno vrsto ali <strong>in</strong>tegralom. Njen obstoj,<br />

to pomeni izračunljivost, ni sam po sebi umeven. Tako na primer,<br />

konvolucija dveh <strong>signalov</strong>, ki imata konstantno vrednost,<br />

divergira oziroma narašča preko vseh mej.<br />

Zadostni pogoji <strong>za</strong> obstoj konvolucije so naslednji:<br />

1. če sta funkciji x ali y ali obe omejenega trajanja,<br />

2. če je evklidska norma <strong>signalov</strong> x <strong>in</strong> y končna:<br />

‖x‖ 2 < ∞, ‖y‖ 2 < ∞<br />

(pri tem ni nujno, da obstaja tudi evklidska norma konvolucije:<br />

‖x ∗ y‖ 2 ), ali<br />

3. če obstaja ‖x‖ 1 <strong>in</strong> ‖y‖ 1 .<br />

Z omejenim trajanjem funkcije povemo, da je funkcija različna<br />

od nič le na omejenem območju def<strong>in</strong>icijske domene. Na primer,<br />

da je različna od nič le nad delom osi T Z pri diskretnih signalih<br />

ali T R pri zveznih signalih, oziroma da obstajata meji<br />

−∞ < M N < ∞, zunaj katerih je funkcija enaka nič (slika<br />

4.6-1).<br />

x( t) x( n)<br />

- M 0 N t<br />

0 n<br />

Slika 4.6-1<br />

Funkciji omejenega trajanja.<br />

Levo: x(t), −M t N,<br />

desno: x(n), n 0.<br />

Iz def<strong>in</strong>icije konvolucije x ∗ y tudi sledi, da je je rezultat konvolucije<br />

lahko različen od nič nad <strong>in</strong>tervalom, katerega dolž<strong>in</strong>a je<br />

vsota <strong>in</strong>tervalov, nad katerimi sta funkciji x <strong>in</strong> y različni od nič.<br />

4.7 Zaključek<br />

V tem poglavju smo podali elementarni pregled sistemov. Bolj<br />

podrobna <strong>in</strong> bolj matematično usmerjena izpeljana znanja so v<br />

učbenikih, ki se ukvarjajo s teorijo sistemov, na primer [17].


104 4. Sistemi<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


Seznam oznak<br />

operacije<br />

〈 , 〉 skalarni ali notranji produkt<br />

x(t)<br />

na primer 〈a, b〉 = ab s<strong>in</strong> α α : kot med a <strong>in</strong> b, a, b iznos a, b.<br />

povprečna vrednost x(t), na primer x(t) = 1 ∫ T /2<br />

T −T /2 x(t) dt = 1 T<br />

spremenljivke<br />

f frekvenca, f = 1/T<br />

h impulzni odziv<br />

h(t) impulzni odziv zveznega l<strong>in</strong>earnega sistema<br />

h(n) impulzni odziv diskretni l<strong>in</strong>earnega sistema<br />

p(t) trenutna moč<br />

r(τ) korelacija<br />

r x (τ): avtokorelacija signala x(t)<br />

r xy (τ): križna korelacija signala x(t) s signalom y(t)<br />

u(t) enotska stopnica, u(t) = 1, t 0<br />

v vhodni signal<br />

x signal<br />

x(t) zvezni časovni signal<br />

x(n) diskretni časovni signal<br />

x q amplitudno diskretni signal<br />

x q (t) amplitudno diskretni zvezni časovni signal<br />

x q (n) amplitudno diskretni časovni signal<br />

∫<br />

T<br />

x(t) dt<br />

105


106 4. Sistemi<br />

x ∗<br />

y<br />

konjugirano kompleksni signal<br />

če je x(t) = a(t) + jb(t), potem x ∗ (t) = a(t) − jb(t)<br />

izhodni signal<br />

E<br />

P<br />

T<br />

δ(t)<br />

∆(n)<br />

φ(t)<br />

ω<br />

A<br />

C<br />

N<br />

R<br />

T<br />

Z<br />

energija<br />

povprečna moč, P = 1 ∫<br />

T T<br />

p(t) dt<br />

perioda (periodičnega) signala, x(t) = x(t + T ), T = 1 f<br />

Diracov impulz<br />

enotski ali Kroneckerjev impulz<br />

bazična funkcija<br />

krožna frekvenca, ω = 2πf = 2π<br />

T<br />

množice<br />

amplitudno območje signala, amplitudni razmah (<strong>in</strong>terval realnih števil)<br />

prostor kompleksnih števil<br />

C N : prostor kompleksnih N-teric<br />

prostor naravnih števil<br />

prostor realnih števil<br />

R N : prostor realnih N-teric<br />

signalna os<br />

prostor celih števil<br />

Z − : nepozitivna (sem<strong>in</strong>egativna) cela števila<br />

Z + : nenegativna (semipozitivna) cela števila<br />

transformacije<br />

T transformacija, preslikava signala<br />

T −1 <strong>in</strong>verzna transformacija<br />

posebni signali<br />

Sa s<strong>in</strong> x/x, ”sample function”<br />

ramp(t) strm<strong>in</strong>a<br />

rect(t) pravokotni pulz<br />

trian(t) trikotni pulz<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328


Literatura<br />

[1] A.V. Oppenheim, Ronald W. Schafer (1989). Discrete-Time<br />

Signal Process<strong>in</strong>g. Prentice Hall Process<strong>in</strong>g Series, ISBN<br />

0-13-216292-X<br />

[2] E.C. Ifeachor, B.W. Jervis (1997). Digital signal proces<strong>in</strong>g, A<br />

practical approach. Addison-Wesley, ISABN 0-201-54413-X<br />

[3] H. S. Carslaw (1930). An <strong>in</strong>troduction to Fourier’s series and<br />

<strong>in</strong>tegrals (third edition). Dover Publications, <strong>in</strong>c. (ponatis 1960)<br />

[4] I.N. Sneddon (1951). Fourier transforms. Dover publications<br />

I<strong>in</strong>c., (ponatis 1995), ISBN 0-486-68522-5 (pbk)<br />

[5] H.F. Davis (1963). Fourier series and orthogonal functions.<br />

Dover publications Inc., IISBN 0-486-65973-9<br />

[6] M. Reed, B. Simon (1975). Fourier analysis, Self-Adja<strong>in</strong>tness.<br />

Academic press Inc., ISBN 0-12-585002-6(v.2)<br />

[7] H. Kwakernaak, R. Sivan (1991). An <strong>in</strong>troduction to the analysis<br />

and process<strong>in</strong>g of signals (third edition). The McMillan Press<br />

LTD., ISBN 0-333-48887-3<br />

[8] M.R. Spiegel (1974). Theory and problems of Fourier analysis<br />

with applications to boundary value problems. Schaum’s outl<strong>in</strong>e<br />

series, McGraw-Hill (18.izdaja 1994). ISBN 0-07-060219-0<br />

[9] M.R. Spiegel (1974). Theory and problems of Lapalace transform.<br />

Schaum’s outl<strong>in</strong>e series, McGraw-Hill (18.izdaja 1994). ISBN<br />

0-07-06231-X<br />

[10] M.E. van Valkeburg (19xx). Network Analysis.<br />

[11] D. Lange (19xx). Methoden der Signal und sistemanalise.<br />

[12] Dietmar Achilles (1985). Die Fourier-Transformation <strong>in</strong> der<br />

Signalverabeitung. Spr<strong>in</strong>ger Verlag<br />

107


108 LITERATURA<br />

[13] Charles K. Chui, Guanrong Chen (1992). Signal Process<strong>in</strong>g and<br />

System Theory (Selected topics). Spr<strong>in</strong>ger Verlag<br />

[14] Paul A. Lynn (1994). An <strong>in</strong>troduction to the analysis and<br />

Process<strong>in</strong>g of signals. MacMillan Press LTD.<br />

[15] I.N. Bronšte<strong>in</strong>, K.A. Semendjajev, G. Musol, H. Mühlig (1997).<br />

Matematični priročnik. Tehniška <strong>za</strong>ložba Ljubljana.<br />

[16] Ludvig Gyergyek (1987). Teorija obdelave <strong>signalov</strong> <strong>in</strong> statistične<br />

metode. Založba FER Ljubljana.<br />

[17] R. Svečko (1994). Teorija l<strong>in</strong>earnih sistemov. Založba FERI<br />

Maribor ISBN 86-435-0076-3<br />

[18] Dali -Donlagić (1992). Teorija obdelave <strong>signalov</strong>. (<strong>za</strong>piski<br />

predavanj) TF - ERI Maribor.<br />

[19] Žarko Čučej, Peter Plan<strong>in</strong>šič (1999). Teorija <strong>signalov</strong>: Uvod v<br />

teorijo, Založba TF - ERI Maribor, ISBN 86-435-0267-7<br />

http://SPaRC.feri.uni-mb/publikacije<br />

[20] Charles L. Phillips, John M. Parr (1995). Signals, systems, and<br />

transformas, Prentice Hall Inc., ISBN 0-13-795253-8<br />

[21] John J. Komo (1987). Signals, systems, and transformas,<br />

Academic pres <strong>in</strong>c., ISBN 0-12-418660-2<br />

Žarko Čučej: Teorija <strong>signalov</strong> revizija 20010328

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!