下面的函数为四种低通滤波器的实现

function [b, a]=analogfile(Wp, Ws, alphamax, alphamin, Wmax, ind)if ind==1,[N,Wn]=buttord(Wp, Ws, alphamax, alphamin, 's');[b,a]=butter(N, Wn, 's');elseif ind==2,[N, Wn]=cheb1ord(Wp, Ws, alphamax, alphamin, 's');[b, a]=cheby1(N, alphamax, Wn, 's');elseif ind==3,[N, Wn]=cheb2ord(Wp, Ws, alphamax, alphamin, 's');[b, a]=cheby2(N, alphamin, Wn, 's');else[N, Wn]=ellipord(Wp, Ws, alphamax, alphamin, 's');[b, a]=ellip(N, alphamax, alphamin, Wn, 's');endW=0:0.001:Wmax;H=freqs(b, a, W);Hm=abs(H);Ha=unwrap(angle(H));N=length(W);alpha1=alphamax * ones(1, N);alpha2 = alphamin * ones(1, N);subplot(221);plot(W, Hm);grid;axis([0  Wmax  0  1.1*max(Hm)]);subplot(222);plot(W, Ha);grid;axis([0  Wmax   1.1*min(Ha)  1.1*max(Ha)]);subplot(223);splane(b, a);figure(1);subplot(224);plot(W, -20*log10(abs(H)));hold on;plot(W, alpha1, 'r', W, alpha2, 'r');grid;axis([0 max(W) -0.1 100]);hold off;
end

其中使用到了:
unwrap 函数,推荐阅读:matlab之unwrap函数
freqs 函数,推荐阅读: matlab freqz-m,Matlab函数freqs和freqz

函数调用

clear all; clf
alphamax=0.1;
alphamin=60;
Wp=10;
Ws=15;
Wmax=25;
ind=4;%指定滤波器的种类
[b, a]=analogfile(Wp, Ws, alphamax, alphamin, Wmax, ind);



matlab 实现低通巴特沃斯滤波器、切比雪夫1型/2型滤波器 和 椭圆滤波器相关推荐

  1. Matlab语音信号去噪程序,使用低通巴特沃斯滤波器

    Matlab语音信号去噪程序,使用低通巴特沃斯滤波器. 1.读取一段歌曲的信号,绘制时域频域图,并播放. 2.添加正弦噪声: 3.设计巴特沃斯低通滤波器: 4.使用滤波器去除噪声,并画出时域频域图,播 ...

  2. 一个设计低通巴特沃斯数字滤波器的实例

    本人本科渣渣一个,前两天导师让我设计一个数字滤波器.由于本人基本没有数字信号处理基础,于是只能依靠百度和matlab,折腾了半天总算是摸索明白了.百度上有一些文章不靠谱,很容易误导别人,故在此发一篇博 ...

  3. java巴特沃斯滤波器编程_EMG信号的低通巴特沃斯滤波器

    使用matlab中自带的randn函数产生一组随机数,作为EMG信号,然后EMG信号的采样率为2048hz.这里随机数产生的随机数种子采用的机遇系统时钟的随机数种子.系统输入有两个,一个是仿真时间,单 ...

  4. 设计低通巴特沃斯数字滤波器(matlab)去除ECG基线漂移实例

    matlab代码************************** M=load ('datafile'); ecg=M(1:15000,1); 读取数据 fs=512; 采样率 fmaxd_1=5 ...

  5. 【图像去噪】均值+中值+高通低通+巴特沃斯+PCA+小波+维纳滤波图像去噪【含GUI Matlab源码 1705期】

    ⛄一.案例说明 三幅经典的图像(lena.peppers.barbara)分别加上了高斯噪声.乘性噪声.椒盐噪声.高斯白噪声:均值为0,方差为0.05.椒盐噪声:噪声密度0.10. 乘性噪声:这里为均 ...

  6. 5.6.2 低通滤波器(理想低通+巴特沃斯低通滤波器)

    1.低通滤波器 低通滤波是将频域图像中的高频部分滤除而通过低频部分.图像的边缘和噪声对应于频域图像中的高频部分,而低通滤波的作用即是减弱这部分的能量,从而达到图像平滑去噪的目的. 2.理想低通滤波器 ...

  7. 巴特沃斯、切比雪夫Ⅰ型、切比雪夫Ⅱ型和椭圆型滤波器的相同和不同之处

    数字滤波器是由数字乘法器.加法器和延时单元组成的一种算法或装置.数字滤波器的功能是对输入离散信号的数字代码进行运算处理,以达到改变信号频谱的目的. 数字滤波器有多种分类方法,每一种方法都从不同的侧面揭 ...

  8. 巴特沃斯、切比雪夫Ⅰ型、切比雪夫Ⅱ型、椭圆滤波器代码

    1. % %% 巴特沃斯带通滤波器 2. % %设计带通滤波器 3. % %巴特沃斯.切比雪夫I型.切比雪夫II型.椭圆型滤波器 4. % %wp和ws分别是通带和阻带的频率(截止频率).当wp和ws ...

  9. 滤波器基础05——巴特沃斯、切比雪夫与贝塞尔滤波器

    滤波器基础系列博客,传送门: 滤波器基础01--滤波器的种类与特性 滤波器基础02--滤波器的传递函数与性能参数 滤波器基础03--Sallen-Key滤波器.多反馈滤波器与Bainter陷波器 滤波 ...

最新文章

  1. PTA数据结构与算法题目集(中文)7-18
  2. SQL server的with的用法(一)
  3. 算法分析中递推式的一般代数解法
  4. 信息论4—无失真信源编码(非延长码,霍夫曼编码)
  5. AtCoder Regular Contest 064
  6. 电商如何“链”动新春,看聚划算助力“原年人”的好日子
  7. Java基础系列--static关键字
  8. javascript 操作服务器 文件,使用javascript读取服务器文件
  9. (云苍穹)表单插件常用代码
  10. vue实例练习和vue样式练习
  11. 如何创建一个微信公众号?
  12. 《写给大家看的设计书》(第四版)分享
  13. 微商卖货怎么引流?微商怎么找客源?
  14. 利用Python打造一个语音合成系统,听说你想要御姐音?
  15. 51单片机,点亮LED灯和闪烁
  16. web一阶段 day14
  17. 第十三届蓝桥杯 EDA 设计与开发科目 模拟试题(详细解读)
  18. 既然Talk is cheap, 那么就用代码教你如何进行正交设计
  19. Unity Steam_VR开发工具插件---VRTK 自带案例分析
  20. oracle查询当前时间前10分钟到当前时间的数据

热门文章

  1. 京东云服务器使用教程
  2. 深入理解JVM(十五)垃圾回收相关概念
  3. VTK:频域处理——低通滤波(理想+巴特沃斯)
  4. Arduino实验二十九 PCF8591数模转换传感器实验
  5. C语言的三种参数传递方式
  6. pyecharts去除小圆圈
  7. C++多线程——CreateThread
  8. 生产企业原材料订购与运输的研究-数据处理课程设计
  9. element之el-scrollbar
  10. Android 多窗口实现