Filter의 주파수 응답 plotter 설계

목차

1. 실험과정
2. 결론

본문내용

matlab에서 제공하는 함수를 사용하지 않고 pole-zero plot을 하기 위해서 우선 필터의 분
모와 분자의 계수들을 입력으로 받아들이는 함수로 설정하였다. 받아들인 계수들을 roots라
함수를 이용하여 근을 구하였고 그 값이 각각 zero와 pole이기 때문에 zeros, poles라는
변수에 저장시켰다. 여기서 roots함수는 벡터의 순서에 따라 차수를 부여하여 근을 계산하
기 때문에 반드시 처음 계수를 입력 받을 때 순서대로 받아야 한다.
그 후 pole-zero plot을 하기 위해서 원을 plot해야 하는데 이를 위해 반지름을 r이라는 변수에 저장하고(단위 원이기 때문에) 전체 각은 0.00001간격으로 쪼개서 연속적인 값으로 보이도록 설정하였다. 이 각의 법위를 이용하여 x,y의 polar form을 각각 x,y에 저장하였고 이 x,y값을 가지고 plot을 하면 원이 나오게 된다(polar form을 한 값을 각각 축으로 삼았기 때문에) 여기에 hold on을 하고 zplane(zeros, poles)이라는 함수를 사용하여 위에서 구한 pole과 zero를 원 위에 겹쳐 그리도록 하면 pole-zero plot을 그릴 수 있다. zplane 함수 자체에 xlabel과 ylabel tatle이 지정되어 있기 때문에 추가적으로 설정할 필요 없이 위와 같은 pole-zero plot을 얻을 수 있었다.
(나)에서 구한 계수를 넣어서 결과를 확인해 보니 zero는 1,-1값을 가지고 pole은 원 내부에 있어서 앞서 말한 조건을 충족한다고 할 수 있다.

내장함수 pzmap과 freqz라는 함수를 사용하여 plot하면 위와 같은 그래프를 얻을 수 있다.
(나)에서 구한 a의 계수와 조건으로 주어진 b의 계수를 vector로 입력하고 agle은 각도로 -pi부터 pi까지 0.0001 간격으로 잘게 쪼개어 거의 연속인 값을 가지도록 설정하였다.
pole zero plot을 표현하기 위해서는 tf함수를 사용하여 변수들을 지정해 줘야 한다. tf함수는 필터의 분자 계수 vector, 분모 계수 vector, angle의 간격, z^-1로 그린다는 정보를 받아 re 값에 저장하게 된다.

참고 자료

없음