一、前言

1. 设计流程

2.系统频率响应

2.1频响图

系统函数 H 是一个复数，其图谱分为：幅度谱、相位谱

幅度谱

x轴：模拟频率f（数字频率w转化来）【单位：赫兹Hz】

y轴：|H1|幅度【一般用:20 * log10|H1|】【单位：分贝dB】

相位谱

x轴：模拟频率f（数字频率w转化来）【单位：赫兹Hz】

y轴：H1 的相位

2.2 各个频率转换关系

$\Omega = 2\pi f\$
$T = \frac{1}{{\mathop f\nolimits_s }}\$ 【 ${\mathop f\nolimits_s }\$ 为采样率】
$w = \Omega T = 2\pi f*\frac{1}{{\mathop f\nolimits_s }}\$
所以可以推出f： $f = \frac{{w\mathop f\nolimits_s }}{{2\pi }}\$

3.巴特沃斯滤波器简介

N：滤波器阶数

$\mathop \Omega \nolimits_{\rm{c}} \$ ：3dB截止频率

3.1特点

最大平坦性

在截止频率前较为平坦，这个平坦也保证了信号的原始值，不会因为滤波被衰减。巴特沃斯低通滤波器的通频带最大扁平效应使通频带的增益得到扁平优化。（由上图可知：N值越大，通频带越平坦）

幅频特性、相频特性都是单调递减

巴特沃斯低通滤波器是全极点滤波器，所有零点都在无穷远处；它的幅度随频率的增加而单调下降，且频率无穷大时它的幅度才衰减完。

3.2设计思路

幅度平方函数：表现了N 阶巴特沃斯低通滤波器的振幅和频率关系

传递函数

$H\left( s \right) = \frac{1}{{\mathop s\nolimits^N + \mathop b\nolimits_{N - 1} \mathop s\nolimits^{N - 1} + \cdots + \mathop b\nolimits_1 s + \mathop b\nolimits_0 }}\$

4.时域、频域、z域下的系统描述

时域

z域

频域

数字滤波器

特性描述方式

单位冲激响应

h(n)

系统函数

H(z)

频率响应

$H(\mathop e\nolimits^{jw} )\$

输出输入关系

y(n)=x(n)*h(n)

Y(z)=X(z)·H(z)

$Y(\mathop e\nolimits^{jw} ) = X(\mathop e\nolimits^{jw} ) \cdot H(\mathop e\nolimits^{jw} )\$

二、具体设计步骤

1.数字滤波器技术指标

先要已知以下指标：

$\mathop w\nolimits_p \$ ：数字通带截止频率

$\mathop w\nolimits_s \$ ：数字阻带截止频率

$\mathop \delta \nolimits_p \$ ：通带最大衰减（单位：dB）

$\mathop \delta \nolimits_s \$ ：阻带最小衰减（单位：dB）

衰减 $\mathop \delta \nolimits_p \$ 、 $\mathop \delta \nolimits_s \$ 的推导：

$\mathop \delta \nolimits_p = 20\lg \frac{{\left| {H\left( {\mathop e\nolimits^{j0} } \right)} \right|}}{{\left| {H\left( {\mathop e\nolimits^{j\mathop w\nolimits_p } } \right)} \right|}} = - 20\lg \left| {H\left( {\mathop e\nolimits^{j\mathop w\nolimits_p } } \right)} \right| = - 20\lg \left( {1 - \mathop \alpha \nolimits_1 } \right)\$

$\mathop \delta \nolimits_s = 20\lg \frac{{\left| {H\left( {\mathop e\nolimits^{j0} } \right)} \right|}}{{\left| {H\left( {\mathop e\nolimits^{j\mathop w\nolimits_s } } \right)} \right|}} = - 20\lg \left| {H\left( {\mathop e\nolimits^{j\mathop w\nolimits_s } } \right)} \right| = - 20\lg \left( {1 - \mathop \alpha \nolimits_2 } \right)\$

