该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

clear;close all;

fp=40;fs=150;Fs=1000;Rp=3;Rs=60;T=1/Fs;

W1p=2*tan(2*pi*fp*T/2)/pi;

W1s=2*tan(2*pi*fs*T/2)/pi;  %求归一化频率

[N,Wn]=buttord(W1p,W1s,Rp,Rs);%确定最小阶数和频率参数

[z,p,k]=buttap(N);%设定模拟低通原型的零极点增益参数

[bp,ap]=zp2tf(z,p,k);%将零极点增益转换成分子分母参数

[bp,ap]=butter(N,1,'s');%直接获取归一化低通原型

[bs,as]=lp2lp(bp,ap,Wn*pi*Fs);%将低通原型转换为模拟低通

[bz,az]=bilinear(bs,as,Fs);

sys=tf(bz,az,T);%给出传输函数

[H,W]=freqz(bz,az,512,Fs);

plot(W,20*log10(abs(H)));

figure(2)

n=0:100;

s=sin(2*pi*500*T)+sin(2*pi*7000*T); %含不同频率余弦信号

subplot(2,1,1);plot(T,s) %混频信号波形

title('输入信号')

xlabel('t');

ylabel('s');

s1=fft(s,100);

mag1=abs(s1);

f1=(0:length(s1)-1)'*fs/length(s1);

subplot(2,1,2);plot(f1,mag1) %混频信号频谱波形

title('输入信号频谱')

xlabel('w');

ylabel('db');

axis([0,0.8,0,40]);

figure(3)

y=filter(bz,az,s); %数字滤波函数输出

subplot(2,1,1);plot(T,y) %通过数字滤波器滤波后输出信号波形

title('输出信号')

xlabel('t');

ylabel('s');

y1=fft(y,100);

mag2=abs(y1);

f2=(0:length(y1)-1)'*fs/length(y1);

subplot(2,1,2);plot(f2,mag2) %输出信号频谱波形

title('输出信号频谱')

xlabel('w');

ylabel('db');

axis([0,0.8,0,40]);

这是一个基于双线性变换的巴特沃斯IIR低通滤波器,在运行时,只有滤波器图形,但输出输入的波形和频谱却没有了,不知哪里出问题,请大家帮忙

butter滤波器是iir吗_IIR低通滤波器求改相关推荐

  1. butter滤波器是iir吗_IIR数字滤波器的设计与滤波

    设计IIR滤波器,实现对存在加性干扰的时域离散信号进行滤波. 已知带加性干扰的信号用x(n)表示,x(n)=xs(n)+η(n),式中xs(n)是有用的信号,是一个0-0.2πrad的带限信号.η(n ...

  2. butter滤波器是iir吗_MATLAB IIR滤波器设计函数buttord与butter

    (4)Ws -阻带截止频率 标量| 二元向量 阻带截止频率,指定为标量或二元向量,值在0到1之间,其中1对应于归一化的奈奎斯特频率πrad / sample. 数据类型:单| 双 (5)Rp-通带纹波 ...

  3. butter滤波器是iir吗_学习随笔之IIR滤波器与FIR滤波器

    学习随笔之IIR滤波器与FIR滤波器 IIR滤波器(Infinite Impulse Response Digital Filter无限冲击响应数字滤波器)与FIR滤波器(Finite Impulse ...

  4. butter滤波器是iir吗_如何快速设计一个IIR滤波器

    在文章如何快速设计一个FIR滤波器(一)以及如何快速设计一个FIR滤波器(二)等文章中,我们讨论了如何设计FIR(Finite Impulse Response Filter),FIR有很多优点,比如 ...

  5. butter滤波器是iir吗_如何快速设计应用一个IIR滤波器

    在小平:如何快速设计应用一个FIR滤波器 中,我们讨论了如何设计一个FIR滤波器,接下来我们介绍IIR滤波器.和设计FIR滤波器一样,我们可以粗略的设计IIR滤波器(幅频响应不精确,设计简单),也可以 ...

  6. 使用matlab设计IIR巴特沃斯低通滤波器

    目的和要求 设计IIR巴特沃斯低通滤波器 将滤波器用于加噪声信号的处理 改变参数指标比较不同 实验结果与分析 (1)滤波器频率特性如图 输入信号特性如图: 经过滤波器之后信号特性如图: (2)改变指标 ...

  7. iir陷波滤波器 matlab,IIR数字滤波器设计50Hz陷波器(MATLAB代码)

    %% IIR陷波器设计 % 目的:设计一个陷波器阻带在50±1.5Hz以内,采样频率为400Hz的滤波器, % 并要求通带最大衰减为0.1dB,阻带最小衰减为60dB. clc; clear;clos ...

  8. fir 滤波参数的含义_详解FIR滤波器与IIR滤波器的具体区别

    姓名:高新诚       学号:17021223282 原文链接:http://mp.weixin.qq.com/s/tC_mdXIeusuggDV1tn3sOw [嵌牛导读]:在大学中学习到了FIR ...

  9. 简谈FIR滤波器和IIR滤波器的区别

    最近总是遇到很多大侠在问滤波器相关的问题,之前对这一方面接触不是很多,最近也是在学习一些这方面的知识,今天先和大侠简单聊一聊FIR滤波器和IIR滤波器的区别,后续等研究的差不多了,再更新有关滤波器的详 ...

最新文章

  1. 学习OGRE - 读取天龙八部scene文件
  2. HashMap 散列初体验
  3. 操作系统学习笔记-01-操作系统的概念(定义),功能和目标
  4. uni-app获取腾讯地图计算两经纬度的实际距离(可批量)
  5. web前端开发技术要求会什么
  6. git 学习1--查看全局配置
  7. html2canvas在iphone7 iphonex都生成不了
  8. python 调用 .netcore api_浅谈Python调用XBee的API来进行通讯
  9. 北大oj-1423- Big Number
  10. 高并发实战之幂等处理
  11. 如何添加、管理和删除PDF高亮文本?
  12. NB-IoT移远BC95使用小结
  13. R可视化绘图二-ggplot2
  14. wifi共享精灵2014最新版 v04.25.001 官方正式版
  15. excel两列数据对比找不同_Excel | 如何核对两列/行数据,找不同?
  16. 太乙超级计算机,从“启明”到“太乙”,南科大的超算发展之路
  17. 海康威视录像机刷linux系统,海康硬盘录像机怎么刷机?详尽图文刷机教程刷机工具亲测100%好用...
  18. 水洼数 DFS 蓝桥杯 java代码
  19. SAP 荣获「中国好公司」头衔
  20. 阿里旺旺自动回复工具开发一

热门文章

  1. 如何在Windows上安装Python
  2. 包管理器_包管理器的演变
  3. 公有云 私有云 混合云_您需要了解的有关混合云的知识
  4. Bootstrap3 模态对话框的方法
  5. Bootstrap3 动态模态对话框
  6. Bootstrap3栅格系统布局实例
  7. 学习经验浅谈——论做学习笔记的重要性
  8. 深度学习笔记(4) 浅层神经网络
  9. zkcli远程连接_高级框架第一天RPC:远程过程调用
  10. java 插入数据 主键_JDBC插入数据返回数据主键代码实例