매트랩을 이용한 BPF 설계하기(bandpass filter) 주어진 스펙을 만족(fir1)

소개글

매트랩을 이용한 BPF 설계하기(bandpass filter) 주어진 스펙을 만족(fir1)

목차

I. 매트랩을 이용한 BPF 설계하기

II. test signal을 통과시켜보기

III. Perfectly sharp filter인가? 왜 그런가?

IV. 반성 및 고찰

본문내용

I. 매트랩을 이용한 BPF 설계하기
II. test signal을 통과시켜보기
<Matlab code 및 설명>
% 먼저 모든 변수를 지우고 시작
clear all
% 가장 먼저 신호를 sampling 해야 하는데, 이는 너무 크게 잡으면 필터설계가 어려워지므로, 샘플링은 Nyquist Rate의 2배정도가 적절하다. 따라서 나는 test signal의 최대주파수 20000의 2배인 40000정도를 그 값으로 했다.
fs=40000; % 샘플링 주파수를 40000으로 정함.

% t의 범위를 정하였다. 이는 test signal의 그것에 맞추어서 정해야 하는데, test signal의 최대주파수를 고려하여 범위를 정하였다.(너무 세밀하게 하면 plot 하는 시간이 너무 느려져서 적당하게 조정하였음)
t = 0:0.00005:0.03; % t=0:0.1:100; %Fig.2-2, Fig.3-2를 plot 할 때 대신 이것을 씀
% test signal
sig=sin(2*pi*100*t)+sin(2*pi*200*t)+sin(2*pi*500*t)+sin(2*pi*1000*t)+sin(2*pi*2000*t)+sin(2*pi*5000*t)+sin(2*pi*10000*t)+sin(2*pi*20000*t);

% Band-pass filter design
% W에서 필터가 cutoff해야 할 주파수를 정한다. 문제에서 주어진 것처럼 1000Hz와 5000Hz를 cutoff 하도록 정했다. M은 필터의 차수를 정하는 것으로 일반적으로 100~200정도로 하면 적절하다.

W = [1000 5000]/fs*2; %W=[850 5150]/fs*2; %Fig 1-1 plot
% 이 cutoff주파수는 6dB를 떨어뜨리는 주파수이므로 W =[850 5150]/fs*2; 정도로 적절히 수정하면 3dB주파수 1kHz와 5kHz를 갖는 필터를 설계할 수 있다. 바꾼 W로 그린 그래프는 아래에 따로 첨부함. (Fig.1-1)

참고 자료

없음