2.模拟滤波器指标

【需要将数字指标转化为模拟指标：冲激响应不变法、双线性变换法的转换公式不同】

$\mathop \Omega \nolimits_p \$ ：模拟通带截止频率

$\mathop \Omega \nolimits_s \$ ：模拟阻带截止频率

$\mathop \delta \nolimits_p \$ ：通带最大衰减（单位：dB）

$\mathop \delta \nolimits_s \$ ：阻带最小衰减（单位：dB）

2.1冲激响应不变法

$\mathop \Omega \nolimits_p = \frac{{\mathop w\nolimits_p }}{T}\$

$\mathop \Omega \nolimits_s = \frac{{\mathop w\nolimits_s }}{T}\$

$\mathop \delta \nolimits_p \$ 、 $\mathop \delta \nolimits_s \$ 不变

2.2双线性变换法

【需要预畸变】

$\mathop \Omega \nolimits_p = c \cdot \tan \left( {\frac{{\mathop w\nolimits_p }}{2}} \right)\$

$\mathop \Omega \nolimits_s = c \cdot \tan \left( {\frac{{\mathop w\nolimits_s }}{2}} \right)\$

（c为任意常数，其中低通滤波器c取 $\frac{2}{T}\$ ）

$\mathop \delta \nolimits_p \$ 、 $\mathop \delta \nolimits_s \$ 不变

3.过渡模拟滤波器

需要计算：

N：阶数

$\mathop \Omega \nolimits_{\rm{c}} \$ ：3dB截止频率

3.1求阶数N

（1）计算 $\mathop k\nolimits_{sp} \$

$\mathop k\nolimits_{sp} = \sqrt {\frac{{\mathop {10}\nolimits^{0.1\mathop \delta \nolimits_s } - 1}}{{\mathop {10}\nolimits^{0.1\mathop \delta \nolimits_p } - 1}}} \$

（2）计算归一化频率 $\mathop \lambda \nolimits_{sp} \$

$\mathop \lambda \nolimits_{sp} = \frac{{\mathop \Omega \nolimits_s }}{{\mathop \Omega \nolimits_p }} = \frac{{\mathop f\nolimits_s }}{{\mathop f\nolimits_p }}\$

（3）得到阶数N

$N = \left\lceil {\frac{{\lg \mathop k\nolimits_{sp} }}{{\lg \mathop \lambda \nolimits_{sp} }}} \right\rceil \$ （向上取整）

3.2 查表确定归一化低通原型系统函数 $\mathop H\nolimits_{an} \left( s \right)\$

3.3求3dB截止频率 $\mathop \Omega \nolimits_{\rm{c}} \$

【冲激响应不变法、双线性变换法公式不同】

冲激响应不变法：应用通带指标 $\mathop \Omega \nolimits_p \$

$\mathop \Omega \nolimits_c = \mathop \Omega \nolimits_p \mathop {\left( {\mathop {10}\nolimits^{0.1\mathop \delta \nolimits_p } - 1} \right)}\nolimits^{ - \frac{1}{{2N}}} \$

（原因：冲激响应不变法需要阻带效果好）

双线性变换法：应用阻带指标 $\mathop \Omega \nolimits_s \$

$\mathop \Omega \nolimits_c = \mathop \Omega \nolimits_s \mathop {\left( {\mathop {10}\nolimits^{0.1\mathop \delta \nolimits_s } - 1} \right)}\nolimits^{ - \frac{1}{{2N}}} \$

3.4去归一化得到模拟系统函数 $\mathop H\nolimits_a (s)\$

$\mathop H\nolimits_a (s) = \mathop H\nolimits_{an} (\frac{s}{{\mathop \Omega \nolimits_c }})\$

4.数字化：得到数字滤波器系统函数H(z)

4.1冲激响应不变法

先将模拟系统函数 $\mathop H\nolimits_a (s)\$ 展开成部分分式形式：

