MATLAB에서 디지털 필터 시뮬레이션을 구현하는 방법

fp = 30;

fs = 35;

Fs = 100;

wp = 2*pi*fp/Fs;

ws = 2*pi*fs/Fs;

wp = tan(wp/2)

ws = tan(ws/2); 대역의 최대 감쇠는 0.5dB이고 정지 대역의 최소 감쇠는 40dB입니다.

);

set(gca,'XTickMode','manual','XTick' ,);

subplot(2, 2, 4);plot(w/pi, grd); title('그룹 지연')

xlabel('w(/pi )') ;ylabel('샘플');

axis();

set(gca,'XTickMode','manual','XTick',[0 0.2 0.3 1 ]); 그리드

작동 결과:

그림 2 버터워스 디지털 저역 통과 필터 진폭 주파수-위상 주파수 특성

1-3 디자인-버터워스 고역 디지털 필터에는 0.6*pi의 통과대역 차단 주파수, 1dB 이하의 통과대역 감쇠, 15DB의 저지대역 감쇠 및 샘플링 T=1이 필요합니다.

Wp=0.6*pi;

Ws=0.4*pi;

Ap=1;

As=15;

[N, wn]=buttord(Wp/pi, Ws/pi, Ap, As); 버터워스 필터 차수 및 차단 주파수 계산

주파수 변환을 통해 버터워스 설계 method Scotch 필터

[db, mag, pha, grd, w]=freqz_m(b, a) 디지털 필터 응답

plot(w, mag);

title('디지털 필터 진폭-주파수 응답|H(ej\Omega)|')

그림 3 버터워스 디지털 고역 통과 필터

2-1 사용 선형 위상 FIR 저역통과 필터를 설계하고 성능 지표를 충족하는 창 함수 방법: 통과대역 경계 주파수

Wp=0.5*pi, 저지대역 경계 주파수 Ws=0.66*pi, 저지대역 감쇠는 아님 40dB 미만, 통과대역 리플은 3dB를 넘지 않습니다. Hanning 창을 선택합니다.

코드:

wp =0.5*pi;

ws=0.66*pi;

wdelta =ws-wp;

N= ceil(8*pi/wdelta)

if rem(N, 2)==0

N=N 1;

끝

);