LTI-järjestelmät ja signaalien suodatus.
LTI-järjestelmät ja signaalien suodatus.
LTI-järjestelmät ja signaalien suodatus.
Transform your PDFs into Flipbooks and boost your revenue!
Leverage SEO-optimized Flipbooks, powerful backlinks, and multimedia content to professionally showcase your products and significantly increase your reach.
Signaalien <strong>suodatus</strong><br />
Signaalinkäsittelyjärjestelmät muokkaavat lähtösignaalista tulosignaalin. Järjestelmät<br />
koostuvat (tai ne voidaan mallintaa) tyypillisesti suotimista. Näin <strong>signaalien</strong> <strong>suodatus</strong><br />
on keskeinen signaalinkäsittelyn operaatio.<br />
Voidaan helposti osoittaa, että tietyin edellytyksin tulo- <strong>ja</strong> lähtösignaalin välillä on<br />
riippuvuus. Kun järjestelmä (suodin) on lineaarinen <strong>ja</strong> aikainvariantti (<strong>LTI</strong>-järjestelmä)<br />
<strong>ja</strong> sen impulssivaste tunnetaan, voidaan lähtösignaali laskea aikatasossa tulon <strong>ja</strong><br />
impulssivasteen konvoluutiona. Fourier-muunnoksen konvoluutioteoreeman<br />
perusteella aikatason konvoluutiota vastaa taajuustason kertolasku, joten<br />
taajuustasossa lähtösignaali voidaan määrittää kertomalla tulo impulssivasteen<br />
Fourier-muunnoksella (eli järjestelmän taajuusvasteella).<br />
Aikatasossa impulssivaste määrittää suotimen ominaisuudet (esim. viiveen,<br />
nousua<strong>ja</strong>n <strong>ja</strong> <strong>suodatus</strong>ominaisuudet). Useimmiten suotimen ominaisuuksia<br />
tarkastellaan taajuustasossa muodostamalla amplitudi- <strong>ja</strong> vaihespektri.<br />
Amplitudispektri saadaan taajuusvasteen itseisarvona <strong>ja</strong> vaihespektri taajuusvasteen<br />
argumenttina. Vaihespektristä voidaan edelleen määrittää ryhmäviive, joka kertoo<br />
järjestelmässä syntyvät viiveet taajuuden funktiona.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 1
vaihespektri<br />
Fouriermuunnos<br />
Aikataso<br />
Taajuustaso<br />
tulo<br />
x(t)<br />
Signaalien <strong>suodatus</strong><br />
<strong>LTI</strong>-järjestelmä<br />
impulssivaste<br />
h(t)<br />
taajuusvaste<br />
(siirtofunktio)<br />
X ( f )<br />
H ( f )<br />
lähtö konvoluutio<br />
y( t)<br />
= h(<br />
t)<br />
∗ x(<br />
t)<br />
Y ( f ) = H ( f ) X ( f<br />
{ X }<br />
{ H } arg{ H ( f )}<br />
konvoluutioteoreema<br />
X ( f ) arg ( f ) H ( f ) arg ( f ) Y ( f )<br />
amplitudi vaihespektri<br />
amplitudi vaihe-<br />
amplitudi<br />
-spektri<br />
-spektri spektri<br />
-spektri<br />
)<br />
G( f )<br />
ryhmäviive<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 2
Signaalien <strong>suodatus</strong><br />
Lineaarinen aikainvariantti (<strong>LTI</strong>) järjestelmä<br />
Järjestelmällä tarkoitetaan tässä mitä tahansa ohjelmaa tai laitetta joka tuottaa<br />
vasteen (lähtösignaali) herätteestä (tulosignaali).<br />
x(t)<br />
Järjestelmä<br />
T<br />
y(t) = T{x(t)}<br />
Järjestelmä on lineaarinen, jos sille on voimassa<br />
1. Additiivisuus<br />
T{x 1<br />
(t) + x 2<br />
(t)} = T{x 1<br />
(t)} + T{x 2<br />
(t)} = y 1<br />
(t) + y 2<br />
(t) kaikille signaaleille x 1<br />
(t) <strong>ja</strong> x 2<br />
(t)<br />
2. Homogeenisuus<br />
T{ax(t)} = aT{x(t)} = ay(t) kaikille signaaleille x(t) <strong>ja</strong> vakiokertoimelle a.<br />
Järjestelmät, jotka eivät toteuta kumpaakin yllä olevista ehdoista ovat epälineaarisia.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 3
Signaalien <strong>suodatus</strong><br />
Lineaarisuuden ohella toinen tärkeä järjestelmien yleinen ominaisuus on<br />
aikainvarianttius (t. siirtoinvarianttius). Järjestelmä on aikainvariantti, jos sille pätee<br />
T{x(t-t 0<br />
)} = y(t-t 0<br />
). Herätteen viivästäminen aiheuttaa aikainvariantin järjestelmän<br />
vasteeseen siis samansuuruisen viiveen.<br />
x(t-t 0<br />
)<br />
Järjestelmä<br />
T<br />
y(t-t 0<br />
)= T{x(t-t 0<br />
)}<br />
Järjestelmiä, jotka ovat sekä lineaarisia että aikainvariantte<strong>ja</strong> kutsutaan <strong>LTI</strong>järjestelmiksi.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 4
Signaalien <strong>suodatus</strong><br />
Konvoluutio<br />
Tarkastellaan lineaarista aikainvarianttia järjestelmää T, jossa mielivaltainen heräte<br />
x(t) tuottaa vasteen y(t):<br />
x(t)<br />
<strong>LTI</strong>järjestelmä<br />
T<br />
y(t)= T{x(t)}<br />
Olkoon h(t) järjestelmän T impulssivaste (eli vaste, kun herätteenä on<br />
yksikköimpulssijono):<br />
{ δ ( t)<br />
} h(<br />
t)<br />
T =<br />
Koska järjestelmä T on aikainvariantti, pätee myös<br />
T<br />
{ δ ( t −τ<br />
)} = h(<br />
t −τ<br />
)<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 5
Signaalien <strong>suodatus</strong><br />
Herätesignaali x(t) voidaan esittää impulssifunktiota käyttäen muodossa<br />
∞<br />
∫<br />
−∞<br />
x ( t)<br />
= x(<br />
τ ) δ ( t −τ<br />
) dτ<br />
Tässä siis poimitaan impulssifunktiolla kaikki herätesignaalin arvot, koska muuttu<strong>ja</strong><br />
τ saa kaikki arvot väliltä -∞…+∞ <strong>ja</strong> impulssifunktion määritelmän mukaisesti<br />
esimerkiksi x(1000)δ(t-1000)=x(1000), koska δ(t-1000) on arvoltaan 1, kun t = 1000<br />
<strong>ja</strong> muulloin nolla. Vastaavasti x(2000)δ(t-2000) tuottaa herätteen arvon a<strong>ja</strong>nhetkellä t<br />
= 2000, jne. Kun kaikki mahdolliset herätesignaalin arvot summataan yhteen<br />
(integrointi) saadaan tuloksena herätesignaali x(t).<br />
Laitetaan järjestelmään T herätteeksi x(t) <strong>ja</strong> hyödynnetään järjestelmän T<br />
lineaarisuus (homogeenisuus <strong>ja</strong> additiivisuus) ominaisuutta:<br />
y(<br />
t)<br />
∞<br />
∞<br />
∞<br />
⎧<br />
⎫<br />
= T ⎨ ∫<br />
⎬ ∫<br />
∫ )<br />
⎩−∞<br />
⎭ −∞<br />
−∞<br />
{ x(<br />
t)<br />
} = T x(<br />
τ ) δ ( t −τ<br />
) dτ<br />
= x(<br />
τ ) T{ δ ( t −τ<br />
)} dτ<br />
= x(<br />
τ ) h(<br />
t −τ<br />
dτ<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 6
Signaalien <strong>suodatus</strong><br />
Tulos kertoo siis, että minkä tahansa lineaarisen aikainvariantin järjestelmän vaste<br />
mielivaltaiseen herätesignaaliin voidaan määrittää edellä olevan kaavan mukaisesti<br />
herätesignaalin <strong>ja</strong> järjestelmän impulssivasteen integraalina. Integraalia kutsutaan<br />
konvoluutiointegraaliksi tai konvoluutioksi, <strong>ja</strong> yleisesti käytetään merkintää<br />
∞<br />
y ( t)<br />
x(<br />
t)<br />
∗ h(<br />
t)<br />
= x(<br />
τ ) h(<br />
t −τ<br />
) dτ<br />
= h(<br />
τ ) x(<br />
t −τ<br />
) dτ<br />
= h(<br />
t)*<br />
x(<br />
t)<br />
= ∫ ∫<br />
−∞<br />
Edellä järjestelmä voi olla esimerkiksi suodin (t. suodatin), joka rajoittaa vasteen<br />
taajuudet jollekin tietylle välille, tai (tiedonsiirto)kanava, joka liittää lähettimen <strong>ja</strong><br />
vastaanottimen toisiinsa. Yleisemmin konvoluutio voidaan laskea minkä tahansa<br />
kahden signaalin välillä.<br />
Esimerkki. Signaalien x 1<br />
(t)=e -αt u(t) <strong>ja</strong> x 2<br />
(t)=e -βt u(t) konvoluutio, kun α>0 <strong>ja</strong> β>0.<br />
∞<br />
−∞<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 7
©Jyrki Laitinen TL5231 Signaaliteoria S2004 8<br />
Signaalien <strong>suodatus</strong><br />
[ ] ( )<br />
t<br />
t<br />
t<br />
t<br />
t<br />
t<br />
t<br />
t<br />
t<br />
t<br />
t<br />
e<br />
e<br />
e<br />
e<br />
e<br />
e<br />
d<br />
e<br />
e<br />
d<br />
e<br />
e<br />
t<br />
y<br />
t<br />
t<br />
u<br />
u<br />
t<br />
t<br />
u<br />
u<br />
d<br />
t<br />
u<br />
e<br />
u<br />
e<br />
d<br />
t<br />
x<br />
x<br />
t<br />
x<br />
t<br />
x<br />
t<br />
y<br />
α<br />
β<br />
β<br />
α<br />
β<br />
τ<br />
β<br />
α<br />
β<br />
τ<br />
β<br />
α<br />
β<br />
τ<br />
β<br />
ατ<br />
τ<br />
β<br />
ατ<br />
β<br />
α<br />
β<br />
α<br />
β<br />
α<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
τ<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
−<br />
∞<br />
−∞<br />
−<br />
−<br />
−<br />
∞<br />
−∞<br />
−<br />
−<br />
=<br />
−<br />
−<br />
= −<br />
−<br />
−<br />
=<br />
=<br />
=<br />
⇒<br />
⎩<br />
⎨<br />
⎧ <<br />
<<br />
=<br />
−<br />
⇒<br />
⎩<br />
⎨<br />
⎧ <<br />
=<br />
−<br />
⎩<br />
⎨<br />
⎧ ><br />
=<br />
−<br />
=<br />
−<br />
=<br />
∗<br />
=<br />
∫<br />
∫<br />
∫<br />
∫<br />
1<br />
1<br />
1<br />
/<br />
)<br />
(<br />
muulloin<br />
0,<br />
0<br />
1,<br />
)<br />
(<br />
)<br />
(<br />
muulloin<br />
0,<br />
1,<br />
)<br />
(<br />
muulloin<br />
0,<br />
0<br />
1,<br />
)<br />
(<br />
)<br />
(<br />
)<br />
(<br />
)<br />
(<br />
)<br />
(<br />
)<br />
(<br />
)<br />
(<br />
)<br />
(<br />
)<br />
(<br />
)<br />
(<br />
0<br />
0<br />
)<br />
(<br />
0<br />
)<br />
(<br />
)<br />
(<br />
2<br />
1<br />
2<br />
1
Signaalien <strong>suodatus</strong><br />
Alla on piirretty signaalit x 1<br />
(t), x 2<br />
(t) <strong>ja</strong> x 1<br />
(t)*x 2<br />
(t), kun α = 0.1 <strong>ja</strong> β=1.<br />
1<br />
0.9<br />
x1(t)<br />
0.8<br />
0.7<br />
x1(t)*x2(t)<br />
0.6<br />
0.5<br />
0.4<br />
0.3<br />
0.2<br />
x2(t)<br />
0.1<br />
0<br />
0 1 2 3 4 5 6 7 8 9 10<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 9
Signaalien <strong>suodatus</strong><br />
Stabiilius<br />
Järjestelmä on stabiili, jos vaste on äärellinen kaikilla herätteen äärellisillä arvoilla.<br />
Oletetaan, että herätesignaali x(t) on rajoitettu siten, että sen itseisarvo on pienempi<br />
tai yhtäsuuri kuin äärellinen reaaliluku M, eli<br />
x ( t)<br />
≤ M ,<br />
− ∞ < t<br />
< ∞<br />
Vaste y(t) voidaan määrittää konvoluutiointegraalina<br />
y(<br />
t)<br />
=<br />
∞<br />
∫<br />
−∞<br />
h(<br />
τ ) x(<br />
t<br />
−τ<br />
) dτ<br />
⇒<br />
y(<br />
t)<br />
≤<br />
∞<br />
∫<br />
−∞<br />
h(<br />
τ ) x(<br />
t<br />
−τ<br />
) dτ<br />
=<br />
∞<br />
∫<br />
−∞<br />
h(<br />
τ )<br />
x(<br />
t<br />
−τ<br />
) dτ<br />
≤<br />
M<br />
∞<br />
∫<br />
−∞<br />
h(<br />
τ ) dτ<br />
Vaste y(t) pysyy siis äärellisenä <strong>ja</strong> järjestelmä stabiilina, jos on voimassa ehto<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 10
Signaalien <strong>suodatus</strong><br />
∞<br />
∫<br />
−∞<br />
h ( t)<br />
dt<br />
< ∞<br />
Kausaalius<br />
Järjestelmä on kausaali, jos se tuottaa vastetta vasta kun heräte on annettu.<br />
Kausaalille <strong>LTI</strong>-järjestelmälle pätee<br />
h( t)<br />
= 0, t <<br />
0<br />
Reaaliaikajärjestelmät ovat yleensä kausaale<strong>ja</strong>. Joissakin sovelluksissa signaale<strong>ja</strong><br />
voidaan tallentaa ennen käsittelyä muistiin, jolloin signaalia muokkaava järjestelmä<br />
voi olla myös ei-kausaali. Esimerkiksi kuvasignaalit käsitellään myös monissa<br />
reaaliaikajärjestelmissä siten, että yksittäiset kuvat talletetaan kokonaisuudessaan<br />
muistiin ennen käsittelyoperaatioita. Tällöin yksittäiset suodatukset voivat olla eikausaale<strong>ja</strong>.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 11
Signaalien <strong>suodatus</strong><br />
Taajuusvaste<br />
Konvoluutioteoreeman mukaisesti aikatasossa tapahtuvaa konvoluutiota vastaa<br />
taajuustasossa kertolasku:<br />
y ( t)<br />
= x(<br />
t)<br />
∗h(<br />
t)<br />
⇔ Y ( ω)<br />
= X ( ω)<br />
H ( ω)<br />
Tässä X(ω) on herätteen x(t), H(ω) impulssivasteen h(t) <strong>ja</strong> Y(ω) vasteen y(t) Fouriermuunnos.<br />
Impulssivasteen Fourier-muunnosta H(ω) kutsutaan järjestelmän<br />
taajuusvasteeksi (t. siirtofunktioksi).<br />
H ( ω)<br />
=<br />
F<br />
{ h(<br />
t)<br />
}<br />
1<br />
H(ω)<br />
=<br />
=<br />
∞<br />
∫<br />
−∞<br />
h(<br />
t)<br />
e<br />
Y ( ω)<br />
X ( ω)<br />
− jωt<br />
dt<br />
δ(t)<br />
x(t)<br />
X(ω)<br />
<strong>LTI</strong>järjestelmä<br />
h(t)<br />
y(t)= x(t)*h(t)<br />
Y(ω)=X(ω)H(ω)<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 12
Signaalien <strong>suodatus</strong><br />
Vasteen Fourier-muunnos saadaan siis kertomalla taajuusvaste H(ω) <strong>ja</strong> herätteen<br />
Fourier-muunnos X(ω) keskenään. Usein on käytännössä konvoluutiota<br />
tehokkaampaa laskea nopealla Fourier-muunnoksella taajuusvaste <strong>ja</strong> herätteen<br />
Fourier-muunnos, jotka kertomalla saadaan vasteen Fourier-muunnos. Vastesignaali<br />
saadaan tällöin käänteisellä Fourier-muunnoksella.<br />
y(<br />
t)<br />
∞<br />
∫<br />
−∞<br />
1<br />
j<br />
= ω ω<br />
ω t<br />
X ( ) H ( ) e dω<br />
2π<br />
Esimerkiksi tietoliikennesovelluksissa tiedonsiirtokanava usein vääristää siirrettävää<br />
signaalia. Jos kanavan impulssi- tai taajuusvaste tunnetaan, voidaan vääristymä<br />
kor<strong>ja</strong>ta <strong>ja</strong>kamalla vastaanotetun signaalin Fourier-muunnos kanavan<br />
taajuusvasteella:<br />
X ( ω)<br />
=<br />
Y ( ω)<br />
H ( ω)<br />
Vääristymätön signaali x(t) saadaan nyt käänteisellä Fourier-muunnoksella<br />
signaalista X(ω). Tässä suoritetaan ns. dekonvoluutio taajuustasossa.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 13
Amplitudi- <strong>ja</strong> vaihevaste<br />
Signaalien <strong>suodatus</strong><br />
Taajuusvaste H(ω) on yleisessä tapauksessa kompleksiarvoinen funktio, joka<br />
voidaan esittää itseisarvonsa <strong>ja</strong> vaiheensa avulla:<br />
H ( ω)<br />
=<br />
H ( ω)<br />
e<br />
j arg<br />
{ H ( ω )}<br />
Itseisarvoa |H(ω)| sanotaan järjestelmän amplitudivasteeksi (t.<br />
amplitudispektriksi) <strong>ja</strong> vaihefunktiota arg{H(ω)} järjestelmän vaihevasteeksi (t.<br />
vaihespektriksi).<br />
Jos impulssivaste h(t) on reaaliarvoinen, on taajuusvaste H(ω) ns.<br />
konjugaattisymmetrinen, jolloin<br />
H ( ω)<br />
arg<br />
= H ( −ω<br />
)<br />
|H(ω)|<br />
{ H ( −ω<br />
)} = −arg{ H ( ω)<br />
}<br />
θ h<br />
(ω)<br />
H ( ω)<br />
= H ( −ω<br />
)<br />
θ ( −ω<br />
) = −θ<br />
( ω)<br />
h<br />
h<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 14
Merkitään<br />
H ( ω)<br />
=<br />
X ( ω)<br />
=<br />
H ( ω)<br />
e<br />
X ( ω)<br />
e<br />
j arg<br />
arg<br />
{ H ( ω )}<br />
{ X ( ω )}<br />
Signaalien <strong>suodatus</strong><br />
Tällöin vastesignaalin taajuustason esitys voidaan kirjoittaa muotoon<br />
Y ( ω)<br />
= Y ( ω)<br />
e<br />
j arg<br />
{ Y ( ω )} j arg{ X ( ω )} j arg{ H ( ω )}<br />
=<br />
=<br />
X ( ω)<br />
e<br />
X ( ω)<br />
H ( ω)<br />
e<br />
j<br />
H ( ω)<br />
e<br />
[ arg{ X ( ω )} + arg{ H ( ω )})<br />
]<br />
Vastesignaalin amplitudispektri saadaan siis kertomalla herätteen <strong>ja</strong> suotimen<br />
amplitudispektrit keskenään <strong>ja</strong> vastesignaalin vaihespektri puolestaan laskemalla<br />
yhteen herätteen <strong>ja</strong> suotimen vaihespektrit.<br />
Y ( ω)<br />
=<br />
arg<br />
X ( ω)<br />
H ( ω)<br />
{ Y ( ω)<br />
} = arg{ X ( ω)<br />
} + arg{ H ( ω)<br />
}<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 15
Signaalien <strong>suodatus</strong><br />
Ryhmäviive<br />
Ryhmäviive määrittää suotimessa herätteeseen syntyvän viiveen taajuuden<br />
funktiona. Ryhmäviive G(ω) määritellään järjestelmän vaihevasteen derivaattaa<br />
käyttäen muodossa<br />
G(<br />
ω)<br />
= −<br />
dθ<br />
( ω)<br />
dω<br />
⇒ G(<br />
f<br />
1 dθ<br />
h<br />
( f<br />
) = − ⋅<br />
2π<br />
df<br />
h<br />
)<br />
Jos suotimen vaihevaste on lineaarinen, saa ryhmäviive vakioarvon <strong>ja</strong> suodatettavat<br />
signaalit viivästyvät taajuudesta riippumatta vakioa<strong>ja</strong>n.<br />
arg{X(ω)}=-kω<br />
-d[arg{X(ω)]/dω = k<br />
ω<br />
ω<br />
Kulmakerroin -k<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 16
Signaalien <strong>suodatus</strong><br />
Esimerkki. Tiedonsiirto epälineaarisen vaihevasteen omaavassa kanavassa.<br />
Siirretään 0.5 Hz:n taajuisella sinisignaalilla moduloitu 5 Hz:n taajuinen kantoaalto<br />
tiedonsiirtokanavassa, jossa ryhmäviive ei ole siirrettävän signaalin taajuusalueella<br />
vakio. Siirrettyyn signaaliin syntyy vääristymiä, koska moduloidun signaalin taajuudet<br />
viivästyvät kanavassa eripituisen a<strong>ja</strong>n.<br />
Kanavan ominaisuudet<br />
1<br />
0.5<br />
0<br />
-0.5<br />
-1<br />
0 2 4 6 8 10<br />
Amplitudi<br />
Vaihe [rad]<br />
Ryhmäviive [s ]<br />
1<br />
0.8<br />
0.6<br />
0.4<br />
0.2<br />
0<br />
0 2.5 5<br />
0<br />
-5<br />
-10<br />
0 2.5 5<br />
1.5<br />
0.5<br />
1<br />
0<br />
0 2.5 5<br />
Taajuus [Hz]<br />
1<br />
0.5<br />
0<br />
-0.5<br />
-1<br />
0 2 4 6 8 10<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 17
Peräkkäiset järjestelmät<br />
Signaalien <strong>suodatus</strong><br />
Tarkastellaan yleisemmin tilannetta, jossa kaksi <strong>LTI</strong>-suodinta on asetettu peräkkäin<br />
eli kaskadiin. Suotimien impulssivasteet ovat h 1<br />
(t) <strong>ja</strong> h 2<br />
(t) <strong>ja</strong> taajuusvasteet<br />
vastaavasti H 1<br />
(ω) <strong>ja</strong> H 2<br />
(ω).<br />
x(t)<br />
h 1<br />
(t)<br />
h 2<br />
(t)<br />
y(t)<br />
Voidaan osoittaa, että nämä on mahdollista yhdistää konvoloimalla impulssivasteet<br />
keskenään<br />
x(t)<br />
h 1<br />
(t) * h 2<br />
(t)<br />
y(t) = [h 1<br />
(t) * h 2<br />
(t)] * x(t)<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 18
Signaalien <strong>suodatus</strong><br />
Taajuustasossa konvoluutiota vastaa kertolasku, joten taajuusvasteet voidaan<br />
yhdistää kertomalla ne keskenään.<br />
X(ω)<br />
H 1<br />
(ω)<br />
H 2<br />
(ω)<br />
y(ω)<br />
X(ω)<br />
H 1<br />
(ω) · H 2<br />
(ω)<br />
Y(ω) = H 1<br />
(ω) · H 2<br />
(ω)· X(ω)<br />
Yhdistetyn järjestelmän amplitudispektri |H(ω)| saadaan nyt kertomalla<br />
yksittäisten järjestelmien amplitudispektrit keskenään <strong>ja</strong> vaihespektri arg{H(ω)}<br />
puolestaan laskemalla yksittäisten järjestelmien vaihespektrit yhteen.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 19
Signaalien <strong>suodatus</strong><br />
H ( ω)<br />
= H<br />
= H ( ω)<br />
e<br />
⇒<br />
H ( ω)<br />
arg<br />
=<br />
1<br />
H<br />
( ω)<br />
⋅ H<br />
1<br />
( ω)<br />
⋅<br />
2<br />
( ω)<br />
=<br />
H<br />
( ω)<br />
( ω)<br />
{ H ( ω )} j arg{ H ( ω )}<br />
{ H ( ω )}<br />
{ H ( ω)<br />
} = arg{ H ( ω)<br />
} + arg{ H ( ω)<br />
}<br />
1<br />
2<br />
=<br />
H<br />
H<br />
1<br />
1<br />
( ω)<br />
⋅<br />
⋅e<br />
j arg<br />
2<br />
j arg<br />
H<br />
2<br />
1<br />
( ω)<br />
e<br />
j<br />
⋅<br />
H<br />
2<br />
( ω)<br />
[ arg{ H ( ω )} + arg{ H ( ω )}]<br />
1<br />
⋅e<br />
2<br />
2<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 20
Rinnakkaiset järjestelmät<br />
Signaalien <strong>suodatus</strong><br />
Tarkastellaan tilannetta, jossa kaksi <strong>LTI</strong>-suodinta on asetettu rinnakkain. Suotimien<br />
impulssivasteet ovat h 1<br />
(t) <strong>ja</strong> h 2<br />
(t) <strong>ja</strong> taajuusvasteet vastaavasti H 1<br />
(ω) <strong>ja</strong> H 2<br />
(ω).<br />
Rinnakkaiset järjestelmät voidaan tarvittaessa yhdistää laskemalla impulssivasteet<br />
tai taajuusvasteet yhteen.<br />
h 1<br />
(t)<br />
x(t)<br />
y(t) = h 1<br />
(t)*x(t) + h 2<br />
(t) *x(t)<br />
h 2<br />
(t)<br />
x(t)<br />
h 1<br />
(t) + h 2<br />
(t)<br />
y(t) = [h 1<br />
(t) + h 2<br />
(t)] * x(t)<br />
X(ω)<br />
H 1<br />
(ω) + H 2<br />
(ω)<br />
Y(ω) = [H 1<br />
(ω) + H 2<br />
(ω)] · X(ω)<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 21
Signaalien <strong>suodatus</strong><br />
Alla on määritetty yhdistetyn suotimen amplitudi <strong>ja</strong> vaihespektri. Huomaa merkinnät<br />
A 1<br />
= |H 1<br />
(ω)|, A 2<br />
= |H 2<br />
(ω)|, θ 1<br />
= arg{H 1<br />
(ω)} <strong>ja</strong> θ 2<br />
= arg{H 2<br />
(ω)}.<br />
H ( ω)<br />
= H<br />
1<br />
( ω)<br />
+ H<br />
= A1<br />
⋅e<br />
= A1<br />
⋅<br />
= A1<br />
⋅<br />
⇒<br />
H ( ω)<br />
=<br />
arg<br />
{ H ( ω )} j arg{ H ( ω )}<br />
{ cosθ1<br />
+ j sinθ1} + A2<br />
⋅{ cosθ<br />
2<br />
+ j sinθ<br />
2}<br />
cosθ<br />
+ A ⋅cosθ<br />
+ j ⋅{ A ⋅sinθ<br />
+ A ⋅sinθ<br />
}<br />
2<br />
( A ⋅cosθ<br />
+ A ⋅cosθ<br />
) + ( A ⋅sinθ<br />
+ A ⋅sinθ<br />
)<br />
1<br />
2<br />
( ω)<br />
1<br />
⎡<br />
⎢<br />
⎣ A<br />
⋅e<br />
A<br />
⎤<br />
1<br />
⋅sinθ1<br />
+ A2<br />
⋅sinθ<br />
2<br />
{ H ( ω)<br />
} = arctan<br />
⎥ ⎦<br />
1<br />
=<br />
2<br />
H<br />
1<br />
( ω)<br />
⋅e<br />
1<br />
jθ<br />
1<br />
2<br />
+<br />
1<br />
⋅cosθ<br />
+ A<br />
2<br />
A<br />
j arg<br />
2<br />
2<br />
1<br />
1<br />
jθ<br />
⋅cosθ<br />
2<br />
2<br />
+<br />
2<br />
1<br />
H<br />
2<br />
( ω)<br />
⋅e<br />
2<br />
1<br />
2<br />
2<br />
1<br />
2<br />
2<br />
2<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 22
Vääristymätön tiedonsiirto<br />
Signaalien <strong>suodatus</strong><br />
Monissa sovelluksissa signaali halutaan siirtää tiedonsiirtokanavaa pitkin<br />
mahdollisimman muuttumattomana. Tällöin puhutaan vääristymättömästä<br />
tiedonsiirrosta, jossa signaali toteuttaa ehdon<br />
y( t)<br />
= kx(<br />
t − td<br />
)<br />
Eli vääristämättömässä tiedonsiirtokanavassa signaalin vahvistus on k <strong>ja</strong> signaali voi<br />
viivästyä a<strong>ja</strong>n t d<br />
verran. Yhtälön Fourier-muunnos on<br />
Y ( ω)<br />
=<br />
⇒ H ( ω)<br />
=<br />
⇒<br />
H ( ω)<br />
θ ( ω)<br />
h<br />
X ( ω)<br />
H ( ω)<br />
=<br />
ke<br />
k<br />
− jωt<br />
= −ωt<br />
d<br />
d<br />
.<br />
= kX ( ω)<br />
e<br />
− jωt<br />
Signaalin x(t) jokainen taajuuskomponentti siis vahvistuu tekijällä k <strong>ja</strong> viivästyy<br />
tekijällä t d<br />
.<br />
d<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 23
Vääristymät tiedonsiirtokanavassa<br />
Signaalien <strong>suodatus</strong><br />
Vääristämättömän <strong>LTI</strong>-tiedonsiirtokanavan amplitudispektri on vakio <strong>ja</strong> vaihespektri<br />
lineaarinen. Jos amplitudi- tai vaihespektri ei täytä näitä ehto<strong>ja</strong>, signaali levenee<br />
aikatasossa kulkiessaan kanavan läpi. Tällaista vääristymistä kutsutaan dispersioksi<br />
<strong>ja</strong> se on erityisen ongelmallista aika<strong>ja</strong>koisessa (TDM) tiedonsiirrossa, jossa<br />
peräkkäisten pulssien <strong>ja</strong> siten myös naapurikanavien välille voi syntyä interferenssiä.<br />
Taajuus<strong>ja</strong>koisessa (FDM) tiedonsiirrossa virheet syntyvät kunkin kanavan sisällä<br />
eivätkä naapurikanavat häiritse toisiaan tällaisten vääristymien seurauksena.<br />
<strong>LTI</strong>-kanavassa syntyvää vääristymää voidaan vähentää liittämällä<br />
tiedonsiirtokanavaan ylimääräinen järjestelmä, joka kompensoi epäideaalisen<br />
kanavan vaikutusta. Tällaista tekniikkaa kutsutaan kanavan ekvalisoinniksi <strong>ja</strong><br />
vastaavaa järjestelmää ekvalisaattoriksi.<br />
tulo<br />
H c (f)<br />
H eq (f)<br />
lähtö<br />
<strong>LTI</strong>-kanava<br />
ekvalisaattori<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 24
Signaalien <strong>suodatus</strong><br />
Ekvalisaattori suunnitellaan siten, että siirtotien kokonaissiirtofunktio vastaa<br />
vääristämättömän kanavan siirtofunktiota, eli<br />
H<br />
c<br />
(<br />
f<br />
) H<br />
eq<br />
(<br />
f<br />
)<br />
=<br />
ke<br />
− j 2πft d<br />
Ekvalisaattorin siirtofunktioksi saadaan tällöin<br />
H<br />
eq<br />
(<br />
f<br />
)<br />
=<br />
ke<br />
H<br />
− j 2πft<br />
c<br />
d<br />
( f )<br />
Käytännössä ekvalisointi toteutetaan siten, että sen siirtofunktio approksimoi<br />
mahdollisimman hyvin ideaalista siirtofunktiota. Tällainen siirtofunktio toteutetaan<br />
yleensä FIR-suodinrakenteella, jossa lähtöarvo muodostetaan summaamalla<br />
suodinkertoimilla painotettu<strong>ja</strong> tulosignaalin arvo<strong>ja</strong> yhteen eli tulosignaalin arvojen<br />
lineaarikombinaationa.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 25
Signaalien <strong>suodatus</strong><br />
Jos tiedonsiirtokanava on epälineaarinen, syntyy lähtösignaaliin kanavassa sellaisia<br />
taajuuksia, joita alkuperäisessä signaalissa ei ole, mikä on ongelma erityisesti<br />
taajuus<strong>ja</strong>koisessa tiedonsiirrossa. Epälineaarista kanavaa voidaan mallintaa<br />
sar<strong>ja</strong>kehitelmänä<br />
2<br />
3<br />
k<br />
y(<br />
t)<br />
= a g(<br />
t)<br />
+ a g ( t)<br />
+ a g ( t)<br />
+ L + a g ( ) +L<br />
1 2<br />
3<br />
k<br />
t<br />
Esimerkiksi termin a 2<br />
g 2 (t) spektri saadaan Fourier-muunnoksen<br />
kertolaskuominaisuuden perusteella spektrin G(f) konvoluutiona itsensä kanssa. Eli,<br />
jos g(t):n kaistanleveys on B, niin g 2 (t):n kaistanleveys on 2B. Tällöin g k (t):n<br />
kaistanleveys on kB. Epälineaarisuudesta syntyvien uusien taajuuksien määrä<br />
riippuu siis epälineaarisuuden asteesta.<br />
Tiedonsiirtoon aiheutuu vääristymiä myös silloin, kun signaali saapuu<br />
vastaanottimeen kahta tai useampaa eri viiveen omaavaa siirtotietä pitkin.<br />
Esimerkiksi kaapelissa, jossa impedanssisovitus on virheellinen, vastaanottopäähän<br />
saapuu alkuperäinen signaali sekä sen hei<strong>ja</strong>stuksia erilaisten viiveiden jälkeen.<br />
Radiosiirrossa eri kohteista syntyy hei<strong>ja</strong>stuksia, jotka saapuvat vastaanottimeen eri<br />
aikana.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 26
Signaalien <strong>suodatus</strong><br />
Edellä on oletettu, että tiedonsiirtokanavan ominaisuudet säilyvät muuttumattomina<br />
a<strong>ja</strong>n hetkestä toiseen. Esimerkiksi radiosiirrossa ilmankehän ominaisuuksien<br />
vaihtelun vuoksi signaalin vaimennus muuttuu a<strong>ja</strong>n funktiona. Tätä voidaan<br />
kompensoida mm. automaattisella tason säädöllä (AGC).<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 27
Signaalien <strong>suodatus</strong><br />
Suodattimet<br />
Suodin on järjestelmä, joka muokkaa herätteen amplitudia- <strong>ja</strong> vaihetta halutulla<br />
tavalla. Suodatuksessa herätteen taajuussisältö yleensä muuttuu. Perussuotimet<br />
ovat alipäästö-,ylipäästö-, kaistanpäästö- <strong>ja</strong> kaistanestosuodin, jotka muokkaavat<br />
amplitude<strong>ja</strong>. Tämän lisäksi tärkeitä suotimia ovat vaiheen muokkaukseen käytettävät<br />
kokopäästösuodin (kulkuaikakor<strong>ja</strong>in) <strong>ja</strong> Hilbert-muunnin.<br />
Suotimet määritellään yleensä antamalla päästö-, esto- <strong>ja</strong> siirtymäkaistan<br />
ra<strong>ja</strong>taajuudet. Suotimelle voidaan taajuusvasteen perusteella määritellä myös<br />
kaistanleveys esimerkiksi taajuusvasteen nollakohtien tai –3 dB:n pisteen<br />
perusteella.<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 28
Signaalien <strong>suodatus</strong><br />
|H(ω)|<br />
1.2<br />
1<br />
Siirtymäkaista<br />
Päästökaista<br />
Estokaista<br />
0.8<br />
0.6<br />
0.707|H(ω)|<br />
0.4<br />
Kaistanleveys ω Β<br />
0.2<br />
0<br />
-0.2<br />
0 ω 1 ω B ω 2<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 29
Signaalien <strong>suodatus</strong><br />
1.2<br />
|H(w)|<br />
1<br />
0.8<br />
0.6<br />
Siirtymäkaistat<br />
Päästökaista<br />
Estokaista<br />
Estokaista<br />
0.707|H(ω)|<br />
0.4<br />
0.2<br />
0<br />
-0.2<br />
Kaistanleveys<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 30
Ideaalinen alipäästösuodin<br />
Signaalien <strong>suodatus</strong><br />
Ideaalisen alipäästösuotimen taajuusvaste saa päästökaistalla arvon 1 <strong>ja</strong><br />
estokaistalla arvon 0. Oletetaan, että suotimen ra<strong>ja</strong>taajuus on B. Taajuusvaste<br />
voidaan tällöin esittää muodossa<br />
H<br />
LPF<br />
− j<br />
⎧e<br />
( ω)<br />
= ⎨<br />
⎩ 0,<br />
ω<br />
t d<br />
,<br />
ω ≤ B<br />
muulloin<br />
Impulssivaste voidaan määrittää käänteisellä Fourier-muunnoksella<br />
h<br />
LPF<br />
∞<br />
1<br />
j t B<br />
( t)<br />
= ∫<br />
ω<br />
H<br />
LPF<br />
( ω)<br />
e dω<br />
= sinc<br />
d<br />
2π<br />
π<br />
−∞<br />
[ B(<br />
t − t )]<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 31
Signaalien <strong>suodatus</strong><br />
h LPF<br />
(t)<br />
B/π<br />
Ideaalisen<br />
alipäästösuotimen<br />
impulssivaste<br />
Ideaalisen<br />
alipäästösuotimen<br />
amplitudi- <strong>ja</strong><br />
vaihevaste<br />
|H LPF<br />
(ω)|<br />
1<br />
t<br />
-B<br />
0 B<br />
θ LPF<br />
(ω)<br />
ω<br />
t d<br />
2π/Β<br />
B<br />
-B 0<br />
ω<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 32
Ideaalinen ylipäästösuodin<br />
Signaalien <strong>suodatus</strong><br />
Ideaalisen ylipäästösuotimen taajuusvaste saa päästökaistalla arvon 1 <strong>ja</strong><br />
estokaistalla arvon 0. Oletetaan, että suotimen ra<strong>ja</strong>taajuus on B. Taajuusvaste<br />
voidaan tällöin esittää muodossa<br />
H<br />
HPF<br />
− jωtd<br />
⎧e<br />
, ω ≥ B − jωtd<br />
( ω)<br />
= ⎨<br />
= e − H<br />
LPF<br />
( ω)<br />
⎩0,<br />
muulloin<br />
|H HPF<br />
(ω)|<br />
θ HPF<br />
(ω)<br />
1<br />
-B<br />
0 B<br />
ω<br />
B<br />
-B 0<br />
ω<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 33
Ideaalinen kaistanpäästösuodin<br />
Signaalien <strong>suodatus</strong><br />
Ideaalisen kaistanpäästösuotimen taajuusvaste voidaan esittää muodossa<br />
H<br />
BPF<br />
( ω)<br />
=<br />
−<br />
⎧e<br />
⎨<br />
⎩0,<br />
jω<br />
t d<br />
,<br />
B 1<br />
≤ ω ≤ B 2<br />
muulloin<br />
|H BPF<br />
(ω)| θ BPF<br />
(ω)<br />
1<br />
ω<br />
-B 2<br />
-B 1<br />
0 B 1<br />
B 2<br />
-B 2<br />
-B 1<br />
0<br />
B 1<br />
B 2<br />
ω<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 34
Signaalien <strong>suodatus</strong><br />
Ideaalinen kaistanestosuodin<br />
Ideaalisen kaistanestosuotimen taajuusvaste voidaan esittää muodossa<br />
jωt<br />
H ( ω)<br />
= e<br />
− d − H ( ω)<br />
BSF<br />
BPF<br />
|H BSF<br />
(ω)|<br />
1<br />
θ BSF<br />
(ω)<br />
ω<br />
-B 2<br />
-B 1<br />
0 B 1<br />
B 2<br />
-B 2<br />
B 1<br />
B 2<br />
ω<br />
-B 1<br />
0<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 35
Signaalien <strong>suodatus</strong><br />
Hilbert-muunnin<br />
Tarkastellaan suodinta, jonka taajuusvaste on muotoa<br />
⎧<br />
⎪e<br />
H ( ω)<br />
= ⎨<br />
⎪<br />
⎩e<br />
π<br />
− j<br />
2<br />
π<br />
+ j<br />
2<br />
= −<br />
= +<br />
j,<br />
j,<br />
ω > 0<br />
ω < 0<br />
= −<br />
j sgn( ω)<br />
Suodin toteuttaa positiivisilla taajuuksilla -π/2:n suuruisen vaihesiirron <strong>ja</strong><br />
vastaavasti negatiivisilla taajuuksilla +π/2:n suuruisen vaihesiirron.<br />
Amplitudivahvistus on sen si<strong>ja</strong>an kaikilla taajuuksilla 1, joten amplitudi ei<br />
suodatuksessa muutu. Tällaista suodinta kutsutaan Hilbert-muuntimeksi <strong>ja</strong> sillä on<br />
keskeinen merkitys tiedonsiirtotekniikassa kapeakaistaisten <strong>signaalien</strong> käsittelyssä.<br />
|H Hilbert<br />
(ω)| θ Hilbert<br />
(ω)<br />
1<br />
+π/2<br />
ω<br />
0 0<br />
-π/2<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 36<br />
ω
Signaalien <strong>suodatus</strong><br />
Hilbert-muuntimen impulssivaste voidaan määrittää taajuusvasteesta käänteisellä<br />
Fourier-muunnoksella. Tulokseksi saadaan<br />
h Hilbert<br />
( t)<br />
=<br />
1<br />
πt<br />
Mielivaltaisen signaalin x(t) Hilbert-muunnos voidaan määrittää Hilbert-muuntimen<br />
impulssivasteen <strong>ja</strong> signaalin konvoluutiona<br />
∞<br />
1 1<br />
xˆ(<br />
t)<br />
= x(<br />
t)<br />
∗ hHilbert ( t)<br />
= ∫ x(<br />
τ )* dτ<br />
= ∫<br />
π ( t −τ<br />
) π<br />
−∞<br />
x(<br />
τ )<br />
dτ<br />
t −τ<br />
Tässä on erityisesti huomattava, että Hilbert-muunnettu signaali on a<strong>ja</strong>n funktio, eli<br />
Hilbert-muunnos ei tuota signaalista taajuustason esitystä. Hilbert-muunnoksen<br />
käänteismuunnos määritellään muodossa<br />
∞<br />
1 1<br />
x(<br />
t)<br />
= −xˆ(<br />
t)<br />
∗ hHilbert ( t)<br />
= − ∫ xˆ(<br />
τ )* dτ<br />
= − ∫<br />
π ( t −τ<br />
) π<br />
−∞<br />
©Jyrki Laitinen TL5231 Signaaliteoria S2004 37<br />
∞<br />
−∞<br />
∞<br />
−∞<br />
xˆ(<br />
τ )<br />
dτ<br />
t −τ