$\mathop H\nolimits_a (s) = \sum\limits_{k = 1}^N {\frac{{\mathop A\nolimits_k }}{{s - \mathop s\nolimits_k }}} \$

再变换成巴特沃斯数字滤波器：

$H(z) = \sum\limits_{k = 1}^N {\frac{{T\mathop A\nolimits_k }}{{1 - \mathop e\nolimits^{\mathop s\nolimits_k T} \mathop z\nolimits^{ - 1} }}} \$

4.2双线性变换法

$H(z) = \mathop {\mathop H\nolimits_a (s)}\nolimits_{|s = \frac{2}{T} \cdot \frac{{1 - \mathop z\nolimits^{ - 1} }}{{1 + \mathop z\nolimits^{ - 1} }}} \$

三、Matlab代码实现

任务要求：

1.冲激响应不变法法

butttord函数

功能：查找巴特沃思模拟滤波器的最小阶数n和截止频率Wc
格式：
```
[N,wc] = buttord(wp,ws,ap,as,’s’)
```
参数：

wp:通带截止频率 ws ：阻带截止频率

ap:通带最大允许衰减 as:阻带最小允许衰减

‘s’：表示模拟滤波器（省略s表示数字滤波器）

N: 滤波器阶数 wc: 3db截止频率

buttap 函数

功能：可设计出N阶巴特沃斯低通滤波器的零、极点
格式：
```
[z，p，k] = buttap(N)
```
参数：

z：N阶巴特沃斯低通滤波器的零点

p：N阶巴特沃斯低通滤波器的极点

k：增益

N：阶数

zp2tf 函数

功能：由零极点和增益确定归一化Han(s)系数
格式：
```
[Bap,Aap]=zp2tf(z,p,k))
```
参数：

z：N阶巴特沃斯低通滤波器的零点

p：N阶巴特沃斯低通滤波器的极点

Bap：归一化系统函数Han(s)分子

Aap：归一化系统函数Han(s)分母

lp2lp 函数（低通使用）

功能：用s/wc 代替s，去归一化Han(s)，得到H(s)
格式：
```
[Bbs,Abs]=lp2lp(Bap,Aap,wac)
```
参数：

Bap：归一化系统函数Han(s)分子

Aap：归一化系统函数Han(s)分母

wac：3dB截止频率

Bbs：数字域系统函数分子

Abs：数字域系统函数分母

impinvar函数（冲激响应不变法）

功能：使用冲击响应不变法进行数字化
格式：
```
[b,a]=impinvar(B,A,fs)
```
参数：

B,A：模拟滤波器系统函数分子分母多项式系数向量

b,a：数字滤波器系统函数分子分母多项式系数向量

fs：采样频率

freqz函数

功能：求频率响应函数
格式：
```
[h,w] = freqz(b,a,n,'whole')
```
参数：

输入：

b,a：数字滤波器系统函数分子分母多项式系数向量

n：计算点数（默认为512）

'whole'：w范围为0-2π

输出：（若默认n，h和w大小都是512）

h：频率响应向量，为幅度【一般单位：dB，用公式20 * log10(abs(h))来转换】

w：数字角频率【画图时要转换为模拟频率f】

代码总览

