FIR滤波器设计源码程序 - matlab通信信号 - 谷速源码
下载频道> 资源分类> matlab源码> 通信信号> FIR滤波器设计源码程序

标题:FIR滤波器设计源码程序
分享到:

所属分类: 通信信号 资源类型: 文件大小: 2.35 KB 上传时间: 2016-01-28 20:45:44 下载次数: 14 资源积分:1分 提 供 者: xiaopeng2 FIR滤波器设计源码程序
内容:
FIR滤波器设计源码程序,程序员在编程的过程中可以参考学习使用,希望对IT程序员有用,此源码程序简单易懂、方便阅读,有很好的学习价值!
% 【例4-7】根据下列技术指标,设计一个FIR数字低通滤波器:
%            =0.3π, =0.4π, =3dB, =50dB
% 选择一个适当的窗函数,确定单位冲激响应,绘出所设计的滤波器的幅度响应。
% 根据窗函数最小阻带衰减的特性表4-2,可采用海明窗和布拉克曼窗可提供大于50dB的衰减。本例选择海明窗,其过渡带为6.6π/N,因此具有较小的阶次
% 例程4-13  利用海明窗设计Ⅱ型数字低通滤波器
clear all;
Wp=0.2*pi;
Ws=0.4*pi;
rp=3;
rs=50;
Rip=10^(-rp/20);
Atn=10^(-rs/20);
 
tr_width=Ws-Wp;                         %过渡带宽度
N=ceil(6.6*pi/tr_width)+1                    %滤波器长度
n=0:1:N-1;
Wc=(Ws+Wp)/2;                         %理想低通滤波器的截止频率
hd=ideal_lp1(Wc,N);                      %理想低通滤波器的单位冲激响应
w_ham=(hamming (N))';                    %海明窗
h=hd.*w_ham;                           %截取得到实际的单位脉冲响应
 
[db,mag,pha,grd,w]=freqz_m2(h,[1]);              %计算实际滤波器的幅度响应
 
delta_w=2*pi/1000;
Ap=-(min(db(1:1:Wp/delta_w+1)))            %实际通带纹波
As=-(max(db(Ws/delta_w+1:1:501)))     %实际阻带纹波
 
figure(1);
subplot(221);
stem(n,hd);
title('理想单位脉冲响应hd(n)');
 
subplot(222);
stem(n,w_ham);
title('海明窗w(n)');
 
subplot(223);
stem(n,h);
title('实际单位脉冲响应hd(n)');
 
subplot(224);
plot(w/pi,db);
title('幅度响应(dB)');
axis([0,1,-100,10]);
 
figure(2);
subplot(221);
plot(w/pi,mag);
title('幅频特性');
xlabel('w(/pi)');
ylabel('|H(jw)|');
axis([0,1,0,1.2]);
set(gca,'XTickMode','manual','XTick',[Wp/pi Ws/pi]);
set(gca,'YTickMode','manual','YTick',[Atn Rip]);
grid
 
subplot(222);
plot(w/pi,db);
title('幅频特性(db)');
xlabel('w(/pi)');
ylabel('dB');
axis([0,1,-100,5]);
set(gca,'XTickMode','manual','XTick',[Wp/pi Ws/pi]);
set(gca,'YTickMode','manual','YTick',[-rs -rp 0]);
grid
 
subplot(223);
plot(w/pi,pha/pi);
title('相频特性');
xlabel('w(/pi)');
ylabel('pha(/pi)');
axis([0,1,-1,1]);
set(gca,'XTickMode','manual','XTick',[Wp/pi Ws/pi]);
grid
 
subplot(224);
plot(w/pi,grd);
title('群延时');
xlabel('w(/pi)');
ylabel('Sample');
axis([0,1,0,50]);
set(gca,'XTickMode','manual','XTick',[Wp/pi Ws/pi]);
grid
 

文件列表(点击上边下载按钮,如果是垃圾文件请在下面评价差评或者投诉):

关键词: 滤波器 源码 程序

Top_arrow
回到顶部
联系方式| 版权声明| 招聘信息| 广告服务| 银行汇款| 法律顾问| 兼职技术| 付款方式| 关于我们|
网站客服网站客服 程序员兼职招聘 程序员兼职招聘
沪ICP备19040327号-3
公安备案号:沪公网安备 31011802003874号
库纳格流体控制系统(上海)有限公司 版权所有
Copyright © 1999-2014, GUSUCODE.COM, All Rights Reserved