matlab 实现低通巴特沃斯滤波器、切比雪夫1型/2型滤波器 和 椭圆滤波器
下面的函数为四种低通滤波器的实现
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型滤波器 和 椭圆滤波器相关推荐
- Matlab语音信号去噪程序,使用低通巴特沃斯滤波器
Matlab语音信号去噪程序,使用低通巴特沃斯滤波器. 1.读取一段歌曲的信号,绘制时域频域图,并播放. 2.添加正弦噪声: 3.设计巴特沃斯低通滤波器: 4.使用滤波器去除噪声,并画出时域频域图,播 ...
- 一个设计低通巴特沃斯数字滤波器的实例
本人本科渣渣一个,前两天导师让我设计一个数字滤波器.由于本人基本没有数字信号处理基础,于是只能依靠百度和matlab,折腾了半天总算是摸索明白了.百度上有一些文章不靠谱,很容易误导别人,故在此发一篇博 ...
- java巴特沃斯滤波器编程_EMG信号的低通巴特沃斯滤波器
使用matlab中自带的randn函数产生一组随机数,作为EMG信号,然后EMG信号的采样率为2048hz.这里随机数产生的随机数种子采用的机遇系统时钟的随机数种子.系统输入有两个,一个是仿真时间,单 ...
- 设计低通巴特沃斯数字滤波器(matlab)去除ECG基线漂移实例
matlab代码************************** M=load ('datafile'); ecg=M(1:15000,1); 读取数据 fs=512; 采样率 fmaxd_1=5 ...
- 【图像去噪】均值+中值+高通低通+巴特沃斯+PCA+小波+维纳滤波图像去噪【含GUI Matlab源码 1705期】
⛄一.案例说明 三幅经典的图像(lena.peppers.barbara)分别加上了高斯噪声.乘性噪声.椒盐噪声.高斯白噪声:均值为0,方差为0.05.椒盐噪声:噪声密度0.10. 乘性噪声:这里为均 ...
- 5.6.2 低通滤波器(理想低通+巴特沃斯低通滤波器)
1.低通滤波器 低通滤波是将频域图像中的高频部分滤除而通过低频部分.图像的边缘和噪声对应于频域图像中的高频部分,而低通滤波的作用即是减弱这部分的能量,从而达到图像平滑去噪的目的. 2.理想低通滤波器 ...
- 巴特沃斯、切比雪夫Ⅰ型、切比雪夫Ⅱ型和椭圆型滤波器的相同和不同之处
数字滤波器是由数字乘法器.加法器和延时单元组成的一种算法或装置.数字滤波器的功能是对输入离散信号的数字代码进行运算处理,以达到改变信号频谱的目的. 数字滤波器有多种分类方法,每一种方法都从不同的侧面揭 ...
- 巴特沃斯、切比雪夫Ⅰ型、切比雪夫Ⅱ型、椭圆滤波器代码
1. % %% 巴特沃斯带通滤波器 2. % %设计带通滤波器 3. % %巴特沃斯.切比雪夫I型.切比雪夫II型.椭圆型滤波器 4. % %wp和ws分别是通带和阻带的频率(截止频率).当wp和ws ...
- 滤波器基础05——巴特沃斯、切比雪夫与贝塞尔滤波器
滤波器基础系列博客,传送门: 滤波器基础01--滤波器的种类与特性 滤波器基础02--滤波器的传递函数与性能参数 滤波器基础03--Sallen-Key滤波器.多反馈滤波器与Bainter陷波器 滤波 ...
最新文章
- PTA数据结构与算法题目集(中文)7-18
- SQL server的with的用法(一)
- 算法分析中递推式的一般代数解法
- 信息论4—无失真信源编码(非延长码,霍夫曼编码)
- AtCoder Regular Contest 064
- 电商如何“链”动新春,看聚划算助力“原年人”的好日子
- Java基础系列--static关键字
- javascript 操作服务器 文件,使用javascript读取服务器文件
- (云苍穹)表单插件常用代码
- vue实例练习和vue样式练习
- 如何创建一个微信公众号?
- 《写给大家看的设计书》(第四版)分享
- 微商卖货怎么引流?微商怎么找客源?
- 利用Python打造一个语音合成系统,听说你想要御姐音?
- 51单片机,点亮LED灯和闪烁
- web一阶段 day14
- 第十三届蓝桥杯 EDA 设计与开发科目 模拟试题(详细解读)
- 既然Talk is cheap, 那么就用代码教你如何进行正交设计
- Unity Steam_VR开发工具插件---VRTK 自带案例分析
- oracle查询当前时间前10分钟到当前时间的数据