基于matlab的FIR加窗滤波器设计
1.确定参数:
Wp 通带截止频率
Ws 阻带截止频率
As 衰减dB
Rp 纹波dB
2.根据参数确定理想 频率响应 (M (Wp + Ws )/2 )
双边响应
3.选窗
过滤带宽
阻带衰减
4.得出最终响应
时域相乘 频域卷积
5.观察各项参数是否达标(采用归一化后的数据)
衰减
幅频特性
相频
群延时 ??
···
Ps:观察频率响应的函数
这里a = 1 因为FIR可由差分方程描述
故分母为1
代码:
wp =0.2*pi;
ws = 0.3*pi;
tr_width = ws -wp;
M = ceil(6.6*pi/tr_width)+1;
n = [0:1:M-1];
wc = ((ws + wp)/2) %ideal LPF cutoff frequency
hd = ideal_lp(wc,M);
w_ham = (hamming(M))';
h = hd .* w_ham;
[db,mag,pha,grd,w] = freqz_m(h,[1]);
delta_w = 2*pi/1000;
Rp = -(min(db(1:1:wp/delta_w+1)));
As = -round(max(db(ws/delta_w+1:1:501)));
%plot
subplot(2,2,1);
stem(n,hd);
title('Ideal Impulse Response');
axis([0 M-1 -0.1 0.3]);
xlabel('n');
ylabel('hd(n)');
subplot(2,2,2);
stem(n,w_ham);
title(' Hamming Window');
axis([0 M-1 0 1.1]);
xlabel('n');
ylabel('w(n)');
subplot(2,2,3);
stem(n,h);
title('Actual Impluse Response');
axis([0 M-1 -0.1 0.3]);
xlabel('n');
ylabel('h(n)');
subplot(2,2,4);
plot(w/pi,db);
title('Magnitude Response in dB');grid
axis([0 1 -100 10]);
xlabel('frequency in pi units');
ylabel('Decibels');
带通:
wp1 =0.35*pi;
ws1 = 0.2*pi;
wp2 =0.65*pi;
ws2 = 0.8*pi;
tr_width = min((wp1 -ws1),(ws2 - wp2));
M = ceil(11*pi/tr_width)+1;
M
n = [0:1:M-1];
wc1 = ((ws1 + wp1)/2);
wc2 = ((ws2 + wp2)/2);
hd = ideal_lp(wc2,M)-ideal_lp(wc1,M);
w_ham = (blackman(M))';
h = hd .* w_ham;
[db,mag,pha,grd,w] = freqz_m(h,[1]);
delta_w = 2*pi/1000;
Rp = -(min(db(wp1/delta_w+1:1:wp2/delta_w+1)));
Rp
As = -round(max(db(ws2/delta_w+1:1:501)));
As
%plot
subplot(2,2,1);
stem(n,hd);
title('Ideal Impulse Response');
axis([0 M-1 -0.5 0.5]);
xlabel('n');
ylabel('hd(n)');
subplot(2,2,2);
stem(n,w_ham);
title(' blackman Window');
axis([0 M-1 0 1.1]);
xlabel('n');
ylabel('w(n)');
subplot(2,2,3);
stem(n,h);
title('Actual Impluse Response');
axis([0 M-1 -0.5 0.5]);
xlabel('n');
ylabel('h(n)');
subplot(2,2,4);
plot(w/pi,db);
title('Magnitude Response in dB');grid
axis([0 1 -150 10]);
xlabel('frequency in pi units');
ylabel('Decibels');
基于matlab的FIR加窗滤波器设计的更多相关文章
- IIR滤波器设计(调用MATLAB IIR函数来实现)
转载请注明文章来源 – http://blog.csdn.net/v_hyx ,请勿用于任何商业用途 对于滤波器设计,以前虽然学过相关的理论(现代数字信号处理和DSP设计),但一直不求 ...
- 基于MATLAB的GUI(Graphical User Interface)音频实时显示设计
摘要:本文章的设计主要讲基于matlab的gui音频实时显示设计,此次设计的gui相当于一个简洁的音乐播放器,界面只有”录音“和”播放“两个控件,哈哈,够简洁吧.通过”录音“按钮可以实现声音从电脑的声 ...
- Matlab滤波器设计(转)
滤波器设计是一个创建满足指定滤波要求的滤波器参数的过程.滤波器的实现包括滤波器结构的选择和滤波器参数的计算.只有完成了滤波器的设计和实现,才能最终完成数据的滤波. 滤波器设计的目标是实现数据序列的频率 ...
- FIR滤波器设计
FIR滤波器的优越性: 相位对应为严格的线性,不存在延迟失真,仅仅有固定的时间延迟: 因为不存在稳定性问题,设计相对简单: 仅仅包括实数算法,不涉及复数算法,不须要递推运算,长度为M,阶数为M-1,计 ...
- 信号、系统与滤波器设计(matlab)
0. 基本概念 AWG:Additive White(zero-mean) Gaussian,可加白噪声: AWGN:Additive White(zero-mean) Gaussian Noise ...
- 基于Matlab的MMSE的语音增强算法的研究
本课题隶属于学校的创新性课题研究项目.2012年就已经做完了,今天一并拿来发表. 目录: --基于谱减法的语音信号增强算法..................................... ...
- 基于MATLAB的多功能语音处理器
一.设计功能 录制音频,保存音频 对录制的语音信号进行频谱分析,确定该段语音的主要频率范围: 利用采样定理,对该段语音信号进行采样,观察不用采样频率(过采样.欠采样.临界采样)对信号的影响: 实现语音 ...
- fdatool的滤波器设计
作者:桂. 时间:2017-08-15 20:28:11 链接:http://www.cnblogs.com/xingshansi/p/7367738.html 前言 本文主要记录滤波器设计的基本流 ...
- 基于MATLAB的离散小波变换
申明,本文非笔者原创,原文转载自: 基于Matlab的离散小波变换 http://blog.sina.com.cn/s/blog_725866260100ryh3.html 简介 在 ...
随机推荐
- CGI/MIME/servlet术语解释
CGI→一种协议, 一种标准, 一种规范 使用CGI协议, 能够让用户访问某些动态资源的时候, 触发web服务器, 让web服务器根据CGI协议能够调用外部(web服务器外部)的程序来执行处理这个动态 ...
- 什么是Ajax和JSON,它们的优缺点
什么是Ajax??? 术语Ajax用来描述一组技术,它使浏览器可以为用户提供更为自然的浏览体验. Ajax它是“Asynchronous JavaScript + XML的简写” 定义Ajax: Aj ...
- netty实现远程调用RPC功能
netty实现远程调用RPC功能 依赖 服务端功能模块编写 客户端功能模块编写 netty实现远程调用RPC功能 PRC的功能一句话说白了,就是远程调用其他电脑的api 依赖 <dependen ...
- Python之二维数组N*N顺时针旋转90度
需求:把一个二维数组顺时针旋转90度,现实数据的替换. 比如把4*4的二维数组顺时针旋转90度 原始数据是一个嵌套列表:[['A', 'B', 'C', 'D'], ['A', 'B', 'C', ' ...
- 15 更多jQuery知识
jqueryUI 官网: https://jqueryui.com/ jqueryUI 中文网: http://www.jqueryui.org.cn/ jquery插件内容包含 官网demo: ht ...
- Jquery改变td内容为1的颜色
Jquery改变td内容为1的颜色<table id="tb" > <tr> <td val="1">1</td> ...
- 从cocos2d-x-2.x到cocos2d-x-3.x: lua项目配置
cocos2dx-x3.0的正式版出来也有一段时间了,现在最新的版本是到了3.2alpha,和2.x系列相比,能够找到的相关资料除了官网上的wiki,其他的也不见得多,遇到的一些和2.x的差异和问题在 ...
- sonarLint 插件配置sonarQube Server
Connected Mode You can bind Eclipse projects to a SonarQube project (supporting SonarQube servers 5. ...
- eclipse egit(分支管理 上)
这一章比较重要,讲述了Git比svn强大的地方,直接转载廖雪峰老师的文字,更好的理解 什么是分支 和 为什么分支git比svn做的更好 分支在实际中有什么用呢?假设你准备开发一个新功能,但是需要两周才 ...
- SqlServer数据导入到ORACLE
ORACLE中执行 select * from SYSTEM."employ_epl"