y1.m

来自「基于matlab的数字信号处理经典实验代码及所有完整实验报告」· M 代码 · 共 9 行

clc
ws1=0.2*pi;wp1=0.35*pi;
wp2=0.65*pi;ws2=0.8*pi;
wd=min((wp1-ws1),(ws2-wp2));
wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;
w_bla=(blackman(N+1))';
N=ceil(11*pi/wd);
hd=ideal_lp(wc2,N+1)-ideal_lp(wc1,N+1);
h=hd.*w_bla