FIR Filter using MATLAB

본문내용

1 .9];
b = fir1(N, stopband,'stop');
subplot(3,1,1); plot(fhz, H,ff, abs(freqz(b, 1, n)));
subplot(3,1,3); plot(t, filter(b, 1, s));
lab5 실험 1~3
fs=20000;
n=256;
ff=fs/(2*n)*(0:n-1);
N=10;
f=[0 .1 .1 .1 1];
H=[1 1 1 0 0 ];
fhz=f*fs/2;
b = fir1(N, 0.1);
t=(1:20000)/fs;
s=sin(2*pi*t*5) + sin(2*pi*t*15) + sin(2*pi*t*30);
subplot(3,1,1);plot(t,s);
subplot(3,1,2); plot(fhz, H, ff, abs(freqz(b, 1, n)));
legend('FIR LPF N=20');
subplot(3,1,3); plot(t, filter(b, 1, s));
실험4
fs=20000;
n=256;
ff=fs/(2*n)*(0:n-1);
N=100;
f=[0 .05 .05 .05 1];
H=[0 0 1 1 1];
fhz=f*fs/2;
b = fir1(N, 0.05, 'high');
t=(1:20000)/fs;
s=sin(2*pi*t*5) + sin(2*pi*t*15) + sin(2*pi*t*30);
subplot(3,1,1);plot(t,s);
subplot(3,1,2); plot(fhz, H, ff, abs(freqz(b, 1, n)));
legend('FIR HPF N=10');
subplot(3,1,3); plot(t, filter(b, 1, s));
실험5
fs=20000;
n=256;
ff=fs/(2*n)*(0:n-1);
N=10;
f=[0 .1 .1 .9 .9 1];
H=[0 0 1 1 0 0];
fhz=f*fs/2;
t=(1:200)/fs;
s=sin(2*pi*t*1000) + sin(2*pi*t*2000) + sin(2*pi*t*5000);
subplot(3,1,1);plot(t,s);
passband=[.1 .9];
b = fir1(N, passband);
subplot(3,1,2); plot(fhz, H,ff, abs(freqz(b, 1, n)));
subplot(3,1,3); plot(t, filter(b, 1, s));
실험6
fs=20000;
n=256;
ff=fs/(2*n)*(0:n-1);
N=10;
f=[0 .1 .1 .9 .9 1];
H=[1 1 0 0 1 1];
fhz=f*fs/2;
t=(1:200)/fs;
s=sin(2*pi*t*1000) + sin(2*pi*t*2000) + sin(2*pi*t*5000);
subplot(3,1,1);plot(t,s);
stopband=[.1 .9];
b = fir1(N, stopband,'stop');
subplot(3,1,2); plot(fhz, H,ff, abs(freqz(b, 1, n)));
subplot(3,1,3); plot(t, filter(b, 1, s));
fs=20000;
n=256;
ff=fs/(2*n)*(0:n-1);
N=4;
f=[0 .1 .1 .1 1];
H=[1 1 1 0 0 ];
fhz=f*fs/2;
b = fir1(N, 0.1);
t=(1:200)/fs;
s=sin(2*pi*t*500)+sin(2*pi*t*1000)+sin(2*pi*t*2000)+sin(2*pi*t*5000)+sin(2*pi*t*9500);
figure; plot(t,s);
subplot(5,2,1); plot(fhz, H, ff, abs(freqz(b, 1, n)));
legend('FIR LPF N=6');
subplot(5,2,2); plot(t, filter(b, 1, s));
N=6;
f=[0 .1 .1 .1 1];
H=[1 1 1 0 0 ];
fhz=f*fs/2;
b = fir1(N, 0.1);
t=(1:200)/fs;
s=sin(2*pi*t*500)+sin(2*pi*t*1000)+sin(2*pi*t*2000)+sin(2*pi*t*5000)+sin(2*pi*t*9500);
subplot(5,2,3); plot(fhz, H, ff, abs(freqz(b, 1, n)));
legend('FIR LPF N=6');
subplot(5,2,4); plot(t, filter(b, 1, s));
N=20;
f=[0 .1 .1 .1 1];
H=[1 1 1 0 0 ];
fhz=f*fs/2;
b = fir1(N, 0.1);
t=(1:200)/fs;
s=sin(2*pi*t*500)+sin(2*pi*t*1000)+sin(2*pi*t*2000)+sin(2*pi*t*5000)+sin(2*pi*t*9500);
subplot(5,2,5); plot(fhz, H, ff, abs(freqz(b, 1, n)));
legend('FIR LPF N=20');
subplot(5,2,6); plot(t, filter(b, 1, s));
N=50;
f=[0 .1 .1 .1 1];
H=[1 1 1 0 0 ];
fhz=f*fs/2;
b = fir1(N, 0.1);
t=(1:200)/fs;
s=sin(2*pi*t*500)+sin(2*pi*t*1000)+sin(2*pi*t*2000)+sin(2*pi*t*5000)+sin(2*pi*t*9500);
subplot(5,2,7); plot(fhz, H, ff, abs(freqz(b, 1, n)));
legend('FIR LPF N=50');
subplot(5,2,8); plot(t, filter(b, 1, s));
N=100;
f=[0 .1 .1 .1 1];
H=[1 1 1 0 0 ];
fhz=f*fs/2;
b = fir1(N, 0.1);
t=(1:200)/fs;
s=sin(2*pi*t*500)+sin(2*pi*t*1000)+sin(2*pi*t*2000)+sin(2*pi*t*5000)+sin(2*pi*t*9500);
subplot(5,2,9); plot(fhz, H, ff, abs(freqz(b, 1, n)));
legend('FIR LPF N=100');
subplot(5,2,10); plot(t, filter(b, 1, s));
실험3

키워드

FIR,   MATLAB,   필터