Project2. Filter Design Using Bilinear Transformation
- 최초 등록일
- 2015.05.01
- 최종 저작일
- 2014.11
- 10페이지/
MS 워드
- 가격 2,500원
![할인쿠폰받기](/images/v4/document/ico_det_coupon.gif)
소개글
디지털신호처리 프로젝트2
목차
1. Butterworth Filter
2. Chebyshev type I Filter
3. Chebyshev type II Filter
4. Elliptic Filter
본문내용
Project2. Filter Design Using Bilinear Transformation
Given Specs
①0.99 ≤ |H(ejω)| ≤ 1.01, 0 ≤ |ω| ≤ 0.4π
②|H(ejω)| ≤ 0.001, 0.6π ≤ |ω| ≤ π
<중략>
wp=0.4*pi;
ws=0.6*pi;
rp=-20*log10(0.99); % 1-dp in dB
rs=-20*log10(0.001); % ds in dB
% (a) Prewarping with Td=1. Note that Omega=2 tan (w/2)/Td
Td=1; % Td=2
Omegap=2*tan(wp/2)/Td;
Omegas=2*tan(ws/2)/Td;
[N ,Omegac]=buttord(Omegap,Omegas,rp,rs,'s'); % Matlab solutions
[z, p, k]=butter(N,Omegac,'s'); % zero-pole-gain
[num ,den]=butter(N,Omegac,'s');
tf2sos(num,den) % H_c(s) in second-order-section
% (c) Bilinear transform
[bz ,az]=bilinear(num,den,1/Td); % default 1 for fs=1/Td
tf2sos(bz,az) % digital filter
% plot the pole zero diagram
figure(1);
zplane(bz,az)
% Other stuff below
[H ,w]=freqz(bz,az,512);
% plot (in dB) the magnitude response
figure(2);
plot(w/pi,20*log10(abs(H)));
axis([0 1 -100 20]);
xlabel('Normalized Frequency (/times/pi)') % dB plot
% plot the magnitude response (in linear scale)
figure(3);
plot(w/pi,abs(H)); % linear plot
xlabel('Normalized Frequency (/times/pi)')
% plot the group delay
figure(4);
grpdelay(bz,az)
참고 자료
없음