|
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 ==
|
|
|
1 .9];
b = fir1(N, stopband,'stop');
subplot(3,1,1); plot(fhz, H,ff, abs(freqz(b, 1, n)));
subplot(3,1,3); plot(t, filter(b, 1, s));
lab5 실험 1~3
fs=20000;
n=256;
ff=fs/(2*n)*(0:n-1);
N=10;
f=[0 .1 .1 .1 1];
H=[1 1 1 0 0 ];
fhz=f*fs/2;
b = fir1(N, 0.1);
t=(1:20000)/fs;
s=sin(2*pi*t*5) + sin(2*pi*t*
|
|
|
freqz의 함수를 이용한 결과와 비교, 분석하시오. (분석, 평가)
⇒ testba.m 파일을 실행한다. 우선 tf함수를 통하여 나온 결과를 확인해 본다.
transfer function이 나타났는데 이것은 주어진 계수들을 통하여 system function의 계수들을 대입한 것과 같은
|
|
|
freqz(b0);
figure(1) % #1 ideal filter, design filter plot
plot(f0/2,a0,w0/(2*pi),abs(h0)), grid;
legend('Ideal','firpm Design')
figure(2) % #2 필터 계수 stem
stem([0:n0],b0), grid;
figure(3) %#3 위상 plot
plot(w0/(2*pi),angle(h0)), grid;
%---------------양자화 에러체크--------%
bl=8;
m=1;
d1=abs(b0
|
|
|
freqz(b,a,128,fs);
진
폭
및
위
상
파
형
시
간
축
파
형
계속
구분
내 용(BSF)
주
파
수
축
파
형
Z-변환
매
트
랩
코
드
Zplane(b,a)
복
소
평
면
(4) LPF/BPF/HPF 사용 sample3.wav 잡음제거 과정 및 결과
- 저역 통과 필터 사용 잡음 제거
구분
내 용(LPF)
매
트
랩
|
|
메뉴펼치기
