|
freqz(b, 1, w);//%freqs, freqz함수는 H(s)와H(z)의 값들을 구하는데 사용된다
subplot(2,2,2);
plot(w, abs(H));
subplot(2,2,4);
plot(w, angle(H));//%함수 abs, angle은 복소함수H(s)와 H(Z)의 위상과 크기를 나타내는데 사용한다.
% Impulse response
b=[1];
a=[1 -0.9];
n=[-20:120];
x=[n ==
|
|
|
r를 convolution하였기 때문이다. 이차이는 그리 크지 않다. 이는 다음 page에서 MATLAB으로 구현된 spectrum에서 확인할 수 있다.
결과를 보다 확실히 하기 위해 주어진 입력 데이터를 MATLAB “freqz\"함수를 이용하여 주파수 성분을 분석하여 보았다.
주
|
|
|
스펙트럼을 나타낸다.
>> figure(4)
>> idx=[(16*up_rate):(32*up_rate)];
>> subplot(2,1,1);plot(t_up(idx),s_up(idx));
>>axis([t_up(min(idx)),t_up(max(idx)),-5,5]);
>> title(\'본래 신호\');
>>subplot(2,1,2);plot(t_up(idx),x_up(idx));
>>axis([t_up(min(idx)),t_up(
|
|
|
스펙트럼을 비교해 본 결과 성공적으로 두 성분이 제거되었음을 눈으로 확인할 수 있었다. 극점들이 모두 단위원 안에 있어 역시나 안정한 시스템임을 알 수 있다.
4) 실습 3에서 설계했던 코드에서 butter함수만 bandpass 필터로 바꿔주고 동일하
|
|
메뉴펼치기
