juxing=boxcar(n);；sanjiao=triang(n);%构造矩形窗；hanming=hamming(n);%构造汉宁；bulaikeman=blackman(n);%；kaize=kaiser(n);%构造凯塞窗；plot(x,juxing,'b.&#3；legend('矩形窗',&#3；运行结果如图2-1所示：；图
juxing=boxcar(n);
sanjiao=triang(n);
%构造矩形窗
%构造三角窗
hanming=hamming(n);
%构造汉宁窗
bulaikeman=blackman(n);
%构造布莱克曼窗
kaize=kaiser(n);
%构造凯塞窗
plot(x,juxing,'b.',x,sanjiao,'gx',x,hanming,'r+',x,bulaikeman,'cd',x,kaize,'k*');
legend('矩形窗','三角窗','汉明窗','布莱克曼窗','凯塞窗');
运行结果如图2-1所示：
图2-1 5种窗函数绘图比较
实际利用窗函数法进行FIR滤波器设计时，依据所给的技术指标一般需要经过以下几个步骤进行设计：
1． 给定理想的频率响应函数Hd(ejw)及技术指标；
2． 求出理想的单位抽样响应hd(n);
3． 根据阻带衰减选择窗函数w(n)
4． 根据过渡带宽度确定N值;
5． 求出所设计的FIR滤波器的单位抽样响应；
6． 计算频率响应，验算指标是否满足要求。
MATLAB作为一款优秀的数值计算软件，本身就内置了丰富的函数，其中便有用于通信仿真的一系列函数，并且MATLAB中还集成了通信设计的工具箱，不管是内置的函数，还是通信工具箱，均有专用于滤波器设计的工具，常用的主要有用函数法设计和用通信工具箱设计，下面分别予以介绍。
依据设计步骤编写
此种方法不依赖MATLAB中的滤波器设计函数，而是依据FIR滤波器的设计步骤自己求解理想滤波器的冲击响应，然后用窗函数对冲击响应进行截短，从而得到FIR滤波器，由于低通滤波器设计较为简单，因而可以通过两个低通的理想冲击响应函数相减得到理想带通的冲击响应，再通过窗函数对其截短，从而得到实际可行的FIR滤波器，这里以用汉明窗进行截短实现带通FIR滤波器，关键代码如下：完整代码在后面
wp1=0.3* ws1=0.45*
wp2=0.65* ws2=0.75*
tr_width=ws1-wp1;
%求过渡带宽度 %求得所需窗函数的长度 M=ceil(6.6*pi/tr_width)+1;
n=[0:1:M-1];
wc1=(ws1+wp1)/2; wc2=(ws2+wp2)/2;
hd=ideal_lp(wc2,M)-ideal_lp(wc1,M);
w_ham=(hamming(M))';
h=hd .* w_
程序运行结果：
%求截止频率 %求得理想带通的冲击响应 %得到长度为M的汉宁窗 %利用窗函数截短
图2-2 汉宁窗带通滤波器
利用MATLAB自带函数设计
利用原理说明介绍的fir1函数进行设计，这种设计方法只需要给出滤波器的阶数，截止频率，窗函数等参数，MATLAB即可自行完成设计，并可通过freqz函数查看滤波器的幅频响应和相频响应，已验证滤波器是否满足设计要求，下面给出利用fir1函数设计的不同窗函数的数字带通滤波器。
1. 利用矩形窗进行设计
%设定采样频率
fp1=4500;fp2=6500;
fs1=3000;fs2=7500;
%第一截止频率 %第二截止频率
%最小阻带衰减
Ws1=(fp1+fs1)/Ws2=(fp2+fs2)/
w=(fp1-fs1)/
M=ceil((As-7.95)/(14.36*w))
juxing=boxcar(M+1);
boxb=fir1(M,[Ws1,Ws2],juxing);
freqz(boxb,1,fs,fs);
%截止频率归一化处理 %求归一化过渡带 %计算所需滤波器的阶数 %生成长度为M+1的矩形窗 %生成矩形窗设计的fir滤波器 %绘制幅频和相频响应曲线
运行结果：
图2-3矩形窗fir滤波器幅频和相频响应曲线
从幅频响应上看，通带基本无波纹，阻带中波纹较大，因而阻带较不理想，相频响应曲线在通带内为直线，效果较好，信号失真小。
2. 利用三角窗进行设计
利用三角窗进行设计时，原理与矩形窗基本相同，只不过生成窗函数时采用triang（）函数生成三角窗，程序运行结果如下：
图2-4三角窗设计的fir滤波器幅频和相频响应曲线
3. 利用汉宁窗进行设计
利用汉宁窗进行设计时，原理与矩形窗基本相同，只不过生成窗函数时采用hamming（）函数生成三角窗，程序运行结果如下：
图2-5汉宁窗设计的fir滤波器幅频和相频响应曲线
4. 利用布莱克曼窗进行设计
利用布莱克曼窗进行设计时，原理与矩形窗基本相同，只不过生成窗函数时采用blackman（）函数生成三角窗，程序运行结果如下：
图2-6布莱克曼窗设计的fir滤波器幅频和相频响应曲线
5. 利用凯塞窗进行设计
利用凯塞窗进行设计时，滤波器的参数可以用以下函数求得：调用函数
[n,wn,bta,ftype]=kaiserord（f，a，dev，fs）
参数： 1）f为对应的归一化频率
2）a为由f指定的各个频带上的幅值向量，一般只有0和1，和
f长度关系为（2*a的长度）―2=（f的长度）
三亿文库包含各类专业文献、中学教育、高等教育、幼儿教育、小学教育、应用写作文书、文学作品欣赏、行业资料、外语学习资料、12MATLAB设计数字带通FIR滤波器的几种窗函数的比较等内容。　
　Matlab结合窗函数法设计数字带通FIR滤波器_IT/计算机_专业资料。matlab结合窗函数设计fir数字滤波器今日推荐 88份文档 2014全国高考状元联手分享状元笔记 ... 　窗函数法设计一个数字带通... 20页 20财富值 基于matlab的FIR数字滤波器......MATLAB 信号工具箱主要提供了以下几种窗函数,如表 1-1 所示:表 1-1 MATLAB... 　利用MATLAB仿真软件系统结合窗函数法设计一个数字带通FIR滤波器_高等教育_教育专区...几种常用的滤波方法: 1. 算术平均值法 2. 中值滤波法 1 《Matlab 课程设计... 　基于窗函数法的数字带通FIR滤波器设计_信息与通信_工程科技_专业资料。利用MATLAB仿真软件系统结合窗函数法设计一个数字带通FIR滤波器性能指标如下:通带下限截止频率... 　Matlab仿真软件系统结合窗函数法设计一个数字带通FIR滤波器_信息与通信_工程科技_专业资料。Matlab仿真软件系统结合窗函数法设计一个数字带通FIR滤波器武汉... 　滤波器的基本原理的理解 2、学习用 MATLAB 语言的窗函数法编写设计 FIR 数字...加窗线性相位 FIR 滤波器设计,可设计出标准的低通、带通、高通和带阻滤波 器... 　MATLAB窗函数法实现FIR的高通,带通和低通滤波器的程序_信息与通信_工程科技_专业...FIR 滤波器简介 FIR 滤波器的窗函数设计法及性能比较 数字滤波器是一种用来... 　利用 MATLAB 仿真软件系统结合窗函数法设计一个数字带通 FIR 滤波器 设计任务与...几种常用的滤波方法: 1. 算术平均值法 2. 中值滤波法 3. 滑动平均值法 4... 　Matlab结合窗函数法设计数字带通FIR滤波器_信息与通信_工程科技_专业资料。武汉理工...几种常用的滤波方法: 1. 算术平均值法 2. 中值滤波法 3 武汉理工大学《...Matlab实现数字FIR的高通_、带通、低通、带阻滤波器的程序
利用汉宁窗设计Ⅰ型数字高通滤波器clear&Wp=0.6*Ws=0.4*tr_width=Wp-Ws;&&&&&&&&&&&&&&&&&&&&&&&&&%过渡带宽度N=ceil(6.2*pi/tr_width)&&&&&&&&&&&&&&&&&&&&&&&&%滤波器长度n=0:1:N-1;Wc=(Ws+Wp)/2;&&&&&&&&&&&&&&&&&&&&&&&&&%理想低通滤波器的截止频率hd=ideal_hp1(Wc,N);&&&&&&&&&&&&&&&&&&&&&&%理想低通滤波器的单位冲激响应w_han=(hanning(N))';&&&&&&&&&&&&&&&&&&&&&%汉宁窗h=hd.*w_&&&&&&&&&&&&&&&&&&&&&&&&&&&&%截取得到实际的单位脉冲响应[db,mag,pha,w]=freqz_m2(h,[1]);&&&&&&&&&&&&&&%计算实际滤波器的幅度响应delta_w=2*pi/1000;Ap=-(min(db(Wp/delta_w+1:1:501)))&&&&&&&&&%实际通带纹波As=-round(max(db(1:1:Ws/delta_w+1)))&&&&&&&&&%实际阻带纹波subplot(221)stem(n,hd)title('理想单位脉冲响应hd(n)')subplot(222)stem(n,w_han)title('汉宁窗w(n)')subplot(223)stem(n,h)title('实际单位脉冲响应hd(n)')subplot(224)plot(w/pi,db)title('幅度响应(dB)')axis([0,1,-100,10])clear&Wp=0.6*Ws=0.4*tr_width=Wp-Ws;&&&&&&&&&&&&&&&&&&&&&&&&&%过渡带宽度N=ceil(6.2*pi/tr_width)&&&&&&&&&&&&&&&&&&&&&&&&%滤波器长度n=0:1:N-1;Wc=(Ws+Wp)/2;&&&&&&&&&&&&&&&&&&&&&&&&&%理想低通滤波器的截止频率hd=ideal_hp1(Wc,N);&&&&&&&&&&&&&&&&&&&&&&%理想低通滤波器的单位冲激响应w_han=(hanning(N))';&&&&&&&&&&&&&&&&&&&&&%汉宁窗h=hd.*w_&&&&&&&&&&&&&&&&&&&&&&&&&&&&%截取得到实际的单位脉冲响应[db,mag,pha,w]=freqz_m2(h,[1]);&&&&&&&&&&&&&&%计算实际滤波器的幅度响应delta_w=2*pi/1000;Ap=-(min(db(Wp/delta_w+1:1:501)))&&&&&&&&&%实际通带纹波As=-round(max(db(1:1:Ws/delta_w+1)))&&&&&&&&&%实际阻带纹波subplot(221)stem(n,hd)title('理想单位脉冲响应hd(n)')subplot(222)stem(n,w_han)title('汉宁窗w(n)')subplot(223)stem(n,h)title('实际单位脉冲响应hd(n)')subplot(224)plot(w/pi,db)title('幅度响应(dB)')axis([0,1,-100,10])基于切比雪夫一致逼近法设计FIR数字低通滤波器clear&&&&&&&&&&&&&&&f=[0&0.6&0.7&1]
分享这篇日志的人也喜欢
进来被骗不怪我
我又来了……
没守护，很绝望呀～
听歌就来看痴梦cc?
直播艰难，只想吃糖?
热门日志推荐
人人最热标签
北京千橡网景科技发展有限公司：
文网文[号··京公网安备号·甲测资字
文化部监督电子邮箱：wlwh@··
文明办网文明上网举报电话： 举报邮箱：&&&&&&&&&&&&
请输入手机号，完成注册
请输入验证码
密码必须由6-20个字符组成
下载人人客户端
品评校花校草，体验校园广场扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
下载作业帮安装包
扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
用matlab设计一个带阻滤波器,信号有5Hz,40Hz,50Hz,要求把40Hz的滤除,程序是这样的,滤波器和产生的信w1=38/100;w2=42/100;[B,A]=butter(1,[w1,w2],'stop');[h,w]=freqz(B,A);f=w/pi*100;figure(1),plot(f,20*log10(abs(h)));axis([0,80,-30,10]);xlabel('频率/Hz')ylabel('幅度/dB')title('频率响应曲线')fs=200;t=0:1/fs:1;x1=sin(2*pi*5*t)+sin(2*pi*50*t)+sin(2*pi*40*t);Y1=fft(x1,128);Y2=abs(Y1);f=(0:64)*100/64;figure(2),subplot(2,1,1),plot(f,Y2(1:65)),title('原始信号');y=filter(B,A,Y2);subplot(2,1,2),plot(y),title('滤波后信号');
有爱的小猴子29
扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
% 代码如下：fs=200;t=0:1/fs:1;x1=sin(2*pi*5*t)+sin(2*pi*50*t);plot(t,x1);title('sin(2*pi*5*t)+sin(2*pi*50*t)波形图');x2=sin(2*pi*5*t)+sin(2*pi*50*t)+sin(2*pi*40*t);plot(t,x2);title('sin(2*pi*5*t)+sin(2*pi*50*t)+sin(2*pi*40*t)波形图');fx = fft(x);fx1 = fftshift(fx);iCenter = 101;for i = 1 :200if abs(i-iCenter) >38 && abs(i - iCenter)
为您推荐：
其他类似问题
扫描下载二维码