41t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 19 行

Rp=3;Rs=20;
N=15;n=[0:1:N-1];%只需要改变N的值
wc=pi/4;
hd=ideal_lp(wc,N);
w_han=(hanning(N))';
h=hd.*w_han;
[db,mag,pha,grd,w]=freqz_m(h,1);
subplot(211);plot(w/pi,db);xlabel('w/pi');ylabel('db');
hold on;
subplot(212);plot(w/pi,pha);xlabel('w/pi');ylabel('pha');
hold on;
N=33;n=[0:1:N-1];%只需要改变N的值
wc=pi/4;
hd=ideal_lp(wc,N);
w_han=(hanning(N))';
h=hd.*w_han;
[db,mag,pha,grd,w]=freqz_m(h,1);
subplot(211);plot(w/pi,db,'r');xlabel('w/pi');ylabel('db');
subplot(212);plot(w/pi,pha,'r');xlabel('w/pi');ylabel('pha');