• Synteza i analiza filtru NOI (górnoprzepustowego)

Mplik do kreślący charakterystyki:

wga=1.3888

fgc=260;

fp=1600;

Tp=1/fp;

wgc=fgc*(2*pi);

B=wga*tan(wgc*Tp/2);

L=[1.43 -2.88 1.43];

M=[(B^2+1.4256*B+1.5162) (2*B^2-2*1.5162) (B^2-1.4256*B+1.5162)];

[H W]=freqz(L,M, 512,fp);

plot(W, abs(H)); %char amp-czestotliwosciowa(widmowa);

grid on;

%char fazowa;

P=angle(H)*180/pi;

figure(2)

plot(W,P);

grid on;

figure(3)

phasedelay(L,M,512,fp); %char opoznienia fazy;

%odpowiedz na skok jednostkowy;

figure(4)

[o ts]= stepz(L,M);

plot(ts,o);

grid on;

%odpowiedz na impuls;

figure(5);

impz(L,M);

grid on;

Charakterystyka widmowa:

Charakterystyka fazowa:

Charakterystyka opóźnienia fazowego:

Chakaterystyka odpowiedzi na skok jednostkowy:

Chakaterystyka odpowiedzi na impuls:

• Analiza filtu SOI

fodc=350;

fp=1100+100*5;

f1=50;

N=fp/f1;

H=[zeros(1,7), exp(-j*(7:15)*pi*(N-1)/N), zeros(1,1),exp(j*(N-(17:25))*pi*(N-1)/N),zeros(1,6)];

figure(15)

stem(0:N-1,abs(H));

Charakterystyka odpowiedzi impulsowej:

impz(L,1);

Charakterystyka widmowa:

M=1;

[H W]=freqz(L,M, 512,fp);

plot(W, abs(H)); %char amp-czestotliwosciowa(widmowa);

grid on;

Charakterystyka widmowa po dodaniu okna blackmana:

M=1;

L=L.*blackman(N)

[H W]=freqz(L,M, 512,fp);

plot(W, abs(H)); %char amp-czestotliwosciowa(widmowa);

grid on;

Charakterystyka częstotliwościowa:

P=angle(H)*180/pi;

plot(W,P);

Opóźnienie fazowe

phasedelay(L,M,512,fp)

Charakterystyka odpowiedzi na skok:

[o ts]= stepz(L,M);

plot(ts,o)

• Filtracja

Sygnał zaszumiony

N=40;

fodc=400;

fp=2000;

t=(0:(1/fp):0.16);

signal=5*sin(100*pi*t)+2.5*sin(2*pi*(fodc+150)*t);

plot(t,signal)

grid on;

Furier

N=40;

fodc=400;

fp=2000;

t=(0:(1/fp):0.16);

signal=5*sin(100*pi*t)+2.5*sin(2*pi*(fodc+150)*t);

Le = length(t)-1;

signal1 = abs(fft(signal))/(Le/2);

ff = 0:(fp./Le):fp;

plot(ff,signal1)

grid on;

Filtr NOI

gr=10;

wga=0.7874;

fp=1000+gr*100;

fgc=200+gr*10;

wgc=2*pi*fgc;

Tp=1/fp;

A=wga*(1/tan(wgc*Tp/2));

L=[1 2 1];

M=[(1+1.73*A+A^2) (-2*A^2+2) (A^2-1.73*A+1)];

N=40;

fodc=400;

fp=2000;

t=(0:(1/fp):0.16);

signal=5*sin(100*pi*t)+2.5*sin(2*pi*(fodc+150)*t);

signal1_f = filter(L,M,signal);

plot(t,signal1_f)

grid on;

gr=10;

wga=0.7874;

fp=1000+gr*100;

fgc=200+gr*10;

wgc=2*pi*fgc;

Tp=1/fp;

A=wga*(1/tan(wgc*Tp/2));

L=[1 2 1];

M=[(1+1.73*A+A^2) (-2*A^2+2) (A^2-1.73*A+1)];

N=40;

fodc=400;

fp=2000;

t=(0:(1/fp):0.16);

signal=5*sin(100*pi*t)+2.5*sin(2*pi*(fodc+150)*t);

signal1_fw = abs(fft(signal1_f))/(Le/2);

plot(ff,signal1_fw)

grid on;

FILTR SOI

fp=2000;

fodc=400;

f1=50;

N=fp/f1 ;

deltaf=fp/N

H=[exp(-i*(0:8)*pi*(N-1)/N),zeros(1,23),exp(i*(N-(32:39))*pi*(N-1)/N)]

L=ifft(H);

signal1_f = filter(L,H,signal);

plot(t,signal1_f)

grid on;

fp=2000;

fodc=400;

f1=50;

N=fp/f1 ;

deltaf=fp/N

H=[exp(-i*(0:8)*pi*(N-1)/N),zeros(1,23),exp(i*(N-(32:39))*pi*(N-1)/N)]

L=ifft(H);

signal1_fw = abs(fft(signal1_f))/(Le/2);

plot(ff,signal1_fw)

grid on;

WNIOSKI