马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
clear;wp1=0.2*pi;ws1=0.25*pi;ws2=0.30*pi;wp2=0.35*pi;<BR>rp=2;rs=50;<BR>tr_width=min((ws1-wp1),(wp2-ws2));<BR>M=ceil(6.6*pi/tr_width)+1<BR>n=[0:1:M-1];<BR>wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;<BR>hd=ideal_lp(wc1,M)+ideal_lp(pi,M)-ideal_lp(wc2,M);<BR>w_bla=(blackman(M));<BR><FONT color=#0909f7>h=hd.*w_bla</FONT>;<BR>[db,mag,pha,grd,w]=freqz_m(h,[1]);<BR>delta_w=2*pi/1000;<BR>rp=-min(db(wp1/delta_w+1:1:wp2/delta_w))<BR>rs=-round(max(db(ws2/delta_w+1:1:ws1/delta_w)))<BR>subplot(2,3,1);stem(n,hd);title('理想脉冲响应')<BR>axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('gd(n)');<BR>subplot(2,3,2);stem(n,w_box);title('矩形窗')<BR>axis([0 M-1 0 1.1]);xlabel('n');ylabel('w(n)');<BR>subplot(2,3,3);stem(n,h);title('实际脉冲响应')<BR>axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('h(n)');<BR>subplot(2,3,4);plot(w/pi,db);axis([0 1 -100 10]);<BR>title('幅度响应(单位:dB)');grid;<BR><BR>有错,但我不知如何改
[ 本帖最后由 eight 于 2008-5-7 15:46 编辑 ] |