MATLAB filter, filtic 함수를 이용한 이산 시스템 응답

MATLAB 을 이용하여 아래와 같은 조건의 시스템에 대한 시스템 응답을 구해 보자.

 

<code>

 

a=[1 -1 1];     % 피드백 루프의 계수

b=[1 0 0];        % 포워드 루프의 계수

n=(0:40)';        % Sample index 값

x=inline('cos(2*pi*n/6).*(n>=0)','n'); % x(n) 을 인라인 함수로 생성

z_i=filtic(b,a,[1 2]); % filter initial condition 설정

y=filter(b,a,x(n),z_i); % filter 함수를 이용한 시스템 응답

stem(n,y), grid on, xlabel('n'), ylabel('y[n]')