采用窗函数法（Kaiser窗）设计一个FIR数字低通滤波器,说我14行的 'ideal'没有定义。
技术指标是：
要求：
采样频率为8kHz；
通带：0Hz~1kHz，带内波动小于5%；
阻带：1.5kHz，带内最小衰减：Rs=40dB
程序是：
% 采样频率为8kHz；
% 通带：0Hz~1kHz，带内波动小于5%； wp=0.、25pi
% 阻带：1.5kHz，带内最小衰减：Rs=40dB。 wst=0.375pi
%
clc
clear
Rs=40;
Wp=0.25*pi; %根据通带：0Hz~1kHz，带内波动小于5%； 得 wp=0.125pi
Wst=0.375*pi; % 阻带：1.5kHz，带内最小衰减：Rs=40dB。 得wst=0.1875pi
dert_w=Wst-Wp;
%N=ceil((Rs-7.95)*2*pi/(14.36*dert_w)+1);
N=ceil((10*pi/dert_w)+1);
beta=0.5842*(Rs-21)^0.4+0.07886*(Rs-21);
hd=ideal((Wst-Wp)/2,N); %滤波器在时域系统的冲击响应
B=kaiser(N,beta); %凯泽窗
h=hd.*(B)'; %加窗后
[H,m]=freqz(h,[1],1024,'whole'); %获取频率响应
mag=abs(H); %幅值
db=20*log10((mag+eps)/max(mag)); %分贝数
pha=angle(H); %相位
%绘图
w=m/pi
figure(1);
subplot(2,2,1);
stem(hd);
xlabel('n');
ylabel('hd');
title('滤波器时域');
subplot(2,2,2);
plot(w,mag);
xlabel('w');
ylabel('h');
title('加窗后幅度响应');
subplot(2,2,3);
plot(w,db);
xlabel('w');
ylabel('db');
title('分贝数');
axis([0 1 -100 0]);
subplot(2,2,4);
plot(w,pha);%实际低通滤波器单位脉冲响应
xlabel('w');
ylabel('相位');
title('相频响应');
axis([0 1 -4 4]);
ideal没定义应该是要新建一个ideal函数。我是matlab小菜鸟，希望哪位大神帮帮忙。
% Ideal LowPass filter computation
% --------------------------------
% [hd] = ideal_lp(wc,M)
%hd= ideal impulse response between 0 to M-1
%wc= cutoff frequency in radians
% M= length of the ideal filter
%
alpha = (M-1)/2; n = 0:1:(M-1);
m=n- alpha; fc = wc/pi; hd = fc*sinc(fc*m);
end
注：此处的wc=(Wst-Wp)/2本回答被网友采纳
窗函数法FIR低通数字滤波
3.窗函数FIR滤波器 式(7-4-3)至式(7-4-8)窗函数都满足关于中心偶对称的线性相位滤波器的约束要求,结合式(7-4-1)至式(7-4-2)可以得到相应窗函数的FIR低通数字滤波器函数(郭志宏,罗锋,等,2007): 航空重力勘探理论方法及应用 航空重力勘探理论方法及应用 用该滤波器窗口对时间域(或空间域)长度为M的数据...
用窗函数法设计FIR数字滤波器
2-1用窗函数法设计一个线性相位FIR低通滤波器，并满足性能指标:通带边界频率 Wp=0.5*pi,阻带边界频率Ws=0.66*pi,阻带衰减不小于40dB,通带波纹不大于3dB。选择汉宁窗。代码：wp =0.5*pi;ws=0.66*pi;wdelta =ws-wp;N= ceil(8*pi\/wdelta)if rem(N,2)==0 N=N+1;end Nw =N;w...
用窗函数法设计一个低通滤波器,阶数为50 截止频率为 [0 0.4 0.5 1...
用窗函数法设计一个低通滤波器,阶数为50 截止频率为 [0 0.4 0.5 1] 相应的幅值为 [1 1 0 0] 5 我来答 分享 微信扫一扫 新浪微博 QQ空间 举报 浏览22 次 可选中1个或多个下面的关键词,搜索相关资料。也可直接点“搜索资料”搜索整个问题。 函数 低通滤波器 截止频率 设计 幅值 搜索资料 本地...
如何利用窗函数设计fir滤波器
Matlab中，函数fir1（）和fir2（）利用加窗傅里叶级数法设计FIR滤波器。函数fir1（）用来设计传统的LP（低通）、HP（高通）、BP（带通）、BS（带阻）和多频带FIR滤波器；而函数fir2（）用来设计具有任意幅度响应的的FIR滤波器。一、函数fir1的各种形式如下：b = fir1(N,Wn)b = fir1(N,Wn,'f...
1、如何用窗函数法设计线性相位低通滤波器?
1、如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器，写出设计步骤。答：将模拟频率转换成数字频率，确定理想滤波器 的特性；由 求出 ；选择适当的窗函数，并根据线性相位条件确定窗函数的长度N；在MATLAB中，可由w=boxcar(N)（矩形窗）、w=hanning(N)（...
急!!!用窗函数法设计FIR滤波器的主要特点是什么?
没有输出到输入的反馈，但有些结构中（例如频率抽样结构）也包含有反馈的递归部分。设FIR滤波器的单位冲激响应h (n)为一个N点序列，0 ≤ n ≤N —1，则滤波器的系统函数为H(z)=∑h(n)*z^-k。就是说，它有（N—1）阶极点在z = 0处，有（N—1）个零点位于有限z平面的任何位置。
用汉宁窗函数法设计一个数字FIR低通滤波器,要求通带边界频率为400Hz...
要求通带边界频率为400Hz，阻带边界频率为500Hz，通带最大衰减1dB，阻带最小衰减40dB，抽样频率为2000Hz，用MATLAB画出幅频特性，画出并分析滤波器传输函数的零极点；信号经过该滤波器，其中300Hz，600Hz，滤波器的输出是什么？用Matlab验证你的结论并给出的图形。
...列出窗函数法设计列线性相位FIR数字低通滤波器的设计步骤?
如果要使用窗函数法来设计列线性相位FIR数字低通滤波器，可以按照以下步骤进行：确定滤波器的通带边界频率Wp和阻带边界频率Ws，以及通带和阻带中的最大响应衰减值Ap和As。使用频率响应规格，利用最小化差值法（Parks-McClellan算法）或拉普拉斯变换，求出滤波器的最小阶数N和对应的系数h。为了实现列线性...
数字滤波器设计
采用窗函数法(分别用汉宁窗、哈明窗、布莱克曼窗函数)设计满足要求的FIR低通滤波器,求出h(n),并画出幅频响应损耗函数曲线.用频率采样法设计满足要求的FIR低通滤波器,求出h(n),并画出幅频响应损耗函数曲线。具体内容如下:(1)设计巴特沃斯模拟低通滤波器,求出Ha(s)的分子、分母多项式系数B和A,并画出幅频...
哪位能发个数字语音滤波器的发展历史和现状啊>高分求助.
fir数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题,设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。在对滤波器实际设计时,整个过程的运算量是很大的。例如利用窗函数法【2】设计m阶fir低通滤波器时,首先要根据(1)式计算出理想低通滤波器的...
