clear;clc;close;
Rp=1;Rs=60;
N=49;n=[0:1:N-1];
wp1=0.35*pi;wp2=0.65*pi;ws1=0.2*pi;ws2=0.8*pi;
wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;
ideal_lp(wc1,N);
ideal_lp(wc2,N);
hd1=ideal_lp(wc1,N);
hd2=ideal_lp(wc2,N);
w_han=(kaiser(N,5.658))';
h=(hd2-hd1).*w_han
[db,mag,pha,grd,w]=freqz_m(h,1);
plot(w/pi,db);xlabel('w/pi');ylabel('db');hold on;
用matlab自带的fdatool吧。。。简单说来，就是在命令行打fdatool，然后回车
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024