plzhi4 发表于 2005-12-21 18:07

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

本帖最后由 wdhd 于 2016-3-14 15:05 编辑

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=;
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;
=freqz_m(h,);
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();xlabel('n');ylabel('gd(n)');
subplot(2,3,2);stem(n,w_box);title('矩形窗')
axis();xlabel('n');ylabel('w(n)');
subplot(2,3,3);stem(n,h);title('实际脉冲响应')
axis();xlabel('n');ylabel('h(n)');
subplot(2,3,4);plot(w/pi,db);axis();
title('幅度响应(单位:dB)');grid;

plzhi4 发表于 2005-12-21 18:09

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

songzy41 发表于 2006-8-12 07:27

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

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=;
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;
=freqz_m(h,);
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();xlabel('n');ylabel('gd(n)');
subplot(2,2,2);stem(n,w_bla);title('Blackman窗')
axis();xlabel('n');ylabel('w(n)');
subplot(2,2,3);stem(n,h);title('实际脉冲响应')
axis();xlabel('n');ylabel('h(n)');
subplot(2,2,4);plot(w/pi,db);axis();
title('幅度响应(单位:dB)');grid;

计算结果有：
页: [1]
查看完整版本: 求助:请帮我看一下这个程序(急)