|
면 고주파에서 불안정한 성능을 보인다.
여기서 한 가지 특이점은 IIR필터의 차수가 20일때이다. 예상대로라면 IIR필터의 차수가 10일때 보다 더 이상적인 필터의 모양과 가깝게 나와야 하나 0~1khz의 구간에서 반듯한 모양이 아닌 불규칙한 파형
|
|
|
IIR Butterworth LPF의 주파수 응답 곡선과 pole-zero diagram을 plot하여 보고서에 첨부하라. (차수를 2, 5, 10, 20로 변화시키면서)
차수(N)이 증가할 수 록 이상적인 LPF에 가까워 짐을 볼 수 있다.
Pole-Zero 다이어그램에서 보이듯이 차수(N)만큼 Pole 개수 Zero
|
|
|
IIR필터는 Elloptic filter를 이용한 설계
[N,wn]=ellipord(wp/pi,ws/pi,Rp,As);
[b,a]=ellip(N,Rp,As,wn);
db=20*log10(abs(freqz(b,a,-pi:dw:pi)));
subplot(2,1,2); plot(-1:dw/pi:1,db,'b-'); title('Optimal IIR (using Elliptic filter)');
Asd=-max(db(wsi+501:1:1001));
axis([0 1 -100 10]); grid on; xlabel('No
|
|
|
IIR Filter Design (fir2와 비슷함)
이외에도 MATLAB Toolbox에서 제공해주는 Filter Design Function들이 다수 있음. I. IIR Filter Order Selection
1. buttord
2. Cheb1ord / Cheb2ord
II. IIR Filter Design (Filter Coefficients - b, a Parameter 추출)
1. butter
2. cheby1 / cheby2
3. 기타
|
|
|
IIR filter>
#include <stdio.h>
#include "volume.h"
#define taps 7 //when BPF, 9
/* Global declarations */
float inp_buffer[BUFSIZE]; /* processing data buffers */
float out_buffer[BUFSIZE];
//edited
float x_buff[taps]={0.0,0.0,0.0,0.0,0.0,0.0,0.0}; //when BPF ,0.0,0.0};
float y_buff[taps]={
|
|
메뉴펼치기
