Data	Wykonawcy	Temat
5-12-95	Gabryjela Kremser Marek Macioszek Janusz Sychla	Obróbka sygnałów

1.Przykłady określenia charakterystyk sygnałów cyfrowych.

-przebieg ćwiczenia z pakietem Matlab:

f=[5,10,3,0,1,6,15,12,13,2,11,8,9,14,7,4]'f = 5 10 3 0 1 6 15 12 13 2 11 8 9 14 7 4 mean(f)ans = 7.5000 fd=dtrend(f)fd = -2.5000 2.5000 -4.5000 -7.5000 -6.5000 -1.5000 7.5000 4.5000 5.5000 -5.5000 3.5000

0.5000 1.5000 6.5000 -0.5000 -3.5000 mean(fd)ans = 0 plot(f) plot(fd)

1.1Funkcja autokorelacji

fc=xcorr(fd) fc=fc/16fc = 0.5469 -0.4687 -0.1094 2.5625 -0.0156 -3.4062 -3.6719 -5.3750 1.4219 2.6562 1.7656 -5.3125 -6.1406 -0.2812 5.2031 21.2500 5.2031 -0.2812 -6.1406 -5.3125 1.7656 2.6562 1.4219 -5.3750 -3.6719 -3.4062 -0.0156 2.5625 -0.1094 -0.4687 0.5469 plot(fc)

Szybka transmisja Fouriera FTT

ff=fft(f)/16ff = 7.5000 -1.2071 + 1.0412i 0.5000 + 1.2071i 0.2071 - 1.8066i -0.5000 - 0.5000i 0.2071 - 0.8066i 0.5000 + 0.2071i -1.2071 + 0.0412i 0.5000 -1.2071 - 0.0412i 0.5000 - 0.2071i 0.2071 + 0.8066i -0.5000 + 0.5000i 0.2071 + 1.8066i 0.5000 - 1.2071i -1.2071 - 1.0412i

Liczby zespolone sprzężone, częstotliwość zmienia się od O do Π , tak zwana częstotliwość Neyquista Π/T_s gdzie T_s =1 1.2 Szum biały e=rand(150,1); plot(e)

1.3 Funkcja autokorelacji e=randn(150,1); ec=xcorr(e)/150; plot(ec)

1.4 Skok jednostkowy u1=ones(150,1); plot(u1)

un=u1+1*e plot(un)

2. Dyskretyzacja obiektu ciągłego

num=2

den=[5 1]

den = 5 1

t=0:.1:15;

h=step(num,den,t),

plot(t,h)

g=impulse(num,den,t)

plot(g)

nyquist(num,den)

bode(num,den)

2.2 Model w przestrzeni stanu.

[a,b,c,d]=tf2ss(num,den)

a = -0.2000

b = 1

c = 0.4000

d = 0

[a,b,c,d]=tf2ss(num,den)

a = -0.2000

b = 1

c = 0.4000

d = 0

[ad,bd]=c2d(a,b,.1)

ad = 0.9802

bd = 0.0990

[num,den]=ss2tf(ad,bd,c,d)

num = 0 0.0396

den = 1.0000 -0.9802

hd=dstep(num,den,150);

plot(hd)

hd=dstep(num,den,15);

plot(hd)

[numd,dend]=ss2tf(ad,bd,c,d)

numd = 0 0.0396

dend = 1.0000 -0.9802

hd=dstep(numd,dend,15);

plot(hd)

---