[매트랩]IIR필터, FIR 필터 설계.

소개글

DSP(Digital signal processing ) 디지털 신호처리 라는 과목에서
배우는 각종 IIR(Infinite impulse response) 필터와 FIR(Finite impulse response)필터를
매트랩으로 구현하였습니다. 레포트 마지막에 매트랩 코드를 첨부하였습니다.

IIR필터:버터워스필터 를 임펄스불변(Impulse invariance)방법과 쌍일차변환(bilinear transform) 방법을 이용하여 구현

FIR필터: 사각(rectengular)윈도우, Hann 윈도우 Hamming 윈도우 Black man 윈도우
Kaiser 윈도우, Optimal 윈도우 등을 통해서 구현

목차

1. 아래의 기준을 갖는 Butterworth 저역통과필터(LPF)를 설계한다. 임펄스 불변(impulse invariance) 방법과 쌍일차 변환 (bilinear transformation)을 이용하여 구현하고, 결과를 비교하라.

<윈도우 함수를 이용한 FIR 필터 설계>
2. 다음의 기준을 만족시키는 FIR 필터를 사각(rectangular) 윈도우, Hann 윈도우, Hamming 윈도우, Blackman 윈도우를 사용하여 설계해보라. 그리고, Kaiser 윈도우를 사용하여 설계해보고 결과를 비교하라.
3. 아래 기준을 만족시키는 FIR 필터를 Kaiser 윈도우를 이용하여 설계하라.

<최적의 FIR 및 IIR 필터 설계>
4. 아래 기준을 만족시키는 가장 효율적인 FIR 및 IIR 저역통과필터를 설계하고, 설계된 두 필터를 비교하라.

본문내용

소스중 일부.

%using bilinear transformation
T=1;
omegap = (2/T)*tan(wp/2);
omegas = (2/T)*tan(ws/2);
N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(omegap/omegas)));
omegac = omegap/((10^(Rp/10)-1)^(1/(2*N))); % cutoff freq.
wn = 2*atan((omegac*T)/2);
wn = wn/pi; % scaling
[b,a]=butter(N,wn);
dw=pi/500;
db=20*log10(abs(freqz(b,a,-pi:dw:pi)));
subplot(2,1,2); plot(-1:dw/pi:1,db,`b-`); title(`Butterworth Filter by bilinear transformation`);
axis([0 1 min(db) 10]); grid on; xlabel(`Freq.in pi units`); ylabel(`|H(e^{jw})| (dB)`)
text(0.9,-100,{[`N=`,num2str(N)]},`FontSize`,9,`BackgroundColor`,[0.8314 0.8157 0.7843],`EdgeColor`,[0 0 0]);

참고 자료

없음