bpf.m

来自「用MATLAB实现的信号处理！ 低通滤波器高通滤波器可以输出方波！」· M 代码 · 共 21 行

fr1=700;
fp1=900;
fp2=1100;
fr2=1300;
fs=1800;
tr_width=min((fp1-fr1),(fr2-fp2))/2;
N= ceil(8/(tr_width*2))+1;
n=0:1:N-1;
wc1=(fr1+fp1)*pi;
wc2=(fr2+fp2)*pi;
hd=ideal_bp1(wc1,wc2,N);
w_ham=(hamming(N))';
h=hd.*w_ham;
[db,mag,pha,w]=freqz_m(h,[1]);
subplot(2,2,1);
stem(n,h);
title('带通滤波器');
subplot(2,2,2);
plot(w/pi,db);
title('带通滤波器幅度响应(dB)');