c语言低通滤波参数调节,低通滤波器的参数设置
低通滤波器是容许低于截止频率的信号通过,但高于截止频率的信号不能通过的电子滤波装置。
低通滤波器参数:Fs=8000,fp=2500,fs=3500,Rp=1dB,As=30dB,其他滤波器可以通过与低通之间的映射关系实现。
%%模拟滤波器
%巴特沃斯——滤波器设计
wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;
[N,wc]=buttord(wp,ws,Rp,As,‘s’)%计算率波器的阶数和3dB截止频率
[B,A]=butter(N,wc,‘s’);%计算滤波器系统函数分子分母多项式
fk=0:800/512:8000;wk=2*pi*fk;
Hk=freqs(B,A,wk);
figure
plot(fk/1000,20*log10(abs(Hk)));
grid on,xlabel(‘频率(kHz)’),ylabel(‘幅度(dB)’)
TItle(‘巴特沃斯模拟滤波器’)
axis([0,4,-35,5])
%%
%切比雪夫I——滤波器设计
wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;
[N1,wp1]=cheb1ord(wp,ws,Rp,As,‘s’)%计算切比雪夫滤波器的阶数和通带边界频率
[B1,A1]=cheby1(N1,Rp,wp1,‘s’);%计算滤波器系统函数分子分母多项式
fk=0:800/512:8000;wk=2*pi*fk;
Hk=freqs(B1,A1,wk);figure,
plot(fk/1000,20*log10(abs(Hk)));
grid on,xlabel(‘频率(kHz)’),ylabel(‘幅度(dB)’)
TItle(‘切比雪夫I模拟滤波器’)
axis([0,4,-35,5])
%%
%切比雪夫II——滤波器设计
wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;
[N2,wso]=cheb2ord(wp,ws,Rp,As,‘s’)%计算切比雪夫滤波器的阶数和通带边界频率
[B2,A2]=cheby2(N1,Rp,wso,‘s’);%计算滤波器系统函数分子分母多项式
fk=0:800/512:8000;wk=2*pi*fk;
Hk=freqs(B1,A1,wk);figure,
plot(fk/1000,20*log10(abs(Hk)));
grid on,xlabel(‘频率(kHz)’),ylabel(‘幅度(dB)’)
TItle(‘切比雪夫II模拟滤波器’)
axis([0,4,-35,5])
%%
%椭圆——滤波器设计
wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30;
[N,wpo]=ellipord(wp,ws,Rp,As,‘s’)%计算滤波器的阶数和通带边界频率
[B,A]=ellip(N,Rp,As,wpo,‘s’);%计算滤波器系统函数分子分母多项式
fk=0:800/512:8000;wk=2*pi*fk;
Hk=freqs(B1,A1,wk);figure,
plot(fk/1000,20*log10(abs(Hk)));
grid on,xlabel(‘频率(kHz)’),ylabel(‘幅度(dB)’)
axis([0,4,-35,5]),TItle(‘椭圆模拟滤波器’)
%%
%数字滤波器
%脉冲响应法滤波器设计
fp=2500;fs=3500;Fs=8000;
wp=2*fp/Fs;ws=2*fs/Fs;%求归一化数字通带截止频率,求归一化数字阻带起始频率
deltaw=ws-wp;%求过渡带宽
N0=ceil(6.6/deltaw);%求窗口长度
N=N0+mod(N0+1,2); %确保窗口长度 N为奇数
n=N-1;%求出滤波器的阶数 n
wn=(ws+wp)/2; %求滤波器的截止频率
b=fir1(n,wn)%利用 fir1 函数求出滤波器的系数
[Hk,w] = freqz(b,1); % 计算频率响应
mag = abs(Hk); % 求幅频特性
db = 20*log10(mag/max(mag)); % 化为分贝值
dw =pi/512; %关于pi归一化
Rp = -(min(db(1:wp*pi/dw+1))) % 检验通带波动
As = -(max(db(ws*pi/dw+1:512))) % 检验最小阻带衰减
figure,plot(0:pi/511:pi,db),grid on
axis([0,4.0,-80,5]),title(‘数字滤波器——脉冲响应法’)
%%
%fir1窗函数法
fp=2500;fs=3500;Fs=8000;rs=30;
wp=2*fp*pi/Fs;ws=2*fs*pi/Fs;%求归一化数字通带截止频率,求归一化数字阻带起始频率
Bt=ws-wp;%求过渡带宽
alpha=0.5842*(rs-21)^0.4+0.07886*(rs-21);%计算kaiser窗的控制参数
M=ceil((rs-8)/2.285/Bt);%求出滤波器的阶数
wc=(ws+wp)/2/pi; %求滤波器的截止频率并关于pi归一化
hk=fir1(M,wc,kaiser(M+1,alpha))%利用 fir1 函数求出滤波器的系数
[Hk,w] = freqz(hk,1); % 计算频率响应
mag = abs(Hk); % 求幅频特性
db = 20*log10(mag/max(mag)); % 化为分贝值
db1=db‘;
figure,plot(0:pi/511:pi,db1),grid on
axis([0,4.0,-80,5]),title(’数字滤波器——fir1窗函数法‘)
%%
%频率采样法
fp=2500;fs=3500;Fs=8000;rs=30;
wp=2*fp*pi/Fs;ws=2*fs*pi/Fs;%求归一化数字通带截止频率,求归一化数字阻带起始频率
Bt=ws-wp;%求过渡带宽
m=1;alpha=0.5842*(rs-21)^0.4+0.07886*(rs-21);%计算kaiser窗的控制参数
N=ceil(m+1)*2*pi/Bt;%求出滤波器的阶数
N=N+mod(N+1,2);
Np=fix(wp/(2*pi/N));
Ns=N-2*Np-1;
Hk=[ones(1,Np+1),zeros(1,Ns),ones(1,Np)];
wc=(ws+wp)/2/pi; %求滤波器的截止频率并关于pi归一化
hk=fir1(M,wc,kaiser(M+1,alpha))%利用 fir1 函数求出滤波器的系数
[Hk,w] = freqz(hk,1); % 计算频率响应
mag = abs(Hk); % 求幅频特性
db = 20*log10(mag/max(mag)); % 化为分贝值
db1=db’;
figure,plot(0:pi/511:pi,db1),grid on
axis([0,4.0,-80,5]),title(‘数字滤波器——频率采样法’)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
%利用等波纹最佳逼近法设计FIR数字滤波器
Fs=8000;f=[2500,3500];m=[1,0];
rp=1;rs=30;
delta1=(10^(rp/20)-1)/(10^(rp/20)+1);delta2=10^(-rs/20);
rip=[delta1,delta2];
[M,fo,mo,w]=remezord(f,m,rip,Fs);%边界频率为模拟频率时必须加入采样频率
M=M+1;%估算的M直达不到要求,家1后满足要求
hn=remez(M,fo,mo,w);
[Hk,w] = freqz(hn,1); % 计算频率响应
mag = abs(Hk); % 求幅频特性
db = 20*log10(mag/max(mag)); % 化为分贝值
db1=db‘;
figure,plot(0:pi/511:pi,db1),grid on
axis([0,4.0,-80,5]),title(’数字滤波器——等波纹最佳逼近法‘)
c语言低通滤波参数调节,低通滤波器的参数设置相关推荐
- c语言低通滤波程序,一阶低通滤波器c语言
1. 一阶滤波算法的原理 一阶滤波,又叫一阶惯性滤波,或一阶低通滤波.是使用软件编程实现普通硬件RC低通滤波器的功能. 一阶低通滤波的算法公式为: Y(n)=αX(n) (1-α)Y(n-1) 式中: ...
- 用理想低通滤波器在频率域实现低通滤波、高通滤波。
文章目录 理想低通滤波 理想高通滤波 总结: 在matlab中,处理频域图像的过程大概为: 1.对图像进行傅里叶变换. 2.为了方便构建滤波器,将频域图像中心化(即将0频显示在图像中心) 3.设计所需 ...
- python低通滤波器_[Python图像处理] 二十三.傅里叶变换之高通滤波和低通滤波
[Python图像处理] 二十三.傅里叶变换之高通滤波和低通滤波 发布时间:2019-04-28 20:35, 浏览次数:2739 , 标签: Python 该系列文章是讲解Python OpenCV ...
- 一阶低通滤波器方程_一阶低通滤波原理详解
在汽车标定中,使用最多的滤波算法即低通滤波,很多朋友可能知道怎么标定,但是不清楚具体的原理,因此本文将介绍一阶低通滤波的原理.算法建模仿真和优缺点:一阶滤波算法的原理 一阶滤波,又叫一阶惯性滤波,或一 ...
- 图像 理想低通滤波_图像处理之滤波(下)
[toc]目录 一.常规滤波 低通 高通 带通 带阻 二.非局部均值滤波 三.维纳滤波 四.卡尔曼滤波 前言 所谓滤波,其实就是从混合在一起的诸多信号中提取出所需要的信号. 信号的分类: 确定型信号, ...
- 中值滤波_图像处理之低通滤波
在图像处理或者计算机视觉应用中,在正式对图像进行分析处理前一般需要一个预处理的过程.预处理是对图像作一些诸如降维.降噪的操作,主要是为后续处理提供一个体积合适的.只包含所需信息的图像.这里通常会用到一 ...
- # MPU6050之2阶低通滤波,滤波效果非常好
自从打算搞无人机之后,最终倒在了姿态的问题上,每次到了调参数PID,飞机都不稳.之后锁定了几个可能的因素,一是参数没做到最佳,二就是电机高转速下,6050受到震动的情况下,姿态乱了,导致PID无法 ...
- 实战低通滤波和卡尔曼滤波
目录 背景 所需硬件与软件 理论基础 一阶低通滤波 卡尔曼滤波 仿真验证 实际验证 总结 背景 滤波这个词对任何一个工科生都不会陌生,尤其是做控制或者信号方面的从业者和学生.我们不仅可以通过硬件滤波也 ...
- opencv 四 Mat的基本操作3(高通滤波、低通滤波、对比度调节)
图像滤波分为高通滤波和低通滤波,高通滤波用于求图形的边缘,低通滤波用于图像去噪.图像模糊化等.这里的频是指变化(相邻像素值的变化),高通滤波是指使变化大也就是图像的边缘)的通过(低通滤波是指使变化小( ...
- 数字图像处理学习笔记5:频率域滤波1(傅里叶频谱图,低通滤波-平滑,高通滤波-锐化)
文章目录 前言 一.傅里叶变换:傅里叶频谱图 二.低通滤波 1.理想低通滤波 2.布特沃斯低通滤波 3.高斯低通滤波 4.小结 三.高通滤波 1.理想高通滤波 2.布特沃斯高通滤波 3.高斯高通滤波 ...
最新文章
- ubuntu下制作u盘启动盘
- Python Scipy 科学计算库
- java垃圾回收之循环引用
- 在python中、处理的一切都是对象_Python 3+ 一切都是对象
- 问题 C: 【例2-3】围圈报数
- D3.js、echar.js 前端必备大数据技能
- Java NIO---1
- 目前服务器操作系统版本,目前服务器操作系统版本号
- ssms 缺少索引信息_搜索引擎整合营销
- 【模拟遍历】剑指offer:替换空格
- 林德物料搬运公司成功案例:基于功能需求开发的软件模型质量保障
- java编写个倒计时_怎么编写一个倒计时java程序?求具体步骤!
- 翻译:swift 5初始化 被忽略的Convenience便捷初始化、Required和继承
- 拼多多电商玩家如何利用软件机器人快速采集平台数据
- 一款强大的网站在线客服聊天系统:whisper搭建教程
- 设计模式六大原则详解
- Cannot resolve method ‘create‘ in ‘String‘
- 电脑误删的文件怎么恢复?分享90%的人都会的这2招
- 火焰检测网络(October, 2018,ICIP)
- 使用伪类来实现类似微信群聊的头像样式
热门文章
- 父级fixed_CSS-应用position:fixed时是否可以保持父元素的宽度?
- 动态数据的表格页面展示
- 示波器的带宽与采样率是什么关系
- 交换机端口mtu值最大_思科交换机基本配置命令全集-mtu值怎么设置
- bom实现方块移动_javascript利用键盘控制小方块的移动
- 如何测量多个变量之间的非线性关系
- 如何粘贴winedit编辑器console中的错误文本信息
- java中floor和ceil_floor() 和 ceil()函数
- python向网页上传excle、text文件
- Java中将List分组到Map中算法(可用于android联系人拼音分组)