求助:请帮我看一下这个程序(急)

plzhi4

clear;wp1=0.2*pi;ws1=0.25*pi;ws2=0.30*pi;wp2=0.35*pi;
rp=2;rs=50;
tr_width=min((ws1-wp1),(wp2-ws2));
M=ceil(6.6*pi/tr_width)+1
n=[0:1:M-1];
wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;
hd=ideal_lp(wc1,M)+ideal_lp(pi,M)-ideal_lp(wc2,M);
w_bla=(blackman(M));
h=hd.*w_bla;
[db,mag,pha,grd,w]=freqz_m(h,[1]);
delta_w=2*pi/1000;
rp=-min(db(wp1/delta_w+1:1:wp2/delta_w))
rs=-round(max(db(ws2/delta_w+1:1:ws1/delta_w)))
subplot(2,3,1);stem(n,hd);title('理想脉冲响应')
axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('gd(n)');
subplot(2,3,2);stem(n,w_box);title('矩形窗')
axis([0 M-1 0 1.1]);xlabel('n');ylabel('w(n)');
subplot(2,3,3);stem(n,h);title('实际脉冲响应')
axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('h(n)');
subplot(2,3,4);plot(w/pi,db);axis([0 1 -100 10]);
title('幅度响应(单位:dB)');grid;

plzhi4

我在线等.谢谢各位了!!

songzy41

我把该程序作了修改如下可能是迟到的回复。

clear;wp1=0.2*pi;ws1=0.25*pi;ws2=0.30*pi;wp2=0.35*pi;
rp=2;rs=50;
tr_width=min((ws1-wp1),(wp2-ws2));
M=ceil(6.6*pi/tr_width)+1
n=[0:1:M-1];
wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;
hd=ideal_lp(wc1,M)+ideal_lp(pi,M)-ideal_lp(wc2,M);
w_bla=(blackman(M));
h=hd'.*w_bla;
[db,mag,pha,grd,w]=freqz_m(h,[1]);
delta_w=2*pi/1000;
rp=-min(db(wp1/delta_w+1:1:wp2/delta_w))
rs=-max(db(ws1/delta_w+1:1:ws2/delta_w))
subplot(2,2,1);stem(n,hd);title('理想脉冲响应')
axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('gd(n)');
subplot(2,2,2);stem(n,w_bla);title('Blackman窗')
axis([0 M-1 0 1.1]);xlabel('n');ylabel('w(n)');
subplot(2,2,3);stem(n,h);title('实际脉冲响应')
axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('h(n)');
subplot(2,2,4);plot(w/pi,db);axis([0 1 -100 10]);
title('幅度响应(单位:dB)');grid;

计算结果有：