clear;close all;clc
%1.数字滤波器的技术指标
ap = 1;%通带最大衰减
as = 30;%阻带最小衰减
fp = 200;%通带截止频率
fs = 400;%阻带截止频率
Fs = 1000;%抽样间隔
T = 1/Fs;%2.将数字指标转化成模拟滤波器技术指标（采用冲激响应不变法）
wap = 2*pi*fp;
was = 2*pi*fs;%3.设计滤波器
[N,wac] = buttord(wap,was,ap,as,'s');%计算阶数N和3dB截止频率wac
[z,p,k] = buttap(N);%创建巴特沃斯低通滤波器 z零点p极点k增益
[Bap,Aap] = zp2tf(z,p,k);%由零极点和增益确定归一化Han(s)系数
[Bbs,Abs] = lp2lp(Bap,Aap,wac);%将s/wc 代替s，去归一化
[B,A] = impinvar(Bbs,Abs,Fs);%模拟域到数字域——冲激响应不变法
[H1,w] = freqz(B,A);%根据H(z)求频率响应%绘制数字滤波器频响幅度谱
figure(1);
f = w*Fs/(2*pi);
subplot(2,1,1);
plot(f,20*log10(abs(H1)));%绘制幅度响应
title('冲激响应不变法——巴特沃斯BLPF(幅度)');
xlabel('频率/Hz');
ylabel('H1幅值/dB');
subplot(2,1,2);
plot(f,unwrap(angle(H1)));% 绘制相位响应
xlabel('频率/Hz');
ylabel('角度/Rad');
title('冲激响应不变法——巴特沃斯BLPF(相位)');

结果展示

2.双线性变换法

【除了数字化用到的函数不同，其它都相同】

bilinear函数

功能：使用双线性变换法进行数字化
格式：
```
[b,a]=bilinear(B,A,fs)
```
参数：

B,A：模拟滤波器系统函数分子分母多项式系数向量

b,a：数字滤波器系统函数分子分母多项式系数向量

fs：采样频率

代码总览

clear,clc,close;
%1.数字滤波器的技术指标要求
ap = 1;%通带最大衰减
as = 30;%阻带最小衰减
fp = 200;%通带截止频率
fs = 400;%阻带截止频率
Fs = 1000;%抽样间隔
T = 1/Fs;%2.将数字指标转化成模拟滤波器指标
wp=(2*pi*fp)/Fs;
ws=(2*pi*fs)/Fs;
% 数字指标转模拟指标 预畸变,前面要× (2/T)
wap=2*Fs*tan(wp/2);
was=2*Fs*tan(ws/2);%3.设计模拟滤波器
[N,wac]=buttord(wap,was,ap,as,'s');% N为阶数，wac为3dB截止频率
[z,p,k]=buttap(N);% 创建巴特沃斯低通滤波器 z零点p极点k增益
[Bap,Aap]=zp2tf(z,p,k);% 由零极点和增益确定归一化Han(s)系数
[Bbs,Abs]=lp2lp(Bap,Aap,wac);% 低通到低通 计算去归一化Ha(s)
[B,A] = bilinear(Bbs,Abs,Fs); % 模拟域到数字域:双线性不变法
[H1,w] = freqz(B,A);% 根据H(z)求频率响应特性%绘制数字滤波器频响幅度谱
figure(2);
f=w*Fs/(2*pi);
subplot(211);
plot(f,20*log10(abs(H1))); % 绘制幅度响应
title('双线性变换法——巴特沃斯BLPF(幅度)');
xlabel('频率/Hz');
ylabel('H1幅值/dB');
subplot(212);
plot(f,unwrap(angle(H1)));% 绘制相位响应
xlabel('频率/Hz');
ylabel('角度/Rad')
title('双线性变换法——巴特沃斯BLPF(相位)');

结果展示

数字信号处理|Matlab设计巴特沃斯低通滤波器（冲激响应不变法和双线性变换法）相关推荐

