ako pronađeš grešku napravi PR na: https://github.com/ciglenecki/fer-oos
2024-02-22: dodan minus predznak kod "Radi preglednosti uvodimo oznaku"
Sinusoida/kosinusoida:
Reprezetnacija sinunsoide:
Eksponencijalna funkcija
Veza između diskretne i kont. eksponencijale:
Veza između Kroneckove delta funkcije i jedinične stepenice: sumacija
od
Signal rastavljamo u njegove gradivne komponente koje se sastoje od
drugih signala (bazne funkcije
Rastav signala na bazne funckije i njihov spektar može sadržavati prebrojivo ili neprebrojivo mnogo baznih funckija.
Najpoznatnija metoda rastava signala je Fourierova transformacija.
Skalarnim umnoškom možemo pronaći spektar
Primjetiti da je
gdje:
Operacija
Operacija
Fouriervoa baza (
Spektar (DTFN):
Signal (IDFTN):
Očitavanje/uzorkovanje/otpikavanje kontinuiranog signala
Signalu
Ako uz zadovoljen uvjet teorema očitavanja uzorci signala jednoznačno opisuju signal onda to znači da mora postojati postupak rekonstrukcije signala iz njegovih uzoraka. Taj postupak zovemo interpolacija.
Međutim, želimo vremenski kont. signal
Mora biti
Očitavanje na mjestima na kojima želimo radimo sa Diracovom delta funkcijom.
Za jednostavniju interpretaciju koristimo Diracov češalj (periodizirane Diracove delta funkcije)
Sad originalan signal
Što je sa spektrom? Kako se on očitava?
Spektar
što je periodizacija spektra
Ako uz zadovoljen (za sada nedefinirani) uvjet teorema očitavanja uzorci signala jednoznačno opisuju signal onda to znači da mora postojati postupak rekonstrukcije signala iz njegovih uzoraka. Taj postupak zovemo interpolacija.
Za rekonstrukciju u spektralnoj domeni potrebno je
zadovoljiti Nyquistov uvjet. Za pojasno ograničene signale maksimalne
frekvencije
Jednostavnije: frekvencije očitavanja (
Ako je Nyquistov uvjet zadovoljen iz pozantog spektra
Taj postupak se zove filtracija a matematički se može zapisati kao:
gdje je spektar
Konvolucija
Konvolucija u vremenskoj domeni rezultira množenjem u spektralnoj domeni.
Sad se mogu dva spektra pomnožiti (spektar se definira kao integral
originalnog signala pomnožen sa
Sad znamo da spektralna filtracija postaje konvolucija u vremenskoj domeni:
Dodatno,
Shannonova interpolacijska formula ili jednostavnije idealnom sinc interpolacijom.
Formula nam u osnovi govori da kontinuirani signal rekonstruiramo tako da na mjesto svakog uzorka postavimo funkciju sinc skaliranu tako da njena amplituda bude upravo jednaka amplitudi uzorka.
Ulazne veličine
Sustav prema tome možemo opisati tako da modeliramo kako ulazne veličine utječu na izlazne veličine.
Sustav
Primjeri sustava:
Malo šire o svojstvima sustava:
Sustav je vremenski nepromjenjiv ako za svaki ulazno-izlazni par
signala
U praksi smjer pomaka prema budućnosti je nemoguć jer to znači da ovisimo o budućim uzorcima signala.
Svojstvo linearnosti sastoji se od:
a sustavi koji su linearni zadovoljavaju svojstvo superpozicije.
Linearno vremenski nepromjenjivi (LTI) sustavi su oni koji zadovoljavaju svojstvo:
Promatramo dva ulazna signala
ako je odziv
Sustav je homogen ako ulazni signal pomnožimo sa
Sustav je linearan ako je aditivan i homogen za svaki ulazno-izlazni par:
Prije svega, želimo predstaviti sustav
Međutim želimo ograničiti broj uzoraka koji je potreban za računanje odziva na neku konačnu vrijednost (računala…)
Dva svojstva koja ograničavaju broj uzoraka su:
Sustav je bezmemorijski ako računamo odziv sustava u koraku
Sustav je memorijski ako za računanje odziva sustava u koraku
Primjer memorijskog sustava prvog reda. Ovaj sustav ovisi o samo
jednom dodatnom uzorku osim
Kauzalnost u općem smislu jest veza koja opisuje kako neki događaj kojeg zovemo uzrok utječe na neki drugi događaj kojeg zovemo posljedica.
U obradbi signala kauzalnost sustava znači da odziv sustava ne smije ovisiti o pobudi sustava koja se još nije dogodila.
Odziv
Formalno (Sustav je kauzalan ako ne postoji budući korak
Nekauzalan sustav je sustav koji nije kauzalan. Dakle, sustav čija
pobuda
Impulsni odziv sustava je skup odziva koji su pomaknuti sa jediničnim
impulsom
Skup svih
Bitno je napomenuti da je
Odziv na
Sad možemo ispustiti
Sad za linearni vremenski nepromjenjivi sustav dobivamo konvoluciju:
Pomoću konvolucije možemo opisati ponašanje svakog linearnog vremenski nepromjenjivog sustava.
Rezime:
Linearan vremenski nepromjenjiv sustav S opisan svojim impulsnim
odzivom
Konvolucija za kauzalne signale za koje vrijedi
Dodatan način reprezentacije LTI sustava
Krećemo od Z transformacija koje djeluje na konvolucije:
presjek njihovih područja konvergencije je određen s
ako su svi signali kauzalni rasta ne većeg od eksponencijalnog onda
Kako Z transformacija konvoluciju pretvara u množak imamo:
Obzirom da impulsni odziv
Sad imamo tri glavne zadaće a mi promatramo uglavnom konvolucijsku zadaću:
Odziv LIT sustava na diskretnu eksponencijalu
gdje je
Ukoliko zamjenimo Z transformaciju sa DTFT-om dobivamo:
DTFT transformacija impulznog odziva
ako RoC prijenosne funkcije H(z) obuhvaća jediničnu kružnicu.
Sustav je (BIBO, bounded-input, bounded-output) stabilan ako za svaku amplitudno ograničenu pobudu daje amplitudno ograničeni odziv.
Signal
Dakle sustav
Neka je
LTI sustav određen impulsnim odzivom
ako je LTI sustav stabilan onda područje konvergencije prijenosne
funckije
Područje konvergencije stabilnog LTI sustava mora sadržavati
jediničnu kružnicu pa je
Vremenski diskretan LTI sustav opisan linearnom diferencijskom jednadžbom sa stalnim koeficijentima je stabilan ako se svi polovi njegove prijenosne funkcije nalaze unutar jedinične kružnice.
Filtracija je postupak u kojem se dizajniraju sustavi koji propuštaju i/ili prigušavaju odredene frekvencije signala. To je postupak uklanjanja neželjenih značajki ili komponenta signala. Najčešće to znači uklanjanje nekih frekvencija ili frekvencijskih pojaseva.
Dizajniranje filtra znači odrediti koeficijente prijenosne funkcije digitalnog filtra koje zadovoljavaju zadanu filtarsku specifikaciju. Tu zadaću pronalaženja filtarskih koeficijenata nazivamo računalnom podržani dizajn filtara. Glavni problem kod dizajniranja filtra jest određivanje koeficijenata prijenosne funkcije:
takvih da frekvencijska karakteristika filtra
za kauzalne vremenski diskretne LTI sustave je
IIR (Infinite Impulse Response) sustavi su oni za koje je impulsni
odziv definiran za beskonačni interval
FIR (Finate Impulse Response) sustavi su oni za koje je impulsni
odziv definiran za konačni interval
Originalan izraz pretvaramo u implicitan oblik preko lienarne diferencijske jednadžbe s konstantnim koeficijentima:
Zbog toga što je postupak izračuna odziva rekurzivan onda IIR nazivamo rekurzivnim sustavima.
za FIR vrijedi da su svi koef.
Sustavi su ujedno i IIR sustavi (beskonačno trajanje impulsnog odziva) koji se sastoje od:
Svi
Sustavi su ujedno i FIR sustavi (konačno trajanje impulsnog odziva) koji se sastoje od:
Svi
Utjecaj položaja polova i nula na frekvencijsku karakteristiku
fun fact: Prematanje faze za bilo koji višekratnik od
Fazno kašnjenje je negativni omjer fazne karakteristike i kružne frekvencije:
Kod filtriranja često želimo mijenjati samo amplitudu pojedinih
komponenti signala bez mijenjanja njihove faze. Taj zahtjev da filtar
utječe samo na amplitudu postavlja ogrnaičenja na željenu
fazno-frekvencijsku karakteristiku. Preciznije, željena fazna
karakteristika
Pošto je fazno kašnjenje konstata, to znači da sve komponente složenog signala kasne podjednako i zbog toga takva konstantna faza ne unosi fazna izobličenja.
Sustav koji ima različita kašnjenja za komponente složenog signala uzrokuje fazno izobličenje signala.
Razramtramo spektar
Fazno-frekvencijska karakteristika
Kod filtriranja često želimo mijenjati samo amplitudu pojedinih
komponenti signala bez mijenjanja njihove faze. Taj zahtjev da filtar
utječe samo na amplitudu postavlja ograničenje na željenu
fazno-frekvencijsku karakteristiku. Preciznije, željena fazna
karakteristika
gdje su
LTI sustave ili filtre čija fazna karakteristika je opisana gore navedenom jednadžbom nazivamo sustavima s generaliziranom linearnom fazom. Glavna prednost takvih filtara je da ne unose fazna izobličenja, no moguće je da uzrokuje grupno kašnjenje.
Možemo pokazati da se taj uvjet generalizirane linearne faze koji je iskazan u frekvencijskoj domeni transformira u uvjet simetričnost ili antisimetričnosti impulsnog odziva u vremenskoj domeni. Prema tome impulsni odziv svakog filtra generalizirane linearne faze mora imati točku oko koje je ili simetričan ili antisimetričan.
Promatramo signal:
Grupno vrijeme kašnjenja jest negativna derivacija fazno-frekvencijske karakteristike:
Idealno grupno vrijeme kašnjenje za filtre savršeno linearne faze je konstantno. Za dizajnirani IIR filtar stoga želimo da grupno vrijeme kašnjenje u području propuštanja bude što bliže konstanti. Grupno kašnjenje u području gušenja nas ne zanima jer se te frekvencije uklanjaju iz signala.
Primijetite sljedeće:
Onaj koji ne utječe na amplitudno frekvencijsku karakteristiku već samo na faznu. Zbog toga se nazivaju i korektori faze. Ovaj sustav ima recipročne polove (unutar jedinične kružnice) i nule (izvan jedinične kružnice).
Kod klasičnog dizajna filtara (i digitalnih i analognih) postupak dizajniranja se razdvaja u dva koraka:
Takvo razdvajanje dizajna u dva odvojena koraka se temelji na činjenici da svaku prijenosnu funkciju stabilnog i kauzalnog filtra možemo jednoznačno rastaviti na kaskadu (produkt prijenosnih funkcija) minimalno-faznog filtra (engl. minimum-phase filter) i svepropusnog filtra (engl. all-pass filter), odnosno
Minimalno-fazni filtar
Svepropusni filtar
Ova mogućnost rastava prijenosne funkcije znači da bilo koji IIR filtar možemo prikazati preko minimalno-faznog dijela koji u potpunosti određuje amplitudnu karakteristiku filtra, te preko svepropusnog dijela koji dodatno određuje faznu karakteristiku.
Četiri tipa amplitudno-selektivih filtara obzirom na oblik željene amplitudne karakteristike:
Odgovarajuće kratice u engleskom jeziku su LP za nisko-propusni filtar, od engl. low-pass, HP za visoko-propusni filtar, od engl. high-pass, BP za pojasno-propusni filtar, od engl. band-pass, i BS za pojasnu branu, od engl. band-stop.
Nisko-propusni filtar propušta frekvencije manje od
granične frekvencije
Visoko-propusni filtar guši frekvencije manje od
granične frekvencije
Pojasno-propusni filtar propušta samo frekvencije
unutar odabranog frekvencijskog pojasa
Pojasni brana gušiti sve frekvencije unutar
odabranog frekvencijskog pojasa
Ostali tipove filtara koji ne odgovaraju nekom od standardna četiri tipa amplitudno-selektivnih filtara možemo dizajnirati klasičnim postupkom samo ako se njihova amplitudna karakteristika može iskazati kao kombinacija standardnih NP, VP, PP i PB filtara, bilo kao njihov paralelni spoj (zbroj prijenosnih funkcija) bilo kao njihov kaskadni spoj (produkt prijenosnih funkcija).
FIR smo već definirali gore.
tldr; tako da impulsni odziv tih filtara zadovoljava svojstvo simetrije ili antisimetrije. To znači da su prvi i zadnji koeficijent jednaki, drugi i predzadnji jednaki, itd.
Svaki digitalni filtar u potpunosti je opisan položajem vlastitih polova i nula. FIR filtri nemaju ne-trivijalnih polova što znači da samo nule utječu na njihovu frekvencijsku karakteristiku. Prema tome kod FIR filtara nas zanima isključivo položaj vlastitih nula.
Znamo da ako je impulsni odziv realan onda su nule digitalnog filtra ili čisto realne ili su sparene u konjugirano-kompleksne parove.
Također, nula filtra postavljena na jediničnu kružnicu
Kod određivanja položaja nula od posebnog interesa su FIR filtri generalizirane linearne faze. Impulsni odziv tih filtara zadovoljava svojstvo simetrije ili antisimetrije.
To znači da su prvi i zadnji koeficijent jednaki, drugi i predzadnji jednaki, itd… Gore navedeno svojstvo impulsnog odziva omogućava nam da prijenosnu funkcija filtra izrazimo na dva različita načina,
Zamjena varijable
Također vrijedi
pa na kraju dobivamo
Polinome
Za takav polinom vrijedi da ako je
To svojstvo proizlazi odmah iz
Navedeno ograničenje položaja nula FIR filtra linearne faze je očekivano obzirom na to da simetrija ili antisimetrija pripadajućeg impulsnog odziva smanjuje broj stupnjeve slobode u optimizacijskom problemu na pola. Prema tome kada bi dizajnirali takve filtre korištenjem reprezentacije preko nula, polova i pojačanja onda je samo dio nula slobodno raspoređen i podložan optimizaciji, a preostale nule su unaprijed ograničene uvjetima koji ih sparuju u konjugirano-kompleksne i u recipročne parove.
FIR prednosti:
FIR Mane:
IIR prenosti:
IIR Mane:
Bazna funckije
Radi preglednosti uvodimo oznaku
Uoči da se “broj” i “nazivnik” kod
Sad spektar možemo zapisati kao:
a
Broj operacija za izrarvno računanje
dakle složenost je
Problem složenosti napadamo sa principom podjeli-pa-vladaj. Neka je
Sad svaki indeks rastavljamo na dva načina: 1. decimacija u vremenu
,
Razlažemo signal
Interpretacija člana
Složenost spajanja dvije
Dakle za rastavljanje imamo
Linearna konvolucija:
možemo ju lakše računati da izračunamo tablicu množenja pa radimo sumacije.
Cirkularna konvolucija (
Teorem o konvoluciji kaže da se konvolucija u vremenskoj domeni preslikava u običan umnožak u frekvencijskoj domeni.
Za velike N možemo koristit FFT da izračunamo cirkularnu konvoluciju
tako da prvo napravimo Fourierovu transformaciju nad oba signala. Sad
imamo 2 spektra. Pomnožimo ih, a zatim ih vratimo u signal sa inverzom
Foruierove transformacije. Složenost te operacije je
Nije komutativna! No vrijedi
Nije komutativna! No vrijedi
Teorem o korelaciji kaže da se korelacija u vremenskoj domeni preslikava u umnožak dva spektra od kojih je drugi konjugiran u frekvencijskoj domeni.
Vrijedi:
Ako djelujemo sa IDTFT operacijom dobivamo sljedeće:
Dakle korelaciju možemo izračunati pomoću DTFT-a na sljedeći način: - primjeni DTFT na oba signala - sad imamo dva spektra - primjeni konjugaciju za drugi spektar - pomnoži dva spektra - koristi IDTFT da se vratiš u vremensku domenu
Složenost za efikasno računanje korelacije uz FFT je
Realizacija filtra mora biti numerički stabila pa ćemo definirati česte numeričke probleme poput preljev, podljev i greška zaokruživanja
Svaki program koji izravno računa sljedeći izraz kažemo da odgovara direktnoj realizaciji
Kažemo da se na početku filtracije javlja prijelazna pojava jer
početne vrijednosti
Možemo smanjtii broj međurezutlata tako da spojimo ulaz
… TODO …
Konačna preciznost određena je broj značajnih znamenki koje pohranjujemo i položaj decimalne točke. Kod prelaska u konačnu preciznost uvijek se javljaju greške zaokruživanja. Te operacije onda označavamo sa zaokruživanjem simbola.
Rezultat približnih operacija u konačnoj preciznosti koje unose grešku mora biti jednak zaokruživanju rezultata kojeg bi dobili korištenjem artimetike beskonačne preciznosti na najbliži broj kojeg možemo zapisati.
Ako je najveći broj kojeg možemo zapisati 9.9 a najmanji -9.9, a najmanji pozitivan 0.1
onda 6.2 + 5.3 = 9.9 (nije 11.5) jer se dogodio preljev.
Također ovisi o vrsti artimetike: - aritmetika zasićenjem: 9.9 (stane na max broju) - modulna artimetika: -8.4 (dolazi do prematanja)
Pri množenju 6.2 * 0.2 = 1.24 ali zapravo dobivamo 1.2 i tu se dogodila greška zaokruživanja.
Također 0.2 * 0.2 = 0.04 ali zapravo dobivamo 0 i to je greška podljeva.
Prelijev se događa u čvoru grafa toka signala čija apsolutna
vrijednost preskoči dinamku
Ako postavimo ograničenje amplitude (ili dinamike) ulaza
onda za svaki čvor kroz analizu najgoreg slučaja možemo odrediti kolika može biti najveća vrijednost u tom čvoru.
Ako znamo sve dinamike čvorova
Dinamiku određujemo u tri koraka: 1. određujemo prijenosnu funkciju
Hv (z) između ulaza i čvora - označimo sve čvorove i iz njih postavljamo
jednadžbe. Na kraju svedemo sve na izraz: