wp=0.4*pi; ws=0.6*pi; deltaw=ws-wp;% 过渡带宽的计算
N0 = ceil(6.6*pi/ deltaw)+1;%按海明窗计算所需的滤波器阶数N0
N=N0+mod(N0+1,2);%为了实现第一类偶对称滤波器,应使其长度N为奇数
wdham = (hamming(N))';%求窗函数
wc=(ws+wp)/2;%截止频率取为两边界频率的平均值
tao=(N-1)/2;
n=[0:N-1];
m=n-tao+eps;
hd=sin(wc*m)./(pi*m);%求理想脉冲响应
h=hd.*wdham; %设计的脉冲响应为理想脉冲响应与窗函数乘积
[H,w]=freqz(h,[1],1000,'whole');
H=(H(1:1:501))'; w=(w(1:1:501))';%对设计结果进行检验
mag=abs(H);
db=20*log10((mag+eps)./max(mag));
pha=angle(H);
grd=grpdelay(h,[1],w);
dw=2*pi/1000;
Rp=-(min(db(1:wp/dw+1)));%检验通带波动
As=-round(max(db(ws/dw+1:501)));%检验最小带阻衰减
%绘图
n=0:N-1;
subplot(2,2,1);stem(n,hd,'.');title('理想脉冲响应')
axis([0 N-1 -0.2 0.3]); ylabel('hd(n)'); text(N+1,-0.1,'n')
subplot(2,2,2);stem(n,wdham,'.');title('海明窗')
axis([0 N-1 0 1.1]); ylabel('w(n)'); text(N+1,0,'n')
subplot(2,2,3); stem(n,h,'.');title('实际脉冲响应')
axis([0 N-1 -0.2 0.3]); xlabel('n'); ylabel('h(n)')
sbplot(2,2,4);plot(w/pi,db);title('幅度响应(单位: db)');grid
axis([0 1 -100 10]);xlabel('频率(单位;pi)'); ylabel('分贝')
set(gca,'XTickMode','manual','XTick',[0,0.4,0.6,1])
set(gca,'XTickMode','manual','XTick',[-50,0])
set(gca,'YTickLabelMode','manual','YTickLabels',['50';' 0'])
set(gcf,'color','w');
运行出现下面的怎么办求大神
??? Error: File: C:\matlab7.0\work\Untitled.m Line: 23 Column: 48
Missing variable or function.
N0 = ceil(6.6*pi/ deltaw)+1;%按海明窗计算所需的滤波器阶数N0
N=N0+mod(N0+1,2);%为了实现第一类偶对称滤波器,应使其长度N为奇数
wdham = (hamming(N))';%求窗函数
wc=(ws+wp)/2;%截止频率取为两边界频率的平均值
tao=(N-1)/2;
n=[0:N-1];
m=n-tao+eps;
hd=sin(wc*m)./(pi*m);%求理想脉冲响应
h=hd.*wdham; %设计的脉冲响应为理想脉冲响应与窗函数乘积
[H,w]=freqz(h,[1],1000,'whole');
H=(H(1:1:501))'; w=(w(1:1:501))';%对设计结果进行检验
mag=abs(H);
db=20*log10((mag+eps)./max(mag));
pha=angle(H);
grd=grpdelay(h,[1],w);
dw=2*pi/1000;
Rp=-(min(db(1:wp/dw+1)));%检验通带波动
As=-round(max(db(ws/dw+1:501)));%检验最小带阻衰减
%绘图
n=0:N-1;
subplot(2,2,1);stem(n,hd,'.');title('理想脉冲响应')
axis([0 N-1 -0.2 0.3]); ylabel('hd(n)'); text(N+1,-0.1,'n')
subplot(2,2,2);stem(n,wdham,'.');title('海明窗')
axis([0 N-1 0 1.1]); ylabel('w(n)'); text(N+1,0,'n')
subplot(2,2,3); stem(n,h,'.');title('实际脉冲响应')
axis([0 N-1 -0.2 0.3]); xlabel('n'); ylabel('h(n)')
sbplot(2,2,4);plot(w/pi,db);title('幅度响应(单位: db)');grid
axis([0 1 -100 10]);xlabel('频率(单位;pi)'); ylabel('分贝')
set(gca,'XTickMode','manual','XTick',[0,0.4,0.6,1])
set(gca,'XTickMode','manual','XTick',[-50,0])
set(gca,'YTickLabelMode','manual','YTickLabels',['50';' 0'])
set(gcf,'color','w');
运行出现下面的怎么办求大神
??? Error: File: C:\matlab7.0\work\Untitled.m Line: 23 Column: 48
Missing variable or function.