巴特沃斯低通滤波器归一化matlab,基于MATLAB设计巴特沃斯低通滤波器
<基于MATLAB设计巴特沃斯低通滤波器>由会员分享,可在线阅读,更多相关<基于MATLAB设计巴特沃斯低通滤波器(6页珍藏版)>请在人人文库网上搜索. 1.基于MATLAB设 ...
MATLAB：巴特沃斯低通滤波器过滤信号
MATLAB:巴特沃斯低通滤波器过滤信号实验内容产生两个不同频率的正弦信号,设计合适的滤波器,输出其中的一个信号,滤除另外的一个信号. 代码: %生成两个不同频率的信号 clc; T=2;%时域长度 ...
用matlab做巴特沃斯低通滤波器
用matlab做巴特沃斯低通滤波器趁着暑假,做一个心电图的matlab实验,遇到了滤波器问题,网上代码比较杂乱,做了一个汇总整理. 主要做了一个简单的低通滤波器并以三角函数为例子进行低通滤波. 基本 ...
【信号处理】数字信号处理Matlab设计附GUI界面和报告
1 简介信号处理是现代信息处理的基本内容,数字信号的处理更是重中之重.数字信号处理的研究内容主要是语音信号和图像信号,而研究语音信号对于现代语音通信有着积极的意义.研究语音信号又分为时域和频域两个 ...
数字信号处理matlab设计滤波器
目录 1 设计目的.内容及要求 1 1.1 设计目的 1 1.2 设计内容 1 1.3 设计思考 2 1.4 设计要求 2 2 设计方案与设计原理 3 2.1 设计思路 3 2.2 采样定理 3 ...
C语言实现巴特沃斯低通滤波器
本文C语言的实现公式和思路参考自: 数字信号处理|Matlab设计巴特沃斯低通滤波器(冲激响应不变法和双线性变换法) 还有北京交通大学陈后金教授和华东理工大学万永菁教授的网课文章目录 ma ...
四旋翼惯导融合之观测传感器滞后问题汇总与巴特沃斯低通滤波器设计（气压计MS5611、GPS模块M8N、超声波、PX4FLOW等）
前文讲到APM的三阶互补方案,之前附的图是从学长博客里面抠的,感觉还不是很详细,于是自己就画了下,顺便重新理一下思路. 上图中下标为O的表示原始量(Origion),C表示矫正后的量(Correcti ...
使用matlab设计IIR巴特沃斯低通滤波器
目的和要求设计IIR巴特沃斯低通滤波器将滤波器用于加噪声信号的处理改变参数指标比较不同实验结果与分析 (1)滤波器频率特性如图输入信号特性如图: 经过滤波器之后信号特性如图: (2)改变指标 ...
MATLAB_数字信号处理_模拟滤波器_设计巴特沃斯滤波器
简介巴特沃斯滤波器-百度百科巴特沃斯滤波器-维基百科简介:巴特沃斯滤波器是一种模拟滤波器,它在频率响应方面具有特殊的属性.它被设计为具有均匀的幅度响应,即在通带内,它对所有频率的增益是相等的,而 ...

数字信号处理|Matlab设计巴特沃斯低通滤波器（冲激响应不变法和双线性变换法）

一、前言

1. 设计流程

2.系统频率响应

2.1频响图

2.2 各个频率转换关系

3.巴特沃斯滤波器简介

3.1特点

3.2设计思路

4.时域、频域、z域下的系统描述

二、具体设计步骤

1.数字滤波器技术指标

2.模拟滤波器指标

2.1冲激响应不变法

2.2双线性变换法

3.过渡模拟滤波器

3.1求阶数N

3.2 查表确定归一化低通原型系统函数 $\mathop H\nolimits_{an} \left( s \right)\$

3.3求3dB截止频率 $\mathop \Omega \nolimits_{\rm{c}} \$

3.4去归一化得到模拟系统函数 $\mathop H\nolimits_a (s)\$

4.数字化：得到数字滤波器系统函数H(z)

4.1冲激响应不变法

4.2双线性变换法

三、Matlab代码实现

1.冲激响应不变法法

butttord函数

buttap 函数

zp2tf 函数

lp2lp 函数（低通使用）

impinvar函数（冲激响应不变法）

freqz函数

代码总览

结果展示

2.双线性变换法

bilinear函数

代码总览

结果展示

数字信号处理|Matlab设计巴特沃斯低通滤波器（冲激响应不变法和双线性变换法）相关推荐

最新文章

热门文章